TurboRVB: A many-body toolkit for ab initio electronic simulations by quantum Monte Carlo.

TurboRVB is a computational package for ab initio Quantum Monte Carlo (QMC) simulations of both molecular and bulk electronic systems. The code implements two types of well established QMC algorithms: Variational Monte Carlo (VMC) and diffusion Monte Carlo in its robust and efficient lattice regularized variant. A key feature of the code is the possibility of using strongly correlated many-body wave functions (WFs), capable of describing several materials with very high accuracy, even when standard mean-field approaches [e.g., density functional theory (DFT)] fail. The electronic WF is obtained by applying a Jastrow factor, which takes into account dynamical correlations, to the most general mean-field ground state, written either as an antisymmetrized geminal power with spin-singlet pairing or as a Pfaffian, including both singlet and triplet correlations. This WF can be viewed as an efficient implementation of the so-called resonating valence bond (RVB) Ansatz, first proposed by Pauling and Anderson in quantum chemistry [L. Pauling, The Nature of the Chemical Bond (Cornell University Press, 1960)] and condensed matter physics [P.W. Anderson, Mat. Res. Bull 8, 153 (1973)], respectively. The RVB Ansatz implemented in TurboRVB has a large variational freedom, including the Jastrow correlated Slater determinant as its simplest, but nontrivial case. Moreover, it has the remarkable advantage of remaining with an affordable computational cost, proportional to the one spent for the evaluation of a single Slater determinant. Therefore, its application to large systems is computationally feasible. The WF is expanded in a localized basis set. Several basis set functions are implemented, such as Gaussian, Slater, and mixed types, with no restriction on the choice of their contraction. The code implements the adjoint algorithmic differentiation that enables a very efficient evaluation of energy derivatives, comprising the ionic forces. Thus, one can perform structural optimizations and molecular dynamics in the canonical NVT ensemble at the VMC level. For the electronic part, a full WF optimization (Jastrow and antisymmetric parts together) is made possible, thanks to state-of-the-art stochastic algorithms for energy minimization. In the optimization procedure, the first guess can be obtained at the mean-field level by a built-in DFT driver. The code has been efficiently parallelized by using a hybrid MPI-OpenMP protocol, which is also an ideal environment for exploiting the computational power of modern Graphics Processing Unit accelerators.

Quantum Monte Carlo Study of the Protonated Water Dimer.

We report an extensive theoretical study of the protonated water dimer H5O2(+) (Zundel ion) by means of the highly correlated variational Monte Carlo and lattice regularized Monte Carlo approaches. This system represents the simplest model for proton transfer (PT), and a correct description of its properties is essential in order to understand the PT mechanism in more complex aqueous systems. Our Jastrow correlated AGP wave function ensures an accurate treatment of electron correlation. By exploiting the advantage of contracting the primitive basis set over atomic hybrid orbitals, we are able to limit dramatically the number of variational parameters with a systematic control on the numerical precision, a crucial ingredient in order to simulate larger systems. For both energetics and geometrical properties, our QMC results are found to be in excellent agreement with state-of-the-art coupled cluster CCSD(T) techniques. A comparison with density functional theory in the PBE approximation points to the crucial role of electron correlation for a correct description of the PT in the dimer. We prove that the QMC framework used in this work is able to resolve the tiny energy differences (∼0.3 kcal/mol) and structural variations involved in proton transfer reactions. Our approach combines these features and a favorable N(4) scaling with the number of particles which paves the way to the simulation of more realistic PT models. A test calculation on a larger protonated water cluster is carried out. The QMC approach used here represents a promising candidate to provide the first high-level ab initio description of PT in water.