pith. sign in

arxiv: 2605.29069 · v1 · pith:EEFJRSH2new · submitted 2026-05-27 · 🌌 astro-ph.GA · astro-ph.CO· gr-qc· math-ph· math.MP

Equilibrium Core and Vortex Solutions of Bose Einstein Condensate Dark Matter around a Black Hole

Pith reviewed 2026-06-29 10:28 UTC · model grok-4.3

classification 🌌 astro-ph.GA astro-ph.COgr-qcmath-phmath.MP
keywords Bose-Einstein condensate dark matterblack holevortex solutionsaxisymmetric equilibriastability analysisself-interactionenthalpy functional
0
0 comments X

The pith

BEC dark matter forms stable cores and vortices around a black hole

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

The paper constructs families of stationary axisymmetric solutions for Bose-Einstein condensate dark matter surrounding a central black hole, covering both ground-state cores and first-winding-number vortex lines. An imaginary-time method solves the coupled equations across ranges of self-interaction strength and black hole mass, yielding density profiles that vary systematically with those parameters. Stability along each family is classified by the turning point criterion applied to the enthalpy functional, separating stable from unstable branches. For attractive self-interaction the solutions are bounded by a maximum mass, which restricts the parameter window in which the configurations can form. This extends earlier results on core solutions acting as attractors in collapse to show that vortex solutions can likewise occupy stable branches compatible with formation scenarios.

Core claim

Stationary axisymmetric solutions of BECDM around a point-mass black hole exist for both the ground-state core and the line vortex with nonzero winding number; families of these solutions are obtained for varying self-interaction and black hole mass, with the turning point criterion on the enthalpy functional identifying stable branches, and a maximum mass appearing when the self-interaction is attractive.

What carries the argument

Imaginary-time relaxation of the axisymmetric stationary Gross-Pitaevskii-Poisson system, combined with the turning point criterion on the enthalpy functional to separate stable and unstable branches.

If this is right

  • Core solutions continue to act as attractors in collapse around a black hole.
  • Vortex solutions with nonzero winding number possess stable branches.
  • Attractive self-interaction imposes a maximum mass that bounds the allowable parameter range.
  • Density profiles change continuously with self-interaction strength and black hole mass.

Where Pith is reading between the lines

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

  • If the stable vortex branches survive in full dynamical simulations they could influence angular-momentum transport near galactic-center black holes.
  • The same numerical construction could be repeated for Kerr black holes to test whether frame-dragging alters the location of the stability turning points.
  • The maximum-mass bound for attractive interactions supplies a concrete target for N-body or fluid simulations of BECDM halo assembly around supermassive black holes.

Load-bearing premise

The turning point criterion applied to the enthalpy functional correctly classifies stable and unstable branches for these axisymmetric, self-gravitating BEC configurations.

What would settle it

A long-term numerical time evolution in which a configuration labeled stable by the turning-point criterion remains intact while one labeled unstable collapses or disperses would confirm or refute the stability classification.

Figures

Figures reproduced from arXiv: 2605.29069 by Francisco S. Guzman, Ivan Alvarez-Rios.

Figure 1
Figure 1. Figure 1: FIG. 1. Radial profiles on the equatorial plane ( [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Invariant mass [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Maximum invariant mass ( [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Critical invariant ( [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
read the original abstract

We present the construction of stationary solutions of Bose-Einstein condensate dark matter (BECDM) around a point-like gravitational source representing a black hole. The problem is formulated for general axisymmetric configurations, and we focus on two cases: the ground-state core solution and the first nonzero winding number configuration corresponding to a line vortex solution. The stationary equations are solved using an imaginary-time approach, which enables the construction of families of solutions across a wide range of self-interaction and black hole masses. We analyze the impact of these parameters on the density distribution and on the stability properties of the solutions, assessing stability through the turning point criterion based on the enthalpy functional, which allows us to identify stable and unstable branches along each family of solutions. It has been shown in the past that spherical core solutions act as attractors in the collapse of BECDM around black holes in the non-interacting case ($g=0$), supporting their astrophysical relevance. In the present work, the existence of a maximum mass for configurations with attractive self-interaction ($g<0$) allows us to infer the parameter range in which such solutions may also arise in this regime. Building on this picture, we show that stable vortex solutions of BECDM can also exist in the presence of a black hole, whose stability properties suggest that these configurations may likewise be compatible with physically relevant formation scenarios.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript constructs stationary axisymmetric solutions of the Gross-Pitaevskii-Poisson system describing BECDM around a central point-mass black hole. It obtains families of ground-state core solutions and first-excited vortex solutions (nonzero winding number) via an imaginary-time method across ranges of self-interaction strength g and black-hole mass. Stability is classified by applying the turning-point criterion to the enthalpy functional, identifying stable and unstable branches; for attractive g<0 a maximum mass is reported, and the existence of stable vortices is argued to be compatible with astrophysical formation scenarios.

Significance. If the stability classification holds, the work extends earlier results on spherical g=0 cores to include vortices and nonzero g, supplying concrete parameter ranges (maximum mass for attractive interactions) that could be tested against formation simulations. The systematic construction of solution families parameterized by g and M_BH is a methodological strength.

major comments (1)
  1. [Abstract and stability analysis] Abstract (final paragraph) and stability analysis: the turning-point criterion on the enthalpy functional is used to classify stable versus unstable vortex branches, yet no derivation or cross-validation is supplied showing that the one-parameter turning-point argument remains equivalent to dynamical stability once nonzero phase winding and the central point-mass singularity are present. This directly underpins the central claim that stable vortex solutions exist and are astrophysically relevant.
minor comments (2)
  1. [Abstract] The abstract would benefit from a concise statement of the governing equations and the explicit definition of the enthalpy functional used for the turning-point test.
  2. [Numerical method] Convergence tests, grid resolution, and residual norms for the imaginary-time solver should be reported to allow assessment of the accuracy of the density profiles and masses entering the stability analysis.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments. Below we respond point-by-point to the single major comment, indicating the revision that will be made to address the concern about the turning-point criterion.

read point-by-point responses
  1. Referee: [Abstract and stability analysis] Abstract (final paragraph) and stability analysis: the turning-point criterion on the enthalpy functional is used to classify stable versus unstable vortex branches, yet no derivation or cross-validation is supplied showing that the one-parameter turning-point argument remains equivalent to dynamical stability once nonzero phase winding and the central point-mass singularity are present. This directly underpins the central claim that stable vortex solutions exist and are astrophysically relevant.

    Authors: We acknowledge that the manuscript applies the turning-point criterion to the enthalpy functional without supplying an explicit derivation or numerical cross-validation tailored to the combination of nonzero winding number and central point-mass singularity. The criterion follows from the variational structure of the GPP energy functional (extremization with respect to density at fixed particle number), which has been used previously for spherical cores and for vortex states in the absence of a central singularity. The addition of a fixed central potential and azimuthal phase winding enters the effective potential but preserves the one-parameter family structure used in the turning-point analysis. Nevertheless, we agree that an explicit statement of applicability strengthens the central claim. In the revised manuscript we will add a concise paragraph (or short appendix) in the stability section that (i) recalls the standard derivation of the turning-point method for the GPP system and (ii) notes why the same argument carries over to axisymmetric configurations with winding and an external 1/r potential. This revision will be made without altering the reported solution families or stability classifications. revision: yes

Circularity Check

0 steps flagged

No significant circularity; numerical construction from GPP system

full rationale

The paper numerically constructs stationary axisymmetric solutions (core and vortex) of the Gross-Pitaevskii-Poisson system with central point-mass potential via imaginary-time relaxation, then applies the standard turning-point criterion to the enthalpy functional along families parameterized by g and black-hole mass. No quoted step reduces a claimed prediction or stability classification to a fitted input, self-defined quantity, or load-bearing self-citation chain; the attractor statement for g=0 cores is external prior support. The derivation remains self-contained against the underlying equations.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

Paper rests on standard numerical assumptions for BECDM; no new entities or heavily fitted parameters are introduced in the abstract.

free parameters (2)
  • self-interaction strength g
    Varied parametrically to generate solution families; not fitted to external data.
  • black hole mass
    Varied parametrically to generate solution families; not fitted to external data.
axioms (2)
  • domain assumption Imaginary-time evolution converges to stationary solutions of the axisymmetric Gross-Pitaevskii equation with Newtonian gravity.
    Invoked to construct the families.
  • domain assumption Turning-point criterion on the enthalpy functional identifies dynamical stability for these configurations.
    Used to label stable and unstable branches.

pith-pipeline@v0.9.1-grok · 5788 in / 1307 out tokens · 38691 ms · 2026-06-29T10:28:01.620847+00:00 · methodology

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

65 extracted references · 26 canonical work pages · 4 internal anchors

  1. [1]

    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)

  2. [2]

    Cold and fuzzy dark matter,

    Wayne Hu, Rennan Barkana, and Andrei Gruzinov, “Cold and fuzzy dark matter,” Phys. Rev. Lett.85, 1158– 1161 (2000)

  3. [3]

    Self-gravitating bose-einstein condensates,

    Pierre-Henri Chavanis, “Self-gravitating bose-einstein condensates,” inQuantum Aspects of Black Holes, edited by Xavier Calmet (Springer International Publishing, Cham, 2015) pp. 151–194. 9

  4. [4]

    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)

  5. [5]

    Ultra-Light Dark Matter,

    Elisa G. M. Ferreira, “Ultra-Light Dark Matter,” arXiv e-prints , arXiv:2005.03254 (2020), arXiv:2005.03254 [astro-ph.CO]

  6. [6]

    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 Physics113, 103787 (2020)

  7. [7]

    Wave dark matter,

    Lam Hui, “Wave dark matter,” Annual Review of As- tronomy and Astrophysics59, 247–289 (2021)

  8. [8]

    Systems of self-gravitating particles in general relativity and the concept of an equa- tion of state,

    R. Ruffini and S. Bonazzola, “Systems of self-gravitating particles in general relativity and the concept of an equa- tion of state,” Phys. Rev.187, 1767–1783 (1969)

  9. [9]

    Evolution of the schr¨ odinger-newton system for a self-gravitating scalar field,

    F. S. Guzm´ an and L. Arturo Ure˜ na L´ opez, “Evolution of the schr¨ odinger-newton system for a self-gravitating scalar field,” Phys. Rev. D69, 124033 (2004)

  10. [10]

    Gravita- tional cooling of self-gravitating bose condensates,

    F. S. Guzm´ an and L. Arturo Ure˜ na L´ opez, “Gravita- tional cooling of self-gravitating bose condensates,” The Astrophysical Journal645, 814D819 (2006)

  11. [11]

    Scalar field dark mat- ter: Nonspherical collapse and late-time behavior,

    Argelia Bernal and F. S. Guzm´ an, “Scalar field dark mat- ter: Nonspherical collapse and late-time behavior,” Phys- ical Review D74(2006), 10.1103/physrevd.74.063504

  12. [12]

    Cosmic Structure as the Quantum Interference of a Coherent Dark Wave

    Hsi-Yu Schive, Tzihong Chiueh, and Tom Broadhurst, “Cosmic Structure as the Quantum Interference of a Co- herent Dark Wave,” Nature Phys.10, 496–499 (2014a), arXiv:1406.6586 [astro-ph.GA]

  13. [13]

    Galaxy Formation with BECDM: I. Turbulence and relaxation of idealised haloes

    Philip Mocz, Mark Vogelsberger, Victor H. Robles, Jes´ us Zavala, Michael Boylan-Kolchin, Anastasia Fialkov, and Lars Hernquist, “Galaxy formation with becdm i. turbu- lence and relaxation of idealized haloes,” Mon. Not. Roy. Astron. Soc.471, 4559–4570 (2017), arXiv:1705.05845 [astro-ph.CO]

  14. [15]

    Deep zoom-in simulation of a fuzzy dark matter galactic halo,

    Bodo Schwabe and Jens C. Niemeyer, “Deep zoom-in simulation of a fuzzy dark matter galactic halo,” Phys. Rev. Lett.128, 181301 (2022)

  15. [16]

    Gravita- tional bose-einstein condensation in the kinetic regime,

    D. G. Levkov, A. G. Panin, and I. I. Tkachev, “Gravita- tional bose-einstein condensation in the kinetic regime,” Phys. Rev. Lett.121, 151301 (2018)

  16. [17]

    Formation and mass growth of axion stars in axion miniclusters,

    Benedikt Eggemeier and Jens C. Niemeyer, “Formation and mass growth of axion stars in axion miniclusters,” Phys. Rev. D100, 063528 (2019)

  17. [18]

    New insights into the formation and growth of boson stars in dark matter ha- los,

    Jiajun Chen, Xiaolong Du, Erik W. Lentz, David J. E. Marsh, and Jens C. Niemeyer, “New insights into the formation and growth of boson stars in dark matter ha- los,” Phys. Rev. D104, 083022 (2021)

  18. [19]

    Core-halo mass relation of ultralight ax- ion dark matter from merger history,

    Xiaolong Du, Christoph Behrens, Jens C. Niemeyer, and Bodo Schwabe, “Core-halo mass relation of ultralight ax- ion dark matter from merger history,” Phys. Rev. D95, 043519 (2017)

  19. [20]

    Soliton formation and the core- halo mass relation: An eigenstate perspective,

    J. Luna Zagorac, Emily Kendall, Nikhil Padmanabhan, and Richard Easther, “Soliton formation and the core- halo mass relation: An eigenstate perspective,” Phys. Rev. D107, 083513 (2023)

  20. [21]

    Effect of boundary conditions on structure for- mation in fuzzy dark matter,

    Iv´ an´Alvarez-Rios, Francisco S. Guzm´ an, and Paul R. Shapiro, “Effect of boundary conditions on structure for- mation in fuzzy dark matter,” Phys. Rev. D107, 123524 (2023)

  21. [22]

    Mass-radius relation of self- gravitating bose-einstein condensates with a central black hole,

    Pierre-Henri Chavanis, “Mass-radius relation of self- gravitating bose-einstein condensates with a central black hole,” Eur. Phys. J. Plus134, 352 (2019)

  22. [23]

    Core mass-halo mass relation of bosonic and fermionic dark matter halos harboring a su- permassive black hole,

    Pierre-Henri Chavanis, “Core mass-halo mass relation of bosonic and fermionic dark matter halos harboring a su- permassive black hole,” Physical Review D101(2020), 10.1103/physrevd.101.063532

  23. [24]

    Station- ary solutions of the schr¨ odinger-poisson-euler system and their stability,

    Iv´ an´Alvarez-Rios and Francisco S. Guzm´ an, “Station- ary solutions of the schr¨ odinger-poisson-euler system and their stability,” Physics Letters B843, 137984 (2023)

  24. [25]

    Construction of ground-state solutions of the gross–pitaevskii–poisson system using genetic algo- rithms,

    Carlos Tena-Contreras, Iv´ an Alvarez-R´ ıos, and Fran- cisco S. Guzm´ an, “Construction of ground-state solutions of the gross–pitaevskii–poisson system using genetic algo- rithms,” Universe10(2024), 10.3390/universe10080309

  25. [26]

    Gravita- tional atoms: General framework for the construction of multistate axially symmetric solutions of the schr¨ odinger- poisson system,

    F. S. Guzm´ an and L. Arturo Ure˜ na L´ opez, “Gravita- tional atoms: General framework for the construction of multistate axially symmetric solutions of the schr¨ odinger- poisson system,” Phys. Rev. D101, 081302 (2020)

  26. [27]

    Scalar field dark matter as an alternative explanation for the anisotropic distribution of satellite galaxies,

    Jordi Sol´ ıs-L´ opez, Francisco S. Guzm´ an, Tonatiuh Matos, Victor H. Robles, and L. Arturo Ure˜ na L´ opez, “Scalar field dark matter as an alternative explanation for the anisotropic distribution of satellite galaxies,” Phys. Rev. D103, 083535 (2021)

  27. [28]

    Vortices in a bose-einstein condensate,

    M. R. Matthews, B. P. Anderson, P. C. Haljan, D. S. Hall, C. E. Wieman, and E. A. Cornell, “Vortices in a bose-einstein condensate,” Phys. Rev. Lett.83, 2498– 2501 (1999)

  28. [29]

    Imprint- ing vortices in a bose-einstein condensate using topolog- ical phases,

    A. E. Leanhardt, A. G¨ orlitz, A. P. Chikkatur, D. Kielpin- ski, Y. Shin, D. E. Pritchard, and W. Ketterle, “Imprint- ing vortices in a bose-einstein condensate using topolog- ical phases,” Phys. Rev. Lett.89, 190403 (2002)

  29. [30]

    Rotating trapped bose-einstein condensates,

    Alexander L. Fetter, “Rotating trapped bose-einstein condensates,” Rev. Mod. Phys.81, 647–691 (2009)

  30. [31]

    C. J. Pethick and H. Smith,Bose–Einstein Condensation in Dilute Gases, 2nd ed. (Cambridge University Press, 2008)

  31. [32]

    Angular mo- mentum and vortex formation in bose-einstein-condensed cold dark matter haloes: Angular momentum in bec-cdm haloes,

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

  32. [33]

    Vortices in bose-einstein condensate dark matter,

    Ben Kain and Hong Y. Ling, “Vortices in bose-einstein condensate dark matter,” Physical Review D82(2010), 10.1103/physrevd.82.064042

  33. [34]

    Dynamical galactic effects induced by solitonic vortex structure in bosonic dark matter,

    K. Korshynska, Y. M. Bidasyuk, E. V. Gorbar, Junji Jia, and A. I. Yakimenko, “Dynamical galactic effects induced by solitonic vortex structure in bosonic dark matter,” The European Physical Journal C83(2023), 10.1140/epjc/s10052-023-11548-1

  34. [36]

    Cosmic filament spin from dark matter vortices,

    Stephon Alexander, Christian Capanelli, Elisa G. M. Fer- reira, and Evan McDonough, “Cosmic filament spin from dark matter vortices,” Phys. Lett. B833, 137298 (2022), arXiv:2111.03061 [astro-ph.CO]

  35. [37]

    Kinematic imprints of vortex lines of bec dark matter on baryonic matter,

    Iv´ an´Alvarez Rios, Carlos Tena-Contreras, and Fran- cisco S. Guzm´ an, “Kinematic imprints of vortex lines of bec dark matter on baryonic matter,” Physical Review D111(2025), 10.1103/x9mt-wprk. 10

  36. [38]

    Barra and P

    K. Korshynska, O. O. Prykhodko, E. V. Gorbar, Junji Jia, and A. I. Yakimenko, “Vortex lines in ultralight bosonic dark matter around rotating supermassive black holes,” Physical Review D111(2025), 10.1103/phys- revd.111.023006

  37. [39]

    Scalar dark matter vortex stabilization with black holes,

    Noah Glennon, Anthony E. Mirasola, Nathan Musoke, Mark C. Neyrinck, and Chanda Prescod-Weinstein, “Scalar dark matter vortex stabilization with black holes,” Journal of Cosmology and Astroparticle Physics 2023, 004 (2023)

  38. [40]

    Measuring Distance and Properties of the Milky Way's Central Supermassive Black Hole with Stellar Orbits

    A. M. Ghez, S. Salim, N. N. Weinberg,et al., “Measur- ing distance and properties of the milky way’s central supermassive black hole with stellar orbits,” Astrophys. J.689, 1044–1062 (2008), arXiv:0808.2870 [astro-ph]

  39. [41]

    Monitoring stellar orbits around the Massive Black Hole in the Galactic Center

    S. Gillessen, F. Eisenhauer, S. Trippe,et al., “Monitor- ing stellar orbits around the massive black hole in the galactic center,” Astrophys. J.692, 1075–1109 (2009), arXiv:0810.4674 [astro-ph]

  40. [42]

    A geometric distance measurement to the galactic center black hole with 0.3% uncertainty,

    GRAVITY Collaboration, R. Abuter, A. Amorim,et al., “A geometric distance measurement to the galactic center black hole with 0.3% uncertainty,” Astron. Astrophys. 625, L10 (2019), arXiv:1904.05721 [astro-ph.GA]

  41. [43]

    Axion star nucleation in dark mini- halos around primordial black holes,

    Mark P. Hertzberg, Enrico D. Schiappacasse, and Tsu- tomu T. Yanagida, “Axion star nucleation in dark mini- halos around primordial black holes,” Phys. Rev. D102, 023013 (2020)

  42. [44]

    Parasitic black holes: The swallowing of a fuzzy dark matter soliton,

    Vitor Cardoso, Taishi Ikeda, Rodrigo Vicente, and Miguel Zilh˜ ao, “Parasitic black holes: The swallowing of a fuzzy dark matter soliton,” Physical Review D106 (2022), 10.1103/physrevd.106.l121302

  43. [45]

    Supersonic friction of a black hole traversing a self- interacting scalar dark matter cloud,

    Alexis Boudon, Philippe Brax, and Patrick Valageas, “Supersonic friction of a black hole traversing a self- interacting scalar dark matter cloud,” Phys. Rev. D108, 103517 (2023)

  44. [46]

    Ac- cretion of self-interacting scalar field dark matter onto a reissner-nordstr¨ om black hole,

    Yuri Ravanal, Gabriel G´ omez, and Normal Cruz, “Ac- cretion of self-interacting scalar field dark matter onto a reissner-nordstr¨ om black hole,” Physical Review D108 (2023), 10.1103/physrevd.108.083004

  45. [47]

    Dynamical friction from ultralight dark matter,

    Yourong Wang and Richard Easther, “Dynamical friction from ultralight dark matter,” Phys. Rev. D105, 063523 (2022)

  46. [48]

    Black holes as condensation points of fuzzy dark matter cores,

    Curicaveri Palomares-Ch´ avez, Iv´ an´Alvarez Rios, and Francisco S. Guzm´ an, “Black holes as condensation points of fuzzy dark matter cores,” Physical Review D 112(2025), 10.1103/fwf5-n21g

  47. [49]

    Black hole merger simulations in wave dark matter environments,

    Jamie Bamber, Josu C. Aurrekoetxea, Katy Clough, and Pedro G. Ferreira, “Black hole merger simulations in wave dark matter environments,” Physical Review D107 (2023), 10.1103/physrevd.107.024035

  48. [50]

    Effect of wave dark matter on equal mass black hole mergers,

    Josu C. Aurrekoetxea, Katy Clough, Jamie Bamber, and Pedro G. Ferreira, “Effect of wave dark matter on equal mass black hole mergers,” Physical Review Letters132 (2024), 10.1103/physrevlett.132.211401

  49. [51]

    Self-interacting scalar dark mat- ter around binary black holes,

    Josu C. Aurrekoetxea, James Marsden, Katy Clough, and Pedro G. Ferreira, “Self-interacting scalar dark mat- ter around binary black holes,” Physical Review D110 (2024), 10.1103/physrevd.110.083011

  50. [52]

    Supermassive black hole binaries in ul- tralight dark matter,

    Benjamin C. Bromley, Pearl Sandick, and Barmak Shams Es Haghi, “Supermassive black hole binaries in ul- tralight dark matter,” Phys. Rev. D110, 023517 (2024)

  51. [53]

    Stationary solu- tions of the wave equation in a medium with nonlinearity saturation,

    N. G. Vakhitov and A. A. Kolokolov, “Stationary solu- tions of the wave equation in a medium with nonlinearity saturation,” Radiophysics and Quantum Electronics16, 783–789 (1973)

  52. [54]

    Stable vortex in bose-einstein condensate dark matter,

    Y. O. Nikolaieva, A. O. Olashyn, Y. I. Kuriatnikov, S. I. Vilchynskii, and A. I. Yakimenko, “Stable vortex in bose-einstein condensate dark matter,” Low Tempera- ture Physics47, 684–692 (2021)

  53. [55]

    Linear stability of nonrelativistic self-interacting boson stars,

    Emmanuel Ch´ avez Nambo, Alberto Diez-Tejedor, Ar- mando A. Roque, and Olivier Sarbach, “Linear stability of nonrelativistic self-interacting boson stars,” Physical Review D109(2024), 10.1103/physrevd.109.104011

  54. [56]

    Scalar dark matter vortex stabilization with black holes,

    Noah Glennon, Anthony E. Mirasola, Nathan Musoke, Mark C. Neyrinck, and Chanda Prescod-Weinstein, “Scalar dark matter vortex stabilization with black holes,” JCAP07, 004 (2023), arXiv:2301.13220 [astro- ph.CO]

  55. [57]

    Numerical solution of partial differ- ential equations using the discrete fourier transform,

    Daniela Estefan´ ıa Rodr´ ıguez Lara, Iv´ an´Alvarez, and Francisco Guzm´ an, “Numerical solution of partial differ- ential equations using the discrete fourier transform,” Re- vista Mexicana de F´ ısica E22(2025), 10.31349/revmex- fise.22.020221

  56. [58]

    Thomas,Numerical Partial Differential Equations: Conservation Laws and Elliptic Equations, Texts in Ap- plied Mathematics (Springer New York, 2013)

    J.W. Thomas,Numerical Partial Differential Equations: Conservation Laws and Elliptic Equations, Texts in Ap- plied Mathematics (Springer New York, 2013)

  57. [59]

    Fermion-boson stars as attractors in fuzzy dark matter and ideal gas dynamics,

    Iv´ an Alvarez-Rios, Francisco S. Guzm´ an, and Jens Niemeyer, “Fermion-boson stars as attractors in fuzzy dark matter and ideal gas dynamics,” Physical Review Letters135(2025), 10.1103/4tkh-7hjs

  58. [60]

    Fuzzy dark matter soliton cores around supermassive black holes,

    Elliot Y Davies and Philip Mocz, “Fuzzy dark matter soliton cores around supermassive black holes,” Monthly Notices of the Royal Astronomical Society492, 5721– 5729 (2020)

  59. [61]

    Dynamical friction in a fuzzy dark matter uni- verse,

    Lachlan Lancaster, Cara Giovanetti, Philip Mocz, Yonatan Kahn, Mariangela Lisanti, and David N. Spergel, “Dynamical friction in a fuzzy dark matter uni- verse,” Journal of Cosmology and Astroparticle Physics 2020, 001–001 (2020)

  60. [62]

    Ejection of supermassive black holes and implica- tions for merger rates in fuzzy dark matter haloes,

    Amr A El-Zant, Zacharias Roupas, and Joseph Silk, “Ejection of supermassive black holes and implica- tions for merger rates in fuzzy dark matter haloes,” Monthly Notices of the Royal Astronomical Society499, 2575–2586 (2020)

  61. [63]

    Self-Interacting Dark Matter Solves the Final Parsec Problem of Supermassive Black Hole Mergers,

    Gonzalo Alonso- ´Alvarez, James M. Cline, and Caitlyn Dewar, “Self-Interacting Dark Matter Solves the Final Parsec Problem of Supermassive Black Hole Mergers,” Phys. Rev. Lett.133, 021401 (2024), arXiv:2401.14450 [astro-ph.CO]

  62. [64]

    https://zenodo.org/records/19561144,

    “https://zenodo.org/records/19561144,” . Appendix A: Regularity of solutions

  63. [65]

    The regularity conditions (23)on the axis and equatorial plane are respectively∂ r⊥ ϕ0(0, z) = 0 and∂ zϕ0(r⊥,0) =

    Casem= 0 Consider the axial equation for the core configuration withm= 0 −1 2 ∂2 r⊥ ϕ0 + 1 r⊥ ∂r⊥ ϕ0 +∂ 2 z ϕ0 +V T (r⊥, z)ϕ0 =µϕ 0. The regularity conditions (23)on the axis and equatorial plane are respectively∂ r⊥ ϕ0(0, z) = 0 and∂ zϕ0(r⊥,0) =

  64. [66]

    Using the chain rule,∂ r⊥ ϕ0 =ϕ ′(r)r ⊥/rand∂ zϕ0 = ϕ′(r)z/r, the regularity conditions are automatically sat- isfied

    Assuming spherical symmetry of the total potential, 11 VT (r⊥, z) =V T (r) withr= p r2 ⊥ +z 2, the ground state inherits this symmetry so thatϕ 0(r⊥, z) =ϕ(r). Using the chain rule,∂ r⊥ ϕ0 =ϕ ′(r)r ⊥/rand∂ zϕ0 = ϕ′(r)z/r, the regularity conditions are automatically sat- isfied. The Laplacian becomesϕ ′′(r) + 2 r ϕ′(r), and the equation reduces to −1 2 ϕ′′...

  65. [67]

    Near the symmetry axis, the dominant contribution is the centrifugal term−m 2ϕm/r2 ⊥, so to leading order the equation reduces to∂ 2 r⊥ ϕm + 1 r⊥ ∂r⊥ ϕm − m2 r2 ⊥ ϕm ≃0

    Regularity of the Axialm >0Vortex Solution For a vortex configuration with winding numberm > 0, the axial stationary equation is −1 2 ∂2ϕm ∂r2 ⊥ + 1 r⊥ ∂ϕm ∂r⊥ − m2 r2 ⊥ ϕm +∂ 2 z ϕm + VT (r⊥, z)ϕm =µϕ m. Near the symmetry axis, the dominant contribution is the centrifugal term−m 2ϕm/r2 ⊥, so to leading order the equation reduces to∂ 2 r⊥ ϕm + 1 r⊥ ∂r⊥ ϕm...