Source author record

Xavier Gonze

Xavier Gonze appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2025arXiv

Low-temperature behavior of density-functional theory for metals based on density-functional perturbation theory and Sommerfeld expansion

The temperature dependence of most solid-state properties is dominated by lattice vibrations, but metals display notable purely electronic effects at low temperature, such as the linear specific heat and the linear entropy, that were derived by Sommerfeld for the non-interacting electron gas via the low-temperature expansion of Fermi-Dirac integrals. Here we treat temperature as a perturbation within density-functional perturbation theory (DFPT). For finite temperature, we show how self-consistency screens the bare, temperature-induced density change obtained in the non-interacting picture: the inverse transpose of the electronic dielectric operator, that includes Adler-Wiser and a term related to the shift in Fermi level, links the self-consistent density response to the bare thermal density change. This approach is implemented in DFTK, and demonstrated by the computation of the second-order derivative of the free energy, and the first-order derivative of entropy for aluminum. Then, we examine the $T\!\to\!0$ limit. The finite temperature formalism contains divergences, that we cure using the Sommerfeld expansion to analyze metallic systems at 0 K. The electronic free energy is quadratic in $T$ provided the Fermi level is not at a Van Hove singularity of the density of states. If the latter happens, another temperature behavior might appear, depending on the type of Van Hove singularity, that we analyze. Our formulation applies to systems periodic in one, two, or three dimensions, and provides a basis for studying temperature-dependent electronic instabilities (e.g., charge-density waves) within density-functional theory and DFPT.

preprint2022arXiv

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

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.

preprint2022arXiv

Spectroscopic signatures of nonpolarons : the case of diamond

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.

preprint2022arXiv

Variational Polaron Equations Applied to the Anisotropic Fröhlich Model

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.

preprint2022arXiv

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

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.

preprint2021arXiv

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

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.

preprint2020arXiv

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

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.

preprint2016arXiv

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

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.

preprint2016arXiv

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

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.

preprint2015arXiv

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

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.

preprint2015arXiv

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

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.

preprint2015arXiv

Temperature dependence of the electronic structure of semiconductors and insulators

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.

preprint2015arXiv

Understanding thermal quenching of photoluminescence from first principles

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.

preprint2014arXiv

The temperature dependence of electronic eigenenergies in the adiabatic harmonic approximation

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.

preprint2013arXiv

Band widths and gaps from the Tran-Blaha functional : Comparison with many-body perturbation theory

For a set of ten crystalline materials (oxides and semiconductors), we compute the electronic band structures using the Tran-Blaha [Phys. Rev. Lett. 102, 226401 (2009)] (TB09) functional. The band widths and gaps are compared with those from the local-density approximation (LDA) functional, many-body perturbation theory (MBPT), and experiments. At the density-functional theory (DFT) level, TB09 leads to band gaps in much better agreement with experiments than LDA. However, we observe that it globally underestimates, often strongly, the valence (and conduction) band widths (more than LDA). MBPT corrections are calculated starting from both LDA and TB09 eigenenergies and wavefunctions. They lead to a much better agreement with experimental data for band widths. The band gaps obtained starting from TB09 are close to those from quasi-particle self-consistent GW calculations, at a much reduced cost. Finally, we explore the possibility to tune one of the semi-empirical parameters of the TB09 functional in order to obtain simultaneously better band gaps and widths. We find that these requirements are conflicting.

preprint2013arXiv

First-principles study of excitonic effects in Raman intensities

The ab initio prediction of Raman intensities for bulk solids usually relies on the hypothesis that the frequency of the incident laser light is much smaller than the band gap. However, when the photon frequency is a sizeable fraction of the energy gap, or higher, resonance effects appear. In the case of silicon, when excitonic effects are neglected, the response of the solid to light increases by nearly three orders of magnitude in the range of frequencies between the static limit and the gap. When excitonic effects are taken into account, an additional tenfold increase in the intensity is observed. We include these effects using a finite-difference scheme applied on the dielectric function obtained by solving the Bethe-Salpeter equation. Our results for the Raman susceptibility of silicon show stronger agreement with experimental data compared with previous theoretical studies. For the sampling of the Brillouin zone, a double-grid technique is proposed, resulting in a significant reduction in computational effort.

preprint2011arXiv

Implementation and testing of Lanczos-based algorithms for Random-Phase Approximation eigenproblems

The treatment of the Random-Phase Approximation Hamiltonians, encountered in different frameworks, like Time-Dependent Density Functional Theory or Bethe-Salpeter equation, is complicated by their non-Hermicity. Compared to their Hermitian Hamiltonian counterparts, computational methods for the treatment of non-Hermitian Hamiltonians are often less efficient and less stable, sometimes leading to the breakdown of the method. Recently [Grüning et al. Nano Lett. {\bf 8}, 2820 (2009)], we have identified that such Hamiltonians are usually pseudo-Hermitian. Exploiting this property, we have implemented an algorithm of the Lanczos type for random-Phase Approximation Hamiltonians that benefits from the same stability and computational load as its Hermitian counterpart, and applied it to the study of the optical response of carbon nanotubes. We present here the related theoretical grounds and technical details, and study the performance of the algorithm for the calculation of the optical absorption of a molecule within the Bethe-Salpeter equation framework.

preprint2008arXiv

Exciton-plasmon states in nanoscale materials: breakdown of the Tamm-Dancoff approximation

Within the Tamm-Dancoff approximation ab initio approaches describe excitons as packets of electron-hole pairs propagating only forward in time. However, we show that in nanoscale materials excitons and plasmons hybridize, creating exciton--plasmon states where the electron-hole pairs oscillate back and forth in time. Then, as exemplified by the trans-azobenzene molecule and carbon nanotubes, the Tamm-Dancoff approximation yields errors as large as the accuracy claimed in ab initio calculations. Instead, we propose a general and efficient approach that avoids the Tamm--Dancoff approximation, and correctly describes excitons, plasmons and exciton-plasmon states.

preprint2004arXiv

First-principles study of the electrooptic effect in ferroelectric oxides

We present a method to compute the electrooptic tensor from first principles, explicitly taking into account the electronic, ionic and piezoelectric contributions. It allows us to study the non-linear optic behavior of three ferroelectric ABO_3 compounds : LiNbO_3, BaTiO_3 and PbTiO_3. Our calculations reveal the dominant contribution of the soft mode to the electrooptic coefficients in LiNbO_3 and BaTiO_3 and identify the coupling between the electric field and the polar atomic displacements along the B-O chains as the origin of the large electrooptic response in these compounds.