A complete OSV-MP2 analytical gradient theory for molecular structure and dynamics simulations

Abstract

We propose an exact algorithm for computing the analytical gradient within the framework of the orbital-specific-virtual (OSV) second-order Møller-Plesset (MP2) theory in resolution-of-identity (RI) approximation. We implement the exact relaxation of perturbed OSVs through the explicit constraints of the perturbed orthonormality, the perturbed diagonality and the perturbed singular value condition. We explicitly show that the rotation of OSVs within the retained OSV subspace makes no contribution to gradients, as long as the iterative solution of the unperturbed Hylleraas residual equation is well converged. The OSV relaxation is solved as the perturbed non-degenerate singular value problem between the retained and discarded OSV subspaces. The detailed derivation and preliminary implementations for gradient working equations are discussed. The coupled-perturbed localization method is implemented for meta-Löwdin localization function. The numerical accuracy of computed OSV-MP2 gradients is demonstrated for the geometries of selected molecules that are often discussed in other theories. Moreover, the OSV-MP2 analytical gradients can generate atomic forces that are utilized to drive the Born-Oppenheimer molecular dynamics (BOMD) simulation for studying structural and vibrational properties with respect to OSV selections. By performing the OSV-MP2 NVE BOMD calculation using the normal OSV selection, the structural and vibrational details of protonated water cations are well reproduced. The 200 picoseconds NVT well-tempered metadynamics at 300 K has been simulated to compute the OSV-MP2 rotational free energy surface of coupled hydroxyl and methyl rotors for ethanol molecule.

References in corpus (4)

• Towards Exact Molecular Dynamics Simulations with Machine-Learned Force Fields
• Tensor factorizations of local second-order Møller Plesset theory
• Molecular Force Fields with Gradient-Domain Machine Learning: Construction and Application to Dynamics of Small Molecules with Coupled Cluster Forces
• Exploration of Reduced Scaling Formulation of Equation of Motion Coupled-Cluster Singles and Doubles Based on State-Averaged Pair Natural Orbitals

Cited by in corpus (3)

• Low-data deep quantum chemical learning for accurate MP2 and coupled-cluster correlations
• Performant Automatic Differentiation of Local Coupled Cluster Theories: Response Properties and Ab Initio Molecular Dynamics
• Third-order many-body expansion of OSV-MP2 wavefunction for low-order scaling analytical gradient computation