Coupled Electron-Ion Monte Carlo simulation of hydrogen molecular crystals

Abstract

We performed simulations for solid molecular hydrogen at high pressures (250GPa$\le$P$\le$500GPa) along two isotherms at T=200 K (phases III and VI) and at T=414 K (phase IV). At T=200K we considered likely candidates for phase III, the C2c and Cmca12 structures, while at T=414K in phase IV we studied the Pc48 structure. We employed both Coupled Electron-Ion Monte Carlo (CEIMC) and Path Integral Molecular Dynamics (PIMD) based on Density Functional Theory (DFT) using the vdW-DF approximation. The comparison between the two methods allows us to address the question of the accuracy of the xc approximation of DFT for thermal and quantum protons without recurring to perturbation theories. In general, we find that atomic and molecular fluctuations in PIMD are larger than in CEIMC which suggests that the potential energy surface from vdW-DF is less structured than the one from Quantum Monte Carlo. We find qualitatively different behaviors for systems prepared in the C2c structure for increasing pressure. Within PIMD the C2c structure is dynamically partially stable for P$\le$250GPa only: it retains the symmetry of the molecular centers but not the molecular orientation; at intermediate pressures it develops layered structures like Pbcn or Ibam and transforms to the metallic Cmca-4 structure at P$\ge$450GPa. Instead, within CEIMC, the C2c structure is found to be dynamically stable at least up to 450GPa; at increasing pressure the molecular bond length increases and the nuclear correlation decreases. For the other two structures the two methods are in qualitative agreement although quantitative differences remain. We discuss various structural properties and the electrical conductivity. We find these structures become conducting around 350GPa but the metallic Drude-like behavior is reached only at around 500GPa, consistent with recent experimental claims.

arXiv comment16 pages, 13 figures, invited article submitted to the special issue on Nuclear Quantum Effects in the Journal of Chemical Physics

References in corpus (23)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• Van der Waals Density Functional for General Geometries
• Observation of the Wigner-Huntington Transition to Solid Metallic Hydrogen
• Accelerating the convergence of path integral dynamics with a generalized Langevin equation
• Nuclear Quantum Effects and Nonlocal Exchange-Correlation Functionals Applied to Liquid Hydrogen at High Pressure
• Density functional theory study of phase IV of solid hydrogen
• A First-order Phase Transition to Metallic Hydrogen
• Dissociation of high-pressure solid molecular hydrogen: Quantum Monte Carlo and anharmonic vibrational study
• Backflow Correlations for the Electron Gas and Metallic Hydrogen
• Quantum Monte Carlo study of the phase diagram of solid molecular hydrogen at extreme pressures
• Benchmark of Exchange-Correlation Functionals for High Pressure Hydrogen using Quantum Monte Carlo
• Towards a Predictive First-Principles Description of Solid Molecular Hydrogen with Density-Functional Theory
• Proton Transfer in Phase IV of Solid Hydrogen and Deuterium
• Bonding, structures, and band gap closure of hydrogen at high pressures
• Nature of the Metallization Transition in Solid Hydrogen
• Electrical conductivity of high-pressure liquid hydrogen by quantum Monte Carlo methods
• Low temperature phase diagram of hydrogen at pressures up to 380 GPa. A possible metallic phase at 360 GPa and 200 K
• Trial wave functions for High-Pressure Metallic Hydrogen
• Classical and quantum ordering of protons in cold solid hydrogen under megabar pressures
• Free energy methods in Coupled Electron Ion Monte Carlo
• New low temperature phase in dense hydrogen: The phase diagram to 421 GPa
• Response to critiques on Observation of the Wigner-Huntington transition to metallic hydrogen
• Energetics of H${}_2$ clusters from density functional and coupled cluster theories

Cited by in corpus (16)

• A Perspective on Conventional High-Temperature Superconductors at High Pressure: Methods and Materials
• Towards high-throughput superconductor discovery via machine learning
• First principles simulations of dense hydrogen
• Energy gap closure of crystalline molecular hydrogen with pressure
• Ab Initio Molecular Dynamics on Quantum Computers
• Local structure in dense hydrogen at the liquid-liquid phase transition by Coupled Electron-Ion Monte Carlo
• Probing anharmonic phonons by quantum correlators: A path integral approach
• Benchmarking vdW-DF first principle predictions against Coupled Electron-Ion Monte Carlo for high pressure liquid hydrogen
• Electronic structure and optical properties of quantum crystals from first principles calculations in the Born-Oppenheimer approximation
• Nuclear magnetic resonance spectroscopy as a dynamical structural probe of high pressure hydrogen
• A structural optimization algorithm with stochastic forces and stresses
• Electron localization properties in high pressure hydrogen at the liquid-liquid phase transition by Coupled Electron-Ion Monte Carlo
• High-pressure II-III phase transition in solid hydrogen: Insights from state-of-the-art ab initio calculations
• Semimetallic molecular hydrogen at pressure above 350 GPa
• Hydrogen phase-IV characterization by full account of quantum anharmonicity
• Experimental phase transition mapping for hydrogen above 300 K up to 300 GPa