We estimate the mass of the intermediate-mass black hole at the heart of the dwarf elliptical galaxy NGC 404 using Atacama Large Millimeter/submillimeter Array (ALMA) observations of the molecular interstellar medium at an unprecedented linear resolution of ~0.5 pc, in combination with existing stellar kinematic information. These ALMA observations reveal a central disc/torus of molecular gas clearly rotating around the black hole. This disc is surrounded by a morphologically and kinematically complex flocculent distribution of molecular clouds, that we resolve in detail. Continuum emission is detected from the central parts of NGC 404, likely arising from the Rayleigh-Jeans tail of emission from dust around the nucleus, and potentially from dusty massive star-forming clumps at discrete locations in the disc. Several dynamical measurements of the black hole mass in this system have been made in the past, but they do not agree. We show here that both the observed molecular gas and stellar kinematics independently require a ~5x10$^5$ Msun black hole once we include the contribution of the molecular gas to the potential. Our best estimate comes from the high-resolution molecular gas kinematics, suggesting the black hole mass of this system is 5.5$^{+4.1}_{-3.8}\times$10$^5$ Msun (at the 99% confidence level), in good agreement with our revised stellar kinematic measurement and broadly consistent with extrapolations from the black hole mass - velocity dispersion and black hole mass - bulge mass relations. This highlights the need to accurately determine the mass and distribution of each dynamically important component around intermediate-mass black holes when attempting to estimate their masses.

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}$.

Rotation curves have traditionally been difficult to trace for early-type galaxies (ETGs) because they often lack a high-density disk of cold gas as in late-type galaxies (LTGs). We derive rotation curves for three lenticular galaxies from the ATLAS3D survey, combining CO data in the inner parts with deep HI data in the outer regions, extending out to 10-20 effective radii. We also use Spitzer photometry at 3.6 um to decompose the rotation curves into the contributions of baryons and dark matter (DM). We find that (1) the rotation-curve shapes of these ETGs are similar to those of LTGs of similar mass and surface brightness; (2) the dynamically$-$inferred stellar mass-to-light ratios are small for quiescent ETGs but similar to those of star-forming LTGs; (3) the DM halos follow the same scaling relations with galaxy luminosity as those of LTGs; (4) one galaxy (NGC 3626) is poorly fitted by cuspy DM profiles, suggesting that DM cores may exist in high-mass galaxies too. Our results indicate that these lenticular galaxies have recently transitioned from LTGs to ETGs without altering their DM halo structure (e.g., via a major merger) and could be faded spirals. We also confirm that ETGs follow the same radial acceleration relation as LTGs, reinforcing the notion that this is a universal law for all galaxy types.

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.

The phenomenology of modified Newtonian dynamics (MOND) on galaxy scales may point to more fundamental theories of either modified gravity (MG) or modified inertia (MI). In this paper, we test the applicability of the global deep-MOND parameter $Q$ which is predicted to vary at the $10\%$ level between MG and MI theories. Using mock-observed analytical models of disk galaxies, we investigate several observational uncertainties, establish a set of quality requirements for actual galaxies, and derive systematic corrections in the determination of $Q$. Implementing our quality requirements to the SPARC database yields $15$ galaxies, which are close enough to the deep-MOND regime as well as having rotation curves that are sufficiently extended and sampled. For these galaxies, the average and median values of $Q$ seem to favor MG theories, albeit both MG and MI predictions are in agreement with the data within $1.5\sigma$. Improved precision in the determination of $Q$ can be obtained by measuring extended and finely-sampled rotation curves for a significant sample of extremely low-surface-brightness galaxies.

Studies of nucleated dwarf galaxies can constrain the scenarios for the formation and evolution of nuclear star clusters (NSC) in low-mass galaxies and give us insights on the origin of ultra compact dwarf galaxies (UCDs). We report the discovery of a NSC in the dwarf galaxy KKs58 and investigate its properties together with those of another NSC in KK197. Both NSCs are hosted by dwarf elliptical galaxies of the Centaurus group. Combining ESO VLT MUSE data with photometry from VLT FORS2, CTIO Blanco DECam, and HST ACS, as well as high-resolution spectroscopy from VLT UVES, we analyse the photometric, kinematic and stellar population properties of the NSCs and their host galaxies. We confirm membership of the NSCs based on their radial velocities and location close to the galaxy centres. We also confirm the membership of two globular clusters (GCs) and detect oblate rotation in the main body of KK197. Based on high signal-to-noise spectra taken with MUSE of the NSCs of both KKs58 and KK197 we measure low metallicities, [Fe/H] = $-1.75 \pm 0.06$ dex and [Fe/H] = $-1.84 \pm 0.05$ dex, and stellar masses of $7.3 \times 10^5 M_\odot$ and $1.0 \times 10^6 M_\odot$, respectively. Both NSCs are more metal-poor than their hosts that have metallicities of $-1.35 \pm 0.23$ dex (KKs58) and $-0.84 \pm 0.12$ dex (KK197). This can be interpreted as NSC formation via the in-spiral of GCs. The masses, sizes and metallicities of the two NSCs place them among other NSCs, but also among the known UCDs of the Centaurus group. This indicates that NSCs might constitute the progenitors of a part of the low-mass UCDs, although their properties are almost indistinguishable from typical GCs.

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 $\Lambda$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 $\Lambda$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 $\Omega_{\rm m,gal} \approx 0.03$ over the probed mass range.

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.

Dwarf galaxy satellite systems are essential probes to test models of structure formation, making it necessary to establish a census of dwarf galaxies outside of our own Local Group. We present deep FORS2 VI band images from the ESO Very Large Telescope for 15 dwarf galaxy candidates in the Centaurus group of galaxies. We confirm nine dwarfs to be members of CenA by measuring their distances using a Bayesian approach to determine the tip of the red giant branch luminosity. We have also fitted theoretical isochrones to measure their mean metallicities. The properties of the new dwarfs are similar to those in the Local Group in terms of their sizes, luminosities, and mean metallicities. Within our photometric precision there is no evidence of a metallicity spread, but we do observe possible extended star formation in several galaxies, as evidenced by a population of asymptotic giant branch stars brighter than the red giant branch tip. The new dwarfs do not show any signs of tidal disruption. Together with the recently reported dwarf galaxies by the complementary PISCeS survey, we study the luminosity function and 3D structure of the group. By comparing the observed luminosity function to the high-resolution cosmological simulation IllustrisTNG, we find agreement within a 90% confidence interval. However, CenA seems to be missing its brightest satellites and has an overabundance of the faintest dwarfs in comparison to its simulated analogs. In terms of the overall 3D distribution of the observed satellites, we find that the whole structure is flattened along the line-of-sight, with an root-mean-square (rms) height of 130 kpc and a rms semi major axis length of 330 kpc. Future distance measurements of the remaining dwarf galaxy candidates are needed to complete the census of dwarf galaxies in the Centaurus group.

The formation of globular clusters remains an open debate. Dwarf starburst galaxies are efficient at forming young massive clusters with similar masses as globular clusters and may hold the key to understanding their formation. We study star cluster formation in a tidal debris - including the vicinity of three tidal dwarf galaxies - in a massive gas dominated collisional ring around NGC~5291. These dwarfs have physical parameters which differ significantly from local starbursting dwarfs. They are gas-rich, highly turbulent, have a gas metallicity already enriched up to half-solar, and are expected to be free of dark matter. The aim is to study massive star cluster formation in this as yet unexplored type of environment. We use imaging from the Hubble Space Telescope using broadband filters covering the wavelength range from the near-ultraviolet to the near-infrared. We determine the masses and ages of the cluster candidates by using the spectral energy distribution-fitting code CIGALE, carefully considering age-extinction degeneracy effects on the estimation of the physical parameters. We find that the tidal dwarf galaxies in the ring of NGC 5291 are forming star clusters with an average efficiency of $\sim40\%$, comparable to blue compact dwarf galaxies. We also find massive star clusters for which the photometry suggests that they were formed at the very birth of the tidal dwarf galaxies and have survived for several hundred million years. Therefore our study shows that extended tidal dwarf galaxies and compact clusters may be formed simultaneously. In the specific case observed here, the young star clusters are not massive enough to survive for a Hubble time. However one may speculate that similar objects at higher redshift, with higher star formation rate, might form some of the long lived globular clusters.

The so-called ultra-diffuse galaxy NGC~1052-DF2 was announced to be a galaxy lacking dark matter based on a spectroscopic study of its constituent globular clusters. Here we present the first spectroscopic analysis of the stellar body of this galaxy using the MUSE integral-field spectrograph at the (ESO) Very Large Telescope. The MUSE datacube simultaneously provides DF2's stellar velocity field and systemic velocities for seven globular clusters (GCs). We further discovered three planetary nebulae (PNe) that are likely part of this galaxy. While five of the clusters had velocities measured in the literature, we were able to confirm the membership of two more candidates through precise radial velocity measurements, which increases the measured specific frequency of GCs in DF2. The mean velocity of the diffuse stellar body, 1792.9$^{-1.8}_{+1.4}$~\kms, is consistent with the mean globular cluster velocity. We detect a weak but significant velocity gradient within the stellar body, with a kinematic axis close to the photometric major-axis, making it a prolate-like rotator. We estimate a velocity dispersion from the clusters and PNe of $\sigma_{\mathrm{int}}=10.6^{+3.9}_{-2.3}$~\kms. The velocity dispersion $\sigma_{\rm{DF2}\star}$(\re) for the stellar body within one effective radius is $10.8^{-4.0}_{+3.2}$~\kms. Considering various sources of systemic uncertainties this central value varies between 5 and 13~\kms, and we conservatively report a 95\% confidence upper limit to the dispersion within one \re\ of 21~\kms. We provide updated mass estimates based on these dispersions corresponding to the different distances to NGC~1052-DF2 that have been reported in the recent literature.

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.

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

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 (\Upsilon_*) 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 \Upsilon_* values as a function of galaxy color with an uncertainty of only 0.05 dex. We provide prescriptions to deduce \Upsilon_* for optical and near-IR bandpasses, with near-IR bandpasses having the least uncertainty (\Upsilon_* from 0.40 to 0.55). We also provide the community with a webtool, with flexible stellar population parameters, to generate their own \Upsilon_* values over the wavelength range for most galaxy surveys.

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 $\Lambda$CDM-motivated priors. We find that the characteristic volume density $\rho_{\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 $\rho_s\cdot r_s$, instead, correlate with galaxy luminosity. These scaling relations provide an empirical benchmark to cosmological simulations of galaxy formation.

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\mu]}$ 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.

In a $\Lambda$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.

Maximum disc decompositions of rotation curves place a dynamical upper limit to the mass attributable to stars in galaxies. The precise definition of this term, however, can be vague and varies in usage. We develop an algorithm to robustly quantify maximum-disc mass models and apply it to 153 galaxies from the SPARC database. Our automatic procedure recovers classic results from manual decompositions. High-mass, high-surface-brightness galaxies have mean maximum-disc mass-to-light ratios of $\sim 0.7 \;{\mathrm{M}_\odot}/{\mathrm{L}_\odot}$ in the Spitzer 3.6 $\mu$m band, which are close to the expectations from stellar population models, suggesting that these galaxies are nearly maximal. Low-mass, low-surface-brightness galaxies have very high maximum-disc mass-to-light ratios (up to 10 $\mathrm{M}_\odot/\mathrm{L}_\odot$), which are unphysical for standard stellar population models, confirming they are sub-maximal. The maximum-disc mass-to-light ratios are more closely correlated with surface brightness than luminosity. The mean ratio between baryonic and observed velocity at the peak of the baryonic contribution is $\mathrm{V}_{\mathrm{bar}}/{\mathrm{V}_\mathrm{p}} \approx 0.88$, but correlates with surface brightness, so it is unwise to use this mean value to define the maximum disc concept. Our algorithm requires no manual intervention and could be applied to large galaxy samples from future HI surveys with Apertif, Askap, and SKA.

We present a multiwavelength study of the AGN-host starburst galaxy PKS 0529-549 at z~2.6. We use (1) new ALMA observations of the dust continuum and of the [CI] 370 um line, tracing molecular gas, (2) SINFONI spectroscopy of the [OIII] 5007 Ang line, tracing ionized gas, and (3) ATCA radio continuum images, tracing synchrotron emission. Both [CI] and [OIII] show regular velocity gradients, but their systemic velocities and position angles differ by ~300 km/s and ~30 degrees, respectively. The [CI] is consistent with a rotating disc, aligned with the dust and stellar continuum, while the [OIII] likely traces an outflow, aligned with two AGN-driven radio lobes. We model the [CI] cube using 3D disc models, which give best-fit rotation velocities V~310 km/s and velocity dispersions sigma<30 km/s. Hence, the [CI] disc has V/sigma>10 and is not particularly turbulent, similar to local galaxy discs. The dynamical mass (~10^11 Msun) is comparable to the baryonic mass within the errors. This suggests that baryons dominate the inner galaxy dynamics, similar to massive galaxies at z=0. Remarkably, PKS 0529-549 lies on the local baryonic Tully-Fisher relation, indicating that at least some massive galaxies are already in place and kinematically relaxed at z~2.6. This work highlights the potential of the [CI] line to trace galaxy dynamics at high z, as well as the importance of multiwavelength data to interpret gas kinematics.

Galaxies follow a tight radial acceleration relation (RAR): the acceleration observed at every radius correlates with that expected from the distribution of baryons. We use the Markov Chain Monte Carlo method to fit the mean RAR to 175 individual galaxies in the SPARC database, marginalizing over stellar mass-to-light ratio ($\Upsilon_{\star}$), galaxy distance, and disk inclination. Acceptable fits with astrophysically reasonable parameters are found for the vast majority of galaxies. The residuals around these fits have an rms scatter of only 0.057 dex ($\sim$13$\%$). This is in agreement with the predictions of modified Newtonian dynamics (MOND). We further consider a generalized version of the RAR that, unlike MOND, permits galaxy-to-galaxy variation in the critical acceleration scale. The fits are not improved with this additional freedom: there is no credible indication of variation in the critical acceleration scale. The data are consistent with the action of a single effective force law. The apparent universality of the acceleration scale and the small residual scatter are key to understanding galaxies.

The Milky Way and Andromeda galaxy are each surrounded by a thin plane of satellite galaxies that may be corotating. Cosmological simulations predict that most satellite galaxy systems are close to isotropic with random motions, so those two well-studied systems are often interpreted as rare statistical outliers. We test this assumption using the kinematics of satellite galaxies around the Centaurus A galaxy. Our statistical analysis reveals evidence for corotation in a narrow plane: of the 16 Centaurus A's satellites with kinematic data, 14 follow a coherent velocity pattern aligned with the long axis of their spatial distribution. In standard cosmology simulations, < 0.5% of Centaurus A-like systems show such behavior. Corotating satellite systems may be common in the Universe, challenging small-scale structure formation in the prevailing cosmological paradigm.

We explore the star forming properties of late type, low surface brightness (LSB) galaxies. The star forming main sequence (SFR-$M_*$) of LSB dwarfs has a steep slope, indistinguishable from unity ($1.04 \pm 0.06$). They form a distinct sequence from more massive spirals, which exhibit a shallower slope. The break occurs around $M_* \approx 10^{10}\;M_{\odot}$, and can also be seen in the gas mass-stellar mass plane. The global Kennicutt-Schmidt law (SFR-$M_g$) has a slope of $1.47 \pm 0.11$ without the break seen in the main sequence. There is an ample supply of gas in LSB galaxies, which have gas depletion times well in excess of a Hubble time, and often tens of Hubble times. Only $\sim 3\%$ of this cold gas need be in the form of molecular gas to sustain the observed star formation. In analogy with the faint, long-lived stars of the lower stellar main sequence, it may be appropriate to consider the main sequence of star forming galaxies to be defined by thriving dwarfs (with $M_* < 10^{10}\;M_{\odot}$) while massive spirals (with $M_* > 10^{10}\;M_{\odot}$) are weary giants that constitute more of a turn-off population.

Verlinde (2016) has recently proposed that spacetime and gravity may emerge from an underlying microscopic theory. In a de Sitter spacetime, such emergent gravity (EG) contains an additional gravitational force due to dark energy, which may explain the mass discrepancies observed in galactic systems without the need of dark matter. For a point mass, EG is equivalent to Modified Newtonian Dynamics (MOND). We show that this equivalence does not hold for finite-size galaxies: there are significant differences between EG and MOND in the inner regions of galaxies. We confront theoretical predictions with the empirical Radial Acceleration Relation (RAR). We find that (i) EG is consistent with the observed RAR only if we substantially decrease the fiducial stellar mass-to-light ratios; the resulting values are in tension with other astronomical estimates; (ii) EG predicts that the residuals around the RAR should correlate with radius; such residual correlation is not observed.

We study the link between baryons and dark matter in 240 galaxies with spatially resolved kinematic data. Our sample spans 9 dex in stellar mass and includes all morphological types. We consider (i) 153 late-type galaxies (LTGs; spirals and irregulars) with gas rotation curves from the SPARC database; (ii) 25 early-type galaxies (ETGs; ellipticals and lenticulars) with stellar and HI data from ATLAS^3D or X-ray data from Chandra; and (iii) 62 dwarf spheroidals (dSphs) with individual-star spectroscopy. We find that LTGs, ETGs, and "classical" dSphs follow the same radial acceleration relation: the observed acceleration (gobs) correlates with that expected from the distribution of baryons (gbar) over 4 dex. The relation coincides with the 1:1 line (no dark matter) at high accelerations but systematically deviates from unity below a critical scale of ~10^-10 m/s^2. The observed scatter is remarkably small (<0.13 dex) and largely driven by observational uncertainties. The residuals do not correlate with any global or local galaxy property (baryonic mass, gas fraction, radius, etc.). The radial acceleration relation is tantamount to a Natural Law: when the baryonic contribution is measured, the rotation curve follows, and vice versa. Including ultrafaint dSphs, the relation may extend by another 2 dex and possibly flatten at gbar<10^-12 m/s^2, but these data are significantly more uncertain. The radial acceleration relation subsumes and generalizes several well-known dynamical properties of galaxies, like the Tully-Fisher and Faber-Jackson relations, the "baryon-halo" conspiracies, and Renzo's rule.

Cosmological $N$-body simulations predict dark matter (DM) haloes with steep central cusps (e.g. NFW, Navarro et al. 1996). This contradicts observations of gas kinematics in low-mass galaxies that imply the existence of shallow DM cores. Baryonic processes such as adiabatic contraction and gas outflows can, in principle, alter the initial DM density profile, yet their relative contributions to the halo transformation remain uncertain. Recent high resolution, cosmological hydrodynamic simulations (Di Cintio et al. 2014, DC14) predict that inner density profiles depend systematically on the ratio of stellar to DM mass (M$_*$/M$_{\text{halo}}$). Using a Markov Chain Monte Carlo approach, we test the NFW and the M$_*$/M$_{\text{halo}}$-dependent DC14 halo models against a sample of 147 galaxy rotation curves from the new {\it Spitzer} Photometry and Accurate Rotation Curves (SPARC) data set. These galaxies all have extended H{\small I} rotation curves from radio interferometry as well as accurate stellar mass density profiles from near-infrared photometry. The DC14 halo profile provides markedly better fits to the data compared to the NFW profile. Unlike NFW, the DC14 halo parameters found in our rotation curve fits naturally fall within two standard deviations of the mass-concentration relation predicted by $\Lambda$CDM and the stellar mass-halo mass relation inferred from abundance matching with few outliers. Halo profiles modified by baryonic processes are therefore more consistent with expectations from $\Lambda$ cold dark matter ($\Lambda$CDM) cosmology and provide better fits to galaxy rotation curves across a wide range of galaxy properties than do halo models that neglect baryonic physics. Our results offer a solution to the decade long cusp-core discrepancy.

We report a correlation between the radial acceleration traced by rotation curves and that predicted by the observed distribution of baryons. The same relation is followed by 2693 points in 153 galaxies with very different morphologies, masses, sizes, and gas fractions. The correlation persists even when dark matter dominates. Consequently, the dark matter contribution is fully specified by that of the baryons. The observed scatter is small and largely dominated by observational uncertainties. This radial acceleration relation is tantamount to a natural law for rotating galaxies.

We use the SPARC (Spitzer Photometry & Accurate Rotation Curves) database to study the relation between the central surface density of stars Sstar and dynamical mass Sdyn in 135 disk galaxies (S0 to dIrr). We find that Sdyn correlates tightly with Sstar over 4 dex. This central density relation can be described by a double power law. High surface brightness galaxies are consistent with a 1:1 relation, suggesting that they are self-gravitating and baryon dominated in the inner parts. Low surface brightness galaxies systematically deviate from the 1:1 line, indicating that the dark matter contribution progressively increases but remains tightly coupled to the stellar one. The observed scatter is small (~0.2 dex) and largely driven by observational uncertainties. The residuals show no correlations with other galaxy properties like stellar mass, size, or gas fraction.

We introduce SPARC (Spitzer Photometry & Accurate Rotation Curves): a sample of 175 nearby galaxies with new surface photometry at 3.6 um and high-quality rotation curves from previous HI/Halpha studies. SPARC spans a broad range of morphologies (S0 to Irr), luminosities (~5 dex), and surface brightnesses (~4 dex). We derive [3.6] surface photometry and study structural relations of stellar and gas disks. We find that both the stellar mass-HI mass relation and the stellar radius-HI radius relation have significant intrinsic scatter, while the HI mass-radius relation is extremely tight. We build detailed mass models and quantify the ratio of baryonic-to-observed velocity (Vbar/Vobs) for different characteristic radii and values of the stellar mass-to-light ratio (M/L) at [3.6]. Assuming M/L=0.5 Msun/Lsun (as suggested by stellar population models) we find that (i) the gas fraction linearly correlates with total luminosity, (ii) the transition from star-dominated to gas-dominated galaxies roughly corresponds to the transition from spiral galaxies to dwarf irregulars in line with density wave theory; and (iii) Vbar/Vobs varies with luminosity and surface brightness: high-mass, high-surface-brightness galaxies are nearly maximal, while low-mass, low-surface-brightness galaxies are submaximal. These basic properties are lost for low values of M/L=0.2 Msun/Lsun as suggested by the DiskMass survey. The mean maximum-disk limit in bright galaxies is M/L=0.7 Msun/Lsun at [3.6]. The SPARC data are publicly available and represent an ideal test-bed for models of galaxy formation.

The physical mechanisms driving star formation (SF) in galaxies are still not fully understood. Tidal dwarf galaxies (TDGs), made of gas ejected during galaxy interactions, seem to be devoid of dark matter and have a near-solar metallicity. The latter makes it possible to study molecular gas and its link to SF using standard tracers (CO, dust) in a peculiar environment. We present a detailed study of a nearby TDG in the Virgo Cluster, VCC 2062, using new high-resolution CO(1--0) data from the Plateau de Bure, deep optical imaging from the Next Generation Virgo Cluster Survey (NGVS), and complementary multiwavelength data. Until now, there was some doubt whether VCC 2062 was a true TDG, but the new deep optical images from the NGVS reveal a stellar bridge between VCC 2062 and its parent galaxy, NGC 4694, which is clear proof of its tidal origin. Several high-resolution tracers (\halpha, UV, 8~\mi, and 24~\mi) of the star formation rate (SFR) are compared to the molecular gas distribution as traced by the CO(1-0). Coupled with the SFR tracers, the NGVS data are used with the CIGALE code to model the stellar populations throughout VCC 2062, yielding a declining SFR in the recent past, consistent with the low \halpha/UV ratio, and a high burst strength. HI emission covers VCC 2062, whereas the CO is concentrated near the HI maxima. The CO peaks correspond to two very distinct regions: one with moderate SF to the NE and one with only slightly weaker CO emission but with nearly no SF. Even where SF is clearly present, the SFR is below the value expected from the surface density of the molecular and the total gas as compared to spiral galaxies and other TDGs. After discussing different possible explanations, we conclude that the low surface brightness is a crucial parameter to understand the low SFR.

In a LCDM cosmology, the baryonic Tully-Fisher relation (BTFR) is expected to show significant intrinsic scatter resulting from the mass-concentration relation of dark matter halos and the baryonic-to-halo mass ratio. We study the BTFR using a sample of 118 disc galaxies (spirals and irregulars) with data of the highest quality: extended HI rotation curves (tracing the outer velocity) and Spitzer photometry at 3.6 $\mu$m (tracing the stellar mass). Assuming that the stellar mass-to-light ratio (M*/L) is nearly constant at 3.6 $\mu$m, we find that the scatter, slope, and normalization of the BTFR systematically vary with the adopted M*/L. The observed scatter is minimized for M*/L > 0.5, corresponding to nearly maximal discs in high-surface-brightness galaxies and BTFR slopes close to ~4. For any reasonable value of M*/L, the intrinsic scatter is ~0.1 dex, below general LCDM expectations. The residuals show no correlations with galaxy structural parameters (radius or surface brightness), contrary to the predictions from some semi-analytic models of galaxy formation. These are fundamental issues for LCDM cosmology.

The mass discrepancy acceleration relation (MDAR) describes the coupling between baryons and dark matter (DM) in galaxies: the ratio of total-to-baryonic mass at a given radius anti-correlates with the acceleration due to baryons. The MDAR has been seen as a challenge to the $\Lambda$CDM galaxy formation model, while it can be explained by Modified Newtonian Dynamics. In this Letter we show that the MDAR arises in a $\Lambda$CDM cosmology once observed galaxy scaling relations are taken into account. We build semi-empirical models based on $\Lambda$CDM haloes, with and without the inclusion of baryonic effects, coupled to empirically motivated structural relations. Our models can reproduce the MDAR: specifically, a mass-dependent density profile for DM haloes can fully account for the observed MDAR shape, while a universal profile shows a discrepancy with the MDAR of dwarf galaxies with $\rm M^{\star}$$<$$\rm10^{9.5}M_{\odot}$, a further indication suggesting the existence of DM cores. Additionally, we reproduce slope and normalization of the baryonic Tully-Fisher relation (BTFR) with 0.17 dex scatter. These results imply that in $\Lambda$CDM (i) the MDAR is driven by structural scaling relations of galaxies and DM density profile shapes, and (ii) the baryonic fractions determined by the BTFR are consistent with those inferred from abundance-matching studies.

Collisional debris around interacting and post-interacting galaxies often display condensations of gas and young stars that can potentially form gravitationally bound objects: Tidal Dwarf Galaxies (TDGs). We summarise recent results on TDGs, which are originally published in Lelli et al. (2015, A&A). We study a sample of six TDGs around three different interacting systems, using high-resolution HI observations from the Very Large Array. We find that the HI emission associated to TDGs can be described by rotating disc models. These discs, however, would have undergone less than one orbit since the time of the TDG formation, raising the question of whether they are in dynamical equilibrium. Assuming that TDGs are in dynamical equilibrium, we find that the ratio of dynamical mass to baryonic mass is consistent with one, implying that TDGs are devoid of dark matter. This is in line with the results of numerical simulations where tidal forces effectively segregate dark matter in the halo from baryonic matter in the disc, which ends up forming tidal tails and TDGs.

Tidal dwarf galaxies (TDGs) are recycled objects that form within the collisional debris of interacting/merging galaxies. They are expected to be devoid of non-baryonic dark matter, since they can form only from dissipative material ejected from the discs of the progenitor galaxies. We investigate the gas dynamics in a sample of six bona-fide TDGs around three interacting and post-interacting systems: NGC 4694, NGC 5291, and NGC 7252 ("Atoms for Peace"). For NGC 4694 and NGC 5291 we analyse existing HI data from the Very Large Array (VLA), while for NGC 7252 we present new HI observations from the Jansky VLA together with long-slit and integral-field optical spectroscopy. For all six TDGs, the HI emission can be described by rotating disc models. These HI discs, however, have undergone less than a full rotation since the time of the interaction/merger event, raising the question of whether they are in dynamical equilibrium. Assuming that these discs are in equilibrium, the inferred dynamical masses are consistent with the observed baryonic masses, implying that TDGs are devoid of dark matter. This puts constraints on putative "dark discs" (either baryonic or non-baryonic) in the progenitor galaxies. Moreover, TDGs seem to systematically deviate from the baryonic Tully-Fisher relation. These results provide a challenging test for alternative theories like MOND.

This is a brief rebuttal to arXiv:1502.03821, which claims to provide the first observational proof of dark matter interior to the solar circle. We point out that this result is not new, and can be traced back at least a quarter century.

Recent studies have shown that starburst dwarf galaxies have steeply rising rotation curves in their inner parts, pointing to a close link between the intense star formation and a centrally concentrated mass distribution (baryons and dark matter). More quiescent dwarf irregulars typically have slowly rising rotation curves, although some "compact" irregulars with steep, inner rotation curves exist. We analyze archival Hubble Space Telescope images of two nearby "compact" irregular galaxies (NGC 4190 and NGC 5204), which were selected solely on the basis of their dynamical properties and their proximity. We derive their recent star-formation histories by fitting color-magnitude diagrams of resolved stellar populations, and find that the star-formation properties of both galaxies are consistent with those of known starburst dwarfs. Despite the small sample, this strongly reinforces the notion that the starburst activity is closely related to the inner shape of the potential well.

Strong bursts of star formation in galaxies may be triggered either by internal or external mechanisms. We study the distribution and kinematics of the HI gas in the outer regions of 18 nearby starburst dwarf galaxies, that have accurate star-formation histories from HST observations of resolved stellar populations. We find that starburst dwarfs show a variety of HI morphologies, ranging from heavily disturbed HI distributions with major asymmetries, long filaments, and/or HI-stellar offsets, to lopsided HI distributions with minor asymmetries. We quantify the outer HI asymmetry for both our sample and a control sample of typical dwarf irregulars. Starburst dwarfs have more asymmetric outer HI morphologies than typical irregulars, suggesting that some external mechanism triggered the starburst. Moreover, galaxies hosting an old burst (>100 Myr) have more symmetric HI morphologies than galaxies hosting a young one (<100 Myr), indicating that the former ones probably had enough time to regularize their outer HI distribution since the onset of the burst. We also investigate the nearby environment of these starburst dwarfs and find that most of them ($\sim$80$\%$) have at least one potential perturber at a projected distance <200 kpc. Our results suggest that the starburst is triggered either by past interactions/mergers between gas-rich dwarfs or by direct gas infall from the IGM.

For disk galaxies (spirals and irregulars), the inner circular-velocity gradient (inner steepness of the rotation curve) correlates with the central surface brightness with a slope of ~0.5. This implies that the central dynamical mass density scales almost linearly with the central baryonic density. Here I show that this empirical relation is consistent with a simple model where the central baryonic fraction f_bar(0) is fixed to 1 (no dark matter) and the observed scatter is due to differences in the baryonic mass-to-light ratio M_bar/L (ranging from 1 to 3 in the R-band) and in the characteristic thickness of the central stellar component dz (ranging from 100 to 500 pc). Models with lower baryonic fractions are possible, although they require some fine-tuning in the values of M_bar/L and dz. Regardless of the actual value of f_bar(0), the fact that different types of galaxies do not show strong variations in f_bar(0) is surprising, and may represent a challenge for models of galaxy formation in a LCDM cosmology.

We investigate the dynamics of starbursting dwarf galaxies, using both new and archival HI observations. We consider 18 nearby galaxies that have been resolved into single stars by HST observations, providing their star formation history and total stellar mass. We find that 9 objects have a regularly-rotating HI disk, 7 have a kinematically disturbed HI disk, and 2 show unsettled HI distributions. Two galaxies (NGC 5253 and UGC 6456) show a velocity gradient along the minor axis of the HI disk, that we interpret as strong radial motions. For galaxies with a regularly rotating disk we derive rotation curves, while for galaxies with a kinematically disturbed disk we estimate the rotation velocities in their outer parts. We derive baryonic fractions within about 3 optical scale lengths and find that, on average, baryons constitute at least 30$\%$ of the total mass. Despite the star formation having injected $\sim$10$^{56}$ ergs in the ISM in the last $\sim$500 Myr, these starbursting dwarfs have both baryonic and gas fractions similar to those of typical dwarf irregulars, suggesting that they did not eject a large amount of gas out of their potential wells.

UGC 4483 is a nearby Blue Compact Dwarf (BCD) galaxy. HST observations have resolved the galaxy into single stars and this has led to the derivation of its star formation history and to a direct estimate of its stellar mass. We have analysed archival VLA observations of the 21 cm line and found that UGC 4483 has a steeply-rising rotation curve which flattens in the outer parts at a velocity of ~20 km/s. Radial motions of ~5 km/s may also be present. As far as we know, UGC 4483 is the lowest-mass galaxy with a differentially rotating HI disk. The steep rise of the rotation curve indicates that there is a strong central concentration of mass. We have built mass models using the HST information on the stellar mass to break the disk-halo degeneracy: old stars contribute ~50$%$ of the observed rotation velocity at 2.2 disk scale-lengths. Baryons (gas and stars) constitute an important fraction of the total dynamical mass. These are striking differences with respect to typical dwarf irregular galaxies (dIrrs), which usually have slowly-rising rotation curves and are thought to be entirely dominated by dark matter. BCDs appear to be different from non-starbursting dIrrs in terms of their HI and stellar distributions and their internal dynamics. To their high central surface brightnesses and high central HI densities correspond strong central rotation-velocity gradients. This implies that the starburst is closely related with the gravitational potential and the concentration of gas. We discuss the implications of our results on the properties of the progenitors/descendants of BCDs.

For a rotating galaxy, the inner circular-velocity gradient d_{R}V(0) provides a direct estimate of the central dynamical mass density, including gas, stars, and dark matter. We consider 60 low-mass galaxies with high-quality HI and/or stellar rotation curves (including starbursting dwarfs, irregulars, and spheroidals), and estimate d_{R}V(0) as V(R_d)/R_d, where R_d is the galaxy scale-length. For gas-rich dwarfs, we find that V(R_d)/R_d correlates with the central surface brightness mu(0), the mean atomic gas surface density Sigma_gas, and the star formation rate surface density Sigma_SFR. Starbursting galaxies, such as blue compact dwarfs (BCDs), generally have higher values of V(R_d)/R_d than dwarf irregulars, suggesting that the starburst is closely related to the inner shape of the potential well. There are, however, some "compact" irregulars with values of V(R_d)/R_d similar to BCDs. Unless a redistribution of mass takes place, BCDs must evolve into compact irregulars. Rotating spheroidals in the Virgo cluster follow the same correlation between V(R_d)/R_d and mu(0) as gas-rich dwarfs. They have values of V(R_d)/R_d comparable to those of BCDs and compact irregulars, pointing at evolutionary links between these types of dwarfs. Finally, we find that, similarly to spiral galaxies and massive starbursts, the star-formation activity in dwarfs can be parametrized as Sigma_SFR = epsilon*Sigma_gas/t_orb, where t_orb is the orbital time and epsilon = 0.02.

For disc galaxies, a close relation exists between the distribution of light and the shape of the rotation curve. We quantify this relation by measuring the inner circular-velocity gradient d_{R}V(0) for spiral and irregular galaxies with high-quality rotation curves. We find that d_{R}V(0) correlates with the central surface brightness mu_{0} over more than two orders of magnitude in d_{R}V(0) and four orders of magnitudes in mu_{0}. This is a scaling-relation for disc galaxies. It shows that the central stellar density of a galaxy closely relates to the inner shape of the potential well, also for low-luminosity and low-surface-brightness galaxies that are expected to be dominated by dark matter.

I Zw 18 is a prototype Blue Compact Dwarf (BCD), characterized by a strong starburst and extremely low metallicity (Z ~ 0.02 Zsun). It has long been considered a candidate young galaxy in the Local Universe, but recent studies indicate the presence of old stars. We analysed archival VLA observations of the 21 cm line and found that the HI associated to the starburst region forms a compact fast-rotating disk. The HI column densities are very high, up to ~50-100 Msun/pc^2 (~0.6-1.2 x 10^22 atoms/cm^2). The rotation curve is flat with a steep rise in the inner parts, indicating the presence of a strong central concentration of mass. Mass models with a dark matter halo show that baryons may dominate the gravitational potential in the inner regions. A radial inflow/outflow motion of ~15 km/s is also present. I Zw 18 appears structurally different from typical dwarf irregulars in terms of gas distribution, stellar distribution and dynamics. It may be considered as a "miniature" high-surface-brightness disk galaxy. These dynamical properties must be tightly related to the starburst. They also shed new light on the question of the descendants of BCDs. There is also extended HI emission towards the outlying stellar complex I Zw 18 C and a ~13.5 kpc HI tail. An interaction/merger between gas-rich dwarfs is the most likely explanation for the starburst.

Giant low surface brightness (GLSB) galaxies are commonly thought to be massive, dark matter dominated systems. However, this conclusion is based on highly uncertain rotation curves. We present here a new study of two prototypical GLSB galaxies: Malin 1 and NGC 7589. We re-analysed existing HI observations and derived new rotation curves, which were used to investigate the distributions of luminous and dark matter in these galaxies. In contrast to previous findings, the rotation curves of both galaxies show a steep rise in the central parts, typical of high surface brightness (HSB) systems. Mass decompositions with a dark matter halo show that baryons may dominate the dynamics of the inner regions. Indeed, a "maximum disk" fit gives stellar mass-to-light ratios in the range of values typically found for HSB galaxies. These results, together with other recent studies, suggest that GLSB galaxies are systems with a double structure: an inner HSB early-type spiral galaxy and an outer extended LSB disk. We also tested the predictions of MOND: the rotation curve of NGC 7589 is reproduced well, whereas Malin 1 represents a challenging test for the theory.