Computer simulations are used to test whether a recently introduced generalization of Rosenfeld's excess-entropy scaling method for estimating transport coefficients in systems obeying molecular dynamics can be extended to predict long-time diffusivities in fluids of particles undergoing Brownian dynamics in the absence of interparticle hydrodynamic forces. Model fluids with inverse-power-law, Gaussian-core, and Hertzian pair interactions are considered. Within the generalized Rosenfeld scaling method, long-time diffusivities of ultrasoft Gaussian-core and Hertzian particle fluids, which display anomalous trends with increasing density, are predicted (to within 20%) based on knowledge of interparticle interactions, excess entropy, and scaling behavior of simpler inverse-power-law fluids.
Dense suspensions of hard-sphere colloids have long served as useful experimental models for exploring how constraints of particle packing affect the properties of condensed phases.1–3 Nonetheless, many technologically relevant suspensions contain high concentrations of “particles” (e.g., micelles, microgel colloids, dendrimers, or star polymers), which interact with one another through considerably softer effective potentials.4–6 Such systems are of fundamental interest because they exhibit novel properties that cannot be readily understood in terms of concepts derived from idealized hard-particle systems such as “excluded volume,” “crowding,” and “collision rate,” etc. Examples highlighted in theoretical and computer simulation studies include re-entrant and cluster crystalline phases,7–17 as well as anomalous dependencies of fluid-phase thermodynamic,18–22 structural,20,22–29 and dynamic14,18,20–27,29 properties on particle concentration.
To date, much of the theoretical work on these systems has been devoted to developing and testing statistical mechanical approaches for predicting their structural and thermodynamic behavior.4,30–32 However, some of the focus is beginning to shift toward predicting dynamics. For example, it has recently been shown that generalized Langevin theories24,33,34 can qualitatively (but not yet quantitatively) predict the anomalous density dependence of the long-time self-diffusivity observed in Brownian dynamics simulations of the Gaussian-core model. Likewise, mode-coupling theory has been able to capture some of the nonmonotonic dynamic trends displayed by fluids of particles that interact via star-polymer-like,23 harmonic,35 Hertzian,29 or Gaussian-core27,36 pair potentials. Kinetic theory has also been used to gain insights into the nontrivial temperature- and density-dependent trends in the long-time molecular dynamics of particles with bounded, penetrable-sphere interactions.37
Another promising, albeit more heuristic, approach for predicting the dynamic properties of soft-particle fluids is a recently proposed generalization21 of the excess-entropy scaling method of Rosenfeld.38,39 Here, excess entropy sex refers to the difference between the entropy per particle of the fluid and that of an ideal gas of particles with the same number density ρ. Excess entropy is a negative quantity, and its magnitude characterizes the extent to which static interparticle correlations—present due to interparticle interactions— reduce the number of microstates accessible to the fluid. Although a rigorous statistical mechanical derivation of excess entropy scalings for dynamics has proven to be elusive,21,39–41 the qualitative expectation—corroborated by data from both computer simulations and experiments (see, e.g., Refs. 39 and 42–48)—is that changes to a fluid that strengthen its interparticle correlations also result in slower dynamic relaxation processes.
In the generalized Rosenfeld (GR) scaling approach,21 information about interparticle interactions is used to recast the species’ long-time tracer diffusivities as dimensionless (GR-reduced) combinations that are, by construction, single-valued functions of sex in the low ρ limit. For fluids of particles that interact solely via an inverse-power-law pair (IPL) repulsion, one can easily show that GR-reduced diffusivities remain single-valued functions of sex for all values of ρ and temperature T. Interestingly, Krekelberg et al. have further discovered that GR-reduced tracer diffusivities obtained from molecular dynamics simulations of more complex model systems—equilibrium fluid mixtures of (additive or nonadditive) hard-sphere or Gaussian-core particles at different values of ρ, T, and composition—approximately follow a similar excess-entropy based scaling relation.21 From a practical viewpoint, this means that knowledge of the interparticle interactions and the thermodynamic excess entropy for these model fluids allows one to predict key aspects of their long-time dynamic behavior.
Can the GR scaling method be extended to treat systems with different types of dynamics (e.g., with equations of motion that incorporate dynamic effects of “solvent” surrounding the particles)? In this work, we take a step toward addressing this question. Specifically, we present simulation data that tests whether the GR scaling method can be recast in a manner useful for predicting long-time tracer diffusivities of suspended particles undergoing Brownian (i.e., overdamped Langevin) dynamics in the absence of interparticle hydrodynamic forces. The model fluids considered here comprise particles that interact via pair potentials
For the case of Brownian dynamics considered here, it is convenient to express the GR-reduced long-time diffusivity as21
where B is the second virial coefficient, β ≡ 1/kBT, kB is Boltzmann's constant, and D2 = −(∂[D/D0]/∂ρ)T|ρ = 0 quantifies how pair interactions at low particle concentration modify the tracer diffusivity D relative to the infinite dilution value D0. As has been shown elsewhere,21 (1 − D/D0)GR = −sex/kB for densities low enough that
Based on Eq. (1), it is evident that four quantities must be computed to test the GR excess-entropy scaling relation for state points across the T-ρ plane: B, D2, D/D0, and sex. The first two depend solely on β and
The physical significance of Q(r) is that it characterizes the extent to which application of a small external force f to a particle in the low-ρ fluid would distort its pair correlation function with its neighbors,
To compute D/D0, we use Brownian dynamics simulations, where particle translations are governed by the Langevin equation, solved in the overdamped limit. The algorithm we employ, wherein the position ri of each particle i ∈ [1, N] updates at each time step according to
The Brownian dynamics simulations track N particles in a periodically-replicated cubic simulation cell whose volume V sets ρ = N/V. We use N = 1000, 2000, and 3000 particles for the IPL, Hertzian, and Gaussian-core fluids, respectively. Interparticle potentials were truncated at r = 1.4 σ– 3.9 σ,62σ, and 3.71 σ, respectively, for these systems. For simulations of the IPL fluids, we set D0 = 0.001(σ2ε/m)1/2 and the Δt = 0.025τB, where τB ≡ mD0/kBT is a characteristic Brownian time scale. For simulations of Gaussian-core and Hertzian particle systems, we set D0 = 0.0001(σ2ε/m)1/2 and Δt = 0.1τB. Long-time tracer diffusivities are computed from the mean-squared displacements via the relation, D = 〈Δr2〉/6t as t → ∞.
We employ free-energy-based simulation techniques to evaluate sex. In the first step, we obtain the ρ dependence of the Helmholtz free energy at high T using grand canonical transition matrix Monte Carlo simulation.56 In the second step, we perform a canonical temperature expanded ensemble simulation57 with a transition matrix Monte Carlo algorithm58 to evaluate the change in Helmholtz free energy with T at constant ρ. Collectively, these simulations provide values for excess Helmholtz free energy and excess energy, and hence sex, at a state point of interest. Interested readers can find additional details in our earlier papers.59,60
We now examine the simulated behavior of six IPL fluids with exponents of μ = 6, 8, 10, 12, 18, and 36. Specifically, Fig. 1(a) shows how (1 − D/D0) and its GR-reduced counterpart (1 − D/D0)GR vary with density ρσ3 and −sex/kB, respectively, for these fluids at T = ε/kB. As should be expected, systems with softer interparticle repulsions (i.e., lower μ) display considerably faster single-particle dynamics at a given ρσ3. For example, the μ = 6 fluid must realize a density twice that of the μ = 36 fluid for both systems to exhibit D/D0 = 0.2. In contrast, (1 − D/D0)GR for these fluids shows only minor systematic variation when plotted versus −sex/kB, with less than 20% difference between the μ = 6 and μ = 36 values.
To put to the GR scaling relationship for Brownian dynamics to a more stringent test, we further analyze the behaviors of model fluids with ultrasoft Hertzian and Gaussian-core interactions [see Figs. 1(b) and 1(c), respectively]. In particular, we explore a broad range of ρσ3 for these two systems along isotherms spanning from the moderately supercooled liquid (T ⩽ 0.004ε/kB) to temperatures high enough (T ⩾ 0.01ε/kB) for the models to exhibit mean-field-like behavior.4,5,7,26 Note that this data set includes state points for which D/D0 increases with increasing ρσ3 in both systems, behavior that is anomalous when compared to the trends exhibited by simpler atomic and colloidal fluids. The main conclusion to be drawn from Figs. 1(b) and 1(c) is that, despite the anomalous relationship between single-particle dynamics and ρ in these systems, the GR excess entropy scaling relation looks strikingly similar to the quasi-universal form of the IPL fluids. In fact, the Hertzian and Gaussian-core data deviate by less than 20% from the least-squares fit to the (1 − D/D0)GR versus sex data from Fig. 1(a).
One consequence of this quasi-universal scaling behavior is that one can estimate Brownian dynamics D/D0 data, even for nontrivial systems with ultrasoft potentials, simply from knowledge of their interparticle interactions, the excess entropy, and, e.g., a least-squares fit of the IPL data of Fig. 1(a). Figure 2 illustrates the quality of the predictions made by this approach for the Hertzian and Gaussian-core model fluids. As can be seen—with no adjustable parameters—the estimated values of diffusivity are accurate for low-to-moderate ρ and remain qualitatively reliable for higher ρ where the anomalous trends emerge. Collectively, these results illustrate a strong apparent link between static interparticle correlations and long-time dynamics for model fluids with either hard or ultrasoft interactions.
Finally, whether the approach outlined here can provide an accurate estimate to the “thermodynamic” component of diffusivity for systems which also have strong hydrodynamic interparticle contributions to dynamics is an open question, which deserves attention in future studies. Furthermore, it will be interesting to test whether the GR scaling for diffusivity in Brownian dynamics also holds for systems with more complex interparticle interactions (e.g., with more than one characteristic length scale).
We thank Professor David S. Dean for suggesting the relaxation method for computing D2. T.M.T. acknowledges support of the Welch Foundation (F-1696) and the David and Lucile Packard Foundation. J.R.E. acknowledges financial support of the National Science Foundation, Grant No. CBET-0828979. M.J.P. acknowledges the support of the Thrust 2000—Harry P. Whitworth Endowed Graduate Fellowship in Engineering. The Texas Advanced Computing Center (TACC), the University at Buffalo Center for Computational Research, and the Rensselaer Polytechnic Institute Computational Center for Nanotechnology Innovations provided computational resources for this study.