Arxiv:1702.07221V1 [Cond-Mat.Soft] 23 Feb 2017
Total Page:16
File Type:pdf, Size:1020Kb
View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Queen Mary Research Online Direct links between dynamical, thermodynamic and structural properties of liquids: modelling results L. Wang1, C. Yang1, M. T. Dove1, Yu. D. Fomin2, V. V. Brazhkin2, and K. Trachenko1 1 School of Physics and Astronomy, Queen Mary University of London, Mile End Road, London, E1 4NS, UK and 2 Institute for High Pressure Physics, RAS, 142190, Moscow, Russia We develop an approach to liquid thermodynamics based on collective modes. We perform exten- sive molecular dynamics simulations of noble, molecular and metallic liquids and provide the direct evidence that liquid energy and specific heat are well-described by the temperature dependence of the Frenkel (hopping) frequency. The agreement between predicted and calculated thermody- namic properties is seen in the notably wide range of temperature spanning tens of thousands of Kelvin. The range includes both subcritical liquids and supercritical fluids. We discuss the struc- tural crossover and inter-relationships between structure, dynamics and thermodynamics of liquids and supercritical fluids. INTRODUCTION because the interactions are both strong and system- dependent, the resulting energy and other thermody- It is an interesting fact that the liquid state has proven namic functions will also be strongly system-dependent, to be difficult to describe by theory throughout the his- precluding their calculation in general form and under- tory of condensed matter research [1{12]. The problem standing using basic principles, in contrast to solids and extends beyond condensed matter and exists in other ar- gases [12]. Consistent with this somewhat pessimistic eas where strong interactions are combined with dynam- view, the discussion of liquid thermodynamic properties ical disorder such as field theory. has remained scarce. Indeed, physics textbooks have very little, if anything, to say about liquid specific heat, in- In a weakly-interacting system such as a dense gas, the cluding textbooks dedicated to liquids [1{12]. potential energy is much smaller than the kinetic energy. As recently reviewed [17], emerging evidence advances These systems are amenable to perturbation treatment our understanding of the thermodynamics of the liquid giving corrections to the non-interacting case [5]. Pertur- state. The start point is the early theoretical idea of J bation approaches have been widely explored to calculate Frenkel [1] who proposed that liquids can be considered liquid thermodynamic properties but have not been able as solids at times smaller than liquid relaxation time, τ, to agree with experiments. For example, the analysis of the average time between two particle rearrangements at tractable models such as van der Waals or hard-spheres one point in space. This implies that phonons in liquids systems returns the gas-like result for the liquid constant- 3 will be similar to those in solids for frequencies above the volume specific heat cv = 2 kB [12{14]. This is in contrast Frenkel frequency !F: to experimental results showing that cv of monatomic liq- uids close to the melting point is nearly identical to the solid-like result, cv = 3kB and decreases to about 2kB 1 at high temperature [15, 16]. As expected on general ! > !F = (1) τ grounds, the perturbation approach does not work for strongly-interacting systems. The above argument predicts that liquids are capa- Strong interactions are successfully treated in solids, ble of supporting shear modes, the property hitherto at- crystals or glasses, where the harmonic model is a good tributable to solids only, but only for frequencies above starting point and gives the most of the vibrational en- !F. ergy. However, this approach requires fixed reference We note that low-frequency modes in liquids, sound points around which the energy expansion can be made. waves, are well-understood in the hydrodynamic regime arXiv:1702.07221v1 [cond-mat.soft] 23 Feb 2017 With small vibrations around mean atomic positions, !τ < 1 [18], however Eq. (1) denotes a distinct, solid- solids meet this requirement but liquids seemingly do not: like elastic regime of wave propagation where !τ > 1. In liquid ability to flow implies that the reference lattice is essence, this suggests the existence of a cutoff frequency non-existent. !F above which particles in the liquid can be described Therefore, liquids seemingly have no simplifying fea- by the same equations of motion as in, for example, solid tures such as small interactions of gases or small displace- glass. Therefore, liquid collective modes include both ments of solids [12]. In other words, liquids have no small longitudinal and transverse modes with frequency above parameter. One might adopt a general approach not re- !F in the solid-like elastic regime and one longitudinal hy- lying on approximations and seek to directly calculate drodynamic mode with frequency below !F (shear mode the liquid energy for a model system where interactions is non-propagating below frequency !F as discussed be- and structure are known. This meets another challenge: low). 2 Recall the earlier textbook assertion [12] that a gen- collective modes in liquids play the same role, in view of eral thermodynamic theory of liquids can not be devel- the earlier Frenkel proposal and recent experimental evi- oped because liquids have no small parameter. How is dence? We have started exploring this question [37] just this fundamental problem addressed here? According to before the high-frequency transverse modes were directly Frenkel's idea, liquids behave like solids with small oscil- measured and subsequently developed it in a number of lating particle displacements serving as a small param- ways [17]. This involves calculating the liquid energy eter. Large-amplitude diffusive particle jumps continue as the phonon energy where transverse modes propagate to play an important role, but do not destroy the exis- above !F in Eq. (1). tence of the small parameter. Instead, the jumps serve to The main aim of this paper is to provide direct com- modify the phonon spectrum: their frequency, !F, sets putational evidence to the phonon theory of liquid ther- the minimal frequency above which the small-parameter modynamics and its predictions. We achieve this by cal- description applies and solid-like modes propagate. culating the liquid energy and !F in extensive molecular It has taken a long time to verify this picture exper- dynamics simulations. In the next chapter, we briefly imentally. The experimental evidence supporting the discuss the main steps involved in calculating the liq- propagation of high-frequency modes in liquids currently uid energy. We then proceed to calculating the liquid includes inelastic X-ray, neutron and Brillouin scattering energy and Frenkel frequency independently from molec- experiments but most important evidence is recent and ular dynamics simulations using several methods which follows the deployment of powerful synchrotron sources agree with each other. We do this for three systems cho- of X-rays [19{36]. sen from different classes of liquids: noble, metallic and Early experiments detected the presence of high- molecular, and find good agreement between predicted frequency longitudinal acoustic propagating modes and and calculated results in the wide range of temperature mapped dispersion curves which were in striking resem- and pressure. The range includes both subcritical liq- blance to those in solids [19]. These and similar results uids and supercritical state below the Frenkel line where were generated at temperature just above the melting. transverse waves propagate. We calculate and analyze The measurements were later extended to high temper- liquid energy and cv using several different methods. Fi- atures considerably above the melting point, confirming nally, we discuss how our results offer insights into inter- the same result. It is now well established that liquids relationships between structure, dynamics and thermo- sustain propagating modes with wavelengths extending dynamics in liquids and supercritical fluids. down towards interatomic separations, comparable to the wave vectors of phonons in crystals at the Brillouin zone boundaries [20{32]. More recently, the same result has PHONON APPROACH TO LIQUID been asserted for supercritical fluids [24, 31, 32]. Impor- THERMODYNAMICS tantly, the propagating modes in liquids include acoustic transverse modes. These were first seen in highly viscous Calculating liquid energy fluids (see, e.g., Refs. [33, 34]), but were then studied in low-viscosity liquids on the basis of positive disper- We summarize the main result of calculation of the liq- sion [20{22, 25] (the presence of high-frequency trans- uid energy on the basis of propagating modes. A detailed verse modes increases sound velocity from the hydrody- discussion can be found in a recent review [17]. namic to the solid-like value). These studies included wa- According to the previous discussion, the propagating ter [35], where it was found that the onset of transverse modes in liquids include two transverse modes propagat- excitations coincides with the inverse of liquid relaxation ing in the solid-like elastic regime with frequency ! > !F. time [36], as predicted by Frenkel [1]. The energy of these modes, together with the energy of More recently, high-frequency transverse modes in liq- the longitudinal mode gives the liquid vibrational energy. uids were directly measured in the form of distinct dis- In addition to vibrations, particles in the liquids undergo persion branches and verified on the basis of computer diffusive jumps between quasi-equilibrium positions as modeling [26{30], and the striking similarity between discussed above. Adding the energy of these jumps to dispersion curves in liquids and their crystalline (poly- the phonon energy in the Debye model gives the total crystalline) counterparts was noted. We note that the energy of thermal motion in the liquid [17, 37]: contribution of high-frequency modes is particularly im- portant for liquid thermodynamics because these modes 3! make the largest contribution to the energy due to !F ET = NT 3 − (2) quadratic density of states.