{"id":"39bb752b-2198-46db-9ce2-0b4684ef26bf","arxiv_id":"2507.08583","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"The authors show that classical fluid analogies for Schrödinger-Newton systems are only consistent in a semi-classical limit and that a pseudo-Reynolds number can be defined for void dynamics.","lead":"This paper asks how closely the Schrödinger-Poisson system, used for fuzzy dark matter and large-scale structure, can be described as a classical fluid with viscosity and pressure. It finds that the analogy works only under restrictive conditions, and proposes a pseudo-Reynolds number for such systems.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (34) is an imposed formal redefinition, not a derived viscous balance; its claimed ℏ→0 limit is internally inconsistent, so the central viscosity analogy is unsupported.","rationale":"The reader's weakest-assumption analysis is exactly right, so my pass agrees rather than moves the verdict. Eq. (34) is the only bridge between the Madelung system and the Navier-Stokes form, and it is imposed to force the match rather than derived from the Schrödinger equation. My read sharpens the problem: the condition is not just 'may fail'; in 3D it requires a nongeneric pointwise alignment of ∇U and ∇²v, and the paper's own scaling is contradictory, with η of order ℏ from Eq. (37) but η→∞ in the ℏ→0 limit claimed in Section 4. A concrete test with an exact 1D solution can determine whether η is a positive finite material coefficient or merely a formal state-dependent ratio. I also confirm the reader's note that combining Eqs. (54)–(59) yields a +4ν/3 ρ0 Ptxx term, not the minus sign in Eq. (60), and Section 6's scaling evidence relies on a figure not present in the manuscript. These are correctable if the authors reframe the paper as a formal analogy and fix the derivations, so the existing CONDITIONAL verdict is appropriate.","tokens_in":8356,"tokens_out":16652,"duration_ms":189277,"concrete_test":"Take a known exact 1D Schrödinger-Poisson solution (for example, a Gaussian wave packet in a static potential or a stable soliton), compute N, φ, v, and U from it, and evaluate η(x,t) from Eqs. (37)–(41) for a sequence of ℏ values, rescaling φ so the classical velocity v stays fixed as ℏ→0. If η changes sign, diverges at zeros of v_xx, or does not approach a finite nonzero value in the semiclassical limit, then Eq. (34) is not a physical viscous balance and the paper's 'η→∞ consistency limit' claim is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is Eq. (34), (1/m)∇U = −η/ρ ∇²v, adopted from Fernández de Córdoba et al. (2016) to force the Madelung momentum equation (32) into the Navier-Stokes form (19). This balance is imposed, not derived from the Schrödinger equation. In three dimensions a single scalar η exists only if ∇U is pointwise parallel to ∇²v; since v = (ℏ/m)∇φ, that requires ∇(∇²N/N) ∥ ∇(∇²φ), a nongeneric constraint on the wave function. The paper never states or verifies this constraint, and its consistency claim is internally contradictory: Eq. (37) gives η of order ℏ, so η→0 as ℏ→0, while the final paragraph of Section 4 asserts η→∞ in that limit. Unless an extra divergence in ∇²v is assumed, both statements cannot hold. The identification is therefore a formal redefinition, not an emergent viscous property of Schrödinger-Newton dynamics; the Section 6 pseudo-Reynolds argument inherits this problem and also cites Fig. 1, which is not present in the manuscript.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper examines whether Schrödinger-Newton / Schrödinger-Poisson systems can be interpreted with classical fluid concepts, in particular viscosity. After a standard Madelung transformation, the authors identify the quantum potential with a viscous term in the momentum equation (Eq. 34), leading to a spatially dependent effective viscosity of order ℏ. They then study attenuation of acoustic waves in such a fluid and define a pseudo-Reynolds number for void dynamics, concluding that the fluid analogy is possible only under restrictions.","tokens_in":8621,"tokens_out":13609,"duration_ms":139638,"significance":"If the derivation were sound, this would be a useful clarification of the limits of fluid analogies for fuzzy dark matter, and it would connect the Madelung quantum-pressure term with viscous dissipation. The paper is concise and draws on a relevant literature. However, the central identification is imposed rather than derived, the ℏ→0 limit is internally inconsistent, and the wave equation contains a sign error; the pseudo-Reynolds-number section is not self-contained. The paper therefore does not yet establish its main conclusions.","major_comments":[{"comment":"Equation (34) is imposed to force Eq. (32) into the Navier-Stokes form Eq. (19); it is not derived from the Schrödinger equation. For a scalar η to exist in three dimensions, ∇U must be pointwise parallel to ∇²v. With v=(ℏ/m)∇φ this requires ∇(∇²N/N) ∥ ∇(∇²φ), a nongeneric constraint that the manuscript neither states nor verifies. The one-dimensional example sidesteps the issue, but the general claim in Section 4 that the system can be described by a Navier-Stokes equation with scalar viscosity is therefore not established.","section":"Section 4, Eq. (34)"},{"comment":"Equation (37) with Eqs. (38)-(40) gives η ∝ ℏ, since A(x) contains ℏ² and B(x) contains ℏ. Hence η→0 as ℏ→0. The final paragraph of Section 4 states the opposite: the Navier-Stokes description 'is only consistent in the limit ℏ → 0 in which case η → ∞'. These statements are mutually contradictory as written. In addition, combining Eq. (34) with the definitions of A and B yields η = -ρ A/B, not η = ρ A/B as printed in Eq. (37); this sign must be resolved before the effective viscosity can be given a physical interpretation.","section":"Section 4, Eqs. (35)-(37) and final paragraph"},{"comment":"Combining Eqs. (54), (56), (58) and (59) gives ν²/4 P_xxxx − ρ0 P_tt + (4ν/3) ρ0 P_txx = 0, not the equation printed with a minus sign before the P_txx term in Eq. (60). The dispersion relation Eq. (63) is therefore derived from a different equation and the acoustic attenuation example must be recomputed. There are also linearization inconsistencies in this section: Eq. (45) gives P ≈ (ν²/4ρ0) δρ_xx to first order, whereas Eq. (46) gives P ≈ (ν²/2ρ0) δρ_xx; Eq. (52) would follow from Eq. (45) only after the denominator ρ0+δρ is replaced by ρ0, and Eq. (55) retains a nonlinear denominator ρ0+2δρ before being rearranged into Eq. (56).","section":"Section 5, Eqs. (54)-(60)"},{"comment":"The pseudo-Reynolds-number conclusion is not self-contained. The section summarizes Gallagher & Coles (2022) and refers to Figure 1, but the manuscript provides no quantitative data, no definition of the ν used in the simulations, and no error analysis; the figure appears only as a caption placeholder in the version under review. Since the abstract advertises the pseudo-Reynolds number as one of the paper's conclusions, this section needs to present the supporting evidence or be explicitly demoted to a citation of prior work.","section":"Section 6"}],"minor_comments":[{"comment":"The second term in the viscous stress tensor should read ∂v_j/∂x_i, not ∂v_j/dx_i, and the index placement should be made consistent.","section":"Eq. (16)"},{"comment":"The phrase 'discussed in this thesis' should be 'discussed in this paper'.","section":"Section 4, paragraph before Eq. (33)"},{"comment":"The hydrostatic condition ∇P/ρ = ∇V/m should be motivated explicitly; with the sign convention of Eq. (6) this identification is not the usual hydrostatic equilibrium condition and the reader is left to reconstruct the convention.","section":"Section 4, Eq. (33)"},{"comment":"The notation for pressure is inconsistent: Eq. (43) introduces P as a vector, while Eqs. (45)-(60) treat P as a scalar. The mixed notation makes the linearization hard to follow.","section":"Section 5, Eqs. (43)-(56)"},{"comment":"The central references Fernández de Córdoba et al. (2016) and Gallagher & Coles (2022) are not summarized beyond the equations borrowed from them; a short statement of their derivations would make the argument more self-contained.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper has a clear scope and the topic is appropriate, but the central derivation currently contains sign errors and an internal contradiction in the ℏ→0 limit. The pseudo-Reynolds section is also too dependent on the companion paper Gallagher & Coles (2022), with no data presented in the current manuscript. These issues are fixable in principle, but they affect the paper's main claims rather than just the presentation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a readable, honest paper, but the one new result has a sign error, and the central consistency argument in Section 4 contradicts itself. The correct response is to send it back for corrections, not to reject it.\n\nWhat is actually new: Section 5 derives a damped wave equation for quantum-pressure fluctuations, Eq. (60), and gives an explicit dispersion relation. The derivation is routine once you accept the Madelung transformation and the equation of state, but I don't think that exact combination appears in earlier work. The paper is also careful to label Eq. (34) as a condition imposed to match the Navier-Stokes form, and it explicitly says a full fluid analogy would need heat transfer and entropy production. That is a fair and useful caveat.\n\nNow the soft spots. The reader's sign analysis is correct: combining Eqs. (54)–(59) yields a plus sign on the (4ν/3)ρ0P_txx term, not the minus in Eq. (60). That flips the sign of the imaginary term in the dispersion relation, so the predicted attenuation in Section 5 comes from their equations only if you make an additional sign convention. Easy to fix, but needed.\n\nMore importantly, Section 4 contains an internal contradiction. From Eq. (37), η is of order ℏ, so η→0 as ℏ→0. But the last paragraph of the section says the Navier-Stokes description is consistent only in the limit ℏ→0, in which case η→∞. Both cannot be true unless the denominator in Eq. (37) diverges as ℏ→0, which the paper neither states nor argues. This is load-bearing because the connection to the adhesion model rests on that limit.\n\nThe stress-test also notes that Eq. (34) is imposed, not derived, and that in three dimensions it requires ∇U ∥ ∇²v pointwise, which is not generic. I agree. The authors are open about using Fernández de Córdoba et al.'s condition, but they never test whether that alignment holds for realistic wave functions. So the viscosity interpretation remains a formal redefinition rather than an emergent property. That doesn't kill the idea, but it does mean the conceptual claim is weaker than the abstract suggests.\n\nSection 6 is a summary of the authors' earlier paper, and Fig. 1 is not present in the manuscript, so the pseudo-Reynolds scaling is not independently checkable here. That is a minor but annoying omission.\n\nWho should read this: people working on fuzzy dark matter simulations and wave-mechanical structure formation. It clarifies what you can and can't call viscosity in the SP system, which is useful even if the details need work. I would send it to a serious referee, because the issues are fixable and the topic has an audience. I would not cite it in its current form.\n\nRecommendation: engage it in peer review, but require the authors to correct the sign in Eq. (60), resolve the ℏ→0 limit contradiction, and either include Fig. 1 or provide the data behind the scaling claim.","headline":"Modest but honest paper; the new wave equation has a sign error and the central ℏ→0 claim is internally inconsistent—worth a referee but needs corrections.","tokens_in":9096,"tokens_out":8525,"would_cite":false,"duration_ms":83921,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that a Schrödinger-Newton system can be represented as a viscous Navier-Stokes fluid only in the semiclassical limit ℏ→0, where the effective viscosity diverges.","keywords":["Schrödinger-Newton system","Schrödinger-Poisson formalism","Madelung transformation","quantum pressure","effective viscosity","Navier-Stokes analogy","pseudo-Reynolds number","fuzzy dark matter"],"falsifier":"Run a full numeric Schrödinger-Poisson simulation of a one-dimensional void or halo, extract the Madelung density and velocity fields, and compare the local vector (1/m)∇U with -(η/ρ)∇²v for any reasonable choice of effective η: the ratio would need to be uniform in space for the balance to hold, and a spatially changing ratio would falsify the identification. Alternatively, measure the damping of small-amplitude density waves in the same simulation and check whether the attenuation matches the dispersion relation k=β-iα derived from Eq. (63); no attenuation would also falsify the viscous analogy.","tokens_in":8169,"feed_emoji":"🌊","tokens_out":8097,"duration_ms":83467,"temperature":0.7,"pith_summary":"This paper asks whether the Schrödinger-Newton (SN) system, and its cosmological Schrödinger-Poisson (SP) limit, can be described as a classical fluid with pressure and viscosity, in the way often assumed for fuzzy dark matter and large-scale structure. It argues that a Navier-Stokes description is possible but only under restrictive conditions: the gradient of the quantum potential must be forced to balance a viscous force term, the effective viscosity becomes spatially varying and of order ℏ, and the description is consistent only in the limit ℏ→0, where viscosity diverges and the flow is nearly incompressible. It then shows that small-amplitude sound waves in such a quantum fluid obey a damped wave equation with a definite attenuation coefficient, and that a pseudo-Reynolds number can be used to scale void-expansion dynamics. A careful reader cares because this sets boundaries on how far wave-mechanical descriptions of dark matter can be interpreted as ordinary fluids.","feed_headline":"Fluid viscosity in Schrödinger-Newton systems works only as ℏ→0","feed_subtitle":"Effective viscosity is spatially varying and of order ℏ; the Navier-Stokes picture survives only as ℏ→0.","key_machinery":"The argument is carried by the Madelung transformation ψ=ψ0 exp(S+iϕ), with density ρ=m|ψ|² and the quantum potential U=-(ℏ²/2m)∇²N/N, together with the imposed identity Eq. (34): (1/m)∇U=-(η/ρ)∇²v. This identity is the bridge that lets the quantum momentum equation be read as the irrotational Navier-Stokes equation, and it is what forces the viscosity to be spatially varying and of order ℏ. The same machinery produces the quantum pressure equation of state P=(ν²/2)∇²(√ρ)/√ρ used in the wave-damping calculation.","core_discovery":"On the paper's own terms, the central discovery is that the Madelung-transformed Schrödinger-Newton equations can be made to match a viscous Navier-Stokes equation only if one imposes the balance condition (1/m)∇U = -(η/ρ)∇²v, together with the hydrostatic condition ∇P/ρ = ∇V/m. This makes the effective viscosity η a spatially dependent quantity of order ℏ, not a constant material parameter. Consequently the fluid analogy is consistent only in the semiclassical limit ℏ→0, in which η→∞ and the wave-function's spatial variations are small, forcing the flow to be almost incompressible and isentropic. The paper also derives an attenuation coefficient for small-amplitude acoustic waves from the resulting viscous wave equation and defends a pseudo-Reynolds number R=ul/ν as a scaling tool for one-dimensional SP void evolution.","pith_inferences":["Going further: because the matching condition is imposed rather than derived, the apparent viscosity of the SP system is better read as a bookkeeping term in the Madelung mapping than as a genuine dissipative transport coefficient; simulations that show viscous damping may be displaying the mapping's structure rather than new physics.","Going further: the same balance condition could be tested in other Madelung-mapped wave equations, such as Bose-Einstein condensates or nonlinear optical systems, where the effective viscosity would predict a measurable relation between amplitude gradients and velocity; a failure there would suggest the SN case is not generic.","Going further: if the pseudo-Reynolds number is a true scaling invariant of void dynamics, it could be used to predict a transition from smooth void expansion to wave-dominated behaviour as the effective ℏ is varied, much as classical Reynolds numbers mark laminar-turbulent transitions.","Going further: because η diverges as ℏ→0, any finite-resolution SP simulation with a non-zero effective Planck constant automatically contains a physically large viscosity, so fuzzy-dark-matter resolution studies should report effective ℏ alongside grid scales."],"forward_implications":["The effective viscosity of a Schrödinger-Newton fluid must be treated as spatially varying and of order ℏ, so constant-viscosity Navier-Stokes codes cannot faithfully represent the full quantum system.","A consistent Navier-Stokes description exists only in the semiclassical limit ℏ→0, where η→∞ and the flow is nearly incompressible and isentropic; this is the same limit used in analytic adhesion-model treatments.","Small-amplitude acoustic waves in a fuzzy dark matter fluid are attenuated with coefficient α from the dispersion relation k=β-iα, suppressing small-scale structure in a way analogous to free-streaming in hot dark matter.","A pseudo-Reynolds number R=ul/ν can be defined for one-dimensional SP void expansion, with scaling relations among peak velocity, length, and ν that make it a useful dimensionless scaling parameter."],"supporting_citations":[{"why":"Supplies the two matching conditions (hydrostatic pressure balance and the quantum-potential/viscosity balance) that the paper borrows to turn the Madelung equations into Navier-Stokes form.","marker":"Fernández de Córdoba et al. (2016)"},{"why":"Prior work that defined the pseudo-Reynolds number and provided the one-dimensional void-expansion scaling test the paper extends and defends.","marker":"Gallagher & Coles (2022)"},{"why":"Established the wave-mechanical form with quantum pressure P and identified ν as having the dimensions of kinematic viscosity.","marker":"Short & Coles (2006a)"},{"why":"Companion treatment of the ν parameter as a viscosity-like term in the Schrödinger-Poisson system.","marker":"Short & Coles (2006b)"},{"why":"Introduced the Schrödinger-Poisson approach in cosmology, the system whose fluid analogy this paper interrogates.","marker":"Widrow & Kaiser (1993)"},{"why":"Source of the Reynolds number construction that the pseudo-Reynolds number adapts for the SP system.","marker":"Reynolds (1883)"}],"fun_headline_variants":["Schrödinger-Newton fluid link only valid as ℏ→0","Quantum fluid analogy for Schrödinger-Newton breaks down","Viscosity in Schrödinger-Newton is ℏ-order, not constant","Madelung fluid picture fails except in semiclassical limit","Schrödinger-Newton waves mimic viscous fluids only classically"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument rests on the assumption that the gradient of the quantum potential, a term built from the wave-function's amplitude, exactly balances the viscous force term -η/ρ∇²v; this balance is imposed by hand, and if it does not hold physically the viscosity interpretation collapses.","fun_headline_variants_meta":{"raw":{"variants":["Schrödinger-Newton fluid link only valid as ℏ→0","Quantum fluid analogy for Schrödinger-Newton breaks down","Viscosity in Schrödinger-Newton is ℏ-order, not constant","Madelung fluid picture fails except in semiclassical limit","Schrödinger-Newton waves mimic viscous fluids only classically"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000237,"raw_usage":{"total_tokens":1448,"prompt_tokens":828,"completion_tokens":620,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":444,"completion_tokens_details":{"reasoning_tokens":531}},"tokens_in":444,"tokens_out":620,"duration_ms":6427,"temperature":1.0,"reasoning_tokens":531,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T18:15:44.260279+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a full numeric Schrödinger-Poisson simulation of a one-dimensional void or halo, extract the Madelung density and velocity fields, and compare the local vector (1/m)∇U with -(η/ρ)∇²v for any reasonable choice of effective η: the ratio would need to be uniform in space for the balance to hold, and a spatially changing ratio would falsify the identification. Alternatively, measure the damping of small-amplitude density waves in the same simulation and check whether the attenuation matches the dispersion relation k=β-iα derived from Eq. (63); no attenuation would also falsify the viscous analogy.","supporting_citations":[],"review_version":1}