The Core-Cusp Problem Revisited: ULDM vs. CDM

TL;DR

This paper investigates whether ultra-light dark matter (ULDM) can resolve the core-cusp problem in galaxy halos by comparing its density profiles to those of cold dark matter (CDM), highlighting the need for more data and simulations.

Contribution

It introduces a parameterization of ULDM halo profiles accounting for environmental variability, analyzing their potential to address the core-cusp discrepancy.

Findings

Feasible ULDM profiles can have lower central densities than NFW at relevant radii.

Both ULDM and NFW profiles can fit observational data reasonably well.

More comprehensive data and baryonic feedback simulations are needed for definitive tests.

Abstract

The core-cusp problem is a widely cited motivation for the exploration of dark matter models beyond standard CDM. One such alternative is ULDM; extremely light scalar particles exhibiting wavelike properties on kiloparsec scales. Astrophysically realistic ULDM halos are expected to consist of inner solitonic cores embedded in NFW-like outer halos. The presence of the solitonic core suggests that ULDM may resolve the core-cusp discrepancy associated with pure NFW halos without recourse to baryonic physics. However, it has been demonstrated that the density of ULDM halos can exceed those of comparable NFW configurations at some radii and halo masses, apparently exacerbating the problem rather than solving it. This situation arises because, although solitonic cores are flat at their centres, they obey an inverse mass-radius scaling relationship. Meanwhile, the mass of the inner soliton increases with the total halo mass, and therefore the inner core becomes more peaked at large halo masses. We describe a parameterisation of the radial density profiles of ULDM halos that allows for environmental variability of the core-halo mass relation in order to investigate this issue in more detail. For halos up to 10^12 solar masses we find feasible ULDM profiles for which the central density is lower than their NFW counterparts at astrophysically accessible radii. However, comparisons to observed profiles do not strongly favour either option; both give reasonable fits to subsets of the data for some parameter choices. Consequently, we find that robust tests of the core-cusp problem in ULDM will require more comprehensive observational data and simulations that include baryonic feedback.

Equations

$${\rho }_{\mathrm{NFW}}(r)=\frac{{\rho }_0}{\frac{r}{R_s}{\left(1+\frac{r}{R_s}\right)}^2}.$$

$$i\dot{\psi }=-\frac{1}{2}{\nabla }^2\psi +\Phi \psi ,$$

$${\nabla }^2\Phi =4\pi |\psi {|}^2,$$

$${\psi }^{\prime }(x)=\gamma \psi (\sqrt{\gamma }x),$$

$$\rho$$

$$r$$

$$\mathcal{L}={\left(\frac{8\pi {\hslash }^2}{3m^2{H}_0^2{\Omega }_{m_0}}\right)}^{\frac{1}{4}}\approx 121{\left(\frac{10^{-23}\mathrm{eV}}{m}\right)}^{\frac{1}{2}}\mathrm{kpc},$$

$$\mathcal{M}$$

$$\approx 7\times 10^7{\left(\frac{10^{-23}\mathrm{eV}}{m}\right)}^{\frac{3}{2}}{M}_{\odot }.$$

$$\rho (r)=\{\begin{array}{ll}{\rho }_{\mathrm{sol}}(r),& 0\le r\le r_{\alpha }\\ {\rho }_{\mathrm{NFW}}(r),& r_{\alpha }\le r\le r_{\mathrm{vir}},\end{array}$$

$$M_c=\alpha {\left(|E|/M\right)}^{1/2},$$

$${\rho }_c=2.94\times 10^6{M}_{\odot }{\mathrm{kpc}}^{-3}{\left(\frac{{\mathrm{M}}_{\mathrm{vir}}}{10^9{\mathrm{M}}_{\odot }}\right)}^{4/3}{{\mathrm{m}}_{22}}^2,$$

$$r_c=1.6\mathrm{kpc}{\left(\frac{{\mathrm{M}}_{\mathrm{vir}}}{10^9{\mathrm{M}}_{\odot }}\right)}^{-1/3}\frac{1}{{\mathrm{m}}_{22}},$$

$$V(r{)}^2=\frac{4\pi G}{r}{\int }_0^r\rho (r^{\prime })r^{\prime 2}d{r}^{\prime },$$

$$V^2=V_{\mathrm{disk}}^2+V_{\mathrm{bulge}}^2+V_{\mathrm{gas}}^2+V_{\mathrm{halo}}^2.$$

$$\mathrm{BIC}=k\ln \mathrm{N}-2\ln \mathcal{L},$$

$$\frac{1}{2}{M}_{tot}{v}^2=\frac{1}{4}\frac{G{M}_{tot}^2}{R_{tot}},$$

$$v^2$$

$$\to R_{core}=\frac{M_{core}{R}_{\mathrm{vir}}}{{\mathrm{M}}_{\mathrm{vir}}}.$$

$$R_{core}^2\propto \frac{R_{\mathrm{vir}}}{{\mathrm{M}}_{\mathrm{vir}}}$$

$$\to R_{core}^2\propto \frac{{{\mathrm{M}}_{\mathrm{vir}}}^{1/3}}{{\mathrm{M}}_{\mathrm{vir}}}$$

$$\to R_{core}\propto {\left({{\mathrm{M}}_{\mathrm{vir}}}^{-2/3}\right)}^{1/2}$$

$$\to R_{core}\propto {{\mathrm{M}}_{\mathrm{vir}}}^{-1/3}.$$

Full text

1 INTRODUCTION

It is widely agreed that non-baryonic dark matter constitutes the majority of the mass of the observable universe, but its precise nature remains an open question. Many dark matter models have been proposed, with particle CDM [Cold Dark Matter] being the most widely studied. This scenario successfully accounts for the large scale structure of the universe [43] and the spectrum of anisotropies in the microwave background [12, 16, 15, 32, 20, 1, 18]. Nevertheless, the so-called “small-scale crisis” remains a challenge [47]. A key issue is the apparent tension between the central density profiles of dark matter halos in simulations containing only gravitationally interacting CDM, and those inferred from observational data. Simulations tend to produce ‘cuspy’ central density profiles [30], which grow as $1/r$ at small radii, but observational data appears to favour flattened central cores [29]. The tension between the two is widely known as the “core-cusp problem” [13, 35, 14].

The seriousness of the core-cusp problem is the subject of ongoing debate, as it has been shown that it may be ameliorated in some cases by adding baryonic matter to CDM simulations [5]. Nevertheless, the wider category of “small-scale” problems in standard CDM, along with tighter constraints from direct-detection experiments [39], motivates the study of alternative dark matter models. One scenario which has gained substantial traction is ultra-light dark matter [ULDM], also known as scalar-field dark matter, $\Psi$ dark matter, BEC dark matter and fuzzy dark matter. Often these models are referred to as ‘axion-like’, given that they describe an extremely light scalar field.

As reviewed by Hui et al. [19], current constraints prefer a ULDM particle mass of $\mathcal{O}(\sim 10^{-22}\mathrm{eV})$, corresponding to a kiloparsec-scale de Broglie wavelength. ULDM thus exhibits novel wave-like behaviour on astrophysically interesting scales, as well as supporting soliton-like gravitationally confined Bose-Einstein condensates. ULDM simulations suggest that realistic astrophysical halos have an inner core consisting of a kiloparsec scale condensate, surrounded by a virialised outer halo [40, 45]. The outer region resembles a standard CDM halo which is well approximated by the NFW profile characteristic of collisionless CDM, and most commonly associated with WIMP dark matter [Weakly Interacting Dark Matter] [30]. Because the (solitonic) profiles of the inner condensates are flat, it has been suggested that ULDM can resolve the core-cusp problem without the need to invoke baryonic astrophysics. However, solitonic density profiles obey an inverse mass-radius scaling law, so the density of the ULDM halo might exceed that of an analogous NFW halo over a finite range of small radii in larger galaxies. In particular, Ref [37] concludes that DM-only NFW profiles may outperform ULDM profiles for galaxies with halo masses $\mathrm{M_{h}}\gtrsim 10^{11}\mathrm{M}_{\odot}$.111 Ref [37] refers to these galaxies as “large dwarfs”, though we note that the upper limit on this category is around $10^{12}\mathrm{M_{\odot}}$, approaching estimates of the Milky Way mass [46].

To further explore the possible worsening of the core-cusp problem in ULDM, we examine the effect of scatter in the core-halo mass scaling relation. Starting from the semi-analytic density profile of Ref. [37], we look at the scatter in the parameters implied by Ref. [38]. We show that the resulting statistical variability may ease concerns that the core-cusp problem is exacerbated for ULDM relative to CDM for “large dwarf” galaxies.

Our analysis also highlights a number of caveats that apply to all such discussions. First, the incoherent outer regions of ULDM halos are subject to strong fluctuations, both temporally and spatially. These are not captured by semi-analytic halo density profiles and we argue that these fluctuations may accentuate the intrinsic scatter in halo parameters. Moreover, baryonic feedback is known to be significant for dwarf galaxies [10, 5] and neither the NFW or ULDM profiles incorporate this effect. Consequently, we caution against attempting to discriminate between ULDM and CDM models based on DM-only simplified theoretical profiles.

Observationally, we find that neither semi-analytic ULDM halos (for ULDM particle mass $0.8-2.5\times 10^{-22}\operatorname{eV}$) nor NFW halos provide a particularly convincing fit to rotation curves of large dwarf galaxies in the SPARC database [21]. Moreover, many rotation curves are extracted from a few data points with significant uncertainties and which only span a small range of radial distances, further complicating attempts to draw robust conclusions. These issues are exacerbated by the relatively large number of free parameters in the theoretical models. For instance, a ULDM mass parameter of $10^{-23}\operatorname{eV}$ seemingly ameliorates the core-cusp problem in galaxies exhibiting a steep decrease in rotation velocity at small radii, but such a small mass is in tension with other constraints. Consequently, the primary conclusion to be drawn from this type of parameter-fitting exercise seems to be that analyses of the core-cusp problem (and potentially other “small scale” anomalies) based on simplified semi-analytic DM-only models cannot meaningfully test these scenarios, especially when observational data is limited and detailed numerical simulations with baryonic feedback have not been performed.

The structure of the paper is as follows. In Section 2, we review the construction of semi-analytic density profiles for both the ULDM and CDM models and briefly discuss aspects of realistic ULDM halos which are not captured by the semi-analytic model. In Section 3 we compare the semi-analytic density profiles for ULDM and CDM halos in the dwarf galaxy mass range $10^{11}-10^{12}\operatorname{M}_{\odot}$, taking into account statistical variation in both the NFW concentration parameter and the ULDM core-halo mass relation. We then compare the radial velocity profiles inferred from these density profiles with astrophysical data from the SPARC database [21]. We conclude in Section 4.

2 Semi-analytic Halos

We begin by looking at the semi-analytic parametrisations of ULDM and CDM halos. The well known NFW profile of CDM [30, 26] is given by

[TABLE]

At small radii the profile is proportional to $1/r$, while at large radii it goes as $1/r^{3}$. The parameters $\rho_{0}$ and $R_{s}$ vary from halo to halo; $\rho_{0}$ can be interpreted as a characteristic density, while $R_{s}$ is the scale radius and determines the distance from the centre at which the transition between the ‘small $r$’ and ‘large $r$’ limits occurs.

The NFW halo is assumed to be radially symmetric, and requires truncation at a finite radius in order to prevent the integrated mass diverging as $r\rightarrow\infty$. The truncation is typically set by the virial radius, which is itself determined approximately via the spherical top-hat collapse model describing the evolution of a uniform spherical overdensity in a smooth expanding background [48, 44, 17]. Gravitational collapse of the overdensity halts when virial equilibrium is reached. In this scenario the corresponding virial radius is the radius at which the mean internal density is $\Delta_{c}\rho_{\mathrm{crit}}(t)$. Here $\rho_{\mathrm{crit}}(t)$ is the critical density of the universe at time $t$. The factor $\Delta_{c}$ is of order $10^{2}$ and while different conventions exist, we make the common choice $\Delta_{c}=200$ [36] in what follows.

Once the virial radius is specified as the outer limit of the halo, Equation 1 completely determines the density profile for given $\rho_{0}$ and $R_{s}$. For any given virial mass, there is a range of corresponding NFW density profiles, with the distributions of $\rho_{0}$ and $R_{s}$ emerging from the mass-concentration-redshift relation seen in N-body simulations and observations [25, 34].

Whereas CDM halos can be described by NFW distributions, a different approach must be taken in the case of ULDM. ULDM dynamics is governed by the Schrödinger-Poisson system of coupled differential equations. In a static background, they take the dimensionless form

[TABLE]

where $\psi$ is the ULDM wavefunction, $\Phi$ is the Newtonian potential, and the density $\rho\propto|\psi|^{2}$. The solitonic ground state profile cannot be written down analytically, but given a numerically computed spherically symmetric profile $\psi$ for $\psi(0)=1$, the full family of solutions is then given by

[TABLE]

where $\gamma$ is a scaling parameter and the dimensionless mass of the soliton is proportional to $\sqrt{\gamma}$, while the dimensionless radius is proportional to $1/\sqrt{\gamma}$. The dimensionless density $|\psi|^{2}$ and dimensionless radius $x$ can be transformed into dimensionful quantities by

[TABLE]

where

[TABLE]

and

[TABLE]

Ref. [37] gives a piecewise parameterization of the generic ULDM profile

[TABLE]

where $\rho_{\mathrm{sol}}(r)$ is the appropriately scaled density profile of the ground state soliton solution. The contribution from the solitonic core and the overall virial mass is predicted to obey a scaling relationship [38, 8] which sets the central density, $\rho_{c}$, of a ULDM halo with virial mass, $\mathrm{M_{vir}}$. This yields an expression relating the core size to the velocity dispersion, and finally to the halo virial mass.222The authors of [38] suggest the following general expression:

$M_{c}=\alpha\left(|E|/M\right)^{1/2},$

(10)

where the core mass $M_{c}$ is determined by the total energy, $E$, and the total mass of the halo, $M$ where $\alpha$ is a constant of order unity. They then explain that the right hand side of the equation represents the halo velocity dispersion, while the left hand side represents the inverse core size due to soliton scaling laws. By invoking the virial condition of the spherical collapse model, the authors then construct the redshift dependent relationship between the solitonic core mass and the halo virial mass for a ULDM halo.

This core-halo mass relation can also be understood simply by matching the virial velocities of the core and the wider halo (see Appendix A for details). At $z=0$ the relationship is found to be [38]

[TABLE]

and

[TABLE]

where $r_{c}$ is the radius at which the density is half of the central value, and $\mathrm{m_{22}}$ is given by $\mathrm{m_{22}}\equiv m/10^{-22}\operatorname{eV}$ where $m$ is the ULDM particle mass.

While the piecewise semi-analytic ULDM profile is a useful tool, one should be mindful of its limitations. For example, while a number of studies have attempted to establish ‘universal’ properties of ULDM halos, many of these analyses generated ULDM halos through the mergers of smaller compact objects [40, 28]. This method of halo assembly is not representative of a realistic structure formation process, however it has the advantage of avoiding the computational difficulty of undertaking large-scale ULDM cosmological simulations. For this reason there is currently limited information from which to draw robust conclusions about the properties of astrophysical ULDM halos. In particular, universal application of the core-halo mass relation cannot be fully justified until more work is done to understand the characteristic timescales associated with the formation of quantum pressure-supported cores in scenarios including condensation from a fluctuating background, gravitational collapse in an expanding background, and mergers of objects with and without stable central cores. Moreover, it is difficult to accurately predict the effect that baryonic feedback will have on the formation of solitonic cores in halos of different masses, which could be significant at small radii in the present context.

Halo substructure is likewise missing from the semi-analytic model presented above. In simulations of soliton mergers the resulting halos have turbulent outer regions, with fluctuations on scales comparable to the core size, as illustrated in Figure 1. In addition to the fluctuations inherent in a large ULDM halo, smaller halos are likely to orbit or interact with larger halos. This substructure is not captured by the semi-analytic model described above, and predictions for tracer velocity profiles may thus not match those of realistic astrophysical objects. Furthermore, temporal fluctuations in the core density are also missing from the semi-analytic model. Realistic halo cores are not exact soliton solutions of the Schrödinger-Poisson equation, they interact non-trivially with the fluctuating NFW-like outer halo, and their central densities can vary with time by as much as a factor of two [45].

Taken together, these limitations suggest that the core-halo mass relation of the semi-analytic model should not be interpreted as an inviolable rule, but as a statement about the averaged characteristics of a statistical distribution. To estimate the variance corresponding to this distribution, we can consider a range of possible central densities for a given virial mass (somewhat analogous to the scatter in NFW concentration parameters [26]). The results of Ref. [38] indicate that a scatter in the core mass $\mathrm{M_{c}}$ of up to $\pm 50\%$ is possible for a given virial mass. Unfortunately, the small sample size and limited halo mass range ($\mathrm{M_{vir}}\approx 10^{8}-10^{11}\operatorname{M}_{\odot}$) found in [38] preclude a detailed analysis of the statistical properties of realistic astrophysical halos, but future simulations (especially those including baryonic feedback) should lead to improved predictions for this distribution.

To partially account for statistical variance in halo properties, one can allow for variation in the radius at which the solitonic profile of the ULDM halo transitions into an NFW profile. This is acknowledged in Ref [37] and is captured by the parameter $\alpha$: the transition radius, $r_{\alpha}$, is given by $r_{\alpha}=\alpha r_{c}$, with $3\leq\alpha\leq 4$. For a given theoretical halo, an adjustment to the transition radius should be accompanied by changes in the parameters of the outer NFW halo, so as to maintain the core-halo mass ratio.

Thus, by taking the central soliton density and transition radius as variable parameters, we can create a range of plausible ULDM halo profiles for a given halo by using the virial mass to predict $\rho_{c}$, and assuming variation of $\pm 50\%$ around this central value. Given specific values for the central density and transition parameter $\alpha$, the solitonic piece of the ULDM profile is then completely specified, and its mass can be calculated. The remainder of the virial mass must be accounted for by the NFW tail of the profile. By matching the densities of the NFW tail to the inner soliton at the transition radius, the values of $R_{s}$ and $\rho_{0}$ for the NFW tail are obtained.

3 ULDM and CDM halos and astrophysical data

We now compare the radial profiles of ULDM halos to NFW halos using the semi-analytic profiles described above, focusing on masses in the range $10^{11}$ and $10^{12}\operatorname{M}_{\odot}$, which may show an apparent worsening of the core-cusp problem [37]. Figure 2 compares such halos; the shaded blue region represents the ULDM halos for which the core-halo mass relation has a scatter $\operatorname{M}_{c}=\operatorname{M}_{\mathrm{cp}}\pm 50\%$ range, where $\operatorname{M}_{\mathrm{cp}}$ is the theoretical prediction for the core mass. Note that because higher central densities correspond to narrower soliton profiles, the shaded region possesses ‘crossover points’ near the transition from the solitonic to NFW profile, appearing somewhat skewed from the median line. Were we to vary more parameters in the model (such as transition radius and axion mass), we would see a broadening of the shaded region, such that the median line would be fully encompassed by the shaded region. Because we are here focusing primarily on the core mass (and therefore central density), we illustrate only the changes in density profile attributable to this, hence the restricted range of profiles shown as the shaded blue region.

The Schrödinger-Poisson soliton scaling relations show that the $\operatorname{M}_{c}=\operatorname{M}_{\mathrm{cp}}\pm 50\%$ mass range corresponds to a range of $\gamma_{p}/4\leq\gamma\leq 9\gamma_{p}/4$, where $\gamma_{p}$ is the theoretical prediction of the square root of the dimensionless central density. Consequently, there is a large variation in the central density and thus widely varying predictions for the ULDM profiles. We fix $\alpha=3.5$ (in the middle of the predicted range) which does not affect the central density as the core lies well within the solitonic region. Changing the value of $\alpha$ will, however, affect the predicted velocity profiles for each halo. We do not attempt to fit this parameter to data in this Section; the previously discussed limitations of the semi-analytic models employed here suggest that this would be unlikely to be a meaningful exercise. The blue ULDM profiles are compared to the red shaded regions of Figure 2, showing the $2\sigma$ variation about the theoretical prediction for the concentration parameter of the NFW halo with the same virial mass [26].

Following Ref [37], we plot to a minimum radius of $r/r_{\mathrm{vir}}=10^{-4}$ and for the same choices of $\mathrm{m_{22}}$. For any $\operatorname{M}_{\mathrm{vir}}$, the NFW halo density will inevitably exceed that of the ULDM halo at very small radii, though the threshold for this transition may be arbitrarily small, and not observationally relevant. However, we note that the apparent worsening of the core-cusp discrepancy does depend on the choice of inner radial cutoff.

From Figure 2 we see that for halo masses of $10^{11}\operatorname{M}_{\odot}$ there is a wide range of $\mathrm{M_{c}}$ for which the ULDM profile is ‘less cuspy’ than its NFW counterpart. For a halo mass of $10^{12}\operatorname{M}_{\odot}$ and a ULDM particle mass $\mathrm{m_{22}}=0.8$ the range of plausible ULDM profiles likewise includes those which are ‘less cuspy’ than the corresponding NFW profile. At higher particle mass ($\mathrm{m_{22}}=2.5$) for $10^{12}\operatorname{M}_{\odot}$ halos, the NFW profiles tend to be less peaked than corresponding ULDM profiles at radial distances in the range $10^{-4}\leq r/r_{\mathrm{vir}}\leq 1$.

To assess the suitability of these semi-analytic profiles, we compare to observations drawn from the SPARC database. Because observations yield the (line of sight) velocity distributions of tracer stars as a function of galactocentric radius rather than the halo density itself, we must first transform our theoretical density profiles into velocity profiles. In so doing, we acknowledge that the effects of non-circular motion and kinematic irregularities constitute a non-trivial source of random error in observed velocities, which should be kept in mind especially when working with limited data sets.

We convert density profiles to velocity distributions [42] via

[TABLE]

where

[TABLE]

The SPARC database contains photometric data for 175 galaxies and rotation curves from $\mathrm{H}_{\mathrm{I}}$/$\mathrm{H}_{\alpha}$ studies. The disk and bulge velocities in the SPARC database are given for $\Upsilon=1\operatorname{M}_{\odot}/\operatorname{L}_{\odot}$ at $3.6\operatorname{\mu m}$. However, the greatest source of uncertainty in mass modelling is the assumed stellar mass-to-light ratio, $\Upsilon_{\star}$ [21]. As in [37], we assume a constant value of $\Upsilon_{\star}=0.2\operatorname{M}_{\odot}/\operatorname{L}_{\odot}$ at $3.6\operatorname{\mu m}$, likewise noting that this constitutes a non-trivial source of uncertainty. Moreover, there is significant uncertainty in the SPARC data itself. Error bars are omitted in the following graphs for ease of viewing, however, they are discussed in AppendixC.

The characteristics of the velocity profiles in the SPARC database vary widely from galaxy to galaxy; however we qualitatively identify two subsets of galaxies; those with maximum tracer velocities $75\leq v<125\operatorname{kms}^{-1}$, and those for which $125\leq v<175\operatorname{kms}^{-1}$. The former group tends to exhibit a strong steepening in the radial velocity profile toward the inner halo, while the profiles for the latter group are comparatively flat333We exclude data for which the velocities calculated according to Equation 14 are inconsistent - this can occur due to the uncertainty in the assumption for $\Upsilon_{\star}$.. We assume that higher asymptotic velocities correspond to a larger halo mass, and consider halo masses in the range $10^{11}-10^{12}\operatorname{M}_{\odot}$, expecting that masses at the top end of the range will give a better match to galaxies with higher asymptotic velocities.

In Figure 3, we see that galaxies with asymptotic velocities at the higher end of the range do not always exhibit a pronounced steepening of the velocity profile at small radii. Indeed in some cases there is simply no data at small radii. From this figure we see that while a halo mass of $10^{12}\mathrm{M}_{\odot}$ with $\mathrm{m_{22}}=2.5$ provides a reasonable fit to the data at radii $>10\mathrm{kpc}$, it is difficult to judge the fit at small radii, where the ULDM and NFW profiles differ most strongly, due to a lack of data. Furthermore, while the data at higher radii seems to be relatively clustered, there are significant deviations within the limited data that exists at small radii. For example, the curves for NGC1090 and NGC6946 are widely disparate at small radii, but seem to converge at larger radii. Attempting to fit such data to a single set of model parameters would be of limited utility without a much more comprehensive data set from which statistical outliers could be properly identified. Furthermore, we note that there are substantial changes in theoretical ULDM velocity profiles under variation in the ULDM particle mass. The scale of these changes is exhibited in B. As such, we remark that analyses of the sort presented here would benefit greatly from tighter constraints on the ULDM particle mass.

By contrast, for galaxies with smaller maximum velocities ($75\leq v<125\operatorname{kms}^{-1}$), there is more data at smaller radii. For such galaxies we see the steepening rotation curves characteristic of cored density profiles, as shown in Figure 4. In this case, choosing parameters such that the theoretical profiles overlap with the data at small radii is easy (in this case $\mathrm{m_{22}}=0.1$, $\mathrm{M}_{\mathrm{vir}}5\times 10^{11}\mathrm{M}_{\odot}$), however it is not clear whether the behaviour of this profile would fit data at larger radii were it available. Furthermore, while the choice $\mathrm{m_{22}}=0.1$ provides a reasonable fit to the data in this case, a ULDM particle mass $\mathrm{m_{22}}=0.1$ is in tension with constraints from the Lyman-$\alpha$ forest, as well as high-redshift UV luminosity function comparisons444$0.8<\mathrm{m_{22}}<2.5$ is preferred by current constraints, as mentioned in [37] [2, 6, 3, 33, 31]. It must be acknowledged, however, that baryonic feedback is expected to have the greatest impact in the innermost regions of realistic halos. As such, agreement between our semi-analytic DM-only model and observational data at small radii should be interpreted cautiously, especially since this is also the region where assumptions regarding the stellar mass to light ratio have the greatest significance.

4 Conclusions

The ULDM model has gained attention in part because it may offer a solution to the CDM core-cusp problem. However, in some cases ULDM profiles can actually have higher densities than their NFW counterparts at observationally relevant radii in the interior of halos with mass $\mathrm{M}\gtrsim 10^{12}\mathrm{M}_{\odot}$, where the central density is determined by the theoretical core-halo mass relation. However, apparent statistical spread in the ULDM core-halo mass relation [38] leads to a sizeable range of plausible central densities for a halo of any given mass. Furthermore, analyses of oscillations of the cores of ULDM halos on timescales much smaller than the relaxation time have demonstrated significant fluctuations in central density [45]. This suggests that theoretical core-halo mass relations should not be interpreted too literally for any individual ULDM halo, and this should be taken into account when performing model-selection analyses. The limited available simulation data means that the exact features of the distribution of halo properties in ULDM are poorly characterised. Nevertheless, it remains apparent that core masses at the lower end of the plausible range could mitigate the apparent worsening of the core-cusp discrepancy for ULDM halos.

When the spread in the theoretical core-halo relation is accounted for, comparisons of theoretical ULDM and NFW profiles to the SPARC database yield inconclusive results as far as interior regions of the halos are concerned. Parameters can be easily chosen to provide a superficial fit to given subsets of data. However, the available data often do not span a large enough range of radial values to assess the relative merits of the UDLM and NFW profiles over the whole profile. In particular, rotation curves in the SPARC database which exhibit a strong steepening at small radii often lack data at large radii, whereas rotation curves which include large radii information often lack measurements at small radii. From our restricted analysis, however, it appears that neither the theoretical ULDM nor CDM model can reliably reproduce the data across a broad range. Wide deviations at small galactocentric radii cannot be accounted for by random measurement errors, and suggest that sophisticated modelling of baryonic physics will be necessary before any conclusions can be reasonably drawn. The lack of baryonic physics in both the semi-analytic CDM and ULDM models is significant, as are other limitations due to poorly characterised statistics and simplistic assumptions about halo modelling.

In principle, one could perform a BIC analysis to determine which of the ULDM or NFW models is more strongly favoured by the data [22]. The model with the lowest BIC, defined as

[TABLE]

is preferred. Here $\mathcal{L}$, $\mathrm{N}$, and $k$ are the maximised likelihood function, the sample size, and the number of model parameters, respectively. The utility of BIC or other model selection tools, however, is hampered by the lack of comprehensive data, the high number of free parameters (the stellar mass to light ratio in the SPARC data, assumed virial mass of the galactic halos, ULDM particle mass, the NFW concentration parameter, the UDLM soliton to NFW transition radius and variation in the ULDM core-halo mass relation) and the omission of baryonic feedback in ULDM simulations. Indeed, BIC analyses are known to be compromised when the sample size of the data is not sufficiently large in comparison to the number of free parameters in the model. Because of this limitation, and large and unquantified systematic biases in both the observational data and theoretical predictions, such analyses are premature at this point. Previous studies such as Ref. [4] of the core-halo mass relation and the fitting of semi-analytical profiles to galaxy data also emphasise the necessarily preliminary and tentative nature of all analyses of ULDM-derived rotation curves.

To summarise, the parameter space describing “typical” ULDM halos is larger than simple semi-analytical models suggest. It is necessary to constrain this parameter space in order to make robust model selection possible. Tightening the constraints on the plausible ULDM particle mass [7, 23, 11] and obtaining additional spectroscopic data with improved uncertainties covering a greater halo mass range and radius would be of tremendous benefit in this regard. Such improved data can be expected from future surveys [41]. In addition, better ULDM cosmological structure formation simulations are needed to improve the understanding of ULDM halo evolution [24, 9, 27] and these should also incorporate baryonic feedback. Thus, the key conclusion to be drawn from this work is that more information from simulations and astrophysical observations is needed, as is more sophisticated incorporation of baryonic effects within semi-analytic models of both ULDM and CDM, before the relative successes of each model can be fairly compared.

Acknowledgements.

We acknowledge invaluable discussions with Jens Niemeyer, Shaun Hotchkiss, and Mateja Gosenca in completing this work. We also acknowledge support from the Marsden Fund of the Royal Society of New Zealand. This research was supported by use of the Nectar Research Cloud, a collaborative Australian research platform supported by the National Collaborative Research Infrastructure Strategy (NCRIS).

Appendix A Core-halo mass relation

The core-halo mass relation can be simply interpreted as the statement that the average internal velocity of a tracer mass in the core must be equal to the virial velocity of a tracer mass in the wider halo. If this were not the case, and instead the average velocity were higher within the core, these higher velocity particles would move outward, resulting in dynamical mass redistribution within the halo. During this process, the halo would not be in equilibrium and would thus not be virialised.

From the virial theorem we have that $E_{K}=-1/2\ E_{P}$, where $E_{K}$ and $E_{P}$ represent kinetic and potential energies, respectively. Alternatively we can write:

[TABLE]

where $G$ is the gravitational constant, $M_{tot}$ and $R_{tot}$ are the total mass and radius, and $v^{2}$ is the mean of the squares of individual tracer velocities. Demanding that $v^{2}$ is the same for the core as for the total virialised halo allows us to then write

[TABLE]

We know from the soliton scaling properties that $R_{core}\propto M_{core}^{-1}$, and since $\mathrm{M_{vir}}=4/3\ \pi R_{\mathrm{vir}}^{3}\bar{\rho}$, we also have $R_{\mathrm{vir}}\propto\mathrm{M_{vir}}^{1/3}$. Hence, Equation A becomes

[TABLE]

With this scaling relation in mind, the constant of proportionality may be determined through analysis of simulated halos.

Appendix B - Impact of ULDM particle mass on halo velocity profiles

Figure 5 demonstrates the scale of the changes to the velocity profiles of theoretical ULDM halos under changes in the ULDM particle mass. In order to perform a meaningful parameter fitting exercise, observational data would be required to span the wide range of radii illustrated here. In this way, the regions of the rotation curves most sensitive to the assumption for the ULDM particle mass could be tested simultaneously. Presently available data, when binned according to e.g. maximum velocity, is likely to yield disparate preferences for the ULDM particle mass, as illustrated in Figures 3 and 4. Further work to constrain the plausible range of the particle mass will make comparisons of the ULDM and CDM models with astrophysical data more effective.

Appendix C Errors in SPARC data

As discussed in Section 3, the uncertainties associated with the SPARC rotation curves for the galaxies studied here make it difficult to draw robust conclusions as to the suitability of one or the other dark matter model. Indeed, sources of error quoted in the SPARC database relate not only to the random error in measured line-of-site velocities, but also to errors on the galaxy distance measurement, inclination, and total luminosity. Furthermore, inaccuracy in the assumed stellar mass-to-light ratio may lead to skewed velocity decompositions, a systematic effect that could exceed the stochastic measurement errors.

Figures 6 and 7 show the error bars associated with the low asymptotic velocity ($75\leq v<125\operatorname{kms}^{-1}$) and high asymptotic velocity ($125\leq v<175\operatorname{kms}^{-1}$) measurements, respectively. In Figure 6, we see large error bars at small radii. It is precisely this regime in which accurate velocity profiles are needed to assess the significance of the core-cusp discrepancy - a key differentiating factor between ULDM and CDM models as illustrated in Figure 4. Hence, more comprehensive, accurate data in this regime would be of tremendous benefit. Furthermore, in Figure 4 we also observe that at higher radii, the ULDM model exhibits a characteristic dip in the radial velocity profile. The observational data does not extend far enough into the high radius regime to reveal whether such features are exist in astrophysical structures. Indeed, the data in this case tapers off at galactocentric radii exceeding around 10kpc. From the limited number of galaxies for which data approaches this regime (in particular UGC03580, UGC01230, NGC3769, NGC1003), there does not seem to be a tendency toward a dip. The error bars are relatively constrained in this high radius region, as shown in Figure 6, so the absence of a dip arguably weakens support for the ULDM model. More data at higher radii is required to make a strong determination on this point.

Meanwhile, in Figure 7, the spread of data at high radius is within the scale of the error bars, but at smaller radii the data is not encompassed by random measurement error alone. This may suggest that grouping galaxies by asymptotic velocities alone is insufficient as a method of characterisation. This spread of data may suggest that grouping galaxies based on asymptotic velocities alone is an insufficient method of characterisation. However, there is very limited data in this sample at small radii, so both larger data sets and comprehensive modelling of baryonic effects in high density inner regions are required to resolve this issue.