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

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.

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.

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.

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.

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.

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.

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.

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.

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 $\Theta_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 ($\rho_{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.

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.

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.

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 dwarf galaxy NGC1052-DF2 has recently been identified as potentially lacking dark matter. If correct, this could be a challenge for MOND, which predicts that low surface brightness galaxies should evince large mass discrepancies. However, the correct prediction of MOND depends on both the internal field of the dwarf and the external field caused by its proximity to the giant elliptical NGC1052. Taking both into consideration under plausible assumptions, we find $\sigma_{\rm MOND} = 13.4^{+4.8}_{-3.7}\;\mathrm{km}\,\mathrm{s}^{-1}$. This is only marginally higher than the claimed 90\% upper limit on the velocity dispersion ($\sigma < 10.5\;\mathrm{km}\,\mathrm{s}^{-1}$), and compares well with the observed root mean square velocity dispersion ($\sigma = 14.3\;\mathrm{km}\,\mathrm{s}^{-1}$). We also discuss a few caveats on both the observational and theoretical side. On the theory side, the internal virialization time in this dwarf may be longer that the time scale of variation of the external field. On the observational side, the paucity of data and their large uncertainties call for further analysis of the velocity dispersion of NGC1052-DF2, to check whether it poses a challenge to MOND or is a success thereof.

The recently reported detection of redshifted 21cm absorption at $z \approx 17$ is a significant advance in the exploration of the cosmic dark ages. The observed signal ($T_{\mathrm{21}} \approx -0.5$ K with the limit $T_{\mathrm{21}} < -0.3$ K at 99\% c.l.) is anomalously strong for $\Lambda$CDM, which predicts at most $-0.24$ K. Here I show that the strong observed signal is expected in a purely baryonic universe.

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.

MHONGOOSE is a deep survey of the neutral hydrogen distribution in a representative sample of 30 nearby disk and dwarf galaxies with HI masses from 10^6 to ~10^{11} M_sun, and luminosities from M_R ~ -12 to M_R ~ -22. The sample is selected to uniformly cover the available range in log(M_HI). Our extremely deep observations, down to HI column density limits of well below 10^{18} cm^{-2} - or a few hundred times fainter than the typical HI disks in galaxies - will directly detect the effects of cold accretion from the intergalactic medium and the links with the cosmic web. These observations will be the first ever to probe the very low-column density neutral gas in galaxies at these high resolutions. Combination with data at other wavelengths, most of it already available, will enable accurate modelling of the properties and evolution of the mass components in these galaxies and link these with the effects of environment, dark matter distribution, and other fundamental properties such as halo mass and angular momentum. MHONGOOSE can already start addressing some of the SKA-1 science goals and will provide a comprehensive inventory of the processes driving the transformation and evolution of galaxies in the nearby universe at high resolution and over 5 orders of magnitude in column density. It will be a Nearby Galaxies Legacy Survey that will be unsurpassed until the advent of the SKA, and can serve as a highly visible, lasting statement of MeerKAT's capabilities.

The existence of a spatially thin, kinematically coherent Disk of Satellites (DoS) around the Milky Way (MW), is a problem that often garners vivacious debate in the literature or at scientific meetings. One of the most recent incarnations of this wrangle occurred with two papers by Maji et al, who argued that these structures "maybe a misinterpretation of the data". These claims are in stark contrast to previous works. Motivated by this and other recent publications on this problem, we discuss necessary considerations to make, observational effects to consider, and pitfalls to avoid when investigating satellite galaxy planes such as the MW's DoS. In particular, we emphasize that conclusions need to have a statistical basis including a determination of the significance of satellite alignments, observational biases must not be ignored, and measurement errors (e.g. for proper motions) need to be considered. We discuss general problems faced by attempts to determine the dynamical stability of the DoS via orbit integrations of MW satellite galaxies, and demonstrate that to interpret simulations, it is helpful to compare them with a null case of isotropically distributed satellite positions and velocities. Based on these criteria, we find that the conclusions of Maji et al. do not hold up to scrutiny, and that their hydrodynamic cosmological simulation of a single host shows no evidence for a significant kinematic coherence among the simulated satellite galaxies, in contrast to the observed MW system.

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.

The Candidate Cluster and Protocluster Catalog (CCPC) contains 218 galaxy overdensities com- posed of more than 2000 galaxies with spectroscopic redshifts spanning the first few Gyrs after the Big Bang (2.0 < z < 6.6). We use Spitzer archival data to track the underlying stellar mass of these overdense regions in various temporal cross sections by building rest-frame near-infrared luminosity functions across the span of redshifts.This exercise maps the stellar growth of protocluster galaxies, as halos in the densest environments should be the most massive from hierarchical accretion. The characteristic apparent magnitude, m*(z), is relatively flat from 2.0 < z < 6.6, consistent with a passive evolution of an old stellar population. This trend maps smoothly to lower redshift results of cluster galaxies from other works. We find no difference in the luminosity

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.

Crater II is an unusual object among the dwarf satellite galaxies of the Local Group in that it has a very large size for its small luminosity. This provides a strong test of MOND, as Crater II should be in the deep MOND regime ($g_{in} \approx 34\;\mathrm{km}^2\,\mathrm{s}^{-2}\,\mathrm{kpc}^{-1} \ll a_0 = 3700\;\mathrm{km}^2\,\mathrm{s}^{-2}\,\mathrm{kpc}^{-1}$). Despite its great distance ($\approx 120$ kpc) from the Milky Way, the external field of the host ($g_{ex} \approx 282\; \mathrm{km}^2\,\mathrm{s}^{-2}\,\mathrm{kpc}^{-1}$) comfortably exceeds the internal field. Consequently, Crater II should be subject to the external field effect, a feature unique to MOND. This leads to the prediction of a very low velocity dispersion: $\sigma_{efe} = 2.1^{+0.9}_{-0.6}\;\mathrm{km}\,\mathrm{s}^{-1}$.

The Candidate Cluster and Protocluster Catalog (CCPC) is a list of objects at redshifts z > 2 composed of galaxies with spectroscopically confirmed redshifts that are coincident on the sky and in redshift. These protoclusters are identified by searching for groups in volumes corresponding to the expected size of the most massive protoclusters at these redshifts. In CCPC1 we identified 43 candidate protoclusters among 14,000 galaxies between 2.74 < z < 3.71. Here we expand our search to more than 40,000 galaxies with spectroscopic redshifts z > 2.00, resulting in an additional 173 candidate structures. The most significant of these are 36 protoclusters with overdensities {\delta} > 7. We also identify three large proto-supercluster candidates containing multiple protoclusters at z = 2.3, 3.5 and z = 6.56. Eight candidates with N > 10 or more galaxies are found at redshifts z > 4.0. The last system in the catalog is the most distant spectroscopic protocluster candidate known to date at z = 6.56.

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.

We have developed a search methodology to identify galaxy protoclusters at $z>2.74$, and implemented it on a sample of $\sim$14,000 galaxies with previously measured redshifts. The results of this search are recorded in the Candidate Cluster and Protocluster Catalog (CCPC). The catalog contains 12 clusters that are highly significant overdensities ($\delta_{gal}>7$), 6 of which are previously known. We also identify another 31 candidate protoclusters (including 4 previously identified structures) of lower overdensity. CCPC systems vary over a wide range of physical sizes and shapes, from small, compact groups to large, extended, and filamentary collections of galaxies. This variety persists over the range from $z=3.71$ to $z=2.74$. These structures exist as galaxy overdensities ($\delta_{gal}$) with a mean value of 2, similar to the values found for other protoclusters in the literature. The median number of galaxies for CCPC systems is 11. Virial mass estimates are large for these redshifts, with thirteen cases apparently having $M > 10^{15}\, M_{\odot}$. If these systems are virialized, such masses would pose a challenge to $\Lambda$CDM.

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 distribution of the Galactic disk is constructed from the terminal velocity curve and the mass discrepancy-acceleration relation. Mass models numerically quantifying the detailed surface density profiles are tabulated. For $R_0 = 8$ kpc, the models have stellar mass $5 < M_* < 6 \times 10^{10}$ M$_{\odot}$, scale length $2.0 \le R_d \le 2.9$ kpc, LSR circular velocity $222 \le \Theta_0 \le 233$ km s$^{-1}$, and solar circle stellar surface density $34 \le \Sigma_d(R_0) \le 61$ M$_{\odot}$ pc$^{-2}$. The present inter-arm location of the solar neighborhood may have a somewhat lower stellar surface density than average for the solar circle. The Milky Way appears to be a normal spiral galaxy that obeys scaling relations like the Tully-Fisher relation, the size-mass relation, and the disk maximality-surface brightness relation. The stellar disk is maximal, and the spiral arms are massive. The bumps and wiggles in the terminal velocity curve correspond to known spiral features (e.g., the Centaurus Arm is a $\sim 50\%$ overdensity). The rotation curve switches between positive and negative over scales of hundreds of parsecs. The rms amplitude $\langle$$|$$dV/dR$$|^2$$\rangle$$^{1/2} \approx 14$ km s$^{-1}$ kpc$^{-1}$, implying that commonly neglected terms in the Jeans equations may be non-negligible. The spherically averaged local dark matter density is $\rho_{0,DM} \approx 0.009$ M$_{\odot}$ pc$^{-3}$ (0.3 GeV cm$^{-3}$). Adiabatic compression of the dark matter halo may help reconcile the Milky Way with the $c$-$V_{200}$ relation expected in $\Lambda$CDM while also helping to mitigate the too big to fail problem, but it remains difficult to reconcile the inner bulge/bar dominated region with a cuspy halo. We note that NGC 3521 is a near twin to the Milky Way, having a similar luminosity, scale length, and rotation curve.

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.

The evidence that stellar systems surrounding the Milky Way (MW) are distributed in a Vast Polar Structure (VPOS) may be observationally biased by satellites detected in surveys of the northern sky. The recent discoveries of more than a dozen new systems in the southern hemisphere thus constitute a critical test of the VPOS phenomenon. We report that the new objects are located close to the original VPOS, with half of the sample having offsets less than 20 kpc. The positions of the new satellite galaxy candidates are so well aligned that the orientation of the revised best-fitting VPOS structure is preserved to within 9 degrees and the VPOS flattening is almost unchanged (31 kpc height). Interestingly, the shortest distance of the VPOS plane from the MW center is now only 2.5 kpc, indicating that the new discoveries balance out the VPOS at the Galactic center. The vast majority of the MW satellites are thus consistent with sharing a similar orbital plane as the Magellanic Clouds, confirming a hypothesis proposed by Kunkel & Demers and Lynden-Bell almost 40 years ago. We predict the absolute proper motions of the new objects assuming they orbit within the VPOS. Independent of the VPOS results we also predict the velocity dispersions of the new systems under three distinct assumptions: that they (i) are dark-matter-free star clusters obeying Newtonian dynamics, (ii) are dwarf satellites lying on empirical scaling relations of galaxies in dark matter halos, and (iii) obey MOND.

We present WFC3 observations of three low surface brightness (LSB) galaxies from the Schombert et. al LSB catalog that are within 11 Mpc of the Milky Way. Deep imaging at F336W, F555W and F814W allow the construction of the V-I color-magnitude diagrams (CMD) to M_I = -2. Overall 1869, 465 and 501 stellar sources are identified in the three LSB galaxies F415-3, F608-1 and F750-V1. The spatial distribution of young blue stars matches the H-alpha maps from ground-based imaging, indicating that star formation in LSB galaxies follows the same style as in other irregular galaxies. Several star complexes are identified, matching regions of higher surface brightness as seen from ground-based imaging. The color-magnitude diagrams for each LSB galaxy has the similar morphology to Local Volume (LV) dwarf galaxies, i.e. a blue main sequence, blue and red He burning branches and asymptotic giant branch (AGB) stars. The LSB CMD's distinguish themselves from nearby dwarf CMD's by having a higher proportional of blue main sequence stars and fewer AGB stars than expected from their mean metallicities. Current [Fe/H] values below -0.6 are deduced from the position of the red helium-burning branch (rHeB) stars in the V-I diagram. The distribution of stars on the blue helium-burning branch (bHeB) and rHeB from the U-V and V-I CMD indicate a history of constant star formation for the last 100 Myrs.

We present near-infrared observations of the environments around three radio-loud sources (MG1 J0442+0202, 3C 068.2, and MS 1426.9+1052) at redshifts z=1.10,1.57, and 1.83 (respectively), that are surrounded by near-infrared galaxy overdensities. Overdensities with respect to field counts were found to be significant up to 19-sigma, with twelve times the expected number of galaxies within the inner regions of the densest proto-cluster. Color-magnitude relations are constructed in K_s, J-K_s, with each candidate cluster exhibiting a feature consistent with the beginnings of a red sequence. Galaxy models based on the redshift of the radio source are used to compare expected color-magnitude relations for a given formation epoch with the observed red sequence of each candidate, and are found to be consistent with an old (z_f > 5) formation epoch for a few bright, red galaxies on the red sequence.

This article explores the agreement between the predictions of Modified Newtonian Dynamics (MOND) and the rotation curves and stellar velocity dispersion profiles measured by the DiskMass Survey. A bulge-disk decomposition was made for each of the thirty published galaxies, and a MOND Poisson solver was used to simultaneously compute, from the baryonic mass distributions, model rotation curves and vertical velocity dispersion profiles, which were compared to the measured values. The two main free parameters, the stellar disk's mass-to-light ratio ($M/L$) and its exponential scale-height ($h_z$), were estimated by Markov Chain Monte Carlo modelling. The average best-fit K-band stellar mass-to-light ratio was $M/L \simeq 0.55 \pm 0.15$. However, to match the DiskMass Survey data, the vertical scale-heights would have to be in the range $h_z=200$ to $400$ pc which is a factor of two lower than those derived from observations of edge-on galaxies with a similar scale-length. The reason is that modified gravity versions of MOND characteristically require a larger $M/L$ to fit the rotation curve in the absence of dark matter and therefore predict a stronger vertical gravitational field than Newtonian models. It was found that changing the MOND acceleration parameter, the shape of the velocity dispersion ellipsoid, the adopted vertical distribution of stars, as well as the galaxy inclination, within any realistic range, all had little impact on these results.

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.

We estimate the stellar masses of disk galaxies with two independent methods: a photometrically self-consistent color$-$mass-to-light ratio relation (CMLR) from population synthesis models, and the Baryonic Tully-Fisher relation (BTFR) calibrated by gas rich galaxies. These two methods give consistent results. The CMLR correctly converts distinct Tully-Fisher relations in different bands into the same BTFR. The BTFR is consistent with $M_b \propto V_f^4$ over nearly six decades in mass, with no hint of a change in slope over that range. The intrinsic scatter in the BTFR is negligible, implying that the IMF of disk galaxies is effectively universal. The gas rich BTFR suggests an absolute calibration of the stellar mass scale that yields nearly constant mass-to-light ratios in the near-infrared (NIR): $0.57\;M_{\odot}/L_{\odot}$ in $K_s$ and $0.45\;M_{\odot}/L_{\odot}$ at $3.6\mu$. There is only modest intrinsic scatter ($\sim 0.12$ dex) about these typical values. There is no discernible variation with color or other properties: the NIR luminosity is a good tracer of stellar mass.

I review the connection between dynamics and the baryonic mass distribution in rotationally supported galaxies. The enclosed dynamical mass-to-light ratio increases with decreasing galaxy luminosity and surface brightness. The correlation with surface brightness appears to be the more fundamental, with the dependence on luminosity following simply from the weaker correlation between luminosity and surface brightness. In addition to this global relation, there is also a local relation between the amplitude of the mass discrepancy and the acceleration predicted by the observed distribution of baryons. I provide an empirical calibration of this mass discrepancy-acceleration relation. The data are consistent with the operation of a singe effective force law in disk galaxies, making this relation tantamount to a natural law. I further provide formulae by which the radial dark matter distribution can be estimated from surface photometry. The form of the dark matter halo depends uniquely on the distribution of baryons in each galaxy, and in general is neither a cusp nor a core. It remains difficult to see how galaxy formation models can reproduce the observed behavior, which is uniquely predicted by MOND.

A series of population models are designed to explore the star formation history of gas-rich, low surface brightness (LSB) galaxies. LSB galaxies are unique in having properties of very blue colors, low H$\alpha$ emission and high gas fractions that indicated a history of constant star formation (versus the declining star formation models used for most spirals and irregulars). The model simulations use an evolving multi-metallicity composite population that follows a chemical enrichment scheme based on Milky Way observations. Color and time sensitive stellar evolution components (i.e., BHB, TP-AGB and blue straggler stars) are included, and model colors are extended into the Spitzer wavelength regions for comparison to new observations. In general, LSB galaxies are well matched to the constant star formation scenario with the variation in color explained by a fourfold increase/decrease in star formation over the last 0.5 Gyrs (i.e., weak bursts). Early-type spirals, from the S$^4$G sample, are better fit by a declining star formation model where star formation has decreased by 40% in the last 12 Gyrs.

We combine Spitzer $3.6\mu$ observations of a sample of disk galaxies spanning over 10 magnitudes in luminosity with optical luminosities and colors to test population synthesis prescriptions for computing stellar mass. Many commonly employed models fail to provide self-consistent results: the stellar mass estimated from the luminosity in one band can differ grossly from that of another band for the same galaxy. Independent models agree closely in the optical ($V$-band), but diverge at longer wavelengths. This effect is particularly pronounced in recent models with substantial contributions from TP-AGB stars. We provide revised color--mass-to-light ratio relations that yield self-consistent stellar masses when applied to real galaxies. The $B-V$ color is a good indicator of the mass-to-light ratio. Some additional information is provided by $V-I$, but neither it nor $J-K_s$ are particularly useful for constraining the mass-to-light ratio on their own. In the near-infrared, the mass-to-light ratio depends weakly on color, with typical values of $0.6\; \mathrm{M}_{\odot}/\mathrm{L}_{\odot}$ in the $K_s$-band and $0.47\; \mathrm{M}_{\odot}/\mathrm{L}_{\odot}$ at $3.6\mu$.

Sub-halos in dark-matter-based cosmological simulations tend to be distributed approximately isotropically around their host. The existence of highly flattened, co-orbiting planes of satellite galaxies has therefore been identified as a possible problem for these cosmological models, but so far studies have not considered the hosts' environments. That satellite planes are now known around both major galaxies in the Local Group raises the question whether they are more likely around paired hosts. In a first attempt to investigate this possibility we focus on the flattening and orbital coherence of the 11 brightest satellite galaxies of the vast polar structure (VPOS) around the Milky Way (MW). We search for VPOS analogs in the ELVIS suite of cosmological simulations, which consist of 24 paired and 24 isolated host halos. We do not find significant differences between the properties of sub-halo distributions around paired and isolated hosts. The observed flattening and the observed orbital alignment are each reproduced by only 0.2 to 2 per cent of paired and isolated systems incorporating the obscuration of satellites by randomly oriented galactic discs. Only one of all 4800 analyzed realizations (0.02 per cent) reproduces both parameters simultaneously, but the average orbital pole of this sub-halo system does not align as well with the normal to the plane fit as observed. That the MW is part of a galaxy pair thus does not help in explaining the existence of the VPOS if the satellite galaxies are identified with sub-halos found in dissipationless simulations.

The concordance model of cosmology, LCDM, provides a satisfactory description of the evolution of the universe and the growth of large scale structure. Despite considerable effort, this model does not at present provide a satisfactory description of small scale structure and the dynamics of bound objects like individual galaxies. In contrast, MOND provides a unique and predictively successful description of galaxy dynamics, but is mute on the subject of cosmology. Here I briefly review these contradictory world views, emphasizing the wealth of distinct, interlocking lines of evidence that went into the development of LCDM while highlighting the practical impossibility that it can provide a satisfactory explanation of the observed MOND phenomenology in galaxy dynamics. I also briefly review the baryon budget in groups and clusters of galaxies where neither paradigm provides an entirely satisfactory description of the data. Relatively little effort has been devoted to the formation of structure in MOND; I review some of what has been done. The amplitude ratio of the first to second peak in the CMB power spectrum was correctly predicted a priori, but the third peak is more natural to LCDM. MOND anticipates that structure forms more quickly than in LCDM. This motivated the prediction that reionization would happen earlier in MOND than originally expected in LCDM, as subsequently observed. This also provides a natural explanation for massive, early clusters of galaxies and large, empty voids. However, it is far from obvious that the mass spectrum of galaxy clusters or the power spectrum of galaxies can be explained in MOND, two things that LCDM does well. Critical outstanding issues are the development of an acceptable relativistic parent theory for MOND, and the reality of the non-baryonic dark matter of LCDM. Do suitable dark matter particles exist, or are they a modern aether?

We present new Spitzer 3.6 micron observations of a sample of disk galaxies spanning over 10 magnitudes in luminosity and ranging in gas fraction from ~10% to over 90%. We use these data to test population synthesis prescriptions for computing stellar mass. Many commonly employed models fail to provide self-consistent stellar masses in the sense that the stellar mass estimated from the optical luminosity typically exceeds that estimated from the near-infrared (NIR) luminosity. This problem is present in models both with and without TP-AGB stars, but is more severe in the former. Self-consistency can be achieved if NIR mass-to-light ratios are approximately constant with a mean value near 0.5 Msun/Lsun at 3.6 microns. We use the Baryonic Tully-Fisher relation calibrated by gas rich galaxies to provide an independent estimate of the color-mass to light ratio relation. This approach also suggests that the typical 3.6 micron mass-to-light ratio is 0.5 (0.65 in the K band) for rotationally supported galaxies. These values are consistent with a Kroupa IMF.

The recently discovered dwarf galaxy Perseus I appears to be associated with the dominant plane of non-satellite galaxies in the Local Group (LG). We predict its velocity dispersion and those of the other isolated dSphs Cetus and Tucana to be 6.5, 8.2, and 5.5 km/s, respectively. The NGC 3109 association, including the recently discovered dwarf galaxy Leo P, aligns with the dwarf galaxy structures in the LG such that all known nearby non-satellite galaxies in the northern Galactic hemisphere lie in a common thin plane (rms height 53 kpc; diameter 1.2 Mpc). This plane has an orientation similar to the preferred orbital plane of the Milky Way (MW) satellites in the vast polar structure. Five of seven of these northern galaxies were identified as possible backsplash objects, even though only about one is expected from cosmological simulations. This may pose a problem, or instead the search for local backsplash galaxies might be identifying ancient tidal dwarf galaxies expelled in a past major galaxy encounter. The NGC 3109 association supports the notion that material preferentially falls towards the MW from the Galactic south and recedes towards the north, as if the MW were moving through a stream of dwarf galaxies.

We compile a sample of dispersion-supported and rotation-supported galaxies for which the dynamical mass enclosed at 500 pc can be measured. We find that this dynamical quantity increases only slowly with stellar mass (M$_{total}$(<500pc) $\propto$ M$_{stars}$$^{0.16}$) over ~9 decades in baryonic mass and ~3 in length scale, with a sudden upturn at the highest masses (M$_{stars}$ > 10$^{10}$M$_{sun}$). This upturn occurs for the earliest type (S0 and Sa) disk galaxies, and is consistent with the "extra" mass within 500 pc being predominantly stellar. This sudden change may be indicative of different bulge formation mechanisms between early and late type disks. We also compare the data to the results of recent galaxy formation simulations. These broadly track the observed trend, with some tendency to put too much mass in the inner 500 pc, depending on the details of each simulation.

We utilize Young's algorithm to model the adiabatic compression of the dark matter haloes of galaxies in the THINGS survey to determine the relationship between the halo fit to the rotation curve and the corresponding primordial halo prior to compression. Young's algorithm conserves radial action and angular momentum, resulting in less halo compression than more widely utilized approximations. We find that estimates of the parameters of the NFW haloes fit to the current dark matter distribution systematically overestimate the concentration and underestimate the virial velocity of the corresponding primordial halo. It is the latter that is predicted by dark matter simulations; so accounting for compression is a necessary step for evaluating whether massive galaxies are consistent with dark matter-only simulations. The inferred primordial haloes broadly follow the c-V200 relation expected in a LCDM cosmogony, but often scatter to lower concentrations. We are unable to obtain fits at all for those galaxies whose current dark matter haloes are poorly described by the NFW form. We thus find a mixed bag: some galaxies are reasonably well described by adiabatic compression within a primordial NFW halo, while others require additional mechanisms that reduce the density of dark matter below the primordial initial condition.

Surface photometry at 3.6$\mu$m is presented for 61 low surface brightness (LSB) galaxies ($\mu_o < 19$ 3.6$\mu$m mag arcsecs$^{-2}$). The sample covers a range of luminosity from $-$11 to $-$22 in $M_{3.6}$ and size from 1 to 25 kpc. The morphologies in the mid-IR are comparable to those in the optical with 3.6$\mu$m imaging reaches similar surface brightness depth as ground-based optical imaging. A majority of the resulting surface brightness profiles are single exponential in shape with very few displaying upward or downward breaks. The mean $V-3.6$ color of LSB is 2.3 with a standard deviation of 0.5. Color-magnitude and two color diagrams are well matched to models of constant star formation, where the spread in color is due to small changes in the star formation rate (SFR) over the last 0.5 Gyrs as also suggested by the specific star formation rate measured by H$\alpha$.

The Lambda-CDM cosmological model is succesful at reproducing various independent sets of observations concerning the large-scale Universe. This model is however currently, and actually in principle, unable to predict the gravitational field of a galaxy from it observed baryons alone. Indeed the gravitational field should depend on the relative contribution of the particle dark matter distribution to the baryonic one, itself depending on the individual assembly history and environment of the galaxy, including a lot of complex feedback mechanisms. However, for the last thirty years, Milgrom's formula, at the heart of the MOND paradigm, has been consistently succesful at predicting rotation curves from baryons alone, and has been resilient to all sorts of observational tests on galaxy scales. We show that the few individual galaxy rotation curves that have been claimed to be highly problematic for the predictions of Milgrom's formula, such as Holmberg II or NGC 3109, are actually false alarms. We argue that the fact that it is actually possible to predict the gravitational field of galaxies from baryons alone presents a challenge to the current Lambda-CDM model, and may indicate a breakdown of our understanding of gravitation and dynamics, and/or that the actual lagrangian of the dark sector is very different and richer than currently assumed. On the other hand, it is obvious that any alternative must also, in fine, reproduce the successes of the Lambda-CDM model on large scales, where this model is so well-tested that it presents by itself a challenge to any such alternative.

We employ recently published measurements of the velocity dispersions in the newly discovered dwarf satellite galaxies of Andromeda to test our previously published predictions of this quantity. The data are in good agreement with our specific predictions for each dwarf made a priori with MOND, with reasonable stellar mass-to-light ratios, and no dark matter, while Newtonian dynamics point to quite large mass discrepancies in these systems. MOND distinguishes between regimes where the internal field of the dwarf, or the external field of the host, dominates. The data appear to recognize this distinction, which is a unique feature of MOND not explicable in LCDM.

The luminosities, colors and H\alpha\ emission for 429 HII regions in 54 LSB galaxies are presented. While the number of HII regions per galaxy is lower in LSB galaxies compared to star-forming irregulars and spirals, there is no indication that the size or luminosity function of HII regions differs from other galaxy types. The lower number of HII regions per galaxy is consistent with their lower total star formation rates. The fraction of total $L_{H\alpha}$ contributed by HII regions varies from 10 to 90% in LSB galaxies (the rest of the H$\alpha$ emission being associated with a diffuse component) with no correlation with galaxy stellar or gas mass. Bright HII regions have bluer colors, similar to the trend in spirals; their number and luminosities are consistent with the hypothesis that they are produced by the same HII luminosity function as spirals. Comparison with stellar population models indicates that the brightest HII regions in LSB galaxies range in cluster mass from a few $10^3 M_{\sun}$ (e.g., $\rho$ Oph) to globular cluster sized systems (e.g., 30 Dor) and that their ages are consistent with clusters from 2 to 15 Myrs old. The faintest HII regions are comparable to those in the LMC powered by a single O or B star. Thus, star formation in LSB galaxies covers the full range of stellar cluster mass.

We examine the formation of clusters of galaxies in numerical simulations of a QUMOND cosmogony with massive sterile neutrinos. Clusters formed in these exploratory simulations develop higher velocities than those found in {\Lambda}CDM simulations. The bulk motions of clusters attain about 1000 km/s by low redshift, comparable to observations whereas {\Lambda}CDM simulated clusters tend to fall short. Similarly, high pairwise velocities are common in cluster-cluster collisions like the Bullet cluster. There is also a propensity for the most massive clusters to be larger in QUMOND and to appear earlier than in {\Lambda}CDM, potentially providing an explanation for 'pink elephants' like El Gordo. However, it is not obvious that the cluster mass function can be recovered.

Spiral galaxies must acquire gas to maintain their observed level of star formation beyond the next few billion years (Leroy et al. 2008). A source of this material may be the gas that resides between galaxies, but our understanding of the state and distribution of this gas is incomplete (Shull et al. 2012). Radio observations (Braun & Thilker 2004) of the Local Group of galaxies have revealed hydrogen gas extending from the disk of the galaxy M31 at least halfway to M33. This feature has been interpreted to be the neutral component of a condensing intergalactic filament (Dav\'e et al. 2001) which would be able to fuel star formation in M31 and M33, but simulations suggest that such a feature could also result from an interaction between both galaxies within the past few billion years (Bekki 2008). Here we report radio observations showing that about 50 per cent percent of this gas is composed of clouds, while the rest is distributed in an extended, diffuse component. The clouds have velocities comparable to those of M31 and M33, and have properties suggesting they are unrelated to other Local Group objects. We conclude that the clouds are likely to be transient condensations of gas embedded in an intergalactic filament and are therefore a potential source of fuel for future star formation in M31 and M33.

We compare the recently published velocity dispersions for seventeen Andromeda dwarf spheroidals with estimates of the MOND predictions, based on the luminosities of these dwarfs, with reasonable stellar M/L values, and no dark matter. We find that the two are consistent within the uncertainties. We further predict the velocity dispersions of another ten dwarfs for which only photometric data are currently available.

The Universe on large scales is well described by the Lambda-CDM cosmological model. There however remain some heavy clouds on our global understanding, especially on galaxy scales, which we review here. While some of these clouds might perhaps disappear through small compensatory adjustments of the model, such as changing the mass of the dark matter particles or accounting better for baryonic physics, others should rather be taken as strong indications that the physics of the dark sector is, at the very least, much richer and complex than currently assumed, and that our understanding of gravity and dynamics might also be at play. For instance, the empirically well-tested MOND phenomenology in galaxies, whatever its final explanation, should be understood in any model of galaxy formation and dynamics. Current alternatives to Lambda-CDM however bring with them many unsolved questions and challenges.

A wealth of astronomical data indicate the presence of mass discrepancies in the Universe. The motions observed in a variety of classes of extragalactic systems exceed what can be explained by the mass visible in stars and gas. Either (i) there is a vast amount of unseen mass in some novel form - dark matter - or (ii) the data indicate a breakdown of our understanding of dynamics on the relevant scales, or (iii) both. Here, we first review a few outstanding challenges for the dark matter interpretation of mass discrepancies in galaxies, purely based on observations and independently of any alternative theoretical framework. We then show that many of these puzzling observations are predicted by one single relation - Milgrom's law - involving an acceleration constant (or a characteristic surface density) of the order of the square-root of the cosmological constant in natural units. This relation can at present most easily be interpreted as the effect of a single universal force law resulting from a modification of Newtonian dynamics (MOND) on galactic scales. We exhaustively review the current observational successes and problems of this alternative paradigm at all astrophysical scales, and summarize the various theoretical attempts (TeVeS, GEA, BIMOND, and others) made to effectively embed this modification of Newtonian dynamics within a relativistic theory of gravity.

Here we present a new particle-mesh galactic N-body code that uses the full multigrid algorithm for solving the modified Poisson equation of the Quasi Linear formulation of Modified Newtonian Dynamics (QUMOND). A novel approach for handling the boundary conditions using a refinement strategy is implemented and the accuracy of the code is compared with analytical solutions of Kuzmin disks. We then employ the code to compute the predicted rotation curves for a sample of five spiral galaxies from the THINGS sample. We generated static N-body realisations of the galaxies according to their stellar and gaseous surface densities and allowed their distances, mass-to-light ratios (M/L) and both the stellar and gas scale-heights to vary in order to estimate the best fit parameters. We found that NGC 3621, NGC 3521 and DDO 154 are well fit by MOND using expected values of the distance and M/L. NGC 2403 required a moderately larger $M/L$ than expected and NGC 2903 required a substantially larger value. The surprising result was that the scale-height of the dominant baryonic component was well constrained by the rotation curves: the gas scale-height for DDO 154 and the stellar scale-height for the others. In fact, if the suggested stellar scale-height (one-fifth the stellar scale-length) was used in the case of NGC 3621 and NGC 3521 it would not be possible to produce a good fit to the inner rotation curve. For each of the four stellar dominated galaxies, we calculated the vertical velocity dispersions which we found to be, on the whole, quite typical compared with observed stellar vertical velocity dispersions of face on spirals. We conclude that modelling the gas scale-heights of the gas rich dwarf spiral galaxies will be vital in order to make precise conclusions about MOND.

The Baryonic Tully-Fisher Relation (BTFR) is an empirical relation between baryonic mass and rotation velocity in disk galaxies. It provides tests of galaxy formation models in LCDM and of alternative theories like MOND. Observations of gas rich galaxies provide a measure of the slope and normalization of the BTFR that is more accurate (if less precise) than that provided by star dominated spirals, as their masses are insensitive to the details of stellar population modeling. Recent independent data for such galaxies are consistent with Mb = AVf^4 with A = 47+/-6 Msun (km/s)^-4. This is equivalent to MOND with $a_0 = 1.3 +/- 0.3 A/s/s. The scatter in the data is consistent with being due entirely to observational uncertainties. It is unclear why the physics of galaxy formation in LCDM happens to pick out the relation predicted by MOND. We introduce a feedback efficacy parameter E to relate halo properties to those of the galaxies they host. E correlates with star formation rate and gas fraction in the sense that galaxies that have experienced the least star formation have been most impacted by feedback.

This paper presents optical and H-alpha imaging for a large sample of LSB galaxies selected from the PSS-II catalogs (Schombert et. al 1992). As noted in previous work, LSB galaxies span a range of luminosities (-10 > M_V > -20) and sizes (0.3 kpc < R_V25 < 10 kpc), although they are consistent in their irregular morphology. Their H-alpha luminosities (L(H-alpha) range from 10^36 to 10^41 ergs s^-1 (corresponding to a range in star formation, using canonical prescriptions, from 10^-5 to 1 M_solar yr^-1). Although their optical colors are at the extreme blue edge for galaxies, they are similar to the colors of dwarf galaxies (van Zee 2001) and gas-rich irregulars (Hunter & Elmegreen 2006). However, their star formation rates per unit stellar mass are a factor of ten less than other galaxies of the same baryonic mass, indicating that they are not simply quiescent versions of more active star forming galaxies. This paper presents the data, reduction techniques and new philosophy of data storage and presentation. Later papers in this series will explore the stellar population and star formation history of LSB galaxies using this dataset.

In recent postings, Foreman & Scott and Gnedin criticize my work on the baryonic Tully-Fisher relation (BTFR) of gas rich galaxies as tests of MOND and LCDM (McGaugh 2011a,b). These criticisms are rather redundant, as they mostly rehash material I have already discussed. Gnedin is concerned with explaining away the problem of the apparently missing baryons with ionized gas. This is one hypothetical possibility that I have previously discussed (McGaugh et al. 2010; McGaugh & Wolf 2010; McGaugh 2011b). Foreman & Scott claim to find a difference between the acceleration scale that they measure and that I measure from the same data, but the apparent difference results simply from their failure to account for the well known fact (Binney & Tremaine 1987) that flattened disks rotate faster than the equivalent spherical mass distribution (as spelled out for this application in McGaugh 2011b). They further argue that the the intrinsic scatter in the BTFR is merely tiny rather than zero. This difference is within the uncertainty of the uncertainties, so confirms rather than refutes my result. Worse, these papers miss the basic point. Why does MOND have any predictions come true? Why does it provide such an effective and economical description of galaxy dynamics? They do nothing to address these deeper questions. One can imagine that the apparently MONDian behavior of galaxies is merely an approximation to some phenomenology that emerges in the context of LCDM, but we have yet to demonstrate how this occurs. Until we confront this issue seriously, and confirm the existence of non-baryonic cold dark matter in the laboratory, it is reasonable to remain skeptical of LCDM as well as of MOND.

The current cosmological paradigm, LCDM, requires that the mass-energy of the universe be dominated by invisible components: dark matter and dark energy. An alternative to these dark components is that the law of gravity be modified on the relevant scales. A test of these ideas is provided by the Baryonic Tully-Fisher Relation (BTFR), an empirical relation between the observed mass of a galaxy and its rotation velocity. Here I report a test using gas rich galaxies for which both axes of the BTFR can be measured independently of the theories being tested and without the systematic uncertainty in stellar mass that affects the same test with star dominated spirals. The data fall precisely where predicted a priori by the modified Newtonian dynamics (MOND). The scatter in the BTFR is attributable entirely to observational uncertainty. This is consistent with the action of a single effective force law but poses a serious fine-tuning problem for LCDM.

Standard cosmology has many successes on large scales, but faces some fundamental difficulties on small, galactic scales. One such difficulty is the cusp/core problem. High resolution observations of the rotation curves for dark matter dominated low surface brightness (LSB) galaxies imply that galactic dark matter halos have a density profile with a flat central core, whereas N-body structure formation simulations predict a divergent (cuspy) density profile at the center. It has been proposed that this problem can be resolved by stellar feedback driving turbulent gas motion that erases the initial cusp. However, strong gravitational lensing prefers a cuspy density profile for galactic halos. In this paper, we use the most recent high resolution observations of the rotation curves of LSB galaxies to fit the core size as a function of halo mass, and compare the resultant lensing probability to the observational results for the well defined combined sample of the Cosmic Lens All-Sky Survey (CLASS) and Jodrell Bank/Very Large Array Astrometric Survey (JVAS). The lensing probabilities based on such density profiles are too low to match the observed lensing in CLASS/JVAS. High baryon densities in the galaxies that dominate the lensing statistics can reconcile this discrepancy, but only if they steepen the mass profile rather than making it more shallow. This places contradictory demands upon the effects of baryons on the central mass profiles of galaxies.

Local Group dwarf spheroidal satellite galaxies are the faintest extragalactic stellar systems known. We examine recent data for these objects in the plane of the Baryonic Tully-Fisher Relation (BTFR). While some dwarf spheroidals adhere to the BTFR, others deviate substantially. We examine the residuals from the BTFR and find that they are not random. The residuals correlate with luminosity, size, metallicity, ellipticity, and susceptibility of the dwarfs to tidal disruption in the sense that fainter, more elliptical, and tidally more susceptible dwarfs deviate farther from the BTFR. These correlations disfavor stochastic processes and suggest a role for tidal effects. We identify a test to distinguish between the {\Lambda}CDM and MOND based on the orbits of the dwarf satellites of the Milky Way and how stars are lost from them.

Dwarf and low surface brightness galaxies are ideal objects to test modified Newtonian dynamics (MOND), because in most of these galaxies the accelerations fall below the threshold below where MOND supposedly applies. We have selected from the literature a sample of 27 dwarf and low surface brightness galaxies. MOND is successful in explaining the general shape of the observed rotation curves for roughly three quarters of the galaxies in the sample presented here. However, for the remaining quarter, MOND does not adequately explain the observed rotation curves. Considering the uncertainties in distances and inclinations for the galaxies in our sample, a small fraction of poor MOND predictions is expected and is not necessarily a problem for MOND. We have also made fits taking the MOND acceleration constant, a_0, as a free parameter in order to identify any systematic trends. We find that there appears to be a correlation between central surface brightness and the best-fit value of a_0, in the sense that lower surface brightness galaxies tend to have lower a_0. However, this correlation depends strongly on a small number of galaxies whose rotation curves might be uncertain due to either bars or warps. Without these galaxies, there is less evidence of a trend, but the average value we find for a_0 ~ 0.7*10^-8 cm s^-2 is somewhat lower than derived from previous studies. Such lower fitted values of a_0 could occur if external gravitational fields are important.

We consider dark masses measured from kinematic tracers at discrete radii in galaxies for which baryonic contributions to overall potentials are either subtracted or negligible. Recent work indicates that rotation curves due to dark matter (DM) halos at intermediate radii in spiral galaxies are remarkably similar, with a mean rotation curve given by $\log_{10}[V_{c,\mathrm{DM}}/(\mathrm{km s^{-1}})]=1.47_{-0.19}^{+0.15}+0.5\log_{10}[r/\mathrm{kpc}]$. Independent studies show that while estimates of the dark mass of a given dwarf spheroidal (dSph) galaxy are robust only near the half-light radius, data from the Milky Way's (MW's) dSph satellites are consistent with a narrow range of mass profiles. Here we combine published constraints on the dark halo masses of spirals and dSphs and include available measurements of low surface brightness galaxies for additional comparison. We find that most measured MW dSphs lie on the extrapolation of the mean rotation curve due to DM in spirals. The union of MW-dSph and spiral data appears to follow a mass-radius relation of the form $M_{\mathrm{DM}}(r)/M_{\odot}=200_{-120}^{+200}(r/\mathrm{pc})^2$, or equivalently a constant acceleration $g_{\mathrm{DM}}=3_{-2}^{+3}\times 10^{-9}\mathrm{cm s^{-2}}$, spanning $0.02\la r \la 75$ kpc. Evaluation at specific radii immediately generates two results from the recent literature: a common mass for MW dSphs at fixed radius and a constant DM central surface density for galaxies ranging from MW dSphs to spirals. However, recent kinematic measurements indicate that M31's dSph satellites are systematically less massive than MW dSphs of similar size. Such deviations from what is otherwise a surprisingly uniform halo relation presumably hold clues to individual formation and evolutionary histories.

Galaxies are composed of baryonic stars and gas embedded in dark matter halos. Here I briefly review two aspects of the connection between baryons and their halos. (1) The observed baryon content of galaxies falls short of the cosmic baryon fraction by an amount that varies systematically with mass. Where these missing baryons now reside is unclear. (2) The characteristic acceleration in disk galaxies correlates strongly with their baryonic mass surface density. This implies a close coupling between the gravitational dynamics, which is presumably dominated by dark matter, and the purely baryonic components of galaxies.

We make an inventory of the baryonic and gravitating mass in structures ranging from the smallest galaxies to rich clusters of galaxies. We find that the fraction of baryons converted to stars reaches a maximum between M500 = 1E12 and 1E13 Msun, suggesting that star formation is most efficient in bright galaxies in groups. The fraction of baryons detected in all forms deviates monotonically from the cosmic baryon fraction as a function of mass. On the largest scales of clusters, most of the expected baryons are detected, while in the smallest dwarf galaxies, fewer than 1% are detected. Where these missing baryons reside is unclear.

We calibrate the Baryonic Tully-Fisher (BTF) Relation using a sample of gas dominated galaxies. These determine the absolute scale of the baryonic mass--rotation speed relation independent of the choice of stellar mass estimator. We find a BTF slope of 3.94 +/- 0.07 (random) +/- 0.08 (systematic) and a zero point of 1.79 +/- 0.26 (random) +/- 0.25 (systematic). We apply this relation to estimate the stellar masses of star dominated galaxies. This procedure reproduces the trend of mass-to-light ratio with color predicted by population synthesis models. The normalization is also correct, consistent with empirical estimates of the IMF used in such models.

We present optical long-slit and SparsePak Integral Field Unit emission line spectroscopy along with optical broadband and near IR images of the edge-on spiral galaxy NGC 2683. We find a multi-valued, figure-of-eight velocity structure in the inner 45 arcsec of the long-slit spectrum and twisted isovelocity contours in the velocity field. We also find, regardless of wavelength, that the galaxy isophotes are boxy. We argue that taken together, these kinematic and photometric features are evidence for the presence of a bar in NGC 2683. We use our data to constrain the orientation and strength of the bar.

The baryonic Tully-Fisher relation (BTF) is a fundamental relation between baryonic mass and maximum rotation velocity. It can be used to estimate distances, as well as to constrain the properties of dark matter and its relation with the visible matter. In this paper, we explore if extremely low-mass dwarf galaxies follow the same BTF relation as high-mass galaxies. We quantify the scatter in the BTF relation and use this to constrain the allowed elongations of dark matter halo potentials. We obtained HI synthesis data of 11 dwarf galaxies and derive several independent estimates for the maximum rotation velocity. Constructing a BTF relation using data from the literature for the high-mass end, and galaxies with detected rotation from our sample for the low-mass end results in a BTF with a scatter of 0.33 mag. This scatter constrains the ellipticities of the potentials in the plane of the disks of the galaxies to an upper limit of 0-0.06 indicating that dwarf galaxies are at most only mildly tri-axial. Our results indicate that the BTF relation is a fundamental relation which all rotationally dominated galaxies seem to follow.

We present Chandra detections of x-ray emission from the AGN in two giant Low Surface Brightness (LSB) galaxies, UGC 2936 and UGC 1455. Their x-ray luminosities are 1.8\times10^{42} ergs/s and 1.1\times10^{40} ergs/s respectively. Of the two galaxies, UGC 2936 is radio loud. Together with another LSB galaxy UGC 6614 (XMM archival data) both appear to lie above the X-ray-Radio fundamental plane and their AGN have black hole masses that are low compared to similar galaxies lying on the correlation. However, the bulges in these galaxies are well developed and we detect diffuse x-ray emission from four of the eight galaxies in our sample. Our results suggest that the bulges of giant LSB galaxies evolve independently of their halo dominated disks which are low in star formation and disk dynamics. The centers follow an evolutionary path similar to that of bulge dominated normal galaxies on the Hubble Sequence but the LSB disks remain unevolved. Thus the bulge and disk evolution are decoupled and so whatever star formation processes produced the bulges did not affect the disks.

We model the NFW potential to determine if, and under what conditions, the NFW halo appears consistent with the observed velocity fields of low surface brightness (LSB) galaxies. We present mock DensePak IFU velocity fields and rotation curves of axisymmetric and non-axisymmetric potentials that are well-matched to the spatial resolution and velocity range of our sample galaxies. We find that the DensePak IFU can accurately reconstruct the velocity field produced by an axisymmetric NFW potential and that a tilted-ring fitting program can successfully recover the corresponding NFW rotation curve. We also find that non-axisymmetric potentials with fixed axis ratios change only the normalization of the mock velocity fields and rotation curves and not their shape. The shape of the modeled NFW rotation curves does not reproduce the data: these potentials are unable to simultaneously bring the mock data at both small and large radii into agreement with observations. Indeed, to match the slow rise of LSB galaxy rotation curves, a specific viewing angle of the non-axisymmetric potential is required. For each of the simulated LSB galaxies, the observer's line-of-sight must be along the minor axis of the potential, an arrangement which is inconsistent with a random distribution of halo orientations on the sky.

Aims:We have studied the bulge and the disk kinematics of the giant low surface brightness galaxy ESO 323-G064 in order to investigate its dynamical properties and the radial mass profile of the dark matter (DM) halo. Methods:We observed the galaxy with integral field spectroscopy (VLT/VIMOS, in IFU configuration), measured the positions of the ionized gas by fitting Gaussian functions to the O[III] and Hbeta emission lines, and fit stellar templates to the galaxy spectra to determine velocity and velocity dispersions. We modeled the stellar kinematics in the bulge with spherical isotropic Jeans models and explored the implications of self consistent and dark matter scenarios for NFW and pseudo isothermal halos. Results:In the bulge-dominated region, r<5", the emission lines show multi-peaked profiles. The disk dominated region of the galaxy, 13"<r<30", exhibits regular rotation, with a flat rotation curve that reaches 248 +/- 6 km/sec. From this we estimate the total barionic mass to be M_bar ~ 1.9 10^11 M_sun and the total DM halo mass to be M_DM ~ 4.8 10^12 M_sun. The stellar velocity and velocity dispersion have been measured only in the innermost ~5" of the bulge, and reveal a regular rotation with an observed amplitude of 140 km/sec and a central dispersion of sigma=180 km/sec. Our simple Jeans modeling shows that dark matter is needed in the central 5" to explain the kinematics of the bulge, for which we estimate a mass of (7 +/- 3) 10^10 M_sun. However, we are not able to disentangle different DM scenarios. The computed central mass density of the bulge of ESO 323-G064 resembles the central mass density of some high surface brightness galaxies, rather than that of low surface brightness galaxies.

Using the Tuorla-Heidelberg model for the mass distribution of the Milky Way, I determine the rotation curve predicted by MOND. The result is in good agreement with the observed terminal velocities interior to the solar radius and with estimates of the Galaxy's rotation curve exterior thereto. There are no fit parameters: given the mass distribution, MOND provides a good match to the rotation curve. The Tuorla-Heidelberg model does allow for a variety of exponential scale lengths; MOND prefers short scale lengths in the range 2.0 to 2.5 kpc. The favored value of scale length depends somewhat on the choice of interpolation function. There is some preference for the `simple' interpolation function as found by Famaey & Binney. I introduce an interpolation function that shares the advantages of the simple function on galaxy scales while having a much smaller impact in the solar system. I also solve the inverse problem, inferring the surface mass density distribution of the Milky Way from the terminal velocities. The result is a Galaxy with `bumps and wiggles' in both its luminosity profile and rotation curve that are reminiscent of those frequently observed in external galaxies.

We present high-resolution optical velocity fields from DensePak integral field spectroscopy, along with derived rotation curves, for a sample of low surface brightness galaxies. In the limit of no baryons, we fit the NFW and pseudoisothermal halo models to the data and find the rotation curve shapes and halo central densities to be better described by the isothermal halo. For those galaxies with photometry, we present halo fits for three assumptions of the stellar mass-to-light ratio. We find that the velocity contribution from the baryons is significant enough in the maximum disk case that maximum disk and the NFW halo are mutually exclusive. We find a substantial cusp mass excess at the centers of the galaxies, with at least two times more mass expected in the cuspy CDM halo than is allowed by the data. We also find that to reconcile the data with LCDM, ~20 km/s noncircular motions are needed and/or the power spectrum has a lower amplitude on the scales we probe.

We studied the bulge and the disk kinematic of the giant Low Surface Brightness (LSB) galaxy ESO 323-G064 in order to investigate its dynamics and its Dark Matter (DM) content. We observed the galaxy with the integral field spectroscopy (VLT/VIMOS, in IFU configuration). Results for the gaseous kinematics (bulge and disk) and stellar kinematics (bulge) are presented, together with a Jeans model for the stellar bulge kinematics.

We consider the orbit of the bullet cluster 1E 0657-56 in both CDM and MOND using accurate mass models appropriate to each case in order to ascertain the maximum plausible collision velocity. Impact velocities consistent with the shock velocity (~ 4700km/s) occur naturally in MOND. CDM can generate collision velocities of at most ~ 3800km/s, and is only consistent with the data provided that the shock velocity has been substantially enhanced by hydrodynamical effects.

The global missing baryon problem - that the sum of observed baryons falls short of the number expected form BBN - is well known. In addition to this, there is also a local missing baryon problem that applies to individual dark matter halos. This halo by halo missing baryon problem is such that the observed mass fraction of baryons in individual galaxies falls short of the cosmic baryon fraction. This deficit is a strong function of circular velocity. I give an empirical estimate of this function, and note the presence of a critical scale of ~900 km/s therein. I also briefly review Omega_b from BBN, highlighting the persistent tension between lithium and the CMB, and discuss some pros and cons of individual galaxies and clusters of galaxies as potential reservoirs of dark baryons.

We present preliminary results of a study of the low frequency radio continuum emission from the nuclei of Giant Low Surface Brightness (LSB) galaxies. We have mapped the emission and searched for extended features such as radio lobes/jets associated with AGN activity. LSB galaxies are poor in star formation and generally less evolved compared to nearby bright spirals. This paper presents low frequency observations of 3 galaxies; PGC 045080 at 1.4 GHz, 610 MHz, 325MHz, UGC 1922 at 610 MHz and UGC 6614 at 610 MHz. The observations were done with the GMRT. Radio cores as well as extended structures were detected and mapped in all three galaxies; the extended emission may be assocated with jets/lobes associated with AGN activity. Our results indicate that although these galaxies are optically dim, their nuclei can host AGN that are bright in the radio domain.

We present radio observations and optical spectroscopy of the giant low surface brightness (LSB) galaxy PGC 045080 (or 1300+0144). PGC 045080 is a moderately distant galaxy having a highly inclined optical disk and massive HI gas content. Radio continuum observations of the galaxy were carried out at 320 MHz, 610 MHz and 1.4 GHz. Continuum emission was detected and mapped in the galaxy. The emission appears extended over the inner disk at all three frequencies. At 1.4 GHz and 610 MHz it appears to have two distinct lobes. We also did optical spectroscopy of the galaxy nucleus; the spectrum did not show any strong emission lines associated with AGN activity but the presence of a weak AGN cannot be ruled out. Furthermore, comparison of the H$\alpha$ flux and radio continuum at 1.4 GHz suggests that a significant fraction of the emission is non-thermal in nature. Hence we conclude that a weak or hidden AGN may be present in PGC 045080. The extended radio emission represents lobes/jets from the AGN. These observations show that although LSB galaxies are metal poor and have very little star formation, their centers can host significant AGN activity. We also mapped the HI gas disk and velocity field in PGC 045080. The HI disk extends well beyond the optical disk and appears warped. In the HI intensity maps, the disk appears distinctly lopsided. The velocity field is disturbed on the lopsided side of the disk but is fairly uniform in the other half. We derived the HI rotation curve for the galaxy from the velocity field. The rotation curve has a flat rotation speed of ~ 190 km/s.

We examine the amplitude of the rotation velocity that can be attributed to the dark matter halos of disk galaxies, focusing on well measured intermediate radii. The data for 60 galaxies spanning a large range of mass and Hubble types, taken together, are consistent with a dark halo velocity log(Vh) = C + B log(r) with C = 1.47 (+0.15, -0.19) and B = 1/2 over the range 1 < r < 74 kpc. The range in C stems from different choices of the stellar mass estimator, from minimum to maximum disk. For all plausible choices of stellar mass, the implied densities of the dark halos are lower than expected from structure formation simulations in LCDM, which anticipate C > 1.6. This problem is not specific to a particular type of galaxy or to the innermost region of the halo (cusp or core); the velocity attributable to dark matter is too low at all radii.

LSB galaxies have low metallicities, diffuse stellar disks, and massive HI disks. We have detected molecular gas in two giant LSB galaxies, UGC 6614 and F568-6. A millimeter continuum source has been detected in UGC 6614 as well. At centimeter wavelengths we have detected and mapped the continuum emission from the giant LSB galaxy 1300+0144. The emission is extended about the nucleus and is most likely originating from the AGN in the galaxy. The HI gas distribution and velocity field in 1300+0144 was also mapped. The HI disk extends well beyond the optical disk and appears lopsided in the intensity maps.

We present BIMA and IRAM CO(1--0) observations of seven low surface brightness (LSB) galaxies, including three large spirals with faint disks but prominent bulges, and four relatively small LSB galaxies with irregular disks. The giant LSB galaxies are UGC 5709, UGC 6614 and F568-6 (Malin2). The smaller LSB galaxies are NGC 5585, UGC 4115, UGC 5209 and F583-1. The galaxies were selected based on their relatively high metallicity and apparent signs of star formation in their disks. The BIMA maps suggested the presence of molecular gas in 2 of the giant LSB galaxies, F568-6 and UGC 6614. Using the 30m IRAM telescope we detected CO (1--0) emission in the disks of both galaxies and in the nucleus of F568-6. The molecular gas in these galaxies is clearly offset from the nucleus and definitely associated with the LSB disk. In addition we also detected a millimeter continuum source in the center of UGC 6614. When compared with VLA 1.5 GHz observations of the galaxy, the emission was found to have a flat spectrum indicating that the millimeter continuum emission is most likely due to an active galactic nucleus (AGN) in the galaxy. Our results show that giant LSB spirals may contain significant quantities of molecular gas in their disks and also harbor radio bright AGN in their centers.

I review the data relating to the appearance of the missing mass problem at a particular acceleration scale. Rotation curves are examined in detail, with emphasis on the empirical connection between baryonic and total mass distributions. Other lines of evidence are also considered, including disk stability, dwarf spheroidals, giant ellipticals, clusters of galaxies, the microwave background, BBN, and the growth of large scale structure.

We present high resolution two-dimensional velocity fields from integral field spectroscopy along with derived rotation curves for eleven low surface brightness galaxies. We fit NFW and pseudo-isothermal halo models to the new data combined with previous long-slit and HI data. In most cases we find the pseudo-isothermal halo to better represent the data than the NFW halo, as the NFW concentrations are often lower than expected for a LambdaCDM cosmology. We also compare our results to previous studies and find that including the new two-dimensional optical data does not significantly alter the halo parameters, but does decrease the uncertainties by roughly a factor of 2.

The rotation curves of spiral galaxies obey strong scaling relations. These include the Tully-Fisher and baryonic Tully-Fisher relations, and the mass discrepancy--acceleration relation. These relations can be used to place constraints on the mass-to-light ratios of stars. Once the stellar mass is constrained, the distribution of dark matter follows. The shape of the dark matter distribution is consistent with the expectations of NFW halos exterior to 1 kpc, but the amplitude is wrong. This is presumably related to the long-standing problem of the normalization of the Tully-Fisher relation and may imply a downturn in the amplitude of the power spectrum at small scales. More fundamentally, the persistent success of MOND remains a troubling fact.

We present high resolution two dimensional velocity fields from integral field spectroscopy along with derived rotation curves for nine low surface brightness galaxies. This is a positive step forward in terms of both data quality and number of objects studied. We fit NFW and pseudo-isothermal halo models to the observations. We find that the pseudo-isothermal halo better represents the data in most cases than the NFW halo, as the resulting concentrations are lower than would be expected for LCDM.

A fine balance between dark and baryonic mass is observed in spiral galaxies. As the contribution of the baryons to the total rotation velocity increases, the contribution of the dark matter decreases by a compensating amount. This poses a fine-tuning problem for \LCDM galaxy formation models, and may point to new physics for dark matter particles or even a modification of gravity.

I investigate the Baryonic Tully-Fisher relation for a sample of galaxies with extended 21 cm rotation curves spanning the range 20 < Vf < 300 km/s. A variety of scalings of the stellar mass-to-light ratio are considered. For each prescription for M/L, I give fits of the form Md = A Vf^x. Presumably, the prescription that comes closest to the correct value will minimize the scatter in the relation. The fit with minimum scatter has A = 50 solar masses/(km/s)^4 and x = 4. This relation holds over five decades in mass. Galaxy color, stellar fraction, and M/L are correlated with each other and with Md, in the sense that more massive galaxies tend to be more evolved. There is a systematic dependence of the degree of maximality of disks on surface brightness. High surface brightness galaxies typically have M/L ~ 3/4 of the maximum disk value, while low surface brightness galaxies typically attain ~ 1/4 of this amount.

The initial radial density profiles of dark matter halos are laid down by gravitational collapse in hierarchical structure formation scenarios and are subject to further compression as baryons cool and settle to the halo centers. We here describe an explicit implementation of the algorithm, originally developed by Young, to calculate changes to the density profile as the result of adiabatic infall in a spherical halo model. Halos with random motion are more resistant to compression than are those in which random motions are neglected, which is a key weakness of the simple method widely employed. Young's algorithm results in density profiles in excellent agreement with those from N-body simulations. We show how the algorithm may be applied to determine the original uncompressed halos of real galaxies, a step which must be computed with care in order to enable a confrontation with theoretical predictions from theories such as LCDM.

We validate the baryonic Tully Fisher (BTF) relation by exploring the Tully Fish er (TF) and BTF properties of optically and HI-selected disk galaxies. The data includes galaxies from: Sakai et al. (2000) calibrator sample; McGaugh et al. (2000: MC2000) I-band sample; and 18 newly acquired HI-selected field dwarf galaxies observed with the ANU 2.3m telescope and the ATNF Parkes telescope from Gurovich's thesis sample (2005). As in MC2000, we re-cast the TF and BTF relations as relationships between baryo n mass and W_{20}. First we report some numerical errors in MC2000. Then, we c alculate weighted bi-variate linear fits to the data, and finally we compare the fits of the intrinsically fainter dwarfs with the brighter galaxies of Sakai et al. (2000). With regards to the local calibrator disk galaxies of Sakai et al. (2000), our results suggest that the BTF relation is indeed tighter than the T F relation and that the slopes of the BTF relations are statistically flatter th an the equivalent TF relations. Further, for the fainter galaxies which include the I-band MCG2000 and HI-selected galaxies of Gurovich's thesis sample, we calc ulate a break from a simple power law model because of what appears to be real c osmic scatter. Not withstanding this point, the BTF models are marginally better models than the equivalent TF ones with slightly smaller reduced chi^2.

I present a model devoid of non-baryonic cold dark matter (CDM) which provides an acceptable fit to the WMAP data for the power spectrum of temperature fluctuations in the cosmic background radiation (CBR). An a priori prediction of such no-CDM models was a first-to-second peak amplitude ratio A1:2 = 2.4. WMAP measures A1:2 = 2.34 +/- 0.09. The baryon content is the dominant factor in fixing this ratio; no-CDM models which are consistent with the WMAP data are also consistent with constraints on the baryon density from the primordial abundances of 2H, 4He, and 7Li. However, in order to match the modest width of the acoustic peaks observed by WMAP, a substantial neutrino mass is implied: m(nu) ~ 1 eV. Even with such a heavy neutrino, structure is expected to form rapidly under the influence of MOND. Consequently, the epoch of reionization should occur earlier than is nominally expected in LCDM. This prediction is realized in the polarization signal measured by WMAP. An outstanding test is in the amplitude of the third acoustic peak. Experiments which probe high-L appear to favor a third peak which is larger than predicted by the no-CDM model.

The mass discrepancy in disk galaxies is shown to be well correlated with acceleration, increasing systematically with decreasing acceleration below a critical scale a0 = 3700 km^2/s^2/kpc = 1.2E-10 m/s/s. For each galaxy, there is an optimal choice of stellar mass-to-light ratio which minimizes the scatter in this mass discrepancy-acceleration relation. The same mass-to-light ratios also minimize the scatter in the baryonic Tully-Fisher relation and are in excellent agreement with the expectations of stellar population synthesis. Once the disk mass is determined in this fashion, the dark matter distribution is specified. The circular velocity attributable to the dark matter can be expressed as a simple equation which depends only on the observed distribution of baryonic mass. It is a challenge to understand how this very fine-tuned coupling between mass and light comes about.

We report the oxygen abundances of the HII regions of a sample of low surface brightness (LSB) galaxies. We provide analytic functions describing the McGaugh (1991) calibration of the R_23 method. We use this and the equivalent width (EW) method to determine oxygen abundances, and also make direct estimates in a few cases where the temperature sensitive [OIII] 4363 line is available. We find LSB galaxies to be metal poor, consistent with the L-Z relation of other galaxies. The large gas mass fractions of these objects provide an interesting test of chemical evolution models. We find no obvious deviation from the closed-box model of galactic chemical evolution. Based on our abundance and gas mass fraction measurements, we infer that LSB galaxies are not fundamentally different than other galaxy types but are perhaps at an early stage of evolution.

Low Surface Brightness (LSB) galaxies are dominated by dark matter, and their rotation curves thus reflect their dark matter distribution. Recent high-resolution rotation curves suggest that their dark matter mass-density distributions are dominated by a constant-density core. This seems inconsistent with the predictions of Cold Dark Matter (CDM) models which produce halos with compact density cusps and steep mass-density profiles. However, the observationally determined mass profiles may be affected by non-circular motions, asymmetries and offsets between optical and dynamical centres, all of which tend to lower the observed slopes. Here we determine the impact of each of these effects on a variety of halo models, and compare the results with observed mass-density profiles. Our simulations suggest that no single systematic effect can reconcile the data with the cuspy CDM halos. The data are best described by a model with a soft core with an inner power-law mass-density slope alpha = -0.2 +/- 0.2. However, no single universal halo profile provides a completely adequate description of the data.

The concentrations of the cuspy dark matter halos predicted by simulations of cold dark matter are related to the cosmology in which the halos form. Observational constraints on halo concentration therefore map into constraints on cosmological parameters. In order to explain the observed concentrations of dark matter dominated low surface brightness galaxies, we require a cosmology with rather little power on galaxy scales. Formally, we require s_8 G_0.6 < 0.23, where G_0.6 is a modified shape parameter appropriate to this problem. Practically, this means either Omega_m < 0.2 or s_8 < 0.8. These limits apply as long as we insist that the cuspy halos found in simulations must describe the halos of low surface brightness galaxies. A low density cosmology helps with the low observed concentrations, but it offers no explanation of the many cases where the shape of the density profile deviates from the predicted cuspy form. These cases must have suffered very extensive mass redistribution if the current halo formation picture is not to fail outright. It is far from clear whether any of the mass redistribution mechanisms which have been suggested (e.g., feedback) are viable.

Modified Newtonian dynamics (MOND) is an empirically motivated modification of Newtonian gravity or inertia suggested by Milgrom as an alternative to cosmic dark matter. The basic idea is that at accelerations below a0 ~ 10^{-8} cm/s^2 ~ cH0/6 the effective gravitational attraction approaches sqrt{gN*a0} where gN is the usual Newtonian acceleration. This simple algorithm yields flat rotation curves for spiral galaxies and a mass-rotation velocity relation of the form M ~ V^4 that forms the basis for the observed luminosity-rotation velocity relation-- the Tully-Fisher law. We review the phenomenological success of MOND on scales ranging from dwarf spheroidal galaxies to superclusters, and demonstrate that the evidence for dark matter can be equally well interpreted as evidence for MOND. We discuss the possible physical basis for an acceleration-based modification of Newtonian dynamics as well as the extension of MOND to cosmology and structure formation.

I review what we currently do and do not know about the masses of disk galaxies and their dark matter halos. The prognosis for disks is good: the asymptotic rotation velocity provides a good indicator of total disk mass. The prognosis for halos is bad: cuspy halos provide a poor description of the data, and the total mass of individual dark matter halos remains ill-constrained.

Stars have mass; halos lack cusps.

We present mass models for a sample of 30 high-resolution rotation curves of low surface brightness (LSB) galaxies. We fit both pseudo-isothermal (core-dominated) and Cold Dark Matter (CDM) (cusp-dominated) halos for a wide variety of assumptions about the stellar mass-to-light ratio. We find that the pseudo-isothermal model provides superior fits. CDM fits show systematic deviations from the data, and often have a small statistical likelihood of being the appropriate model. The distribution of concentration parameters is too broad, and has too low a mean, to be explained by $\Lambda$CDM. This failing becomes more severe as increasing allowance is made for stellar mass: NFW fits require uncomfortably low mass-to-light ratios. In contrast, the maximum disk procedure does often succeed in predicting the inner shape of the rotation curves, but requires uncomfortably large stellar mass-to- light ratios. The data do admit reasonable stellar population mass-to-light ratios if halos have cores rather than cusps.

We present long slit Halpha observations of 50 low surface brightness galaxies. Of these, 36 are of sufficient quality to form rotation curves. These data provide a large increase in the number of low surface brightness galaxies for which accurate rotation curves are available. They also represent an order of magnitude improvement in spatial resolution over previous 21 cm studies (1" to 2" instead of 13" to 45"). The improved resolution and accuracy of the data extend and strengthen the scientific conclusions previously inferred from 21 cm data.

We derive the mass density profiles of dark matter halos that are implied by high spatial resolution rotation curves of low surface brightness galaxies. We find that at small radii, the mass density distribution is dominated by a nearly constant density core with a core radius of a few kpc. For rho(r) ~ r^a, the distribution of inner slopes a is strongly peaked around a = -0.2. This is significantly shallower than the cuspy a < -1 halos found in CDM simulations. While the observed distribution of alpha does have a tail towards such extreme values, the derived value of alpha is found to depend on the spatial resolution of the rotation curves: a ~ -1 is found only for the least well resolved galaxies. Even for these galaxies, our data are also consistent with constant density cores (a = 0) of modest (~ 1 kpc) core radius, which can give the illusion of steep cusps when insufficiently resolved. Consequently, there is no clear evidence for a cuspy halo in any of the low surface brightness galaxies observed.

The optical and HI properties for a sample of low surface brightness (LSB) dwarf galaxies, cataloged from the Second Palomar Sky Survey, is presented. Gas mass fractions for LSB dwarfs reach the highest levels of any know galaxy type (f_g=95%) confirming that their low stellar densities are due to inefficient conversion of gas mass into stellar mass. Comparison with star formation models indicates that the blue optical colors of LSB dwarfs is not due to low metallicity or recent star formation and can only be explained by a dominant stellar population that is less than 5 Gyrs in mean age. If star formation occurs in OB complexes, similar to normal galaxies, then LSB dwarfs must undergo weak bursts traveling over the extent of the galaxy to maintain their LSB nature, which contributes to their irregular morphological appearance.

The amplitude of the second peak in the angular power spectrum of the cosmic microwave background radiation is constrained to be small by recent experiments like Boomerang. This is surprising in the context of the LambdaCDM model, which predicted a large second peak. However, this result is expected if CDM does not exist. The observed shape of the power spectrum was accurately predicted (before the fact) by a model motivated by the surprising recent successes of the modified dynamics (MOND) hypothesized by Milgrom.

I describe the disk mass-rotation velocity relation which underpins the familiar luminosity-linewidth relation. Continuity of this relation favors nearly maximal stellar mass-to-light ratios. This contradicts the low mass-to-light ratios implied by the lack of surface brightness dependence in the same relation.

The amplitudes of peaks in the angular power spectrum of anisotropies in the microwave background radiation depend on the mass content of the universe. The second peak should be prominent when cold dark matter is dominant, but is depressed when baryons dominate. Recent microwave background data are consistent with a purely baryonic universe with Omega(matter) = Omega(baryon) ~ 0.03 and Omega(Lambda) ~ 1.

We explore the Tully-Fisher relation over five decades in stellar mass in galaxies with circular velocities ranging over 30 < Vc < 300 km/s. We find a clear break in the optical Tully-Fisher relation: field galaxies with Vc < 90 km/s fall below the relation defined by brighter galaxies. These faint galaxies are however very gas rich; adding in the gas mass and plotting baryonic disk mass Md = M* + Mg in place of luminosity restores a single linear relation. The Tully-Fisher relation thus appears fundamentally to be a relation between rotation velocity and total baryonic mass of the form Md = A Vc^4.

We present VLA H I observations of the low surface brightness galaxy UGC 12695 and its two companions, UGC 12687 and a newly discovered dwarf galaxy 2333+1234. UGC 12695 shows solid body rotation but has a very lopsided morphology of the H I disk, with the majority of the H I lying in the southern arm of the galaxy. The H I column density distribution of this very blue, LSB galaxy coincides in detail with its light distribution. Comparing the H I column density of UGC 12695 with the empirical (but not well understood) value of Sigma_c = 10E21 atoms/cm^2 found in, i.e., Skillman's 1986 paper shows the star formation to be a local affair, occurring only in those regions where the column density is above this star formation threshold. The low surface brightness nature of this galaxy could thus be attributed to an insufficient gas surface density, inhibiting star formation on a more global scale. Significantly, though, the Toomre criterion places a much lower critical density on the galaxy (+/-10E20 atoms/cm^2), which is shown by the galaxy's low SFR to not be applicable. Within a projected distance of 300kpc/30kms of UGC 12695 lie two companion galaxies - UGC 12687, a high surface brightness barred spiral galaxy, and 2333+1234, a dwarf galaxy discovered during this investigation. The close proximity of the three galaxies, combined with UGC 12695's extremely blue color and regions of localized starburst and UGC 12687's UV excess bring to mind mutually induced star formation through tidal activity.

The rotation curves of disk galaxies exhibit a number of striking regularities. The amplitude of the rotation is correlated with luminosity (Tully-Fisher), the shape of the rotation curve is well predicted by the luminous mass distribution, and the magnitude of the mass discrepancy increases systematically with decreasing centripetal acceleration. These properties indicate a tight connection between light and mass, and impose strong constraints on theories of galaxy formation.

Two hypothesized solutions of the mass discrepancy problem are cold dark matter (CDM) and modified Newtonian dynamics (MOND). The virtues and vices of these very different hypotheses are largely disjoint, making the process of distinguishing between them very dependent on how we weigh disparate lines of evidence. One clear difference is the nature of the principal mass constituent of the universe (CDM or baryons). This difference in the baryon fraction ($f_b \approx 0.1$ vs. 1) should leave a distinctive signature in the spectrum of fluctuations in the cosmic microwave background. Here I discuss some of the signatures which should be detectable in the near future. The most promising appears to be the ratio of the amplitudes of the first two peaks relative to the intervening trough.

I review current theories of disk galaxy formation. All fail in the light of recent data for low surface brightness galaxies. The difficulty stems from the presence of a unique acceleration scale in the data.

I explore some consequences of Milgrom's modified dynamics for cosmology. There appear to be two promising tests for distinguishing MOND from CDM: (1) the rate of growth of structure and (2) the baryon fraction. These should be testable with observations of clusters at high redshift and the microwave background, respectively.

We present models for the interstellar medium in disk galaxies. In particular, we investigate whether the ISM in low surface brightness galaxies can support a significant fraction of molecular gas given their low metallicity and surface density. It is found that the abundance and line brightness of CO in LSB galaxies is small and typically below current observational limits. Still, depending on physical details of the ISM, the fraction of gas in the form of molecular hydrogen can be significant in the inner few kpc of a low surface brightness galaxy. This molecular gas would be at temperatures of $\sim$ 30 -- 50 K, rather higher than in high surface brightness galaxies. These results may help explain the star forming properties and inferred evolutionary history of LSB galaxies.

Our view of the properties of galaxies is strongly affected by the way in which we survey for them. I discuss some aspects of selection effects and methods to compensate for them. One result is an estimate of the surface brightness distribution. I believe this is progress, but considerable uncertainty remains.

Utilizing the F814W and F300W filters, short exposure HST WFPC-2 images were taken of UGC 12695 a nearby (z ~ 0.021) low surface brightness disk galaxy. UGC 12695 has an unusual morphology, consisting of a Y-shaped nucleus surrounded by a faint spiral arm with a number of bright H II regions interspersed throughout the galaxy. Surface photometry indicates the majority of recent star formation in this galaxy occurred in these very localized regions, most of which have a radius of < 2''. This uneven stellar distribution combined with the galaxy's overall extremely blue color and low metallicity indicates UGC 12695 is an unevolved galaxy. Some of the structural peculiarities of this galaxy arise because a number of background galaxies, previously thought to be morphological components of this galaxy, are showing through both the outer nucleus and spiral arms of UGC 12695. Surface photometry of these galaxies shows them to be fairly small (alpha < 1.8'') disk galaxies with total magnitudes ranging from 19.6 through 24.2, and central surface brightnesses between 20.2 mag/arcsec^2 < mu(0) < 23.1 mag/arcsec^2. When possible, the U, B, V, and I colors of these galaxies were measured using ground-based images, which show the galaxies to be fairly red indicatingthey are likely at redshifts z > 0.5.Inclusion of them in the photometry of UGC 12695 makes the galaxyappear significantly redder. When these galaxies are masked out, the resultant U-I color of UGC 12695 is -0.2 +/- 0.1, making it perhaps the bluest galaxy every measured in this color system and confirming its nature as a very unevolved galaxy at low redshift.Spectroscopy of these background galaxies, through the transparent disk of UGC 12695, may help to clarify its chemical evolution and current heavy element content.

We present MOND (Modified Newtonian Dynamics) fits to 15 rotation curves of LSB galaxies. Good fits are readily found, although for a few galaxies minor adjustments to the inclination are needed. Reasonable values for the stellar mass-to-light ratios are found, as well as an approximately constant value for the total (gas and stars) mass-to-light ratio. We show that the LSB galaxies investigated here lie on the one, unique Tully-Fisher relation, as predicted by MOND. The scatter on the Tully-Fisher relation can be completely explained by the observed scatter in the total mass-to-light ratio. We address the question of whether MOND can fit any arbitrary rotation curve by constructing a plausible fake model galaxy. While MOND is unable to fit this hypothetical galaxy, a normal dark halo fit is readily found, showing that dark matter fits are much less selective in producing fits. The good fits to rotation curves of LSB galaxies support MOND, especially as these are galaxies with large mass discrepancies deep in the MOND regime.

The severity of the mass discrepancy in spiral galaxies is strongly correlated with the central surface brightness of the disk. Progressively lower surface brightness galaxies have ever larger mass discrepancies. No other parameter (luminosity, size, velocity, morphology) is so well correlated with the magnitude of the mass deficit. The rotation curves of low surface brightness disks thus provide a unique data set with which to probe the dark matter distribution in galaxies. The mass discrepancy is apparent from $R = 0$ giving a nearly direct map of the halo mass distribution. The luminous mass is insignificant. Interpreting the data in terms of dark matter leads to troublesome fine-tuning problems. Different observations require contradictory amounts of dark matter. Structure formation theories are as yet far from able to explain the observations.

The rotation curves of low surface brightness galaxies provide a unique data set with which to test alternative theories of gravitation over a large dynamic range in size, mass, surface density, and acceleration. Many clearly fail, including any in which the mass discrepancy appears at a particular length-scale. One hypothesis, MOND [Milgrom 1983, ApJ, 270, 371], is consistent with the data. Indeed, it accurately predicts the observed behavior. We find no evidence on any scale which clearly contradicts MOND, and a good deal which supports it.

A number of observations strongly suggest that the baryon fraction is not a universal constant. One obvious interpretation is that there is some distribution of Fb, and the different observations sample different portions of the distribution. However, the small intrinsic scatter in the Tully-Fisher relation requires that the baryon fraction be very nearly universal. It is not easy to resolve this paradox in the framework of the standard picture.

We present mass models of a sample of 19 low surface brightness (LSB) galaxies and compare the properties of their constituent mass components with those of a sample of high surface brightness (HSB) galaxies.We find that LSB galaxies are dark matter dominated. Their halo parameters are only slightly affected by assumptions on stellar mass-to-light ratios. Comparing LSB and HSB galaxies we find that mass models derived using the maximum disk hypothesis result in the disks of LSB galaxies having systematically higher stellar mass-to-light ratios than HSB galaxies of similar rotation velocity. This is inconsistent with all other available evidence on the evolution of LSB galaxies. We argue therefore that the maximum disk hypothesis does not provide a representative description of the LSB galaxies and their evolution. Mass models with stellar mass-to-light ratios determined by the colors and stellar velocity dispersions of galactic disks imply that LSB galaxies have dark matter halos that are more extended and less dense than those of HSB galaxies. Surface brightness is thus related to the halo properties. LSB galaxies are slowly evolving, low density and dark matter dominated galaxies.

In twenty years, low surface brightness (LSB) galaxies have evolved from being an idiosyncratic notion to being one of the major baryonic repositories in the Universe. The story of their discovery and the characterization of their properties is told here. Their recovery from the noise of the night sky background is a strong testament to the severity of surface brightness selection effects. LSB galaxies have a number of remarkable properties which distinguish them from the more familiar Hubble Sequence of spirals. The two most important are 1) they evolve at a significantly slower rate and may well experience star formation outside of the molecular cloud environment, 2) they are embedded in dark matter halos which are of lower density and more extended than the halos around high surface brightness (HSB) disk galaxies. Compared to HSB disks, LSB disks are strongly dark matter dominated at all radii and show a systematic increase in M/L with central surface brightness. In addition, the recognition that large numbers of LSB galaxies actually exist has changed the form of the galaxy luminosity function and has clearly increased the space density of galaxies at z = 0. Recent CCD surveys have uncovered a population of red LSB disks that may be related to the excess of faint blue galaxies detected at moderate redshifts. LSB galaxies offer us a new window into galaxy evolution and formation which is every bit as important as those processes which have produced easy to detect galaxies. Indeed, the apparent youth of some LSB galaxies suggest that galaxy formation is a greatly extended process. While the discovery of LSB galaxies have lead to new insights, it remains unwise to presume that we now have a representative sample which encompasses all galaxy types and forms.

Using analytic stability criteria, we demonstrate that, due to their low surface mass density and large dark matter content, LSB disks are quite stable against the growth of global nonaxisymmetric modes such as bars. However, depending on their (poorly constrained) stellar velocity dispersions, they may be only marginally stable against local instabilities. We simulate a collision between an LSB and HSB galaxy and find that, while the HSB galaxy forms a strong bar, the response of the LSB disk is milder, manifesting weaker rings and spiral features. The lack of sufficient disk self-gravity to amplify dynamical instabilities naturally explains the rarity of bars in LSB disks. The stability of LSB disks may also inhibit interaction-driven gas inflow and starburst activity in these galaxies.

We show that the gas mass fraction of spiral galaxies is strongly correlated with luminosity and surface brightness. It is not correlated with linear size. Gas fraction varies with luminosity and surface brightness at the same rate, indicating evolution at fixed size. Dim galaxies are clearly less evolved than bright ones, having consumed only $\sim 1/2$ of their gas. This resolves the gas consumption paradox, since there exist many galaxies with large gas reservoirs. These gas rich galaxies must have formed the bulk of their stellar populations in the last half of a Hubble time. The existence of such immature galaxies at $z = 0$ indicates that either galaxy formation is a lengthy or even ongoing process, or the onset of significant star formation can be delayed for arbitrary periods in tenuous gas disks.

The gas mass fraction (fg) of spiral galaxies is strongly correlated with the central surface brightness of their disks. There exist many dim galaxies with long gas consumption time scales and fg > 0.5. This resolves the gas consumption paradox. The surface density of gas follows the optical surface brightness, but does not vary by as large a factor. This is the signature of a critical density threshold for star formation. Such a mechanism seems to be responsible for the slow evolution of dim galaxies.

We demonstrate that, due to their low surface mass density and large dark matter content, LSB disks are quite stable against the growth of global bar modes. However, they may be only marginally stable against local disk instabilities. We simulate a collision between an LSB and HSB galaxy and find that, while the HSB galaxy forms a strong bar, the response of the LSB disk is milder, in the form of spiral features and an oval distortion. Unlike its HSB counterpart, the LSB disk does not suffer strong inflow of gas into the central regions. The lack of sufficient disk self-gravity to amplify dynamical instabilities makes it difficult to explain strong interaction-driven starbursts in LSB galaxies without invoking mergers.

Low Surface Brightness (LSB) galaxies form a large population of disc galaxies that extend the Hubble sequence towards extreme late-types. They are only slowly evolving, and still in an early evolutionary state. The Tully-Fisher relation and rotation curves of LSB galaxies both show that LSB galaxies are very dark-matter dominated with respect to ``normal'' high surface brightness (HSB) galaxies. Mass models derived from the rotation curves of LSB and HSB galaxies show that LSB galaxies inhabit less dense and more extended halos. Mass density, which changes with surface brightness, is as important in determining the evolution of a galaxy as total mass is.

Suto et al. (1996, astro-ph/9602061 and ApJ, 461, L33) have suggested that an X-ray halo in the Local Group might explain both the observed low-energy excess in the X-ray background and the quadrupole anisotropy in the cosmic microwave background. Recent observations of poor groups of galaxies by the ROSAT PSPC set reasonable limits on how extensive and dense such a halo could be. The poor groups most similar to the Local Group do not have a detectable halo, and the upper limits of these observations suggest that any Local Group halo would be nearly two orders of magnitude too tenuous to produce the effects Suto et al. (1996) discuss. In particular, the Sunyaev-Zel'dovich effect cannot contribute significantly to the quadrupole anisotropy measured by COBE.

We compare the dynamical properties of two galaxies at identical positions on the Tully-Fisher relation, but with different surface brightnesses. We find that the low surface brightness galaxy UGC 128 has a higher mass-to-light ratio, and yet has lower mass densities than the high surface brightness galaxy NGC 2403. This is true for the gas surface density, the stellar surface density, and the total mass density.

We present Very Large Array ({\sc vla}) and Westerbork Synthesis Radio Telescope ({\sc wsrt}) 21-cm H{\sc i} observations of 19 late-type low surface brightness (LSB) galaxies. Our main findings are that these galaxies, as well as having low surface brightnesses, have low H{\sc i} surface densities, about a factor of $\sim 3$ lower than in normal late-type galaxies. We show that LSB galaxies in some respects resemble the outer parts of late-type normal galaxies, but may be less evolved. LSB galaxies are more gas-rich than their high surface brightness counterparts. The rotation curves of LSB galaxies rise more slowly than those of HSB galaxies of the same luminosity, with amplitudes between 50 and 120~km~s$^{-1}$, and are often still increasing at the outermost measured point. The shape of the rotation curves suggests that LSB galaxies have low matter surface densities. We use the average total mass surface density of a galaxy as a measure for the evolutionary state, and show that LSB galaxies are among the least compact, least evolved galaxies. We show that both $M_{\rm HI}/L_B$ and $M_{\rm dyn}/L_B$ depend strongly on central surface brightness, consistent with the surface brightness--mass-to-light ratio relation required by the Tully-Fisher relation. LSB galaxies are therefore slowly evolving galaxies, and may well be low surface density systems in all respects.

I give analytic expressions for the relative number, luminosity, and mass density of disc galaxies as a function of surface brightness. These surface brightness distributions are asymmetric, with long tails to lower surface brightnesses. This asymmetry induces systematic errors in most determinations of the galaxy luminosity function. Galaxies of low surface brightness exist in large numbers, but the additional contribution to the integrated luminosity density is modest, probably 10 -- 30%.

I review the evolutionary status of the field populations of starbursting dwarf galaxies and low surface brightness spirals. Particular attention is paid to the long standing problem of the progenitors for star burst galaxies, which may be related to a population of small, low surface brightness dwarfs. Larger low surface brightness spirals are known to evolve slowly, converting little of their gas mass into stars and metals in a Hubble time. Their slow evolution is controlled by their low densities, which arise naturally from formation by late collapse of low amplitude perturbations in the initial density field.

Optical surveys for galaxies are biased against the inclusion of low surface brightness (LSB) galaxies. Disney (1976) suggested that the constancy of disk central surface brightness noticed by Freeman (1970) was not a physical result, but instead was an artifact of sample selection. Since LSB galaxies do exist, the pertinent and still controversial issue is if these newly discovered galaxies constitute a significant percentage of the general galaxy population. In this paper, we address this issue by determining the space density of galaxies as a function of disk central surface brightness. Using the physically reasonable assumption (which is motivated by the data) that central surface brightness is independent of disk scale length, we arrive at a distribution which is roughly flat (i.e. approximately equal numbers of galaxies at each surface brightness) faintwards of the Freeman (1970) value. Brightwards of this, we find a sharp decline in the distribution which is analogous to the turn down in the luminosity function at $L^*$. An intrinsically sharply peaked ``Freeman Law'' distribution can be completely ruled out, and no Gaussian distribution can fit the data. Low surface brightness galaxies (those with central surface brightnesses fainter than $22~B\sb$) comprise $\gtrsim 1/2$ the general galaxy population, so a representative sample of galaxies at $z = 0$ does not really exist at present since past surveys have been insensitive to this component of the general galaxy population.

We present the B band Tully-Fisher relation for Low Surface Brightness (LSB) galaxies. These LSB galaxies follow the same Tully-Fisher relation as normal spiral galaxies. This implies that the mass-to-light ratio (M/L) of LSB galaxies is typically a factor of 2 larger than that of normal galaxies of the same total luminosity and morphological type. Since the dynamical mass of a galaxy is related to the rotation velocity and scale length via M \propto V^2 h, at fixed linewidth LSB galaxies must be twice as large as normal galaxies. This is confirmed by examining the relation between scale length and linewidth for LSB and normal galaxies. The universal nature of the Tully-Fisher relation can be understood if LSB galaxies are galaxies with low mass surface density, \sigma. The mass surface density apparently controls the luminosity evolution of a galaxy such as to keep the product \sigma M/L constant.

We present $UBVI$ and H$\alpha$ images of a sample of Low Surface Brightness (LSB) disk galaxies. These galaxies are generally late types, if they can be sensibly classified at all. However, they are not dwarfs, being intrinsically large and luminous. The morphology of LSB galaxies is discussed in terms of the physical interpretation of the Hubble sequence. Galaxies with high contrast relative to the sky background are subject to being more finely typed than those which appear merely as fuzzy blobs on photographic plates. This causes the stages of the Hubble sequence to be nonlinear in the sense that large morphological type distinctions are made between high surface brightness spirals when only small physical differences exist, and small morphological distinctions are made between low surface brightness galaxies even when large physical differences exist. Many LSB galaxies lack the old red disk conspicuous in higher surface brightness spirals. Their morphology is strikingly similar in all bands from $U$ to $I$, suggesting farily homogeneous stellar populations lacking a well developed giant branch. These properties, together with their very blue colors, suggest that LSB galaxies are relatively younger than their high surface brightness counterparts. A few of these LSB galaxies appear to be very young ($\simlt 1$~Gyr), and as such may represent local examples of protogalaxies.

Low Surface-Brightness (LSB) galaxies are severely underrepresented in surveys used to define the local luminosity function (LF), but could it in principle be detected in deep surveys. To explore the possible contribution of such objects to faint galaxy counts, we construct catalogs of simulated non-evolving galaxies drawn from a multivariate distribution of galaxy luminosities, central surface brightnesses, bulge/disk ratios and spectral-energy distributions. We compare two models dominated by LSB galaxies to a standard non-evolving model. Model galaxies are convolved with seeing and selected in a manner that closely matches real surveys. For each model we compute the local LF, HI mass function, number counts, redshift and color distributions. We find it possible to include a large population of LSB galaxies and incorporate a steep faint-end slope of the LF in our simulations without violating the constraints on the local LF or HI mass function. For $q_0 = 0.5$, the most favorable model matches the counts to B=23, but falls short of the observations at fainter magnitudes. The colors and redshift distributions remain roughly consistent with observations to B=24. The most serious discrepancy with observations is in the distribution of $r_e$ at faint magnitudes, suggesting that the model contains too many LSB galaxies. Nevertheless, the results suggest that LSB galaxies could be a significant contributor to faint-galaxy counts, reducing the need for such extreme models of galaxy evolution as rapid merging or bursting dwarf galaxies.

OBSERVATIONS of galaxies to very faint magnitudes have revealed a population of blue galaxies at intermediate redshift$^{1-5}$. These galaxies represent a significant excess over the expectation of standard cosmological models for reasonable amounts of evolution of the locally observed galaxy population. However, the surveys which define the local galaxy population are strongly biased against objects of low surface brightness$^{6-9}$. Low surface brightness galaxies have properties very similar to those of the excess blue population$^{10,11}$, and recent work suggests that they are comparable in abundance to the more readily detected normal galaxies$^{9,12}$. I show that the very deep surveys which reveal the excess population can easily detect low surface brightness galaxies to large redshifts, but that local surveys will miss them because they are not comparably sensitive. This suggests that the excess faint galaxies {\it are\/} low surface brightness galaxies. No alteration of standard cosmology is required, but it is necessary to reconsider the way in which the galaxy distribution function is specified.

The oxygen abundances in the \HII regions of a sample of low surface brightness (LSB) disk galaxies are presented. In general, LSB galaxies are found to be metal poor ($Z < {1 \over 3} Z\solar$). Indeed, some LSB galaxies rival the lowest abundance extragalactic objects known. These low metallicities indicate that LSB galaxies evolve slowly, forming relatively few stars over a Hubble time. The low metallicities of LSB galaxies occur even though many are comparable in size and mass to the prominent spirals which define the Hubble sequence. As well as being low in surface brightness, these galaxies tend to be isolated. This suggests that surface mass density and environment are more relevant to galaxy evolution than gross size. Despite the low surface brightness of the disks, massive ($M > 60 M\solar$) stars are inferred to be present and no abnormality of the IMF is indicated. Many low excitation \HII regions exist at low metallicity in LSB galaxies, and the ionization parameter is not tightly correlated with metallicity. However, there does seem to be a significant envelope of maximum ionization at a given metallicity.

Using redshift samples, we calculate the cross-correlation functions of LSB galaxies with normal galaxies in complete samples (i.e. CfA and IRAS), which enables us to compare directly the amplitudes and shapes of the correlation functions. For pair separations $r\gs 2\mpc$, we find that the shape is in agreement with that of the correlation functions for other galaxies. The amplitudes ($A$) of $\xiab (r)$ are lower than those of the autocorrelation functions for the CfA and IRAS samples, with $A_{\rm LSB-CfA}:A_{\rm CfA-CfA}\approx 0.4$ and $A_{\rm LSB-IRAS}:A_{\rm IRAS-IRAS}\approx 0.6$. These results suggest that LSB galaxies are imbedded in the same large scale structure as other galaxies, but are less strongly clustered. This offers the hope that LSB galaxies may be unbiased tracers of the mass density on large scales. For $r\ls 2\mpc$, the cross-correlation functions are significantly lower than that expected from the extrapolation of $\xiab$ on larger scales, showing that the formation and survival of LSB galaxies may be inhibited by interaction with neighboring galaxies. The model which matches these observations suggests that strong luminosity segregation in galaxy clustering is not a necessary consequence of biased galaxy formation, unless the effect of surface brightness is taken in to account. It is also implies a significant mass density in LSB galaxies.

We present UBVI surface photometry of a sample of low surface brightness (LSB) disk galaxies. LSB disk galaxies are fairly well described as exponential disks with no preferred value for either scale length, central surface brightness, or rotational velocity. Indeed, the distribution of scale lengths is indistinguishable from that of high surface brightness spirals, indicating that dynamically similar galaxies (e.g., those with comparable Rv^2) exist over a large range in surface density. These LSB galaxies are strikingly blue. The complete lack of correlation between central surface brightness and color rules out any fading scenario. Similarly, the oxygen abundances inferred from HII region spectra are uncorrelated with color so the low metallicities are not the primary cause of the blue colors. While these are difficult to interpret in the absence of significant star formation, the most plausible scenario is a stellar population with a young mean age stemming from late formation and subsequent slow evolution. These properties suggest that LSB disks formed from low initial overdensities with correspondingly late collapse times.