Local-Density-Derived Semiempirical Nonlocal Pseudopotentials for InP with Applications to Large Quantum Dots

In the same way that atomic calculations have been used previously to extract bare ionic pseudopotentials, self-consistent bulk calculations can be used to construct screened atomic pseudopotentials. We use such a method to construct screened nonlocal atomic pseudopotentials for InP. A series of bulk, local-density-approximation (LDA) calculations are performed on a few InP crystal structures,covering a range of unit-cell volumes, to produce bulk potentials {VLDA(G)}. By solving a set of linear equations, we extract from these crystalline potentials the corresponding screened atomic 'spherical LDA' (SLDA) potentials VaSLDA (|q|) for sites a=In or P. In combination with the nonlocal part of the usual LDA pseudopotentials, these SLDA potentials give band structures and wave functionsthat are virtually indistinguishable from the self-consistent LDA results for bulk InP. In the next step, we apply linear changes to the local SLDA potentials (while keeping the nonlocal potentials at their LDA values), to fit the band structures to experiment. Interestingly, this removal of LDA eigenvalue errors requires only small and subtle changes in the potential--mostly an upshift in theregion near the cation core, with nearly no change at the bond center. Furthermore, the linear changes to the SLDA potentials result mostly in an upshift of the conduction bands with little effect on the valence bands. Because only small changes in the potential suffice to fit the bands to experimental results, the wave functions remain virtually unchanged relative to those in the original LDAcalculation. Hence, we obtain semiempirical pseudopotentials which can produce ab initio LDA-quality wave functions with experimentally measured band structures, effective masses, and deformation potentials. The potentials obtained here were deposited on a FTP site and can be used by interested readers. Since the resulting pseudopotentials are 'soft' (with small high-momentum components), theycan be applied within a plane-wave basis in combination with a Gaussian correction to large systems for which LDA calculations are prohibitively expensive. As an illustration, we apply our InP screened atomic pseudopotentials to calculate quantum size effects on the band gaps of InP dots with sizes up to 700 atoms. Good agreement is found between the theoretical and the experimental band gaps.Fitting the calculated band gaps E8 (in unit of eV) versus the effect dot sizes D (in unit of .ANG.) gives E8=1.45+37.295/D1.16. This prediction differs significantly from the quadratic size dependence D.hiven.2.0 expected from simple effective-mass theory.