Nicoletta Gnan

Pressure-energy correlations in liquids. I. Results from computer simulations

We show that a number of model liquids at fixed volume exhibit strong correlations between equilibrium fluctuations of the configurational parts of (instantaneous) pressure and energy. We present detailed results for 13 systems, showing in which systems these correlations are significant. These include Lennard-Jones liquids (both single- and two-component) and several other simple liquids, neither hydrogen-bonding liquids such as methanol and water, nor the Dzugutov liquid, which has significant contributions to pressure at the second nearest neighbor distance. The pressure-energy correlations, which for the Lennard-Jones case are shown to also be present in the crystal and glass phases, reflect an effective inverse power-law potential dominating fluctuations, even at zero and slightly negative pressure. An exception to the inverse power-law explanation is a liquid with hard-sphere repulsion and a square-well attractive part, where a strong correlation is observed, but only after time averaging. The companion paper [N. P. Bailey et al., J. Chem. Phys. 129, 184508 (2008)] gives a thorough analysis of the correlations, with a focus on the Lennard-Jones liquid, and a discussion of some experimental and theoretical consequences.

Pressure-energy correlations in liquids. II. Analysis and consequences

We present a detailed analysis and discuss consequences of the strong correlations of the configurational parts of pressure and energy in their equilibrium fluctuations at fixed volume reported for simulations of several liquids in the previous paper [N. P. Bailey et al., J. Chem. Phys. 129, 184507 (2008)]. The analysis concentrates specifically on the single-component Lennard-Jones system. We demonstrate that the potential may be replaced, at fixed volume, by an effective power law but not simply because only short-distance encounters dominate the fluctuations. Indeed, contributions to the fluctuations are associated with the whole first peak of the radial distribution function, as we demonstrate by an eigenvector analysis of the spatially resolved covariance matrix. The reason the effective power law works so well depends crucially on going beyond single-pair effects and on the constraint of fixed volume. In particular, a better approximation to the potential includes a linear term, which contributes to the mean values of potential energy and virial, but little to their fluctuations, for density fluctuations which conserve volume. We also study in detail the zero temperature limit of the (classical) crystalline phase, where the correlation coefficient becomes very close, but not equal, to unity, in more than one dimension; in one dimension the limiting value is exactly unity. In the second half of the paper we consider four consequences of strong pressure-energy correlations: (1) analyzing experimental data for supercritical argon we find 96% correlation; (2) we discuss the particular significance acquired by the correlations for viscous van der Waals liquids approaching the glass transition: For strongly correlating viscous liquids knowledge of just one of the eight frequency-dependent thermoviscoelastic response functions basically implies knowledge of them all; (3) we reinterpret aging simulations of ortho-terphenyl carried out by Mossa et al. [Eur. Phys. J. B 30, 351 (2002)], showing their conclusions follow from the strongly correlating property; and (4) we briefly discuss the presence of the correlations (after appropriate time averaging) in model biomembranes, showing that significant correlations may be present even in quite complex systems.

Pressure-energy correlations in liquids. III. Statistical mechanics and thermodynamics of liquids with hidden scale invariance

Computer simulations recently revealed that several liquids exhibit strong correlations between virial and potential energy equilibrium fluctuations in the NVT ensemble (U. R. Pedersen et al., Phys. Rev. Lett. 100, 015701 (2008)). In order to investigate whether these correlations are present also far from equilibrium constant-volume aging following a temperature down jump from equilibrium was simulated for two strongly correlating liquids, an asymmetric dumbbell model and Lewis-Wahnstrom OTP, as well as for SPC water that is not strongly correlating. For the two strongly correlating liquids virial and potential energy follow each other closely during the aging towards equilibrium. For SPC water, on the other hand, virial and potential energy vary with little correlation as the system ages towards equilibrium. Further proof that strong pressure-energy correlations express a configuration space property comes from monitoring pressure and energy during the crystallization (reported here for the first time) of supercooled Lewis-Wahnstrom OTP at constant temperature.

Multidimensional stationary probability distribution for interacting active particles

We derive the stationary probability distribution for a non-equilibrium system composed by an arbitrary number of degrees of freedom that are subject to Gaussian colored noise and a conservative potential. This is based on a multidimensional version of the Unified Colored Noise Approximation. By comparing theory with numerical simulations we demonstrate that the theoretical probability density quantitatively describes the accumulation of active particles around repulsive obstacles. In particular, for two particles with repulsive interactions, the probability of close contact decreases when one of the two particle is pinned. Moreover, in the case of isotropic confining potentials, the radial density profile shows a non trivial scaling with radius. Finally we show that the theory well approximates the “pressure” generated by the active particles allowing to derive an equation of state for a system of non-interacting colored noise-driven particles.

Predicting the Effective Temperature of a Glass

We explain the findings by Di Leonardo et al. [Phys. Rev. Lett. 84, 6054 (2000)10.1103/PhysRevLett.84.6054] that the effective temperature of a Lennard-Jones glass depends only on the final density in the volume and/or temperature jump that produces the glass. This is not only a property of the Lennard-Jones liquid, but a feature of all strongly correlating liquids. For such liquids data from a single quench simulation provide enough information to predict the effective temperature of any glass produced by jumping from an equilibrium state. This prediction is validated by simulations of the Kob-Andersen binary Lennard-Jones liquid and shown not to apply for the nonstrongly correlating monatomic Lennard-Jones Gaussian liquid.

Velocity distribution in active particles systems.

We derive an analytic expression for the distribution of velocities of multiple interacting active particles which we test by numerical simulations. In clear contrast with equilibrium we find that the velocities are coupled to positions. Our model shows that, even for two particles only, the individual velocities display a variance depending on the interparticle separation and the emergence of correlations between the velocities of the particles. When considering systems composed of many particles we find an analytic expression connecting the overall velocity variance to density, at the mean-field level, and to the pair distribution function valid in the limit of small noise correlation times. Finally we discuss the intriguing analogies and main differences between our effective free energy functional and the theoretical scenario proposed so far for phase-separating active particles.

Erratum: “Pressure-energy correlations in liquids. I. Results from computer simulations” [J. Chem. Phys. 129, 184507 (2008)]

1 On page 1 Ref. 1 three lines above Eq. 1 , the phrase in parentheses “ the parts in addition to the ideal gas terms ” should be moved up one line and set as end to the previous sentence, which should be as follows: “Specifically, the fluctuations that are in many cases strongly correlated are those of the configurational parts of pressure and energy the parts in addition to the ideal gas terms .” 2 On page 2, the left hand side of Eq. 3 should be the potential energy not the virial, i.e., the term “Wpair” should be replaced by “Upair.”

Discontinous change from thermally- to geometrically-dominated effective interactions in colloidal solutions

We report numerical results for the effective potential arising between two colloids immersed in a self-assembling cosolute which forms reversible clusters. The potential is evaluated at cosolute state points with different densities and temperatures but with the same connectivity properties. We find that the range of the resulting effective potential is controlled only by the cosolute thermal correlation length rather than by its connectivity length. We discuss the significant differences from previous results focusing on the cosolute forming irreversible clusters and we show that the irreversible bond case represents a singular limit which cannot be accessed in equilibrium by continuously increasing the bond lifetime.

Corrigendum: Velocity distribution in active particles systems

Scientific Reports 6: Article number: 2329710.1038/srep23297; published online: March222016; updated: May312016 In this Article, Roberto Di Leonardo is incorrectly affiliated with ‘CNR-IMIP, UOS Roma, Dipartimento di Fisica Universita Sapienza, I-00185, Roma, Italy’. The correct affiliation is listed below: NANOTEC-CNR, Institute of Nanotechnology, Soft and Living Matter Laboratory, Roma, Italy.

Modelling realistic microgels in an explicit solvent

Thermoresponsive microgels are polymeric colloidal networks that can change their size in response to a temperature variation. This peculiar feature is driven by the nature of the solvent-polymer interactions, which triggers the so-called volume phase transition from a swollen to a collapsed state above a characteristic temperature. Recently, an advanced modelling protocol to assemble realistic, disordered microgels has been shown to reproduce experimental swelling behavior and form factors. In the original framework, the solvent was taken into account in an implicit way, condensing solvent-polymer interactions in an effective attraction between monomers. To go one step further, in this work we perform simulations of realistic microgels in an explicit solvent. We identify a suitable model which fully captures the main features of the implicit model and further provides information on the solvent uptake by the interior of the microgel network and on its role in the collapse kinetics. These results pave the way for addressing problems where solvent effects are dominant, such as the case of microgels at liquid-liquid interfaces.