Quantum Wavepacket Ab Initio Molecular Dynamics for Extended Systems

Rare Events Sampling Methods for Quantum and Classical Ab Initio Molecular Dynamics

We provide an approach to sample rare events during classical ab initio molecular dynamics and quantum wavepacket dynamics. For classical AIMD, a set of fictitious degrees of freedom are introduced that may harmonically interact with the electronic and nuclear degrees of freedom to steer the dynamics in a conservative fashion toward energetically forbidden regions. A similar approach when introduced for quantum wavepacket dynamics has the effect of biasing the trajectory of the wavepacket centroid toward the regions of the potential surface that are difficult to sample. The approach is demonstrated for a phenol-amine system, which is a prototypical problem for condensed phase-proton transfer, and for model potentials undergoing wavepacket dynamics. In all cases, the approach yields trajectories that conserve energy while sampling rare events.

Capturing Weak Interactions in Surface Adsorbate Systems at Coupled Cluster Accuracy: A Graph-Theoretic Molecular Fragmentation Approach Improved through Machine Learning

The accurate and efficient study of the interactions of organic matter with the surface of water is critical to a wide range of applications. For example, environmental studies have found that acidic polyfluorinated alkyl substances, especially perfluorooctanoic acid (PFOA), have spread throughout the environment and bioaccumulate into human populations residing near contaminated watersheds, leading to many systemic maladies. Thus, the study of the interactions of PFOA with water surfaces became important for the mitigation of their activity as pollutants and threats to public health. However, theoretical study of the interactions of such organic adsorbates on the surface of water, and their bulk concerted properties, often necessitates the use of ab initio methods to properly incorporate the long-range electronic properties that govern these extended systems. Notable theoretical treatments of "on-water" reactions thus far have employed hybrid DFT and semilocal DFT, but the interactions involved are weak interactions that may be best described using post-Hartree-Fock theory. Here, we aim to demonstrate the utility of a graph-theoretic approach to molecular fragmentation that accurately captures the critical "weak" interactions while maintaining an efficient ab initio treatment of the long-range periodic interactions that underpin the physics of extended systems. We apply this graph-theoretical treatment to study PFOA on the surface of water as a model system for the study of weak interactions seen in the wide range of surface interactions and reactions. The approach divides a system into a set of vertices, that are then connected through edges, faces, and higher order graph theoretic objects known as simplexes, to represent a collection of locally interacting subsystems. These subsystems are then used to construct ab initio molecular dynamics simulations and for computing multidimensional potential energy surfaces. To further improve the computational efficiency of our graph theoretic fragmentation method, we use a recently developed transfer learning protocol to construct the full system potential energy from a family of neural networks each designed to accurately model the behavior of individual simplexes. We use a unique multidimensional clustering algorithm, based on the k-means clustering methodology, to define our training space for each separate simplex. These models are used to extrapolate the energies for molecular dynamics trajectories at PFOA water interfaces, at less than one-tenth the cost as compared to a regular molecular fragmentation-based dynamics calculation with excellent agreement with couple cluster level of full system potential energies.

Born-Oppenheimer Molecular Dynamics Study on Proton Dynamics of Strong Hydrogen Bonds in Aspirin Crystals, with Emphasis on Differences between Two Crystal Forms

In this study, the proton dynamics of hydrogen bonds for two forms of crystalline aspirin was investigated by the Born-Oppenheimer molecular dynamics (BOMD) method. Analysis of the geometrical parameters of hydrogen bonds using BOMD reveals significant differences in hydrogen bonding between the two crystalline forms of aspirin, Form I and Form II. Analysis of the trajectory for Form I shows spontaneous proton transfer in cyclic dimers, which is absent in Form II. Quantization of the O-H stretching modes allows a detailed discussion on the strength of hydrogen-bonding interactions. The focal point of our study is examination of the hydrogen bond characteristics in the crystal structure and clarification of the influence of hydrogen bonding on the presence of the two crystalline forms of aspirin. In the BOMD method, thermal motions were taken into account. Solving the Schrödinger equation for the snapshots of 2D proton potentials, extracted from MD, gives the best agreement with IR spectra. The character of medium-strong hydrogen bonds in Form I of aspirin was compared with that of weaker hydrogen bonds in aspirin Form II. Two proton minima are present in the potential function for the hydrogen bonds in Form I. The band contours, calculated by using one- and two-dimensional O-H quantization, reflect the differences in the hydrogen bond strengths between the two crystalline forms of aspirin, as well as the strong hydrogen bonding in the cyclic dimers of Form I and the medium-strong hydrogen bonding in Form II.

Car-Parrinello simulation of the vibrational spectrum of a medium strong hydrogen bond by two-dimensional quantization of the nuclear motion: application to 2-hydroxy-5-nitrobenzamide

The nature of medium strong intra- and intermolecular hydrogen bonding in 2-hydroxy-5-nitrobenzamide in the crystal phase was examined by infrared spectroscopy and Car-Parrinello molecular dynamics simulation. The focal point of our study was the part of the infrared spectra associated with the O-H and N-H stretching modes that are very sensitive to the strength of hydrogen bonding. For spectra calculations we used an isolated dimer and the fully periodic crystal environment. We calculated the spectra by using harmonic approximation, the time course of the dipole moment function as obtained from the Car-Parrinello simulation, and the quantization of the nuclear motion of the proton for an instantaneous snapshot of the structures in one and two dimensions. Although quantitative assessment of the agreement between the computed and experimental band contour is difficult due to the fact that the experimental band is very broad, we feel that the most reasonable qualitative agreement with the experiment is obtained from snapshot structures and two-dimensional quantization of the proton motion. We have also critically examined the methods of constructing the one-dimensional proton potential. Perspectives are given for the treatment of nuclear quantum effects in biocatalysis.