-
Quantum effects and temperature dependence of the free energy barrier and shear frequency in bilayer graphene
Authors:
Jean Paul Nery,
Lorenzo Monacelli,
Francesco Mauri
Abstract:
Little is known about the quantum effects of friction, and first-principles studies are virtually non-existent. In multilayer graphene, the order of the stacking of the layers plays a crucial role in the electronic properties and the manifestation of superconductivity, and friction affects how layers slide to change the stacking. Here, we study the effects of friction in bilayer graphene by focusi…
▽ More
Little is known about the quantum effects of friction, and first-principles studies are virtually non-existent. In multilayer graphene, the order of the stacking of the layers plays a crucial role in the electronic properties and the manifestation of superconductivity, and friction affects how layers slide to change the stacking. Here, we study the effects of friction in bilayer graphene by focusing on quantum fluctuations in the free energy barrier that layers have to overcome to slide with respect to each other. We evaluate the free energy at the barrier in bilayer graphene through the self-consistent harmonic approximation, accounting for quantum and thermal ionic fluctuations. We obtain a very large reduction of the barrier of more than 30%. Temperature effects are weaker, although high temperatures induce a significant reduction of the barrier as well. Our approach is general and paves the way for systematically accounting for nuclear quantum and thermal effects in free energy barriers of other macroscopic systems.
△ Less
Submitted 5 July, 2024;
originally announced July 2024.
-
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.
-
Non-perturbative Green's function method to determine the electronic spectral function due to electron-phonon interactions: Application to a graphene model from weak to strong coupling
Authors:
Jean Paul Nery,
Francesco Mauri
Abstract:
In solid state physics, the electron-phonon interaction (EPI) is central to many phenomena. The theory of the renormalization of electronic properties due to EPIs became well established with the theory of Allen-Heine-Cardona, usually applied to second order in perturbation theory (P2). However, this is only valid in the weak coupling regime, while strong EPIs have been reported in many materials.…
▽ More
In solid state physics, the electron-phonon interaction (EPI) is central to many phenomena. The theory of the renormalization of electronic properties due to EPIs became well established with the theory of Allen-Heine-Cardona, usually applied to second order in perturbation theory (P2). However, this is only valid in the weak coupling regime, while strong EPIs have been reported in many materials. Although non-perturbative (NP) methods have started to arise in the last years, they are usually not well justified, and it is not clear to what degree they reproduce the exact theory. To address this issue, we present a stochastic approach for the evaluation of the non-perturbative interacting Green's function in the adiabatic limit, and show it is equivalent to the Feynman expansion to all orders in the perturbation. Also, by defining a self-energy, we can reduce the effect of broadening needed in numerical calculations, improving convergence in the supercell size. In addition, we clarify whether it is better to average the Green's function or self-energy. Then we apply the method to a graphene tight-binding model, and obtain several interesting results: (i) The Debye-Waller term, which is normally neglected, does affect the change of the Fermi velocity. (ii) The P2 and NP self-energies differ even at room temperature for some k-points, raising the question of how well P2 works in other materials. (iii) Close to the Dirac point, positive and negative energy peaks merge. (iv) In the strong coupling regime, a peak appears at energy E=0, which is consistent with previous works on disorder and localization in graphene. (v) The spectral function becomes more asymmetric at stronger coupling and higher temperatures. Finally, in the Appendix we show that the method has better convergent properties when the coupling is strong relative to when it is weak, and discuss other technical aspects.
△ Less
Submitted 24 May, 2022; v1 submitted 21 March, 2022;
originally announced March 2022.
-
Ab-initio energetics of graphite and multilayer graphene: stability of Bernal versus rhombohedral stacking
Authors:
Jean Paul Nery,
Matteo Calandra,
Francesco Mauri
Abstract:
There has been a lot of excitement around the observation of superconductivity in twisted bilayer graphene, associated to flat bands close to the Fermi level. Such correlated electronic states also occur in multilayer rhombohedral stacked graphene (RG), which has been receiving increasing attention in the last years. In both natural and artificial samples however, multilayer stacked Bernal graphen…
▽ More
There has been a lot of excitement around the observation of superconductivity in twisted bilayer graphene, associated to flat bands close to the Fermi level. Such correlated electronic states also occur in multilayer rhombohedral stacked graphene (RG), which has been receiving increasing attention in the last years. In both natural and artificial samples however, multilayer stacked Bernal graphene (BG) occurs more frequently, making it desirable to determine what is their relative stability and under which conditions RG might be favored. Here, we study the energetics of BG and RG in bulk and also multilayer stacked graphene using first-principles calculations. It is shown that the electronic temperature, not accounted for in previous studies, plays a crucial role in determining which phase is preferred. We also show that the low energy states at room temperature consist of BG, RG and mixed BG-RG systems with a particular type of interface. Energies of all stacking sequences (SSs) are calculated for N = 12 layers, and an Ising model is used to fit them, which can be used for larger N as well. In this way, the ordering of low energy SSs can be determined and analyzed in terms of a few parameters. Our work clarifies inconsistent results in the literature, and sets the basis to studying the effect of external factors on the stability of multilayer graphene systems in first principles calculations.
△ Less
Submitted 24 February, 2021; v1 submitted 18 November, 2020;
originally announced November 2020.
-
Long-range rhombohedral-stacked graphene through shear
Authors:
Jean Paul Nery,
Matteo Calandra,
Francesco Mauri
Abstract:
The discovery of superconductivity and correlated electronic states in the flat bands of twisted bilayer graphene has raised a lot of excitement. Flat bands also occur in multilayer graphene flakes that present rhombohedral (ABC) stacking order on many consecutive layers. Although Bernal-stacked (AB) graphene is more stable, long-range ABC-ordered flakes involving up to 50 layers have been surpris…
▽ More
The discovery of superconductivity and correlated electronic states in the flat bands of twisted bilayer graphene has raised a lot of excitement. Flat bands also occur in multilayer graphene flakes that present rhombohedral (ABC) stacking order on many consecutive layers. Although Bernal-stacked (AB) graphene is more stable, long-range ABC-ordered flakes involving up to 50 layers have been surprisingly observed in natural samples. Here we present a microscopic atomistic model, based on first-principles density functional theory calculations, that demonstrates how shear stress can produce long-range ABC order. A stress-angle phase diagram shows under which conditions ABC-stacked graphene can be obtained, providing an experimental guide for its synthesis.
△ Less
Submitted 18 February, 2020;
originally announced February 2020.
-
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.
-
Low temperature-semiconductor band gap thermal shifts: T^4 shifts from ordinary acoustic and T^2 from piezoacoustic coupling
Authors:
Philip B. Allen,
Jean Paul Nery
Abstract:
At low temperature T, the experimental gap of silicon decreases as E_g(T)=E_g(0)-AT^4. The main reason is electron-phonon renormalization. The physics behind the T^4-power law is more complex than has been realized. Renormalization by intraband scattering requires a careful non-adiabatic treatment in order to correctly include acoustic phonons and avoid divergences from piezoacoustic phonon intera…
▽ More
At low temperature T, the experimental gap of silicon decreases as E_g(T)=E_g(0)-AT^4. The main reason is electron-phonon renormalization. The physics behind the T^4-power law is more complex than has been realized. Renormalization by intraband scattering requires a careful non-adiabatic treatment in order to correctly include acoustic phonons and avoid divergences from piezoacoustic phonon interactions. The result is an unexpected low T term E_g(0)+A' T^p with positive coefficient A', and power p=4 for non-piezoelectric materials, and power p=2 for piezoelectric materials. The acoustic phonons in piezoelectric semiconductors generate a piezoelectric field, modifying the electron-phonon coupling. However, at higher T, when thermal acoustic phonons of energy hbar v_s q acquire energies comparable to the electronic intermediate state (higher than the band-edge state by hbar^2 q^2 /2m*), the low q and higher q intraband contributions to T^p rapidly cancel, giving little thermal effect. But there is an additional T-dependence from interband effects of acoustic phonons. This turns out to have power law T^4 for both non-piezoelectric and piezoelectric semiconductors. This term can have either sign, but usually reduces the size of gaps as T increases. It arises after cancellation of the T^2 terms that appear separately in Debye-Waller and Fan parts of the acoustic phonon interband renormalization. The cancellation occurs because of the acoustic sum rule.
△ Less
Submitted 26 January, 2017; v1 submitted 6 May, 2016;
originally announced May 2016.
-
Influence of Fröhlich polaron coupling on renormalized electron bands in polar semiconductors. Results for zincblende GaN
Authors:
Jean Paul Nery,
Philip B. Allen
Abstract:
\ni We develop a simple method to study the zero-point and thermally renormalized electron energy $\varepsilon_{\mathbf{k}n}(T)$ for $\mathbf{k}n$ the conduction band minimum or valence maximum in polar semiconductors. We use the adiabatic approximation, including an imaginary broadening parameter $iδ$ to supress noise in the density-functional integrations. Fröhlich polaron methods provide analyt…
▽ More
\ni We develop a simple method to study the zero-point and thermally renormalized electron energy $\varepsilon_{\mathbf{k}n}(T)$ for $\mathbf{k}n$ the conduction band minimum or valence maximum in polar semiconductors. We use the adiabatic approximation, including an imaginary broadening parameter $iδ$ to supress noise in the density-functional integrations. Fröhlich polaron methods provide analytic expressions for the contribution of the problematic optical phonon mode. We use this to correct the renormalization obtained from the adiabatic approximation. Test calculations are done for zincblende GaN for an 18x18x18 integration grid. The Fröhlich correction is of order -0.02 eV for the zero-point energy shift of the conduction band minimum, and +0.03 eV for the valence band maximum; the correction to renormalization of the 3.28 eV gap is -0.05 eV, a significant fraction of the total zero point renormalization of -0.15 eV.
△ Less
Submitted 21 July, 2016; v1 submitted 14 March, 2016;
originally announced March 2016.
-
Boundary divergences in vacuum self-energies and quantum field theory in curved spacetime
Authors:
Francisco D. Mazzitelli,
Jean Paul Nery,
Alejandro Satz
Abstract:
It is well known that boundary conditions on quantum fields produce divergences in the renormalized energy-momentum tensor near the boundaries. Although irrelevant for the computation of Casimir forces between different bodies, the self-energy couples to gravity, and the divergences may, in principle, generate large gravitational effects. We present an analysis of the problem in the context of qua…
▽ More
It is well known that boundary conditions on quantum fields produce divergences in the renormalized energy-momentum tensor near the boundaries. Although irrelevant for the computation of Casimir forces between different bodies, the self-energy couples to gravity, and the divergences may, in principle, generate large gravitational effects. We present an analysis of the problem in the context of quantum field theory in curved spaces. Our model consists of a quantum scalar field coupled to a classical field that, in a certain limit, imposes Dirichlet boundary conditions on the quantum field. We show that the model is renormalizable and that the divergences in the renormalized energy-momentum tensor disappear for sufficiently smooth interfaces.
△ Less
Submitted 16 October, 2011;
originally announced October 2011.
