1
|
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
|
2
|
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
|
3
|
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
|
4
|
Ranieri U, Formisano F, Gorelli FA, Santoro M, Koza MM, De Francesco A, Bove LE. Crossover from gas-like to liquid-like molecular diffusion in a simple supercritical fluid. Nat Commun 2024; 15:4142. [PMID: 38755136 PMCID: PMC11099187 DOI: 10.1038/s41467-024-47961-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2023] [Accepted: 04/15/2024] [Indexed: 05/18/2024] Open
Abstract
According to textbooks, no physical observable can be discerned allowing to distinguish a liquid from a gas beyond the critical point. Yet, several proposals have been put forward challenging this view and various transition boundaries between a gas-like and a liquid-like behaviour, including the so-called Widom and Frenkel lines, and percolation line, have been suggested to delineate the supercritical state space. Here we report observation of a crossover from gas-like (Gaussian) to liquid-like (Lorentzian) self-dynamic structure factor by incoherent quasi-elastic neutron scattering measurements on supercritical fluid methane as a function of pressure, along the 200 K isotherm. The molecular self-diffusion coefficient was derived from the best Gaussian (at low pressures) or Lorentzian (at high pressures) fits to the neutron spectra. The Gaussian-to-Lorentzian crossover is progressive and takes place at about the Widom line intercept (59 bar). At considerably higher pressures, a liquid-like jump diffusion mechanism properly describes the supercritical fluid on both sides of the Frenkel line. The present observation of a gas-like to liquid-like crossover in the self dynamics of a simple supercritical fluid confirms emerging views on the unexpectedly complex physics of the supercritical state, and could have planet-wide implications and possible industrial applications in green chemistry.
Collapse
Affiliation(s)
- Umbertoluca Ranieri
- Dipartimento di Fisica, Università di Roma La Sapienza, Piazzale Aldo Moro 5, Roma, 00187, Italy
- Centre for Science at Extreme Conditions and School of Physics and Astronomy, University of Edinburgh, Edinburgh, EH9 3FD, UK
| | - Ferdinando Formisano
- CNR - Istituto Officina dei Materiali (IOM), Grenoble, INSIDE@ILL, 71 Avenue des Martyrs, Grenoble, Cedex 9, France.
- Institut Laue-Langevin, 71 Avenue des Martyrs, Grenoble, Cedex 9, France.
| | - Federico A Gorelli
- Center for High Pressure Science and Technology Advanced Research (HPSTAR), 1690 Cailun Road, Shanghai, 201203, China.
- Shanghai Advanced Research in Physical Sciences (SHARPS), Pudong, Shanghai, 201203, China.
- Consiglio Nazionale delle Ricerche, Istituto Nazionale di Ottica, CNR-INO, Via Nello Carrara 1, Sesto Fiorentino (FI), 50019, Italy.
| | - Mario Santoro
- Consiglio Nazionale delle Ricerche, Istituto Nazionale di Ottica, CNR-INO, Via Nello Carrara 1, Sesto Fiorentino (FI), 50019, Italy
- European Laboratory for Nonlinear Spectroscopy, LENS, Via Nello Carrara 1, Sesto Fiorentino (FI), 50019, Italy
| | - Michael Marek Koza
- Institut Laue-Langevin, 71 Avenue des Martyrs, Grenoble, Cedex 9, France
| | - Alessio De Francesco
- CNR - Istituto Officina dei Materiali (IOM), Grenoble, INSIDE@ILL, 71 Avenue des Martyrs, Grenoble, Cedex 9, France
- Institut Laue-Langevin, 71 Avenue des Martyrs, Grenoble, Cedex 9, France
| | - Livia E Bove
- Dipartimento di Fisica, Università di Roma La Sapienza, Piazzale Aldo Moro 5, Roma, 00187, Italy
- Laboratory of Quantum Magnetism, Institute of Physics, École Polytechnique Fédérale de Lausanne, Lausanne, CH-1015, Switzerland
- Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Sorbonne Université, UMR CNRS 7590, 5 Place Jussieu, Paris, 75005, France
| |
Collapse
|
5
|
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
|
6
|
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
|
7
|
Kondratyuk N, Ryltsev R, Ankudinov V, Chtchelkatchev N. First-principles calculations of the viscosity in multicomponent metallic melts: Al-Cu-Ni as a test case. J Mol Liq 2023. [DOI: 10.1016/j.molliq.2023.121751] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/03/2023]
|
8
|
Jäger A, Steinberg L, Mickoleit E, Thol M. Residual Entropy Scaling for Long-Chain Linear Alkanes and Isomers of Alkanes. Ind Eng Chem Res 2023. [DOI: 10.1021/acs.iecr.2c04238] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/19/2023]
Affiliation(s)
- A. Jäger
- Institute of Power Engineering, Faculty of Mechanical Science and Engineering, Technische Universität Dresden, Helmholtzstraße 14, 01069 Dresden, Germany
| | - L. Steinberg
- Thermodynamics, Ruhr University Bochum, 44780 Bochum, Germany
| | - E. Mickoleit
- Institute of Power Engineering, Faculty of Mechanical Science and Engineering, Technische Universität Dresden, Helmholtzstraße 14, 01069 Dresden, Germany
| | - M. Thol
- Thermodynamics, Ruhr University Bochum, 44780 Bochum, Germany
| |
Collapse
|
9
|
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
|
10
|
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]
|
11
|
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
|
12
|
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
|
13
|
Herrero C, Pauletti M, Tocci G, Iannuzzi M, Joly L. Connection between water's dynamical and structural properties: Insights from ab initio simulations. Proc Natl Acad Sci U S A 2022; 119:e2121641119. [PMID: 35588447 PMCID: PMC9173753 DOI: 10.1073/pnas.2121641119] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2021] [Accepted: 04/12/2022] [Indexed: 01/25/2023] Open
Abstract
SignificanceFirst-principles calculations, which explicitly account for the electronic structure of matter, can shed light on the molecular structure and dynamics of water in its supercooled state. In this work, we use density functional theory, which relies on a functional to describe electronic exchange and correlations, to evaluate which functional best describes the temperature evolution of bulk water transport coefficients. We also assess the validity of the Stokes-Einstein relation for all the functionals in the temperature range studied, and explore the link between structure and dynamics. Based on these results, we show how transport coefficients can be computed from structural descriptors, which require shorter simulation times to converge, and we point toward strategies to develop better functionals.
Collapse
Affiliation(s)
- Cecilia Herrero
- Univ Lyon, Univ Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| | - Michela Pauletti
- Department of Chemistry, Universität Zürich, 8057 Zürich, Switzerland
| | - Gabriele Tocci
- Department of Chemistry, Universität Zürich, 8057 Zürich, Switzerland
| | - Marcella Iannuzzi
- Department of Chemistry, Universität Zürich, 8057 Zürich, Switzerland
| | - Laurent Joly
- Univ Lyon, Univ Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
- Institut Universitaire de France (IUF), 75005 Paris, France
| |
Collapse
|
14
|
Philippi F, Rauber D, Eliasen KL, Bouscharain N, Niss K, Kay CWM, Welton T. Pressing matter: why are ionic liquids so viscous? Chem Sci 2022; 13:2735-2743. [PMID: 35340854 PMCID: PMC8890108 DOI: 10.1039/d1sc06857a] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2021] [Accepted: 02/06/2022] [Indexed: 01/09/2023] Open
Abstract
Room temperature ionic liquids are considered to have huge potential for practical applications such as batteries. However, their high viscosity presents a significant challenge to their use changing from niche to ubiquitous. The modelling and prediction of viscosity in ionic liquids is the subject of an ongoing debate involving two competing hypotheses: molecular and local mechanisms versus collective and long-range mechanisms. To distinguish between these two theories, we compared an ionic liquid with its uncharged, isoelectronic, isostructural molecular mimic. We measured the viscosity of the molecular mimic at high pressure to emulate the high densities in ionic liquids, which result from the Coulomb interactions in the latter. We were thus able to reveal that the relative contributions of coulombic compaction and the charge network interactions are of similar magnitude. We therefore suggest that the optimisation of the viscosity in room temperature ionic liquids must follow a dual approach.
Collapse
Affiliation(s)
- Frederik Philippi
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London White City Campus London W12 0BZ UK
| | - Daniel Rauber
- Department of Chemistry, Saarland University Campus B2.2 Saarbrücken Germany
| | - Kira Lieberkind Eliasen
- "Glass and Time", IMFUFA, Department of Science and Environment, Roskilde University P.O. Box 260 DK-4000 Roskilde Denmark
| | | | - Kristine Niss
- "Glass and Time", IMFUFA, Department of Science and Environment, Roskilde University P.O. Box 260 DK-4000 Roskilde Denmark
| | - Christopher W M Kay
- Department of Chemistry, Saarland University Campus B2.2 Saarbrücken Germany.,London Centre for Nanotechnology, University College London 17-19 Gordon Street London WC1H 0AH UK
| | - Tom Welton
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London White City Campus London W12 0BZ UK
| |
Collapse
|
15
|
Mitra S, Marín-Aguilar S, Sastry S, Smallenburg F, Foffi G. Correlation between plastic rearrangements and local structure in a cyclically driven glass. J Chem Phys 2022; 156:074503. [DOI: 10.1063/5.0077851] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | | | - Srikanth Sastry
- Theoretical Sciences Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, India
| | | | - Giuseppe Foffi
- Laboratoire de Physique des Solides, Laboratoire de Physique des Solides, France
| |
Collapse
|
16
|
Li N, Wang XH, Gao N, Chen GM. Simple Direct Relationship between Scaled Viscosity and a Dimensionless Calorimetric Parameter for Saturated Liquids. Ind Eng Chem Res 2022. [DOI: 10.1021/acs.iecr.1c03508] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- N. Li
- Ningbo Research Institute, Zhejiang University, Ningbo 315100, China
- NingboTech University, Ningbo 315100, China
| | - X. H. Wang
- Fluids and Thermal Engineering Research Group, Faculty of Engineering, University of Nottingham, Nottingham NG7 2RD, United Kingdom
| | - N. Gao
- NingboTech University, Ningbo 315100, China
| | - G. M. Chen
- Ningbo Research Institute, Zhejiang University, Ningbo 315100, China
- NingboTech University, Ningbo 315100, China
| |
Collapse
|
17
|
Nandi UK, Patel P, Moid M, Nandi MK, Sengupta S, Karmakar S, Maiti PK, Dasgupta C, Maitra Bhattacharyya S. Thermodynamics and its correlation with dynamics in a mean-field model and pinned systems: A comparative study using two different methods of entropy calculation. J Chem Phys 2022; 156:014503. [PMID: 34998317 DOI: 10.1063/5.0065668] [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 recent study introduced a novel mean-field model system where each particle over and above the interaction with its regular neighbors interacts with k extra pseudo-neighbors. Here, we present an extensive study of thermodynamics and its correlation with the dynamics of this system. We surprisingly find that the well-known thermodynamic integration (TI) method of calculating the entropy provides unphysical results. It predicts vanishing of the configurational entropy at temperatures close to the onset temperature of the system and negative values of the configurational entropy at lower temperatures. Interestingly, well below the temperature at which the configurational entropy vanishes, both the collective and the single-particle dynamics of the system show complete relaxation. Negative values of the configurational entropy are unphysical, and complete relaxation when the configurational entropy is zero violates the prediction of the random first-order transition theory (RFOT). However, the entropy calculated using the two-phase thermodynamics (2PT) method remains positive at all temperatures for which we can equilibrate the system, and its values are consistent with RFOT predictions. We find that with an increase in k, the difference in the entropy computed using the two methods increases. A similar effect is also observed for a system where a randomly selected fraction of the particles are pinned in their positions in the equilibrated liquid. We show that the difference in entropy calculated via the 2PT and TI methods increases with pinning density.
Collapse
Affiliation(s)
- Ujjwal Kumar Nandi
- Polymer Science and Engineering Division, CSIR-National Chemical Laboratory, Pune 411008, India
| | - Palak Patel
- Polymer Science and Engineering Division, CSIR-National Chemical Laboratory, Pune 411008, India
| | - Mohd Moid
- Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India
| | - Manoj Kumar Nandi
- Department of Engineering, University of Campania "Luigi Vanvitelli", 81031 Aversa (CE), Italy
| | - Shiladitya Sengupta
- Department of Physics, Indian Institute of Technology, Roorkee 247667, India
| | - Smarajit Karmakar
- Centre for Interdisciplinary Sciences, Tata Institute of Fundamental Research, 36/P, Gopanpally Village, Serilingampally Mandal, RR District, Hyderabad 500019, India
| | - Prabal K Maiti
- Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India
| | - Chandan Dasgupta
- Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India
| | | |
Collapse
|
18
|
Knudsen S, Todd BD, Dyre JC, Hansen JS. Generalized hydrodynamics of the Lennard-Jones liquid in view of hidden scale invariance. Phys Rev E 2021; 104:054126. [PMID: 34942805 DOI: 10.1103/physreve.104.054126] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2021] [Accepted: 10/12/2021] [Indexed: 11/07/2022]
Abstract
In recent years lines along which structure and dynamics are invariant to a good approximation, so-called isomorphs, have been identified in the thermodynamic phase diagrams of several model liquids and solids. This paper reports computer simulation data of the transverse and longitudinal collective dynamics at different length scales along an isomorph of the Lennard-Jones system. Our findings are compared to corresponding results along an isotherm and an isochore. Confirming the theoretical prediction, the reduced-unit dynamics of the transverse momentum density is invariant to a good approximation along the isomorph on all time and length scales. Likewise, the wave-vector dependent shear-stress autocorrelation function is found to be isomorph invariant (with minor deviations at very short times). A similar invariance is not seen along the isotherm or the isochore. Using a spatially nonlocal hydrodynamic model for the transverse momentum-density time-autocorrelation function, the macroscopic shear viscosity and its wave dependence are determined, demonstrating that the shear viscosity is isomorphic invariant on all length scales studied. This analysis implies the existence of a length scale that is isomorph invariant in reduced units, i.e., which characterizes each isomorph. The transverse sound-wave velocity, the Maxwell relaxation time, and the rigidity shear modulus are also isomorph invariant. In contrast to the isomorph invariance of all aspects of the transverse dynamics, the reduced-unit dynamics of the mass density is not invariant on length scales longer than the interparticle distance. By fitting to a generalized hydrodynamic model, we extract values for the wave-vector-dependent thermal diffusion coefficient, sound attenuation coefficient, and adiabatic sound velocity. The isomorph variation of these quantities in reduced units on long length scales can be eliminated by scaling with the density-scaling exponent, a fundamental quantity in the isomorph theory framework; this is an empirical observation that remains to be explained theoretically.
Collapse
Affiliation(s)
- Solvej Knudsen
- "Glass and Time", IMFUFA, Department of Science and Environment, Roskilde University, Postbox 260, DK-4000 Roskilde, Denmark.,Department of Mathematics, School of Science, Computing and Engineering Technologies, Swinburne University of Technology, Hawtorn, Victoria 3122, Australia
| | - B D Todd
- Department of Mathematics, School of Science, Computing and Engineering Technologies, Swinburne University of Technology, Hawtorn, Victoria 3122, Australia
| | - Jeppe C Dyre
- "Glass and Time", IMFUFA, Department of Science and Environment, Roskilde University, Postbox 260, DK-4000 Roskilde, Denmark
| | - J S Hansen
- "Glass and Time", IMFUFA, Department of Science and Environment, Roskilde University, Postbox 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
19
|
Saliou A, Jarry P, Jakse N. Excess entropy scaling law: A potential energy landscape view. Phys Rev E 2021; 104:044128. [PMID: 34781503 DOI: 10.1103/physreve.104.044128] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2021] [Accepted: 09/28/2021] [Indexed: 11/07/2022]
Abstract
The relationship between excess entropy and diffusion is revisited by means of large-scale computer simulation combined to supervised learning approach to determine the excess entropy for the Lennard-Jones potential. Results reveal a strong correlation with the properties of the potential energy landscape (PEL). In particular the exponential law holding in the liquid is seen to be linked with the landscape-influenced regime of the PEL whereas the fluidlike power-law corresponds to the free diffusion regime.
Collapse
Affiliation(s)
- Anthony Saliou
- Université Grenoble Alpes, CNRS, Grenoble INP, SIMaP, F-38000 Grenoble, France
| | - Philippe Jarry
- C-TEC, Parc Economique Centr'alp, 725 rue Aristide Bergès, CS10027, Voreppe 38341 cedex, France
| | - Noel Jakse
- Université Grenoble Alpes, CNRS, Grenoble INP, SIMaP, F-38000 Grenoble, France
| |
Collapse
|
20
|
Bell IH, Delage-Santacreu S, Hoang H, Galliero G. Dynamic Crossover in Fluids: From Hard Spheres to Molecules. J Phys Chem Lett 2021; 12:6411-6417. [PMID: 34232673 DOI: 10.1021/acs.jpclett.1c01594] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
We propose a simple and generic definition of a demarcation reconciling structural and dynamic frameworks when combined with the entropy scaling framework. This crossover line between gas- and liquid-like behaviors is defined as the curve for which an individual property, the contribution to viscosity due to molecules' translation, is exactly equal to a collective property, the contribution to viscosity due to molecular interactions. Such a definition is shown to be consistent with the one based on the minima of the kinematic viscosity. For the hard sphere, this is shown to be an exact solution. For Lennard-Jones spheres and dimers and for some simple real fluids, this relation holds very well. This crossover line passes nearby the critical point, and for all studied fluids, it is well captured by the critical excess entropy curve for atomic fluids, emphasizing the link between transport properties and local structure.
Collapse
Affiliation(s)
- Ian H Bell
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, United States
| | - Stéphanie Delage-Santacreu
- Université de Pau et des Pays de l'Adour, e2s UPPA, Laboratoire de Mathematiques et de leurs Applications de Pau (IPRA, CNRS UMR5142), Pau 64000, France
| | - Hai Hoang
- Institute of Fundamental and Applied Sciences, Duy Tan University, 10C Tran Nhat Duat Street, District 1, Ho Chi Minh City 700000, Vietnam
- Faculty of Natural Sciences, Duy Tan University, Da Nang 550000, Vietnam
| | - Guillaume Galliero
- Université de Pau et des Pays de l'Adour, e2s UPPA, TOTAL, CNRS, LFCR, UMR 5150, Laboratoire des fluides complexes et leurs reservoirs, Pau 64000, France
| |
Collapse
|
21
|
Hultmark S, Cravcenco A, Kushwaha K, Mallick S, Erhart P, Börjesson K, Müller C. Vitrification of octonary perylene mixtures with ultralow fragility. SCIENCE ADVANCES 2021; 7:7/29/eabi4659. [PMID: 34272241 PMCID: PMC8284888 DOI: 10.1126/sciadv.abi4659] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/11/2021] [Accepted: 06/02/2021] [Indexed: 05/08/2023]
Abstract
Strong glass formers with a low fragility are highly sought-after because of the technological importance of vitrification. In the case of organic molecules and polymers, the lowest fragility values have been reported for single-component materials. Here, we establish that mixing of organic molecules can result in a marked reduction in fragility. Individual bay-substituted perylene derivatives display a high fragility of more than 70. Instead, slowly cooled perylene mixtures with more than three components undergo a liquid-liquid transition and turn into a strong glass former. Octonary perylene mixtures display a fragility of 13 ± 2, which not only is a record low value for organic molecules but also lies below values reported for the strongest known inorganic glass formers. Our work opens an avenue for the design of ultrastrong organic glass formers, which can be anticipated to find use in pharmaceutical science and organic electronics.
Collapse
Affiliation(s)
- Sandra Hultmark
- Department of Chemistry and Chemical Engineering, Chalmers University of Technology, 41296 Göteborg, Sweden
| | - Alex Cravcenco
- Department of Chemistry and Molecular Biology, University of Gothenburg, Kemigården 4, 41296 Göteborg, Sweden
| | - Khushbu Kushwaha
- Department of Chemistry and Molecular Biology, University of Gothenburg, Kemigården 4, 41296 Göteborg, Sweden
| | - Suman Mallick
- Department of Chemistry and Molecular Biology, University of Gothenburg, Kemigården 4, 41296 Göteborg, Sweden
| | - Paul Erhart
- Department of Physics, Chalmers University of Technology, 41296 Göteborg, Sweden
| | - Karl Börjesson
- Department of Chemistry and Molecular Biology, University of Gothenburg, Kemigården 4, 41296 Göteborg, Sweden
| | - Christian Müller
- Department of Chemistry and Chemical Engineering, Chalmers University of Technology, 41296 Göteborg, Sweden.
| |
Collapse
|
22
|
Zhang M, Chen Y, Dai LH. Universal Scaling in the Temperature-Dependent Viscous Dynamics of Metallic Glasses. J Phys Chem B 2021; 125:3419-3425. [PMID: 33764771 DOI: 10.1021/acs.jpcb.1c00034] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The essential query about glass formation is how to understand the sheer temperature dependence of viscous dynamics of glass-forming liquids near the liquid-to-glass-transition temperature Tg. In this work, we report a universal scaling in the temperature-dependent viscous dynamics of metallic glasses (MGs) in the form of the Williams-Landel-Ferry equation on the basis of compiled data on the temperature-dependent viscosity and structural relaxation times of 89 MGs ever-reported in the past decades. Implications of this universal scaling are illustrated in the framework of the Adam-Gibbs relation, suggesting a universal vitrification mechanism in MGs mediated by configurational entropy wherein configurational entropy vanishes universally for all supercooled metallic liquids after a further decrease in temperature of ∼170.7 K (whereas with a relatively large error of ±150 K) below Tg. This result corroborates the thermodynamic origin of glass formation and suggests that MGs are an ideal research subject for understanding in depth the nature of glass transition for their relatively simple molecular structures.
Collapse
Affiliation(s)
- Meng Zhang
- Institute of Advanced Wear & Corrosion Resistant and Functional Materials, Jinan University, Guangzhou 510632, China.,State Key Laboratory of Nonlinear Mechanics, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China
| | - Yan Chen
- State Key Laboratory of Nonlinear Mechanics, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China.,School of Engineering Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Lan-Hong Dai
- State Key Laboratory of Nonlinear Mechanics, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China.,School of Engineering Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
23
|
Scheiner B, Yoon TJ. Calculation of self-diffusion coefficients in supercritical carbon dioxide using mean force kinetic theory. J Chem Phys 2021; 154:134101. [PMID: 33832259 DOI: 10.1063/5.0045211] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This paper presents an application of mean force kinetic theory (MFT) to the calculation of the self-diffusivity of CO2 in the supercritical fluid regime. Two modifications to the typical application of MFT are employed to allow its application to a system of molecular species. The first is the assumption that the inter-particle potential of mean force can be obtained from the molecule center-of-mass pair correlation function, which in the case of CO2 is the C-C pair correlation function. The second is a new definition of the Enskog factor that describes the effect of correlations at the surface of the collision volume. The new definition retains the physical picture that this quantity represents a local density increase, resulting from particle correlations, relative to that in the zero density homogeneous fluid limit. These calculations are facilitated by the calculation of pair correlation functions from molecular dynamics (MD) simulations using the FEPM2 molecular CO2 model. The self-diffusivity calculated from theory is in good agreement with that from MD simulations up to and slightly beyond the density at the location of the Frenkel line. The calculation is compared with and is found to perform similarly well to other commonly used models but has a greater potential for application to systems of mixed species and to systems of particles with long range interatomic potentials due to electrostatic interactions.
Collapse
Affiliation(s)
- Brett Scheiner
- Los Alamos National Laboratory, Los Alamos, New Mexico 87544, USA
| | - Tae Jun Yoon
- Los Alamos National Laboratory, Los Alamos, New Mexico 87544, USA
| |
Collapse
|
24
|
Singh AN, Dyre JC, Pedersen UR. Solid–liquid coexistence of neon, argon, krypton, and xenon studied by simulations. J Chem Phys 2021; 154:134501. [DOI: 10.1063/5.0045398] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Affiliation(s)
- Aditya N. Singh
- Department of Chemistry, University of California, Berkeley, California 94720, USA and Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53703, USA
| | - 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
|
25
|
Diffusion in dense supercritical methane from quasi-elastic neutron scattering measurements. Nat Commun 2021; 12:1958. [PMID: 33785748 PMCID: PMC8009954 DOI: 10.1038/s41467-021-22182-4] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2020] [Accepted: 02/26/2021] [Indexed: 11/08/2022] Open
Abstract
Methane, the principal component of natural gas, is an important energy source and raw material for chemical reactions. It also plays a significant role in planetary physics, being one of the major constituents of giant planets. Here, we report measurements of the molecular self-diffusion coefficient of dense supercritical CH4 reaching the freezing pressure. We find that the high-pressure behaviour of the self-diffusion coefficient measured by quasi-elastic neutron scattering at 300 K departs from that expected for a dense fluid of hard spheres and suggests a density-dependent molecular diameter. Breakdown of the Stokes-Einstein-Sutherland relation is observed and the experimental results suggest the existence of another scaling between self-diffusion coefficient D and shear viscosity η, in such a way that Dη/ρ=constant at constant temperature, with ρ the density. These findings underpin the lack of a simple model for dense fluids including the pressure dependence of their transport properties.
Collapse
|
26
|
Mehri S, Ingebrigtsen TS, Dyre JC. Single-parameter aging in a binary Lennard-Jones system. J Chem Phys 2021; 154:094504. [PMID: 33685153 DOI: 10.1063/5.0039250] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This paper studies physical aging by computer simulations of a 2:1 Kob-Andersen binary Lennard-Jones mixture, a system that is less prone to crystallization than the standard 4:1 composition. Starting from thermal-equilibrium states, the time evolution of the following four quantities is monitored by following up and down jumps in temperature: potential energy, virial, average squared force, and the Laplacian of the potential energy. Despite the fact that significantly larger temperature jumps are studied here than in typical similar experiments, to a good approximation, all four quantities conform to the single-parameter-aging scenario derived and validated for small jumps in experiments [T. Hecksher, N. B. Olsen, and J. C. Dyre, J. Chem. Phys. 142, 241103 (2015)]. As a further confirmation of single-parameter aging with a common material time for the four different quantities monitored, their relaxing parts are found to be almost identical for all temperature jumps.
Collapse
Affiliation(s)
- Saeed Mehri
- 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
| | - Jeppe C Dyre
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
27
|
Lucco Castello F, Tolias P. Theoretical Estimate of the Glass Transition Line of Yukawa One-Component Plasmas. Molecules 2021; 26:molecules26030669. [PMID: 33525346 PMCID: PMC7865523 DOI: 10.3390/molecules26030669] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2020] [Revised: 01/07/2021] [Accepted: 01/25/2021] [Indexed: 12/04/2022] Open
Abstract
The mode coupling theory of supercooled liquids is combined with advanced closures to the integral equation theory of liquids in order to estimate the glass transition line of Yukawa one-component plasmas from the unscreened Coulomb limit up to the strong screening regime. The present predictions constitute a major improvement over the current literature predictions. The calculations confirm the validity of an existing analytical parameterization of the glass transition line. It is verified that the glass transition line is an approximate isomorphic curve and the value of the corresponding reduced excess entropy is estimated. Capitalizing on the isomorphic nature of the glass transition line, two structural vitrification indicators are identified that allow a rough estimate of the glass transition point only through simple curve metrics of the static properties of supercooled liquids. The vitrification indicators are demonstrated to be quasi-universal by an investigation of hard sphere and inverse power law supercooled liquids. The straightforward extension of the present results to bi-Yukawa systems is also discussed.
Collapse
|
28
|
Ingebrigtsen TS, Schrøder TB, Dyre JC. Hidden Scale Invariance in Polydisperse Mixtures of Exponential Repulsive Particles. J Phys Chem B 2021; 125:317-327. [PMID: 33369412 DOI: 10.1021/acs.jpcb.0c09726] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Polydisperse systems of particles interacting by the purely repulsive exponential (EXP) pair potential are studied in regard to how structure and dynamics vary along isotherms, isochores, and isomorphs. The sizable size polydispersities of 23%, 29%, 35%, and 40%, as well as energy polydispersity 35%, were considered. For each system an isomorph was traced out covering about one decade in density. For all systems studied, the structure and dynamics vary significantly along the isotherms and isochores but are invariant to a good approximation along the isomorphs. We conclude that the single-component EXP system's hidden scale invariance (implying isomorph invariance of structure and dynamics) is maintained even when a sizable polydispersity is introduced into the system.
Collapse
Affiliation(s)
- Trond S Ingebrigtsen
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Thomas B Schrøder
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Jeppe C Dyre
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|