{"id":"f00be33d-5c64-46e3-b82f-330d6a340c56","arxiv_id":"2501.14244","paper_version":7,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":4,"one_line_summary":"SPH kernel smoothing is proposed as an inverse LES filter that lets a classical two-fluid model mimic microscopic quantum fluctuations in liquid helium-4.","lead":"This paper argues that a classical two-fluid model of liquid helium-4, solved with smoothed particle hydrodynamics, can mimic microscopic quantum fluctuations because SPH kernel smoothing acts like an inverse large-eddy-simulation filter. The relevance is a possible bridge between centimeter-scale helium flow simulation and quantum vortex physics, though the authors repeatedly label the analogy speculative.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"SPH truncation noise is asserted to substitute for physical microscopic fluctuations, but the only test (Eq. 29, Fig. 6) assumes white-noise third derivatives and never measures actual SPH field statistics; the paper itself calls the analogy speculative.","rationale":"The reader's weakest-assumption identification points to Eq. (29), and that is the same primary concern here: the white-noise character of the third derivative is assumed, not established for actual SPH helium-4 fields. I add two refinements. First, even the synthetic Monte Carlo test does not establish literal whiteness for finite h: the PSD of Gaussian-smoothed white noise is exp(−h²k²), flat only at low wavenumbers, so the Fig. 6 claim 'irrespective of h' is too strong. Second, the decisive missing evidence is not the formal Taylor identity but a quantitative comparison of f^SPH_epsilon from real simulations with the physical fluctuation spectrum of helium-4; without that, the substitution claim remains an analogy. The paper is internally consistent as a formal exercise and explicitly labels the physical interpretation speculative, so I do not move the verdict from CONDITIONAL to REJECT. The Eq. (39) concern about Biot–Savart being built into the definition of ω0 is real but secondary: it affects how strongly one can claim the Condiff term is an SGS model, not the internal mathematics. The Monte Carlo and simulation results are genuine evidence for the formal claims, but they do not close the statistical gap in the central physical assertion.","tokens_in":27743,"tokens_out":9748,"duration_ms":93021,"concrete_test":"Run the Fig. 3 counterflow SPH setup at kernel radii h∈{0.01,0.05,0.1,0.2}; from the particle data compute the kernel-smoothed field ⟨f⟩_W and a high-resolution reference f_ref (e.g., the same problem at h→0 or a well-resolved grid solution), form the error field ε=f_ref−⟨f⟩_W, and estimate its autocorrelation and PSD. Also compute the PSD of d³f/dx³ from the SPH fields. If the ACF is not delta-like at resolved scales, the PSD is not flat up to the SPH cutoff, or the error variance differs by more than an order of magnitude from the physical helium-4 fluctuation level (from GP/vortex-filament simulations or dynamic structure-factor data), then the §3.1 substitution claim is quantitatively unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the substitution claim in §3.1: SPH truncation error f^SPH_epsilon can serve as a proxy for physical subscale fluctuation f^LES_epsilon. For this, conditions (i)–(iii) are required; only condition (ii) is tested, and only for a synthetic input. Equation (29) assumes d³f/dx³ is spatially uncorrelated Gaussian white noise, and the Monte Carlo convolution of that noise with a Gaussian kernel is then shown to be \"white\" (Fig. 6). Mathematically, however, smoothing white noise gives PSD |φ̂(k)|² = exp(−h²k²), which is flat only for k ≪ 1/h and is strongly colored for h=0.2 at high wavenumbers. More importantly, real SPH fields of helium-4 contain vortices, interfaces, and particle disorder; their third derivatives are never measured in the paper. No quantitative estimate is given of f^SPH_epsilon in the actual counterflow/vortex-lattice simulations, nor a comparison with the physical fluctuation spectrum of helium-4 at the corresponding scales. The paper explicitly concedes in the abstract and §3.1 that the resemblance \"lacks first-principle justification and should be viewed as a speculative analogy.\" The formal Taylor identities in §3.1–3.2 are internally consistent, but the physical substitution claim rests on an untested statistical-equivalence assumption. A secondary formal issue is Eq. (39), where the 1/|r| convolution for ω0 is postulated, so the Biot–Savart form in Eq. (43) is to some degree built in rather than derived.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper argues that classical and quantum hydrodynamic two-fluid models of liquid helium-4 can be connected by scale transformations: LES filtering maps microscopic to macroscopic scales, while SPH smoothing provides an inverse transformation. It claims that SPH truncation errors can, under specified conditions, stand in for physical microscopic fluctuations, and that the Condiff rotational-viscosity term can be interpreted as a subgrid-scale model containing Biot-Savart vortex interactions. The manuscript reviews the author's earlier SPH simulations of vortex lattices and counterflow, derives the formal Taylor-expansion relation between kernel smoothing and LES filtering, presents a Monte Carlo test of the whiteness of the SPH truncation error, and concludes with a speculative picture of the microscopic composition of liquid helium-4.","tokens_in":28119,"tokens_out":8601,"duration_ms":81580,"significance":"If the substitution claim were quantitatively established, the paper would offer a practical multiscale bridge between classical CFD simulations and quantum two-fluid behavior. The formal kernel/filter correspondence in Section 3.1 is standard and useful, and the paper is commendably explicit about the speculative character of the analogy. However, the load-bearing statistical equivalence is not demonstrated: the whiteness test assumes the very property it claims to verify, and the Biot-Savart form in Section 3.2 is largely built into the postulated convolution. The paper therefore has value as a clearly framed formal analogy and as a research program statement, but not yet as an established physical connection.","major_comments":[{"comment":"The Monte Carlo test does not establish condition (ii). The quantity fε^(2)(x) defined in Eq. (30) is a convolution of the assumed white noise with the deterministic kernel r^3 W(r,h), so its power spectral density is proportional to |K̂(k)|^2 and is not flat for wavenumbers k ≳ 1/h when h = 0.2. In addition, the Taylor remainder contains an integration point θ = θ(x,r) that depends on x and r, so treating it as a fixed convolution is not mathematically justified. The paper also never measures the third-derivative statistics of actual SPH fields in the helium-4 simulations, so the white-noise property of the SPH truncation error is assumed rather than demonstrated.","section":"§3.1, Eq. (30) and Fig. 6"},{"comment":"The central assertion that SPH truncation errors can substitute for physical microscopic fluctuations is not supported by quantitative evidence. No estimate of f_SPH^epsilon is given for the counterflow or vortex-lattice simulations, and no comparison is made with the expected spectrum of thermal or quantum fluctuations in helium-II at the corresponding scales. The manuscript itself states in the abstract, §3, and §4 that the resemblance 'lacks first-principle justification' and should be viewed as a 'speculative analogy.' The paper should either supply a quantitative comparison or consistently present the result as a formal analogy rather than as a demonstrated connection.","section":"Abstract and §3.1, conditions (i)–(iii)"},{"comment":"The 1/|r| convolution defining the macroscopic spin field ω0 is postulated rather than derived, and it introduces an unspecified constant C_q. Because Eq. (43) follows directly from inserting the filament vorticity (42) into this postulated convolution, the Biot-Savart form is built into the construction. To make the SGS identification substantive, the convolution should be derived from a filtering operation on the microscopic vortex field, and C_q should be fixed by a physical normalization such as the quantum of circulation.","section":"§3.2, Eq. (39)"},{"comment":"The assertion that the second term F[2η_r ∇×ω0] vanishes for quantum vortices is not justified for singular filament vorticity. The curl of a delta-distributed vorticity field, as in Eq. (42), does not automatically vanish in a distributional sense, so the recurrence argument leading to Eq. (43) requires a careful treatment of the filament singularity. Without that treatment, Eq. (43) is an algebraic restatement of Eq. (39) rather than a derivation.","section":"§3.2, Eq. (41)"},{"comment":"The counterflow velocity profiles are reproduced by adjusting the initial distribution of viscous particles to match the target vortex-line-density profile, so the tail-flattened profile is a prescribed input rather than an emergent prediction. This does not invalidate the formal analogy, but the results of Fig. 3 cannot be cited as independent evidence for the inverse coarse-graining mechanism, as is done in §3.3.","section":"§2, Fig. 3 and §3.3"}],"minor_comments":[{"comment":"The notation does not visually distinguish the true value, the kernel-filtered value, and the LES-filtered value, because the overbars are lost or used inconsistently; please define explicit symbols such as f, f^W, and f^G throughout.","section":"Eqs. (20)–(28)"},{"comment":"The prefactor, sign, and θ dependence of the Taylor remainder term fε^(2) are not derived; please show the intermediate steps leading from Eq. (25) to Eq. (30).","section":"§3.1, Eq. (30)"},{"comment":"The Monte Carlo estimators for the autocorrelation function and power spectral density are not fully specified; please state the normalization, frequency range, sample size, and error bars so the whiteness claim can be checked.","section":"Fig. 6"},{"comment":"The phrase 'local field ionization' is used without definition or supporting discussion; please clarify what is meant and how it connects to the cited helium-3 and ion experiments.","section":"§3.3"},{"comment":"No code or data availability statement is provided for the simulations shown; please add a reproducibility statement, even if only as a link to prior work or a data repository.","section":"Figs. 1–4"}],"recommendation":"major_revision","confidential_remarks":"The paper leans heavily on the author's prior publications (Refs. [27]–[30]) without shipping code or data, and the abstract claims more than the caveats allow. If the journal is willing to publish speculative theory with clearly stated limitations, a major revision that reframes the central claim as a formal analogy and adds the requested quantitative checks could make this acceptable. If the journal requires a validated physical result, the current manuscript would be a better fit as a perspective or research-program article."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the quick read. This is a speculative, interpretive paper, not a derivation, and the author knows it. The genuinely new move is reading the Condiff rotational viscosity term in the classical two-fluid model as a subgrid-scale closure whose 1/|r| convolution form routes directly into a Biot-Savart vortex interaction term. I haven't seen that identification in the literature the paper cites, and the vector algebra leading to Eq. (43) is internally consistent. The paper also frames SPH kernel smoothing as an 'inverse coarse-graining' operation, the reverse of the LES filter. The Taylor-expansion mapping between the two is standard, and the author is careful to say the equality f_SPH_eps = -f_LES_eps is formal and only holds when particle approximation and numerical errors are ignored.\n\nWhat the paper does well is its level of candor. The abstract, §3.1, and §3.3 all say the resemblance between SPH truncation noise and physical microscopic fluctuations lacks first-principle justification and should be read as a speculative analogy. That is not a paper overclaiming.\n\nThe soft spot is the load-bearing statistical claim. The substitution requires conditions (i)–(iii), and the only test is for condition (ii), inside a Monte Carlo setup where the third derivative of f is assumed to be white noise. Smoothing white noise with a Gaussian kernel leaves the PSD roughly flat only for k ≪ 1/h; at h=0.2 the high-wavenumber end is already colored. More importantly, real helium-4 SPH fields contain vortices, interfaces, and particle disorder, and the paper never measures f_SPH_eps in its own counterflow or vortex-lattice simulations, nor compares it to the physical fluctuation spectrum of helium-4. So the proposed mechanism is plausible, but it is not demonstrated. The Biot-Savart result is also partially built in: Eq. (39) postulates a 1/|r| convolution for the macroscopic spin field, so the vortex-interaction form follows from that choice. The counterflow profile reproductions rely on adjusted initial particle distributions, and the simulations themselves live in prior papers without shipped code or data. Those are real empirical holes, though the paper is upfront about most of them.\n\nWho should read it: people working on SPH as a multiscale method, or on classical models of helium-4, will find a useful review and a concrete proposal to test. But the central claims need quantitative tests: measure SPH truncation-error statistics in actual simulations, compare with helium-4 spectra, and show that the Biot-Savart SGS term changes the resolved dynamics in a way that matches experiments. That is referee work, not desk-reject work. Send it to peer review and ask for those tests.","headline":"A self-aware, interpretive paper that reads SPH truncation error as inverse LES and the Condiff rotational viscosity as a Biot-Savart SGS closure, but the central statistical equivalence is never tested on real helium-4 fields.","tokens_in":28648,"tokens_out":3607,"would_cite":false,"duration_ms":31512,"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 claims that smoothed particle hydrodynamics, run on the classical two-fluid model, produces kernel-truncation noise that can stand in for the microscopic fluctuations of quantum helium-4, and that the spin-conserving viscosity…","keywords":["liquid helium-4","two-fluid model","smoothed particle hydrodynamics","large eddy simulation","subgrid-scale model","spin angular momentum conservation","quantum vortices","inverse coarse-graining"],"falsifier":"Compute the third derivative of a physical field from an actual SPH simulation of the two-fluid model and measure its spatial autocorrelation and power spectrum: if the autocorrelation is not delta-like or the spectrum is not flat across the relevant wavenumbers, the white-noise premise of Eq. (29) fails, so the truncation error cannot serve as a statistically faithful substitute for microscopic fluctuations. A complementary check is to compare the velocity-fluctuation statistics produced by varying kernel radius $h$ against experimental or high-resolution quantum simulation data for superfluid helium-4.","tokens_in":27475,"feed_emoji":"🌀","tokens_out":5742,"duration_ms":48794,"temperature":0.7,"pith_summary":"The paper argues that the classical and quantum two-fluid descriptions of cryogenic helium-4 are connected by scale transformations: large-eddy-simulation filtering projects the microscopic quantum description onto the macroscopic scale, while the smoothing step of smoothed particle hydrodynamics performs the inverse transformation, projecting macroscopic fields back into microscopic-looking fluctuations. The load-bearing consequence is that SPH kernel truncation error, which grows with kernel radius, can serve as a stand-in for the subgrid fluctuations that a quantum treatment would supply, so classical SPH simulations can reproduce quantum-looking vortex lattices and counterflow profiles. The paper further claims that the spin-angular-momentum-conserving rotational viscosity term is formally a subgrid-scale model, and that under a point-vortex approximation it reduces to Biot-Savart vortex interactions, embedding quantum vortex physics into a classical two-fluid solver. The author is explicit that the resemblance between numerical noise and physical fluctuations is a speculative analogy rather than a first-principle derivation.","feed_headline":"SPH noise can substitute for quantum fluctuations in helium-4","feed_subtitle":"A scale-transformation argument makes classical two-fluid SPH reproduce quantum-like vortex and counterflow behavior.","key_machinery":"The central object is the SPH kernel-approximation error $f^{\\rm SPH}_\\epsilon$, defined by $f = \\bar f + f^{\\rm SPH}_\\epsilon$ after replacing the Dirac delta function with a Gaussian kernel. The Taylor expansion of this replacement shows that $f^{\\rm SPH}_\\epsilon$ is a sum of even powers of the kernel width $h$, so larger smoothing widths produce larger apparent fluctuations. A second machinery element is the white-noise assumption of Eq. (29), that the third derivative of the field is Gaussian white noise; Monte Carlo convolution with the kernel then yields white-noise autocorrelation and flat power spectra. The third element is the rotational viscosity term $2\\eta_r\\nabla\\times\\omega_0$ in the Condiff viscosity model, which, when $\\omega_0$ is defined via a $1/|r|$ convolution, is rearranged into a Biot-Savart integral over quantum vortex filaments, turning the term into a subgrid-scale closure.","core_discovery":"On the paper's own terms, the discovery is that SPH discretization of the classical two-fluid model is an inverse coarse-graining operation: the same Taylor expansion that defines the LES filtering error defines the SPH kernel approximation error, with opposite sign when the Gaussian kernel is used for both. Thus $f^{\\rm SPH}_\\epsilon$ corresponds to $-f^{\\rm LES}_\\epsilon$, and the truncation error behaves like a white-noise subgrid fluctuation when the third derivative of the field is spatially uncorrelated Gaussian white noise. This makes the kernel radius $h$ the control knob for fluctuation amplitude, which explains why the author's simulations reproduce microscopic-scale fluctuations at macroscopic scales. The spin-angular-momentum-conserving term in the viscous fluid equation is shown to correspond to the divergence of the SGS stress tensor, and with the macroscopic spin field defined by a $1/|r|$ convolution, the rotational viscosity term becomes a Biot-Savart integral over quantum vortex filaments, giving a concrete SGS closure for the two-fluid model.","pith_inferences":["The white-noise equivalence suggests a design principle: SPH truncation noise could be deliberately calibrated, rather than suppressed, as a cheap stochastic subgrid model for weakly compressible flows beyond helium-4.","The same scale-transformation argument should apply to other particle methods, such as MPS, where density fluctuations are also a few percent; if their truncation errors share the white-noise property, the analogy generalizes.","A testable extension is to run two SPH simulations of the same macroscopic helium-4 flow with different kernel radii and compare fluctuation amplitude to the even-power scaling in $h$ predicted by the Taylor expansion; a clean power law would support the mechanism, and a failure would expose missing corrections.","Because the Biot-Savart result follows from choosing a $1/|r|$ convolution for the spin field, a different, physically motivated coarse-graining kernel would produce a different closure; checking which form best matches mutual-friction data would distinguish the model from alternatives."],"forward_implications":["Classical hydrodynamic two-fluid SPH simulations can reproduce quantum-looking vortex lattices and counterflow profiles because kernel smoothing supplies the missing subgrid fluctuation statistics, with amplitude controlled by the kernel radius $h$.","The spin-angular-momentum-conserving viscosity term can be read as an SGS model: it transfers small-scale vortex motion into the macroscopic equations, and under point-vortex quantization it reduces to Biot-Savart interactions.","The normal fluid should be viewed as a mixture of inviscid and viscous fluid particles; molecular viscosity acts at microscopic scales, while large-scale effective viscosity is dominated by eddy viscosity, so in laminar regimes the normal fluid can be treated as inviscid.","Larger kernel radii amplify microscopic fluctuations but require higher particle density to maintain the regularity condition, so reproducing quantum-like fluctuations carries a computational cost."],"supporting_citations":[{"why":"Supplies the polar-fluid Navier-Stokes form containing the rotational viscosity term $2\\eta_r\\nabla\\times\\omega_0$ that is later identified as the SGS closure.","marker":"[38]"},{"why":"Gives the SPH discretization of the spin-angular-momentum-conserving equation that the two-fluid model uses.","marker":"[39]"},{"why":"Establishes the prior result that SPH smoothing can be viewed as LES filtering, which the inverse coarse-graining argument builds on.","marker":"[73]"},{"why":"Reports the vortex-lattice reproduction that motivates the claim that quantum phenomena can emerge from the classical SPH model.","marker":"[27]"},{"why":"Gives the counterflow simulations and the formal bridge between the quantum and classical equations that this paper extends.","marker":"[28]"},{"why":"Provides the derivation of Eq. (14) and the physical meaning of $\\omega_0$ used in the rotational-viscosity analysis.","marker":"[29]"},{"why":"Supplies the Taylor expansion of the kernel approximation that yields the even-power series in $h$ for the truncation error.","marker":"[70]"},{"why":"Provides the regularity condition that limits how large the kernel radius can be while the truncation-error expansion remains valid.","marker":"[75]"}],"fun_headline_variants":["SPH noise mimics quantum fluctuations in helium-4","Classical SPH reproduces quantum-like fluctuations in superfluid","Helium-4: SPH truncation error as quantum proxy","Inverse coarse-graining: SPH noise as quantum stand-in","SPH kernel error mimics quantum vortex dynamics in He-4"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quantitative bridge rests on assuming that the third derivative of the physical field is spatially uncorrelated Gaussian white noise; if real SPH fields in cryogenic helium-4 do not have white-noise third derivatives, the claim that kernel truncation error reproduces microscopic fluctuations loses its basis, and the Biot-Savart result additionally depends on choosing a $1/|r|$ convolution for the spin field.","fun_headline_variants_meta":{"raw":{"variants":["SPH noise mimics quantum fluctuations in helium-4","Classical SPH reproduces quantum-like fluctuations in superfluid","Helium-4: SPH truncation error as quantum proxy","Inverse coarse-graining: SPH noise as quantum stand-in","SPH kernel error mimics quantum vortex dynamics in He-4"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000222,"raw_usage":{"total_tokens":1512,"prompt_tokens":1065,"completion_tokens":447,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":681,"completion_tokens_details":{"reasoning_tokens":361}},"tokens_in":681,"tokens_out":447,"duration_ms":4015,"temperature":1.0,"reasoning_tokens":361,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T15:14:52.160840+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the third derivative of a physical field from an actual SPH simulation of the two-fluid model and measure its spatial autocorrelation and power spectrum: if the autocorrelation is not delta-like or the spectrum is not flat across the relevant wavenumbers, the white-noise premise of Eq. (29) fails, so the truncation error cannot serve as a statistically faithful substitute for microscopic fluctuations. A complementary check is to compare the velocity-fluctuation statistics produced by varying kernel radius $h$ against experimental or high-resolution quantum simulation data for superfluid helium-4.","supporting_citations":[{"cited_title":"M ¨uller, D","cited_arxiv_id":null,"evidence_quote":"Gives the SPH discretization of the spin-angular-momentum-conserving equation that the two-fluid model uses."},{"cited_title":"Di Mascio, M","cited_arxiv_id":null,"evidence_quote":"Establishes the prior result that SPH smoothing can be viewed as LES filtering, which the inverse coarse-graining argument builds on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the vortex-lattice reproduction that motivates the claim that quantum phenomena can emerge from the classical SPH model."},{"cited_title":"Tsuzuki, Theoretical framework bridging classical and quantum mechanics for the dynamics of cryogenic liquid helium-4 using smoothed- particle hydrodynamics, Phys","cited_arxiv_id":null,"evidence_quote":"Gives the counterflow simulations and the formal bridge between the quantum and classical equations that this paper extends."},{"cited_title":"Tsuzuki, A hydrodynamic approach to reproduce multiple spinning vortices in horizontally rotating three-dimensional liquid helium-4, Phys","cited_arxiv_id":null,"evidence_quote":"Provides the derivation of Eq. (14) and the physical meaning of $\\omega_0$ used in the rotational-viscosity analysis."},{"cited_title":"Stranex, S","cited_arxiv_id":null,"evidence_quote":"Supplies the Taylor expansion of the kernel approximation that yields the even-power series in $h$ for the truncation error."},{"cited_title":"Imoto, Truncation error estimates of approximate operators in a generalized particle method, JJIAM 37 (2) (2020) 565–598","cited_arxiv_id":null,"evidence_quote":"Provides the regularity condition that limits how large the kernel radius can be while the truncation-error expansion remains valid."}],"review_version":1}