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

Disentangling structural information from core-level excitation spectra

Core-level spectra of liquids can be difficult to interpret due to the presence of a range of local environments. We present computational methods for investigating core-level spectra based on the idea that both local structural parameters and the x-ray spectra behave as functions of the local atomic configuration around the absorbing site. We identify correlations between structural parameters and spectral intensities in defined regions of interest, using the oxygen K-edge excitation spectrum of liquid water as a test case. Our results show that this kind of analysis can find the main structure-spectral relationships of ice, liquid water, and supercritical water.

Enhancing the catalytic activity of hydronium ions through constrained environments

The dehydration of alcohols is involved in many organic conversions but has to overcome high free-energy barriers in water. Here we demonstrate that hydronium ions confined in the nanopores of zeolite HBEA catalyse aqueous phase dehydration of cyclohexanol at a rate significantly higher than hydronium ions in water. This rate enhancement is not related to a shift in mechanism; for both cases, the dehydration of cyclohexanol occurs via an E1 mechanism with the cleavage of C_β–H bond being rate determining. The higher activity of hydronium ions in zeolites is caused by the enhanced association between the hydronium ion and the alcohol, as well as a higher intrinsic rate constant in the constrained environments compared with water. The higher rate constant is caused by a greater entropy of activation rather than a lower enthalpy of activation. These insights should allow us to understand and predict similar processes in confined spaces.

Alcohol dehydration can be challenging in aqueous phase. Here the authors show that hydronium ions confined with zeolite pores catalyse alcohol dehydration at a significantly increased rate relative to aqueous phase hydronium ions, driven by an increased association between the ion and alcohol and a greater entropy of activation.

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.