No Arabic abstract
This paper presents a numerical implementation of a first-principles envelope-function theory derived recently by the author [B. A. Foreman, Phys. Rev. B 72, 165345 (2005)]. The examples studied deal with the valence subband structure of GaAs/AlAs, GaAs/Al(0.2)Ga(0.8)As, and In(0.53)Ga(0.47)As/InP (001) superlattices calculated using the local density approximation to density-functional theory and norm-conserving pseudopotentials without spin-orbit coupling. The heterostructure Hamiltonian is approximated using quadratic response theory, with the heterostructure treated as a perturbation of a bulk reference crystal. The valence subband structure is reproduced accurately over a wide energy range by a multiband envelope-function Hamiltonian with linear renormalization of the momentum and mass parameters. Good results are also obtained over a more limited energy range from a single-band model with quadratic renormalization. The effective kinetic-energy operator ordering derived here is more complicated than in many previous studies, consisting in general of a linear combination of all possible operator orderings. In some cases the valence-band Rashba coupling differs significantly from the bulk magnetic Luttinger parameter. The splitting of the quasidegenerate ground state of no-common-atom superlattices has non-negligible contributions from both short-range interface mixing and long-range dipole terms in the quadratic density response.
In this paper a multi-band envelope-function Hamiltonian for lattice-matched semiconductor heterostructures is derived from first-principles norm-conserving pseudopotentials. The theory is applicable to isovalent or heterovalent heterostructures with macroscopically neutral interfaces and no spontaneous bulk polarization. The key assumption -- proved in earlier numerical studies -- is that the heterostructure can be treated as a weak perturbation with respect to some periodic reference crystal, with the nonlinear response small in comparison to the linear response. Quadratic response theory is then used in conjunction with k.p perturbation theory to develop a multi-band effective-mass Hamiltonian (for slowly varying envelope functions) in which all interface band-mixing effects are determined by the linear response. To within terms of the same order as the position dependence of the effective mass, the quadratic response contributes only a bulk band offset term and an interface dipole term, both of which are diagonal in the effective-mass Hamiltonian. Long-range multipole Coulomb fields arise in quantum wires or dots, but have no qualitative effect in two-dimensional systems beyond a dipole contribution to the band offsets.
The valence band offsets (VBO) for the b{eta}-type A3B6 layered compounds depending on the thickness of the crystals have been investigated from the first principles, based on the density functional theory. To simulate the structure of a given thickness the periodic slab model was used. Two adjacent crystal slabs consisting of several layers were separated by a vacuum region of two-layer width. It is shown that at the crystal thickness more than 12 layers, photothreshold practically becomes independent on the thickness of the crystal.
We study nonideal mixing effects in the regime of warm dense matter (WDM) by computing the shock Hugoniot curves of BN, MgO, and MgSiO_3. First, we derive these curves from the equations of state (EOS) of the fully interacting systems, which were obtained using a combination of path integral Monte Carlo calculations at high temperature and density functional molecular dynamics simulations at lower temperatures. We then use the ideal mixing approximation at constant pressure and temperature to rederive these Hugoniot curves from the EOS tables of the individual elements. We find that the linear mixing approximation works remarkably well at temperatures above ~2*10^5 K, where the shock compression ratio exceeds ~3.2. The shape of the Hugoniot curve of each compound is well reproduced. Regions of increased shock compression, that emerge because of the ionization of L and K shell electrons, are well represented and the maximum compression ratio on the Hugoniot curves is reproduced with high precision. Some deviations are seen near the onset of the L shell ionization regime, where ionization equilibrium in the fully interacting system cannot be well reproduced by the ideal mixing approximation. This approximation also breaks down at lower temperatures, where chemical bonds play an increasingly import role. However, the results imply that equilibrium properties of binary and ternary mixtures in the regime of WDM can be derived from the EOS tables of the individual elements. This significantly simplifies the characterization of binary and ternary mixtures in the WDM and plasma phases, which otherwise requires large numbers of more computationally expensive first-principles computer simulations.
We use a recently developed self-consistent GW approximation to present first principles calculations of the conduction band spin splitting in GaAs under [110] strain. The spin orbit interaction is taken into account as a perturbation to the scalar relativistic hamiltonian. These are the first calculations of conduction band spin splitting under deformation based on a quasiparticle approach; and because the self-consistent GW scheme accurately reproduces the relevant band parameters, it is expected to be a reliable predictor of spin splittings. We also discuss the spin relaxation time under [110] strain and show that it exhibits an in-plane anisotropy, which can be exploited to obtain the magnitude and sign of the conduction band spin splitting experimentally.
The band offsets between crystalline and hydrogenated amorphous silicon (a-Si:H) are key parameters governing the charge transport in modern silicon hetrojunction solar cells. They are an important input for macroscopic simulators that are used to further optimize the solar cell. Past experimental studies, using X-ray photoelectron spectroscopy (XPS) and capacitance-voltage measurements, have yielded conflicting results on the band offset. Here we present a computational study on the band offsets. It is based on atomistic models and density-functional theory (DFT). The amorphous part of the interface is obtained by relatively long DFT first-principles molecular-dynamics (MD) runs at an elevated temperature on 30 statistically independent samples. In order to obtain a realistic conduction band position the electronic structure of the interface is calculated with a hybrid functional. We find a slight asymmetry in the band offsets, where the offset in the valence band (0.30 eV) is larger than in the conduction band (0.17 eV). Our results are in agreement with the latest XPS measurements that report a valence band offset of 0.3 eV [M. Liebhaber et al., Appl. Phys. Lett. 106, 031601 (2015)].