Analysis of the basis set superposition error in molecular dynamics of hydrogen-bonded liquids: application to methanol.

An efficient protocol is presented to compensate for the basis set superposition error (BSSE) in DFT molecular dynamics (MD) simulations using localized Gaussian basis sets. We propose a classical correction term that can be added a posteriori to account for BSSE. It is tested to what extension this term will improve radial distribution functions (RDFs). The proposed term is pairwise between certain atoms in different molecules and was calibrated by fitting reference BSSE data points computed with the counterpoise method. It is verified that the proposed exponential decaying functional form of the model is valid. This work focuses on hydrogen-bonded liquids, i.e., methanol, and more specific on the intermolecular hydrogen bond, but in principle the method is generally applicable on any type of interaction where BSSE is significant. We evaluated the relative importance of the Grimme-dispersion versus BSSE and found that they are of the same order of magnitude, but with an opposite sign. Upon introduction of the correction, the relevant RDFs, obtained from MD, have amplitudes equal to experiment.

Critical analysis of the accuracy of models predicting or extracting liquid structure information.

This work aims at a critical assessment of properties predicting or extracting information on the density and structure of liquids. State-of-the-art NVT and NpT molecular dynamics (MD) simulations have been performed on five liquids: methanol, chloroform, acetonitrile, tetrahydrofuran, and ethanol. These simulations allow the computation of properties based on first principles, including the equilibrium density and radial distribution functions (RDFs), characterizing the liquid structure. Refinements have been incorporated in the MD simulations by taking into account basis set superposition errors (BSSE). An extended BSSE model for an instantaneous evaluation of the BSSE corrections has been proposed, and their impact on the liquid properties has been assessed. If available, the theoretical RDFs have been compared with the experimentally derived RDFs. For some liquids, significant discrepancies have been observed, and a profound but critical investigation is presented to unravel the origin of these deficiencies. This discussion is focused on tetrahydrofuran where the experiment reveals some prominent peaks completely missing in any MD simulation. Experiments providing information on liquid structure consist mainly of neutron diffraction measurements offering total structure factors as the primary observables. The splitting of these factors in reciprocal space into intra- and intermolecular contributions is extensively discussed, together with their sensitivity in reproducing correct RDFs in coordinate space.

Atomic Velocity Projection Method: A New Analysis Method for Vibrational Spectra in Terms of Internal Coordinates for a Better Understanding of Zeolite Nanogrowth.

An efficient protocol is presented to identify signals in vibrational spectra of silica oligomers based on theoretical molecular dynamics (MD) simulations. The method is based on the projection of the atomic velocity vectors on the tangential directions of the trajectories belonging to a predefined set of internal coordinates. In this way only contributions of atomic motions along these internal coordinates are taken into consideration. The new methodology is applied to the spectra of oligomers and rings, which play an important role in zeolite synthesis. A suitable selection of the relevant internal coordinates makes the protocol very efficient but relies on intuition and theoretical insight. The simulation data necessary to compute vibrational spectra of relevant silica species are obtained through MD using proper force fields. The new methodology-the so-called velocity projection method-makes a detailed analysis of vibrational spectra possible by establishing a one-to-one correspondence between a spectral signal and a proper internal coordinate. It offers valuable perspectives in understanding the elementary steps in silica organization during zeolite nanogrowth. The so-called velocity projection method is generally applicable on data obtained from all types of MD and is a highly valuable alternative to normal-mode analysis which has its limitations due to the presence of many local minima on the potential energy surface. In this work the method is exclusively applied to inelastic neutron scattering, but extension to the infrared power spectrum is apparent.

MD-tracks: a productive solution for the advanced analysis of molecular dynamics and Monte Carlo simulations.

In this paper, we present MD-TRACKS, an advanced statistical analysis toolkit for Molecular Dynamics and Monte Carlo simulations. The program is compatible with different molecular simulation codes, and the analysis results can be loaded into spreadsheet software and plotting tools. The analysis is performed with commands that operate on a binary trajectory database. These commands process not only plain trajectory data but also the output of other MD-TRACKS commands, which enables complex analysis work flows that are easily programmed in shell scripts. The applicability, capabilities, and ease of use of MD-TRACKS are illustrated by means of examples, that is, the construction of vibrational spectra and radial distribution functions from a molecular dynamics run is discussed in the case of tetrahydrofuran. These properties are compared with the experimental data available in the literature. MD-TRACKS is open-source software distributed at http://molmod.ugent.be/code/ .