1
|
Nitzke I, Stephan S, Vrabec J. Topology of thermodynamic potentials using physical models: Helmholtz, Gibbs, Grand, and Null. J Chem Phys 2024; 160:214104. [PMID: 38828811 DOI: 10.1063/5.0207592] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2024] [Accepted: 05/17/2024] [Indexed: 06/05/2024] Open
Abstract
Thermodynamic potentials play a substantial role in numerous scientific disciplines and serve as basic constructs for describing the behavior of matter. Despite their significance, comprehensive investigations of their topological characteristics and their connections to molecular interactions have eluded exploration due to experimental inaccessibility issues. This study addresses this gap by analyzing the topology of the Helmholtz energy, Gibbs energy, Grand potential, and Null potential that are associated with different isothermal boundary conditions. By employing Monte Carlo simulations in the NVT, NpT, and μVT ensembles and a molecular-based equation of state, methane, ethane, nitrogen, and methanol are investigated over a broad range of thermodynamic conditions. The predictions from the two independent methods are overall in very good agreement. Although distinct quantitative differences among the fluids are observed, the overall topology of the individual thermodynamic potentials remains unaffected by the molecular architecture, which is in line with the corresponding states principle-as expected. Furthermore, a comparative analysis reveals significant differences between the total potentials and their residual contributions.
Collapse
Affiliation(s)
- Isabel Nitzke
- Thermodynamics, Technische Universität Berlin, 10587 Berlin, Germany
| | - Simon Stephan
- Laboratory of Engineering Thermodynamics, RPTU Kaiserslautern, 67663 Kaiserslautern, Germany
| | - Jadran Vrabec
- Thermodynamics, Technische Universität Berlin, 10587 Berlin, Germany
| |
Collapse
|
2
|
Schwerdtfeger P, Wales DJ. 100 Years of the Lennard-Jones Potential. J Chem Theory Comput 2024; 20:3379-3405. [PMID: 38669689 DOI: 10.1021/acs.jctc.4c00135] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/28/2024]
Abstract
It is now 100 years since Lennard-Jones published his first paper introducing the now famous potential that bears his name. It is therefore timely to reflect on the many achievements, as well as the limitations, of this potential in the theory of atomic and molecular interactions, where applications range from descriptions of intermolecular forces to molecules, clusters, and condensed matter.
Collapse
Affiliation(s)
- Peter Schwerdtfeger
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study, Massey University Auckland, Private Bag 102904, Auckland 0745, New Zealand
| | - David J Wales
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K
| |
Collapse
|
3
|
Paterson T, Bannerman MN, Lue L. Using the Zeno line to assess and refine molecular models. J Chem Phys 2024; 160:154503. [PMID: 38624125 DOI: 10.1063/5.0192770] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2023] [Accepted: 04/01/2024] [Indexed: 04/17/2024] Open
Abstract
The Zeno line is the locus of points on the temperature-density plane where the compressibility factor of the fluid is equal to one. It has been observed to be straight for a broad variety of real fluids, although the underlying reasons for this are still unclear. In this work, a detailed study of the Zeno line and its relation to the vapor-liquid coexistence curve is performed for two simple model pair-potential fluids: attractive square-well fluids with varying well-widths λ and Mie n-6 fluids with different repulsive exponents n. Interestingly, the Zeno lines of these fluids are curved, regardless of the value of λ or n. We find that for square-well fluids, λ ≈ 1.8 presents a Zeno line, which is the most linear over the largest temperature range. For Mie n-6 fluids, we find that the straightest Zeno line occurs for n between 8 and 10. Additionally, the square-well and Mie fluids with the straightest Zeno line showed the closest quantitative agreement with the vapor-liquid coexistence curve for experimental fluids that follow the principle of corresponding states (e.g., argon, xenon, krypton, methane, nitrogen, and oxygen). These results suggest that the Zeno line can provide a useful additional feature, in complement to other properties, such as the phase envelope, to evaluate molecular models.
Collapse
Affiliation(s)
- Thomas Paterson
- Department of Chemical and Process Engineering, University of Strathclyde, James Weir Building, 75 Montrose Street, Glasgow G1 1XJ, United Kingdom
| | - Marcus N Bannerman
- School of Engineering, University of Aberdeen, Aberdeen AB24 3UE, United Kingdom
| | - Leo Lue
- Department of Chemical and Process Engineering, University of Strathclyde, James Weir Building, 75 Montrose Street, Glasgow G1 1XJ, United Kingdom
| |
Collapse
|
4
|
Korneva M, Zhilyaev P. Solid-liquid phase transition inside van der Waals nanobubbles: an atomistic perspective. Phys Chem Chem Phys 2023. [PMID: 37432424 DOI: 10.1039/d3cp01285a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/12/2023]
Abstract
The liquid-solid phase transition during the confinement of a van der Waals bubble is studied using molecular dynamics simulations. In particular, argon is considered inside a graphene bubble, where the outer membrane is a sheet of graphene, and the substrate is atomically flat graphite. A methodology to avoid metastable states of argon is developed and implemented to derive a melting curve of trapped argon. It is found that in the confinement, the melting curve of argon shifts toward higher temperatures, and the temperature shift is about 10-30 K. The ratio of the height to the radius of the GNB (H/R) decreases with increasing temperature. It also most likely undergoes an abrupt change through the liquid-crystal phase transition. The semi-liquid state of argon was detected in the transition region. At this state, the argon structure stays layered, but the atoms travel distances of several lattice constants.
Collapse
Affiliation(s)
- Mariia Korneva
- Center for Materials Technologies, Skolkovo Institute of Science and Technology, Skolkovo Innovation Center, Building 3, Moscow, 143026, Russia.
| | - Petr Zhilyaev
- Center for Materials Technologies, Skolkovo Institute of Science and Technology, Skolkovo Innovation Center, Building 3, Moscow, 143026, Russia.
| |
Collapse
|
5
|
Antolović I, Staubach J, Stephan S, Vrabec J. Phase equilibria of symmetric Lennard-Jones mixtures and a look at the transport properties near the upper critical solution temperature. Phys Chem Chem Phys 2023. [PMID: 37365979 DOI: 10.1039/d3cp01434g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/28/2023]
Abstract
This study investigates phase equilibria and transport properties of five symmetric binary Lennard-Jones mixtures using molecular simulation and equation of state models. The mixtures are selected for their representation of different types of phase behavior and the research contributes to the development of simulation techniques, mixture theories and understanding of thermophysical mixture properties. A novel method is introduced for determining the critical end point (CEP) and critical azeotropic end point (CAEP) by molecular simulation. The van der Waals one-fluid theory is assessed for its performance in conjunction with Lennard-Jones equation of state models, while addressing different phase equilibrium types simultaneously. An empirical correlation is introduced to account for deviations between the equation of state and simulation that arise when using the same binary interaction parameter. This study also investigates the influence of the liquid-liquid critical point on thermophysical properties, which are found to exhibit no significant anomalies or singularities. System-size effects of diffusion coefficients are addressed by extrapolating simulation data to the thermodynamic limit and applying analytical finite-size corrections.
Collapse
Affiliation(s)
- Ivan Antolović
- Thermodynamics, Technical University Berlin, 10587 Berlin, Germany.
| | - Jens Staubach
- Laboratory of Engineering Thermodynamics (LTD), RPTU Kaiserslautern, Kaiserslautern, Germany
| | - Simon Stephan
- Laboratory of Engineering Thermodynamics (LTD), RPTU Kaiserslautern, Kaiserslautern, Germany
| | - Jadran Vrabec
- Thermodynamics, Technical University Berlin, 10587 Berlin, Germany.
| |
Collapse
|
6
|
Chaparro G, Müller EA. Development of thermodynamically consistent machine-learning equations of state: Application to the Mie fluid. J Chem Phys 2023; 158:2890032. [PMID: 37161943 DOI: 10.1063/5.0146634] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2023] [Accepted: 04/24/2023] [Indexed: 05/11/2023] Open
Abstract
A procedure for deriving thermodynamically consistent data-driven equations of state (EoS) for fluids is presented. The method is based on fitting the Helmholtz free energy using artificial neural networks to obtain a closed-form relationship between the thermophysical properties of fluids (FE-ANN EoS). As a proof-of-concept, an FE-ANN EoS is developed for the Mie fluids, starting from a database obtained by classical molecular dynamics simulations. The FE-ANN EoS is trained using first- (pressure and internal energy) and second-order (e.g., heat capacities, Joule-Thomson coefficients) derivative data. Additional constraints ensure that the data-driven model fulfills thermodynamically consistent limits and behavior. The results for the FE-ANN EoS are shown to be as accurate as the best available analytical model while being developed in a fraction of the time. The robustness of the "digital" equation of state is exemplified by computing physical behavior it has not been trained on, for example, fluid phase equilibria. Furthermore, the model's internal consistency is successfully assessed using Brown's characteristic curves.
Collapse
Affiliation(s)
- Gustavo Chaparro
- Department of Chemical Engineering, Sargent Centre for Process Systems Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - Erich A Müller
- Department of Chemical Engineering, Sargent Centre for Process Systems Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| |
Collapse
|
7
|
Urschel M, Stephan S. Determining Brown's Characteristic Curves Using Molecular Simulation. J Chem Theory Comput 2023; 19:1537-1552. [PMID: 36802608 DOI: 10.1021/acs.jctc.2c01102] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/23/2023]
Abstract
Brown's characteristic curves define lines on the thermodynamic surface where special thermodynamic conditions hold. These curves are an important tool for the development of thermodynamic models of fluids. Yet, practically no experimental data for Brown's characteristic curves is available. In this work, a rigorous and generalized method for determining Brown's characteristic curves based on molecular simulation was developed. As multiple thermodynamic equivalent definitions apply for the characteristic curves, different simulation routes were compared. Based on this systematic approach, the most favorable route for determining each characteristic curve was identified. The computational procedure developed in this work combines molecular simulation, molecular-based equation of state, and the evaluation of the second virial coefficient. The new method was tested on a simple model system (the classical Lennard-Jones fluid) and different types of real substances (toluene, methane, ethane, propane, and ethanol). It is thereby shown that the method is robust and yields accurate results. Moreover, a computer code implementation of the method is presented.
Collapse
Affiliation(s)
- Maximilian Urschel
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern 67663, Germany
| | - Simon Stephan
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern 67663, Germany
| |
Collapse
|
8
|
Pohl S, Fingerhut R, Thol M, Vrabec J, Span R. Equation of state for the Mie (λ r,6) fluid with a repulsive exponent from 11 to 13. J Chem Phys 2023; 158:084506. [PMID: 36859099 DOI: 10.1063/5.0133412] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/03/2023] Open
Abstract
An empirical multi-parameter equation of state in terms of the reduced Helmholtz energy is presented for the Mie (λr-6) fluid with a repulsive exponent λr from 11 to 13. The equation is fitted to an extensive dataset from molecular dynamics simulation as well as the second and third thermal virial coefficients. It is comprehensively compared with the SAFT-VR model and is a more accurate description of the considered fluid class. The equation is valid for reduced temperatures T/Tc from 0.55 to 4.5 and for reduced pressures of up to p/pc = 265. A good extrapolation behavior and the occurrence of a single Maxwell loop down to the vicinity of the triple point temperature are realized.
Collapse
Affiliation(s)
- Sven Pohl
- Thermodynamics, Ruhr-University Bochum, 44801 Bochum, Germany
| | - Robin Fingerhut
- Thermodynamics, Technical University of Berlin, 10587 Berlin, Germany
| | - Monika Thol
- Thermodynamics, Ruhr-University Bochum, 44801 Bochum, Germany
| | - Jadran Vrabec
- Thermodynamics, Technical University of Berlin, 10587 Berlin, Germany
| | - Roland Span
- Thermodynamics, Ruhr-University Bochum, 44801 Bochum, Germany
| |
Collapse
|
9
|
Fertig D, Stephan S. Influence of dispersive long-range interactions on transport and excess properties of simple mixtures. Mol Phys 2023. [DOI: 10.1080/00268976.2022.2162993] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- David Fertig
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| | - Simon Stephan
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| |
Collapse
|
10
|
Alsaifi NM, Elliott JR. Avoiding Artifacts in Noncubic Equations of State. Ind Eng Chem Res 2022. [DOI: 10.1021/acs.iecr.2c01923] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Nayef M. Alsaifi
- Center for Refining and Advanced Chemicals, Chemical Engineering Department, King Fahd University of Petroleum & Minerals, Dhahran31261, Saudi Arabia
| | - J. Richard Elliott
- Chemical and Biomolecular Engineering Department, The University of Akron, Akron, Ohio44325-3906, United States
| |
Collapse
|
11
|
Apfelbaum EM. The Line of Ideal Isothermal Compressibility. J Phys Chem B 2022; 126:2912-2920. [PMID: 35389650 DOI: 10.1021/acs.jpcb.2c01247] [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
We have considered the line along which the values of the isothermal compressibility of a system are the same as they would be for an ideal gas. It was called the κT line. Various substances and models have been studied with the use of the multi-parameter equations of states, implemented in REFPROP (Lemmon, E. W.; Bell, I. H.; Huber, M. L.; McLinden, M. O. NIST standard reference database 23: Reference fluid thermodynamic and transport properties-REFPROP, version 10.0 (National Institute of Standards and Technology, 2018)), and numerical simulations, respectively. It is shown that in the reduced variables, all κT lines can be described by the same relation following from the van der Waals equation. This similarity appears to be valid even for the systems with non-straight Zeno lines. Deviations within ∼5% are observed only near the binodals.
Collapse
Affiliation(s)
- E M Apfelbaum
- Joint Institute for High Temperatures of the Russian Academy of Sciences, Izhorskaya Street 13 Bldg. 2, Moscow 125412, Russia
| |
Collapse
|
12
|
Stephan S, Dyga M, Alabd Alhafez I, Lenhard J, Urbassek HM, Hasse H. Reproducibility of atomistic friction computer experiments: a molecular dynamics simulation study. MOLECULAR SIMULATION 2021. [DOI: 10.1080/08927022.2021.1987430] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Simon Stephan
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| | - Maximilian Dyga
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| | - Iyad Alabd Alhafez
- Physics Department and Research Center (OPTIMAS), TU Kaiserslautern, Kaiserslautern, Germany
| | - Johannes Lenhard
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| | - Herbert M. Urbassek
- Physics Department and Research Center (OPTIMAS), TU Kaiserslautern, Kaiserslautern, Germany
| | - Hans Hasse
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| |
Collapse
|