{"id":"abb86376-0ca4-4c07-b459-67307d9dc603","arxiv_id":"2607.23127","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"In distinguishable helium-4, cooling through a crossover produces a gas-like liquid (HQDL) with ultralow viscosity and superdiffusion, driven by nuclear quantum delocalization alone.","lead":"Simulations of helium-4 atoms treated as distinguishable—no Bose exchange, no superfluidity—show the liquid turning gas-like as it is cooled, with ultralow viscosity, superdiffusive motion, and a breakdown of the Stokes–Einstein relation, all driven by quantum spreading of atomic positions. The result offers an exchange-free benchmark for what nuclear quantum effects alone can do in a quantum liquid.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Centroid approximation for nonlinear stress/energy-current correlations is unvalidated in HQDL; ultralow viscosity and gas-like VAF could be artifacts.","rationale":"The paper's central claim is that transport in distinguishable 4He becomes gas-like at low temperature solely through nuclear quantum effects. The observable evidence for this — ultralow shear viscosity, superdiffusion, monotonic VAF, SE breakdown — is generated by centroid molecular dynamics combined with the centroid approximation for Green–Kubo correlation functions. The reader identified this as the weakest assumption, and the manuscript itself acknowledges the lack of rigorous justification. My stress test agrees with that assessment and sharpens it: the approximation is applied precisely in the regime where it is least tested. The very fast stress relaxation (τ_s ~10^-3 ps) and the anomalous VAF shape are the quantities that would change the verdict if they are artifacts. At the same time, the paper deserves credit for explicitly flagging the limitation and for noting that the static delocalization is independent of the dynamical approximation, which supports the qualitative crossover picture even if the precise transport coefficients shift. Therefore the honest verdict remains CONDITIONAL: the qualitative claim is plausible and internally consistent, but it cannot be accepted as established until an independent dynamical method validates the extreme HQDL points. I do not see a basis for rejection — no internal inconsistency or undeniable artifact is identified — but I also do not see a basis for upgrading to ACCEPT without the proposed cross-check.","tokens_in":32262,"tokens_out":4851,"duration_ms":54274,"concrete_test":"For the HQDL state at 0.15 K, 5 bar (and, if feasible, 0.1 K, 1 bar), recompute ηs, λ, and C_v(t) using RPMD at the same N=256, N_b=500, same Aziz potential, or alternatively with CMD at N_b=1000 and N=2048. If the RPMD/alternative ηs differs from the CMD value by more than 50% in either direction, or if the VAF develops oscillations, the centroid approximation is not robust in HQDL and the gas-like classification is method-dependent. Agreement within ~10% would support the current conclusion.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Every quantitative transport coefficient in the paper (ηs, λ, the fractional SE exponent, and the second Frenkel line) flows through Eq. (14), which replaces the canonical correlation function of a nonlinear operator — the stress tensor (Eq. 18) and energy current (Eq. 23) — by the classical time correlation of centroid variables. The equality is exact only for operators linear in positions/momenta, and the authors explicitly concede in Sec. V F that 'a rigorous theoretical justification for the centroid approximation in the Green–Kubo formula, Eq. (14), is generally not available.' Their empirical validation is limited to He I above Tλ and para-H2, warm regimes where centroids are nearly classical. At 0.1–0.3 K in HQDL, λ_quantum reaches 6.6 Å and necklaces overlap; there is no benchmark showing that the centroid stress autocorrelation (which decays in ~10^-3 ps at 0.1 K) or the centroid VAF (which becomes monotonic and superdiffusive) matches the true Kubo dynamics. The paper also reports no finite-size (N=256 only) or Trotter (N_b=500 only) convergence for the dynamical quantities, and the thermal conductivity from the same scheme is overestimated relative to experiment. Since the 'one of the lowest viscosities' (7.83×10^-7 Pa·s) and the 'second Frenkel line' are defined by these centroid correlations, a failure of Eq. (14) in this regime would collapse the central claim. This is a known limitation, not a hidden error, but it is still the point on which the argument hinges.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports extensive path-integral centroid molecular dynamics (CMD) simulations of distinguishable (Boltzmann) 4He over 0.1–3.3 K and 1–60 bar, with transport coefficients obtained from centroid-approximated Green–Kubo formulas. It identifies two liquid states, LQDL and HQDL, and claims that upon cooling the system crosses from a conventional liquid to a gas-like liquid driven solely by nuclear quantum effects. Evidence includes a monotonic centroid VAF in HQDL, MSD exponent γ>1, shear viscosity minimum 7.83×10^-7 Pa·s, fractional Stokes–Einstein exponent ξ≈0.23–0.30, Prandtl number >1, and transport minima in ηs, λ, ν, and α. The authors introduce a 'second Frenkel line' in the subcritical region, distinguished from the supercritical Frenkel line by the driver being quantum rather than thermal fluctuations.","tokens_in":32625,"tokens_out":8845,"duration_ms":90882,"significance":"If the central claims are correct, the paper provides a striking qualitative result: a liquid can become gas-like upon cooling because of nuclear quantum delocalization, without Bose exchange or superfluidity. This would establish an exchange-free reference for interpreting transport anomalies in real and confined 4He. The study is systematic—115 state points, multiple independent observables, error bars from integration-window variation, and favorable comparison to He I in the high-temperature regime—which are genuine strengths. However, the quantitative claims and the very existence of the HQDL transport crossover rest on the centroid approximation for nonlinear correlation functions, which the authors concede is not rigorously justified. The significance is therefore conditional until that approximation is validated in the strongly quantum regime.","major_comments":[{"comment":"All collective transport coefficients (ηs, λ, ν, α, Pr, SE exponent) and the VAF-shape-based second Frenkel line flow through Eq. (14), which replaces canonical Kubo correlations of nonlinear operators—the stress tensor (Eq. 18) and energy current (Eq. 23)—by classical centroid correlations. The authors state in Sec. V F that rigorous justification is unavailable, and exactness for linear operators [13] does not apply here. The cited empirical validation (He I above Tλ [27], para-H2 [25,26]) is in the nearly classical regime, not at 0.1–0.3 K where λ_quantum reaches 6.6 Å and necklaces overlap. No benchmark is provided for the HQDL regime. Moreover, Sec. III H reports that the same scheme overestimates λ relative to experiment even in the validated He I regime, showing the approximation has quantitative error. I request a direct validation at representative HQDL state points, e.g., finit","section":"Eq. (14), Secs. II B, V F, III G–H"},{"comment":"No finite-size or Trotter convergence checks are reported for dynamical quantities; N=256 and Nb=500 are used exclusively. At 0.1 K and 1 bar, λ_quantum=6.6 Å, and the simulation box is only about 23 Å across, so necklace overlap and periodic boundary effects could be severe. The superdiffusive MSD exponent γ is fitted over t=20–60 ps (Sec. III E), and the VAF relaxation time τ_v≈4.4 ps (Sec. IV C) is comparable to the time needed to traverse a significant fraction of the box. The authors should report box lengths and provide N- and Nb-convergence data for the MSD, VAF, and stress autocorrelation at least at the 0.1 K, 1 bar and 5 bar state points.","section":"Sec. II C, Sec. III E, Sec. IV C"},{"comment":"The second Frenkel line is defined only by a qualitative change in the VAF from oscillatory to monotonic decay. No numerical criterion, uncertainty, or reproducible algorithm is given for T_F2, yet the line is drawn in Fig. 2 and used as a central diagnostic. The authors should specify a quantitative operational definition—e.g., the first zero of C_v(t), a threshold on the depth of the first minimum of the normalized VAF, or a fit parameter—and show the crossover points and their scatter. Without this, the existence and location of the second Frenkel line cannot be independently assessed.","section":"Sec. V B, Table III, Fig. 2"},{"comment":"The NVE simulations are labeled by the nominal pressure of the preceding NPT run, but Sec. III A explicitly states that the average pressure in the NVE runs deviates significantly for some state points. The isobaric temperature dependences in Figs. 4, 5, 7, 9, 10 and the extracted minima and crossover temperatures therefore mix different actual pressures. The authors should report the actual NVE pressures (or densities) for each state point and quantify the deviation, or replot the data against density. As it stands, the minima claimed at '1 bar' or '5 bar' may be shifted by the pressure drift.","section":"Sec. III A, Figs. 4–10"}],"minor_comments":[{"comment":"The partition function includes a 1/N! factor while the text repeatedly calls the particles 'distinguishable.' If the 1/N! is retained as the Gibbs-correction factor for Boltzmann statistics, this should be stated explicitly, since for truly distinguishable particles the factor should be absent (it cancels in the transport averages but is conceptually confusing).","section":"Eq. (1)"},{"comment":"The notation 'P r' and 'V AF' appears with extra spaces; it should read 'Pr' and 'VAF' consistently. Also, Table I uses 'P r <1' without space.","section":"Throughout"},{"comment":"The thermal diffusivity α uses C_P obtained from fitted enthalpy–temperature curves (Figs. S6–S7), but the fitting procedure and its uncertainty are not described. Since Pr = ηs C_P / λ depends on C_P, the error in C_P should be propagated into Pr and the stated crossover.","section":"Sec. IV A"},{"comment":"The comparison with Nakayama et al. [49] is plausible but under-supported: only two NVT temperatures and one NPT gas-like run are reported, and the original simulation conditions of [49] are not fully analyzed. The authors should provide a clearer demonstration that the Nakayama conditions lie on the gas side of the coexistence boundary.","section":"Sec. III B"},{"comment":"The 'Frenkel line' in the supercritical region is reproduced from Ref. 15 and not computed here; the caption should state this more prominently, since the reader may otherwise mistake it for a new result of this paper.","section":"Fig. 2(b)"},{"comment":"The order-of-magnitude estimate of the entropy from the KSS bound is clearly flagged as speculative, but it is not central to the paper and could be moved to the Supplementary Material to sharpen the main text.","section":"Sec. V A"}],"recommendation":"major_revision","confidential_remarks":"The paper is technically interesting and the simulation campaign is substantial, but the central quantitative claims depend on the centroid approximation in a regime where it has not been validated. This is a fixable issue—additional convergence and benchmark tests at representative HQDL points would either substantiate or refute the claims. I do not see grounds for rejection, but the current version is too conditional for acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The qualitative picture is probably right: distinguishable helium-4, without any Bose exchange, develops a gas-like transport regime at low temperature because quantum delocalization smooths the potential and kills the oscillatory backscattering that defines ordinary liquid dynamics. The evidence is consistent across several independent observables—monotonic VAF, superdiffusive MSD, fractional Stokes–Einstein exponent near 0.25, Prandtl number crossing above one, and transport minima in viscosity, conductivity, and diffusivity. That consistency is the paper's real strength, and the authors deserve credit for laying it out clearly and for flagging their main methodological weakness in Sec. V F rather than burying it.\n\nThe soft spot is exactly where the stress-test note lands: all quantitative transport coefficients flow from Eq. (14), replacing canonical correlation functions of nonlinear operators (stress, energy current) with centroid-variable correlations. For linear operators that replacement is exact; for these nonlinear ones it is not, and the authors say outright that rigorous justification is unavailable. Their validation comes from warm regimes—He I above Tλ and para-H2—where centroids are nearly classical. At 0.1–0.3 K in HQDL the necklaces overlap heavily, and the stress autocorrelation decays on a 10^-3 ps timescale; there is no benchmark showing the centroid version matches true Kubo dynamics in that regime. The viscosity minimum at 7.83×10^-7 Pa·s, the steep SE breakdown, and the second Frenkel line all depend on this approximation holding. That is not a hidden error, but it is load-bearing.\n\nOther issues are secondary but worth listing. N=256 with no finite-size check, Nb=500 with no Trotter-convergence check for dynamical quantities. The ξ exponents have no error bars. Thermal conductivity is overestimated against experiment, which the authors acknowledge and plausibly attribute to the centroid kinetic energy term. The comparison with Nakayama et al.'s monotonic VAF at 1.18 K is explained by ensemble sampling, but the explanation is speculative without more direct evidence. No code or trajectory data is released, only \"available upon reasonable request,\" so the central numbers cannot be independently checked.\n\nFor a reader working on quantum liquids or path-integral dynamics, this is worth a careful read. It is a serious computational study with a bold, clearly-stated claim, and the authors engage honestly with the prior literature. But the right verdict is conditional: the qualitative crossover is well-supported by static and dynamic trends, while the quantitative magnitudes—especially the \"one of the lowest viscosities\" language—need validation of the centroid approximation in the extreme regime, ideally via RPMD or an independent method, plus finite-size and Trotter checks. I would send it to peer review and ask for those tests before accepting the second Frenkel line as more than a plausible interpretation.","headline":"Plausible and interesting qualitative crossover, but every exciting number rides on an unvalidated centroid approximation, so treat the magnitudes as provisional.","tokens_in":33149,"tokens_out":1839,"would_cite":true,"duration_ms":22767,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82D15","82C31"],"pacs":[],"model":"deepseek-v4-flash","headline":"In distinguishable helium-4, nuclear quantum fluctuations alone turn the liquid gas-like on cooling—superdiffusion, ultralow viscosity, and a second Frenkel line, with no superfluidity.","keywords":["distinguishable helium-4","nuclear quantum effects","transport crossover","second Frenkel line","centroid molecular dynamics","quantum liquid without superfluidity","Stokes-Einstein breakdown","shear viscosity minimum"],"falsifier":"Perform a benchmark simulation of distinguishable 4He at 0.1–0.3 K and 1–10 bar using a method that does not rely on the centroid approximation for nonlinear operators—for example, a ring-polymer approach with an energy-current estimator based on the intermediate scattering function—and check whether the VAF remains monotonic, γ > 1, and η_s ≈ 10⁻⁷ Pa·s. Alternatively, measure diffusion and viscosity of 4He confined in extremely narrow nanopores where exchange is suppressed; if the predicted superdiffusion and viscosity minimum do not appear, the centroid approximation is the likely culprit.","tokens_in":32108,"feed_emoji":"🧊","tokens_out":4764,"duration_ms":41320,"temperature":0.7,"pith_summary":"The paper simulates helium-4 whose atoms obey Boltzmann statistics (distinguishable, no exchange) and finds two liquid states separated by a continuous crossover. At higher temperatures, the low quantum-dispersion liquid (LQDL) behaves like an ordinary liquid: oscillatory velocity autocorrelation, Stokes-Einstein relation, Prandtl number below one. On cooling below roughly 0.5 K, nuclear quantum delocalization strengthens and the liquid enters the high quantum-dispersion liquid (HQDL), which paradoxically behaves like a gas: superdiffusive motion, monotonic velocity autocorrelation decay, one of the lowest shear viscosities reported for any atomic liquid, and breakdown of the Stokes-Einstein relation. The crossover shows up as minima in shear and kinematic viscosity, thermal conductivity, and thermal diffusivity, and defines a second Frenkel line in the subcritical regime, distinct from the supercritical one driven by thermal motion. The claim is that quantum fluctuations alone—not Bose exchange or superfluidity—can produce gas-like fluidity.","feed_headline":"Cooling helium-4 turns its liquid gas-like—no superfluidity","feed_subtitle":"Simulations of distinguishable helium-4 find a second Frenkel line and record-low viscosity as quantum delocalization grows.","key_machinery":"The central object is the atomic necklace of path-integral centroid molecular dynamics: each atom is a ring polymer whose spatial delocalization is quantified by the quantum wavelength λ_quantum ≈ 2R_g, with expansion factor α_λ = λ_quantum/λ_dB. The transport coefficients (self-diffusion, shear viscosity, thermal conductivity, and derived kinematic viscosity, thermal diffusivity, Prandtl number) are computed via the centroid approximation to the Green–Kubo formula, where canonical correlation functions of stress and energy current are replaced by centroid-variable correlations. The crossover is diagnosed by the shape of the normalized velocity autocorrelation function (oscillatory vs. monot","core_discovery":"The central discovery is that distinguishable helium-4, modeled without Bose exchange, exhibits a transport crossover upon cooling rather than simply becoming more solid-like. In the HQDL state below about 0.5 K, the velocity autocorrelation function decays monotonically instead of oscillating, the mean-square displacement exponent exceeds unity (superdiffusion), the shear viscosity reaches 7.83×10⁻⁷ Pa·s at 0.15 K and 5 bar—below the experimental minimum of the normal-fluid component of He II—and the fractional Stokes-Einstein exponent drops to ξ ≈ 0.23–0.30. These gas-like transport properties arise from nuclear quantum delocalization, which widens the atomic 'necklaces' to quantum wavelen","pith_inferences":["The paper's mechanism—quantum delocalization smoothing the effective potential—suggests a testable prediction for other light-atom liquids (e.g., hydrogen or neon isotopes): if zero-point motion is large enough to prevent freezing, they too should exhibit a low-temperature gas-like transport crossover and a second Frenkel line.","The fractional Stokes-Einstein exponent ξ≈0.2–0.3 with no dynamical heterogeneity challenges the usual attribution of SE breakdown to heterogeneity; if confirmed, percolated necklace interpenetration would be a new structural route to SE breakdown, possibly relevant to metallic liquids and polymer networks.","Because the centroid approximation's validity is empirical, the quantitative claims (e.g., the exact viscosity minimum) could be tested by recomputing with a different dynamical approximation that avoids nonlinear centroid currents; if the monotonic VAF and superdiffusion persist, the crossover is robust.","If Bose statistics were switched on, the paper's estimate that the exchange timescale (≈h/(k_B T) ≈ 100 ps at 0.1 K) exceeds all relaxation times (τ_v ≈ 4.4 ps) suggests that transport in real He II below 0.1 K might be viewed as the same quantum-fluctuation-driven fluidity, with superfluidity adding a coherent component rather than being the sole origin."],"forward_implications":["If distinguishable helium-4 is realized (e.g., in nanopore confinement that suppresses exchange), its low-temperature liquid should show the gas-like signatures predicted here: superdiffusion, monotonic VAF, fractional Stokes-Einstein exponent ≈0.2–0.3, and Prandtl number >1.","The transport minima found in η_s, λ, ν, and α in the subcritical regime are a new class of minima, distinct from the well-known supercritical ones, and should appear as a 'reentrant gas-like' regime in any sufficiently quantum liquid that does not freeze.","The identification of a second Frenkel line extends the liquid–gas dynamical crossover concept to the subcritical low-temperature region, where the driving fluctuation switches from thermal to nuclear quantum.","For glassy states, HQDA should be a 'glassy yet fluidic' state: suppressed self-diffusion but viscosity and thermal conductivity close to the liquid values, providing a new criterion for vitrification.","A direct comparison with bosonic 4He would isolate the role of Bose exchange: the same ingredients without exchange should produce no superfluid transition and no λ divergence, but retain the low-viscosity HQDL behavior."],"fun_headline_variants":["Quantum jitter turns helium liquid gas-like","Chilling helium-4 creates gas-like liquid—no superfluidity","Helium's liquid acts like a gas when cooled—quantum twist","Quantum delocalization flips helium-4's transport to gas-like","Second Frenkel line marks helium-4's gas-like liquid on cooling"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The key load-bearing assumption is that the centroid approximation to the Green–Kubo formula, which replaces the true canonical correlation functions by centroid-variable correlations for nonlinear stress and energy-current operators, accurately reproduces the long-time decay of these correlation functions at 0.1–0.3 K; the authors state that no rigorous justification exists and validation is only for He I above T_λ and para-hydrogen.","fun_headline_variants_meta":{"raw":{"variants":["Quantum jitter turns helium liquid gas-like","Chilling helium-4 creates gas-like liquid—no superfluidity","Helium's liquid acts like a gas when cooled—quantum twist","Quantum delocalization flips helium-4's transport to gas-like","Second Frenkel line marks helium-4's gas-like liquid on cooling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000362,"raw_usage":{"total_tokens":1855,"prompt_tokens":870,"completion_tokens":985,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":614,"completion_tokens_details":{"reasoning_tokens":910}},"tokens_in":614,"tokens_out":985,"duration_ms":10416,"temperature":1.0,"reasoning_tokens":910,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T03:31:51.659593+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Perform a benchmark simulation of distinguishable 4He at 0.1–0.3 K and 1–10 bar using a method that does not rely on the centroid approximation for nonlinear operators—for example, a ring-polymer approach with an energy-current estimator based on the intermediate scattering function—and check whether the VAF remains monotonic, γ > 1, and η_s ≈ 10⁻⁷ Pa·s. Alternatively, measure diffusion and viscosity of 4He confined in extremely narrow nanopores where exchange is suppressed; if the predicted superdiffusion and viscosity minimum do not appear, the centroid approximation is the likely culprit.","supporting_citations":[],"review_version":1}