Showing 1–50 of 179 results for author: McGaugh, S

Accelerated Structure Formation: the Early Emergence of Massive Galaxies and Clusters of Galaxies

Authors: Stacy S. McGaugh, James M. Schombert, Federico Lelli, Jay Franck

Abstract: Galaxies in the early universe appear to have grown too big too fast, assembling into massive, monolithic objects more rapidly than anticipated in the hierarchical $Λ$CDM structure formation paradigm. The available data are consistent with there being a population of massive galaxies that form early ($z \gtrsim 10$) and follow an approximately exponential star formation history with a short (… ▽ More Galaxies in the early universe appear to have grown too big too fast, assembling into massive, monolithic objects more rapidly than anticipated in the hierarchical $Λ$CDM structure formation paradigm. The available data are consistent with there being a population of massive galaxies that form early ($z \gtrsim 10$) and follow an approximately exponential star formation history with a short ($\lesssim 1$ Gyr) e-folding timescale on the way to becoming massive ($M_* \approx 10^{11}\;\mathrm{M}_{\odot}$) galaxies by $z = 0$, consistent with the traditional picture for the evolution of giant elliptical galaxies. Observations of the kinematics of spiral galaxies as a function of redshift similarly show that massive disks and their scaling relations were in place at early times, indicating a genuine effect in mass that cannot be explained as a quirk of luminosity evolution. That massive galaxies could form by $z = 10$ was explicitly predicted in advance by MOND. We discuss some further predictions of MOND, such as the early emergence of clusters of galaxies and the cosmic web. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: 23 pages, 11 figures. Submitted to ApJ

Indefinitely Flat Circular Velocities and the Baryonic Tully-Fisher Relation from Weak Lensing

Authors: Tobias Mistele, Stacy McGaugh, Federico Lelli, James Schombert, Pengfei Li

Abstract: We use a new deprojection formula to infer the gravitational potential around isolated galaxies from weak gravitational lensing. The results imply circular velocity curves that remain flat for hundreds of kpc, greatly extending the classic result from 21 cm observations. Indeed, there is no clear hint of a decline out to 1 Mpc, well beyond the expected virial radii of dark matter halos. Binning th… ▽ More We use a new deprojection formula to infer the gravitational potential around isolated galaxies from weak gravitational lensing. The results imply circular velocity curves that remain flat for hundreds of kpc, greatly extending the classic result from 21 cm observations. Indeed, there is no clear hint of a decline out to 1 Mpc, well beyond the expected virial radii of dark matter halos. Binning the data by mass reveals a correlation with the flat circular speed that closely agrees with the Baryonic Tully-Fisher Relation known from kinematic data. These results apply to both early and late type galaxies, indicating a common universal behavior. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: 14 pages, 4 figures; accepted for publication in ApJL

Journal ref: ApJL 969 L3 (2024)

MHONGOOSE discovery of a gas-rich low-surface brightness galaxy in the Dorado Group

Authors: F. M. Maccagni, W. J. G. de Blok, P. E. Mancera Piña, R. Ragusa, E. Iodice, M. Spavone, S. McGaugh, K. A. Oman, T. A. Oosterloo, B. S. Koribalski, M. Kim, E. A. K. Adams, P. Amram, A. Bosma, F. Bigiel, E. Brinks, L. Chemin, F. Combes, B. Gibson, J. Healy, B. W. Holwerda, G. I. G. Józsa, P. Kamphuis, D. Kleiner, S. Kurapati , et al. (6 additional authors not shown)

Abstract: We present the discovery of a low-mass gas-rich low-surface brightness galaxy in the Dorado Group, at a distance of 17.7 Mpc. Combining deep MeerKAT 21-cm observations from the MeerKAT HI Observations of Nearby Galactic Objects: Observing Southern Emitters (MHONGOOSE) survey with deep photometric images from the VST Early-type Galaxy Survey (VEGAS) we find a stellar and neutral atomic hydrogen (HI… ▽ More We present the discovery of a low-mass gas-rich low-surface brightness galaxy in the Dorado Group, at a distance of 17.7 Mpc. Combining deep MeerKAT 21-cm observations from the MeerKAT HI Observations of Nearby Galactic Objects: Observing Southern Emitters (MHONGOOSE) survey with deep photometric images from the VST Early-type Galaxy Survey (VEGAS) we find a stellar and neutral atomic hydrogen (HI) gas mass of $M_\star = 2.23\times10^6$ M$_\odot$ and $M_{\rm HI}=1.68\times10^6$ M$_\odot$, respectively. This low-surface brightness galaxy is the lowest mass HI detection found in a group beyond the Local Universe ($D\gtrsim 10$ Mpc). The dwarf galaxy has the typical overall properties of gas-rich low surface brightness galaxies in the Local group, but with some striking differences. Namely, the MHONGOOSE observations reveal a very low column density ($\sim 10^{18-19}$ cm$^{-2}$) HI disk with asymmetrical morphology possibly supported by rotation and higher velocity dispersion in the centre. There, deep optical photometry and UV-observations suggest a recent enhancement of the star formation. Found at galactocentric distances where in the Local Group dwarf galaxies are depleted of cold gas (at $390$ projected-kpc distance from the group centre), this galaxy is likely on its first orbit within the Dorado group. We discuss the possible environmental effects that may have caused the formation of the HI disk and the enhancement of star formation, highlighting the short-lived phase (a few hundreds of Myr) of the gaseous disk, before either SF or hydrodynamical forces will deplete the gas of the galaxy. △ Less

Submitted 27 May, 2024; originally announced May 2024.

Comments: 18 pages, 14 figures, 5 tables; accepted for publication in A&A

Satellite group infall into the Milky Way: exploring the Crater-Leo case with new HST proper motions

Authors: Mariana P. Júlio, Marcel S. Pawlowski, Sangmo Tony Sohn, Salvatore Taibi, Roeland P. van der Marel, Stacy S. McGaugh

Abstract: Context. Within $Λ$ Cold Dark Matter ($Λ$CDM) simulations, Milky Way-like galaxies accrete some of their satellite galaxies in groups of 3-5 members rather than individually. It was also suggested that this might be the reason behind the origin of satellite planes. Objects accreted in groups are expected to share similar specific total energy and angular momentum, and also identical orbital planes… ▽ More Context. Within $Λ$ Cold Dark Matter ($Λ$CDM) simulations, Milky Way-like galaxies accrete some of their satellite galaxies in groups of 3-5 members rather than individually. It was also suggested that this might be the reason behind the origin of satellite planes. Objects accreted in groups are expected to share similar specific total energy and angular momentum, and also identical orbital planes and directions. Aims. Looking at observations of Milky Way satellites, the dwarf galaxies Leo II, IV, V, Crater II, and the star cluster Crater 1 were proposed to be a vestige of group infall. The suggested "Crater-Leo group" shows a coherent distance gradient and all these objects align along a great circle on the sky. We use proper motion data to investigate whether the phase-space distribution of the members of the proposed group are indeed consistent with group infall. Methods. To further investigate this possibility, we use Gaia Data Release 3 (DR3) and new Hubble Space Telescope (HST) proper motions, $(μ_{α*}, μ_δ) = (-0.1921 \pm 0.0514, -0.0686 \pm 0.0523)$ mas yr$^{-1}$ for Leo IV and $(μ_{α*}, μ_δ) = (0.1186 \pm 0.1943, -0.1183 \pm 0.1704)$ mas yr$^{-1}$ for Leo V, to derive accurate orbital properties for the proposed group objects. In addition, we explore other possible members of this putative association. Results. Leo II, Leo IV, and Crater 1 show orbital properties consistent with those we predict from assuming group infall. However, our results suggest that Crater II was not accreted with the rest of the objects. If confirmed with increasingly accurate proper motions in the future, the Crater-Leo objects appear to constitute the first identified case of a cosmologically expected, typical group infall event, as opposed to the highly hierarchical Magellanic Cloud system. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: Accepted for publication in A&A on 23-04-2024. 20 pages, 18 figures, 5 tables

The Baryonic Tully-Fisher Relation I: WISE/Spitzer Photometry

Authors: Francis Duey, James Schombert, Stacy McGaugh, Federico Lelli

Abstract: We present WISE W1 photometry of the SPARC (Spitzer Photometry and Accurate Rotation Curves) sample. The baseline of near-IR fluxes is established for use by stellar mass models, a key component to the baryonic Tully-Fisher relation and other kinematic galaxies scaling relations. We focus this paper on determination of the characteristics of the W1 fluxes compared to IRAC 3.6 fluxes, internal accu… ▽ More We present WISE W1 photometry of the SPARC (Spitzer Photometry and Accurate Rotation Curves) sample. The baseline of near-IR fluxes is established for use by stellar mass models, a key component to the baryonic Tully-Fisher relation and other kinematic galaxies scaling relations. We focus this paper on determination of the characteristics of the W1 fluxes compared to IRAC 3.6 fluxes, internal accuracy limitations from photometric techniques, external accuracy by comparison to other work in the literature and the range of W1 to IRAC 3.6 colors. We outline the behavior of SDSS g, W1 and IRAC 3.6 colors with respect to underlying SED features. We also note a previously unknown correlation between WISE colors and the central surface brightness, probably related to the low metallicity of low surface brightness dwarfs. △ Less

Submitted 2 April, 2024; originally announced April 2024.

Comments: Accepted to AJ, 19 pages, 10 figures

MHONGOOSE -- A MeerKAT Nearby Galaxy HI Survey

Authors: W. J. G. de Blok, J. Healy, F. M. Maccagni, D. J. Pisano, A. Bosma, J. English, T. Jarrett, A. Marasco, G. R. Meurer, S. Veronese, F. Bigiel, L. Chemin, F. Fraternali, B. W. Holwerda, P. Kamphuis, H. R. Klöckner, D. Kleiner, A. K. Leroy, M. Mogotsi, K. A. Oman, E. Schinnerer, L. Verdes-Montenegro, T. Westmeier, O. I. Wong, N. Zabel , et al. (35 additional authors not shown)

Abstract: The MHONGOOSE (MeerKAT HI Observations of Nearby Galactic Objects: Observing Southern Emitters) survey maps the distribution and kinematics of the neutral atomic hydrogen (HI) gas in and around 30 nearby star-forming spiral and dwarf galaxies to extremely low HI column densities. The HI column density sensitivity (3 sigma over 16 km/s) ranges from ~ 5 x 10^{17} cm^{-2} at 90'' resolution to ~4 x 1… ▽ More The MHONGOOSE (MeerKAT HI Observations of Nearby Galactic Objects: Observing Southern Emitters) survey maps the distribution and kinematics of the neutral atomic hydrogen (HI) gas in and around 30 nearby star-forming spiral and dwarf galaxies to extremely low HI column densities. The HI column density sensitivity (3 sigma over 16 km/s) ranges from ~ 5 x 10^{17} cm^{-2} at 90'' resolution to ~4 x 10^{19} cm^{-2} at the highest resolution of 7''. The HI mass sensitivity (3 sigma over 50 km/s) is ~5.5 X 10^5 M_sun at a distance of 10 Mpc (the median distance of the sample galaxies). The velocity resolution of the data is 1.4 km/s. One of the main science goals of the survey is the detection of cold, accreting gas in the outskirts of the sample galaxies. The sample was selected to cover a range in HI masses, from 10^7 M_sun to almost 10^{11} M_sun, to optimally sample possible accretion scenarios and environments. The distance to the sample galaxies ranges from 3 to 23 Mpc. In this paper, we present the sample selection, survey design, and observation and reduction procedures. We compare the integrated HI fluxes based on the MeerKAT data with those derived from single-dish measurement and find good agreement, indicating that our MeerKAT observations are recovering all flux. We present HI moment maps of the entire sample based on the first ten percent of the survey data, and find that a comparison of the zeroth- and second-moment values shows a clear separation between the physical properties of the HI in areas with star formation and areas without, related to the formation of a cold neutral medium. Finally, we give an overview of the HI-detected companion and satellite galaxies in the 30 fields, five of which have not previously been catalogued. We find a clear relation between the number of companion galaxies and the mass of the main target galaxy. △ Less

Submitted 6 June, 2024; v1 submitted 2 April, 2024; originally announced April 2024.

Comments: Accepted for publication in Astronomy & Astrophysics

A Distinct Radial Acceleration Relation across Brightest Cluster Galaxies and Galaxy Clusters

Authors: Yong Tian, Chung-Ming Ko, Pengfei Li, Stacy McGaugh, Shemile L. Poblete

Abstract: Recent studies reveal a radial acceleration relation (RAR) in galaxies, which illustrates a tight empirical correlation connecting the observational acceleration and the baryonic acceleration with a characteristic acceleration scale. However, a distinct RAR has been revealed on BCG-cluster scales with a seventeen times larger acceleration scale by the gravitational lensing effect. In this work, we… ▽ More Recent studies reveal a radial acceleration relation (RAR) in galaxies, which illustrates a tight empirical correlation connecting the observational acceleration and the baryonic acceleration with a characteristic acceleration scale. However, a distinct RAR has been revealed on BCG-cluster scales with a seventeen times larger acceleration scale by the gravitational lensing effect. In this work, we systematically explored the acceleration and mass correlations between dynamical and baryonic components in 50 Brightest Cluster Galaxies (BCGs). To investigate the dynamical RAR in BCGs, we derived their dynamical accelerations from the stellar kinematics using the Jeans equation through Abel inversion and adopted the baryonic mass from the SDSS photometry. We explored the spatially resolved kinematic profiles with the largest integral field spectroscopy (IFS) data mounted by the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. Our results demonstrate that the dynamical RAR in BCGs is consistent with the lensing RAR on BCG-cluster scales as well as a larger acceleration scale. This finding may imply that BCGs and galaxy clusters have fundamental differences from field galaxies. We also find a mass correlation, but it is less tight than the acceleration correlation. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Comments: 10 pages, 6 figures

Journal ref: A&A, 684, A180 (2024)

BST1047+1156: A (Failing) Ultradiffuse Tidal Dwarf in the Leo I Group

Authors: J. Christopher Mihos, Patrick R. Durrell, Aaron E. Watkins, Stacy S. McGaugh, John J. Feldmeier

Abstract: We use deep Hubble Space Telescope imaging to study the resolved stellar populations in BST1047+1156, a gas-rich, ultradiffuse dwarf galaxy found in the intragroup environment of the Leo I galaxy group. While our imaging reaches approximately two magnitudes below the tip of the red giant branch at the Leo I distance of 11 Mpc, we find no evidence for an old red giant sequence that would signal an… ▽ More We use deep Hubble Space Telescope imaging to study the resolved stellar populations in BST1047+1156, a gas-rich, ultradiffuse dwarf galaxy found in the intragroup environment of the Leo I galaxy group. While our imaging reaches approximately two magnitudes below the tip of the red giant branch at the Leo I distance of 11 Mpc, we find no evidence for an old red giant sequence that would signal an extended star formation history for the object. Instead, we clearly detect the red and blue helium burning sequences of its stellar populations, as well as the fainter blue main sequence, all indicative of a recent burst of star formation having taken place over the past 50--250 Myr. Comparing to isochrones for young metal-poor stellar populations, we infer this post-starburst population to be moderately metal poor, with metallicity [M/H] in the range -1 to -1.5. The combination of a young, moderately metal-poor post starburst population and no old stars motivates a scenario in which BST1047 was recently formed during a weak burst of star formation in gas that was tidally stripped from the outskirts of the neighboring massive spiral M96. BST1047's extremely diffuse nature, lack of ongoing star formation, and disturbed HI morphology all argue that it is a transitory object, a "failing tidal dwarf" in the process of being disrupted by interactions within the Leo I group. Finally, in the environment surrounding BST1047, our imaging also reveals the old, metal-poor ([M/H]=-1.3 +/- 0.2) stellar halo of M96 at a projected radius of 50 kpc. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: 15 pages, 9 figures, accepted for publication in The Astrophysical Journal

Discord in Concordance Cosmology and Anomalously Massive Early Galaxies

Authors: Stacy McGaugh

Abstract: Cosmological parameters are constrained by a wide variety of observations. We examine the concordance diagram for modern measurements of the Hubble constant, the shape parameter from large scale structure, the cluster baryon fraction, and the age of the universe, all from non-CMB data. There is good agreement for $H_0 = 73.24 \pm 0.38\;\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{Mpc}^{-1}$ and… ▽ More Cosmological parameters are constrained by a wide variety of observations. We examine the concordance diagram for modern measurements of the Hubble constant, the shape parameter from large scale structure, the cluster baryon fraction, and the age of the universe, all from non-CMB data. There is good agreement for $H_0 = 73.24 \pm 0.38\;\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{Mpc}^{-1}$ and $Ω_m = 0.237 \pm 0.015$. This concordance value is indistinguishable from the WMAP3 cosmology but is not consistent with that of Planck: there is a tension in $Ω_m$ as well as $H_0$. These tensions have emerged as progressively higher multipoles have been incorporated into CMB fits. This temporal evolution is suggestive of a systematic effect in the analysis of CMB data at fine angular scales, and may be related to the observation of unexpectedly massive galaxies at high redshift. These are overabundant relative to $Λ$CDM predictions by an order of magnitude at $z > 7$. Such massive objects are anomalous and could cause gravitational lensing of the surface of last scattering in excess of the standard calculation made in CMB fits, potentially skewing the best-fit cosmological parameters and contributing to the Hubble tension. △ Less

Submitted 16 January, 2024; v1 submitted 5 December, 2023; originally announced December 2023.

Comments: 13 pages, 6 figures. Expanded figures at the suggestion of referee. Universe, in press

Journal ref: 2024, Universe, 10, 48

Radial acceleration relation of galaxies with joint kinematic and weak-lensing data

Authors: Tobias Mistele, Stacy McGaugh, Federico Lelli, James Schombert, Pengfei Li

Abstract: We combine kinematic and gravitational lensing data to construct the Radial Acceleration Relation (RAR) of galaxies over a large dynamic range. We improve on previous weak-lensing studies in two ways. First, we compute stellar masses using the same stellar population model as for the kinematic data. Second, we introduce a new method for converting excess surface density profiles to radial accelera… ▽ More We combine kinematic and gravitational lensing data to construct the Radial Acceleration Relation (RAR) of galaxies over a large dynamic range. We improve on previous weak-lensing studies in two ways. First, we compute stellar masses using the same stellar population model as for the kinematic data. Second, we introduce a new method for converting excess surface density profiles to radial accelerations. This method is based on a new deprojection formula which is exact, computationally efficient, and gives smaller systematic uncertainties than previous methods. We find that the RAR inferred from weak-lensing data smoothly continues that inferred from kinematic data by about $2.5\,\mathrm{dex}$ in acceleration. Contrary to previous studies, we find that early- and late-type galaxies lie on the same joint RAR when a sufficiently strict isolation criterion is adopted and their stellar and gas masses are estimated consistently with the kinematic RAR. △ Less

Submitted 8 April, 2024; v1 submitted 23 October, 2023; originally announced October 2023.

Comments: 41 pages, 16 figures, published version

Journal ref: JCAP04(2024)020

Molecular and Ionized Gas in Tidal Dwarf Galaxies: The Spatially Resolved Star-Formation Relation

Authors: Navyasree Kovakkuni, Federico Lelli, Pierre-alain Duc, Médéric Boquien, Jonathan Braine, Elias Brinks, Vassilis Charmandaris, Francoise Combes, Jeremy Fensch, Ute Lisenfeld, Stacy McGaugh, J. Chris Mihos, Marcel. S. Pawlowski, Yves. Revaz, Peter. M. Weilbacher

Abstract: Tidal dwarf galaxies (TDGs) are low-mass objects that form within tidal and/or collisional debris ejected from more massive interacting galaxies. We use CO($1-0$) observations from ALMA and integral-field spectroscopy from MUSE to study molecular and ionized gas in three TDGs: two around the collisional galaxy NGC 5291 and one in the late-stage merger NGC 7252. The CO and H$α$ emission is more com… ▽ More Tidal dwarf galaxies (TDGs) are low-mass objects that form within tidal and/or collisional debris ejected from more massive interacting galaxies. We use CO($1-0$) observations from ALMA and integral-field spectroscopy from MUSE to study molecular and ionized gas in three TDGs: two around the collisional galaxy NGC 5291 and one in the late-stage merger NGC 7252. The CO and H$α$ emission is more compact than the HI emission and displaced from the HI dynamical center, so these gas phases cannot be used to study the internal dynamics of TDGs. We use CO, HI, and H$α$ data to measure the surface densities of molecular gas ($Σ_{\rm mol}$), atomic gas ($Σ_{\rm atom}$) and star-formation rate ($Σ_{\rm SFR}$), respectively. We confirm that TDGs follow the same spatially integrated $Σ_{\rm SFR}-Σ_{\rm gas}$ relation of regular galaxies, where $Σ_{\rm gas} = Σ_{\rm mol} + Σ_{\rm atom}$, even though they are HI dominated. We find a more complex behaviour in terms of the spatially resolved $Σ_{\rm SFR}-Σ_{\rm mol}$ relation on sub-kpc scales. The majority ($\sim$60$\%$) of SF regions in TDGs lie on the same $Σ_{\rm SFR}-Σ_{\rm mol}$ relation of normal spiral galaxies but show a higher dispersion around the mean. The remaining fraction of SF regions ($\sim$40$\%$) lie in the starburst region and are associated with the formation of massive super star clusters, as shown by Hubble Space Telescope images. We conclude that the local SF activity in TDGs proceeds in a hybrid fashion, with some regions comparable to normal spiral galaxies and others to extreme starbursts. △ Less

Submitted 12 September, 2023; originally announced September 2023.

Comments: 11 pages, 6 figures, 3 tables, Accepted for publication in MNRAS

The oblateness of dark matter halos of nearby galaxies and its correlation with gas mass fractions

Authors: Mousumi Das, Roger Ianjamasimanana, Stacy McGaugh, James Schombert, K. S. Dwarakanath

Abstract: We present a method to measure the the oblateness parameter q of the dark matter halos of gas rich galaxies that have extended HI disks. We have applied our model to a sample of 20 nearby galaxies that are gas rich and close to face-on, of which 6 are large disk galaxies, 8 have moderate stellar masses and 6 are low surface brightness (LSB) dwarf galaxies. We have used the stacked HI velocity disp… ▽ More We present a method to measure the the oblateness parameter q of the dark matter halos of gas rich galaxies that have extended HI disks. We have applied our model to a sample of 20 nearby galaxies that are gas rich and close to face-on, of which 6 are large disk galaxies, 8 have moderate stellar masses and 6 are low surface brightness (LSB) dwarf galaxies. We have used the stacked HI velocity dispersion and HI surface densities to derive q in the outer disk regions. Our most important result is that gas dominated galaxies (such as LSB dwarfs) that have M(gas)/M(baryons)>0.5 have oblate halos (q<0.55), whereas stellar dominated galaxies have a range of q values from 0.2 to 1.3. We also find a significant positive correlation between q and stellar mass, which indicates that galaxies with massive stellar disks have a higher probability of having halos that are spherical or slightly prolate, whereas low mass galaxies preferably have oblate halos. We briefly also discuss how the halo shape affects the disks of galaxies, especially the oblate halos. △ Less

Submitted 3 August, 2023; originally announced August 2023.

Comments: 6 pages, 3 figures. Appeared in the proceedings of Dynamical Masses of Local Group Galaxies: IAU Symposium 379, held in Potsdam, March 2023

Local Group Galaxies from an External Perspective

Authors: Stacy McGaugh

Abstract: I discuss Local Group galaxies from the perspective of external galaxies that define benchmark scaling relations. Making use of this information leads to a model for the Milky Way that includes bumps and wiggles due to spiral arms. This model reconciles the terminal velocities observed in the interstellar medium with the rotation curve derived from stars, correctly predicts the gradual decline of… ▽ More I discuss Local Group galaxies from the perspective of external galaxies that define benchmark scaling relations. Making use of this information leads to a model for the Milky Way that includes bumps and wiggles due to spiral arms. This model reconciles the terminal velocities observed in the interstellar medium with the rotation curve derived from stars, correctly predicts the gradual decline of the outer rotation curve ($dV/dR = -1.7\;\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{kpc}^{-1}$), and extrapolates well out to 50 kpc. Rotationally supported Local Group galaxies are in excellent agreement with the baryonic Tully-Fisher relation. Pressure supported dwarfs that are the most likely to be in dynamical equilibrium also align with this relation. Local Group galaxies thus appear to be normal members of the low redshift galaxy population. There is, however, a serious tension between the dynamical masses of the Milky Way and M31 ($M_{200} \approx 1.4$ and $1.6 \times 10^{12}\;\mathrm{M}_{\odot}$, respectively) and those expected from the stellar mass-halo mass relation of abundance matching ($M_{200} \approx 3$ and $20 \times 10^{12}\;\mathrm{M}_{\odot}$, respectively). △ Less

Submitted 1 May, 2023; originally announced May 2023.

Comments: 6 pages, 4 figures. Contribution to IAU Symposium 379: Dynamical Masses of Local Group Galaxies

Measuring galaxy cluster mass profiles into the low acceleration regime with galaxy kinematics

Authors: Pengfei Li, Yong Tian, Mariana P. Júlio, Marcel S. Pawlowski, Federico Lelli, Stacy S. McGaugh, James M. Schombert, Justin I. Read, Po-Chieh Yu, Chung-Ming Ko

Abstract: We probe the dynamical mass profiles of 10 galaxy clusters from the HIghest X-ray FLUx Galaxy Cluster Sample (HIFLUGCS) using galaxy kinematics. We numerically solve the spherical Jeans equation, and parameterize the dynamical mass profile and the galaxy velocity anisotropy profile using two general functions to ensure that our results are not biased towards any specific model. The mass-velocity a… ▽ More We probe the dynamical mass profiles of 10 galaxy clusters from the HIghest X-ray FLUx Galaxy Cluster Sample (HIFLUGCS) using galaxy kinematics. We numerically solve the spherical Jeans equation, and parameterize the dynamical mass profile and the galaxy velocity anisotropy profile using two general functions to ensure that our results are not biased towards any specific model. The mass-velocity anisotropy degeneracy is ameliorated by using two "virial shape parameters" that depend on the fourth moment of velocity distribution. The resulting velocity anisotropy estimates consistently show a nearly isotropic distribution in the inner regions, with an increasing radial anisotropy towards large radii. We compare our derived dynamical masses with those calculated from X-ray gas data assuming hydrostatic equilibrium, finding that massive and rich relaxed clusters generally present consistent mass measurements, while unrelaxed or low-richness clusters have systematically larger total mass than hydrostatic mass by an average of 50\%. This might help alleviate current tensions in the measurement of $σ_8$, but it also leads to cluster baryon fractions below the cosmic value. Finally, our approach probes accelerations as low as $10^{-11}$ m s$^{-2}$, comparable to the outskirts of individual late-type galaxies. We confirm that galaxy clusters deviate from the radial acceleration relation defined by galaxies. △ Less

Submitted 14 June, 2023; v1 submitted 16 March, 2023; originally announced March 2023.

Comments: 16 pages, 8 figures, 1 table. A&A accepted version

Journal ref: A&A 677, A24 (2023)

Superfluid dark matter in tension with weak gravitational lensing data

Authors: Tobias Mistele, Stacy McGaugh, Sabine Hossenfelder

Abstract: Superfluid dark matter (SFDM) is a model that promises to reproduce the successes of both particle dark matter on cosmological scales and those of Modified Newtonian Dynamics (MOND) on galactic scales. SFDM reproduces MOND only up to a certain distance from the galactic center, and only for kinematic observables: It does not affect trajectories of light. We test whether this is consistent with a r… ▽ More Superfluid dark matter (SFDM) is a model that promises to reproduce the successes of both particle dark matter on cosmological scales and those of Modified Newtonian Dynamics (MOND) on galactic scales. SFDM reproduces MOND only up to a certain distance from the galactic center, and only for kinematic observables: It does not affect trajectories of light. We test whether this is consistent with a recent analysis of weak gravitational lensing that has probed accelerations around galaxies to unprecedentedly large radii. This analysis found the data to be close to the prediction of MOND, suggesting they might be difficult to fit with SFDM. To investigate this matter, we solved the equations of motion of the model and compared the result to observational data. Our results show that the SFDM model is incompatible with the weak-lensing observations, at least in its current form. △ Less

Submitted 5 September, 2023; v1 submitted 15 March, 2023; originally announced March 2023.

Comments: 20 pages, 12 figures

Journal ref: JCAP09(2023)004

Estimating the Oblateness Of Dark Matter Halos Using Neutral Hydrogen Velocity Dispersion

Authors: Mousumi Das, Roger Ianjamasimanana, Stacy S. McGaugh, James Schombert, K. S. Dwarakanath

Abstract: We derive the oblateness parameter q of the dark matter halo of a sample of gas rich, face-on disk galaxies. We have assumed that the halos are triaxial in shape but their axes in the disk plane (a and b) are equal, so that q=c/a measures the halo flattening. We have used the HI velocity dispersion, derived from the stacked HI emission lines and the disk surface density to determine the disk poten… ▽ More We derive the oblateness parameter q of the dark matter halo of a sample of gas rich, face-on disk galaxies. We have assumed that the halos are triaxial in shape but their axes in the disk plane (a and b) are equal, so that q=c/a measures the halo flattening. We have used the HI velocity dispersion, derived from the stacked HI emission lines and the disk surface density to determine the disk potential and the halo shape at the R25 and 1.5R25 radii. We have applied our model to 20 nearby galaxies, of which 6 are large disk galaxies with M(stellar)>10^10 solar mass, 8 have moderate stellar masses and 6 are low surface brightness dwarf galaxies. Our most important result is that gas rich galaxies that have M(gas)/M(baryons)>0.5 have oblate halos (q < 0.55), whereas stellar dominated galaxies have a range of q values, ranging from 0.21+-0.07 in NGC4190 to 1.27+-0.61 in NGC5194. Our results also suggest a positive correlation between the stellar mass and the halo oblateness q, which indicates that galaxies with massive stellar disks have a higher probability of having halos that are spherical or slightly prolate, whereas low mass galaxies have oblate halos (q < 0.55). △ Less

Submitted 9 March, 2023; originally announced March 2023.

Comments: 9 pages, 4 figures. Accepted for publication in ApJL

Cold gas disks in main-sequence galaxies at cosmic noon: Low turbulence, flat rotation curves, and disk-halo degeneracy

Authors: Federico Lelli, Zhi-Yu Zhang, Thomas G. Bisbas, Lingrui Lin, Padelis Papadopoulos, James M. Schombert, Enrico Di Teodoro, Antonino Marasco, Stacy S. McGaugh

Abstract: We study the dynamics of cold molecular gas in two main-sequence galaxies at cosmic noon (zC-488879 at $z\simeq1.47$ and zC-400569 at $z\simeq2.24$) using new high-resolution ALMA observations of multiple $^{12}$CO transitions. For zC-400569 we also re-analyze high-quality H$α$ data from the SINS/zC-SINF survey. We find that (1) Both galaxies have regularly rotating CO disks and their rotation cur… ▽ More We study the dynamics of cold molecular gas in two main-sequence galaxies at cosmic noon (zC-488879 at $z\simeq1.47$ and zC-400569 at $z\simeq2.24$) using new high-resolution ALMA observations of multiple $^{12}$CO transitions. For zC-400569 we also re-analyze high-quality H$α$ data from the SINS/zC-SINF survey. We find that (1) Both galaxies have regularly rotating CO disks and their rotation curves are flat out to $\sim$8 kpc contrary to previous results pointing to outer declines in the rotation speed $V_{\rm rot}$; (2) The intrinsic velocity dispersions are low ($σ_{\rm CO}\lesssim15$ km/s for CO and $σ_{\rm Hα}\lesssim37$ km/s for H$α$) and imply $V_{\rm rot}/σ_{\rm CO}\gtrsim17-22$ yielding no significant pressure support; (3) Mass models using HST images display a severe disk-halo degeneracy: models with inner baryon dominance and models with "cuspy" dark matter halos can fit the rotation curves equally well due to the uncertainties on stellar and gas masses; (4) Milgromian dynamics (MOND) can successfully fit the rotation curves with the same acceleration scale $a_0$ measured at $z\simeq0$. The question of the amount and distribution of dark matter in high-$z$ galaxies remains unsettled due to the limited spatial extent of the available kinematic data; we discuss the suitability of various emission lines to trace extended rotation curves at high $z$. Nevertheless, the properties of these two high-$z$ galaxies (high $V_{\rm rot}/σ_{\rm V}$ ratios, inner rotation curve shapes, bulge-to-total mass ratios) are remarkably similar to those of massive spirals at $z\simeq0$, suggesting weak dynamical evolution over more than 10 Gyr of the Universe's lifetime. △ Less

Submitted 25 May, 2023; v1 submitted 31 January, 2023; originally announced February 2023.

Comments: 18 pages, 11 figures, 4 tables, 2 appendices. Minor revisions to match the published version in Astronomy and Astrophysics

Journal ref: A&A 672, A106 (2023)

Aether scalar tensor theory confronted with weak lensing data at small accelerations

Authors: Tobias Mistele, Stacy McGaugh, Sabine Hossenfelder

Abstract: The recently proposed aether scalar tensor (AeST) model reproduces both the successes of particle dark matter on cosmological scales and those of modified Newtonian dynamics (MOND) on galactic scales. But the AeST model reproduces MOND only up to a certain maximum galactocentric radius. Since MOND is known to fit very well to observations at these scales, this raises the question of whether the Ae… ▽ More The recently proposed aether scalar tensor (AeST) model reproduces both the successes of particle dark matter on cosmological scales and those of modified Newtonian dynamics (MOND) on galactic scales. But the AeST model reproduces MOND only up to a certain maximum galactocentric radius. Since MOND is known to fit very well to observations at these scales, this raises the question of whether the AeST model comes into tension with data. We tested whether or not the AeST model is in conflict with observations using a recent analysis of data for weak gravitational lensing. We solved the equations of motion of the AeST model, analyzed the solutions' behavior, and compared the results to observational data. The AeST model shows some deviations from MOND at the radii probed by weak gravitational lensing. The data show no clear indication of these predicted deviations. △ Less

Submitted 4 July, 2023; v1 submitted 9 January, 2023; originally announced January 2023.

Comments: 18 pages, 14 figures, forthcoming in A&A

Journal ref: A&A 676, A100 (2023)

Deep Narrowband Photometry of the M101 Group: Strong-Line Abundances of 720 HII Regions

Authors: Ray Garner III, J. Christopher Mihos, Paul Harding, Aaron E. Watkins, Stacy S. McGaugh

Abstract: We present deep, narrowband imaging of the nearby spiral galaxy M101 and its satellites to analyze the oxygen abundances of their HII regions. Using CWRU's Burrell Schmidt telescope, we add to the narrowband dataset of the M101 Group, consisting of H$α$, H$β$, and [OIII] emission lines, the blue [OII]$λ$3727 emission line for the first time. This allows for complete spatial coverage of the oxygen… ▽ More We present deep, narrowband imaging of the nearby spiral galaxy M101 and its satellites to analyze the oxygen abundances of their HII regions. Using CWRU's Burrell Schmidt telescope, we add to the narrowband dataset of the M101 Group, consisting of H$α$, H$β$, and [OIII] emission lines, the blue [OII]$λ$3727 emission line for the first time. This allows for complete spatial coverage of the oxygen abundance of the entire M101 Group. We used the strong-line ratio $R_{23}$ to estimate oxygen abundances for the HII regions in our sample, utilizing three different calibration techniques to provide a baseline estimate of the oxygen abundances. This results in ~650 HII regions for M101, 10 HII regions for NGC 5477, and ~60 HII regions for NGC 5474, the largest sample for this Group to date. M101 shows a strong abundance gradient while the satellite galaxies present little or no gradient. There is some evidence for a flattening of the gradient in M101 beyond $R \sim 14 \text{ kpc}$. Additionally, M101 shows signs of azimuthal abundance variations to the west and southwest. The radial and azimuthal abundance variations in M101 are likely explained by an interaction it had with its most massive satellite NGC 5474 ~300 Myr ago combined with internal dynamical effects such as corotation. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: 24 pages, 14 figures, 5 tables. Accepted to ApJ

Testing Modified Gravity Theories with Numerical Solutions of the External Field Effect in Rotationally Supported Galaxies

Authors: Kyu-Hyun Chae, Federico Lelli, Harry Desmond, Stacy S. McGaugh, James M. Schombert

Abstract: The strong equivalence principle is violated by gravity theories of Milgromian dynamics (MOND) through the action of the external field effect. We test two different Lagrangian theories AQUAL and QUMOND based on their numerical solutions of the external field effect, by comparing two independent estimates of the mean external field strength of the nearby universe: a theory-deduced value from fitti… ▽ More The strong equivalence principle is violated by gravity theories of Milgromian dynamics (MOND) through the action of the external field effect. We test two different Lagrangian theories AQUAL and QUMOND based on their numerical solutions of the external field effect, by comparing two independent estimates of the mean external field strength of the nearby universe: a theory-deduced value from fitting the outer rotation curves of 114 galaxies and an empirical value from the large-scale distribution of cosmic baryons. The AQUAL-deduced external field strength from rotation curves agrees with that from the large-scale cosmic environment, while QUMOND-deduced value is somewhat higher. This suggests that AQUAL is likely to be preferred over QUMOND as an effective non-relativistic limit of a potential relativistic modified gravity theory. △ Less

Submitted 18 October, 2022; v1 submitted 15 September, 2022; originally announced September 2022.

Comments: 6 pages, 2 figures, 1 table (results revised, conclusions strengthened)

Incorporating baryon-driven contraction of dark matter halos in rotation curve fits

Authors: Pengfei Li, Stacy S. McGaugh, Federico Lelli, James M. Schombert, Marcel S. Pawlowski

Abstract: The condensation of baryons within a dark matter (DM) halo during galaxy formation should result in some contraction of the halo as the combined system settles into equilibrium. We quantify this effect on the cuspy primordial halos predicted by DM-only simulations for the baryon distributions observed in the galaxies of the SPARC database. We find that the DM halos of high surface brightness galax… ▽ More The condensation of baryons within a dark matter (DM) halo during galaxy formation should result in some contraction of the halo as the combined system settles into equilibrium. We quantify this effect on the cuspy primordial halos predicted by DM-only simulations for the baryon distributions observed in the galaxies of the SPARC database. We find that the DM halos of high surface brightness galaxies (with $Σ_{\rm eff}\gtrsim100$ $L_\odot$ pc$^{-2}$ at 3.6 $μ$m) experience strong contraction. Halos become more cuspy as a result of compression: the inner DM density slope increases with the baryonic surface mass density. We iteratively fit rotation curves to find the balance between initial halo parameters (constrained by abundance matching), compression, and stellar mass-to-light ratio. The resulting fits often require lower stellar masses than expected for stellar populations, particularly in galaxies with bulges: stellar mass must be reduced to make room for the DM it compresses. This trade off between dark and luminous mass is reminiscent of the cusp-core problem in dwarf galaxies, but occurs in more massive systems: the present-epoch DM halos cannot follow from cuspy primordial halos unless (1) the stellar mass-to-light ratios are systematically smaller than expected from standard stellar population synthesis models, and/or (2) there is a net outward mass redistribution from the initial cusp, even in massive galaxies widely considered to be immune from such effects. △ Less

Submitted 8 August, 2022; originally announced August 2022.

Comments: 13 pages, 9 figures, 1 table. Accepted for publication in A&A

Journal ref: A&A 665, A143 (2022)

The Effect of Adiabatic Compression on Dark Matter Halos and the Radial Acceleration Relation

Authors: Pengfei Li, Stacy S. McGaugh, Federico Lelli, Yong Tian, James M. Schombert, Chung-Ming Ko

Abstract: We use a semi-empirical model to investigate the radial acceleration relation (RAR) in a cold dark matter (CDM) framework. Specifically, we build 80 model galaxies covering the same parameter space as the observed galaxies in the SPARC database, assigning them to dark matter halos using abundance matching and halo mass-concentration relations. We consider several abundance matching relations, find… ▽ More We use a semi-empirical model to investigate the radial acceleration relation (RAR) in a cold dark matter (CDM) framework. Specifically, we build 80 model galaxies covering the same parameter space as the observed galaxies in the SPARC database, assigning them to dark matter halos using abundance matching and halo mass-concentration relations. We consider several abundance matching relations, finding some to be a better match to the kinematic data than others. We compute the unavoidable gravitational interactions between baryons and their dark matter halos, leading to an overall compression of the original NFW halos. Before halo compression, high-mass galaxies approximately lie on the observed RAR whereas low-mass galaxies display up-bending "hooks" at small radii due to DM cusps, making them deviate systematically from the observed relation. After halo compression, the initial NFW halos become more concentrated at small radii, making larger contributions to rotation curves. This increases the total accelerations, moving all model galaxies away from the observed relation. These systematic deviations suggest that the CDM model with abundance matching alone cannot explain the observed RAR. Further effects (e.g., feedback) would need to counteract the compression with precisely the right amount of halo expansion, even in high mass galaxies with deep potential wells where such effects are generally predicted to be negligible. △ Less

Submitted 7 February, 2022; originally announced February 2022.

Comments: 12 pages, 6 figures. Accepted for publication in ApJ

Stellar Mass-to-light Ratios: Composite Bulge+Disk Models and the Baryonic Tully-Fisher Relation

Authors: James Schombert, Stacy McGaugh, Federico Lelli

Abstract: We present stellar population models to calculate the mass-to-light ratio ($Υ_*$) based on galaxy's colors ranging from $GALEX$ FUV to Spitzer IRAC1 at 3.6$μ$m. We present a new composite bulge+disk $Υ_*$ model that considers the varying contribution from bulges and disks based on their optical and near-IR colors. Using these colors, we build plausible star formation histories and chemical enrichm… ▽ More We present stellar population models to calculate the mass-to-light ratio ($Υ_*$) based on galaxy's colors ranging from $GALEX$ FUV to Spitzer IRAC1 at 3.6$μ$m. We present a new composite bulge+disk $Υ_*$ model that considers the varying contribution from bulges and disks based on their optical and near-IR colors. Using these colors, we build plausible star formation histories and chemical enrichment scenarios based on the star formation rate-stellar mass and mass-metallicity correlations for star-forming galaxies. The most accurate prescription is to use the actual colors for the bulge and disk components to constrain $Υ_*$; however, a reasonable bulge+disk model plus total color only introduces 5% more uncertainty. Full bulge+disk $Υ_*$ prescriptions applied to the baryonic TF relation improves the linearity of the correlation, increases the slope and reduces the total scatter by 4%. △ Less

Submitted 4 February, 2022; originally announced February 2022.

Comments: 22 pages, six figures, accepted by AJ

Galactic Mass-to-Light Ratios With Superfluid Dark Matter

Authors: Tobias Mistele, Stacy McGaugh, Sabine Hossenfelder

Abstract: We make rotation curve fits to test the superfluid dark matter model. 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 resem… ▽ More We make rotation curve fits to test the superfluid dark matter model. 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). We fitted the superfluid dark matter model to the rotation curves of 169 galaxies in the SPARC sample. 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. 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. △ Less

Submitted 3 August, 2022; v1 submitted 18 January, 2022; originally announced January 2022.

Comments: 29 pages, 29 figures, published version

Journal ref: A&A 664, A40 (2022)

Testing the Strong Equivalence Principle. II. Relating the External Field Effect in Galaxy Rotation Curves to the Large-Scale Structure of the Universe

Authors: Kyu-Hyun Chae, Harry Desmond, Federico Lelli, Stacy S. McGaugh, James M. Schombert

Abstract: Theories of modified gravity generically violate the strong equivalence principle, so that the internal dynamics of a self-gravitating system in free fall depends on the strength of the external gravitational field (the external field effect). We fit rotation curves (RCs) from the SPARC database with a model inspired by Milgromian dynamics (MOND), which relates the outer shape of a RC to the exter… ▽ More Theories of modified gravity generically violate the strong equivalence principle, so that the internal dynamics of a self-gravitating system in free fall depends on the strength of the external gravitational field (the external field effect). We fit rotation curves (RCs) from the SPARC database with a model inspired by Milgromian dynamics (MOND), which relates the outer shape of a RC to the external Newtonian field from the large-scale baryonic matter distribution through a dimensionless parameter $e_{\rm N}$. We obtain a $>4σ$ statistical detection of the external field effect (i.e. $e_{\rm N}>0$ on average), confirming previous results. We then locate the SPARC galaxies in the cosmic web of the nearby Universe and find a striking contrast in the fitted $e_{\rm N}$ {values} for galaxies in underdense versus overdense regions. Galaxies in an underdense region between 22 and 45 Mpc from the celestial axis in the northern sky have RC fits consistent with $e_{\rm N}\simeq0$, while those in overdense regions adjacent to the CfA2 great wall and the Perseus-Pisces supercluster return $e_{\rm N}$ that are a factor of two larger than the median for SPARC galaxies. We also calculate independent estimates of $e_{\rm N}$ from galaxy survey data and find that they agree with the $e_{\rm N}$ inferred from the RCs within the uncertainties, the chief uncertainty being the spatial distribution of baryons not contained in galaxies or clusters. △ Less

Submitted 10 September, 2021; originally announced September 2021.

Comments: Astrophysical Journal, in press, 13 figures, 3 tables

The Baryonic Tully-Fisher Relation in the Local Group and the Equivalent Circular Velocity of Pressure Supported Dwarfs

Authors: Stacy McGaugh, Federico Lelli, Jim Schombert, Pengfei Li, Tiffany Visgaitis, Kaelee Parker, Marcel Pawlowski

Abstract: We explore the Baryonic Tully-Fisher Relation in the Local Group. Rotationally supported Local Group galaxies adhere precisely to the relation defined by more distant galaxies. For pressure supported dwarf galaxies, we determine the scaling factor $β_c$ that relates their observed velocity dispersion to the equivalent circular velocity of rotationally supported galaxies of the same mass such that… ▽ More We explore the Baryonic Tully-Fisher Relation in the Local Group. Rotationally supported Local Group galaxies adhere precisely to the relation defined by more distant galaxies. For pressure supported dwarf galaxies, we determine the scaling factor $β_c$ that relates their observed velocity dispersion to the equivalent circular velocity of rotationally supported galaxies of the same mass such that $V_o = β_c σ_*$. For a typical mass-to-light ratio $Υ_* = 2\;\mathrm{M}_{\odot}/\mathrm{L}_{\odot}$ in the $V$-band, we find that $β_c = 2$. More generally, $\log β_c = 0.25 \log Υ_* +0.226$. This provides a common kinematic scale relating pressure and rotationally supported dwarf galaxies. △ Less

Submitted 7 September, 2021; originally announced September 2021.

Comments: 14 pages, 4 figures. Accepted for publication in the Astronomical Journal

Mass-Velocity Dispersion Relation in MaNGA Brightest Cluster Galaxies

Authors: Yong Tian, Han Cheng, Stacy S. McGaugh, Chung-Ming Ko, Yun-Hsin Hsu

Abstract: We investigate a kinematic scaling relation between the baryonic mass and the flat velocity dispersion, i.e. mass-velocity dispersion relation (MVDR), from the brightest cluster galaxies (BCGs) to the galaxy clusters. In our studies, the baryonic mass of BCGs is mainly estimated by photometry. The velocity dispersion profiles are explored with the integrated field unit (IFU) by Mapping Nearby Gala… ▽ More We investigate a kinematic scaling relation between the baryonic mass and the flat velocity dispersion, i.e. mass-velocity dispersion relation (MVDR), from the brightest cluster galaxies (BCGs) to the galaxy clusters. In our studies, the baryonic mass of BCGs is mainly estimated by photometry. The velocity dispersion profiles are explored with the integrated field unit (IFU) by Mapping Nearby Galaxies at Apache Point Observatory (MaNGA). For the first time, we reveal two significant results with 54 MaNGA BCGs: (1) the flat velocity dispersion profiles; (2) a tight empirical relation on the BCG-cluster scale together with cluster samples, i.e., MVDR, $\log(M_\mathrm{bar}/M_\odot)=4.1^{+0.1}_{-0.1}\log(σ_{\mathrm{los}}/\mathrm{km}\,\mathrm{s}^{-1})+1.6^{+0.3}_{-0.3}$, with a tiny lognormal intrinsic scatter of $10^{+2}_{-1}\%$. This slope is identical to the acceleration relation in galaxy clusters, which is reminiscent of the spiral galaxies, albeit at a larger characteristic acceleration scale. The residuals of the MVDR represent a Gaussian distribution, displaying no correlations with four properties: baryonic mass, scale length, surface density, and redshift. Notably, the MVDR on the BCG-cluster scale provides a strict test, which disfavors the general prediction of the slope of three in the dark matter model. △ Less

Submitted 19 August, 2021; originally announced August 2021.

Comments: 6 pages, 3 figures, 1 table

Journal ref: The Astrophysical Journal Letters, 917:L24 (6pp), 2021 August 20

Comparing the Inner and Outer Star Forming Complexes in the Nearby Spiral Galaxies NGC 628, NGC 5457 and NGC 6946 using UVIT Observations

Authors: Jyoti Yadav, Mousumi Das, Narendra Nath Patra, K. S. Dwarakanath, P. T. Rahna, Stacy S. McGaugh, James Schombert, Jayant Murthy

Abstract: We present a far-UV (FUV) study of the star-forming complexes (SFCs) in three nearby galaxies using the Ultraviolet Imaging Telescope (UVIT). The galaxies are close to face-on and show significant outer disk star formation. Two of them are isolated (NGC 628, NGC 6946), and one is interacting with distant companions (NGC 5457). We compared the properties of the SFCs inside and outside the optical r… ▽ More We present a far-UV (FUV) study of the star-forming complexes (SFCs) in three nearby galaxies using the Ultraviolet Imaging Telescope (UVIT). The galaxies are close to face-on and show significant outer disk star formation. Two of them are isolated (NGC 628, NGC 6946), and one is interacting with distant companions (NGC 5457). We compared the properties of the SFCs inside and outside the optical radius (R$_{25}$). We estimated the sizes, star formation rates (SFRs), metallicities, and the Toomre Q parameter of the SFCs. We find that the outer disk SFCs are at least ten times smaller in area than those in the inner disk. The SFR per unit area ($Σ_{SFR}$) in both regions have similar mean values, but the outer SFCs have a much smaller range of $Σ_{SFR}$. They are also metal-poor compared to the inner disk SFCs. The FUV emission is well correlated with the neutral hydrogen gas (\HI) distribution and is detected within and near several \HI~holes. Our estimation of the Q parameter in the outer disks of the two isolated galaxies suggests that their outer disks are stable (Q$>$1). However, their FUV images indicate that there is ongoing star formation in these regions. This suggests that there may be some non-luminous mass or dark matter in their outer disks, which increases the disk surface density and supports the formation of local gravitational instabilities. In the interacting galaxy, NGC 5457, the baryonic surface density is sufficient (Q$<$1) to trigger local disk instabilities in the outer disk. △ Less

Submitted 31 March, 2021; originally announced March 2021.

Comments: 24 pages, 8 figures , Accepted in ApJ

Testing Galaxy Formation and Dark Matter with Low Surface Brightness Galaxies

Authors: Stacy McGaugh

Abstract: Galaxies are the basic structural element of the universe; galaxy formation theory seeks to explain how these structures came to be. I trace some of the foundational ideas in galaxy formation, with emphasis on the need for non-baryonic cold dark matter. Many elements of early theory did not survive contact with observations of low surface brightness galaxies, leading to the need for auxiliary hypo… ▽ More Galaxies are the basic structural element of the universe; galaxy formation theory seeks to explain how these structures came to be. I trace some of the foundational ideas in galaxy formation, with emphasis on the need for non-baryonic cold dark matter. Many elements of early theory did not survive contact with observations of low surface brightness galaxies, leading to the need for auxiliary hypotheses like feedback. The failure points often trace to the surprising predictive successes of an alternative to dark matter, the Modified Newtonian Dynamics (MOND). While dark matter models are flexible in accommodating observations, they do not provide the predictive capacity of MOND. If the universe is made of cold dark matter, why does MOND get any predictions right? △ Less

Submitted 8 March, 2021; originally announced March 2021.

Comments: 29 pages (plus 20 pages of references), 8 figures. Invited review for a conference on Dark Matter and Modified Gravity, part of the project `LHC and Gravity'. Includes an introduction accessible to a broad audience and some historical narrative leading to more advanced material

A cautionary tale in fitting galaxy rotation curves with Bayesian techniques: does Newton's constant vary from galaxy to galaxy?

Authors: Pengfei Li, Federico Lelli, Stacy McGaugh, James Schombert, Kyu-Hyun Chae

Abstract: The application of Bayesian techniques to astronomical data is generally non-trivial because the fitting parameters can be strongly degenerated and the formal uncertainties are themselves uncertain. An example is provided by the contradictory claims over the presence or absence of a universal acceleration scale (g$_\dagger$) in galaxies based on Bayesian fits to rotation curves. To illustrate the… ▽ More The application of Bayesian techniques to astronomical data is generally non-trivial because the fitting parameters can be strongly degenerated and the formal uncertainties are themselves uncertain. An example is provided by the contradictory claims over the presence or absence of a universal acceleration scale (g$_\dagger$) in galaxies based on Bayesian fits to rotation curves. To illustrate the situation, we present an analysis in which the Newtonian gravitational constant $G_N$ is allowed to vary from galaxy to galaxy when fitting rotation curves from the SPARC database, in analogy to $g_{\dagger}$ in the recently debated Bayesian analyses. When imposing flat priors on $G_N$, we obtain a wide distribution of $G_N$ which, taken at face value, would rule out $G_N$ as a universal constant with high statistical confidence. However, imposing an empirically motivated log-normal prior returns a virtually constant $G_N$ with no sacrifice in fit quality. This implies that the inference of a variable $G_N$ (or g$_{\dagger}$) is the result of the combined effect of parameter degeneracies and unavoidable uncertainties in the error model. When these effects are taken into account, the SPARC data are consistent with a constant $G_{\rm N}$ (and constant $g_\dagger$). △ Less

Submitted 28 January, 2021; v1 submitted 27 January, 2021; originally announced January 2021.

Comments: Accepted for publication in A&A Letters

Journal ref: A&A 646, L13 (2021)

Anomalous Stellar Populations in LSB Galaxies

Authors: James Schombert, Stacy McGaugh

Abstract: We present new HST WFC3 near-IR observations of the CMD's in two LSB galaxies, F575-3 and F615-1, notable for having no current star formation based on a lack of H$α$ emission. Key features of the near-IR CMD's are resolved, such as the red giant branch (RGB), the asymptotic giant branch (AGB) region and the top of the blue main sequence (bMS). F575-3 has the bluest RGB of any CMD in the literatur… ▽ More We present new HST WFC3 near-IR observations of the CMD's in two LSB galaxies, F575-3 and F615-1, notable for having no current star formation based on a lack of H$α$ emission. Key features of the near-IR CMD's are resolved, such as the red giant branch (RGB), the asymptotic giant branch (AGB) region and the top of the blue main sequence (bMS). F575-3 has the bluest RGB of any CMD in the literature, indicating an extremely low mean metallicity. F615-1 has unusually wide RGB and AGB sequences suggesting multiple episodes of star formation from metal-poor gas, possibly infalling material. Both galaxies have an unusual population of stars to the red of the RGB and lower in luminosity than typical AGB stars. These stars have normal optical colors but abnormal near-IR colors. We suggest that this population of stars might be analogous to local peculiar stars like Be stars with strong near-IR excesses owing to a surrounding disk of hot gas. △ Less

Submitted 17 December, 2020; originally announced December 2020.

Comments: 29 pages, 11 figures, accepted by AJ

Mass-Velocity Dispersion Relation in HIFLUGCS Galaxy Clusters

Authors: Yong Tian, Po-Chieh Yu, Pengfei Li, Stacy S. McGaugh, Chung-Ming Ko

Abstract: We investigate the mass-velocity dispersion relation (MVDR) in 29 galaxy clusters in the HIghest X-ray FLUx Galaxy Cluster Sample (HIFLUGCS). We measure the spatially resolved, line-of-sight velocity dispersion profiles of these clusters, which we find to be mostly flat at large radii, reminiscent of the rotation curves of galaxies. We discover a tight empirical relation between the baryonic mass… ▽ More We investigate the mass-velocity dispersion relation (MVDR) in 29 galaxy clusters in the HIghest X-ray FLUx Galaxy Cluster Sample (HIFLUGCS). We measure the spatially resolved, line-of-sight velocity dispersion profiles of these clusters, which we find to be mostly flat at large radii, reminiscent of the rotation curves of galaxies. We discover a tight empirical relation between the baryonic mass $M_\mathrm{bar}$ and the flat velocity dispersion $σ$ of the member galaxies, i.e. MVDR, $\log(M_\mathrm{bar}/M_\odot)=4.1^{+0.4}_{-0.4}\,\log(σ/\mathrm{km}\,\mathrm{s}^{-1})+1.6^{+1.0}_{-1.3}$, with the lognormal intrinsic scatter of $12^{+3}_{-3}\%$. The residuals of the MVDR are uncorrelated with other cluster properties like temperature, cluster radius, baryonic mass surface density, and redshift. These characteristics are reminiscent of the MVDR for individual galaxies, albeit about ten times larger characteristic acceleration scale. The cluster baryon fraction falls short of the cosmic value, exposing a problem: the discrepancy increases systematically for clusters of lower mass and lower baryonic acceleration. △ Less

Submitted 9 February, 2021; v1 submitted 2 October, 2020; originally announced October 2020.

Comments: Version accepted by ApJ; 10 pages, 5 figures, 2 Tables

Testing the Strong Equivalence Principle: Detection of the External Field Effect in Rotationally Supported Galaxies

Authors: Kyu-Hyun Chae, Federico Lelli, Harry Desmond, Stacy S. McGaugh, Pengfei Li, James M. Schombert

Abstract: The strong equivalence principle (SEP) distinguishes General Relativity from other viable theories of gravity. The SEP demands that the internal dynamics of a self-gravitating system under free-fall in an external gravitational field should not depend on the external field strength. We test the SEP by investigating the external field effect (EFE) in Milgromian dynamics (MOND), proposed as an alter… ▽ More The strong equivalence principle (SEP) distinguishes General Relativity from other viable theories of gravity. The SEP demands that the internal dynamics of a self-gravitating system under free-fall in an external gravitational field should not depend on the external field strength. We test the SEP by investigating the external field effect (EFE) in Milgromian dynamics (MOND), proposed as an alternative to dark matter in interpreting galactic kinematics. We report a detection of this EFE using galaxies from the Spitzer Photometry and Accurate Rotation Curves (SPARC) sample together with estimates of the large-scale external gravitational field from an all-sky galaxy catalog. Our detection is threefold: (1) the EFE is individually detected at $8σ$ to $11σ$ in "golden" galaxies subjected to exceptionally strong external fields, while it is not detected in exceptionally isolated galaxies, (2) the EFE is statistically detected at more than $4σ$ from a blind test of 153 SPARC rotating galaxies, giving a mean value of the external field consistent with an independent estimate from the galaxies' environments, and (3) we detect a systematic downward trend in the weak gravity part of the radial acceleration relation at the right acceleration predicted by the EFE of the MOND modified gravity. Tidal effects from neighboring galaxies in the $Λ$CDM context are not strong enough to explain these phenomena. They are not predicted by existing $Λ$CDM models of galaxy formation and evolution, adding a new small-scale challenge to the $Λ$CDM paradigm. Our results point to a breakdown of the SEP, supporting modified gravity theories beyond General Relativity. △ Less

Submitted 4 March, 2021; v1 submitted 24 September, 2020; originally announced September 2020.

Comments: ApJ (published), 14 figures, 2 tables (several environmental gravitational fields corrected, Table 2 indexing error fixed, and Figs 5 & 6 revised from an Erratum accepted for publication)

Self-consistent Color-Mass-to-Light-Ratio relations for Low Surface Brightness Galaxies

Authors: Wei Du, Stacy S. McGaugh

Abstract: The color - stellar mass-to-light ratio relation (CMLR) is a widely accepted tool to estimate the stellar mass (M*) of a galaxy. However, an individual CMLR tends to give distinct M* for a same galaxy when it is applied in different bands. Examining five representative CMLRs from literature, we find that the difference in M* predicted in different bands from optical to near-infrared by a CMLR is 0… ▽ More The color - stellar mass-to-light ratio relation (CMLR) is a widely accepted tool to estimate the stellar mass (M*) of a galaxy. However, an individual CMLR tends to give distinct M* for a same galaxy when it is applied in different bands. Examining five representative CMLRs from literature, we find that the difference in M* predicted in different bands from optical to near-infrared by a CMLR is 0.1-0.3 dex. Therefore, based on a sample of low surface brightness galaxies (LSBG) that covers a wide range of color and luminosity, we re-calibrated each original CMLR in r, i, z, J, H, and K bands to give internally self-consistent M* for a same galaxy. The g-r is the primary color indicator in the re-calibrated relations which show little dependenceon red (r - z) or near-infrared (J - K) colors.Additionally, the external discrepancies in the originally predicted stellar mass-to-light ratio (M*/L) by the five independent CMLRs have been greatly reduced after re-calibration, especially in near-infrared bands, implying that the near-infrared luminosities are more robust to predict M*/L. For each CMLR, the re-calibrated relations provided in this work could produce internally self-consistent M* from divergent photometric bands, and are extensions of the re-calibrations from Johnson-Cousin filter system by the pioneering work of McGaugh & Schombert (2014) to SDSS filter system. △ Less

Submitted 16 July, 2020; originally announced July 2020.

Comments: 19 pages,11 figures, 7 tables, accepted for publication in The Astronomical Journal

Using The Baryonic Tully-Fisher Relation to Measure $H_o$

Authors: James Schombert, Stacy McGaugh, Federico Lelli

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 tip magnitude of the red giant branch, we calibrate the bTFR on… ▽ More 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 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}$. △ Less

Submitted 1 July, 2020; v1 submitted 15 June, 2020; originally announced June 2020.

Comments: 21 pages, 3 figures, accepted by AJ

Can Bars Erode Cuspy Halos ?

Authors: Sandeep Kumar Kataria, Mousumi Das, Stacy McGaugh

Abstract: One of the major and widely known small scale problem with the Lambda CDM model of cosmology is the core-cusp problem. In this study we investigate whether this problem can be resolved using bar instabilities. We see that all the initial bars are thin (b/a $<$ 0.3) in our simulations and the bar becomes thick (b/a $>$ 0.3$) faster in the high resolution simulations. By increasing the resolution, w… ▽ More One of the major and widely known small scale problem with the Lambda CDM model of cosmology is the core-cusp problem. In this study we investigate whether this problem can be resolved using bar instabilities. We see that all the initial bars are thin (b/a $<$ 0.3) in our simulations and the bar becomes thick (b/a $>$ 0.3$) faster in the high resolution simulations. By increasing the resolution, we mean a larger number of disk particles. The thicker bars in the high resolution simulations transfer less angular momentum to the halo. Hence, we find that in the high resolution simulations it takes around 7 Gyr for the bar to remove inner dark matter cusp which is too long to be meaningful in galaxy evolution timescales. Physically, the reason is that as the resolution increases, the bar buckles faster and becomes thicker much earlier on. △ Less

Submitted 8 May, 2020; originally announced May 2020.

Comments: 2 pages, 2 figures, Galactic Dynamics in the Era of Large Surveys Proceedings IAU Symposium No 353, 2019

Predictions and Outcomes for the Dynamics of Rotating Galaxies

Authors: Stacy McGaugh

Abstract: A review is given of a priori predictions made for the dynamics of rotating galaxies. One theory - MOND - has had many predictions corroborated by subsequent observations. While it is sometimes possible to offer post hoc explanations for these observations in terms of dark matter, it is seldom possible to use dark matter to predict the same phenomena. A review is given of a priori predictions made for the dynamics of rotating galaxies. One theory - MOND - has had many predictions corroborated by subsequent observations. While it is sometimes possible to offer post hoc explanations for these observations in terms of dark matter, it is seldom possible to use dark matter to predict the same phenomena. △ Less

Submitted 29 April, 2020; originally announced April 2020.

Comments: 36 pages (10 are references), 9 figures. Invited review for the Galaxies special Issue "Debate on the Physics of Galactic Rotation and the Existence of Dark Matter." Provides test cases for the importance of prior predictions in the application of the scientific method

Journal ref: Galaxies 2020, 8(2), 35

A comprehensive catalog of dark matter halo models for SPARC galaxies

Authors: Pengfei Li, Federico Lelli, Stacy McGaugh, James Schombert

Abstract: We present rotation curve fits to 175 late-type galaxies from the Spitzer Photometry & Accurate Rotation Curves (SPARC) database using seven dark matter (DM) halo profiles: pseudo-isothermal (pISO), Burkert, Navarro-Frenk-White (NFW), Einasto, Di Cintio (2014, DC14), coreNFW, and a new semi-empirical profile named Lucky13. We marginalize over stellar mass-to-light ratio, galaxy distance, disk incl… ▽ More We present rotation curve fits to 175 late-type galaxies from the Spitzer Photometry & Accurate Rotation Curves (SPARC) database using seven dark matter (DM) halo profiles: pseudo-isothermal (pISO), Burkert, Navarro-Frenk-White (NFW), Einasto, Di Cintio (2014, DC14), coreNFW, and a new semi-empirical profile named Lucky13. We marginalize over stellar mass-to-light ratio, galaxy distance, disk inclination, halo concentration and halo mass (and an additional shape parameter for Einasto) using a Markov Chain Monte Carlo method. We find that cored halo models such as the DC14 and Burkert profiles generally provide better fits to rotation curves than the cuspy NFW profile. The stellar mass-halo mass relation from abundance matching is recovered by all halo profiles once imposed as a Bayesian prior, whereas the halo mass-concentration relation is not reproduced in detail by any halo model. We provide an extensive set of figures as well as best-fit parameters in machine-readable tables to facilitate model comparison and the exploration of DM halo properties. △ Less

Submitted 30 January, 2020; v1 submitted 28 January, 2020; originally announced January 2020.

Comments: 11 pages, 6 figures, 1 table. Accepted for publication in ApJS

Tracing the Dynamical Mass in Galaxy Disks Using HI Velocity Dispersion and its Implications for the Dark Matter Distribution in Galaxies

Authors: Mousumi Das, Stacy McGaugh, Roger Ianjamasimanana, James Schombert, K. S. Dwarakanath

Abstract: We present a method to derive the dynamical mass of face-on galaxy disks using their neutral hydrogen (HI) velocity dispersion. We have applied the method to nearby, gas rich galaxies that have extended HI gas disks and have low inclinations. The galaxy sample includes 4 large disk galaxies; NGC628, NGC6496, NGC3184, NGC4214 and 3 dwarf galaxies DDO46, DDO63 and DDO187. We have used archival data… ▽ More We present a method to derive the dynamical mass of face-on galaxy disks using their neutral hydrogen (HI) velocity dispersion. We have applied the method to nearby, gas rich galaxies that have extended HI gas disks and have low inclinations. The galaxy sample includes 4 large disk galaxies; NGC628, NGC6496, NGC3184, NGC4214 and 3 dwarf galaxies DDO46, DDO63 and DDO187. We have used archival data from the THINGS and LITTLE THINGS surveys to derive the HI gas distributions and SPITZER mid-infrared images to determine the stellar disk mass distributions. We examine the disk dynamical and baryonic mass ratios in the extreme outer disks where there is HI gas but no visible stellar disk. We find that for the large galaxies the disk dynamical and Hi gas mass surface densities are comparable in the outer disks. But in the smaller dwarf galaxies, for which the total HI gas mass dominates the baryonic mass i.e. M(HI)>M(stars), the disk dynamical mass is much larger than the baryonic mass. For these galaxies there must either be a very low luminosity stellar disk which provides the vertical support for the HI gas disk or there is halo dark matter associated with their disks, which is possible if the halo has an oblate shape so that the inner part of the dark matter halo is concentrated around the disk. Our results are important for explaining the equilibrium of HI disks in the absence of stellar disks, and is especially important for gas rich, dwarf galaxies that appear to have significant dark matter masses associated with their disks. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: 13 pages, 4 figures and 2 tables. To appear in ApJ

The halo mass function of late-type galaxies from HI kinematics

Authors: Pengfei Li, Federico Lelli, Stacy McGaugh, Marcel S. Pawlowski, Martin A. Zwaan, James Schombert

Abstract: We present an empirical method to measure the halo mass function (HMF) of galaxies. We determine the relation between the \hi\ line-width from single-dish observations and the dark matter halo mass ($M_{200}$) inferred from rotation curve fits in the SPARC database, then we apply this relation to galaxies from the \hi\ Parkes All Sky Survey (HIPASS) to derive the HMF. This empirical HMF is well fi… ▽ More We present an empirical method to measure the halo mass function (HMF) of galaxies. We determine the relation between the \hi\ line-width from single-dish observations and the dark matter halo mass ($M_{200}$) inferred from rotation curve fits in the SPARC database, then we apply this relation to galaxies from the \hi\ Parkes All Sky Survey (HIPASS) to derive the HMF. This empirical HMF is well fit by a Schecther function, and matches that expected in $Λ$CDM over the range $10^{10.5} < M_{200} < 10^{12}\;\mathrm{M}_{\odot}$. More massive halos must be poor in neutral gas to maintain consistency with the power law predicted by $Λ$CDM. We detect no discrepancy at low masses. The lowest halo mass probed by HIPASS, however, is just greater than the mass scale where the Local Group missing satellite problem sets in. The integrated mass density associated with the dark matter halos of \hi-detected galaxies sums to $Ω_{\rm m,gal} \approx 0.03$ over the probed mass range. △ Less

Submitted 1 November, 2019; originally announced November 2019.

Comments: 7 pages, 3 figures, 1 table. Accepted for publication in ApJ Letters

The Imprint of Spiral Arms on the Galactic Rotation Curve

Authors: Stacy McGaugh

Abstract: We discuss a model for the Milky Way obtained by fitting the observed terminal velocities with the radial acceleration relation. The resulting stellar surface density profile departs from a smooth exponential disk, having bumps and wiggles that correspond to massive spiral arms. These features are used to estimate the term for the logarithmic density gradient in the Jeans equation, which turn out… ▽ More We discuss a model for the Milky Way obtained by fitting the observed terminal velocities with the radial acceleration relation. The resulting stellar surface density profile departs from a smooth exponential disk, having bumps and wiggles that correspond to massive spiral arms. These features are used to estimate the term for the logarithmic density gradient in the Jeans equation, which turn out to have exactly the right location and amplitude to reconcile the apparent discrepancy between the stellar rotation curve and that of the interstellar gas. This model also predicts a gradually declining rotation curve outside the solar circle with slope $-1.7\;\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{kpc}^{-1}$, as subsequently observed. △ Less

Submitted 24 September, 2019; originally announced September 2019.

Comments: Accepted for publication in the Astrophysical Journal. 16 pages, 4 figures, 4 tables

Journal ref: 2019, ApJ, 885, 87

Dynamical Regularities in Galaxies

Authors: Stacy McGaugh, Federico Lelli, Pengfei Li, Jim Schombert

Abstract: Galaxies are observed to obey a strict set of dynamical scaling relations. We review these relations for rotationally supported disk galaxies spanning many decades in mass, surface brightness, and gas content. The behavior of these widely varied systems can be summarized with a handful of empirical laws connected by a common acceleration scale. Galaxies are observed to obey a strict set of dynamical scaling relations. We review these relations for rotationally supported disk galaxies spanning many decades in mass, surface brightness, and gas content. The behavior of these widely varied systems can be summarized with a handful of empirical laws connected by a common acceleration scale. △ Less

Submitted 4 September, 2019; originally announced September 2019.

Comments: 8 pages, 6 figures; invited review for IAU Symposium 353: Galactic Dynamics in the Era of Large Surveys

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

Authors: Federico Lelli, Stacy S. McGaugh, James M. Schombert, Harry Desmond, Harley Katz

Abstract: 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… ▽ More 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. △ Less

Submitted 17 January, 2019; originally announced January 2019.

Comments: 12 pages, 6 figures, accepted for publication in MNRAS

Journal ref: MNRAS 2019 484 (3): 3267-3278

The Mass-to-light Ratios and the Star Formation Histories of Disk Galaxies

Authors: J. Schombert, S. McGaugh, F. Lelli

Abstract: We combine new data from the main sequence (M_* versus SFR) of star-forming galaxies and galaxy colors (from GALEX to Spitzer) with a flexible stellar population scheme to deduce the mass-to-light ratio (Υ_*) of star-forming galaxies from the SPARC and S^4G samples. We find that the main sequence for galaxies, particular the low-mass end, combined with the locus of galaxy colors, constrains the po… ▽ More We combine new data from the main sequence (M_* versus SFR) of star-forming galaxies and galaxy colors (from GALEX to Spitzer) with a flexible stellar population scheme to deduce the mass-to-light ratio (Υ_*) of star-forming galaxies from the SPARC and S^4G samples. We find that the main sequence for galaxies, particular the low-mass end, combined with the locus of galaxy colors, constrains the possible star formation histories of disk and dwarf galaxies to a similar shape found by Speagle et al. (2014). Combining the deduced star formation history with stellar population models in the literature produces reliable Υ_* values as a function of galaxy color with an uncertainty of only 0.05 dex. We provide prescriptions to deduce Υ_* for optical and near-IR bandpasses, with near-IR bandpasses having the least uncertainty (Υ_* from 0.40 to 0.55). We also provide the community with a webtool, with flexible stellar population parameters, to generate their own Υ_* values over the wavelength range for most galaxy surveys. △ Less

Submitted 26 November, 2018; originally announced November 2018.

Comments: 17 pages, 12 figures, accepted by MNRAS

A constant characteristic volume density of dark matter haloes from SPARC rotation curve fits

Authors: Pengfei Li, Federico Lelli, Stacy S. McGaugh, Nathaniel Starkman, James M. Schombert

Abstract: We study the scaling relations between dark matter (DM) haloes and galaxy discs using 175 galaxies from the SPARC database. We explore two cosmologically motivated DM halo profiles: the Einasto profile from DM-only simulations and the DC14 profile from hydrodynamic simulations. We fit the observed rotation curves using a Markov Chain Monte Carlo method and break the disc-halo degeneracy using near… ▽ More We study the scaling relations between dark matter (DM) haloes and galaxy discs using 175 galaxies from the SPARC database. We explore two cosmologically motivated DM halo profiles: the Einasto profile from DM-only simulations and the DC14 profile from hydrodynamic simulations. We fit the observed rotation curves using a Markov Chain Monte Carlo method and break the disc-halo degeneracy using near-infrared photometry and $Λ$CDM-motivated priors. We find that the characteristic volume density $ρ_{\rm s}$ of DM haloes is nearly constant over $\sim$5 decades in galaxy luminosity. The scale radius $r_s$ and the characteristic surface density $ρ_s\cdot r_s$, instead, correlate with galaxy luminosity. These scaling relations provide an empirical benchmark to cosmological simulations of galaxy formation. △ Less

Submitted 1 November, 2018; originally announced November 2018.

Comments: 19 pages, 11 figures, 3 tables. Accepted for publication in MNRAS

The Tight Empirical Relation between Dark Matter Halo Mass and Flat Rotation Velocity for Late-Type Galaxies

Authors: Harley Katz, Harry Desmond, Stacy McGaugh, Federico Lelli

Abstract: We present a new empirical relation between galaxy dark matter halo mass (${\rm M_{halo}}$) and the velocity along the flat portion of the rotation curve (${\rm V_{flat}}$), derived from 120 late-type galaxies from the SPARC database. The orthogonal scatter in this relation is comparable to the observed scatter in the baryonic Tully-Fisher relation (BTFR), indicating a tight coupling between total… ▽ More We present a new empirical relation between galaxy dark matter halo mass (${\rm M_{halo}}$) and the velocity along the flat portion of the rotation curve (${\rm V_{flat}}$), derived from 120 late-type galaxies from the SPARC database. The orthogonal scatter in this relation is comparable to the observed scatter in the baryonic Tully-Fisher relation (BTFR), indicating a tight coupling between total halo mass and galaxy kinematics at $r\ll R_{\rm vir}$. The small vertical scatter in the relation makes it an extremely competitive estimator of total halo mass. We demonstrate that this conclusion holds true for different priors on $M_*/L_{[3.6μ]}$ that give a tight BTFR, but requires that the halo density profile follows DC14 rather than NFW. We provide additional relations between ${\rm M_{halo}}$ and other velocity definitions at smaller galactic radii (i.e. ${\rm V_{2.2}}$, ${\rm V_{eff}}$, and ${\rm V_{max}}$) which can be useful for estimating halo masses from kinematic surveys, providing an alternative to abundance matching. Furthermore, we constrain the dark matter analog of the Radial Acceleration Relation and also find its scatter to be small, demonstrating the fine balance between baryons and dark matter in their contribution to galaxy kinematics. △ Less

Submitted 29 October, 2018; originally announced October 2018.

Comments: 6 pages, 4 figures, Accepted to MNRAS Letters

A Precise Milky Way Rotation Curve Model for an Accurate Galactocentric Distance

Authors: Stacy McGaugh

Abstract: I provide a model rotation curve for the Milky Way that matches the details of the terminal velocity curve normalized to the Galactocentric distance $R_0 = 8.122$ kpc obtained by the GRAVITY collaboration and the corresponding circular speed of the LSR $Θ_0 = 233.3$ km/s. The model provides a numerical representation of the azimuthally averaged radial run of the gravitational potential of each mas… ▽ More I provide a model rotation curve for the Milky Way that matches the details of the terminal velocity curve normalized to the Galactocentric distance $R_0 = 8.122$ kpc obtained by the GRAVITY collaboration and the corresponding circular speed of the LSR $Θ_0 = 233.3$ km/s. The model provides a numerical representation of the azimuthally averaged radial run of the gravitational potential of each mass component of the Galaxy (bulge-bar, stellar disk, gas disk, and dark matter) as represented by the rotation curve of each. It provides precise estimates of quantities like the stellar mass of the Galaxy ($6.16 \pm 0.31 \times 10^{10}\;\mathrm{M}_{\odot}$) and the local density of dark matter ($ρ_{DM}(R_0) = 6.76^{+0.08}_{-0.14} \times 10^{-3}\; \mathrm{M}_{\odot}\,\mathrm{pc}^{-3} = 0.257^{+0.003}_{-0.005}\; \mathrm{GeV}\,\mathrm{cm}^{-3}$). The dark matter density implied by the radial force is less than that found in many studies of the vertical force, perhaps indicating that the usual assumption of a spherical dark matter halo is no longer adequate. △ Less

Submitted 28 August, 2018; originally announced August 2018.

Comments: To be published in RNAAS. 2 pages, 1 figure. Numerical data available at http://astroweb.case.edu/ssm/models/

Predictions for the sky-averaged depth of the 21cm absorption signal at high redshift in cosmologies with and without non-baryonic cold dark matter

Authors: Stacy McGaugh

Abstract: We consider the 21cm absorption signal expected at high redshift in cosmologies with and without non-baryonic cold dark matter. The expansion of the early universe decelerates strongly with dark matter, but approximately coasts without it. This results in a different path length across the epochs when absorption is expected, with the consequence that the absorption is predicted to be a factor of… ▽ More We consider the 21cm absorption signal expected at high redshift in cosmologies with and without non-baryonic cold dark matter. The expansion of the early universe decelerates strongly with dark matter, but approximately coasts without it. This results in a different path length across the epochs when absorption is expected, with the consequence that the absorption is predicted to be a factor of $\sim 2$ greater without dark matter than with it. Observation of such a signal would motivate consideration of extended theories of gravity in lieu of dark matter. △ Less

Submitted 7 August, 2018; originally announced August 2018.

Comments: Accepted for publication in Phys. Rev. Letters. 9 pages (4 of which are bibliography), 1 figure. Model data available at http://astroweb.case.edu/ssm/models/

Stellar Feedback and the Energy Budget of Late-Type Galaxies: Missing Baryons and Core Creation

Authors: Harley Katz, Harry Desmond, Federico Lelli, Stacy McGaugh, Arianna Di Cintio, Chris Brook, James Schombert

Abstract: In a $Λ$CDM cosmology, galaxy formation is a globally inefficient process: it is often the case that far fewer baryons are observed in galaxy disks than expected from the cosmic baryon fraction. The location of these "missing baryons" is unclear. By fitting halo profiles to the rotation curves of galaxies in the SPARC data set, we measure the "missing baryon" mass for individual late-type systems.… ▽ More In a $Λ$CDM cosmology, galaxy formation is a globally inefficient process: it is often the case that far fewer baryons are observed in galaxy disks than expected from the cosmic baryon fraction. The location of these "missing baryons" is unclear. By fitting halo profiles to the rotation curves of galaxies in the SPARC data set, we measure the "missing baryon" mass for individual late-type systems. Assuming that haloes initially accrete the cosmological baryon fraction, we show that the maximum energy available from supernovae is typically not enough to completely eject these "missing baryons" from a halo, but it is often sufficient to heat them to the virial temperature. The energy available from supernovae has the same scaling with galaxy mass as the energy needed to heat or eject the "missing baryons", indicating that the coupling efficiency of the feedback to the ISM may be constant with galaxy virial mass. We further find that the energy available from supernova feedback is always enough to convert a primordial cusp into a core and has magnitude consistent with what is required to heat the "missing baryons" to the virial temperature. Taking a census of the baryon content of galaxies with ${\rm 10^9<M_{vir}/M_{\odot}<10^{12}}$ reveals that $\sim86\%$ of baryons are likely to be in a hot phase surrounding the galaxies and possibly observable in the X-ray, $\sim7\%$ are in the form of cold gas, and $\sim7\%$ are in stars. △ Less

Submitted 2 August, 2018; originally announced August 2018.

Comments: 16 pages, 8 figures, MNRAS accepted

Uncorrelated velocity and size residuals across galaxy rotation curves

Authors: Harry Desmond, Harley Katz, Federico Lelli, Stacy McGaugh

Abstract: 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 mo… ▽ More 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 $ΔV_r$ and mass-size residuals $ΔR$ as a function of the radius $r$ at which the velocity is measured, and thus present the $ΔV_r-Δ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 $ΔV_r-Δ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. △ Less

Submitted 11 January, 2019; v1 submitted 1 August, 2018; originally announced August 2018.

Comments: 7 pages, 4 figures; revised to match MNRAS published version

Journal ref: MNRAS 2019 484 (1): 239-244