Perspective: Excess-entropy scaling

Entropy Scaling of Viscosity - II: Predictive Scheme for Normal Alkanes

In this work, a residual entropy value 6/10 of the way between the critical point and a value of -2/3 of Boltzmann's constant is shown to collapse the scaled viscosity for the family of normal alkanes. Based on this approach, a nearly universal correlation is proposed that can reproduce 95% of the experimental data for normal alkanes within ±18% (without removal of clearly erroneous data). This universal correlation has no new fluid-specific empirical parameters and is based on experimentally accessible values. This collapse is shown to be valid to a residual entropy half way between the critical point and the triple point, beyond which the macroscopically-scaled viscosity has a super-exponential dependence on residual entropy, terminating at the triple point. A key outcome of this study is a better understanding of entropy scaling for fluids with intramolecular degrees of freedom. A study of the transport and thermodynamic properties at the triple point rounds out the analysis.

Isomorph theory beyond thermal equilibrium

This paper generalizes isomorph theory to systems that are not in thermal equilibrium. The systems are assumed to be R-simple, i.e., to have a potential energy that as a function of all particle coordinates R obeys the hidden-scale-invariance condition U(R_a) < U(R_b) ⇒ U(λR_a) < U(λR_b). "Systemic isomorphs" are introduced as lines of constant excess entropy in the phase diagram defined by density and systemic temperature, which is the temperature of the equilibrium state point with the average potential energy equal to U(R). The dynamics is invariant along a systemic isomorph if there is a constant ratio between the systemic and the bath temperature. In thermal equilibrium, the systemic temperature is equal to the bath temperature and the original isomorph formalism is recovered. The new approach rationalizes within a consistent framework previously published observations of isomorph invariance in simulations involving nonlinear steady-state shear flows, zero-temperature plastic flows, and glass-state isomorphs. This paper relates briefly to granular media, physical aging, and active matter. Finally, we discuss the possibility that the energy unit defining the reduced quantities should be based on the systemic rather than the bath temperature.

Thermodynamic or density scaling of the thermal conductivity of liquids

Thermodynamic or density scaling is applied to thermal conductivity (λ) data from the literature for the model Lennard-Jones (12-6) fluid; the noble gases neon to xenon; nitrogen, ethene, and carbon dioxide as examples of linear molecules; the quasi-spherical molecules methane and carbon tetrachloride; the flexible chain molecules n-hexane and n-octane; the planar toluene and m-xylene; the cyclic methylcyclohexane; the polar R132a and chlorobenzene; and ammonia and methanol as H-bonded fluids. Only data expressed as Rosenfeld reduced properties could be scaled successfully. Two different methods were used to obtain the scaling parameter γ, one based on polynomial fits to the group (TV^γ) and the other based on the Avramov equation. The two methods agree well, except for λ of CCl₄. γ for the thermal conductivity is similar to those for the viscosity and self-diffusion coefficient for the smaller molecules. It is significantly larger for the Lennard-Jones fluid, possibly due to a different dependence on packing fraction, and much larger for polyatomic molecules where heat transfer through internal modes may have an additional effect. Methanol and ammonia, where energy can be transmitted through intermolecular hydrogen bonding, could not be scaled. This work is intended as a practical attempt to examine thermodynamic scaling of the thermal conductivity of real fluids. The divergence of the scaling parameters for different properties is unexpected, suggesting that refinement of theory is required to rationalize this result. For the Lennard-Jones fluid, the Ohtori-Iishi version of the Stokes-Einstein-Sutherland relation applies at high densities in the liquid and supercritical region.

Coarse-graining strategy for modeling effective, highly diffusive fluids with reduced polydispersity: A dynamical study

We present a coarse-graining strategy for reducing the number of particle species in mixtures to achieve a simpler system with higher diffusion while preserving the total particle number and characteristic dynamic features. As a system of application, we chose the bidisperse Lennard-Jones-like mixture, discovered by Kob and Andersen [Phys. Rev. Lett. 73, 1376 (1994)], possessing a slow dynamics due to the fluid's multi-component character with its apparently unconventional choice for the pair potential of the type-A-type-B arrangement. We further established in a so-formed coarse-grained and temperature-independent monodisperse system an equilibrium structure with a radial distribution function resembling its mixture counterpart. This one-component system further possesses similar dynamic features such as glass transition temperature and critical exponents while subjected to Newtonian mechanics. This strategy may finally lead to the manufacturing of new nanoparticle/colloidal fluids by experimentally modeling only the outcoming effective pair potential(s) and no other macroscopic quantity.

Excess-entropy scaling in supercooled binary mixtures

Transport coefficients, such as viscosity or diffusion coefficient, show significant dependence on density or temperature near the glass transition. Although several theories have been proposed for explaining this dynamical slowdown, the origin remains to date elusive. We apply here an excess-entropy scaling strategy using molecular dynamics computer simulations and find a quasiuniversal, almost composition-independent, relation for binary mixtures, extending eight orders of magnitude in viscosity or diffusion coefficient. Metallic alloys are also well captured by this relation. The excess-entropy scaling predicts a quasiuniversal breakdown of the Stokes-Einstein relation between viscosity and diffusion coefficient in the supercooled regime. Additionally, we find evidence that quasiuniversality extends beyond binary mixtures, and that the origin is difficult to explain using existing arguments for single-component quasiuniversality.

Thermodynamics-Structure-Dynamics Correlations and Nonuniversal Effects in the Elastically Collective Activated Hopping Theory of Glass-Forming Liquids

We employ the microscopic Elastically Collective Nonlinear Langevin Equation (ECNLE) theory of activated dynamics in combination with crystal-avoiding simulations to study four inter-related questions for metastable monodisperse hard sphere fluids. The first is how significantly improved integral equation theory structural input (Modified-Verlet (MV) closure) changes the dynamical predictions of ECNLE theory. The main consequence is a modest enhancement of the importance of the collective elastic barrier relative to its local cage contribution, which increases the alpha relaxation time and fragility relative to prior results based on the Percus-Yevick closure. Second, ECNLE-MV theory predictions for the alpha time and self-diffusion constant in the metastable regime are quantitatively compared to our new simulations. The small adjustment of a numerical prefactor that enters the collective elastic barrier leads to quantitative agreement over three decades. Third, using the more accurate MV structural input, ECNLE theory is shown to predict thermodynamics-structure-dynamics "correlations" based on various long and short wavelength scalar properties all related to static two-point collective density fluctuations. The logarithm of the alpha relaxation time scales as a power law with these scalar metrics with an exponent that is significantly lower in the less dense noncooperative activated regime compared to the very dense highly cooperative regime. However, the discovered correlation of activated relaxation with a thermodynamic property (dimensionless compressibility) is not causal in ECNLE theory, but rather reflects a strong connection between the local structural quantities that quantify kinetic constraints in the theory with the amplitude of long wavelength density fluctuations. Fourth, the consequences of chemically specific nonuniversalities associated with the onset condition and relative importance of collective elasticity are studied. The predicted thermodynamics-structure-dynamics correlations are found to be robust, albeit with nontrivial shifts of the onset condition.

An entropy scaling demarcation of gas- and liquid-like fluid behaviors

In this work, we propose a generic and simple definition of a line separating gas-like and liquid-like fluid behaviors from the standpoint of shear viscosity. This definition is valid even for fluids such as the hard sphere and the inverse power law that exhibit a unique fluid phase. We argue that this line is defined by the location of the minimum of the macroscopically scaled viscosity when plotted as a function of the excess entropy, which differs from the popular Widom lines. For hard sphere, Lennard-Jones, and inverse-power-law fluids, such a line is located at an excess entropy approximately equal to -2/3 times Boltzmann's constant and corresponds to points in the thermodynamic phase diagram for which the kinetic contribution to viscosity is approximately half of the total viscosity. For flexible Lennard-Jones chains, the excess entropy at the minimum is a linear function of the chain length. This definition opens a straightforward route to classify the dynamical behavior of fluids from a single thermodynamic quantity obtainable from high-accuracy thermodynamic models.

Scaling of relaxation and excess entropy in plastically deformed amorphous solids

When stressed sufficiently, solid materials yield and deform plastically via reorganization of microscopic constituents. Indeed, it is possible to alter the microstructure of materials by judicious application of stress, an empirical process utilized in practice to enhance the mechanical properties of metals. Understanding the interdependence of plastic flow and microscopic structure in these nonequilibrium states, however, remains a major challenge. Here, we experimentally investigate this relationship, between the relaxation dynamics and microscopic structure of disordered colloidal solids during plastic deformation. We apply oscillatory shear to solid colloidal monolayers and study their particle trajectories as a function of shear rate in the plastic regime. Under these circumstances, the strain rate, the relaxation rate associated with plastic flow, and the sample microscopic structure oscillate together, but with different phases. Interestingly, the experiments reveal that the relaxation rate associated with plastic flow at time t is correlated with the strain rate and sample microscopic structure measured at earlier and later times, respectively. The relaxation rate, in this nonstationary condition, exhibits power-law, shear-thinning behavior and scales exponentially with sample excess entropy. Thus, measurement of sample static structure (excess entropy) provides insight about both strain rate and constituent rearrangement dynamics in the sample at earlier times.

Effective hardness of interaction from thermodynamics and viscosity in dilute gases

The hardness of the effective inverse power law (IPL) potential, which can be obtained from thermodynamics or collision integrals, can be used to demonstrate similarities between thermodynamic and transport properties. This link is investigated for systems of increasing complexity (i.e., the EXP, square-well, Lennard-Jones, and Stockmayer potentials; ab initio results for small molecules; and rigid linear chains of Lennard-Jones sites). These results show that while the two approaches do not yield precisely the same values of effective IPL exponent, their qualitative behavior is intriguingly similar, offering a new way of understanding the effective interactions between molecules, especially at high temperatures. In both approaches, the effective hardness is obtained from a double-logarithmic temperature derivative of an effective area.

Residual entropy model for predicting the viscosities of dense fluid mixtures

In this work, we have investigated the mono-variant relationship between the reduced viscosity and residual entropy in pure fluids and in binary mixtures of hydrocarbons and hydrocarbons with dissolved carbon dioxide. The mixtures considered were octane + dodecane, decane + carbon dioxide, and 1,3-dimethylbenzene (m-xylene) + carbon dioxide. The reduced viscosity was calculated according to the definition of Bell, while the residual entropy was calculated from accurate multi-parameter Helmholtz-energy equations of state and, for mixtures, the multi-fluid Helmholtz energy approximation. The mono-variant dependence of reduced viscosity upon residual molar entropy was observed for the pure fluids investigated, and by incorporating two scaling factors (one for reduced viscosity and the other for residual molar entropy), the data were represented by a single universal curve. To apply this method to mixtures, the scaling factors were determined from a mole-fraction weighted sum of the pure-component values. This simple model was found to work well for the systems investigated. The average absolute relative deviation (AARD) was observed to be between 1% and 2% for pure components and a mixture of similar hydrocarbons. Larger deviations, with AARDs of up to 15%, were observed for the asymmetric mixtures, but this compares favorably with other methods for predicting the viscosity of such systems. We conclude that the residual-entropy concept can be used to estimate the viscosity of mixtures of similar molecules with high reliability and that it offers a useful engineering approximation even for asymmetric mixtures.

Viscous properties of nickel-containing binary metal melts

The paper presents the results of molecular dynamics study of the viscosity of nickel-containing binary metal melts for a wide range of temperatures, including the region of the equilibrium liquid phase and supercooled melt. It is shown that the temperature dependencies of the viscosity of binary metal melts are described by the Kelton's quasi-universal model. Based on the analysis of the viscosity coefficient of the binary melt composition within the framework of the Rosenfeld's scale transformations, it has been established that to correctly describe the viscosity of binary/multicomponent metal melts within the framework of entropy models, it is necessary to use a more complex representation of the excess entropy S _ex than in the approximation of pair correlation entropy S ₂.

Temperature and pressure dependence of the alpha relaxation in ortho-terphenyl

Molecular dynamics (MD) simulations of ortho-terphenyl using an all-atom model with the optimized potentials for liquid simulations (OPLS) force field were performed both in the high temperature Arrhenian region and at lower temperatures that include the onset of the super-Arrhenian region. From the MD simulations, the internal energy of both the equilibrium liquid and crystal was determined from 300 K to 600 K and at pressures from 0.1 MPa to 1 GPa. The translational and rotational diffusivities were also determined at these temperatures and pressures for the equilibrium liquid. It is shown that within a small offset, the excess internal energy Ū_x from the MD simulations is consistent with the experimentally determined excess internal energy reported earlier [Caruthers and Medvedev, Phys. Rev. Mater. 2, 055604, (2018)]. The MD mobility data {including extremely long-time 1 atm simulations from the study by Eastwood et al. [J. Phys. Chem. B 117, 12898, (2013)]} were combined with experimental data to form a unified dataset, where it was shown that in both the high temperature Arrhenian region and the lower temperature super-Arrhenian region, the mobility is a linear function of 1/Ū_x(T,p), albeit with different proportionality constants. The transition between the Arrhenian and super-Arrhenian regions is relatively sharp at a critical internal energy Ū_x ^α. The 1/Ū_x(T,p) model is able to describe the mobility data over nearly 16 orders-of-magnitude. Other excess thermodynamic properties such as excess enthalpy and excess entropy (i.e., the Adam-Gibbs model) are unable to unify the pressure dependence of the mobility.

A comprehensive study of the thermal conductivity of the hard sphere fluid and solid by molecular dynamics simulation

This work reports a new set of hard sphere (HS) thermal conductivity coefficient, λ, data obtained by Molecular Dynamics (MD) computer simulation, over a density range covering the dilute fluid to near the close-packed solid, and for a large number of particles (up to N = 13 1072) and long simulation times. The N-dependence of the thermal conductivity is shown to be proportional to N-2/3 to a good approximation over a wide range of system sizes, which enabled λ values in the thermodynamic limit to be predicted accurately. The fluid and solid λ can be represented well by the Enskog theory (ET) formula, λE, times a density-dependent correction term, which is close to unity for the fluid and practically constant for the solid. The convergence of the MD λ data back towards ET in the metastable fluid starts just above the freezing density. For the HS solid and dense fluid it was found that the thermal conductivity is nearly linear in pressure, as has been observed experimentally for a number of solids. Simple excess entropy scaling over the higher density fluid phase region was found, and Rosenfeld's exponential relationship can be fitted to the simulation data for the solid to a high degree of accuracy. The simulation analysis has revealed a number of new trends in the behaviour of the HS thermal conductivity which could be useful in building more accurate models for heat conduction in experimental systems.

Relationship between the Transport Coefficients of Polar Substances and Entropy

An expression is proposed that relates the transport properties of polar substances (diffusion coefficient, viscosity coefficient, and thermal conductivity coefficient) with entropy. To calculate the entropy, an equation of state with a good description of the properties in a wide region of the state is used. Comparison of calculations based on the proposed expressions with experimental data showed good agreement. A deviation exceeding 20% is observed only in the region near the critical point as well as at high pressures.

Recent Progress towards Chemically-Specific Coarse-Grained Simulation Models with Consistent Dynamical Properties

Coarse-grained (CG) models can provide computationally efficient and conceptually simple characterizations of soft matter systems. While generic models probe the underlying physics governing an entire family of free-energy landscapes, bottom-up CG models are systematically constructed from a higher-resolution model to retain a high level of chemical specificity. The removal of degrees of freedom from the system modifies the relationship between the relative time scales of distinct dynamical processes through both a loss of friction and a “smoothing” of the free-energy landscape. While these effects typically result in faster dynamics, decreasing the computational expense of the model, they also obscure the connection to the true dynamics of the system. The lack of consistent dynamics is a serious limitation for CG models, which not only prevents quantitatively accurate predictions of dynamical observables but can also lead to qualitatively incorrect descriptions of the characteristic dynamical processes. With many methods available for optimizing the structural and thermodynamic properties of chemically-specific CG models, recent years have seen a stark increase in investigations addressing the accurate description of dynamical properties generated from CG simulations. In this review, we present an overview of these efforts, ranging from bottom-up parameterizations of generalized Langevin equations to refinements of the CG force field based on a Markov state modeling framework. We aim to make connections between seemingly disparate approaches, while laying out some of the major challenges as well as potential directions for future efforts.

Modified Entropy Scaling of the Transport Properties of the Lennard-Jones Fluid

Rosenfeld proposed two different scaling approaches to model the transport properties of fluids, separated by 22 years, one valid in the dilute gas, and another in the liquid phase. In this work, we demonstrate that these two limiting cases can be connected through the use of a novel approach to scaling transport properties and a bridging function. This approach, which is empirical and not derived from theory, is used to generate reference correlations for the transport properties of the Lennard-Jones 12-6 fluid of viscosity, thermal conductivity, and self-diffusion. This approach, with a very simple functional form, allows for the reproduction of the most accurate simulation data to within nearly their statistical uncertainty. The correlations are used to confirm that for the Lennard-Jones fluid the appropriately scaled transport properties are nearly monovariate functions of the excess entropy from low-density gases into the supercooled phase and up to extreme temperatures. This study represents the most comprehensive metastudy of the transport properties of the Lennard-Jones fluid to date.

Entropy and the Tolman Parameter in Nucleation Theory

Thermodynamic aspects of the theory of nucleation are commonly considered employing Gibbs’ theory of interfacial phenomena and its generalizations. Utilizing Gibbs’ theory, the bulk parameters of the critical clusters governing nucleation can be uniquely determined for any metastable state of the ambient phase. As a rule, they turn out in such treatment to be widely similar to the properties of the newly-evolving macroscopic phases. Consequently, the major tool to resolve problems concerning the accuracy of theoretical predictions of nucleation rates and related characteristics of the nucleation process consists of an approach with the introduction of the size or curvature dependence of the surface tension. In the description of crystallization, this quantity has been expressed frequently via changes of entropy (or enthalpy) in crystallization, i.e., via the latent heat of melting or crystallization. Such a correlation between the capillarity phenomena and entropy changes was originally advanced by Stefan considering condensation and evaporation. It is known in the application to crystal nucleation as the Skapski–Turnbull relation. This relation, by mentioned reasons more correctly denoted as the Stefan–Skapski–Turnbull rule, was expanded by some of us quite recently to the description of the surface tension not only for phase equilibrium at planar interfaces, but to the description of the surface tension of critical clusters and its size or curvature dependence. This dependence is frequently expressed by a relation derived by Tolman. As shown by us, the Tolman equation can be employed for the description of the surface tension not only for condensation and boiling in one-component systems caused by variations of pressure (analyzed by Gibbs and Tolman), but generally also for phase formation caused by variations of temperature. Beyond this particular application, it can be utilized for multi-component systems provided the composition of the ambient phase is kept constant and variations of either pressure or temperature do not result in variations of the composition of the critical clusters. The latter requirement is one of the basic assumptions of classical nucleation theory. For this reason, it is only natural to use it also for the specification of the size dependence of the surface tension. Our method, relying on the Stefan–Skapski–Turnbull rule, allows one to determine the dependence of the surface tension on pressure and temperature or, alternatively, the Tolman parameter in his equation. In the present paper, we expand this approach and compare it with alternative methods of the description of the size-dependence of the surface tension and, as far as it is possible to use the Tolman equation, of the specification of the Tolman parameter. Applying these ideas to condensation and boiling, we derive a relation for the curvature dependence of the surface tension covering the whole range of metastable initial states from the binodal curve to the spinodal curve.

Density fluctuations and random walks in an overdamped and supercooled simple liquid

In this work, the short-time dynamics of simple liquid is explored both analytically and numerically with the focus on the interplay between the density fluctuations in a volume surrounding a chosen particle and its random walk motion. The particles interact via the Lennard-Jones potential with parameters corresponding to liquid argon. For large times, analytical calculations based on the fluctuation theory provides an explicit expression reproducing isothermal change of the self-diffusion coefficient in liquid argon corresponding to the experimental data. These results lead to the conclusion that such behavior is based on the reduced mobility of particles reflected in their density fluctuations that can be equivalently achieved in the cases of either low temperatures and pressures (supercooling) or moderate temperatures and high pressures (overdamping).

Universal Scaling Law for Colloidal Diffusion in Complex Media

Using video microscopy and simulations, we study the diffusion of probe particles in a wide range of complex backgrounds, both crystalline and disordered, in quasi-2D colloidal systems. The dimensionless diffusion coefficients D^{*} from different systems collapse to a single master curve when plotted as a function of the structural entropy of the backgrounds, confirming the universal relation between diffusion dynamics and the structure of the medium. A new scaling equation is proposed with consideration for the viscous friction from the solvent, which is absent in previous theoretical models. This new universal law quantitatively predicts the diffusion coefficients from different systems over several orders of magnitude of D^{*}, with a single common fitting parameter.

Excess entropy and long-time diffusion in colloidal fluids with short-range interparticle attraction

Liquid structure and dynamics are experimentally investigated in colloidal suspensions with short-range depletion attraction. The colloidal fluid samples consist of hard-sphere colloidal particles suspended along with rodlike depletants based on surfactant micelles. The spheres have a range of surface chemistries, diameters, and packing fractions, and the rodlike micelle length depends on the temperature. Thus, the combination of hard-spheres and depletants generates a sample wherein short-range interparticle attraction can be temperature-tuned in situ. Video optical microscopy and particle tracking techniques are employed to measure particle trajectories from which structural and dynamical quantities are derived, including the particle pair correlation function [g(r)], mean square displacement, long-time diffusion coefficient, and the sample two-body excess entropy (S₂). The samples with stronger short-range attractions exhibit more order, as characterized by g(r) and S₂. The stronger short-range attractions are also observed to lead to slower long-time diffusion and more heterogeneous dynamics at intermediate time scales. Finally, the excess entropy scaling law prediction, i.e., the exponential relationship between two-body excess entropy and long-time diffusivity, is observed across the full range of samples.

Thermodynamic and dynamical properties of the hard sphere system revisited by molecular dynamics simulation

Revised thermodynamic and dynamical properties of the hard sphere (HS) system are obtained from extensive molecular dynamics calculations carried out with large system sizes (number of particles, N) and long times. Accurate formulas for the compressibility factor of the HS solid and fluid branches are proposed, which represent the metastable region and take into account its divergence at close packing. Some basic second-order thermodynamic properties are obtained and a maximum in some of their derivatives in the metastable fluid region is found. The thermodynamic parameters associated with the melting-freezing transition have been determined to four digit accuracy, which generates accurate new values for the coexistence properties of the HS system. For the self-diffusion coefficient, D, it is shown that relatively large systems (N > 104) are required to achieve an accurate linear extrapolation of D to the infinite size limit with a D vs. N-1/3 plot. Moreover, it is found that there is a density dependence of the value of the slope in the linear regime. The density dependent correction becomes practically insignificant at higher densities and the hydrodynamic formula found in the literature is still accurate. However, with decreasing density the density dependence of the size correction cannot be neglected, which indicates that other sources of N-dependence, apart from those derived on purely hydrodynamic grounds, may also be important (and as yet unaccounted for). A detailed analytic representation of the density dependence of the HS self-diffusion coefficient and the HS viscosity, η, is given. It is shown that the HS viscosity near freezing and in the metastable region can be described well by the Krieger-Dougherty equation. Both D and η start to scale at high densities and in the metastable region in such a way that Dηp = const, where p ≃ 0.97, and D → 0 and η → ∞ at a packing fraction of 0.58, which coincides with some previous predictions of the HS glass transition density.

Protein Solvent Shell Structure Provides Rapid Analysis of Hydration Dynamics

The solvation layer surrounding a protein is clearly an intrinsic part of protein structure-dynamics-function, and our understanding of how the hydration dynamics influences protein function is emerging. We have recently reported simulations indicating a correlation between regional hydration dynamics and the structure of the solvation layer around different regions of the enzyme Candida antarctica lipase B, wherein the radial distribution function (RDF) was used to calculate the pairwise entropy, providing a link between dynamics (diffusion) and thermodynamics (excess entropy) known as Rosenfeld scaling. Regions with higher RDF values/peaks in the hydration layer (the first peak, within 6 Å of the protein surface) have faster diffusion in the hydration layer. The finding thus hinted at a handle for rapid evaluation of hydration dynamics at different regions on the protein surface in molecular dynamics simulations. Such an approach may move the analysis of hydration dynamics from a specialized venture to routine analysis, enabling an informatics approach to evaluate the role of hydration dynamics in biomolecular function. This paper first confirms that the correlation between regional diffusive dynamics and hydration layer structure (via water center of mass around protein side-chain atom RDF) is observed as a general relationship across a set of proteins. Second, it seeks to devise an approach for rapid analysis of hydration dynamics, determining the minimum amount of information and computational effort required to get a reliable value of hydration dynamics from structural data in MD simulations based on the protein-water RDF. A linear regression model using the integral of the hydration layer in the water-protein RDF was found to provide statistically equivalent apparent diffusion coefficients at the 95% confidence level for a set of 92 regions within five different proteins. In summary, RDF analysis of 10 ns of data after simulation convergence is sufficient to accurately map regions of fast and slow hydration dynamics around a protein surface. Additionally, it is anticipated that a quick look at protein-water RDFs, comparing peak heights, will be useful to provide a qualitative ranking of regions of faster and slower hydration dynamics at the protein surface for rapid analysis when investigating the role of solvent dynamics in protein function.

Probing the link between residual entropy and viscosity of molecular fluids and model potentials

This work investigates the link between residual entropy and viscosity based on wide-ranging, highly accurate experimental and simulation data. This link was originally postulated by Rosenfeld in 1977 [Rosenfeld Y (1977) Phys Rev A 15:2545-2549], and it is shown that this scaling results in an approximately monovariate relationship between residual entropy and reduced viscosity for a wide range of molecular fluids [argon, methane, [Formula: see text], [Formula: see text], refrigerant R-134a (1,1,1,2-tetrafluoroethane), refrigerant R-125 (pentafluoroethane), methanol, and water] and a range of model potentials (hard sphere, inverse power, Lennard-Jones, and Weeks-Chandler-Andersen). While the proposed "universal" correlation of Rosenfeld is shown to be far from universal, when used with the appropriate density scaling for molecular fluids, the viscosity of nonassociating molecular fluids can be mapped onto the model potentials. This mapping results in a length scale that is proportional to the cube root of experimentally measurable liquid volume values.

Experimental Evidence for a State-Point-Dependent Density-Scaling Exponent of Liquid Dynamics

A large class of liquids obey density scaling characterized by an exponent, which quantifies the relative roles of temperature and density for the dynamics. We present experimental evidence that the density-scaling exponent γ is state-point dependent for the glass formers tetramethyl-tetraphenyl-trisiloxane (DC704) and 5-polyphenyl ether (5PPE). A method is proposed that from dynamic and thermodynamic properties at equilibrium estimates the value of γ. The method applies at any state point of the pressure-temperature plane, both in the supercooled and the normal liquid regimes. We find that γ is generally state-point dependent, which is confirmed by reanalyzing data for 20 metallic liquids and two model liquids.