Effects of magnetic field on the electronic structure of wurtzite quantum dots: Calculations using effective-mass envelope function theory

The Hamiltonian of the wurtzite quantum dots in the presence of an external homogeneous magnetic field is given. The electronic structure and optical properties are studied in the framework of effective-mass envelope function theory. The energy levels have new characteristics, such as parabolic property, antisymmtric splitting, and so on, different from the Zeeman splitting. With the crystal field splitting energy Delta(c)=25 meV, the dark excitons appear when the radius is smaller than 25.85 A in the absence of external magnetic field. This result is more consistent with the experimental results reported by Efros [Phys. Rev. B 54, 4843 (1996)]. It is found that dark excitons become bright under appropriate magnetic field depending on the radius of dots. The circular polarization factors of the optical transitions of randomly oriented dots are zero in the absence of external magnetic field and increase with the increase of magnetic field, in agreement with the experimental results. The circular polarization factors of single dots change from nearly 0 to about 1 as the orientation of the magnetic field changes from the x axis of the crystal structure to the z axis, which can be used to determine the orientation of the z axis of the crystal structure of individual dots. The antisymmetric Hamiltonian is very important to the effects of magnetic field on the circular polarization of the optical transition of quantum dots.

Giant Zeeman splitting in manganese-doped ZnSe quantum spheres: Eight-band effective-mass calculations

We deduce the eight-band effective-mass Hamiltonian model for a manganese-doped ZnSe quantum sphere in the presence of the magnetic field, including the interaction between the conduction and valence bands, the spin-orbit coupling within the valence bands, the intrinsic spin Zeeman splitting, and the sp-d exchange interaction between the carriers and magnetic ion in the mean-field approximation. The size dependence of the electron and hole energy levels as well as the giant Zeeman splitting energies are studied theoretically. We find that the hole giant Zeeman splitting energies decrease with the increasing radius, smaller than that in the bulk material, and are different for different J(z) states, which are caused by the quantum confinement effect. Because the quantum sphere restrains the excited Landau states and exciton states, in the experiments we can observe directly the Zeeman splitting of basic states. At low magnetic field, the total Zeeman splitting energy increases linearly with the increasing magnetic field and saturates at modest field which is in agreement with recent experimental results. Comparing to the undoped case, the Zeeman splitting energy is 445 times larger which provides us with wide freedom to tailor the electronic structure of DMS nanocrystals for technological applications.

Design of shallow acceptors in ZnO: First-principles band-structure calculations

p-type doping is a great challenge for the full utilization of ZnO as short-wavelength optoelectronic material. Due to a large electronegative characteristic of oxygen, the ionization energy of acceptors in ZnO is usually too high. By analyzing the defect wave-function character, we propose several approaches to lower the acceptor ionization energy by codoping acceptors with donor or isovalent atoms. Using the first-principles band-structure method, we show that the acceptor transition energies of V-Zn-O-O can be reduced by introducing F-O next to V-Zn to reduce electronic potential, whereas the acceptor transition energy of N-O-nZn(Zn) (n=1-4) can be reduced if we replace Zn by isovalent Mg or Be to reduce the anion and cation kinetic p-d repulsion, as well as the electronic potential.

Comparative study for colloidal quantum dot conduction band state calculations

By comparing the results of some well-controlled calculation methods, we analyze the relative importance of bulk band structure, multi-bulk-band coupling, and boundary conditions in determining colloidal quantum dot conduction band eigenenergies. We find that while the bulk band structure and correct boundary conditions are important, the effects of multi-bulk-band coupling are small.

Variational calculations of ionized-donor-bound excitons in GaAs-AlxGa1-xAs quantum wells

Using a two-parameter wave function, we calculate variationally the binding energy of an exciton bound to an ionized donor impurity (D+,X) in GaAs-AlxGa1-xAs quantum wells for the values of the well width from 10 to 300 Angstrom, when the dopant is located in the center of the well and at the edge of the well. The theoretical results confirm that the previous experimental speculation proposed by Reynolds tit al. [Phys. Rev. B 40, 6210 (1989)] is the binding energy of D+,X for the dopant at the edge of the well. in addition, we also calculate the center-of-mass wave function of the exciton and the average interparticle distances. The results are discussed in detail.

Electronic structure and magnetic coupling properties of Gd-doped AlN: first-principles calculations

In this work, the electronic structure and magnetic coupling properties of Gd doped AlN have been investigated using first-principles method. We found that in the AlN:Gd system, due to the s-f coupling allowed by the symmetry, the exchange splitting of the conduction band is much larger than that of the valence band, which makes the electron-mediated ferromagnetism possible in this material. This property is also confirmed by the energy differences between anti-ferromagnetic and ferromagnetic phase for Al14Gd2N16 with different concentrations of electrons (holes), as well as by the calculated exchange constants. The result indicates that Gd-doped AlN is a promising candidate for the applications in future spintronic devices.

Band crossing in isovalent semiconductor alloys with large size mismatch: First-principles calculations of the electronic structure of Bi and N incorporated GaAs

For large size- and chemical-mismatched isovalent semiconductor alloys, such as N and Bi substitution on As sites in GaAs, isovalent defect levels or defect bands are introduced. The evolution of the defect states as a function of the alloy concentration is usually described by the popular phenomenological band anticrossing (BAC) model. Using first-principles band-structure calculations we show that at the impurity limit the N-(Bi)-induced impurity level is above (below) the conduction- (valence-) band edge of GaAs. These trends reverse at high concentration, i.e., the conduction-band edge of GaAs1-xNx becomes an N-derived state and the valence-band edge of GaAs1-xBix becomes a Bi-derived state, as expected from their band characters. We show that this band crossing phenomenon cannot be described by the popular BAC model but can be naturally explained by a simple band broadening picture.

INTERVALLEY-GAMMA-CHI DEFORMATION POTENTIALS IN III-V ZINCBLENDE AND SI SEMICONDUCTORS BY AB-INITIO PSEUDOPOTENTIAL CALCULATIONS

Intervalley GAMMA - X deformation potential constants (IVDP's) have been calculated by first principle pseudopotential method for the III-V zincblende semiconductors AlP, AlAs, AlSb, GaP, GaAs, GaSb, InP, InAs and InSb. As a prototype crystal we have also carried out calculations on Si. When comparing the calculated IVDP's of LA phonon for GaP, InP and InAs and LO phonon for AlAs, AlSb, GaAs, GaSb and InSb with a previous calculation by EPM in rigid approximation, good agreements are found. However, our ab initio pseudopotential results of LA phonon for AlAs, AlSb, GaAs, GaSb and InSb and LO phonon for GaP, InP and InAs are about one order of magnitude smaller than those obtained by EPM calculations, which indicate that the electron redistributions upon the phonon deformations may be important in affecting GAMMA - X intervalley shatterings for these phonon modes when the anions are being displaced. In our calculations the phonon modes of LA and LO at X point have been evaluated in frozen phonon approximation. We have obtained, at the same time, the LAX and LOX phonon frequencies for these materials from total energy calculations. The calculated phonon frequencies agree very well with experimental values for these semiconductors.

INTERVALLEY GAMMA-CHI DEFORMATION POTENTIALS IN III-V ZINCBLENDE SEMICONDUCTORS BY ABINITIO PSEUDOPOTENTIAL CALCULATIONS

Intervalley GAMMA-X deformation-potential constants (IVDP's) have been calculated by use of a first-principles pseudopotential method for the III-V zinc-blende semiconductors AlP, AlAs, AlSb, GaP, GaAs, GaSb, InP, InAs, and InSb. When the calculated IVDP's of LA phonons for GaP, InP, and InAs and of LO phonons for AlAs, AlSb, GaAs, GaSb, and InSb are compared with results of a previous calculation that used the empirical pseudopotential method (EPM) and a rigid-ion approximation, good agreement is found. However, our ab initio pseudopotential results on IVDP's of LA phonons for AlAs, AlSb, GaAs, GaSb, and InSb and of LO phonons for GaP, InP, and InAs are about one order of magnitude smaller than those obtained by use of EPM calculations, indicating that the electron redistribution accompanying crystal-lattice deformation has a significant effect on GAMMA-X intervalley scattering for these phonon modes when the anions are being displaced. In our calculations the LA- and LO-phonon modes at the X point have been evaluated in the frozen-phonon approximation. We have also obtained the LAX- and LOX-phonon frequencies for these materials from total-energy calculations, which agree very well with experimental values for these semiconductors. We have also calculated GAMMA-X hole-phonon scattering matrix elements for the top valence bands in these nine semiconductors, from which the GAMMA-X IVDP's of the top valence bands for the longitudinal phonons and transverse phonons are evaluated, respectively.

ENERGY-LEVEL CALCULATIONS OF THE RECONSTRUCTED 90-DEGREES PARTIAL DISLOCATION IN SILICON

An LCAO scheme taking into account 10 atomic orbitals (s-, p-, and d-type) applied to a supercell containing 256 atoms is used to calculate the bound states of the reconstructed 90-degrees partial dislocation in Si. The results differ significantly from our earlier calculations on the unreconstructed 90-degrees partial using the same method. We find two bands separate from each other in the entire Brillouin zone and the upper band penetrates deep into the indirect band gap which is in contradiction with the general opinion that core reconstruction clears the band gap of dislocation states.

1ST PRINCIPLE CALCULATIONS OF QUASI-PARTICLE ENERGIES FOR BAND STRUCTURES OF SI AND GAAS

We have applied the Green function theory in GW approximation to calculate the quasiparticle energies for semiconductors Si and GaAs. Good agreements of the calculated excitation energies and fundamental energy gaps with the experimental band structures were achieved. We obtained the calculated fundamental gaps of Si and GaAs to be 1.22 and 1.42 eV in comparison to the experimental values of 1.17 and 1.52 eV, respectively. Ab initio pseudopotential method has been used to generate basis wavefunctions and charge densities for calculating dielectric matrix elements and electron self-energies.

1ST PRINCIPLES PSEUDOPOTENTIAL CALCULATIONS OF ZONE BOUNDARY PHONON FREQUENCIES OF III-V-ZINCBLENDE SEMICONDUCTORS

Longitudinal zone boundary X phonon frequencies have been calculated by a first principles pseudopotential method for III-V zincblende semiconductors AlP, AlAs, AlSb, GaP, GaAs, GaSb, InP, InAs and InSb. The phonon frequencies have been evaluated from total energy calculations in the frozen phonon approximation. The calculated phonon frequencies agree very well with the experimental values.

1ST-PRINCIPLES CALCULATIONS FOR QUASI-PARTICLE ENERGIES OF GAP AND GAAS

We have applied the Green-function method in the GW approximation to calculate quasiparticle energies for the semiconductors GaP and GaAs. Good agreement between the calculated excitation energies and the experimental results was achieved. We obtained calculated direct band gaps of GaP and GaAs of 2.93 and 1.42 eV, respectively, in comparison with the experimental values of 2.90 and 1.52 eV, respectively. An ab initio pseudopotential method has been used to generate basis wave functions and charge densities for calculating the dielectric matrix elements and self-enegies. To evaluate the dynamical effects of the screened interaction, the generalized-plasma-pole model has been utilized to extend the dielectric matrix elements from static results to finite frequencies. We presen the calculated quasiparticle energies at various high-symmetry points of the Brillouin zone and compare them with the experimental results and other calculations.

1ST PRINCIPLE CALCULATIONS OF QUASI-PARTICLE ENERGIES FOR BAND STRUCTURES OF SEMICONDUCTORS

We successfully applied the Green function theory in GW approximation to calculate the quasiparticle energies for semiconductors Si and GaAs. Ab initio pseudopotential method was adopted to generate basis wavefunctions and charge densities for calculating dielectric matrix elements and electron self-energies. To evaluate dynamical effects of screened interaction, GPP model was utilized to extend dieletric matrix elements from static results to finite frequencies. We give a full account of the theoretical background and the technical details for the first principle pseudopotential calculations of quasiparticle energies in semiconductors and insulators. Careful analyses are given for the effective and accurate evaluations of dielectric matrix elements and quasiparticle self-energies by using the symmetry properties of basis wavefunctions and eigenenergies. Good agreements between the calculated excitation energies and fundamental energy gaps and the experimental band structures were achieved.

2-DELTA-FUNCTION GENERALIZED PLASMA POLE MODEL FOR 1ST PRINCIPLE CALCULATIONS OF QUASI-PARTICLE ENERGIES IN MANY-ELECTRON SYSTEMS

To evaluate the dynamical effects of the screened interaction in the calculations of quasiparticle energies in many-electron systems a two-delta-function generalized plasma pole model (GPP) is introduced to simulate the dynamical dielectric function. The usual single delta-function GPP model has the drawback of over simplifications and for the crystals without the center of symmetry is inappropriate to describe the finite frequency behavior for dielectric function matrices. The discrete frequency summation method requires too much computation to achieve converged results since ab initio calculations of dielectric function matrices are to be carried out for many different frequencies. The two-delta GPP model is an optimization of the two approaches. We analyze the two-delta GPP model and propose a method to determine from the first principle calculations the amplitudes and effective frequencies of these delta-functions. Analytical solutions are found for the second order equations for the parameter matrices entering the model. This enables realistic applications of the method to the first principle quasiparticle calculations and makes the calculations truly adjustable parameter free.