Fragmentation Method for Computing Quantum Mechanics and Molecular Mechanics Gradients for Force Matching: Validation with Hydration Free Energy Predictions Using Adaptive Force Matching.

Abstract

A fragmentation approach referred to as a simple overlapping region method for force matching (SORForM) is presented. SORForM is designed to enable efficient computation of quantum mechanical (QM) forces for large molecules and is validated in the framework of adaptive force matching (AFM) to develop solute models in water. The SORForM method divides a molecule into overlapping QM regions with each region containing a gradient zone and a buffer zone. The buffer zone ensures that the atoms in the gradient zone have their surroundings unchanged with fragmentation. The performance of the method is validated with mefenamic acid and linalyl acetate by comparing the hydration free energies of AFM models developed with and without SORForM. The AFM hydration free energies are also compared with that of the experiments. The models developed with B3LYP-D3(BJ) and def2-TZVP are in excellent agreement with experiments. Our work shows that PBE-D3(BJ) provides less satisfactory results when compared to B3LYP-D3(BJ). The def2-TZVP basis set is found to greatly improve the agreement with experiments when compared to a double-zeta quality basis set.