1
|
Busch J, Paschek D. Computing Accurate True Self-Diffusion Coefficients and Shear Viscosities Using the OrthoBoXY Approach. J Phys Chem B 2024; 128:1040-1052. [PMID: 38240259 DOI: 10.1021/acs.jpcb.3c07540] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2024]
Abstract
In a recent paper [Busch, J.; Paschek, D. J. Phys. Chem. B 2023, 127, 7983-7987], we have shown that for molecular dynamics (MD) simulations using orthorhombic periodic boundary conditions with "magic" box length ratios of Lz/Lx = Lz/Ly = 2.7933596497, the self-diffusion coefficients Dx and Dy in x- and y-directions are independent of the system size. They both represent the true self-diffusion coefficient D0 = (Dx + Dy)/2, while the shear viscosity can be calculated from diffusion coefficients in x-, y-, and z-directions, using η = kBT·8.1711245653/[3πLz(Dx + Dy - 2Dz)]. In this contribution, we test this "OrthoBoXY" approach by its application to a variety of different systems: liquid water, dimethyl ether, methanol, triglyme, water/methanol mixtures, water/triglyme mixtures, and imidazolium-based ionic liquids. The chosen systems range from small-sized molecular liquids to complex mixtures and ionic liquids, while spanning a viscosity range of almost 3 orders of magnitude. We assess the efficiency of the method for computing true self-diffusion and viscosity data and provide simple formulas for estimating the required MD simulation lengths and sizes for delivering reliable data with targeted uncertainty levels. Our analysis of the system size dependence of statistical uncertainties for both the viscosity and the self-diffusion coefficient leads us to the conclusion that it is preferable to extend the simulation length instead of increasing the system size. MD simulations consisting of 768 molecules or ion pairs seem to be perfectly adequate.
Collapse
Affiliation(s)
- Johanna Busch
- Institut für Chemie, Abteilung Physikalische und Theoretische Chemie, Universität Rostock, Albert-Einstein-Str. 27, D-18059 Rostock, Germany
| | - Dietmar Paschek
- Institut für Chemie, Abteilung Physikalische und Theoretische Chemie, Universität Rostock, Albert-Einstein-Str. 27, D-18059 Rostock, Germany
| |
Collapse
|
2
|
Maffioli L, Smith ER, Ewen JP, Daivis PJ, Dini D, Todd BD. Slip and stress from low shear rate nonequilibrium molecular dynamics: The transient-time correlation function technique. J Chem Phys 2022; 156:184111. [PMID: 35568555 DOI: 10.1063/5.0088127] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We derive the transient-time correlation function (TTCF) expression for the computation of phase variables of inhomogenous confined atomistic fluids undergoing boundary-driven planar shear (Couette) flow at constant pressure. Using nonequilibrium molecular dynamics simulations, we then apply the TTCF formalism to the computation of the shear stress and the slip velocity for atomistic fluids at realistic low shear rates, in systems under constant pressure and constant volume. We show that, compared to direct averaging of multiple trajectories, the TTCF method dramatically improves the accuracy of the results at low shear rates and that it is suitable to investigate the tribology and rheology of atomistically detailed confined fluids at realistic flow rates.
Collapse
Affiliation(s)
- Luca Maffioli
- Department of Mathematics, School of Science, Computing and Engineering Technologies, Swinburne University of Technology, P.O. Box 218, Hawthorn, Victoria 3122, Australia
| | - Edward R Smith
- Mechanical and Aerospace Engineering, Brunel University London, Kingston Lane, Uxbridge UB8 3PH, United Kingdom
| | - James P Ewen
- Department of Mechanical Engineering, Imperial College London, South Kensington Campus, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Peter J Daivis
- School of Science, RMIT University, GPO Box 2476, Victoria 3001, Australia
| | - Daniele Dini
- Department of Mechanical Engineering, Imperial College London, South Kensington Campus, Exhibition Road, London SW7 2AZ, United Kingdom
| | - B D Todd
- Department of Mathematics, School of Science, Computing and Engineering Technologies, Swinburne University of Technology, P.O. Box 218, Hawthorn, Victoria 3122, Australia
| |
Collapse
|
3
|
Desgranges C, Delhommelle J. Unraveling liquid polymorphism in silicon driven out-of-equilibrium. J Chem Phys 2020; 153:054502. [DOI: 10.1063/5.0015417] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Caroline Desgranges
- Department of Chemistry, New York University, New York, New York 10003, USA and Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| | - Jerome Delhommelle
- Department of Chemistry, New York University, New York, New York 10003, USA and Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| |
Collapse
|
4
|
Celebi AT, Nguyen CT, Hartkamp R, Beskok A. The role of water models on the prediction of slip length of water in graphene nanochannels. J Chem Phys 2019; 151:174705. [DOI: 10.1063/1.5123713] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Alper Tunga Celebi
- Process and Energy Department, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | - Chinh Thanh Nguyen
- Lyle School of Engineering, Southern Methodist University, 3101 Dyer Street, Dallas, Texas 75205, USA
| | - Remco Hartkamp
- Process and Energy Department, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | - Ali Beskok
- Lyle School of Engineering, Southern Methodist University, 3101 Dyer Street, Dallas, Texas 75205, USA
| |
Collapse
|
5
|
Jamali SH, Hartkamp R, Bardas C, Söhl J, Vlugt TJH, Moultos OA. Shear Viscosity Computed from the Finite-Size Effects of Self-Diffusivity in Equilibrium Molecular Dynamics. J Chem Theory Comput 2018; 14:5959-5968. [PMID: 30296092 PMCID: PMC6236468 DOI: 10.1021/acs.jctc.8b00625] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
A method is proposed for calculating
the shear viscosity of a liquid
from finite-size effects of self-diffusion coefficients in Molecular
Dynamics simulations. This method uses the difference in the self-diffusivities,
computed from at least two system sizes, and an analytic equation
to calculate the shear viscosity. To enable the efficient use of this
method, a set of guidelines is developed. The most efficient number
of system sizes is two and the large system is at least four times
the small system. The number of independent simulations for each system
size should be assigned in such a way that 50%–70% of the total
available computational resources are allocated to the large system.
We verified the method for
250 binary and 26 ternary Lennard-Jones systems, pure water, and an
ionic liquid ([Bmim][Tf2N]). The computed shear viscosities
are in good agreement with viscosities obtained from equilibrium Molecular
Dynamics simulations for all liquid systems far from the critical
point. Our results indicate that the proposed method is suitable for
multicomponent mixtures and highly viscous liquids. This may enable
the systematic screening of the viscosities of ionic liquids and deep
eutectic solvents.
Collapse
Affiliation(s)
- Seyed Hossein Jamali
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft University of Technology , Leeghwaterstraat 39 , 2628CB Delft , The Netherlands
| | - Remco Hartkamp
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft University of Technology , Leeghwaterstraat 39 , 2628CB Delft , The Netherlands
| | - Christos Bardas
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft University of Technology , Leeghwaterstraat 39 , 2628CB Delft , The Netherlands
| | - Jakob Söhl
- Delft Institute of Applied Mathematics , Delft University of Technology , van Mourik Broekmanweg 6 , 2628XE Delft , The Netherlands
| | - Thijs J H Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft University of Technology , Leeghwaterstraat 39 , 2628CB Delft , The Netherlands
| | - Othonas A Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft University of Technology , Leeghwaterstraat 39 , 2628CB Delft , The Netherlands
| |
Collapse
|
6
|
Ciccotti G, Ferrario M. Non-equilibrium by molecular dynamics: a dynamical approach. MOLECULAR SIMULATION 2016. [DOI: 10.1080/08927022.2015.1121543] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
7
|
|
8
|
Jain A, Sasmal C, Hartkamp R, Todd B, Prakash JR. Brownian dynamics simulations of planar mixed flows of polymer solutions at finite concentrations. Chem Eng Sci 2015. [DOI: 10.1016/j.ces.2014.09.035] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
9
|
Hartkamp R, Coasne B. Structure and transport of aqueous electrolytes: From simple halides to radionuclide ions. J Chem Phys 2014; 141:124508. [DOI: 10.1063/1.4896380] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
10
|
|
11
|
Hartkamp R, Todd BD, Luding S. A constitutive framework for the non-Newtonian pressure tensor of a simple fluid under planar flows. J Chem Phys 2013; 138:244508. [DOI: 10.1063/1.4810746] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
12
|
Costa D, Sergi A, Ferrario M. Transient behavior of a model fluid under applied shear. J Chem Phys 2013; 138:184501. [PMID: 23676050 DOI: 10.1063/1.4803147] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We study the transient behavior of a model fluid composed by soft repulsive spheres subjected to a planar uniform shear. To this aim, we use a dynamical non-equilibrium molecular dynamics method originally developed by Ciccotti and Jacucci [Phys. Rev. Lett. 35, 789 (1975)] and recently applied to the study of the transient regimes in various fluid systems. We show that the dynamical method allows one to study the transient behavior of the viscous time-dependent response over a wide range of applied shear rates, provided that a temperature control is enforced on the system. In this study, we adopt in particular the configurational thermostat of Braga and Travis [J. Chem. Phys. 123, 134101 (2005)]. The initial behavior of the dynamical response to a θ-like perturbation is characterized by a rapid increase, culminating in a pronounced peak, later relaxing to a plateau value. The latter positively reproduces the values of the viscosity observed in standard steady-state non-equilibrium molecular dynamics.
Collapse
Affiliation(s)
- Dino Costa
- Dipartimento di Fisica e di Scienze della Terra, Università degli Studi di Messina Viale F. Stagno d'Alcontres 31, 98166 Messina, Italy
| | | | | |
Collapse
|
13
|
Bernardi S, Brookes SJ, Searles DJ, Evans DJ. Response theory for confined systems. J Chem Phys 2012; 137:074114. [DOI: 10.1063/1.4746121] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|