Prospects for release-node quantum Monte Carlo

Stochastic nodal surfaces in quantum Monte Carlo calculations

Treating the fermionic ground state-problem as a constrained stochastic optimization problem, a formalism for fermionic quantum Monte Carlo is developed that makes no reference to a trial wave function. Exchange symmetry is enforced by nonlocal terms appearing in the Green's function corresponding to an additional walker propagation channel. Complemented by a treatment of diffusion that encourages the formation of a stochastic nodal surface, we find that an approximate long-range extension of walker cancellations can be employed without introducing significant bias, reducing the number of walkers required for a stable calculation. A proof-of-concept implementation is shown to give a stable fermionic ground state for simple harmonic and atomic systems.

Fermion sign problem in imaginary-time projection continuum quantum Monte Carlo with local interaction

We use the shadow wave function formalism as a convenient model to study the fermion sign problem affecting all projector quantum Monte Carlo methods in continuum space. We demonstrate that the efficiency of imaginary-time projection algorithms decays exponentially with increasing number of particles and/or imaginary-time propagation. Moreover, we derive an analytical expression that connects the localization of the system with the magnitude of the sign problem, illustrating this behavior through numerical results. Finally, we discuss the computational complexity of the fermion sign problem and methods for alleviating its severity.

Molecular-Atomic Transition along the Deuterium Hugoniot Curve with Coupled Electron-Ion Monte Carlo Simulations

We have performed simulations of the principal deuterium Hugoniot curve using coupled electron-ion Monte Carlo calculations. Using highly accurate quantum Monte Carlo methods for the electrons, we study the region of maximum compression along the Hugoniot, where the system undergoes a continuous transition from a molecular fluid to a monatomic fluid. We include all relevant physical corrections so that a direct comparison to experiment can be made. Around 50 GPa we find a maximum compression of 4.85. This compression is approximately 5.5% higher than previous theoretical predictions and 15% higher than the most accurate experimental data. Thus first-principles simulations encompassing the most advanced techniques are in disagreement with the results of the best experiments.

A survey on pure sampling in quantum Monte Carlo methods

The most commonly employed diffusion Monte Carlo algorithm and some of its variants afford a way to sample configuration space from a so-called “mixed distribution”, the product of an input trial solution to the Schrödinger equation for the ground state and its unknown exact solution. This mixed distribution is sufficient to compute the ground state energy and other properties represented by operators that commute with the Hamiltonian. These energy-related properties are exact, save for a small bias introduced by the input trial function’s incorrect exchange nodes, the so-called “fixed-node error”. However, properties represented by operators that commute with the position operator are also of interest. When calculated by sampling from the mixed distribution, these properties are much more strongly biased by the input trial function. Our objective is to review methods that allow sampling from the desired “pure” distribution, one that is unbiased except for the exchange node error. Thereby, one accurately calculates physical properties such as the dipole and other electrical moments, electrical response properties of molecules, and particle distribution functions for clusters. We survey the results of calculations that employ pure-sampling methods through what has been published in year 2012. Our review also touches on truly exact sampling methods.