SPARC Sample Research Articles

On the time dependency of [formula omitted

In this paper, we test one of the prediction of the Scale Invariant Vacuum (SIV) theory on MOND. According to that theory, MOND’s acceleration a0 is not a universal quantity but depends from several parameters, and from time. To check this prediction, we compare the dependency of a0 from redshift, using the values of a0 obtained by Marra et al. (2020) carrying a Bayesian inference for 153 galaxies of the SPARC sample. Since this sample does not contain galaxies at large redshift, we estimated a0 from the Nestor Shacha et al. (2023) data. The SPARC sample, for small values of the redshift, in the redshift range 0.00032–0.032 gives a correlation with z, and the Nestor Shacha et al. (2023) data, in the higher redshift range 0.5–2.5 gives an anti-correlation with the redshift z. Both samples show a dependency of a0 from z. The different behavior at low and high redshift can be related to a change of the a0(z) with redshift, or to the lower precision with which the high redshift value of a0 are known.

CASCO: Cosmological and AStrophysical parameters from Cosmological simulations and Observations – I. Constraining physical processes in local star-forming galaxies

ABSTRACT We compare the structural properties and dark matter content of star-forming galaxies taken from the Cosmology and Astrophysics with MachinE Learning Simulations (camels) cosmological simulations to the observed trends derived from the SPARC sample in the stellar mass range $[10^{9}, 10^{11}]\, \textrm {M}_{\odot }$, to provide constraints on the value of cosmological and astrophysical (supernova- and active galactic nucleus-related) parameters. We consider the size–, internal DM fraction–, internal DM mass–, and total–stellar mass relations for all the 1065 simulations, all having different cosmological and astrophysical parameters, from the IllustrisTNG, SIMBA, and ASTRID suites of camels, and search for the parameters that minimize the χ2 with respect to the observations. For the IllustrisTNG suite, we find the following constraints for the cosmological parameters: $\Omega _{\textrm {m}} = 0.27_{-0.05}^{+0.01}$, $\sigma _{8} = 0.83_{-0.11}^{+0.08}$, and $S_{8} = 0.78_{-0.09}^{+0.03}$, which are consistent within 1σ with the results from the 9-yr WMAP observations. SN feedback-related astrophysical parameters, which describe the departure of outflow wind energy per unit star formation rate and wind velocity from the reference IllustrisTNG simulations, assume the following values: $A_{\textrm {SN1}} = 0.48_{-0.16}^{+0.25}$ and $A_{\textrm {SN2}} = 1.21_{-0.34}^{+0.03}$, respectively. Therefore, simulations with a lower value of outflow wind energy per unit star formation rate with respect to the reference illustrisTNG simulation better reproduce the observations. Variation of active galactic nucleus feedback parameters, on the other hand, shows negligible effects on the scaling relation trends in the mass range probed. Simulations based on SIMBA and ASTRID suites predict central dark matter masses substantially larger than those observed in real galaxies, which can be reconciled with observations only by requiring values of Ωm inconsistent with cosmological constraints for SIMBA, or simulations characterized by unrealistic galaxy mass distributions for ASTRID.

A test of linearity of the ratio of dark matter to baryonic matter in galaxy clusters

We search for a linearity in the ratio of dark matter to baryonic matter as a function of radius for galaxy clusters, motivated by a recent result by Lovas (arXiv:2206.11431), who has discovered such a linearity for a diverse suite of galaxies in the SPARC sample. For our analysis, we used a sample of 54 non-cool core clusters from the HIFLUGCS sample. We do not find any evidence for a linear trend in the aforementioned ratio as a function of radius for individual clusters. We then repeat this analysis for the stacked sample, which also does not show this linearity. Therefore, the linear scaling found by Lovas is not a universal property of dark matter haloes at all scales.

Galactic mass-to-light ratios with superfluid dark matter

Context. We make rotation curve fits to test the superfluid dark matter model. Aims. In addition to verifying that the resulting fits match the rotation curve data reasonably well, we aim to evaluate how satisfactory they are with respect to two criteria, namely, how reasonable the resulting stellar mass-to-light ratios are and whether the fits end up in the regime of superfluid dark matter where the model resembles modified Newtonian dynamics (MOND). Methods. We fitted the superfluid dark matter model to the rotation curves of 169 galaxies in the SPARC sample. Results. We found that the mass-to-light ratios obtained with superfluid dark matter are generally acceptable in terms of stellar populations. However, the best-fit mass-to-light ratios have an unnatural dependence on the size of the galaxy in that giant galaxies have systematically lower mass-to-light ratios than dwarf galaxies. A second finding is that the superfluid often fits the rotation curves best in the regime where the superfluid’s force cannot resemble that of MOND without adjusting a boundary condition separately for each galaxy. In that case, we can no longer expect superfluid dark matter to reproduce the phenomenologically observed scaling relations that make MOND appealing. If, on the other hand, we consider only solutions whose force approximates MOND well, then the total mass of the superfluid is in tension with gravitational lensing data. Conclusions. We conclude that even the best fits with superfluid dark matter are still unsatisfactory for two reasons. First, the resulting stellar mass-to-light ratios show an unnatural trend with galaxy size. Second, the fits do not end up in the regime that automatically resembles MOND, and if we force the fits to do so, the total dark matter mass is in tension with strong lensing data.

MIGHTEE-H i: the baryonic Tully–Fisher relation over the last billion years

ABSTRACT Using a sample of 67 galaxies from the MeerKAT International GigaHertz Tiered Extragalactic Exploration Survey Early Science data, we study the H i-based baryonic Tully–Fisher relation (bTFr), covering a period of ∼1 billion years (0 ≤ z ≤ 0.081). We consider the bTFr based on two different rotational velocity measures: the width of the global H i profile and Vout, measured as the outermost rotational velocity from the resolved H i rotation curves. Both relations exhibit very low intrinsic scatter orthogonal to the best-fitting relation (σ⊥ = 0.07 ± 0.01), comparable to the SPARC sample at z ≃ 0. The slopes of the relations are similar and consistent with the z ≃ 0 studies ($3.66^{+0.35}_{-0.29}$ for W50 and $3.47^{+0.37}_{-0.30}$ for Vout). We find no evidence that the bTFr has evolved over the last billion years, and all galaxies in our sample are consistent with the same relation independent of redshift and the rotational velocity measure. Our results set-up a reference for all future studies of the H i-based bTFr as a function of redshift that will be conducted with the ongoing deep SKA pathfinders surveys.

Using the Baryonic Tully–Fisher Relation to Measure Ho

Abstract We explore the use of the baryonic Tully–Fisher relation (bTFR) as a new distance indicator. Advances in near-IR imaging and stellar population models, plus precise rotation curves, have reduced the scatter in the bTFR such that distance is the dominant source of uncertainty. Using 50 galaxies with accurate distances from Cepheids or the tip magnitude of the red giant branch, we calibrate the bTFR on a scale independent of H o . We then apply this calibrated bTFR to 95 independent galaxies from the SPARC sample, using CosmicFlows-3 velocities, to deduce the local value of H o . We find H o = 75.1 ± 2.3 (stat) ±1.5 (sys) km s−1 Mpc−1.

Dark matter cores and cusps in spiral galaxies and their explanations

We compare proposed solutions to the core vs cusp issue of spiral galaxies, which has also been framed as a diversity problem, and demonstrate that the cuspiness of dark matter halos is correlated with the stellar surface brightness. We compare the rotation curve fits to the SPARC sample from a self-interacting dark matter (SIDM) model, which self-consistently includes the impact of baryons on the halo profile, and hydrodynamical N-body simulations with cold dark matter (CDM) . The SIDM model predicts a strong correlation between the core size and the stellar surface density, and it provides the best global fit to the data. The CDM simulations without strong baryonic feedback effects fail to explain the large dark matter cores seen in low surface brightness galaxies. On the other hand, with strong feedback, CDM simulations do not produce galaxy analogs with high stellar and dark matter densities, and therefore they have trouble in explaining the rotation curves of high surface brightness galaxies. This implies that current feedback implementations need to be modified. We also explicitly show how the concentration-mass and stellar-to-halo mass relations together lead to a radial acceleration relation (RAR) in an averaged sense, and reiterate the point that the RAR does not capture the diversity of galaxy rotation curves in the inner regions. These results make a strong case for SIDM as the explanation for the cores and cusps of field galaxies.

Does dark matter admixed pulsar exist?

In this paper, we have considered a twofluid model assuming that the pulsars are made of ordinary matter admixed with dark matter.Contribution of dark matter comes from the fitting of the rotation curves of the SPARC sample of galaxies[95]. For this we have investigated the dark matter based on the Singular Isothermal Sphere (SIS) dark matter density profile in the galactic halo region. Considering this twofluid model, we have studied the physical features of the pulsars present in different galaxy in details. Here, we compute the probable radii, compactness (u) and surface red-shift (Zs) of the four pulsars namely : PSR J1748-2021B in NGC 6440B galaxy, PSR J1911-5958A in NGC 6752 galaxy, PSR B1802-07 in NGC 6539 galaxy and PSR J1750-37A in NGC 6441 galaxy.

On the absence of a universal surface density, and a maximum Newtonian acceleration in dark matter haloes: Consequences for MOND

Abstract We study the dark matter (DM) surface density using the SPARC sample and compare it to Donato et al. (2009) result. By means of MCMC method, we infer the best-fitting parameters for each galaxy. We reobtain the scaling relation between the surface density and luminosity, and several other scaling laws relating the dark matter halo properties to that of the galactic disc properties. We conclude, in contrast with Donato et al. (2009), that the dark matter surface density is not a universal (constant) quantity but correlates with the luminosity as well as with other galactic disc properties. A derived posterior probability distribution of ρ 0 r 0 shows that the null hypothesis of constancy is rejected at a very high confidence level. These results leave little room for the claimed universality of dark matter surface density. Since MOND has strong prediction on the surface density (Milgrom, 2009), we compared our result with those predictions, finding that MOND predictions are violated by data. To strengthen the previous result, we compared our results to another prediction of MOND (Milgrom and Sanders, 2005), the existence of a maximum Newtonian dark matter acceleration in the halo. Also in this case, MOND predictions are in contradiction with data. The dark matter Newtonian acceleration correlates with all the previously presented galactic disc properties, and data are distributed outside the bound predicted by Milgrom and Sanders (2005). We also find that the null hypothesis (constancy of DM Newtonian acceleration) is rejected at a very high confidence level.

Reconciling the Diversity and Uniformity of Galactic Rotation Curves with Self-Interacting Dark Matter

Galactic rotation curves exhibit diverse behavior in the inner regions, while obeying an organizing principle, i.e., they can be approximately described by a radial acceleration relation or the Modified Newtonian Dynamics phenomenology. We analyze the rotation curve data from the SPARC sample, and explicitly demonstrate that both the diversity and uniformity are naturally reproduced in a hierarchical structure formation model with the addition of dark matter self-interactions. The required concentrations of the dark matter halos are fully consistent with the concentration-mass relation predicted by the Planck cosmological model. The inferred stellar mass-to-light ($3.6 \mu m$) ratios scatter around $0.5 M_\odot/L_\odot$, as expected from population synthesis models, leading to a tight radial acceleration relation and baryonic Tully-Fisher relation. The inferred stellar-halo mass relation is consistent with the expectations from abundance matching. These results indicate that the inner dark matter halos of galaxies are thermalized due to the self-interactions of dark matter particles.

Constraints on chameleon f(R)-gravity from galaxy rotation curves of the SPARC sample

ABSTRACT In chameleon f(R)-gravity, the fifth force will lead to ‘upturns’ in galaxy rotation curves near the screening radius. The location of the upturn depends on the cosmic background value of the scalar field $\bar{f}_\mathrm{R0}$, as well as the mass, size, and environment of the galaxy. We search for this signature of modified gravity in the SPARC sample of measured rotation curves, using an MCMC technique to derive constraints on $\bar{f}_\mathrm{R0}$. Assuming NFW dark matter haloes and with $\bar{f}_\mathrm{R0}$ freely varying for each galaxy, most galaxies prefer f(R) gravity to ΛCDM, but there is a large spread of inferred $\bar{f}_\mathrm{R0}$ values, inconsistent with a single global value. Requiring instead a consistent $\bar{f}_\mathrm{R0}$ value for the whole sample, models with $\log _{10}|\bar{f}_\mathrm{R0}|\gt -6.1$ are excluded. On the other hand, models in the range $-7.5\lt \log _{10}|\bar{f}_\mathrm{R0}|\lt -6.5$ seem to be favoured with respect to ΛCDM, with a significant peak at −7. However, this signal is largely a result of galaxies for which the f(R) signal is degenerate with the core/cusp problem, and when the NFW profile is replaced with a cored halo profile, ΛCDM gives better fits than any given f(R) model. Thus, we find no convincing evidence of f(R) gravity down to the level of $|\bar{f}_\mathrm{R0}|\sim 6 \times 10^{-8}$, with the caveat that if cored halo density profiles cannot ultimately be explained within ΛCDM, a screened modified gravity theory could possibly provide an alternative solution for the core/cusp problem. However, the f(R) models studied here fall short of achieving this.

Modified Gravity (MOG) fits to observed radial acceleration of SPARC galaxies

The equation of motion in the generally covariant modified gravity (MOG) theory leads, for weak gravitational fields and non-relativistic motion, to a modification of Newton's gravitational acceleration law. In addition to the metric $g_{\mu\nu}$, MOG has a vector field $\phi_\mu$ that couples with gravitational strength to all baryonic matter. The gravitational coupling strength is determined by the MOG parameter $\alpha$, while parameter $\mu$ is the small effective mass of $\phi_\mu$. The MOG acceleration law has been demonstrated to fit a wide range of galaxies, galaxy clusters and the Bullet Cluster and Train Wreck Cluster mergers. For the SPARC sample of rotationally supported spiral and irregular galaxies, McGaugh et al. [24] (MLS) have found a radial acceleration relation (RAR) that relates accelerations derived from galaxy rotation curves to Newtonian accelerations derived from galaxy mass models. Using the same SPARC galaxy data, mass models independently derived from that data, and MOG parameters $\alpha$ and $\mu$ that run with galaxy mass, we demonstrate that adjusting galaxy parameters within $\pm 1$-sigma bounds can yield MOG predictions consistent with the given rotational velocity data. Moreover, the same adjusted parameters yield a good fit to the RAR of MLS, with the RAR parameter $a_0=(5.4\pm .3)\times 10^{-11}\,{\rm m/s^2}$.

The radial acceleration relation and dark baryons in MOND

Recent observations of rotationally supported galaxies show a tight correlation between the observed radial acceleration at every radius and the Newtonian acceleration generated by the baryonic mass distribution, the so-called radial acceleration relation (RAR). The rotation curves (RCs) of the SPARC sample of disk galaxies with different morphologies, masses, sizes and gas fractions are investigated in the context of modified Newtonian dynamics (MOND). We include the effect of cold dark baryons by scaling the measured mass in the atomic form by a factor of $c$ in the mass budget of galaxies. In addition to the standard interpolating function, we also fit the RCs and the RAR with the empirical RAR-inspired interpolating function. Slightly better fits for about $47\%$ of galaxies in our sample are achieved in the presence of dark baryons ($c>1$) with the mean value of $c = 2.4\pm 1.3$. Although the MOND fits are not significantly improved by including dark baryons, it results in a decrease in the characteristic acceleration $g_\dag$ by $40\%$. We find no correlation between the MOND critical acceleration $a_0$ and the central surface brightness of the stellar disk, $\mu_{3.6}$. This supports $a_0$ being a universal constant for all galaxies.

The baryonic Tully–Fisher relation for different velocity definitions and implications for galaxy angular momentum

We study the baryonic Tully-Fisher relation (BTFR) at z=0 using 153 galaxies from the SPARC sample. We consider different definitions of the characteristic velocity from HI and H-alpha rotation curves, as well as HI line-widths from single-dish observations. We reach the following results: (1) The tightest BTFR is given by the mean velocity along the flat part of the rotation curve. The orthogonal intrinsic scatter is extremely small (6%) and the best-fit slope is 3.85+/-0.09, but systematic uncertainties may drive the slope from 3.5 to 4.0. Other velocity definitions lead to BTFRs with systematically higher scatters and shallower slopes. (2) We provide statistical relations to infer the flat rotation velocity from HI line-widths or less extended rotation curves (like H-alpha and CO data). These can be useful to study the BTFR from large HI surveys or the BTFR at high redshifts. (3) The BTFR is more fundamental than the relation between angular momentum and galaxy mass (the Fall relation). The Fall relation has about 7 times more scatter than the BTFR, which is merely driven by the scatter in the mass-size relation of galaxies. The BTFR is already the "fundamental plane" of galaxy discs: no value is added with a radial variable as a third parameter.

Uncorrelated velocity and size residuals across galaxy rotation curves

The mass--velocity--size relation of late-type galaxies decouples into independent correlations between mass and velocity (the Tully-Fisher relation), and between mass and size. This behaviour is different to early-type galaxies which lie on a Fundamental Plane. We study the coupling of the Tully-Fisher and mass-size relations in observations (the SPARC sample) and in empirical galaxy formation models based on halo abundance matching, and rotation curve fits with a hydrodynamically motivated halo profile. We systematically investigate the correlation coefficient between the Tully-Fisher residuals $\Delta V_r$ and mass-size residuals $\Delta R$ as a function of the radius $r$ at which the velocity is measured, and thus present the $\Delta V_r-\Delta R$ relation across rotation curves. We find no significant correlation in either the data or models for any $r$, aside from $r \ll R_\text{eff}$ where baryonic mass dominates. We show that this implies an anticorrelation between galaxy size and halo concentration (or halo mass) at fixed baryonic mass, and provides evidence against the hypothesis that galaxy and halo specific angular momentum are proportional. Finally, we study the $\Delta V_r-\Delta R$ relations produced by the baryons and dark matter separately by fitting halo profiles to the rotation curves. The balance between these components illustrates the "disk-halo conspiracy" required for no overall correlation.

On the Dwarf Galaxy Rotation Curve Diversity Problem

In this paper, we show how baryonic physics can solve the problem of the striking diversity in dwarf galaxy rotation curves shapes. To this aim, we compare the distribution of galaxies of the SPARC sample, in the plane V2kpc-VRlast (V2kpc being the galaxy rotation velocity at 2 kpc, and VRlast being the outermost one), with that of galaxies that we simulated, taking account of baryonic effects. The scatter in the rotation curves in the V2kpc-VRlast plane, as well as the trend of the SPARC sample, and our simulated galaxy distribution is in good agreement. The solution of the “diversity” problem lies in the ability of the baryonic process to produce non-self-similar haloes, contrary to DM-only simulations. We show also that baryonic effects can reproduce the rotation curves of galaxies such as IC2574, which are characterized by a slow rise in radius. A solution to the diversity problem can be obtained taking the baryon physics effects appropriately into account.