Showing 1–50 of 86 results for author: Gonze, X

Polarons in the Cubic Generalized Fröhlich Model: Spontaneous Symmetry Breaking

Authors: Vasilii Vasilchenko, Xavier Gonze

Abstract: Within the variational polaron equation framework, the Fröhlich model for cubic systems with three-fold degenerate electronic bands is numerically solved in the strong coupling regime, for a wide range of its input parameters. By comparing the results to the previously reported ones obtained with the Gaussian Ansatz approach, the inadequacy of the latter is uncovered, especially when degenerate ba… ▽ More Within the variational polaron equation framework, the Fröhlich model for cubic systems with three-fold degenerate electronic bands is numerically solved in the strong coupling regime, for a wide range of its input parameters. By comparing the results to the previously reported ones obtained with the Gaussian Ansatz approach, the inadequacy of the latter is uncovered, especially when degenerate bands are present in a system. Moreover, the symmetry groups of polaronic solutions in the cubic generalized Fröhlich model without spin-orbit coupling are investigated: we provide and discuss a phase diagram of symmetry groups of ground-state polarons, showing spontaneous symmetry breaking. While the cubic symmetry of the three-band degenerate model Hamiltonian corresponds to the full octahedral group $O_h$, lowest-energy polarons possess either $D_{4h}$ or $D_{3d}$ point groups. This phase diagram bears some similarities but differs nevertheless from the one that is obtained by the straight analysis of the band effective masses. The obtained results will provide a firm ground for further exploration of the generalized Fröhlich model and will likely be applicable beyond the model's inherent approximations. △ Less

Submitted 2 May, 2024; originally announced May 2024.

Validation of the GreenX library time-frequency component for efficient GW and RPA calculations

Authors: Maryam Azizi, Jan Wilhelm, Dorothea Golze, Francisco A. Delesma, Ramón L. Panadés-Barrueta, Patrick Rinke, Matteo Giantomassi, Xavier Gonze

Abstract: Electronic structure calculations based on many-body perturbation theory (e.g. GW or the random-phase approximation (RPA)) require function evaluations in the complex time and frequency domain, for example inhomogeneous Fourier transforms or analytic continuation from the imaginary axis to the real axis. For inhomogeneous Fourier transforms, the time-frequency component of the GreenX library provi… ▽ More Electronic structure calculations based on many-body perturbation theory (e.g. GW or the random-phase approximation (RPA)) require function evaluations in the complex time and frequency domain, for example inhomogeneous Fourier transforms or analytic continuation from the imaginary axis to the real axis. For inhomogeneous Fourier transforms, the time-frequency component of the GreenX library provides time-frequency grids that can be utilized in low-scaling RPA and GW implementations. In addition, the adoption of the compact frequency grids provided by our library also reduces the computational overhead in RPA implementations with conventional scaling. In this work, we present low-scaling GW and conventional RPA benchmark calculations using the GreenX grids with different codes (FHI-aims, CP2K and ABINIT) for molecules, two-dimensional materials and solids. Very small integration errors are observed when using 30 time-frequency points for our test cases, namely $<10^{-8}$ eV/electron for the RPA correlation energies, and 10 meV for the GW quasiparticle energies. △ Less

Submitted 12 March, 2024; v1 submitted 11 March, 2024; originally announced March 2024.

Comments: 15 pages, 6 figures

Variational Density Functional Perturbation Theory for Metals

Authors: Xavier Gonze, Samare Rostami, Christian Tantardini

Abstract: Density functional perturbation theory is a well-established method to study responses of molecules and solids, especially responses to atomic displacements or to different perturbing fields (electric, magnetic). Like for density functional theory, the treatment of metals is delicate, due to the Fermi-Dirac statistics and electronic bands crossing the Fermi energy. At zero temperature, there is an… ▽ More Density functional perturbation theory is a well-established method to study responses of molecules and solids, especially responses to atomic displacements or to different perturbing fields (electric, magnetic). Like for density functional theory, the treatment of metals is delicate, due to the Fermi-Dirac statistics and electronic bands crossing the Fermi energy. At zero temperature, there is an abrupt transition from occupied states to unoccupied ones, usually addressed with smearing schemes. Also, at finite temperature, fractional occupations are present, and the occupation numbers may vary in response to the perturbation. The present work establishes the characteristics of density functional perturbation theory stemming from the underlying variational principle, in the case of metals. After briefly reviewing variational density functional theory for metals, the convexity of the entropy function of the occupation number is analyzed, and, at finite temperature, the benefit of resmearing the Fermi-Dirac broadening with the Methfessel-Paxton one is highlighted. Then the variational expressions for the second-order derivative of the free energy are detailed, exposing the different possible gauge choices. The influence of the inaccuracies in the unperturbed wavefunctions from the prior density functional theory calculation is studied. The whole formalism is implemented in the ABINIT software package. △ Less

Submitted 15 February, 2024; originally announced February 2024.

Journal ref: Phys. Rev. B 109, 014317 (2024/01/30)

Surface Enhanced Infrared Absorption mechanism and modification of the plasmonic response

Authors: Tanguy Colleu, Adam Fekete, Xavier Gonze, Alexandre Cloots, Vincent Liégeois, Gian-Marco Rignanese, Luc Henrard

Abstract: Surface Enhanced Infrared Absorption (SEIRA) is an experimental method where trace amount of a compound can be detected with high sensibility. This high detection sensibility is the result of the interaction of the molecules with a localized plasmon, usually from a metallic nano-particle. In this study we numerically investigate by discrete dipole approximation the origin of the Fano-like response… ▽ More Surface Enhanced Infrared Absorption (SEIRA) is an experimental method where trace amount of a compound can be detected with high sensibility. This high detection sensibility is the result of the interaction of the molecules with a localized plasmon, usually from a metallic nano-particle. In this study we numerically investigate by discrete dipole approximation the origin of the Fano-like response of the system, including the induced transparency when the plasmon resonance and the molecular vibrational mode coincide. The detailed analysis of the localization of the absorption show that the modification of the absorption cross-section when the molecule is present comes from a change of the plasmonic resonance, not from the direct molecular response which is negligible. This sheds a new light on the SEIRA mechanism. In particular, it demonstrates that the sensibility is associated with the influence of the molecule on the plasmon resonance rather than with the local field enhancement itself. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: 12 pages, 5 figures

Generating and grading 34 Optimized Norm-Conserving Vanderbilt Pseudopotentials for Actinides and Super Heavy Elements in the PseudoDojo

Authors: Christian Tantardini, Miroslav Iliaš, Matteo Giantomassi, Alexander G. Kvashnin, Valeria Pershina, Xavier Gonze

Abstract: In the last decades, material discovery has been a very active research field driven by the need to find new materials for many different applications. This has also included materials with heavy elements, beyond the stable isotopes of lead, as most actinides exhibit unique properties that make them useful in various applications. Furthermore, new heavy elements beyond actinides, collectively refe… ▽ More In the last decades, material discovery has been a very active research field driven by the need to find new materials for many different applications. This has also included materials with heavy elements, beyond the stable isotopes of lead, as most actinides exhibit unique properties that make them useful in various applications. Furthermore, new heavy elements beyond actinides, collectively referred to as super-heavy elements (SHEs), have been synthesized, filling previously empty space of Mendeleev periodic table. Their chemical bonding behavior, of academic interest at present, would also benefit of state-of-the-art modeling approaches. In particular, in order to perform first-principles calculations with planewave basis sets, one needs corresponding pseudopotentials. In this work, we present a series of scalar- and fully-relativistic optimized norm-conserving Vanderbilt pseudopotentials (ONCVPs) for thirty-four actinides and super-heavy elements, for three different exchange-correlation functionals (PBE, PBEsol and LDA). The scalar-relativistic version of these ONCVPs is tested by comparing equations of states for crystals, obtained with \textsc{abinit} 9.6, with those obtained by all-electron zeroth-order regular approximation (ZORA) calculations, without spin-orbit coupling, performed with the Amsterdam Modeling Suite \textsc{band} code. $Δ$-Gauge and $Δ_1$-Gauge indicators are used to validate these pseudopotentials. This work is a contribution to the PseudoDojo project, in which pseudopotentials for the whole periodic table are developed and systematically tested. The pseudopotential files are available on the PseudoDojo web-interface pseudo-dojo.org in psp8 and UPF2 formats, both suitable for \textsc{abinit}, the latter being also suitable for Quantum ESPRESSO. △ Less

Submitted 7 November, 2023; v1 submitted 6 September, 2023; originally announced September 2023.

Journal ref: Computer Physics Communications, 295, 109002 (2024)

Effect of spin-orbit coupling on the zero-point renormalization of the electronic band gap in cubic materials: First-principles calculations and generalized Fröhlich model

Authors: Véronique Brousseau-Couture, Xavier Gonze, Michel Côté

Abstract: The electronic structure of semiconductors and insulators is affected by ionic motion through electron-phonon interaction, yielding temperature-dependent band gap energies and zero-point renormalization (ZPR) at absolute zero temperature. For polar materials, the most significant contribution to the band gap ZPR can be understood in terms of the Fröhlich model, which focuses on the nonadiabatic in… ▽ More The electronic structure of semiconductors and insulators is affected by ionic motion through electron-phonon interaction, yielding temperature-dependent band gap energies and zero-point renormalization (ZPR) at absolute zero temperature. For polar materials, the most significant contribution to the band gap ZPR can be understood in terms of the Fröhlich model, which focuses on the nonadiabatic interaction between an electron and the macroscopic electrical polarization created by a long-wavelength optical longitudinal phonon mode. On the other hand, spin-orbit interaction (SOC) modifies the bare electronic structure, which will, in turn, affect the electron-phonon interaction and the ZPR. We present a comparative investigation of the effect of SOC on the band gap ZPR of twenty semiconductors and insulators with cubic symmetry using first-principles calculations. We observe a SOC-induced decrease of the ZPR, up to 30%, driven by the valence band edge, which almost entirely originates from the modification of the bare electronic eigenenergies and the decrease of the hole effective masses near the $Γ$ point. We also incorporate SOC into a generalized Fröhlich model, addressing the Dresselhaus splitting which occurs in noncentrosymmetric materials, and confirm that the predominance of nonadiabatic effects on the band gap ZPR of polar materials is unchanged when including SOC. Our generalized Fröhlich model with SOC provides a reliable estimate of the SOC-induced decrease of the polaron formation energy obtained from first principles and brings to light some fundamental subtleties in the numerical evaluation of the effective masses with SOC for noncentrosymmetric materials. We finally warn about a possible breakdown of the parabolic approximation within the physically relevant energy range of the Fröhlich interaction for materials with high phonon frequencies treated with SOC. △ Less

Submitted 11 July, 2023; originally announced July 2023.

Comments: 34 pages including Supplemental Material; 10 figures, 1 table, 3 Supplementary figures, 9 Supplementary tables

Journal ref: Phys. Rev. B 107, 115173 (2023)

A First-Principles Explanation of the Luminescent Line Shape of SrLiAl$_3$N$_4$:Eu$^{2+}$ Phosphor for Light-Emitting Diode Applications

Authors: Julien Bouquiaux, Samuel Poncé, Yongchao Jia, Anna Miglio, Masayoshi Mikami, Xavier Gonze

Abstract: White light-emitting diodes are gaining popularity and are set to become the most common light source in the U.S. by 2025. However, their performance is still limited by the lack of an efficient red-emitting component with a narrow band emission. The red phosphor SrLiAl$_3$N$_4$:Eu$^{2+}$ is among the first promising phosphors with a small bandwidth for next-generation lighting, but the microscopi… ▽ More White light-emitting diodes are gaining popularity and are set to become the most common light source in the U.S. by 2025. However, their performance is still limited by the lack of an efficient red-emitting component with a narrow band emission. The red phosphor SrLiAl$_3$N$_4$:Eu$^{2+}$ is among the first promising phosphors with a small bandwidth for next-generation lighting, but the microscopic origin of this narrow emission remains elusive. In the present work, density functional theory, the $Δ$SCF-constrained occupation method, and a generalized Huang-Rhys theory are used to provide an accurate description of the vibronic processes occurring at the two Sr$^{2+}$ sites that the Eu$^{2+}$ activator can occupy. The emission band shape of Eu(Sr1), with a zero-phonon line at 1.906 eV and a high luminescence intensity, is shown to be controlled by the coupling between the 5d$_{z^2}$-4f electronic transition and the low-frequency phonon modes associated with the Sr and Eu displacements along the Sr channel. The good agreement between our computations and experimental results allows us to provide a structural assignment of the observed total spectrum. By computing explicitly the effect of the thermal expansion on zero-phonon line energies, the agreement is extended to the temperature-dependent spectrum. These results provide insight into the electron-phonon coupling that accompanies the 5d-4f transition in similar UCr$_4$C$_4$-type phosphors. Furthermore, these results highlight the importance of the Sr channel in shaping the narrow emission of SrLiAl$_3$N$_4$:Eu$^{2+}$, and they shed new light on the structure-property relations of such phosphors. △ Less

Submitted 7 July, 2023; v1 submitted 15 June, 2023; originally announced June 2023.

Comments: 10 pages, 6 figures

Journal ref: Chem. Mater. 2023

Roadmap on Electronic Structure Codes in the Exascale Era

Authors: Vikram Gavini, Stefano Baroni, Volker Blum, David R. Bowler, Alexander Buccheri, James R. Chelikowsky, Sambit Das, William Dawson, Pietro Delugas, Mehmet Dogan, Claudia Draxl, Giulia Galli, Luigi Genovese, Paolo Giannozzi, Matteo Giantomassi, Xavier Gonze, Marco Govoni, Andris Gulans, François Gygi, John M. Herbert, Sebastian Kokott, Thomas D. Kühne, Kai-Hsin Liou, Tsuyoshi Miyazaki, Phani Motamarri , et al. (16 additional authors not shown)

Abstract: Electronic structure calculations have been instrumental in providing many important insights into a range of physical and chemical properties of various molecular and solid-state systems. Their importance to various fields, including materials science, chemical sciences, computational chemistry and device physics, is underscored by the large fraction of available public supercomputing resources d… ▽ More Electronic structure calculations have been instrumental in providing many important insights into a range of physical and chemical properties of various molecular and solid-state systems. Their importance to various fields, including materials science, chemical sciences, computational chemistry and device physics, is underscored by the large fraction of available public supercomputing resources devoted to these calculations. As we enter the exascale era, exciting new opportunities to increase simulation numbers, sizes, and accuracies present themselves. In order to realize these promises, the community of electronic structure software developers will however first have to tackle a number of challenges pertaining to the efficient use of new architectures that will rely heavily on massive parallelism and hardware accelerators. This roadmap provides a broad overview of the state-of-the-art in electronic structure calculations and of the various new directions being pursued by the community. It covers 14 electronic structure codes, presenting their current status, their development priorities over the next five years, and their plans towards tackling the challenges and leveraging the opportunities presented by the advent of exascale computing. △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: Submitted as a roadmap article to Modelling and Simulation in Materials Science and Engineering; Address any correspondence to Vikram Gavini (vikramg@umich.edu) and Danny Perez (danny_perez@lanl.gov)

Zero-point lattice expansion and band gap renormalization: Grüneisen approach versus free energy minimization

Authors: Véronique Brousseau-Couture, Émile Godbout, Xavier Gonze, Michel Côté

Abstract: The zero-point lattice expansion (ZPLE) is a small variation of the lattice parameters induced by the presence of phonons in a material compared to the static lattice picture. It contributes significantly to the zero-point renormalization (ZPR) of the band gap energy, but its consequences have not been investigated as thoroughly as those stemming from electron-phonon interactions. In the usual fir… ▽ More The zero-point lattice expansion (ZPLE) is a small variation of the lattice parameters induced by the presence of phonons in a material compared to the static lattice picture. It contributes significantly to the zero-point renormalization (ZPR) of the band gap energy, but its consequences have not been investigated as thoroughly as those stemming from electron-phonon interactions. In the usual first-principles approach, one evaluates the ZPLE by minimizing the $T=0$ K Helmholtz free energy. In this work, we show that the formalism based on the Grüneisen parameters, which commonly neglects zero-point effects, can be efficiently used to compute ZPLE for both isotropic and anisotropic materials at much lower computational cost. We systematically test this formalism on 22 cubic and wurtzite materials and obtain excellent agreement with free energy minimization results for both the ZPLE and the resulting band gap ZPR. We use our results to validate an empirical expression estimating the ZPLE-induced ZPR and unveil its sensitivity to the temperature range involved in estimating the ZPLE from experimental data. Our findings finally reveal that the ZPLE contribution to the band gap ZPR can reach 20% to more than 80% of the electron-phonon interaction contribution for heavier or more ionic materials, including materials containing light atoms. Considering both contributions on an equal footing is thus essential should one attempt to compare theoretical ZPR results with experimental data. △ Less

Submitted 25 August, 2022; v1 submitted 6 July, 2022; originally announced July 2022.

Comments: 30 pages, 5 figures, 14 tables. Sec. II D expanded for internal consistency, typo corrected in Fig.4., clarified Tables captions, minor corrections to the text; results and conclusions unchanged

Journal ref: Phys. Rev. B 106, 085137 (2022)

High-throughput analysis of Fröhlich-type polaron models

Authors: Pedro Miguel M. C. de Melo, Joao C. de Abreu, Bogdan Guster, Matteo Giantomassi, Zeila Zanolli, Xavier Gonze, Matthieu J. Verstraete

Abstract: The electronic structure of condensed matter can be significantly affected by the electron-phonon interaction, leading to important phenomena such as electrical resistance, superconductivity or the formation of polarons. This interaction is often neglected in band structure calculations but can have a strong impact on band gaps or optical spectra. Commonly used frameworks for electron-phonon energ… ▽ More The electronic structure of condensed matter can be significantly affected by the electron-phonon interaction, leading to important phenomena such as electrical resistance, superconductivity or the formation of polarons. This interaction is often neglected in band structure calculations but can have a strong impact on band gaps or optical spectra. Commonly used frameworks for electron-phonon energy corrections are the Allen-Heine-Cardona theory and the Fröhlich model. While the latter shows qualitative agreement with experiment for many polar materials, its simplicity should bring hard limits to its applicability in real materials. Improvements can be made by introducing a generalized version of the model, which considers anisotropic and degenerate electronic bands, and multiple phonon branches. In this work, we search for trends and outliers on over a thousand materials in existing databases of phonon and electron band structures. We use our results to identify the limits of applicability of the standard Frölich model by comparing to the generalized version, and by testing its basic hypothesis of a large radius for the polaronic wavefunction and the corresponding atomic displacement cloud. Among our extended set of materials, most exhibit large polaron behavior as well as validity of the perturbative treatment. For the valence band, there is also a significant fraction of the materials for which the perturbative treatment cannot be applied and/or for which the size of the self-trapping region is close to the atomic repetition distance. We find a large variety of behaviors, and employ much more accurate, fully ab initio Allen-Heine-Cardona calculations to understand extreme cases, where the Fröhlich model should fail and unusually large zero-point renormalization energies occur. △ Less

Submitted 1 July, 2022; originally announced July 2022.

Variational Polaron Equations Applied to the Anisotropic Fröhlich Model

Authors: Vasilii Vasilchenko, Andriy Zhugayevych, Xavier Gonze

Abstract: Starting from recent advances in the first-principles modeling of polarons, variational polaron equations in the strong-coupling adiabatic approximation are formulated in Bloch space. In this framework, polaron formation energy as well as individual electron, phonon and electron-phonon contributions are obtained. We suggest an efficient gradient-based optimization algorithm and apply these equatio… ▽ More Starting from recent advances in the first-principles modeling of polarons, variational polaron equations in the strong-coupling adiabatic approximation are formulated in Bloch space. In this framework, polaron formation energy as well as individual electron, phonon and electron-phonon contributions are obtained. We suggest an efficient gradient-based optimization algorithm and apply these equations to the generalized Fröhlich model with anisotropic non-degenerate electronic bands, both in two- and three-dimensional cases. The effect of the divergence of Fröhlich electron-phonon matrix elements at $Γ$-point is treated analytically, improving the convergence with respect to the sampling in reciprocal space. The whole methodology is validated by obtaining the known asymptotic solution of the standard Fröhlich model in isotropic scenario and also by comparing our results with the Gaussian ansatz approach, showing the difference between the numerically exact and Gaussian trial wavefunctions. Additionally, decomposition of the energy into individual terms allows one to recover the Pekar's 1:2:3:4 theorem, which is shown to be valid even in the anisotropic case. We expect that the improvements in the formalism and numerical implementation will be applicable beyond the large polaron hypothesis inherent to Fröhlich model. △ Less

Submitted 2 June, 2022; originally announced June 2022.

Comments: 14 pages, 5 figures

Spectroscopic signatures of nonpolarons : the case of diamond

Authors: Joao C. de Abreu, Jean Paul Nery, Matteo Giantomassi, Xavier Gonze, Matthieu J. Verstraete

Abstract: Polarons are quasi-particles made from electrons interacting with vibrations in crystal lattices. They derive their name from the strong electron-vibration polar interaction in ionic systems, that induces associated spectroscopic and optical signatures of such quasi-particles in these materials. In this paper, we focus on diamond, a non-polar crystal with inversion symmetry which nevertheless show… ▽ More Polarons are quasi-particles made from electrons interacting with vibrations in crystal lattices. They derive their name from the strong electron-vibration polar interaction in ionic systems, that induces associated spectroscopic and optical signatures of such quasi-particles in these materials. In this paper, we focus on diamond, a non-polar crystal with inversion symmetry which nevertheless shows characteristic signatures of polarons, better denoted "nonpolarons" in this case. The polaronic effects are produced by short-range crystal fields with only a small influence of long-range quadrupoles. The many-body spectral function has a characteristic energy dependence, showing a plateau structure that is similar to but distinct from the satellites observed in the polar Fröhlich case. The temperature-dependent spectral function of diamond is determined by two methods: the standard Dyson-Migdal approach, which calculates electron-phonon interactions within the lowest-order expansion of the self-energy, and the cumulant expansion, which includes higher orders of electron-phonon interactions. The latter corrects the nonpolaron energies and broadening, providing a more realistic spectral function, which we examine in detail for both conduction and valence band edges. △ Less

Submitted 25 March, 2022; originally announced March 2022.

Comments: 10 pages 16 Figures + Subfigures 2 Tables 82 References Supplementary Information (SI) at the end: SI 10 pages SI 14 Figures + Subfigures SI 1 Table

Sr-Doped Molecular Hydrogen: Synthesis and Properties of SrH$_{22}$

Authors: Dmitrii V. Semenok, Wuhao Chen, Xiaoli Huang, Di Zhou, Ivan A. Kruglov, Arslan B. Mazitov, Michele Galasso, Christian Tantardini, Xavier Gonze, Alexander G. Kvashnin, Artem R. Oganov, Tian Cui

Abstract: Recently, several research groups announced reaching the point of metallization of hydrogen above 400 GPa. Following the mainstream of extensive investigations of compressed polyhydrides, in this work we demonstrate that small (4 atom %) doping of molecular hydrogen by strontium leads to a dramatic reduction in the metallization pressure to about 200 GPa. Studying the high-pressure chemistry of th… ▽ More Recently, several research groups announced reaching the point of metallization of hydrogen above 400 GPa. Following the mainstream of extensive investigations of compressed polyhydrides, in this work we demonstrate that small (4 atom %) doping of molecular hydrogen by strontium leads to a dramatic reduction in the metallization pressure to about 200 GPa. Studying the high-pressure chemistry of the Sr-H system at 56-180 GPa, we observed the formation of several previously unknown compounds: C2/m-Sr$_3$H$_{13}$, pseudocubic SrH$_6$, SrH$_9$ with cubic F-43m Sr sublattice, and pseudotetragonal P1-SrH$_{22}$, the metal hydride with the highest hydrogen content discovered so far. Unlike Ca and Y, strontium forms molecular semiconducting polyhydrides, whereas calcium and yttrium polyhydrides are high-Tc superconductors with an atomic H sublattice. The latter phase, SrH$_{22}$ or Sr$_{0.04}$H$_{0.96}$, may be considered as a convenient model of the consistent bandgap closure and metallization of hydrogen. Using the impedance measurements in diamond anvil cells at 300-440 K, we estimated the direct bandgap of the Pm-3n-like compound P1-SrH$_6$ to be 0.44-0.51 eV at 150 GPa, and its metallization pressure to be 220 GPa. Together with the machine learning interatomic potentials, the impedance spectroscopy allowed us to estimate the diffusion coefficients of hydrogen D$_H$ = 1.0-2.8 E-10 m$^2$/s in SrH$_6$ and 1.2-2.1 E-9 m$^2$/s in P1-SrH$_{22}$ at 500-600 K. △ Less

Submitted 5 November, 2021; v1 submitted 29 October, 2021; originally announced October 2021.

Comments: Supporting information was compressed and reduced by 2 times to 36 figures

Journal ref: Adv Mater . 2022 Apr 21; e2200924

Fröhlich polaron effective mass and localization length in cubic materials: degenerate and anisotropic electronic bands

Authors: Bogdan Guster, Pedro Melo, Bradley A. A. Martin, Véronique Brousseau-Couture, Anna Miglio, Matteo Giantomassi, Michel Côté, Jarvist M. Frost, Matthieu J. Verstraete, Xavier Gonze

Abstract: Polarons, that is, charge carriers correlated with lattice deformations, are ubiquitous quasiparticles in semiconductors, and play an important role in electrical conductivity. To date most theoretical studies of so-called large polarons, in which the lattice can be considered as a continuum, have focused on the original Fröhlich model: a simple (non-degenerate) parabolic isotropic electronic band… ▽ More Polarons, that is, charge carriers correlated with lattice deformations, are ubiquitous quasiparticles in semiconductors, and play an important role in electrical conductivity. To date most theoretical studies of so-called large polarons, in which the lattice can be considered as a continuum, have focused on the original Fröhlich model: a simple (non-degenerate) parabolic isotropic electronic band coupled to one dispersionless longitudinal optical phonon branch. The Fröhlich model allows one to understand characteristics such as polaron formation energy, radius, effective mass and mobility. Real cubic materials, instead, have electronic band extrema that are often degenerate or anisotropic and present several phonon modes. In the present work, we address such issues. We keep the continuum hypothesis inherent to the large polaron Fröhlich model, but waive the isotropic and non-degeneracy hypotheses, and also include multiple phonon branches. For polaron effective masses, working at the lowest order of perturbation theory, we provide analytical results for the case of anisotropic electronic energy dispersion, with two distinct effective masses (uniaxial) and numerical simulations for the degenerate 3-band case, typical of III-V and II-VI semiconductor valence bands. We also deal with the strong-coupling limit, using a variational treatment: we propose trial wavefunctions for the above-mentioned cases, providing polaron radii and energies. Then, we evaluate the polaron formation energies, effective masses and localisation lengths using parameters representative of a dozen II-VI, III-V and oxide semiconductors, for both electron and hole polarons...In the non-degenerate case, we compare the perturbative approach with the Feynman path integral approach in characterisizing polarons in the weak coupling limit... △ Less

Submitted 26 September, 2021; originally announced September 2021.

Comments: 18 pages, 10 figures in main text; 11 pages, 1 figure in supplementary material

Design rule for the emission linewidth of Eu$^{2+}$-activated phosphors

Authors: Yongchao Jia, Samuel Poncé, Anna Miglio, Masayoshi Mikami, Xavier Gonze

Abstract: We study from first principles the emission linewidth of Eu$^{2+}$-doped LED phosphors. Based on the one-dimensional configuration coordinate model, an analysis of first principles data obtained for fifteen compounds show that, at working temperature, the linewidth of Eu$^{2+}$ emission band in solids is negligibly affected by quantum effects, and can be extracted from the Franck-Condon energy shi… ▽ More We study from first principles the emission linewidth of Eu$^{2+}$-doped LED phosphors. Based on the one-dimensional configuration coordinate model, an analysis of first principles data obtained for fifteen compounds show that, at working temperature, the linewidth of Eu$^{2+}$ emission band in solids is negligibly affected by quantum effects, and can be extracted from the Franck-Condon energy shifts. For a fixed Stokes shift, the difference of Franck-Condon energy shifts in the excited and ground states is the key factor for the FWHM determination. Narrow emission Eu$^{2+}$-doped LED phosphors are expected for the case with large positive value of such difference. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: 6 pages, 4 figures

CZTS Raman spectra beyond kesterite: a first-principles study

Authors: S. P. Ramkumar, H. P. C. Miranda, X. Gonze, G. -M. Rignanese

Abstract: 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… ▽ More 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. △ Less

Submitted 12 January, 2021; originally announced January 2021.

Comments: 16 pages, 5 figures (+7 figures in the appendix)

Predominance of non-adiabatic effects in zero-point renormalization of the electronic band gap

Authors: Anna Miglio, Véronique Brousseau-Couture, Emile Godbout, Gabriel Antonius, Yang-Hao Chan, Steven G. Louie, Michel Côté, Matteo Giantomassi, Xavier Gonze

Abstract: Electronic and optical properties of materials are affected by atomic motion through the electron-phonon interaction: not only band gaps change with temperature, but even at absolute zero temperature, zero-point motion causes band-gap renormalization. We present a large-scale first-principles evaluation of the zero-point renormalization of band edges beyond the adiabatic approximation. For materia… ▽ More Electronic and optical properties of materials are affected by atomic motion through the electron-phonon interaction: not only band gaps change with temperature, but even at absolute zero temperature, zero-point motion causes band-gap renormalization. We present a large-scale first-principles evaluation of the zero-point renormalization of band edges beyond the adiabatic approximation. For materials with light elements, the band gap renormalization is often larger than 0.3 eV, and up to 0.7 eV. This effect cannot be ignored if accurate band gaps are sought. For infrared-active materials, global agreement with available experimental data is obtained only when non-adiabatic effects are taken into account. They even dominate zero-point renormalization for many materials, as shown by a generalized Fröhlich model that includes multiple phonon branches, anisotropic and degenerate electronic extrema, whose range of validity is established by comparison with first-principles results. △ Less

Submitted 25 November, 2020; originally announced November 2020.

Comments: Note that the supplementary material, included in the main file, is structured differently than in the published version, with a different numbering scheme

Journal ref: npj Computational Materials 6, 167 (2020)

Importance of long-range channel Sr displacements for the narrow emission in Sr[Li$_2$Al$_2$O$_2$N$_2$]:Eu$^{2+}$ phosphor

Authors: Julien Bouquiaux, Samuel Poncé, Yongchao Jia, Anna Miglio, Masayochi Mikami, Xavier Gonze

Abstract: The recently discovered Sr[Li$_2$Al$_2$O$_2$N$_2$]:Eu$^{2+}$ red phosphor, candidate for the next generation of eco-efficient white light-emitting diodes, exhibits excellent emission spectral position and exceptionally small linewidth. It belongs to the UCr$_4$C$_4$-structure family of phosphors, containing many potential candidates commercial phosphors, whose small linewidth, tentatively ascribed… ▽ More The recently discovered Sr[Li$_2$Al$_2$O$_2$N$_2$]:Eu$^{2+}$ red phosphor, candidate for the next generation of eco-efficient white light-emitting diodes, exhibits excellent emission spectral position and exceptionally small linewidth. It belongs to the UCr$_4$C$_4$-structure family of phosphors, containing many potential candidates commercial phosphors, whose small linewidth, tentatively ascribed to the high-symmetry cuboid environment of the doping site, has drawn the attention of researchers in the last five years. We use density functional theory, $Δ$SCF method and configuration coordinate models (CCM) to provide a complete characterization of this material. Using a multi-dimensional CCM, an accurate description of the coupling of the vibronic structure with the electronic 5d$\rightarrow$4f transition is obtained, including the partial Huang-Rhys factors and frequency of the dominant modes. We show that, in addition to the first-coordination shell cuboid deformation mode, low-frequency phonon modes involving chains of strontium atoms along the tetragonal axis shape the emission linewidth in Sr[Li$_2$Al$_2$O$_2$N$_2$]:Eu$^{2+}$. This finding sheds new light on the emission properties of UCr$_4$C$_4$-structure phosphors, possessing similar Ca/Sr/Ba channel. Our approach provides a robust theoretical framework to systematically study the emission spectra of such Eu-doped phosphors, and predict candidates with expected similar or even sharper linewidth. △ Less

Submitted 9 September, 2021; v1 submitted 1 October, 2020; originally announced October 2020.

Comments: Published version

Journal ref: Adv. Optical Mater. 2021, 2100649

Band gap renormalization, carrier mobilities, and the electron-phonon self-energy in crystalline naphthalene

Authors: Florian Brown-Altvater, Gabriel Antonius, Tonatiuh Rangel, Matteo Giantomassi, Claudia Draxl, Xavier Gonze, Steven G. Louie, Jeffrey B. Neaton

Abstract: Organic molecular crystals are expected to feature appreciable electron-phonon interactions that influence their electronic properties at zero and finite temperature. In this work, we report first-principles calculations and an analysis of the electron-phonon self-energy in naphthalene crystals. We compute the zero-point renormalization and temperature dependence of the fundamental band gap, and t… ▽ More Organic molecular crystals are expected to feature appreciable electron-phonon interactions that influence their electronic properties at zero and finite temperature. In this work, we report first-principles calculations and an analysis of the electron-phonon self-energy in naphthalene crystals. We compute the zero-point renormalization and temperature dependence of the fundamental band gap, and the resulting scattering lifetimes of electronic states near the valence- and conduction-band edges employing density functional theory. Further, our calculated phonon renormalization of the $GW$-corrected quasiparticle band structure predicts a fundamental band gap of 5 eV for naphthalene at room temperature, in good agreement with experiments. From our calculated phonon-induced electron lifetimes, we obtain the temperature-dependent mobilities of electrons and holes in good agreement with experimental measurements at room temperatures. Finally, we show that an approximate energy self-consistent computational scheme for the electron-phonon self-energy leads to the prediction of strong satellite bands in the electronic band structure. We find that a single calculation of the self-energy can reproduce the self-consistent results of the band gap renormalization and electrical mobilities for naphthalene, provided that the on-the-mass-shell approximation is used, i.e., if the self-energy is evaluated at the bare eigenvalues. △ Less

Submitted 19 July, 2020; originally announced July 2020.

Comments: 12 pages, 7 figures, 3 tables

Journal ref: Phys. Rev. B 101 (2020) 165102

Phonon-limited electron mobility in Si, GaAs and GaP with exact treatment of dynamical quadrupoles

Authors: Guillaume Brunin, Henrique Pereira Coutada Miranda, Matteo Giantomassi, Miquel Royo, Massimiliano Stengel, Matthieu J. Verstraete, Xavier Gonze, Gian-Marco Rignanese, Geoffroy Hautier

Abstract: We describe a new approach to compute the electron-phonon self-energy and carrier mobilities in semiconductors. Our implementation does not require a localized basis set to interpolate the electron-phonon matrix elements, with the advantage that computations can be easily automated. Scattering potentials are interpolated on dense $\mathbf{q}$ meshes using Fourier transforms and ab initio models to… ▽ More We describe a new approach to compute the electron-phonon self-energy and carrier mobilities in semiconductors. Our implementation does not require a localized basis set to interpolate the electron-phonon matrix elements, with the advantage that computations can be easily automated. Scattering potentials are interpolated on dense $\mathbf{q}$ meshes using Fourier transforms and ab initio models to describe the long-range potentials generated by dipoles and quadrupoles. To reduce significantly the computational cost, we take advantage of crystal symmetries and employ the linear tetrahedron method and double-grid integration schemes, in conjunction with filtering techniques in the Brillouin zone. We report results for the electron mobility in Si, GaAs, and GaP obtained with this new methodology. △ Less

Submitted 21 September, 2020; v1 submitted 3 February, 2020; originally announced February 2020.

Journal ref: Phys. Rev. B 102, 094308 (2020)

Electron-Phonon Beyond Fröhlich: Dynamical Quadrupoles in Polar and Covalent Solids

Authors: Guillaume Brunin, Henrique Pereira Coutada Miranda, Matteo Giantomassi, Miquel Royo, Massimiliano Stengel, Matthieu J. Verstraete, Xavier Gonze, Gian-Marco Rignanese, Geoffroy Hautier

Abstract: We include the treatment of quadrupolar fields beyond the Fröhlich interaction in the first-principles electron-phonon vertex in semiconductors. Such quadrupolar fields induce long-range interactions that have to be taken into account for accurate physical results. We apply our formalism to Si (nonpolar), GaAs, and GaP (polar) and demonstrate that electron mobilities show large errors if dynamical… ▽ More We include the treatment of quadrupolar fields beyond the Fröhlich interaction in the first-principles electron-phonon vertex in semiconductors. Such quadrupolar fields induce long-range interactions that have to be taken into account for accurate physical results. We apply our formalism to Si (nonpolar), GaAs, and GaP (polar) and demonstrate that electron mobilities show large errors if dynamical quadrupoles are not properly treated. △ Less

Submitted 21 September, 2020; v1 submitted 3 February, 2020; originally announced February 2020.

Journal ref: Phys. Rev. Lett. 125, 136601 (2020)

First-principles prediction of lattice coherency in van der Waals heterostructures

Authors: Benoit Van Troeye, Aurélien Lherbier, Simon M. -M. Dubois, Jean-Christophe Charlier, Xavier Gonze

Abstract: The emergence of superconductivity in slightly-misaligned graphene bilayer [1] and moiré excitons in MoSe$_2$-WSe$_2$ van der Waals (vdW) heterostructures [2] is intimately related to the formation of a 2D superlattice in those systems. At variance, perfect primitive lattice matching of the constituent layers has also been reported in some vdW-heterostructures [3-5], highlighting the richness of i… ▽ More The emergence of superconductivity in slightly-misaligned graphene bilayer [1] and moiré excitons in MoSe$_2$-WSe$_2$ van der Waals (vdW) heterostructures [2] is intimately related to the formation of a 2D superlattice in those systems. At variance, perfect primitive lattice matching of the constituent layers has also been reported in some vdW-heterostructures [3-5], highlighting the richness of interfaces in the 2D world. In this work, the determination of the nature of such interface, from first principles, is demonstrated. To do so, an extension of the Frenkel-Kontorova (FK) model [6] is presented, linked to first-principles calculations, and used to predict lattice coherency for a set of 56 vdW-heterostructures. Computational predictions agree with experiments, when available. New superlattices as well as perfectly-matching interfaces are predicted. △ Less

Submitted 15 November, 2019; originally announced November 2019.

Comments: 16 pages, 3 figures

Beyond the one-dimensional configuration coordinate model of photoluminescence

Authors: Yongchao Jia, Samuel Poncé, Anna Miglio, Masayoshi Mikami, Xavier Gonze

Abstract: The one-dimensional configuration coordinate model (1D-CCM) is widely used for the analysis of photoluminescence in molecules and doped solids, and relies on a linear combination of the equilibrium nuclear configurations of ground and excited states. It delivers an estimation of the energy barrier at which ground and excited state curves cross, semi-classically linked to non-radiative transition r… ▽ More The one-dimensional configuration coordinate model (1D-CCM) is widely used for the analysis of photoluminescence in molecules and doped solids, and relies on a linear combination of the equilibrium nuclear configurations of ground and excited states. It delivers an estimation of the energy barrier at which ground and excited state curves cross, semi-classically linked to non-radiative transition rate and thermal quenching. To assess its predictive power for the latter properties, we propose a new \textit{optimized configuration path (OCP) method in which} the ground-state and excited-state forces are mixed instead of their configurations. We also define another one-parameter model thanks a double energy parabola hypothesis (DEPH). We compare the OCP method and the DEPH reference with the 1D-CCM for three paradigmatic 4f-5d phosphors Y$_3$Al$_5$O$_{12}$:Ce, Lu$_2$SiO$_5$:Ce, and YAlO$_3$:Ce. We find that the OCP and DEPH methods yield similar results with geometries that have significantly lower ground-state energies than the 1D-CCM for the same 4f-5d energy difference. However, the OCP method suffers from the appearance of multiple local minima, rendering the clear determination of the optimal geometry very difficult in practice. Still the OCP method allows one to quantify the deviations from the 1D-CCM, therefore increasing confidence in the lower bound obtained from the DEPH for the 4f-5d crossing barrier, and its comparison with the energy of the auto-ionization thermal quenching mechanism. We expect the OCP approach to be applicable to other luminescent materials or molecules. △ Less

Submitted 20 September, 2019; originally announced September 2019.

Comments: 12 pages, 6 figure (Accepted by Phys. Rev. B)

Vibrational and dielectric properties of monolayer transition metal dichalcogenides

Authors: Nicholas A. Pike, Antoine Dewandre, Benoit Van Troeye, Xavier Gonze, Matthieu J. Verstraete

Abstract: First-principles studies of two-dimensional transition metal dichalcogenides have contributed considerably to the understanding of their dielectric, optical, elastic, and vibrational properties. The majority of works to date focus on a single material or physical property. Here we use a single first-principles methodology on the whole family of systems, to investigate in depth the relationships be… ▽ More First-principles studies of two-dimensional transition metal dichalcogenides have contributed considerably to the understanding of their dielectric, optical, elastic, and vibrational properties. The majority of works to date focus on a single material or physical property. Here we use a single first-principles methodology on the whole family of systems, to investigate in depth the relationships between different physical properties, the underlying symmetry and the composition of these materials, and observe trends. We compare to bulk counterparts to show strong interlayer effects in triclinic compounds. Previously unobserved relationships between these monolayer compounds become apparent. These trends can then be exploited by the materials science, nanoscience, and chemistry communities to better design devices and heterostructures for specific functionalities. △ Less

Submitted 25 March, 2019; originally announced March 2019.

Comments: 4 figures, 11 pages

Journal ref: Phys. Rev. Materials 3, 074009 (2019)

Ab initio Study of Luminescence in Ce-doped Lu$_2$SiO$_5$: The Role of Oxygen Vacancies on Emission Color and Thermal Quenching Behavior

Authors: Yongchao Jia, Anna Miglio, Masayoshi Mikami, Xavier Gonze

Abstract: We study from first principles the luminescence of Lu$_2$SiO$_5$:Ce$^{3+}$ (LSO:Ce), a scintillator widely used in medical imaging applications, and establish the crucial role of oxygen vacancies (V$_O$) in the generated spectrum. The excitation energy, emission energy and Stokes shift of its luminescent centers are simulated through a constrained density-functional theory method coupled with a… ▽ More We study from first principles the luminescence of Lu$_2$SiO$_5$:Ce$^{3+}$ (LSO:Ce), a scintillator widely used in medical imaging applications, and establish the crucial role of oxygen vacancies (V$_O$) in the generated spectrum. The excitation energy, emission energy and Stokes shift of its luminescent centers are simulated through a constrained density-functional theory method coupled with a $Δ$SCF analysis of total energies, and compared with experimental spectra. We show that the high-energy emission band comes from a single Ce-based luminescent center, while the large experimental spread of the low-energy emission band originates from a whole set of different Ce-V$_O$ complexes together with the other Ce-based luminescent center. Further, the luminescence thermal quenching behavior is analyzed. The $4f-5d$ crossover mechanism is found to be very unlikely, with a large crossing energy barrier (E$_{fd}$) in the one-dimensional model. The alternative mechanism usually considered, namely the electron auto-ionization, is also shown to be unlikely. In this respect, we introduce a new methodology in which the time-consuming accurate computation of the band gap for such models is bypassed. We emphasize the usually overlooked role of the differing geometry relaxation in the excited neutral electronic state Ce$^{3+,*}$ and in the ionized electronic state Ce$^{4+}$. The results indicate that such electron auto-ionization cannot explain the thermal stability difference between the high- and low-energy emission bands. Finally, a hole auto-ionization process is proposed as a plausible alternative. With the already well-established excited state characterization methodology, the approach to color center identification and thermal quenching analysis proposed here can be applied to other luminescent materials in the presence of intrinsic defects. △ Less

Submitted 16 November, 2018; originally announced November 2018.

Comments: 13 pages, 8 figures, accepted by Phys. Rev. Materials

Journal ref: Physical Review MATERIALS 2, 125202 (2018)

Vibrational and dielectric properties of the bulk transition metal dichalcogenides

Authors: Nicholas A. Pike, Antoine Dewandre, Benoit Van Troeye, Xavier Gonze, Matthieu J. Verstraete

Abstract: Interest in the bulk transition metal dichalcogenides for their electronic, photovoltaic, and optical properties has grown and led to their use in many technological applications. We present a systematic investigation of their interlinked vibrational and dielectric properties, using density functional theory and density functional perturbation theory, studying the effects of the spin-orbit interac… ▽ More Interest in the bulk transition metal dichalcogenides for their electronic, photovoltaic, and optical properties has grown and led to their use in many technological applications. We present a systematic investigation of their interlinked vibrational and dielectric properties, using density functional theory and density functional perturbation theory, studying the effects of the spin-orbit interaction and of the long-range e$^-$- e$^-$ correlation as part of our investigation. This study confirms that the spin-orbit interaction plays a small role in these physical properties, while the direct contribution of dispersion corrections is of crucial importance in the description of the interatomic force constants. Here, our analysis of the structural and vibrational properties, including the Raman spectra, compare well to experimental measurement. Three materials with different point groups are showcased and data trends on the full set of fifteen existing hexagonal, trigonal, and triclinic materials are demonstrated. This overall picture will enable the modeling of devices composed of these materials for novel applications. △ Less

Submitted 6 April, 2018; originally announced April 2018.

Comments: 11 pages, 6 figures

Journal ref: Phys. Rev. Materials 2, 063608 (2018)

First-principle study of paraelectric and ferroelectric CsH$_2$PO$_4$ including dispersion forces: stability and related vibrational, dielectric and elastic properties

Authors: Benoit Van Troeye, Michiel Jan van Setten, Matteo Giantomassi, Marc Torrent, Gian-Marco Rignanese, Xavier Gonze

Abstract: Using density functional theory (DFT) and density functional perturbation theory (DFPT), we investigate the stability and response functions of CsH$_2$PO$_4$, a ferroelectric material at low temperature. This material cannot be described properly by the usual (semi-)local approximations within DFT. The long-range e$^-$-e$^-$ correlation needs to be properly taken into account, using, for instance,… ▽ More Using density functional theory (DFT) and density functional perturbation theory (DFPT), we investigate the stability and response functions of CsH$_2$PO$_4$, a ferroelectric material at low temperature. This material cannot be described properly by the usual (semi-)local approximations within DFT. The long-range e$^-$-e$^-$ correlation needs to be properly taken into account, using, for instance, Grimme's DFT-D methods, as investigated in this work. We find that DFT-D3(BJ) performs the best for the members of the dihydrogenated alkali phosphate family (KH$_2$PO$_4$, RbH$_2$PO$_4$, CsH$_2$PO$_4$), leading to experimental lattice parameters reproduced with an average deviation of 0.5 %. With these DFT-D methods, the structural, dielectric, vibrational and mechanical properties of CsH$_2$PO$_4$ are globally in excellent agreement with the available experiments ($<$ 2% MAPE for Raman-active phonons). Our study suggests the possible existence of a new low-temperature phase for CsH$_2$PO$_4$, not yet reported experimentally. Finally, we report the implementation of DFT-D contributions to elastic constants within DFPT. △ Less

Submitted 26 January, 2018; originally announced January 2018.

Comments: This paper was published in Physical Review B the 25 January 2017 (21 pages, 4 figures)

Journal ref: B. Van Troeye, M. J. van Setten, M. Giantomassi, M. Torrent, G.-M. Rignanese, and X. Gonze, Phys. Rev. B 95, 024112 (2017)

Interatomic force constants including the DFT-D dispersion contribution

Authors: Benoit Van Troeye, Marc Torrent, Xavier Gonze

Abstract: Grimme's DFT-D dispersion contribution to interatomic forces constants, required for the computation of the phonon band structures in density-functional perturbation theory, has been derived analytically. The implementation has then been validated with respect to frozen phonons, and applied on materials where weak cohesive forces play a major role i.e. argon, graphite, benzene. We show that these… ▽ More Grimme's DFT-D dispersion contribution to interatomic forces constants, required for the computation of the phonon band structures in density-functional perturbation theory, has been derived analytically. The implementation has then been validated with respect to frozen phonons, and applied on materials where weak cohesive forces play a major role i.e. argon, graphite, benzene. We show that these dispersive contributions have to be considered to properly reproduce the experimental vibrational properties of these materials, although the lattice parameter change, coming from the ground-state relaxation with the proper functional, induces the most important change with respect to a treatment without dispersion corrections. In the current implementation, the contribution of these dispersion corrections to the dynamical matrices (with a number of elements that is proportional to the square of the number of atoms) has only a cubic scaling with the number of atoms. In practice the overload with respect to density-functional calculations is small, making this methodology promising to study vibrational properties of large dispersive systems. △ Less

Submitted 26 January, 2018; originally announced January 2018.

Comments: This paper was published in Physical Review B the 20 April 2016 (22 pages, 7 figures)

Journal ref: Phys. Rev. B 93, 144304 (2016)

Large phosphorene in-plane contraction induced by interlayer interactions in graphene-phosphorene heterostructures

Authors: Benoit Van Troeye, Aurélien Lherbier, Jean-Christophe Charlier, Xavier Gonze

Abstract: Intralayer deformation in van der Waals (vdW) heterostructures is generally assumed to be negligible due to the weak nature of the interactions between the layers, especially when the interfaces are found incoherent. In the present work, graphene-phosphorene vdW-heterostructures are investigated with the Density Functional Theory (DFT). The challenge of treating nearly incommensurate (very large)… ▽ More Intralayer deformation in van der Waals (vdW) heterostructures is generally assumed to be negligible due to the weak nature of the interactions between the layers, especially when the interfaces are found incoherent. In the present work, graphene-phosphorene vdW-heterostructures are investigated with the Density Functional Theory (DFT). The challenge of treating nearly incommensurate (very large) supercell in DFT is bypassed by considering different energetic quantities in the grand canonical ensemble, alternative to the formation energy, in order to take into account the mismatch elastic contribution of the different layers. In the investigated heterostructures, it is found that phosphorene contracts by ~4% in the armchair direction when compared to its free-standing form. This large contraction leads to important changes in term of electronic properties, with the direct electronic optical transition of phosphorene becoming indirect in specific vdW-heterostructures. More generally, such a contraction indicates strong substrate effects in supported or encapsulated phosphorene -neglected hitherto- and paves the way to substrate-controlled stress- tronic in such 2D crystal. In addition, the stability of these vdW-heterostructures are investigated as a function of the rotation angle between the layers and as a function of the stacking composition. The alignment of the specific crystalline directions of graphene and phosphorene is found energetically favored. In parallel, several several models based on DFT-estimated quantities are presented; they allow notably a better understanding of the global mutual accommodation of 2D materials in their corresponding interfaces, that is predicted to be non-negligible even in the case of incommensurate interfaces. △ Less

Submitted 23 January, 2018; originally announced January 2018.

Comments: 33 pages, 6 figures

Journal ref: Phys. Rev. Materials 2, 074001 (2018)

Incipient Metals: Functional Materials with a Unique Bonding Mechanism

Authors: Matthias Wuttig, Volker L. Deringer, Xavier Gonze, Christophe Bichara, Jean-Yves Raty

Abstract: While solid-state materials are commonly classified as covalent, ionic, or metallic, there are cases that defy these iconic bonding mechanisms. A prominent example is given by phase-change materials (PCMs) for data storage or photonics, which have recently been argued to show 'resonant' bonding; a clear definition of this mechanism, however, has been lacking until the present day. Here we show tha… ▽ More While solid-state materials are commonly classified as covalent, ionic, or metallic, there are cases that defy these iconic bonding mechanisms. A prominent example is given by phase-change materials (PCMs) for data storage or photonics, which have recently been argued to show 'resonant' bonding; a clear definition of this mechanism, however, has been lacking until the present day. Here we show that these solids are clearly different from resonant bonding in the pi-orbital systems of benzene and graphene. Instead, they exhibit a unique mechanism between covalent and metallic bonding, which we call 'metavalent' bonding. The materials are on the verge of electron delocalization, which explains their exceptional property portfolio, and we therefore argue that they represent 'incipient metals'. This yields deeper, fundamental insight into the bonding nature of solid-state materials, and is expected to accelerate the discovery and design of new functional materials including PCMs and thermoelectrics. △ Less

Submitted 10 December, 2017; originally announced December 2017.

The PseudoDojo: Training and grading a 85 element optimized norm-conserving pseudopotential table

Authors: M. J. van Setten, M. Giantomassi, E. Bousquet, M. J. Verstraete, D. R. Hamann, X. Gonze, G. -M. Rignanese

Abstract: 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 appro… ▽ More 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. ... △ Less

Submitted 30 October, 2017; v1 submitted 27 October, 2017; originally announced October 2017.

Comments: abstract truncated, 17 pages, 25 figures, 8 tables

Journal ref: Computer Physics Communications 226, 39-54, 2018

Quasiparticles and phonon satellites in spectral functions of semiconductors and insulators: Cumulants applied to full first principles theory and Fröhlich polaron

Authors: Jean Paul Nery, Philip B. Allen, Gabriel Antonius, Lucia Reining, Anna Miglio, Xavier Gonze

Abstract: The electron-phonon interaction causes thermal and zero-point motion shifts of electron quasiparticle (QP) energies $ε_k(T)$. Other consequences of interactions, visible in angle-resolved photoemission spectroscopy (ARPES) experiments, are broadening of QP peaks and appearance of sidebands, contained in the electron spectral function $A(k,ω)=-{\Im m}G_R(k,ω) /π$, where $G_R$ is the retarded Gree… ▽ More The electron-phonon interaction causes thermal and zero-point motion shifts of electron quasiparticle (QP) energies $ε_k(T)$. Other consequences of interactions, visible in angle-resolved photoemission spectroscopy (ARPES) experiments, are broadening of QP peaks and appearance of sidebands, contained in the electron spectral function $A(k,ω)=-{\Im m}G_R(k,ω) /π$, where $G_R$ is the retarded Green's function. Electronic structure codes (e.g. using density-functional theory) are now available that compute the shifts and start to address broadening and sidebands. Here we consider MgO and LiF, and determine their nonadiabatic Migdal self energy. The spectral function obtained from the Dyson equation makes errors in the weight and energy of the QP peak and the position and weight of the phonon-induced sidebands. Only one phonon satellite appears, with an unphysically large energy difference (larger than the highest phonon energy) with respect to the QP peak. By contrast, the spectral function from a cumulant treatment of the same self energy is physically better, giving a quite accurate QP energy and several satellites approximately spaced by the LO phonon energy. In particular, the positions of the QP peak and first satellite agree closely with those found for the Fröhlich Hamiltonian by Mishchenko $\textit{et al.}$ (2000) using diagrammatic Monte Carlo. We provide a detailed comparison between the first-principles MgO and LiF results and those of the Fröhlich Hamiltonian. Such an analysis applies widely to materials with infra-red active phonons. We also compare the retarded and time-ordered cumulant treatments: they are equivalent for the Fröhlich Hamiltonian, and only slightly differ in first-principles electron-phonon results for wide-band gap materials. △ Less

Submitted 20 October, 2017; originally announced October 2017.

Comments: 21 pages, 19 figures

Journal ref: Phys. Rev. B 97, 115145 (2018)

Convergence and pitfalls of density functional perturbation theory phonons calculations from a high-throughput perspective

Authors: Guido Petretto, Xavier Gonze, Geoffroy Hautier, Gian-Marco Rignanese

Abstract: The diffusion of large databases collecting different kind of material properties from high-throughput density functional theory calculations has opened new paths in the study of materials science thanks to data mining and machine learning techniques. Phonon calculations have already been employed successfully to predict materials properties and interpret experimental data, e.g. phase stability, f… ▽ More The diffusion of large databases collecting different kind of material properties from high-throughput density functional theory calculations has opened new paths in the study of materials science thanks to data mining and machine learning techniques. Phonon calculations have already been employed successfully to predict materials properties and interpret experimental data, e.g. phase stability, ferroelectricity and Raman spectra, so their availability for a large set of materials will further increase the analytical and predictive power at hand. Moving to a larger scale with density functional perturbation calculations, however, requires the presence of a robust framework to handle this challenging task. In light of this, we automatized the phonon calculation and applied the result to the analysis of the convergence trends for several materials. This allowed to identify and tackle some common problems emerging in this kind of simulations and to lay out the basis to obtain reliable phonon band structures from high-throughput calculations, as well as optimizing the approach to standard phonon simulations. △ Less

Submitted 3 January, 2018; v1 submitted 16 October, 2017; originally announced October 2017.

Journal ref: Computational Materials Science Volume 144, March 2018, Pages 331-337

Automation methodologies and large-scale validation for $GW$, towards high-throughput $GW$ calculations

Authors: M. J. van Setten, M. Giantomassi, X. Gonze, G. -M. Rignanese, G. Hautier

Abstract: 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-functiona… ▽ More 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. △ Less

Submitted 25 September, 2017; v1 submitted 21 September, 2017; originally announced September 2017.

Comments: 12 pages, 11 figures

Journal ref: Phys. Rev. B 96, 155207 (2017)

First-principles Study of the Luminescence of Eu2+-doped Phosphors

Authors: Yongchao Jia, Anna Miglio, Samuel Poncé, Masayoshi Mikami, Xavier Gonze

Abstract: The luminescence of fifteen representative Eu$^{2+}$-doped phosphors used for white-LED and scintillation applications is studied through a Constrained Density Functional Theory. Transition energies and Stokes shift are deduced from differences of total energies between the ground and excited states of the systems, in the absorption and emission geometries. The general applicability of such method… ▽ More The luminescence of fifteen representative Eu$^{2+}$-doped phosphors used for white-LED and scintillation applications is studied through a Constrained Density Functional Theory. Transition energies and Stokes shift are deduced from differences of total energies between the ground and excited states of the systems, in the absorption and emission geometries. The general applicability of such methodology is first assessed: for this representative set, the calculated absolute error with respect to experiment on absorption and emission energies is within 0.3 eV. This set of compounds covers a wide range of transition energies that extents from 1.7 to 3.5 eV. The information gained from the relaxed geometries and total energies is further used to evaluate the thermal barrier for the $4f-5d$ crossover, the full width at half-maximum of the emission spectrum and the temperature shift of the emission peak, using a one-dimensional configuration-coordinate model. The former results indicate that the $4f-5d$ crossover cannot be the dominant mechanism for the thermal quenching behavior of Eu$^{2+}$-doped phosphors and the latter results are compared to available experimental data and yield a 30$\%$ mean absolute relative error. Finally, a semi-empirical model used previously for Ce$^{3+}$-doped hosts is adapted to Eu$^{2+}$-doped hosts and gives the absorption and emission energies within 0.9 eV of experiment, underperforming compared to the first-principles calculation. △ Less

Submitted 1 September, 2017; originally announced September 2017.

Comments: 17 pages, 13 figures, (Phys. Rev. B 2017 Accept)

Journal ref: Phys. Rev. B 96, 125132 (2017)

Assessment of First-Principles and Semiempirical Methodologies for Absorption and Emission Energies of Ce$^{3+}$-Doped Luminescent Materials

Authors: Yongchao Jia, Samuel Poncé, Anna Miglio, Masayoshi Mikami, Xavier Gonze

Abstract: In search of a reliable methodology for the prediction of light absorption and emission of Ce$^{3+}$-doped luminescent materials, 13 representative materials are studied with first-principles and semiempirical approaches. In the first-principles approach, that combines constrained density-functional theory and $Δ$SCF, the atomic positions are obtained for both ground and excited states of the Ce… ▽ More In search of a reliable methodology for the prediction of light absorption and emission of Ce$^{3+}$-doped luminescent materials, 13 representative materials are studied with first-principles and semiempirical approaches. In the first-principles approach, that combines constrained density-functional theory and $Δ$SCF, the atomic positions are obtained for both ground and excited states of the Ce$^{3+}$ ion. The structural information is fed into Dorenbos' semiempirical model. Absorption and emission energies are calculated with both methods and compared with experiment. The first-principles approach matches experiment within 0.3 eV, with two exceptions at 0.5 eV. In contrast, the semiempirical approach does not perform as well (usually more than 0.5 eV error). The general applicability of the present first-principles scheme, with an encouraging predictive power, opens a novel avenue for crystal site engineering and high-throughput search for new phosphors and scintillators. △ Less

Submitted 21 August, 2017; originally announced August 2017.

Comments: 12 pages, 3 figures

Journal ref: Advanced Optical Materials 5, 1600997 (2017)

First-principles study of Ce$^{3+}$ doped lanthanum silicate nitride phosphors: Neutral excitation, Stokes shift, and luminescent center identification

Authors: Yongchao Jia, Anna Miglio, Samuel Poncé, Xavier Gonze, Masayoshi Mikami

Abstract: We study from first principles two lanthanum silicate nitride compounds, LaSi$_{3}$N$_{5}$ and La$_{3}$Si$_{6}$N$_{11}$, pristine as well as doped with Ce$^{3+}$ ion, in view of explaining their different emission color, and characterising the luminescent center. The electronic structures of the two undoped hosts are similar, and do not give a hint to quantitatively describe such difference. The… ▽ More We study from first principles two lanthanum silicate nitride compounds, LaSi$_{3}$N$_{5}$ and La$_{3}$Si$_{6}$N$_{11}$, pristine as well as doped with Ce$^{3+}$ ion, in view of explaining their different emission color, and characterising the luminescent center. The electronic structures of the two undoped hosts are similar, and do not give a hint to quantitatively describe such difference. The $4f\rightarrow 5d$ neutral excitation of the Ce$^{3+}$ ions is simulated through a constrained density-functional theory method coupled with a $Δ$SCF analysis of total energies, yielding absorption energies. Afterwards, atomic positions in the excited state are relaxed, yielding the emission energies and Stokes shifts. Based on these results, the luminescent centers in LaSi$_{3}$N$_{5}$:Ce and La$_{3}$Si$_{6}$N$_{11}$:Ce are identified. The agreement with the experimental data for the computed quantities is quite reasonable and explains the different color of the emitted light. Also, the Stokes shifts are obtained within 20\% difference relative to experimental data. △ Less

Submitted 21 August, 2017; originally announced August 2017.

Comments: 12 pages, 10 figures

Journal ref: Physical Review B 93, 155111 (2016)

Precise effective masses from density functional perturbation theory

Authors: Jonathan Laflamme Janssen, Yannick Gillet, Samuel Poncé, Alexandre Martin, Marc Torrent, Xavier Gonze

Abstract: The knowledge of effective masses is a key ingredient to analyze numerous properties of semiconductors, like carrier mobilities, (magneto-)transport properties, or band extrema characteristics yielding carrier densities and density of states. Currently, these masses are usually calculated using finite-difference estimation of density functional theory (DFT) electronic band curvatures. However, fin… ▽ More The knowledge of effective masses is a key ingredient to analyze numerous properties of semiconductors, like carrier mobilities, (magneto-)transport properties, or band extrema characteristics yielding carrier densities and density of states. Currently, these masses are usually calculated using finite-difference estimation of density functional theory (DFT) electronic band curvatures. However, finite differences require an additional convergence study and are prone to numerical noise. Moreover, the concept of effective mass breaks down at degenerate band extrema. We assess the former limitation by developing a method that allows to obtain the Hessian of DFT bands directly, using density functional perturbation theory (DFPT). Then, we solve the latter issue by adapting the concept of `transport equivalent effective mass' to the $\vec{k} \cdot \hat{\vec{p}}$ framework. The numerical noise inherent to finite-difference methods is thus eliminated, along with the associated convergence study. The resulting method is therefore more general, more robust and simpler to use, which makes it especially appropriate for high-throughput computing. After validating the developed techniques, we apply them to the study of silicon, graphane, and arsenic. The formalism is implemented into the ABINIT software and supports the norm-conserving pseudopotential approach, the projector augmented-wave method, and the inclusion of spin-orbit coupling. The derived expressions also apply to the ultrasoft pseudopotential method. △ Less

Submitted 19 August, 2017; originally announced August 2017.

Comments: 22 pages, 6 figures

Journal ref: Phys. Rev. B 93, 205147 (2016)

Ab Initio Approach to Second-order Resonant Raman Scattering Including Exciton-Phonon Interaction

Authors: Yannick Gillet, Stefan Kontur, Matteo Giantomassi, Claudia Draxl, Xavier Gonze

Abstract: Raman spectra obtained by the inelastic scattering of light by crystalline solids contain contributions from first-order vibrational processes (e.g. the emission or absorption of one phonon, a quantum of vibration) as well as higher-order processes with at least two phonons being involved. At second order, coupling with the entire phonon spectrum induces a response that may strongly depend on the… ▽ More Raman spectra obtained by the inelastic scattering of light by crystalline solids contain contributions from first-order vibrational processes (e.g. the emission or absorption of one phonon, a quantum of vibration) as well as higher-order processes with at least two phonons being involved. At second order, coupling with the entire phonon spectrum induces a response that may strongly depend on the excitation energy, and reflects complex processes more difficult to interpret. In particular, excitons (i.e. bound electron-hole pairs) may enhance the absorption and emission of light, and couple strongly with phonons in resonance conditions. We design and implement a first-principles methodology to compute second-order Raman scattering, incorporating dielectric responses and phonon eigenstates obtained from density-functional theory and many-body theory. We demonstrate our approach for the case of silicon, relating frequency-dependent relative Raman intensities, that are in excellent agreement with experiment, to different vibrations and regions of the Brillouin zone. We show that exciton-phonon coupling, computed from first principles, indeed strongly affect the spectrum in resonance conditions. The ability to analyze second-order Raman spectra thus provides direct insight into this interaction. △ Less

Submitted 19 August, 2017; originally announced August 2017.

Comments: 10 pages, 8 figures

Journal ref: Scientific Reports 7, 7344 (2017)

Origin of the Counterintuitive Dynamic Charge in the Transition-Metal Dichalcogenides

Authors: Nicholas A. Pike, Benoit Van Troeye, Antoine Dewandre, Xavier Gonze, Matthieu J. Verstraete

Abstract: We investigate the chemical bonding characteristics of the transition metal dichalcogenides based on their static and dynamical atomic charges within Density Functional Theory. The dynamical charges of the trigonal transition metal dichalcogenides are anomalously large, while in their hexagonal counterparts, their sign is even counterintuitive i.e. the transition metal takes the negative charge. T… ▽ More We investigate the chemical bonding characteristics of the transition metal dichalcogenides based on their static and dynamical atomic charges within Density Functional Theory. The dynamical charges of the trigonal transition metal dichalcogenides are anomalously large, while in their hexagonal counterparts, their sign is even counterintuitive i.e. the transition metal takes the negative charge. This phenomenon cannot be understood simply in terms of a change in the static atomic charge as it results from a local change of polarization. We present our theoretical understanding of these phenomena based on the perturbative response of the system to a static electric field and by investigating the hybridization of the molecular orbitals near the Fermi level. Furthermore, we establish a link between the sign of the Born effective charge and the $π$-backbonding in organic chemistry and propose an experimental procedure to verify the calculated sign of the dynamical charge in the transition metal dichalcogenides. △ Less

Submitted 12 February, 2017; v1 submitted 6 January, 2017; originally announced January 2017.

Comments: 3 figures, 6 pages

Journal ref: Phys. Rev. B 95, 201106 (2017)

First-principles investigation of the structural, dynamical and dielectric properties of kesterite, stannite and PMCA phases of Cu2ZnSnS4

Authors: S. Poyyapakkam Ramkumar, Y. Gillet, A. Miglio, M. J. van Setten, X. Gonze, G. -M. Rignanese

Abstract: 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 perturbatio… ▽ More 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. △ Less

Submitted 18 October, 2016; originally announced October 2016.

Journal ref: Phys. Rev. B 94, 224302 (2016)

Efficient On-the-Fly Interpolation Technique for Bethe-Salpeter Calculations of Optical Spectra

Authors: Yannick Gillet, Matteo Giantomassi, Xavier Gonze

Abstract: The Bethe-Salpeter formalism represents the most accurate method available nowadays for computing neutral excitation energies and optical spectra of crystalline systems from first principles. Bethe-Salpeter calculations yield very good agreement with experiment but are notoriously difficult to converge with respect to the sampling of the electronic wavevectors. Well-converged spectra therefore req… ▽ More The Bethe-Salpeter formalism represents the most accurate method available nowadays for computing neutral excitation energies and optical spectra of crystalline systems from first principles. Bethe-Salpeter calculations yield very good agreement with experiment but are notoriously difficult to converge with respect to the sampling of the electronic wavevectors. Well-converged spectra therefore require significant computational and memory resources, even by today's standards. These bottlenecks hinder the investigation of systems of great technological interest. They are also barriers to the study of derived quantities like piezoreflectance, thermoreflectance or resonant Raman intensities. We present a new methodology that decreases the workload needed to reach a given accuracy. It is based on a double-grid on-the-fly interpolation within the Brillouin zone, combined with the Lanczos algorithm. It achieves significant speed-up and reduction of memory requirements. The technique is benchmarked in terms of accuracy on silicon, gallium arsenide and lithium fluoride. The scaling of the performance of the method as a function of the Brillouin Zone point density is much better than a conventional implementation. We also compare our method with other similar techniques proposed in the literature. △ Less

Submitted 4 February, 2016; originally announced February 2016.

Comments: 25 pages, 7 figures, Submitted to Computer Physics Communications

Temperature evolution of the band-gap in BiFeO3 traced by resonant Raman scattering

Authors: Mads Christof Weber, Mael Guennou, Constance Toulouse, Maximilien Cazayous, Yannick Gillet, Xavier Gonze, Jens Kreisel

Abstract: Knowledge of the electronic band structure of multiferroic oxides, crucial for the understanding and tuning of photo-induced effects, remains very limited even in the model and thoroughly studied BiFeO3. Here, we investigate the electronic band structure of BiFeO3 using Raman scattering with twelve different excitation wavelengths ranging from the blue to the near infrared. We show that resonant R… ▽ More Knowledge of the electronic band structure of multiferroic oxides, crucial for the understanding and tuning of photo-induced effects, remains very limited even in the model and thoroughly studied BiFeO3. Here, we investigate the electronic band structure of BiFeO3 using Raman scattering with twelve different excitation wavelengths ranging from the blue to the near infrared. We show that resonant Raman signatures can be assigned to direct and indirect electronic transitions, as well as in-gap electronic levels, most likely associated to oxygen vacancies. Their temperature evolution establishes that the remarkable and intriguing variation of the optical band-gap can be related to the shrinking of an indirect electronic band-gap, while the energies for direct electronic transitions remains nearly temperature independent. △ Less

Submitted 8 March, 2016; v1 submitted 13 October, 2015; originally announced October 2015.

Comments: 5 pages, 4 figures

Journal ref: Phys. Rev. B 93, 125204 (2016)

Understanding thermal quenching of photoluminescence from first principles

Authors: Samuel Ponce, Yongchao Jia, Matteo Giantomassi, Masayoshi Mikami, Xavier Gonze

Abstract: Understanding the physical mechanisms behind thermal effects in phosphors is crucial for white light-emitting diodes (WLEDs) applications, as thermal quenching of their photoluminescence might render them useless. The two chemically close Eu-doped \Hosta and \Hostb crystals are typical phosphors studied for WLEDs. The first one sustains efficient light emission at 100$^{\circ}$C while the second o… ▽ More Understanding the physical mechanisms behind thermal effects in phosphors is crucial for white light-emitting diodes (WLEDs) applications, as thermal quenching of their photoluminescence might render them useless. The two chemically close Eu-doped \Hosta and \Hostb crystals are typical phosphors studied for WLEDs. The first one sustains efficient light emission at 100$^{\circ}$C while the second one emits very little light at that temperature. Herein, we analyze from first principles their electronic structure and atomic geometry, before and after absorption/emission of light. Our results, in which the Eu-5d levels are obtained inside the band gap thanks to the removal of an electron from the 4f$^7$ shell, attributes the above-mentioned experimental difference to an auto-ionization model of the thermal quenching, based on the energy difference between Eu$_{\text{5d}}$ and the conduction band minimum. For both Eu-doped phosphors, we identify the luminescent center, and we show that the atomic relaxation in their excited state is of crucial importance for a realistic description of the emission characteristics. △ Less

Submitted 25 August, 2015; originally announced August 2015.

Comments: 5 pages and 3 figures

Dynamical and anharmonic effects on the electron-phonon coupling and the zero-point renormalization of the electronic structure

Authors: Gabriel Antonius, Samuel Poncé, Étienne Lantagne-Hurtubise, Gabriel Auclair, Xavier Gonze, Michel Côté

Abstract: The renormalization of the band structure at zero temperature due to electron-phonon coupling is investigated in diamond, BN, LiF and MgO crystals. We implement a dynamical scheme to compute the frequency-dependent self-energy and the resulting quasiparticle electronic structure. Our calculations reveal the presence of a satellite band below the Fermi level of LiF and MgO. We show that the renorma… ▽ More The renormalization of the band structure at zero temperature due to electron-phonon coupling is investigated in diamond, BN, LiF and MgO crystals. We implement a dynamical scheme to compute the frequency-dependent self-energy and the resulting quasiparticle electronic structure. Our calculations reveal the presence of a satellite band below the Fermi level of LiF and MgO. We show that the renormalization factor (Z), which is neglected in the adiabatic approximation, can reduce the zero-point renormalization (ZPR) by as much as 40%. Anharmonic effects in the renormalized eigenvalues at finite atomic displacements are explored with the frozen-phonon method. We use a non-perturbative expression for the ZPR, going beyond the Allen-Heine-Cardona theory. Our results indicate that high-order electron-phonon coupling terms contribute significantly to the zero-point renormalization for certain materials. △ Less

Submitted 21 August, 2015; v1 submitted 28 May, 2015; originally announced May 2015.

Journal ref: Physical Review B 92, 085137 (2015)

Temperature dependence of the electronic structure of semiconductors and insulators

Authors: Samuel Poncé, Yannick Gillet, Jonathan Laflamme Janssen, Andrea Marini, Matthieu Verstraete, Xavier Gonze

Abstract: The renormalization of electronic eigenenergies due to electron-phonon coupling is sizable in many materials with light atoms. This effect, often neglected in ab-initio calculations, can be computed using the perturbation-based Allen-Heine-Cardona theory in the adiabatic or non-adiabatic harmonic approximation. After a short description of the numerous recent progresses in this field, and a brief… ▽ More The renormalization of electronic eigenenergies due to electron-phonon coupling is sizable in many materials with light atoms. This effect, often neglected in ab-initio calculations, can be computed using the perturbation-based Allen-Heine-Cardona theory in the adiabatic or non-adiabatic harmonic approximation. After a short description of the numerous recent progresses in this field, and a brief overview of the theory, we focus on the issue of phonon wavevector sampling convergence, until now poorly understood. Indeed, the renormalization is obtained numerically through a q-point sampling inside the BZ. For q-points close to G, we show that a divergence due to non-zero Born effective charge appears in the electron-phonon matrix elements, leading to a divergence of the integral over the BZ for band extrema. Although it should vanish for non-polar materials, unphysical residual Born effective charges are usually present in ab-initio calculations. Here, we propose a solution that improves the coupled q-point convergence dramatically. For polar materials, the problem is more severe: the divergence of the integral does not disappear in the adiabatic harmonic approximation, but only in the non-adiabatic harmonic approximation. In all cases, we study in detail the convergence behavior of the renormalization as the q-point sampling goes to infinity and the imaginary broadening parameter goes to zero. This allows extrapolation, thus enabling a systematic way to converge the renormalization for both polar and non-polar materials. Finally, the adiabatic and non-adiabatic theory, with corrections for the divergence problem, are applied to the study of five semiconductors and insulators: a-AlN, b-AlN, BN, diamond and silicon. For these five materials, we present the zero-point renormalization, temperature dependence, phonon-induced lifetime broadening and the renormalized electronic bandstructure. △ Less

Submitted 22 April, 2015; originally announced April 2015.

Comments: 27 pages and 26 figures

A many-body perturbation theory approach to the electron-phonon interaction with density-functional theory as a starting point

Authors: Andrea Marini, Samuel Poncé, Xavier Gonze

Abstract: The electron-phonon interaction plays a crucial role in many fields of physics and chemistry. Nevertheless, its actual calculation by means of modern many-body perturbation theory is weakened by the use of model Hamiltonians that are based on parameters difficult to extract from the experiments. Such shortcoming can be bypassed by using density-functional theory to evaluate the electron-phonon sca… ▽ More The electron-phonon interaction plays a crucial role in many fields of physics and chemistry. Nevertheless, its actual calculation by means of modern many-body perturbation theory is weakened by the use of model Hamiltonians that are based on parameters difficult to extract from the experiments. Such shortcoming can be bypassed by using density-functional theory to evaluate the electron-phonon scattering amplitudes, phonon frequencies and electronic bare energies. In this work, we discuss how a consistent many-body diagrammatic expansion can be constructed on top of density-functional theory. In that context, the role played by screening and self-consistency when all the components of the electron-nucleus and nucleus-nucleus interactions are taken into account is paramount. A way to avoid over-screening is notably presented. Finally, we derive cancellations rules as well as internal consistency constraints in order to draw a clear, sound and practical scheme to merge many-body perturbation and density-functional theory. △ Less

Submitted 14 May, 2015; v1 submitted 2 March, 2015; originally announced March 2015.

Comments: 25 pages, 13 figures

High-Mobility Bismuth-based Transparent P-Type Oxide from High-throughput Material Screening

Authors: Amit Bathia, Geoffroy Hautier, Tan Nilgianskul, Anna Miglio, Gian-Marco Rignanese, Xavier Gonze, Jin Suntivich

Abstract: Transparent oxides are essential building blocks to many technologies, ranging from components in transparent electronics, transparent conductors, to absorbers and protection layers in photovoltaics and photoelectrochemical devices. However, thus far, it has been difficult to develop p-type oxides with wide band gap and high hole mobility; current state-of-art transparent p-type oxides have hole m… ▽ More Transparent oxides are essential building blocks to many technologies, ranging from components in transparent electronics, transparent conductors, to absorbers and protection layers in photovoltaics and photoelectrochemical devices. However, thus far, it has been difficult to develop p-type oxides with wide band gap and high hole mobility; current state-of-art transparent p-type oxides have hole mobility in the range of < 10 cm$^2$/Vs, much lower than their n-type counterparts. Using high-throughput computational screening to guide the discovery of novel oxides with wide band gap and high hole mobility, we report the computational identification and the experimental verification of a bismuth-based double-perovskite oxide that meets these requirements. Our identified candidate, Ba$_2$BiTaO$_6$, has an optical band gap larger than 4 eV and a Hall hole mobility above 30 cm$^2$/Vs. We rationalize this finding with molecular orbital intuitions; Bi$^{3+}$ with filled s-orbitals strongly overlap with the oxygen p, increasing the extent of the metal-oxygen covalency and effectively reducing the valence effective mass, while Ta$^{5+}$ forms a conduction band with low electronegativity, leading to a high band gap beyond the visible range. Our concerted theory-experiment effort points to the growing utility of a data-driven materials discovery and the combination of both informatics and chemical intuitions as a way to discover future technological materials. △ Less

Submitted 14 December, 2014; originally announced December 2014.

Journal ref: Chem. Mater., 2016, 28, pp 30-34

The temperature dependence of electronic eigenenergies in the adiabatic harmonic approximation

Authors: Samuel Poncé, Gabriel Antonius, Yannick Gillet, Paul Boulanger, Jonathan Laflamme Janssen, Andrea Marini, Michel Côté, Xavier Gonze

Abstract: The renormalization of electronic eigenenergies due to electron-phonon interactions (temperature dependence and zero-point motion effect) is important in many materials. We address it in the adiabatic harmonic approximation, based on first principles (e.g. Density-Functional Theory), from different points of view: directly from atomic position fluctuations or, alternatively, from Janak's theorem g… ▽ More The renormalization of electronic eigenenergies due to electron-phonon interactions (temperature dependence and zero-point motion effect) is important in many materials. We address it in the adiabatic harmonic approximation, based on first principles (e.g. Density-Functional Theory), from different points of view: directly from atomic position fluctuations or, alternatively, from Janak's theorem generalized to the case where the Helmholtz free energy, including the vibrational entropy, is used. We prove their equivalence, based on the usual form of Janak's theorem and on the dynamical equation. We then also place the Allen-Heine-Cardona (AHC) theory of the renormalization in a first-principle context. The AHC theory relies on the rigid-ion approximation, and naturally leads to a self-energy (Fan) contribution and a Debye-Waller contribution. Such a splitting can also be done for the complete harmonic adiabatic expression, in which the rigid-ion approximation is not required. A numerical study within the Density-Functional Perturbation theory framework allows us to compare the AHC theory with frozen-phonon calculations, with or without the rigid-ion terms. For the two different numerical approaches without rigid-ion terms, the agreement is better than 7 $μ$eV in the case of diamond, which represent an agreement to 5 significant digits. The magnitude of the non rigid-ion terms in this case is also presented, distinguishing specific phonon modes contributions to different electronic eigenenergies. △ Less

Submitted 12 August, 2014; originally announced August 2014.

Structural, electronic, vibrational and dielectric properties of LaBGeO$_5$ from first principles

Authors: R. Shaltaf, H. K. Juwhari, B. Hamad, J. Khalifeh, G. -M. Rignanese, X. Gonze

Abstract: 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 fiv… ▽ More 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. △ Less

Submitted 29 November, 2013; originally announced November 2013.