SciPost Submission Page
Superconductivity in RbH$_{12}$ at low pressures: an \emph{ab initio} study
by Đorđe Dangić, Yuewen Fang, Ion Errea
Submission summary
| Authors (as registered SciPost users): | Dorde Dangic |
| Submission information | |
|---|---|
| Preprint Link: | scipost_202507_00044v1 (pdf) |
| Date submitted: | July 16, 2025, 2:20 p.m. |
| Submitted by: | Dorde Dangic |
| Submitted to: | SciPost Physics |
| Ontological classification | |
|---|---|
| Academic field: | Physics |
| Specialties: |
|
| Approaches: | Theoretical, Computational |
Abstract
High-pressure polyhydrides are leading contenders for room temperature superconductivity. The next frontier lies in stabilizing them at ambient pressure, which would allow their practical applications. In this first-principles computational study, we investigate the potential for record-low pressure stabilization of binary superhydrides within the RbH$_{12}$ system, including lattice quantum anharmonic effects in the calculations. We identify five competing phases for the pressure range between 0 and 100 GPa. Incorporating anharmonic and quantum effects on ion dynamics, we find the $Immm$ and $P6_3/mmc$ phases to be the most probable, potentially metastable even at pressures as low as 10 GPa. Notably, all phases exhibit metallic properties, with critical temperatures between 50 and 100 K, within the pressure range where they are dynamically stable. These findings have the potential to inspire future experimental exploration of high-temperature superconductivity at low pressures in Rb-H binary compounds.
Author indications on fulfilling journal expectations
- Provide a novel and synergetic link between different research areas.
- Open a new pathway in an existing or a new research direction, with clear potential for multi-pronged follow-up work
- Detail a groundbreaking theoretical/experimental/computational discovery
- Present a breakthrough on a previously-identified and long-standing research stumbling block
Current status:
Reports on this Submission
Strengths
-
This is a very detailed, careful and accurate calculation of potential superconducting properties in a binary hydride with a potentially "low-pressure" high Tc
-
Both zero point fluctuations and anharmonicity, often overlooked in calculations of this kind, have been accounted for.
-
The paper is extremely well written and easy to read. Because of this, I would expect it to gain a wide audience, not only for its specific content (RbH12) but for the clear presentation of methodology and explanation of results which are generalizable to other compounds
Weaknesses
- The results themselves are somewhat underwhelming as they don't vary from earlier results even with all the extra work and detail supplied by SSCHA
Report
Requested changes
One extremely trivial matter:
In the last paragraph of the introduction the phrase "leaving *to* Immm and Cmcm phases to emerge" occurs. The "to" is odd. Is it supposed to be "two"? Or maybe "the"?
Recommendation
Publish (meets expectations and criteria for this Journal)
Strengths
- The authors employ an advanced method to study the superconducting properties of the RbH12 polyhydride. The calculations are sound.
- The authors employ a clever trick to estimate the Raman activity of phonon modes in the system despite it being a metal.
Weaknesses
- The article does not propose a new material, nor it improves substantially the understanding of an existing one.
- There are a few worrying methodological aspects. They are detailed in the requested changes. One important point is that the authors stretch the metastability range of this hydride, dismissing imaginary phonon frequencies as an interpolation error. The fact that SSCHA suggests that this system should be dynamically stable at a significantly lower pressure is one of the core points of the paper. It appears to be rather weak, in light of the fact that the interpolated phonon frequencies are imaginary.
Report
From the side of the importance, it is hard to evaluate this paper. Rubidium hydrides have been studied already from the computational and experimental side. This paper proposes an improved study, in the sense that it employs SSCHA to compute the phonon properties. On the other hand, the authors do not study the experimentally published phases, but rather study a phase predicted in a previous high-throughput study. The main effect to SSCHA is to lower the stabilization pressure of this known phase, and may motivate further experimental studies.
Overall, I think this article satisfies the general acceptance criteria, but I believe it does not meet the more stringent expectations of SciPost physics. There is no breakthrough, novel link, or groundbreaking discovery.
Requested changes
1. In the abstract, and later in the results, the authors write that the Immm and P63/mmc phases will be dynamically stable down to a certain pressure even though their phonon dispersion exhibit imaginary phonons. The reason, verbatim, is “these are probably interpolation issues”. They refer to the fact that a point for which the SSCHA matrices were computed has become stable, while another for which they were interpolated remained unstable. It is very hard to consider this argument as reliable evidence to support the claim that the structure is metastable at this pressure. I understand that the a larger supercell may be too expensive to compute, but this is not a valid argument to claim stability at that pressure.
2. The authors estimate the Raman activity of the phonon modes in a rather unconventional way. Since in metals the polarizability is not defined, they take a random tensor, symmetrize it, and compute some intensities. This is a clever trick, but I do not think that showing a simulated spectrum that is wrong (by the authors’ own admission) is a good idea. I suggest they just show the peak centers as vertical lines with different colors to mark which one is active or not (by symmetry).
3. The smearing on the DOS in Fig. 5 appears to be the 0.02 Ry used to compute the charge density. The result unfortunately is a DOS which is all smeared out. I suggest they employ the 0.008 Ry smearing which they used for the matrix elements, which also gives an indication if the chosen k-grid is reasonable.
4. There is an inconsistency between the y axis label of Fig. 1 and what the authors write (“Fig. 1 shows the Gibbs free energy”). I think it’s because they take it at T = 0, but distinguish it from enthalpy because there is the zero point energy. However this is all very misleading. At T = 0 it makes no sense to call this a free energy. The zero-point energy is not part of the free energy, even though one may include it in the expression for Fh.
5. Concerning Fig. 4, it is also rather worrying that the lines are so irregular, and that they include so few points, especially since these are done with plain DFT calculations which should be very cheap. Moreover, there should also be a “RbH+11H” line, or something similar, i.e. a line showing how much higher are the proposed structure from the most important points in the convex hull. Showing only the enthalpy (or free energy) compared between RbH12 phases may mislead the reader into thinking that there is one RbH12 phase on the hull even at ambient pressure, which is not the case.
6. A minor point is that the authors write that “a first ab initio prediction of high-temperature superconductivity was done for H3S, followed immediately by the experimental confirmation”. The real story is slightly more complicated, as the group of Eremets had been working on hydrogen sulfide independently, and tried multiple times to double-check the validity of the experiment. In fact, in the time between Duan’s and Eremets’ arXiv papers there are only a few days, which would not have been enough to do all those experiments. This was written by Eremets himself in Physics Reports 856, 1-78 (2020). I think it would be right to give the late Eremets credit for that.
7. The authors used a mu* of 0.1 to solve the isotropic Eliashberg equations. In Ref. Nature Reviews Physics 6 509–523 (2024) it is clearly argued that this is inappropriate if the goal is to compare the results with the McMillan formula, as that formula assumes a mu* with a different cutoff on the Matsubara frequencies. Please take it into account.
Recommendation
Accept in alternative Journal (see Report)
