{"id":"ae6aec31-d135-49eb-9248-c3b83853f453","arxiv_id":"1909.02614","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A self-interacting scalar dark matter soliton around a supermassive black hole settles into a unique critical accretion flow with flux ~ r_s^2 m^4 / λ4 and survives for many Hubble times.","lead":"This paper calculates how fast a dark matter soliton made of a self-interacting scalar field falls into the supermassive black hole at a galaxy's center. It finds a unique critical accretion rate that is small enough for the soliton to survive far longer than the age of the universe.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The unique flux F_c ≈ 0.66 F_s hinges on a branch-switching argument that is illustrated but not rigorously established, and no numerical KG solution is provided; a targeted simulation would confirm or refute it.","rationale":"The paper presents a coherent analytic derivation: the ansatz (51), the constant-flux condition (84), and the two-branch structure yield a unique flux via a standard transonic argument. I checked the flux expression and found no algebraic inconsistency (Eq. (58) appears to have a typo with (1-2k^2) in the numerator, but Eq. (84) uses the denominator version required by Eq. (75)). The main gap is that the uniqueness proof is not fully rigorous (unimodality is unproven, boundary conditions are heuristic), and there is no numerical or experimental check. The paper itself acknowledges the approximate matching and calls for simulations in Sec. VI. This does not undermine the internal consistency, but it justifies a conditional verdict. I agree with the reader's conditional assessment, though I locate the weakest assumption less in the near-horizon WKB breakdown and more in the branch-switching/uniqueness argument and the absence of numerical verification.","tokens_in":22614,"tokens_out":39758,"duration_ms":394254,"concrete_test":"Run a spherically symmetric numerical evolution of the Klein-Gordon equation (49) on the fixed Schwarzschild background in Eddington coordinates, for m rs = 10 and m rs = 100 (satisfying the large-mass bound (40)), with quartic coupling and an initial profile matching the soliton core. After the transient, measure the time-averaged flux F = √(f h³) r² ⟨T^r_0⟩ at r ≈ 2.43 rs (x ≈ x*) and compare with F_c = 0.66 |F_s|. Agreeing fluxes within ~20% would validate the branch-switching selection; a substantially different or time-varying flux would falsify the uniqueness claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that a steady-state infall onto the BH selects the unique critical flux F_c = F_* F_s with F_* ≈ 0.66 (Eqs. (86), (91))—rests on the two-branch structure of the constant-flux equation (84). The argument in Sec. IV D assumes that F(k,x) has a single maximum in k at every radius, so that for F < F_c there are exactly two solutions k1(x) < k2(x), and that a continuous profile cannot switch branches. This is illustrated in Figs. 1-3 but not proven analytically; if F(k,x) had a different global shape, the uniqueness would fail. The branch assignment itself (k1 near the horizon, k2 at large radii) is justified by physical expectation rather than rigorous matching: the matching to the static soliton is explicitly approximate (Sec. IV D 2), leaving a nonzero velocity. The paper's own text in Sec. IV D 2 states that near the horizon k is 'close to zero', which is inconsistent with Eq. (98), where kc(1/4) = k_s ≈ 0.54; this indicates that the near-horizon branch assignment is not cleanly controlled. No numerical solution of the full nonlinear Klein-Gordon equation (49) is given, and Sec. VI calls for dedicated simulations. Therefore the quantitative value F_* ≈ 0.66 and the uniqueness statement lack independent verification. A simulation would either confirm the flux or reveal that the leading-order WKB solution does not realize the predicted steady state.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper analyzes steady-state accretion of a massive scalar field with a repulsive quartic self-interaction onto a Schwarzschild black hole, in the large scalar-mass limit. The authors introduce a Jacobi-elliptic ansatz, Eq. (51), reduce the nonlinear Klein-Gordon equation to algebraic conditions, and derive the flux function F(k,x), Eq. (84). They argue that the boundary conditions—matching to a static soliton at large radii and to an ingoing free-fall solution near the horizon—force the modulus k(x) to switch branches at the minimum of the maximum-flux curve, thereby selecting a unique critical flux Fc = F* Fs with F* ≈ 0.66, Eq. (91). They derive the density profile <ρ_φ> ∝ r^{-1} in the black-hole-dominated region, Eq. (110), and estimate the soliton lifetime tc ≫ t_H, Eqs. (118)–(119). The free-field case is treated separately and gives an unconstrained flux. The paper closes by noting that dedicated numerical simulations would be useful to confirm the nonlinear dynamics.","tokens_in":22895,"tokens_out":5768,"duration_ms":63062,"significance":"Should the central claim hold, the paper provides a concrete, essentially parameter-free prediction for scalar dark-matter soliton infall onto supermassive black holes: a unique flux scaling as F_c ∼ -r_s^2 m^4/λ4, a universal density slope in the BH-dominated region, and survival of the soliton for many Hubble times. The derivation is explicit and transparent, with the flux obtained from the equations of motion rather than fitted to observables. The comparison with the free-field case and with the hydrodynamic transonic solutions gives useful physical context. The main caveats are that the uniqueness argument is heuristic and the near-horizon validity of the large-mass ansatz is not fully controlled; the quantitative value F_* ≈ 0.66 should therefore be regarded as a well-motivated but unverified prediction until a numerical benchmark is provided.","major_comments":[{"comment":"The uniqueness of the critical flux rests entirely on the assumptions that F(k,x) has a single maximum in k for every radius and that a continuous steady profile cannot switch branches except at the peak. These facts are illustrated numerically but not proven analytically. If F(k,x) had additional extrema for some x, or if the steady solution were allowed to jump between branches, the selection argument would fail. Since the value F_* ≈ 0.66 is the central quantitative output, this gap should be closed by an analytic argument or, preferably, by a targeted numerical solution of the full Klein-Gordon equation (49). The paper itself calls for such simulations in Sec. VI, which reinforces the need for this check.","section":"IV D 1, Eq. (84), Figs. 1–3"},{"comment":"The near-horizon branch choice is not fully controlled. The text states that for x → 1/4 the modulus is on the lower branch and 'close to zero', yet Eq. (98) gives k_c(1/4) = k_s ≈ 0.54. The critical solution at the horizon is therefore not in the small-k regime used to justify why the flux is small. This inconsistency should be resolved, as it directly bears on the identification of k1 with the near-horizon branch and hence on the uniqueness argument.","section":"IV D 2 and Eq. (98)"},{"comment":"The large-mass ansatz (51) is assumed valid down to the Schwarzschild radius, but the phase β diverges logarithmically, Eq. (104), and β' diverges as (r − r_s/4)^{-1}. The paper checks density gradients, Eq. (28), but does not demonstrate that all subleading radial-gradient terms in the Klein-Gordon equation are negligible in this regime. If those terms become important near the horizon, the derived critical flux F_c could change. A controlled asymptotic expansion in powers of (m r_s)^{-1}, or a numerical integration, is needed to justify the leading-order solution all the way to the horizon.","section":"II B, IV E"},{"comment":"The large-radius matching to the static soliton is explicitly approximate: the text notes that 'there remains a nonzero velocity β′' and a nonzero flux. This is acknowledged, but the uniqueness of the outer branch k2 and the quantitative value of the transition radius rsg depend on the smallness of this residual velocity. No estimate of the matching error is given. The authors should quantify how small the residual velocity is compared with the free-fall velocity at the matching radius, or explain why the precise value does not affect the selected flux.","section":"IV D 2, Eq. (95)"}],"minor_comments":[{"comment":"The title and abstract contain typographical artifacts ('gal actic'); the manuscript should be proofread before resubmission.","section":"Title and abstract"},{"comment":"The notation ⟨ρφ⟩ used in Eqs. (46) and (107) is not defined before its first appearance; please define the average over the fast oscillation period when it is introduced.","section":"II B, Eq. (16)"},{"comment":"The caption is confusing: it describes both the F = Fc/3 curves and the F = Fc curves, but the phrase 'the inner dotted curves that meet at x⋆' could be read as referring to the dashed curves. Please make the distinction between dashed and dotted curves explicit.","section":"Fig. 3 caption"},{"comment":"The discussion calls the study 'fully nonrelativistic' while also referring to 'relativistic infall at small radii'; this wording is internally inconsistent and should be clarified.","section":"VI"},{"comment":"The identification v_r = πβ′/(2m) is introduced without explaining the factor π/2; a short justification after Eq. (59) would improve readability.","section":"IV D 3, Eq. (100)"}],"recommendation":"major_revision","confidential_remarks":"The paper is well written and the derivation is transparent, but the central quantitative claim F_* ≈ 0.66 is not yet independently verified. The branch-switching argument and the near-horizon validity of the WKB-type ansatz are the two load-bearing points where the result could fail. I would recommend that the editor obtain an assessment from someone with expertise in numerical solutions of the Klein-Gordon equation, since a targeted simulation would settle both issues."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper you want to know about: it derives a unique steady-state accretion flux for a quartic self-interacting scalar soliton falling into a Schwarzschild black hole, and finds the soliton survives many Hubble times. This is genuinely new in its specifics, and the derivation is explicit enough to follow line by line. The main caveat is that the uniqueness argument is visually compelling but not a proof, and there is no numerical solution of the full Klein-Gordon equation backing the leading-order WKB ansatz.\n\nWhat it does well: the free-field case reduces cleanly to independent-particle infall, a good sanity check. The interacting case uses a Duffing/elliptic-function ansatz, and the critical flux F_c = F_* F_s with F_* ≈ 0.66 and F_s ≈ r_s^2 m^4/λ4 is computed from elliptic integrals, not fitted. The density profile in the BH-dominated region scales as r^{-1}, distinct from the free-field r^{-3/2}, and the consequences for stellar dynamics and the soliton lifetime are worked out. The paper is honest about its approximations and explicitly calls for dedicated simulations.\n\nSoft spots: the branch-switching argument in Sec. IV.D is the load-bearing part. The authors show that F(k,x) has a single maximum at the relevant radii and that the boundary conditions put the solution on different branches at small and large radii, so only the critical flux permits a continuous k(r). This is the same logic as Bondi's transonic solution, and I find it plausible, but the 'single maximum' claim is supported by figures rather than analysis. A different global shape would break the uniqueness. The near-horizon branch assignment is argued from physical expectation—large velocity means small k—but the actual value at the horizon is k_s ≈ 0.54, which is not tiny; the text's 'close to zero' is a slip, not an error. The divergent phase near the horizon is a coordinate artifact and disappears in Eddington coordinates, so the ansatz is not obviously broken, though subleading radial gradients are not fully controlled.\n\nCitation pattern is fine: the earlier soliton result [85] is used as an outer boundary condition, and the critical flux is not fitted to anything.\n\nBottom line: this is a careful, honest physics paper with a concrete, falsifiable prediction. I would send it to peer review without hesitation. The constant 0.66 may or may not survive a full simulation, but the analytic framework deserves a serious referee and citation.","headline":"A careful analytic construction of a unique accretion flux for self-interacting scalar solitons; the branch-selection argument is the main soft spot, but it deserves a serious referee.","tokens_in":23450,"tokens_out":7066,"would_cite":true,"duration_ms":75896,"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":"A self-interacting scalar dark-matter soliton falling into a Schwarzschild black hole reaches a unique steady state with a critical flux of order $r_s^2 m^4/\\lambda_4$, and in that state it survives many Hubble times.","keywords":["scalar dark matter","soliton","supermassive black holes","quartic self-interaction","Jacobi elliptic function","constant flux selection","transonic accretion","soliton lifetime"],"falsifier":"A direct check would be to solve the full nonlinear Klein-Gordon equation for the quartic potential on a Schwarzschild background with a soliton-like outer boundary condition and see whether the steady infall approaches the predicted critical flux $F_c\\simeq 0.66\\, r_s^2 m^4/\\lambda_4$ and a $1/r$ density profile down to $r_s/4$; a complementary check is to carry the large-mass expansion to next order and see whether the minimum of the peak flux shifts.","tokens_in":22345,"feed_emoji":"🕳️","tokens_out":7830,"duration_ms":79163,"temperature":0.7,"pith_summary":"The paper asks what happens to a scalar-field dark-matter soliton when a supermassive black hole sits at its centre and steadily pulls the field inward. The paper's central claim is that, for a repulsive quartic self-interaction in the large-mass limit, the infall cannot be arbitrary: the requirement of a constant flux and the boundary conditions at the soliton and at the horizon select one unique critical flux, analogous to the transonic accretion solution in hydrodynamics. This flux is small, scaling as $F_c \\simeq 0.66\\, r_s^2 m^4/\\lambda_4$, so the soliton loses mass only slowly. The result matters because it determines whether such solitons can survive to the present day inside galaxies; the paper finds typical lifetimes $t_c \\gg t_H$, so scalar solitons around supermassive black holes are not destroyed by accretion. It also predicts a specific density profile, $\\langle\\rho_\\varphi\\rangle \\propto r^{-1}$ in the black-hole-dominated region, which is testable with stellar-dynamics constraints near galactic centres.","feed_headline":"Solitons outlive the universe around supermassive black holes","feed_subtitle":"Critical accretion rate is tiny: a quartic scalar soliton loses mass slowly enough to survive many Hubble times.","key_machinery":"The central object is the leading-order large-mass ansatz $\\varphi=\\varphi_0(r)\\,\\mathrm{cn}[\\,\\omega(r)t-K(r)\\beta(r),k(r)]$, in which $\\mathrm{cn}$ is the Jacobi elliptic function and $K(k)$ the complete elliptic integral of the first kind; it generalises the free-field cosine $\\cos(mt-s)$ to the anharmonic oscillator generated by the quartic potential. Substituting this ansatz into the nonlinear Klein-Gordon equation turns the radial problem into two algebraic equations: one fixes the amplitude $\\varphi_0$ in terms of the modulus $k$, and the other fixes the radial velocity $\\beta'$. The argument is carried by the flux function $F(k,x)$ that comes from the averaged energy-momentum conservation; its peak value as a function of radius has a unique minimum $F_c$, and that minimum is the only flux at which the solution can switch from the outer soliton branch to the inner horizon branch, in direct analogy with the transonic accretion solution of spherical hydrodynamics.","core_discovery":"Treating the scalar field as a coherent wave $\\varphi = \\varphi_0(r)\\,\\mathrm{cn}[\\,\\omega(r)t-K(r)\\beta(r),k(r)]$ and keeping only the leading-order radial gradients in the large-mass limit, the authors reduce the relativistic Klein-Gordon equation to two algebraic conditions for the amplitude and the modulus $k(r)$. Imposing a time-independent averaged flux $F$ through each spherical shell gives a flux function $F(k,x)$ whose maximum over $k$ has a global minimum $F_* \\simeq 0.66$ at $x_* \\simeq 2.43$, where $x=r/r_s$. Matching to the static soliton at large radii (upper branch $k_2$) and to the free-fall horizon condition at small radii (lower branch $k_1$) forces the flux to take exactly this critical value $F_c$; no other constant flux admits a continuous branch-switching profile. The resulting density profile falls as $1/r$ in the region dominated by the black hole, and the associated soliton depletion time exceeds the Hubble time for the parameters of interest.","pith_inferences":["The paper does not explore attractive quartic self-interactions; if the same branch-switching selection applies there, the negative pressure might allow a different critical flux or no steady transonic-type solution, possibly giving much shorter soliton lifetimes.","The predicted $1/r$ density profile could be probed indirectly through precise orbital precession of the closest stars or through light-bending effects near Sgr A*, since the dark-matter distribution changes the effective potential at the few-percent level.","The analytic method should extend to higher-order self-interactions such as $\\phi^6$: the elliptic ansatz would be replaced by a more general periodic solution, but the same minimum-of-the-peak-flux argument should select the accretion rate whenever a single branch switch exists."],"forward_implications":["Galactic scalar solitons with $\\rho_a\\sim 1\\,\\mathrm{eV}^4$ are not swallowed by their central supermassive black holes: the critical accretion flux gives $t_c\\sim 10^3\\, t_H (\\rho_s/\\bar\\rho_c)(\\rho_a/1\\,\\mathrm{eV}^4)^{-5/2}(M/10^8M_\\odot)^{-2}\\gg t_H$.","In the black-hole-dominated region the time-averaged energy density and the radial velocity both decay as $1/r$, instead of the $r^{-3/2}$ free-infall profile, so the dark-matter mass inside radius $r$ grows as $M_\\varphi(