Pith. sign in

REVIEW 3 major objections 4 minor 48 references

A Schwarzschild black hole's absorption of a massive scalar packet is governed by the ratio of its horizon radius to the field's reduced Compton wavelength, with a sharp behavioral change when the two are equal.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-01 13:59 UTC pith:GNSEMCAA

load-bearing objection Honest, reproducible numerical maps of mode-by-mode Noether-charge accretion; the low-k0 floor in the Rs>lambdaC regime is real for their data, but its interpretation as mass-dominated slow transport is not established because the initial data is not an on-shell massive ingoing packet. the 3 major comments →

arxiv 2607.18920 v1 pith:GNSEMCAA submitted 2026-07-21 gr-qc astro-ph.GA

Accretion of multipolar massive complex scalar field packets by a Schwarzschild black hole

Flavio Rosales-Infante , Ivan Alvarez-Rios , Francisco S. Guzman This is my paper
classification gr-qc astro-ph.GA MSC 83C57 PACS 04.70.-s04.25.dg
keywords massive scalar fieldsSchwarzschild black holeNoether charge accretionmultipolar decompositioneffective potential barrierCompton wavelengthfuzzy dark matterKlein-Gordon equation
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper studies a localized, nearly monochromatic packet of a complex massive scalar field falling onto a Schwarzschild black hole, and asks which multipole components are absorbed and which remain outside. It answers with accretion-efficiency maps in the (k0, ℓ) plane, where each pixel is the fraction of that mode's conserved Noether charge that crosses the horizon by a fixed final time. The paper's central claim is that these maps are organized by the ratio μ̃ = Rs/λC: when the horizon is smaller than the Compton wavelength, a broad transition follows the classical effective-potential threshold; when the horizon is larger, the transition sharpens and a partial-accretion floor appears at low k0. A sympathetic reader would care because this provides a time-domain, charge-based classification of what a black hole accretes from a fuzzy, wave-like dark matter distribution, and which high-multipole remnant it leaves behind.

Core claim

The central discovery is that mode-resolved accretion of massive scalar wave packets separates into regimes determined by Rs/λC, with the transition from inefficient to efficient accretion tracked by the classical transmission threshold k0^tr = sqrt(Vmax − μ̃²) for Rs ≲ λC, and a distinct low-k0 partial-accretion floor for Rs > λC. The paper shows this with time-domain evolutions of the Klein–Gordon equation decomposed into independent multipoles, measuring the flux of the Noether current through the horizon. Its concluding statement is that this change across Rs = λC is the main result encoded in the parameter-space maps.

What carries the argument

The argument rests on the massive-scalar effective potential in tortoise coordinates, V_eff(r; ℓ, μ̃) = (1 − Rs/r)(μ̃² + ℓ(ℓ+1)/r² + Rs/r³), whose maximum defines the classical transmission threshold via the on-shell dispersion ω² = k0² + μ̃². Around this barrier the paper organizes the mode-by-mode accretion fractions, computed as accumulated Noether-current flux through the horizon from 1+1D multipolar evolutions, into two-dimensional response maps.

Load-bearing premise

The maps are interpreted as if each packet has the massive on-shell frequency ω² = k0² + μ̃² and enough time to reach the horizon, but the initial data sets the momentum with a massless-like relation and the evolution runs only for 200 Rs/c; when μ̃ > k0, which is exactly the regime where the new floor appears, both assumptions are strained.

What would settle it

Compute the same (k0, ℓ) accretion maps with the initial momentum set from the massive on-shell dispersion (πℓ built from ω = sqrt(k0² + μ̃²)) and with evolution time extended to roughly 10^4 Rs/c. If the low-k0 partial-accretion floor for Rs > λC and the sharpened transition persist, the paper's central claim is supported; if the floor disappears or the transition reverts to the single classical threshold, the claimed regime change is an artifact of the chosen initial data and finite time window.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • If the central claim is correct, the same effective-potential landscape that supports long-lived massive-scalar configurations also decides which multipoles survive a close passage: low-ℓ sectors are absorbed, high-ℓ sectors are scattered back.
  • For Rs ≲ λC, the boson mass barely changes the response maps; the angular barrier is the controlling factor.
  • For Rs > λC, low-k0 packets no longer saturate to complete accretion even at low ℓ, so the hole becomes a selective filter rather than a total absorber.
  • The integrated accretion efficiency over representative sectors increases monotonically with k0 and rises steeply when many multipoles cross the threshold together.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Because the evolution time is 200 Rs/c for packets launched at 80 Rs, and on-shell massive packets with k0 ≲ μ̃ arrive later, part of the low-k0 floor may be an arrival-time cutoff. Extending the runs to thousands of Rs/c would separate genuine scattering from packets that simply have not reached the horizon.
  • The initial momentum is set by πℓ = sqrt(γrr)ψℓ (Eqs. 17–18), a massless-like choice, while the interpretation assumes ω² = k0² + μ̃². Reinitializing with the massive on-shell frequency is a direct test of whether the Rs > λC transition persists.
  • If mode filtering is real, fuzzy-dark-matter remnants near supermassive black holes should be systematically depleted of low multipoles, so scalar structures near the horizon would be dominated by high-ℓ angular structure.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper studies the finite-time accretion of complex massive scalar wave packets by a Schwarzschild black hole in the test-field limit. After a spherical-harmonic decomposition, the authors evolve nearly monochromatic Gaussian packets mode by mode in a 1+1 first-order formulation on an Eddington–Finkelstein background, and use the conserved Noether current flux through the horizon as an accretion diagnostic. For a grid of radial wavenumbers k0 and multipoles ℓ, they construct accretion-efficiency maps, find a transition between inefficient, partial, and efficient accretion, and compare the transition with the classical effective-potential threshold. Their central claim is that the response is controlled by μ~=Rs/λC: for Rs≲λC the transition is broad and angular-barrier dominated, while for Rs>λC it sharpens and a partial-accretion floor appears at low k0. They interpret this floor as mass-dominated slow radial transport and summarize this as the main new result.

Significance. If established, the paper would provide a useful time-domain, Noether-charge-based classification of massive-scalar absorption by black holes, complementing the stationary scattering framework. The numerical core is a real strength: second-order self-convergence, a vanishing conservation defect under refinement, an inner-boundary-location test, and open data are reported. The effective-potential threshold Eq. (21) contains no fitted parameters. However, the main new feature — the low-k0 partial-accretion floor in the Rs>λC regime — is not yet established, because the initial data actually evolved are not the on-shell massive packets assumed in the interpretation, and the finite evolution time is too short to distinguish genuine reflection from arrival-time truncation. The maps remain valid numerical results for the specified initial-data family, but the physical interpretation of the central new feature requires additional work.

major comments (3)
  1. [Sec. V-A, Eqs. (16)-(18); Sec. V-C, Eq. (21)] The initial condition πℓ=√γ^rr ψℓ selects a massless-like frequency content (ω≈k0) in the approximately flat initial region, not the on-shell massive relation ω=√(k0²+μ~²) used to define the threshold Eq. (21) and used in the interpretation of the floor. For a purely ingoing massive positive-frequency packet with carrier k0 one would need π≈(ω/k0)ψ; π=ψ corresponds to ω≈k0. In the regime μ~≫k0 the prescribed data contains a substantial outgoing component — for μ~=10k0 the relative amplitude is ≈(ω−k0)/(ω+k0)≈0.82. The paper's own check in Sec. V-A, Qℓ(0)=μ~²+ℓ(ℓ+1)/r², follows from Eq. (14) only if ω²≈γ^rr k²≈k0², confirming the off-shell nature of the data. The low-k0 floor in Fig. 7 may therefore be an artifact of this mixed initial data rather than mass-dominated slow transport. Please evolve genuinely massive ingoing data, or provide a Fourier decomposition of the initial data and a
  2. [Sec. V-A (r0=80Rs, σ=10Rs, T=200Rs/c) and Sec. V-C (ηℓ,acc final values)] For an on-shell massive packet with small k0 and μ~>1, the group velocity is k0/√(k0²+μ~²), so the packet moves much more slowly than c. For example, μ~=10 and k0=0.1 gives a travel time from r0=80Rs to the horizon of order 8000Rs/c, far exceeding T=200Rs/c. The statement that T is 'long enough to capture the accretion stage' is therefore not justified for the low-k0 floor points. The final-time values of ηℓ,acc may not be saturated, and the floor could be an arrival-time truncation effect. Please show time series of ηℓ,acc(t) for representative low-k0 floor points, or otherwise demonstrate that the floor is a final-state property rather than a finite-integration-time artifact.
  3. [Sec. V-C, Eq. (21) and Figs. 6-7] The classical threshold curve ktr0=√(Vmax−μ~²) is derived from the stationary effective potential using the asymptotic on-shell dispersion relation. Since the initial data in Sec. V-A are not on-shell massive modes, the comparison of the maps against this curve is not fully coherent in the Rs>λC regime where μ~ is large and the mismatch is largest. For Rs≲λC the threshold is a meaningful benchmark, but the sharp-transition and floor claims concern exactly the regime where the initial data are far off-shell. A threshold based on the actual frequency content of the initial data, or a modified initial data set that is on-shell, would be needed to support the claimed interpretation.
minor comments (4)
  1. [Sec. V-C] Typographical issues: 'dependient' should be 'dependent', 'bsoson' should be 'boson', and 'counterpat' should be 'counterpart'.
  2. [Sec. IV] The sign convention for the shift βr and the identification of the characteristic fields e± with 'ingoing'/'outgoing' directions is not stated precisely; clarifying this would make the claim that Eq. (18) selects 'ingoing' initial data unambiguous.
  3. [References] References [23] and [38] are the same paper, and [25] and [41] are also the same paper; please consolidate to avoid duplication.
  4. [Sec. V-A] The sentence 'This relation shows that ingoing data are initially compatible with the modal Hamilton–Jacobi balance equation' is misleading: as noted above, the relation actually demonstrates that the initial data satisfy a massless-like dispersion balance, not the massive on-shell relation used later. Please rephrase or correct.

Circularity Check

0 steps flagged

No significant circularity: the threshold curve is an external benchmark and the maps are direct initial-value outputs.

full rationale

The derivation chain is not circular. The accretion maps in Figs. 6–7 are direct numerical outputs of the initial-value problem defined by Eqs. (13) and (16)–(18); no parameter in the map construction is fitted to the threshold curve. The only derived benchmark, Eq. (21), is computed from the standard massive-scalar effective potential with no fitted constants and is used as an external comparison, not as an input that generates the maps. The statement that μ̃ = Rs/λC controls the response is dimensional: after the rescaling leading to Eq. (10), μ̃ is the only dimensionless parameter, so this is a scaling statement rather than a self-derived result. Self-citations such as Refs. [10–12, 42] are motivational and contextual; they are not the load-bearing justification for the central numerical results, and the validation in Appendix A is self-contained. The skeptic’s concern that Eq. (18) does not produce a purely ingoing on-shell massive packet, and that the finite window T = 200 Rs/c may truncate slow low-k0 transport, is a physical-interpretation and correctness risk, not a circular reduction: the low-k0 floor is an output of the specified Cauchy data, and no equation in the paper is constructed from the quantity it is used to explain.

Axiom & Free-Parameter Ledger

5 free parameters · 5 axioms · 0 invented entities

The central claims rest on the standard KG-in-Schwarzschild framework plus a set of modeling choices. The only physical input is μ̃ = Rs/λC; the numerical parameters (r0, σ, T, ℓmax, contour levels) are hand-chosen and shape the quantitative maps. The main unexamined assumption is that the initial packet is well described by the massive dispersion relation, which its actual momentum content does not satisfy.

free parameters (5)
  • initial packet placement r0 = 80 Rs
    Chosen to be far enough to be quasi-Minkowskian but inside the grid; affects travel time and therefore the finite-time accretion fractions.
  • initial packet width σ = 10 Rs
    Sets the spectral width of the packet; affects how 'nearly monochromatic' the carrier is.
  • final integration time T = 200 Rs/c
    Determines the accretion fractions; the low-k0 'partial-accretion floor' for Rs > λC is partly a finite-time effect since slow packets need ≳1000 Rs/c to travel from r0 to the horizon.
  • regime classification thresholds = η = 0.01, 0.5, 0.99 contours
    The inefficient/partial/efficient classification boundaries are arbitrary contour levels, not derived quantities.
  • multipole truncation ℓmax = 80
    Truncation of the modal sum; the integrated efficiency ηacc depends on it.
axioms (5)
  • domain assumption Test-field regime: the scalar field does not backreact on the spacetime
    The entire study fixes a Schwarzschild background (Sec. II-D); accretion diagnostics ignore the self-gravity of the field.
  • domain assumption Asymptotic dispersion relation ω² = k0² + μ̃² applies to the initial packet
    Used to define the transmission threshold ktr0 in Eq. (21) and to interpret the partial-accretion floor (Sec. V-C). The initial data (17)–(18) actually fixes π ≈ +ψ, i.e., ω ≈ k0 near flat space, so the on-shell massive relation is not satisfied by the initial data when μ̃ > k0.
  • ad hoc to paper Initial data is 'ingoing' and compatible with the modal Hamilton–Jacobi balance
    Sec. V-A claims Qℓ(t=0) = μ̃² + ℓ(ℓ+1)/r²; direct evaluation of Eq. (14) with S = −ℏk0r gives an additional k0²(r+2Rs)/(r+Rs) term, so the stated compatibility appears to be an approximation or error.
  • standard math m-degeneracy of the reduced equations
    Eq. (13) has no m dependence, so m=0 represents each ℓ; standard for spherically symmetric backgrounds.
  • domain assumption Classical barrier transmission condition ω² > Vmax organizes the maps
    The effective potential Veff and threshold (21) are the standard stationary-scattering picture; its agreement with the maps is asserted qualitatively for the Rs ≤ λC regime.

pith-pipeline@v1.3.0-alltime-deepseek · 19253 in / 33179 out tokens · 299324 ms · 2026-08-01T13:59:17.439113+00:00 · methodology

0 comments
read the original abstract

We study the finite-time accretion of complex massive scalar wave packets by a Schwarzschild black hole in the test-field regime, with parameters motivated by ultralight fuzzy dark matter around supermassive black holes. Our goal is to determine how the scalar content of a localized configuration is redistributed after interacting with the black hole, and which spectral and multipolar components are more efficiently absorbed. We decompose the Klein--Gordon field into independent multipolar sectors and evolve nearly monochromatic Gaussian packets mode by mode, reducing the problem to a set of 1+1 dimensional evolutions. Accretion is quantified with the flux of the conserved Noether current through the horizon surface, providing a direct measure of the scalar charge absorbed by the black hole. For a carrier radial wavenumber $k_0$ and multipole index $\ell$, we construct accretion-efficiency maps in the $(k_0,\ell)$ plane that contain the fraction of accreted modal charge. These maps exhibit a transition between inefficient, partial, and efficient accretion regimes, which we relate to the structure of an effective potential. We show that the process is controlled by the ratio between the Schwarzschild radius $R_s$ and the reduced Compton wavelength $\lambdabar_C$. For $R_s \lesssim \lambdabar_C$, the transition is broad and dominated by the angular momentum barrier, while for $R_s > \lambdabar_C$ it sharpens across a narrower range of $k_0$ and a partial-accretion floor emerges at low $k_0$. These results provide a time-domain, Noether-charge-based classification of black hole accretion for massive scalar wave packets.

Figures

Figures reproduced from arXiv: 2607.18920 by Flavio Rosales-Infante, Francisco S. Guzman, Ivan Alvarez-Rios.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6 [PITH_FULL_IMAGE:figures/full_fig_p012_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7 [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8 [PITH_FULL_IMAGE:figures/full_fig_p014_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: shows the resulting time-dependent exponent p(t). After a short initial transient, p(t) remains close to p ' 2 over most of the evolution, with moderate time variability associated with the wave-packet dynam￾ics. This behavior is consistent with the nominal second￾order accuracy of the radial finite difference discretiza￾tion, that dominates over the third order of the evolu- [PITH_FULL_IMAGE:figures/full… view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10 [PITH_FULL_IMAGE:figures/full_fig_p018_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11 [PITH_FULL_IMAGE:figures/full_fig_p019_11.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

48 extracted references · 1 canonical work pages

  1. [1]

    Scattering of scalar waves from a Schwarzschild black hole,

    Norma G. Sanchez, “Scattering of scalar waves from a Schwarzschild black hole,” Journal of Mathematical Physics 17, 688–692 (1976)

    1976

  2. [2]

    It represents the cumulative contribution of many modes with moderate accretion efficiency. Thus, the summed curve should be read together with the maps, that is, a nonzero value of ηacc at small k0 indicates par- tial charge absorption across several representative sec- tors, while the approach to ηacc ' 1 indicates that al- most all evolved sectors carry...

    2025

  3. [3]

    Absorption cross section of small black holes,

    W. G. Unruh, “Absorption cross section of small black holes,” Physical Review D14, 3251–3259 (1976)

    1976

  4. [4]

    Scattering and absorption of electromagnetic waves by a Schwarzschild black hole,

    R. Fabbri, “Scattering and absorption of electromagnetic waves by a Schwarzschild black hole,”Physical Review D 12, 933–942 (1975)

    1975

  5. [5]

    Total transmission through the Schwarzschild black-hole potential barrier,

    Nils Andersson, “Total transmission through the Schwarzschild black-hole potential barrier,”Classical and Quantum Gravity11, L39–L44 (1994)

    1994

  6. [6]

    Universality of low energy absorption cross sections for black holes,

    Sumit R. Das, Gary W. Gibbons, and Samir D. Mathur, “Universality of low energy absorption cross sections for black holes,” Physical Review Letters 78, 417–419 (1997), arXiv:hep-th/9609052

    Pith/arXiv arXiv 1997

  7. [7]

    In these maps, the color scale represents ηℓ,acc(T ;k0) 2 [0, 1]. Values close to zero indicate that most of the modal charge remains in the exterior after the interaction, val- ues close to unity indicate efficient accretion, and inter- mediate values correspond to partial accretion. In this way, the maps provide a response diagram for the black hole, the...

  8. [8]

    Absorption of massive scalar field by a charged black hole,

    T. Nakamura and H. Sato, “Absorption of massive scalar field by a charged black hole,”Physics Letters B61, 371– 374 (1976)

    1976

  9. [9]

    Absorptionofamassive scalar field by a charged black hole,

    Carolina L. Benone, Ednilton S. de Oliveira, Sam R. Dolan, andLuísC.B.Crispino,“Absorptionofamassive scalar field by a charged black hole,”Physical Review D 89, 104053 (2014), arXiv:1404.0687 [gr-qc]

    Pith/arXiv arXiv 2014

  10. [10]

    Scalar absorption by charged rotat- ing black holes,

    Luiz C. S. Leite, Carolina L. Benone, and Luís C. B. Crispino, “Scalar absorption by charged rotat- ing black holes,”Physical Review D 96, 044043 (2017), arXiv:1708.03370 [gr-qc]

    Pith/arXiv arXiv 2017

  11. [11]

    Scalar absorption: Black holes ver- sus wormholes,

    Haroldo C. D. Lima Junior, Carolina L. Benone, and Luís C. B. Crispino, “Scalar absorption: Black holes ver- sus wormholes,”Physical Review D101, 124009 (2020), arXiv:2006.03967 [gr-qc]

    Pith/arXiv arXiv 2020

  12. [12]

    Scalar field dark matter: behavior around black holes,

    A. Cruz-Osorio, F. S. Guzmán, and F. D. Lora-Clavijo, “Scalar field dark matter: behavior around black holes,” Journal of Cosmology and Astroparticle Physics2011, 029 (2011)

    2011

  13. [13]

    Spherical nonlin- ear absorption of cosmological scalar fields onto a black hole,

    F. S. Guzmán and F. D. Lora-Clavijo, “Spherical nonlin- ear absorption of cosmological scalar fields onto a black hole,” Phys. Rev. D85, 024036 (2012)

    2012

  14. [14]

    Fre- quency shift of light emitted from growing and shrinking black holes,

    F.G.Guzmán, I.Alvarez-Ríos, andJ.A.González,“Fre- quency shift of light emitted from growing and shrinking black holes,”Phys. Rev. D104, 084014 (2021)

    2021

  15. [15]

    A Fur- ther analysis of a cosmological model of quintessence and scalar dark matter,

    Tonatiuh Matos and L. Arturo Urena-Lopez, “A Fur- ther analysis of a cosmological model of quintessence and scalar dark matter,”Phys. Rev. D63, 063506 (2001)

    2001

  16. [16]

    Ultralight scalars as cosmological dark matter,

    Lam Hui, Jeremiah P. Ostriker, Scott Tremaine, and Edward Witten, “Ultralight scalars as cosmological dark matter,” Phys. Rev. D95, 043541 (2017)

    2017

  17. [17]

    Small-scale structure of fuzzy and axion-like dark matter,

    Jens C. Niemeyer, “Small-scale structure of fuzzy and axion-like dark matter,”Progress in Particle and Nuclear Physics 113, 103787 (2020)

    2020

  18. [18]

    Wald, General Relativity (The University of Chicago Press, Chicago, 1984)

    Robert M. Wald, General Relativity (The University of Chicago Press, Chicago, 1984)

    1984

  19. [19]

    N. D. Birrell and P. C. W. Davies,Quantum Fields in Curved Space (Cambridge University Press, Cambridge, 1982)

    1982

  20. [20]

    Miguel Alcubierre, Introduction to 3+1 Numerical Rela- tivity (Oxford University Press, Oxford, 2008)

    2008

  21. [21]

    Peskin and Daniel V

    Michael E. Peskin and Daniel V. Schroeder,An Introduc- tion to Quantum Field Theory (Addison-Wesley, Read- ing, Massachusetts, 1995)

    1995

  22. [22]

    Quantentheorie in hydrodynamischer form,

    Erwin Madelung, “Quantentheorie in hydrodynamischer form,” Zeitschrift für Physik40, 322–326 (1926)

    1926

  23. [23]

    Asuggestedinterpretationofthequantum theory in terms of “hidden

    DavidBohm,“Asuggestedinterpretationofthequantum theory in terms of “hidden” variables.i,”Physical Review 85, 166–179 (1952)

    1952

  24. [24]

    Klein-gordon equation in hydrody- namical form,

    Cheuk-Yin Wong, “Klein-gordon equation in hydrody- namical form,” Journal of Mathematical Physics 51, 122304 (2010)

    2010

  25. [25]

    Are black holes a serious threat to scalar field dark matter models?

    Juan Barranco, Argelia Bernal, Juan Carlos Degollado, Alberto Diez-Tejedor, Miguel Megevand, Miguel Alcu- bierre, Darío Núñez, and Olivier Sarbach, “Are black holes a serious threat to scalar field dark matter models?” Physical Review D 84, 083008 (2011), arXiv:1108.0931 [gr-qc]

    Pith/arXiv arXiv 2011

  26. [26]

    Horizon growth of supermassive black hole seeds fed with collisional dark matter,

    F. D. Lora-Clavijo, M. Gracia-Linares, and F. S. Guz- man, “Horizon growth of supermassive black hole seeds fed with collisional dark matter,” Monthly Notices of the Royal Astronomical Society443, 2242–2251 (2014), arXiv:1406.7233 [astro-ph.GA]

    Pith/arXiv arXiv 2014

  27. [27]

    Qua- sistationary solutions of scalar fields around accreting black holes,

    Nicolas Sanchis-Gual, Juan Carlos Degollado, Paula Izquierdo, José A. Font, and Pedro J. Montero, “Qua- sistationary solutions of scalar fields around accreting black holes,” Physical Review D 94, 043004 (2016), arXiv:1606.05146 [gr-qc]

    Pith/arXiv arXiv 2016

  28. [28]

    Grav- itational atoms beyond the test field limit: the case of sgr a* and ultralight dark matter,

    Miguel Alcubierre, Juan Barranco, Argelia Bernal, Juan Carlos Degollado, Alberto Diez-Tejedor, Miguel Megevand, Darío Núñez, and Olivier Sarbach, “Grav- itational atoms beyond the test field limit: the case of sgr a* and ultralight dark matter,”Classical and Quan- tum Gravity 42, 21LT01 (2025)

    2025

  29. [29]

    First m87 event horizon telescope results. vi. the shadow and mass ofthecentralblackhole,

    The Event Horizon Telescope Collaboration, “First m87 event horizon telescope results. vi. the shadow and mass ofthecentralblackhole,” TheAstrophysicalJournalLet- ters 875, L6 (2019), arXiv:1906.11243 [astro-ph.GA]

    Pith/arXiv arXiv 2019

  30. [30]

    First constraints on Fuzzy Dark Matter from Lyman-α Forest data and hy- drodynamical simulations,

    Vid Iršič, Matteo Viel, Martin G. Haehnelt, James S. Bolton, and George D. Becker, “First constraints on Fuzzy Dark Matter from Lyman-α Forest data and hy- drodynamical simulations,”Physical Review Letters119 (2017), 10.1103/physrevlett.119.031302

    doi:10.1103/physrevlett.119.031302 2017

  31. [31]

    Constraining the mass of light bosonic dark matter using SDSS Lyman-α forest,

    Eric Armengaud, Nathalie Palanque-Delabrouille, Christophe Yèche, David J. E. Marsh, and Julien Baur, “Constraining the mass of light bosonic dark matter using SDSS Lyman-α forest,” Monthly Notices of the Royal Astronomical Society471, 4606–4614 (2017). 16

    2017

  32. [32]

    A search for ultralight axions using pre- cision cosmological data,

    Renée Hlozek, Daniel Grin, David J. E. Marsh, and Pe- dro G. Ferreira, “A search for ultralight axions using pre- cision cosmological data,”Physical Review D91 (2015), 10.1103/physrevd.91.103512

    doi:10.1103/physrevd.91.103512 2015

  33. [33]

    Galaxy uv- luminosity function and reionization constraints on axion dark matter,

    Brandon Bozek, David J. E. Marsh, Joseph Silk, and Rosemary F. G. Wyse, “Galaxy uv- luminosity function and reionization constraints on axion dark matter,” Monthly Notices of the Royal Astronomical Society 450, 209–222 (2015), https://academic.oup.com/mnras/article- pdf/450/1/209/18503963/stv624.pdf

    2015

  34. [34]

    Theeffects of the small-scale dm power on the cosmological neutral hydrogen (hi) distribution at high redshifts,

    Abir Sarkar, Rajesh Mondal, Subinoy Das, Shiv.K. Sethi, SomnathBharadwaj, andDavidJ.E.Marsh,“Theeffects of the small-scale dm power on the cosmological neutral hydrogen (hi) distribution at high redshifts,”Journal of Cosmology and Astroparticle Physics2016, 012 (2016)

    2016

  35. [35]

    Angular mo- mentumand vortexformation in bose-einstein-condensed cold dark matter haloes: Angular momentum in bec-cdm haloes,

    Tanja Rindler-Daller and Paul R. Shapiro, “Angular mo- mentumand vortexformation in bose-einstein-condensed cold dark matter haloes: Angular momentum in bec-cdm haloes,” Monthly Notices of the Royal Astronomical So- ciety 422, 135–161 (2012)

    2012

  36. [36]

    Exact solution to finite temperature sfdm: Natural cores without feedback,

    Victor H. Robles and T. Matos, “Exact solution to finite temperature sfdm: Natural cores without feedback,”The Astrophysical Journal 763, 19 (2012)

    2012

  37. [37]

    A fast spherical har- monics transform algorithm,

    Reiji Suda and Masayasu Takami, “A fast spherical har- monics transform algorithm,”Math. Comput. 71, 703– 715 (2002)

    2002

  38. [38]

    Black hole excision for dynamic black holes,

    Miguel Alcubierre, Bernd Brügmann, Denis Pollney, Ed- ward Seidel, and Ryoji Takahashi, “Black hole excision for dynamic black holes,”Phys. Rev. D 64, 061501(R) (2001)

    2001

  39. [39]

    thesis, University of Wales, Cardiff, Cardiff, UK (1994)

    Miguel Alcubierre, Investigations in Numerical Relativ- ity, Ph.D. thesis, University of Wales, Cardiff, Cardiff, UK (1994)

    1994

  40. [40]

    Are black holes a serious threat to scalar field dark matter models?

    Juan Barranco, Argelia Bernal, Juan Carlos Degol- lado, Alberto Diez-Tejedor, Miguel Megevand, Miguel Alcubierre, Darío Núñez, and Olivier Sarbach, “Are black holes a serious threat to scalar field dark matter models?” Physical Review D 84 (2011), 10.1103/phys- revd.84.083008

    doi:10.1103/phys- 2011

  41. [41]

    Schwarzschild black holes can wear scalar wigs,

    Juan Barranco, Argelia Bernal, Juan Carlos Degol- lado, Alberto Diez-Tejedor, Miguel Megevand, Miguel Alcubierre, Darío Núñez, and Olivier Sarbach, “Schwarzschild black holes can wear scalar wigs,”Physi- cal Review Letters109, 081102 (2012), arXiv:1207.2153 [gr-qc]

    Pith/arXiv arXiv 2012

  42. [42]

    Schwarzschild scalar wigs: Spectral analysis and late time behavior,

    Juan Barranco, Argelia Bernal, Juan Carlos Degol- lado, Alberto Diez-Tejedor, Miguel Megevand, Miguel Alcubierre, Darío Núñez, and Olivier Sarbach, “Schwarzschild scalar wigs: Spectral analysis and late time behavior,” Physical Review D 89, 083006 (2014), arXiv:1312.5808 [gr-qc]

    Pith/arXiv arXiv 2014

  43. [43]

    Quasis- tationary solutions of scalar fields around accreting black holes,

    Nicolas Sanchis-Gual, Juan Carlos Degollado, Paula Izquierdo, José A. Font, and Pedro J. Montero, “Quasis- tationary solutions of scalar fields around accreting black holes,” Physical Review D94, 043004 (2016)

    2016

  44. [44]

    Black holes as condensation points of fuzzy dark matter cores,

    Curicaveri Palomares-Chavez, Ivan Alvarez-Rios, and Francisco S. Guzman, “Black holes as condensation points of fuzzy dark matter cores,” Physical Review D 112, 023545 (2025),arXiv:2412.15465 [astro-ph.GA]

    Pith/arXiv arXiv 2025

  45. [45]

    Dataset for

    Flavio Rosales-Infante, Francisco S. Guzman, and Iván Álvarez Rios,“Dataset for ”accretion of multipolar mas- sive complex scalar field packets by a schwarzschild black hole”,” (2026), 10.5281/zenodo.20130371. 0 25 50 75 100 125 150 175 200 t 1.95 2.00 2.05 2.10 2.15 2.20 p ( t ) Time evolution of the Richardson auto-convergence exponent FIG. 9. Time evol...

    work page doi:10.5281/zenodo.20130371 2026

  46. [46]

    Procedure for estimation and reporting of uncertainty due to discretization in cfd applications,

    Ismail B. Celik, Urs Ghia, Patrick J. Roache, Christo- pher J. Freitas, Hector Coleman, and Paul E. Raad, “Procedure for estimation and reporting of uncertainty due to discretization in cfd applications,”Journal of Flu- ids Engineering 130, 078001 (2008)

    2008

  47. [47]

    Guzmánn, Numerical Methods for Initial Value Problems in Physics (Springer, 2023)

    Francisco S. Guzmánn, Numerical Methods for Initial Value Problems in Physics (Springer, 2023)

    2023

  48. [48]

    Verifi- cation and validation in computational fluid dynamics,

    William L. Oberkampf and Timothy G. Trucano, “Verifi- cation and validation in computational fluid dynamics,” Progress in Aerospace Sciences38, 209–272 (2002). Appendix A: Numerical validation This appendix presents a convergence test of the evo- lution scheme. For this we use the wave-packet in Eqs. ( 16)–(18), including the same physically motivated cho...

    2002