S. V. Buldyrev

Liquid-liquid phase transitions for soft-core attractive potentials.

Using event-driven molecular dynamics simulations, we study a three-dimensional one-component system of spherical particles interacting via a discontinuous potential combining a repulsive square soft core and an attractive square well. In the case of a narrow attractive well, it has been shown that this potential has two metastable gas-liquid critical points. Here we systematically investigate how the changes of the parameters of this potential affect the phase diagram of the system. We find a broad range of potential parameters for which the system has both a gas-liquid critical point C1 and a liquid-liquid critical point C2. For the liquid-gas critical point we find that the derivatives of the critical temperature and pressure, with respect to the parameters of the potential, have the same signs: they are positive for increasing width of the attractive well and negative for increasing width and repulsive energy of the soft core. This result resembles the behavior of the liquid-gas critical point for standard liquids. In contrast, for the liquid-liquid critical point the critical pressure decreases as the critical temperature increases. As a consequence, the liquid-liquid critical point exists at positive pressures only in a finite range of parameters. We present a modified van der Waals equation which qualitatively reproduces the behavior of both critical points within some range of parameters, and gives us insight on the mechanisms ruling the dependence of the two critical points on the potentials parameters. The soft-core potential studied here resembles model potentials used for colloids, proteins, and potentials that have been related to liquid metals, raising an interesting possibility that a liquid-liquid phase transition may be present in some systems where it has not yet been observed.

Nanoscale Dynamics of Phase Flipping in Water near its Hypothesized Liquid-Liquid Critical Point

One hypothesized explanation for waters anomalies imagines the existence of a liquid-liquid (LL) phase transition line separating two liquid phases and terminating at a LL critical point. We simulate the classic ST2 model of water for times up to 1000 ns and system size up to N = 729. We find that for state points near the LL transition line, the entire system flips rapidly between liquid states of high and low density. Our finite-size scaling analysis accurately locates both the LL transition line and its associated LL critical point. We test the stability of the two liquids with respect to the crystal and find that of the 350 systems simulated, only 3 of them crystallize and these 3 for the relatively small system size N = 343 while for all other simulations the incipient crystallites vanish on a time scales smaller than ≈ 100 ns.

Metastable liquid-liquid phase transition in a single-component system with only one crystal phase and no density anomaly

We investigate the phase behavior of a single-component system in three dimensions with spherically-symmetric, pairwise-additive, soft-core interactions with an attractive well at a long distance, a repulsive soft-core shoulder at an intermediate distance, and a hard-core repulsion at a short distance, similar to potentials used to describe liquid systems such as colloids, protein solutions, or liquid metals. We showed [Nature (London) 409, 692 (2001)] that, even with no evidence of the density anomaly, the phase diagram has two first-order fluid-fluid phase transitions, one ending in a gas-low-density-liquid (LDL) critical point, and the other in a gas-high-density-liquid (HDL) critical point, with a LDL-HDL phase transition at low temperatures. Here we use integral equation calculations to explore the three-parameter space of the soft-core potential and perform molecular dynamics simulations in the interesting region of parameters. For the equilibrium phase diagram, we analyze the structure of the crystal phase and find that, within the considered range of densities, the structure is independent of the density. Then, we analyze in detail the fluid metastable phases and, by explicit thermodynamic calculation in the supercooled phase, we show the absence of the density anomaly. We suggest that this absence is related to the presence of only one stable crystal structure.

Discrete molecular dynamics simulations of peptide aggregation

We study the aggregation of peptides using the discrete molecular dynamics simulations. Specifically, at temperatures above the alpha-helix melting temperature of a single peptide, the model peptides aggregate into a multilayer parallel beta-sheet structure. This structure has an interstrand distance of 4.8 A and an intersheet distance of 10 A, which agree with experimental observations. Our model explains these results as follows: hydrogen-bond interactions give rise to the interstrand spacing in beta sheets, while Gō interactions between side chains make beta strands parallel to each other and allow beta sheets to pack into layers. An important feature of our results is that the aggregates contain free edges, which may allow for further aggregation of model peptides to form elongated fibrils.

Models for a liquid–liquid phase transition

We use molecular dynamics simulations to study two- and three-dimensional models with the isotropic double-step potential which in addition to the hard core has a repulsive soft core of larger radius. Our results indicate that the presence of two characteristic repulsive distances (hard core and soft core) is sufficient to explain liquid anomalies and a liquid–liquid phase transition, but these two phenomena may occur independently. Thus liquid–liquid transitions may exist in systems like liquid metals, regardless of the presence of the density anomaly. For 2D, we propose a model with a specific set of hard core and soft core parameters, that qualitatively reproduces the phase diagram and anomalies of liquid water. We identify two solid phases: a square crystal (high density phase), and a triangular crystal (low density phase) and discuss the relation between the anomalies of liquid and the polymorphism of the solid. Similarly to real water, our 2D system may have the second critical point in the metastable liquid phase beyond the freezing line. In 3D, we find several sets of parameters for which two fluid–fluid phase transition lines exist: the first line between gas and liquid and the second line between high-density liquid (HDL) and low-density liquid (LDL). In all cases, the LDL phase shows no density anomaly in 3D. We relate the absence of the density anomaly with the positive slope of the LDL–HDL phase transition line.

The puzzling behavior of water at very low temperature

Although H2O has been the topic of considerable research since the beginning of the century, most of its anomalous physical properties are still not well understood. First we discuss some of the anomalies of this ‘‘complex fluid. ’’ Then we describe a qualitative interpretation in terms of percolation concepts. Finally, we discuss recent experiments and simulations relating to the liquid–liquid phase transition hypothesis that, in addition to the known critical point in water, there may exist a ‘‘second ’’ critical point at low temperatures. In particular, we discuss very recent measurements at Tsukuba of the compression-induced melting and decompression-induced melting lines of high-pressure forms of ice. We show how knowledge of these lines enables one to obtain an approximation for the Gibbs potential G(P,T) and the equation of state V(P,T) for water, both of which are consistent with the possible continuity of liquid water and the amorphous forms of solid water.

Scaling in nature: from DNA through heartbeats to weather

The purpose of this report is to describe some recent progress in applying scaling concepts to various systems in nature. We review several systems characterized by scaling laws such as DNA sequences, heartbeat rates and weather variations. We discuss the finding that the exponent alpha quantifying the scaling in DNA in smaller for coding than for noncoding sequences. We also discuss the application of fractal scaling analysis to the dynamics of heartbeat regulation, and report the recent finding that the scaling exponent alpha is smaller during sleep periods compared to wake periods. We also discuss the recent findings that suggest a universal scaling exponent characterizing the weather fluctuations.

Evidence for an unusual dynamical-arrest scenario in short-ranged colloidal systems

Extensive molecular dynamics simulation studies of particles interacting via a short-ranged attractive square-well potential are reported. The calculated loci of constant diffusion coefficient D in the temperature-packing fraction plane show a reentrant behavior, i.e., an increase of diffusivity on cooling, confirming an important part of the high volume-fraction dynamical-arrest scenario earlier predicted by theory for particles with short-ranged potentials. The more efficient localization mechanism induced by the short-range bonding provides, on average, additional free volume as compared to the hard-sphere case and results in faster dynamics.

Average time spent by Levy flights and walks on an interval with absorbing boundaries

We consider a Lévy flyer of order alpha that starts from a point x(0) on an interval [O,L] with absorbing boundaries. We find a closed-form expression for the average number of flights the flyer takes and the total length of the flights it travels before it is absorbed. These two quantities are equivalent to the mean first passage times for Lévy flights and Lévy walks, respectively. Using fractional differential equations with a Riesz kernel, we find exact analytical expressions for both quantities in the continuous limit. We show that numerical solutions for the discrete Lévy processes converge to the continuous approximations in all cases except the case of alpha-->2, and the cases of x(0)-->0 and x(0)-->L. For alpha>2, when the second moment of the flight length distribution exists, our result is replaced by known results of classical diffusion. We show that if x(0) is placed in the vicinity of absorbing boundaries, the average total length has a minimum at alpha=1, corresponding to the Cauchy distribution. We discuss the relevance of this result to the problem of foraging, which has received recent attention in the statistical physics literature.

Dynamics of Plaque Formation in Alzheimer's Disease

Plaques that form in the brains of Alzheimer patients are made of deposits of the amyloid-beta peptide. We analyze the time evolution of amyloid-beta deposition in immunostained brain slices from transgenic mice. We find that amyloid-beta deposits appear in clusters whose characteristic size increases from 14 microm in 8-month-old mice to 22 microm in 12-month-old mice. We show that the clustering has implications for the biological growth of amyloid-beta by presenting a growth model that accounts for the experimentally observed structure of individual deposits and predicts the formation of clusters of deposits and their time evolution.