PHYSICAL REVIEW E 85, 031203 (2012)

Two liquid states of matter: A dynamic line on a phase diagram

V. V. Brazhkin,1,* Yu. D. Fomin,1 A. G. Lyapin,1 V. N. Ryzhov,1 and K. Trachenko2
1
Institute for High Pressure Physics RAS, 142190 Troitsk Moscow Region, Russia
2
South East Physics Network and School of Physics, Queen Mary University of London, Mile End Road, London E1 4NS, United Kingdom
(Received 15 August 2011; published 30 March 2012)
It is generally agreed that the supercritical region of a liquid consists of one single state (supercritical fluid).
On the other hand, we show here that liquids in this region exist in two qualitatively different states: “rigid”
and “nonrigid” liquids. Rigid to nonrigid transition corresponds to the condition τ ≈ τ0 , where τ is the liquid
relaxation time and τ0 is the minimal period of transverse quasiharmonic waves. This condition defines a new
dynamic crossover line on the phase diagram and corresponds to the loss of shear stiffness of a liquid at all
available frequencies and, consequently, to the qualitative change in many important liquid properties. We
analyze this line theoretically as well as in real and model fluids and show that the transition corresponds to
the disappearance of high-frequency sound, to the disappearance of roton minima, qualitative changes in the
temperature dependencies of sound velocity, diffusion, viscous flow, and thermal conductivity, an increase in
particle thermal speed to half the speed of sound, and a reduction in the constant volume specific heat to 2kB per
particle. In contrast to the Widom line that exists near the critical point only, the new dynamic line is universal:
It separates two liquid states at arbitrarily high pressure and temperature and exists in systems where liquid-gas
transition and the critical point are absent altogether. We propose to call the new dynamic line on the phase
diagram “Frenkel line”.

DOI: 10.1103/PhysRevE.85.031203 PACS number(s): [Link], 62.10.+s, 63.50.−x, [Link]

I. INTRODUCTION We begin our discussion with the paper of Frenkel [1], who
provided a microscopic description of Maxwell phenomeno-
Our current understanding and discussion of basic states of
logical viscoelastic theory of liquid flow [2] by introducing
matter such as solid, liquid, and gas is illustrated by “tempera-
liquid relaxation time τ : τ is the average time between two
ture, pressure” (T ,P ) or “temperature, density” (T ,ρ) phase di-
consecutive atomic jumps in a liquid at one point in space.
agrams. Crossing a line on such diagrams corresponds to ther-
Each jump can approximately be viewed as a jump of an atom
modynamic phase transitions, leading to qualitative changes in
from its neighboring cage into a new equilibrium position with
the physical behavior of the system. Below, we show that, for
subsequent cage relaxation. These atomic jumps give a liquid
one basic state of matter, the liquid phase, an equally important
its ability to flow. The relaxation time τ is a fundamental
qualitative change in system behavior exists, which is related
flow property of a liquid, and it defines liquid viscosity η and
to the change in its dynamics rather than thermodynamics.
diffusion coefficient D.
Consequently, we propose that all liquids have two qualita-
The above picture implies that the motion of an atom in a
tively different states, and that a new additional crossover line
liquid consists of two types: quasiharmonic vibrational motion
(narrow zone) should be added to the phase diagram which
around an equilibrium position as in a solid and diffusive
separates the two states. This line does not correspond to any
motion between two neighboring positions where typical
thermodynamic phase transition, i.e., it is a crossover, and
diffusion distances exceed vibrational distances by about a
there is no discontinuity of any physical properties along the
factor of 5–10 [ Figs. 2(a)–2(c)]. Therefore, atomic motion
line. We called it the “Frenkel line” as discussed below.
in a liquid combines both elements of the short-amplitude
More specifically, a typical (T ,P ) diagram [Fig. 1(a)]
vibrational motion as in a solid and the large-amplitude
implies that a liquid is separated from a gas by the boiling
ballistic-collisional motion as in a gas. One should mention
line ending at the critical point. The diagram further implies
that Frenkel’s ideas were discussed and used for the past 20
that only one single state (frequently called “supercritical
years by Wallace [3] and Chisolm and Wallace [4] to calculate
fluid”) exists for all pressures and temperatures above the
the thermodynamic and dynamic properties of a liquid.
critical point. On the other hand, we propose that an important
In this paper, we propose that the point at which the solid-
qualitative change in a fluid behavior takes place on crossing
like oscillatory motion ceases, leaving only the gaslike motion,
our new line. Importantly, this dynamical crossover line
marks the change in the most important properties of a fluid.
extends for arbitrary values of pressure and temperature above
the critical point [Fig. 1(a)]. In addition, the new line is not
related to the critical point from the physical perspective and,
therefore, exists in systems where the liquid-gas transition is II. THE FRENKEL LINE ON THE PHASE DIAGRAM
absent altogether as is the case in some soft-matter subjects The value of τ decreases with temperature increase,
with short-range attractive forces as well as in the model spanning many orders of magnitude. On the other hand, the
soft-sphere system [Fig. 1(b)]. minimal (Debye) vibration period, τ0 (τ0 ≈ 0.1–1 ps), is weakly
temperature dependent and is mostly defined by interactions
in a given system. At certain high temperatures, the solidlike
*
brazhkin@[Link] vibration character ceases [Figs. 2(a), 2(d), and 2(e)]. This

1539-3755/2012/85(3)/031203(12) 031203-1 ©2012 American Physical Society

BRAZHKIN, FOMIN, LYAPIN, RYZHOV, AND TRACHENKO PHYSICAL REVIEW E 85, 031203 (2012)

(a) 40
Nonrigid liquid ρ =1, T=1.325 ρ =2, T=70
ρ =0.5, T=4.5 ρ =1, T=100 (a)
Dynamic line 30
Temperature

Critical
Rigid liquid 20

Displacement
point
Tm
Gas
10

Solid 4
Triple
point
2
Pressure
(b) 0

Nonrigid liquid -2
Dynamic line
Temperature

0 2 4 6 8 10 12
Rigid liquid Time
Tm -1.4
ρ =1, T=1.325 3.5 ρ =2, T=70
-1.6
Solid
-1.8
3.0
Pressure Displacement -2.0

-2.2
(b) (c)
FIG. 1. (Color online) Pressure-temperature phase diagrams (in
arbitrary units) of (a) an ordinary substance and (b) a system without 2 3 4 5 11.2 11.6
the boiling line and liquid-gas critical point. In both cases, there exists 20
a dynamic line separating rigid and nonrigid liquids. ρ =0.5, T=4.5 ρ =1, T=100
8 18

point is reached when τ becomes comparable to τ0 ,

16
τ ≈ τ0 . (1) 6
(d) 14 (e)
The sign of the approximate equality means that the values
2 3 4 5 6 1.0 1.5
can differ by tens of percent. In the following discussion,
we consider τ as the average time it takes an atom to move Time
the average interparticle distance a. Then, τ quantifies the
motion envisaged by Frenkel where an atom jumps distance FIG. 2. (Color online) Examples of particle trajectories (x coor-
a during time τ between two equilibrium positions at low dinate) for the LJ liquid at different conditions, presented in LJ units
(ρc = 0.314, Tc = 1.31). Panels (b)–(e) show selected fragments from
temperatures as well as the motion at high temperatures where
(a). Cases (b) and (c) correspond to the rigid state where vibrations
two equilibrium positions are absent altogether and the motion
are present (thin lines); (d) and (e) correspond to collisional motion in
between collisions is ballistic as in a gas.
the nonrigid state (thick lines). Time and displacement are in standard
We note that, for τ > τ0 where dynamics is mostly LJ units.
vibrational, the atomic jumps take place by activation over
the barrier created by the potential energy of interaction
[1]. Therefore, the transition from solidlike vibrations to as witnessed by the change in its elastic, structural, dynamic,
continuous gaslike ballistic motion takes place when kinetic diffusive, and thermodynamic properties.
energy K of a particle becomes comparable to the potential
energy of their interaction. Hence, condition (1) implies that
approximately, III. SOLIDLIKE ELASTIC AND STRUCTURAL
PROPERTIES OF A LIQUID
3kB T /2 ≈ Epot . (2)
An important insight from the introduction of relaxation
For most substances, the ratio of kinetic to potential energy time is that, if observation time is smaller than τ , the local
K/Epot at the melting temperature is significantly lower than structure of a liquid does not change and is the same as that
1. In this case, the particles are primarily within the range of of a solid. This enabled Frenkel to predict that a liquid should
the action of the quasiharmonic potential and, after melting, maintain solidlike shear waves at all frequencies ω > 2π /τ
relatively long-lived regions with a well-defined vibration [1]. This prediction was later confirmed experimentally for
spectrum and short-range order remain. We will elaborate on different kinds of liquids [5–8].
this below. The maximum oscillation frequency available in the system
Condition (1) is achieved on a certain new crossover line is ω0 = 2π/τ0 . Therefore, solidlike shear waves exist in the
on the phase diagram. Below, we show that crossing this line range of 2π /τ < ω < 2π/τ0. Consequently, condition (1)
leads to important qualitative changes in the system behavior (τ ≈ τ0 ) corresponds to the complete loss of shear waves and,

031203-2
TWO LIQUID STATES OF MATTER: A DYNAMIC LINE . . . PHYSICAL REVIEW E 85, 031203 (2012)

therefore, to the loss of shear resistance, or rigidity, at all is defined from the dispersion relation ω = Vs k. Using linear
frequencies existing in the system. Debye approximation and taking maximum frequency ω as
We note that the condition of the crossover (1) and Debye frequency ω0 = 2π /τ0 and kmax = π /a gives Vs =
other conditions considered below correspond to approximate 2a/τ0 . Let us now approach the fluid from high temperatures
equality. Moreover, there are wide distributions of the τ and τ0 where the oscillatory motion is lost and recall that τ is the
values, and we consider their average values. Nevertheless, in time between two consecutive collisions over distance a. Then,
all liquids, there are definite T ,P conditions corresponding to Vth ≈ a/τ , where Vth is the particle thermal velocity. Therefore,
the loss of transverselike vibrations in the excitation spectra, condition (1) implies
hence, we can speak about the definite line or, at least, about
Vs ≈ 2Vth . (3)
the narrow crossover instead of a wide crossover.
Usually, the sharp transition in the macroscopic system In condensed phases, such as solids and liquids, the speed
is associated with some divergent value. The condition (1) of sound is primarily determined by the interactions between
defines the crossover in the dynamics of a system. To propose atoms: The sound velocities are given by elastic moduli. In
a more or less strict criterion of the transition, one can dense liquids, moduli vary insignificantly with temperature
consider the analogy with percolation. Let us consider the along isochors [12], and sound velocities are weakly temper-
high-temperature state of a fluid. If, for time τ , the projection ature dependent. On the other hand, the thermal velocity of
of the particle momentum on its total displacement changes a classical particle increases with temperature without bound.
sign at least once, then this particle is considered an oscillating Therefore, a temperature range must exist where the speed of
one. A small fraction of these particles will be present sound and thermal velocity become comparable in magnitude.
even in the low-density gas state. The proportion of such The physical meaning of condition (3) is that particles cease
particles increases with the temperature decrease, and at a to feel elastic resistance of the medium and start moving in a
certain temperature, the infinite percolation cluster consisting ballistic way.
of these particles appears. We can assume that this temperature The condition for thermal velocity of the particles becoming
corresponds to the appearance of the first high-frequency half the speed of sound as a condition for the change in atomic
transverse-like excitations in the fluid. Of course, depending dynamics is quite natural. The factor of 2 arises from the fact
on the shape of the distributions of the τ and τ0 values, that the neighboring particles move in opposite phases at the
this percolation occurs at τ ≈ τ0 where the proportionality highest frequency, and their relative thermal velocity equals
coefficient can differ from 1 by a dozen percent. the speed of sound.
The ability of liquids to flow is often associated with The speed of sound in the nonrigid liquid is defined by
zero rigidity, or shear resistance, that markedly distinguishes thermal velocity, such as in a gas, and this speed should
liquids from solids. However, this implies zero rigidity at increase with the temperature on isobars. On the other hand, the
low frequencies only, whereas, at higher frequencies, a liquid speed of sound in the rigid liquid decreases with temperature
supports shear stress. On the other hand, condition (1) marks due to the decrease in the bulk modulus on isobars. Thus, on
the qualitative change, from the physical point of view, in isobars near the line defined by the conditions (1) and (3), the
system elastic properties: shear resistance is lost completely at qualitative change in the temperature dependence of the sound
all frequencies available in the liquid. Therefore, condition (1) velocity should take place.
marks the crossover between a “rigid” liquid where rigidity We note that Eq. (3) is based on the same physical
exists in a certain frequency range and a “nonrigid” liquid grounds as Eq. (2). Indeed, the speed of sound is governed
which cannot sustain rigidity at any frequency. by the elastic moduli, which are, in turn, proportional to the
The crossover (1) has important implications for the ability potential energy of the system per unit volume. However,
of liquids to undergo phase transitions, an emerging and the proportionality coefficient in Eq. (2) may be significantly
fast-developing area of research [9–11]. When τ > τ0 at (several times) different from 1. Indeed, the elastic moduli
low temperatures, local structure or definite short-range order are the second derivatives of the potential energy, i.e., the
structure (SROS) of the liquid remains unchanged during times nature of the interparticle interaction is very important for the
shorter than τ . In this case, pressure and temperature changes proportionality coefficient. For example, in Ref. [13], it was
can induce a phase transition in a liquid, accompanied by shown that a van der Waals fluid near the critical point holds
the change in the SROS and dynamics. On the other hand, kB T ≈ Epot , while a Coulomb plasma with a compensating
when τ < τ0 at high temperatures, only the random packing background near the critical point holds 3kB T ≈ Epot .
structure type can exist because atoms are in the continuous Another interesting consequence of Eq. (1) is related
gaslike state of motion. In this case, pressure can only induce to the phenomenon of “high-frequency sound” or positive
a change in density but not a phase transition with a change in dispersion of sound velocity, which is the increase in the
SROS. Therefore, the Frenkel line demarcates the regions on speed of sound at high frequencies. Frenkel predicted [1]
the phase diagram where phase transitions in a liquid operate. that this effect should exist for frequencies ω > 2π /τ . If,
as he argued, shear waves kick in at frequency 2π /τ , the
speed of sound increases from (B/ρ)1/2 to [(B + 4G/3)/ρ]1/2 ,
IV. DYNAMICS AND SOUND PROPAGATION
where B is the bulk modulus and G is the shear modulus of
The qualitative change in atomic dynamics, defined by a liquid because G becomes nonzero at this frequency. For
Eq. (1), has important consequences for sound propagation and viscous liquids, such behavior has been indeed observed many
diffusion. Let us approach the liquid from low temperatures times in the past 50 years (see, for example, Ref. [8]). The
where τ > τ0 . The speed of sound in a liquid or a solid Vs viscoelastic model was later developed in detail, including

031203-3
BRAZHKIN, FOMIN, LYAPIN, RYZHOV, AND TRACHENKO PHYSICAL REVIEW E 85, 031203 (2012)

to a collisional one [23,24], hence, the near equality of D ∗ and

Dc is not surprising in our picture. Thus, for a rather crude
1 estimation of the crossover (1) at moderate pressures, we can
use the condition,
D ≈ Dc . (5)
2
In this sense, the proposed dynamic line starts not far from
ω

3 the end of the liquid-gas transition. However, it is important

to stress that all basic conditions above, Eqs. (1)–(4), are not
1 related to the liquid-gas transition and to the existence of a
2 critical point from the physical point of view and continue to
operate in systems where the critical point is absent altogether.
3
We will further comment on this below.
0
Importantly, condition (1) corresponds to the crossover
Q between two different qualitative temperature dependencies
of diffusion D and viscosity η. Indeed, when τ > τ0 at
FIG. 3. (Color online) Schematic dispersion curves for longitude
low temperatures, τ ∼ exp(U /T ), where U is the effective
excitations in a liquid at different temperatures and fixed pressure.
activation barrier. Then, D ∼ a 2 /τ ∼ exp(−U /T ). On the
Curves 1–3 correspond to temperatures below, in the vicinity, and
other hand, when τ < τ0 at high temperatures, τ quantifies
above the dynamical line, respectively (T1 < T2 < T3 ). The inset
shows the initial parts of the same curves where the black straight
thermal motion as discussed above: τ ∼ 1/Vth ∼ 1/T 1/2 , giving
dashed lines indicate initial slopes of the dispersion curves. D ∼ T 1/2 for a low-density gas or D ∼ T α , where α is almost
constant for a dense fluid [25]. Therefore, condition (1) gives
memory function formalism, nonlocal mode coupling theory, the crossover of D from exponential to power-law temperature
etc. [14–18]. Following the prediction, the high-frequency dependence.
sound was recently observed in numerous simple nonviscous Similar to D, temperature behavior of viscosity η changes
liquids, receiving particular attention since the development of at the crossover (1). Indeed, when τ > τ0 , η decreases with
inelastic x-ray techniques [19–22]. We now observe that the temperature almost exponentially, which can be seen most
proposed crossover (1) marks the point at which the positive easily by applying the Maxwell relation η = G∞ τ , where G∞
dispersion disappears completely because, as discussed above, is the instantaneous shear modulus having weak temperature
this crossover corresponds to the complete loss of shear dependence in comparison with the exponential decrease in
waves that can exist in a liquid and corresponding loss of τ . On the other hand, when τ < τ0 , η ∼ T 1/2 . This follows
the high-frequency sound (see Fig. 3). from applying either the Stokes-Einstein-Debye relationship
Another important feature of the liquid excitation spectra is η ∼ T /D or the Maxwell relationship η = G∞ τ , where τ ∼
the presence of a rotonlike minimum of the dispersion curve 1/T 1/2 from above and recalling that G∞ is proportional to
in a second pseudo-Brillouin zone [19] (see Fig. 3). At the the kinetic ∼T term in this regime [26]. We, therefore, con-
crossover, defined by the condition (1), the contributions of clude that condition (1) corresponds to the qualitative change
jumping “ballistic” atoms and “vibrating” atoms become com- in the temperature behavior of viscosity as it crosses over
parable to each other, and the rotonlike minimum disappears from the exponential decrease at τ > τ0 to the power-law
(see Fig. 3 for the illustration of the general trend). increase at τ < τ0 .
Thermal conductivity κ is another important physical
property that should have a minimum close to the crossover
V. DIFFUSION, VISCOSITY, AND THERMAL (1). Indeed, usually, thermal conductivity behaves similar to
CONDUCTIVITY viscosity. The value κ can be estimated as κ ∼ cV V L, where
The change in the character of atomic diffusion in the liquid cV is the constant-volume specific heat, V ≈ Vs for the rigid
at the crossover (1) occurs at a particular value of diffusion liquid, V ≈ Vth for the nonrigid liquid, and L is the free path
constant D ∗ . D can be estimated as D = a 2 /6τ . When τ ≈ τ0 length of the phonon or atom, correspondingly. When τ > τ0 , κ
at the crossover (1), we have decreases with temperature because L and cV decrease. On the
other hand, when τ < τ0 , κ ∼ Vth ∼ T 1/2 since L and cV values
D ≈ D ∗ ≈ a 2 /6τ0 . (4) have a weak temperature dependence in the nonrigid liquid.
We, therefore, conclude that three important physical
Taking a ≈ 1 Å and τ0 ≈ 0.1 ps gives D ∗ ≈ 10−8 m2 /s. The values—viscosity, thermal conductivity, and sound velocity
condition D = D ∗ provides a good estimation of the dynamic should have minima on the isobars in the vicinity of the
line defined by condition (1) because it is fairly insensitive to crossover line.
the increases in pressure and temperature. Indeed, both a 2 and
τ0 decrease with pressure only slightly, and their ratio becomes
VI. THERMAL ENERGY AND SPECIFIC HEAT
even less sensitive to pressure and temperature. We note that
the D ∗ value is consistent with the experimental values of Experimentally, constant-volume specific heat of liquids at
diffusion near critical point Dc [12]. Indeed, not far from the ambient pressure decreases from about 3kB per particle around
critical point, liquids are known to lose their elastic properties, the melting temperature to about 2kB at high temperatures
and the relaxation process changes its nature from an activation [12,27]. A further decrease in specific heat with a temperature

031203-4
TWO LIQUID STATES OF MATTER: A DYNAMIC LINE . . . PHYSICAL REVIEW E 85, 031203 (2012)

increase up to gaslike values 3kB /2 is observed at high 100

pressures in the supercritical region [12]. This behavior T Critical Liquid
Liquid-gas isochore
was quantitatively explained on the basis of decreasing the ρ =ρ (a)
contribution of transverselike modes to liquid energy with the τ =τ II I
10 D=D
temperature [28]. In this model, the liquid thermal energy per

T/Tc
c =2k
atom is V =2V
  3  Solid
E τ0 1
= kB T 3 − . (6)
N τ
0.1 1 10
P/Pc 100 1000 10 000
According to Eq. (6), when τ considerably exceeds τ0 at Calculation
low temperatures, liquid energy is close to 3kB T per atom, T Critical
isochore
giving the Dulong-Petit value of specific heat of 3kB . When ρ =ρ (b) II
τ =τ
τ approaches τ0 at high temperatures, liquid energy becomes

T/Tc
10 Experiment
2kB T per atom, giving the specific heat of 2kB , consistent with Ne, D=D
the experimental results. At this temperature, shear waves are Ar, D=D I
Ar, c =2k
completely lost at all frequencies, and longitudinal modes only
Ar, V ~2V
contribute to the heat capacity. Therefore, the crossover from N , V ~2V
a rigid to a nonrigid liquid at τ ≈ τ0 is accompanied by the Liquid
decrease in the specific heat from its solid-state value to the
Solid
1
value of 2kB , 1 10 P/Pc 100

cV ≈ 2kB . (7) The crossover lines II

10
according to five criteria (c) Liquid
Equation (7) gives the minimal value of specific heat that for the simulated LJ liquid
and the experimental
T/Tc

a rigid liquid can possibly have. A further decrease in heat lines for Ar, Ne, and N I
capacity corresponds to the loss of longitudinal modes as the
temperature is increased until the dense gas state (nonrigid 1
liquid) is reached with cV = 3kB /2. Solid
Having discussed the main physical properties that change 1 10 P/Pc 100 1000
at the crossover (1), we now provide numerical and experi- T
Liquid
mental evidence supporting our proposal. 100 Liquid-gas
II I
ρ =ρ (d)
τ =τ
D=D
VII. EVIDENCE FROM MOLECULAR DYNAMICS
T/Tc

10 c =2k
SIMULATIONS AND EXPERIMENTAL DATA V =2V
N , V ~2V
We have drawn the lines determined by conditions (1),
1
(2), (3), (5), and (7), for real substances (Ar, Ne, and N2 ) Solid
and model particle systems with Lennard-Jones (LJ) and soft- 1 ρ /ρc 10
sphere potentials (SSps) (Figs. 4 and 5). Details of molecular
dynamics (MD) simulations are given in the Appendix.
FIG. 4. (Color online) (a)–(c) (T ,P ) phase diagram and (d) (T ,ρ)
It is well known that the LJ potential adequately represents
phase diagram of the LJ liquid in the relative critical coordinates.
the behavior of many molecular and rare gas liquids, whereas,
Panels (a) and (d) present calculated lines defined by different
SSp systems describe a behavior of many real substances at
criteria (see the text). Panel (b) presents some experimental data from
very high pressures where an attractive term of the interparticle Ref. [12]. Stars in panel (b) correspond to known experimental points
potential can be neglected. We have calculated points on the where the liquid loses shear waves and positive dispersion (open
phase diagram that correspond to strict equalities τ = τ0 , symbol for Ar [20] and solid symbol for N2 [22]). Experimental
3kT /2 = Epot , Vs = 2Vth , D = Dc , and cV = 2kB , hence, critical isochors are also shown in panel (b) (dashed olive line for Ne
moderate differences between the lines are not unexpected as and dashed-dotted olive line for Ar). Experimental data for criterion
discussed above. (3) for nitrogen and points (solid circles) from the simulation of the
According to Fig. 4, the position of the crossover line LJ system in Fig. 2 (with the same color) are shown in panel (d). Panel
[condition (1)] for the Lennard-Jones system agrees with the (c) illustrates a narrow crossover zone according to different criteria
line determined by condition (5) at moderate pressures P < for the LJ system and real liquids. In all cases, number I corresponds
10 − 102 Pc and with the lines determined by conditions to the rigid liquid, and number II corresponds to the nonrigid one.
(3) and (7) at high pressures P > 10Pc . At low pressures, In panel (d),the density from the liquid side is shown for the melting
the lines determined by conditions (3) and (7) shift from curve.
the crossover line (1) due to critical point anomalies and
loss of Debye approximation at low densities. As discussed (n = 6), 3kT /2 ≈ 0.3Epot ]. For a soft-sphere system, the lines
above, condition (2) is not directly based on condition (1), determined by conditions (1), (3), and (7) match well over the
and the proportionality coefficient in Eq. (2) differs from 1 entire pressure range (Fig. 5) as no critical point and associated
significantly [for LJ particles, 3kT /2 ≈ 5Epot, for SSp particles anomalies exist for this system.

031203-5
BRAZHKIN, FOMIN, LYAPIN, RYZHOV, AND TRACHENKO PHYSICAL REVIEW E 85, 031203 (2012)

100 1400
(a) II Neon

Sound velocity (m/s)

Temperature

1200

1000
10 Liquid T

I τ =τ 800
c =2k
V V 600
1 Solid K=E P=10Pc
η=η 400
P=30Pc
10 100 1000 200 P=100Pc
Pressure
1.0x10
-4 P=10Pc f~exp(Ea /T)
100 α
(b) II P=30Pc f~T

Viscosity (Pa s)
Temperature

-5
8.0x10 P=100Pc
Liquid
10 T 6.0x10
-5

τ =τ
I c =2k -5
4.0x10
1 V =2V
K=E -5
Solid 2.0x10
η=η

1 2 0.0
Density

Thermoconductivity (W m-1 K-1)

0.20
100 P=10Pc
P=30Pc
(c) II
Temperature

0.15 P=100Pc
10
Liquid
T 0.10
1 I τ =τ
c =2k
Solid K=E 0.05
0.1
10 100 0 100 200 300 400 500 600 700

Pressure Temperature (K)

100
FIG. 6. (Color online) Experimental [12] isobaric temperature
(d) II
Temperature

dependences of sound velocity, viscosity, and thermal conductivity

10 of neon. Asymptotes for viscosity at low and high temperatures
Liquid correspond to exponents of inverse temperature and power law,
1 T respectively. Arrows indicated the corresponding temperatures of
I τ =τ condition (1), τ = τ0 . The exponent α = 0.59 at P = 10Pc , 0.52
c =2k at P = 30Pc , and 0.32 at P = 100Pc .
0.1 Solid K=E

1 2 3 4
that are significantly above the melting line, implying that
Density the new line separates two distinct liquid states. We also
FIG. 5. (Color online) (T ,P ) and (T ,ρ) phase diagrams of the note that this Frenkel line is located significantly below the
simulated soft-sphere systems with (a) and (b) n = 12 and (c) and isochore, in contrast to the Widom line discussed below in
(d) n = 6. This figure presents calculated lines defined by different more detail. Importantly, the lines of all conditions [with the
criteria including those for the minimum of viscosity η along the exception of condition (5), which is approximate and is valid
isochor in panel (b). In panels (b) and (d), the density from the liquid only at moderate pressures] lie almost parallel to each other
side is shown for the melting curve. Pressure, density, and temperature in the whole range of pressures and temperatures. Notably,
are in standard soft-sphere units. pressure in this range varies by almost 4 orders of magnitude.
We further observe that, by only a slight variation in the
In addition to model systems, we find good agreement proportionality coefficients in criteria (1), (3), and (5), namely,
between the theoretical predictions and experimental data for by a factor of 1.3–1.4, the lines defined by these conditions
liquid Ar, Ne, and N2 (Fig. 4). In particular, we observe a good practically coincide, which testifies to the commonality of
match between the region of the disappearance of the positive physical mechanisms underlying the above criteria. The line
dispersion of sound velocity in liquid Ar [20] and N2 [22] and of criterion (7) coincides with the main line [criterion (1)] for
the dynamic line (Fig. 4). cV = 1.9kB , which differs from cV = 2kB by 5% only.
Despite the approximate way in which conditions (2), (3), The qualitative change in the temperature dependence of
(5), and (7) correspond to condition (1), all the lines defined sound velocity Vsl , viscosity η, and thermal conductivity κ
by different conditions are located within quite a narrow for real substances (Ne, N2 ) also occurs near the crossover
range. We note that these lines are located at temperatures line (Figs. 6 and 7). In fact, condition (1) holds close to the

031203-6
TWO LIQUID STATES OF MATTER: A DYNAMIC LINE . . . PHYSICAL REVIEW E 85, 031203 (2012)

2.8
LJ liquid, isochoreρ =1
2.6 Ar, isobar P =50P

2.4

cV /kB
2.2

2.0

1.8

1.6
0 1 2 3 4 5 10 15 20 25 30
T/Tc

FIG. 8. (Color online) The dependencies of heat capacity versus

temperature for the LJ liquid along the isochor (for ρ = 1) (orange cir-
cles with interpolation lines) and for Ar along isobar P = 50Pc (blue
solid line). The arrows indicate the points of intersection of the cor-
responding curves by the Frenkel line according to the condition (1).
-1
-1

VIII. THE FRENKEL LINE AND THE MELTING LINE

Interestingly, all crossover lines correspond to a density
increase with increasing temperature: The relation ρ ∼ T k
is met, where k ≈ 0.25 for liquid N2 , Lennard-Jones, and
soft-sphere systems with n = 12 (Figs. 4 and 5). We recall
that, for any system of particles with a uniform potential, there
are scaling relations for physical values [29]. In particular,
ρ n/3 /T = const along the melting line for the soft-sphere
1/4
system, giving ρ ∼ Tm for the melting temperature Tm for
n = 12 [29]. The similarity in the dependence of ρ(T ) for the
FIG. 7. (Color online) Experimental [12] isobaric temperature melting and dynamic lines implies similar scaling relations
dependences of sound velocity, viscosity, and thermal conductivity for the dynamic line. This point will be discussed in detail
of nitrogen. Asymptotes for viscosity at low and high temperatures elsewhere. As a result, the region of the rigid liquid does not,
correspond to exponents of inverse temperature and power law, under any pressures, disappear (Figs. 4 and 5). Consequently,
respectively. Arrows indicated the corresponding temperatures of the Frenkel line continues for arbitrarily high pressures and
condition (1), τ = τ0 . The exponent α = 0.59 at P = 10Pc , 0.53 temperatures.
at P = 30Pc , and 0.31 at P = 100Pc . It should be mentioned that the “softer” the repulsion
potential, the wider the region of existence of a rigid liquid
temperature where the deviation from the exponential viscosity (Fig. 5). We have calculated the ratio between the kinetic
decrease starts. It means that activated behavior operates even energy of the particles and the potential energy of the particle
if a particle exhibits one oscillation before jumping to another interaction K/Epot along the melting curve for the soft-sphere
position. system with different repulsion coefficients n (see Fig. 9). This
The temperature dependencies of the constant-volume heat ratio is almost constant along the melting curve and varies from
capacity for Ar (along isobars) and for the LJ liquid (along 0 for n = 3 to infinity for n → ∞ (hard spheres). The value
isochors) are presented in Fig. 8. The dynamic crossover of K/Epot near the melting curve is considerably larger than
in both cases is located at condition (7) when cV = 1 for only very large coefficients n > 30. For even larger
1.9kB . values of n, the dynamic line can hide behind the melting
One can mention that, in most of real metallic, covalent, curve. Particularly, the condition (7), cV = 2kB , is fulfilled
and ionic liquids, the Frenkel line lies at extreme experi- near the melting curve with n ≈ 60 (see Fig. 10). Thus, for
mental conditions, e.g., ∼10 GPa and ∼104 K, and only n > 50–60, the dynamic line falls within the region of crystal
shock-wave experiments can be used for its study. However, stability and separates a low-temperature almost harmonic
for molecular and rare-gas liquids, the line is situated at state of the crystal and a high-temperature strongly anharmonic
the “static pressure” experimental conditions. For example, state. In the high-temperature state, a particle moves, most of
condition (1) for Ne at P ∼ 3 GPa should take place at the time, in a ballistic way outside the region of the action
T ∼ 1100 K (approximately five to six times higher than of the potential, similar to the behavior of the particles of a
the melting temperature). For many soft-matter systems, such hard-sphere crystal for which the potential energy is equal to
as colloidal systems, macromolecules, etc., the dynamic line zero and cV = 3kB /2. In this case, the condition (1), τ ≈ τ0 ,
lies in the accessible range of pressures and temperatures as practically implies that a particle in the crystal spends similar
well. time moving (almost) harmonically and ballistically.

031203-7
BRAZHKIN, FOMIN, LYAPIN, RYZHOV, AND TRACHENKO PHYSICAL REVIEW E 85, 031203 (2012)

are not independent quantities, the structure corresponds to the

MD calculation random packing of spheres, and all transverselike excitations
4 Linear fit are damped at wavelengths on the order of one wavelength.
Therefore, it is appropriate to treat the hard-sphere liquid as a
nonrigid liquid at all temperatures.
3 The above results are important for understanding the
K/Epot

difference between the liquid-glass transition and a jamming

0.40
transition, the question that has been widely discussed and
2 0.35 debated. For very large values, n > 100, the glass transition
in an undercooled metastable melt occurs in a nonrigid
K/Epot
0.30
n=12 liquid state and, like the hard-sphere glass transition, is
1 governed by geometrical confinement conditions, and there
0.25
T=Tm
is no temperature region with activated behavior.
0.20
10 100 1000 It is known that the softer the repulsion potential, the lower
Pressure
0 the glass transition temperature or Kauzmann temperature
0 20 40 60 80 100 120 [30,31]. Thus, the region of existence of the rigid liquid for the
n soft potential is wider both in the stable state over the melting
temperature as well as in the metastable undercooled liquid
FIG. 9. (Color online) Kinetic-to-potential energy ratio for the
state.
SSp liquid at the melting curve (it is pressure independent) as a
function of the soft-sphere potential exponent n. The inset shows
It is known that, for particle systems with a narrow region of
the calculated pressure dependence of the kinetic-to-potential energy the action of the potential, the liquid-gas equilibrium line in the
ratio at the melting curve for the SSp liquid with n = 12 (the equilibrium phase diagram is absent. Instead, there is a solid
uncertainty of the calculation is on the order of the symbol size). state crystal-crystal isostructural transition ending at the criti-
Temperature and density are in standard soft-sphere units. cal point [32–34]. In this case, the discussed Frenkel line, sep-
arating the harmonic and anharmonic states of the crystal, is, in
Formally, we note that particle motion for the hard-sphere fact, a continuation of the above isostructural transition into the
system in the undercooled region can also be tentatively split supercritical region. Thus, for some colloidal and macromolec-
into two types: a small amplitude motion inside the cages ular systems for which the boiling line is absent, the Frenkel
and rare jumps over large distances. Despite the similarity of line can also lie in the region of stability of a solid phase.
this picture to the behavior of particle trajectories in a normal
liquid, there are considerable distinctions. The hard-sphere IX. THE FRENKEL LINE AGAINST THE WIDOM LINE
motion inside the cages is of a purely collisional absolutely
nonharmonic character, and the jumps only involve geometry It is interesting to note the recent attempts [20,21] to link
restrictions rather than overcoming the activation barrier. In the change in the excitation spectrum to the “thermodynamic”
this case of cV = 3/2kB , the speed of sound and thermal speed continuation of the boiling curve, the so-called Widom line, the
line of the maxima of thermodynamic properties in the vicinity
of the critical point [35]. In Ref. [20], the Widom line was
10
n=50 cv=2.0 discussed for heat capacity only and, notably, experimental
Temperature

data were extrapolated at very high pressures to give an

extrapolated thermodynamic line.
From the physical point of view, this extrapolated line is
T=Tm qualitatively different from the Frenkel line proposed here.
Indeed, the proposed crossover line is not related to the
1 extrapolation of the boiling curve and exists in systems where
liquid-gas transition and the correspondent Widom line are
6 absent altogether, including, in the model soft-sphere system,
5
n=60 T=Tm some colloidal systems, macromolecules, and so on [32–34].
4
Temperature

3 In addition, there are several other important differences. First,

cv=2.0 the maxima of heat capacity cP become smeared at T /Tc >
2 2.5 and P /Pc > 15 [12], and therefore, the extrapolation
into the high pressure range (P /Pc ≈ 100) performed in
1 Ref. [20] is not physically meaningful. Second, we have
recently calculated compressibility βT , expansion coefficient
0.94 0.96 0.98 1 1.02 αP , heat capacity cP , and density fluctuations ζ along the
Density isotherms for the Lennard-Jones particle system [36]. The
results, together with experimental data for Ar and Ne, are
FIG. 10. (Color online) Relation between melting line (blue solid summed up in Fig. 11. We observe that the thermodynamic
line) and line corresponding to the condition cV = 2.0 (orange continuation of the boiling curve gives a single line within
symbols correspond to calculation, and orange dashed lines are 10% departure in temperature from the critical point only.
approximations) for the SSp liquids with n = 50 and 60. Upon further departure, it represents a rapidly widening

031203-8
TWO LIQUID STATES OF MATTER: A DYNAMIC LINE . . . PHYSICAL REVIEW E 85, 031203 (2012)

6 critical
microscopic basis of what is now known as the “viscoelastic”
Calculation
max of α
isochore picture of liquids [1].
4 max of ζ
II (a)
We note that the past few decades have seen considerable re-
max of c
I search into the undercooled metastable state of liquids related
2 to the problem of glass transition. Several dynamic crossovers
have been found that correspond to the mode-coupling
T/Tc

1
Liquid Experiment temperature, the breakdown of the Stokes-Einstein relation,
Ar, max of ζ
Ar, max of c
the crossover from exponential to nonexponential relaxation,
Ne, max of ζ between the Arrhenius and the Vogel-Fulcher-Tammann laws,
Solid Ne, max of c and so on (see, e.g., Ref. [37] and references therein).
1 10 100 1000 On the other hand, there have been no attempts to introduce
P/Pc a dynamic crossover at high temperatures that corresponds to
a fundamental change in the character of atomic motion in a
4
liquid. We note that the line partially related to our condition
(b) (2) was considered earlier, at least, at low pressures [38,39].
3 II I This line corresponds to percolation of “physical clusters.”
T/Tc

A cluster, by definition, is made of particles linked pairwise

2
where a pair of particles is linked if the sum of their relative
Liquid Solid kinetic and interaction energies is less than zero. This line
may be not far from the Frenkel line only in the vicinity of
1 max of α the critical point. Rosenfeld also discussed a possible change
max of ζ in the nature of the dynamics of a simple liquid [40]. The
max of c
0 corresponding crossover in Ref. [40] is reflected in a change
0 1 2 3 4 from power law to exponential scaling of reduced transport
ρ /ρc coefficients when plotted versus excess entropy. In particular,
Rosenfeld predicted a minimum of viscosity as a function
FIG. 11. (Color online) (T ,P ) and (T ,ρ) phase diagrams of the of temperature for all inverse-power potentials. The crossover
simulated LJ liquid (the basic notations are the same as in Fig. 4 and discussed in Ref. [40] has been obtained only at the low-density
the Frenkel line τ ≈ τ0 and calculated detectable points of maximums region where a weakly coupled liquid is considered to be close
for thermal expansion α, fluctuations ζ , and isobaric thermal capacity to a diluted gas with taking into account only first and second
cP . The similar experimental data for Ar and Ne [12] are shown in virial coefficients. Nevertheless, for SSp systems, at least at a
panel (a). The maximum of β becomes smeared at T > 1.1Tc , and its small density region, the line of the crossover considered by
line is not shown. Rosenfeld lies close to the Frenkel line. We also note the earlier
paper of Hiwatari et al. [41], who found qualitatively different
bunch of lines instead of one single line, implying further behaviors of the velocity autocorrelation function C(t) depend-
that the extrapolation of the thermodynamic Widom line ing on the temperature and density for the SSp fluids [41].
is not unique. Third, we observe that the lines of maxima For these systems, Hiwatari et al. assumed the existence of
of all calculated properties rapidly decrease in magnitude the critical reduced density that separated the liquidlike and
and become smeared at T > (2–2.5)Tc and P > (10–15)Pc gaslike dynamical regimes. For SSp with n = 12, Hiwatari
and, therefore, cannot be meaningfully extrapolated to higher et al. estimated the critical reduced density ρ ∗ ≈ 0.54 close
pressures and temperatures. This is in contrast to the proposed to 0.58 found in our simulations for the same system. Unlike
Frenkel line, which, as discussed above, exists for arbitrarily in the paper of Hiwatari et al., the main parameter controlling
high pressures and temperatures. Finally, we note that, apart the crossover in our theory was not density but relaxation time
from the line of the maxima of the heat capacity lying close to as given by our Eq. (1). In different systems, Eq. (1) gives
the critical isochore, the lines of thermodynamic anomalies different values of reduced density at the crossover.
correspond to the decrease in density with a temperature As mentioned above, the concept of the Frenkel line may be
increase, in strong contrast to the Frenkel line (Fig. 11). applied to the colloidal liquids and macromolecular systems.
In such complex systems, different time scales are present.
There is some evidence of a mapping between long-time
molecular and Brownian dynamics [42]. It is obvious that the
X. OTHER CROSSOVER LINES AND CONCLUSIONS
dynamic rigid to non-rigid liquid transition should exist for the
To summarize, we introduced the new line separating Brownian dynamics as well.
two liquid states with qualitatively different properties. We In conclusion, we state that the proposed Frenkel line
proposed to call the line, defined by Eq. (1), the Frenkel line separates a rigid liquid where solidlike shear waves exist
to honor the contribution of Frenkel to the area of liquid and the diffusion regime is jumplike and is activated as in
dynamics. The contribution started from the microscopic a solid from a nonrigid liquid where no shear modes exist and
definition and discussion of liquid relaxation time τ . On the diffusion is collisional as in a gas. This line can be mapped
basis of this property, Frenkel made a number of important in future experiments using several conditions for liquid
predictions regarding flow, relaxation, as well as elastic and properties that we discussed, including the disappearance
phonon properties of liquids that subsequently formed the of SROS peculiarities, the disappearance of the positive

031203-9
BRAZHKIN, FOMIN, LYAPIN, RYZHOV, AND TRACHENKO PHYSICAL REVIEW E 85, 031203 (2012)

dispersion of sound velocity, as well as specific values of from the diffusion value both at the low-temperature and at the
diffusion constants, speed of sound, and specific heat: D ≈ high-temperature limit where there is a ballistic-collisional
Dc , Vs /Vth ≈ 2, and cV (T ) ≈ 2kB . regime. It makes necessary to compute the values of τ0
and τ along the same isochors at different temperatures and
ACKNOWLEDGMENTS extrapolate the data to see the cross points. The accuracy of
the estimation of the temperature of crossover τ = τ0 is about
The authors wish to thank S. M. Stishov, G. G. Simeoni, ±20%–30%; the errors in the calculations of other values are
F. Gorelli, and H. E. Stanley for valuable discussions. This less than 10%.
work has been supported by the RFBR (Grants No. 11-02- The values of τ0 and τ as well as the temperature of
00303, No. 11-02-00341, and No. 10-02-01407) and by the disappearance of solidlike vibration motion under heating also
Programs of the Presidium of RAS. K.T. is grateful to EPSRC. were extracted from the analysis of the self-intermediate scat-
tering function Fs (q,t) and velocity autocorrelation function
APPENDIX
We have studied the LJ liquid and two soft-sphere (SSp) liq-
uids with n = 12 and n = 6 in a very wide range of parameters, q=0.79 ρ =1.0
from temperature T = 0.6 (well below the critical point) to T =
100.0 and densities from 0.1 to 2.7 for the LJ liquid, T = 0
0.1–100 and density in the range of 0.1–3.33 for the SSp liquid 1.0

2
with n = 12, and T = 1.0–10 and density in the range of 0.5–4.5

dFs (q,t)/dt
for the SSp liquid with n = 6. We also have studied SSp liquids 0.8
T=10
with high values of n up to 100 for the calculation of Epot and T=20
0.6
cV . An essential property of soft spheres is that the phase di- -5 T=30
2

Fs(q,t)
agram corresponds to the equationγ = ρσ 3 ( kBεT )3/n = const. 0.4
T=50
T=70
The parameters γ for n = 12 and n = 6 systems were taken as T=100
γ12 = 1.15 and γ6 =1.22 [30]. System size in the simulations 0.2

varied depending on the density reaching 4000 particles at the 0.0 (a)
0.01 0.1 1 10
highest density. The cutoff radius was set at 2.5σ for LJ and -10 t
SSp with n = 12 and half the box size for SSp with n = 6. 0.1 1
The equations of state were integrated by the velocity Verlet t
algorithm. The temperature was kept constant during the equi- 1.0
libration by velocities rescaling. When the equilibrium was ρ =1.0
reached, the system was simulated in NVE ensemble. The usual
equilibration period was 1.5 × 106 steps, and the production 0.06
run was 0.5 × 106 steps where the time step was dt = 0.001 LJ
units. The soft-sphere system was simulated in NVE ensemble.
〈v(0)v(t)〉

The system consisted of 1000 particles, and the time step 0.5
was 0.0005. The equilibration and production periods were 0.00
3.5 × 106 and 0.5 × 106 steps, respectively. The simulations T=2
and computation of properties were performed in the same way T=5
T=10
as for the LJ system. We also use the data of the potential energy T=20
0.1 0.2 0.3
and equation of state for the LJ fluid from Ref. [43] and the T=30
melting line of the LJ system from Ref. [44]. Everywhere in the 0.0 T=50
T=70
text and in the figures, densities, temperatures, and pressures T=100 (b)
of the LJ liquid are given in the units ρ/ρc , T /Tc , and P /Pc .
The following critical parameters, averaged from literature 0.01 0.1
sources, were used for the LJ system: ρc = 0.314, Tc = 1.31. t
To calculate the Debye period τ0 , we analyzed many
particle trajectories. At low temperatures, τ0 was obtained FIG. 12. (Color online) Self-intermediate scattering functions (in
as the average time of oscillatory motion at one equilibrium the inset) and (a) their second-order time derivatives and (b) velocity
position of an atom before jumping to another position. At high autocorrelation functions of the LJ liquid along the isochor ρ = 1.
Temperature T is pointed out in LJ units. One can see the maxima
temperatures, the oscillations disappear, and the motion of the
of the second-order derivatives in panel (a). The maxima at shorter
particle consists of collisional movement only. The value of
times correspond to the presence of vibrations, whereas, the maxima
τ was calculated as the time it takes a particle to move the
at larger times correspond to ballisticlike jumps. The disappearance
distance close to the average interparticle separation. Then, of the maxima at the short time region occurs at T between 30 and
−2/3
it can be estimated as τ = ρ6D , where D is the diffusion 50, and it corresponds to the rigid liquid nonrigid liquid transition.
coefficient. D was calculated from the long-time limit of mean Similarly, there are qualitative changes in the velocity autocorrelation
square displacement using the Einstein relation. We note that function behavior in the same temperature region. Oscillations in
τ0 can be obtained in the low-temperature limit where the the rigid liquid state take place, whereas, monotonic decay of
oscillations are well pronounced, whereas, τ can be calculated autocorrelation functions in the nonrigid liquid state is observed.

031203-10
TWO LIQUID STATES OF MATTER: A DYNAMIC LINE . . . PHYSICAL REVIEW E 85, 031203 (2012)

C(t) [15], Infinite-frequency shear modulus G∞ was calculated as in

 N  Ref. [26] where it was shown that G∞ for the pair-potential
1  systems can be obtained from the radial distribution function.
Fs (q,t) = exp{i q · [rj (t) − rj (0)]} , ∂ρ −1
N j =1
The bulk modulus B = β −1 = ρ( ∂P )T was obtained from
the equation of state. The longitudinal and transverse sound
C(t) = v(0)v(t) ,
velocities were calculated as Vsl = (B/ρ)1/2 and Vst =
at different temperatures (see Fig. 12). Usually, the function (G∞ /ρ)1/2 . The heat capacities at constant volume were
Fs (q,t) is used to analyze the dynamics of glass-forming obtained by differentiating the internal energy at isochors
undercooled liquids [45,46]. However, more detailed analysis cV = ( ∂U ) .
∂T V
demonstrates that oscillating and ballistic-collisional regimes The temperature corresponding to cV = 2.0 was deter-
also can be seen from the behavior of these functions [47]. mined from the dependence of the isochoric heat capacity
From Fig. 12, one can see qualitative changes in the self- on the temperature along the isochors by linear interpolation
intermediate scattering function and the velocity autocorrela- of the data. The same method was applied to construct other
tion function at the same temperature interval around T ≈ 30 lines presented in our paper. Shear viscosity was computed by
for ρ = 1. integrating the shear stress autocorrelation function.
The kinetic energy is obtained as K/N = 3kB T /2, and The diffusion coefficients have been calculated from the
the potential energy is obtained as the total energy of the viscosity data using the Stokes-Einstein relationship. To prove
interaction. In the case of the LJ liquid, there is a complication condition (3) for real liquids, we have taken the relationship Vsl
associated with the attractive part and the sign change in = 2.3Vth [where the coefficient 2.3 for the longitudinal sound
the potential energy at compression. That is why, in this velocity instead of 2 for the shear sound velocity was taken to
case, we estimated the potential energy as Epot = ELJ (V ) − match with condition (5) for the same systems at P ∼ 100Pc ].
ELJ (V0 )+P0 (V0 − V ) where V0 corresponds to the volume of Thermodynamic and dynamic data for real fluids (Ar, Ne, and
minimal possible (negative) pressure P0 at zero temperature. N2 ) were taken from Ref. [12].

[1] J. Frenkel, Kinetic Theory of Liquids (Oxford University Press, [18] W-C. Pilgrim and C. Morcel, J. Phys.: Condens. Matter 18, R585
New York, 1946). (2006).
[2] J. C. Maxwell, Philos. Trans. R. Soc. London 157, 49 [19] T. Scopigno, G. Ruocco, and F. Sette, Rev. Mod. Phys. 77, 881
(1867). (2005).
[3] D. C. Wallace, Phys. Rev. E 56, 4179 (1997). [20] G. G. Simeoni et al., Nat. Phys. 6, 503 (2010).
[4] E. D. Chisolm and D. C. Wallace, J. Phys.: Condens. Matter 13, [21] F. Gorelli, M. Santoro, T. Scopigno, M. Krisch, and G. Ruocco,
R739 (2001). Phys. Rev. Lett. 97, 245702 (2006).
[5] M. Grimsditch, R. Bhadra, and L. M. Torell, Phys. Rev. Lett. 62, [22] F. Bencivenga et al., Europhys. Lett. 75, 70 (2006).
2616 (1989). [23] F. Bencivenga, A. Cunsolo, M. Krisch, G. Monaco, L. Orsingher,
[6] T. Pezeril, C. Klieber, S. Andrieu, and K. A. Nelson, Phys. Rev. G. Ruocco, F. Sette, and A. Vispa, Phys. Rev. Lett. 98, 085501
Lett. 102, 107402 (2009). (2007).
[7] S. Hosokawa et al., Phys. Rev. Lett. 102, 105502 (2009). [24] F. Bencivenga et al., J. Chem. Phys. 130, 064501 (2009).
[8] Y. Jeong, S. R. Nagel, and S. Bhattacharya, Phys. Rev. A 34, [25] C. R. Reid and T. K. Sherwood, in The Properties of Gases and
602 (1986). Liquids, 2nd ed. (McGraw-Hill, New York, 1966).
[9] C. J. Roberts, A. Z. Panagiotopoulos, and P. G. Debenedetti, [26] R. Zwanzig and R. D. Mountain, J. Chem. Phys. 43, 4464 (1965).
Phys. Rev. Lett. 77, 4386 (1996). [27] G. Grimvall, Phys. Scr. 11, 81 (1975).
[10] New Kinds of Phase Transitions: Transformations in Disordered [28] K. Trachenko, Phys. Rev. B 78, 104201 (2008).
Substances, Proceedings of NATO Advanced Research Work- [29] W. G. Hoover, S. G. Gray, and K. W. Johnson, J. Chem. Phys.
shop, Volga River, edited by V. V. Brazhkin, S. V. Buldyrev, 55, 1128 (1971).
V. N. Ryzhov, and H. [Link] (Kluwer, Dordrecht, [30] D. M. Heyes and A. C. Branka, Mol. Phys. 107, 309 (2009).
2002). [31] Z. Shi, P. G. Debenedetti, F. H. Stillinger, and P. Ginart, J. Chem.
[11] P. F. McMillan, J. Mater. Chem. 14, 1506 (2004). Phys. 135, 084513 (2011).
[12] NIST Chemistry WebBook [[Link] [32] G. Malescio, J. Phys.: Condens. Matter 19, 073101 (2007).
[13] S. M. Stishov, JETP Lett. 57, 196 (1993). [33] M. H. J. Hagen et al., Nature (London) 365, 425 (1993).
[14] J. P. Boon and S. Yip, Molecular Hydrodynamics (McGraw-Hill [34] A. P. Gast and W. B. Russel, Phys. Today 51(12), 24 (1998).
International Book Company, New York, 1980). [35] L. Xu et al., Proc. Natl. Acad. Sci. USA 102, 16558 (2005).
[15] U. Balucani and M. Zoppi, Dynamics of the Liquid State [36] V. V. Brazhkin et al., J. Phys. Chem. B 115, 14112 (2011).
(Clarendon, Oxford, 1994). [37] F. Mallamace et al., Proc. Natl. Acad. Sci. USA 107, 22457
[16] E. Pontecorvo, M. Krisch, A. Cunsolo, G. Monaco, A. Mermet, (2010).
R. Verbeni, F. Sette, and G. Ruocco, Phys. Rev. E 71, 011501 [38] A. Coniglio, U. De Angelis, and A. Forlani, J. Phys. A 10, 1123
(2005). (1977).
[17] T. Bryk et al., J. Chem. Phys. 133, 024502 (2010). [39] X. Campi, H. Krivine, and N. Sator, Physica A 296, 24 (2001).

031203-11
BRAZHKIN, FOMIN, LYAPIN, RYZHOV, AND TRACHENKO PHYSICAL REVIEW E 85, 031203 (2012)

[40] Y. Rosenfeld, J. Phys.: Condens. Matter 11, 5415 (1999). [44] R. Agrawal and D. A. Kofke, Mol. Phys. 85, 43 (1995).
[41] Y. Hiwatari et al., Progr. Theor. Phys. 52, 1105 (1974). [45] W. Kob and H. C. Andersen, Phys. Rev. Lett. 73, 1376 (1994).
[42] M. J. Pond, J. R. Errigton, and T. M. Truskett, Soft Matter 7, [46] T. B. Schroder, S. Sastry, J. C. Dyre, and S. C. Glotzer, J. Chem.
9859 (2011). Phys. 112, 9834 (2000).
[43] K. Meier, Ph.D. thesis, University of the Federal Armed Forces [47] G. De Lorenzi-Venneri, E. D. Chisolm, and D. C. Wallace, Phys.
Hamburg, 2002. Rev. E 78, 041205 (2008).

031203-12