Related papers: Ab-initio molecular dynamics for high-pressure liq…
The article is devoted to numerical studies of atomic (metal) hydrogen with Path Integral Monte Carlo (PIMC) technique. The research is focused on the range of temperatures and densities where quantum statistics effects are crucial for…
New hybrid Molecular Dynamics-Monte Carlo methods are proposed to increase the efficiency of constant-pressure simulations. Two variations of the isobaric Molecular Dynamics component of the algorithms are considered. In the first, we use…
We present a set of new numerical methods that are relevant to calculating radiation pressure terms in hydrodynamics calculations, with a particular focus on massive star formation. The radiation force is determined from a Monte Carlo…
We construct an effective low-energy Hamiltonian from the classical action via Monte Carlo with importance sampling. We use Monte Carlo (i) to compute matrix elements of the transition amplitude and (ii) to construct stochastically a basis.…
We introduce a general technique to compute finite temperature electronic properties by a novel covariant formulation of the electronic partition function. By using a rigorous variational upper bound to the free energy we are led to the…
We present a new Monte Carlo method which couples Path Integral for finite temperature protons with Quantum Monte Carlo for ground state electrons, and we apply it to metallic hydrogen for pressures beyond molecular dissociation. This…
Identifying the atomic structure and properties of solid hydrogen under high pressures is a long-standing problem of high-pressure physics with far-reaching significance in planetary and materials science. Determining the…
We describe an efficient algorithm to compute forces in quantum Monte Carlo using adjoint algorithmic differentiation. This allows us to apply the space warp coordinate transformation in differential form, and compute all the 3M force…
The static and dynamical properties of heavy water have been studied at ambient conditions with extensive Car-Parrinello molecular-dynamics simulations in the canonical ensemble, with temperatures ranging between 325 K and 400 K.…
Fermionic path integral Monte Carlo simulations have been applied to study the equilibrium properties of the hydrogen and deuterium in the density and temperature range of 1.6 < rs < 14.0 and 5000K < T < 167000K. We use this technique to…
We present extensive molecular dynamics (MD) simulations investigating numerous candidate crystal structures for hydrogen in conditions around the present experimental frontier (400GPa). Spontaneous phase transitions in the simulations…
A Monte Carlo algorithm for computing quantum mechanical expectation values of coordinate operators in many body problems is presented. The algorithm, that relies on the forward walking method, fits naturally in a Green's Function Monte…
We present a technique for extracting Raman intensities from ab initio molecular dynamics (MD) simulations at high temperature. The method is applied to the highly anharmonic case of dense hydrogen up to 500 K for pressures ranging from 180…
The metallic regime of the hydrogen plasma is studied by ab initio Molecular Dynamics simulations, for classical protons and fully degenerate electrons, in the strong coupling regime of the protons. The breakdown of linear screening…
First-principle modeling of dense hydrogen is crucial in materials and planetary sciences. Despite its apparent simplicity, predicting the ionic and electronic structure of hydrogen is a formidable challenge, and it is connected with the…
We present simple and practical strategies to reduce the variance of Monte Carlo estimators. Our focus is on variational Monte Carlo calculations of atomic forces and pressure in electronic systems, although we show that the underlying…
We have mapped the molecular-atomic transition in liquid hydrogen using first principles molecular dynamics. We predict that a molecular phase with short-range orientational order exists at pressures above 100 GPa. The presence of this…
In this work we develop tools that enable the study of non-adiabatic effects with variational and diffusion Monte Carlo methods. We introduce a highly accurate wave function ansatz for electron-ion systems that can involve a combination of…
We study a partially ionized hydrogen plasma by means of quantum molecular dynamics, which is based on wave packets. We introduce a new model which distinguishes between free and bound electrons. The free electrons are modelled as Gaussian…
We develop a formalism to accurately account for the renormalization of electronic structure due to quantum and thermal nuclear motions within the Born-Oppenheimer approximation. We focus on the fundamental energy gap obtained from…