Source author record

Matteo Calandra

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

43works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

43 published item(s)

preprint2022arXiv

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

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

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

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

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

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

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

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

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

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

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

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 phonons in few layer MoS$_2$: Raman spectroscopy of ultra low energy compression and shear modes

Molybdenum disulfide (MoS$_2$) is a promising material for making two-dimensional crystals and flexible electronic and optoelectronic devices at the nanoscale. MoS$_2$ flakes can show high mobilities and have even been integrated in nanocircuits . A fundamental requirement for such use is efficient thermal transport. Electronic transport generates heat which needs to be evacuated, more crucially so in nanostructures. Anharmonic phonon-phonon scattering is the dominant intrinsic limitation to thermal transport in insulators. Here, using appropriate samples, ultra-low energy Raman spectroscopy and first principles calculations, we provide a full experimental and theoretical description of compression and shear modes of few-layer (FL) MoS$_2$. We demonstrate that the compression modes are strongly anharmonic with a marked enhancement of phonon-phonon scattering as the number of layers is reduced, most likely a general feature of nanolayered materials with weak interlayer coupling.

preprint2013arXiv

Chemically-exfoliated single-layer MoS$_2$ : stability, lattice dynamics and catalytic adsorption from first principles

Chemically and mechanically exfoliated MoS$_2$ single-layer samples have substantially different properties. While mechanically exfoliated single-layers are mono-phase (1H polytype with Mo in trigonal prismatic coordination), the chemically exfoliated samples show coexistence of three different phases, 1H, 1T (Mo in octahedral coordination) and 1T$^{'}$ (a distorted $2\times 1$ 1T-superstructure). By using first-principles calculations, we investigate the energetics and the dynamical stability of the three phases. We show that the 1H phase is the most stable one, while the metallic 1T phase, strongly unstable, undergoes a phase transition towards a metastable and insulating 1T$^{'}$ structure composed of separated zig-zag chains. We calculate electronic structure, phonon dispersion, Raman frequencies and intensities for the 1T$^{'}$ structure. We provide a microscopical description of the J$_1$, J$_2$ and J$_3$ Raman features first detected more then $20$ years ago, but unexplained up to now. Finally, we show that H adsorbates, that are naturally present at the end of the chemical exfoliation process, stabilize the 1T$^{\prime}$ over the 1H one.

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

Projector augmented wave calculation of x-ray absorption spectra at the L2,3 edges

We develop a technique based on density functional theory and the projector augmented wave method in order to obtain the x-ray absorption cross section at a general edge, both in the electric dipole and quadrupole approximations. The method is a generalization of Taillefumier et al., PRB 66, 195107 (2002). We apply the method to the calculation of the Cu L2,3 edges in fcc copper and cuprite (Cu2O), and to the S L2,3 edges in molybdenite (2H-MoS2). The role of core-hole effects, modeled in a supercell approach, as well as the decomposition of the spectrum into different angular momentum channels are studied in detail. In copper we find that the best agreement with experimental data is obtained when core-hole effects are neglected. On the contrary, core-hole effects need to be included both in Cu2O and 2H-MoS2. Finally we show that a non-negligible component of S L2,3 edges in 2H-MoS2 involves transition to states with s character at all energy scales. The inclusion of this angular momentum channel is mandatory to correctly describe the angular dependence of the measured spectra. We believe that transitions to s character states are quantitatively significant at the L2,3 edges of third row elements from Al to Ar.

preprint2012arXiv

Anharmonic suppression of Charge density wave in 2H-NbS$_2$

The temperature dependence of the phonon spectrum in the superconducting transition metal dichalcogenide 2H-NbS$_2$ is measured by diffuse and inelastic x-ray scattering. A deep, wide and strongly temperature dependent softening, of the two lowest energy longitudinal phonons bands, appears along the $\mathrm{ΓM}$ symmetry line in reciprocal space. In sharp contrast to the iso-electronic compounds 2H-NbSe$_2$, the soft phonons energies are finite, even at very low temperature, and no charge density wave instability occurs, in disagreement with harmonic ab-initio calculations. We show that 2H-NbS$_2$ is at the verge of the charge density wave transition and its occurrence is only suppressed by the large anharmonic effects. Moreover, the anharmonicity and the electron phonon coupling both show a strong in-plane anisotropy.

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.

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.

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

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.

preprint2009arXiv

Angular dependence of core hole screening in LiCoO2: A DFT+U calculation of the oxygen and cobalt K-edge x-ray absorption spectra

Angular dependent core-hole screening effects have been found in the cobalt K-edge x-ray absorption spectrum of LiCoO2, using high-resolution data and parameter-free GGA+U calculations. The Co 1s core-hole on the absorber causes strong local attraction. The core-hole screening on the nearest neighbours cobalt induces a 2 eV shift in the density of states with respect to the on-site 1s-3d transitions, as detected in the Co K pre-edge spectrum. Our DFT+U calculations reveal that the off-site screening is different in the out-of-plane direction, where a 3 eV shift is visible in both calculations and experiment. The detailed analysis of the inclusion of the core-hole potential and the Hubbard parameter U shows that the core-hole is essential for the off-site screening, while U improves the description of the angular dependent screening effects. In the case of oxygen K-edge, both the core-hole potential and the Hubbard parameter improve the relative positions of the spectral features.

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.