The Coupled Electronic-Ionic Monte Carlo Simulation Method

Equation of state of metallic hydrogen from coupled electron-ion Monte Carlo simulations

We present a study of hydrogen at pressures higher than molecular dissociation using the coupled electron-ion Monte Carlo method. These calculations use the accurate reptation quantum Monte Carlo method to estimate the electronic energy and pressure while doing a Monte Carlo simulation of the protons. In addition to presenting simulation results for the equation of state over a large region of the phase diagram, we report the free energy obtained by thermodynamic integration. We find very good agreement with density-functional theory based molecular-dynamics calculations for pressures beyond 600 GPa and densities above rho=1.4 g/cm(3) , both for thermodynamic and structural properties. This agreement provides a strong support to the different approximations employed in the density-functional treatment of the system, specifically the approximate exchange-correlation potential and the use of pseudopotentials for the range of densities considered. We find disagreement with chemical models, which suggests that a reinvestigation of planetary models--previously constructed using the Saumon-Chabrier-Van Horn equations of state--might be needed.

Backflow correlations for the electron gas and metallic hydrogen

We justify and evaluate backflow three-body wave functions for a two-component system of electrons and protons. Based on the generalized Feynman-Kacs formula, many-body perturbation theory, and band structure calculations, we analyze the use and the analytical form of the backflow function from different points of view. The resulting wave functions are used in variational and diffusion Monte Carlo calculations of the electron gas and of solid and liquid metallic hydrogen. For the electron gas, the purely analytic backflow and three-body form gives lower energies than those of previous calculations. For bcc hydrogen, analytical and optimized backflow-three-body wave functions lead to energies nearly as low as those from using local density approximation orbitals in the trial wave function. However, compared to wave functions constructed from density functional solutions, backflow wave functions have the advantage of only few parameters to estimate, the ability to include easily and accurately electron-electron correlations, and that they can be directly generalized from the crystal to a disordered liquid of protons.