Source author record

Ion Errea

Ion Errea 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

20works
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

20 published item(s)

preprint2026arXiv

The Impact of Ionic Anharmonicity on Superconductivity in Metal-Stuffed B-C Clathrates

Metal-stuffed B$-$C compounds with sodalite clathrate structure have captured increasing attention due to their predicted exceptional superconductivity above liquid nitrogen temperature at ambient pressure. However, by neglecting the quantum lattice anharmonicity, the existing studies may result in an incomplete understanding of such a lightweight system. Here, using state-of-the-art ab initio methods incorporating quantum effects and machine learning potentials, we revisit the properties of a series of $XY$$\text{B}_{6}\text{C}_{6}$ clathrates where $X$ and $Y$ are metals. Our findings show that ionic quantum and anharmonic effects can harden the $E_g$ and $E_u$ vibrational modes, enabling the dynamical stability of 15 materials previously considered unstable in the harmonic approximation, including materials with previously unreported ($XY$)$^{1+}$ state, which is demonstrated here to be crucial to reach high critical temperatures. Further calculations based on the anisotropic Migdal-Eliashberg equation demonstrate that the $T_\text{c}$ values for KRb$\text{B}_{6}\text{C}_{6}$ and Rb$\text{B}_{3}\text{C}_{3}$ among these stabilized compounds are 102 and 115 K at 0 and 15 GPa, respectively, both being higher than $T_\text{c}$ of 92 K of KPb$\text{B}_{6}\text{C}_{6}$ at the anharmonic level. These record-high $T_\text{c}$ values, surpassing liquid nitrogen temperatures, emphasize the importance of anharmonic effects in stabilizing B-C clathrates with large electron-phonon coupling strength and advancing the search for high-$T_\text{c}$ superconductivity at (near) ambient pressure.

preprint2022arXiv

Impact of ionic quantum fluctuations on the thermodynamic stability and superconductivity of LaBH$_8$

The recent prediction of a metastable high-symmetry Fm$\bar{3}$m phase of LaBH$_8$ gives hopes to reach high superconducting critical temperatures at affordable pressures among ternary hydrogen-rich compounds. Making use of first-principles calculations within density functional theory and the stochastic self-consistent harmonic approximation, we determine that ionic quantum fluctuations drive the system dynamically unstable below 77 GPa, a much higher pressure than the 45 GPa expected classically. Quantum anharmonic effects stretch the covalent B-H bond in the BH$_8$ units of the structure and, consequently, soften all hydrogen-character modes. Above 77 GPa Fm$\bar{3}$m LaBH$_8$ remains metastable and, interestingly, its superconducting critical temperature is largely enhanced by quantum anharmonic effects, reaching critical temperatures around 170 K at the verge of the dynamical instability. Our results suggest that low pressure metastable phases with covalently bonded symmetric XH$_8$ units will be destabilized by ionic quantum fluctuations.

preprint2022arXiv

Superconducting hydrides on a quantum landscape

Reaching superconductivity at ambient conditions is one of the biggest scientific dreams. The discoveries in the last few years at high pressures place hydrogen-based compounds as the best candidates for making it true. As the recent history shows, first-principles calculations are expected to continue guiding the experimental quest in the right track in the coming years. Considering that ionic quantum fluctuations largely affect the structural and vibrational properties of superconducting hydrides, in many cases making them stable at much lower pressures than expected, it will be crucial to include such effects on the future ab initio predictions. The prospects for low-pressure high critical-temperature compounds are wide open, even at ambient pressure.

preprint2021arXiv

Quantum anharmonic enhancement of superconductivity in $P6_3/mmc$ ScH$_6$ at high pressures: a first-principles study

Making use of first-principles calculations, we analyze the effect of quantum ionic fluctuations and lattice anharmonicity on the crystal structure and superconductivity of P63/mmc ScH6 in the 100-160 GPa pressure range within the stochastic self-consistent harmonic approximation. We predict a strong correction to the crystal structure, the phonon spectra, and the superconducting critical temperatures, which have been estimated in previous calculations without considering ionic fluctuations on the crystal structure and assuming the harmonic approximation for the lattice dynamics. Quantum ionic fluctuations have a large impact on the H2 molecular-like units present in the crystal by increasing the hydrogen-hydrogen distance about a 5%. According to our anharmonic phonon spectra, this structure will be dynamically stable at least above 85 GPa, which is 45 GPa lower than the pressure given by the harmonic approximation. Contrary to many superconducting hydrogen-rich compounds, where quantum ionic effects and the consequent anharmonicity tend to lower the superconducting critical temperature, our results show that it can be enhanced in P63/mmc ScH6 by approximately a 15%. We attribute the enhancement of the critical temperature to the stretching of the H$_2$ molecular-like units and the associated increase of the electron-phonon interaction. Our results suggest that quantum ionic effects increase the superconducting critical temperature in hydrogen-rich materials with H2 units by increasing the hydrogen-hydrogen distance and, consequently, the electron-phonon interaction.

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

Infrared permittivity of the biaxial van der Waals semiconductor $α$-MoO$_3$ from near- and far-field correlative studies

The biaxial van der Waals semiconductor $α$-phase molybdenum trioxide ($α$-MoO$_3$) has recently received significant attention due to its ability to support highly anisotropic phonon polaritons (PhPs) -infrared (IR) light coupled to lattice vibrations in polar materials-, offering an unprecedented platform for controlling the flow of energy at the nanoscale. However, to fully exploit the extraordinary IR response of this material, an accurate dielectric function is required. Here, we report the accurate IR dielectric function of $α$-MoO$_3$ by modelling far-field, polarized IR reflectance spectra acquired on a single thick flake of this material. Unique to our work, the far-field model is refined by contrasting the experimental dispersion and damping of PhPs, revealed by polariton interferometry using scattering-type scanning near-field optical microscopy (s-SNOM) on thin flakes of $α$-MoO$_3$, with analytical and transfer-matrix calculations, as well as full-wave simulations. Through these correlative efforts, exceptional quantitative agreement is attained to both far- and near-field properties for multiple flakes, thus providing strong verification of the accuracy of our model, while offering a novel approach to extracting dielectric functions of nanomaterials, usually too small or inhomogeneous for establishing accurate models only from standard far-field methods. In addition, by employing density functional theory (DFT), we provide insights into the various vibrational states dictating our dielectric function model and the intriguing optical properties of $α$-MoO$_3$.

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.

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

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.

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.

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

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.

preprint2011arXiv

Anharmonic stabilization of the high-pressure simple cubic phase of calcium

The phonon spectrum of the high-pressure simple cubic phase of calcium, in the harmonic approx- imation, shows imaginary branches that make it mechanically unstable. In this letter, the phonon spectrum is recalculated using density-functional theory (DFT) ab initio methods fully including anharmonic effects up to fourth order at 50 GPa. Considering that perturbation theory cannot be employed with imaginary harmonic frequencies, a variational procedure based on the Gibbs- Bogoliubov inequality is used to estimate the renormalized phonon frequencies. The results show that strong quantum anharmonic effects make the imaginary phonons become positive even at zero temperature so that the simple cubic phase becomes mechanically stable, as experiments suggest. Moreover, our calculations find a superconducting Tc in agreement with experiments and predict an anomalous behavior of the specific heat.