QM/MM simulation of liquid water with an adaptive quantum region.

The simulation of complex chemical systems often requires a multi-level description, in which a region of special interest is treated using a computationally expensive quantum mechanical (QM) model while its environment is described by a faster, simpler molecular mechanical (MM) model. Furthermore, studying dynamic effects in solvated systems or bio-molecules requires a variable definition of the two regions, so that atoms or molecules can be dynamically re-assigned between the QM and MM descriptions during the course of the simulation. Such reassignments pose a problem for traditional QM/MM schemes by exacerbating the errors that stem from switching the model at the boundary. Here we show that stable, long adaptive simulations can be carried out using density functional theory with the BLYP exchange-correlation functional for the QM model and a flexible TIP3P force field for the MM model without requiring adjustments of either. Using a primary benchmark system of pure water, we investigate the convergence of the liquid structure with the size of the QM region, and demonstrate that by using a sufficiently large QM region (with radius 6 Å) it is possible to obtain radial and angular distributions that, in the QM region, match the results of fully quantum mechanical calculations with periodic boundary conditions, and, after a smooth transition, also agree with fully MM calculations in the MM region. The key ingredient is the accurate evaluation of forces in the QM subsystem which we achieve by including an extended buffer region in the QM calculations. We also show that our buffered-force QM/MM scheme is transferable by simulating the solvated Cl(-) ion.

Gaussian Approximation Potential: an interatomic potential derived from first principles Quantum Mechanics

Simulation of materials at the atomistic level is an important tool in studying microscopic structure and processes. The atomic interactions necessary for the simulation are correctly described by Quantum Mechanics. However, the computational resources required to solve the quantum mechanical equations limits the use of Quantum Mechanics at most to a few hundreds of atoms and only to a small fraction of the available configurational space. This thesis presents the results of my research on the development of a new interatomic potential generation scheme, which we refer to as Gaussian Approximation Potentials. In our framework, the quantum mechanical potential energy surface is interpolated between a set of predetermined values at different points in atomic configurational space by a non-linear, non-parametric regression method, the Gaussian Process. To perform the fitting, we represent the atomic environments by the bispectrum, which is invariant to permutations of the atoms in the neighbourhood and to global rotations. The result is a general scheme, that allows one to generate interatomic potentials based on arbitrary quantum mechanical data. We built a series of Gaussian Approximation Potentials using data obtained from Density Functional Theory and tested the capabilities of the method. We showed that our models reproduce the quantum mechanical potential energy surface remarkably well for the group IV semiconductors, iron and gallium nitride. Our potentials, while maintaining quantum mechanical accuracy, are several orders of magnitude faster than Quantum Mechanical methods.

Communication: energy benchmarking with quantum Monte Carlo for water nano-droplets and bulk liquid water.

We show the feasibility of using quantum Monte Carlo (QMC) to compute benchmark energies for configuration samples of thermal-equilibrium water clusters and the bulk liquid containing up to 64 molecules. Evidence that the accuracy of these benchmarks approaches that of basis-set converged coupled-cluster calculations is noted. We illustrate the usefulness of the benchmarks by using them to analyze the errors of the popular BLYP approximation of density functional theory (DFT). The results indicate the possibility of using QMC as a routine tool for analyzing DFT errors for non-covalent bonding in many types of condensed-phase molecular system.

Quantum Confinement Effects and Electronic Properties of SnO2 Quantum Wires and Dots

On the basis of the density functional theory (DFT) within local density approximations (LDA) approach, we calculate the band gaps for different size SnO2 quantum wire (QWs) and quantum dots (QDs). A model is proposed to passivate the surface atoms of SnO2 QWs and QDs. We find that the band gap increases between QWs and bulk evolve as Delta E-g(wire) = 1.74/d(1.20) as the effective diameter d decreases, while being Delta E-g(dot) = 2.84/d(1.26) for the QDs. Though the similar to d(1.2) scale is significantly different from similar to d(2) of the effective mass result, the ratio of band gap increases between SnO2 QWs and QDs is 0.609, very close to the effective mass prediction. We also confirm, although the LDS calculations underestimate the band gap, that they give the trend of band gap shift as much as that obtained by the hybrid functional (PBE0) with a rational mixing of 25% Fock exchange and 75% of the conventional Perdew-Burke-Ernzerhof (PBE) exchange functional for the SnO2 QWs and QDs. The relative deviation of the LDA calculated band gap difference Lambda E-g compared with the corresponding PBE0 results is only within 5%. Additionally, it is found the states of valence band maximum (VBM) and conduction band minimum (CBM) of SnO2 QWs or QDs have a mostly p- and s-like envelope function symmetry, respectively, from both LDA and PBE0 calculations.

Rashba electron transport in one-dimensional quantum waveguides

The properties of Rashba wave function in the planar one-dimensional waveguide are studied, and the following results are obtained. Due to the Rashba effect, the plane waves of electron with the energy E divide into two kinds of waves with the wave vectors k(1)=k(0)+k(delta) and k(2)=k(0)-k(delta), where k(delta) is proportional to the Rashba coefficient, and their spin orientations are +pi/2 (spin up) and -pi/2 (spin down) with respect to the circuit, respectively. If there is gate or ferromagnetic contact in the circuit, the Rashba wave function becomes standing wave form exp(+/- ik(delta)l)sin[k(0)(l-L)], where L is the position coordinate of the gate or contact. Unlike the electron without considering the spin, the phase of the Rashba plane or standing wave function depends on the direction angle theta of the circuit. The travel velocity of the Rashba waves with the wave vector k(1) or k(2) are the same hk(0)/m*. The boundary conditions of the Rashba wave functions at the intersection of circuits are given from the continuity of wave functions and the conservation of current density. Using the boundary conditions of Rashba wave functions we study the transmission and reflection probabilities of Rashba electron moving in several structures, and find the interference effects of the two Rashba waves with different wave vectors caused by ferromagnetic contact or the gate. Lastly we derive the general theory of multiple branches structure. The theory can be used to design various spin polarized devices.

Electrically controllable RKKY interaction in semiconductor quantum wires

We demonstrate in theory that it is possible to all-electrically manipulate the RKKY interaction in a quasi-one-dimensional electron gas embedded in a semiconductor heterostructure, in the presence of Rashba and Dresselhaus spin-orbit interaction. In an undoped semiconductor quantum wire where intermediate excitations are gapped, the interaction becomes the short-ranged Bloembergen-Rowland superexchange interaction. Owing to the interplay of different types of spin-orbit interaction, the interaction can be controlled to realize various spin models, e.g., isotropic and anisotropic Heisenberg-like models, Ising-like models with additional Dzyaloshinsky-Moriya terms, by tuning the external electric field and designing the crystallographic directions. Such controllable interaction forms a basis for quantum computing with localized spins and quantum matters in spin lattices.

Photoluminescence properties of tensile-strained GaAsP/GaInP single quantum wells grown by metal organic chemical vapor deposition

Tensile-strained GaAsP/GaInP single quantum well (QW) laser diode (I-D) structures have been grown by low-pressure metal organic chemical vapor deposition (LP-MOCVD) and related photoluminescence (PL) properties have been investigated in detail. The samples have the same well thickness of 16 nm but different P compositions in a GaAsP QW. Two peaks in room temperature (RT) PL spectra are observed for samples with a composition larger than 0.10. Temperature and excitation-power-dependent PL spectra have been measured for a sample with it P composition of 0.15. It is found that the two peaks have a 35 meV energy separation independent of temperature and only the low-energy peak exists below 85 K. Additionally, both peak intensities exhibit a monotonous increase as excitation power increases. Analyses indicate that the two peaks arise from the intrinsic-exciton recombination mechanisms of electron-heavy hole (e-hh) and electron-light hole (e-hh). A theoretical calculation based oil model-solid theory, taking, into account the spin-orbit splitting energy, shows good agreement with our experimental results. The temperature dependence of PL intensity ratio is well explained using the spontaneous emission theory for e-hh and e-hh transitions. front which the ratio can be characterized mainly by the energy separation between the fill and Ill states.

Temperature-dependent transport and spin accumulation in a quantum wire with Rashba spin-orbit interaction

We study the theory of temperature-dependent electron transport, spin polarization, and spin accumulation in a Rashba spin-orbit interaction (RSOI) quantum wire connected nonadiabatically to two normal conductor electrode leads. The influence of both the wire-lead connection and the RSOI on the electron transport is treated analytically by means of a scattering matrix technique and by using an effective free-electron approximation. Through analytical analysis and numerical examples, we demonstrate a simple way to design a sensitive spin-transfer switch that operates without applying any external magnetic fields or attaching ferromagnetic contacts. We also demonstrate that the antisymmetry of the spin accumulation can be destroyed slightly by the coupling between the leads and the wire. Moreover, temperature can weaken the polarization and smear out the oscillations in the spin accumulation.

Electronic states of InAs/GaAs tyre-shape quantum ring

In the framework of the effective mass theory, this paper calculates the electron energy levels of an InAs/GaAs tyre-shape quantum ring (TSQR) by using the plane wave basis. The results show that the electron energy levels are sensitively dependent on the TSQR's section thickness d, and insensitively dependent on TSQR's section inner radius R-1 and TSQR's inner radius R-2. The model and results provide useful information for the design and fabrication of InAs/GaAs TSQRs.

Electronic energy levels in an asymmetric quantum-dots-in-a-well structure for infrared photodetectors

Theoretical calculation of electronic energy levels of an asymmetric InAs/InGaAS/GaAS quantum-dots-in-a-well (DWELL) structure for infrared photodetectors is performed in the framework of effective-mass envelope-function theory. Our calculated results show that the electronic energy levels in quantum dots (QDs) increase when the asymmetry increases and the ground state energy increases faster than the excited state energies. Furthermore, the results also show that the electronic energy levels in QDs decrease as the size of QDs and the width of quantum well (QW) in the asymmetric DWELL structure increase. Additionally, the effects of asymmetry, the size of QDs and the width of QW on the response peak of asymmetry DWELL photodetectors are also discussed.

Confined LO-phonon-assisted magnetotunneling in a parabolic quantum well with double barriers

The confined longitudinal-optical phonon-assisted tunneling through a parabolic quantum well with double barriers in a magnetic field perpendicular to the interfaces is studied theoretically based on a dielectric continuum model. The numerical results show that the applied magnetic field sharpens and heightens the phonon-assisted tunneling peaks in agreement with experimental observation. Furthermore, the phonon-assisted magnetotunneling peaks shift towards the higher biases as the magnetic field increases. In contrast to the results for a rectangular quantum well, the ratio of peak to valley of the phonon-assisted tunneling is larger for the wider well case. It also indicates that the phonon-assisted tunneling current peaks can be easily observed for a wider parabolic quantum well. (C) 2008 Published by Elsevier B.V.

Transmission properties of electron in quantum rings

We investigated the transmission probability of a single electron transmission through a quantum ring device based on the single-band effective mass approximation method and transfer matrix theory. The time-dependent Schrodinger equation is applied on a Gaussian wave packet passing through the quantum ring system. The electron tunneling resonance peaks split when the electron transmits through a double quantum ring. The splitting energy increases as the distance between the two quantum rings decreases. We studied the tunneling time through the single electron transmission quantum ring from the temporal evolution of the Gaussian wave packet. The electron probability density is sensitive to the thickness of the barrier between the two quantum rings. (C) 2008 American Institute of Physics.

Theory of huge tunneling magnetoresistance in graphene

We investigate theoretically the spin-independent tunneling magnetoresistance effect in a graphene monolayer modulated by two parallel ferromagnets deposited on a dielectric layer. For the parallel magnetization configuration, Klein tunneling can be observed in the transmission spectrum but at specific oblique incident angles. For the antiparallel magnetization configuration, the transmission can be blocked by the magneticelectric barrier provided by the ferromagnets. Such a transmission discrepancy results in a tremendous magnetoresistance ratio and can be tuned by the inclusion of an electric barrier.

Evolution of Ge and SiGe quantum dots under excimer laser annealing

We present different relaxation mechanisms of Ge and SiGe quantum dots under excimer laser annealing. Investigation of the coarsening and relaxation of the dots shows that the strain in Ge dots on Ge films is relaxed by dislocation since there is no interface between the Ge dots and the Ge layer, while the SiGe dots on Si0.77Ge0.23 film relax by lattice distortion to coherent clots, which results from the obvious interface between the SiGe clots and the Si0.77Ge0.23 film. The results are suggested and sustained by Vanderbilt and Wickham's theory, and also demonstrate that no bulk diffusion oGeurs during the excimer laser annealing.

Spin-orbit splitting of a hydrogenic donor impurity in GaAs/GaAlAs quantum wells

The electronic states of a hydrogenic donor impurity in GaAs/GaAlAs quantum wells are investigated theoretically in the framework of effective-mass envelope function theory, including the effect of Rashba spin-orbit coupling. The splits of electron energy levels are calculated. The results show that (1) the split energy of the excited state is larger than that of the ground state; (2) the split energy peak appears as the GaAs well width increases from zero; and (3) the maximum split energy reaches about 1.6 meV. Our results are useful for the application of Rashba spin-orbit coupling to photoelectric devices. (c) 2008 American Institute of Physics.