Phase behavior of model confined fluids. Influence of substrate-fluid interaction strength

Molecular dynamics of an embedded-charge model of lysozyme aqueous solutions

The onset of liquid-vapor separation in an interaction site model of a lysozyme aqueous solution is investigated by means of molecular dynamics (MD). Calculations are performed for a soft-core version of a potential early introduced by Carlsson et al. (J. Phys. Chem. B 2001, 105, 9040; 2001, 105, 12189.) whose liquid-vapor coexistence was studied by Rosch and Errington (J. Phys. Chem. B 2007, 111, 12591.); our modified model preserves the tailoring onto the experimental lysozyme solution properties embodied by those descriptions. We first show that the structural results obtained by Carlsson et al. at ambient conditions are quite well reproduced by our approach. Then, we cool the system along an isochoric path by monitoring the structural and dynamical properties at various temperatures. We thus show that a fluid-fluid separation takes place at a temperature 15% below the presumed binodal; in particular, we observe that density inhomogeneities develop rather early in the MD run and evolve over tens of nanoseconds into two dense aggregates that eventually merge, after approximately 100 ns more, into a single liquid phase separated from a vapor region by a well-defined planar interface. The densities of the two coexisting fluids are compatible with previous determinations of the binodal line. The connections of this work to the overall scenario of phase stability investigations in protein solutions, as well as possible developments based on the use of more refined models, are discussed.

Excess-entropy scaling of dynamics for a confined fluid of dumbbell-shaped particles

We use molecular simulation to study the ability of excess entropy scaling relationships to describe the kinetic properties of a confined molecular system. We examine a model for a confined fluid consisting of dumbbell-shaped molecules that interact with atomistically detailed pore walls via a Lennard-Jones potential. We obtain kinetic, thermodynamic, and structural properties of the system at three wall-fluid interaction strengths and over a temperature range that includes sub- and super-critical conditions. Four dynamic properties are considered: translational and rotational diffusivities, a characteristic relaxation time for rotational motion, and a collective relaxation time stemming from analysis of the coherent intermediate scattering function. We carefully consider the reference state used to define the excess entropy of a confined fluid. Three ideal-gas reference states are considered, with the cases differentiated by the extent to which one-body spatial and orientational correlations are accounted for in the reference state. Our results indicate that a version of the excess entropy that includes information related to the one-body correlations in a confined fluid serves as the best scaling variable for dynamic properties. When adopting such a definition for the reference state, to a very good approximation, bulk and confined data for a specified dynamic property at a given temperature collapse onto a common curve when plotted against the excess entropy.