Researcher profile

Francesco Mauri

Francesco Mauri contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

13 published item(s)

preprint2022arXiv

Many-body Green's function approach to lattice thermal transport

Recent progress in understanding thermal transport in complex crystals has highlighted the prominent role of heat conduction mediated by interband tunneling processes, which emerge between overlapping phonon bands (i.e. with energy differences smaller than their broadenings). These processes have recently been described in different ways, relying on the Wigner or Green-Kubo formalism, leading to apparently different results which question the definition of the heat-current operator. Here, we implement a full quantum approach based on the Kubo formula, elucidating analogies and differences with the recently introduced Wigner or Green-Kubo formulations, and extending the description of thermal transport to the overdamped regime of atomic vibrations, where the phonon quasiparticle picture breaks down. We rely on first-principles calculations on complex crystals with ultralow conductivity to compare numerically the thermal conductivity obtained within the aforementioned approaches, showing that at least in the quasiparticle regime the differences are negligible for practical applications.

preprint2022arXiv

Non-perturbative Green's function method to determine the electronic spectral function due to electron-phonon interactions: Application to a graphene model from weak to strong coupling

In solid state physics, the electron-phonon interaction (EPI) is central to many phenomena. The theory of the renormalization of electronic properties due to EPIs became well established with the theory of Allen-Heine-Cardona, usually applied to second order in perturbation theory (P2). However, this is only valid in the weak coupling regime, while strong EPIs have been reported in many materials. Although non-perturbative (NP) methods have started to arise in the last years, they are usually not well justified, and it is not clear to what degree they reproduce the exact theory. To address this issue, we present a stochastic approach for the evaluation of the non-perturbative interacting Green's function in the adiabatic limit, and show it is equivalent to the Feynman expansion to all orders in the perturbation. Also, by defining a self-energy, we can reduce the effect of broadening needed in numerical calculations, improving convergence in the supercell size. In addition, we clarify whether it is better to average the Green's function or self-energy. Then we apply the method to a graphene tight-binding model, and obtain several interesting results: (i) The Debye-Waller term, which is normally neglected, does affect the change of the Fermi velocity. (ii) The P2 and NP self-energies differ even at room temperature for some k-points, raising the question of how well P2 works in other materials. (iii) Close to the Dirac point, positive and negative energy peaks merge. (iv) In the strong coupling regime, a peak appears at energy E=0, which is consistent with previous works on disorder and localization in graphene. (v) The spectral function becomes more asymmetric at stronger coupling and higher temperatures. Finally, in the Appendix we show that the method has better convergent properties when the coupling is strong relative to when it is weak, and discuss other technical aspects.

preprint2022arXiv

The microscopic origin of anomalous properties of ice relies on the strong quantum anharmonic regime of atomic vibrations

Water ice is a unique material presenting intriguing physical properties, like negative thermal expansion and anomalous volume isotope effect (VIE). They arise from the interplay between weak hydrogen bonds and nuclear quantum fluctuations, making theoretical calculations challenging. Here, we employ the stochastic self-consistent harmonic approximation (SSCHA) to investigate how thermal and quantum fluctuations affect the physical properties of ice XI ab initio. Regarding the anomalous VIE, our work reveals that quantum effects on hydrogen are so strong to be in a nonlinear regime: when progressively increasing the mass of hydrogen from protium to infinity (classical limit), the volume firstly expands and then contracts, with a maximum slightly above the mass of tritium. We observe an anharmonic renormalization of about 10% in the bending and stretching phonon frequencies probed in IR and Raman experiments. For the first time, we report an accurate comparison of the low energy phonon dispersion with the experimental data, possible only thanks to high-level accuracy in the electronic correlation and nuclear quantum and thermal fluctuations, paving the way for the study of thermal transport in ice from first principles and the simulation of ice under pressure.

preprint2022arXiv

Wigner formulation of thermal transport in solids

Two different heat-transport mechanisms are discussed in solids: in crystals, heat carriers propagate and scatter particle-like as described by Peierls' formulation of the Boltzmann transport equation for phonon wavepackets. In glasses, instead, carriers behave wave-like, diffusing via a Zener-like tunneling between quasi-degenerate vibrational eigenstates, as described by the Allen-Feldman equation. Recently, it has been shown that these two conduction mechanisms emerge from a Wigner transport equation, which unifies and extends the Peierls-Boltzmann and Allen-Feldman formulations, allowing to describe also complex crystals where particle-like and wave-like conduction mechanisms coexist. Here, we discuss the theoretical foundations of such transport equation as is derived from the Wigner phase-space formulation of quantum mechanics, elucidating how the interplay between disorder, anharmonicity, and the quantum Bose-Einstein statistics of atomic vibrations determines thermal conductivity. This Wigner formulation argues for a preferential phase convention for the dynamical matrix in the reciprocal Bloch representation and related off-diagonal velocity operator's elements; such convention is the only one yielding a conductivity which is invariant with respect to the non-unique choice of the crystal's unit cell and is size-consistent. We rationalize the conditions determining the crossover from particle-like to wave-like heat conduction, showing that phonons below the Ioffe-Regel limit (i.e. with a mean free path shorter than the interatomic spacing) contribute to heat transport due to their wave-like capability to interfere and tunnel. Finally, we show that the present approach overcomes the failures of the Peierls-Boltzmann formulation for materials with ultralow or glass-like thermal conductivity.

preprint2021arXiv

Ab-initio energetics of graphite and multilayer graphene: stability of Bernal versus rhombohedral stacking

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.

preprint2021arXiv

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

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.

preprint2020arXiv

Black metal hydrogen above 360 GPa driven by proton quantum fluctuations

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.

preprint2020arXiv

Hybrid-functional electronic structure of multilayer graphene

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.

preprint2020arXiv

Long-range rhombohedral-stacked graphene through shear

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2019arXiv

Quantum Crystal Structure in the 250 K Superconducting Lanthanum Hydride

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