Transferable Anisotropic United-Atom Force Field Based on the Mie Potential for Phase Equilibrium Calculations: n-Alkanes and n-Olefins

TAMie Force Field for Alkanethiols: Multifidelity Gaussian Processes for Dealing with Scarce Experimental Data

This study extends the transferable anisotropic Mie potential (TAMie) to alkanethiols. The force field parameters are optimized by using an analytic equation of state as a surrogate model. Given the lack of experimental density data at elevated temperatures where Monte Carlo simulations have high statistical precision, the equation of state is supplemented by a linear multifidelity Gaussian process approach to bridge the temperature gap. Force field parameters are adjusted by minimizing squared deviations of calculated vapor pressures and liquid densities from experimental data of 1-propanethiol, 1-butanethiol and 1-pentanethiol leading to small mean absolute relative deviations in liquid densities and vapor pressures. The force field is transferable to higher 1-thiols, as shown for 1-hexanethiol and 1-octanethiol. Individual parameter sets are provided for methanethiol and ethanethiol. The shear viscosity of pure substances is predicted in fair agreement with experimental data, considering that it is not included in the parametrization. Further, the phase behavior of binary mixtures of alkanethiols with alkanes is studied, and predictions of the TAMie model are found in excellent agreement with experimental data.

Transferable Anisotropic Mie Potential Force Field for Alkanediols

The development of force fields for polyfunctional molecules, such as alkanediols, requires a careful account of different average intramolecular conformations for gas states compared to dense liquid states, where intra- and intermolecular hydrogen bonds compete. In the present work, the transferable anisotropic Mie (TAMie) potential is extended to 1,n-alkanediols. Using the convention that intramolecular nonbonded interactions up to and including the third neighbor are excluded, all force field parameters developed previously for 1-alcohols were transferred to 1,5-pentanediol and beyond, with good agreement with experimental phase equilibrium data. To obtain trans-gauche ratios of 1,2-ethanediol and 1,3-propanediol that are consistent with experimental results, the propensities for intra- and intermolecular hydrogen bonds had to be balanced. This was achieved by parameterizing the intramolecular dihedral energy functions governing the O-C-C-O and O-C-C-C angles while intramolecular charge-charge interactions were active. All partial charges belonging to a functional group are collected in a charge group and all interactions among two charge groups are evaluated even if they are separated by less than three bonds. With this approach, it is possible to apply the nonbonded parameters from 1-alcohols to alkanediols without further refinement. The agreement with experimental phase equilibrium and shear viscosity data is of similar quality as for the 1-alcohols and the trans-gauche ratio agrees with literature results from spectroscopic measurements and ab initio calculations.

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.

Molecular modelling of the thermophysical properties of fluids: expectations, limitations, gaps and opportunities

This manuscript provides an overview of the current state of the art in terms of the molecular modelling of the thermophysical properties of fluids. It is intended to manage the expectations and serve as guidance to practising physical chemists, chemical physicists and engineers in terms of the scope and accuracy of the more commonly available intermolecular potentials along with the peculiarities of the software and methods employed in molecular simulations while providing insights on the gaps and opportunities available in this field. The discussion is focused around case studies which showcase both the precision and the limitations of frequently used workflows.

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.

Comparison of the United- and All-Atom Representations of (Halo)alkanes Based on Two Condensed-Phase Force Fields Optimized against the Same Experimental Data Set

The level of accuracy that can be achieved by a force field is influenced by choices made in the interaction-function representation and in the relevant simulation parameters. These choices, referred to here as functional-form variants (FFVs), include for example the model resolution, the charge-derivation procedure, the van der Waals combination rules, the cutoff distance, and the treatment of the long-range interactions. Ideally, assessing the effect of a given FFV on the intrinsic accuracy of the force-field representation requires that only the specific FFV is changed and that this change is performed at an optimal level of parametrization, a requirement that may prove extremely challenging to achieve in practice. Here, we present a first attempt at such a comparison for one specific FFV, namely the choice of a united-atom (UA) versus an all-atom (AA) resolution in a force field for saturated acyclic (halo)alkanes. Two force-field versions (UA vs AA) are optimized in an automated way using the CombiFF approach against 961 experimental values for the pure-liquid densities ρ_liq and vaporization enthalpies ΔH_vap of 591 compounds. For the AA force field, the torsional and third-neighbor Lennard-Jones parameters are also refined based on quantum-mechanical rotational-energy profiles. The comparison between the UA and AA resolutions is also extended to properties that have not been included as parameterization targets, namely the surface-tension coefficient γ, the isothermal compressibility κ_T, the isobaric thermal-expansion coefficient α_P, the isobaric heat capacity c_P, the static relative dielectric permittivity ϵ, the self-diffusion coefficient D, the shear viscosity η, the hydration free energy ΔG_wat, and the free energy of solvation ΔG_che in cyclohexane. For the target properties ρ_liq and ΔH_vap, the UA and AA resolutions reach very similar levels of accuracy after optimization. For the nine other properties, the AA representation leads to more accurate results in terms of η; comparably accurate results in terms of γ, κ_T, α_P, ϵ, D, and ΔG_che; and less accurate results in terms of c_P and ΔG_wat. This work also represents a first step toward the calibration of a GROMOS-compatible force field at the AA resolution.

Equation of state and Helmholtz energy functional for fused heterosegmented hard chains

Modern equations of state for real nonspherical molecules are often based on Wertheim's first-order thermodynamic perturbation theory (TPT1). A major drawback of TPT1 is that it assumes tangentially bonded spheres. In this work, we develop a Helmholtz energy functional for systems comprising hard heterosegmented chains with arbitrary bond lengths. This is achieved by using hard-sphere fragments (i.e., hard spheres with spherical caps removed at the intersection to their neighbors) as monomers as opposed to full hard spheres. The model is written as a Helmholtz energy functional for inhomogeneous systems and the equation of state for a homogeneous system is determined as a special case. We thereby obtain an equation of state that can be used as a reference to develop statistical associating fluid theory models that more accurately describe the thermodynamic properties of nonspherical molecules. The model is validated against molecular simulation results of bulk pressures and density profiles in slit pores. For the bulk pressures, we show that the equation of state is in excellent agreement with results from molecular simulation for dimers, trimers, and chains of up to 20 segments. The density profiles of individual segments of the chains are regarded in slit pores. Some deviations of the theory from results of molecular simulations are observed for strongly fused chains. Overall, however, good agreement is found for inhomogeneous systems.

Implanting Polypyrrole in Metal-Porphyrin MOFs: Enhanced Electrocatalytic Performance for CO2RR

Metal-organic frameworks (MOFs) with plenty of active sites and high porosity have been considered as an excellent platform for the electroreduction of CO₂, yet they are still restricted by the low conductivity or low efficiency. Herein, we insert the electron-conductive polypyrrole (PPy) molecule into the channel of MOFs through the in situ polymerization of pyrrole in the pore of MOF-545-Co to increase the electron-transfer ability of MOF-545-Co and the obtained hybrid materials present excellent electrocatalytic CO₂RR performance. For example, FE_CO of PPy@MOF-545-Co can reach up to 98% at -0.8 V, almost 2 times higher than that of bare MOF-545-Co. The high performance might be attributed to the incorporation of PPy that can serve as electric cables in the channel of MOF to facilitate electron transfer during the CO₂RR process. This attempt might provide new insights to improve the electrocatalytic performance of MOFs for CO₂RR.

Molecular dynamics simulation of octacosane for phase diagrams and properties via the united-atom scheme

We used the united-atom scheme to build three types of crystalline structures for octacosane (C₂₈H₅₈) and carried out molecular dynamics simulations to investigate their phase properties. By gradually heating the three polymorphs, we managed to reproduce the sequence of experimentally reported crystalline phases and rotator phases. By studying the system density, molecule morphology, chain tilt angle and cell anisotropy, we hypothesized three mechanisms behind the observed system deformations and phase transformations during the annealing process. Furthermore, our model successfully predicted the melting temperature and heat of fusion. We also reproduced the characteristics of the rotator phases and the liquid phase, validating the transferability of the united-atom scheme among the different condensed phases of octacosane. Our methodology represents an effective and efficient means of numerical study for octacosane and may be used for other members of the n-alkane family.

Effect of the range of particle cohesion on the phase behavior and thermodynamic properties of fluids

Molecular simulations are performed for the (m + 1, m) potential to systematically investigate the effect of changing the range of particle cohesion on both vapor-liquid equilibria and thermodynamic properties of fluids. The results are reported for m = 4-11, which represent a progressive narrowing of the potential energy well. The conventional Lennard-Jones potential is used as a reference point for normal fluid behavior. Small values of m result in a broadening of the phase envelope compared with the Lennard-Jones potential, whereas a contraction is observed in other cases. The critical properties are reported, and a relationship between the critical temperature and the Boyle temperature is determined. The low values of the critical compressibility factor when m < 6 reflect the behavior observed for real fluids such as n-alkanes. The results for supercritical thermodynamic properties are much more varied. Properties such as pressure, potential energy, isochoric thermal pressure coefficient, and thermal expansion coefficient vary consistently with m, whereas other properties such as the Joule-Thomson coefficient exhibit much more nuanced behavior. Maximum and minimum values are reported for both the isochoric heat capacity and isothermal compressibility. A minimum in the speed of sound is also observed.

Prediction of Adsorption Isotherms and Selectivities: Comparison between Classical Density Functional Theory Based on the Perturbed-Chain Statistical Associating Fluid Theory Equation of State and Ideal Adsorbed Solution Theory

This study gives an assessment of the predictive capability of classical density functional theory (DFT) for adsorption processes of pure substances and mixtures of spherical and nonspherical molecular species. A Helmholtz energy functional based on the perturbed-chain statistical associating fluid theory (PC-SAFT) is applied to calculate isotherms and selectivities of multicomponent adsorption. In order to unambiguously assess the accuracy of the DFT model, we conduct molecular simulations. Monte Carlo (MC) simulations are performed in the grand canonical ensemble using the transition matrix. Two types of systems are studied: a model system, where fluid-fluid and solid-fluid interactions are defined as (single-site) Lennard-Jones interactions, and a more realistic methane-n-butane mixture in a graphite-like pore. Differences between a slit-shaped and a cylindrical pore geometry are examined for the model system. Adsorption isotherms and selectivities obtained from DFT calculations and MC simulations are found in very good agreement, particularly at high pressures. Capillary condensation observed along adsorption isotherms containing n-butane was accurately predicted, both, in equilibrium pressure and in density-increase. Comparisons with results from the ideal adsorbed solution theory are presented, confirming powerful predictions of the DFT approach.

A Modified Shifted Force Approach to the Wolf Summation

The Wolf method for calculation of electrostatic interactions in molecular simulations is known to describe the energy well, whereas the forces have discontinuities. For a more reliable description of the forces this method can be extended with a shifted force approach. This leads to a good description of the forces and precise molecular dynamics simulation, but the description of the energy becomes poorer. In this study we propose a modification of a shifted force extension to describe the energy as well as the forces in better agreement to reference data as determined from the Ewald summation. We show that vapor-liquid phase equilibria (VLE) calculated with Monte Carlo simulations in the grand canonical ensemble and dynamic properties calculated with molecular dynamics simulations can be calculated reliably using this modification to describe the electrostatic interactions.

Prediction of Contact Angles and Density Profiles of Sessile Droplets Using Classical Density Functional Theory Based on the PCP-SAFT Equation of State

This study demonstrates the capability of the density functional theory (DFT) formalism to predict contact angles and density profiles of model fluids and of real substances in good quantitative agreement with molecular simulations and experimental data. The DFT problem is written in cylindrical coordinates, and the solid-fluid interactions are defined as external potentials toward the fluid phase. Monte Carlo (MC) molecular simulations are conducted in order to assess the density profiles resulting from the Helmholtz energy functional used in the DFT formalism. Good quantitative agreement between DFT predictions and MC results for Lennard-Jones and ethane nanodroplets is observed, both for density profiles and for contact angles. That comparison suggests, first, that the Helmholtz energy functional proposed in a previous study [ Sauer , E. ; Gross , J. Ind. Eng. Chem. Res. 56 , 2017 , 4119 - 4135 ] is suitable for three-phase contact lines and, second, that Lagrange multipliers can be used to constrain the number of molecules, similar to a canonical ensemble. Experiments of sessile droplets on solid surfaces are performed to assess whether a real solid with its microscopic roughness can be described through a simple model potential. Comparison of DFT results to experimental data is done for a Teflon surface because Teflon can be regarded as a substrate exhibiting only attractive interactions of van der Waals type. It is shown that the real solid can be described as a perfectly planar solid with effective solvent-to-solid interactions, defined through a single adjustable parameter for the solid. Subsequent predictions for the contact angle of eight solvents, including polar components such as water, are found in very good agreement to experimental data using simple Berthelot-Lorentz combining rules. For the eight investigated solvents, we find mean absolute deviations of 3.77°.

Uncertainty quantification confirms unreliable extrapolation toward high pressures for united-atom Mie λ-6 force field

Molecular simulation results at extreme temperatures and pressures can supplement experimental data when developing fundamental equations of state. Since most force fields are optimized to agree with vapor-liquid equilibria (VLE) properties, however, the reliability of the molecular simulation results depends on the validity/transferability of the force field at higher temperatures and pressures. As demonstrated in this study, although state-of-the-art united-atom Mie λ-6 potentials for normal and branched alkanes provide accurate estimates for VLE, they tend to over-predict pressures for dense supercritical fluids and compressed liquids. The physical explanation for this observation is that the repulsive barrier is too steep for the "optimal" united-atom Mie λ-6 potential parameterized with VLE properties. Bayesian inference confirms that no feasible combination of non-bonded parameters (ϵ, σ, and λ) is capable of simultaneously predicting saturated vapor pressures, saturated liquid densities, and pressures at high temperatures and densities. This conclusion has both practical and theoretical ramifications, as more realistic non-bonded potentials may be required for accurate extrapolation to high pressures of industrial interest.