{"id":"64cfc79f-7c39-488f-9b22-3885f4bb236b","arxiv_id":"2412.07663","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"Shear relaxation times of Lennard-Jones, Yukawa, soft-sphere, and hard-sphere fluids collapse onto a quasi-universal curve, with freezing-point values near 0.18 in microscopic reduced units.","lead":"The authors combine published viscosity and elasticity data for four model fluids and find that, when time is measured in microscopic units, the shear relaxation time follows the same curve: it dips at intermediate density and rises again toward freezing. If the pattern holds, it gives a simple route to estimate relaxation times in liquids, dusty plasmas, and colloids from only density and temperature.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Hard-sphere anchor rests on an untested approximate modulus; direct MD stress-autocorrelation data would settle whether the 0.18±0.04 freezing band holds.","rationale":"I agree with the reader that the hard-sphere shear modulus is the most load-bearing fragile input. The central claim is a quantitative band at freezing: τ*M ≈ 0.18 ± 0.04. The HS system is one of four data points, and unlike the LJ, Yukawa, and SS moduli—which derive from well-tested equations of state or explicitly established fits—the HS value rests on a non-standard route (Miller [66] plus Tao-Song-Mason [67]) that the paper itself notes is required because the usual Zwanzig-Mountain expression diverges. The Miller derivation and the TSM derivative are both approximate in ways that are not cross-checked within the paper. If the HS G* were, say, 30% lower, τ*M would rise to about 0.24, pushing the spread of values from 0.14–0.21 to 0.14–0.24, still quasi-universal but with a noticeably wider band; if G* were 30% higher, τ*M would drop to about 0.13, overlapping the LJ value and artificially narrowing the band. Thus the precise claim of a tight band is sensitive to this input. Other possible concerns—such as the soft-sphere minimum appearing as '0.9' in Section II E, which is almost certainly a typo for '0.09' given the freezing value of 0.16 and the stated 'factor of two' increase—are less substantive. The qualitative minimum-then-increase trend is robust across all four systems regardless of the HS modulus details, so the appropriate action is to maintain the CONDITIONAL verdict pending a direct MD test of the HS modulus, not to reject the paper. My recommendation is therefore UNCHANGED: the reader's conditional assessment is appropriate and my concern reinforces the same condition.","tokens_in":13110,"tokens_out":9787,"duration_ms":84641,"concrete_test":"Run event-driven molecular dynamics for a hard-sphere fluid at the freezing packing fraction (φ ≈ 0.494, the value used in the paper) and compute the shear stress autocorrelation function C(t). The instantaneous shear modulus is C(0) in the appropriate reduced units, and the Maxwell relaxation time is τ_M = η/G∞ = [∫₀^∞ C(t) dt] / C(0). Compare the resulting G*∞ and τ*M with the Miller/TSM values used in Table I (G* ≈ 40, τ* ≈ 0.17). If the MD result differs from the Miller/TSM G* by more than about 10%, the hard-sphere anchor moves and the claimed 0.18 ± 0.04 freezing band must be re-evaluated. As a secondary check, vary the Tao-Song-Mason g'(σ) within its known uncertainty and report the induced shift in the HS freezing τ*M.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that reduced Maxwell relaxation times at freezing form a narrow quasi-universal band, τ*M ≈ 0.18 ± 0.04, depends on four anchor values. Three of these are obtained from established equations of state or direct Zwanzig-Mountain integrations: the LJ modulus from the Thol EoS, the Yukawa modulus from published fits consistent with [48], and the SS modulus from the exact inverse-power relation plus an 8th-order virial EoS. The hard-sphere modulus is the exception. Section II D computes G∞ using Miller's derivation [66] with the derivative g'(σ) taken from the Tao-Song-Mason approximation [67]. This is the only modulus in the paper that cannot be checked against a known Zwanzig-Mountain integral because that integral diverges in the HS limit. The HS freezing point is set by this route: with η*≈6.8 and G*≈40 it yields τ*M ≈ 0.17. If the Miller/TSM estimate is inaccurate, the HS anchor shifts and the four-system band either widens or narrows. Since the paper provides no direct stress-autocorrelation data or sensitivity analysis for this input, the load-bearing 'all comparable' statement is not yet independently secured for the one system where the standard derivation fails. The concern is not that the Miller route is obviously wrong, but that its accuracy is untested in the present context, and the HS point is exactly the one most in need of such a test.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript examines the Maxwell shear relaxation time τ_M = η/G∞ for four simple monatomic model fluids (Lennard-Jones, Yukawa, soft-sphere, hard-sphere) using the same reduced units based on the interparticle spacing Δ = ρ^{-1/3} and thermal velocity v_T = (T/m)^{1/2}. Using published equations of state, viscosity fits, and elastic-modulus formulas, the authors compute τ_M^* over the entire fluid range up to freezing. They report a common qualitative trend: τ_M^* first decreases with density, reaches a minimum, and then increases toward the freezing point. They further claim that the reduced relaxation times at both the minima and the freezing point are quasi-universal across the four systems, with freezing values τ_M^* ≈ 0.18 ± 0.04. The paper then derives implications for the transverse-mode k-gap, the Stokes-Einstein relation, and the hierarchy of relaxation times in dense fluids.","tokens_in":13455,"tokens_out":6641,"duration_ms":55294,"significance":"If the quasi-universal band is correct, it provides a robust, unit-independent estimate of the Maxwell relaxation time near freezing for simple monatomic fluids, and it connects the rheological Maxwell time to the k-gap and the vibrational model of transport. The paper's compilation is useful because it expresses previously scattered results in a common normalization and identifies the hard-sphere modulus as a delicate case. However, the central evidence is not fully independent: two of the four systems are built from existing quasi-universal fits, and the hard-sphere branch relies on an unvalidated approximate modulus. The paper would be strengthened by a direct test or uncertainty analysis of the hard-sphere input, and by correcting an evident inconsistency in the soft-sphere minimum. The claimed universality at freezing may well survive these tests, but as it stands the evidence is suggestive rather than conclusive.","major_comments":[{"comment":"The hard-sphere shear modulus is the only input not obtained from a direct Zwanzig-Mountain integration or an established equation of state: it relies on Miller's derivation [66] with the contact derivative g'(σ) from the Tao-Song-Mason approximation [67], because Eq. (4) diverges in the HS limit. Since the HS freezing point (τ_M^* ≈ 0.17) is one of the four anchors of the claimed band τ_M^* ≈ 0.18 ± 0.04, the central claim is sensitive to this approximation. The paper provides no sensitivity analysis and no comparison against direct stress-autocorrelation MD data for G∞ or τ_M in hard spheres. I request such a test, or at least an estimate of the uncertainty of the TSM route, before the quasi-universal band can be considered secure.","section":"II.D (hard-sphere fluid)"},{"comment":"The values of the minima listed for the soft-sphere fluid are inconsistent with the stated trend: the text gives τ_M^* ≈ 0.9 at the minimum, while Table I gives τ_M^* ≈ 0.16 at freezing for the same system. If 0.9 is the true minimum, then τ_M^* decreases from 0.9 to 0.16 on approaching freezing, contradicting the claim that it increases toward freezing and that minima are comparable across systems (LJ 0.07-0.08, Yukawa 0.12, HS 0.15). This appears to be a typo for 0.09, but as written it undermines the summary and must be corrected.","section":"II.E (Summary)"},{"comment":"The quasi-universal collapse for the LJ and Yukawa systems is in large part inherited from the empirical inputs: the LJ viscosity is generated by the modified excess entropy scaling of Ref. [22] (already a freezing-density-scaling collapse) and the Yukawa branch uses the practical viscosity formula of Ref. [44] and the quasi-universal G∞ of Ref. [48]. Thus, for these two systems, the near-constancy of τ_M^* across the dense-fluid regime is not a fully independent test of quasi-universality. The strongest independent evidence comes from the soft-sphere and hard-sphere branches, which use different sources. The paper would be more complete if it stated this explicitly and discussed how much of the spread in Figs. 3-4 reflects the uncertainty of the underlying fits.","section":"II.A and II.B"}],"minor_comments":[{"comment":"The phrase 'The dotes denote the original calculation' contains a typo: 'dotes' should be 'dots'.","section":"Fig. 2 caption"},{"comment":"The sentence 'This is not surprising, because the crossover is considered, and thus there is no \"exact demarcation line\"' is unclear; it should be reworded, for example as 'the crossover is not sharp, and thus there is no exact demarcation line'.","section":"II.E (Summary)"},{"comment":"The text states 'The freezing packing fraction is tabulated in the same work' without giving the numerical value; including the freezing packing fraction used would help the reader reproduce Fig. 5.","section":"II.C (Soft-sphere fluid)"},{"comment":"The table would be easier to interpret if each row also reported the specific state point (e.g., T* for LJ, Γ/Γ_fr for Yukawa, ρ/ρ_fr for SS and HS) at which the tabulated values were evaluated.","section":"Table I"}],"recommendation":"major_revision","confidential_remarks":"This is a well-written paper from a group that has made many contributions to the transport properties of simple fluids. The cross-system comparison is new and of interest, but the paper leans heavily on prior fits from the same group, and the hard-sphere modulus is the least tested input. The editor might consider asking a reviewer with specific expertise in hard-sphere rheology to assess the validity of the Miller/TSM route. The soft-sphere minimum value in Section II.E appears to be a typographical error that should be corrected before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this paper does something simple and useful: it puts the Maxwell relaxation time for LJ, Yukawa, soft-sphere, and hard-sphere fluids on a common reduced scale (Delta/vT) and shows they all follow the same trend, with freezing-point values clustered around tau*M ~ 0.18 +/- 0.04. The new content is the systematic cross-system comparison and the hard-sphere calculation, which nobody had done cleanly because the standard Zwanzig-Mountain modulus diverges in the HS limit. That part is legitimately new, and the qualitative picture--decrease at low density, minimum, then rise toward freezing--holds up across all four systems. The paper is also honest. It explicitly flags the HS modulus issue, uses established equations of state and simulation fits, and does not oversell the scatter. The implications section, especially the k-gap estimate and the time-scale hierarchy, follows naturally and gives the result practical value for people working on dusty plasmas or soft matter. The soft spots are real but not fatal. First, the HS anchor relies on Miller's derivation plus the Tao-Song-Mason approximation for g'(sigma), and there is no direct stress-autocorrelation simulation in the paper to check it. If that modulus is off, the 0.18 +/- 0.04 band loses one of its supports. Second, part of the apparent universality is inherited: the LJ and Yukawa curves use the authors' own prior freezing-density and excess-entropy scaling fits, so the comparison flatters the conclusion. Third, there are no error bars anywhere, and the freezing-point band is essentially a read-off of four numbers. These do not kill the claim--the qualitative behavior is consistent across independent sources--but they do mean the quantitative band should be treated as provisional, especially for HS. I would send this to peer review. It is a well-organized, clearly written paper that gives the community a convenient rule of thumb and identifies exactly where further verification is needed. For the right reader--someone working on viscosity, relaxation, or transport in simple fluids--it is a solid reference. A good referee should push for a direct HS stress-tau simulation, or at least a sensitivity analysis of the Miller/TSM input, before trusting the 0.18 +/- 0.04 value as more than an estimate.","headline":"A useful, honest compilation showing quasi-universal reduced Maxwell relaxation times across four simple fluids, with the hard-sphere anchor as the one genuinely fragile input.","tokens_in":685,"tokens_out":796,"would_cite":true,"duration_ms":17392,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper argues that the Maxwell shear relaxation time, reduced by interparticle spacing and thermal velocity, follows a single quasi-universal curve in simple fluids from dilute gas to freezing, with a common value near 0.18 at freezing.","keywords":["shear relaxation time","Maxwell relaxation time","quasi-universality","simple fluids","freezing density scaling","shear viscosity","instantaneous shear modulus","hard-sphere fluid"],"falsifier":"Run a direct molecular-dynamics simulation of the hard-sphere fluid at packing fractions up to freezing and compute the instantaneous shear modulus from stress fluctuations, a route that stays finite, then form $\\tau_M = \\eta/G_\\infty$ with the Green-Kubo viscosity. If the resulting $\\tau_M^*$ at freezing lands outside $0.18 \\pm 0.04$, the hard-sphere anchor, and with it the quasi-universal band, collapses.","tokens_in":12886,"feed_emoji":"⏱️","tokens_out":11780,"duration_ms":97810,"temperature":0.7,"pith_summary":"The Maxwell shear relaxation time $\\tau_M = \\eta/G_\\infty$, the ratio of shear viscosity to instantaneous shear modulus, spans more than ten orders of magnitude in absolute terms across simple fluids. The paper asks whether this spread is due to the interaction potential or merely to different length and time scales. It computes $\\tau_M$ for Lennard-Jones, Yukawa, soft-sphere, and hard-sphere fluids in common reduced units $\\Delta = \\rho^{-1/3}$ and $v_T = \\sqrt{T/m}$, and finds that $\\tau_M^* = \\tau_M v_T/\\Delta$ behaves the same way in all four: it decreases at low density, reaches a minimum at intermediate density, and then rises toward freezing. At freezing the values cluster around $\\tau_M^* \\simeq 0.18 \\pm 0.04$ despite the enormous absolute spread. The upshot is that near their freezing point, simple monatomic melts share a universal viscoelastic time scale set mainly by density relative to freezing, not by interaction details.","feed_headline":"Near freezing, four simple fluids share one relaxation time","feed_subtitle":"Four simple fluids agree within 30 percent at freezing despite hugely different absolute times.","key_machinery":"The load-bearing identity is $\\tau_M = \\eta/G_\\infty$, evaluated in the system-independent reduced units $\\Delta = \\rho^{-1/3}$ and $v_T = \\sqrt{T/m}$, so that $\\tau_M^* = \\tau_M v_T/\\Delta$. The viscosity coefficients come from established fits and molecular-dynamics data for each fluid; the instantaneous shear modulus $G_\\infty$ comes from the standard high-frequency elastic-modulus formula (a kinetic term plus an integral over the radial distribution function) for the Lennard-Jones, Yukawa, and soft-sphere fluids, and from a finite hard-sphere derivation that bypasses the divergence of that standard formula in the hard-sphere limit, using an approximate contact derivative of the radial distribution function. Normalizing density by the freezing density is what brings the four curves together.","core_discovery":"On the paper's own terms, the central discovery is that the reduced Maxwell relaxation time is a quasi-universal function of density normalized by its freezing value. For the four model fluids, the same qualitative density dependence appears—decrease, minimum, increase—and the numerical values at both the minimum and the freezing point are comparable: $\\tau_M^* \\simeq 0.07$ to $0.15$ at the minima and $\\tau_M^* \\simeq 0.18 \\pm 0.04$ at freezing. The authors read this as evidence that interaction softness or the presence of long-range attraction plays no systematic role over the investigated range. They then use the quasi-universal freezing value to derive further near-freezing regularities: a common reduced cutoff wave number $k_{\\rm gap}^* \\simeq 0.50 \\pm 0.06$ for transverse shear waves, a large separation between the diffusion time and the Maxwell time ($\\tau_D/\\tau_M \\simeq 27$ to $44$), and the time-scale ordering $1/\\Omega_E < \\tau_M \\ll \\tau_D$ consistent with a vibrational picture of transport in dense fluids.","pith_inferences":["If the band survives closer scrutiny, the Maxwell relaxation time near freezing becomes a purely thermodynamic estimate, $\\tau_M \\simeq 0.18\\Delta/v_T$, applicable to complex plasmas and colloidal suspensions where the interaction potential is poorly characterized.","The apparent universality suggests a corresponding-states principle for viscoelasticity: simple fluids at equal $\\rho/\\rho_{\\rm fr}$ have equal reduced $\\tau_M$; testing molecular liquids with anisotropic or bounded potentials would show whether the principle extends beyond monatomic pairwise-additive models.","The hard-sphere result is the least secure anchor; an independent stress-fluctuation simulation would settle whether the band is real physics or an artifact of the approximate contact derivative used for the hard-sphere modulus."],"forward_implications":["Near freezing, any simple monatomic fluid has $\\tau_M^* \\simeq 0.18 \\pm 0.04$; given temperature and number density, this yields an absolute relaxation time estimate without knowing the interaction potential.","The reduced cutoff wave number for transverse collective modes is quasi-universal, $k_{\\rm gap}^* \\simeq 0.50 \\pm 0.06$, so the Maxwell time directly fixes where shear waves begin to propagate.","The large separation $\\tau_D/\\tau_M \\simeq 27$ to $44$ at freezing supports the vibrational model of dense-liquid transport and justifies treating atomic oscillations as temporarily solid-like.","The minimum of $\\tau_M^*$ marks the gas-like to liquid-like dynamical crossover but sits deeper in the dense regime than other crossover indicators such as extrema of transport coefficients.","For the potentials studied, no systematic dependence on interaction softness or long-range attraction appears, so the same quasi-universal behaviour should hold for other simple monatomic fluids."],"supporting_citations":[{"why":"Supplies the high-frequency elastic-modulus formula used for the Lennard-Jones, Yukawa, and soft-sphere fluids.","marker":"[18]"},{"why":"Provides the modified excess entropy scaling used to obtain the Lennard-Jones shear viscosity along isotherms.","marker":"[22]"},{"why":"Equation of state for the Lennard-Jones fluid used to evaluate excess energy and pressure in the modulus calculation.","marker":"[29]"},{"why":"Practical formula for the Yukawa shear viscosity that produces the Yukawa relaxation-time curve.","marker":"[44]"},{"why":"Demonstrates the quasi-universal scaling of the Yukawa transverse sound speed used for the Yukawa shear modulus.","marker":"[48]"},{"why":"Supplies the fitting formula for the soft-sphere shear viscosity.","marker":"[53]"},{"why":"Provides the eighth-order virial expansion for the soft-sphere compressibility and the freezing packing fraction used for normalization.","marker":"[56]"},{"why":"Provides the molecular-dynamics hard-sphere shear viscosity data used directly in the calculation.","marker":"[62]"},{"why":"Supplies the finite hard-sphere derivation of elastic moduli used to bypass the divergent standard formula.","marker":"[66]"},{"why":"Approximation for the derivative of the hard-sphere radial distribution function at contact, required in the hard-sphere modulus calculation.","marker":"[67]"}],"fun_headline_variants":["Quasi-universal shear relaxation time in simple fluids","Four simple fluids, one reduced relaxation time near freezing","Shear relaxation time collapses to universal curve for four fluids","Simple fluids share quasi-universal relaxation time at freezing","Reduced shear relaxation time universal across four model fluids"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The hard-sphere shear modulus near freezing is the fragile input: because the standard formula diverges in the hard-sphere limit, the calculation relies on a specific finite-derivation route and on an approximate value for the derivative of the radial distribution function at contact; if that approximate value is off, the hard-sphere anchor of the universal band shifts.","fun_headline_variants_meta":{"raw":{"variants":["Quasi-universal shear relaxation time in simple fluids","Four simple fluids, one reduced relaxation time near freezing","Shear relaxation time collapses to universal curve for four fluids","Simple fluids share quasi-universal relaxation time at freezing","Reduced shear relaxation time universal across four model fluids"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000687,"raw_usage":{"total_tokens":3083,"prompt_tokens":880,"completion_tokens":2203,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":496,"completion_tokens_details":{"reasoning_tokens":2124}},"tokens_in":496,"tokens_out":2203,"duration_ms":15544,"temperature":1.0,"reasoning_tokens":2124,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T18:38:00.765431+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a direct molecular-dynamics simulation of the hard-sphere fluid at packing fractions up to freezing and compute the instantaneous shear modulus from stress fluctuations, a route that stays finite, then form $\\tau_M = \\eta/G_\\infty$ with the Green-Kubo viscosity. If the resulting $\\tau_M^*$ at freezing lands outside $0.18 \\pm 0.04$, the hard-sphere anchor, and with it the quasi-universal band, collapses.","supporting_citations":[{"cited_title":"High-frequency elastic moduli of simple fluids,","cited_arxiv_id":null,"evidence_quote":"Supplies the high-frequency elastic-modulus formula used for the Lennard-Jones, Yukawa, and soft-sphere fluids."},{"cited_title":"Modified entropy scaling of the transport properties of the Lennard-Jones fluid,","cited_arxiv_id":null,"evidence_quote":"Provides the modified excess entropy scaling used to obtain the Lennard-Jones shear viscosity along isotherms."},{"cited_title":"Universal scaling of transverse sound speed and its iso- morphic property in Yukawa fluids,","cited_arxiv_id":null,"evidence_quote":"Demonstrates the quasi-universal scaling of the Yukawa transverse sound speed used for the Yukawa shear modulus."},{"cited_title":"Probing the link between residual entropy and viscosity of molecular fluids and model potentials,","cited_arxiv_id":null,"evidence_quote":"Supplies the fitting formula for the soft-sphere shear viscosity."},{"cited_title":"Thermo- dynamic properties and entropy scaling law for diffusivity in soft spheres,","cited_arxiv_id":null,"evidence_quote":"Provides the eighth-order virial expansion for the soft-sphere compressibility and the freezing packing fraction used for normalization."},{"cited_title":"Elastic moduli of a fluid of rigid spheres,","cited_arxiv_id":null,"evidence_quote":"Supplies the finite hard-sphere derivation of elastic moduli used to bypass the divergent standard formula."},{"cited_title":"Derivative of the hard-sphere radial distribution function at contact,","cited_arxiv_id":null,"evidence_quote":"Approximation for the derivative of the hard-sphere radial distribution function at contact, required in the hard-sphere modulus calculation."}],"review_version":1}