1
|
Grzybowski A, Koperwas K, Paluch M. Role of anisotropy in understanding the molecular grounds for density scaling in dynamics of glass-forming liquids. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2024; 87:084501. [PMID: 38861964 DOI: 10.1088/1361-6633/ad569d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/14/2024] [Accepted: 06/11/2024] [Indexed: 06/13/2024]
Abstract
Molecular Dynamics (MD) simulations of glass-forming liquids play a pivotal role in uncovering the molecular nature of the liquid vitrification process. In particular, much focus was given to elucidating the interplay between the character of intermolecular potential and molecular dynamics behaviour. This has been tried to achieve by simulating the spherical particles interacting via isotropic potential. However, when simulation and experimental data are analysed in the same way by using the density scaling approaches, serious inconsistency is revealed between them. Similar scaling exponent values are determined by analysing the relaxation times and pVT data obtained from computer simulations. In contrast, these values differ significantly when the same analysis is carried out in the case of experimental data. As discussed thoroughly herein, the coherence between results of simulation and experiment can be achieved if anisotropy of intermolecular interactions is introduced to MD simulations. In practice, it has been realized in two different ways: (1) by using the anisotropic potential of the Gay-Berne type or (2) by replacing the spherical particles with quasi-real polyatomic anisotropic molecules interacting through isotropic Lenard-Jones potential. In particular, the last strategy has the potential to be used to explore the relationship between molecular architecture and molecular dynamics behaviour. Finally, we hope that the results presented in this review will also encourage others to explore how 'anisotropy' affects remaining aspects related to liquid-glass transition, like heterogeneity, glass transition temperature, glass forming ability, etc.
Collapse
Affiliation(s)
- A Grzybowski
- Institute of Physics, University of Silesia in Katowice, 41-500 Chorzów, Poland
| | - K Koperwas
- Institute of Physics, University of Silesia in Katowice, 41-500 Chorzów, Poland
| | - M Paluch
- Institute of Physics, University of Silesia in Katowice, 41-500 Chorzów, Poland
| |
Collapse
|
2
|
Yu N, Huang D, Feng Y. Melting curve of two-dimensional Yukawa systems predicted by isomorph theory. Phys Rev E 2024; 109:065212. [PMID: 39020935 DOI: 10.1103/physreve.109.065212] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2024] [Accepted: 06/04/2024] [Indexed: 07/20/2024]
Abstract
The analytical expression for the conditions of the solid-fluid phase transition, i.e., the melting curve, for two-dimensional (2D) Yukawa systems is derived theoretically from the isomorph theory. To demonstrate that the isomorph theory is applicable to 2D Yukawa systems, molecular dynamical simulations are performed under various conditions. Based on the isomorph theory, the analytical isomorphic curves of 2D Yukawa systems are derived using the local effective power-law exponent of the Yukawa potential. From the obtained analytical isomorphic curves, the melting curve of 2D Yukawa systems is directly determined using only two known melting points. The determined melting curve of 2D Yukawa systems well agrees with the previous obtained melting results using completely different approaches.
Collapse
Affiliation(s)
- Nichen Yu
- Institute of Plasma Physics and Technology, Jiangsu Key Laboratory of Frontier Material Physics and Devices, School of Physical Science and Technology, Soochow University, Suzhou 215006, China
| | - Dong Huang
- Institute of Plasma Physics and Technology, Jiangsu Key Laboratory of Frontier Material Physics and Devices, School of Physical Science and Technology, Soochow University, Suzhou 215006, China
| | - Yan Feng
- Institute of Plasma Physics and Technology, Jiangsu Key Laboratory of Frontier Material Physics and Devices, School of Physical Science and Technology, Soochow University, Suzhou 215006, China
- National Laboratory of Solid State Microstructures, Nanjing University, Nanjing 210093, China
| |
Collapse
|
3
|
Yu N, Huang D, Lu S, Khrapak S, Feng Y. Universal scaling of transverse sound speed and its isomorphic property in Yukawa fluids. Phys Rev E 2024; 109:035202. [PMID: 38632806 DOI: 10.1103/physreve.109.035202] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2024] [Accepted: 02/13/2024] [Indexed: 04/19/2024]
Abstract
Molecular dynamical simulations are performed to investigate the scaling of the transverse sound speed in two-dimensional (2D) and 3D Yukawa fluids. From the calculated diagnostics of the radial distribution function, the mean-squared displacement, and the Pearson correlation coefficient, the approximate isomorphic curves for 2D and 3D liquidlike Yukawa systems are obtained. It is found that the structure and dynamics of 2D and 3D liquidlike Yukawa systems exhibit the isomorphic property under the conditions of the same relative coupling parameter Γ/Γ_{m}=const. It is demonstrated that the reduced transverse sound speed, i.e., the ratio of the transverse sound speed to the thermal speed, is an isomorph invariant, which is a quasiuniversal function of Γ/Γ_{m}. The obtained isomorph invariant of the reduced transverse sound speed can be useful to estimate the transverse sound speed, or determine the coupling strength, with applications to dusty (complex) plasma or colloidal systems.
Collapse
Affiliation(s)
- Nichen Yu
- Institute of Plasma Physics and Technology, School of Physical Science and Technology, Jiangsu Key Laboratory of Thin Films, Soochow University, Suzhou 215006, China
| | - Dong Huang
- Institute of Plasma Physics and Technology, School of Physical Science and Technology, Jiangsu Key Laboratory of Thin Films, Soochow University, Suzhou 215006, China
| | - Shaoyu Lu
- Institute of Plasma Physics and Technology, School of Physical Science and Technology, Jiangsu Key Laboratory of Thin Films, Soochow University, Suzhou 215006, China
| | - Sergey Khrapak
- Joint Institute for High Temperatures, Russian Academy of Sciences, 125412 Moscow, Russia
| | - Yan Feng
- Institute of Plasma Physics and Technology, School of Physical Science and Technology, Jiangsu Key Laboratory of Thin Films, Soochow University, Suzhou 215006, China
| |
Collapse
|
4
|
Jin J, Hwang J, Voth GA. Gaussian representation of coarse-grained interactions of liquids: Theory, parametrization, and transferability. J Chem Phys 2023; 159:184105. [PMID: 37942867 DOI: 10.1063/5.0160567] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2023] [Accepted: 10/06/2023] [Indexed: 11/10/2023] Open
Abstract
Coarse-grained (CG) interactions determined via bottom-up methodologies can faithfully reproduce the structural correlations observed in fine-grained (atomistic resolution) systems, yet they can suffer from limited extensibility due to complex many-body correlations. As part of an ongoing effort to understand and improve the applicability of bottom-up CG models, we propose an alternative approach to address both accuracy and transferability. Our main idea draws from classical perturbation theory to partition the hard sphere repulsive term from effective CG interactions. We then introduce Gaussian basis functions corresponding to the system's characteristic length by linking these Gaussian sub-interactions to the local particle densities at each coordination shell. The remaining perturbative long-range interaction can be treated as a collective solvation interaction, which we show exhibits a Gaussian form derived from integral equation theories. By applying this numerical parametrization protocol to CG liquid systems, our microscopic theory elucidates the emergence of Gaussian interactions in common phenomenological CG models. To facilitate transferability for these reduced descriptions, we further infer equations of state to determine the sub-interaction parameter as a function of the system variables. The reduced models exhibit excellent transferability across the thermodynamic state points. Furthermore, we propose a new strategy to design the cross-interactions between distinct CG sites in liquid mixtures. This involves combining each Gaussian in the proper radial domain, yielding accurate CG potentials of mean force and structural correlations for multi-component systems. Overall, our findings establish a solid foundation for constructing transferable bottom-up CG models of liquids with enhanced extensibility.
Collapse
Affiliation(s)
- Jaehyeok Jin
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, 5735 S. Ellis Ave., Chicago, Illinois 60637, USA
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - Jisung Hwang
- Department of Statistics, The University of Chicago, 5747 S. Ellis Ave., Chicago, Illinois 60637, USA
| | - Gregory A Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, 5735 S. Ellis Ave., Chicago, Illinois 60637, USA
| |
Collapse
|
5
|
Kaśkosz F, Koperwas K, Grzybowski A, Paluch M. The origin of the density scaling exponent for polyatomic molecules and the estimation of its value from the liquid structure. J Chem Phys 2023; 158:144503. [PMID: 37061492 DOI: 10.1063/5.0141975] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2023] Open
Abstract
In this article, we unravel the problem of interpreting the density scaling exponent for the polyatomic molecules representing the real van der Waals liquids. Our studies show that the density scaling exponent is a weighted average of the exponents of the repulsive terms of all interatomic interactions that occur between molecules, where the potential energy of a given interaction represents its weight. It implies that potential energy is a key quantity required to calculate the density scaling exponent value for real molecules. Finally, we use the well-known method for potential energy estimation and show that the density scaling exponent could be successfully predicted from the liquid structure for fair representatives of the real systems.
Collapse
Affiliation(s)
- F Kaśkosz
- University of Silesia in Katowice, Institute of Physics, 75 Pułku Piechoty 1, 41-500 Chorzów, Poland
| | - K Koperwas
- University of Silesia in Katowice, Institute of Physics, 75 Pułku Piechoty 1, 41-500 Chorzów, Poland
| | - A Grzybowski
- University of Silesia in Katowice, Institute of Physics, 75 Pułku Piechoty 1, 41-500 Chorzów, Poland
| | - M Paluch
- University of Silesia in Katowice, Institute of Physics, 75 Pułku Piechoty 1, 41-500 Chorzów, Poland
| |
Collapse
|
6
|
Heyes DM, Dini D, Pieprzyk S, Brańka AC. Departures from perfect isomorph behavior in Lennard-Jones fluids and solids. J Chem Phys 2023; 158:134502. [PMID: 37031156 DOI: 10.1063/5.0143651] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/05/2023] Open
Abstract
Isomorphs are lines on a fluid or solid phase diagram along which the microstructure is invariant on affine density scaling of the molecular coordinates. Only inverse power (IP) and hard sphere potential systems are perfectly isomorphic. This work provides new theoretical tools and criteria to determine the extent of deviation from perfect isomorphicity for other pair potentials using the Lennard-Jones (LJ) system as a test case. A simple prescription for predicting isomorphs in the fluid range using the freezing line as a reference is shown to be quite accurate for the LJ system. The shear viscosity and self-diffusion coefficient scale well are calculated using this method, which enables comments on the physical significance of the correlations found previously in the literature to be made. The virial–potential energy fluctuation and the concept of an effective IPL system and exponent, n′, are investigated, particularly with reference to the LJ freezing and melting lines. It is shown that the exponent, n′, converges to the value 12 at a high temperature as ∼ T−1/2, where T is the temperature. Analytic expressions are derived for the density, temperature, and radius derivatives of the radial distribution function along an isomorph that can be used in molecular simulation. The variance of the radial distribution function and radial fluctuation function are shown to be isomorph invariant.
Collapse
Affiliation(s)
- D. M. Heyes
- Department of Mechanical Engineering, Imperial College London, Exhibition Road, South Kensington, London SW7 2AZ, United Kingdom
| | - D. Dini
- Department of Mechanical Engineering, Imperial College London, Exhibition Road, South Kensington, London SW7 2AZ, United Kingdom
| | - S. Pieprzyk
- Institute of Molecular Physics, Polish Academy of Sciences, M. Smoluchowskiego 17, 60-179 Poznań, Poland
| | - A. C. Brańka
- Institute of Molecular Physics, Polish Academy of Sciences, M. Smoluchowskiego 17, 60-179 Poznań, Poland
| |
Collapse
|
7
|
Saw S, Costigliola L, Dyre JC. Configurational temperature in active matter. I. Lines of invariant physics in the phase diagram of the Ornstein-Uhlenbeck model. Phys Rev E 2023; 107:024609. [PMID: 36932558 DOI: 10.1103/physreve.107.024609] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2022] [Accepted: 01/23/2023] [Indexed: 06/18/2023]
Abstract
This paper shows that the configurational temperature of liquid-state theory, T_{conf}, defines an energy scale, which can be used for adjusting model parameters of active Ornstein-Uhlenbeck particle (AOUP) models in order to achieve approximately invariant structure and dynamics upon a density change. The required parameter changes are calculated from the variation of a single configuration's T_{conf} for a uniform scaling of all particle coordinates. The resulting equations are justified theoretically for models involving a potential-energy function with hidden scale invariance. The validity of the procedure is illustrated by computer simulations of the Kob-Andersen binary Lennard-Jones AOUP model, showing the existence of lines of approximate invariance of the reduced-unit radial distribution function and time-dependent mean-square displacement.
Collapse
Affiliation(s)
- Shibu Saw
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Lorenzo Costigliola
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Jeppe C Dyre
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
8
|
Jiang Y, Weeks ER, Bailey NP. Isomorphs in sheared binary Lennard-Jones glass: Transient response. Phys Rev E 2023; 107:014610. [PMID: 36797950 DOI: 10.1103/physreve.107.014610] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2022] [Accepted: 01/06/2023] [Indexed: 01/31/2023]
Abstract
We have studied shear deformation of binary Lennard-Jones glasses to investigate the extent to which the transient part of the stress strain curves is invariant when the thermodynamic state point is varied along an isomorph. Shear deformations were carried out on glass samples of varying stability, determined by cooling rate, and at varying strain rates, at state points deep in the glass. Density changes up to and exceeding a factor of two were made. We investigated several different methods for generating isomorphs but none of the previously developed methods could generate sufficiently precise isomorphs given the large density changes and nonequilibrium situation. Instead, the temperatures for these higher densities were chosen to give state points isomorphic to the starting state point by requiring the steady-state flow stress for isomorphic state points to be invariant in reduced units. In contrast to the steady-state flow stress, we find that the peak stress on the stress strain curve is not invariant. The peak stress decreases by a few percent for each ten percent increase in density, although the differences decrease with increasing density. Analysis of strain profiles and nonaffine motion during the transient phase suggests that the root of the changes in peak stress is a varying tendency to form shear bands, with the largest tendency occurring at the lowest densities. We suggest that this reflects the effective steepness of the potential; a higher effective steepness gives a greater tendency to form shear bands.
Collapse
Affiliation(s)
- Yonglun Jiang
- Department of Physics, Emory University, 400 Dowman Drive, Atlanta, Georgia 30322, USA
| | - Eric R Weeks
- Department of Physics, Emory University, 400 Dowman Drive, Atlanta, Georgia 30322, USA
| | - Nicholas P Bailey
- "Glass and Time," IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
9
|
Kaśkosz F, Koperwas K, Paluch M. The role of the excluded volume in the molecular dynamics for molecular systems revealed by the direct computational approach. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120321] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/31/2022]
|
10
|
Knudsen S, Todd BD, Dyre JC, Hansen JS. Generalized hydrodynamics of the Lennard-Jones liquid in view of hidden scale invariance. Phys Rev E 2021; 104:054126. [PMID: 34942805 DOI: 10.1103/physreve.104.054126] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2021] [Accepted: 10/12/2021] [Indexed: 11/07/2022]
Abstract
In recent years lines along which structure and dynamics are invariant to a good approximation, so-called isomorphs, have been identified in the thermodynamic phase diagrams of several model liquids and solids. This paper reports computer simulation data of the transverse and longitudinal collective dynamics at different length scales along an isomorph of the Lennard-Jones system. Our findings are compared to corresponding results along an isotherm and an isochore. Confirming the theoretical prediction, the reduced-unit dynamics of the transverse momentum density is invariant to a good approximation along the isomorph on all time and length scales. Likewise, the wave-vector dependent shear-stress autocorrelation function is found to be isomorph invariant (with minor deviations at very short times). A similar invariance is not seen along the isotherm or the isochore. Using a spatially nonlocal hydrodynamic model for the transverse momentum-density time-autocorrelation function, the macroscopic shear viscosity and its wave dependence are determined, demonstrating that the shear viscosity is isomorphic invariant on all length scales studied. This analysis implies the existence of a length scale that is isomorph invariant in reduced units, i.e., which characterizes each isomorph. The transverse sound-wave velocity, the Maxwell relaxation time, and the rigidity shear modulus are also isomorph invariant. In contrast to the isomorph invariance of all aspects of the transverse dynamics, the reduced-unit dynamics of the mass density is not invariant on length scales longer than the interparticle distance. By fitting to a generalized hydrodynamic model, we extract values for the wave-vector-dependent thermal diffusion coefficient, sound attenuation coefficient, and adiabatic sound velocity. The isomorph variation of these quantities in reduced units on long length scales can be eliminated by scaling with the density-scaling exponent, a fundamental quantity in the isomorph theory framework; this is an empirical observation that remains to be explained theoretically.
Collapse
Affiliation(s)
- Solvej Knudsen
- "Glass and Time", IMFUFA, Department of Science and Environment, Roskilde University, Postbox 260, DK-4000 Roskilde, Denmark.,Department of Mathematics, School of Science, Computing and Engineering Technologies, Swinburne University of Technology, Hawtorn, Victoria 3122, Australia
| | - B D Todd
- Department of Mathematics, School of Science, Computing and Engineering Technologies, Swinburne University of Technology, Hawtorn, Victoria 3122, Australia
| | - Jeppe C Dyre
- "Glass and Time", IMFUFA, Department of Science and Environment, Roskilde University, Postbox 260, DK-4000 Roskilde, Denmark
| | - J S Hansen
- "Glass and Time", IMFUFA, Department of Science and Environment, Roskilde University, Postbox 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
11
|
Carter BMGD, Royall CP, Dyre JC, Ingebrigtsen TS. Isomorphs in nanoconfined liquids. SOFT MATTER 2021; 17:8662-8677. [PMID: 34515711 PMCID: PMC8494272 DOI: 10.1039/d1sm00233c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/13/2021] [Accepted: 08/03/2021] [Indexed: 06/13/2023]
Abstract
We study in this paper the possible existence of Roskilde-simple liquids and their isomorphs in a rough-wall nanoconfinement. Isomorphs are curves in the thermodynamic phase diagram along which structure and dynamics are invariant in suitable nondimensionalized units. Two model liquids using molecular dynamics computer simulations are considered: the single-component Lennard-Jones (LJ) liquid and the Kob-Andersen binary LJ mixture, both of which in the bulk phases are known to have good isomorphs. Nanoconfinement is implemented by adopting a slit-pore geometry with fcc crystalline walls; this implies inhomogenous density profiles both parallel and perpendicular to the confining walls. Despite this fact and consistent with an earlier study [Ingebrigtsen et al., Phys. Rev. Lett., 2013, 111, 235901] we find that these two nanoconfined liquids have isomorphs to a good approximation. More specifically, we show good invariance along the isomorphs of inhomogenous density profiles, mean-square displacements, and higher-order structures probed using the topological cluster classification algorithm. Our study thus provides an alternative framework for understanding nanoconfined liquids.
Collapse
Affiliation(s)
- Benjamin M G D Carter
- H.H. Wills Physics Laboratory, Tyndall Avenue, Bristol, BS8 1TL, UK
- Bristol Centre for Functional Nanomaterials, Tyndall Avenue, Bristol, BS8 1TL, UK
| | - C Patrick Royall
- H.H. Wills Physics Laboratory, Tyndall Avenue, Bristol, BS8 1TL, UK
- School of Chemistry, University of Bristol, Cantock's Close, Bristol, BS8 1TS, UK
- Centre for Nanoscience and Quantum Information, Tyndall Avenue, Bristol, BS8 1FD, UK
| | - Jeppe C Dyre
- Department of Science and Environment, Roskilde University, Postbox 260, DK-4000 Roskilde, Denmark.
| | - Trond S Ingebrigtsen
- Department of Science and Environment, Roskilde University, Postbox 260, DK-4000 Roskilde, Denmark.
| |
Collapse
|
12
|
Yang X, Kim D, May EF, Bell IH. Entropy Scaling of Thermal Conductivity: Application to Refrigerants and Their Mixtures. Ind Eng Chem Res 2021. [DOI: 10.1021/acs.iecr.1c02154] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Xiaoxian Yang
- Fluid Science & Resources Division, Department of Chemical Engineering, University of Western Australia, Crawley, Western Australia 6009, Australia
| | - Dongchan Kim
- Fluid Science & Resources Division, Department of Chemical Engineering, University of Western Australia, Crawley, Western Australia 6009, Australia
| | - Eric F. May
- Fluid Science & Resources Division, Department of Chemical Engineering, University of Western Australia, Crawley, Western Australia 6009, Australia
| | - Ian H. Bell
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, United States
| |
Collapse
|
13
|
Attia E, Dyre JC, Pedersen UR. Extreme case of density scaling: The Weeks-Chandler-Andersen system at low temperatures. Phys Rev E 2021; 103:062140. [PMID: 34271644 DOI: 10.1103/physreve.103.062140] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2021] [Accepted: 05/28/2021] [Indexed: 11/07/2022]
Abstract
This paper studies numerically the Weeks-Chandler-Andersen system, which is shown to obey hidden scale invariance with a density-scaling exponent that varies from below 5 to above 500. This unprecedented variation makes it advantageous to use the fourth-order Runge-Kutta algorithm for tracing out isomorphs. Good isomorph invariance of structure and dynamics is observed over more than three orders of magnitude temperature variation. For all state points studied, the virial potential-energy correlation coefficient and the density-scaling exponent are controlled mainly by the temperature. Based on the assumption of statistically independent pair interactions, a mean-field theory is developed that rationalizes this finding and provides an excellent fit to data at low temperatures.
Collapse
Affiliation(s)
- Eman Attia
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, 4000 Roskilde, Denmark
| | - Jeppe C Dyre
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, 4000 Roskilde, Denmark
| | - Ulf R Pedersen
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, 4000 Roskilde, Denmark
| |
Collapse
|
14
|
Bell IH, Delage-Santacreu S, Hoang H, Galliero G. Dynamic Crossover in Fluids: From Hard Spheres to Molecules. J Phys Chem Lett 2021; 12:6411-6417. [PMID: 34232673 DOI: 10.1021/acs.jpclett.1c01594] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
We propose a simple and generic definition of a demarcation reconciling structural and dynamic frameworks when combined with the entropy scaling framework. This crossover line between gas- and liquid-like behaviors is defined as the curve for which an individual property, the contribution to viscosity due to molecules' translation, is exactly equal to a collective property, the contribution to viscosity due to molecular interactions. Such a definition is shown to be consistent with the one based on the minima of the kinematic viscosity. For the hard sphere, this is shown to be an exact solution. For Lennard-Jones spheres and dimers and for some simple real fluids, this relation holds very well. This crossover line passes nearby the critical point, and for all studied fluids, it is well captured by the critical excess entropy curve for atomic fluids, emphasizing the link between transport properties and local structure.
Collapse
Affiliation(s)
- Ian H Bell
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, United States
| | - Stéphanie Delage-Santacreu
- Université de Pau et des Pays de l'Adour, e2s UPPA, Laboratoire de Mathematiques et de leurs Applications de Pau (IPRA, CNRS UMR5142), Pau 64000, France
| | - Hai Hoang
- Institute of Fundamental and Applied Sciences, Duy Tan University, 10C Tran Nhat Duat Street, District 1, Ho Chi Minh City 700000, Vietnam
- Faculty of Natural Sciences, Duy Tan University, Da Nang 550000, Vietnam
| | - Guillaume Galliero
- Université de Pau et des Pays de l'Adour, e2s UPPA, TOTAL, CNRS, LFCR, UMR 5150, Laboratoire des fluides complexes et leurs reservoirs, Pau 64000, France
| |
Collapse
|
15
|
Nandi MK, Bhattacharyya SM. Microscopic Theory of Softness in Supercooled Liquids. PHYSICAL REVIEW LETTERS 2021; 126:208001. [PMID: 34110221 DOI: 10.1103/physrevlett.126.208001] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/04/2020] [Accepted: 04/19/2021] [Indexed: 06/12/2023]
Abstract
We introduce a new measure of the structure of a liquid which is the softness of the mean-field potential developed by us earlier. We find that this softness is sensitive to small changes in the structure. Then, we study its correlation with the supercooled liquid dynamics. The study involves a wide range of liquids (fragile, strong, attractive, repulsive, and active) and predicts some universal behaviors such as the softness being linearly proportional to the temperature and inversely proportional to the activation barrier of the dynamics with system dependent proportionality constants. We establish a master equation between the dynamics and the softness parameter and show that, indeed, the dynamics, when scaled by the temperature and system dependent parameters, show a data collapse when plotted against softness. The dynamics of fragile liquids show a strong softness dependence, whereas that of strong liquids show a much weaker softness dependence. We also connect the present study with the earlier studies of softness involving machine learning (ML), thus, providing a theoretical framework for understanding the ML results.
Collapse
Affiliation(s)
- Manoj Kumar Nandi
- Polymer Science and Engineering Division, CSIR-National Chemical Laboratory, Pune-411008, India
| | | |
Collapse
|
16
|
Rahman M, Carter BMGD, Saw S, Douglass IM, Costigliola L, Ingebrigtsen TS, Schrøder TB, Pedersen UR, Dyre JC. Isomorph Invariance of Higher-Order Structural Measures in Four Lennard-Jones Systems. Molecules 2021; 26:molecules26061746. [PMID: 33804670 PMCID: PMC8003765 DOI: 10.3390/molecules26061746] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2021] [Revised: 03/07/2021] [Accepted: 03/09/2021] [Indexed: 11/16/2022] Open
Abstract
In the condensed liquid phase, both single- and multicomponent Lennard-Jones (LJ) systems obey the "hidden-scale-invariance" symmetry to a good approximation. Defining an isomorph as a line of constant excess entropy in the thermodynamic phase diagram, the consequent approximate isomorph invariance of structure and dynamics in appropriate units is well documented. However, although all measures of the structure are predicted to be isomorph invariant, with few exceptions only the radial distribution function (RDF) has been investigated. This paper studies the variation along isomorphs of the nearest-neighbor geometry quantified by the occurrence of Voronoi structures, Frank-Kasper bonds, icosahedral local order, and bond-orientational order. Data are presented for the standard LJ system and for three binary LJ mixtures (Kob-Andersen, Wahnström, NiY2). We find that, while the nearest-neighbor geometry generally varies significantly throughout the phase diagram, good invariance is observed along the isomorphs. We conclude that higher-order structural correlations are no less isomorph invariant than is the RDF.
Collapse
Affiliation(s)
- Mahajabin Rahman
- Department of Physics, Emory University, Atlanta, GA 30322, USA;
| | | | - Shibu Saw
- “Glass and Time”, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark; (S.S.); (I.M.D.); (L.C.); (T.S.I.); (T.B.S.); (U.R.P.)
| | - Ian M. Douglass
- “Glass and Time”, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark; (S.S.); (I.M.D.); (L.C.); (T.S.I.); (T.B.S.); (U.R.P.)
| | - Lorenzo Costigliola
- “Glass and Time”, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark; (S.S.); (I.M.D.); (L.C.); (T.S.I.); (T.B.S.); (U.R.P.)
| | - Trond S. Ingebrigtsen
- “Glass and Time”, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark; (S.S.); (I.M.D.); (L.C.); (T.S.I.); (T.B.S.); (U.R.P.)
| | - Thomas B. Schrøder
- “Glass and Time”, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark; (S.S.); (I.M.D.); (L.C.); (T.S.I.); (T.B.S.); (U.R.P.)
| | - Ulf R. Pedersen
- “Glass and Time”, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark; (S.S.); (I.M.D.); (L.C.); (T.S.I.); (T.B.S.); (U.R.P.)
| | - Jeppe C. Dyre
- “Glass and Time”, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark; (S.S.); (I.M.D.); (L.C.); (T.S.I.); (T.B.S.); (U.R.P.)
- Correspondence:
| |
Collapse
|
17
|
Singh A, Bhattacharyya SM, Singh Y. Emergence of cooperatively reorganizing cluster and super-Arrhenius dynamics of fragile supercooled liquids. Phys Rev E 2021; 103:032611. [PMID: 33862818 DOI: 10.1103/physreve.103.032611] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2020] [Accepted: 03/12/2021] [Indexed: 06/12/2023]
Abstract
In this paper, we develop a theory to calculate the structural relaxation time τ_{α} of fragile supercooled liquids. Using the information of the configurational entropy and structure, we calculate the number of dynamically free, metastable, and stable neighbors around a central particle. In supercooled liquids, the cooperatively reorganizing clusters (CRCs) in which the stable neighbors form "stable" nonchemical bonds with the central particle emerge. For an event of relaxation to take place, these bonds have to reorganize irreversibly; the energy involved in the processes is the effective activation energy of relaxation. The theory brings forth a temperature T_{a} and a temperature-dependent parameter ψ(T) which characterize slowing down of dynamics on cooling. It is shown that the value of ψ(T) is equal to 1 for T>T_{a}, indicating that the underlying microscopic mechanism of relaxation is dominated by the entropy-driven processes, while for T<T_{a}, ψ(T) decreases on cooling, indicating the emergence of the energy-driven processes. This crossover of ψ(T) from high to low temperatures explains the crossover seen in τ_{α}. The dynamics of systems that may have similar static structure but very different dynamics can be understood in terms of ψ(T). We present results for the Kob-Anderson model for three densities and show that the calculated values of τ_{α} are in excellent agreement with simulation values for all densities. We also show that when ψ(T), τ_{α}, and other quantities are plotted as a function of T/T_{a} (or T_{a}/T), the data collapse on master curves.
Collapse
Affiliation(s)
- Ankit Singh
- Department of Physics, Banaras Hindu University, Varanasi 221 005, India
| | | | - Yashwant Singh
- Department of Physics, Banaras Hindu University, Varanasi 221 005, India
| |
Collapse
|
18
|
Mairhofer J. A Residual Entropy Scaling Approach for Viscosity Based on the GERG-2008 Equation of State. Ind Eng Chem Res 2021. [DOI: 10.1021/acs.iecr.0c04938] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
19
|
Lucco Castello F, Tolias P, Dyre JC. Testing the isomorph invariance of the bridge functions of Yukawa one-component plasmas. J Chem Phys 2021; 154:034501. [PMID: 33499616 DOI: 10.1063/5.0036226] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
It has been recently conjectured that bridge functions remain nearly invariant along phase diagram lines of constant excess entropy for the broad class of R-simple liquids. To test this hypothesis, the bridge functions of Yukawa systems are computed outside the correlation void with the Ornstein-Zernike inversion method employing structural input from ultra-accurate molecular dynamics simulations and inside the correlation void with the cavity distribution method employing structural input from ultra-long specially designed molecular dynamics simulations featuring a tagged particle pair. Yukawa bridge functions are revealed to be isomorph invariant to a very high degree. The observed invariance is not exact, however, since isomorphic deviations exceed the overall uncertainties.
Collapse
Affiliation(s)
- F Lucco Castello
- Space and Plasma Physics, Royal Institute of Technology, Stockholm SE-100 44, Sweden
| | - P Tolias
- Space and Plasma Physics, Royal Institute of Technology, Stockholm SE-100 44, Sweden
| | - J C Dyre
- Glass and Time, IMFUFA, Roskilde University, Roskilde DK-4000, Denmark
| |
Collapse
|
20
|
Koperwas K, Grzybowski A, Paluch M. Virial-potential-energy correlation and its relation to density scaling for quasireal model systems. Phys Rev E 2021; 102:062140. [PMID: 33466035 DOI: 10.1103/physreve.102.062140] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2020] [Accepted: 11/18/2020] [Indexed: 11/07/2022]
Abstract
In this paper, we examine the virial- and the potential-energy correlation for quasireal model systems. This correlation constitutes the framework of the theory of the isomorph in the liquid phase diagram commonly examined using simple liquids. Interestingly, our results show that for the systems characterized by structural anisotropy and flexible bonds, the instantaneous values of total virial and total potential energy are entirely uncorrelated. It is due to the presence of the intramolecular interactions because the contributions to the virial and potential energy resulting from the intermolecular interactions still exhibit strong linear dependence. Interestingly, in contrast to the results reported for simple liquids, the slope of the mentioned linear dependence is different than the values of the density scaling exponent. However, our findings show that for quasireal materials, the slope of dependence between the virial and potential energy (resulting from the intermolecular interactions) strongly depends on the interval of intermolecular distances that are taken into account. Consequently, the value of the slope of the discussed relationship, which enables satisfactory density scaling, can be obtained. Interestingly, this conclusion is supported by the results obtained for analogous systems without intermolecular attraction, for which the value the slope of the virial-potential-energy correlation is independent of considered intermolecular distances, directly corresponds to the exponent of the intermolecular repulsion, and finally leads to accurate density scaling.
Collapse
Affiliation(s)
- K Koperwas
- University of Silesia in Katowice, Institute of Physics, 75 Pułku Piechoty 1, 41-500 Chorzów, Poland and Silesian Center for Education and Interdisciplinary Research SMCEBI, 75 Pułku Piechoty 1a, 41-500 Chorzów, Poland
| | - A Grzybowski
- University of Silesia in Katowice, Institute of Physics, 75 Pułku Piechoty 1, 41-500 Chorzów, Poland and Silesian Center for Education and Interdisciplinary Research SMCEBI, 75 Pułku Piechoty 1a, 41-500 Chorzów, Poland
| | - M Paluch
- University of Silesia in Katowice, Institute of Physics, 75 Pułku Piechoty 1, 41-500 Chorzów, Poland and Silesian Center for Education and Interdisciplinary Research SMCEBI, 75 Pułku Piechoty 1a, 41-500 Chorzów, Poland
| |
Collapse
|
21
|
Saw S, Dyre JC. Structure of the Lennard-Jones liquid estimated from a single simulation. Phys Rev E 2021; 103:012110. [PMID: 33601502 DOI: 10.1103/physreve.103.012110] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2020] [Accepted: 12/08/2020] [Indexed: 11/07/2022]
Abstract
Combining the recent Piskulich-Thompson approach [Z. A. Piskulich and W. H. Thompson, J. Chem. Phys. 152, 011102 (2020)JCPSA60021-960610.1063/1.5135932] with isomorph theory, from a single simulation the structure of a single-component Lennard-Jones (LJ) system is obtained at an arbitrary state point in almost the whole liquid region of the temperature-density phase diagram. The LJ system exhibits two temperature ranges where the van't Hoff assumption that energetic and entropic forces are temperature independent is valid to a good approximation. A method to evaluate the structure at an arbitrary state point along an isochore from the knowledge of structures at two temperatures on the isochore is also discussed. We argue that, in general, the structure of any hidden scale-invariant system obeying the van't Hoff assumption in the whole range of temperatures can be determined in the whole liquid region of the phase diagram from a single simulation.
Collapse
Affiliation(s)
- Shibu Saw
- "Glass and Time," IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Jeppe C Dyre
- "Glass and Time," IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
22
|
Bell IH. Entropy Scaling of Viscosity - II: Predictive Scheme for Normal Alkanes. JOURNAL OF CHEMICAL AND ENGINEERING DATA 2020; 65:10.1021/acs.jced.0c00749. [PMID: 34121765 PMCID: PMC8191377 DOI: 10.1021/acs.jced.0c00749] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
In this work, a residual entropy value 6/10 of the way between the critical point and a value of -2/3 of Boltzmann's constant is shown to collapse the scaled viscosity for the family of normal alkanes. Based on this approach, a nearly universal correlation is proposed that can reproduce 95% of the experimental data for normal alkanes within ±18% (without removal of clearly erroneous data). This universal correlation has no new fluid-specific empirical parameters and is based on experimentally accessible values. This collapse is shown to be valid to a residual entropy half way between the critical point and the triple point, beyond which the macroscopically-scaled viscosity has a super-exponential dependence on residual entropy, terminating at the triple point. A key outcome of this study is a better understanding of entropy scaling for fluids with intramolecular degrees of freedom. A study of the transport and thermodynamic properties at the triple point rounds out the analysis.
Collapse
Affiliation(s)
- Ian H Bell
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, CO 80305
| |
Collapse
|
23
|
Abstract
This paper generalizes isomorph theory to systems that are not in thermal equilibrium. The systems are assumed to be R-simple, i.e., to have a potential energy that as a function of all particle coordinates R obeys the hidden-scale-invariance condition U(Ra) < U(Rb) ⇒ U(λRa) < U(λRb). "Systemic isomorphs" are introduced as lines of constant excess entropy in the phase diagram defined by density and systemic temperature, which is the temperature of the equilibrium state point with the average potential energy equal to U(R). The dynamics is invariant along a systemic isomorph if there is a constant ratio between the systemic and the bath temperature. In thermal equilibrium, the systemic temperature is equal to the bath temperature and the original isomorph formalism is recovered. The new approach rationalizes within a consistent framework previously published observations of isomorph invariance in simulations involving nonlinear steady-state shear flows, zero-temperature plastic flows, and glass-state isomorphs. This paper relates briefly to granular media, physical aging, and active matter. Finally, we discuss the possibility that the energy unit defining the reduced quantities should be based on the systemic rather than the bath temperature.
Collapse
Affiliation(s)
- Jeppe C Dyre
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
24
|
Bell IH, Dyre JC, Ingebrigtsen TS. Excess-entropy scaling in supercooled binary mixtures. Nat Commun 2020; 11:4300. [PMID: 32855393 PMCID: PMC7453028 DOI: 10.1038/s41467-020-17948-1] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2020] [Accepted: 07/22/2020] [Indexed: 11/09/2022] Open
Abstract
Transport coefficients, such as viscosity or diffusion coefficient, show significant dependence on density or temperature near the glass transition. Although several theories have been proposed for explaining this dynamical slowdown, the origin remains to date elusive. We apply here an excess-entropy scaling strategy using molecular dynamics computer simulations and find a quasiuniversal, almost composition-independent, relation for binary mixtures, extending eight orders of magnitude in viscosity or diffusion coefficient. Metallic alloys are also well captured by this relation. The excess-entropy scaling predicts a quasiuniversal breakdown of the Stokes-Einstein relation between viscosity and diffusion coefficient in the supercooled regime. Additionally, we find evidence that quasiuniversality extends beyond binary mixtures, and that the origin is difficult to explain using existing arguments for single-component quasiuniversality.
Collapse
Affiliation(s)
- Ian H Bell
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, CO, 80305, USA
| | - Jeppe C Dyre
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, Postbox 260, Roskilde, DK-4000, Denmark
| | - Trond S Ingebrigtsen
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, Postbox 260, Roskilde, DK-4000, Denmark.
| |
Collapse
|
25
|
Bell IH, Galliero G, Delage-Santacreu S, Costigliola L. An entropy scaling demarcation of gas- and liquid-like fluid behaviors. J Chem Phys 2020; 152:191102. [PMID: 33687260 DOI: 10.1063/1.5143854] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we propose a generic and simple definition of a line separating gas-like and liquid-like fluid behaviors from the standpoint of shear viscosity. This definition is valid even for fluids such as the hard sphere and the inverse power law that exhibit a unique fluid phase. We argue that this line is defined by the location of the minimum of the macroscopically scaled viscosity when plotted as a function of the excess entropy, which differs from the popular Widom lines. For hard sphere, Lennard-Jones, and inverse-power-law fluids, such a line is located at an excess entropy approximately equal to -2/3 times Boltzmann's constant and corresponds to points in the thermodynamic phase diagram for which the kinetic contribution to viscosity is approximately half of the total viscosity. For flexible Lennard-Jones chains, the excess entropy at the minimum is a linear function of the chain length. This definition opens a straightforward route to classify the dynamical behavior of fluids from a single thermodynamic quantity obtainable from high-accuracy thermodynamic models.
Collapse
Affiliation(s)
- Ian H Bell
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, USA
| | - Guillaume Galliero
- Universite de Pau et des Pays de l'Adour, e2s UPPA, TOTAL, CNRS, LFCR, UMR 5150, Laboratoire des fluides complexes et leurs reservoirs, Pau, France
| | - Stéphanie Delage-Santacreu
- Universite de Pau et des Pays de l'Adour, e2s UPPA, Laboratoire de Mathematiques et de leurs Applications de Pau (IPRA, CNRS UMR5142), Pau, France
| | - Lorenzo Costigliola
- Department of Science and Environment, Roskilde University, Postbox 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
26
|
Bell IH. Effective hardness of interaction from thermodynamics and viscosity in dilute gases. J Chem Phys 2020; 152:164508. [PMID: 32357769 DOI: 10.1063/5.0007583] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The hardness of the effective inverse power law (IPL) potential, which can be obtained from thermodynamics or collision integrals, can be used to demonstrate similarities between thermodynamic and transport properties. This link is investigated for systems of increasing complexity (i.e., the EXP, square-well, Lennard-Jones, and Stockmayer potentials; ab initio results for small molecules; and rigid linear chains of Lennard-Jones sites). These results show that while the two approaches do not yield precisely the same values of effective IPL exponent, their qualitative behavior is intriguingly similar, offering a new way of understanding the effective interactions between molecules, especially at high temperatures. In both approaches, the effective hardness is obtained from a double-logarithmic temperature derivative of an effective area.
Collapse
Affiliation(s)
- Ian H Bell
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, USA
| |
Collapse
|
27
|
Bacher AK, Pedersen UR, Schrøder TB, Dyre JC. The EXP pair-potential system. IV. Isotherms, isochores, and isomorphs in the two crystalline phases. J Chem Phys 2020; 152:094505. [DOI: 10.1063/1.5144871] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Andreas Kvist Bacher
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Ulf R. Pedersen
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Thomas B. Schrøder
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Jeppe C. Dyre
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
28
|
Bell IH. Entropy Scaling of Viscosity - I: A Case Study of Propane. JOURNAL OF CHEMICAL AND ENGINEERING DATA 2020; 65:10.1021/acs.jced.0c00209. [PMID: 33364635 PMCID: PMC7754705 DOI: 10.1021/acs.jced.0c00209] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
In this work, a broadly-applicable and simple approach for building high accuracy viscosity correlations is demonstrated for propane. The approach is based on the combination of a number of recent insights related to the use of residual entropy scaling, especially a new way of scaling the viscosity for consistency with the dilute-gas limit. With three adjustable parameters in the dense phase, the primary viscosity data for propane are predicted with a mean absolute relative deviation of 1.38%, and 95% of the primary data are predicted within a relative error band of less than 5%. The dimensionality of the dense-phase contribution is reduced from the conventional two dimensional approach (temperature and density) to a one-dimensional correlation with residual entropy as the independent variable. The simplicity of the model formulation ensures smooth extrapolation behavior (barring errors in the equation of state itself). The approach proposed here should be applicable to a wide range of chemical species. The supporting information includes the relevant data in tabular form and a Python implementation of the model.
Collapse
Affiliation(s)
- Ian H Bell
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, CO 80305
| |
Collapse
|
29
|
Heyes DM, Dini D, Costigliola L, Dyre JC. Transport coefficients of the Lennard-Jones fluid close to the freezing line. J Chem Phys 2019; 151:204502. [DOI: 10.1063/1.5128707] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- D. M. Heyes
- Department of Mechanical Engineering, Imperial College London, Exhibition Road, South Kensington, London SW7 2AZ, United Kingdom
| | - D. Dini
- Department of Mechanical Engineering, Imperial College London, Exhibition Road, South Kensington, London SW7 2AZ, United Kingdom
| | - L. Costigliola
- “Glass and Time,” IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - J. C. Dyre
- “Glass and Time,” IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
30
|
Jiang Y, Weeks ER, Bailey NP. Isomorph invariance of dynamics of sheared glassy systems. Phys Rev E 2019; 100:053005. [PMID: 31869994 DOI: 10.1103/physreve.100.053005] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2019] [Indexed: 06/10/2023]
Abstract
We study hidden scale invariance in the glassy phase of the Kob-Andersen binary Lennard-Jones system. After cooling below the glass transition, we generate a so-called isomorph from the fluctuations of potential energy and virial in the NVT ensemble: a set of density, temperature pairs for which structure and dynamics are identical when expressed in appropriate reduced units. To access dynamical features, we shear the system using the SLLOD algorithm coupled with Lees-Edwards boundary conditions and study the statistics of stress fluctuations and the particle displacements transverse to the shearing direction. We find good collapse of the statistical data, showing that isomorph theory works well in this regime. The analysis of stress fluctuations, in particular the distribution of stress changes over a given strain interval, allows us to identify a clear signature of avalanche behavior in the form of an exponential tail on the negative side. This feature is also isomorph invariant. The implications of isomorphs for theories of plasticity are discussed briefly.
Collapse
Affiliation(s)
- Yonglun Jiang
- Department of Physics, Emory University, 400 Dowman Drive, Atlanta, Georgia 30322, USA
| | - Eric R Weeks
- Department of Physics, Emory University, 400 Dowman Drive, Atlanta, Georgia 30322, USA
| | - Nicholas P Bailey
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, 4000 Roskilde, Denmark
| |
Collapse
|
31
|
Bell IH, Messerly R, Thol M, Costigliola L, Dyre JC. Modified Entropy Scaling of the Transport Properties of the Lennard-Jones Fluid. J Phys Chem B 2019; 123:6345-6363. [PMID: 31241958 DOI: 10.1021/acs.jpcb.9b05808] [Citation(s) in RCA: 53] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Rosenfeld proposed two different scaling approaches to model the transport properties of fluids, separated by 22 years, one valid in the dilute gas, and another in the liquid phase. In this work, we demonstrate that these two limiting cases can be connected through the use of a novel approach to scaling transport properties and a bridging function. This approach, which is empirical and not derived from theory, is used to generate reference correlations for the transport properties of the Lennard-Jones 12-6 fluid of viscosity, thermal conductivity, and self-diffusion. This approach, with a very simple functional form, allows for the reproduction of the most accurate simulation data to within nearly their statistical uncertainty. The correlations are used to confirm that for the Lennard-Jones fluid the appropriately scaled transport properties are nearly monovariate functions of the excess entropy from low-density gases into the supercooled phase and up to extreme temperatures. This study represents the most comprehensive metastudy of the transport properties of the Lennard-Jones fluid to date.
Collapse
Affiliation(s)
- Ian H Bell
- Applied Chemicals and Materials Division , National Institute of Standards and Technology , Boulder , Colorado 80305 , United States
| | - Richard Messerly
- Applied Chemicals and Materials Division , National Institute of Standards and Technology , Boulder , Colorado 80305 , United States
| | - Monika Thol
- Thermodynamics , Ruhr-Universität Bochum , Universitätsstraße 150 , 44801 Bochum , Germany
| | - Lorenzo Costigliola
- DNRF Centre "Glass and Time," IMFUFA, Department of Science and Environment , Roskilde University , Postbox 260, DK-4000 Roskilde , Denmark
| | - Jeppe C Dyre
- DNRF Centre "Glass and Time," IMFUFA, Department of Science and Environment , Roskilde University , Postbox 260, DK-4000 Roskilde , Denmark
| |
Collapse
|
32
|
Fragiadakis D, Roland C. Intermolecular distance and density scaling of dynamics in molecular liquids. J Chem Phys 2019; 150:204501. [DOI: 10.1063/1.5098455] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- D. Fragiadakis
- Naval Research Laboratory, Chemistry Division, Washington, District of Columbia 20375-5342, USA
| | - C.M. Roland
- Naval Research Laboratory, Chemistry Division, Washington, District of Columbia 20375-5342, USA
| |
Collapse
|
33
|
Nikiteas I, Heyes DM. Bounded inverse power potentials: Isomorphism and isosbestic points. J Chem Phys 2019; 150:144504. [DOI: 10.1063/1.5089491] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- I. Nikiteas
- Department of Physics, Royal Holloway, University of London, Egham, Surrey TW20 0EX, United Kingdom
| | - D. M. Heyes
- Department of Physics, Royal Holloway, University of London, Egham, Surrey TW20 0EX, United Kingdom
| |
Collapse
|
34
|
Friedeheim L, Dyre JC, Bailey NP. Hidden scale invariance at high pressures in gold and five other face-centered-cubic metal crystals. Phys Rev E 2019; 99:022142. [PMID: 30934297 DOI: 10.1103/physreve.99.022142] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2018] [Indexed: 11/07/2022]
Abstract
Recent density functional theory simulations showed that metals have a hitherto overlooked symmetry termed "hidden scale invariance" [Hummel et al., Phys. Rev. B 92, 174116 (2015)PRBMDO1098-012110.1103/PhysRevB.92.174116]. This scaling property implies the existence of lines in the thermodynamic phase diagram, so-called isomorphs, along which structure and dynamics are invariant to a good approximation when given in properly reduced units. This means that the phase diagram becomes effectively one-dimensional with regard to several physical properties. This paper investigates consequences and implications of the isomorph theory in six metallic crystals: Au, Ni, Cu, Pd, Ag, and Pt. The data are obtained from molecular dynamics simulations employing many-body effective medium theory (EMT) to model the atomic interactions realistically. We test the predictions from isomorph theory for structure and dynamics by means of the radial distribution and the velocity autocorrelation functions, as well as the prediction of instantaneous equilibration after a jump between two isomorphic state points. Many properties of crystals tend to be dominated by defects, and many of the properties associated with these defects are expected to be isomorph invariant as well. This is investigated in this paper for the case of vacancy diffusion. In regard to the perfect crystal properties, we find the predicted invariance of structure and also, though less perfectly, of dynamics. We show results on the variation of the density-scaling exponent γ, which can be related to the Grüneisen parameter, for all six metals. We consider large density changes up to a factor of two, corresponding to very high pressures. Unlike systems modeled using the Lennard-Jones potential where the density-scaling exponent γ is almost constant, this quantity varies substantially when using the EMT potential and is also strongly material dependent.
Collapse
Affiliation(s)
- Laura Friedeheim
- "Glass and Time," IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Jeppe C Dyre
- "Glass and Time," IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Nicholas P Bailey
- "Glass and Time," IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
35
|
Probing the link between residual entropy and viscosity of molecular fluids and model potentials. Proc Natl Acad Sci U S A 2019; 116:4070-4079. [PMID: 30770449 DOI: 10.1073/pnas.1815943116] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
This work investigates the link between residual entropy and viscosity based on wide-ranging, highly accurate experimental and simulation data. This link was originally postulated by Rosenfeld in 1977 [Rosenfeld Y (1977) Phys Rev A 15:2545-2549], and it is shown that this scaling results in an approximately monovariate relationship between residual entropy and reduced viscosity for a wide range of molecular fluids [argon, methane, [Formula: see text], [Formula: see text], refrigerant R-134a (1,1,1,2-tetrafluoroethane), refrigerant R-125 (pentafluoroethane), methanol, and water] and a range of model potentials (hard sphere, inverse power, Lennard-Jones, and Weeks-Chandler-Andersen). While the proposed "universal" correlation of Rosenfeld is shown to be far from universal, when used with the appropriate density scaling for molecular fluids, the viscosity of nonassociating molecular fluids can be mapped onto the model potentials. This mapping results in a length scale that is proportional to the cube root of experimentally measurable liquid volume values.
Collapse
|
36
|
Abstract
This article gives an overview of excess-entropy scaling, the 1977 discovery by Rosenfeld that entropy determines properties of liquids like viscosity, diffusion constant, and heat conductivity. We give examples from computer simulations confirming this intriguing connection between dynamics and thermodynamics, counterexamples, and experimental validations. Recent uses in application-related contexts are reviewed, and theories proposed for the origin of excess-entropy scaling are briefly summarized. It is shown that if two thermodynamic state points of a liquid have the same microscopic dynamics, they must have the same excess entropy. In this case, the potential-energy function exhibits a symmetry termed hidden scale invariance, stating that the ordering of the potential energies of configurations is maintained if these are scaled uniformly to a different density. This property leads to the isomorph theory, which provides a general framework for excess-entropy scaling and illuminates, in particular, why this does not apply rigorously and universally. It remains an open question whether all aspects of excess-entropy scaling and related regularities reflect hidden scale invariance in one form or other.
Collapse
Affiliation(s)
- Jeppe C Dyre
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
37
|
Bacher AK, Schrøder TB, Dyre JC. The EXP pair-potential system. I. Fluid phase isotherms, isochores, and quasiuniversality. J Chem Phys 2018; 149:114501. [DOI: 10.1063/1.5043546] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Andreas Kvist Bacher
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Thomas B. Schrøder
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Jeppe C. Dyre
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
38
|
Costigliola L, Pedersen UR, Heyes DM, Schrøder TB, Dyre JC. Communication: Simple liquids' high-density viscosity. J Chem Phys 2018; 148:081101. [PMID: 29495764 DOI: 10.1063/1.5022058] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This paper argues that the viscosity of simple fluids at densities above that of the triple point is a specific function of temperature relative to the freezing temperature at the density in question. The proposed viscosity expression, which is arrived at in part by reference to the isomorph theory of systems with hidden scale invariance, describes computer simulations of the Lennard-Jones system as well as argon and methane experimental data and simulation results for an effective-pair-potential model of liquid sodium.
Collapse
Affiliation(s)
- Lorenzo Costigliola
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Ulf R Pedersen
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - David M Heyes
- Department of Mechanical Engineering, Imperial College London, Exhibition Road, South Kensington, London SW7 2AZ, United Kingdom
| | - Thomas B Schrøder
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Jeppe C Dyre
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
39
|
Abstract
This paper derives and discusses the configuration-space Langevin equation describing a physically aging R-simple system and the corresponding Smoluchowski equation. Externally controlled thermodynamic variables like temperature, density, and pressure enter the description via the single parameter Ts/T, in which T is the bath temperature and Ts is the "systemic" temperature defined at any time t as the thermodynamic equilibrium temperature of the state point with density ρ(t) and potential energy U(t). In equilibrium, Ts ≅ T with fluctuations that vanish in the thermodynamic limit. In contrast to Tool's fictive temperature and other effective temperatures in glass science, the systemic temperature is defined for any configuration with a well-defined density, even if it is not close to equilibrium. Density and systemic temperature define an aging phase diagram, in which the aging system traces out a curve. Predictions are discussed for aging following various density-temperature and pressure-temperature jumps from one equilibrium state to another, as well as for a few other scenarios. The proposed theory implies that R-simple glass-forming liquids are characterized by the dynamic Prigogine-Defay ratio being equal to unity.
Collapse
Affiliation(s)
- Jeppe C Dyre
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
40
|
Pedersen UR, Schrøder TB, Dyre JC. Phase Diagram of Kob-Andersen-Type Binary Lennard-Jones Mixtures. PHYSICAL REVIEW LETTERS 2018; 120:165501. [PMID: 29756931 DOI: 10.1103/physrevlett.120.165501] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/04/2017] [Revised: 01/19/2018] [Indexed: 06/08/2023]
Abstract
The binary Kob-Andersen (KA) Lennard-Jones mixture is the standard model for computational studies of viscous liquids and the glass transition. For very long simulations, the viscous KA system crystallizes, however, by phase separating into a pure A particle phase forming a fcc crystal. We present the thermodynamic phase diagram for KA-type mixtures consisting of up to 50% small (B) particles showing, in particular, that the melting temperature of the standard KA system at liquid density 1.2 is 1.028(3) in A particle Lennard-Jones units. At large B particle concentrations, the system crystallizes into the CsCl crystal structure. The eutectic corresponding to the fcc and CsCl structures is cutoff in a narrow interval of B particle concentrations around 26% at which the bipyramidal orthorhombic PuBr_{3} structure is the thermodynamically stable phase. The melting temperature's variation with B particle concentration at two constant pressures, as well as at the constant density 1.2, is estimated from simulations at pressure 10.19 using isomorph theory. Our data demonstrate approximate identity between the melting temperature and the onset temperature below which viscous dynamics appears. Finally, the nature of the solid-liquid interface is briefly discussed.
Collapse
Affiliation(s)
- Ulf R Pedersen
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Thomas B Schrøder
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Jeppe C Dyre
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
41
|
Pieprzyk S, Brańka AC, Maćkowiak S, Heyes DM. Comprehensive representation of the Lennard-Jones equation of state based on molecular dynamics simulation data. J Chem Phys 2018; 148:114505. [DOI: 10.1063/1.5021560] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- S. Pieprzyk
- Institute of Molecular Physics, Polish Academy of Sciences, M. Smoluchowskiego 17, 60-179 Poznań, Poland
| | - A. C. Brańka
- Institute of Molecular Physics, Polish Academy of Sciences, M. Smoluchowskiego 17, 60-179 Poznań, Poland
| | - Sz. Maćkowiak
- Institute of Physics, Poznań University of Technology, Piotrowo 3, 60-965 Poznań, Poland
| | - D. M. Heyes
- Department of Physics, Royal Holloway, University of London, Egham, Surrey TW20 0EX, United Kingdom
| |
Collapse
|
42
|
Lötgering-Lin O, Fischer M, Hopp M, Gross J. Pure Substance and Mixture Viscosities Based on Entropy Scaling and an Analytic Equation of State. Ind Eng Chem Res 2018. [DOI: 10.1021/acs.iecr.7b04871] [Citation(s) in RCA: 55] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Oliver Lötgering-Lin
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| | - Matthias Fischer
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| | - Madlen Hopp
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| | - Joachim Gross
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| |
Collapse
|
43
|
Ransom TC, Oliver WF. Glass Transition Temperature and Density Scaling in Cumene at Very High Pressure. PHYSICAL REVIEW LETTERS 2017; 119:025702. [PMID: 28753339 DOI: 10.1103/physrevlett.119.025702] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/02/2016] [Indexed: 06/07/2023]
Abstract
We present a new method that allows direct measurements of the glass transition temperature T_{g} at pressures up to 4.55 GPa in the glass-forming liquid cumene (isopropylbenzene). This new method uses a diamond anvil cell and can measure T_{g} at pressures of 10 GPa or greater. Measuring T_{g} at the glass→liquid transition involves monitoring the disappearance of pressure gradients initially present in the glass, but also takes advantage of the large increase in the volume expansion coefficient α_{p} at T_{g} as the supercooled or superpressed liquid is entered. Accurate T_{g}(P) values in cumene allow us to show that density scaling holds along this isochronous line up to pressures much higher than any previous study, corresponding to a density increase of 29%. Our results for cumene over this huge compression range yield ρ^{γ}/T=C, where C is a constant and where γ=4.77±0.02 for this nonassociated glass-forming system. Finally, high-pressure cumene viscosity data from the literature taken at much lower pressures and at several different temperatures, corresponding to a large dynamic range of nearly 13 orders of magnitude, are shown to superimpose on a plot of η vs ρ^{γ}/T for the same value of γ.
Collapse
Affiliation(s)
- T C Ransom
- Naval Research Laboratory, Chemistry Division, Code 6100, Washington DC 20375-5342, USA
- Physics Department, University of Arkansas, Fayetteville, Arkansas 72701, USA
| | - W F Oliver
- Physics Department, University of Arkansas, Fayetteville, Arkansas 72701, USA
| |
Collapse
|
44
|
Friisberg IM, Costigliola L, Dyre JC. Density-scaling exponents and virial potential-energy correlation coefficients for the (2n, n) Lennard-Jones system. J CHEM SCI 2017. [DOI: 10.1007/s12039-017-1307-1] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
45
|
Density dependence of relaxation dynamics in glass formers, and the dependence of their fragility on the softness of inter-particle interactions. J CHEM SCI 2017. [DOI: 10.1007/s12039-017-1314-2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
46
|
Ingebrigtsen TS, Tanaka H. Effect of Energy Polydispersity on the Nature of Lennard-Jones Liquids. J Phys Chem B 2016; 120:7704-13. [PMID: 27434103 DOI: 10.1021/acs.jpcb.6b05486] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In the companion paper [ Ingebrigtsen , T. S. ; Tanaka , H. J. Phys. Chem. B 2015 , 119 , 11052 ] the effect of size polydispersity on the nature of Lennard-Jones (LJ) liquids, which represent most molecular liquids without hydrogen bonds, was studied. More specifically, it was shown that even highly size polydisperse LJ liquids are Roskilde-simple (RS) liquids. RS liquids are liquids with strong correlation between constant volume equilibrium fluctuations of virial and potential energy and are simpler than other types of liquids. Moreover, it was shown that size polydisperse LJ liquids have isomorphs to a good approximation. Isomorphs are curves in the phase diagram of RS liquids along which structure, dynamics, and some thermodynamic quantities are invariant in dimensionless (reduced) units. In this paper, we study the effect of energy polydispersity on the nature of LJ liquids. We show that energy polydisperse LJ liquids are RS liquids. However, a tendency of particle segregation, which increases with the degree of polydispersity, leads to a loss of strong virial-potential energy correlation but is mitigated by increasing temperature and/or density. Isomorphs are a good approximation also for energy polydisperse LJ liquids, although particle-resolved quantities display a somewhat poorer scaling compared to the mean quantities along the isomorph.
Collapse
Affiliation(s)
- Trond S Ingebrigtsen
- Institute of Industrial Science, University of Tokyo , 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Hajime Tanaka
- Institute of Industrial Science, University of Tokyo , 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| |
Collapse
|
47
|
Costigliola L, Schrøder TB, Dyre JC. Freezing and melting line invariants of the Lennard-Jones system. Phys Chem Chem Phys 2016; 18:14678-90. [PMID: 27186598 DOI: 10.1039/c5cp06363a] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The invariance of several structural and dynamical properties of the Lennard-Jones (LJ) system along the freezing and melting lines is interpreted in terms of isomorph theory. First the freezing/melting lines of the LJ system are shown to be approximated by isomorphs. Then we show that the invariants observed along the freezing and melting isomorphs are also observed on other isomorphs in the liquid and crystalline phases. The structure is probed by the radial distribution function and the structure factor and dynamics are probed by the mean-square displacement, the intermediate scattering function, and the shear viscosity. Studying these properties with reference to isomorph theory explains why the known single-phase melting criteria hold, e.g., the Hansen-Verlet and the Lindemann criteria, and why the Andrade equation for the viscosity at freezing applies, e.g., for most liquid metals. Our conclusion is that these empirical rules and invariants can all be understood from isomorph theory and that the invariants are not peculiar to the freezing and melting lines, but hold along all isomorphs.
Collapse
Affiliation(s)
- Lorenzo Costigliola
- "Glass and Time", IMFUFA, Department of Sciences, Roskilde University, Postbox 260, DK-4000 Roskilde, Denmark.
| | | | | |
Collapse
|
48
|
Perakis F, De Marco L, Shalit A, Tang F, Kann ZR, Kühne TD, Torre R, Bonn M, Nagata Y. Vibrational Spectroscopy and Dynamics of Water. Chem Rev 2016; 116:7590-607. [DOI: 10.1021/acs.chemrev.5b00640] [Citation(s) in RCA: 225] [Impact Index Per Article: 28.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Affiliation(s)
- Fivos Perakis
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Luigi De Marco
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
- Department of Chemistry and James Franck Institute, University of Chicago, 929 East 57th Street, Chicago, Illinois 60637, United States
| | - Andrey Shalit
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Fujie Tang
- International Center for Quantum Materials, Peking University, 5 Yiheyuan Road, Haidian, Beijing 100871, China
| | - Zachary R. Kann
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, United States,
| | - Thomas D. Kühne
- Department of Chemistry, University of Paderborn, Warburger Strasse 100, D-33098 Paderborn, Germany
| | - Renato Torre
- European Lab for Nonlinear Spectroscopy and Dipartimento di Fisica e Astronomia, Università di Firenze, Via Nello Carrara 1, Sesto Fiorentino (Firenze) I-50019, Italy
| | - Mischa Bonn
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Yuki Nagata
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| |
Collapse
|
49
|
Veldhorst AA, Dyre JC, Schrøder TB. Scaling of the dynamics of flexible Lennard-Jones chains: Effects of harmonic bonds. J Chem Phys 2016; 143:194503. [PMID: 26590538 DOI: 10.1063/1.4934973] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The previous paper [A. A. Veldhorst et al., J. Chem. Phys. 141, 054904 (2014)] demonstrated that the isomorph theory explains the scaling properties of a liquid of flexible chains consisting of ten Lennard-Jones particles connected by rigid bonds. We here investigate the same model with harmonic bonds. The introduction of harmonic bonds almost completely destroys the correlations in the equilibrium fluctuations of the potential energy and the virial. According to the isomorph theory, if these correlations are strong a system has isomorphs, curves in the phase diagram along which structure, dynamics, and the excess entropy are invariant. The Lennard-Jones chain liquid with harmonic bonds does have curves in the phase diagram along which the structure and dynamics are invariant. The excess entropy is not invariant on these curves, which we refer to as "pseudoisomorphs." In particular, this means that Rosenfeld's excess-entropy scaling (the dynamics being a function of excess entropy only) does not apply for the Lennard-Jones chain with harmonic bonds.
Collapse
Affiliation(s)
- Arno A Veldhorst
- Department of Sciences, DNRF Center "Glass and Time," IMFUFA, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Jeppe C Dyre
- Department of Sciences, DNRF Center "Glass and Time," IMFUFA, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Thomas B Schrøder
- Department of Sciences, DNRF Center "Glass and Time," IMFUFA, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
50
|
Abstract
The location of the melting line (ML) of the Lennard-Jones (LJ) system and its associated physical properties are investigated using molecular dynamics computer simulation. The radial distribution function and the behavior of the repulsive and attractive parts of the potential energy indicate that the ML is not a single isomorph, but the isomorphic state evolves gradually with temperature, i.e., it is only "locally isomorphic." The state point dependence of the unitless isomorphic number, X̃, for a range of static and dynamical properties of the LJ system in the solid and fluid states, and for fluid argon, are also reported. The quantity X̃ typically varies most with state point in the vicinity of the triple point and approaches a plateau in the high density (temperature) limit along the ML.
Collapse
Affiliation(s)
- D M Heyes
- Department of Physics, Royal Holloway, University of London, Egham, Surrey TW20 0EX, United Kingdom
| | - A C Brańka
- Institute of Molecular Physics, Polish Academy of Sciences, M. Smoluchowskiego 17, 60-179 Poznań, Poland
| |
Collapse
|