Although the evidence is still limited, a growing body of research suggests music may have beneficial effects for diseases such as Parkinson’s.
Image credit: Shutterstock/agsandrew.
New Research In
Physical Sciences
Featured Portals
Articles by Topic
Biological Sciences
Featured Portals
Articles by Topic
Contributed by Ho-Kwang Mao, December 15, 2009 (sent for review November 28, 2009)
Abstract
The long-standing prediction that hydrogen can assume a metallic state under high pressure, combined with arguments put forward more recently that this state might even be superconducting up to high temperatures, continues to spur tremendous research activities toward the experimental realization of metallic hydrogen. These efforts have however so far been impeded by the enormous challenges associated with the exceedingly large required pressure. Hydrogen-dense materials, of the MH4 form (where M can be, e.g., Si, Ge, or Sn) or of the MH3 form (with M being, e.g., Al, Sc, Y, or La), allow for the rather exciting opportunity to carry out a proxy study of metallic hydrogen and associated high-temperature superconductivity at pressures within the reach of current techniques. At least one experimental report indicates that a superconducting state might have been observed already in SiH4, and several theoretical studies have predicted superconductivity in pressurized hydrogen-rich materials; however, no systematic dependence on the applied pressure has yet been identified so far. In the present work, we have used first-principles methods in an attempt to predict the superconducting critical temperature (Tc) as a function of pressure (P) for three metal-hydride systems of the MH3 form, namely ScH3, YH3, and LaH3. By comparing the obtained results, we are able to point out a general trend in the Tc-dependence on P. These gained insights presented here are likely to stimulate further theoretical studies of metallic phases of hydrogen-dense materials and should lead to new experimental investigations of their superconducting properties.
Hydrogen, while isoelectronic to the alkali metals, is in fact insulating under ambient conditions, which is a consequence of the comparatively high binding energy of its electron. Enormous pressure would be required to reach a hypothetical metallic phase of hydrogen (1), which remains so far elusive to the available experimental techniques. Most intriguingly, it was proposed about a decade ago that metallic hydrogen could also be superconducting up to room temperature (2), inspiring numerous studies in the direction toward realizing this goal. Hydrogen-dense materials have been recognized as a potential route to investigate such phenomena, because in these metal alloys, the combination of chemical precompression (3) from the metal host and of experimentally applied pressures could facilitate hydrogen densities in a range where metallization can occur. For example, metallization in silane (SiH4) has been experimentally observed (4, 5) and was analyzed from a theoretical perspective by us (6).
The emergence of high-temperature superconducting phases in several pressurized hydrogen-dense materials has been predicted (7–12), but, until now, no systematic study had been carried out on how the critical temperature (Tc) of superconductivity evolves under the applied pressure. Here, we demonstrate on the examples of ScH3, YH3, and LaH3, that the respective pressure dependencies of Tc exhibit remarkably similar behavior, hinting at a general rule for the MH3 class of hydrogen-dense materials.
Many metal-hydride alloys investigated for metallization and superconductivity possess either the stoichiometry MH3 or MH4. The mechanism by which they are expected to conduct is markedly different: In the MH3-type systems, M and H form a 3-D network, and both participate in the conduction process, whereas in the MH4-type systems, M and H separate into individual layers, and conduction occurs only through the H-layers. The three isoelectronic trihydride systems that we have concentrated on in the present work, crystallize in the fcc structure when sufficiently pressurized, with the metal located at the fcc lattice sites, and two distinct groups of hydrogen atoms positioned at the octahedral (HO) and tetrahedral interstitial sites (HT), respectively. The M - HT bonds form a MH2-network that acts as a cage structure, encapsulating the HO.
Results and Discussion
The main result of our study is presented in Fig. 1A, where we provide a summary of our calculated Tc values for ScH3, YH3, and LaH3 as a function of pressure. It is important to note that several other theoretical works have predicted isolated Tc values for several tetrahydrides (7–9) and for one trihydride (AIH3) (11), the latter having only recently been found to possess a metallic phase (13–15). Typically, these studies predicted Tc only at very few pressure values, making it difficult to discern a particular trend. For the sake of facilitating a more direct comparison, these data are put side by side in Fig. 1B with our own data. It becomes noticeable that for the MH3-type systems investigated by us, the pressure at which superconductivity sets in (10–20 GPa) is nearly an order of magnitude lower than that for the MH4 compounds. Around 10 GPa, fcc-LaH3 is dynamically stabilized and possesses its highest Tc of 20 K. Slightly below 20 GPa, fcc-ScH3 and fcc-YH3 are stabilized and mark their highest Tc, 18 K and 40 K, respectively. This phonon-mediated superconductivity is seen to become reduced under increasing pressure, and eventually, all three MH3 systems become normal metallic (i.e., not superconducting) around 35 GPa. In our calculations, we checked the dynamical stability of the three metal hydrides up to much higher pressure and for ScH3 and YH3 we did not observe any dynamical instability up to 80 GPa, but it is worth noting that LaH3 is dynamically destabilized above 65 GPa, which implies that another high pressure phase of LaH3 could exist in that pressure regime. We will return to a more detailed discussion of these results below. First however, let us look at some other calculated properties of the three MH3 systems.
Fig. 1.
Calculated critical temperature of superconductivity (Tc) as a function of pressure. (A) Our predictions for ScH3, YH3 (12), and LaH3. (B) The data for GeH4 (8), SiH4 (9, 10), SnH4 (7), and AIH3 (11) (Hollow Symbols), as predicted by various groups, compared with our data for the three trihydride systems (Filled Symbols). The approximate grouping of the MH3 and MH4 Tc - P data pairs is respectively highlighted by the yellow and green colored areas.
Fig. 2 shows the phonon dispersion curves that we obtained from our lattice dynamics calculations for ScH3, YH3, and LaH3. For all three materials, the dispersion relation was calculated at a respective pressure that lies just slightly above the stability threshold for the corresponding metal hydride. This ensures that the phonon dispersions are captured near the onset of dynamical stability. For the case of ScH3, dynamical stability sets in for pressures > 18 GPa. As it can be seen from the violet dashed line (drawn only for the acoustic mode below the transition pressure in the destabilized dispersion regime) in the corresponding plot of the phonon dispersion for ScH3 in Fig. 2A, the instability is driven by the softening at both the Γ-point and the X-point. In contrast to that, one can note from Fig. 2B and C that for YH3 and LaH3, the instability due to softening occurs at or near the K-point. It becomes apparent that the phonon dispersion curves for YH3 and LaH3 are highly similar to each other, whereas they differ substantially from ScH3. One of the main disparities lies in the much stronger dispersion at the Γ-point in ScH3 than in YH3 and LaH3, which can be seen to even reach down to the acoustic phonon bands.
Fig. 2.
Phonon dispersion curves at the lowest possible stabilizing pressure plotted along high symmetry lines for (A) ScH3 at 18 GPa, (B) YH3 at 18 GPa, and (C) LaH3 at 11 GPa. In all three panels, the black solid lines show the stabilized phonon dispersion relation, whereas the dashed violet line illustrates the destabilized phonon dispersion relation (shown only for the acoustic mode) at a slightly lower pressure in the dynamically unstable regime. The colored filled circles highlight points of strong electron–phonon coupling, with the radius proportional to the respective coupling strengths.
Whereas there exists no comprehensive experimental study for the low-temperature phase diagram of these trihydrides, the calculated phonon softening near the K-point in YH3 and LaH3 indicates equivalent structural transformations in both. For ScH3, the situation regarding the evolution of structural changes with pressure is less clear, and it could very well be that another crystal structure may be stabilized in the pressure range below 18 GPa.
The filled circles in Fig. 2 mark the electron–phonon coupling with the radius proportional to their respective normalized strength. We see that only in the case of ScH3, the electron–phonon coupling is very strong around the Γ-point, whereas for YH3 and LaH3, the strongest coupling occurs around the L-point.
As it could already be seen above from the phonon dispersion relations (Fig. 2), the behavior of ScH3 is markedly different from YH3 and LaH3. A similar conclusion can also be reached from Fig. 3, where we plot the calculated spectral functions for all three systems. We note that the acoustic band and optical band are merged together in the case of ScH3 (Fig. 3A). Thus, Sc–Sc vibrations and Sc–H vibrations are seen to contribute equally to the spectral function. This is in contrast to the situation for YH3 (Fig. 3B) and LaH3 (Fig. 3C) that once again appear rather similar to each other. For a higher pressure level, the spectral function is seen to go to zero for ScH3 and LaH3, whereas it reappears with a finite value only for YH3. This helps to explain the evolution of Tc with pressure as shown in Fig. 1.
Fig. 3.
The spectral function α2F is shown as colored shaded areas together with the phonon density of states (Black Curves) as a function of frequency for (A) ScH3 at 18 GPa, (B) YH3 at 18 GPa, and (C) LaH3 at 11 GPa.
Let us now return to the discussion of the Tc vs. P data of Fig. 1A. Our data for the trihydrides clearly show that Tc is always highest in the vicinity of a structural transformation. This is a strong indication that the superconducting mechanism is driven by electron-phonon coupling. As we recently discussed just for the case of YH3 (12), at the lowest possible stabilizing pressure, the electron–phonon coupling is strong for Y and Y-HO vibrations, which couple with a high density of electronic states near the Fermi level belonging to the atomic-like HO. With rising pressure, the coupling diminishes and the superconducting state eventually ceases to exist. Combined with the new results that we obtained for ScH3 and LaH3, these data suggest the rather intriguing conclusion that this general trend does in fact hold true for all MH3 systems studied by us.
We note that for higher pressures (> 50 GPa), the superconducting phase reemerges in YH3, but not so in ScH3 and LaH3. This second superconductivity region can be mostly attributed to Y-HO vibrations coupling with the electronic states near the Fermi level (12). Because Sc, Y, and La are isoelectronic, the only essential difference between the three trihydrides is the mass of the metal. It thus appears that the mass of Y happens to fall into just the right range to generate vibrational energies of the Y-HO mode that can give rise to this second superconducting regime. Sc and La, on the other hand, seem to have masses either too low or too high that lead to nonmatching vibrational energies, thus preventing the possibility of a reemerging second superconducting regime at higher pressure.
Conclusions
In summary, based on our first-principles calculations, we were able to identify a general trend of pressure-induced superconductivity for the metal trihydrides ScH3, YH3, and LaH3, and drew comparison with other hydrogen-dense materials, including those of the MH4 form. We found that the three transition metal hydrides studied by us possess the highest Tc at the vicinity of a structural transformation to the cubic (fcc) phase and that Tc always drops as the pressure increases, until they turn into a normal metallic state (i.e., not superconducting). A secondary superconductivity regime is seen to arise only for YH3.
As with any other theoretical prediction, we also naturally hope that our work will not merely lead to further computational studies of related metallic hydrogen-dense materials, but actually inspire experimental measurements to test our data. Certainly, the pressure range of 10–80 GPa, for which we have predicted here superconducting states in ScH3, YH3, and LaH3, is easily accessible with contemporary diamond anvil cell techniques. Rather, the challenge might actually lie in preparing metal-hydride samples of the proper stoichiometry. Even slight deviations of the perfect 1∶3 ratio of metal to hydrogen could prevent the superconducting state from emerging. Nevertheless, we hope that experimentalists will embrace this challenge to test our predictions, because we strongly believe that there is tremendous merit for the understanding of the underlying physics in metallic hydrogen-rich systems if such measurements could indeed be carried out.
Methods
Density Functional Theory and Lattice Dynamics Calculations.
Electronic structure calculations are based on the generalized gradient approximation (16) with the Perdew–Burke–Ernzerhof parametrization (17) for the exchange-correlation functional to density functional theory (18, 19). Phonon dispersion and electron–phonon calculations were performed with density functional perturbation theory (20) using the program Quantum ESPRESSO with a kinetic energy cutoff of 60 Ryd and a 24 × 24 × 24 Monkhorst–Pack mesh (21) for Brillouin zone integration. Electron–phonon coupling matrix elements were computed in the first Brillouin zone on a 8 × 8 × 8 q-mesh.
Tc Calculations.
We calculated the critical temperature for superconductivity (Tc) using the Allen–Dynes equation (22) 
[1]In terms of the normalized weighting function that arises in Eliashberg theory, 
[2]The parameter λ is a dimensionless measure of the strength of α2F 
[3]and the logarithmic average frequency ωlog 
[4]
Acknowledgments
We thank the Swedish Foundation for International Cooperation in Research and Higher Education, the Swedish Research Council, Futura Foundation, and Wenner–Gren Foundations for financial support. The Swedish National Infrastructure for Computing and the Uppsala Multidisciplinary Center for Advanced Computational Science are gratefully acknowledged for providing computing time. This material is based upon work supported as part of the EFree, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences Award DE-SC0001057.
Footnotes
- 1To whom correspondence may be addressed. E-mail: rajeev{at}fysik.uu.se or hmao{at}ciw.edu.
Author contributions: D.Y.K., H.-K.M., and R.A. designed research; D.Y.K. performed research; T.W.K. contributed new reagents/analytic tools; D.Y.K., R.H.S., H.-K.M., T.W.K., and R.A. analyzed data; and D.Y.K., R.H.S., H.-K.M., and R.A. wrote the paper; .
The authors declare no conflict of interest.
References
General trend for pressurized superconducting hydrogen-dense materials
Sign up for Article Alerts
Article Classifications
- Physical Sciences
- Physics
Jump to section
You May Also be Interested in
Biomedical communities and journals need to standardize nomenclature of gene products to enhance accuracy in scientific and public communication.
Image credit: Shutterstock/greenbutterfly.
Shapeshifting designs could have wide-ranging pharmaceutical and biomedical applications in coming years.
Image credit: Acacia Dishman/Medical College of Wisconsin.
Amanda Rodewald, Ivan Rudik, and Catherine Kling talk about the hazards of ozone pollution to birds.
CRISPR-Cas9 gene editing can improve the effectiveness of spermatogonial stem cell transplantation in mice and livestock, a study finds.
Image credit: Jon M. Oatley.