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)

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.

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.

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

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$., 8 pages, 7 figures, ref. 19 corrected

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

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 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

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., 11 pages, 1 table, LaTeX

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., 25 pages, 6 Figures, 10 Tables, LaTeX

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)