On the Possibility of Metastable Metallic Hydrogen
On the Possibility of Metastable Metallic Hydrogen Jeffrey M. McMahon
Department of Physics & Astronomy
March 14, 2017
Jeffrey M. McMahon (WSU)
March 14, 2017
1 / 19
The Problem of Metallic Hydrogen In 1935, Wigner and Huntington predicted1 that sufficient pressure would dissociate hydrogen molecules; and any Bravais lattice of such atoms would be metallic. Initial interest was primarily related to astrophysical problems ... ... later to fundamental physics; also with practical applications: • It is expected to have remarkable properties: I I I
High-temperature superconductivity2,3 Novel types of quantum fluids3 Powerful rocket fuel
• Theoretical capabilities have progressed4 ; properties can be accurately calculated from first principles
• ... As have experimental ones; it is now possible to create in the laboratory
• It may be metastable6 1
E. Wigner and H. B. Huntington, J. Chem. Phys. 3, 764–770 (1935) N. W. Ashcroft, Phys. Rev. Lett. 21, 1748–1749 (1968) solid 3 E. Babaev, A. Sudbo, and N. W. Ashcroft, Nature 431, 666-668 (2004) 4 J. M. McMahon, M. A. Morales, C. Pierleoni, and D. M. Ceperley, Rev. Mod. Phys. 84, 1607–1653 (2012) 5 R. P. Dias and I. F. Silvera, Science (2017) 6 E. G. Brovman, Yu. Kagan, and A. Kholas, JETP 34, 1300–1315 (1972) 2
Jeffrey M. McMahon (WSU)
Introduction
metallic hydrogen5
March 14, 2017
2 / 19
Metastable Metallic Hydrogen There is little doubt that hydrogen becomes metallic at high pressures. For practical (terrestrial) applications, ... ... the significant, outstanding1,2 question is whether it is metastable. Note: The answer to this presupposes solutions to the following: (1) Determination of the minimum-energy crystal structure(s), ... ... and proof that it lies at a stationary point (2) Proof of (dynamic) stability (3) Analysis of the relation between the ground- and metastable-state structures of metallic hydrogen, and the molecular phase (4) Determination of the lifetime of the metastable state
It is the purpose of this talk to answer this question. 1 2
E. G. Brovman, Yu. Kagan, and A. Kholas, JETP 34, 1300–1315 (1972) H. Y. Geng, H. X. Song, J. F. Li, and Q. Wu, J. Appl. Phys. 111, 063510 (2012) Jeffrey M. McMahon (WSU)
Introduction
March 14, 2017
3 / 19
Metastable Metallic Hydrogen (Continued)
Brief answer:
Metallic hydrogen is metastable ...
... but, only to (approximately) 250 GPa.
Jeffrey M. McMahon (WSU)
Introduction
March 14, 2017
4 / 19
Outline
This talk is outlined as follows: • Phase diagram of hydrogen • Crystal structure(s) of metallic hydrogen I Ground-state structures of atomic metallic hydrogen I Searching for (new) metastable metallic structures • Dynamic stabilities • Concluding remarks
Jeffrey M. McMahon (WSU)
Introduction
March 14, 2017
5 / 19
Phase Diagram of Hydrogen
Temperature (K)
1500
liquid H2 liquid H
1000
500
I
I'? IV V III
II 0
10
VI H2-PRE
100
1000
solid H 10000
Pressure (GPa) Figure: Combined experimental (lines) and theoretical (lines with points) phase diagram of hydrogen. Data compiled from several sources Jeffrey M. McMahon (WSU)
Phase Diagram
March 14, 2017
6 / 19
Ground-state Structures of Atomic Metallic Hydrogen We earlier found1 several candidate structures of atomic metallic hydrogen:
(Not shown are:
• High-pressure molecular phases, since found2 • Ultrahigh-pressure phases, since found2 • Nearly-degenerate, lower-symmetry distortions3 , not originally reported in Ref. [1])
m
m
Figure: Ground-state enthalpies of the crystal structures of atomic metallic hydrogen relative to fcc ∆H, not including proton zero-point energy. The inset shows an expanded view of the ultrahigh-pressure region.
1
J. M. McMahon and D. M. Ceperley, Phys. Rev. Lett. 106, 165302 (2011) H. Liu, H. Wang, and Y. Ma, J. Phys. Chem. C 116, 9221–9226 (2012) 3 H. Y. Geng, H. X. Song, J. F. Li, and Q. Wu, J. Appl. Phys. 111, 063510 (2012)
2
Jeffrey M. McMahon (WSU)
Atomic Metallic Hydrogen
March 14, 2017
7 / 19
Tetragonal Family of Structures The most promising candidates structures for higher-pressure, metastable metallic hydrogen belong to the tetragonal family. Examples:
x y c/a Diamond
β-Sn
Cs-IV
Fd-3m (c/a = 1.414)
I41 /amd (c/a ∼ 0.55 < 1)
I41 /amd (c/a ∼ 3.73 > 1)
(These structures differ only in their c/a ratio.) Jeffrey M. McMahon (WSU)
Tetragonal Structures
March 14, 2017
8 / 19
0 K Phase Diagram of the Tetragonal Structures The 0 K phase diagram of the tetragonal structures reveals general, qualitative features of metastable metallic hydrogen: Regions of distinct qualitative behavior are indicated with dotted lines:
200
∆H (meV/proton)
100
• > 300 GPa: The (relative) enthalpies re-
0 -100
•
-200 c/a > 1
-300 -400 0
100
200
300
400
500
P (GPa)
•
Figure: 0 K phase diagram of the tetragonal structures of metallic hydrogen.
Jeffrey M. McMahon (WSU)
Tetragonal Structures
main relatively flat (with the exception of diamond). 200–300 GPa: The enthalpies of the βSn, diamond, and Cs-IV structures become nearly degenerate; those of c/a 1 and 1 begin to decrease and increase, respectively. < 200 GPa: The enthalpies of the β-Sn, diamond, and Cs-IV structures becomes degenerate; those of c/a 1 and 1 decrease and increase sharply.
March 14, 2017
9 / 19
Energies of the Tetragonal Structures This behavior is understood by considering the energies of the tetragonal structures over all c/a: Notice that:
• At 500 GPa: The diamond structure forms
E (eV/proton)
-13
•
-14
• -15
-16
3
1.15017 Å (500 GPa) 3 1.40361 Å (300 GPa) 3 1.62966 Å (200 GPa) 3 2.08067 Å (100 GPa) 3 4.72262 Å (0 GPa)
0
2
4
6
8
General feature: The collapse of energy barriers, between 200–300 GPa. 10
c/a
Figure: Energies of the tetragonal structures of metallic hydrogen, as a function of c/a. Fixed volumes are shown, calculated as an average of each structure at the pressure indicated in parenthesis. The c/a 1, β-Sn, diamond, Cs-IV, and 1 structures are indicated with arrows (in order of increasing c/a). Jeffrey M. McMahon (WSU)
the energy barrier (a maximum) between the β-Sn and Cs-IV structures. ... Between 300–200 GPa: The energy barrier collapses (appreciably). ... By 100 GPa: The energy barrier has collapsed, leaving a (single) minimum.
• The increase in energy as c/a → 0 and •
Tetragonal Structures
→ ∞ are due to the coming together of Brillouin planes with like charges. ... This explains the sharp decrease in the c/a 1 structure (nearest-neighbor distance: 0.9900Å): General feature: The tendency to form molecules, below 200 GPa. March 14, 2017
10 / 19
Energy of Metallic Hydrogen There is no guarantee that the high-pressure structures are representative of those (potentially) metastable at lower pressures. Consider the energy E (in the adiabatic approximation): E = Estatic + EZPE • Estatic (static energy) favors anisotropic structures (in metallic hydrogen). • EZPE (zero-point energy) favors symmetric structures. • At low pressures: The energy is determined largely by Estatic . • With increasing pressure: The importance of EZPE increases.
We performed searches1 for metallic (atomic and/or mixed atomic/molecular) structures with the lowest Estatic at 0 GPa. 1
C. J. Pickard and R. J. Needs, J. Phys. Condens. Matter 23, 053201 (2011) Jeffrey M. McMahon (WSU)
Energy of Metallic Hydrogen
March 14, 2017
11 / 19
Relative Enthalpies of Candidate Structures Our searches revealed several candidate structures of (potentially) metastable metallic hydrogen. The results are (qualitatively) understood, and further support the prior results, by considering only the most stable structures. Example:
∆H (meV/proton)
100 0
Regions of distinct behavior (dotted lines):
-100
• < 200 GPa: Molecular structures become
-200 -300 Cs-IV 1H 3H 5H 7H 9H
-400 -500 -600 -700
0
100
200
300
400
• 500
P (GPa)
Figure: Relative enthalpies ∆H (to Cs-IV) of the most stable metallic structures, from searches over different numbers of atoms in the unit cell. Jeffrey M. McMahon (WSU)
Candidate Structures
•
favored. Indicated by sharply decreasing enthalpies, with a magnitude proportional to the number of atoms in the search. 200–300 GPa: Transition region. Energy barriers between (atomic) metallic phases form. Indicated by crossovers. > 300 GPa: Energy barriers between metallic phases exist. Structures found at 0 GPa are also no longer representative. Indicated by lack of trend.
March 14, 2017
12 / 19
The Tendency Towards Molecules Below 200 GPa: All structures found show a tendency towards the formation of molecules. Example: As the number of atoms increases, there is a greater proportion of molecules (to atoms) formed:
(a) C2 (5)
(b) C2/m (7)
(c) Pm (9)
Figure: Most stable candidate structures of metastable metallic hydrogen at ambient conditions, for searches with different numbers of atoms (shown in parenthesis). Jeffrey M. McMahon (WSU)
Candidate Structures
March 14, 2017
13 / 19
The Tendency Towards Molecules (Continued) Example: This tendency is seen even in the (completely) atomic structures (found in searches over a lower number of atoms):
(a) Pmmm
(b) Cmmm
(c) Immm
(d) P6/mmm
Figure: Most stable (nearly degenerate) candidate structures of metastable atomic metallic hydrogen at ambient conditions.
(Notice the proximity of atoms; distances are on the order of 0.9914 Å.)
These results suggest that the atomic hydrogen at lower pressures has no region(s) of stability. Jeffrey M. McMahon (WSU)
Candidate Structures
March 14, 2017
14 / 19
Dynamic Stability A crystal is dynamically (mechanically) stable, if it executes a stable oscillatory motion about its equilibrium configuration. By solving the equation of motions for the normal modes of vibration, ... ... this is equivalent to the condition that their frequencies ω satisfy: ωj (q)2 ≥ 0 for all phonon branches j and wavevectors q. (For an imaginary frequency means that the system, subject to a small displacement, will disrupt exponentially with time.)
Dynamic stabilities were determined by calculating the phonon density of states F(ω) for each structure. M. Born and K. Huang, “Dynamical Theory of Crystal Lattices” (1954) Jeffrey M. McMahon (WSU)
Dynamic Stabilities
March 14, 2017
15 / 19
Dynamic Stabilities of the Tetragonal Structures The Cs-IV structure is dynamically stable below molecular dissociation: 0.004 500 GPa 350 GPa 250 GPa
0.008
F(ω) / proton
F(ω) / proton
0.003
0.002
0.001
0
200 GPa 100 GPa 50 GPa 0 GPa
0.006
0.004
0.002
0
500
1000
1500
2000
2500
3000
0 -2000
-1000
0
1000
2000
3000
-1
-1
ω (cm )
ω (cm )
(a) stable region
(b) unstable region
Figure: Phonon density of states F(ω) of Cs-IV, at several pressures. Pressure regions of (a) stability and (b) instability are shown. The dotted line in (b) is used to separate the stable from unstable frequencies. The arrow also in this plot indicates the growing instability, with decreasing pressure.
... but not to ambient conditions. (All other tetragonal structures are dynamically unstable.) Jeffrey M. McMahon (WSU)
Dynamic Stabilities
March 14, 2017
16 / 19
Dynamic Stabilities of the Candidate Structures The candidate structures are dynamically unstable, at all pressures: 0.020 0.0005
Pmmm Cmmm Immm P6/mmm
0.0020
0.015
Cmmm Immm P6/mmm
0.010 0.0002
F(ω) / proton
0.0003
F(ω) / proton
0.0004
0.0015
0.0010
0.0001
0.005 0.0005 0
-6000
-4000
-2000
0
2000
0.000 4000
0.0000 -2000
-1000
-1
0
1000
2000
3000
4000
-1
ω (cm )
ω (cm )
(a) 0 GPa
(b) 200 GPa
Figure: Phonon density of states F(ω) of candidate structures of metastable metallic hydrogen, at two pressures. An inset of the imaginary-frequency region is shown in (a). The dotted line in (b) is used to separate the stable from unstable phonon frequencies.
(The imaginary frequencies in (a) are not anomalous1 .) 6
E. G. Brovman, Yu. Kagan, and A. Kholas, JETP 34, 1300–1315 (1972) Jeffrey M. McMahon (WSU)
Dynamic Stabilities
March 14, 2017
17 / 19
Summary and Open Questions Summary • The Cs-IV structure of atomic metallic hydrogen is metastable, to approximately 250 GPa • Candidate structures for metastable metallic hydrogen were predicted • ... None were dynamically stable (below 250 GPa) • The processes involved with these phases, as the pressure is changed, was analyzed
Open Questions • Properties (e.g., superconductivity) • Quantitative corrections (calculation uncertainties, anharmonic effects, etc.) • Thermal effects and lifetimes (considering zero-point energy) C. M. Tenney, K. L. Sharkey, and J. M. McMahon, In Preparation (2017) Jeffrey M. McMahon (WSU)
Concluding Remarks
March 14, 2017
18 / 19
Acknowledgments Members of the McMahon Research Group on the hydrogen project: This project:
Other projects: • Jeevake Attapattu • Zachary Croft
Craig M. Tenney
Keeper L. Sharkey
Start-up support:
Department of Physics & Astronomy Jeffrey M. McMahon (WSU)
Concluding Remarks
March 14, 2017
19 / 19
Department of Physics & Astronomy
March 14, 2017
Jeffrey M. McMahon (WSU)
March 14, 2017
1 / 19
The Problem of Metallic Hydrogen In 1935, Wigner and Huntington predicted1 that sufficient pressure would dissociate hydrogen molecules; and any Bravais lattice of such atoms would be metallic. Initial interest was primarily related to astrophysical problems ... ... later to fundamental physics; also with practical applications: • It is expected to have remarkable properties: I I I
High-temperature superconductivity2,3 Novel types of quantum fluids3 Powerful rocket fuel
• Theoretical capabilities have progressed4 ; properties can be accurately calculated from first principles
• ... As have experimental ones; it is now possible to create in the laboratory
• It may be metastable6 1
E. Wigner and H. B. Huntington, J. Chem. Phys. 3, 764–770 (1935) N. W. Ashcroft, Phys. Rev. Lett. 21, 1748–1749 (1968) solid 3 E. Babaev, A. Sudbo, and N. W. Ashcroft, Nature 431, 666-668 (2004) 4 J. M. McMahon, M. A. Morales, C. Pierleoni, and D. M. Ceperley, Rev. Mod. Phys. 84, 1607–1653 (2012) 5 R. P. Dias and I. F. Silvera, Science (2017) 6 E. G. Brovman, Yu. Kagan, and A. Kholas, JETP 34, 1300–1315 (1972) 2
Jeffrey M. McMahon (WSU)
Introduction
metallic hydrogen5
March 14, 2017
2 / 19
Metastable Metallic Hydrogen There is little doubt that hydrogen becomes metallic at high pressures. For practical (terrestrial) applications, ... ... the significant, outstanding1,2 question is whether it is metastable. Note: The answer to this presupposes solutions to the following: (1) Determination of the minimum-energy crystal structure(s), ... ... and proof that it lies at a stationary point (2) Proof of (dynamic) stability (3) Analysis of the relation between the ground- and metastable-state structures of metallic hydrogen, and the molecular phase (4) Determination of the lifetime of the metastable state
It is the purpose of this talk to answer this question. 1 2
E. G. Brovman, Yu. Kagan, and A. Kholas, JETP 34, 1300–1315 (1972) H. Y. Geng, H. X. Song, J. F. Li, and Q. Wu, J. Appl. Phys. 111, 063510 (2012) Jeffrey M. McMahon (WSU)
Introduction
March 14, 2017
3 / 19
Metastable Metallic Hydrogen (Continued)
Brief answer:
Metallic hydrogen is metastable ...
... but, only to (approximately) 250 GPa.
Jeffrey M. McMahon (WSU)
Introduction
March 14, 2017
4 / 19
Outline
This talk is outlined as follows: • Phase diagram of hydrogen • Crystal structure(s) of metallic hydrogen I Ground-state structures of atomic metallic hydrogen I Searching for (new) metastable metallic structures • Dynamic stabilities • Concluding remarks
Jeffrey M. McMahon (WSU)
Introduction
March 14, 2017
5 / 19
Phase Diagram of Hydrogen
Temperature (K)
1500
liquid H2 liquid H
1000
500
I
I'? IV V III
II 0
10
VI H2-PRE
100
1000
solid H 10000
Pressure (GPa) Figure: Combined experimental (lines) and theoretical (lines with points) phase diagram of hydrogen. Data compiled from several sources Jeffrey M. McMahon (WSU)
Phase Diagram
March 14, 2017
6 / 19
Ground-state Structures of Atomic Metallic Hydrogen We earlier found1 several candidate structures of atomic metallic hydrogen:
(Not shown are:
• High-pressure molecular phases, since found2 • Ultrahigh-pressure phases, since found2 • Nearly-degenerate, lower-symmetry distortions3 , not originally reported in Ref. [1])
m
m
Figure: Ground-state enthalpies of the crystal structures of atomic metallic hydrogen relative to fcc ∆H, not including proton zero-point energy. The inset shows an expanded view of the ultrahigh-pressure region.
1
J. M. McMahon and D. M. Ceperley, Phys. Rev. Lett. 106, 165302 (2011) H. Liu, H. Wang, and Y. Ma, J. Phys. Chem. C 116, 9221–9226 (2012) 3 H. Y. Geng, H. X. Song, J. F. Li, and Q. Wu, J. Appl. Phys. 111, 063510 (2012)
2
Jeffrey M. McMahon (WSU)
Atomic Metallic Hydrogen
March 14, 2017
7 / 19
Tetragonal Family of Structures The most promising candidates structures for higher-pressure, metastable metallic hydrogen belong to the tetragonal family. Examples:
x y c/a Diamond
β-Sn
Cs-IV
Fd-3m (c/a = 1.414)
I41 /amd (c/a ∼ 0.55 < 1)
I41 /amd (c/a ∼ 3.73 > 1)
(These structures differ only in their c/a ratio.) Jeffrey M. McMahon (WSU)
Tetragonal Structures
March 14, 2017
8 / 19
0 K Phase Diagram of the Tetragonal Structures The 0 K phase diagram of the tetragonal structures reveals general, qualitative features of metastable metallic hydrogen: Regions of distinct qualitative behavior are indicated with dotted lines:
200
∆H (meV/proton)
100
• > 300 GPa: The (relative) enthalpies re-
0 -100
•
-200 c/a > 1
-300 -400 0
100
200
300
400
500
P (GPa)
•
Figure: 0 K phase diagram of the tetragonal structures of metallic hydrogen.
Jeffrey M. McMahon (WSU)
Tetragonal Structures
main relatively flat (with the exception of diamond). 200–300 GPa: The enthalpies of the βSn, diamond, and Cs-IV structures become nearly degenerate; those of c/a 1 and 1 begin to decrease and increase, respectively. < 200 GPa: The enthalpies of the β-Sn, diamond, and Cs-IV structures becomes degenerate; those of c/a 1 and 1 decrease and increase sharply.
March 14, 2017
9 / 19
Energies of the Tetragonal Structures This behavior is understood by considering the energies of the tetragonal structures over all c/a: Notice that:
• At 500 GPa: The diamond structure forms
E (eV/proton)
-13
•
-14
• -15
-16
3
1.15017 Å (500 GPa) 3 1.40361 Å (300 GPa) 3 1.62966 Å (200 GPa) 3 2.08067 Å (100 GPa) 3 4.72262 Å (0 GPa)
0
2
4
6
8
General feature: The collapse of energy barriers, between 200–300 GPa. 10
c/a
Figure: Energies of the tetragonal structures of metallic hydrogen, as a function of c/a. Fixed volumes are shown, calculated as an average of each structure at the pressure indicated in parenthesis. The c/a 1, β-Sn, diamond, Cs-IV, and 1 structures are indicated with arrows (in order of increasing c/a). Jeffrey M. McMahon (WSU)
the energy barrier (a maximum) between the β-Sn and Cs-IV structures. ... Between 300–200 GPa: The energy barrier collapses (appreciably). ... By 100 GPa: The energy barrier has collapsed, leaving a (single) minimum.
• The increase in energy as c/a → 0 and •
Tetragonal Structures
→ ∞ are due to the coming together of Brillouin planes with like charges. ... This explains the sharp decrease in the c/a 1 structure (nearest-neighbor distance: 0.9900Å): General feature: The tendency to form molecules, below 200 GPa. March 14, 2017
10 / 19
Energy of Metallic Hydrogen There is no guarantee that the high-pressure structures are representative of those (potentially) metastable at lower pressures. Consider the energy E (in the adiabatic approximation): E = Estatic + EZPE • Estatic (static energy) favors anisotropic structures (in metallic hydrogen). • EZPE (zero-point energy) favors symmetric structures. • At low pressures: The energy is determined largely by Estatic . • With increasing pressure: The importance of EZPE increases.
We performed searches1 for metallic (atomic and/or mixed atomic/molecular) structures with the lowest Estatic at 0 GPa. 1
C. J. Pickard and R. J. Needs, J. Phys. Condens. Matter 23, 053201 (2011) Jeffrey M. McMahon (WSU)
Energy of Metallic Hydrogen
March 14, 2017
11 / 19
Relative Enthalpies of Candidate Structures Our searches revealed several candidate structures of (potentially) metastable metallic hydrogen. The results are (qualitatively) understood, and further support the prior results, by considering only the most stable structures. Example:
∆H (meV/proton)
100 0
Regions of distinct behavior (dotted lines):
-100
• < 200 GPa: Molecular structures become
-200 -300 Cs-IV 1H 3H 5H 7H 9H
-400 -500 -600 -700
0
100
200
300
400
• 500
P (GPa)
Figure: Relative enthalpies ∆H (to Cs-IV) of the most stable metallic structures, from searches over different numbers of atoms in the unit cell. Jeffrey M. McMahon (WSU)
Candidate Structures
•
favored. Indicated by sharply decreasing enthalpies, with a magnitude proportional to the number of atoms in the search. 200–300 GPa: Transition region. Energy barriers between (atomic) metallic phases form. Indicated by crossovers. > 300 GPa: Energy barriers between metallic phases exist. Structures found at 0 GPa are also no longer representative. Indicated by lack of trend.
March 14, 2017
12 / 19
The Tendency Towards Molecules Below 200 GPa: All structures found show a tendency towards the formation of molecules. Example: As the number of atoms increases, there is a greater proportion of molecules (to atoms) formed:
(a) C2 (5)
(b) C2/m (7)
(c) Pm (9)
Figure: Most stable candidate structures of metastable metallic hydrogen at ambient conditions, for searches with different numbers of atoms (shown in parenthesis). Jeffrey M. McMahon (WSU)
Candidate Structures
March 14, 2017
13 / 19
The Tendency Towards Molecules (Continued) Example: This tendency is seen even in the (completely) atomic structures (found in searches over a lower number of atoms):
(a) Pmmm
(b) Cmmm
(c) Immm
(d) P6/mmm
Figure: Most stable (nearly degenerate) candidate structures of metastable atomic metallic hydrogen at ambient conditions.
(Notice the proximity of atoms; distances are on the order of 0.9914 Å.)
These results suggest that the atomic hydrogen at lower pressures has no region(s) of stability. Jeffrey M. McMahon (WSU)
Candidate Structures
March 14, 2017
14 / 19
Dynamic Stability A crystal is dynamically (mechanically) stable, if it executes a stable oscillatory motion about its equilibrium configuration. By solving the equation of motions for the normal modes of vibration, ... ... this is equivalent to the condition that their frequencies ω satisfy: ωj (q)2 ≥ 0 for all phonon branches j and wavevectors q. (For an imaginary frequency means that the system, subject to a small displacement, will disrupt exponentially with time.)
Dynamic stabilities were determined by calculating the phonon density of states F(ω) for each structure. M. Born and K. Huang, “Dynamical Theory of Crystal Lattices” (1954) Jeffrey M. McMahon (WSU)
Dynamic Stabilities
March 14, 2017
15 / 19
Dynamic Stabilities of the Tetragonal Structures The Cs-IV structure is dynamically stable below molecular dissociation: 0.004 500 GPa 350 GPa 250 GPa
0.008
F(ω) / proton
F(ω) / proton
0.003
0.002
0.001
0
200 GPa 100 GPa 50 GPa 0 GPa
0.006
0.004
0.002
0
500
1000
1500
2000
2500
3000
0 -2000
-1000
0
1000
2000
3000
-1
-1
ω (cm )
ω (cm )
(a) stable region
(b) unstable region
Figure: Phonon density of states F(ω) of Cs-IV, at several pressures. Pressure regions of (a) stability and (b) instability are shown. The dotted line in (b) is used to separate the stable from unstable frequencies. The arrow also in this plot indicates the growing instability, with decreasing pressure.
... but not to ambient conditions. (All other tetragonal structures are dynamically unstable.) Jeffrey M. McMahon (WSU)
Dynamic Stabilities
March 14, 2017
16 / 19
Dynamic Stabilities of the Candidate Structures The candidate structures are dynamically unstable, at all pressures: 0.020 0.0005
Pmmm Cmmm Immm P6/mmm
0.0020
0.015
Cmmm Immm P6/mmm
0.010 0.0002
F(ω) / proton
0.0003
F(ω) / proton
0.0004
0.0015
0.0010
0.0001
0.005 0.0005 0
-6000
-4000
-2000
0
2000
0.000 4000
0.0000 -2000
-1000
-1
0
1000
2000
3000
4000
-1
ω (cm )
ω (cm )
(a) 0 GPa
(b) 200 GPa
Figure: Phonon density of states F(ω) of candidate structures of metastable metallic hydrogen, at two pressures. An inset of the imaginary-frequency region is shown in (a). The dotted line in (b) is used to separate the stable from unstable phonon frequencies.
(The imaginary frequencies in (a) are not anomalous1 .) 6
E. G. Brovman, Yu. Kagan, and A. Kholas, JETP 34, 1300–1315 (1972) Jeffrey M. McMahon (WSU)
Dynamic Stabilities
March 14, 2017
17 / 19
Summary and Open Questions Summary • The Cs-IV structure of atomic metallic hydrogen is metastable, to approximately 250 GPa • Candidate structures for metastable metallic hydrogen were predicted • ... None were dynamically stable (below 250 GPa) • The processes involved with these phases, as the pressure is changed, was analyzed
Open Questions • Properties (e.g., superconductivity) • Quantitative corrections (calculation uncertainties, anharmonic effects, etc.) • Thermal effects and lifetimes (considering zero-point energy) C. M. Tenney, K. L. Sharkey, and J. M. McMahon, In Preparation (2017) Jeffrey M. McMahon (WSU)
Concluding Remarks
March 14, 2017
18 / 19
Acknowledgments Members of the McMahon Research Group on the hydrogen project: This project:
Other projects: • Jeevake Attapattu • Zachary Croft
Craig M. Tenney
Keeper L. Sharkey
Start-up support:
Department of Physics & Astronomy Jeffrey M. McMahon (WSU)
Concluding Remarks
March 14, 2017
19 / 19
