Full Temperature-Dependent Potential and Anharmonicity in Metallic Hydrogen: Colossal NQE and the Consequences

The temperature-dependent effective potential (TDEP) method for anharmonic phonon dispersion is generalized to the full potential case by combining with path integral formalism. This extension naturally resolves the intrinsic difficulty in the original TDEP at low temperature. The new method is applied to solid metallic hydrogen at high pressure. A colossal nuclear quantum effect (NQE) and subsequent anharmonicity are discovered, which not only leads to unexpectedly large drift of protons, but also slows down the convergence rate substantially when computing the phonon dispersions. By employing direct ab initio path integral molecular dynamics simulations as the benchmark, a possible breakdown of phonon picture in metallic hydrogen due to colossal NQE is indicated, implying novel lattice dynamical phenomena might exist. Inspired by this observation, a general theoretical formalism for quantum lattice dynamics beyond phonon is sketched, with the main features being discussed.

[1]  H. Geng.  Public Debate on Metallic Hydrogen to Boost High Pressure Research , 2017, 1803.11418.

[2]  H. Geng,et al.  Prediction of a Mobile Solid State in Dense Hydrogen under High Pressures. , 2017, The journal of physical chemistry letters.

[3]  H. Geng,et al.  Predicted reentrant melting of dense hydrogen at ultra-high pressures , 2016, Scientific Reports.

[4]  H. Geng.  Accelerating ab initio path integral molecular dynamics with multilevel sampling of potential surface , 2014, J. Comput. Phys..

[5]  Y. Gohda,et al.  Anharmonic force constants extracted from first-principles molecular dynamics: applications to heat transfer simulations , 2014, Journal of physics. Condensed matter : an Institute of Physics journal.

[6]  H. Geng,et al.  High-pressure behavior of dense hydrogen up to 3.5 TPa from density functional theory calculations , 2010, 1010.3392.

[7]  Dario Alfè,et al.  PHON: A program to calculate phonons using the small displacement method , 2009, Comput. Phys. Commun..

[8]  Kim Lefmann,et al.  Avoided crossing of rattler modes in thermoelectric materials. , 2008, Nature materials.

[9]  H. Sato,et al.  Low-Lying Optical Phonon Modes in the Filled Skutterudite CeRu4Sb12 , 2006, cond-mat/0611038.

[10]  N. Ashcroft,et al.  A superconductor to superfluid phase transition in liquid metallic hydrogen , 2004, Nature.

[11]  G. Kresse,et al.  Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set , 1996 .