Energy benchmarks for methane-water systems from quantum Monte Carlo and second-order Møller-Plesset calculations

The quantum Monte Carlo (QMC) method is used to generate accurate energy benchmarks for methane-water clusters containing a single methane monomer and up to 20 water monomers. The benchmarks for each type of cluster are computed for a set of geometries drawn from molecular dynamics simulations. The...

Full description

Bibliographic Details
Published in:The Journal of Chemical Physics
Main Authors: Gillan, M. J., Alfè, D., Manby, F. R.
Language:unknown
Published: 2023
Subjects:
Online Access:http://www.osti.gov/servlets/purl/1565342
https://www.osti.gov/biblio/1565342
https://doi.org/10.1063/1.4926444
Description
Summary:The quantum Monte Carlo (QMC) method is used to generate accurate energy benchmarks for methane-water clusters containing a single methane monomer and up to 20 water monomers. The benchmarks for each type of cluster are computed for a set of geometries drawn from molecular dynamics simulations. The accuracy of QMC is anticipated to be comparable with that of coupled-cluster calculations, and this is confirmed by comparisons for the CH 4 -H 2 O dimer. The benchmarks are used to assess the accuracy of the second-order Møller-Plesset (MP2) approximation close to the complete basis-set limit. A recently created embedded many-body technique is shown to give an efficient procedure for computing basis-set converged MP2 energies for the large clusters. It is found that MP2 values for the methane binding energies and the cohesive energies of the water clusters without methane are in close agreement with the QMC benchmarks, but the agreement is aided by partial cancelation between 2-body and beyond-2-body errors of MP2. The embedding approach allows MP2 to be applied without loss of accuracy to the methane hydrate crystal, and it is shown that the resulting methane binding energy and the cohesive energy of the water lattice agree almost exactly with recently reported QMC values