No Arabic abstract
We derive first- and second-order piezoelectric coefficients for the zinc-blende III-V semiconductors, {Al,Ga,In}-{N,P,As,Sb}. The results are obtained within the Heyd-Scuseria-Ernzerhof hybrid-functional approach in the framework of density functional theory and the Berry-phase theory of electric polarization. To achieve a meaningful interpretation of the results, we build an intuitive phenomenological model based on the description of internal strain and the dynamics of the electronic charge centers. We discuss in detail first- and second-order internal strain effects, together with strain-induced changes in ionicity. This analysis reveals that the relatively large importance in the III-Vs of non-linear piezoelectric effects compared to the linear ones arises because of a delicate balance between the ionic polarization contribution due to internal strain relaxation effects, and the contribution due to the electronic charge redistribution induced by macroscopic and internal strain.
The direct calculation of the elastic and piezoelectric tensors of solids can be accomplished by treating homogeneous strain within the framework of density-functional perturbation theory. By formulating the energy functional in reduced coordinates, we show that the strain perturbation enters only through metric tensors, and can be treated in a manner exactly paralleling the treatment of other perturbations. We present an analysis of the strain perturbation of the plane-wave pseudopotential functional, including the internal strain terms necessary to treat the atomic-relaxation contributions. Procedures for computationally verifying these expressions by comparison with numerical derivatives of ground-state calculations are described and illustrated.
We demonstrate that Daubechies wavelets can be used to construct a minimal set of optimized localized contracted basis functions in which the Kohn-Sham orbitals can be represented with an arbitrarily high, controllable precision. Ground state energies and the forces acting on the ions can be calculated in this basis with the same accuracy as if they were calculated directly in a Daubechies wavelets basis, provided that the amplitude of these contracted basis functions is sufficiently small on the surface of the localization region, which is guaranteed by the optimization procedure described in this work. This approach reduces the computational costs of DFT calculations, and can be combined with sparse matrix algebra to obtain linear scaling with respect to the number of electrons in the system. Calculations on systems of 10,000 atoms or more thus become feasible in a systematic basis set with moderate computational resources. Further computational savings can be achieved by exploiting the similarity of the contracted basis functions for closely related environments, e.g. in geometry optimizations or combined calculations of neutral and charged systems.
We obtain parameters for non-orthogonal and orthogonal TB models from two-atomic molecules for all combinations of elements of period 1 to 6 and group 3 to 18 of the periodic table. The TB bond parameters for 1711 homoatomic and heteroatomic dimers show clear chemical trends. In particular, using our parameters we compare to the rectangular d-band model, the reduced sp TB model as well as canonical TB models for sp- and d-valent systems which have long been used to gain qualitative insight into the interatomic bond. The transferability of our dimer-based TB bond parameters to bulk systems is discussed exemplarily for the bulk ground-state structures of Mo and Si. Our dimer-based TB bond parameters provide a well-defined and promising starting point for developing refined TB parameterizations and for making the insight of TB available for guiding materials design across the periodic table.
We present an implementation of time-dependent density-functional theory (TDDFT) in the linear response formalism enabling the calculation of low energy optical absorption spectra for large molecules and nanostructures. The method avoids any explicit reference to canonical representations of either occupied or virtual Kohn-Sham states and thus achieves linear-scaling computational effort with system size. In contrast to conventional localised orbital formulations, where a single set of localised functions is used to span the occupied and unoccupied state manifold, we make use of two sets of in situ optimised localised orbitals, one for the occupied and one for the unoccupied space. This double representation approach avoids known problems of spanning the space of unoccupied Kohn-Sham states with a minimal set of localised orbitals optimised for the occupied space, while the in situ optimisation procedure allows for efficient calculations with a minimal number of functions. The method is applied to a number of medium sized organic molecules and a good agreement with traditional TDDFT methods is observed. Furthermore, linear scaling of computational cost with system size is demonstrated on a system of carbon nanotubes.
Exchange interactions are a manifestation of the quantum mechanical nature of the electrons and play a key role in predicting the properties of materials from first principles. In density functional theory (DFT), a widely used approximation to the exchange energy combines fractions of density-based and Hartree-Fock (exact) exchange. This so-called hybrid DFT scheme is accurate in many materials, for reasons that are not fully understood. Here we show that a 1/4 fraction of exact exchange plus a 3/4 fraction of density-based exchange is compatible with a correct quantum mechanical treatment of the exchange energy of an electron pair in the unpolarized electron gas. We also show that the 1/4 exact-exchange fraction mimics a correlation interaction between doubly-excited electronic configurations. The relation between our results and trends observed in hybrid DFT calculations is discussed, along with other implications.