ORIGINAL RESEARCH article
Local Density Correlations in Liquids
- 1Shull-Wollan Center – Joint Institute for Neutron Sciences, Department of Materials Science and Engineering, University of Tennessee, Knoxville, Knoxville, TN, United States
- 2Department of Physics and Astronomy, University of Tennessee, Knoxville, Knoxville, TN, United States
- 3Oak Ridge National Laboratory, Oak Ridge, TN, United States
The structure and dynamics of liquids and glasses are usually discussed in terms of the short-range atomic structure and the nature of atomic bonds. However, many properties of liquids are rather similar and only weakly dependent on chemistry. Here we suggest that these properties are controlled not so much by the short-range atomic correlations but mainly by the medium-range density correlations which are similar among various liquids and supersede chemical details of the system. Taking metallic alloy liquids as examples we discuss the origin of the medium-range density correlations and argue they are only weakly dependent on the detailed nature of the atomic bond. We also discuss some consequences of the medium-range correlations on atomic transport. By extrapolating the coherence of density correlation to infinity we can define an ideal glass state which is very diverse in the local atomic structure but is characterized by long-range density correlations. This ideal state is very different from those previously proposed based on the coherence in the atomic structure.
The atomic structure of liquid and glass appears random, but unlike gas they are condensed matters held together by cohesive force. The density of a liquid is comparable, or even higher as in the case of water, than that of a solid of the same composition. Therefore, atoms are strongly correlated in position and momentum, and the structure and dynamics of liquid should be described by not where the atoms are (structure), but in terms of the correlations. Actually the “structure of liquid” is an oxymoron, because liquid is intrinsically dynamic. In solids, the time-scales of the structure, diffusion and phonon are all very different, so that they can be characterized by different scattering methods; diffraction (elastic scattering) for structure, inelastic scattering for phonon and quasielastic scattering for diffusion. However, in liquids the time-scales for them are all in the same order of magnitude, about a pico-second (ps). Consequently liquid shows no elastic scattering , and phonons are overdamped and often localized. Thus, not only the method of characterization needs to be reconsidered, but also the terminology needs to be readdressed; structure, diffusion and phonon are all mixed up and are deeply entangled. We need a new experimental approach and a new description to characterize such mixed dynamics. The same difficulty is found for a large family of soft-matter including colloids and biological matters. Actually to call them soft-matter is missing the point, because softness is relative, and is not the most important characteristic of their nature. They may be more adeptly called dynamic aperiodic matter (DAM), rather than soft-matter.
A powerful approach to describe the dynamics of DAM is the time-dependent pair-distribution function (PDF) called the Van Hove function (VHF) . The VHF is often spelled the van Hove function, but he is not Dutch but Belgian; we follow the original spelling. The VHF has been known for a long time, but only recently it became possible to determine it by inelastic scattering experiments as shown in Figure 1 . At the same time the behavior of the VHF was examined by molecular dynamics (MD) simulation, which gave a new insight on the interpretation of the VHF . In particular, it became clear that the dynamics of the medium-range correlations beyond the first nearest neighbors is quite distinct from that of the short-range structure represented by the first neighbor shell. Whereas, the short-range structure reflects the local chemistry, the medium-range part of the VHF describes collective density fluctuations which are more universal and less dependent on chemistry.
Figure 1. The Van Hove function of water at room temperature determined by inelastic x-ray scattering . The cut at t = 0 corresponds to the snapshot PDF. As time proceeds the nearest neighbors move away, while the second neighbors move in, showing dynamic correlations.
Van Hove Function
The structure of liquid and glass is usually described by the atomic pair-distribution function (PDF),
where ri is the position of the i-th atom at time t, ρ0 is the atomic density of N atoms, and <…> denotes thermal average. However, for liquid this is misleading because ri varies with time. The PDF is merely the same-time (snapshot) correlation function. The time-dependent, or time-delayed, correlation is described by the VHF,
where ri(t) is the position of the i-th atom at time t. At t = 0, G(r, 0), is just the g(r). The VHF describes how g(r) decays with time, and is obtained from the dynamic structure factor measured by inelastic scattering, S(Q, ω) ,
through the double Fourier-transformation,
where Q = |Q|. Sometimes S(Q, ω) is Fourier-transformed to the time domain by
which is called the intermediate scattering function and is widely used for the study of relaxation phenomena . The VHF is obtained by the Fourier-transformation of F(Q, t) to real space. We should also note that the energy integral of S(Q, ω),
yields the snapshot PDF, g(r), upon Fourier-transformation. Therefore, the regular x-ray diffraction measurement with low energy resolution produces the snapshot PDF. For a long time this has been the only way to assess the structural nature of liquid. That is why g(r) is called the “structure of liquid.” Now that the VHF is experimentally accessible, we should focus on the VHF. To perform the double Fourier-transformation by Equation (4) requires the knowledge of S(Q, ω) over a wide Q-ω space, which has been difficult for a long time. But recent progress in the source and instrumentation [6, 7] made it possible to determine the VHF in a reasonable time.
Medium-range Density Correlation
At first it was assumed that the VHF decays uniformly for all distances, controlled by diffusion . Then the decay rate of F(Q, t) should vary as 1/Q2 as expected for diffusion. However, De Gennes showed that the decay rate of F(Q, t) depends non-monotonically on Q . In particular, the quasi-elastic width, Δω, of S(Q, ω) shows strong narrowing at the first peak of S(Q), now known as the De Gennes narrowing. This effect is generally considered to be the evidence of collective excitations, but our recent results have shown otherwise. As shown in Figure 2 the decay time of the VHF increases linearly with distance, with the rate related to the Maxwell relaxation time . In general the r-dependent decay time depends on r as,
where a is the nearest neighbor distance. The power exponent χ depends on the spatial dimension, D, as . The argument for this dimensional dependence is as follows. The number of atoms in the shell from r to r + dr, Ndr, is proportional to rD−1. Then the rate of fluctuation is proportional to , thus to r(D−1)/2. The De Gennes narrowing occurs as a natural consequence of this r-dependence, and is not necessarily indicative of collective excitation. Whereas, the first peak of the VHF describes the dynamics of the nearest neighbors, the higher order peaks describe the dynamics of many atoms. The number of atoms involved in the first peak is the coordination number, NC, which is 12–14 for close-packed systems. Then the number of atoms in the n-th shell is of the order of n2NC, quite a large number as n increases. The width of the peaks of g(r) beyond the first peak is of the order of a/2, whereas the uncertainty in the atomic position, for instance the phonon amplitude, is smaller by an order of magnitude. That is why at large distances the PDF peaks do not reflect detailed atomic positions, but a coarse-grained density correlations.
Figure 2. Linear dependence of the decay time of the Van Hove function on distance, simulated for liquid Fe .
At temperatures above the viscosity crossover temperature, TA , which is about twice the glass transition temperature for metallic liquids , the phonon propagation length, λp = vtτM, where vt is the transverse sound velocity, becomes smaller than the interatomic distance . In other words the structural relaxation time, τM, becomes shorter than the time-scale of phonon, , where ωD is the Debye frequency. Then atoms separated by a distance beyond the range of the interatomic potential cannot communicate, and should show no correlation. The reason why the PDF shows oscillation beyond the first peak even at high temperatures is that these peaks represent not the detailed structural correlations of atoms, but coarse-grained medium-range density correlations. These density correlations interact via long wave longitudinal waves which survive well above TA  and modulate only the local density, not the detailed structure. For this reason the medium-range correlations represent the coarse-grained density rather than the atomic structure.
More than a century ago Ornstein and Zernike published a paper with a remarkable insight . They pointed out that g(r) has two components,
where c(r) describes short-range correlation, and m(r) medium-range correlation. The c(r) is directly influenced by the interatomic potential, whereas m(r) is a medium-range consequence of c(r), and as discussed above reflects coarse-grained density correlations. Then, they proposed a mean-field closure equation,
where h(r) = c(r) + m(r). Then, as a result of integration the medium-range part of h(r) reflects the effect of the interatomic potential only indirectly, usually weakly. We now move to Q space by the Fourier-transformation. Equation (9) becomes ,
c(r) has a relatively sharp peak at the interatomic distance, a, represented by
where f (r) is the peak shape function. Then,
where F(Q) is the Fourier-transform of f (r). Thus c(Q) has the first peak at Q = Q1 ~ 5π/2a. The medium-range part of h(r), m(r), is determined primarily by the first peak of S(Q) , so that only the first peak of c(Q) at Q1 is of interest for m(r). Note that the effects of the interatomic potential and composition are included solely in f (r), thus in F(Q), which is a slowly varying function of Q. Consequently the first peak of c(Q) is nearly universal, only weakly dependent on composition and the nature of bonding. As a result the behavior of m(r) is remarkably general. This makes sense because m(r) represents medium-range density correlations which are almost independent of atomic details. For oxide and polymeric glasses if we use the first sharp diffraction peak (FSDP) , similar results should be obtained.
We expand the first peak of c(Q) around Q1
Because c(r) has a relatively sharp peak F(Q) is a slowly varying function. Then,
Thus, the medium-range coherence depends on both the physical density and the strength of the nearest neighbor correlation, c(Q1). Now the first peak of S(Q) should be Lorentzian by Equation (16). Indeed it is Lorentzian for simple liquids . For some liquids it is more Gaussian, and S(Q1) is relatively small, as discussed elsewhere (Ryu et al., unpublished). The Fourier-transform of the Lorentzian is the screened Coulomb (or Yukawa) potential,
As shown in Figure 3 G(r) = 4πrρ0[g(r) – 1] ~ 4πrh(r) shows an exponential decay with sin(Q1r) oscillation. This form is widely found, even when the shape of the first peak of S(Q) is more Gaussian.
Figure 3. Exponential decay of |G(r)|, where G(r) = 4πrρ0[g(r) – 1], for Pd42.5Ni7.5Cu30P20 liquid at 600 K (Tg = 573K) .
The coherence length, ξs, varies with temperature, through the temperature dependence of c(Q1) by Equation (17). At temperatures above TA τM is shorter than the phonon time-scale, . Thus atoms cannot communicate through phonons, and the transport is controlled by action of cutting the bond, anankeons [12, 19]. Then, the atomic transport is controlled by the nearest neighbors alone, and medium-range order has no direct relevance. The snapshot structure represented by S(Q) has little or no bearing to the transport properties. However, at temperatures well below TA the dynamics is more cooperative, and is controlled by the density correlations represented by ξs(T). Also the time-scale of structural change is much slower than phonons. In this case the snapshot structure plays the usual role of the structure in determining the dynamic properties such as transport and mechanical deformation. The activation energy for viscosity, Ea(T),
Near the glass transition Ea(T) for Pd42.5Ni7.5Cu30P20 liquid depends on the coherence volume as,
as shown in Figure 4. The coherence volume normalized by the atomic volume indicates the number of atoms involved in the coherent volume. Equation (20) means that the activation energy for viscosity depends linearly on the number of atoms involved in the local shear deformation event . Through this relationship the temperature variation of ξs is related to liquid fragility . Similar arguments were advanced for the height of the first peak of S(Q) for metallic liquids [21, 22] and the width of the FSDP of complex liquids .
Figure 4. Dependence of the activation energy for viscosity, Ea(T), on the coherence length, ξs(T), of density correlation .
Various theories, notably the mode-coupling theory , use the first peak of S(Q) as a parameter to represent the atomic structure. Ironically their success originates from the fact that the first peak of S(Q) does not describe the detailed atomic structure but represent the coarse grained density fluctuations which are more general and do not depend on detailed chemistry. Conversely they indeed catch the universal aspects of the behavior of liquid and glass.
Ideal Glass State
The Equation (17) is quite enlightening. If
ξs is negligibly small, and the liquid structure has no correlation. When Z = 2/3, ξs = a. As Z approaches unity the coherence length ξs quickly increases, and diverges when Z = 1. Then h(Q) has a Bragg peak at Q1. In other words,
is the critical condition for the coherence length, ξs, to diverge, attaining the ideal glass state. By taking the limit of ξs → ∞ with Equation (18), its medium-range order is given by
Then, G0(r) = 4πrρ0[g0(r) – 1] has long-range oscillation. This ideal glass state has long-range density correlation without positional order. Quasicrystal  was the first example of a state with long-range order without positional order. It is a crystal in six-dimensions, characterized by six order vectors. The ideal glass is a crystal in the infinite dimensions. Therefore the magnitude of the coherence length at the glass transition is a good measure of how close the structure is to the ideal glass state (Ryu et al., unpublished).
The S(Q) of the ideal glass state has a Bragg sphere with the radius of Q1. The single atom density is given by
where Ω is the solid angle of the vector Q1, and the phase factor δ(Ω) varies randomly with Ω to avoid pile-up of atomic density at the same place. Each density wave has the same wavelength, λ1 = 2π/Q1 = (4/5)a, thus it is fairly coarse-grained. In order to specify the atomic position with higher precision higher order Bragg peaks are needed, just as in crystalline solids. For this reason the density correlation established in the ideal glass does not represent the order in the atomic structure, such as the icosahedral order. Indeed the model of ideal glass obtained by the reverse Monte-Carlo simulation to fit to G0(r) has very diverse local structures. The local topology by Voronoi polyhedral analysis  shows there is no dominant topology with a wide distribution of local structures. For instance icosahedra, which dominate the simple liquid [25, 26], has the population of only 0.7% in this model of ideal glass. A liquid is stabilized by configurational entropy. This diversity in the local structure is an excellent signature of such a high entropy state. The high entropy state is achieved either by increasing temperature (usual way), or mixing a large number of elements. The discovery of bulk metallic glasses with many elements [27, 28] is consistent with the second approach. On the other hand the structure dominated by a single local structural motif, such as an icosahedron, has a low entropy and is far from ideal. It easily forms an ordered state including the crystalline state.
The idea that the growth of icosahedral clusters which are incompatible with periodicity is the origin of glass formation has been widely discussed [29–31]. The current discussion on the ideal glass assumes it is a low enthalpy , low entropy state dominated by one or two local structures . In contrast, the ideal glass proposed here is characterized by the diversity in the local structure and high configurational entropy, which is the hallmark of the liquid state. The dominance by a particular structural motif, such as icosahedron, is seen in a single component system which is a poor glass former, or in glasses with particular compositions . In our view it is a sign of nano-crystallization rather than glass formation.
Usually the structure of liquid is described in terms of the same-time two body correlation function, the snapshot pair-distribution function (PDF). However, liquid does not have a structure in the usual sense because it is inherently dynamic. It is much preferred to characterize the liquid dynamics through the dynamic two-body correlation function, the Van Hove function (VHF). The recent studies of the VHF by experiment and simulation revealed the fundamental difference in nature between the short-range order and the medium-range order. Whereas, the short-range order, represented by the first peak of the PDF, describes the atomic environment of an atom and reflects the nature of the atomic bond and some details of the interatomic potential, the medium-range order, represented by the PDF peaks beyond the first peak, describes more coarse-grained density correlations which supersedes local chemistry. Thus, various liquids with different bonding and chemistry share very similar medium-range order, characterized by the Ornstein-Zernike form and the coherence in density correlation. This density correlation is coarse-grained and does not represent detailed atom-atom correlation. Therefore, it is nearly independent of the details of local structure. The ideal glass state obtained by extrapolation has long-range correlations in density, but its local structure is very diverse and it has no structural coherence at the atomic level. In our view many properties of liquids, such as viscosity just above the glass transition, depend on such density correlations rather than atomic structural correlations. This explains why many properties of liquids depend only weakly upon chemistry, and are largely common for liquids with varying nature of atomic bonds. The efforts to seek the origin of these properties directly to local atomic structure may be misguided.
Data Availability Statement
The datasets for this study can be obtained from the author (firstname.lastname@example.org).
This work was conceived, executed, and authored by TE.
This work was supported by the U.S. Department of Energy, Office of Science, Basic Energy Sciences, Materials Sciences and Engineering Division.
Conflict of Interest
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
6. Mason TE, Abernathy D, Anderson I, Ankner J, Egami T, Ehlers G, et al. The spallation neutron source in Oak Ridge: a powerful tool for materials research. Phys B. (2006) 385–86:955–60. doi: 10.1016/j.physb.2006.05.281
7. Baron AQR. High-resolution inelastic X-ray scattering I: Context, spectrometers, samples, and superconductors. In: Jaeschke E, Khan S, Schneider J, Hastings J, editors. Synchrotron Light Sources and Free-Electron Lasers. Cham: Springer International Publishing (2015). p. 1–68. doi: 10.1007/978-3-319-04507-8_41-1
17. Voylov DN, Griffin PJ, Mercado B, Keum JK, Nakanishi M, Novikov VN, et al. Correlation between temperature variations of static and dynamic properties in glass-forming liquids. Phys Rev E. (2016) 94:060603. doi: 10.1103/PhysRevE.94.060603
32. Swallen SF, Kearns KL, Mapes MK, Kim YS, McMahon RJ, Ediger MD, et al. Organic glasses with exceptional thermodynamic and kinetic stability. Science. (2007) 315:353–56. doi: 10.1126/science.1135795
Keywords: liquid, density correlation, Van Hove function, liquid dynamics, medium-range correlation
Citation: Egami T (2020) Local Density Correlations in Liquids. Front. Phys. 8:50. doi: 10.3389/fphy.2020.00050
Received: 16 December 2019; Accepted: 20 February 2020;
Published: 28 February 2020.
Edited by:Kia Ling Ngai, Institute for Physical Chemical Processes (CNR), Italy
Reviewed by:Evan Ma, Johns Hopkins University, United States
Hai-Bin Yu, Huazhong University of Science and Technology, China
Copyright © 2020 Egami. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Takeshi Egami, email@example.com