Source author record

Francesco Mauri

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

ResearcherUnclaimed source record

cond-mat.mtrl-sci cond-mat.mes-hall cond-mat.supr-con cond-mat.str-el mtrl-th physics.chem-ph cond-mat cond-mat.dis-nn

Catalog footprint

What is connected

63works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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

Building this map preview

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

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

preprint2016arXiv

Anharmonic effects in atomic hydrogen: superconductivity and lattice dynamical stability

We present first-principles calculations of metallic atomic hydrogen in the 400-600 GPa pressure range in a tetragonal structure with space group $I4_1/amd$, which is predicted to be its first atomic phase. Our calculations show a band structure close to the free-electron-like limit due to the high electronic kinetic energy induced by pressure. Bands are properly described even in the independent electron approximation fully neglecting the electron-electron interaction. Linear-response harmonic calculations show a dynamically stable phonon spectrum with marked Kohn anomalies. Even if the electron-electron interaction has a minor role in the electronic bands, the inclusion of electronic ex- change and correlation in the density response is essential to obtain a dynamically stable structure. Anharmonic effects, which are calculated within the stochastic self-consistent harmonic approxima- tion, harden high-energy optical modes and soften transverse acoustic modes up to a 20% in energy. Despite a large impact of anharmonicity has been predicted in several high-pressure hydrides, here the superconducting critical temperature is barely affected by anharmonicity, as it is lowered from its harmonic 318 K value only to 300 K at 500 GPa. We atribute the small impact of anharmoncity on superconductivity to the absence of softened optical modes and the fairly uniform distribution of the electron-phonon coupling among the vibrational modes.

preprint2016arXiv

Determination of scattering time and of valley occupation in transition-metal dichalcogenides doped by field effect

The transition-metal dichalcogenides have attracted a lot of attention as a possible stepping-stone toward atomically thin and flexible field-effect transistors. One key parameter to describe the charge transport is the time between two successive scattering events - the transport scattering time. In a recent report, we have shown that it is possible to use density functional theory to obtain the band structure of two-dimensional semiconductors in presence of field effect doping. Here, we report a simple method to extract the scattering time from the experimental conductivity and from the knowledge of the band structure. We apply our approach to monolayers and multilayers of MoS$_2$, MoSe$_2$, MoTe$_2$, WS$_2$, and WSe$_2$ in presence of a gate. In WS$_2$, for which accurate measurements of mobility have been published, we find that the scattering time is inversely proportional to the density of states at the Fermi level. Finally, we show that it is possible to identify the critical doping at which different valleys start to be occupied from the doping-dependence of the conductivity.

preprint2016arXiv

Spin susceptibility and electron-phonon coupling of two-dimensional materials by range-separated hybrid density functionals: Case study of Li$_x$ZrNCl

We investigate the capability of density functional theory (DFT) to appropriately describe the spin susceptibility, $χ_s$, and the intervalley electron-phonon coupling in Li$_x$ZrNCl. At low doping, Li$_x$ZrNCl behaves as a two-dimensional two-valley electron gas, with parabolic bands. In such a system, $χ_s$ increases with decreasing doping because of the electron-electron interaction. We show that DFT with local functionals (LDA/GGA) is not capable of reproducing this behavior. The use of exact exchange in Hartree-Fock (HF) or in DFT hybrid functionals enhances $χ_s$. HF, B3LYP, and PBE0 approaches overestimate $χ_s$, whereas the range-separated HSE06 functional leads to results similar to those obtained in the random phase approximation (RPA) applied to a two-valley two-spin electron gas. Within HF, Li$_x$ZrNCl is even unstable towards a ferromagnetic state for $x<0.16$. The intervalley phonons induce an imbalance in the valley occupation that can be viewed as the effect of a pseudomagnetic field. Thus, similarly to what happens for $χ_s$, the electron-phonon coupling of intervalley phonons is enhanced by the electron-electron interaction. Only hybrid DFT functionals capture such an enhancement and the HSE06 functional reproduces the RPA results presented in M. Calandra {\it et al.} [Phys. Rev. Lett. {\bf 114}, 077001 (2015)]. These results imply that the description of the susceptibility and electron-phonon coupling with a range-separated hybrid functional would be important also in other two-dimensional weakly doped semiconductors, such as transition-metal dichalcogenides and graphene.

preprint2016arXiv

Two-dimensional Fröhlich interaction in transition-metal dichalcogenide monolayers: Theoretical modeling and first-principles calculations

We perform ab initio calculations of the coupling between electrons and small-momentum polar-optical phonons in monolayer transition metal dichalcogenides of the 2H type: MoS2, MoSe2, MoTe2, WS2, and WSe2. The so-called Fröhlich interaction is fundamentally affected by the dimensionality of the system. In a plane-wave framework with periodic boundary conditions, this coupling is affected by the spurious interaction between the 2D material and its periodic images. To overcome this, we perform density functional perturbation theory calculations with a truncated Coulomb interaction in the out-of-plane direction. We show that the 2D Fröhlich interaction is much stronger than assumed in previous ab initio studies. We provide analytical models depending on the effective charges and dielectric properties of the materials to interpret our ab initio calculations. Screening is shown to play a fundamental role in the phonon-momentum dependency of the polar-optical coupling, with a crossover between two regimes depending on the dielectric properties of the material relative to its environment. The Fröhlich interaction is screened by the dielectric environment in the limit of small phonon momenta and sharply decreases due to stronger screening by the monolayer at finite momenta. The small-momentum regime of the ab initio Fröhlich interaction is reproduced by a simple analytical model, for which we provide the necessary parameters. At larger momenta, however, direct ab initio calculations of electron-phonon interactions are necessary to capture band-specific effects. We compute and compare the carrier relaxation times associated to the scattering by both LO and A1 phonon modes. While both modes are capable of relaxing carriers on timescales under the picosecond at room temperature, their absolute and relative importance vary strongly depending on the material, the band, and the substrate.

preprint2015arXiv

Density-functional calculation of static screening in 2D materials: the long-wavelength dielectric function of graphene

We calculate the long-wavelength static screening properties of both neutral and doped graphene in the framework of density-functional theory. We use a plane-wave approach with periodic images in the third dimension and truncate the Coulomb interactions to eliminate spurious interlayer screening. We carefully address the issue of extracting two dimensional dielectric properties from simulated three-dimensional potentials. We compare this method with analytical expressions derived for two dimensional massless Dirac fermions in the random phase approximation. We evaluate the contributions of the deviation from conical bands, exchange-correlation and local-fields. For momenta smaller than twice the Fermi wavevector, the static screening of graphene within the density-functional perturbative approach agrees with the results for conical bands within random phase approximation and neglecting local fields. For larger momenta, we find that the analytical model underestimates the static dielectric function by $\approx 10%$, mainly due to the conical band approximation.

preprint2015arXiv

Dissociation products and structures of solid H2S at strong compression

Hydrogen sulfides have recently received a great deal of interest due to the record high superconducting temperatures of up to 203 K observed on strong compression of dihydrogen sulfide (H2S). A joint theoretical and experimental study is presented in which decomposition products and structures of compressed H2S are characterized, and their superconducting properties are calculated. In addition to the experimentally known H2S and H3S phases, our first-principles structure searches have identified several energetically competitive stoichiometries that have not been reported previously; H2S3, H3S2, and H4S3. In particular, H4S3 is predicted to be thermodynamically stable within a large pressure range of 25-113 GPa. High-pressure room-temperature X-ray diffraction measurements confirm the presence of H3S and H4S3 through decomposition of H2S that emerge at 27 GPa and coexist with residual H2S, at least up to the highest pressure studied in our experiments of 140 GPa. Electron-phonon coupling calculations show that H4S3 has a small Tc of below 2 K, and that H2S is mainly responsible for the observed superconductivity of samples prepared at low temperature (<100K).

preprint2015arXiv

Electronic and vibrational properties of $\textup{TiSe}_2$ in the charge-density wave phase from first principles

We study the charge-density wave phase in $\textup{TiSe}_2$ by using first principle calculations. We show that, regardless of the local functional used and as long as the cell parameters are in agreement with the experiment, density-functional calculations are able to reproduce not only the structural instability of $\textup{TiSe}_2$, but also the effective distortion observed in the experiments. We study the electronic structure evolution of the system under the charge-density wave deformation. In particular, we show that the energy bands for the distorted superstructure, unfolded into the original Brillouin zone, are in reasonable agreement with angle-resolved photoemission spectroscopy (ARPES) data taken at low temperature. On the contrary, the energy bands for the undistorted structure are not in good agreement with ARPES at high temperature. Motivated by these results, we investigate the effect of the correlation on the electrons of the localized Ti-$d$ orbitals by using the LDA+$U$ method. We show that within this approximation the electronic bands for both the undistorted and distorted structure are in very good agreement with ARPES. On the other hand, the $U$ eliminates the phonon instability of the system. Some possible explanations for this counter intuitive result are proposed. Particularly, the possibility of taking into account the dependence of the parameter $U$ from the atomic positions is suggested.

preprint2015arXiv

First-principles theory of field-effect doping in transition-metal dichalcogenides: Structural properties, electronic structure, Hall coefficient, and electrical conductivity

We investigate how field-effect doping affects the structural properties, the electronic structure and the Hall coefficient of few-layers transition metal dichalcogenides by using density-functional theory. We consider mono-, bi-, and trilayers of the H polytype of MoS2, MoSe2, MoTe2, WS2, and WSe2 and provide a full database of electronic structures and Hall coefficients for hole and electron doping. We find that, for both electron and hole doping, the electronic structure depends on the number of layers and cannot be described by a rigid band shift. Furthermore, it is important to relax the structure under the asymmetric electric field. Interestingly, while the width of the conducting channel depends on the doping, the number of occupied bands at each given k point is almost uncorrelated with the thickness of the doping-charge distribution. Finally, we calculate within the constant-relaxation-time approximation the electrical conductivity and the inverse Hall coefficient. We demonstrate that in some cases the charge determined by Hall-effect measurements can deviate from the real charge by up to 50%. For hole-doped MoTe2 the Hall charge has even the wrong polarity at low temperature. We provide the mapping between the doping charge and the Hall coefficient. In the appendix we present more than 250 band structures for all doping levels of the transition-metal dichalcogenides considered within this work.

preprint2015arXiv

Hydrogen sulphide at high pressure: a strongly-anharmonic phonon-mediated superconductor

We use first principles calculations to study structural, vibrational and superconducting properties of H$_2$S at pressures $P\ge 200$ GPa. The inclusion of zero point energy leads to two different possible dissociations of H$_2$S, namely 3H$_2$S $\to$ 2H$_3$S + S and 5H$_2$S $\to$ 3H$_3$S + HS$_2$, where both H$_3$S and HS$_2$ are metallic. For H$_3$S, we perform non-perturbative calculations of anharmonic effects within the self-consistent harmonic approximation and show that the harmonic approximation strongly overestimates the electron-phonon interaction ($λ\approx 2.64$ at 200 GPa) and T$_c$. Anharmonicity hardens HS bond-stretching modes and softens H--S bond-bending modes. As a result, the electron-phonon coupling is suppressed by $30\%$ ($λ\approx 1.84$ at 200 GPa). Moreover, while at the harmonic level T$_c$ decreases with increasing pressure, the inclusion of anharmonicity leads to a T$_c$ that is almost independent of pressure. High pressure hydrogen sulfide is a strongly anharmonic superconductor.

preprint2015arXiv

Phonon effects on x-ray absorption and nuclear magnetic resonance spectroscopies

In material sciences, spectroscopic approaches combining ab initio calculations with experiments are commonly used to accurately analyze the experimental spectral data. Most state-of-the-art first-principle calculations are usually performed assuming an equilibrium static lattice. Yet, nuclear motion affects spectra even when reduced to the zero-point motion at 0 K. We propose a framework based on Density-Functional Theory that includes quantum thermal fluctuations in theoretical X- ray Absorption Near-Edge Structure (XANES) and solid-state Nuclear Magnetic Resonance (NMR) spectroscopies and allows to well describe temperature effects observed experimentally. Within the Born-Oppenheimer and quasi-harmonic approximations, we incorporate the nuclear motion by generating several non-equilibrium configurations from the dynamical matrix. The averaged calculated XANES and NMR spectral data have been compared to experiments in MgO, proof-of-principle compound. The good agreement obtained between experiments and calculations validates the developed approach, which suggests that calculating the XANES spectra at finite temperature by averaging individual non-equilibrium configurations is a suitable approximation. This study high- lights the relevance of phonon renormalization and the relative contributions of thermal expansion and nuclear dynamics on NMR and XANES spectra on a wide range of temperatures.

preprint2015arXiv

Quantum Hydrogen-Bond Symmetrization and High-Temperature Superconductivity in Hydrogen Sulfide

Hydrogen compounds are peculiar as the quantum nature of the proton can crucially affect their structural and physical properties. A remarkable example are the high-pressure phases of H$_2$O, where quantum proton fluctuations favor the symmetrization of the H bond and lower by 30 GPa the boundary between the asymmetric structure and the symmetric one. Here we show that an analogous quantum symmetrization occurs in the recently discovered sulfur hydride superconductor with the record superconducting critical temperature $T_c=203$ K at 155 GPa. In this system, according to classical theory, superconductivity occurs via formation of a structure of stoichiometry H$_3$S with S atoms arranged on a body-centered-cubic (bcc) lattice. For $P \gtrsim 175$ GPa, the H atoms are predicted to sit midway between two S atoms, in a structure with $Im\bar3m$ symmetry. At lower pressures the H atoms move to an off-center position forming a short H$-$S covalent bond and a longer H$\cdots$S hydrogen bond, in a structure with $R3m$ symmetry. X-ray diffraction experiments confirmed the H$_3$S stoichiometry and the S lattice sites, but were unable to discriminate between the two phases. Our present ab initio density-functional theory (DFT) calculations show that the quantum nuclear motion lowers the symmetrization pressure by 72 GPa. Consequently, we predict that the $Im\bar3m$ phase is stable over the whole pressure range within which a high $T_c$ was measured. The observed pressure-dependence of $T_c$ is closely reproduced in our calculations for the $Im\bar3m$ phase, but not for the $R3m$ phase. Thus, the quantum nature of the proton completely rules the superconducting phase diagram of H$_3$S.

preprint2015arXiv

Strong anharmonicity induces quantum melting of charge density wave in 2H-NbSe$\_2$ under pressure

The interplay between charge density wave (CDW) order and superconductivity has attracted much attention. This is the central issue of along standing debate in simple transition metal dichalcogenides without strong electronic correlations, such as 2H-NbSe$\_2$, in which twosuch phases coexist. The importance of anisotropic electron-phonon interaction has been recently highlighted from both theoretical and experimental point of view, and explains some of the key features of the formation of the CDW in this system. On the other hand, other aspects, such as the effects of anharmonicity, remain poorly understood despite their manifest importance in such soft-phonon driven phase transition. At the theoretical level in particular, their prohibitive computational price usually prevents their investigation within conventional perturbative approaches.Here, we address this issue using a combination of high resolution inelastic X-ray scattering measurements of the phonon dispersion, as afunction of temperature and pressure, with state of the art ab initio calculations. By explicitly taking into account anharmonic effects, we obtain an accurate, quantitative, description of the (P,T) dependence of the phonon spectrum, accounting for the rapid destruction of the CDW under pressure by zero mode vibrations - or quantum fluctuations - of the lattice. The low-energy longitudinal acoustic mode that drives the CDW transition barely contributes to superconductivity, explaining the insensitivity of the superconducting critical temperature to the CDW transition.

preprint2015arXiv

Universal enhancement of superconductivity in two dimensional semiconductors at low doping by electron-electron interaction

In two-dimensional multivalley semiconductors, at low doping, even a moderate electron-electron interaction enhances the response to any perturbation inducing a valley polarization. If the valley polarization is due to the electron-phonon coupling, the electron-electron interaction results in an enhancement of the superconducting critical temperature. By performing first principles calculations beyond density functional theory, we prove that this effect accounts for the unconventional doping-dependence of the superconducting transition-temperature (T$_c$) and of the magnetic susceptibility measured in Li$_x$ZrNCl. By finding the conditions for a maximal T$_c$ enhancement, we show how weakly-doped two-dimensional semiconductors provide a route towards high T$_c$ superconductivity.

preprint2014arXiv

Downfolding electron-phonon Hamiltonians from ab-initio calculations: application to K$_3$Picene

We propose an electron-phonon parameterization which reliably reproduces the geometry and harmonic frequencies of a real system. With respect to standard electron-phonon models, it adds a "double-counting" correction, which takes into account the lattice deformation as the system is dressed by low-energy electron-phonon processes. We show the importance of this correction by studying potassium-doped picene (K$_3$Picene), recently claimed to be a superconductor with a $T_c$ of up to 18 K. The Hamiltonian parameters are derived from ab-initio density functional theory, and the lattice model is solved by dynamical mean-field theory. Our calculations include the effects of electron-electron interactions and local electron-phonon couplings. Even with the inclusion of a strongly coupled molecular phonon, the Hubbard repulsion prevails and the system is an insulator with a small Mott gap of $\approx$ 0.2 eV.

preprint2014arXiv

Electrochemical doping of few layer ZrNCl from first-principles: electronic and structural properties in field-effect configuration

We develop a first-principles theoretical approach to doping in field-effect devices. The method allows for calculation of the electronic structure as well as complete structural relaxation in field-effect configuration using density-functional theory. We apply our approach to ionic-liquid-based field-effect doping of monolayer, bilayer, and trilayer ZrNCl and analyze in detail the structural changes induced by the electric field. We show that, contrary to what is assumed in previous experimental works, only one ZrNCl layer is electrochemically doped and that this induces large structural changes within the layer. Surprisingly, despite these structural and electronic changes, the density of states at the Fermi energy is independent of the doping. Our findings imply a substantial revision of the phase diagram of electrochemically doped ZrNCl and elucidate crucial differences with superconductivity in Li intercalated bulk ZrNCl.

preprint2014arXiv

Electron-Phonon Interactions and the Intrinsic Electrical Resistivity of Graphene

We present a first-principles study of the temperature- and density-dependent intrinsic electrical resistivity of graphene. We use density-functional theory and density-functional perturbation theory together with very accurate Wannier interpolations to compute all electronic and vibrational properties and electron-phonon coupling matrix elements; the phonon-limited resistivity is then calculated within a Boltzmann-transport approach. An effective tight-binding model, validated against first-principles results, is also used to study the role of electron-electron interactions at the level of many-body perturbation theory. The results found are in excellent agreement with recent experimental data on graphene samples at high carrier densities and elucidate the role of the different phonon modes in limiting electron mobility. Moreover, we find that the resistivity arising from scattering with transverse acoustic phonons is 2.5 times higher than that from longitudinal acoustic phonons. Last, high-energy, optical, and zone-boundary phonons contribute as much as acoustic phonons to the intrinsic electrical resistivity even at room temperature and become dominant at higher temperatures.

preprint2014arXiv

First-principles calculations of phonon frequencies, lifetimes and spectral functions from weak to strong anharmonicity: the example of palladium hydrides

The variational stochastic self-consistent harmonic approximation is combined with the calculation of third-order anharmonic coefficients within density-functional perturbation theory and the "$2n+1$" theorem to calculate anharmonic properties of crystals. It is demonstrated that in the perturbative limit the combination of these two methods yields the perturbative phonon linewidth and frequency shift in a very efficient way, avoiding the explicit calculation of fourth-order anharmonic coefficients. Moreover, it also allows calculating phonon lifetimes and inelastic neutron scattering spectra in solids where the harmonic approximation breaks down and a non-perturbative approach is required to deal with anharmonicity. To validate our approach, we calculate the anharmonic phonon linewidth in the strongly anharmonic palladium hydrides. We show that due to the large anharmonicity of hydrogen optical modes the inelastic neutron scattering spectra are not characterized by a Lorentzian line-shape, but by a complex structure including satellite peaks.

preprint2014arXiv

Phonon-limited resistivity of graphene by first-principle calculations: electron-phonon interactions, strain-induced gauge field and Boltzmann equation

Electron-phonon coupling in graphene is extensively modeled and simulated from first principles. We find that using an accurate model for the polarizations of the acoustic phonon modes is crucial to obtain correct numerical results. The interactions between electrons and acoustic phonon modes, the gauge field and deformation potential, are calculated at the DFT level in the framework of linear response. The zero-momentum limit of acoustic phonons is interpreted as a strain pattern, allowing the calculation of the acoustic gauge field parameter in the GW approximation. The role of electronic screening on the electron-phonon matrix elements is investigated. We then solve the Boltzmann equation semi-analytically in graphene, including both acoustic and optical phonon scattering. We show that, in the Bloch-Grüneisen and equipartition regimes, the electronic transport is mainly ruled by the unscreened acoustic gauge field, while the contribution due to the deformation potential is negligible and strongly screened. By comparing with experimental data, we show that the contribution of acoustic phonons to resistivity is doping- and substrate-independent. The DFT+GW approach underestimates this contribution to resistivity by about 30 %. Above 270K, the calculated resistivity underestimates the experimental one more severely, the underestimation being larger at lower doping. We show that, beside remote phonon scattering, a possible explanation for this disagreement is the electron-electron interaction that strongly renormalizes the coupling to intrinsic optical-phonon modes. Finally, after discussing the validity of the Matthiessen rule in graphene, we derive simplified analytical solutions of the Boltzmann equation to extract the coupling to acoustic phonons, related to the strain-induced gauge field, directly from experimental data.

preprint2014arXiv

Two-dimensional analysis of the double-resonant 2D Raman mode in bilayer graphene

By computing the double-resonant Raman scattering cross-section completely from first principles and including electron-electron interaction at the $GW$ level, we unravel the dominant contributions for the double-resonant 2D-mode in bilayer graphene. We show that, in contrast to previous works, the so-called inner processes are dominant and that the 2D-mode lineshape is described by three dominant resonances around the $K$ point. We show that the splitting of the TO phonon branch in $Γ-K$ direction, as large as 12 cm$^{-1}$ in $GW$ approximation, is of great importance for a thorough description of the 2D-mode lineshape. Finally, we present a method to extract the TO phonon splitting and the splitting of the electronic bands from experimental data.

preprint2013arXiv

Ab initio variational approach for evaluating lattice thermal conductivity

We present a first-principles theoretical approach for evaluating the lattice thermal conductivity based on the exact solution of the Boltzmann transport equation. We use the variational principle and the conjugate gradient scheme, which provide us with an algorithm faster than the one previously used in literature and able to always converge to the exact solution. Three-phonon normal and umklapp collision, isotope scattering and border effects are rigorously treated in the calculation. Good agreement with experimental data for diamond is found. Moreover we show that by growing more enriched diamond samples it is possible to achieve values of thermal conductivity up to three times larger than the commonly observed in isotopically enriched diamond samples with 99.93% C12 and 0.07 C13.

preprint2013arXiv

Anharmonic free energies and phonon dispersions from the stochastic self-consistent harmonic approximation: application to platinum and palladium hydrides

Harmonic calculations based on density-functional theory are generally the method of choice for the description of phonon spectra of metals and insulators. The inclusion of anharmonic effects is, however, delicate as it relies on perturbation theory requiring a considerable amount of computer time, fast increasing with the cell size. Furthermore, perturbation theory breaks down when the harmonic solution is dynamically unstable or the anharmonic correction of the phonon energies is larger than the harmonic frequencies themselves.We present a stochastic implementation of the self-consistent harmonic approximation valid to treat anharmonicity at any temperature in the non-perturbative regime. The method is based on the minimization of the free energy with respect to a trial density matrix described by an arbitrary harmonic Hamiltonian. The minimization is performed with respect to all the free parameters in the trial harmonic Hamiltonian, namely, equilibrium positions, phonon frequencies and polarization vectors. The gradient of the free energy is calculated following a stochastic procedure. The method can be used to calculate thermodynamic properties, dynamical properties and anharmonic corrections to the Eliashberg function of the electron-phonon coupling. The scaling with the system size is greatly improved with respect to perturbation theory. The validity of the method is demonstrated in the strongly anharmonic palladium and platinum hydrides. In both cases we predict a strong anharmonic correction to the harmonic phonon spectra, far beyond the perturbative limit. In palladium hydrides we calculate thermodynamic properties beyond the quasiharmonic approximation, while in PtH we demonstrate that the high superconducting critical temperatures at 100 GPa predicted in previous calculations based on the harmonic approximation are strongly suppressed when anharmonic effects are included.

preprint2013arXiv

Anharmonic properties from a generalized third order ab~initio approach: theory and applications to graphite and graphene

We have implemented a generic method, based on the 2n+1 theorem within density functional perturbation theory, to calculate the anharmonic scattering coefficients among three phonons with arbitrary wavevectors. The method is used to study the phonon broadening in graphite and graphene mono- and bi-layer. The broadening of the high-energy optical branches is highly nonuniform and presents a series of sudden steps and spikes. At finite temperature, the two linearly dispersive acoustic branches TA and LA of graphene have nonzero broadening for small wavevectors. The broadening in graphite and bi-layer graphene is, overall, very similar to the graphene one, the most remarkable feature being the broadening of the quasi acoustical ZO' branch. Finally, we study the intrinsic anharmonic contribution to the thermal conductivity of the three systems, within the single mode relaxation time approximation. We find the conductance to be in good agreement with experimental data for the out-of-plane direction but to underestimate it by a factor 2 in-plane.

preprint2013arXiv

Extended Czjzek model applied to NMR parameter distributions in sodium metaphosphate glass

The Extended Czjzek Model (ECM) is applied to the distribution of NMR parameters of a simple glass model (sodium metaphosphate, $\mathrm{NaPO_3}$) obtained by Molecular Dynamics (MD) simulations. Accurate NMR tensors, Electric Field Gradient (EFG) and Chemical Shift Anisotropy (CSA), are calculated from Density Functional Theory (DFT) within the well-established PAW/GIPAW framework. Theoretical results are compared to experimental high-resolution solid-state NMR data and are used to validate the considered structural model. The distributions of the calculated coupling constant $C_Q\propto |V_{zz}|$ and of the asymmetry parameter $η_Q$ that characterize the quadrupolar interaction are discussed in terms of structural considerations with the help of a simple point charge model. Finally, the ECM analysis is shown to be relevant for studying the distribution of CSA tensor parameters and gives new insight into the structural characterization of disordered systems by solid-state NMR.

preprint2013arXiv

First-Principles Theory of Anharmonicity and the Inverse Isotope Effect in Superconducting Palladium-Hydride Compounds

Palladium hydrides display the largest isotope effect anomaly known in literature. Replacement of hydrogen with the heavier isotopes leads to higher superconducting temperatures, a behavior inconsistent with harmonic theory. Solving the self-consistent harmonic approximation by a stochastic approach, we obtain the anharmonic free energy, the thermal expansion and the superconducting properties fully ab initio. We find that the phonon spectra are strongly renormalized by anharmonicity far beyond the perturbative regime. Superconductivity is phonon mediated, but the harmonic approximation largely overestimates the superconducting critical temperatures. We explain the inverse isotope effect, obtaining a -0.38 value for the isotope coefficient in good agreement with experiments, hydrogen anharmonicity being the main responsible for the isotope anomaly.

preprint2013arXiv

Quantum Monte Carlo study of the protonated water dimer

We report an extensive theoretical study of the protonated water dimer (Zundel ion) by means of the highly correlated variational Monte Carlo and lattice regularized Monte Carlo approaches. This system represents the simplest model for proton transfer (PT) and a correct description of its properties is essential in order to understand the PT mechanism in more complex acqueous systems. Our Jastrow correlated AGP wave function ensures an accurate treatment of electron correlations. Exploiting the advantages of contracting the primitive basis set over atomic hybrid orbitals, we are able to limit dramatically the number of variational parameters with a systematic control on the numerical precision, crucial in order to simulate larger systems. We investigate energetics and geometrical properties of the Zundel ion as a function of the oxygen-oxygen distance, taken as reaction coordinate. In both cases, our QMC results are found in excellent agreement with coupled cluster CCSD(T) technique, the quantum chemistry "gold standard". Calculations on proton transfer static barriers and dissociation energies display the same agreement. A comparison with density functional theory (DFT) results in the PBE approximation points out the crucial role of electron correlations for a correct description of the PT in the dimer. Our method is able to resolve the tiny energy differences (~ 0.1 Kcal/mol) and the corresponding structural variations at which the proton transfer takes place; it combines these features with a N^3-N^4 scaling with number of particles (favorable with respect to other post-DFT methods), as proven by a benchmark calculation on a larger protonated water cluster. The QMC approach used here represents a promising candidate to provide the first high-level ab initio description of PT in water.

preprint2012arXiv

Hidden polymorphs drive the vitrification in B2O3

Understanding the conditions which favor crystallisation or vitrification of liquids has been a long-standing scientific problem. Another connected, and not yet well understood question is the relationship between the glassy and the various possible crystalline forms a system may adopt. In this context, B2O3 is a puzzling case of study since i) it is one of the best glass-forming systems despite an apparent lack of low-pressure polymorphism ii) it vitrifies in a glassy form abnormally different from the only known crystalline phase at ambient pressure iii) it never crystallises from the melt unless pressure is applied, an intriguing behaviour known as the crystallisation anomaly. Here, by means of ab-initio calculations, we discover the existence of novel B2O3 crystalline polymorphs with structural properties similar to the glass and formation energies comparable to the known ambient crystal. The resulting configurational degeneracy drives the system vitrification at ambient pressure. The degeneracy is lifted under pressure, unveiling the origin of the crystallisation anomaly. This work reconciles the behaviour of B2O3 with that from other glassy systems and reaffirms the role played by polymorphism in a system's ability to vitrify. Some of the predicted crystals are cage-like materials entirely made of three-fold rings, opening new perspectives for the synthesis of boron-based nanoporous materials.

preprint2012arXiv

Local and non-local electron-phonon couplings in K3Picene and the effect of metallic screening

We analyze the properties of electron-phonon couplings in K3Picene by exploiting a molecular orbital representation derived in the maximally localized Wannier function formalism. This allows us to go beyond the analysis done in Phys. Rev. Lett. 107, 137006 (2011), and separate not only the intra- and intermolecular phonon contributions but also the local and non-local electronic states in the electron-phonon matrix elements. Despite the molecular nature of the crystal, we find that the purely molecular contributions (Holstein-like couplings where the local deformation potential is coupled to intramolecular phonons) account for only 20% of the total electron-phonon interaction lambda. In particular, the Holstein-like contributions to lambda in K3Picene are four times smaller than those computed for an isolated neutral molecule, as they are strongly screened by the metallic bands of the doped crystal. Our findings invalidate the use of molecular electron-phonon calculations to estimate the total electron-phonon coupling in metallic picene, and possibly in other doped metallic molecular crystals. The major contribution (80%) to lambda in K3Picene comes from non-local couplings due to phonon modulated hoppings. We show that the crystal geometry together with the molecular picene structure leads to a strong 1D spatial anisotropy of the non-local couplings. Finally, based on the parameters derived from our density functional theory calculations, we propose a lattice modelization of the electron-phonon couplings in K3Picene which gives 90% of ab-initio lambda.

preprint2012arXiv

Zener tunneling in the electrical transport of quasi-metallic carbon nanotubes

We study theoretically the impact of Zener tunneling on the charge-transport properties of quasi-metallic (Qm) carbon nanotubes (characterized by forbidden band gaps of few tens of meV). We also analyze the interplay between Zener tunneling and elastic scattering on defects. To this purpose we use a model based on the master equation for the density matrix, that takes into account the inter-band Zener transitions induced by the electric field (a quantum mechanical effect), the electron-defect scattering and the electron-phonon scattering. In presence of Zener tunnelling the Qm tubes support an electrical current even when the Fermi energy lies in the forbidden band gap. In absence of elastic scattering (in high quality samples), the small size of the band gap of Qm tubes enables Zener tunnelling for realistic values of the the electric field (above $\sim$ 1 V/μm). The presence of a strong elastic scattering (in low quality samples) further decreases the values of the field required to observe Zener tunnelling. Indeed, for elastic-scattering lengths of the order of 50 nm, Zener tunnelling affects the current/voltage characteristic already in the linear regime. In other words, in quasi-metallic tubes, Zener tunneling is made more visible by defects.

preprint2011arXiv

Charge density wave and superconducting dome in TiSe2 from electron-phonon interaction

At low temperature TiSe2 undergoes a charge density wave instability. Superconductivity is stabilized either by pressure or by Cu intercalation. We show that the pressure phase diagram of TiSe2 is well described by first-principles calculations. At pressures smaller than 4 GPa charge density wave ordering occurs, in agreement with experiments. At larger pressures the disappearing of the charge density wave is due to a stiffening of the short-range force-constants and not to the variation of nesting with pressure. Finally we show that the behavior of Tc as a function of pressure is entirely determined by the electron-phonon interaction without need of invoking excitonic mechanisms. Our work demonstrates that phase-diagrams with competing orders and a superconducting dome are also obtained in the framework of the electron-phonon interaction.

preprint2011arXiv

How to make graphene superconducting

Graphene is the physical realization of many fundamental concepts and phenomena in solid state-physics, but in the long list of graphene remarkable properties, a fundamental block is missing: superconductivity. Making graphene superconducting is relevant as the easy manipulation of this material by nanolytographic techniques paves the way to nanosquids, one-electron superconductor-quantum dot devices, superconducting transistors at the nano-scale and cryogenic solid-state coolers. Here we explore the doping of graphene by adatoms coverage. We show that the occurrence of superconductivity depends on the adatom in analogy with graphite intercalated compounds (GICs). However, most surprisingly, and contrary to the GIC case, Li covered graphene is superconducting at much higher temperature with respect to Ca covered graphene.

preprint2011arXiv

Intercalant and intermolecular phonon assisted superconductivity in K-doped picene

K$_3$-picene is a superconducting molecular crystal with critical temperature T$_c=7$K or 18K, depending on preparation conditions. Using density functional theory we show that electron-phonon interaction accounts for T$_c ~ 3-8$ K. The average electron-phonon coupling, calculated by including the phonon energy scale in the electron-phonon scattering, is $λ=0.73$ and $ω_{\rm log}=18.0\,$ meV. Intercalant and intermolecular phonon-modes contribute substantially (40%) to $λ$ as also shown by the isotope exponents of potassium (0.19) and carbon (0.31). The relevance of these modes makes superconductivity in K-doped picene peculiar and different from that of fullerenes.

preprint2011arXiv

Theory of double-resonant Raman spectra in graphene: intensity and line shape of defect-induced and two-phonon bands

We calculate the double resonant (DR) Raman spectrum of graphene, and determine the lines associated to both phonon-defect processes, and two-phonons ones. Phonon and electronic dispersions reproduce calculations based on density functional theory corrected with GW. Electron-light, -phonon, and -defect scattering matrix elements and the electronic linewidth are explicitly calculated. Defect-induced processes are simulated by considering different kind of idealized defects. For an excitation energy of $ε_L=2.4$ eV, the agreement with measurements is very good and calculations reproduce: the relative intensities among phonon-defect or among two-phonon lines; the measured small widths of the D, $D'$, 2D and $2D'$ lines; the line shapes; the presence of small intensity lines in the 1800, 2000 cm$^{-1}$ range. We determine how the spectra depend on the excitation energy, on the light polarization, on the electronic linewidth, on the kind of defects and on their concentration. According to the present findings, the intensity ratio between the $2D'$ and 2D lines can be used to determine experimentally the electronic linewidth. The intensity ratio between the $D$ and $D'$ lines depends on the kind of model defect, suggesting that this ratio could possibly be used to identify the kind of defects present in actual samples. Charged impurities outside the graphene plane provide an almost undetectable contribution to the Raman signal.

preprint2010arXiv

Clar's Theory, STM Images, and Geometry of Graphene Nanoribbons

We show that Clar's theory of the aromatic sextet is a simple and powerful tool to predict the stability, the π-electron distribution, the geometry, the electronic/magnetic structure of graphene nanoribbons with different hydrogen edge terminations. We use density functional theory to obtain the equilibrium atomic positions, simulated scanning tunneling microscopy (STM) images, edge energies, band gaps, and edge-induced strains of graphene ribbons that we analyze in terms of Clar formulas. Based on their Clar representation, we propose a classification scheme for graphene ribbons that groups configurations with similar bond length alternations, STM patterns, and Raman spectra. Our simulations show how STM images and Raman spectra can be used to identify the type of edge termination.

preprint2010arXiv

Current-voltage characteristics of graphene devices: interplay between Zener-Klein tunneling and defects

We report a theoretical/experimental study of current-voltage characteristics (I-V) of graphene devices near the Dirac point. The I-V can be described by a power law (I \propto V^α, with 1< α<= 1.5). The exponent is higher when the mobility is lower. This superlinear I-V is interpreted in terms of the interplay between Zener-Klein transport, that is tunneling between different energy bands, and defect scattering. Surprisingly, the Zener-Klein tunneling is made visible by the presence of defects.

preprint2010arXiv

Doped graphene as tunable electron-phonon coupling material

We present a new way to tune the electron-phonon coupling (EPC) in graphene by changing the deformation potential with electron/hole doping. We show the EPC for highest optical branch at the high symmetry point K, acquires a strong dependency on the doping level due to electron-electron correlation not accounted in mean-field approaches. Such a dependency influences the dispersion (with respect to the laser energy) of the Raman D and 2D lines and the splitting of the 2D peak in multi-layer graphene. Finally this doping dependence opens the possibility to construct tunable electronic devices through the external control of the EPC.

preprint2010arXiv

Effects of magnetism and doping on the electron-phonon coupling in BaFe$_{2}$As$_{2}$

We calculate the effect of local magnetic moments on the electron-phonon coupling in BaFe$_{2}$As$_{2}+δ$ using the density functional perturbation theory. We show that the magnetism enhances the total electron-phonon coupling by $\sim 50%$, up to $λ\lesssim 0.35$, still not enough to explain the high critical temperature, but strong enough to have a non-negligible effect on superconductivity, for instance, by frustrating the coupling with spin fluctuations and inducing order parameter nodes. The enhancement comes mostly from a renormalization of the electron-phonon matrix elements. We also investigate, in the rigid band approximation, the effect of doping, and find that $λ$ versus doping does not mirror the behavior of the density of states; while the latter decreases upon electron doping, the former does not, and even increases slightly.

preprint2010arXiv

First-principles theory of the orbital magnetization

We compute the orbital magnetization in real materials by evaluating a recently discovered formula for periodic systems, within density functional theory. We obtain improved values of the orbital magnetization in the ferromagnetic metals Fe, Co, and Ni, by taking into account the contribution of the interstitial regions neglected so far in literature. We also use the orbital magnetization to compute the EPR $g$-tensor in molecules and solids. The present approach reproduces the $g$-tensor obtained by linear response (LR), when the spin-orbit can be treated as a perturbation. However, it can also be applied to radicals and defects with an orbital-degenerate ground-state or containing heavy atoms, that can not be properly described by LR.

preprint2010arXiv

Phonon dispersion and low energy anomaly in CaC$_6$

We report measurements of phonon dispersion in CaC$_6$ using inelastic X-ray and neutron scattering. We find good overall agreement, particularly in the 50 meV energy region, between experimental data and first-principles density-functional-theory calculations. However, on the longitudinal dispersion along the $(1 1 1)$ axis of the rhombohedral representation, we find an unexpected anti-crossing with an additional longitudinal mode, at about 11 meV. At a comparable energy, we observe also unexpected intensity on the in-plane direction. These results resolve the previous incorrect assignment of a longitudinal phonon mode to a transverse mode in the same energy range. By calculating the electron susceptibility from first principles we show that this longitudinal excitation is unlikely to be due to a plasmon and consequently can probably be due to defects or vacancies present in the sample.

preprint2010arXiv

Structure and stability of graphene nanoribbons in oxygen, carbon dioxide, water, and ammonia

We determine, by means of density functional theory, the stability and the structure of graphene nanoribbon (GNR) edges in presence of molecules such as oxygen, water, ammonia, and carbon dioxide. As in the case of hydrogen-terminated nanoribbons, we find that the most stable armchair and zigzag configurations are characterized by a non-metallic/non-magnetic nature, and are compatible with Clar's sextet rules, well known in organic chemistry. In particular, we predict that, at thermodynamic equilibrium, neutral GNRs in oxygen-rich atmosphere should preferentially be along the armchair direction, while water-saturated GNRs should present zigzag edges. Our results promise to be particularly useful to GNRs synthesis, since the most recent and advanced experimental routes are most effective in water and/or ammonia-containing solutions.

preprint2009arXiv

Ab-initio study of gap opening and screening effects in gated bilayer graphene

The electronic properties of doped bilayer graphene in presence of bottom and top gates have been studied and characterized by means of Density Functional Theory calculations. Varying independently the bottom and top gates it is possible to control separately the total doping charge on the sample, and the average external electric field acting on the bilayer. We show that, at fixed doping level, the band gap at the K point in the Brillouin zone depends linearly on the average electric field, whereas the corresponding proportionality coefficient has a non-monotonic dependence on doping. We find that the DFT-calculated band gap at K, for small doping levels, is roughly half of the band gap obtained with standard Tight Binding approach. We show that this discrepancy arises from an underestimate, in the TB model, of the screening of the system to the external electric field. In particular, on the basis of our DFT results we observe that, when the bilayer graphene is in presence of an external electric field, both an interlayer and an intralayer screening occur. Only the interlayer screening is included in TB calculations, while both screenings are fundamental for the description of the band gap opening. We finally provide a general scheme to obtain the full band structure of the gated bilayer graphene, for an arbitrary value of the external electric field and of the doping.

preprint2009arXiv

Effect of dimensionality on the charge-density-wave in few-layers 2H-NbSe$_2$

We investigate the charge density wave (CDW) instability in single and double layers, as well as in the bulk 2H-NbSe$_{2}$. We demonstrate that the density functional theory correctly describes the metallic CDW state in the bulk 2H-NbSe$_{2}$. We predict that both mono- and bilayer NbSe$_{2}$ undergo a CDW instability. However, while in the bulk the instability occurs at a momentum $\mathbf{q}_{CDW}\approx{2/3}\mathbf{ΓM}$, in free-standing layers it occurs at $\mathbf{q}_{CDW}\approx{1/2}\mathbf{ΓM}$. Furthermore, while in the bulk the CDW leads to a metallic state, in a monolayer the ground state becomes semimetallic, in agreement with recent experimental data. We elucidate the key role that an enhancement of the electron-phonon matrix element at $\mathbf{q}\approx\mathbf{q}_{CDW}$ plays in forming the CDW ground state.

preprint2009arXiv

Reducing the thermal conductivity of carbon nanotubes below the random isotope limit

We find that introducing segmented isotopic disorder patterns may considerably reduce the thermal conductivity of pristine carbon nanotubes below the uncorrelated disorder value. This is a result of the interplay between different length scales in the phonon scattering process. We use ab-initio atomistic Green's function calculations to quantify the effect of various types of segmentation similar to that experimentally produced by coalescence of isotope-engineered fullerenes.

preprint2009arXiv

Transport properties of graphene in the high-current limit

We present a detailed study of the high-current transport properties of graphene devices patterned in a four-point configuration. The current tends to saturate as the voltage across graphene is increased but never reaches the complete saturation as in metallic nanotubes. Measurements are compared to a model based on the Boltzmann equation, which includes electron scattering processes due to charged and neutral impurities, and graphene optical-phonons. The saturation is incomplete because of the competition between disorder and optical-phonon scattering.

preprint2006arXiv

Born-Oppenheimer Breakdown in Graphene

The Born-Oppenheimer approximation (BO) has proven effective for the accurate determination of chemical reactions, molecular dynamics and phonon frequencies in a wide range of metallic systems. Graphene, recently discovered in the free state, is a zero band-gap semiconductor, which becomes a metal if the Fermi energy is tuned applying a gate-voltage Vg. Graphene electrons near the Fermi energy have twodimensional massless dispersions, described by Dirac cones. Here we show that a change in Vg induces a stiffening of the Raman G peak (i.e. the zone-center E2g optical phonon), which cannot be described within BO. Indeed, the E2g vibrations cause rigid oscillations of the Dirac-cones in the reciprocal space. If the electrons followed adiabatically the Dirac-cone oscillations, no change in the phonon frequency would be observed. Instead, since the electron-momentum relaxation near the Fermi level is much slower than the phonon motion, the electrons do not follow the Dirac-cone displacements. This invalidates BO and results in the observed phonon stiffening. This spectacular failure of BO is quite significant since BO has been the fundamental paradigm to determine crystal vibrations from the early days of quantum mechanics.

preprint2005arXiv

Raman spectra of BN-nanotubes: Ab-initio and bond-polarizability model calculations

We present it ab-initio calculations of the non-resonant Raman spectra of zigzag and armchair BN nanotubes. In comparison, we implement a generalized bond-polarizability model where the parameters are extracted from first-principles calculations of the polarizability tensor of a BN sheet. For light-polarization along the tube-axis, the agreement between model and it ab-initio spectra is almost perfect. For perpendicular polarization, depolarization effects have to be included in the model in order to reproduce the it ab-initio Raman intensities.

preprint2005arXiv

Structural relaxations around Ti, Cr and Fe impurities in alpha-Al2O3 probed by x-ray absorption near edge structure combined with first-principles calculations

We determine the structural relaxations around paramagnetic impurities (Ti, Cr, Fe) in corundum (alpha-Al2O3), by combining x-ray absorption near edge structure (XANES) experiments and ab initio calculations. The structural relaxations are found to be very local. We then show that XANES is sensitive to small variations in interatomic distances within the coordination shell of the absorbing atom. The experiments were carried out on single crystals of ruby and sapphires. Linear dichroic signals are essential to characterize the geometry of the impurity site. The calculations were performed within a self-consistent ``non muffin-tin'' framework, that uses pseudopotentials, plane-wave basis set, and the continued fraction for the absorption cross section.

preprint2001arXiv

All-electron magnetic response with pseudopotentials: NMR chemical shifts

A theory for the ab initio calculation of all-electron NMR chemical shifts in insulators using pseudopotentials is presented. It is formulated for both finite and infinitely periodic systems and is based on an extension to the Projector Augmented Wave approach of Bloechl [P. E. Bloechl, Phys. Rev. B 50, 17953 (1994)] and the method of Mauri et al [F. Mauri, B.G. Pfrommer, and S.G. Louie, Phys. Rev. Lett. 77, 5300 (1996)]. The theory is successfully validated for molecules by comparison with a selection of quantum chemical results, and in periodic systems by comparison with plane-wave all-electron results for diamond.

preprint1995arXiv

First-principle study of excitonic self-trapping in diamond

We present a first-principles study of excitonic self-trapping in diamond. Our calculation provides evidence for self-trapping of the 1s core exciton and gives a coherent interpretation of recent experimental X-ray absorption and emission data. Self-trapping does not occur in the case of a single valence exciton. We predict, however, that self-trapping should occur in the case of a valence biexciton. This process is accompanied by a large local relaxation of the lattice which could be observed experimentally.

preprint1994arXiv

Large scale quantum simulations: C_60 impacts on a semiconducting surface

We present tight binding molecular dynamics simulations of C_60 collisions on the reconstructed diamond(111) surface, carried out with an O(N) method and with cells containing 1140 atoms. The results of our simulations are in very good agreement with experiments performed under the same impact conditions. Furthermore our calculations provide a detailed characterization of the microscopic processes occuring during the collision, and allow the identification of three impact regimes, as a function of the fullerene incident energy. Finally, the study of the reactivity between the cluster and the surface gives insight into the deposition mechanisms of C_60 on semiconducting substrates.

preprint1994arXiv

Total energy global optimizations using non orthogonal localized orbitals

An energy functional for orbital based $O(N)$ calculations is proposed, which depends on a number of non orthogonal, localized orbitals larger than the number of occupied states in the system, and on a parameter, the electronic chemical potential, determining the number of electrons. We show that the minimization of the functional with respect to overlapping localized orbitals can be performed so as to attain directly the ground state energy, without being trapped at local minima. The present approach overcomes the multiple minima problem present within the original formulation of orbital based $O(N)$ methods; it therefore makes it possible to perform $O(N)$ calculations for an arbitrary system, without including any information about the system bonding properties in the construction of the input wavefunctions. Furthermore, while retaining the same computational cost as the original approach, our formulation allows one to improve the variational estimate of the ground state energy, and the energy conservation during a molecular dynamics run. Several numerical examples for surfaces, bulk systems and clusters are presented and discussed.

Francesco Mauri

What is connected

Connect this record

See the researcher in context

Building this map preview

63 published item(s)

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

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

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

Wigner formulation of thermal transport in solids

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

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

Black metal hydrogen above 360 GPa driven by proton quantum fluctuations

Hybrid-functional electronic structure of multilayer graphene

Long-range rhombohedral-stacked graphene through shear

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

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

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

Quantum Crystal Structure in the 250 K Superconducting Lanthanum Hydride

Anharmonic effects in atomic hydrogen: superconductivity and lattice dynamical stability

Determination of scattering time and of valley occupation in transition-metal dichalcogenides doped by field effect

Spin susceptibility and electron-phonon coupling of two-dimensional materials by range-separated hybrid density functionals: Case study of Li$_x$ZrNCl

Two-dimensional Fröhlich interaction in transition-metal dichalcogenide monolayers: Theoretical modeling and first-principles calculations

Density-functional calculation of static screening in 2D materials: the long-wavelength dielectric function of graphene

Dissociation products and structures of solid H2S at strong compression

Electronic and vibrational properties of $\textup{TiSe}_2$ in the charge-density wave phase from first principles

First-principles theory of field-effect doping in transition-metal dichalcogenides: Structural properties, electronic structure, Hall coefficient, and electrical conductivity

Hydrogen sulphide at high pressure: a strongly-anharmonic phonon-mediated superconductor

Phonon effects on x-ray absorption and nuclear magnetic resonance spectroscopies

Quantum Hydrogen-Bond Symmetrization and High-Temperature Superconductivity in Hydrogen Sulfide

Strong anharmonicity induces quantum melting of charge density wave in 2H-NbSe$\_2$ under pressure

Universal enhancement of superconductivity in two dimensional semiconductors at low doping by electron-electron interaction

Downfolding electron-phonon Hamiltonians from ab-initio calculations: application to K$_3$Picene

Electrochemical doping of few layer ZrNCl from first-principles: electronic and structural properties in field-effect configuration

Electron-Phonon Interactions and the Intrinsic Electrical Resistivity of Graphene

First-principles calculations of phonon frequencies, lifetimes and spectral functions from weak to strong anharmonicity: the example of palladium hydrides

Phonon-limited resistivity of graphene by first-principle calculations: electron-phonon interactions, strain-induced gauge field and Boltzmann equation

Two-dimensional analysis of the double-resonant 2D Raman mode in bilayer graphene

Ab initio variational approach for evaluating lattice thermal conductivity

Anharmonic free energies and phonon dispersions from the stochastic self-consistent harmonic approximation: application to platinum and palladium hydrides

Anharmonic properties from a generalized third order ab~initio approach: theory and applications to graphite and graphene

Extended Czjzek model applied to NMR parameter distributions in sodium metaphosphate glass

First-Principles Theory of Anharmonicity and the Inverse Isotope Effect in Superconducting Palladium-Hydride Compounds

Quantum Monte Carlo study of the protonated water dimer

Hidden polymorphs drive the vitrification in B2O3

Local and non-local electron-phonon couplings in K3Picene and the effect of metallic screening

Zener tunneling in the electrical transport of quasi-metallic carbon nanotubes

Charge density wave and superconducting dome in TiSe2 from electron-phonon interaction

How to make graphene superconducting

Intercalant and intermolecular phonon assisted superconductivity in K-doped picene

Theory of double-resonant Raman spectra in graphene: intensity and line shape of defect-induced and two-phonon bands

Clar's Theory, STM Images, and Geometry of Graphene Nanoribbons

Current-voltage characteristics of graphene devices: interplay between Zener-Klein tunneling and defects

Doped graphene as tunable electron-phonon coupling material

Effects of magnetism and doping on the electron-phonon coupling in BaFe$_{2}$As$_{2}$

First-principles theory of the orbital magnetization

Phonon dispersion and low energy anomaly in CaC$_6$

Structure and stability of graphene nanoribbons in oxygen, carbon dioxide, water, and ammonia

Ab-initio study of gap opening and screening effects in gated bilayer graphene

Effect of dimensionality on the charge-density-wave in few-layers 2H-NbSe$_2$

Reducing the thermal conductivity of carbon nanotubes below the random isotope limit

Transport properties of graphene in the high-current limit

Born-Oppenheimer Breakdown in Graphene

Raman spectra of BN-nanotubes: Ab-initio and bond-polarizability model calculations

Structural relaxations around Ti, Cr and Fe impurities in alpha-Al2O3 probed by x-ray absorption near edge structure combined with first-principles calculations

All-electron magnetic response with pseudopotentials: NMR chemical shifts

First-principle study of excitonic self-trapping in diamond

Large scale quantum simulations: C_60 impacts on a semiconducting surface

Total energy global optimizations using non orthogonal localized orbitals