Chemical accuracy from quantum Monte Carlo for the Benzene Dimer

We report an accurate study of interactions between Benzene molecules using variational quantum Monte Carlo (VMC) and diffusion quantum Monte Carlo (DMC) methods. We compare these results with density functional theory (DFT) using different van der Waals (vdW) functionals. In our QMC calculations, we use accurate correlated trial wave functions including three-body Jastrow factors, and backflow transformations. We consider two benzene molecules in the parallel displaced (PD) geometry, and find that by highly optimizing the wave function and introducing more dynamical correlation into the wave function, we compute the weak chemical binding energy between aromatic rings accurately. We find optimal VMC and DMC binding energies of -2.3(4) and -2.7(3) kcal/mol, respectively. The best estimate of the CCSD(T)/CBS limit is -2.65(2) kcal/mol [E. Miliordos et al, J. Phys. Chem. A 118, 7568 (2014)]. Our results indicate that QMC methods give chemical accuracy for weakly bound van der Waals molecular interactions, comparable to results from the best quantum chemistry methods.Comment: Accepted for publication in the Journal of Chemical Physics, Vol. 143, Issue 11, 201

Resonating Valence Bond Quantum Monte Carlo: Application to the ozone molecule

We study the potential energy surface of the ozone molecule by means of Quantum Monte Carlo simulations based on the resonating valence bond concept. The trial wave function consists of an antisymmetrized geminal power arranged in a single-determinant that is multiplied by a Jastrow correlation factor. Whereas the determinantal part incorporates static correlation effects, the augmented real-space correlation factor accounts for the dynamics electron correlation. The accuracy of this approach is demonstrated by computing the potential energy surface for the ozone molecule in three vibrational states: symmetric, asymmetric and scissoring. We find that the employed wave function provides a detailed description of rather strongly-correlated multi-reference systems, which is in quantitative agreement with experiment.Comment: 5 page, 3 figure

Low-pressure phase diagram of crystalline benzene from quantum Monte Carlo

We study the low-pressure (0 to 10 GPa) phase diagram of crystalline benzene using quantum Monte Carlo (QMC) and density functional theory (DFT) methods. We consider the $Pbca$, $P4_32_12$, and $P2_1/c$ structures as the best candidates for phase I and phase II. We perform diffusion quantum Monte Carlo (DMC) calculations to obtain accurate static phase diagrams as benchmarks for modern van der Waals density functionals. We use density functional perturbation theory to compute phonon contribution in the free-energy calculations. Our DFT enthalpy-pressure phase diagram indicates that the $Pbca$ and $P2_1/c$ structures are the most stable phases within the studied pressure range. The DMC Gibbs free-energy calculations predict that the room temperature $Pbca$ to $P2_1/c$ phase transition occurs at 2.1(1) GPa. This prediction is consistent with available experimental results at room temperature. Our DMC calculations show an estimate of 50.6$\pm$0.5 kJ/mol for crystalline benzene lattice energy

Systematic study of finite-size effects in quantum Monte Carlo calculations of real metallic systems

We present a systematic and comprehensive study of finite-size effects in diffusion quantum Monte Carlo calculations of metals. Several previously introduced schemes for correcting finite-size errors are compared for accuracy and efficiency, and practical improvements are introduced. In particular, we test a simple but efficient method of finite-size correction based on an accurate combination of twist averaging and density functional theory. Our diffusion quantum Monte Carlo results for lithium and aluminum, as examples of metallic systems, demonstrate excellent agreement between all of the approaches considered

A consistent description of the iron dimer spectrum with a correlated single-determinant wave function

We study the iron dimer by using an accurate ansatz for quantum chemical calculations based on a simple variational wave function, defined by a single geminal expanded in molecular orbitals and combined with a real space correlation factor. By means of this approach we predict that, contrary to previous expectations, the neutral ground state is $^7 \Delta$ while the ground state of the anion is $^8 \Sigma_g^-$, hence explaining in a simple way a long standing controversy in the interpretation of the experiments. Moreover, we characterize consistently the states seen in the photoemission spectroscopy by Leopold \emph{et al.}. It is shown that the non-dynamical correlations included in the geminal expansion are relevant to correctly reproduce the energy ordering of the low-lying spin states.Comment: 5 pages, submitted to the Chemical Physics Letter

Dissociation of high-pressure solid molecular hydrogen: Quantum Monte Carlo and anharmonic vibrational study

A theoretical study is reported of the molecular-to-atomic transition in solid hydrogen at high pressure. We use the diffusion quantum Monte Carlo method to calculate the static lattice energies of the competing phases and a density-functional-theory-based vibrational self-consistent field method to calculate anharmonic vibrational properties. We find a small but significant contribution to the vibrational energy from anharmonicity. A transition from the molecular Cmca-12 direct to the atomic I4_1/amd phase is found at 374 GPa. The vibrational contribution lowers the transition pressure by 91 GPa. The dissociation pressure is not very sensitive to the isotopic composition. Our results suggest that quantum melting occurs at finite temperature.Comment: Accepted for publication by Phys. Rev. Let

Nature of the metallization transition in solid hydrogen

We present an accurate study of the static-nucleus electronic energy band gap of solid molecular hydrogen at high pressure. The excitonic and quasiparticle gaps of the C 2 / c , P c , P b c n , and P 6 3 / m structures at pressures of 250, 300, and 350 GPa are calculated using the fixed-node diffusion quantum Monte Carlo (DMC) method. The difference between the mean-field and many-body band gaps at the same density is found to be almost independent of system size and can therefore be applied as a scissor correction to the mean-field gap of an infinite system to obtain an estimate of the many-body gap in the thermodynamic limit. By comparing our static-nucleus DMC energy gaps with available experimental results, we demonstrate the important role played by nuclear quantum effects in the electronic structure of solid hydrogen

Nonthermal solid-solid phase transition in ferromagnetic iron

We posit the existence of a nonthermal phase transition in iron, driven by a loss of ferromagnetic ordering on ultrafast timescales with increasing electron temperature. The transition corresponds to a solid-solid BCC to FCC phase transformation and takes place at an electron temperature of 0.62 eV while the ion lattice remains near room temperature. The BCC structure initially undergoes phonon softening during the magnetic transformation, followed by a solid-solid phase transition to the FCC structure, and a subsequent hardening of phonon modes. We present a detailed physical picture of the process, supported by finite-temperature density functional theory simulations of the phonon dispersion curves, electronic density of states, and thermodynamic free energy

Correlation energy of the spin-polarized electron liquid studied using quantum Monte Carlo simulations

Variational and diffusion quantum Monte Carlo (VMC and DMC) methods with Slater-Jastrowbackflow trial wave functions are used to study the spin-polarized three-dimensional uniform electron fluid. We report ground state VMC and DMC energies in the density range 0.5 ≤ rs ≤ 20. Finite-size errors are corrected using canonical-ensemble twist-averaged boundary conditions and extrapolation of the twist-averaged energy per particle calculated at three system sizes (N =113, 259, and 387) to the thermodynamic limit of infinite system size. The DMC energies in the thermodynamic limit are used to parameterize a local spin density approximation correlation function for inhomogeneous electron systems