Thermophysical properties of simple molecular liquid mixtures: On the limitations of some force fields

CO2 and SO2 Capture by Cryptophane-111 Porous Liquid: Insights from Molecular Dynamics Simulations and Computational Chemistry

A computational study of the encapsulation of a gaseous mixture of CO2 and SO2 in a Type II porous liquid is performed under different conditions. The system is composed of cryptophane-111 molecules dispersed in dichloromethane, and it is described using classic molecular dynamics at atomistic resolution. Gaseous CO2 tends to almost fully occupy cryptophane-111's cavities during the first phases of simulation, and, afterwards, it is surpassed by SO2's tendency for occupation. Calculations are performed at five different temperatures in the range of 273 K-310 K, finding a positive correlation between SO2 adsorption and temperature. An evaluation of the radial distribution function of SO2 and CO2 with respect to cryptophane-111 is employed to quantify the number of captured molecules, and an energy study using Density Functional Theory methods is also performed to evaluate the relative stability of the two gases inside the porous liquid.

Extension of the MolMod Database to Transferable Force Fields

MolMod, a web-based database for classical force fields for molecular simulations of fluids [Mol. Sim. 45, 10 (2019), 806-814], was extended to transferable force fields. Eight transferable force fields, including all-atom and united-atom type force fields, were implemented in the MolMod database: OPLS-UA, OPLS-AA, COMPASS, CHARMM, GROMOS, TraPPE, Potoff, and TAMie. These transferable force fields cover a large variety of chemical substance classes. The system is designed such that new transferable force fields can be readily integrated. A graphical user interface was implemented that enables the construction of molecules. The MolMod database compiles the force field for the specified component and force field type and provides the corresponding data and meta data as well as ready-to-use input files for the molecule for different simulation engines. This helps the user to flexibly choose molecular models and integrate them swiftly in their individual workflows, reducing risks of input errors in molecular simulations.

Comparison of Force Fields for the Prediction of Thermophysical Properties of Long Linear and Branched Alkanes

The prediction of thermophysical properties at extreme conditions is an important application of molecular simulations. The quality of these predictions primarily depends on the quality of the employed force field. In this work, a systematic comparison of classical transferable force fields for the prediction of different thermophysical properties of alkanes at extreme conditions, as they are encountered in tribological applications, was carried out using molecular dynamics simulations. Nine transferable force fields from three different classes were considered (all-atom, united-atom, and coarse-grained force fields). Three linear alkanes (n-decane, n-icosane, and n-triacontane) and two branched alkanes (1-decene trimer and squalane) were studied. Simulations were carried out in a pressure range between 0.1 and 400 MPa at 373.15 K. For each state point, density, viscosity, and self-diffusion coefficient were sampled, and the results were compared to experimental data. The Potoff force field yielded the best results.

Benchmarking of Molecular Dynamics Force Fields for Solid-Liquid and Solid-Solid Phase Transitions in Alkanes

Accurate prediction of alkane phase transitions involving solids is needed to prevent catastrophic pipeline blockages, improve lubrication formulations, smart insulation, and energy storage, as well as bring fundamental understanding to processes such as artificial morphogenesis. However, simulation of these transitions is challenging and therefore often omitted in force field development. Here, we perform a series of benchmarks on seven representative molecular dynamics models (TraPPE, PYS, CHARMM36, L-OPLS, COMPASS, Williams, and the newly optimized Williams 7B), comparing with experimental data for liquid properties, liquid-solid, and solid-solid phase transitions of two prototypical alkanes, n-pentadecane (C₁₅) and n-hexadecane (C₁₆). We find that existing models overestimate the melting points by up to 34 K, with PYS and Williams 7B yielding the most accurate results deviating only 2 and 3 K from the experiment. We specially design order parameters to identify crystal-rotator phase transitions in alkanes. United-atom models could only produce a rotator phase with complete rotational disorder, whereas all-atom models using a 12-6 Lennard-Jones potential show no rotator phase even when superheated above the melting point. In contrast, Williams (Buckingham potential) and COMPASS (9-6 Lennard-Jones) reproduce the crystal-to-rotator phase transition, with the optimized Williams 7B model having the most accurate crystal-rotator transition temperature of C₁₅.

Liquid dibromomethane under pressure: a computational study

Molecular dynamics simulations have been performed on liquid dibromomethane at thermodynamic states corresponding to temperature in the range 268-328 K and pressure varying from 1 bar to 3000 bar. The interaction model is a simple effective two-body pair potential with atom-atom Coulomb and Lennard-Jones interactions and molecules are rigid. Thermodynamic properties have been studied, including the isobaric thermal expansion coefficient, the isothermal compressibility, the heat capacities and the speed of sound. The simulation results exhibit a crossing of the isotherms of the isobaric thermal expansion coefficient at about 800 bar in very good agreement with the prediction of an isothermal fluctuation equation of state predicting such a crossing in the pressure range 650-900 bar, though experimental results up to 1000 bar do not find any crossing.