REVIEW 3 major objections 6 minor 46 references
Computing Shear Viscosities from Molecular Dynamics Simulation: Comparing the OrthoBoXY Approach with the Green-Kubo Method
T0 review · 3 major / 6 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read Shear viscosities from the OrthoBoXY approach agree with Green-Kubo results for 15 molecular liquids, and the block length can be shortened by factors of 1/8 to 1/4 for viscous systems, cutting computational cost by up to 24-fold.
desk verdict Genuinely useful OrthoBoXY validation with reproducible data and a practical recipe refinement, but the 'very good agreement' headline overstates the high-viscosity end and the 1/8 block-length scaling rests on a single viscous system. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the OrthoBoXY geometry: an orthorhombic simulation box with $L_z/L_x = L_z/L_y = 2.7933\ldots$, for which the direction-dependent Yeh-Hummer correction to self-diffusion vanishes in the $x$ and $y$ directions ($\zeta_{xx}=\zeta_{yy}=0$). This makes $D_0 = (D_{PBC,xx}+D_{PBC,yy})/2$ a direct estimate of the true self-diffusion coefficient, while the residual $z$-direction correction gives the viscosity through $\eta = k_B T \zeta_{zz} / (6\pi L_z(D_0 - D_{PBC,zz}))$ with $\zeta_{zz}=8.1711\ldots$. The statistical argument is carried by the error-propagation formula $\hat{\sigma}_\eta/\eta = \sqrt{\hat{\sigma}_{D_0}^2+\hat{\sigma}_{D_{PBC,zz}}^2}/|D_0-D_{PBC,zz}|$, whose denominator scales as $N^{-1/3}$ and cancels the $N^{-1/3}$ scaling of the numerator, making the viscosity error nearly independent of system size. The block-length recipe $\tau_{\mathrm{block}} = (8.5s)^2/(6D)$, where $s=(V/N)^{1/3}$, is the target of the paper's cost-saving refinement.
What would settle it
Run OrthoBoXY at $\tau_{\mathrm{block}}/8$ for a viscous liquid outside the four test systems, for instance an ionic liquid or polyol with $\tau_{\mathrm{block}}>100$ ns, and compare with the full-length recipe and with Green-Kubo; a disagreement beyond the standard error would disprove the universal scaling recommendation.
Extended reading notes
Core claim
The central claim is that the OrthoBoXY method is a valid and cheaper alternative to Green-Kubo for computing shear viscosities of molecular liquids. For 15 OPLS-modeled liquids with viscosities from about 0.1 to 230 mPa s, the OrthoBoXY results agree with Green-Kubo values within the reported errors, with the largest deviations appearing for the most viscous systems where the two methods diverge noticeably and it is unclear which method is at fault. The paper also claims that finite-size effects are absent down to 250 molecules, and that the standard error of the viscosity is nearly independent of system size because the $N^{-1/3}$ growth of the weighting factor $|D_0-D_{PBC,zz}|$ in the error formula compensates the $N^{-1/3}$ scaling of the self-diffusion coefficient errors. On this basis it recommends running small systems with long trajectories, and shortens the block-length recipe to $ au_{\mathrm{block}}/8$ for $ au_{\mathrm{block}}>100$ ns and $ au_{\mathrm{block}}/4$ for $1<\tau_{\mathrm{block}}<100$ ns, while keeping the full length for $ au_{\mathrm{block}}<1$ ns.
Load-bearing premise
The shortening of $\tau_{\mathrm{block}}$ by factors of 1/8 and 1/4 is validated on only four liquids, with formaldehyde already failing to give a constant viscosity at short blocks, and the paper assumes these results transfer to other molecular liquids in the same viscosity ranges.
Editorial extensions
If this is right
- OrthoBoXY can serve as a drop-in alternative to Green-Kubo for viscosity screens, giving comparable values from a single equilibrium run.
- Simulations with 250 molecules can replace 1000-molecule boxes without introducing finite-size bias, so computing time can be invested in longer trajectories instead.
- The standard error of the viscosity being nearly independent of system size means that increasing box size is not an effective way to improve precision; longer runs are.
- The refined block-length rules let users cut $\tau_{\mathrm{block}}$ by 1/8 for viscous liquids and 1/4 for moderately viscous ones, which with the smaller box yields up to 24-fold lower cost.
- The averaging variant that averages both $D_0$ and $D_{PBC,zz}$ first avoids the skewed distributions and occasional negative viscosities produced by block-wise averaging.
Reading between the lines
- The 1/8 and 1/4 scaling factors are drawn from only four liquids and should be re-validated on other viscous and moderately viscous systems before being adopted as a universal recipe.
- If the compensation effect holds generally, the optimal allocation of a fixed CPU budget between block length and number of blocks may be derivable from the reported error formula, allowing simulation planning without pilot runs.
- Because OrthoBoXY yields both the true self-diffusion coefficient and the viscosity from one run, combining the shortened blocks with 250-molecule boxes would make large-scale viscosity screening of force-field libraries much cheaper than Green-Kubo-based approaches.
- For highly fluid systems the paper advises against shortening; a direct test of where the breakdown begins, for example by scanning $\tau_{\mathrm{block}}$ continuously for a fluid like acetone, would sharpen the boundary between safe and unsafe scaling.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports equilibrium molecular dynamics simulations of 15 OPLS-AA molecular liquids in the OrthoBoXY geometry and compares the resulting shear viscosities with Green-Kubo values taken from Smith and Sega. It also examines three averaging variants for the OrthoBoXY equation, studies system-size effects (250–1000 molecules) for three liquids, and investigates block-length scaling for four liquids. On this basis it recommends reducing the OrthoBoXY block length by factors of 1/8 for viscous and 1/4 for moderately viscous systems, claiming up to a 24-fold reduction in computational cost without loss of accuracy. The paper includes a useful practical discussion of numerical pitfalls and provides data availability links for topologies and input files.
Significance. If the central agreement claim held uniformly, the paper would be a valuable practical validation of the OrthoBoXY method and a useful resource for practitioners, especially the finite-size analysis and the averaging-variant discussion. The compensation argument connecting self-diffusion error scaling to the system-size independence of the viscosity error is interesting and testable. However, the largest deviations between OrthoBoXY and Green-Kubo occur for precisely the systems where the recommended block-length reduction is applied, and the block-length recommendation rests on a small and partly non-constant sample. These issues limit the significance of the cost-saving claim as it stands, but they are addressable with additional analysis or with appropriately qualified conclusions.
major comments (3)
- [IV.A, Table II] The abstract and Section IV.A state that OrthoBoXY and Green-Kubo viscosities 'agree very well', but this is not supported for the three most viscous systems. In Table II, glycerol is 230±24 vs 160±2 mPa s (OrthoBoXY 44% higher), 3-methylphenol is 70±8 vs 40±2 mPa s (75% higher), and acetamide is 10.0±0.8 vs 15±1 mPa s (33% lower). Each deviation is several times the reported standard error. These are also the systems used to justify the 1/8 block-length scaling in Section IV.D, so the validation is weakest exactly where the recommended cost-saving measure is applied. The authors should either perform Green-Kubo calculations under the same protocol (thermostat, run length, system size) for these systems, or substantially temper the agreement claim and discuss possible systematic bias in this regime.
- [IV.D, Fig. 5] The recommendation to scale tau_block by 1/8 for viscous and 1/4 for moderately viscous systems is based on only four liquids, and formaldehyde, the most fluid of the four, does not show constant viscosities at shortened blocks. It is therefore not established that the scaling transfers to other liquids in the stated tau_block ranges. The paper should provide a physically motivated criterion (for example, a check that the mean-squared displacement is in the diffusive regime) or restrict the recommendation to systems similar to those tested.
- [IV.C, Fig. 4] The claim that the standard error of eta is 'nearly independent' of system size and that the N^{-1/3} weighting 'exactly cancels' the error scaling is stronger than the data support. Only three systems are shown, with large scatter, and no quantitative fit statistics are reported for the N^{-1/3} trends shown in Figs. 4(a) and 4(b). Please report fit parameters and uncertainties, and soften the cancellation statement to 'approximately compensates' unless more data are provided.
minor comments (6)
- [IV.B] There is a typo in the variant headings: 'V ariant' should be 'Variant'.
- [V] In the Conclusions, 'less viscous systemts' should read 'less viscous systems'.
- [IV.C] The word 'indepdent' in 'not completely indepdent from one another' should be 'independent'.
- [Fig. 4] The legend entries 'N^{-1/3} Fit' are not accompanied by visible fit curves or fit parameters; please clarify what is plotted or cite fitted values in the text.
- [III] The choice of the MSD linear-fit window from 0.03×tau_block to 0.4×tau_block is stated without justification; a brief explanation or citation would help readers assess the sensitivity of the reported diffusion coefficients and viscosities.
- [Table I and Table II] The Green-Kubo reference values are attributed to Smith and Sega Ref. 9 in Table II but to Refs. 9 and 10 elsewhere; please harmonize the citation usage.
Circularity Check
No material circularity: OrthoBoXY viscosities are computed from an Ewald-derived geometric relation and benchmarked against independent Green-Kubo data; the block-length scaling is an empirical recommendation, not a fitted target.
full rationale
The central comparison is not circular. Equation (5) computes eta from the simulated difference D0 - DPBC,zz and the geometric constant zeta_zz = 8.1711..., which is obtained by Ewald summation and does not depend on the target viscosities. The Green-Kubo reference values are taken from Smith and Sega (Refs. 9-10), an external data set, so the agreement claim is not constructed from the inputs. The self-citations (Refs. 6-8, 21) supply the OrthoBoXY geometry, the zeta values, and the 8.5-times-diameter block-length heuristic; the heuristic is explicitly revisited and modified, and the zeta values are independently computable constants, so these citations are not load-bearing in a circular sense. The system-size error-compensation argument is a derived prediction from Eq. (13) combined with empirically observed error scalings, and it is tested in Fig. 4 rather than assumed. The block-length scaling factors (1/8, 1/4) are empirical fits to four systems, with formaldehyde explicitly failing the short-block test; this is a generalizability and robustness concern about the cost-saving recommendation, not a circular derivation. Large deviations at the high-viscosity end (glycerol, 3-methylphenol, acetamide in Table II) weaken the headline agreement claim but are benchmark-failure or correctness issues, not evidence that a prediction reduces to its own input. Overall, no step in the claimed derivation chain is equivalent to its inputs by construction.
Assumptions & free parameters
free parameters (3)
- Block-length scaling factors (1/8, 1/4, 1) =
1/8 for tau_block > 100 ns, 1/4 for 1 ns < tau_block < 100 ns, 1 for tau_block < 1 ns
- 8.5 displacement factor in tau_block recipe =
8.5 times the linear diameter s
- MSD linear-fit window =
0.03 tau_block to 0.4 tau_block
assumptions (5)
- domain assumption The Yeh-Hummer equation (Eq. 3) accurately describes the periodic-boundary-condition correction to self-diffusion in orthorhombic simulation boxes.
- standard math For the OrthoBoXY box ratio L_z/L_x = L_z/L_y = 2.7933, zeta_xx = zeta_yy = 0 and zeta_zz = 8.1711.
- domain assumption The standard error of the self-diffusion coefficient for molecular liquids scales as N^-1/3.
- domain assumption The OPLS-AA force field gives realistic viscosities for the 15 investigated liquids.
- domain assumption Blocks of length tau_block are statistically independent.
Cite this review
Pith. "Pith review of Computing Shear Viscosities from Molecular Dynamics Simulation: Comparing the OrthoBoXY Approach with the Green-Kubo Method." pith.science (2026). https://pith.science/paper/43NDTLFE
@misc{pith2026260805845,
author = {Pith},
title = {Pith review of: Computing Shear Viscosities from Molecular Dynamics Simulation: Comparing the OrthoBoXY Approach with the Green-Kubo Method},
year = {2026},
howpublished = {\url{https://pith.science/paper/43NDTLFE}},
note = {Machine review of arXiv:2608.05845}
}
abstract
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.
Figures
Reference graph
Works this paper leans on
-
[1]
D. Paschek , title =. 2026 , publisher =. doi:10.5281/zenodo.21807092 , url =
-
[2]
I. D. Smith and M. Sega , title =. J. Chem. Phys. , year =
- [3]
-
[4]
M. S. Green , title =. J. Chem. Phys. , year =
-
[5]
A. T. Celebi and S. H. Jamali and A. Bardow and T. J. H. Vlugt and O. A. Moultos , title =. Mol. Sim. , year =
-
[6]
S. H. Jamali and A. Bardow and T. J. H. Vlugt and O. A. Moultos , title =. J. Chem. Theory Comput. , year =
-
[7]
Transient hydrodynamic finite-size effects in simulations under periodic boundary conditions , author =. Phys. Rev. E , volume =. 2017 , month =
work page 2017
- [8]
Show all 46 references
-
[9]
M. V. Divergent Diffusion Coefficients in Simulations of Fluids and Lipid Membranes , journal =. 2016 , volume =
2016
-
[10]
D. S. Viswanath and T. K. Ghosh and D. H. L. Prasad and N. V. K. Dutt and K. Y. Rani , title =. 2007 , address =
2007
-
[11]
L. D. Landau and E. M. Lifshitz , title =. 1959 , address =
1959
-
[12]
M. P. Allen and D. J. Tildesley , title =. 1987 , address =
1987
-
[13]
Busch and D
J. Busch and D. Paschek , title =. Phys. Chem. Chem. Phys. , year =
-
[14]
Insights from Virtual Chemistry:
Smith, Imogen Daisy and Sega, Marcello , year = 2025, month = mar, journal =. Insights from Virtual Chemistry:. doi:10.1063/5.0251585 , urldate =
2025 doi
-
[15]
and Messerly, Richard A
Maginn, Edward J. and Messerly, Richard A. and Carlson, Daniel J. and Roe, Daniel R. and Elliot, J. Richard , year = 2018, journal =. Best. doi:10.33011/livecoms.1.1.6324 , urldate =
2018 doi
-
[16]
doi:10.1021/acs.jpcb.3c04492 , urldate =
Busch, Johanna and Paschek, Dietmar , year = 2023, month = sep, journal =. doi:10.1021/acs.jpcb.3c04492 , urldate =
2023 doi
-
[17]
Molecular Dynamics Simulation of a Polymer Chain in Solution , author =. J. Chem. Phys. , volume =. doi:10.1063/1.465445 , urldate =
-
[18]
Yeh, In-Chul and Hummer, Gerhard , year = 2004, month = oct, journal =. System-. doi:10.1021/jp0477147 , urldate =
2004 doi
-
[19]
Kikugawa, Gota and Ando, Shotaro and Suzuki, Jo and Naruke, Yoichi and Nakano, Takeo and Ohara, Taku , title =. J. Chem. Phys. , volume =. 2015 , doi =
2015
-
[20]
Hydrodynamic Consideration of the Finite Size Effect on the Self-Diffusion Coefficient in a Periodic Rectangular Parallelepiped System , author =. J. Chem. Phys. , volume =. doi:10.1063/1.4926841 , urldate =
-
[21]
Jamali, Seyed Hossein and Hartkamp, Remco and Bardas, Christos and S. Shear. J. Chem. Theory Comput. , volume =. doi:10.1021/acs.jctc.8b00625 , urldate =
-
[22]
and Maxwell, David S
Jorgensen, William L. and Maxwell, David S. and. Development and. J. Am. Chem. Soc. , volume =. doi:10.1021/ja9621760 , urldate =
-
[23]
Caleman, Carl and. Force. J. Chem. Theory Comput. , volume =. doi:10.1021/ct200731v , urldate =
-
[24]
Mart. J. Comput. Chem. , volume =. doi:10.1002/jcc.21224 , urldate =
-
[25]
SoftwareX , volume =
Abraham, Mark James and Murtola, Teemu and Schulz, Roland and P. SoftwareX , volume =. doi:10.1016/j.softx.2015.06.001 , urldate =
2015 doi
-
[26]
and Abraham, M
Lindahl, E. and Abraham, M. J. and Hess, B. and. doi:10.5281/zenodo.3685922 , urldate =
-
[27]
A Molecular Dynamics Method for Simulations in the Canonical Ensemble , author =. Mol. Phys. , volume =. doi:10.1080/00268978400101201 , urldate =
-
[28]
, year = 1985, month = mar, journal =
Hoover, William G. , year = 1985, month = mar, journal =. Canonical Dynamics:. doi:10.1103/PhysRevA.31.1695 , urldate =
1985 doi
-
[29]
Canonical Sampling through Velocity Rescaling , author =. J. Chem. Phys. , volume =. doi:10.1063/1.2408420 , urldate =
-
[30]
J. Comput. Chem. , volume =. doi:10.1002/jcc.21787 , urldate =
-
[31]
and Linke, Max and Barnoud, Jonathan and Reddy, Tyler J
Gowers, Richard J. and Linke, Max and Barnoud, Jonathan and Reddy, Tyler J. E. and Melo, Manuel N. and Seyler, Sean L. and Doma. Proceedings of the 15th Python in Science Conference , pages =. doi:10.25080/Majora-629e541a-00e , urldate =
-
[32]
2003 , publisher =
Yaws' Handbook of Thermodynamic and Physical Properties of Chemical Compounds , author =. 2003 , publisher =
2003
-
[33]
2016 , publisher =
CRC Handbook of Chemistry and Physics , editor =. 2016 , publisher =
2016
-
[34]
Prediction of Liquid Viscosity for Organic Compounds by a Quantitative Structure--Property Relationship , author =. J. Phys. Org. Chem. , volume =. doi:10.1002/(SICI)1099-1395(200001)13:1<80::AID-POC179>3.0.CO;2-8 , urldate =
-
[35]
Viscosity of
Wu, Jiangtao and Liu, Zhigang and Bi, Shengshan and Meng, Xianyang , year = 2003, month = mar, journal =. Viscosity of. doi:10.1021/je0256232 , urldate =
2003 doi
-
[36]
The Viscosity of Glycerol , author =. J. Chem. Thermodyn. , volume =. doi:10.1016/j.jct.2017.05.042 , urldate =
2017 doi
-
[37]
Fluid Ph
Viscosities and Densities for Binary Mixtures of Cresols , author =. Fluid Ph. Equilib. , volume =. doi:10.1016/S0378-3812(03)00076-1 , urldate =
-
[38]
Density,
Yasmin, Maimoona and Gupta, Manisha , year = 2011, month = sep, journal =. Density,. doi:10.1007/s10953-011-9731-1 , urldate =
2011 doi
-
[39]
Density and
Yang, Changsheng and Lai, Hexi and Liu, Zhanguang and Ma, Peisheng , year = 2006, month = jul, journal =. Density and. doi:10.1021/je0600808 , urldate =
2006 doi
-
[40]
Guinda, L. M. and Santaf. Viscosity Measurements on Aniline, p-Toludine and p-Anisidine + Ph\'enol in the Temperature Range of 303.15-343.15. J. Chim. Phys. , volume =. doi:10.1051/jcp/1986830631 , urldate =
-
[41]
System-size corrections for self-diffusion coefficients calculated from molecular dynamics simulations:. J. Chem. Phys. , author =. 2016 , pages =. doi:10.1063/1.4960776 , number =
2016 doi
-
[42]
Nature of Intrinsic Uncertainties in Equilibrium Molecular Dynamics Estimation of Shear Viscosity for Simple and Complex Fluids , author =. J. Chem. Phys. , volume =. doi:10.1063/1.5035119 , urldate =
-
[43]
and Vlugt, Thijs J
Hulikal Chakrapani, Thejas and Hajibeygi, Hadi and Moultos, Othonas A. and Vlugt, Thijs J. H. , year = 2025, month = nov, journal =. Impact of Finite-Size Effects on Computed Transport Properties: A Molecular Dynamics Study of Dilute Systems , shorttitle =. doi:10.1080/0026897...
2025
-
[44]
, year = 2015, month = aug, journal =
Zhang, Yong and Otani, Akihito and Maginn, Edward J. , year = 2015, month = aug, journal =. Reliable. doi:10.1021/acs.jctc.5b00351 , urldate =
2015 doi
-
[45]
, year = 1969, month = jun, journal =
Zwanzig, Robert and Ailawadi, Narinder K. , year = 1969, month = jun, journal =. Statistical. doi:10.1103/PhysRev.182.280 , urldate =
1969 doi
-
[46]
Computing
Busch, Johanna and Paschek, Dietmar , year = 2024, month = feb, journal =. Computing. doi:10.1021/acs.jpcb.3c07540 , urldate =
2024 doi
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.