Search

Cu$_2$ZnSnS$_4$ is an earth-abundant photovoltaic absorber material predicted to provide a sustainable solution for commercial solar applications. One of the main limitations restricting its commercialization is the issue of cation disorder. Raman spectroscopy has been a sought after technique to characterize disorder in CZTS while a clear consensus between theoretical and experimental results is yet to be achieved. In the present study, via the virtual crystal approximation, we take into account the progressive nature of Cu-Zn disorder in CZTS: we obtain the phonon frequencies at zone-center within the density functional perturbation theory formalism, and further compute the Raman spectra for the disordered phases, achieving a consensus between theory and experiment. These calculations confirm the presence of complete disorder in Cu-Zn 2$a$, 2$c$ and 2$d$ Wyckoff sites. They also show that the Raman intensities of two prominent $A$ phonon modes characterized by motion of S atoms, also known to be experimentally significant, play a key role in understanding the nature of disorder in CZTS., Comment: 16 pages, 5 figures (+7 figures in the appendix)

First-principles calculations in crystalline structures are often performed with a planewave basis set. To make the number of basis functions tractable two approximations are usually introduced: core electrons are frozen and the diverging Coulomb potential near the nucleus is replaced by a smoother expression. The norm-conserving pseudopotential was the first successful method to apply these approximations in a fully ab initio way. Later on, more efficient and more exact approaches were developed based on the ultrasoft and the projector augmented wave formalisms. These formalisms are however more complex and developing new features in these frameworks is usually more difficult than in the norm-conserving framework. Most of the existing tables of norm- conserving pseudopotentials, generated long ago, do not include the latest developments, are not systematically tested or are not designed primarily for high accuracy. In this paper, we present our PseudoDojo framework for developing and testing full tables of pseudopotentials, and demonstrate it with a new table generated with the ONCVPSP approach. The PseudoDojo is an open source project, building on the AbiPy package, for developing and systematically testing pseudopotentials. At present it contains 7 different batteries of tests executed with ABINIT, which are performed as a function of the energy cutoff. The results of these tests are then used to provide hints for the energy cutoff for actual production calculations. Our final set contains 141 pseudopotentials split into a standard and a stringent accuracy table. In total around 70.000 calculations were performed to test the pseudopotentials. The process of developing the final table led to new insights into the effects of both the core-valence partitioning and the non-linear core corrections on the stability, convergence, and transferability of norm-conserving pseudopotentials. ..., Comment: abstract truncated, 17 pages, 25 figures, 8 tables

The search for new materials, based on computational screening, relies on methods that accurately predict, in an automatic manner, total energy, atomic-scale geometries, and other fundamental characteristics of materials. Many technologically important material properties directly stem from the electronic structure of a material, but the usual workhorse for total energies, namely density-functional theory, is plagued by fundamental shortcomings and errors from approximate exchange-correlation functionals in its prediction of the electronic structure. At variance, the $GW$ method is currently the state-of-the-art {\em ab initio} approach for accurate electronic structure. It is mostly used to perturbatively correct density-functional theory results, but is however computationally demanding and also requires expert knowledge to give accurate results. Accordingly, it is not presently used in high-throughput screening: fully automatized algorithms for setting up the calculations and determining convergence are lacking. In this work we develop such a method and, as a first application, use it to validate the accuracy of $G_0W_0$ using the PBE starting point, and the Godby-Needs plasmon pole model ($G_0W_0^\textrm{GN}$@PBE), on a set of about 80 solids. The results of the automatic convergence study utilized provides valuable insights. Indeed, we find correlations between computational parameters that can be used to further improve the automatization of $GW$ calculations. Moreover, we find that $G_0W_0^\textrm{GN}$@PBE shows a correlation between the PBE and the $G_0W_0^\textrm{GN}$@PBE gaps that is much stronger than that between $GW$ and experimental gaps. However, the $G_0W_0^\textrm{GN}$@PBE gaps still describe the experimental gaps more accurately than a linear model based on the PBE gaps., Comment: 12 pages, 11 figures

Cu2ZnSnS4 (CZTS) is a promising material as an absorber in photovoltaic applications. The measured efficiency, however, is far from the theoretically predicted value for the known CZTS phases. To improve the understanding of this discrepancy we investigate the structural, dynamical, and dielectric of the three main phases of CZTS (kesterite, stannite, and PMCA) using density functional perturbation theory (DFPT). The effect of the exchange-correlation functional on the computed properties is analyzed. A qualitative agreement of the theoretical Raman spectrum with measurements is observed. However, none of the phases correspond to the experimental spectrum within the error bar that is usually to be expected for DFPT. This corroborates the need to consider cation disorder and other lattice defects extensively in this material.

Structural, electronic, vibrational and dielectric properties of LaBGeO$_5$ with the stillwellite structure are determined based on \textit{ab initio} density functional theory. The theoretically relaxed structure is found to agree well with the existing experimental data with a deviation of less than $0.2\%$. Both the density of states and the electronic band structure are calculated, showing five distinct groups of valence bands. Furthermore, the Born effective charge, the dielectric permittivity tensors, and the vibrational frequencies at the center of the Brillouin zone are all obtained. Compared to existing model calculations, the vibrational frequencies are found in much better agreement with the published experimental infrared and Raman data, with absolute and relative rms values of 6.04 cm$^{-1}$, and $1.81\%$, respectively. Consequently, numerical values for both the parallel and perpendicular components of the permittivity tensor are established as 3.55 and 3.71 (10.34 and 12.28), respectively, for the high-(low-)frequency limit.

With the ever-increasing sophistication of codes, the verification of the implementation of advanced theoretical formalisms becomes critical. In particular, cross comparison between different codes provides a strong hint in favor of the correctness of the implementations, and a measure of the (hopefully small) possible numerical differences. We lead a rigorous and careful study of the quantities that enter in the calculation of the zero-point motion renormalization of the direct band gap of diamond due to electron-phonon coupling, starting from the total energy, and going through the computation of phonon frequencies and electron-phonon matrix elements. We rely on two independent implementations : Quantum Espresso + Yambo and ABINIT. We provide the order of magnitude of the numerical discrepancies between the codes, that are present for the different quantities: less than $10^{-5}$ Hartree per atom on the total energy (-5.722 Ha/at), less than 0.07 cm$^{-1}$ on the $\Gamma,L,X$ phonon frequencies (555 to 1330 cm$^{-1}$), less than 0.5% on the square of the electron-phonon matrix elements and less than 4 meV on the zero-point motion renormalization of each eigenenergies (44 to 264 meV). Within our approximations, the DFT converged direct band gap renormalization in diamond due to the electron-phonon coupling is -0.409 eV (reduction of the band gap)., Comment: 11 pages, 5 figures

Combining LSDA+$U$ and an analysis of superexchange interactions beyond DFT, we describe the magnetic ground states in rutile and anatase Cr-doped TiO$_2$. In parallel, we correct our LSDA+$U$ ground state through GW corrections ($GW$@LSDA+$U$) that reproduce the position of impurity states and the band gaps in satisfying agreement with experiments. Because of the different topological coordinations of Cr-Cr bonds in the ground states of rutile and anatase, superexchange interactions induce either ferromagnetic or antiferromagnetic couplings of Cr ions. In Cr-doped anatase, this interaction leads to a new mechanism which stabilizes a ferromagnetic ground state, in keeping with experimental evidence, without the need to invoke F-center exchange., Comment: 5

The theoretical treatment of homogeneous static magnetic fields in periodic systems is challenging, as the corresponding vector potential breaks the translational invariance of the Hamiltonian. Based on density operators and perturbation theory, we propose, for insulators, a periodic framework for the treatment of magnetic fields up to arbitrary order of perturbation, similar to widely used schemes for electric fields. The second-order term delivers a new, remarkably simple, formulation of the macroscopic orbital magnetic susceptibility for periodic insulators. We validate the latter expression using a tight-binding model, analytically from the present theory and numerically from the large-size limit of a finite cluster, with excellent numerical agreement., Comment: 5 pages including 2 figures; accepted for publication in Phys. Rev. B

We extend recent \textit{ab initio} calculations of the electronic band structure and the phonon dispersion relations of rhombohedral GeTe to calculations of the density of phonon states and the temperature dependent specific heat. The results are compared with measurements of the specific heat. It is discovered that the specific heat depends on hole concentration, not only in the very low temperature region (Sommerfeld term) but also at the maximum of $C_p/T^3$ (around 16 K). To explain this phenomenon, we have performed \textit{ab initio} lattice dynamical calculations for GeTe rendered metallic through the presence of a heavy hole concentration ($p$ $\sim$ 2$\times$ 10$^{21}$ cm$^{-3}$). They account for the increase observed in the maximum of $C_p/T^3$., Comment: 8 pages, 7 figures, ref. 19 corrected

During the past five years the low temperature heat capacity of simple semiconductors and insulators has received renewed attention. Of particular interest has been its dependence on isotopic masses and the effect of spin- orbit coupling in ab initio calculations. Here we concentrate on the lead chalcogenides PbS, PbSe and PbTe. These materials, with rock salt structure, have different natural isotopes for both cations and anions, a fact that allows a systematic experimental and theoretical study of isotopic effects e.g. on the specific heat. Also, the large spin-orbit splitting of the 6p electrons of Pb and the 5p of Te allows, using a computer code which includes spin-orbit interaction, an investigation of the effect of this interaction on the phonon dispersion relations and the temperature dependence of the specific heat and on the lattice parameter. It is shown that agreement between measurements and calculations significantly improves when spin-orbit interaction is included., Comment: 25 pages, 12 Figures, 1 table, submitted to PRB

The dynamical, dielectric and elastic properties of GeTe, a ferroelectric material in its low temperature rhombohedral phase, have been investigated using first-principles density functional theory. We report the electronic energy bands, phonon dispersion curves, electronic and low frequency dielectric tensors, infra-red reflectivity, Born effective charges, elastic and piezoelectric tensors and compare them with the existing theoretical and experimental results, as well as with similar quantities available for other ferroelectric materials, when appropriate.

We use many-body perturbation theory, the state-of-the-art method for band gap calculations, to compute the band offsets at the Si/SiO$_2$ interface. We examine the adequacy of the usual approximations in this context. We show that (i) the separate treatment of band-structure and potential lineup contributions, the latter being evaluated within density-functional theory, is justified, (ii) most plasmon-pole models lead to inaccuracies in the absolute quasiparticle corrections, (iii) vertex corrections can be neglected, (iv) eigenenergy self-consistency is adequate. Our theoretical offsets agree with the experimental ones within 0.3 eV.

In order to allow different software applications, in constant evolution, to interact and exchange data, flexible file formats are needed. A file format specification for different types of content has been elaborated to allow communication of data for the software developed within the European Network of Excellence "NANOQUANTA", focusing on first-principles calculations of materials and nanosystems. It might be used by other software as well, and is described here in detail. The format relies on the NetCDF binary input/output library, already used in many different scientific communities, that provides flexibility as well as portability accross languages and platforms. Thanks to NetCDF, the content can be accessed by keywords, ensuring the file format is extensible and backward compatible.

In recent years, there has been increasing interest in the specific heat $C$ of insulators and semiconductors because of the availability of samples with different isotopic masses and the possibility of performing \textit{ab initio} calculations of its temperature dependence $C(T)$ using as a starting point the electronic band structure. Most of the crystals investigated are elemental (e.g., germanium) or binary (e.g., gallium nitride) semiconductors. The initial electronic calculations were performed in the local density approximation and did not include spin-orbit interaction. Agreement between experimental and calculated results was usually found to be good, except for crystals containing heavy atoms (e.g., PbS) for which discrepancies of the order of 20% existed at the low temperature maximum found for $C/T^3$. It has been conjectured that this discrepancies result from the neglect of spin-orbit interaction which is large for heavy atoms ($\Delta_0\sim$1.3eV for the $p$ valence electrons of atomic lead). Here we discuss measurements and \textit{ab initio} calculations of $C(T)$ for crystalline bismuth ($\Delta_0\sim$1.7 eV), strictly speaking a semimetal but in the temperature region accessible to us ($T >$ 2K) acting as a semiconductor. We extend experimental data available in the literature and notice that the \textit{ab initio} calculations without spin-orbit interaction exhibit a maximum at $\sim$8K, about 20% lower than the measured one. Inclusion of spin-orbit interaction decreases the discrepancy markedly: The maximum of $C(T)$ is now only 7% larger than the measured one. Exact agreement is obtained if the spin-orbit hamiltonian is reduced by a factor of $\sim$0.8., Comment: 4 pages, 3 figures

We show that density functional theory within the RPA (random phase approximation for the exchange-correlation energy) provides a correct description of bond dissociation in H$_2$ in a spin-restricted Kohn-Sham formalism, i.e. without artificial symmetry breaking. We present accurate adiabatic connection curves both at equilibrium and beyond the Coulson-Fisher point. The strong curvature at large bond length implies important static (left-right) correlation, justifying modern hybrid functional constructions but also demonstrating their limitations. Although exact at infinite and accurate around the equilibrium bond length, the RPA dissociation curve displays unphysical repulsion at larger but finite bond lengths. Going beyond the RPA by including the exact exchange kernel (RPA+X), we find a similar repulsion. We argue that this deficiency is due to the absence of double excitations in adiabatic linear response theory. Further analyzing the H$_2$ dissociation limit we show that the RPA+X is not size-consistent, in contrast to the RPA., Comment: 15 pages, 5 figures

The non-linear response of infinite periodic solids to homogenous electric fields and collective atomic displacements is discussed in the framework of density functional perturbation theory. The approach is based on the 2n + 1 theorem applied to an electric-field-dependent energy functional. We report the expressions for the calculation of the non-linear optical susceptibilities, Raman scattering efficiencies and electrooptic coefficients. Different formulations of third-order energy derivatives are examined and their convergence with respect to the k-point sampling is discussed. We apply our method to a few simple cases and compare our results to those obtained with distinct techniques. Finally, we discuss the effect of a scissors correction on the EO coefficients and non-linear optical susceptibilities.

Using a plane wave pseudopotential approach to density functional theory we investigate the electron localization length in various oxides. For this purpose, we first set up a theory of the band-by-band decomposition of this quantity, more complex than the decomposition of the spontaneous polarization (a related concept), because of the interband coupling. We show its interpretation in terms of Wannier functions and clarify the effect of the pseudopotential approximation. We treat the case of different oxides: BaO, $\alpha$-PbO, BaTiO$_3$ and PbTiO$_3$. We also investigate the variation of the localization tensor during the ferroelectric phase transitions of BaTiO$_3$ as well as its relationship with the Born effective charges.

The Born effective charge, Z*, that describes the polarization created by collective atomic displacements, can be computed from first-principles following different techniques. We establish the connections existing between these different formulations, and analyze the related band-by-band decompositions. We show that unlike for the full Z*, the different band-by-band values are not equal, and emphasize that one of them has a natural physical meaning in terms of Wannier functions., Comment: 9 pages

Using first-principles calculations based on a variational density functional perturbation theory, we investigate the lattice dynamics of solid solutions of barium and strontium titanates. Averaging the information available for the related pure compounds yields results equivalent to those obtained within the virtual crystal approximation, providing frequencies which are a good approximation to those computed for a (111) ordered supercell. Using the same averaging technique we report the evolution of the ferroelectric and antiferrodistortive instabilities with composition., Comment: 9 pages, 2 figures, Proceedings for Fundamental Physics of Ferroelectrics, Aspen (CO), Feb. 13-20, 2000

Specific matrix elements of exchange and correlation kernels in time-dependent density-functional theory are computed. The knowledge of these matrix elements not only constraints approximate time-dependent functionals, but also allows to link different practical approaches to excited states, either based on density-functional theory, or on many-body perturbation theory, despite the approximations that have been performed to derive them., Comment: Submitted to Phys. Rev. Lett. (February 4, 1999). Other related publications can be found at http://www.fhi-berlin.mpg.de/th/paper.html

Based on recent first-principles computations in perovskite compounds, especially BaTiO3, we examine the significance of the Born effective charge concept and contrast it with other atomic charge definitions, either static (Mulliken, Bader...) or dynamical (Callen, Szigeti...). It is shown that static and dynamical charges are not driven by the same underlying parameters. A unified treatment of dynamical charges in periodic solids and large clusters is proposed. The origin of the difference between static and dynamical charges is discussed in terms of local polarizability and delocalized transfers of charge: local models succeed in reproducing anomalous effective charges thanks to large atomic polarizabilities but, in ABO3 compounds, ab initio calculations favor the physical picture based upon transfer of charges. Various results concerning barium and strontium titanates are presented. The origin of anomalous Born effective charges is discussed thanks to a band-by-band decomposition which allows to identify the displacement of the Wannier center of separated bands induced by an atomic displacement. The sensitivity of the Born effective charges to microscopic and macroscopic strains is examined. Finally, we estimate the spontaneous polarization in the four phases of barium titanate., Comment: 25 pages, 6 Figures, 10 Tables, LaTeX

The Born effective charge, Z*, that describes the polarization created by collective atomic displacements, can be computed from first-principles by different techniques. Its band-by-band decomposition appeared recently as a powerful tool for the microscopic characterisation of the bonding in solids. We describe the connections between the different expressions used for the computation of Z*, and analyze the possible associated band-by-band decompositions. We show that unlike for the full Z*, the different band-by-band values are not equal, and emphasize that one of them has an interesting physical meaning., Comment: 11 pages, 1 table, LaTeX

We consider two perturbative schemes to calculate excitation energies, each employing the Kohn-Sham Hamiltonian as the unperturbed system. Using accurate exchange-correlation potentials generated from essentially exact densities and their exchange components determined by a recently proposed method, we evaluate energy differences between the ground state and excited states in first-order perturbation theory for the Helium, ionized Lithium and Beryllium atoms. It was recently observed that the zeroth-order excitations energies, simply given by the difference of the Kohn-Sham eigenvalues, almost always lie between the singlet and triplet experimental excitations energies, corrected for relativistic and finite nuclear mass effects. The first-order corrections provide about a factor of two improvement in one of the perturbative schemes but not in the other. The excitation energies within perturbation theory are compared to the excitations obtained within $\Delta$SCF and time-dependent density functional theory. We also calculate the excitation energies in perturbation theory using approximate functionals such as the local density approximation and the optimized effective potential method with and without the Colle-Salvetti correlation contribution.

The phonon dispersion curves of cubic BaTiO_3 have been computed within a first-principles approach and the results compared to the experimental data. The curves obtained are very similar to those reported for KNbO_3 by Yu and Krakauer [Phys. Rev. Lett. 74, 4067 (1995)]. They reveal that correlated atomic displacements along <100> chains are at the origin of the ferroelectric instability. A simplified model illustrates that spontaneous collective displacements will occur when a dozen of aligned atoms are coupled. The longitudinal interatomic force constant between nearest neighbour Ti and O atoms is relatively weak in comparison to that between Ti atoms in adjacent cells. The small coupling between Ti and O displacements seems however necessary to reproduce a ferroelectric instability., Comment: 12 pages, 4 figures

The polarization-dependence of the exchange-correlation (XC) energy functional of periodic insulators within Kohn-Sham (KS) density-functional theory requires a ${\cal O} (1/q^2)$ divergence in the XC kernel for small vectors q. This behaviour, exemplified for a one-dimensional model semiconductor, is also observed when an insulator happens to be described as a KS metal, or vice-versa. Although it can occur in the exchange-only kernel, it is not found in the usual local, semi-local or even non-local approximations to KS theory. We also show that the test-charge and electronic definitions of the macroscopic dielectric constant differ from one another in exact KS theory, but are equivalent in the above-mentioned approximations., Comment: 7 pages, 1 figure

We examine the density-functional theory of macroscopic insulators, obtained in the large-cluster limit or under periodic boundary conditions. For polar crystals, we find that the two procedures are not equivalent. In a large-cluster case, the exact exchange-correlation potential acquires a homogeneous ``electric field'' which is absent from the usual local approximations, and the Kohn-Sham electronic system becomes metallic. With periodic boundary conditions, such a field is forbidden, and the polarization deduced from Kohn-Sham wavefunctions is incorrect even if the exact functional is used., Comment: 9 pages, 3 figures, to appear in Phys. Rev. Lett

Using first-principles calculations, the phonon frequencies at the $\Gamma$ point and the dielectric tensor are determined and analysed for the cubic and rhombohedral phases of BaTiO$_{3}$. The dipole-dipole interaction is then separated \`a la Cochran from the remaining short-range forces, in order to investigate their respective influence on lattice dynamics. This analysis highlights the delicate balance of forces leading to an unstable phonon in the cubic phase and demonstrates the extreme sensitivity of this close compensation to minute effective charge changes. Within our decomposition, the stabilization of the unstable mode in the rhombohedral phase or under isotropic pressure has a different origin., Comment: 9 pages, 4 tables, 1 figure

Parallel algorithms for ab initio calculations of vibrations modes of solids are presented and implemented under PVM. Load balancing and communication problems are dealt with in order to increase parallelism efficiency. For accurate time measurements, synchronisation of processes must be achieved. The results obtained by working with 1,2,4 and 8 processors of a Convex MetaSeries computer are presented, showing a very good efficiency. Further parallelisation of the codes (involving a three-dimensional Fast Fourier Transform) with compiler directives on a Convex Exemplar, is discussed., Comment: Published in the Proceedings of the 7th European Convex Users' Conference (ECUC '95), 17 pages (figures included)

Motivated by the negative thermal expansion observed for silicon between 20 K and 120 K, we present first an ab initio study of the volume dependence of interatomic force constants, phonon frequencies of TA(X) and TA(L) modes, and of the associated mode Gruneisen parameters. The influence of successive nearest neighbors shells is analysed. Analytical formulas, taking into account interactions up to second nearest neighbors, are developped for phonon frequencies of TA(X) and TA(L) modes and the corresponding mode Gruneisen parameters. We also analyze the volume and pressure dependence of various thermodynamic properties (specific heat, bulk modulus, thermal expansion), and point out the effect of the negative mode Gruneisen parameters of the acoustic branches on these properties. Finally, we present the evolution of the mean square atomic displacement and of the atomic temperature factor with the temperature for different volumes, for which the anomalous effects are even greater., Comment: 24 pages, Revtex 3.0, 11 figures, accepted for publication in Phys. Rev. B

We present a new technique aimed at preventing plane-wave based total energy and stress calculations from the effect of abrupt changes in basis set size. This s cheme relies on the interpolation of energy as a function of the number of plane waves, and on a scaling hypothesis that allows to perform the interpolation for a unique reference volume. From a theoretical point of view, the new method is compared to those already pr oposed in the literature, and its more rigorous derivation is emphasized. From a practical point of view, we illustrate the importance of the correction o n different materials (Si, BaTiO3, and He) corresponding to different types of b onding, and to different k-point samplings and cut-off energies. Then, we compare the different approaches for the calculation of the lattice par ameter, the bulk modulus, and its derivative versus pressure in bulk silicon., Comment: 41 pages, Revtex 3.0, 10 figures

The response of an infinite, periodic, insulating, solid to an infinitesimally small electric field is investigated in the framework of Density Functional Theory. We find that the applied perturbing potential is not a unique functional of the periodic density change~: it depends also on the change in the macroscopic {\em polarization}. Moreover, the dependence of the exchange-correlation energy on polarization induces an exchange-correlation electric field. These effects are exhibited for a model semiconductor. We also show that the scissor-operator technique is an approximate way of bypassing this polarization dependence., Comment: 11 pages, 1 Fig.

The Born effective charge tensors of Barium Titanate have been calculated for each of its 4 phases. Large effective charges of Ti and O, also predicted by shell model calculations and made plausible by a simplified model, reflect the partial covalent character of the chemical bond. A band by band decomposition confirms that orbital hybridization is not restricted to Ti and O atoms but also involves Ba which appears more covalent than generally assumed. Our calculations reveal a strong dependence of the effective charges on the atomic positions contrasting with a relative insensitivity on isotropic volume changes., Comment: 13 pages