{"id":"913b0a35-5947-4953-ae60-93b567a0754b","arxiv_id":"2412.15428","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"For circular orbits in ultralight dark matter, the dynamical friction on a Plummer sphere is the point-probe result multiplied by the cluster's momentum-space form factor, with deviations growing at large size-to-orbit ratios and high Mach numbers.","lead":"Globular clusters moving through ultralight dark matter feel a drag force that depends on their finite size. This paper derives analytic formulas for that drag for clusters modeled as Plummer spheres, and shows when the finite-size correction becomes important for circular orbits.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The point-probe limit of Eq. (15) is ill-defined because lmax=pi*r0/lp diverges as lp->0, so the 'point probe' reference curves must use some other undisclosed cutoff; the claimed finite-size deviations may be a cutoff artifact.","rationale":"The paper's formal derivation of the total force in Eqs. (10)-(17) is clean: replacing the point source by the Plummer form factor in the linear-response kernel and in the force density is a natural Fourier-space construction, and the point-probe expressions are formally recovered in the lp -> 0 limit of the integrand. The reader's weakest_assumption rightly flagged the inherited angular cutoff and the neglected self-gravity term. My stress test sharpens the cutoff issue into a concrete internal inconsistency: the point-probe limit of Eq. (15) is not attainable by setting lp = 0 because lmax = pi*r0/lp diverges. Since the point-probe curves in Figs. 2-5 must use another cutoff, the comparison is not controlled. This is load-bearing because the paper's headline conclusion is precisely that finite size matters only above lp/r0 ~ 10^-2 and strongly at large Mach number; both statements are statements about the difference between two curves that use different regulators. The proposed test settles the matter by recomputing both curves with a shared, physically motivated cutoff. If the deviations persist with a shared cutoff, the central claim survives and the paper can be accepted after clarifying the cutoff. If they do not, the claimed finite-size effect is an artifact and the conclusion would need substantial revision. I therefore maintain CONDITIONAL rather than ACCEPT or REJECT: the issue is specific and fixable, but the numerical evidence as presented is not yet trustworthy.","tokens_in":10566,"tokens_out":9617,"duration_ms":89253,"concrete_test":"For the parameters of Fig. 3 (r0 = 3*hbar/(m*c_s), lp/r0 = 0.1, Mach M = 2), recompute Eq. (15) for both the point probe and the Plummer sphere using identical cutoffs: first lmax = pi*r0/xi with xi = hbar/(m*c_s) (ULDM healing length), then lmax = 200, and also lmax = 1000 for the point probe to check convergence. If the ratio F_Plummer/F_point changes by more than about 10% between the two cutoff choices, the reported deviations are cutoff artifacts. Also rerun the lp/r0 = 10^-2 case with the same shared cutoff to see whether the claimed threshold survives.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central numerical comparison depends on the angular-momentum cutoff lmax = pi*r0/lp introduced in Eq. (15), where lp is the Plummer radius. The text claims the point-probe limit is recovered by setting lp -> 0 in Eqs. (20)-(21), but in that limit lmax diverges, so Eq. (15) has no well-defined point-probe limit. The point-probe reference curves in Figs. 2-5 must therefore have been computed with a different, unspecified cutoff, presumably inherited from Ref. [24] where the analogous lp denotes the ULDM healing/de Broglie scale rather than a body radius. For supersonic motion the imaginary-part sum in Eq. (20) is logarithmically sensitive to this cutoff, while for subsonic motion the real-part sum also depends on how many high-l modes are kept. The Plummer form factor rho_Pl(k) in Eq. (17) already suppresses high-k modes, so using lmax = pi*r0/lp in addition is an independent regulator, not the point-probe limit of the same expression. Consequently, the deviations between Plummer and point-probe forces claimed in Figs. 3-5 and the threshold lp/r0 ~ 10^-2 are not isolated finite-size effects; they conflate the physical form factor with an ad hoc truncation. The paper's own limitation section confirms other idealizations (homogeneous ULDM, linear response, spherical symmetry), but those are secondary; the cutoff inconsistency directly affects the headline comparison.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives analytic formulas for the dynamical friction force acting on a Plummer sphere moving on a circular orbit in a homogeneous ultralight dark matter (ULDM) superfluid. Working in linear response and neglecting the self-gravity term in the density-perturbation equation, the authors show that the force is the point-probe expression with an additional momentum-space form factor [ρ_Pl(k lp)/M]^2 inserted in the k-integral, leading to analytic expressions for the imaginary part and numerically evaluated real parts of the relevant mode sums (Eqs. (14)-(21)). The paper then computes radial and tangential force components for several Plummer radii and Mach numbers, concluding that the force is essentially unchanged from the point-probe result when lp/r0 is below about 10^-2, and that visible deviations appear for lp/r0 = 5×10^-2 and 10^-1 and at large Mach numbers.","tokens_in":10850,"tokens_out":8667,"duration_ms":70621,"significance":"If the central comparison is robust, the paper provides a useful semi-analytic tool for estimating finite-size corrections to dynamical friction for globular clusters in ULDM models, complementing the point-probe analyses of Refs. [10,23,24]. Strengths include the clean Fourier transform of the Plummer profile in Eq. (2), the explicit reduction of the integrand to the point-probe expression in the limit lp→0, and the absence of fitted parameters; the derivation from Eq. (3) to Eq. (20) is coherent and reproducible. However, the numerical comparison against the point-probe reference is compromised by an undefined cutoff procedure, as detailed in the major comments, so the headline finite-size threshold and the deviations reported in Figs. 3-5 are not yet established. The paper also transparently lists several idealizations (spherical symmetry, homogeneity, linear response) in its conclusions, though it omits the neglect of the self-gravity term from that list.","major_comments":[{"comment":"The point-probe limit is not well-defined in the manuscript's own formalism. The sum in Eq. (15) extends to ℓmax = π r0 / lp, which diverges as lp → 0. However, the text states that Fig. 2 and the point-probe curves in Figs. 3-5 are obtained by 'setting lp = 0 in Eqs. (20) and (21) and using Eq. (15)'. Since ℓmax would be infinite in that limit, the point-probe reference curves must have been computed with some other, undisclosed cutoff, most plausibly the healing-length cutoff of Ref. [24] where lp has a different meaning. As written, the comparison between the Plummer sphere and the point probe does not use the same cutoff prescription, so the deviations in Figs. 3-5 may be an artifact of the mismatch rather than a genuine finite-size effect. Please specify exactly which cutoff is used for the point-probe reference and demonstrate that the comparison is made consistently.","section":"Sec. 3, Eq. (15) and the paragraph before Fig. 2"},{"comment":"The angular cutoff ℓmax = π r0 / lp is inherited from the point-probe analysis of Ref. [24], where the analogous length scale is the healing length of the superfluid, not the radius of the moving body. For the Plummer sphere, the form factor ρ_Pl(k lp) in Eq. (17) already suppresses large-k contributions exponentially for k lp ≳ 1, so it is not self-evident that an additional angular truncation at ℓmax = π r0 / lp is needed or physically justified. No convergence test is provided: for a given lp/r0, one does not know whether the sums in Eq. (15) have converged before ℓmax is reached. Without such a test, the reported threshold lp/r0 ≃ 10^-2 and the deviations in Figs. 3-5 could be consequences of the arbitrary truncation rather than of the extended mass profile. I ask the authors to (i) show the convergence of Eqs. (15)-(17) as ℓmax is increased beyond π r0 / lp for several lp/r0 values, or (ii) justify the cutoff from the physics of the Plummer source, or (iii) remove the cutoff dependence from the central comparison by using a converged ℓmax for the Plummer case and an appropriate independent prescription for the point-probe case.","section":"Sec. 3, Eq. (17) and the definition of ℓmax after it"},{"comment":"The self-gravity term −4πGρ0α is neglected in the linearized equation (3), and the introduction correctly notes that this term is responsible for a non-zero dynamical friction at subsonic speeds in the linear-motion case (Ref. [38]). In the present circular-orbit results, both the radial and tangential components vanish as M → 0 (e.g., Fig. 3), which is consistent with the neglect of this term. The conclusions acknowledge homogeneity and linear-response idealizations but do not explicitly list the neglect of self-gravity. Since the abstract and conclusions make general statements about 'dynamical friction acting on circularly moving globular clusters', the omission of this physical effect and its expected impact on the subsonic regime should be stated explicitly, and the claims restricted accordingly or extended to include the self-gravity term.","section":"Sec. 2, Eq. (3) and the Introduction (discussion of Ref. [38])"}],"minor_comments":[{"comment":"The phrase 'light-hand side' should be 'left-hand side'.","section":"Sec. 2, after Eq. (3)"},{"comment":"The expressions for k1,2 appear garbled in the typeset text ('±mcsif + ml' should presumably be '±i m c_s f^+_ml' with the plus as a superscript). Also the footnote about restoring ℏ is confusing because ℏ is omitted from the main equations except in Eq. (21); please make the units explicit throughout.","section":"Eq. (19)"},{"comment":"The dimensionless force is denoted \\(\\vec{F}\\) in Eq. (14) while the total dimensional force is \\(F_{fr}\\); this reuse of notation is confusing. Please use distinct symbols, e.g., \\(\\vec{\\mathcal{F}}\\) for the dimensionless force.","section":"Sec. 3, Eq. (14)"},{"comment":"The caption does not state the angular cutoff used for the point-probe reference curves. This is directly related to the first major comment; the figure caption should specify the value or prescription of ℓmax.","section":"Caption of Fig. 2"},{"comment":"The paper calls the results 'analytic expressions' but the real part of S^ml_ℓ,ℓ−1 in Eq. (21) is evaluated numerically. Suggest using 'semi-analytic' or 'analytic up to a one-dimensional integral' to describe the force components.","section":"Abstract and Sec. 3"},{"comment":"The notation lp is used for the Plummer radius, while the same symbol denotes the healing length in Ref. [24]. This difference should be stated explicitly to avoid confusion when adopting the cutoff ℓmax from that reference.","section":"Sec. 3, after Eq. (17)"}],"recommendation":"major_revision","confidential_remarks":"The derivation from the Plummer Fourier transform through the Sokhotski-Plemelj evaluation is sound and the paper is within scope for an astrophysics journal. The main issue is the uncontrolled cutoff in the point-probe comparison, which is fixable with a convergence study and an explicit cutoff prescription. I do not see grounds for rejection, but the present version does not yet establish the headline finite-size threshold."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this paper gives the first analytic expressions for the radial and tangential dynamical friction components on a Plummer sphere moving on a circular orbit in a homogeneous ULDM superfluid. That is a real, useful extension of the point-probe framework from Berezhiani et al. The derivation is transparent: Fourier transform of the Plummer profile, linear-response equation, Sokhotski-Plemelj for the imaginary part, and the formal lp→0 limit recovers the point-probe kernel. No fitted parameters, no circularity. The limitation section is honest about the idealizations: homogeneous background, linear response, spherical symmetry, and neglect of the self-gravity term.\n\nNow the soft spots, in proportion. The stress-test concern lands. The cutoff lmax = πr0/lp in Eq. (15) diverges as lp→0, so the point-probe curves in Figs. 2–5 cannot be obtained by setting lp=0 in the same expression. They must have been computed with a different, undisclosed cutoff, presumably the healing length used in Ref. [24], where lp means something else. That means the deviations between Plummer and point-probe forces shown in Figs. 3–5 conflate the physical form factor with a change of regulator. The threshold lp/r0 ~ 10^-2 is therefore not as robust as claimed. The finite-lp formula itself is better defined, because the form factor suppresses high multipoles, but the comparison baseline is ambiguous.\n\nA second issue: the paper cites the numerical simulations by Lancaster et al. and Glennon et al. as motivation, stating they found weaker drag for Plummer spheres, but it never compares its own results to those simulations. For realistic size ratios it concludes the finite size hardly matters, which seems to be in tension with the motivation. That should be addressed.\n\nMinor notes: the real part of S is computed only numerically, which is acceptable but should be documented more carefully. The self-citation to [39] is contextual and harmless.\n\nWho is this for? Anyone modeling dynamical friction of globular clusters or extended objects in ULDM cores will want the finite-size formulas. But they should treat the point-probe comparison with caution until the cutoff question is resolved.\n\nRecommendation: send to peer review with a request for revision. The new analytic result is worth publishing, but the authors need to specify the cutoff used for the point-probe limit, check whether their deviations survive a consistent treatment, and compare directly with the numerical simulations they cite.","headline":"Useful analytic extension of point-probe dynamical friction to Plummer spheres in ULDM, but the headline comparison to the point-probe limit is muddied by an undefined cutoff and needs rework.","tokens_in":11429,"tokens_out":4300,"would_cite":true,"duration_ms":41581,"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":"For a Plummer sphere in ultralight dark matter, the dynamical friction force is the point-probe force with each momentum mode weighted by $[k\\ell_p K_1(k\\ell_p)]^2$: identical to a point mass for compact clusters, measurably different for…","keywords":["ultralight dark matter","Bose-Einstein condensate","dynamical friction","Plummer sphere","globular clusters","superfluid dark matter","linear response","orbital decay"],"falsifier":"Run a fully nonlinear numerical simulation of the ultralight scalar field with self-gravity for a Plummer-sphere source with radius one-tenth of the orbital radius on a circular orbit at Mach number 2, and compare the measured radial and tangential drag forces with Eqs. (14)-(17); agreement would support the form-factor claim, while a mismatch at these parameters would falsify it.","tokens_in":10330,"feed_emoji":"🌌","tokens_out":9086,"duration_ms":77265,"temperature":0.7,"pith_summary":"The paper derives analytic formulas for the radial and tangential components of the dynamical friction force on a globular cluster modelled as a Plummer sphere moving on a circular orbit through a homogeneous superfluid of ultralight bosonic dark matter. The central result is that the drag equals the known point-probe force with the squared momentum-space density profile of the sphere inserted into the momentum integrals, Eqs. (14) with (15) and (17). When the Plummer radius is below about one hundredth of the orbital radius, the extended body feels essentially the same force as a point mass of the same mass; at larger radii, or at large Mach numbers, the two components deviate visibly. The paper also finds that a spherically symmetric cluster receives no internal torque from the drag, so dynamical friction does not spin it up. The formulas matter because globular clusters are usually treated as point particles in dark-matter halo studies, and these equations state when that treatment is safe.","feed_headline":"Same dark-matter drag for compact globular clusters, different for large ones","feed_subtitle":"New analytic radial and tangential friction formulas for Plummer spheres show when treating a cluster as a point mass fails.","key_machinery":"The load-bearing object is the momentum-space form factor of the Plummer sphere, $\\rho_{\\mathrm{Pl}}(k\\ell_p)=M\\,k\\ell_p\\,K_1(k\\ell_p)$. In the linear response equation the moving sphere enters only through this factor, so the entire finite-size effect on the force is packaged into the squared ratio $[k\\ell_p K_1(k\\ell_p)]^2$ multiplying the point-probe integrand. The spherical-harmonic decomposition with the cutoff $\\ell_{\\max}=\\pi r_0/\\ell_p$ and the real positive zero $k_3$ of the dispersion relation turn the force into the two-component formulas of Eqs. (15)-(21).","core_discovery":"On the paper's own terms, the discovery is a structural extension of the point-probe result: for a circularly moving Plummer sphere the total dynamical friction force is obtained from the point-probe expression by the replacement $\\rho_p(k)\\to\\rho_{\\mathrm{Pl}}(k\\ell_p)$, where $\\rho_{\\mathrm{Pl}}(k\\ell_p)=M\\,k\\ell_p\\,K_1(k\\ell_p)$ and $K_1$ is the modified Bessel function of the second kind. The only modification is the factor $[\\rho_{\\mathrm{Pl}}(k\\ell_p)/M]^2=[k\\ell_p K_1(k\\ell_p)]^2$ inside the momentum integral. Evaluating the integral with a distribution identity for the imaginary part gives a closed formula for the tangential component, Eq. (20), and a principal-value integral for the radial component, Eq. (21), with the harmonic sum cut off at $\\ell_{\\max}=\\pi r_0/\\ell_p$. The quantitative content is that the point-probe approximation holds for $\\ell_p/r_0\\lesssim 10^{-2}$ and fails at larger radii or large Mach numbers, precisely where the extended density suppresses modes with $k\\sim 1/\\ell_p$.","pith_inferences":["The same squared-form-factor replacement should extend the point-probe formulas to other spherically symmetric mass models, so the momentum-space density ratio is a general handle for finite-size corrections to dynamical friction.","Because the suppression sets in at $k\\sim 1/\\ell_p$, large globular clusters and dwarf galaxies in ultralight-dark-matter halos could sink more slowly than point-probe estimates suggest; the paper illustrates this but does not develop the orbital-decay consequence.","The vanishing inner torque is a symmetry test: a tidally deformed, aspherical cluster would acquire spin through dynamical friction, so measuring cluster rotation could simultaneously probe shape and dark-matter response.","A direct numerical simulation of the fully nonlinear ultralight scalar field with self-gravity, for $\\ell_p/r_0=0.1$ at $M=2$, would test whether the linear homogeneous form-factor formula survives nonlinear wake effects."],"forward_implications":["For globular clusters with $\\ell_p/r_0\\lesssim 10^{-2}$, the point-mass idealization reproduces both components of the drag, so existing decay-time estimates built on point-probe formulas remain valid in that regime.","For clusters with radius a tenth of the orbital radius or larger, and at Mach numbers near the force maximum, the finite size changes both components, so orbital-decay modelling must use the extended-body expressions.","The total dynamical-friction torque on a spherical cluster is simply $\\mathbf{r}_{\\mathrm{CM}}\\times\\mathbf{F}_{\\mathrm{fr}}$; a spherical Plummer sphere does not spin up from dynamical friction.","The radial and tangential components depend on orbital radius in qualitatively different ways in the sampled cases, so fitting both components separately carries more information about the dark-matter parameters than the total drag alone."],"supporting_citations":[{"why":"Supplies the point-probe dynamical friction result in dark-matter superfluids whose momentum-space integral the Plummer form factor multiplies; the paper's Eq. (14) is built directly on this result.","marker":"[24]"},{"why":"Provides the analytic framework for the dynamical friction force on circularly moving perturbers that the paper follows for the geometry and the response equations.","marker":"[10]"},{"why":"Numerical fuzzy-dark-matter simulations establishing the linear response regime and showing Plummer spheres experience weaker friction than point masses, motivating the analytic calculation.","marker":"[17]"},{"why":"Derives the ultralight-scalar density-response equations used to set up the linear perturbation problem for a moving perturber.","marker":"[16]"},{"why":"Extends the analytic circular-orbit treatment to fuzzy dark matter with quantum pressure, providing the spherical-harmonic expansion and integration procedure.","marker":"[23]"},{"why":"Identifies the self-gravity term that the paper neglects, defining the boundary of validity of the approximation.","marker":"[38]"},{"why":"Defines the Plummer sphere density profile used to model the globular cluster.","marker":"[35]"}],"fun_headline_variants":["Plummer sphere drag in ultralight dark matter deviates from point mass","New analytic friction formulas for globular clusters in ultralight dark matter","Cluster size affects dynamical friction in ultralight dark matter","Point mass fails for large globular clusters in ultralight dark matter","Dynamical friction on Plummer spheres: point mass fails at large radii"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation assumes the dark-matter response stays linear and spatially uniform, with the gravitational back-reaction of the disturbed dark matter (the term $-4\\pi G\\rho_0\\alpha$) neglected, so the cluster's size enters only as a momentum-space filter; if nonlinear wakes, halo granularity, or that back-reaction matter for real globular clusters, the given force formulas do not apply.","fun_headline_variants_meta":{"raw":{"variants":["Plummer sphere drag in ultralight dark matter deviates from point mass","New analytic friction formulas for globular clusters in ultralight dark matter","Cluster size affects dynamical friction in ultralight dark matter","Point mass fails for large globular clusters in ultralight dark matter","Dynamical friction on Plummer spheres: point mass fails at large radii"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000705,"raw_usage":{"total_tokens":3153,"prompt_tokens":893,"completion_tokens":2260,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":509,"completion_tokens_details":{"reasoning_tokens":2165}},"tokens_in":509,"tokens_out":2260,"duration_ms":13735,"temperature":1.0,"reasoning_tokens":2165,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T11:26:26.874741+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a fully nonlinear numerical simulation of the ultralight scalar field with self-gravity for a Plummer-sphere source with radius one-tenth of the orbital radius on a circular orbit at Mach number 2, and compare the measured radial and tangential drag forces with Eqs. (14)-(17); agreement would support the form-factor claim, while a mismatch at these parameters would falsify it.","supporting_citations":[{"cited_title":"Desjacques, A","cited_arxiv_id":null,"evidence_quote":"Provides the analytic framework for the dynamical friction force on circularly moving perturbers that the paper follows for the geometry and the response equations."},{"cited_title":"Plummer, On the problem of distribution in globular star clusters , Mon","cited_arxiv_id":null,"evidence_quote":"Defines the Plummer sphere density profile used to model the globular cluster."}],"review_version":1}