1
|
Grzybowski A, Koperwas K, Paluch M. Role of anisotropy in understanding the molecular grounds for density scaling in dynamics of glass-forming liquids. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2024; 87:084501. [PMID: 38861964 DOI: 10.1088/1361-6633/ad569d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/14/2024] [Accepted: 06/11/2024] [Indexed: 06/13/2024]
Abstract
Molecular Dynamics (MD) simulations of glass-forming liquids play a pivotal role in uncovering the molecular nature of the liquid vitrification process. In particular, much focus was given to elucidating the interplay between the character of intermolecular potential and molecular dynamics behaviour. This has been tried to achieve by simulating the spherical particles interacting via isotropic potential. However, when simulation and experimental data are analysed in the same way by using the density scaling approaches, serious inconsistency is revealed between them. Similar scaling exponent values are determined by analysing the relaxation times and pVT data obtained from computer simulations. In contrast, these values differ significantly when the same analysis is carried out in the case of experimental data. As discussed thoroughly herein, the coherence between results of simulation and experiment can be achieved if anisotropy of intermolecular interactions is introduced to MD simulations. In practice, it has been realized in two different ways: (1) by using the anisotropic potential of the Gay-Berne type or (2) by replacing the spherical particles with quasi-real polyatomic anisotropic molecules interacting through isotropic Lenard-Jones potential. In particular, the last strategy has the potential to be used to explore the relationship between molecular architecture and molecular dynamics behaviour. Finally, we hope that the results presented in this review will also encourage others to explore how 'anisotropy' affects remaining aspects related to liquid-glass transition, like heterogeneity, glass transition temperature, glass forming ability, etc.
Collapse
Affiliation(s)
- A Grzybowski
- Institute of Physics, University of Silesia in Katowice, 41-500 Chorzów, Poland
| | - K Koperwas
- Institute of Physics, University of Silesia in Katowice, 41-500 Chorzów, Poland
| | - M Paluch
- Institute of Physics, University of Silesia in Katowice, 41-500 Chorzów, Poland
| |
Collapse
|
2
|
Jin J, Lee EK, Voth GA. Understanding dynamics in coarse-grained models. III. Roles of rotational motion and translation-rotation coupling in coarse-grained dynamics. J Chem Phys 2023; 159:164102. [PMID: 37870140 DOI: 10.1063/5.0167158] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2023] [Accepted: 09/29/2023] [Indexed: 10/24/2023] Open
Abstract
This paper series aims to establish a complete correspondence between fine-grained (FG) and coarse-grained (CG) dynamics by way of excess entropy scaling (introduced in Paper I). While Paper II successfully captured translational motions in CG systems using a hard sphere mapping, the absence of rotational motions in single-site CG models introduces differences between FG and CG dynamics. In this third paper, our objective is to faithfully recover atomistic diffusion coefficients from CG dynamics by incorporating rotational dynamics. By extracting FG rotational diffusion, we unravel, for the first time reported to our knowledge, a universality in excess entropy scaling between the rotational and translational diffusion. Once the missing rotational dynamics are integrated into the CG translational dynamics, an effective translation-rotation coupling becomes essential. We propose two different approaches for estimating this coupling parameter: the rough hard sphere theory with acentric factor (temperature-independent) or the rough Lennard-Jones model with CG attractions (temperature-dependent). Altogether, we demonstrate that FG diffusion coefficients can be recovered from CG diffusion coefficients by (1) incorporating "entropy-free" rotational diffusion with translation-rotation coupling and (2) recapturing the missing entropy. Our findings shed light on the fundamental relationship between FG and CG dynamics in molecular fluids.
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
| | - Eok Kyun Lee
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, South Korea
| | - 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
|
3
|
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
|
4
|
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]
|
5
|
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
|
6
|
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
|
7
|
Dehlouz A, Privat R, Galliero G, Bonnissel M, Jaubert JN. Revisiting the Entropy-Scaling Concept for Shear-Viscosity Estimation from Cubic and SAFT Equations of State: Application to Pure Fluids in Gas, Liquid and Supercritical States. Ind Eng Chem Res 2021. [DOI: 10.1021/acs.iecr.1c01386] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Aghilas Dehlouz
- École Nationale Supérieure des Industries Chimiques, Laboratoire Réactions et Génie des Procédés (UMR CNRS 7274), Université de Lorraine, 1 rue Grandville, 54000 Nancy, France
- Gaztransport & Technigaz (GTT), 1 route de Versailles, 78470 Saint-Rémy-lès-Chevreuse, France
| | - Romain Privat
- École Nationale Supérieure des Industries Chimiques, Laboratoire Réactions et Génie des Procédés (UMR CNRS 7274), Université de Lorraine, 1 rue Grandville, 54000 Nancy, France
| | - Guillaume Galliero
- E2S UPPA, CNRS Total Energies, LFCR UMR 5150, Université de Pau et des Pays de l’Adour 64000 Pau, France
| | - Marc Bonnissel
- Gaztransport & Technigaz (GTT), 1 route de Versailles, 78470 Saint-Rémy-lès-Chevreuse, France
| | - Jean-Noël Jaubert
- École Nationale Supérieure des Industries Chimiques, Laboratoire Réactions et Génie des Procédés (UMR CNRS 7274), Université de Lorraine, 1 rue Grandville, 54000 Nancy, France
| |
Collapse
|
8
|
Xu WS, Douglas JF, Sun ZY. Polymer Glass Formation: Role of Activation Free Energy, Configurational Entropy, and Collective Motion. Macromolecules 2021. [DOI: 10.1021/acs.macromol.0c02740] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- 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
| | - Jack F. Douglas
- Materials Science and Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Zhao-Yan Sun
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, P. R. China
| |
Collapse
|
9
|
Bell IH, Dyre JC, Ingebrigtsen TS. Excess-entropy scaling in supercooled binary mixtures. Nat Commun 2020; 11:4300. [PMID: 32855393 PMCID: PMC7453028 DOI: 10.1038/s41467-020-17948-1] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2020] [Accepted: 07/22/2020] [Indexed: 11/09/2022] Open
Abstract
Transport coefficients, such as viscosity or diffusion coefficient, show significant dependence on density or temperature near the glass transition. Although several theories have been proposed for explaining this dynamical slowdown, the origin remains to date elusive. We apply here an excess-entropy scaling strategy using molecular dynamics computer simulations and find a quasiuniversal, almost composition-independent, relation for binary mixtures, extending eight orders of magnitude in viscosity or diffusion coefficient. Metallic alloys are also well captured by this relation. The excess-entropy scaling predicts a quasiuniversal breakdown of the Stokes-Einstein relation between viscosity and diffusion coefficient in the supercooled regime. Additionally, we find evidence that quasiuniversality extends beyond binary mixtures, and that the origin is difficult to explain using existing arguments for single-component quasiuniversality.
Collapse
Affiliation(s)
- Ian H Bell
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, CO, 80305, USA
| | - Jeppe C Dyre
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, Postbox 260, Roskilde, DK-4000, Denmark
| | - Trond S Ingebrigtsen
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, Postbox 260, Roskilde, DK-4000, Denmark.
| |
Collapse
|
10
|
Binti Mohd Taib M, Trusler JPM. Residual entropy model for predicting the viscosities of dense fluid mixtures. J Chem Phys 2020; 152:164104. [PMID: 32357798 DOI: 10.1063/5.0002242] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
In this work, we have investigated the mono-variant relationship between the reduced viscosity and residual entropy in pure fluids and in binary mixtures of hydrocarbons and hydrocarbons with dissolved carbon dioxide. The mixtures considered were octane + dodecane, decane + carbon dioxide, and 1,3-dimethylbenzene (m-xylene) + carbon dioxide. The reduced viscosity was calculated according to the definition of Bell, while the residual entropy was calculated from accurate multi-parameter Helmholtz-energy equations of state and, for mixtures, the multi-fluid Helmholtz energy approximation. The mono-variant dependence of reduced viscosity upon residual molar entropy was observed for the pure fluids investigated, and by incorporating two scaling factors (one for reduced viscosity and the other for residual molar entropy), the data were represented by a single universal curve. To apply this method to mixtures, the scaling factors were determined from a mole-fraction weighted sum of the pure-component values. This simple model was found to work well for the systems investigated. The average absolute relative deviation (AARD) was observed to be between 1% and 2% for pure components and a mixture of similar hydrocarbons. Larger deviations, with AARDs of up to 15%, were observed for the asymmetric mixtures, but this compares favorably with other methods for predicting the viscosity of such systems. We conclude that the residual-entropy concept can be used to estimate the viscosity of mixtures of similar molecules with high reliability and that it offers a useful engineering approximation even for asymmetric mixtures.
Collapse
Affiliation(s)
- Malyanah Binti Mohd Taib
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - J P Martin Trusler
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| |
Collapse
|
11
|
Yungbluth JC, Medvedev GA, Savoie BM, Caruthers JM. Temperature and pressure dependence of the alpha relaxation in ortho-terphenyl. J Chem Phys 2020; 152:094504. [PMID: 33480716 DOI: 10.1063/1.5144283] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Molecular dynamics (MD) simulations of ortho-terphenyl using an all-atom model with the optimized potentials for liquid simulations (OPLS) force field were performed both in the high temperature Arrhenian region and at lower temperatures that include the onset of the super-Arrhenian region. From the MD simulations, the internal energy of both the equilibrium liquid and crystal was determined from 300 K to 600 K and at pressures from 0.1 MPa to 1 GPa. The translational and rotational diffusivities were also determined at these temperatures and pressures for the equilibrium liquid. It is shown that within a small offset, the excess internal energy Ūx from the MD simulations is consistent with the experimentally determined excess internal energy reported earlier [Caruthers and Medvedev, Phys. Rev. Mater. 2, 055604, (2018)]. The MD mobility data {including extremely long-time 1 atm simulations from the study by Eastwood et al. [J. Phys. Chem. B 117, 12898, (2013)]} were combined with experimental data to form a unified dataset, where it was shown that in both the high temperature Arrhenian region and the lower temperature super-Arrhenian region, the mobility is a linear function of 1/Ūx(T,p), albeit with different proportionality constants. The transition between the Arrhenian and super-Arrhenian regions is relatively sharp at a critical internal energy Ūx α. The 1/Ūx(T,p) model is able to describe the mobility data over nearly 16 orders-of-magnitude. Other excess thermodynamic properties such as excess enthalpy and excess entropy (i.e., the Adam-Gibbs model) are unable to unify the pressure dependence of the mobility.
Collapse
Affiliation(s)
- Jack C Yungbluth
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, Indiana 47907, USA
| | - Grigori A Medvedev
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, Indiana 47907, USA
| | - Brett M Savoie
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, Indiana 47907, USA
| | - James M Caruthers
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, Indiana 47907, USA
| |
Collapse
|
12
|
Hopp M, Mele J, Hellmann R, Gross J. Thermal Conductivity via Entropy Scaling: An Approach That Captures the Effect of Intramolecular Degrees of Freedom. Ind Eng Chem Res 2019. [DOI: 10.1021/acs.iecr.9b03998] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Madlen Hopp
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, 70569 Stuttgart, Germany
| | - Julia Mele
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, 70569 Stuttgart, Germany
| | - Robert Hellmann
- Institute of Chemistry, University of Rostock, 18055 Rostock, Germany
| | - Joachim Gross
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, 70569 Stuttgart, Germany
| |
Collapse
|
13
|
Dahanayake JN, Shahryari E, Roberts KM, Heikes ME, Kasireddy C, Mitchell-Koch KR. Protein Solvent Shell Structure Provides Rapid Analysis of Hydration Dynamics. J Chem Inf Model 2019; 59:2407-2422. [PMID: 30865440 DOI: 10.1021/acs.jcim.9b00009] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
The solvation layer surrounding a protein is clearly an intrinsic part of protein structure-dynamics-function, and our understanding of how the hydration dynamics influences protein function is emerging. We have recently reported simulations indicating a correlation between regional hydration dynamics and the structure of the solvation layer around different regions of the enzyme Candida antarctica lipase B, wherein the radial distribution function (RDF) was used to calculate the pairwise entropy, providing a link between dynamics (diffusion) and thermodynamics (excess entropy) known as Rosenfeld scaling. Regions with higher RDF values/peaks in the hydration layer (the first peak, within 6 Å of the protein surface) have faster diffusion in the hydration layer. The finding thus hinted at a handle for rapid evaluation of hydration dynamics at different regions on the protein surface in molecular dynamics simulations. Such an approach may move the analysis of hydration dynamics from a specialized venture to routine analysis, enabling an informatics approach to evaluate the role of hydration dynamics in biomolecular function. This paper first confirms that the correlation between regional diffusive dynamics and hydration layer structure (via water center of mass around protein side-chain atom RDF) is observed as a general relationship across a set of proteins. Second, it seeks to devise an approach for rapid analysis of hydration dynamics, determining the minimum amount of information and computational effort required to get a reliable value of hydration dynamics from structural data in MD simulations based on the protein-water RDF. A linear regression model using the integral of the hydration layer in the water-protein RDF was found to provide statistically equivalent apparent diffusion coefficients at the 95% confidence level for a set of 92 regions within five different proteins. In summary, RDF analysis of 10 ns of data after simulation convergence is sufficient to accurately map regions of fast and slow hydration dynamics around a protein surface. Additionally, it is anticipated that a quick look at protein-water RDFs, comparing peak heights, will be useful to provide a qualitative ranking of regions of faster and slower hydration dynamics at the protein surface for rapid analysis when investigating the role of solvent dynamics in protein function.
Collapse
Affiliation(s)
- Jayangika N Dahanayake
- Department of Chemistry , Wichita State University , 1845 Fairmount Street , Wichita , Kansas 67260-0051 , United States
| | - Elaheh Shahryari
- Department of Chemistry , Wichita State University , 1845 Fairmount Street , Wichita , Kansas 67260-0051 , United States
| | - Kirsten M Roberts
- Department of Chemistry , Wichita State University , 1845 Fairmount Street , Wichita , Kansas 67260-0051 , United States
| | - Micah E Heikes
- Department of Chemistry , Wichita State University , 1845 Fairmount Street , Wichita , Kansas 67260-0051 , United States
| | - Chandana Kasireddy
- Department of Chemistry , Wichita State University , 1845 Fairmount Street , Wichita , Kansas 67260-0051 , United States
| | - Katie R Mitchell-Koch
- Department of Chemistry , Wichita State University , 1845 Fairmount Street , Wichita , Kansas 67260-0051 , United States
| |
Collapse
|
14
|
Abstract
This article gives an overview of excess-entropy scaling, the 1977 discovery by Rosenfeld that entropy determines properties of liquids like viscosity, diffusion constant, and heat conductivity. We give examples from computer simulations confirming this intriguing connection between dynamics and thermodynamics, counterexamples, and experimental validations. Recent uses in application-related contexts are reviewed, and theories proposed for the origin of excess-entropy scaling are briefly summarized. It is shown that if two thermodynamic state points of a liquid have the same microscopic dynamics, they must have the same excess entropy. In this case, the potential-energy function exhibits a symmetry termed hidden scale invariance, stating that the ordering of the potential energies of configurations is maintained if these are scaled uniformly to a different density. This property leads to the isomorph theory, which provides a general framework for excess-entropy scaling and illuminates, in particular, why this does not apply rigorously and universally. It remains an open question whether all aspects of excess-entropy scaling and related regularities reflect hidden scale invariance in one form or other.
Collapse
Affiliation(s)
- Jeppe C Dyre
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
15
|
Witman M, Mahynski NA, Smit B. Flat-Histogram Monte Carlo as an Efficient Tool To Evaluate Adsorption Processes Involving Rigid and Deformable Molecules. J Chem Theory Comput 2018; 14:6149-6158. [PMID: 30296088 DOI: 10.1021/acs.jctc.8b00534] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Monte Carlo simulations are the foundational technique for predicting thermodynamic properties of open systems where the process of interest involves the exchange of particles. Thus, they have been used extensively to computationally evaluate the adsorption properties of nanoporous materials and are critical for the in silico identification of promising materials for a variety of gas storage and chemical separation applications. In this work we demonstrate that a well-known biasing technique, known as "flat-histogram" sampling, can be combined with temperature extrapolation of the free energy landscape to efficiently provide significantly more useful thermodynamic information than standard open ensemble MC simulations. Namely, we can accurately compute the isosteric heat of adsorption and number of particles adsorbed for various adsorbates over an extremely wide range of temperatures and pressures from a set of simulations at just one temperature. We extend this derivation of the temperature extrapolation to adsorbates with intramolecular degrees of freedom when Rosenbluth sampling is employed. Consequently, the working capacity and isosteric heat can be computed for any given combined temperature/pressure swing adsorption process for a large range of operating conditions with both rigid and deformable adsorbates. Continuous thermodynamic properties can be computed with this technique at very moderate computational cost, thereby providing a strong case for its application to the in silico identification of promising nanoporous adsorbents.
Collapse
Affiliation(s)
- Matthew Witman
- Department of Chemical and Biomolecular Engineering , University of California , Berkeley 94720 , United States.,Laboratory of Molecular Simulation (LSMO) , Institut des Sciences et Ingénierie Chimiques, Valais, École Polytechnique Fédérale de Lausanne (EPFL) , Rue de l'Industrie 17 , CH-1951 Sion , Switzerland
| | - Nathan A Mahynski
- Chemical Sciences Division , National Institute of Standards and Technology , Gaithersburg , Maryland 20899-8320 , United States
| | - Berend Smit
- Department of Chemical and Biomolecular Engineering , University of California , Berkeley 94720 , United States.,Laboratory of Molecular Simulation (LSMO) , Institut des Sciences et Ingénierie Chimiques, Valais, École Polytechnique Fédérale de Lausanne (EPFL) , Rue de l'Industrie 17 , CH-1951 Sion , Switzerland
| |
Collapse
|
16
|
Tsimpanogiannis IN, Moultos OA, Franco LFM, Spera MBDM, Erdős M, Economou IG. Self-diffusion coefficient of bulk and confined water: a critical review of classical molecular simulation studies. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1511903] [Citation(s) in RCA: 58] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Affiliation(s)
- Ioannis N. Tsimpanogiannis
- Environmental Research Laboratory, National Center for Scientific Research “Demokritos”, Aghia Paraskevi Attikis, Greece
- Institute of Nanoscience and Nanotechnology, National Center for Scientific Research “Demokritos”, Aghia Paraskevi Attikis, Greece
| | - Othonas A. Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| | - Luís F. M. Franco
- School of Chemical Engineering, University of Campinas, Campinas, Brazil
| | | | - Máté Erdős
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| | - Ioannis G. Economou
- Institute of Nanoscience and Nanotechnology, National Center for Scientific Research “Demokritos”, Aghia Paraskevi Attikis, Greece
- Chemical Engineering Program, Texas A&M University at Qatar, Doha, Qatar
| |
Collapse
|
17
|
Mahynski NA, Errington JR, Shen VK. Temperature extrapolation of multicomponent grand canonical free energy landscapes. J Chem Phys 2018; 147:054105. [PMID: 28789543 DOI: 10.1063/1.4996759] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We derive a method for extrapolating the grand canonical free energy landscape of a multicomponent fluid system from one temperature to another. Previously, we introduced this statistical mechanical framework for the case where kinetic energy contributions to the classical partition function were neglected for simplicity [N. A. Mahynski et al., J. Chem. Phys. 146, 074101 (2017)]. Here, we generalize the derivation to admit these contributions in order to explicitly illustrate the differences that result. Specifically, we show how factoring out kinetic energy effects a priori, in order to consider only the configurational partition function, leads to simpler mathematical expressions that tend to produce more accurate extrapolations than when these effects are included. We demonstrate this by comparing and contrasting these two approaches for the simple cases of an ideal gas and a non-ideal, square-well fluid.
Collapse
Affiliation(s)
- Nathan A Mahynski
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8320, USA
| | - Jeffrey R Errington
- Department of Chemical and Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200, USA
| | - Vincent K Shen
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8320, USA
| |
Collapse
|
18
|
Lötgering-Lin O, Fischer M, Hopp M, Gross J. Pure Substance and Mixture Viscosities Based on Entropy Scaling and an Analytic Equation of State. Ind Eng Chem Res 2018. [DOI: 10.1021/acs.iecr.7b04871] [Citation(s) in RCA: 55] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Oliver Lötgering-Lin
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| | - Matthias Fischer
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| | - Madlen Hopp
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| | - Joachim Gross
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| |
Collapse
|
19
|
Yadav HOS, Chakravarty C. Thiolated gold nanoparticle solvation in near-critical fluids: The role of density, temperature, and topology. J Chem Phys 2017; 146:174902. [DOI: 10.1063/1.4982755] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Hari O. S. Yadav
- Department of Chemistry, Indian Institute of Technology-Delhi, New Delhi 110016, India
| | - Charusita Chakravarty
- Department of Chemistry, Indian Institute of Technology-Delhi, New Delhi 110016, India
| |
Collapse
|
20
|
Hopp M, Gross J. Thermal Conductivity of Real Substances from Excess Entropy Scaling Using PCP-SAFT. Ind Eng Chem Res 2017. [DOI: 10.1021/acs.iecr.6b04289] [Citation(s) in RCA: 63] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Madlen Hopp
- Institute of Thermodynamics
and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring
9, 70569 Stuttgart, Germany
| | - Joachim Gross
- Institute of Thermodynamics
and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring
9, 70569 Stuttgart, Germany
| |
Collapse
|
21
|
Shrivastav G, Agarwal M, Chakravarty C, Kashyap HK. Thermodynamic regimes over which homologous alkane fluids can be treated as simple liquids. J Mol Liq 2017. [DOI: 10.1016/j.molliq.2017.01.085] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
22
|
Olsen AE, Dyre JC, Schrøder TB. Communication: Pseudoisomorphs in liquids with intramolecular degrees of freedom. J Chem Phys 2016; 145:241103. [DOI: 10.1063/1.4972860] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Andreas Elmerdahl Olsen
- 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
|
23
|
Dhabal D, Chakravarty C, Molinero V, Kashyap HK. Comparison of liquid-state anomalies in Stillinger-Weber models of water, silicon, and germanium. J Chem Phys 2016; 145:214502. [DOI: 10.1063/1.4967939] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- Debdas Dhabal
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Charusita Chakravarty
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, USA
| | - Hemant K. Kashyap
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| |
Collapse
|
24
|
Lötgering-Lin O, Schöniger A, Nowak W, Gross J. Bayesian Model Selection Helps To Choose Objectively between Thermodynamic Models: A Demonstration of Selecting a Viscosity Model Based on Entropy Scaling. Ind Eng Chem Res 2016. [DOI: 10.1021/acs.iecr.6b02671] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Oliver Lötgering-Lin
- Institute
of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| | - Anneli Schöniger
- Center
for Applied Geoscience, University of Tübingen, Hölderlinstraße 12, 72074 Tübingen, Germany
| | - Wolfgang Nowak
- Institute
for Modelling Hydraulic and Environmental Systems (LS3)/SimTech, University of Stuttgart, Pfaffenwaldring 5a, 70569 Stuttgart, Germany
| | - Joachim Gross
- Institute
of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| |
Collapse
|
25
|
Yadav HOS, Shrivastav G, Agarwal M, Chakravarty C. Effective interactions between nanoparticles: Creating temperature-independent solvation environments for self-assembly. J Chem Phys 2016; 144:244901. [DOI: 10.1063/1.4954325] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
26
|
Jeong C, Douglas JF. Mass dependence of the activation enthalpy and entropy of unentangled linear alkane chains. J Chem Phys 2015; 143:144905. [DOI: 10.1063/1.4932601] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
|
27
|
Lötgering-Lin O, Gross J. Group Contribution Method for Viscosities Based on Entropy Scaling Using the Perturbed-Chain Polar Statistical Associating Fluid Theory. Ind Eng Chem Res 2015. [DOI: 10.1021/acs.iecr.5b01698] [Citation(s) in RCA: 88] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Oliver Lötgering-Lin
- Institute of Thermodynamics
and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring
9, 70569 Stuttgart, Germany
| | - Joachim Gross
- Institute of Thermodynamics
and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring
9, 70569 Stuttgart, Germany
| |
Collapse
|
28
|
Grzybowska K, Grzybowski A, Pawlus S, Pionteck J, Paluch M. Role of entropy in the thermodynamic evolution of the time scale of molecular dynamics near the glass transition. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:062305. [PMID: 26172717 DOI: 10.1103/physreve.91.062305] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/13/2015] [Indexed: 06/04/2023]
Abstract
In this paper, we investigate how changes in the system entropy influence the characteristic time scale of the system molecular dynamics near the glass transition. Independently of any model of thermodynamic evolution of the time scale, against some previous suppositions, we show that the system entropy S is not sufficient to govern the time scale defined by structural relaxation time τ. In the density scaling regime, we argue that the decoupling between τ and S is a consequence of different values of the scaling exponents γ and γ(S) in the density scaling laws, τ=f(ρ(γ)/T) and S=h(ρ(γ(S))/T), where ρ and T denote density and temperature, respectively. It implies that the proper relation between τ and S requires supplementing with a density factor, u(ρ), i.e., τ=g(u(ρ)w(S)). This meaningful finding additionally demonstrates that the density scaling idea can be successfully used to separate physically relevant contributions to the time scale of molecular dynamics near the glass transition. The relation reported by us between τ and S constitutes a general pattern based on nonconfigurational quantities for describing the thermodynamic evolution of the characteristic time scale of molecular dynamics near the glass transition in the density scaling regime, which is a promising alternative to the approaches based as the Adam-Gibbs model on the configurational entropy that is difficult to evaluate in the entire thermodynamic space. As an example, we revise the Avramov entropic model of the dependence τ(T,ρ), giving evidence that its entropic basis has to be extended by the density dependence of the maximal energy barrier for structural relaxation. We also discuss the excess entropy S(ex), the density scaling of which is found to mimic the density scaling of the total system entropy S.
Collapse
Affiliation(s)
- K Grzybowska
- Institute of Physics, University of Silesia, Uniwersytecka 4, 40-007 Katowice, Poland
- Silesian Center for Education and Interdisciplinary Research, 75 Pułku Piechoty 1A, 41-500 Chorzów, Poland
| | - A Grzybowski
- Institute of Physics, University of Silesia, Uniwersytecka 4, 40-007 Katowice, Poland
- Silesian Center for Education and Interdisciplinary Research, 75 Pułku Piechoty 1A, 41-500 Chorzów, Poland
| | - S Pawlus
- Institute of Physics, University of Silesia, Uniwersytecka 4, 40-007 Katowice, Poland
- Silesian Center for Education and Interdisciplinary Research, 75 Pułku Piechoty 1A, 41-500 Chorzów, Poland
| | - J Pionteck
- Leibniz Institute of Polymer Research Dresden, Hohe Strasse 6, D-01069 Dresden, Germany
| | - M Paluch
- Institute of Physics, University of Silesia, Uniwersytecka 4, 40-007 Katowice, Poland
- Silesian Center for Education and Interdisciplinary Research, 75 Pułku Piechoty 1A, 41-500 Chorzów, Poland
| |
Collapse
|
29
|
Abstract
Recent developments show that many liquids and solids have an approximate "hidden" scale invariance that implies the existence of lines in the thermodynamic phase diagram, so-called isomorphs, along which structure and dynamics in properly reduced units are invariant to a good approximation. This means that the phase diagram becomes effectively one-dimensional with regard to several physical properties. Liquids and solids with isomorphs include most or all van der Waals bonded systems and metals, as well as weakly ionic or dipolar systems. On the other hand, systems with directional bonding (hydrogen bonds or covalent bonds) or strong Coulomb forces generally do not exhibit hidden scale invariance. The article reviews the theory behind this picture of condensed matter and the evidence for it coming from computer simulations and experiments.
Collapse
Affiliation(s)
- Jeppe C Dyre
- DNRF Center "Glass and Time", IMFUFA, Department of Sciences, Roskilde University , P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
30
|
Prasad S, Chakravarty C. Onset of simple liquid behaviour in modified water models. J Chem Phys 2014; 140:164501. [DOI: 10.1063/1.4870823] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
31
|
Voyiatzis E, Müller-Plathe F, Böhm MC. Excess entropy scaling for the segmental and global dynamics of polyethylene melts. Phys Chem Chem Phys 2014; 16:24301-11. [DOI: 10.1039/c4cp03559c] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The range of validity of the Rosenfeld and Dzugutov excess entropy scaling laws is analyzed for unentangled linear polyethylene chains.
Collapse
Affiliation(s)
- Evangelos Voyiatzis
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie and Center of Smart Interfaces
- Technische Universität Darmstadt
- D-64287 Darmstadt, Germany
| | - Florian Müller-Plathe
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie and Center of Smart Interfaces
- Technische Universität Darmstadt
- D-64287 Darmstadt, Germany
| | - Michael C. Böhm
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie and Center of Smart Interfaces
- Technische Universität Darmstadt
- D-64287 Darmstadt, Germany
| |
Collapse
|
32
|
Voyiatzis E, Müller-Plathe F, Böhm MC. Do Transport Properties of Entangled Linear Polymers Scale with Excess Entropy? Macromolecules 2013. [DOI: 10.1021/ma401617z] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Evangelos Voyiatzis
- Eduard-Zintl-Institut fur
Anorganische und Physikalische Chemie and Center of Smart Interfaces, Technische Universität Darmstadt, Petersenstrasse 20, D-64287Darmstadt, Germany
| | - Florian Müller-Plathe
- Eduard-Zintl-Institut fur
Anorganische und Physikalische Chemie and Center of Smart Interfaces, Technische Universität Darmstadt, Petersenstrasse 20, D-64287Darmstadt, Germany
| | - Michael C. Böhm
- Eduard-Zintl-Institut fur
Anorganische und Physikalische Chemie and Center of Smart Interfaces, Technische Universität Darmstadt, Petersenstrasse 20, D-64287Darmstadt, Germany
| |
Collapse
|
33
|
Shell MS. Systematic coarse-graining of potential energy landscapes and dynamics in liquids. J Chem Phys 2013; 137:084503. [PMID: 22938246 DOI: 10.1063/1.4746391] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Recent efforts have shown that the dynamic properties of a wide class of liquids can be mapped onto semi-universal scaling laws and constitutive relations that are motivated by thermodynamic analyses of much simpler models. In particular, it has been found that many systems exhibit dynamics whose behavior in state space closely follows that of soft-sphere particles interacting through an inverse power repulsion. In the present work, we show that a recently developed coarse-graining theory provides a natural way to understand how arbitrary liquids can be mapped onto effective soft-sphere models and hence how one might potentially be able to extract underlying dynamical scaling laws. The theory is based on the relative entropy, an information metric that quantifies how well a soft-sphere approximation to a liquid's multidimensional potential energy landscape performs. We show that optimization of the relative entropy not only enables one to extract effective soft-sphere potentials that suggest an inherent scaling of thermodynamic and dynamic properties in temperature-density space, but that also has rather interesting connections to excess entropy based theories of liquid dynamics. We apply the approach to a binary mixture of Lennard-Jones particles, and show that it gives effective soft-sphere scaling laws that well-describe the behavior of the diffusion constants. Our results suggest that the relative entropy formalism may be useful for "perturbative" type theories of dynamics, offering a general strategy for systematically connecting complex energy landscapes to simpler reference ones with better understood dynamic behavior.
Collapse
Affiliation(s)
- M Scott Shell
- Department of Chemical Engineering, University of California Santa Barbara, Santa Barbara, California 93106-5080, USA.
| |
Collapse
|
34
|
Vaz RV, Magalhães AL, Fernandes DL, Silva CM. Universal correlation of self-diffusion coefficients of model and real fluids based on residual entropy scaling law. Chem Eng Sci 2012. [DOI: 10.1016/j.ces.2012.05.006] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
35
|
Abstract
Isomorphs are curves in the phase diagram along which a number of static and dynamic quantities are invariant in reduced units (Gnan, N.; et al. J. Chem. Phys.2009, 131, 234504). A liquid has good isomorphs if and only if it is strongly correlating, i.e., if the equilibrium virial/potential energy fluctuations are more than 90% correlated in the NVT ensemble. Isomorphs were previously discussed with a focus on atomic systems. This paper generalizes isomorphs to liquids composed of rigid molecules and study the isomorphs of systems of small rigid molecules: the asymmetric dumbbell model, a symmetric inverse power-law dumbbell, and the Lewis-Wahnström o-terphenyl (OTP) model. For all model systems, the following quantities are found to a good approximation to be invariant along an isomorph: the isochoric heat capacity, the excess entropy, the reduced molecular center-of-mass self-part of the intermediate scattering function, and the reduced molecular center-of-mass radial distribution function. In agreement with theory, we also find that an instantaneous change of temperature and density from an equilibrated state point to an isomorphic state point leads to no relaxation. The isomorphs of the Lewis-Wahnström OTP model were found to be more approximative than those of the asymmetric dumbbell model; this is consistent with the OTP model being less strongly correlating. The asymmetric dumbbell and Lewis-Wahnström OTP models each have a "master isomorph"; i.e., the isomorphs have identical shape in the virial/potential energy phase diagram.
Collapse
Affiliation(s)
- Trond S Ingebrigtsen
- DNRF Centre Glass and Time, IMFUFA, Department of Sciences, Roskilde University, Postbox 260, DK-4000 Roskilde, Denmark
| | | | | |
Collapse
|
36
|
Pond MJ, Errington JR, Truskett TM. Communication: Generalizing Rosenfeld's excess-entropy scaling to predict long-time diffusivity in dense fluids of Brownian particles: from hard to ultrasoft interactions. J Chem Phys 2011; 134:081101. [PMID: 21361518 DOI: 10.1063/1.3559676] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Computer simulations are used to test whether a recently introduced generalization of Rosenfeld's excess-entropy scaling method for estimating transport coefficients in systems obeying molecular dynamics can be extended to predict long-time diffusivities in fluids of particles undergoing Brownian dynamics in the absence of interparticle hydrodynamic forces. Model fluids with inverse-power-law, Gaussian-core, and Hertzian pair interactions are considered. Within the generalized Rosenfeld scaling method, long-time diffusivities of ultrasoft Gaussian-core and Hertzian particle fluids, which display anomalous trends with increasing density, are predicted (to within 20%) based on knowledge of interparticle interactions, excess entropy, and scaling behavior of simpler inverse-power-law fluids.
Collapse
Affiliation(s)
- Mark J Pond
- Department of Chemical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA
| | | | | |
Collapse
|