In the most recent two decades, the development of impurity controls with ultra-high precision in semiconductors motivates people to put more and more attentions on the solotronic devices, whose properties depend on one or a few dopants. One of the most promising applications of solotronic device is the qubit in quantum computing. In the procedure of exploring qubit candidates, the most straightforward challenges we need face include that the qubit must be highly isolated and can be initialized/manipulated efficiently with high fidelities. It has been proved that qubits based on single defects have excellent performances as quits. For instance, the NV center in diamond forms a ground spin triplet which can be manipulated at room temperature with electromagnetic fields. This work focuses on searching for new single defects as qubit candidates with density functional theory. Lanthanides element possesses excellent optical characteristics and extremely long nuclear coherence time. Therefore, combining it into the diamond platform can be possible design for integrated quantum information processing devices in the future. To investigate the stability of lanthanides dopants in the diamond matrix, the formation energies of charge states of complexes are calculated. The broadening of Eu(III) peak in the photoluminescence spectrum can be verified according to the existence of more than stable configuration and steady 4f electron occupation. In the case of transition-metal dopant in the silicon carbide, it is found that both silicon and carbon substituted nickels in 3C-SiC shows a magnetic-antimagnetic transition under applied strains. The virtual hopping rate of electrons strongly depends on the distance between the spin pair residing in the nickel and dangling bonds. Therefore, the Heisenberg exchange coupling between them can be adjusted subtly by controlling the external strain. According to the the spin Hamiltonian of the defect, the spin state can be manipulated universally with strain and electromagnetic fields. In contrast, the nickel dopant in 4H-SiC exhibits a very stable magnetic property. Other than that, the electronic structure of Cr in 4H-SiC implies that optical manipulations of spin states might be realized in the excited state.

The problem of spin manipulation via the spin-orbit interaction in nonmagnetic semiconductors in the absence of magnetic fields is investigated in this work. We begin with a review of the literature on spin dynamics in semiconductors, then discuss the semi-empirical k ⋅ p method of calculating direct-gap semiconductor properties, which we use to estimate material parameters significant for manipulation of spin even in the absence of a magnetic field. The total effective magnetic fields and precession lengths are calculated for a variety of quantum well orientations, and a class of devices are proposed that will allow for all-electric arbitrary manipulation of spin orientations. The strain- and momentum-dependent spin splitting coefficient C3 has been calculated using a fourteen band Kane k⋅p model for a variety of III-V semiconductors as well as ZnSe and CdSe. It is observed that the spin-splitting parameters C3 and γ, corresponding to the strain-induced spin-orbit interaction and Dresselhaus coefficient, are sensitive to the value of the inter-band spin-orbit coupling Δ− between the p valence and p̄ second conduction band in all cases. The value of Δ− has therefore been recalculated in these materials using a tight-binding model and modern experimental values of the valence and second conduction band spin-orbit splittings. The total effective magnetic field and precession length of spins in strained quantum wells in the (001), (110), and (111) planes are derived with consideration for all known effective magnetic fields except those due to interface effects in non- common-atom heterostructures (native inversion asymmetry). The orientation of the k-linear Dresselhaus field and the strain-dependent fields vary strongly with the growth axis of the quantum well. The precession length in the (110) and (001) cases can achieve infinite anisotropy, while the precession length of (111) quantum wells is always isotropic. We find that the electronic spin rotation induced by drift transport around a closed path in a wide variety of nonmagnetic semiconductors at zero magnetic field depends solely on the physical path taken. Physical paths that produce any possible spin rotation due to transport around a closed path are constructed for electrons experiencing strain or electric fields in (001), (110), or (111)-grown zinc blende semiconductor quantum wells. Spin decoherence due to travel along the path is negligible compared to the background spin decoherence rate. The small size of the designed paths (< 100 nm scale in GaAs) may lead to applications in nanoscale spintronic circuits.

The exact and phaseless variants of auxiliary-field quantum Monte Carlo (AFQMC) have been shown to be capable of producing accurate ground-state energies for a wide variety of systems including those which exhibit substantial electron correlation effects. The first chapter of this thesis will provide an overview of the relevant electronic structure problem, and the phaseless AFQMC (ph-AFQMC) methodology. The computational cost of performing these calculations has to date been relatively high, impeding many important applications of these approaches. In Chapter 2 we present a correlated sampling methodology for AFQMC which relies on error cancellation to dramatically accelerate the calculation of energy differences of relevance to chemical transformations. In particular, we show that our correlated sampling-based ph-AFQMC approach is capable of calculating redox properties, deprotonation free energies, and hydrogen abstraction energies in an efficient manner without sacrificing accuracy. We validate the computational protocol by calculating the ionization potentials and electron affinities of the atoms contained in the G2 test set and then proceed to utilize a composite method, which treats fixed-geometry processes with correlated sampling-based AFQMC and relaxation energies via MP2, to compute the ionization potential, deprotonation free energy, and the O-H bond dissociation energy of methanol, all to within chemical accuracy. We show that the efficiency of correlated sampling relative to uncorrelated calculations increases with system and basis set size and that correlated sampling greatly reduces the required number of random walkers to achieve a target statistical error. This translates to reductions in wall-times by factors of 55, 25, and 24 for the ionization potential of the K atom, the deprotonation of methanol, and hydrogen abstraction from the O-H bond of methanol, respectively. In Chapter 3 we present an implementation of ph-AFQMC utilizing graphical processing units (GPUs). The AFQMC method is recast in terms of matrix operations which are spread across thousands of processing cores and are executed in batches using custom Compute Unified Device Architecture kernels and the hardware-optimized cuBLAS matrix library. Algorithmic advances include a batched Sherman-Morrison-Woodbury algorithm to quickly update matrix determinants and inverses, density-fitting of the two-electron integrals, an energy algorithm involving a high-dimensional precomputed tensor, and the use of single-precision floating point arithmetic. These strategies result in dramatic reductions in wall-times for both single- and multi-determinant trial wavefunctions. For typical calculations we find speed-ups of roughly two orders of magnitude using just a single GPU card. Furthermore, we achieve near-unity parallel efficiency using 8 GPU cards on a single node, and can reach moderate system sizes via a local memory-slicing approach. We illustrate the robustness of our implementation on hydrogen chains of increasing length, and through the calculation of all-electron ionization potentials of the first-row transition metal atoms. We compare long imaginary-time calculations utilizing a population control algorithm with our previously published correlated sampling approach, and show that the latter improves not only the efficiency but also the accuracy of the computed ionization potentials. Taken together, the GPU implementation combined with correlated sampling provides a compelling computational method that will broaden the application of ph-AFQMC to the description of realistic correlated electronic systems. In Chapter 4 the bond dissociation energies of a set of 44 3d transition metal-containing diatomics are computed with ph-AFQMC utilizing the correlated sampling technique. We investigate molecules with H, N, O, F, Cl, and S ligands, including those in the 3dMLBE20 database first compiled by Truhlar and co-workers with calculated and experimental values that have since been revised by various groups. In order to make a direct comparison of the accuracy of our ph-AFQMC calculations with previously published results from 10 DFT functionals, CCSD(T), and icMR-CCSD(T), we establish an objective selection protocol which utilizes the most recent experimental results except for a few cases with well-specified discrepancies. With the remaining set of 41 molecules, we find that ph-AFQMC gives robust agreement with experiment superior to that of all other methods, with a mean absolute error (MAE) of 1.4(4) kcal/mol and maximum error of 3(3) kcal/mol (parenthesis account for reported experimental uncertainties and the statistical errors of our ph-AFQMC calculations). In comparison, CCSD(T) and B97, the best performing DFT functional considered here, have MAEs of 2.8 and 3.7 kcal/mol, respectively, and maximum errors in excess of 17 kcal/mol (for the CoS diatomic). While a larger and more diverse data set would be required to demonstrate that ph-AFQMC is truly a benchmark method for transition metal systems, our results indicate that the method has tremendous potential, exhibiting unprecedented consistency and accuracy compared to other approximate quantum chemical approaches. The energy gap between the lowest-lying singlet and triplet states is an important quantity in chemical photocatalysis, with relevant applications ranging from triplet fusion in optical upconversion to the design of organic light-emitting devices. The ab initio prediction of singlet-triplet (ST) gaps is challenging due to the potentially biradical nature of the involved states, combined with the potentially large size of relevant molecules. In Chapter 5, we show that ph-AFQMC can accurately predict ST gaps for chemical systems with singlet states of highly biradical nature, including a set of 13 small molecules and the ortho-, meta-, and para- isomers of benzyne. With respect to gas-phase experiments, ph-AFQMC using CASSCF trial wavefunctions achieves a mean averaged error of ~1 kcal/mol. Furthermore, we find that in the context of a spin-projection technique, ph-AFQMC using unrestricted single-determinant trial wavefunctions, which can be readily obtained for even very large systems, produces equivalently high accuracy. We proceed to show that this scalable methodology is capable of yielding accurate ST gaps for all linear polyacenes for which experimental measurements exist, i.e. naphthalene, anthracene, tetracene, and pentacene. Our results suggest a protocol for selecting either unrestricted Hartree-Fock or Kohn-Sham orbitals for the single-determinant trial wavefunction, based on the extent of spin-contamination. These findings provide a reliable computational tool with which to investigate specific photochemical processes involving large molecules that may have substantial biradical character. We compute the ST gaps for a set of anthracene derivatives which are potential triplet-triplet annihilators for optical upconversion, and compare our ph-AFQMC predictions with those from DFT and CCSD(T) methods. We conclude with a discussion of ongoing projects, further methodological improvements on the horizon, and future applications of ph-AFQMC to chemical systems of interest in the fields of biology, drug-discovery, catalysis, and condensed matter physics.

<p>A parameter free approach for the calculation of the crystal field splitting of the lowest Russel-Saunders <i>J</i>-multiplet in <i>f</i>-electron systems has been developed and applied to selected compounds. The developed theory is applicable to general symmetries and is based on symmetry constrained density functional theory calculations in the local density or in the generalised gradient approximation.</p><p>The magnetocrystalline anisotropy of Gd has been analysed. It has been shown that the peculiar orientation of the easy axis of magnetisation is consistent with an <i>S</i>-ground state. Further, the temperature dependence of the easy axis of magnetisation has been investigated and it has been shown that the temperature driven reduction of the effective magnetisation is the principal mechanism responsible for it.</p><p>A new method has been developed that allows for theoretical studies of the electronic structure and total energy of elements and compounds in an intermediate valence regime. The method combines model and first-principles band structure calculations, therefore being accurate and computationally efficient. It has been applied to Yb metal under pressure obtaining a remarkable agreement with experimental observations for the equation of state and the x-ray absorption spectroscopy.</p>

The structural, electronic, and superconducting properties of one dimensional materials are calculated from first principles, using the density functional theory. Nanotubes and nanowires are important building blocks in nanotechnology, in particular for nanoelectronics. In this manuscript, the growth of carbon nanotubes is studied through the interaction between carbon and the transition metal atoms used as growth catalysts. The accepted model for a new phase of nanotube-like molybdenum disulfide is critically examined using comparisons of energetic stability and types of chemical bonding in different candidate structures which have similar compositions. The epitaxial growth of diamond carbon on (100) iridium is exceptionally favorable. The differences between various substrates used experimentally are studied, and the specificity of Ir is shown. Finally, the characteristics of the electron-phonon interaction in aluminium nanowires are determined. The structural instabilities and the differences in the electron-phonon coupling are calculated for straight monoatomic wires, zigzag wires, and thicker straight wires. The constrained geometry of the wires generates a coupling which can be very strong or almost vanish, depending on the structural details, but which is concentrated in the longitudinal high-frequency phonons.

This Ph.D. dissertation is focused on the investigation the structure of nitrile anion containing molecules and how the structure and reactivity of those molecules are affected by solvation and counter ion. A systematic approach was employed in this investigation, beginning with an evaluation of the accuracy of three commonly used model chemistries (Hartree-Fock (HF), Second-order Møller-Plesset perturbation theory (MP2), the Becke three-parameter exchange functional coupled with the nonlocal correlation functional of Lee, Yang, and Parr (B3LYP), all paired with the 6-31+G(d) basis set). A series of complexes of various cations with a number of explicit molecules of tetrahydrofuran (THF) and dimethyl ether (DME) were studied with these model chemistries and the results were compared, where possible, with experimental results. From this work, it was determined that the B3LYP models gave the most accurate results for the complexes in question. This work was then extended to acetonitrile anion containing complexes of solvent and cation. Based on the results of that extension, it was determined that cation size and charge density on the cation were critical factors in determining the structure of the acetonitrile anion molecule and in determining if the anion was metalated at the nitrogen or alpha-carbon position, with larger cations favoring carbon metalation and more significant deformation of the alpha-carbon from the expected sp2 hybridization. The final aspect of this dissertation was the determination of reaction coordinate energy profiles for a pair of substitution reactions involving nitrile anion containing cycloaliphatic molecules. The results of this study showed that, due to steric and kinetic factors, the axial products and transitions states associated with these reactions were favored, and that the degree of preference was kinetically controlled. / Bayer School of Natural and Environmental Sciences / Chemistry and Biochemistry / PhD / Dissertation

The structural, electronic, and superconducting properties of one dimensional materials are calculated from first principles, using the density functional theory. Nanotubes and nanowires are important building blocks in nanotechnology, in particular for nanoelectronics. In this manuscript, the growth of carbon nanotubes is studied through the interaction between carbon and the transition metal atoms used as growth catalysts. The accepted model for a new phase of nanotube-like molybdenum disulfide is critically examined using comparisons of energetic stability and types of chemical bonding in different candidate structures which have similar compositions. The epitaxial growth of diamond carbon on (100) iridium is exceptionally favorable. The differences between various substrates used experimentally are studied, and the specificity of Ir is shown. Finally, the characteristics of the electron-phonon interaction in aluminium nanowires are determined. The structural instabilities and the differences in the electron-phonon coupling are calculated for straight monoatomic wires, zigzag wires, and thicker straight wires. The constrained geometry of the wires generates a coupling which can be very strong or almost vanish, depending on the structural details, but which is concentrated in the longitudinal high-frequency phonons.

Clusters are aggregates of a finite number of atoms or molecules. In the present work, free clusters out of metallic parent materials have been created and studied by synchrotron radiation-based photoelectron spectroscopy. The clusters have been formed and studied in a beam and the electronic structure of the clusters has been investigated. Conclusions have been drawn about the spatial distribution of atoms of different elements in bi-component clusters, about the development of metallicity in small clusters, and about the excitation of plasmons. Bi-component alloy clusters of sodium and potassium and of copper and silver have been produced. The site-sensitivity of the photoelectron spectroscopy technique has allowed us to probe the geometric distribution of the atoms of the constituent elements by comparing the responses from the bulk and surface of the clusters. In both cases, we have found evidence for a surface-segregated structure, with the element with the largest atoms and lowest cohesive energy (potassium and silver, correspondingly) dominating the surface and with a mixed bulk. Small clusters of tin and lead have been probed to investigate the development of metallicity. The difference in screening efficiency between metals and non-metals has been utilized to determine in what size range an aggregate of atoms of these metallic parent materials stops to be metallic. For tin this has been found to occur below ~40 atoms while for lead it happened somewhere below 20-30 atoms. The excitation of bulk and surface plasmons has been studied in clusters of sodium, potassium, magnesium and aluminium, with radii in the nanometer range. The excitation energies have been found to be close to those of the corresponding macroscopic solids. We have also observed spectral features corresponding to multi-quantum plasmon excitation in clusters of Na and K. Such features have in macroscopic solids been interpreted as due to harmonic plasmon excitation. Our observations of features corresponding to the excitation of one bulk and one surface plasmon however suggest the presence of sequential excitation in clusters.

A parameter free approach for the calculation of the crystal field splitting of the lowest Russel-Saunders J-multiplet in f-electron systems has been developed and applied to selected compounds. The developed theory is applicable to general symmetries and is based on symmetry constrained density functional theory calculations in the local density or in the generalised gradient approximation. The magnetocrystalline anisotropy of Gd has been analysed. It has been shown that the peculiar orientation of the easy axis of magnetisation is consistent with an S-ground state. Further, the temperature dependence of the easy axis of magnetisation has been investigated and it has been shown that the temperature driven reduction of the effective magnetisation is the principal mechanism responsible for it. A new method has been developed that allows for theoretical studies of the electronic structure and total energy of elements and compounds in an intermediate valence regime. The method combines model and first-principles band structure calculations, therefore being accurate and computationally efficient. It has been applied to Yb metal under pressure obtaining a remarkable agreement with experimental observations for the equation of state and the x-ray absorption spectroscopy.

This thesis is the result of investigations regarding the processes in InAs III-V semiconductor surfaces induced by additional charge incorporated by adsorbates. The aim of the project is to study the development of the accumulation layer on the metal/InAs(111)A/B surfaces and its electronic structure. InAs(111)A is indium-terminated and InAs(111)B is arsenic-terminated. In addition, InAs(100) is also studied. These three substrates are different; InAs(111)A has a (2x2)-termination, explained by an indium vacancy model, and the clean surface exhibits a two-dimensional electron gas (2DEG). InAs(111)B(1x1) is bulk-truncated and unreconstructed and does not host a 2DEG. InAs(100)(4x2)/c(8x2) exhibits a more covalent character of the surface bonds compared to InAs(111)A/B, and the surface is terminated by a complex reconstruction. Photoelectron spectroscopy and LEED (low energy electron diffraction) have been used as the main tools to study surfaces with sub-monolayer to monolayer amounts of adsorbates. A photoemission peak related to a two-dimensional electron gas appears close to the Fermi level. This 2DEG has in most cases InAs bulk properties, since it is located in the InAs conduction band. A systematic study of core levels and valence bands reveals that the appearance of the 2DEGs is a complex process connected to the surface order. Adsorption of lead, tin or bismuth on InAs(111)B(1x1) induces emission from a 2DEG, but only at monolayer coverage and when the surface is ordered. Cobalt reacts strongly with InAs forming InCo islands and no accumulation is observed. Examination of Cs/InAs(111)B does not reveal any 2DEG and the surface reaction is strongly related to the clean surface stabilization process. Examination of the In-terminated InAs(111)A(2x2) surface shows that In reacts strongly with cobalt and tin adatoms and with oxygen in cases of large exposure, which decreases the 2DEG intensity, while adatoms of cesium and small doses of oxygen enhance the emission from the 2DEG. InAs(100) is terminated with one kind of atom - the InAs(100)(4x2)/c(8x2) is indium terminated. Bismuth creates dimers on the surface and a 2DEG is observed. More generally, this thesis describes some of the general physical background applied to surface science and 2DEG. The first part is a general overview of the processes on the surface. The second part concentrates on the methods related to preparation of samples, and the third part on the measurement methods. The photoelectron spectroscopy part concerns the theory used in mapping electronic structure. The inserted figures are taken from different experiments, including results for InAs(111)A not previously published. / QC 20100910