{"id":"dd59162f-f090-4ac8-bafe-3a134e0aae68","arxiv_id":"2608.05845","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"OrthoBoXY and Green-Kubo viscosities agree for 15 molecular liquids, and smaller systems with shorter blocks can cut the cost of OrthoBoXY viscosity runs by up to 24-fold.","lead":"This paper tests a cheaper molecular simulation route to shear viscosity, the OrthoBoXY method, against the standard Green-Kubo method for 15 liquids and finds overall agreement. It also proposes shorter simulation blocks and smaller systems, cutting computational cost by up to 24-fold without losing accuracy.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Claimed 'very good agreement' is not supported at the high-viscosity end: glycerol, 3-methylphenol, and acetamide deviate from Green-Kubo by 33–75%, and these are the same systems used to justify the 1/8 block-length scaling.","rationale":"The central claim in the abstract and Section IV is that OrthoBoXY and Green-Kubo 'agree very well' for 15 liquids spanning three orders of magnitude, and that block lengths can be shortened by up to 1/8 for viscous systems. That claim has two load-bearing parts: (i) the method is validated across the viscosity range, and (ii) the shortened blocks give the same accuracy. My concern targets part (i) at the high-viscosity end. The numbers in Table II are not consistent with the qualitative claim. The three largest deviations are all at high viscosity and are large relative to the stated uncertainties. The paper acknowledges but does not resolve this. Moreover, the block-length recommendation in Section IV.D is based on systems that include glycerol and acetamide, the very systems showing the largest disagreement with Green-Kubo. Demonstrating that shortened blocks reproduce the OrthoBoXY values of those systems only shows internal consistency of the method, not agreement with the reference method that anchors the validation. The reader's concern about generalizing from four systems is real, but it is secondary: even for the four tested systems, the accuracy anchor is compromised. A same-trajectory Green-Kubo calculation would separate a true method bias from a protocol mismatch with the external reference data. This does not change the overall CONDITIONAL verdict; it sharpens the condition.","tokens_in":11022,"tokens_out":5544,"duration_ms":51960,"concrete_test":"Compute Green-Kubo viscosities from the authors' own stored trajectories for glycerol, 3-methylphenol, and acetamide, using the same 40-block decomposition and the same OPLS-AA force field (e.g., integrating the off-diagonal stress autocorrelation function with the same run lengths and thermostat settings). If the in-house Green-Kubo values agree with the OrthoBoXY values within the reported standard errors, the discrepancies in Table II come from external reference conditions and the agreement claim is restored; if the discrepancy exceeds twice the combined standard error, the central 'very good agreement' claim fails at the viscous end and the 1/8 block-scaling recommendation loses its accuracy benchmark.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section IV.A and Table II show the three most viscous systems deviate far beyond the reported 1σ error bars: glycerol 230±24 vs 160±2 (OrthoBoXY 44% high), 3-methylphenol 70±8 vs 40±2 (75% high), and acetamide 10.0±0.8 vs 15±1 (33% low). The abstract's 'agree very well' overstates this. These deviations are not small compared to the three-orders-of-magnitude claim; they occur at the upper end of the range and are the systems where the paper recommends the 1/8 block-length reduction (Section IV.D). If OrthoBoXY is systematically biased in this regime, then the block-length validation—which compares shortened runs only against OrthoBoXY's own longer-block values—does not establish numerical accuracy in the sense that matters for the central claim. The paper's own caveat that 'it can not be said which method is responsible' (Section IV.A) makes the agreement claim conditional rather than established. The cross-check against Green-Kubo is not a same-trajectory comparison: the reference values are taken from Smith and Sega's independent simulations with possibly different thermostatting and run lengths, so the discrepancy could be a simulation-protocol artifact rather than a method error. This is the load-bearing weak point: the headline validation rests on a benchmark that is not met where the method is most needed.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":11333,"tokens_out":3299,"duration_ms":32472,"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":[{"comment":"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.","section":"IV.A, Table II"},{"comment":"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.","section":"IV.D, Fig. 5"},{"comment":"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.","section":"IV.C, Fig. 4"}],"minor_comments":[{"comment":"There is a typo in the variant headings: 'V ariant' should be 'Variant'.","section":"IV.B"},{"comment":"In the Conclusions, 'less viscous systemts' should read 'less viscous systems'.","section":"V"},{"comment":"The word 'indepdent' in 'not completely indepdent from one another' should be 'independent'.","section":"IV.C"},{"comment":"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.","section":"Fig. 4"},{"comment":"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.","section":"III"},{"comment":"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.","section":"Table I and Table II"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the scope of a computational methods journal and the underlying simulations appear competently executed. My main concern is that the headline 'agree very well' claim is not supported at the high-viscosity end, and the block-length scaling recommendation is drawn from a small sample with one clear failure. I believe these issues can be fixed within the manuscript's scope by adding same-protocol Green-Kubo comparisons for the discrepant systems, quantifying the error-scaling fits, and qualifying the transferability of the block-length recipe. I therefore recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a solid methods paper, and the abstract is both too strong and too weak — too strong about 'very good agreement', too weak about what the new analysis actually gives you. The genuinely new pieces are the 15-liquid OrthoBoXY vs. Green-Kubo comparison against Smith and Sega's independent data, the three averaging-variant analysis that identifies a real pitfall (variant A produces skewed viscosity distributions and occasional negative values when the D0 and Dz distributions overlap), and the error-compensation argument for the near-independence of the standard error on system size. Topologies, inputs, and the analysis code are archived, so the results are checkable. That counts for a lot.\n\nThe comparison mostly holds: 12 of 15 systems agree within or near the stated error bars across three orders of magnitude. The exception is the high-viscosity corner, and it is not small: glycerol 44% high (230±24 vs 160±2), 3-methylphenol 75% high (70±8 vs 40±2), acetamide 33% low (10.0±0.8 vs 15±1). The authors honestly say it cannot be determined which method is at fault; Green-Kubo is known to struggle above ~20 mPa s, and for 3-methylphenol both methods sit an order of magnitude above the experimental value (6.6 mPa s), so the OPLS model is in trouble there too. The reference values also come from independent simulations with their own thermostatting and run lengths, so protocol differences could be part of the story. Any of these explanations is plausible; the paper just does not pin it down.\n\nThe block-length recommendation is the second soft spot. The 1/8 scaling for 'highly viscous' systems is effectively supported by glycerol alone — the only tested system in the tau_block > 100 ns class — and the shortened runs are compared only to OrthoBoXY's own longer-block values. That establishes internal consistency, which is the right check for a scaling recommendation, but it lands in the same regime where external agreement with Green-Kubo is worst. The 24-fold cost claim also stitches together two separately tested reductions rather than one tested protocol. Both points are fixable with wording, but they should be fixed.\n\nThe finite-size error-independence argument rests on three systems with large scatter, and the N^-1/3 cancellation is asserted more cleanly than the data resolve. Minor overreach, not a flaw in the practical conclusion.\n\nVerdict: worth a serious referee. Reproducible artifacts, a useful recipe, an honest limitations section, and a headline that needs softening. I would send it to review expecting moderate revision. It is a good reading-group session on what internal consistency can and cannot validate.","headline":"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.","tokens_in":11876,"tokens_out":9025,"would_cite":true,"duration_ms":76387,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82B80","82D15"],"pacs":["02.70.Ns","66.20.-d"],"model":"deepseek-v4-flash","headline":"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.","keywords":["shear viscosity","molecular dynamics","OrthoBoXY","Green-Kubo method","self-diffusion coefficient","finite-size effects","block length recipe","OPLS force field"],"falsifier":"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.","tokens_in":10795,"feed_emoji":"🧪","tokens_out":7740,"duration_ms":65169,"temperature":0.7,"pith_summary":"This paper establishes that the OrthoBoXY approach, which extracts the shear viscosity from the direction-dependent finite-size correction to self-diffusion in a single molecular dynamics run, produces viscosities that agree with the standard Green-Kubo method for 15 neat molecular liquids covering three orders of magnitude in viscosity. It shows that neither the viscosity nor its statistical error depends on the simulated system size down to 250 molecules, because the system-size weighting in the OrthoBoXY equation cancels the size dependence of the diffusion-coefficient errors. It then argues that the previously recommended block length can be shortened by a factor of 1/8 for highly viscous systems and by 1/4 for moderately viscous ones, while fluid systems should keep the full block length. Combined with a 250-molecule system, this reduces the computational cost by up to 24-fold without sacrificing accuracy. A sympathetic reader would care because these recommendations turn a costly transport-coefficient calculation into a cheaper routine one with a clear accuracy target.","feed_headline":"Viscosity from one MD run matches Green-Kubo across 15 liquids","feed_subtitle":"Smaller boxes and shorter blocks cut simulation cost up to 24-fold without losing accuracy.","key_machinery":"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.","core_discovery":"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 $\tau_{\\mathrm{block}}/8$ for $\tau_{\\mathrm{block}}>100$ ns and $\tau_{\\mathrm{block}}/4$ for $1<\\tau_{\\mathrm{block}}<100$ ns, while keeping the full length for $\tau_{\\mathrm{block}}<1$ ns.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the Green-Kubo relation that is the paper's reference method for computing shear viscosity.","marker":"3,4"},{"why":"Introduces the OrthoBoXY geometry and the direction-dependent Yeh-Hummer correction that the method is built on.","marker":"6"},{"why":"Presents the block-length recipe that this paper refines and shortens.","marker":"8"},{"why":"Provides the Green-Kubo viscosities, OPLS densities, and the set of 146 liquids from which the 15 test systems are drawn.","marker":"9"},{"why":"Gives the Yeh-Hummer finite-size formula for self-diffusion, which OrthoBoXY rearranges to obtain viscosity.","marker":"12"},{"why":"Establishes the idea of extracting viscosity from the system-size dependence of diffusion coefficients.","marker":"22"},{"why":"Documents where the Green-Kubo method is considered reliable, used when interpreting deviations at high viscosity.","marker":"5"}],"fun_headline_variants":["OrthoBoXY matches Green-Kubo, cuts viscosity cost 24x","Small MD boxes, short blocks: same viscosity, 24x less cost","Viscosity via OrthoBoXY: 24x cheaper, matches Green-Kubo","OrthoBoXY gives Green-Kubo viscosity at 1/24th the cost"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["OrthoBoXY matches Green-Kubo, cuts viscosity cost 24x","Small MD boxes, short blocks: same viscosity, 24x less cost","Viscosity via OrthoBoXY: 24x cheaper, matches Green-Kubo","OrthoBoXY gives Green-Kubo viscosity at 1/24th the cost"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001707,"raw_usage":{"total_tokens":6866,"prompt_tokens":1160,"completion_tokens":5706,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":776,"completion_tokens_details":{"reasoning_tokens":5618}},"tokens_in":776,"tokens_out":5706,"duration_ms":38296,"temperature":1.0,"reasoning_tokens":5618,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T22:25:39.835302+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the OrthoBoXY geometry and the direction-dependent Yeh-Hummer correction that the method is built on."},{"cited_title":"Botan and V","cited_arxiv_id":null,"evidence_quote":"Presents the block-length recipe that this paper refines and shortens."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Green-Kubo viscosities, OPLS densities, and the set of 146 liquids from which the 15 test systems are drawn."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents where the Green-Kubo method is considered reliable, used when interpreting deviations at high viscosity."}],"review_version":1}