Showing 1–50 of 145 results for author: Calandra, M

Anomalous amplitude mode dynamics below the expected charge-density-wave transition in 1$T$-VSe$_2$

Authors: Charles J. Sayers, Giovanni Marini, Matteo Calandra, Hamoon Hedayat, Xuanbo Feng, Erik van Heumen, Christoph Gadermaier, Stefano Dal Conte, Giulio Cerullo

Abstract: A charge-density-wave (CDW) is characterized by a dynamical order parameter consisting of a time-dependent amplitude and phase, which manifest as optically-active collective modes of the CDW phase. Studying the behaviour of such collective modes in the time-domain, and their coupling with electronic and lattice order, provides important insight into the underlying mechanisms behind CDW formation.… ▽ More A charge-density-wave (CDW) is characterized by a dynamical order parameter consisting of a time-dependent amplitude and phase, which manifest as optically-active collective modes of the CDW phase. Studying the behaviour of such collective modes in the time-domain, and their coupling with electronic and lattice order, provides important insight into the underlying mechanisms behind CDW formation. In this work, we report on femtosecond broadband transient reflectivity experiments on bulk 1$T$-VSe$_2$ using near-infrared excitation. At low temperature, we observe coherent oscillations associated with the CDW amplitude mode and phonons of the distorted lattice. Across the expected transition temperature at 110 K, we confirm signatures of a rearrangement of the electronic structure evident in the quasiparticle dynamics. However, we find that the amplitude mode instead softens to zero frequency at 80 K, possibly indicating an additional phase transition at this temperature. In addition, we demonstrate photoinduced CDW melting, associated with a collapse of the electronic and lattice order, which occurs at moderate excitation densities, consistent with a dominant electron-phonon CDW mechanism. △ Less

Submitted 29 August, 2024; originally announced August 2024.

Comments: 13 pages, 6 figures, supporting information

Ultrafast creation of a light induced semimetallic state in strongly excited 1T-TiSe$_2$

Authors: Maximilian Huber, Yi Lin, Giovanni Marini, Luca Moreschini, Chris Jozwiak, Aaron Bostwick, Matteo Calandra, Alessandra Lanzara

Abstract: Screening, a ubiquitous phenomenon associated with the shielding of electric fields by surrounding charges, has been widely adopted as a means to modify a material's properties. While so far most studies have relied on static changes of screening through doping or gating, here we demonstrate that screening can also drive the onset of distinct quantum states on the ultrafast timescale. By using tim… ▽ More Screening, a ubiquitous phenomenon associated with the shielding of electric fields by surrounding charges, has been widely adopted as a means to modify a material's properties. While so far most studies have relied on static changes of screening through doping or gating, here we demonstrate that screening can also drive the onset of distinct quantum states on the ultrafast timescale. By using time and angle-resolved photoemission spectroscopy we show that intense optical excitation can drive 1T-TiSe$_2$, a prototypical charge density wave material, almost instantly from a gapped into a semimetallic state. By systematically comparing changes in bandstructure over time and excitation strength with theoretical calculations we find that the appearance of this state is likely caused by a dramatic reduction of the screening length. In summary, this work showcases how optical excitation enables the screening driven design of a non-equilibrium semimetallic phase in TiSe$_2$, possibly providing a general pathway into highly screened phases in other strongly correlated materials. △ Less

Submitted 16 August, 2024; originally announced August 2024.

Journal ref: Science Advances 10 (2024)

NbSe$_{2}$'s charge density wave collapse in the (LaSe)$_{1.14}$(NbSe$_{2}$)$_{2}$ misfit layer compound

Authors: Ludovica Zullo, Grégory Setnikar, Amit Pawbake, Tristan Cren, Christophe Brun, Justine Cordiez, Shunsuke Sasaki, Laurent Cario, Giovanni Marini, Matteo Calandra, Marie-Aude Méasson

Abstract: Misfit layer compounds, heterostructures composed by a regular alternating stacking of rocksalt monochalcogenides bilayers and few-layer transition metal dichalchogenides, are an emergent platform to investigate highly doped transition metal dichalcogenides. Among them, (LaSe)$_{1.14}$(NbSe$_2$)$_2$ displays Ising superconductivity, while the presence of a charge density wave (CDW) in the material… ▽ More Misfit layer compounds, heterostructures composed by a regular alternating stacking of rocksalt monochalcogenides bilayers and few-layer transition metal dichalchogenides, are an emergent platform to investigate highly doped transition metal dichalcogenides. Among them, (LaSe)$_{1.14}$(NbSe$_2$)$_2$ displays Ising superconductivity, while the presence of a charge density wave (CDW) in the material is still under debate. Here, by using polarized Raman spectroscopy and first-principles calculations, we show that NbSe$_2$ undergoes a doping-driven collapse of the CDW ordering within the misfit, and no signature of the CDW is detected down to 8~K. We provide a complete experimental and theoretical description of the lattice dynamics of this misfit compound. We show that the vibrational properties are obtained from those of the two subunits, namely the LaSe unit and the NbSe$_2$ bilayer, in the presence of a suitable field-effect doping, and then highlight the 2D nature of the lattice dynamics of NbSe$_2$ within the (LaSe)$_{1.14}$(NbSe$_2$)$_2$ 3D structure. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: First version of the manuscript (before revisions)

First-order rhombohedral to cubic phase transition in photoexcited GeTe

Authors: Matteo Furci, Giovanni Marini, Matteo Calandra

Abstract: Photoexcited GeTe undergoes a non-thermal phase transition from a rhombohedral to a rocksalt crystalline phase. The microscopic mechanism and the nature of the transition are unclear. By using constrained density functional perturbation theory and by accounting for quantum anharmonicity within the stochastic self-consistent harmonic approximation, we show that the non-thermal phase transition is s… ▽ More Photoexcited GeTe undergoes a non-thermal phase transition from a rhombohedral to a rocksalt crystalline phase. The microscopic mechanism and the nature of the transition are unclear. By using constrained density functional perturbation theory and by accounting for quantum anharmonicity within the stochastic self-consistent harmonic approximation, we show that the non-thermal phase transition is strongly first order and does not involve phonon softening, at odd with the thermal one. The transition is driven by the closure of the single particle gap in the photoexcited rhombohedral phase. Finally, we show that ultrafast X-ray diffraction data are consistent with a coexistence of the two phases, as expected in a first order transition. Our results are relevant for the understanding of phase transitions and bonding in phase change materials. △ Less

Submitted 30 April, 2024; originally announced April 2024.

Comments: 5 pages, 6 figures To be published in Phys. Rev. Lett

Optical absorption and photoluminescence of single layer boron nitride from a first principles cumulant approach

Authors: Giovanni Marini, Matteo Calandra, Pierluigi Cudazzo

Abstract: The photoluminescence spectrum of a single-layer boron nitride remains elusive, marked by enigmatic satellites that hint at a significant but unidentified exciton-phonon coupling.Here, by employing a first principles approach based on the many-body cumulant expansion of the charge response, we calculate the optical absorption and photoluminescence of a single layer boron nitride. We identify the s… ▽ More The photoluminescence spectrum of a single-layer boron nitride remains elusive, marked by enigmatic satellites that hint at a significant but unidentified exciton-phonon coupling.Here, by employing a first principles approach based on the many-body cumulant expansion of the charge response, we calculate the optical absorption and photoluminescence of a single layer boron nitride. We identify the specific exciton-phonon scattering channels and unravel their impact on the optical absorption and photoluminescence spectra thereby providing an interpretation of the experimental features. Finally, we show that, even in a strongly polar material such as h-BN monolayer, the electron-hole interaction responsible for the excitonic effect results in the cancellation of the Frölich interaction at small phonon momenta. This effect is captured only if the invariance of the exciton-phonon matrix elements under unitary transformations in the Bloch functions manifold is preserved in the calculation. △ Less

Submitted 6 February, 2024; originally announced February 2024.

Electrical tuning of the magnetic properties of 2D magnets: the case of ${\rm Cr}_2{\rm Ge}_2{\rm Te}_6$

Authors: Guido Menichetti, Matteo Calandra, Marco Polini

Abstract: Motivated by growing interest in atomically-thin van der Waals magnetic materials, we present an {\it ab initio} theoretical study of the dependence of their magnetic properties on the electron/hole density $ρ$ induced via the electrical field effect. By focusing on the case of monolayer ${\rm Cr}_2{\rm Ge}_2{\rm Te}_6$ (a prototypical 2D Ising ferromagnet) and employing a hybrid functional, we fi… ▽ More Motivated by growing interest in atomically-thin van der Waals magnetic materials, we present an {\it ab initio} theoretical study of the dependence of their magnetic properties on the electron/hole density $ρ$ induced via the electrical field effect. By focusing on the case of monolayer ${\rm Cr}_2{\rm Ge}_2{\rm Te}_6$ (a prototypical 2D Ising ferromagnet) and employing a hybrid functional, we first study the dependence of the gap and effective mass on the carrier concentration $ρ$. We then investigate the robustness of magnetism by studying the dependencies of the exchange couplings and magneto-crystalline anisotropy energy (MAE) on $ρ$. In agreement with experimental results, we find that magnetism displays a bipolar electrically-tunable character, which is, however, much more robust for hole ($ρ>0$) rather than electron ($ρ<0$) doping. Indeed, the MAE vanishes for an electron density $ρ\approx - 7.5 \times 10^{13}~{\rm e} \times {\rm cm}^{-2}$, signalling the failure of a localized description based on a Heisenberg-type anisotropic spin Hamiltonian. This is in agreement with the rapid increase of the coupling between fourth-neighbor atoms with increasing electron density. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: 13 pages, 11 figures

Journal ref: Phys. Rev. B 109, 205139 (2024)

Theory of infrared double-resonance Raman spectrum in graphene: the role of the zone-boundary electron-phonon enhancement

Authors: Lorenzo Graziotto, Francesco Macheda, Thibault Sohier, Matteo Calandra, Francesco Mauri

Abstract: We theoretically investigate the double-resonance Raman spectrum of monolayer graphene down to infrared laser excitation energies. By using first-principles density functional theory calculations, we improve upon previous theoretical predictions based on conical models or tight-binding approximations, and rigorously justify the evaluation of the electron-phonon enhancement found in Ref. [Venanzi,… ▽ More We theoretically investigate the double-resonance Raman spectrum of monolayer graphene down to infrared laser excitation energies. By using first-principles density functional theory calculations, we improve upon previous theoretical predictions based on conical models or tight-binding approximations, and rigorously justify the evaluation of the electron-phonon enhancement found in Ref. [Venanzi, T., Graziotto, L. et al., Phys. Rev. Lett. 130, 256901 (2023)]. We proceed to discuss the effects of such enhancement on the room temperature graphene resistivity, hinting towards a possible reconciliation of theoretical and experimental discrepancies. △ Less

Submitted 19 January, 2024; v1 submitted 13 October, 2023; originally announced October 2023.

Comments: 19 pages, 18 figures

Light-Induced Nonthermal Phase Transition to the Topological Crystalline Insulator State in SnSe

Authors: Stefano Mocatti, Giovanni Marini, Matteo Calandra

Abstract: Femtosecond pulses have been used to reveal hidden broken symmetry states and induce transitions to metastable states. However, these states are mostly transient and disappear after laser removal. Photoinduced phase transitions towards crystalline metastable states with a change of topological order are rare and difficult to predict and realize experimentally. Here, by using constrained density fu… ▽ More Femtosecond pulses have been used to reveal hidden broken symmetry states and induce transitions to metastable states. However, these states are mostly transient and disappear after laser removal. Photoinduced phase transitions towards crystalline metastable states with a change of topological order are rare and difficult to predict and realize experimentally. Here, by using constrained density functional perturbation theory and accounting for light-induced quantum anharmonicity, we show that ultra-fast lasers can permanently transform the topologically-trivial orthorhombic structure of SnSe into the topological crystalline insulating rocksalt phase via a first-order non-thermal phase transition. We describe the reaction path and evaluate the critical fluence and the possible decay channels after photoexcitation. Our simulations of the photoexcited structural and vibrational properties are in excellent agreement with recent pump-probe data in the intermediate fluence regime below the transition with an error on the curvature of the quantum free energy of the photoexcited state that is smaller than 2%. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 17 pages, 4 figures, The data underlying this study are openly available in Zenodo at https://doi.org/10.5281/zenodo.8413390

Journal ref: J. Phys. Chem. Lett. 2023, 14, XXX, 9329-9334

Infrared resonance Raman of bilayer graphene: signatures of massive fermions and band structure on the 2D peak

Authors: Lorenzo Graziotto, Francesco Macheda, Tommaso Venanzi, Guglielmo Marchese, Simone Sotgiu, Taoufiq Ouaj, Elena Stellino, Claudia Fasolato, Paolo Postorino, Marvin Metzelaars, Paul Kögerler, Bernd Beschoten, Matteo Calandra, Michele Ortolani, Christoph Stampfer, Francesco Mauri, Leonetta Baldassarre

Abstract: Few-layer graphene possesses low-energy carriers which behave as massive fermions, exhibiting intriguing properties in both transport and light scattering experiments. Lowering the excitation energy of resonance Raman spectroscopy down to 1.17 eV we target these massive quasiparticles in the split bands close to the K point. The low excitation energy weakens some of the Raman processes which are r… ▽ More Few-layer graphene possesses low-energy carriers which behave as massive fermions, exhibiting intriguing properties in both transport and light scattering experiments. Lowering the excitation energy of resonance Raman spectroscopy down to 1.17 eV we target these massive quasiparticles in the split bands close to the K point. The low excitation energy weakens some of the Raman processes which are resonant in the visible, and induces a clearer frequency-separation of the sub-structures of the resonance 2D peak in bi- and trilayer samples. We follow the excitation-energy dependence of the intensity of each sub-structure and, comparing experimental measurements on bilayer graphene with ab initio theoretical calculations, we trace back such modifications on the joint effects of probing the electronic dispersion close to the band splitting and enhancement of electron-phonon matrix elements. △ Less

Submitted 19 January, 2024; v1 submitted 6 October, 2023; originally announced October 2023.

Comments: 22 pages, 12 figures

Electronic structure and lattice dynamics of 1T-VSe$_2$: origin of the 3D-CDW

Authors: Josu Diego, D. Subires, A. H. Said, D. A. Chaney, A. Korshunov, G. Garbarino, F. Diekmann, K. Mahatha, V. Pardo, J. Strempfer, Pablo J. Bereciartua Perez, S. Francoual, C. Popescu, M. Tallarida, J. Dai, Raffaello Bianco, Lorenzo Monacelli, Matteo Calandra, A. Bosak, Francesco Mauri, K. Rossnagel, Adolfo O. Fumega, Ion Errea, S. Blanco-Canosa

Abstract: In order to characterize in detail the charge density wave (CDW) transition of 1$T$-VSe$_2$, its electronic structure and lattice dynamics are comprehensively studied by means of x-ray diffraction, angle resolved photoemission (ARPES), diffuse and inelastic x-ray scattering (IXS), and state-of-the-art first principles density functional theory calculations. Resonant elastic x-ray scattering (REXS)… ▽ More In order to characterize in detail the charge density wave (CDW) transition of 1$T$-VSe$_2$, its electronic structure and lattice dynamics are comprehensively studied by means of x-ray diffraction, angle resolved photoemission (ARPES), diffuse and inelastic x-ray scattering (IXS), and state-of-the-art first principles density functional theory calculations. Resonant elastic x-ray scattering (REXS) does not show any resonant enhancement at either V or Se K-edges, indicating that the CDW peak describes a purely structural modulation of the electronic ordering. ARPES identifies (i) a pseudogap at T$>$T$_{CDW}$, which leads to a depletion of the density of states in the $ML-M'L'$ plane at T$<$T$_{CDW}$, and (ii) anomalies in the electronic dispersion reflecting a sizable impact of phonons on it. A diffuse scattering precursor, characteristic of soft phonons, is observed at room temperature (RT) and leads to the full collapse of the low-energy phonon ($ω_1$) with propagation vector (0.25 0 -0.3) r.l.u. We show that the frequency and linewidth of this mode are anisotropic in momentum space, reflecting the momentum dependence of the electron-phonon interaction (EPI), hence demonstrating that the origin of the CDW is, to a much larger extent, due to the momentum dependence EPI with a small contribution from nesting. The pressure dependence of the $ω_1$ soft mode remains nearly constant up to 13 GPa at RT, with only a modest softening before the transition to the high-pressure monoclinic $C2/m$ phase. The wide set of experimental data are well captured by our state-of-the art first-principles anharmonic calculations with the inclusion of van der Waals (vdW) corrections in the exchange-correlation functional. The description of the electronics and dynamics of VSe$_2$ reported here adds important pieces of information to the understanding of the electronic modulations of TMDs. △ Less

Submitted 28 July, 2023; originally announced July 2023.

Electron-phonon driven charge density wave in CuTe

Authors: Marco Campetella, Giovanni Marini, Jianqiang Sky Zhou, Matteo Calandra

Abstract: The compound CuTe (vulcanite) undergoes a quasi one dimensional charge density wave (CDW) at $T< T_{\mathrm{CDW}}=335$ K with a $5\times1\times2$ periodicity. The mechanism at its origin is debated. Several theoretical works claimed that semilocal functionals are unable to describe its occurrence and ascribed its formation only to strong electron-electron interaction. Moreover, the possible role o… ▽ More The compound CuTe (vulcanite) undergoes a quasi one dimensional charge density wave (CDW) at $T< T_{\mathrm{CDW}}=335$ K with a $5\times1\times2$ periodicity. The mechanism at its origin is debated. Several theoretical works claimed that semilocal functionals are unable to describe its occurrence and ascribed its formation only to strong electron-electron interaction. Moreover, the possible role of quantum anharmonicity has not been addressed. Here, by performing quantum anharmonic calculations, we show that semilocal functionals correctly describe the occurrence of a CDW in CuTe if ultradense electron momentum grids allowing for small electronic temperatures are used. The distortion is driven by the perfect nesting among 1D Fermi surface sheets extending in the $k_y$ direction. Quantum anharmonic effects are important and tend to suppress both the distortion and $T_{\mathrm{CDW}}$. The quantum anharmonic structural minimization of the CDW phase in the generalized gradient approximation leads, however, to distorted Te-Te bond lengths in the low temperature phase that are $21\%$ of the experimental ones at $T=20$ K. This suggests that, even if the electron-electron interaction is not crucial for the mechanism of CDW formation, it is relevant to accurately describe the structural data for the low-T phase. We assess the effect of correlation on the CDW by using the DFT+U+V approximation with parameters calculated from first principles. We find that correlation enhances the Te-Te distortion, T$_{CDW}$ and the total energy gain by the distortion. △ Less

Submitted 9 July, 2023; originally announced July 2023.

Comments: 9 pages, 8 figures, to appear on Phys. Rev. B

Misfit layer compounds as ultra-tunable field effect transistors: from charge transfer control to emergent superconductivity

Authors: Ludovica Zullo, Giovanni Marini, Tristan Cren, Matteo Calandra

Abstract: Misfit layer compounds are heterostructures composed of rocksalt units stacked with few layers transition metal dichalcogenides. They host Ising superconductivity, charge density waves and good thermoelectricity. The design of misfits emergent properties is, however, hindered by the lack of a global understanding of the electronic transfer among the constituents. Here, by performing first principl… ▽ More Misfit layer compounds are heterostructures composed of rocksalt units stacked with few layers transition metal dichalcogenides. They host Ising superconductivity, charge density waves and good thermoelectricity. The design of misfits emergent properties is, however, hindered by the lack of a global understanding of the electronic transfer among the constituents. Here, by performing first principles calculations, we unveil the mechanism controlling the charge transfer and demonstrate that rocksalt units are always donor and dichalcogenides acceptors. We show that misfits behave as a periodic arrangement of ultra-tunable field effect transistors where a charging as large as 6\times10^{14} e^-cm^{-2} can be reached and controlled efficiently by the La-Pb alloying in the rocksalt. Finally, we identify a strategy to design emergent superconductivity and demonstrate its applicability in (LaSe)_{1.27}(SnSe_2)_2. Our work paves the way to the design synthesis of misfit compounds with tailored physical properties. △ Less

Submitted 6 July, 2023; originally announced July 2023.

Comments: First version of the manuscript (before revisions)

EPIq : an open-source software for the calculation of electron-phonon interaction related properties

Authors: Giovanni Marini, Guglielmo Marchese, Gianni Profeta, Jelena Sjakste, Francesco Macheda, Nathalie Vast, Francesco Mauri, Matteo Calandra

Abstract: EPIq (Electron-Phonon wannier Interpolation over k and q-points) is an open-source software for the calculation of electron-phonon interaction related properties from first principles.Acting as a post-processing tool for a density-functional perturbation theory code ( Quantum ESPRESSO ) and wannier90, EPIq exploits the localization of the deformation potential in the Wannier function basis and the… ▽ More EPIq (Electron-Phonon wannier Interpolation over k and q-points) is an open-source software for the calculation of electron-phonon interaction related properties from first principles.Acting as a post-processing tool for a density-functional perturbation theory code ( Quantum ESPRESSO ) and wannier90, EPIq exploits the localization of the deformation potential in the Wannier function basis and the stationary properties of a force-constant functional with respect to the first-order perturbation of the electronic charge density to calculate many electron-phonon related properties with high accuracy and free from convergence issues related to Brillouin zone sampling. EPIq features includes: the adiabatic and non-adiabatic phonon dispersion, superconducting properties (including the superconducting band gap in the Migdal-Eliashberg formulation), double-resonant Raman spectra and lifetime of excited carriers. The possibility to customize most of its input makes EPIq a versatile and interoperable tool. Particularly relevant is the interaction with the Stochastic Self-Consistent Harmonic Approximation (SSCHA) allowing anharmonic effects to be included in the calculation of electron-properties. The scalability offered by the Wannier representation combined with a straightforward workflow and easy-to-read input and output files make EPIq accessible to the wide condensed matter and material science communities. △ Less

Submitted 1 July, 2023; v1 submitted 27 June, 2023; originally announced June 2023.

Protection of Ising spin-orbit coupling in bulk misfit superconductors

Authors: Tomas Samuely, Darshana Wickramaratne, Martin Gmitra, Thomas Jaouen, Ondrej Šofranko, Dominik Volavka, Marek Kuzmiak, Jozef Haniš, Pavol Szabó, Claude Monney, Geoffroy Kremer, Patrick Le Fèvre, François Bertran, Tristan Cren, Shunsuke Sasaki, Laurent Cario, Matteo Calandra, Igor I. Mazin, Peter Samuely

Abstract: Low-dimensional materials have remarkable properties that are distinct from their bulk counterparts. A paradigmatic example is Ising superconductivity that occurs in monolayer materials such as NbSe2 which show a strong violation of the Pauli limit. In monolayers, this occurs due to a combination of broken inversion symmetry and spin-orbit coupling that locks the spins of the electrons out-of-plan… ▽ More Low-dimensional materials have remarkable properties that are distinct from their bulk counterparts. A paradigmatic example is Ising superconductivity that occurs in monolayer materials such as NbSe2 which show a strong violation of the Pauli limit. In monolayers, this occurs due to a combination of broken inversion symmetry and spin-orbit coupling that locks the spins of the electrons out-of-plane. Bulk NbSe2 is centrosymmetric and is therefore not an Ising superconductor. We show that bulk misfit compound superconductors, (LaSe)1.14(NbSe2) and (LaSe)1.14(NbSe2)2, comprised of monolayers and bilayers of NbSe2, exhibit unexpected Ising protection with a Pauli-limit violation comparable to monolayer NbSe2, despite formally having inversion symmetry. We study these misfit compounds using complementary experimental methods in combination with first-principles calculations. We propose theoretical mechanisms of how the Ising protection can survive in bulk materials. We show how some of these mechanisms operate in these bulk compounds due to a concerted effect of charge-transfer, defects, reduction of interlayer hopping, and stacking. This highlights how Ising superconductivity can, unexpectedly, arise in bulk materials, and possibly enable the design of bulk superconductors that are resilient to magnetic fields. △ Less

Submitted 6 April, 2023; originally announced April 2023.

Comments: 14 pages, 4 figures, 4 extended data items

Born effective charges and vibrational spectra in super and bad conducting metals

Authors: Guglielmo Marchese, Francesco Macheda, Luca Binci, Matteo Calandra, Paolo Barone, Francesco Mauri

Abstract: Interactions mediated by electron-phonon coupling are responsible for important cooperative phenomena in metals such as superconductivity and charge-density waves. The same interaction mechanisms produce strong collision rates in the normal phase of correlated metals, causing sizeable reductions of the dc conductivity and reflectivity. As a consequence, low-energy excitations like phonons, which a… ▽ More Interactions mediated by electron-phonon coupling are responsible for important cooperative phenomena in metals such as superconductivity and charge-density waves. The same interaction mechanisms produce strong collision rates in the normal phase of correlated metals, causing sizeable reductions of the dc conductivity and reflectivity. As a consequence, low-energy excitations like phonons, which are crucial for materials characterization, become visible in optical infrared spectra. A quantitative assessment of vibrational resonances requires the evaluation of dynamical Born effective charges, which quantify the coupling between macroscopic electric fields and lattice deformations. We show that the Born effective charges of metals crucially depend on the collision regime of conducting electrons. In particular, we describe, within a first principles framework, the impact of electron scattering on the infrared vibrational resonances, from the undamped, collisionless regime to the overdamped, collision-dominated limit. Our approach enables the interpretation of vibrational reflectance measurements of both super and bad conducting metals, as we illustrate for the case of strongly electron-phonon coupled superhydride H$_3$S. △ Less

Submitted 18 October, 2023; v1 submitted 1 March, 2023; originally announced March 2023.

Probing enhanced electron-phonon coupling in graphene by infrared resonance Raman spectroscopy

Authors: Tommaso Venanzi, Lorenzo Graziotto, Francesco Macheda, Simone Sotgiu, Taoufiq Ouaj, Elena Stellino, Claudia Fasolato, Paolo Postorino, Vaidotas Mišeikis, Marvin Metzelaars, Paul Kögerler, Bernd Beschoten, Camilla Coletti, Stefano Roddaro, Matteo Calandra, Michele Ortolani, Christoph Stampfer, Francesco Mauri, Leonetta Baldassarre

Abstract: We report on resonance Raman spectroscopy measurements with excitation photon energy down to 1.16 eV on graphene, to study how low-energy carriers interact with lattice vibrations. Thanks to the excitation energy close to the Dirac point at $\mathbf{K}$, we unveil a giant increase of the intensity ratio between the double-resonant 2D and 2D$^\prime$ peaks with respect to that measured in graphite.… ▽ More We report on resonance Raman spectroscopy measurements with excitation photon energy down to 1.16 eV on graphene, to study how low-energy carriers interact with lattice vibrations. Thanks to the excitation energy close to the Dirac point at $\mathbf{K}$, we unveil a giant increase of the intensity ratio between the double-resonant 2D and 2D$^\prime$ peaks with respect to that measured in graphite. Comparing with fully \textit{ab initio} theoretical calculations, we conclude that the observation is explained by an enhanced, momentum-dependent coupling between electrons and Brillouin zone-boundary optical phonons. This finding applies to two dimensional Dirac systems and has important consequences for the modeling of transport in graphene devices operating at room temperature. △ Less

Submitted 22 May, 2023; v1 submitted 2 December, 2022; originally announced December 2022.

Comments: 6 pages, 3 figures

Temperature-dependent anharmonic phonons in quantum paraelectric KTaO$_3$ by first principles and machine-learned force fields

Authors: Luigi Ranalli, Carla Verdi, Lorenzo Monacelli, Matteo Calandra, Georg Kresse, Cesare Franchini

Abstract: Understanding collective phenomena in quantum materials from first principles is a promising route toward engineering materials properties on demand and designing new functionalities. This work examines the quantum paraelectric state, an elusive state of matter characterized by the smooth saturation of the ferroelectric instability at low temperature due to quantum fluctuations associated with anh… ▽ More Understanding collective phenomena in quantum materials from first principles is a promising route toward engineering materials properties on demand and designing new functionalities. This work examines the quantum paraelectric state, an elusive state of matter characterized by the smooth saturation of the ferroelectric instability at low temperature due to quantum fluctuations associated with anharmonic phonon effects. The temperature-dependent evolution of the soft ferroelectric phonon mode in the quantum paraelectric KTaO$_3$ in the range 0-300 K is modelled by combining density functional theory (DFT) calculations with the stochastic self-consistent harmonic approximation assisted by an on-the-fly machine-learned force field. The calculated data show that including anharmonic terms is essential to stabilize the spurious imaginary ferroelectric phonon predicted by DFT, in agreement with experiments. Augmenting the DFT workflow with machine-learned force fields allows for efficient stochastic sampling of the configurational space using large supercells in a broad and dense temperature range, inaccessible by conventional ab initio protocols. This work proposes a robust computational workflow capable of accounting for collective behaviors involving different degrees of freedom and occurring at large time/length scales, paving the way for precise modeling and control of quantum effects in materials. △ Less

Submitted 24 September, 2022; originally announced September 2022.

Comments: 10 pages, 6 figures. To be published in Advanced Quantum Technologies. Editor: Luigi Ranalli

On the importance of measuring accurately LDOS maps using scanning tunneling spectroscopy in materials presenting atom-dependent charge order: the case of the correlated Pb/Si(111) single atomic layer

Authors: C. Tresca, T. Bilgeri, T. Bilgeri, G. Menard, V. Cherkez, R. Federicci, D. Longo, M. Hervé, F. Debontridder, P. David, D. Roditchev, G. Profeta, T. Cren, M. Calandra, C. Brun

Abstract: We show how to properly extract the local charge order in two-dimensional materials from scanning tunneling microscopy/spectroscopy (STM/STS) measurements. When the charge order presents spatial variations at the atomic scale inside the unit cell and is energy dependent, particular care should be taken. In such cases the use of the lock-in technique, while acquiring an STM topography in closed fee… ▽ More We show how to properly extract the local charge order in two-dimensional materials from scanning tunneling microscopy/spectroscopy (STM/STS) measurements. When the charge order presents spatial variations at the atomic scale inside the unit cell and is energy dependent, particular care should be taken. In such cases the use of the lock-in technique, while acquiring an STM topography in closed feedback loop, leads to systematically incorrect dI/dV measurements giving a false local charge order. A correct method is either to perform a constant height measurement or to perform a full grid of dI/dV(V) spectroscopies, using a bias voltage setpoint outside the material bandwidth where the local density-of-states (LDOS) is spatially homogeneous. We take as a paradigmatic example of two-dimensional material the 1/3 single-layer Pb/Si(111). As large areas of this phase cannot be grown, charge ordering in this system is not accessible to angular resolved photoemission or grazing x-ray diffraction. Previous investigations by STM/STS supplemented by {\it ab initio} Density Functional Theory (DFT) calculations concluded that this material undergoes a phase transition to a low-temperature $3\times 3$ reconstruction where one Pb atom moves up, the two remaining Pb atoms shifting down. A third STM/STS study by Adler {\it et al.} [PRL 123, 086401 (2019)] came to the opposite conclusion, i.e. that two Pb atoms move up, while one Pb atom shifts down. This latter erroneous conclusion comes from a misuse of the lock-in technique. In contrast, using a full grid of dI/dV(V) spectroscopy measurements, we show that the energy-dependent LDOS maps agree very well with state-of-the-art DFT calculations confirming the one-up two-down charge ordering. This structural and charge re-ordering in the $3\times 3$ unit cell is equally driven by electron-electron interactions and the coupling to the substrate. △ Less

Submitted 13 September, 2022; originally announced September 2022.

Comments: 11 pages, 3 figures

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

Competition between phonon mediated superconductivity and carriers localization in field-effect doped molybdenum dichalcogenides

Authors: Giovanni Marini, Matteo Calandra

Abstract: Superconductivity occurs in electrochemically doped molybdenum dichalcogenides samples thicker than four layers. While the critical temperature (Tc) strongly depends on the field effect geometry (single or double gate) and on the sample (MoS2 or MoSe2), Tc always saturates at high doping. The pairing mechanism and the complicate dependence of Tc on doping, samples and field-effect geometry are cur… ▽ More Superconductivity occurs in electrochemically doped molybdenum dichalcogenides samples thicker than four layers. While the critical temperature (Tc) strongly depends on the field effect geometry (single or double gate) and on the sample (MoS2 or MoSe2), Tc always saturates at high doping. The pairing mechanism and the complicate dependence of Tc on doping, samples and field-effect geometry are currently not understood. Previous theoretical works assumed homogeneous doping of a single layer and attributed the Tc saturation to a charge density wave instability, however the calculated values of the electron-phonon coupling in the harmonic approximation were one order of magnitude larger than the experimental estimates based on transport data. Here, by performing fully relativistic first principles calculations accounting for the sample thickness, the field-effect geometry and anharmonicity, we rule out the occurrence of charge density waves in the experimental doping range and demonstrate a suppression of one order of magnitude in the electron-phonon coupling, now in excellent agreement with transport data. By solving the anisotropic Migdal-Eliashberg equations, we explain the behaviour of Tc in different systems and geometries. From an analysis of mobility data, we propose that the Tc saturation is due to carriers localization and disorder. △ Less

Submitted 5 September, 2022; originally announced September 2022.

High pressure tuning of magnon-polarons in the layered antiferromagnet FePS$_3$

Authors: Amit Pawbake, Thomas Pelini, Alex Delhomme, Davide Romanin, Diana Vaclavkova, Gerard Martinez, Matteo Calandra, Marie-Aude Measson, Martin Veis, Marek Potemski, Milan Orlita, Clement Faugeras

Abstract: Magnetic layered materials have emerged recently as promising systems to introduce magnetism in structures based on two-dimensional (2D) materials and to investigate exotic magnetic ground states in the 2D limit. In this work, we apply high hydrostatic pressures up to P = 8.7 GPa to the bulk layered antiferromagnet FePS$_3$ to tune the collective lattice excitations (phonons) in resonance with mag… ▽ More Magnetic layered materials have emerged recently as promising systems to introduce magnetism in structures based on two-dimensional (2D) materials and to investigate exotic magnetic ground states in the 2D limit. In this work, we apply high hydrostatic pressures up to P = 8.7 GPa to the bulk layered antiferromagnet FePS$_3$ to tune the collective lattice excitations (phonons) in resonance with magnetic excitations (magnons). Close to P = 4 GPa, the magnon-phonon resonance is achieved and the strong coupling between these collective modes leads to the formation of new quasi-particles, the magnon-polarons, evidenced in our low temperature Raman scattering experiments by a particular avoided crossing behavior between the phonon and the doubly degenerate antiferromagnetic magnon. At the pressure-induced magnon-phonon resonance, three distinct coupled modes emerge. As it is mainly defined by intralayer properties, we show that the energy of the magnon is nearly pressure independent. We additionally apply high magnetic fields up to B = 30 T to fully identify and characterize the magnon excitations, and to explore the different magnon-polaron regimes for which the phonon has an energy lower-, equal to-, or higher- than the magnon energy. The description of our experimental data requires introducing a phonon-phonon coupling not taken into account in actual calculations. △ Less

Submitted 23 June, 2022; originally announced June 2022.

Comments: 11 pages + 6 pages of supplementary materials

Theory of ultrafast magnetization of non-magnetic semiconductors with localized conduction bands

Authors: Giovanni Marini, Matteo Calandra

Abstract: The magnetization of a non-magnetic semiconductor by femtosecond light pulses is crucial to achieve an all-optical control of the spin dynamics in materials and to develop faster memory devices. However, the conditions for its detection are largely unknown. In this work we identify the criteria for the observation of ultrafast magnetization and critically discuss the difficulties hindering its exp… ▽ More The magnetization of a non-magnetic semiconductor by femtosecond light pulses is crucial to achieve an all-optical control of the spin dynamics in materials and to develop faster memory devices. However, the conditions for its detection are largely unknown. In this work we identify the criteria for the observation of ultrafast magnetization and critically discuss the difficulties hindering its experimental detection. We show that ultrafast magnetization of a non magnetic semiconductor can be observed in compounds with very localized conduction band states and more delocalized valence bands, such as in the case of a p-d charge transfer gap. By using constrained and time dependent density functional theory simulations, we demonstrate that a transient ferrimagnetic state can be induced in diamagnetic semiconductor V2O5 via ultrafast pulses at realistic fluences. The ferrimagnetic state has opposite magnetic moments on vanadium (conduction) and oxygen (valence) states. Our methodology outruns the case of V2O5 as it identifies the key requirements for a computational screening of ultrafast magnetism in non-magnetic semiconductors. △ Less

Submitted 17 June, 2022; originally announced June 2022.

Comments: https://journals.aps.org/prb/accepted/d907dYc6Zfb16a70e88534261f79f72c0bb0e0d75

Highly tunable optics across a topological transition in organic polymers

Authors: Davide Romanin, Matteo Calandra, Alex W. Chin

Abstract: Controllable topological phase transitions are appealing as they allow for tunable single particle electronic properties. Here, by using state-of-the-art manybody perturbation theory techniques, we show that the topological $Z_2$ phase transition occurring in the single particle spectrum of the recently synthetized ethynylene bridged polyacene polymers is accompanied by a topological excitonic pha… ▽ More Controllable topological phase transitions are appealing as they allow for tunable single particle electronic properties. Here, by using state-of-the-art manybody perturbation theory techniques, we show that the topological $Z_2$ phase transition occurring in the single particle spectrum of the recently synthetized ethynylene bridged polyacene polymers is accompanied by a topological excitonic phase transition: the band inversion in the non-trivial phase yields real-space exciton wave functions in which electrons and holes exchange orbital characters with respect to the trivial phase. The topological excitonic phase transition results in a broad tunability of the singlet-triplet splittings, opening appealing perspectives for the occurrence of singlet fission. Finally, the flatness of the single-particle electronic structure in the topological non trivial phase leads to negatively dispersing triplet excitons in a large portion of the Brillouin zone, opening a route for spontaneously coherent energy transport at room temperature. △ Less

Submitted 1 April, 2022; originally announced April 2022.

Why Mercury is a superconductor

Authors: Cesare Tresca, Gianni Profeta, Giovanni Marini, Giovanni B. Bachelet, Antonio Sanna, Matteo Calandra, Lilia Boeri

Abstract: Despite being the oldest known superconductor, solid mercury is mysteriously absent from all current computational databases of superconductors. In this work, we present a critical study of its superconducting properties based on state-of-the-art superconducting density-functional theory. Our calculations reveal numerous anomalies in electronic and lattice properties, which can mostly be handled,… ▽ More Despite being the oldest known superconductor, solid mercury is mysteriously absent from all current computational databases of superconductors. In this work, we present a critical study of its superconducting properties based on state-of-the-art superconducting density-functional theory. Our calculations reveal numerous anomalies in electronic and lattice properties, which can mostly be handled, with due care, by modern ab-initio techniques. In particular, we highlight an anomalous role of (i) electron-electron correlations on structural properties (ii) spin-orbit coupling on the dynamical stability, and (iii) semicore $d$ levels on the effective Coulomb interaction and, ultimately, the critical temperature. △ Less

Submitted 27 November, 2021; originally announced November 2021.

Light-tunable charge density wave orders in MoTe2 and WTe2 single layers

Authors: Giovanni Marini, Matteo Calandra

Abstract: By using constrained density functional theory modeling, we demonstrate that ultrafast optical pumping unveils hidden charge orders in group VI monolayer transition metal ditellurides. We show that irradiation of the insulating 2H phases stabilizes multiple transient charge density wave orders with light-tunable distortion, periodicity, electronic structure and bandgap. Moreover, optical pumping o… ▽ More By using constrained density functional theory modeling, we demonstrate that ultrafast optical pumping unveils hidden charge orders in group VI monolayer transition metal ditellurides. We show that irradiation of the insulating 2H phases stabilizes multiple transient charge density wave orders with light-tunable distortion, periodicity, electronic structure and bandgap. Moreover, optical pumping of the semimetallic 1T' phases generates a transient charge ordered metallic phase composed of 2D diamond clusters. For each transient phase we identify the critical fluence at which it is observed and the specific optical and Raman fingerprints to directly compare with future ultrafast pump-probe experiments. Our work demonstrates that it is possible to stabilize charge density waves even in insulating 2D transition metal dichalcogenides by ultrafast irradiation. △ Less

Submitted 23 November, 2021; originally announced November 2021.

Comments: Accepted for publication in Physical Review Letters

Charge density wave in single-layer Pb/Ge(111) driven by Pb-substrate exchange interaction

Authors: Cesare Tresca, Matteo Calandra

Abstract: Single layer Pb on top of (111) surfaces of group IV semiconductors hosts charge density wave and superconductivity depending on the coverage and on the substrate. These systems are normally considered to be experimental realizations of single band Hubbard models and their properties are mostly investigated using lattice models with frozen structural degrees of freedom, although the reliability of… ▽ More Single layer Pb on top of (111) surfaces of group IV semiconductors hosts charge density wave and superconductivity depending on the coverage and on the substrate. These systems are normally considered to be experimental realizations of single band Hubbard models and their properties are mostly investigated using lattice models with frozen structural degrees of freedom, although the reliability of this approximation is unclear. Here, we consider the case of Pb/Ge(111) at 1/3 coverage, for which surface X-ray diffraction and ARPES data are available. By performing first principles calculations, we demonstrate that the non-local exchange between Pb and the substrate drives the system into a $3\times 3$ charge density wave. The electronic structure of this charge ordered phase is mainly determined by two effects: the magnitude of the Pb distortion and the large spin-orbit coupling. Finally, we show that the effect applies also to the $3\times 3$ phase of Pb/Si(111) where the Pb-substrate exchange interaction increases the bandwidth by more than a factor 1.5 with respect to DFT+U, in better agreement with STS data. The delicate interplay between substrate, structural and electronic degrees of freedom invalidates the widespread interpretation available in literature considering these compounds as physical realizations of single band Hubbard models. △ Less

Submitted 1 July, 2021; originally announced July 2021.

Comments: To appear in Phys. Rev. B

Lattice dynamics of photoexcited insulators from constrained density-functional perturbation theory

Authors: Giovanni Marini, Matteo Calandra

Abstract: We present a constrained density functional perturbation theory scheme for the calculation of structural and harmonic vibrational properties of insulators in the presence of an excited and thermalized electron-hole plasma. The method is ideal to tame ultrafast light induced structural transitions in the regime where the photocarriers thermalize faster than the lattice, the electron-hole recombinat… ▽ More We present a constrained density functional perturbation theory scheme for the calculation of structural and harmonic vibrational properties of insulators in the presence of an excited and thermalized electron-hole plasma. The method is ideal to tame ultrafast light induced structural transitions in the regime where the photocarriers thermalize faster than the lattice, the electron-hole recombination time is longer than the phonon period and the photocarrier concentration is large enough to be approximated by an electron-hole plasma. The complete derivation presented here includes total energy, forces and stress tensor, variable cell structural optimization, harmonic vibrational properties and the electron-phonon interaction. We discuss in detail the case of zone center optical phonons not conserving the number of electrons and inducing a Fermi shift in the photo-electron and hole distributions. We validate our implementation by comparing with finite differences in Te and VSe2. By calculating the evolution of the phonon spectrum of Te, Si and GaAs as a function of the fluence of the incoming laser light, we demonstrate that even at low fluences, corresponding to approximately 0.1 photocarriers per cell, the phonon spectrum is substantially modified with respect to the ground state one with new Kohn anomalies appearing and a substantial softening of zone center optical phonons. Our implementation can be efficiently used to detect reversible transient phases and irreversible structural transition induced by ultrafast light absorption. △ Less

Submitted 22 June, 2021; originally announced June 2021.

Journal ref: Phys. Rev. B 104, 144103 (2021)

Intrinsic vibrational angular momentum from non-adiabatic effects in non-collinear magnetic molecules

Authors: Oliviero Bistoni, Francesco Mauri, Matteo Calandra

Abstract: We show that in non-collinear magnetic molecules, non-adiabatic (dynamical) effects due to the electron-vibron coupling are time-reversal symmetry breaking interactions for the vibrational field. As in these systems the electronic wavefunction can not be chosen as real, a nonzero geometric vector potential (Berry connection) arises. As a result, an intrinsic nonzero vibrational angular momentum oc… ▽ More We show that in non-collinear magnetic molecules, non-adiabatic (dynamical) effects due to the electron-vibron coupling are time-reversal symmetry breaking interactions for the vibrational field. As in these systems the electronic wavefunction can not be chosen as real, a nonzero geometric vector potential (Berry connection) arises. As a result, an intrinsic nonzero vibrational angular momentum occurs even for non-degenerate modes and in the absence of external probes. The vibronic modes can then be seen as elementary quantum particles carrying a sizeable angular momentum. As a proof of concept, we demonstrate the magnitude of this topological effect by performing non-adiabatic first principles calculations on platinum clusters and by showing that these molecules host sizeable intrinsic phonon angular momenta comparable to the orbital electronic ones in itinerant ferromagnets. △ Less

Submitted 22 April, 2021; originally announced April 2021.

Comments: 5 pages, 3 figures, to appear on Physical Review Letters

The Stochastic Self-Consistent Harmonic Approximation: Calculating Vibrational Properties of Materials with Full Quantum and Anharmonic Effects

Authors: Lorenzo Monacelli, Raffaello Bianco, Marco Cherubini, Matteo Calandra, Ion Errea, Francesco Mauri

Abstract: The efficient and accurate calculation of how ionic quantum and thermal fluctuations impact the free energy of a crystal, its atomic structure, and phonon spectrum is one of the main challenges of solid state physics, especially when strong anharmonicy invalidates any perturbative approach. To tackle this problem, we present the implementation on a modular Python code of the stochastic self-consis… ▽ More The efficient and accurate calculation of how ionic quantum and thermal fluctuations impact the free energy of a crystal, its atomic structure, and phonon spectrum is one of the main challenges of solid state physics, especially when strong anharmonicy invalidates any perturbative approach. To tackle this problem, we present the implementation on a modular Python code of the stochastic self-consistent harmonic approximation method. This technique rigorously describes the full thermodyamics of crystals accounting for nuclear quantum and thermal anharmonic fluctuations. The approach requires the evaluation of the Born-Oppenheimer energy, as well as its derivatives with respect to ionic positions (forces) and cell parameters (stress tensor) in supercells, which can be provided, for instance, by first principles density-functional-theory codes. The method performs crystal geometry relaxation on the quantum free energy landscape, optimizing the free energy with respect to all degrees of freedom of the crystal structure. It can be used to determine the phase diagram of any crystal at finite temperature. It enables the calculation of phase boundaries for both first-order and second-order phase transitions from the Hessian of the free energy. Finally, the code can also compute the anharmonic phonon spectra, including the phonon linewidths, as well as phonon spectral functions. We review the theoretical framework of the stochastic self-consistent harmonic approximation and its dynamical extension, making particular emphasis on the physical interpretation of the variables present in the theory that can enlighten the comparison with any other anharmonic theory. A modular and flexible Python environment is used for the implementation, which allows for a clean interaction with other packages. We briefly present a toy-model calculation to illustrate the potential of the code. △ Less

Submitted 5 March, 2021; originally announced March 2021.

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.

Comments: Fixed typo in the Acknowledgements

Misfit layer compounds: a platform for heavily-doped two-dimensional transition metal dichalcogenides

Authors: Raphaël T. Leriche, Alexandra Palacio-Morales, Marco Campetella, Cesare Tresca, Shunsuke Sasaki, Christophe Brun, François Debontridder, Pascal David, Imad Arfaoui, Ondrej Šofranko, Tomas Samuely, Geoffroy Kremer, Claude Monney, Thomas Jaouen, Laurent Cario, Matteo Calandra, Tristan Cren

Abstract: Transition metal dichalcogenides (TMDs) display a rich variety of instabilities such as spin and charge orders, Ising superconductivity and topological properties. Their physical properties can be controlled by doping in electric double-layer field-effect transistors (FET). However, for the case of single layer NbSe$_2$, FET doping is limited to $\approx 1\times 10^{14}$ cm$^{-2}$, while a somewha… ▽ More Transition metal dichalcogenides (TMDs) display a rich variety of instabilities such as spin and charge orders, Ising superconductivity and topological properties. Their physical properties can be controlled by doping in electric double-layer field-effect transistors (FET). However, for the case of single layer NbSe$_2$, FET doping is limited to $\approx 1\times 10^{14}$ cm$^{-2}$, while a somewhat larger charge injection can be obtained via deposition of K atoms. Here, by performing ARPES, STM, quasiparticle interference measurements, and first principles calculations we show that a misfit compound formed by sandwiching NbSe$_2$ and LaSe layers behaves as a NbSe$_2$ single layer with a rigid doping of $0.55-0.6$ electrons per Nb atom or $\approx 6\times 10^{14}$ cm$^{-2}$. Due to this huge doping, the $3\times3$ charge density wave is replaced by a $2\times2$ order with very short coherence length. As a tremendous number of different misfit compounds can be obtained by sandwiching TMDs layers with rock salt or other layers, our work paves the way to the exploration of heavily doped 2D TMDs over an unprecedented wide range of doping. △ Less

Submitted 23 October, 2020; originally announced October 2020.

Phonon collapse and van der Waals melting of the 3D charge density wave of VSe$_2$

Authors: Josu Diego, A. H. Said, S. K. Mahatha, Raffaello Bianco, Lorenzo Monacelli, Matteo Calandra, Francesco Mauri, K. Rossnagel, Ion Errea, S. Blanco-Canosa

Abstract: Among transition metal dichalcogenides (TMDs), VSe$_2$ is considered to develop a purely 3-dimensional (3D) charge-density wave (CDW) at T$_{CDW}$=110 K. Here, by means of high resolution inelastic x-ray scattering (IXS), we show that the CDW transition is driven by the collapse of an acoustic mode at the critical wavevector \textit{q}$_{CDW}$= (2.25 0 0.7) r.l.u. and critical temperature T… ▽ More Among transition metal dichalcogenides (TMDs), VSe$_2$ is considered to develop a purely 3-dimensional (3D) charge-density wave (CDW) at T$_{CDW}$=110 K. Here, by means of high resolution inelastic x-ray scattering (IXS), we show that the CDW transition is driven by the collapse of an acoustic mode at the critical wavevector \textit{q}$_{CDW}$= (2.25 0 0.7) r.l.u. and critical temperature T$_{CDW}$=110 K. The softening of this mode starts to be pronounced for temperatures below 2$\times$ T$_{CDW}$ and expands over a rather wide region of the Brillouin zone, suggesting a large contribution of the electron-phonon interaction to the CDW formation. This interpretation is supported by our first principles calculations that determine a large momentum-dependence of the electron-phonon interaction, peaking at the CDW wavevector, in the presence of nesting. Fully anharmonic {\it ab initio} calculations confirm the softening of one acoustic branch at \textit{q}$_{CDW}$ as responsible for the CDW formation and show that van der Waals interactions are crucial to melt the CDW. Our work also highlights the important role of out-of-plane interactions to describe 3D CDWs in TMDs. △ Less

Submitted 16 July, 2020; originally announced July 2020.

Journal ref: Nat Commun 12, 598 (2021)

Theory of the thickness dependence of the charge density wave transition in 1T-TiTe$_2$

Authors: Jianqiang Sky Zhou, Raffaello Bianco, Lorenzo Monacelli, Ion Errea, Francesco Mauri, Matteo Calandra

Abstract: Most metallic transition metal dichalcogenides undergo charge density wave (CDW) instabilities with similar or identical ordering vectors in bulk and in single layer, albeit with different critical temperatures. Metallic 1T-TiTe$_2$ is a remarkable exception as it shows no evidence of charge density wave formation in bulk, but it displays a stable $2\times2$ reconstruction in single-layer form. Th… ▽ More Most metallic transition metal dichalcogenides undergo charge density wave (CDW) instabilities with similar or identical ordering vectors in bulk and in single layer, albeit with different critical temperatures. Metallic 1T-TiTe$_2$ is a remarkable exception as it shows no evidence of charge density wave formation in bulk, but it displays a stable $2\times2$ reconstruction in single-layer form. The mechanism for this 3D-2D crossover of the transition is still unclear, although strain from the substrate and the exchange interaction have been pointed out as possible formation mechanisms. Here, by performing non-perturbative anharmonic calculations with gradient corrected and hybrid functionals, we explain the thickness behaviour of the transition in 1T-TiTe$_2$. We demonstrate that the occurrence of the CDW in single-layer TiTe$_2$ occurs from the interplay of non-perturbative anharmonicity and an exchange enhancement of the electron-phonon interaction, larger in the single layer than in the bulk. Finally, we study the electronic and structural properties of the single-layer CDW phase and provide a complete description of its electronic structure, phonon dispersion as well as infrared and Raman active phonon modes. △ Less

Submitted 13 July, 2020; originally announced July 2020.

Comments: 11 pages, 10 pictures

Electronic structure of TiSe$_2$ from a quasi-self-consistent G$_0$W$_0$ approach

Authors: Maria Hellgren, Lucas Baguet, Matteo Calandra, Francesco Mauri, Ludger Wirtz

Abstract: In a previous work it was shown that the inclusion of exact exchange is essential for a first principles description of both the electronic- and the vibrational properties of TiSe$_2$, M. Hellgren et al. [Phys. Rev. Lett. 119, 176401 (2017)]. The $GW$ approximation provides a parameter-free description of screened exchange but is usually employed perturbatively ($G_0W_0$) making results more or le… ▽ More In a previous work it was shown that the inclusion of exact exchange is essential for a first principles description of both the electronic- and the vibrational properties of TiSe$_2$, M. Hellgren et al. [Phys. Rev. Lett. 119, 176401 (2017)]. The $GW$ approximation provides a parameter-free description of screened exchange but is usually employed perturbatively ($G_0W_0$) making results more or less dependent on the starting point. In this work, we develop a quasi-self-consistent extension of $G_0W_0$ based on the random phase approximation (RPA) and the optimized effective potential of hybrid density functional theory. This approach generates an optimal $G_0W_0$ starting-point and a hybrid exchange parameter consistent with the RPA. While self-consistency plays a minor role for systems such as Ar, BN and ScN, it is shown to be crucial for TiS$_2$ and TiSe$_2$. We find the high-temperature phase of TiSe$_2$ to be a semi-metal with a band structure in good agreement with experiment. Furthermore, the optimized hybrid functional agrees well with our previous estimate and therefore accurately reproduces the low-temperature charge density wave phase. △ Less

Submitted 2 February, 2021; v1 submitted 29 May, 2020; originally announced May 2020.

Comments: 11 pages, 5 figures

Journal ref: Phys. Rev. B 103, 075101 (2021)

Bending rigidity, sound propagation and ripples in flat graphene

Authors: Unai Aseginolaza, Josu Diego, Tommaso Cea, Raffaello Bianco, Lorenzo Monacelli, Francesco Libbi, Matteo Calandra, Aitor Bergara, Francesco Mauri, Ion Errea

Abstract: Despite many of the applications of graphene rely on its uneven stiffness and high thermal conductivity, the mechanical properties of graphene, and in general of all 2D materials, are still elusive. The harmonic theory predicts a quadratic dispersion for the flexural acoustic vibrational mode, which leads the unphysical result that long wavelength in-plane acoustic modes decay before vibrating one… ▽ More Despite many of the applications of graphene rely on its uneven stiffness and high thermal conductivity, the mechanical properties of graphene, and in general of all 2D materials, are still elusive. The harmonic theory predicts a quadratic dispersion for the flexural acoustic vibrational mode, which leads the unphysical result that long wavelength in-plane acoustic modes decay before vibrating one period, preventing the propagation of sound. The robustness of the quadratic dispersion has been questioned by arguing that the anharmonic phonon-phonon interaction linearizes it. However, this implies a divergent bending rigidity in the long wavelength limit not reproduced experimentally. Here we show that rotational symmetry protects the quadratic flexural dispersion against phonon-phonon interactions and that, consequently, the bending stiffness is non-divergent irrespective of the temperature. Our non-perturbative anharmonic calculations also determine that sound propagation coexists with a quadratic dispersion. We also show that the temperature dependence of the height fluctuations of the membrane, known as ripples, is fully determined by thermal or quantum fluctuations, but without the anharmonic suppression of their amplitude previously assumed. The universality of our conclusions reconcile experimental evidence and theory not just in graphene, but all 2D materials. △ Less

Submitted 7 June, 2023; v1 submitted 25 May, 2020; originally announced May 2020.

Journal ref: Nature Physics (2024)

Weak Dimensionality Dependence and Dominant Role of Ionic Fluctuations in the Charge-Density-Wave Transition of NbSe$_2$

Authors: Raffaello Bianco, Lorenzo Monacelli, Matteo Calandra, Francesco Mauri, Ion Errea

Abstract: Contradictory experiments have been reported about the dimensionality effect on the charge-density-wave transition in 2H NbSe$_2$. While scanning tunnelling experiments on single layers grown by molecular beam epitaxy measure a charge-density-wave transition temperature in the monolayer similar to the bulk, around 33 K, Raman experiments on exfoliated samples observe a large enhancement of the tra… ▽ More Contradictory experiments have been reported about the dimensionality effect on the charge-density-wave transition in 2H NbSe$_2$. While scanning tunnelling experiments on single layers grown by molecular beam epitaxy measure a charge-density-wave transition temperature in the monolayer similar to the bulk, around 33 K, Raman experiments on exfoliated samples observe a large enhancement of the transition temperature up to 145 K. By employing a non-perturbative approach to deal with anharmonicity, we calculate from first principles the temperature dependence of the phonon spectra both for bulk and monolayer. In both cases, the charge-density-wave transition temperature is estimated as the temperature at which the phonon energy of the mode driving the structural instability vanishes. The obtained transition temperature in the bulk is around 59 K, in rather good agreement with experiments, and it is just slightly increased in the single-layer limit to 73 K, showing the weak dependence of the transition on dimensionality. Environmental factors could motivate the disagreement between the transition temperatures reported by experiments. Our analysis also demonstrates the predominance of ionic fluctuations over electronic ones in the melting of the charge-density-wave order. △ Less

Submitted 1 September, 2020; v1 submitted 17 April, 2020; originally announced April 2020.

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

Electronic structure of pristine and Ni-substituted LaFeO$_3$ from near edge x-ray absorption fine structure experiments and first-principles simulations

Authors: Iurii Timrov, Piyush Agrawal, Xinyu Zhang, Selma Erat, Riping Liu, Artur Braun, Matteo Cococcioni, Matteo Calandra, Nicola Marzari, Daniele Passerone

Abstract: We present a joint theoretical and experimental study of the oxygen $K$-edge spectra for LaFeO$_3$ and homovalent Ni-substituted LaFeO$_3$ (LaFe$_{0.75}$Ni$_{0.25}$O$_3$), using first-principles simulations based on density-functional theory with extended Hubbard functionals and x-ray absorption near edge structure (XANES) measurements. Ground-state and excited-state XANES calculations employ Hubb… ▽ More We present a joint theoretical and experimental study of the oxygen $K$-edge spectra for LaFeO$_3$ and homovalent Ni-substituted LaFeO$_3$ (LaFe$_{0.75}$Ni$_{0.25}$O$_3$), using first-principles simulations based on density-functional theory with extended Hubbard functionals and x-ray absorption near edge structure (XANES) measurements. Ground-state and excited-state XANES calculations employ Hubbard on-site $U$ and inter-site $V$ parameters determined from first principles and the Lanczos recursive method to obtain absorption cross sections, which allows for a reliable description of XANES spectra in transition-metal compounds in a very broad energy range, with an accuracy comparable to that of hybrid functionals but at a substantially lower cost. We show that standard gradient-corrected exchange-correlation functionals fail in capturing accurately the electronic properties of both materials. In particular, for LaFe$_{0.75}$Ni$_{0.25}$O$_3$ they do not reproduce its semiconducting behaviour and provide a poor description of the pre-edge features at the O $K$ edge. The inclusion of Hubbard interactions leads to a drastic improvement, accounting for the semiconducting ground state of LaFe$_{0.75}$Ni$_{0.25}$O$_3$ and for a good agreement between calculated and measured XANES spectra. We show that the partial substitution of Fe for Ni affects the conduction-band bottom by generating a strongly hybridized O($2p$)-Ni($3d$) minority-spin empty electronic state. The present work, based on a consistent correction of self-interaction errors, outlines the crucial role of extended Hubbard functionals to describe the electronic structure of complex transition-metal oxides such as LaFeO$_3$ and LaFe$_{0.75}$Ni$_{0.25}$O$_3$ and paves the way to future studies on similar systems. △ Less

Submitted 24 August, 2020; v1 submitted 8 April, 2020; originally announced April 2020.

Journal ref: Phys. Rev. Research 2, 033265 (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.

Hybrid-functional electronic structure of multilayer graphene

Authors: Marco Campetella, Nguyen Minh Nguyen, Jacopo Baima, Lorenzo Maschio, Francesco Mauri, Matteo Calandra

Abstract: Multilayer graphene with rhombohedral and Bernal stacking are supposed to be metallic, as predicted by density functional theory calculations using semi-local functionals. However recent angular resolved photoemission and transport data have questioned this point of view. In particular, rhombohedral flakes are suggested to be magnetic insulators. Bernal flakes composed of an even number of layers… ▽ More Multilayer graphene with rhombohedral and Bernal stacking are supposed to be metallic, as predicted by density functional theory calculations using semi-local functionals. However recent angular resolved photoemission and transport data have questioned this point of view. In particular, rhombohedral flakes are suggested to be magnetic insulators. Bernal flakes composed of an even number of layers are insulating, while those composed of an odd number of layers are pseudogapped. Here, by systematically benchmarking with plane waves codes, we develop very accurate all-electron Gaussian basis sets for graphene multilayers. We find that, in agreement with our previous calculations, rhombohedral stacked multilayer graphene are gapped for and magnetic. However, the valence band curvature and the details of the electronic structure depend crucially on the basis set. Only substantially extended basis sets are able to correctly reproduce the effective mass of the valence band top at the K point, while the popular POB-TZVP basis set leads to a severe overestimation. In the case of Bernal stacking, we show that exact exchange gaps the flakes composed by four layers and opens pseudogaps for N = 3, 6, 7, 8. However, the gap or pseudogap size and its behaviour as a function of thickness are not compatible with experimental data. Moreover, hybrid functionals lead to a metallic solution for 5 layers and a magnetic ground state for 5, 6 and 8 layers. Magnetism is very weak with practically no effect on the electronic structure and the magnetic moments are mostly concentrated in the central layers. Our hybrid functional calculations on trilayer Bernal graphene multilayers are in excellent agreement with non-magnetic GW calculations. For thicker multilayers, our calculations are a benchmark for manybody theoretical modeling of the low energy electronic structure. △ Less

Submitted 17 January, 2020; originally announced January 2020.

Comments: 13 pages, 12 figures

Journal ref: Phys. Rev. B 101, 165437 (2020)

Thermal conductivity of Bi$_2$Se$_3$ from bulk to thin films: theory and experiment

Authors: Lorenzo Paulatto, Danièle Fournier, Massimiliano Marangolo, Mahmoud Eddrief, Paola Atkinson, Matteo Calandra

Abstract: We calculate the lattice-driven in-plane $(κ_{\parallel})$ and out-of-plane $(κ_{\perp})$ thermal conductivities of Bi$_2$Se$_3$ bulk, and of films of different thicknesses, using the Boltzmann equation with phonon scattering times obtained from anharmonic third order density functional perturbation theory. We compare our results for the lattice component of the thermal conductivity with publish… ▽ More We calculate the lattice-driven in-plane $(κ_{\parallel})$ and out-of-plane $(κ_{\perp})$ thermal conductivities of Bi$_2$Se$_3$ bulk, and of films of different thicknesses, using the Boltzmann equation with phonon scattering times obtained from anharmonic third order density functional perturbation theory. We compare our results for the lattice component of the thermal conductivity with published data for $κ_{\parallel}$ on bulk samples and with our room-temperature thermoreflectance measurements of $κ_{\perp}$ on films of thickness (L) ranging from 18~nm to 191~nm, where the lattice component has been extracted via the Wiedemann-Franz law. Ab-initio theoretical calculations on bulk samples, including an effective model to account for finite sample thickness and defect scattering, compare favorably both for the bulk case (from literature) and thin films (new measurements). In the low-T limit the theoretical in-plane lattice thermal conductivity of bulk Bi$_2$Se$_3$ agrees with previous measurements by assuming the occurrence of intercalated Bi$_2$ layer defects. The measured thermal conductivity monotonically decreases by reducing $L$, its value is $κ_{\perp}\approx 0.39\pm 0.08$~W/m$\cdot$K for $L=18$ nm and $κ_{\perp}=0.68\pm0.14$~W/m$\cdot$K for $L=191$ nm. We show that the decrease of room-temperature $κ_{\perp}$ in Bi$_2$Se$_3$ thin films as a function of sample thickness can be explained by the incoherent scattering of out-of-plane momentum phonons with the film surface. Our work outlines the crucial role of sample thinning in reducing the out-of-plane thermal conductivity. △ Less

Submitted 23 March, 2020; v1 submitted 13 December, 2019; originally announced December 2019.

Black metal hydrogen above 360 GPa driven by proton quantum fluctuations

Authors: Lorenzo Monacelli, Ion Errea, Matteo Calandra, Francesco Mauri

Abstract: Production of metallic hydrogen is one of the top three open quests of physics. Recent low-temperature experiments report different metallization pressures, varying from 360GPa to 490GPa. In this work, we simulate structural properties, vibrational Raman, IR and optical spectra of hydrogen phase III accounting for proton quantum effects. We demonstrate that nuclear quantum fluctuations downshift… ▽ More Production of metallic hydrogen is one of the top three open quests of physics. Recent low-temperature experiments report different metallization pressures, varying from 360GPa to 490GPa. In this work, we simulate structural properties, vibrational Raman, IR and optical spectra of hydrogen phase III accounting for proton quantum effects. We demonstrate that nuclear quantum fluctuations downshift the vibron frequencies by 25%, introduce a broad line-shape in the Raman spectra, and reduce the optical gap by 3eV. We show that hydrogen metallization occurs at 380GPa in phase III due to band overlap, in good agreement with transport data. By simulating the optical properties, we predict this state to be a peculiar black metal, transparent in the IR. The transparent window closes at 450GPa, but the reflectivity remains low, discarding phase III as the shiny metal observed at 490GPa. We predict the conductivity onset to increase by 70GPa and the transparent window to increase by 1.3eV when replacing hydrogen by deuterium at 0K, underlining that metallization is driven by quantum fluctuations and is thus isotope dependent. We show how hydrogen acquires metallic features (conductivity and brightness) at different pressures, explaining the apparent contradictions in existing experimental scenarios. △ Less

Submitted 28 January, 2020; v1 submitted 11 December, 2019; originally announced December 2019.

Gate Tunable Magnetism and Giant Magnetoresistance in ABC-stacked Few-Layer Graphene

Authors: Yongjin Lee, Shi Che, Jairo Velasco Jr., David Tran, Jacopo Baima, Francesco Mauri, Matteo Calandra, Marc Bockrath, Chun Ning Lau

Abstract: Magnetism is a prototypical phenomenon of quantum collective state, and has found ubiquitous applications in semiconductor technologies such as dynamic random access memory (DRAM). In conventional materials, it typically arises from the strong exchange interaction among the magnetic moments of d- or f-shell electrons. Magnetism, however, can also emerge in perfect lattices from non-magnetic elemen… ▽ More Magnetism is a prototypical phenomenon of quantum collective state, and has found ubiquitous applications in semiconductor technologies such as dynamic random access memory (DRAM). In conventional materials, it typically arises from the strong exchange interaction among the magnetic moments of d- or f-shell electrons. Magnetism, however, can also emerge in perfect lattices from non-magnetic elements. For instance, flat band systems with high density of states (DOS) may develop spontaneous magnetic ordering, as exemplified by the Stoner criterion. Here we report tunable magnetism in rhombohedral-stacked few-layer graphene (r-FLG). At small but finite doping (n~10^11 cm-2), we observe prominent conductance hysteresis and giant magnetoconductance that exceeds 1000% as a function of magnetic fields. Both phenomena are tunable by density and temperature, and disappears for n>10^12 cm-2 or T>5K. These results are confirmed by first principles calculations, which indicate the formation of a half-metallic state in doped r-FLG, in which the magnetization is tunable by electric field. Our combined experimental and theoretical work demonstrate that magnetism and spin polarization, arising from the strong electronic interactions in flat bands, emerge in a system composed entirely of carbon atoms. The electric field tunability of magnetism provides promise for spintronics and low energy device engineering. △ Less

Submitted 11 November, 2019; originally announced November 2019.

Anharmonic melting of the charge density wave in single-layer TiSe2

Authors: Jianqiang Sky Zhou, Lorenzo Monacelli, Raffaello Bianco, Ion Errea, Francesco Mauri, Matteo Calandra

Abstract: Low dimensional systems with a vanishing band-gap and a large electron-hole interaction have been proposed to be unstable towards exciton formation. As the exciton binding energy increases in low dimension, conventional wisdom suggests that excitonic insulators should be more stable in 2D than in 3D. Here we study the effects of the electron-hole interaction and anharmonicity in single-layer TiSe2… ▽ More Low dimensional systems with a vanishing band-gap and a large electron-hole interaction have been proposed to be unstable towards exciton formation. As the exciton binding energy increases in low dimension, conventional wisdom suggests that excitonic insulators should be more stable in 2D than in 3D. Here we study the effects of the electron-hole interaction and anharmonicity in single-layer TiSe2. We find that, contrary to the bulk case and to the generally accepted picture, the electron-hole exchange interaction is much smaller in 2D than in 3D and it has negligible effects on phonon spectra. By calculating anharmonic phonon spectra within the stochastic self-consistent harmonic approximation, we obtain TCDW = 440K for an isolated and undoped single-layer and TCDW = 364K for an electron-doping n = 4.6 x 10^13 cm^{-2} , close to the experimental result of 200-280K on supported samples. Our work demonstrates that anharmonicity and doping melt the charge density wave in single-layer TiSe2. △ Less

Submitted 28 October, 2019; originally announced October 2019.

Comments: 5 pages, 4 figures

The stability in the most external region of the Oort Cloud: The evolution of the ejected comets

Authors: J. A. Correa-Otto, M. F. Calandra

Abstract: In this paper, we present a study about the dynamical effects of the Galaxy on the external region of the Oort Cloud. The aims of this paper are: i) to determine an outer limit for the Oort Cloud; and ii) to analyse the dynamical behaviour of the most external objects of the Cloud and how they are ejected from the Solar System. This is undertaken by following the temporal evolution of massless tes… ▽ More In this paper, we present a study about the dynamical effects of the Galaxy on the external region of the Oort Cloud. The aims of this paper are: i) to determine an outer limit for the Oort Cloud; and ii) to analyse the dynamical behaviour of the most external objects of the Cloud and how they are ejected from the Solar System. This is undertaken by following the temporal evolution of massless test particles in the Galactic environment of the solar neighbourhood. Here we show that the effect of the perturbations from the Galactic tide in the particles is similar to that find for the evolution of wide binary stars population. Moreover, in the Oort Cloud we found a dynamical structure around 10 5 au conformed by objects unbound of the Sun. This structure allows us to define a transition region of stability and an outer boundary for the Oort Cloud, and it is also in agreement with previous results about the disruption of wide binary stars. △ Less

Submitted 23 September, 2019; originally announced September 2019.

Comments: Accepted for publishing in MNRAS. 13 pages, 9 figures, 2 tables

Electric field exfoliation and high-T$_{\text{C}}$ superconductivity in field-effect hole-doped hydrogenated diamond (111)

Authors: Davide Romanin, Thibault Sohier, Dario Daghero, Francesco Mauri, Renato S. Gonnelli, Matteo Calandra

Abstract: We investigate the possible occurrence of field-effect induced superconductivity in the hydrogenated $(111)$ diamond surface by first-principles calculations. By computing the band alignment between bulk diamond and the hydrogenated surface we show that the electric field exfoliates the sample, separating the electronic states at the valence band top from the bulk projected ones. At the hole dopin… ▽ More We investigate the possible occurrence of field-effect induced superconductivity in the hydrogenated $(111)$ diamond surface by first-principles calculations. By computing the band alignment between bulk diamond and the hydrogenated surface we show that the electric field exfoliates the sample, separating the electronic states at the valence band top from the bulk projected ones. At the hole doping values considered here, ranging from $n=2.84\times10^{13}$cm$^{-2}$ to $n=6\times 10^{14} $ cm$^{-2}$, the valence band top is composed of up to three electronic bands hosting holes with different effective masses. These bands resemble those of the undoped surface, but they are heavily modified by the electric field and differ substantially from a rigid doping picture. We calculate superconducting properties by including the effects of charging of the slab and of the electric field on the structural properties, electronic structure, phonon dispersion and electron-phonon coupling. We find that at doping as large as $n=6\times 10^{14} $ cm$^{-2}$, the electron-phonon interaction is $λ=0.81$ and superconductivity emerges with $T_{\text{C}}\approx 29-36$K. Superconductivity is mostly supported by in-plane diamond phonon vibrations and to a lesser extent by some out-of-plane vibrations. The relevant electron-phonon scattering processes involve both intra and interband scattering so that superconductivity is multiband in nature. △ Less

Submitted 5 August, 2019; originally announced August 2019.

Anomalous high-temperature superconductivity in YH$_6$

Authors: Ivan A. Troyan, Dmitrii V. Semenok, Alexander G. Kvashnin, Andrey V. Sadakov, Oleg A. Sobolevskiy, Vladimir M. Pudalov, Anna G. Ivanova, Vitali B. Prakapenka, Eran Greenberg, Alexander G. Gavriliuk, Viktor V. Struzhkin, Aitor Bergara, Ion Errea, Raffaello Bianco, Matteo Calandra, Francesco Mauri, Lorenzo Monacelli, Ryosuke Akashi, Artem R. Oganov

Abstract: Pressure-stabilized hydrides are a new rapidly growing class of high-temperature superconductors which is believed to be described within the conventional phonon-mediated mechanism of coupling. Here we report the synthesis of yttrium hexahydride Im3m-YH$_6$ that demonstrates the superconducting transition with T$_c$ = 224 K at 166 GPa, much lower than the theoretically predicted (>270 K). The meas… ▽ More Pressure-stabilized hydrides are a new rapidly growing class of high-temperature superconductors which is believed to be described within the conventional phonon-mediated mechanism of coupling. Here we report the synthesis of yttrium hexahydride Im3m-YH$_6$ that demonstrates the superconducting transition with T$_c$ = 224 K at 166 GPa, much lower than the theoretically predicted (>270 K). The measured upper critical magnetic field B$_c$$_2$(0) of YH$_6$ was found to be 116-158 T, which is 2-2.5 times larger than the calculated value. A pronounced shift of T$_c$ in yttrium deuteride YD$_6$ with the isotope coefficient 0.4 supports the phonon-assisted superconductivity. Current-voltage measurements showed that the critical current I$_c$ and its density J$_c$ may exceed 1.75 A and 3500 A/mm$^2$ at 0 K, respectively, which is comparable with the parameters of commercial superconductors, such as NbTi and YBCO. The superconducting density functional theory (SCDFT) and anharmonic calculations suggest unusually large impact of the Coulomb repulsion in this compound. The results indicate notable departures of the superconducting properties of the discovered YH$_6$ from the conventional Migdal-Eliashberg and Bardeen-Cooper-Schrieffer theories. △ Less

Submitted 12 June, 2020; v1 submitted 5 August, 2019; originally announced August 2019.

Comments: arXiv admin note: text overlap with arXiv:1902.10206

Journal ref: Advanced Materials 2021

Quantum Crystal Structure in the 250 K Superconducting Lanthanum Hydride

Authors: Ion Errea, Francesco Belli, Lorenzo Monacelli, Antonio Sanna, Takashi Koretsune, Terumasa Tadano, Raffaello Bianco, Matteo Calandra, Ryotaro Arita, Francesco Mauri, José A. Flores-Livas

Abstract: The discovery of superconductivity at 200 K in the hydrogen sulfide system at large pressures [1] was a clear demonstration that hydrogen-rich materials can be high-temperature superconductors. The recent synthesis of LaH$_{10}$ with a superconducting critical temperature (T$_{\text{c}}$) of 250 K [2,3] places these materials at the verge of reaching the long-dreamed room-temperature superconducti… ▽ More The discovery of superconductivity at 200 K in the hydrogen sulfide system at large pressures [1] was a clear demonstration that hydrogen-rich materials can be high-temperature superconductors. The recent synthesis of LaH$_{10}$ with a superconducting critical temperature (T$_{\text{c}}$) of 250 K [2,3] places these materials at the verge of reaching the long-dreamed room-temperature superconductivity. Electrical and x-ray diffraction measurements determined a weakly pressure-dependent T$_{\text{c}}$ for LaH$_{10}$ between 137 and 218 gigapascals in a structure with a face-centered cubic (fcc) arrangement of La atoms [3]. Here we show that quantum atomic fluctuations stabilize in all this pressure range a high-symmetry Fm-3m crystal structure consistent with experiments, which has a colossal electron-phonon coupling of $λ\sim3.5$. Even if ab initio classical calculations neglecting quantum atomic vibrations predict this structure to distort below 230 GPa yielding a complex energy landscape with many local minima, the inclusion of quantum effects simplifies the energy landscape evidencing the Fm-3m as the true ground state. The agreement between the calculated and experimental T$_{\text{c}}$ values further supports this phase as responsible for the 250 K superconductivity. The relevance of quantum fluctuations in the energy landscape found here questions many of the crystal structure predictions made for hydrides within a classical approach that at the moment guide the experimental quest for room-temperature superconductivity [4,5,6]. Furthermore, quantum effects reveal crucial to sustain solids with extraordinary electron-phonon coupling that may otherwise be unstable [7]. △ Less

Submitted 27 July, 2019; originally announced July 2019.

Comments: 20 Pages, 14 Figures, 4 Tables

Strong anharmonicity and high thermoelectric efficiency in high temperature SnS from first-principles

Authors: Unai Aseginolaza, Raffaello Bianco, Lorenzo Monacelli, Lorenzo Paulatto, Matteo Calandra, Francesco Mauri, Aitor Bergara, Ion Errea

Abstract: SnS and SnSe are isoelectronic materials with a common phase diagram. Recently, SnSe was found to be the most efficient intrinsic thermoelectric material in its high-temperature Cmcm phase above 800 K. Making use of first-principles calculations, here we show that the electronic and vibrational properties of both materials are very similar in this phase and, consequently, SnS is also expected to h… ▽ More SnS and SnSe are isoelectronic materials with a common phase diagram. Recently, SnSe was found to be the most efficient intrinsic thermoelectric material in its high-temperature Cmcm phase above 800 K. Making use of first-principles calculations, here we show that the electronic and vibrational properties of both materials are very similar in this phase and, consequently, SnS is also expected to have a high thermoelectric figure of merit at high temperature in its Cmcm phase. In fact, the electronic power factor and lattice thermal conductivity are comparable for both materials, which ensures a similar figure of merit. As in the case of SnSe, the vibrational properties of SnS in the Cmcm phase are far from trivial and are dominated by huge anharmonic effects. Its phonon spectra are strongly renormalized by anharmonicity and the spectral functions of some particular in-plane modes depict anomalous non-lorentzian profiles. Finally, we show that non-perturbative anharmonic effects in the third-order force-constants are crucial in the calculation of the lattice thermal conductivity. Our results motivate new experiments in the high temperature regime to measure the figure of merit of SnS. △ Less

Submitted 5 June, 2019; originally announced June 2019.

Journal ref: Phys. Rev. B 100, 214307 (2019)

Quantum effects in muon spin spectroscopy within the stochastic self-consistent harmonic approximation

Authors: Ifeanyi John Onuorah, Pietro Bonfà, Roberto De Renzi, Lorenzo Monacelli, Francesco Mauri, Matteo Calandra, Ion Errea

Abstract: In muon spin rotation experiments the positive implanted muon vibrates with large zero point amplitude by virtue of its light mass. Quantum mechanical calculations of the host material usually treat the muon as a point impurity, ignoring this large vibrational amplitude. As a first order correction, the muon zero point motion is usually described within the harmonic approximation, despite the larg… ▽ More In muon spin rotation experiments the positive implanted muon vibrates with large zero point amplitude by virtue of its light mass. Quantum mechanical calculations of the host material usually treat the muon as a point impurity, ignoring this large vibrational amplitude. As a first order correction, the muon zero point motion is usually described within the harmonic approximation, despite the large anharmonicity of the crystal potential. Here we apply the stochastic self-consistent harmonic approximation, a quantum variational method devised to include strong anharmonic effects in total energy and vibrational frequency calculations, in order to overcome these limitations and provide an accurate ab initio description of the quantum nature of the muon. We applied this full quantum treatment to the calculation of the muon contact hyperfine field in textbook-case metallic systems, such as Fe, Ni, Co including MnSi and MnGe, significantly improving agreement with experiments. Our results show that muon vibrational frequencies are strongly renormalized by anharmonicity. Finally, in contrast to the harmonic approximation, we show that including quantum anharmonic fluctuations, the muon stabilizes at the octahedral site in bcc Fe. △ Less

Submitted 19 July, 2019; v1 submitted 26 April, 2019; originally announced April 2019.

Comments: 10 pages

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

Electronic structure and magnetic properties of few-layer Cr$_2$Ge$_2$Te$_6$: the key role of nonlocal electron-electron interaction effects

Authors: Guido Menichetti, Matteo Calandra, Marco Polini

Abstract: Atomically-thin magnetic crystals have been recently isolated experimentally, greatly expanding the family of two-dimensional materials. In this Article we present an extensive comparative analysis of the electronic and magnetic properties of ${\rm Cr}_2{\rm Ge}_2{\rm Te}_6$, based on density functional theory (DFT). We first show that the often-used ${\rm DFT}+U$ approaches fail in predicting the… ▽ More Atomically-thin magnetic crystals have been recently isolated experimentally, greatly expanding the family of two-dimensional materials. In this Article we present an extensive comparative analysis of the electronic and magnetic properties of ${\rm Cr}_2{\rm Ge}_2{\rm Te}_6$, based on density functional theory (DFT). We first show that the often-used ${\rm DFT}+U$ approaches fail in predicting the ground-state properties of this material in both its monolayer and bilayer forms, and even more spectacularly in its bulk form. In the latter case, the fundamental gap {\it decreases} by increasing the Hubbard-$U$ parameter, eventually leading to a metallic ground state for physically relevant values of $U$, in stark contrast with experimental data. On the contrary, the use of hybrid functionals, which naturally take into account nonlocal exchange interactions between all orbitals, yields good account of the available ARPES experimental data. We then calculate all the relevant exchange couplings (and the magneto-crystalline anisotropy energy) for monolayer, bilayer, and bulk ${\rm Cr}_2{\rm Ge}_2{\rm Te}_6$ with a hybrid functional, with super-cells containing up to $270$ atoms, commenting on existing calculations with much smaller super-cell sizes. In the case of bilayer ${\rm Cr}_2{\rm Ge}_2{\rm Te}_6$, we show that two distinct intra-layer second-neighbor exchange couplings emerge, a result which, to the best of our knowledge, has not been noticed in the literature. △ Less

Submitted 24 April, 2019; originally announced April 2019.

Comments: 13 pages, 6 figures, 3 tables

Journal ref: 2D Mater. 6 045042 (2019)

Competition between exchange-driven dimerization and magnetism in diamond(111)

Authors: Betül Pamuk, Matteo Calandra

Abstract: Strong electron-electron interaction in ultraflat edge states can be responsible for correlated phases of matter, such as magnetism, charge density wave or superconductivity. Here we consider the diamond(111) surface that, after Pandey reconstruction, presents zig-zag carbon chains, generating a flat surface band. By performing full structural optimization with hybrid functionals and neglecting sp… ▽ More Strong electron-electron interaction in ultraflat edge states can be responsible for correlated phases of matter, such as magnetism, charge density wave or superconductivity. Here we consider the diamond(111) surface that, after Pandey reconstruction, presents zig-zag carbon chains, generating a flat surface band. By performing full structural optimization with hybrid functionals and neglecting spin polarization, we find that a substantial dimerization ($0.090$ Å / $0.076$ Å bond disproportionation in the PBE0/HSE06) occurs on the chains; a structural effect absent in calculations based on the LDA/GGA functionals. This dimerization is the primary mechanism for the opening of an insulating gap in the absence of spin polarization. The single-particle direct gap is $1.7$ eV ($1.0$ eV) in the PBE0 (HSE06), comparable with the experimental optical gap of $1.47$ eV, and on the larger(smaller) side of the estimated experimental single particle gap window of 1.57-1.87 eV, after inclusion of excitonic effects. However, by including spin polarization in the calculation, we find that the exchange interaction stabilizes a different ground state, undimerized, with no net magnetization and ferrimagnetic along the Pandey $π$-chains with magnetic moments as large as $0.2-0.3~μ_B$ in the PBE0. The direct single-particle band gap in the equal spin-channel is approximately $2.2$ eV ($1.5$ eV) with the PBE0 (HSE06) functional. Our work is relevant for systems with flat bands in general and wherever the interplay between structural, electronic and magnetic degrees of freedom is crucial, as in twisted bilayer graphene, IVB atoms on IVB(111) surfaces such as Pb/Si(111) or molecular crystals. △ Less

Submitted 21 March, 2019; originally announced March 2019.

Comments: 8 pages, 6 figures