1
|
Wang C, Tian W, Zhou K. Ab Initio Simulation of Liquid Water without Artificial High Temperature. J Chem Theory Comput 2024. [PMID: 39219067 DOI: 10.1021/acs.jctc.4c00650] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
Comprehending the structure and dynamics of water is crucial in various fields, such as water desalination, ion separation, electrocatalysis, and biochemical processes. While reported works show that the ab initio molecular dynamics (AIMD) can accurately portray water's structure, the artificial high temperature (AHT) from 120 to 30 K is needed to mimic the quantum nature of hydrogen-bond network from GGA, metaGGA to hybrid functionals. The AHT proves to be an inadequate approach for systems involving aqueous multiphase mixtures, such as water-solid interfaces and aqueous solutions. This is due to the activation of additional phonons in other phases, which can lead to an overestimation of the dynamics of nearby water molecules. In this work, we find that the regularized SCAN (rSCAN) functional effectively captures both the structure and dynamics of liquid water at ambient conditions without AHT. Moreover, rSCAN closely matches experimental results for the hydration structures of alkali, alkali earth, and halide ions. We anticipate that the versatile and accurate rSCAN functional will emerge as a key tool based on ab initio simulation for investigating chemical processes in aqueous environments.
Collapse
Affiliation(s)
- Chenyu Wang
- College of Energy, Soochow Institute for Energy and Materials InnovationS (SIEMIS), Jiangsu Provincial Key Laboratory for Advanced Carbon Materials and Wearable Energy Technologies, Soochow University, Suzhou 215006, China
| | - Wei Tian
- College of Energy, Soochow Institute for Energy and Materials InnovationS (SIEMIS), Jiangsu Provincial Key Laboratory for Advanced Carbon Materials and Wearable Energy Technologies, Soochow University, Suzhou 215006, China
| | - Ke Zhou
- College of Energy, Soochow Institute for Energy and Materials InnovationS (SIEMIS), Jiangsu Provincial Key Laboratory for Advanced Carbon Materials and Wearable Energy Technologies, Soochow University, Suzhou 215006, China
| |
Collapse
|
2
|
Heyes DM, Dini D, Pieprzyk S, Brańka AC, Costigliola L. Models to predict configurational adiabats of Lennard-Jones fluids and their transport coefficients. J Chem Phys 2024; 161:084502. [PMID: 39193943 DOI: 10.1063/5.0225650] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2024] [Accepted: 08/08/2024] [Indexed: 08/29/2024] Open
Abstract
A comparison is made between three simple approximate formulas for the configurational adiabat (i.e., constant excess entropy, sex) lines in a Lennard-Jones (LJ) fluid, one of which is an analytic formula based on a harmonic approximation, which was derived by Heyes et al. [J. Chem. Phys. 159, 224504 (2023)] (analytic isomorph line, AIL). Another is where the density is normalized by the freezing density at that temperature (freezing isomorph line, FIL). It is found that the AIL formula and the average of the freezing density and the melting density ("FMIL") are configurational adiabats at all densities essentially down to the liquid-vapor binodal. The FIL approximation departs from a configurational adiabat in the vicinity of the liquid-vapor binodal close to the freezing line. The self-diffusion coefficient, D, shear viscosity, ηs, and thermal conductivity, λ, in macroscopic reduced units are essentially constant along the AIL and FMIL at all fluid densities and temperatures, but departures from this trend are found along the FIL at high liquid state densities near the liquid-vapor binodal. This supports growing evidence that for simple model systems with no or few internal degrees of freedom, isodynes are lines of constant excess entropy. It is shown that for the LJ fluid, ηs and D can be predicted accurately by an essentially analytic procedure from the high temperature limiting inverse power fluid values (apart from at very low densities), and this is demonstrated quite well also for the experimental argon viscosity.
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
| | - L Costigliola
- Department of Mechanical Engineering, Imperial College London, Exhibition Road, South Kensington, London SW7 2AZ, United Kingdom
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
3
|
Proctor JE, Trachenko K. Generally applicable physics-based equation of state for liquids. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2024; 87:098001. [PMID: 39094591 DOI: 10.1088/1361-6633/ad6a80] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/05/2024] [Accepted: 08/02/2024] [Indexed: 08/04/2024]
Abstract
Physics-based first-principles pressure-volume-temperature equations of state (EOS) exist for solids and gases but not for liquids due to the long-standing fundamental problems involved in liquid theory. Current EOS models that are applicable to liquids and supercritical fluids at liquid-like density under conditions relevant to planetary interiors and industrial processes are complex empirical models with many physically meaningless adjustable parameters. Here, we develop a generally applicable physics-based (GAP) EOS for liquids including supercritical fluids at liquid-like density. The GAP equation is explicit in the internal energy, and hence links the most fundamental macroscopic static property of fluids, the pressure-volume-temperature EOS, to their key microscopic property: the molecular hopping frequency or liquid relaxation time, from which the internal energy can be obtained. We test our GAP equation against available experimental data in several different ways and find good agreement. Our GAP equation, unavoidably and similarly to solid EOS, contains a semi-empirical term giving the energy of the static sample as a function of volume only (EST(V)). Our testing includes studies along isochores, in order to examine the validity of the GAP equation independently of the validity of any function we may choose to utilize forEST(V). The only other adjustable parameter in the equation is the Grüneisen parameter for the fluid. We observe that the GAP equation is similar to the Mie-Grüneisen solid EOS in a wide range of the liquid phase diagram. This similarity is ultimately related to the condensed state of these two phases. On the other hand, the differences between the GAP equation and EOS for gases are fundamental. Finally, we identify the key gaps in the experimental data that need to be filled in to proceed further with the liquid EOS.
Collapse
Affiliation(s)
- J E Proctor
- Materials and Physics Research Group, University of Salford, Manchester M5 4WT, United Kingdom
| | - Kostya Trachenko
- School of Physical and Chemical Sciences, Queen Mary University of London, Mile End Road, London E1 4NS, United Kingdom
| |
Collapse
|
4
|
Sheydaafar Z, Dyre JC, Schrøder TB. Scaling Properties of Liquid Dynamics Predicted from a Single Configuration: Pseudoisomorphs for Harmonic-Bonded Molecules. J Phys Chem B 2024; 128:8054-8064. [PMID: 39110776 DOI: 10.1021/acs.jpcb.4c03982] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2024]
Abstract
Isomorphs are curves in the thermodynamic phase diagram of invariant excess entropy, structure, and dynamics, while pseudoisomorphs are curves of invariant structure and dynamics, but not of the excess entropy. The latter curves have been shown to exist in molecular models with flexible bonds [Olsen, A. E. J. Chem. Phys. 2016, 145, 241103]. We here present three force-based methods to trace out pseudoisomorphs based on a single configuration and test them on the asymmetric dumbbell and 10-bead Lennard-Jones chain models with bonds modeled as harmonic springs. The three methods are based on requiring that particle forces, center-of-mass forces, and torques, respectively, are invariant in reduced units. For each of the two investigated models we identify a method that works well for tracing out pseudoisomorphs, but these methods are not the same. Overall, we find that the more internal degrees of freedom there are in the molecule studied, the less they appear to affect the gross dynamical behavior. Moreover, the "internal" degrees of freedom (including rotation) do not significantly affect the scaling behavior of the dynamical/transport coefficients provided some 'quenching' is performed.
Collapse
Affiliation(s)
- Zahraa Sheydaafar
- 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
| | - Thomas B Schrøder
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
5
|
Douglass IM, Dyre JC, Costigliola L. Complexity Scaling of Liquid Dynamics. PHYSICAL REVIEW LETTERS 2024; 133:068001. [PMID: 39178431 DOI: 10.1103/physrevlett.133.068001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/30/2024] [Accepted: 07/08/2024] [Indexed: 08/25/2024]
Abstract
According to excess-entropy scaling, dynamic properties of liquids like viscosity and diffusion coefficient are determined by the entropy. This link between dynamics and thermodynamics is increasingly studied and of interest also for industrial applications, but hampered by the challenge of calculating entropy efficiently. Utilizing the fact that entropy is basically the Kolmogorov complexity, which can be estimated from optimal compression algorithms [Avinery et al., Phys. Rev. Lett. 123, 178102 (2019)PRLTAO0031-900710.1103/PhysRevLett.123.178102; Martiniani et al., Phys. Rev. X 9, 011031 (2019)PRXHAE2160-330810.1103/PhysRevX.9.011031], we here demonstrate that the diffusion coefficients of four simple liquids follow a quasiuniversal exponential function of the optimal compression length of a single equilibrium configuration. We conclude that "complexity scaling" has the potential to become a useful tool for estimating dynamic properties of any liquid from a single configuration.
Collapse
|
6
|
Jin J, Voth GA. Understanding dynamics in coarse-grained models. IV. Connection of fine-grained and coarse-grained dynamics with the Stokes-Einstein and Stokes-Einstein-Debye relations. J Chem Phys 2024; 161:034114. [PMID: 39012809 DOI: 10.1063/5.0212973] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2024] [Accepted: 06/28/2024] [Indexed: 07/18/2024] Open
Abstract
Applying an excess entropy scaling formalism to the coarse-grained (CG) dynamics of liquids, we discovered that missing rotational motions during the CG process are responsible for artificially accelerated CG dynamics. In the context of the dynamic representability between the fine-grained (FG) and CG dynamics, this work introduces the well-known Stokes-Einstein and Stokes-Einstein-Debye relations to unravel the rotational dynamics underlying FG trajectories, thereby allowing for an indirect evaluation of the effective rotations based only on the translational information at the reduced CG resolution. Since the representability issue in CG modeling limits a direct evaluation of the shear stress appearing in the Stokes-Einstein and Stokes-Einstein-Debye relations, we introduce a translational relaxation time as a proxy to employ these relations, and we demonstrate that these relations hold for the ambient conditions studied in our series of work. Additional theoretical links to our previous work are also established. First, we demonstrate that the effective hard sphere radius determined by the classical perturbation theory can approximate the complex hydrodynamic radius value reasonably well. Furthermore, we present a simple derivation of an excess entropy scaling relationship for viscosity by estimating the elliptical integral of molecules. In turn, since the translational and rotational motions at the FG level are correlated to each other, we conclude that the "entropy-free" CG diffusion only depends on the shape of the reference molecule. Our results and analyses impart an alternative way of recovering the FG diffusion from the CG description by coupling the translational and rotational motions at the hydrodynamic level.
Collapse
Affiliation(s)
- Jaehyeok Jin
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | - Gregory A Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
7
|
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
|
8
|
Hulikal Chakrapani T, Hajibeygi H, Moultos OA, Vlugt TJH. Mutual Diffusivities of Mixtures of Carbon Dioxide and Hydrogen and Their Solubilities in Brine: Insight from Molecular Simulations. Ind Eng Chem Res 2024; 63:10456-10481. [PMID: 38882502 PMCID: PMC11177264 DOI: 10.1021/acs.iecr.4c01078] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2024] [Revised: 05/14/2024] [Accepted: 05/20/2024] [Indexed: 06/18/2024]
Abstract
H2-CO2 mixtures find wide-ranging applications, including their growing significance as synthetic fuels in the transportation industry, relevance in capture technologies for carbon capture and storage, occurrence in subsurface storage of hydrogen, and hydrogenation of carbon dioxide to form hydrocarbons and alcohols. Here, we focus on the thermodynamic properties of H2-CO2 mixtures pertinent to underground hydrogen storage in depleted gas reservoirs. Molecular dynamics simulations are used to compute mutual (Fick) diffusivities for a wide range of pressures (5 to 50 MPa), temperatures (323.15 to 423.15 K), and mixture compositions (hydrogen mole fraction from 0 to 1). At 5 MPa, the computed mutual diffusivities agree within 5% with the kinetic theory of Chapman and Enskog at 423.15 K, albeit exhibiting deviations of up to 25% between 323.15 and 373.15 K. Even at 50 MPa, kinetic theory predictions match computed diffusivities within 15% for mixtures comprising over 80% H2 due to the ideal-gas-like behavior. In mixtures with higher concentrations of CO2, the Moggridge correlation emerges as a dependable substitute for the kinetic theory. Specifically, when the CO2 content reaches 50%, the Moggridge correlation achieves predictions within 10% of the computed Fick diffusivities. Phase equilibria of ternary mixtures involving CO2-H2-NaCl were explored using Gibbs Ensemble (GE) simulations with the Continuous Fractional Component Monte Carlo (CFCMC) technique. The computed solubilities of CO2 and H2 in NaCl brine increased with the fugacity of the respective component but decreased with NaCl concentration (salting out effect). While the solubility of CO2 in NaCl brine decreased in the ternary system compared to the binary CO2-NaCl brine system, the solubility of H2 in NaCl brine increased less in the ternary system compared to the binary H2-NaCl brine system. The cooperative effect of H2-CO2 enhances the H2 solubility while suppressing the CO2 solubility. The water content in the gas phase was found to be intermediate between H2-NaCl brine and CO2-NaCl brine systems. Our findings have implications for hydrogen storage and chemical technologies dealing with CO2-H2 mixtures, particularly where experimental data are lacking, emphasizing the need for reliable thermodynamic data on H2-CO2 mixtures.
Collapse
Affiliation(s)
- Thejas Hulikal Chakrapani
- Reservoir Engineering, Geoscience and Engineering Department, Faculty of Civil Engineering and Geosciences, Delft University of Technology, Delft 2628 CN, The Netherlands
| | - Hadi Hajibeygi
- Reservoir Engineering, Geoscience and Engineering Department, Faculty of Civil Engineering and Geosciences, Delft University of Technology, Delft 2628 CN, The Netherlands
| | - Othonas A Moultos
- Engineering Thermodynamics, Process and Energy Department, Faculty of Mechanical Engineering, Delft University of Technology, Delft 2628 CB, The Netherlands
| | - Thijs J H Vlugt
- Engineering Thermodynamics, Process and Energy Department, Faculty of Mechanical Engineering, Delft University of Technology, Delft 2628 CB, The Netherlands
| |
Collapse
|
9
|
Khrapak SA, Khrapak AG. Freezing density scaling of transport coefficients in the Weeks-Chandler-Andersen fluid. J Chem Phys 2024; 160:134504. [PMID: 38557849 DOI: 10.1063/5.0199310] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2024] [Accepted: 03/13/2024] [Indexed: 04/04/2024] Open
Abstract
It is shown that the transport coefficients (self-diffusion, shear viscosity, and thermal conductivity) of the Weeks-Chandler-Andersen (WCA) fluid along isotherms exhibit a freezing density scaling (FDS). The functional form of this FDS is essentially the same or closely related to those in the Lennard-Jones fluid, hard-sphere fluid, and some liquefied noble gases. This proves that this FDS represents a quasi-universal corresponding state principle for simple classical fluids with steep interactions. Some related aspects, such as a Stokes-Einstein relation without a hydrodynamic diameter and gas-to-liquid dynamical crossover, are briefly discussed. Simple fitting formulas for the transport coefficients of the dense WCA fluid are suggested.
Collapse
Affiliation(s)
- S A Khrapak
- Joint Institute for High Temperatures, Russian Academy of Sciences, 125412 Moscow, Russia
| | - A G Khrapak
- Joint Institute for High Temperatures, Russian Academy of Sciences, 125412 Moscow, Russia
| |
Collapse
|
10
|
Zhang XZ, Shi R, Lu ZY, Qian HJ. Chemically Specific Systematic Coarse-Grained Polymer Model with Both Consistently Structural and Dynamical Properties. JACS AU 2024; 4:1018-1030. [PMID: 38559727 PMCID: PMC10976574 DOI: 10.1021/jacsau.3c00756] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/30/2023] [Revised: 02/23/2024] [Accepted: 02/28/2024] [Indexed: 04/04/2024]
Abstract
The coarse-grained (CG) model serves as a powerful tool for the simulation of polymer systems; its reliability depends on the accurate representation of both structural and dynamical properties. However, strong correlations between structural and dynamical properties on different scales and also a strong memory effect, enforced by chain connectivity between monomers in polymer systems, render developing a chemically specific systematic CG model a formidable task. In this study, we report a systematic CG approach that combines the iterative Boltzmann inversion (IBI) method and the generalized Langevin equation (GLE) dynamics. Structural properties are ensured by using conservative CG potentials derived from the IBI method. To retrieve the correct dynamical properties in the system, we demonstrate that using a combination of a Rouse-type delta function and a time-dependent short-time kernel in the GLE simulation is practically efficient. The former can be used to adjust the long-time diffusion dynamics, and the latter can be reconstructed from an iterative procedure according to the velocity autocorrelation function (ACF) from all-atomistic (AA) simulations. Taking the polystyrene as an example, we show that not only structural properties of radial distribution function, intramolecular bond, and angle distributions can be reproduced but also dynamical properties of mean-square displacement, velocity ACF, and force ACF resulted from our CG model have quantitative agreement with the reference AA model. In addition, reasonable agreements are observed in other collective properties between our GLE-CG model and the AA simulations as well.
Collapse
Affiliation(s)
| | | | - Zhong-Yuan Lu
- State Key Laboratory of Supramolecular
Structure and Materials, Institute of Theoretical Chemistry, College
of Chemistry, Jilin University, Changchun 130021, China
| | - Hu-Jun Qian
- State Key Laboratory of Supramolecular
Structure and Materials, Institute of Theoretical Chemistry, College
of Chemistry, Jilin University, Changchun 130021, China
| |
Collapse
|
11
|
Saric D, Bell IH, Guevara-Carrion G, Vrabec J. Influence of repulsion on entropy scaling and density scaling of monatomic fluids. J Chem Phys 2024; 160:104503. [PMID: 38456532 DOI: 10.1063/5.0196592] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2024] [Accepted: 02/18/2024] [Indexed: 03/09/2024] Open
Abstract
Entropy scaling is applied to the shear viscosity, self-diffusion coefficient, and thermal conductivity of simple monatomic fluids. An extensive molecular dynamics simulation series is performed to obtain these transport properties and the residual entropy of three potential model classes with variable repulsive exponents: n, 6 Mie (n = 9, 12, 15, and 18), Buckingham's exponential-six (α = 12, 14, 18, and 30), and Tang-Toennies (αT = 4.051, 4.275, and 4.600). A wide range of liquid and supercritical gas- and liquid-like states is covered with a total of 1120 state points. Comparisons to equations of state, literature data, and transport property correlations are made. Although the absolute transport property values within a given potential model class may strongly depend on the repulsive exponent, it is found that the repulsive steepness plays a negligible role when entropy scaling is applied. Hence, the plus-scaled transport properties of n, 6 Mie, exponential-six, and Tang-Toennies fluids lie basically on one master curve, which closely corresponds with entropy scaling correlations for the Lennard-Jones fluid. This trend is confirmed by literature data of n, 6 Mie, and exponential-six fluids. Furthermore, entropy scaling holds for state points where the Pearson correlation coefficient R is well below 0.9. The condition R > 0.9 for strongly correlating liquids is thus not necessary for the successful application of entropy scaling, pointing out that isomorph theory may be a part of a more general framework that is behind the success of entropy scaling. Density scaling reveals a strong influence of the repulsive exponent on this particular approach.
Collapse
Affiliation(s)
- Denis Saric
- Thermodynamics, Technical University of Berlin, Ernst-Reuter-Platz 1, 10587 Berlin, Germany
| | - Ian H Bell
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, USA
| | | | - Jadran Vrabec
- Thermodynamics, Technical University of Berlin, Ernst-Reuter-Platz 1, 10587 Berlin, Germany
| |
Collapse
|
12
|
Huang Y, Widom M. Entropy approximations for simple fluids. Phys Rev E 2024; 109:034130. [PMID: 38632758 DOI: 10.1103/physreve.109.034130] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2023] [Accepted: 02/29/2024] [Indexed: 04/19/2024]
Abstract
Liquid state entropy formulas based on configurational probability distributions are examined for Lennard-Jones fluids across a range temperatures and densities. These formulas are based on expansions of the entropy in a series of n-body distribution functions. We focus on two special cases. One, which we term the "perfect gas" series, starts with the entropy of an ideal gas; the other, which we term the "dense liquid" series, removes a many-body contribution from the ideal gas entropy and reallocates it among the subsequent n-body terms. We show that the perfect gas series is most accurate at low density, while the dense liquid series is most accurate at high density. We propose empirical interpolation methods that are capable of connecting the two series and giving consistent predictions in most situations.
Collapse
Affiliation(s)
- Yang Huang
- Physics Department, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA; University of Science and Technology of China, Hefei 230026, China; and Suzhou Institute for Advanced Research, University of Science and Technology of China, Suzhou 215213, China
| | - Michael Widom
- Physics Department, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| |
Collapse
|
13
|
Shahzad A, Yang F, Steffen J, Neiss C, Panchenko A, Goetz K, Vogel C, Weisser M, Embs JP, Petry W, Lohstroh W, Görling A, Goychuk I, Unruh T. Atomic diffusion in liquid gallium and gallium-nickel alloys probed by quasielastic neutron scattering and molecular dynamic simulations. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2024; 36:175403. [PMID: 38224622 DOI: 10.1088/1361-648x/ad1e9f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/08/2023] [Accepted: 01/15/2024] [Indexed: 01/17/2024]
Abstract
The atomic mobility in liquid pure gallium and a gallium-nickel alloy with 2 at% of nickel is studied experimentally by incoherent quasielastic neutron scattering. The integral diffusion coefficients for all-atom diffusion are derived from the experimental data at different temperatures. DFT-basedab-initiomolecular dynamics (MD) is used to find numerically the diffusion coefficient of liquid gallium at different temperatures, and numerical theory results well agree with the experimental findings at temperatures below 500 K. Machine learning force fields derived fromab-initiomolecular dynamics (AIMD) overestimate within a small 6% error the diffusion coefficient of pure gallium within the genuine AIMD. However, they better agree with experiment for pure gallium and enable the numerical finding of the diffusion coefficient of nickel in the considered melted alloy along with the diffusion coefficient of gallium and integral diffusion coefficient, that agrees with the corresponding experimental values within the error bars. The temperature dependence of the gallium diffusion coefficientDGa(T)follows the Arrhenius law experimentally for all studied temperatures and below 500 K also in the numerical simulations. However,DGa(T)can be well described alternatively by an Einstein-Stokes dependence with the metallic liquid viscosity following the Arrhenius law, especially for the MD simulation results at all studied temperatures. Moreover, a novel variant of the excess entropy scaling theory rationalized our findings for gallium diffusion. Obtained values of the Arrhenius activation energies are profoundly different in the competing theoretical descriptions, which is explained by different temperature-dependent prefactors in the corresponding theories. The diffusion coefficient of gallium is significantly reduced (at the same temperature) in a melted alloy with natural nickel, even at a tiny 2 at% concentration of nickel, as compared with its pure gallium value. This highly surprising behavior contradicts the existing excess entropy scaling theories and opens a venue for further research.
Collapse
Affiliation(s)
- A Shahzad
- Institute for Crystallography and Structural Physics, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Staudtstraße 3, Erlangen 91058, Germany
- Interdisciplinary Center for Nanostructured Films (IZNF) and Center for Nanoanalysis and Electron Microscopy (CENEM), Cauerstraße 3, Erlangen 91058, Germany
- Institute for Material Science, University of Stuttgart, Heisenbergstr. 3, 70569 Stuttgart, Germany
| | - F Yang
- Institut für Materialphysik im Weltraum, Deutsches Zentrum für Luft- und Raumfahrt (DLR), 51170 Köln, Germany
| | - J Steffen
- Chair of Theoretical Chemistry, FAU, 91058 Erlangen, Germany
| | - C Neiss
- Chair of Theoretical Chemistry, FAU, 91058 Erlangen, Germany
| | - A Panchenko
- Institute for Crystallography and Structural Physics, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Staudtstraße 3, Erlangen 91058, Germany
- Interdisciplinary Center for Nanostructured Films (IZNF) and Center for Nanoanalysis and Electron Microscopy (CENEM), Cauerstraße 3, Erlangen 91058, Germany
| | - K Goetz
- Institute for Crystallography and Structural Physics, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Staudtstraße 3, Erlangen 91058, Germany
- Interdisciplinary Center for Nanostructured Films (IZNF) and Center for Nanoanalysis and Electron Microscopy (CENEM), Cauerstraße 3, Erlangen 91058, Germany
| | - C Vogel
- Institute for Crystallography and Structural Physics, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Staudtstraße 3, Erlangen 91058, Germany
- Interdisciplinary Center for Nanostructured Films (IZNF) and Center for Nanoanalysis and Electron Microscopy (CENEM), Cauerstraße 3, Erlangen 91058, Germany
| | - M Weisser
- Institute for Crystallography and Structural Physics, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Staudtstraße 3, Erlangen 91058, Germany
- Interdisciplinary Center for Nanostructured Films (IZNF) and Center for Nanoanalysis and Electron Microscopy (CENEM), Cauerstraße 3, Erlangen 91058, Germany
| | - J P Embs
- Laboratory for Neutron Scattering and Imaging, Paul Scherrer Institut (PSI), CH-5232 Villigen, Switzerland
| | - W Petry
- Physics Department, Technical University of Munich, James-Franck-Str. 1, 85747 Garching, Germany
| | - W Lohstroh
- Research Neutron Source Heinz Maier-Leibnitz (FRM II), Technical University of Munich, Lichtenbergstr. 1, 85748 Garching, Germany
| | - A Görling
- Chair of Theoretical Chemistry, FAU, 91058 Erlangen, Germany
| | - I Goychuk
- Institute for Crystallography and Structural Physics, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Staudtstraße 3, Erlangen 91058, Germany
- Interdisciplinary Center for Nanostructured Films (IZNF) and Center for Nanoanalysis and Electron Microscopy (CENEM), Cauerstraße 3, Erlangen 91058, Germany
| | - T Unruh
- Institute for Crystallography and Structural Physics, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Staudtstraße 3, Erlangen 91058, Germany
- Interdisciplinary Center for Nanostructured Films (IZNF) and Center for Nanoanalysis and Electron Microscopy (CENEM), Cauerstraße 3, Erlangen 91058, Germany
| |
Collapse
|
14
|
Knudsen PA, Heyes DM, Niss K, Dini D, Bailey NP. Invariant dynamics in a united-atom model of an ionic liquid. J Chem Phys 2024; 160:034503. [PMID: 38230811 DOI: 10.1063/5.0177373] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2023] [Accepted: 12/11/2023] [Indexed: 01/18/2024] Open
Abstract
We study a united-atom model of the ionic liquid 1-butyl-1-methylpyrrolidinium bis(trifluoromethyl)sulfonylamide to determine to what extent there exist curves in the phase diagram along which the microscopic dynamics are invariant when expressed in dimensionless, or reduced, form. The initial identification of these curves, termed isodynes, is made by noting that contours of reduced shear viscosity and reduced self-diffusion coefficient coincide to a good approximation. Choosing specifically the contours of reduced viscosity as nominal isodynes, further simulations were carried out for state points on these, and other aspects of dynamics were investigated to study their degree of invariance. These include the mean-squared displacement, shear-stress autocorrelation function, and various rotational correlation functions. These were invariant to a good approximation, with the main exception being rotations of the anion about its long axis. The dynamical features that are invariant have in common that they are aspects that would be relevant for a coarse-grained description of the system; specifically, removing the most microscopic degrees of freedom in principle leads to a simplification of the potential energy landscape, which allows for the existence of isodynes.
Collapse
Affiliation(s)
- Peter A Knudsen
- "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, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - Kristine Niss
- "Glass and Time," IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Daniele Dini
- Department of Mechanical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - Nicholas P Bailey
- "Glass and Time," IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
15
|
Chaparro G, Müller EA. Simulation and Data-Driven Modeling of the Transport Properties of the Mie Fluid. J Phys Chem B 2024; 128:551-566. [PMID: 38181201 PMCID: PMC10801693 DOI: 10.1021/acs.jpcb.3c06813] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2023] [Revised: 12/14/2023] [Accepted: 12/15/2023] [Indexed: 01/07/2024]
Abstract
This work reports the computation and modeling of the self-diffusivity (D*), shear viscosity (η*), and thermal conductivity (κ*) of the Mie fluid. The transport properties were computed using equilibrium molecular dynamics simulations for the Mie fluid with repulsive exponents (λr) ranging from 7 to 34 and at a fixed attractive exponent (λa) of 6 over the whole fluid density (ρ*) range and over a wide temperature (T*) range. The computed database consists of 17,212, 14,288, and 13,099 data points for self-diffusivity, shear viscosity, and thermal conductivity, respectively. The database is successfully validated against published simulation data. The above-mentioned transport properties are correlated using artificial neural networks (ANNs). Two modeling approaches were tested: a semiempirical formulation based on entropy scaling and an empirical formulation based on density and temperature as input variables. For the former, it was found that a unique formulation based on entropy scaling does not yield satisfactory results over the entire density range due to a divergent and incorrect scaling of the transport properties at low densities. For the latter empirical modeling approach, it was found that regularizing the data, e.g., modeling ρ*D* instead of D*, ln η* instead of η*, and ln κ* instead of κ*, as well as using the inverse of the temperature as an input feature, helps to ease the interpolation efforts of the artificial neural networks. The trained ANNs can model seen and unseen data over a wide range of density and temperature. Ultimately, the ANNs can be used alongside equations of state to regress effective force field parameters from volumetric and transport data.
Collapse
Affiliation(s)
- Gustavo Chaparro
- Department of Chemical Engineering,
Sargent Centre for Process Systems Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, U.K.
| | - Erich A. Müller
- Department of Chemical Engineering,
Sargent Centre for Process Systems Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, U.K.
| |
Collapse
|
16
|
Angelis D, Sofos F, Karakasidis TE. Reassessing the transport properties of fluids: A symbolic regression approach. Phys Rev E 2024; 109:015105. [PMID: 38366535 DOI: 10.1103/physreve.109.015105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2023] [Accepted: 01/04/2024] [Indexed: 02/18/2024]
Abstract
The viscosity and thermal conductivity coefficients of the Lennard-Jones fluid are extracted through symbolic regression (SR) techniques from data derived from simulations at the atomic scale. This data-oriented approach provides closed form relations that achieve fine accuracy when compared to well-established theoretical, empirical, or approximate equations, fully transparent, with small complexity and high interpretability. The novelty is further outlined by suggesting analytical expressions for estimating fluid transport properties across the whole phase space, from a dilute gas to a dense liquid, by considering only two macroscopic properties (density and temperature). In such expressions, the underlying physical mechanisms are reflected, while, at the same time, it can be a computationally efficient alternative to costly in time and size first principle and/or molecular dynamics simulations.
Collapse
Affiliation(s)
- Dimitrios Angelis
- Condensed Matter Physics Laboratory, Department of Physics, University of Thessaly, Lamia 35100, Greece
| | - Filippos Sofos
- Condensed Matter Physics Laboratory, Department of Physics, University of Thessaly, Lamia 35100, Greece
| | - Theodoros E Karakasidis
- Condensed Matter Physics Laboratory, Department of Physics, University of Thessaly, Lamia 35100, Greece
| |
Collapse
|
17
|
Jawed A, Khan MN, Khan NA, Hakeem MA, Khan P. Thermodynamic Insights into Variation in Thermomechanical and Physical Properties of Isotactic Polypropylene: Effect of Shear and Cooling Rates. ACS OMEGA 2023; 8:36775-36788. [PMID: 37841128 PMCID: PMC10569009 DOI: 10.1021/acsomega.3c03378] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/15/2023] [Accepted: 08/24/2023] [Indexed: 10/17/2023]
Abstract
In order to elucidate the effect of shear and cooling process on structural, thermomechanical, and physical properties of polymer melt, excess entropy, a thermodynamic quantity is calculated from radial distribution function generated from equilibrated parts of the molecular simulation trajectories. The structural properties are calculated, which includes the density of polypropylene melt, end to end distance, radius of gyration of the polypropylene polymer chain, and monomer-monomer radial distribution function. Non-equilibrium molecular dynamics simulation was employed to investigate the role of the applied shear rate on the properties of polypropylene. Furthermore, a range of cooling rates were employed to cool the melt. Thermomechanical properties, such as Young's modulus, and physical properties, such as glass transition temperature, were determined for different cases. Results showed that slow cooling and high shear substantially improved the Young's modulus and glass transition temperature of the i-PP. Furthermore, a two-body contribution to the excess entropy was used to elucidate the structure-property relationships in the polymer melt as well as the glassy state and the dependence of shear and cooling rate on these properties. We have used the Rosenfeld excess entropy-viscosity relationship to calculate the viscous behavior of the polymer under a steady shear condition.
Collapse
Affiliation(s)
- Ahmad
S. Jawed
- Department of Chemical Engineering, Aligarh Muslim University, Aligarh 202002, India
| | - Mohd Nasir Khan
- Department of Chemical Engineering, Aligarh Muslim University, Aligarh 202002, India
| | - Naseem A. Khan
- Department of Chemical Engineering, Aligarh Muslim University, Aligarh 202002, India
| | - Mohammed A. Hakeem
- Department of Chemical Engineering, Aligarh Muslim University, Aligarh 202002, India
| | - Parvez Khan
- Department of Chemical Engineering, Aligarh Muslim University, Aligarh 202002, India
| |
Collapse
|
18
|
Sorkin B, Diamant H, Ariel G. Universal Relation between Entropy and Kinetics. PHYSICAL REVIEW LETTERS 2023; 131:147101. [PMID: 37862659 DOI: 10.1103/physrevlett.131.147101] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/15/2023] [Accepted: 08/25/2023] [Indexed: 10/22/2023]
Abstract
Relating thermodynamic and kinetic properties is a conceptual challenge with many practical benefits. Here, based on first principles, we derive a rigorous inequality relating the entropy and the dynamic propagator of particle configurations. It is universal and applicable to steady states arbitrarily far from thermodynamic equilibrium. Applying the general relation to diffusive dynamics yields a relation between the entropy and the (normal or anomalous) diffusion coefficient. The relation can be used to obtain useful bounds for the late-time diffusion coefficient from the calculated steady-state entropy or, conversely, to estimate the entropy based on measured diffusion coefficients. We demonstrate the validity and usefulness of the relation through several examples and discuss its broad range of applications, in particular, for systems far from equilibrium.
Collapse
Affiliation(s)
- Benjamin Sorkin
- School of Chemistry and Center for Physics and Chemistry of Living Systems, Tel Aviv University, 69978 Tel Aviv, Israel
| | - Haim Diamant
- School of Chemistry and Center for Physics and Chemistry of Living Systems, Tel Aviv University, 69978 Tel Aviv, Israel
| | - Gil Ariel
- Department of Mathematics, Bar-Ilan University, 52000 Ramat Gan, Israel
| |
Collapse
|
19
|
Liszka K, Grzybowski A, Grzybowska K, Koperwas K, Paluch M. Entropy Scaling of Molecular Dynamics in a Prototypical Anisotropic Model near the Glass Transition. J Phys Chem B 2023. [PMID: 37257018 DOI: 10.1021/acs.jpcb.3c02429] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/02/2023]
Abstract
Dynamics and thermodynamics of molecular systems in the vicinity of the boundary between thermodynamically nonequilibrium glassy and metastable supercooled liquid states are still incompletely explored and their theoretical and simulation models are imperfect despite many previous efforts. Among them, the role of total system entropy, configurational entropy, and excess entropy in the temperature-pressure or temperature-density dependence of global molecular dynamics (MD) timescale relevant to the glass transition is an essential topic intensively studied for over half of a century. By exploiting a well-known simple ellipsoidal model recently successfully applied to simulate the supercooled liquid state and the glass transition, we gain a new insight into the different views on the relationship between entropy and relaxation dynamics of glass formers, showing the molecular grounds for the entropy scaling of global MD timescale. Our simulations in the anisotropic model of supercooled liquid, which involves only translational and rotational degrees of freedom, give evidence that the total system entropy is sufficient to scale global MD timescale. It complies with the scaling effect on relaxation dynamics exerted by the configurational entropy defined as the total entropy diminished by vibrational contributions, which was earlier discovered for measurement data collected near the glass transition. Moreover, we argue that such a scaling behavior is not possible to achieve by using the excess entropy that is in excess of the ideal gas entropy, which is contrary to the results earlier suggested within the framework of simple isotropic models of supercooled liquids. Thus, our findings also warn against an excessive reliance on isotropic models in theoretical interpretations of molecular phenomena, despite their simplicity and popularity, because they may reflect improperly various physicochemical properties of glass formers.
Collapse
Affiliation(s)
- Karol Liszka
- Institute of Physics, University of Silesia in Katowice, ul. 75 Pułku Piechoty 1, 41-500 Chorzów, Poland
| | - Andrzej Grzybowski
- Institute of Physics, University of Silesia in Katowice, ul. 75 Pułku Piechoty 1, 41-500 Chorzów, Poland
| | - Katarzyna Grzybowska
- Institute of Physics, University of Silesia in Katowice, ul. 75 Pułku Piechoty 1, 41-500 Chorzów, Poland
| | - Kajetan Koperwas
- Institute of Physics, University of Silesia in Katowice, ul. 75 Pułku Piechoty 1, 41-500 Chorzów, Poland
| | - Marian Paluch
- Institute of Physics, University of Silesia in Katowice, ul. 75 Pułku Piechoty 1, 41-500 Chorzów, Poland
| |
Collapse
|
20
|
Farshad M, DelloStritto MJ, Suma A, Carnevale V. Detecting Liquid-Liquid Phase Separations Using Molecular Dynamics Simulations and Spectral Clustering. J Phys Chem B 2023; 127:3682-3689. [PMID: 37053472 DOI: 10.1021/acs.jpcb.3c00805] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/15/2023]
Abstract
A stringent test of the accuracy of empirical force fields is reproducing the phase diagram of bulk phases and mixtures. Exploring the phase diagram of mixtures requires the detection of phase boundaries and critical points. In contrast to most solid-liquid transitions, in which a global order parameter (average density) can be used to discriminate between two phases, some demixing transitions entail relatively subtle changes in the local environment of each molecule. In such cases, finite sampling errors and finite-size effects make the identification of trends in local order parameters extremely challenging. Here we analyze one such example, namely a methanol/hexane mixture, and compute several local and global structural properties. We simulate the system at various temperatures and study the structural changes associated with demixing. We show that despite a seemingly continuous transformation between mixed and demixed states, the topological properties of the H-bond network change abruptly as the system crosses the demixing line. In particular, by using spectral clustering, we show that the distribution of cluster sizes develops a fat tail (as expected from percolation theory) in the vicinity of the critical point. We illustrate a simple criterion to identify this behavior, which results from the emergence of large system-spanning clusters from a collection of aggregates. We further tested the spectral clustering analysis on a Lennard-Jones system as a standard example of a system with no H-bonds, and also, in this case, we were able to detect the demixing transition.
Collapse
Affiliation(s)
- Mohsen Farshad
- Institute for Computational Molecular Science, Temple University, Philadelphia, Pennsylvania 19122, United States
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Mark J DelloStritto
- Institute for Computational Molecular Science, Temple University, Philadelphia, Pennsylvania 19122, United States
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Antonio Suma
- Institute for Computational Molecular Science, Temple University, Philadelphia, Pennsylvania 19122, United States
- Dipartimento di Fisica, Università di Bari, 70121 Bari, Italy
| | - Vincenzo Carnevale
- Institute for Computational Molecular Science, Temple University, Philadelphia, Pennsylvania 19122, United States
- Institute for Genomics and Evolutionary Medicine, Temple University, Philadelphia, Pennsylvania 19122, United States
- Department of Biology, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
21
|
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
|
22
|
Khrapak S. Bridgman formula for the thermal conductivity of atomic and molecular liquids. J Mol Liq 2023. [DOI: 10.1016/j.molliq.2023.121786] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/08/2023]
|
23
|
Mehri S, Dyre JC, Ingebrigtsen TS. Hidden scale invariance in the Gay-Berne model. II. Smectic-B phase. Phys Rev E 2023; 107:044702. [PMID: 37198818 DOI: 10.1103/physreve.107.044702] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2023] [Accepted: 04/10/2023] [Indexed: 05/19/2023]
Abstract
This paper complements a previous study of the isotropic and nematic phases of the Gay-Berne liquid-crystal model [Mehri et al., Phys. Rev. E 105, 064703 (2022)2470-004510.1103/PhysRevE.105.064703] with a study of its smectic-B phase found at high density and low temperatures. We find also in this phase strong correlations between the virial and potential-energy thermal fluctuations, reflecting hidden scale invariance and implying the existence of isomorphs. The predicted approximate isomorph invariance of the physics is confirmed by simulations of the standard and orientational radial distribution functions, the mean-square displacement as a function of time, and the force, torque, velocity, angular velocity, and orientational time-autocorrelation functions. The regions of the Gay-Berne model that are relevant for liquid-crystal experiments can thus fully be simplified via the isomorph theory.
Collapse
Affiliation(s)
- Saeed Mehri
- 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
| | - Trond S Ingebrigtsen
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
24
|
Wang Q, Zhao C, Yao Z, Wang J, Wu F, Kumar SGH, Ganapathy S, Eustace S, Bai X, Li B, Lu J, Wagemaker M. Entropy-Driven Liquid Electrolytes for Lithium Batteries. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2023; 35:e2210677. [PMID: 36718916 DOI: 10.1002/adma.202210677] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/17/2022] [Revised: 01/23/2023] [Indexed: 05/17/2023]
Abstract
Developing liquid electrolytes with higher kinetics and enhanced interphase stability is one of the key challenges for lithium batteries. However, the poor solubility of lithium salts in solvents sets constraints that compromises the electrolyte properties. Here, it is shown that introducing multiple salts to form a high-entropy solution, alters the solvation structure, which can be used to raise the solubility of specific salts and stabilize electrode-electrolyte interphases. The prepared high-entropy electrolytes significantly enhance the cycling and rate performance of lithium batteries. For lithium-metal anodes the reversibility exceeds 99%, which extends the cycle life of batteries even under aggressive cycling conditions. For commercial batteries, combining a graphite anode with a LiNi0.8 Co0.1 Mn0.1 O2 cathode, more than 1000 charge-discharge cycles are achieved while maintaining a capacity retention of more than 90%. These performance improvements with respect to regular electrolytes are rationalized by the unique features of the solvation structure in high-entropy electrolytes. The weaker solvation interaction induced by the higher disorder results in improved lithium-ion kinetics, and the altered solvation composition leads to stabilized interphases. Finally, the high-entropy, induced by the presence of multiple salts, enables a decrease in melting temperature of the electrolytes and thus enables lower battery operation temperatures without changing the solvents.
Collapse
Affiliation(s)
- Qidi Wang
- Department of Radiation Science and Technology, Delft University of Technology, Delft, 2629JB, The Netherlands
| | - Chenglong Zhao
- Department of Radiation Science and Technology, Delft University of Technology, Delft, 2629JB, The Netherlands
| | - Zhenpeng Yao
- The State Key Laboratory of Metal Matrix Composites, School of Materials Science and Engineering, Center of Hydrogen Science, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Jianlin Wang
- State Key Laboratory for Surface Physics, Institute of Physics, Chinese Academy of Sciences, Beijing, 100190, China
| | - Fangting Wu
- Shenzhen Key Laboratory on Power Battery Safety and Shenzhen Geim Graphene Center, School of Shenzhen International Graduate, Tsinghua University, Guangdong, 518055, China
| | - Sai Govind Hari Kumar
- Department of Chemistry and Computer Science, University of Toronto, Ontario, M5S 3H6, Canada
| | - Swapna Ganapathy
- Department of Radiation Science and Technology, Delft University of Technology, Delft, 2629JB, The Netherlands
| | - Stephen Eustace
- Department of Biotechnology, Delft University of Technology, Delft, 2629HZ, The Netherlands
| | - Xuedong Bai
- State Key Laboratory for Surface Physics, Institute of Physics, Chinese Academy of Sciences, Beijing, 100190, China
| | - Baohua Li
- Shenzhen Key Laboratory on Power Battery Safety and Shenzhen Geim Graphene Center, School of Shenzhen International Graduate, Tsinghua University, Guangdong, 518055, China
| | - Jun Lu
- College of Chemical and Biological Engineering, Zhejiang University, Hangzhou, 310027, China
| | - Marnix Wagemaker
- Department of Radiation Science and Technology, Delft University of Technology, Delft, 2629JB, The Netherlands
| |
Collapse
|
25
|
Xu WS, Sun ZY. A Thermodynamic Perspective on Polymer Glass Formation. CHINESE JOURNAL OF POLYMER SCIENCE 2023. [DOI: 10.1007/s10118-023-2951-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/29/2023]
|
26
|
Ingebrigtsen TS, Dyre JC. Even Strong Energy Polydispersity Does Not Affect the Average Structure and Dynamics of Simple Liquids. J Phys Chem B 2023; 127:2837-2846. [PMID: 36926946 DOI: 10.1021/acs.jpcb.3c00346] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/18/2023]
Abstract
Size-polydisperse liquids have become standard models for avoiding crystallization, thereby enabling studies of supercooled liquids and glasses formed, e.g., by colloidal systems. Purely energy-polydisperse liquids have been studied much less, but provide an interesting alternative. We here study numerically the difference in structure and dynamics obtained by introducing these two kinds of polydispersity into systems of particles interacting via the Lennard-Jones and EXP pair potentials. To a very good approximation, the average pair structure and dynamics are unchanged even for strong energy polydispersity, which is not the case for size-polydisperse systems. When the system at extreme energy polydispersity undergoes a continuous phase separation into lower and higher particle-energy regions whose structure and dynamics are different from the average, the average structure and dynamics are still virtually the same as for the monodisperse system. Our findings are consistent with the fact that the distribution of forces on the individual particles do not change when energy polydispersity is introduced, while they do change in the case of size polydispersity. A theoretical explanation remains to be found, however.
Collapse
Affiliation(s)
- Trond S Ingebrigtsen
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, Postbox 260, DK-4000 Roskilde, Denmark
| | - Jeppe C Dyre
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, Postbox 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
27
|
Dmitryuk NA, Mistryukova LA, Kryuchkov NP, Khrapak SA, Yurchenko SO. Diffusion mobility increases linearly on liquid binodals above triple point. Sci Rep 2023; 13:2815. [PMID: 36797382 PMCID: PMC9935557 DOI: 10.1038/s41598-022-26390-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2022] [Accepted: 12/14/2022] [Indexed: 02/18/2023] Open
Abstract
Self-diffusion in fluids has been thoroughly studied numerically, but even for simple liquids just a few scaling relationships are known. Relations between diffusion, excitation spectra, and character of the interparticle interactions remain poorly understood. Here, we show that diffusion mobility of particles in simple fluids increases linearly on the liquid branch of the liquid-gas binodal, from the triple point almost up to the critical point. With molecular dynamics simulations, we considered bulk systems of particles interacting via a generalised Lennard-Jones potential, as well as ethane. Using a two-oscillator model for the analysis of excitations, we observed that the mobility (inverse diffusion) coefficient on the liquid-gas binodal increases linearly above the triple point until the dispersion of high-frequency spectra has a solid-like (oscillating) shape. In terms of a separate mode analysis (of longitudinal and transverse modes), this corresponds to crossed modes in the intermediate range of wavenumbers q, between the hydrodynamic regime (small q) and the regime of individual particle motion (large q). The results should be interesting for a broad community in physics and chemistry of fluids, since self-diffusion is among the most fundamental transport phenomena, important for prospective chemical technologies, micro-, nanofluidics, and biotechnologies.
Collapse
Affiliation(s)
- Nikita A. Dmitryuk
- grid.61569.3d0000 0001 0405 5955Bauman Moscow State Technical University, 2nd Baumanskaya Street 5, Moscow, Russia 105005
| | - Lucia A. Mistryukova
- grid.61569.3d0000 0001 0405 5955Bauman Moscow State Technical University, 2nd Baumanskaya Street 5, Moscow, Russia 105005
| | - Nikita P. Kryuchkov
- grid.61569.3d0000 0001 0405 5955Bauman Moscow State Technical University, 2nd Baumanskaya Street 5, Moscow, Russia 105005
| | - Sergey A. Khrapak
- grid.61569.3d0000 0001 0405 5955Bauman Moscow State Technical University, 2nd Baumanskaya Street 5, Moscow, Russia 105005
| | - Stanislav O. Yurchenko
- grid.61569.3d0000 0001 0405 5955Bauman Moscow State Technical University, 2nd Baumanskaya Street 5, Moscow, Russia 105005
| |
Collapse
|
28
|
Saw S, Costigliola L, Dyre JC. Configurational temperature in active matter. II. Quantifying the deviation from thermal equilibrium. Phys Rev E 2023; 107:024610. [PMID: 36932493 DOI: 10.1103/physreve.107.024610] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2022] [Accepted: 01/23/2023] [Indexed: 06/18/2023]
Abstract
This paper proposes using the configurational temperature T_{conf} for quantifying how far an active-matter system is from thermal equilibrium. We measure this "distance" by the ratio of the systemic temperature T_{s} to T_{conf}, where T_{s} is the canonical-ensemble temperature for which the average potential energy is equal to that of the active-matter system. T_{conf} is "local" in the sense that it is the average of a function, which depends only on how the potential energy varies in the vicinity of a given configuration. In contrast, T_{s} is a global quantity. The quantity T_{s}/T_{conf} is straightforward to evaluate in a computer simulation; equilibrium simulations in conjunction with a single steady-state active-matter configuration are enough to determine T_{s}/T_{conf}. We validate the suggestion that T_{s}/T_{conf} quantifies the deviation from thermal equilibrium by data for the radial distribution function of the 3D Kob-Andersen and 2D Yukawa active-matter models with active Ornstein-Uhlenbeck and active Brownian Particle dynamics. Moreover, we show that T_{s}/T_{conf}, structure, and dynamics of the homogeneous phase are all approximately invariant along the motility-induced phase separation boundary in the phase diagram of the 2D Yukawa model. The measure T_{s}/T_{conf} is not limited to active matter and can be used for quantifying how far any system involving a potential-energy function, e.g., a driven Hamiltonian system, is from thermal equilibrium.
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
|
29
|
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
|
30
|
Abstract
High-entropy alloys/compounds have large configurational entropy by introducing multiple components, showing improved functional properties that exceed those of conventional materials. However, how increasing entropy impacts the thermodynamic/kinetic properties in liquids that are ambiguous. Here we show this strategy in liquid electrolytes for rechargeable lithium batteries, demonstrating the substantial impact of raising the entropy of electrolytes by introducing multiple salts. Unlike all liquid electrolytes so far reported, the participation of several anionic groups in this electrolyte induces a larger diversity in solvation structures, unexpectedly decreasing solvation strengths between lithium ions and solvents/anions, facilitating lithium-ion diffusivity and the formation of stable interphase passivation layers. In comparison to the single-salt electrolytes, a low-concentration dimethyl ether electrolyte with four salts shows an enhanced cycling stability and rate capability. These findings, rationalized by the fundamental relationship between entropy-dominated solvation structures and ion transport, bring forward high-entropy electrolytes as a composition-rich and unexplored space for lithium batteries and beyond.
Collapse
|
31
|
Jin J, Schweizer KS, Voth GA. Understanding dynamics in coarse-grained models. I. Universal excess entropy scaling relationship. J Chem Phys 2023; 158:034103. [PMID: 36681649 DOI: 10.1063/5.0116299] [Citation(s) in RCA: 15] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Coarse-grained (CG) models facilitate an efficient exploration of complex systems by reducing the unnecessary degrees of freedom of the fine-grained (FG) system while recapitulating major structural correlations. Unlike structural properties, assessing dynamic properties in CG modeling is often unfeasible due to the accelerated dynamics of the CG models, which allows for more efficient structural sampling. Therefore, the ultimate goal of the present series of articles is to establish a better correspondence between the FG and CG dynamics. To assess and compare dynamical properties in the FG and the corresponding CG models, we utilize the excess entropy scaling relationship. For Paper I of this series, we provide evidence that the FG and the corresponding CG counterpart follow the same universal scaling relationship. By carefully reviewing and examining the literature, we develop a new theory to calculate excess entropies for the FG and CG systems while accounting for entropy representability. We demonstrate that the excess entropy scaling idea can be readily applied to liquid water and methanol systems at both the FG and CG resolutions. For both liquids, we reveal that the scaling exponents remain unchanged from the coarse-graining process, indicating that the scaling behavior is universal for the same underlying molecular systems. Combining this finding with the concept of mapping entropy in CG models, we show that the missing entropy plays an important role in accelerating the CG dynamics.
Collapse
Affiliation(s)
- Jaehyeok Jin
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - Kenneth S Schweizer
- Department of Material Science, Department of Chemistry, Department of Chemical and Biomolecular Engineering, and Materials Research Laboratory, University of Illinois, Urbana, Illinois 61801, USA
| | - Gregory A Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
32
|
Young JM, Bell IH, Harvey AH. Entropy scaling of viscosity for molecular models of molten salts. J Chem Phys 2023; 158:024502. [PMID: 36641388 DOI: 10.1063/5.0127250] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Entropy scaling relates dynamic and thermodynamic properties by reducing the viscosity to a function of only the residual entropy. Molecular simulations are used to investigate the entropy scaling of the viscosity of three models of sodium chloride and five monovalent salts. Even though the correlation between the potential energy and the virial is weak, entropy scaling applies at liquid densities for all models and salts investigated. At lower densities, entropy scaling breaks down due to the formation of ion pairs and chains. Entropy scaling can be used to develop more extendable correlations for the dynamic properties of molten salts.
Collapse
Affiliation(s)
- Jeffrey M Young
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, USA
| | - Ian H Bell
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, USA
| | - Allan H Harvey
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, USA
| |
Collapse
|
33
|
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
|
34
|
Schrøder TB. Predicting Scaling Properties from a Single Fluid Configuration. PHYSICAL REVIEW LETTERS 2022; 129:245501. [PMID: 36563245 DOI: 10.1103/physrevlett.129.245501] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/14/2021] [Revised: 04/25/2022] [Accepted: 11/01/2022] [Indexed: 06/17/2023]
Abstract
Time-dependent dynamical properties of a fluid cannot be estimated directly from a single configuration without performing a simulation. Here, however, we present a method that predicts the scaling properties of both structure and dynamics from a single configuration. The method is demonstrated to work well for the Lennard-Jones fluid as well as the viscous Kob-Andersen Lennard-Jones mixture, both in and out of equilibrium. The method is conceptually simple and easy to implement and, thus, should become a standard tool in the study of scaling properties of fluids and liquids.
Collapse
Affiliation(s)
- Thomas B Schrøder
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
35
|
Thermal conductivity prediction of pure refrigerants and mixtures based on entropy-scaling concept. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120568] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
36
|
Nikitiuk B, Salikova D, Kondratyuk N, Pisarev V. Pair entropy and universal viscosity scaling for molecular systems via molecular dynamics simulations. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120714] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
37
|
Zeng F, Wan R, Xiao Y, Song F, Peng C, Liu H. Predicting the Self-Diffusion Coefficient of Liquids Based on Backpropagation Artificial Neural Network: A Quantitative Structure–Property Relationship Study. Ind Eng Chem Res 2022. [DOI: 10.1021/acs.iecr.2c03342] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Fazhan Zeng
- School of Chemistry & Molecular Engineering, East China University of Science and Technology, Shanghai200237, China
| | - Ren Wan
- School of Chemistry & Molecular Engineering, East China University of Science and Technology, Shanghai200237, China
| | - Yongjun Xiao
- School of Chemistry & Molecular Engineering, East China University of Science and Technology, Shanghai200237, China
| | - Fan Song
- School of Chemistry & Molecular Engineering, East China University of Science and Technology, Shanghai200237, China
| | - Changjun Peng
- School of Chemistry & Molecular Engineering, East China University of Science and Technology, Shanghai200237, China
| | - Honglai Liu
- School of Chemistry & Molecular Engineering, East China University of Science and Technology, Shanghai200237, China
| |
Collapse
|
38
|
Crossover Residual Entropy Scaling of the Viscosity and Thermal Conductivity of Carbon Dioxide. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120799] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
39
|
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]
|
40
|
Li N, Wang X, Chen G, Gao N. Linking Thermal Conductivity and Self-Diffusion Coefficient with a Simple Dimensionless Calorimetric Parameter for Saturated Liquids. Ind Eng Chem Res 2022. [DOI: 10.1021/acs.iecr.2c02447] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Nian Li
- Ningbo Research Institute, Zhejiang University, Ningbo315100, China
- NingboTech University, Ningbo315100, China
- Zhejiang Engineering Research Center for Intelligent Marine Ranch Equipment, Ningbo315100, China
| | - Xuehui Wang
- Department of Chemical Engineering, Imperial College London, LondonSW7 2AZ, U.K
| | - Guangming Chen
- Ningbo Research Institute, Zhejiang University, Ningbo315100, China
- NingboTech University, Ningbo315100, China
| | - Neng Gao
- NingboTech University, Ningbo315100, China
- Zhejiang Engineering Research Center for Intelligent Marine Ranch Equipment, Ningbo315100, China
| |
Collapse
|
41
|
Dehlouz A, Jaubert JN, Galliero G, Bonnissel M, Privat R. Entropy Scaling-Based Correlation for Estimating the Self-Diffusion Coefficients of Pure Fluids. Ind Eng Chem Res 2022. [DOI: 10.1021/acs.iecr.2c01086] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Aghilas Dehlouz
- Université de Lorraine, École Nationale Supérieure des Industries Chimiques, Laboratoire Réactions et Génie des Procédés (UMR CNRS 7274), 1 rue Grandville, 54000 Nancy, France
- Gaztransport & Technigaz (GTT), 1 route de Versailles, 78470 Saint-Rémy-lès-Chevreuse, France
| | - Jean-Noël Jaubert
- Université de Lorraine, École Nationale Supérieure des Industries Chimiques, Laboratoire Réactions et Génie des Procédés (UMR CNRS 7274), 1 rue Grandville, 54000 Nancy, France
| | - Guillaume Galliero
- Université de Pau et des Pays de l’Adour, E2S UPPA, CNRS TotalEnergies, LFCR UMR 5150, 64013 Pau, France
| | - Marc Bonnissel
- Gaztransport & Technigaz (GTT), 1 route de Versailles, 78470 Saint-Rémy-lès-Chevreuse, France
| | - Romain Privat
- Université de Lorraine, École Nationale Supérieure des Industries Chimiques, Laboratoire Réactions et Génie des Procédés (UMR CNRS 7274), 1 rue Grandville, 54000 Nancy, France
| |
Collapse
|
42
|
Xu X, Douglas JF, Xu WS. Thermodynamic–Dynamic Interrelations in Glass-Forming Polymer Fluids. Macromolecules 2022. [DOI: 10.1021/acs.macromol.2c01511] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Xiaolei Xu
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, P. R. China
| | - Jack F. Douglas
- Materials Science and Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Wen-Sheng Xu
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, P. R. China
| |
Collapse
|
43
|
Thermal Conductivity of Ionic Liquids: Recent Challenges Facing Theory and Experiment. J SOLUTION CHEM 2022. [DOI: 10.1007/s10953-022-01205-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/14/2022]
|
44
|
Li X, Kang K, Gu Y, Wang X. Viscosity prediction of pure refrigerants applying the residual entropy scaling theory coupled with a “Generalized Chart” parametrization method for the Statistical Associating Fluid Theory. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120479] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/14/2022]
|
45
|
Attia E, Dyre JC, Pedersen UR. Comparing four hard-sphere approximations for the low-temperature WCA melting line. J Chem Phys 2022; 157:034502. [DOI: 10.1063/5.0097593] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
By combining interface-pinning simulations with numerical integration of the Clausius–Clapeyron equation, we accurately determine the melting-line coexistence pressure and fluid/crystal densities of the Weeks–Chandler–Andersen system, covering four decades of temperature. The data are used for comparing the melting-line predictions of the Boltzmann, Andersen–Weeks–Chandler, Barker–Henderson, and Stillinger hard-sphere approximations. The Andersen–Weeks–Chandler and Barker–Henderson theories give the most accurate predictions, and they both work excellently in the zero-temperature limit for which analytical expressions are derived here.
Collapse
Affiliation(s)
- Eman Attia
- 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
| | - Ulf R. Pedersen
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
46
|
Bell I, Fingerhut R, Vrabec J, Costigliola L. Connecting Entropy Scaling and Density Scaling. J Chem Phys 2022; 157:074501. [DOI: 10.1063/5.0097088] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
It is shown that the residual entropy (entropy minus that of the ideal gas at the same temperature and density) is mostly synonymous with the independent variable of density scaling, identifying a direct link between these two approaches. The residual entropy and the effective hardness of interaction (itself a derivative at constant residual entropy) are studied for the Lennard-Jones monomer and dimer as well as a range of rigid molecular models for carbon dioxide. It is observed that the density scaling exponent appears to be related to the two-body interactions in the dilute-gas limit.
Collapse
Affiliation(s)
- Ian Bell
- National Institute of Standards and Technology Applied Chemicals and Materials Division, United States of America
| | | | - Jadran Vrabec
- Process Engineering, Technical University of Berlin, Germany
| | | |
Collapse
|
47
|
Shirai K, Watanabe K, Momida H. First-principles study of the specific heat of glass at the glass transition with a case study on glycerol. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2022; 34:375902. [PMID: 35785777 DOI: 10.1088/1361-648x/ac7e12] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/19/2022] [Accepted: 07/04/2022] [Indexed: 06/15/2023]
Abstract
The standard method to determine the transition temperature (Tg) of glasses is the jump in the specific heat,ΔCp. Despite its importance, standard theory for this jump is lacking. The difficulties include lack of proper treatment of the specific heat of liquids, hysteresis, and the timescale issue. The first part of this paper provides a non-empirical method for calculating the specific heat in the glass transition. The method consists of molecular dynamics (MD) simulations based on density-functional theory (DFT) and thermodynamics methods. Calculation of the total energy, which is the heart of DFT, is the most general method for obtaining specific heat for any state of matters. The influence of energy dissipation processes on specific heat is treated by adiabatic MD simulations. The problems of hysteresis and the timescale are alleviated by restricting the scope of calculations to equilibrium states only. The second part of this paper demonstrates the validity and usefulness of the methods by applying to the specific-heat jump of glycerol. By decomposingΔCpinto contributions of the structural, phonon, and thermal expansion energies, an appropriate interpretation for the specific-heat jump has been established: the major contribution toΔCpis the change in the structural energy. From this, a neat energy diagram about the glass transition is obtained. An outcome of this study is verification of the empirical relationship between the fragility and the specific-heat jump. These two quantities scale to the ratiok=Tg/ΔTg, whereΔTgis the width of the transition, through which the two quantities are interrelated.
Collapse
Affiliation(s)
- Koun Shirai
- The Institute of Scientific and Industrial Research (SANKEN), Osaka University, 8-1 Mihogaoka, Ibaraki, Osaka 567-0047, Japan
| | - Kota Watanabe
- The Institute of Scientific and Industrial Research (SANKEN), Osaka University, 8-1 Mihogaoka, Ibaraki, Osaka 567-0047, Japan
| | - Hiroyoshi Momida
- The Institute of Scientific and Industrial Research (SANKEN), Osaka University, 8-1 Mihogaoka, Ibaraki, Osaka 567-0047, Japan
| |
Collapse
|
48
|
Harris KR, Kanakubo M. Effect of pressure on the transport properties of 1-hexyl-3-methylimidazolium bis(trifluoromethylsulfonyl)amide and 1-hexyl-3-methylimidazolium tetrafluoroborate. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.119109] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
49
|
Abstract
Active-matter systems feature discrete particles that can convert stored or ambient free energy into motion. To realize the engineering potential of active matter, there is a strong need for predictive and theoretically grounded techniques for describing transport in these systems. In this work, we perform molecular-dynamics (MD) simulations of a model active-matter system, in which we vary the total fraction of active particles (0.01 ≤ ϕ ≤ 0.5) as well as the degree of activity of the active particles. These simulations reveal a fascinating array of transport phenomena, including activity-enhanced diffusion coefficients. By adapting an existing result for binary (inactive) fluids, we demonstrate the existence of an excess entropy scaling relation in an active system. This relationship is well supported by our MD results and establishes a new connection between transport (dynamics) and structure (statics) in active matter, a promising step for predictive and generalizable models of other transport phenomena in such systems.
Collapse
Affiliation(s)
- S Arman Ghaffarizadeh
- Department of Mechanical Engineering, Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, Pennsylvania 15213, United States
| | - Gerald J Wang
- Department of Civil and Environmental Engineering, Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, Pennsylvania 15213, United States
| |
Collapse
|
50
|
Khrapak SA, Khrapak A. Freezing density scaling of fluid transport properties: Application to liquified noble gases. J Chem Phys 2022; 157:014501. [DOI: 10.1063/5.0096947] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A freezing density scaling of transport properties of the Lennard-Jones fluid is rationalized in terms of the Rosenfeld's excess entropy scaling and isomorph theory of Roskilde-simple systems. Then, it is demonstrated that the freezing density scaling operates reasonably well for viscosity and thermal conductivity coefficients of liquid argon, krypton, and xenon. Quasi-universality of the reduced transport coefficients at their minima and at freezing conditions is discussed. The magnitude of the thermal conductivity coefficient at the freezing point is shown to agree remarkably well with the prediction of the vibrational model of thermal transport in dense fluids.
Collapse
Affiliation(s)
- Sergey A. Khrapak
- Complex Plasma, FSBSI Joint Institute for High Temperatures of the Russian Academy of Sciences, Russia
| | - Alexey Khrapak
- Theoretical Department, Joint Institute for High Temperatures RAS, Russia
| |
Collapse
|