Direct evaluation of the force constant matrix in quantum Monte Carlo
We develop a formalism to directly evaluate the matrix of force constants within a Quantum Monte Carlo calculation. We utilize the matrix of force constants to accurately relax the positions of atoms in molecules and determine their vibrational modes, using a combination of Variational and Diffusion Monte Carlo. The computed bond lengths differ by less than 0.007Å from the experimental results for all four tested molecules. For hydrogen and hydrogen chloride, we obtain fundamental vibrational frequencies within 0.1% of experimental results and ~10 times more accurate than leading computational methods. For carbon dioxide and methane, the vibrational frequency obtained is on average within 1.1% of the experimental result, which is at least 3 times closer than results using Restricted Hartree-Fock and Density Functional Theory with a Perdew-Burke-Ernzerhof (PBE) functional and comparable or better than Density Functional Theory with a semi-empirical functional.