Geometry optimization using Quantum Monte Carlo

Daniel Nissenbaum, B. Barbiellini, Arun Bansil

Monika Dash, Claudia Filippi

Vladimir Iglovikov

Charles A. Schuetz, Michael Frenklach, Alexander C. Kollias +1 more

The solution mapping methodology is presented as a method of geometry optimization of molecules in the quantum Monte Carlo method. Applied to formaldehyde as a test system, this approach is found to yield optimized bond lengths and bond angle in the diffusion Monte Carlo method that lie within experimental error. The variational Monte Carlo optimized geometry also lies within experimental error, with the exception of the CH bond length, which is slightly underestimated. Additionally, the resulting quadratic representation of the potential energy surface in the region of the minimum is used to calculate three of the force constants and harmonic frequencies.

Shih‐I Lu

This article describes an approach in determination of equilibrium geometries and harmonic frequencies of molecules by the Ornstein-Uhlenbeck diffusion quantum Monte Carlo method based on the floating spherical Gaussians. In conjunction with a projected and renormalized Hellmann-Feynman gradient and an electronic energy at variational Monte Carlo and diffusion quantum Monte Carlo, respectively, the quasi-Newton algorithm implemented with the Broyden-Fletcher-Goldfarb-Shanno updated Hessian was used to find the optimized molecular geometry. We applied this approach to N2 and H2O molecules. The geometry and harmonic frequencies calculated were consistent with some sophisticated ab initio calculated values within reasonable statistical uncertainty.

Sam Azadi, R. E. Cohen

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 using different van der Waals functionals. In our quantum Monte Carlo (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 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 coupled-cluster theory through perturbative triplets/complete basis set limit is -2.65(2) kcal/mol [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.