Thermodynamics of liquids: standard molar entropies and heat capacities of common solvents from 2PT molecular dynamics

Abstract

We validate here the Two-Phase Thermodynamics (2PT) method for calculating the standard molar entropies and heat capacities of common liquids. In 2PT, the thermodynamics of the system is related to the total density of states (DoS), obtained from the Fourier Transform of the velocity autocorrelation function. For liquids this DoS is partitioned into a diffusional component modeled as diffusion of a hard sphere gas plus a solid component for which the DoS(u) -0 as u -0 as for a Debye solid. Thermodynamic observables are obtained by integrating the DoS with the appropriate weighting functions. In the 2PT method, two parameters are extracted from the DoS self-consistently to describe diffusional contributions: the fraction of diffusional modes, f, and DoS(0). This allows 2PT to be applied consistently and without re-parameterization to simulations of arbitrary liquids. We find that the absolute entropy of the liquid can be determined accurately from a single short MD trajectory (20 ps) after the system is equilibrated, making it orders of magnitude more efficient than commonly used perturbation and umbrella sampling methods. Here, we present the predicted standard molar entropies for fifteen common solvents evaluated from molecular dynamics simulations using the AMBER, GAFF, OPLS AA/L and Dreiding II forcefields. Overall, we find that all forcefields lead to good agreement with experimental and previous theoretical values for the entropy and very good agreement in the heat capacities. These results validate 2PT as a robust and efficient method for evaluating the thermodynamics of liquid phase systems. Indeed 2PT might provide a practical scheme to improve the intermolecular terms in forcefields by comparing directly to thermodynamic properties.

Overstretching of B-DNA with various pulling protocols: Appearance of structural polymorphism and S-DNA

“…To understand the contribution of entropy (and equivalently the enthalpy) to free energy, we employed two-phase thermodynamic (2PT) method developed by Lin et al [74][75][76] and calculate the entropy of DNA during overstretching processes. Fig.…”

Synthesis of single-component metallic glasses by thermal spray of nanodroplets on amorphous substrates

“…To extract the entropy evolution of the nanodroplets during the dynamic impact process, we used the LinBlanco-Goddard two phase thermodynamic (2PT) model 15 to obtain the entropy of the nanodroplets impacting on the substrates. 2PT has been applied to various systems [15][16][17][18] to obtain accurate entropies.…”

Recovery Act: Molecular Simulation of Dissolved Inorganic Carbons for Underground Brine CO2 Sequestration

“…From the associated molecular trajectories of our ReaxFF-DIC simulations, we can access the system thermodynamics using the Two-Phase Thermodynamics method [9][10][11] developed at Caltech, which has been recently applied to studying the stability of hydrocarbons at deep Earth pressures and temperatures [12]. This method allows for a complete description of the entropy and enthalpy effects that determine the stability of the DIC species in water, under varying conditions of brine concentration, temperature and pressure.…”