Abstract:In a recent paper [J. Phys. Chem.B 127, 7983-7987 (2023)] , we have shown that for molecular dynamics (MD) simulations using orthorhombic periodic boundary conditions with "magic" box length ratios of $L_z/L_x=L_z/L_y=2.7933596497$, the self-diffusion coefficients $D_x$ and $D_y$ in $x$- and $y$-direction are independent from the system size. They both represent the true self-diffusion coefficient $D_0=(D_x+D_y)/2$, while the shear viscosity can be calculated from diffusion coefficients in $x$-,$y$- and $z$-direction, using $\eta=k_\mathrm{B}T\cdot 8.1711245653/[3\pi L_z(D_{x}+D_{y}-2D_z)]$. 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 three orders of magnitude. We assess the efficiency of the method for computing true self-diffusion and viscosity data and provide simple formulae 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.

What problem does this paper attempt to address?

### Problems the Paper Aims to Solve The paper aims to address the issue of system size dependence encountered when calculating the true self-diffusion coefficient and shear viscosity of fluids through molecular dynamics (MD) simulations. Specifically, the paper proposes a method called "OrthoBoXY," which uses orthogonal periodic boundary conditions (PBCs) with a specific aspect ratio (Lz/Lx = Lz/Ly ≈ 2.7933596497) to ensure that the self-diffusion coefficients in the x and y directions are not affected by the system size, thereby allowing accurate calculation of the true self-diffusion coefficient $ D_0 $ and shear viscosity $ \eta $. ### Background and Motivation 1. **System Size Dependence of Self-Diffusion Coefficient**: - In traditional MD simulations, the self-diffusion coefficient is often affected by the system size due to hydrodynamic interactions under periodic boundary conditions, leading to inaccurate results. - This dependence manifests as the self-diffusion coefficient varying with the inverse of the system size (L^(-1)), where L is the edge length of the cubic MD cell. 2. **Calculation of Shear Viscosity**: - Calculating shear viscosity typically requires knowledge of the self-diffusion coefficient or a series of MD simulations with different box lengths. - Traditional methods like the Green-Kubo and Einstein relations are effective but require very long simulation times to obtain sufficient statistical samples. ### Proposed Method 1. **OrthoBoXY Method**: - By selecting a specific aspect ratio (Lz/Lx = Lz/Ly ≈ 2.7933596497), the hydrodynamic interactions in the x and y directions can be eliminated, making the self-diffusion coefficients $ D_x $ and $ D_y $ represent the true self-diffusion coefficient $ D_0 $. - The true self-diffusion coefficient $ D_0 $ can be calculated as $ D_0 = (D_x + D_y) / 2 $. - The shear viscosity $ \eta $ can be calculated as $ \eta = k_B T \cdot 8.1711245653 / [3\pi L_z (D_x + D_y - 2D_z)] $, where $ k_B $ is the Boltzmann constant and $ T $ is the temperature. 2. **Applicability**: - This method is applicable to various systems, including liquid water, dimethyl ether, methanol, triethylene glycol, water/methanol mixtures, water/triethylene glycol mixtures, and imidazolium-based ionic liquids. - The viscosity of these systems ranges from a few millipascal-seconds to several hundred millipascal-seconds. ### Main Contributions 1. **Improved Computational Efficiency**: - Reliable data can be effectively calculated using systems with 768 molecules or ion pairs without increasing the system size. - Extending the simulation time is more effective than increasing the system size in improving data reliability. 2. **Simple Formulas Provided**: - The paper provides simple formulas to estimate the required MD simulation time and system size to achieve a target uncertainty level. - By analyzing the impact of system size on statistical uncertainty, it concludes that extending the simulation time is more effective than increasing the system size. ### Conclusion The paper successfully addresses the system size dependence issue of self-diffusion coefficient and shear viscosity in MD simulations through the "OrthoBoXY" method, providing an efficient and accurate approach for calculating these important transport properties.

Computing Accurate True Self-Diffusion Coefficients and Shear Viscosities Using the OrthoBoXY-Approach

Computing Accurate True Self-Diffusion Coefficients and Shear Viscosities Using the OrthoBoXY Approach

OrthoBoXY: A Simple Way to Compute True Self-Diffusion Coefficients from MD Simulations with Periodic Boundary Conditions Without Prior Knowledge of the Viscosity

An OrthoBoXY-Method for Various Alternative Box Geometries

Benchmarking the Diffusion and Viscosity of H-He Mixtures in Warm Dense Matter Regime by Quantum Molecular Dynamics Simulations

Diffusivity of Liquid Mixture Confined in Silicate Pore by Molecular Dynamic Simulation

Self-diffusion and shear viscosity of pure 1-alkanol unary system: molecular dynamics simulation and review of experimental data

Shear stress relaxation and diffusion in simple liquids by molecular dynamics simulations: Analytic expressions and paths to viscosity

Self-diffusion coefficient and viscosity of methane and carbon dioxide via molecular dynamics simulations based on new ab initio-derived force fields

A cluster approach to rationalize shear thinning: Application to polymer solutions and suspension fluids

Assessment of Approaches towards the Relative Permittivity of Mixtures

MD2D: a python module for accurate determination of diffusion coefficient from molecular dynamics

Accurate diffusion coefficients of the excess proton and hydroxide in water via extensive ab initio simulations with different schemes

Efficient single-run implementation of generalized Einstein relation to compute transport coefficients: A binary-based time sampling

Computation of Shear Viscosity by a Consistent Method in Equilibrium Molecular Dynamics Simulations: Applications to 1-Decene Oligomers

Optimal estimates of diffusion coefficients from molecular dynamics simulations

Application of Molecular Dynamics Simulations in Molecular Property Prediction II: Diffusion Coefficient.

A New Single Equation of State to describe the Dynamic Viscosity and Self-Diffusion Coefficient for all Fluid Phases of Water from 200 K to 1800 K based on a New Original Microscopic Model

Accurate Estimation of Transport Coefficients Using Model-free Time Correlation Functions in Equilibrium Simulations

Parametric dependence on shear viscosity of SPC/E water from equilibrium and non-equilibrium molecular dynamics simulations

Parametrization of coarse grained force fields for dynamic property of ethylene glycol oligomers/water binary mixtures