How short-range attractions impact the structural order, self-diffusivity, and viscosity of a fluid

Self-Diffusion in Liquid Copper, Silver, and Gold

The recently developed by us semi-analytical representation of the mean spherical approximation in conjunction with the linear trajectory approximation is applied to the quantitative study of self-diffusivities in liquid Cu, Ag and Au at different temperatures. The square-well model is employed for the description of the interatomic pair interactions in metals under study. It is found that our theoretical results are in good agreement with available experimental and computer-simulation data and can be considered as a prediction when such data are absent.

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.

Self-diffusion coefficient of the square-well fluid from molecular dynamics simulations within the constant force approach

We present a systematic study of the self-diffusion coefficient for a fluid of particles interacting via the square-well pair potential by means of molecular dynamics simulations in the canonical (N, V, T) ensemble. The discrete nature of the interaction potential is modeled by the constant force approximation, and the self-diffusion coefficient is determined for several fluid densities at supercritical thermodynamic states. The dependence of the self-diffusion coefficient on the potential range λ is analyzed in the range of 1.1 ≤ λ ≤ 1.5. The obtained simulation results are in agreement with the self-diffusion coefficient predicted by the Enskog method. Additionally, we show that the diffusion coefficient is very sensitive to the potential range λ. Our results for the self-diffusion coefficient times density extrapolate well to the values in the zero-density limit obtained from the Chapman-Enskog theory for dilute gases. The constant force approximation used in this work to model the discrete pair potentials has shown to be an excellent scheme to compute the transport properties of square-well fluids using molecular dynamics simulations. Finally, the simulation results presented here are useful for improving theoretical approaches, such as the Enskog method.

Relationship between structural order and water-like anomalies in metastable liquid silicon: Ab initio molecular dynamics

The relationship between structural order and water-like anomalies in tetrahedral liquids is still open. Here, first-principle molecular dynamics are performed to study it in metastable liquid Si. It is found that in T-P phase diagram, there indeed exists a structural anomaly region, which encloses density anomaly but not diffusivity anomaly. This is consistent with that of SW Si and BKS SiO₂ but different from that of SPC/E water. Two-body excess entropy anomaly can neither capture the diffusivity, structural, and density anomalies, as it can in a two-scale potential fluid. In structural anomaly region, tetrahedrality order q_tetra (measuring the extent to which an atom and its four nearest neighbours adopt tetrahedral arrangement) and translational order t_trans (measuring the tendency of two atoms to adopt preferential separation) are not perfectly correlated, which is different from that in SW Si and renders it impossible to use the isotaxis line to quantify the degree of structural order needed for water-like anomalies to occur. Along the isotherm of critical temperature T_c, t_trans/q_tetra is approximately linear with pressure. With decreasing pressure along the isotherm below T_c, t_trans/q_tetra departs downward from the line, while it is the opposite case above T_c.

Impact of solvent granularity and layering on tracer hydrodynamics in confinement

Classic hydrodynamic arguments establish that when a spherical tracer particle is suspended between parallel walls, tracer-wall coupling mediated by the solvent will cause the tracer to exhibit position-dependent diffusivity. We investigate how the diffusivity profiles of confined tracers are impacted by the diameter size-ratio of the tracer to solvent: starting from the classic limit of infinite size-ratio (i.e., continuum solvent), we consider size-ratios of four or less to examine how hydrodynamic predictions are disrupted for systems where the tracer and solvent are of similar scale. We use computer simulations and techniques based on the Fokker-Planck formalism to calculate the diffusivity profiles of hard-sphere tracer particles in hard-sphere solvents, focusing on the dynamics perpendicular to the walls. Given wall separations of several tracer diameters, we first consider confinement between hard walls, where anisotropic structuring at the solvent lengthscale generates inhomogeneity in the tracer free-energy landscape and undermines hydrodynamic predictions locally. We then introduce confining planes that we term transparent walls, which restrict tracer and solvent center-accessibilities while completely eliminating static anisotropy, and reveal position-dependent signatures in tracer diffusivity solely attributable to confinement. With or without suppressing static heterogeneity, we find that tracer diffusivity increasingly deviates on a local basis from hydrodynamic predictions at smaller size-ratios. However, hydrodynamic theory still approximately captures spatially-averaged dynamics across the pores even for very small tracer-solvent size-ratios over a wide range of solvent densities and wall separations.

Simple liquids' quasiuniversality and the hard-sphere paradigm

This topical review discusses the quasiuniversality of simple liquids' structure and dynamics and two possible justifications of it. The traditional one is based on the van der Waals picture of liquids in which the hard-sphere system reflects the basic physics. An alternative explanation argues that all quasiuniversal liquids to a good approximation conform to the same equation of motion, referring to the exponentially repulsive pair-potential system as the basic reference system. The paper, which is aimed at non-experts, ends by listing a number of open problems in the field.

How Local and Average Particle Diffusivities of Inhomogeneous Fluids Depend on Microscopic Dynamics

Computer simulations and a stochastic Fokker-Planck equation based approach are used to compare the single-particle diffusion coefficients of equilibrium hard-sphere fluids exhibiting identical inhomogeneous static structure and governed by either Brownian (i.e., overdamped Langevin) or Newtonian microscopic dynamics. The physics of inhomogeneity is explored via the imposition of one-dimensional sinusoidal density profiles of different wavelengths and amplitudes. When imposed density variations are small in magnitude for distances on the scale of a particle diameter, bulk-like average correlations between local structure and mobility are observed. In contrast, when density variations are significant on that length scale, qualitatively different structure-mobility correlations emerge that are sensitive to the governing microscopic dynamics. Correspondingly, a previously proposed scaling between long-time diffusivities for bulk isotropic fluids of particles exhibiting Brownian versus Newtonian dynamics [Pond et al. Soft Matter 2011, 7, 9859-9862] cannot be generalized to describe the position-dependent behaviors of strongly inhomogeneous fluids. While average diffusivities in the inhomogeneous and homogeneous directions are coupled, their qualitative dependencies on inhomogeneity wavelength are sensitive to the details of the microscopic dynamics. Nonetheless, average diffusivities of the inhomogeneous fluids can be approximately predicted for either type of dynamics based on knowledge of bulk isotropic fluid behavior and how inhomogeneity modifies the distribution of available volume. Analogous predictions for average diffusivities of experimental, inhomogeneous colloidal dispersions (based on known bulk behavior) suggest that they will exhibit qualitatively different trends than those predicted by models governed by overdamped Langevin dynamics that do not account for hydrodynamic interactions.

Predicting the structure of fluids with piecewise constant interactions: Comparing the accuracy of five efficient integral equation theories

We use molecular dynamics simulations to test integral equation theory predictions for the structure of fluids of spherical particles with eight different piecewise-constant pair-interaction forms comprising a hard core and a combination of two shoulders and/or wells. Since model pair potentials like these are of interest for discretized or coarse-grained representations of effective interactions in complex fluids (e.g., for computationally intensive inverse optimization problems), we focus here on assessing how accurately their properties can be predicted by analytical or simple numerical closures including Percus-Yevick, hypernetted-chain, and reference hypernetted-chain closures and first-order mean spherical and modified first-order mean spherical approximations. To make quantitative comparisons between the predicted and simulated radial distribution functions, we introduce a cumulative structural error metric. For equilibrium fluid state points of these models, we find that the reference hypernetted-chain closure is the most accurate of the tested approximations as characterized by this metric or related thermodynamic quantities.

Flow enhancement of water-based nanoparticle dispersion through microscale sedimentary rocks

Understanding and controlling fluids flow at the microscale is a matter of growing scientific and technological interest. Flow enhancements of water-based nanoparticle dispersions through microscale porous media are investigated through twelve hydrophilic sedimentary rocks with pore-throat radius between 1.2 and 10 μm, which are quantitatively explained with a simple model with slip length correction for Darcy flow. Both as wetting phase, water exhibited no-slip Darcy flow in all cores; however, flow enhancement of nanoparticle dispersions can be up to 5.7 times larger than that of water, and it increases with the decreasing of pore-throat radius. The experimental data reveals characteristic slip lengths are of order 500 and 1000 nm for 3M® and HNPs-1 nanoparticles, respectively, independent of the lithology or nanoparticle concentration or shear rate. Meanwhile, the phenomenon of flow degradation is observed for HNPs-2 nanoparticles. These results explore the feasible application of using nanoparticle dispersions to control flow at the microscale.

The impact range for smooth wall-liquid interactions in nanoconfined liquids

Bulk and nanoconfined liquids have very different physics; for instance, nanoconfined liquids show stratification and position-dependent relaxation processes. A number of similarities between bulk and nanoconfined liquids have nevertheless been reported in computer simulations during the last decade. Inspired by these observations, we present results from molecular dynamics computer simulations of four nanoconfined liquids (the single-component Lennard-Jones liquid, the Kob-Andersen binary Lennard-Jones mixture, an asymmetric dumbbell model, and the Dzugutov liquid) demonstrating also a microscopic similarity between bulk and nanoconfined liquids. The results show that the interaction range for the wall-liquid and liquid-liquid interactions of the nanoconfined liquid is identical to that of the bulk liquid if the liquid is "Roskilde simple" in bulk as well as nanoconfinement, i.e., exhibits strong correlations between virial and potential-energy equilibrium fluctuations in the NVT ensemble. Under this condition, interactions beyond the first coordination shell may be ignored, in particular for the wall-liquid interactions. This is shown not to be the case for non-Roskilde-simple liquids.

Aggregation dynamics, structure, and mechanical properties of bigels

Recently we have introduced bigels, inter-penetrating gels made of two different colloidal species. Even if particles with simple short-range isotropic potential are employed, the selective interactions enable the tunability of the self-assembly, leading to the formation of complex structures. In the present paper, we explore the non-equilibrium dynamics and the phenomenology underlying the kinetic arrest under quench and the formation of bigels. We demonstrate that the peculiar bigel kinetics can be described through an arrested spinodal decomposition driven by demixing of the colloidal species. The role played by the presence of a second colloidal species on the phase diagram, as expanded to account for the increased number of parameters, is clarified both via extensive numerical simulations and experiments. We provide details on the realisation of bigels, by means of DNA-coated colloids (DNACCs), and the consequent imaging techniques. Moreover we evidence, by comparison with the usual one-component gel formation, the emergence of controllable timescales in the aggregation of the bigels, whose final stages are also experimentally studied to provide morphological details. Finally, we use numerical models to simulate the bigel response to mechanical strain, highlighting how such a new material can bear significantly higher stress compared to the usual one-component gel. We conclude by discussing possible technological uses and by providing insights on the viable research steps to undertake for more complex and yet tuneable multi-component colloidal systems.

Arrested demixing opens route to bigels

Understanding and, ultimately, controlling the properties of amorphous materials is one of the key goals of material science. Among the different amorphous structures, a very important role is played by colloidal gels. It has been only recently understood that colloidal gels are the result of the interplay between phase separation and arrest. When short-ranged attractive colloids are quenched into the phase-separating region, density fluctuations are arrested and this results in ramified amorphous space-spanning structures that are capable of sustaining mechanical stress. We present a mechanism of aggregation through arrested demixing in binary colloidal mixtures, which leads to the formation of a yet unexplored class of materials--bigels. This material is obtained by tuning interspecies interactions. Using a computer model, we investigate the phase behavior and the structural properties of these bigels. We show the topological similarities and the geometrical differences between these binary, interpenetrating, arrested structures and their well-known monodisperse counterparts, colloidal gels. Our findings are supported by confocal microscopy experiments performed on mixtures of DNA-coated colloids. The mechanism of bigel formation is a generalization of arrested phase separation and is therefore universal.

Transport coefficients in silicate melts from structural data via a structure-thermodynamics-dynamics relationship

The viscosity and diffusivities of silicate melts under high-pressure, high-temperature conditions are difficult to obtain experimentally. Estimation and extrapolation of transport coefficients are further complicated by their extreme sensitivity to melt composition. Our molecular-dynamics simulations show that, over a broad range of melt composition, temperature, and pressure, the diffusivities correlate with the excess entropy; approximations to the latter can be obtained from the knowledge of the radial distribution function. Using this structure-thermodynamics-dynamics relationship, we show that transport properties of silicate melts can be estimated quantitatively using static structure factor data from experiments.

Anomalous change in the dynamics of a supercritical fluid

We perform molecular dynamics simulations to investigate dynamical properties of a supercritical Lennard-Jones fluid. We find that in the supercritical region there is a short-ranged deviation in dynamic character. We further find that this anomalous change is associated with the presence of the Widom line, the locus of specific heat maxima, of the liquid-vapor phase transition. The salient change in dynamics is consistent with a crossover in the correlation of the diffusion coefficient with the excess entropy. Our results lead to an interpretation that, even though a supercritical fluid excludes a singularity, its dynamical properties can be significantly affected by the existence of thermodynamic response maxima.

Testing "microscopic" theories of glass-forming liquids

We assess the validity of "microscopic" approaches of glass-forming liquids based on the sole knowledge of the static pair density correlations. To do so, we apply them to a benchmark provided by two liquid models that share very similar static pair density correlation functions while displaying distinct temperature evolutions of their relaxation times. We find that the approaches are unsuccessful in describing the difference in the dynamical behavior of the two models. Our study is not exhaustive, and we have not tested the effect of adding corrections by including, for instance, three-body density correlations. Yet, our results appear strong enough to challenge the claim that the slowdown of relaxation in glass-forming liquids, for which it is well established that the changes of the static structure factor with temperature are small, can be explained by "microscopic" approaches only requiring the static pair density correlations as nontrivial input.

Excess entropy scaling of dynamic quantities for fluids of dumbbell-shaped particles

We use molecular simulation to study the ability of entropy scaling relationships to describe the kinetic properties of two Lennard-Jones dumbbell models. We begin by examining the excess entropy, the key quantity used to correlate dynamic properties within entropy scaling strategies. We compute the thermodynamic excess entropy as well as contributions to the two-body excess entropy stemming from translational and orientational intermolecular correlations. Our results indicate that the total two-body contribution accounts for more than 70% of the thermodynamic excess entropy at all state conditions explored. For the two dumbbell models studied here, the orientational component of the two-body excess entropy dominates at moderate and high fluid densities. We next investigate the relationships between kinetic properties and various contributions to the excess entropy. 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 find that the thermodynamic excess entropy provides the best metric for describing kinetic properties. For each of the dynamic properties considered, reduced data collapse onto a common curve when expressed as a function of the thermodynamic excess entropy. The likelihood of a two-body contribution to the excess entropy serving as a reliable scaling variable is linked to the extent to which it correlates with the thermodynamic excess entropy. The total two-body term contributes significantly to the excess entropy, and therefore this quantity generally serves as a suitable scaling variable.

Generalized Rosenfeld scalings for tracer diffusivities in not-so-simple fluids: mixtures and soft particles

Rosenfeld [Phys. Rev. A 15, 2545 (1977)] originally noticed that casting the transport coefficients of simple monatomic equilibrium fluids in a specific dimensionless form makes them approximately single-valued functions of excess entropy. This observation has predictive value because, while the transport coefficients of dense fluids can be difficult to estimate from first principles, the excess entropy can often be accurately predicted from liquid-state theory. In this work, we use molecular simulations to investigate whether Rosenfeld's observation is a special case of a more general scaling law relating the tracer diffusivities of particles in mixtures to the excess entropy. Specifically, we study the tracer diffusivities, static structure, and thermodynamic properties of a variety of one- and two-component model fluid systems with either additive or nonadditive interactions of the hard-sphere or Gaussian-core form. The results of the simulations demonstrate that the effects of mixture concentration and composition, particle-size asymmetry and additivity, and strength of the interparticle interactions in these fluids are consistent with an empirical scaling law relating the excess entropy to a dimensionless (generalized Rosenfeld) form of tracer diffusivity, which we introduce here. The dimensionless form of the tracer diffusivity follows from knowledge of the intermolecular potential and the transport/thermodynamic behavior of fluids in the dilute limit. The generalized Rosenfeld scaling requires less information and provides more accurate predictions than either Enskog theory or scalings based on the pair-correlation contribution to the excess entropy. As we show, however, it also suffers from some limitations especially for systems that exhibit significant decoupling of individual component tracer diffusivities.

Phase diagram of Hertzian spheres

We report the phase diagram of interpenetrating Hertzian spheres. The Hertz potential is purely repulsive, bounded at zero separation, and decreases monotonically as a power law with exponent 5/2, vanishing at the overlapping threshold. This simple functional describes the elastic interaction of weakly deformable bodies and, therefore, it is a reliable physical model of soft macromolecules, like star polymers and globular micelles. Using thermodynamic integration and extensive Monte Carlo simulations, we computed accurate free energies of the fluid phase and a large number of crystal structures. For this, we defined a general primitive unit cell that allows for the simulation of any lattice. We found multiple re-entrant melting and first-order transitions between crystals with cubic, trigonal, tetragonal, and hexagonal symmetries.

Anomalous structure and dynamics of the Gaussian-core fluid

It is known that there are thermodynamic states for which the Gaussian-core fluid displays anomalous properties such as expansion upon isobaric cooling (density anomaly) and increased single-particle mobility upon isothermal compression (self-diffusivity anomaly). Here, we investigate how temperature and density affect its short-range translational structural order, as characterized by the two-body excess entropy. We find that there is a wide range of conditions for which the short-range translational order of the Gaussian-core fluid decreases upon isothermal compression (structural order anomaly). As we show, the origin of the structural anomaly is qualitatively similar to that of other anomalous fluids (e.g., water or colloids with short-range attractions) and is connected to how compression affects static correlations at different length scales. Interestingly, we find that the self-diffusivity of the Gaussian-core fluid obeys a scaling relationship with the two-body excess entropy that is very similar to the one observed for a variety of simple liquids. One consequence of this relationship is that the state points for which structural, self-diffusivity, and density anomalies of the Gaussian-core fluid occur appear as cascading regions on the temperature-density plane; a phenomenon observed earlier for models of waterlike fluids. There are, however, key differences between the anomalies of Gaussian-core and waterlike fluids, and we discuss how those can be qualitatively understood by considering the respective interparticle potentials of these models. Finally, we note that the self-diffusivity of the Gaussian-core fluid obeys different scaling laws depending on whether the two-body or total excess entropy is considered. This finding, which deserves more comprehensive future study, appears to underscore the significance of higher-body correlations for the behavior of fluids with bounded interactions.

Relation of water anomalies to the excess entropy

Using the five-site transferable intermolecular potential (TIP5P) we perform molecular dynamics simulations to investigate the relationship between the excess entropy and anomalies of water. We find that the two-body excess entropy is an ideal quantity to predict the regions of structural, dynamic, and thermodynamic anomalies of water in its pressure-temperature and density-temperature phase diagrams. From the excess entropy we can also predict the location of the Widom line, associated with the hypothesized liquid-liquid critical point in supercooled water.

Shear-rate-dependent structural order and viscosity of a fluid with short-range attractions

We study a model short-range attractive fluid under shear. For this system, the strength of interparticle attractions strongly influences the equilibrium structural order. We find that shear monotonically decreases structural order regardless of the strength of the attractions. There is a strong correlation between shear-rate-dependent viscosity and a structural order metric, suggesting a structurally based constitutive equation. This correlation also holds for the Lennard-Jones fluid.

Large-amplitude jumps and non-Gaussian dynamics in highly concentrated hard sphere fluids

Our microscopic stochastic nonlinear Langevin equation theory of activated dynamics has been employed to study the real-space van Hove function of dense hard sphere fluids and suspensions. At very short times, the van Hove function is a narrow Gaussian. At sufficiently high volume fractions, such that the entropic barrier to relaxation is greater than the thermal energy, its functional form evolves with time to include a rapidly decaying component at small displacements and a long-range exponential tail. The "jump" or decay length scale associated with the tail increases with time (or particle root-mean-square displacement) at fixed volume fraction, and with volume fraction at the mean alpha relaxation time. The jump length at the alpha relaxation time is predicted to be proportional to a measure of the decoupling of self-diffusion and structural relaxation. At long times corresponding to mean displacements of order a particle diameter, the volume fraction dependence of the decay length disappears. A good superposition of the exponential tail feature based on the jump length as a scaling variable is predicted at high volume fractions. Overall, the theoretical results are in good accord with recent simulations and experiments. The basic aspects of the theory are also compared with a classic jump model and a dynamically facilitated continuous time random-walk model. Decoupling of the time scales of different parts of the relaxation process predicted by the theory is qualitatively similar to facilitated dynamics models based on the concept of persistence and exchange times if the elementary event is assumed to be associated with transport on a length scale significantly smaller than the particle size.

Structural anomalies of fluids: origins in second and higher coordination shells

Compressing or cooling a fluid typically enhances its static interparticle correlations. However, there are notable exceptions. Isothermal compression can reduce the translational order of fluids that exhibit anomalous waterlike trends in their thermodynamic and transport properties, while isochoric cooling (or strengthening of attractive interactions) can have a similar effect on fluids of particles with short-range attractions. Recent simulation studies by Yan [Phys. Rev. E 76, 051201 (2007)] on the former type of system and Krekelberg [J. Chem. Phys. 127, 044502 (2007)] on the latter provide examples where such structural anomalies can be related to specific changes in second and more distant coordination shells of the radial distribution function. Here, we confirm the generality of this microscopic picture through analysis, via molecular simulation and integral equation theory, of coordination shell contributions to the two-body excess entropy for several related model fluids which incorporate different levels of molecular resolution. The results suggest that integral equation theory can be an effective and computationally inexpensive tool for assessing, based on the pair potential alone, whether new model systems are good candidates for exhibiting structural (and hence thermodynamic and transport) anomalies.

Tuning density profiles and mobility of inhomogeneous fluids

Density profiles are the most common measure of inhomogeneous structure in confined fluids, but their connection to transport coefficients is poorly understood. We explore via simulation how tuning particle-wall interactions to flatten or enhance the particle layering of a model confined fluid impacts its self-diffusivity, viscosity, and entropy. Interestingly, interactions that eliminate particle layering significantly reduce confined fluid mobility, whereas those that enhance layering can have the opposite effect. Excess entropy helps to understand and predict these trends.

Glassy Arrest in Colloidal Fluids with Size Polydispersity

Short-range attraction between colloidal particles such as proteins can drive a glass-like structural arrest. For monodisperse systems, mode-coupling theory affords a simple asymptotic prediction of the transition. Here, using a depletion mapping framework, we extend this result to incorporate size polydispersity. For comparison, we also give an energy landscape formulation of the transition. We comment on the relevance to subcellular crowding, recombinant protein expression, and osmotic stress in microbial organisms.