Predicting structural properties of fluids by thermodynamic extrapolation

A structure-dynamics relationship enables prediction of the water hydrogen bond exchange activation energy from experimental data

It has long been understood that the structural features of water are determined by hydrogen bonding (H-bonding) and that the exchange of, or "jumps" between, H-bond partners underlies many of the dynamical processes in water. Despite the importance of H-bond exchanges there is, as yet, no direct method for experimentally measuring the timescale of the process or its associated activation energy. Here, we identify and exploit relationships between water's structural and dynamical properties that provide an indirect route for determining the H-bond exchange activation energy from experimental data. Specifically, we show that the enthalpy and entropy determining the radial distribution function in liquid water are linearly correlated with the activation energies for H-bond jumps, OH reorientation, and diffusion. Using temperature-dependent measurements of the radial distribution function from the literature, we demonstrate how these correlations allow us to infer the value of the jump activation energy, Ea,0, from experimental results. This analysis gives Ea,0 = 3.43 kcal mol-1, which is in good agreement with that predicted by the TIP4P/2005 water model. We also illustrate other approaches for estimating this activation energy consistent with these estimates.

Departures from perfect isomorph behavior in Lennard-Jones fluids and solids

Isomorphs are lines on a fluid or solid phase diagram along which the microstructure is invariant on affine density scaling of the molecular coordinates. Only inverse power (IP) and hard sphere potential systems are perfectly isomorphic. This work provides new theoretical tools and criteria to determine the extent of deviation from perfect isomorphicity for other pair potentials using the Lennard-Jones (LJ) system as a test case. A simple prescription for predicting isomorphs in the fluid range using the freezing line as a reference is shown to be quite accurate for the LJ system. The shear viscosity and self-diffusion coefficient scale well are calculated using this method, which enables comments on the physical significance of the correlations found previously in the literature to be made. The virial–potential energy fluctuation and the concept of an effective IPL system and exponent, n′, are investigated, particularly with reference to the LJ freezing and melting lines. It is shown that the exponent, n′, converges to the value 12 at a high temperature as ∼ T−1/2, where T is the temperature. Analytic expressions are derived for the density, temperature, and radius derivatives of the radial distribution function along an isomorph that can be used in molecular simulation. The variance of the radial distribution function and radial fluctuation function are shown to be isomorph invariant.

Temperature Extrapolation of Henry's Law Constants and the Isosteric Heat of Adsorption

Computational screening of adsorbent materials often uses the Henry's law constant (K_H) (at a particular temperature) as a first discriminator metric due to its relative ease of calculation. The isosteric heat of adsorption in the limit of zero pressure (q_st^∞) is often calculated along with the Henry's law constant, and both properties are informative metrics of adsorbent material performance at low-pressure conditions. In this article, we introduce a method for extrapolating K_H as a function of temperature, using series-expansion coefficients that are easily computed at the same time as K_H itself; the extrapolation function also yields q_st^∞. The extrapolation is highly accurate over a wide range of temperatures when the basis temperature is sufficiently high, for a wide range of adsorbent materials and adsorbate gases. Various results suggest that the extrapolation is accurate when the extrapolation range in inverse-temperature space is limited to |β - β₀ | < 0.5 mol/kJ. Application of the extrapolation to a large set of materials is shown to be successful provided that K_H is not extremely large and/or the extrapolation coefficients converge satisfactorily. The extrapolation is also able to predict q_st^∞ for a system that shows an unusually large temperature dependence. The work provides a robust method for predicting K_H and q_st^∞ over a wide range of industrially relevant temperatures with minimal effort beyond that necessary to compute those properties at a single temperature, which facilitates the addition of practical operating (or processing) conditions to computational screening exercises.

Temperature Dependence of Peptide Conformational Equilibria from Simulations at a Single Temperature

Understanding the structure of proteins is key to unraveling their function in biological processes. Thus, significant attention has been paid to the calculation of conformational free energies. In this paper, we demonstrate a simple extension of fluctuation theory that permits the calculation of the temperature derivative of the conformational free energy, and hence the internal energy and entropy, from single-temperature simulations. The method further enables the decomposition into the contribution of different interactions present in the system to the internal energy surface. We illustrate the method for the canonical test system of alanine dipeptide in aqueous solution, for which we examine the free energy as a function of two dihedral angles. This system, like many, is most effectively treated using accelerated sampling methods and we show how the present approach is compatible with an important class of these, those that introduce a bias potential, by implementing it within metadynamics.

Structure of the Lennard-Jones liquid estimated from a single simulation

Combining the recent Piskulich-Thompson approach [Z. A. Piskulich and W. H. Thompson, J. Chem. Phys. 152, 011102 (2020)JCPSA60021-960610.1063/1.5135932] with isomorph theory, from a single simulation the structure of a single-component Lennard-Jones (LJ) system is obtained at an arbitrary state point in almost the whole liquid region of the temperature-density phase diagram. The LJ system exhibits two temperature ranges where the van't Hoff assumption that energetic and entropic forces are temperature independent is valid to a good approximation. A method to evaluate the structure at an arbitrary state point along an isochore from the knowledge of structures at two temperatures on the isochore is also discussed. We argue that, in general, the structure of any hidden scale-invariant system obeying the van't Hoff assumption in the whole range of temperatures can be determined in the whole liquid region of the phase diagram from a single simulation.

Extrapolation and interpolation strategies for efficiently estimating structural observables as a function of temperature and density

Thermodynamic extrapolation has previously been used to predict arbitrary structural observables in molecular simulations at temperatures (or relative chemical potentials in open-system mixtures) different from those at which the simulation was performed. This greatly reduces the computational cost in mapping out phase and structural transitions. In this work, we explore the limitations and accuracy of thermodynamic extrapolation applied to water, where qualitative shifts from anomalous to simple-fluid-like behavior are manifested through shifts in the liquid structure that occur as a function of both temperature and density. We present formulas for extrapolating in volume for canonical ensembles and demonstrate that linear extrapolations of water's structural properties are only accurate over a limited density range. On the other hand, linear extrapolation in temperature can be accurate across the entire liquid state. We contrast these extrapolations with classical perturbation theory techniques, which are more conservative and slowly converging. Indeed, we show that such behavior is expected by demonstrating exact relationships between extrapolation of free energies and well-known techniques to predict free energy differences. An ideal gas in an external field is also studied to more clearly explain these results for a toy system with fully analytical solutions. We also present a recursive interpolation strategy for predicting arbitrary structural properties of molecular fluids over a predefined range of state conditions, demonstrating its success in mapping qualitative shifts in water structure with density.