Search papers, labs, and topics across Lattice.
This study computes shear viscosities for 15 neat molecular liquids using both the OrthoBoXY approach and the Green-Kubo method, revealing strong agreement between the two. The authors identify optimal strategies for using the OrthoBoXY method, demonstrating that smaller systems (as few as 250 molecules) can yield accurate viscosity results without significant finite size effects. Notably, they propose scaling the block-length parameter in the OrthoBoXY method to reduce computational costs by up to 24 times for highly viscous systems, while maintaining accuracy.
Small molecular systems can deliver accurate shear viscosity results at a fraction of the computational cost, transforming simulation strategies in fluid dynamics.
We calculated shear viscosities of 15 neat molecular liquids from equilibrium molecular dynamics (MD) simulations using the OrthoBoXY approach and compare them to viscosities calculated via the Green-Kubo method. Data from both methods agree very well. Here, we show how to avoid pitfalls while computing the OrthoBoXY-data to obtain optimal results. From simulations of multiple system sizes, we could verify that the viscosity of molecular liquids is not influenced by finite size effects down to systems as small as 250 molecules. Moreover, we demonstrate that also the standard error of the viscosity is nearly independent of the system size. This is shown to be a consequence of a compensation effect of an increasing accuracy of the self-diffusion coefficients with increasing systems-size and the system-size dependent weighting according to the OrthoBoXY-equation. As a consequence, we suggest that it is preferable to run simulations of smaller systems with longer simulation times rather than larger systems with shorter simulation runs. In addition, we discuss a refinement of the recently introduced"recipe"for OrthoBoXY simulations block-lengths $\tau_\mathrm{block}$. Based on data from simulations with varying run-lengths, we suggest the following modification: for highly viscous systems, the value of $\tau_\mathrm{block}$ might safely be scaled by a factor of $1/8$, significantly reducing the computational resources needed. For less viscous systems, the value of $\tau_\mathrm{block}$ might safely be scaled by a factor of $1/4$. For systems with high fluidity, the value of $\tau_\mathrm{block}$ should not be scaled down in order to achieve reliable results. When using a smaller system size of 250 molecules, these refinements are leading up to a 24-fold reduction in computational cost compared to the previous recommended set-up without sacrificing numerical accuracy.
