Jeffrey R. Errington

Direct calculation of liquid–vapor phase equilibria from transition matrix Monte Carlo simulation

An approach for directly determining the liquid–vapor phase equilibrium of a model system at any temperature along the coexistence line is described. The method relies on transition matrix Monte Carlo ideas developed by Fitzgerald, Picard, and Silver [Europhys. Lett. 46, 282 (1999)]. During a Monte Carlo simulation attempted transitions between states along the Markov chain are monitored as opposed to tracking the number of times the chain visits a given state as is done in conventional simulations. Data collection is highly efficient and very precise results are obtained. The method is implemented in both the grand canonical and isothermal–isobaric ensemble. The main result from a simulation conducted at a given temperature is a density probability distribution for a range of densities that includes both liquid and vapor states. Vapor pressures and coexisting densities are calculated in a straightforward manner from the probability distribution. The approach is demonstrated with the Lennard-Jones fluid. ...

Excess-entropy-based anomalies for a waterlike fluid

Many thermodynamic and dynamic properties of water display unusual behavior at low enough temperatures. In a recent study, Yan et al. [Phys. Rev. Lett. 95, 130604 (2005)] identified a spherically symmetric two-scale potential that displays many of the same anomalous properties as water. More specifically, for select parametrizations of the potential, one finds that the regions where isothermal compression anomalously (i) decreases the fluids structural order, (ii) increases its translational self-diffusivity, and (iii) increases its entropy form nested domes in the temperature-density plane. These property relationships are similar to those found for more realistic models of water. In this work, the authors provide evidence that suggests that the anomalous regions specified above can all be linked through knowledge of the excess entropy. Specifically, the authors show how entropy scaling relationships developed by Rosenfeld [Phys. Rev. A 15, 2545 (1977)] can be used to describe the region of diffusivity anomalies and to predict the state conditions for which anomalous viscosity and thermal conductivity behavior might be found.

Surface tension and vapor–liquid phase coexistence of the square-well fluid

Vapor–liquid interfacial tension of square-well ~SW! fluids is calculated using three different methods viz., molecular dynamics ~MD! with collision-based virial evaluation, Monte Carlo with virial computed by volume perturbation, and Binder’s density-distribution method in conjunction with grand-canonical transition-matrix Monte Carlo ~GC-TMMC!. Three values of the SW attractive well range parameter were studied: l51.5, 1.75, and 2.0, respectively. The results from MD and GC-TMMC methods are in very good mutual agreement, while the volume-perturbation method yields data of unacceptable quality. The results are compared with predictions from the statistical associating fluid theory ~SAFT!, and SAFT is shown to give a good estimate for the systems studied. Liquid and vapor coexistence densities and saturation pressure are determined from analysis of GC-TMMC data and the results are found to agree very well with the established literature data.

A computational study of hydration, solution structure, and dynamics in dilute carbohydrate solutions.

We report results from a molecular simulation study of the structure and dynamics of water near single carbohydrate molecules (glucose, trehalose, and sucrose) at 0 and 30 degrees C. The presence of a carbohydrate molecule has a number of significant effects on the microscopic water structure and dynamics. All three carbohydrates disrupt the tetrahedral arrangement of proximal water molecules and restrict their translational and rotational mobility. These destructuring effects and slow dynamics are the result of steric constraints imposed by the carbohydrate molecule and of the ability of a carbohydrate to form stable H bonds with water, respectively. The carbohydrates induce a pronounced decoupling between translational and rotational motions of proximal water molecules.

Evaluating surface tension using grand-canonical transition-matrix Monte Carlo simulation and finite-size scaling.

This Brief Report describes an approach for determining the surface tension of a model system that is applicable over the entire liquid-vapor coexistence region. At the heart of the method is a technique for determining coexistence properties that utilize transition probabilities of attempted Monte Carlo moves during a grand canonical simulation. Finite-size scaling techniques are implemented to determine the infinite system surface tension from a series of finite-size simulations. To demonstrate the method, the surface tension of the Lennard-Jones fluid is determined at temperatures ranging from the triple point to the critical point.

Quantification of order in the Lennard-Jones system

We conduct a numerical investigation of structural order in the shifted-force Lennard-Jones system by calculating metrics of translational and bond-orientational order along various paths in the phase diagram covering equilibrium solid, liquid, and vapor states. A series of nonequilibrium configurations generated through isochoric quenches, isothermal compressions, and energy minimizations are also considered. Simulation results are analyzed using an ordering map representation [Torquato et al., Phys. Rev. Lett. 84, 2064 (2000); Truskett et al., Phys. Rev. E 62, 993 (2000)] that assigns both equilibrium and nonequilibrium states coordinates in an order metric plane. Our results show that bond-orientational order and translational order are not independent for simple spherically symmetric systems at equilibrium. We also demonstrate quantitatively that the Lennard-Jones and hard sphere systems sample the same configuration space at supercritical densities. Finally, we relate the structural order found in fast-quenched and minimum-energy configurations (inherent structures).

Calculation of surface tension via area sampling.

We examine the performance of several molecular simulation techniques aimed at evaluation of the surface tension through its thermodynamic definition. For all methods explored, the surface tension is calculated by approximating the change in Helmholtz free energy associated with a change in interfacial area through simulation of a liquid slab at constant particle number, volume, and temperature. The methods explored fall within three general classes: free-energy perturbation, the Bennett acceptance-ratio scheme, and the expanded ensemble technique. Calculations are performed for both the truncated Lennard-Jones and square-well fluids at select temperatures spaced along their respective liquid-vapor saturation lines. Overall, we find that Bennett and expanded ensemble approaches provide the best combination of accuracy and precision. All of the methods, when applied using sufficiently small area perturbation, generate equivalent results for the Lennard-Jones fluid. However, single-stage free-energy-perturbation methods and the closely related test-area technique recently introduced by Gloor et al. [J. Chem. Phys. 123, 134703 (2005)] generate surface tension values for the square-well fluid that are not consistent with those obtained from the more robust expanded ensemble and Bennett approaches, regardless of the size of the area perturbation. Single-stage perturbation methods fail also for the Lennard-Jones system when applied using large area perturbations. Here an analysis of phase-space overlap produces a quantitative explanation of the observed inaccuracy and shows that the satisfactory results obtained in these cases from the test-area method arise from a cancellation of errors that cannot be expected in general. We also briefly analyze the variation in method performance with respect to the adjustable parameters inherent to the techniques.

Layering and Position-Dependent Diffusive Dynamics of Confined Fluids

We study the diffusive dynamics of a hard-sphere fluid confined between parallel smooth hard walls. The position-dependent diffusion coefficient normal to the walls is larger in regions of high local packing density. High density regions also have the largest available volume, consistent with the fast local diffusivity. Indeed, local and global diffusivities as a function of the Widom insertion probability approximately collapse onto a master curve. Parallel and average normal diffusivities are strongly coupled at high densities and deviate from bulk fluid behavior.

Thermodynamics predicts how confinement modifies the dynamics of the equilibrium hard-sphere fluid.

We study how confining the equilibrium hard-sphere fluid to restrictive one- and two-dimensional channels with smooth interacting walls modifies its structure, dynamics, and entropy using molecular dynamics and transition-matrix Monte Carlo simulations. Although confinement strongly affects local structuring, the relationships between self-diffusivity, excess entropy, and average fluid density are, to an excellent approximation, independent of channel width or particle-wall interactions. Thus, thermodynamics can be used to predict how confinement impacts dynamics.

Relationship between thermodynamics and dynamics of supercooled liquids

Diffusivity, a measure for how rapidly a fluid self-mixes, shows an intimate, but seemingly fragmented, connection to thermodynamics. On one hand, the “configurational” contribution to entropy (related to the number of mechanically stable configurations that fluid molecules can adopt) has long been considered key for predicting supercooled liquid dynamics near the glass transition. On the other hand, the excess entropy (relative to ideal gas) provides a robust scaling for the diffusivity of fluids above the freezing point. Here we provide, to our knowledge, the first evidence that excess entropy also captures how supercooling a fluid modifies its diffusivity, suggesting that dynamics, from ideal gas to glass, is related to a single, standard thermodynamic quantity.