We present a comprehensive test of the relation between stellar and total mass in galaxies as predicted by popular models based on abundance matching (AM) techniques. We use the ``Spectroscopy and H-band Imaging of Virgo cluster galaxies'' (SHIVir) survey with photometric and dynamical profiles for 190 Virgo cluster galaxies to establish a relation between the stellar and dynamical masses measured within the isophotal radius $r_{23.5}$. Various dark matter and galaxy scaling relations are combined with results from the NIHAO suite of hydrodynamical simulations to recast AM predictions in terms of these observed quantities. Our results are quite insensitive to the exact choice of dark matter profile and halo response to baryon collapse. We find that theoretical models reproduce the slope and normalization of the observed stellar-to-halo mass relation (SHMR) over more than three orders of magnitude in stellar mass $(10^8 < M_*/M_{\odot} < 2 \times 10^{11})$. However, the scatter of the observed SHMR exceeds that of AM predictions by a factor $\sim$5. For systems with stellar masses exceeding $5 \times 10^{10}~M_{\odot}$, AM overpredicts the observed stellar masses for a given dynamical mass. The latter offset may support previous indications of a different stellar initial mass function in these massive galaxies. Overall our results support the validity of AM predictions on a wide dynamical range.

Using the NIHAO suite of cosmological zoom-in simulations, we explore a parametric function for the dark-matter halo density profile subject to baryonic effects. This profile has variable inner slope and concentration parameter, and it allows analytic expressions for the gravitational potential, velocity dispersion, and lensing properties. We show that this profile provides better fits to simulated profiles than the Einasto profile and the generalized NFW profile with variable inner slope, in particular towards the halo centers, and that its parameters are correlated with the stellar-to-halo mass ratio $M_{\rm star}/M_{\rm vir}$. This defines a mass-dependent density profile describing the dark matter profiles in all galaxies, which can be directly applied to observed rotation curves of galaxies, gravitational lenses, and semi-analytic models of galaxy formation and satellite-galaxy evolution. The effect of baryons manifests itself in the profile parameters by a sharp drop of the inner density slope and a 20% decrease of the concentration parameter for $\log(M_{\rm star}/M_{\rm vir})$ between -3.5 and -2, which corresponds to $M_{\rm star}$ between $10^7$ and $10^{10}~\rm M_\odot$. The accuracy by which this profile fits simulated galaxies is similar to certain multi-parameter, mass-dependent profiles, but its fewer parameters and analytic nature make it most desirable for many purposes.

Simulating thin and extended galactic disks has long been a challenge in computational astrophysics. We introduce the NIHAO-UHD suite of cosmological hydrodynamical simulations of Milky Way mass galaxies and study stellar disk properties such as stellar mass, size and rotation velocity which agree well with observations of the Milky Way and local galaxies. In particular, the simulations reproduce the age-velocity dispersion relation and a multi-component stellar disk as observed for the Milky Way. Half of our galaxies show a double exponential vertical profile, while the others are well described by a single exponential model which we link to the disk merger history. In all cases, mono-age populations follow a single exponential whose scale height varies monotonically with stellar age and radius. The scale length decreases with stellar age while the scale height increases. The general structure of the stellar disks is already set at time of birth as a result of the inside-out and upside-down formation. Subsequent evolution modifies this structure by increasing both the scale length and height of all mono-age populations. Thus, our results put tight constraints on how much dynamical memory stellar disks can retain over cosmological timescales. Our simulations demonstrate that it is possible to form thin galactic disks in cosmological simulations provided there are no significant stellar mergers at low redshifts. Most of the stellar mass is formed in-situ with only a few percent ($\lesssim5\%$) brought in by merging satellites at early times. Redshift zero snapshots and halo catalogues are publicly available.

We present a simple model for the response of a dissipationless spherical system to an instantaneous mass change at its center, describing the formation of flat cores in dark matter haloes and ultra-diffuse galaxies (UDGs) from feedback-driven outflow episodes in a specific mass range. This model generalizes an earlier simplified analysis of an isolated shell into a system with continuous density, velocity and potential profiles. The response is divided into an instantaneous change of potential at constant velocities due to a given mass loss or gain, followed by energy-conserving relaxation to a new Jeans equilibrium. The halo profile is modeled by a two-parameter function with a variable inner slope and an analytic potential profile (Dekel et al. 2017), which enables determining the associated kinetic energy at equilibrium. The model is tested against NIHAO cosmological zoom-in simulations, where it successfully predicts the evolution of the inner dark-matter profile between successive snapshots in about 75% of the cases, failing mainly in merger situations. This model provides a simple understanding of the formation of dark-matter halo cores and UDGs by supernova-driven outflows, and a useful analytic tool for studying such processes.

This paper is one in a series that explores the importance of radius as a second parameter in galaxy evolution. The topic investigated here is the relationship between star formation rate (SFR) and galaxy radius ($R_{\rm e}$) for main-sequence star-forming galaxies. The key observational result is that, over a wide range of stellar mass and redshift in both CANDELS and SDSS, there is little trend between SFR and $R_{\rm e}$ at fixed stellar mass. The Kennicutt-Schmidt law, or any similar density-related star formation law, then implies that smaller galaxies must have lower gas fractions than larger galaxies (at fixed $M_{\ast}$), and this is supported by observations of local star-forming galaxies. We investigate the implication by adopting the equilibrium "bathtub" model: the ISM gas mass is assumed to be constant over time and the net star formation rate is the difference between the accretion rate of gas onto the galaxy from the halo and the outflow rate due to winds. To match the observed null correlation between SFR and radius, the bathtub model requires that smaller galaxies at fixed mass have weaker galactic winds. Our hypothesis is that galaxies are a 2-dimensional family whose properties are set mainly by halo mass and concentration. Galaxy radius and accretion rate plausibly both depend on halo concentration, which predicts how wind strength should vary with $R_{\rm e}$ and SFR.

We introduce algorithms for black hole physics, i.e., black hole formation, accretion and feedback, into the NIHAO (Numerical Investigation of a Hundred Astrophysical Objects) project of galaxy simulations. This enables us to study high mass, elliptical galaxies, where feedback from the central black hole is generally thought to have a significant effect on their evolution. We furthermore extend the NIHAO suite by 45 simulations that encompass $z=0$ halo masses from $1 \times 10^{12}$ to $4 \times 10^{13}\,\mathrm{M}_{\odot}$, and resimulate five galaxies from the original NIHAO sample with black hole physics, which have $z=0$ halo masses from $8 \times 10^{11}$ to $3 \times 10^{12}\,\mathrm{M}_{\odot}$. Now NIHAO contains 144 different galaxies and thus has the largest sample of zoom-in simulations of galaxies, spanning $z=0$ halo masses from $9 \times 10^{8}$ to $4 \times 10^{13}\,\mathrm{M}_{\odot}$. In this paper we focus on testing the algorithms and calibrating their free parameters against the stellar mass versus halo mass relation and the black hole mass versus stellar mass relation. We also investigate the scatter of these relations, which we find is a decreasing function with time and thus in agreement with observations. For our fiducial choice of parameters we successfully quench star formation in objects above a $z=0$ halo mass of $10^{12}\,\mathrm{M}_{\odot}$, thus transforming them into red and dead galaxies.

We implement an optically thin approximation for the effects of the local radiation field from stars and hot gas on the gas heating and cooling in the N-body SPH code GASOLINE2. We resimulate three galaxies from the NIHAO project: one dwarf, one Milky Way-like and one massive spiral, and study what are the local radiation field effects on various galaxy properties. We also study the effects of varying the Ultra Violet Background (UVB) model, by running the same galaxies with two different UVBs. Galaxy properties at $z=0$ like stellar mass, stellar effective mass radius, HI mass, and radial extent of the HI disc, show significant changes between the models with and without the local radiation field, and smaller differences between the two UVB models. The intrinsic effect of the local radiation field through cosmic time is to increase the equilibrium temperature at the interface between the galaxies and their circumgalactic media (CGM), moving this boundary inwards, while leaving relatively unchanged the gas inflow rate. Consequently, the temperature of the inflow increases when considering the local radiation sources. This temperature increase is a function of total galaxy mass, with a median CGM temperature difference of one order of magnitude for the massive spiral. The local radiation field suppresses the stellar mass growth by $\sim$20 per cent by $z=0$ for all three galaxies, while the HI mass is roughly halfed. The differences in the gas phase diagrams, significantly impact the HI column densities, shifting their peaks in the distributions towards lower $N_{\rm HI}$.

We present catalogs of calibrated photometry and spectroscopic redshifts in the Extended Groth Strip, intended for studies of photometric redshifts (photo-z's). The data includes ugriz photometry from CFHTLS and Y-band photometry from the Subaru Suprime camera, as well as spectroscopic redshifts from the DEEP2, DEEP3 and 3D-HST surveys. These catalogs incorporate corrections to produce effectively matched-aperture photometry across all bands, based upon object size information available in the catalog and Moffat profile point spread function fits. We test this catalog with a simple machine learning-based photometric redshift algorithm based upon Random Forest regression, and find that the corrected aperture photometry leads to significant improvement in photo-z accuracy compared to the original SExtractor catalogs from CFHTLS and Subaru. The deep ugrizY photometry and spectroscopic redshifts are well-suited for empirical tests of photometric redshift algorithms for LSST. The resulting catalogs are publicly available. We include a basic summary of the strategy of the DEEP3 Galaxy Redshift Survey to accompany the recent public release of DEEP3 data.

We use 25 simulated galaxies from the NIHAO project to define and characterize a variety of kinematic stellar structures: thin and thick discs, large scale single discs, classical and pseudo bulges, spheroids, inner discs, and stellar haloes. These structures have masses, spins, shapes and rotational support in good agreement with theoretical expectations and observational data. Above a dark matter halo mass of $2.5\times10^{\rm~11}M_{\rm\odot}$, all galaxies have a classical bulge and 70\% have a thin and thick disc. The kinematic (thin) discs follow a power-law relation between angular momentum and stellar mass $J_{\rm *}=3.4M_{\rm *}^{\rm1.26\pm0.06}$, in very good agreement with the prediction based on the empirical stellar-to-halo mass relation in the same mass range, and show a strong correlation between maximum `observed' rotation velocity and dark matter halo circular velocity $v_{\rm c}=6.4v_{\rm max}^{0.64\pm0.04}$. Tracing back in time these structures' progenitors, we find all to lose a fraction $1-f_j$ of their maximum angular momentum. Thin discs are significantly better at retaining their high-redshift spins ($f_j\sim0.70$) than thick ones ($f_j\sim0.40$). Stellar haloes have their progenitor baryons assembled the latest ($z_{\rm~1/2}\sim1.1$) and over the longest timescales ($\tau\sim6.2$~Gyr), and have the smallest fraction of stars born in-situ ($f_{\rm in-situ}=0.35\pm0.14$). All other structures have $1.5\lesssim z_{\rm1/2}\lesssim3$, $\tau=4\pm2$~Gyr and $f_{\rm in-situ}\gtrsim0.9$.

We use the NIHAO simulations to investigate the effects of baryonic physics on the time evolution of Dark Matter central density profiles. The sample is made of $\approx 70$ independent high resolution hydrodynamical simulations of galaxy formation and covers a wide mass range: 1e10< Mhalo <1e12, i.e., from dwarfs to L* . We confirm previous results on the dependence of the inner dark matter density slope, $\alpha$, on the ratio between stellar-to-halo mass. We show that this relation holds approximately at all redshifts (with an intrinsic scatter of ~0.18 in $\alpha$). This implies that in practically all haloes the shape of their inner density profile changes quite substantially over cosmic time, as they grow in stellar and total mass. Thus, depending on their final stellar-to-halo mass ratio, haloes can either form and keep a substantial density core (size~1 kpc), or form and then destroy the core and re-contract the halo, going back to a cuspy profile, which is even steeper than CDM predictions for massive galaxies (~1e12 Msun). We show that results from the NIHAO suite are in good agreement with recent observational measurements of $\alpha$ in dwarf galaxies. Overall our results suggest that the notion of a universal density profile for dark matter haloes is no longer valid in the presence of galaxy formation.

We use thirty-eight high-resolution simulations of galaxy formation between redshift 10 and 5 to study the impact of a 3 keV warm dark matter (WDM) candidate on the high-redshift Universe. We focus our attention on the stellar mass function and the global star formation rate and consider the consequences for reionization, namely the neutral hydrogen fraction evolution and the electron scattering optical depth. We find that three different effects contribute to differentiate warm and cold dark matter (CDM) predictions: WDM suppresses the number of haloes with mass less than few $10^9$ M$_{\odot}$; at a fixed halo mass, WDM produces fewer stars than CDM; and finally at halo masses below $10^9$ M$_{\odot}$, WDM has a larger fraction of dark haloes than CDM post-reionization. These three effects combine to produce a lower stellar mass function in WDM for galaxies with stellar masses at and below $\sim 10^7$ M$_{\odot}$. For $z > 7$, the global star formation density is lower by a factor of two in the WDM scenario, and for a fixed escape fraction, the fraction of neutral hydrogen is higher by 0.3 at $z \sim 6$. This latter quantity can be partially reconciled with CDM and observations only by increasing the escape fraction from 23 per cent to 34 per cent. Overall, our study shows that galaxy formation simulations at high redshift are a key tool to differentiate between dark matter candidates given a model for baryonic physics.

State-of-the-art cosmological hydrodynamical simulations of galaxy formation have reached the point at which their outcomes result in galaxies with ever more realism. Still, the employed sub-grid models include several free parameters such as the density threshold, $n$, to localize the star-forming gas. In this work, we investigate the possibilities to utilize the observed clustered nature of star formation (SF) in order to refine SF prescriptions and constrain the density threshold parameter. To this end, we measure the clustering strength, correlation length and power-law index of the two-point correlation function of young ($\tau<50$ Myr) stellar particles and compare our results to observations from the HST Legacy Extragalactic UV Survey (LEGUS). Our simulations reveal a clear trend of larger clustering signal and power-law index and lower correlation length as the SF threshold increases with only mild dependence on galaxy properties such as stellar mass or specific star formation rate. In conclusion, we find that the observed clustering of SF is inconsistent with a low threshold for SF ($n<1$ cm$^{-3}$) and strongly favours a high value for the density threshold of SF ($n>10$ cm$^{-3}$), as for example employed in the NIHAO project.

The existence of galaxies with a surface brightness $\mu$ lower than the night sky has been known since three decades. Yet, their formation mechanism and emergence within a $\rm\Lambda CDM$ universe has remained largely undetermined. For the first time, we investigated the origin of Low Surface Brightness (LSB) galaxies with M$_{\star}$$\sim$10$^{9.5-10}$M$_{\odot}$, which we are able to reproduce within hydrodynamical cosmological simulations from the NIHAO suite. The simulated and observed LSBs share similar properties, having large HI reservoir, extended star formation histories and effective radii, low S\'{e}rsic index and slowly rising rotation curves. The formation mechanism of these objects is explored: simulated LSBs form as a result of co-planar co-rotating mergers and aligned accretion of gas at early times, while perpendicular mergers and mis-aligned gas accretion result in higher $\mu$ galaxies by $z$=0. The larger the merger, the stronger the correlation between merger orbital configuration and final $\mu$. While the halo spin parameter is consistently high in simulated LSB galaxies, the impact of halo concentration, feedback-driven gas outflows and merger time only plays a minor-to-no role in determining $\mu$. Interestingly, the formation scenario of such `classical' LSBs differs from the one of less massive, M$_{\star}$$\sim$10$^{7-9}$M$_{\odot}$, Ultra-Diffuse Galaxies, the latter resulting from the effects of SNae driven gas outflows: a M$_{\star}$ of $\sim$10$^9$M$_{\odot}$ thus represents the transition regime between a feedback-dominated to an angular momentum-dominated formation scenario in the LSB realm. Observational predictions are offered regarding spatially resolved star formation rates through LSB discs: these, together with upcoming surveys, can be used to verify the proposed emergence scenario of LSB galaxies.

We use cosmological hydrodynamical galaxy formation simulations from the NIHAO project to investigate the impact of the threshold for star formation on the response of the dark matter (DM) halo to baryonic processes. The fiducial NIHAO threshold, $n=10\, {\rm cm}^{-3}$, results in strong expansion of the DM halo in galaxies with stellar masses in the range $10^{7.5} < M_{star} < 10^{9.5} M_{\odot}$. We find that lower thresholds such as $n=0.1$ (as employed by the EAGLE/APOSTLE and Illustris/AURIGA projects) do not result in significant halo expansion at any mass scale. Halo expansion driven by supernova feedback requires significant fluctuations in the local gas fraction on sub-dynamical times (i.e., < 50 Myr at galaxy half-light radii), which are themselves caused by variability in the star formation rate. At one per cent of the virial radius, simulations with $n=10$ have gas fractions of $\simeq 0.2$ and variations of $\simeq 0.1$, while $n=0.1$ simulations have order of magnitude lower gas fractions and hence do not expand the halo. The observed DM circular velocities of nearby dwarf galaxies are inconsistent with CDM simulations with $n=0.1$ and $n=1$, but in reasonable agreement with $n=10$. Star formation rates are more variable for higher $n$, lower galaxy masses, and when star formation is measured on shorter time scales. For example, simulations with $n=10$ have up to 0.4 dex higher scatter in specific star formation rates than simulations with $n=0.1$. Thus observationally constraining the sub-grid model for star formation, and hence the nature of DM, should be possible in the near future.

We derive the mass weighted total density slopes within the effective (half-light) radius, $\gamma'$, for more than 2000 nearby galaxies from the SDSS-IV MaNGA survey using Jeans-anisotropic-models applied to IFU observations. Our galaxies span a wide range of the stellar mass ($10^9$ $M_{\rm \odot}< M_* < 10^{12}$ M$_{\odot}$) and the velocity dispersion (30 km/s $< \sigma_v <$ 300 km/s). We find that for galaxies with velocity dispersion $\sigma_v>100$ km/s, the density slope has a mean value $\langle \gamma^{\prime} \rangle = 2.24$ and a dispersion $\sigma_{\gamma}=0.22$, almost independent of velocity dispersion. A clear turn over in the $\gamma'-\sigma_v$ relation is present at $\sigma\sim 100$ km/s, below which the density slope decreases rapidly with $\sigma_v$. Our analysis shows that a large fraction of dwarf galaxies (below $M_* = 10^{10}$ M$_{\odot}$) have total density slopes shallower than 1, which implies that they may reside in cold dark matter halos with shallow density slopes. We compare our results with that of galaxies in hydrodynamical simulations of EAGLE, Illustris and IllustrisTNG projects, and find all simulations predict shallower density slopes for massive galaxies with high $\sigma_v$. Finally, we explore the dependence of $\gamma'$ on the positions of galaxies in halos, namely centrals vs. satellites, and find that for the same velocity dispersion, the amplitude of $\gamma'$ is higher for satellite galaxies by about 0.1.

The stellar populations in the inner kiloparsecs of the Milky Way (MW) show complex kinematical and chemical structures. The origin and evolution of these structures is still under debate. Here we study the central region of a fully cosmological hydrodynamical simulation of a disk galaxy that reproduces key properties of the inner kiloparsecs of the MW: it has a boxy morphology and shows an overall rotation and dispersion profile in agreement with observations. We use a clustering algorithm on stellar kinematics to identify a number of discrete kinematic components: a high- and low-spin disk, a stellar halo and two bulge components; one fast rotating and one slow-rotating. We focus on the two bulge components and show that the slow rotating one is spherically symmetric while the fast rotating component shows a boxy/peanut morphology. Although the two bulge components are kinematically discrete populations at present-day, they are both mostly formed over similar time scales, from disk material. We find that stellar particles with lower initial birth angular momentum (most likely thick disc stars) end up in the slow-rotating low-spin bulge, while stars with higher birth angular momentum (most likely thin disc stars) are found in the high-spin bulge. This has the important consequence that a bulge population with a spheroidal morphology does not necessarily indicate a merger origin. In fact, we do find that only $\sim2.3$\% of the stars in the bulge components are ex-situ stars brought in by accreted dwarf galaxies early on. We identify these ex-situ stars as the oldest and most metal-poor stars on highly radial orbits with large vertical excursions from the disk.

The phenomenological basis for Modified Newtonian Dynamics (MOND) is the radial-acceleration-relation (RAR) between the observed acceleration, $a=V^2_{rot}(r)/r$, and the acceleration accounted for by the observed baryons (stars and cold gas), $a_{bar}=V_{bar}^2(r)/r$. We show that the RAR arises naturally in the NIHAO sample of 89 high-resolution LCDM cosmological galaxy formation simulations. The overall scatter from NIHAO is just 0.079 dex, consistent with observational constraints. However, we show that the scatter depends on stellar mass. At high masses ($10^9 <M_{star} <10^{11}$ Msun) the simulated scatter is just $\simeq 0.04$ dex, increasing to $\simeq 0.11$ dex at low masses ($10^7 < M_{star} <10^{9}$Msun). Observations show a similar dependence for the intrinsic scatter. At high masses the intrinsic scatter is consistent with the zero scatter assumed by MOND, but at low masses the intrinsic scatter is non-zero, strongly disfavoring MOND. Applying MOND to our simulations yields remarkably good fits to most of the circular velocity profiles. In cases of mild disagreement the stellar mass-to-light ratio and/or "distance" can be tuned to yield acceptable fits, as is often done in observational mass models. In dwarf galaxies with $M_{star}\sim10^6$Msun MOND breaks down, predicting lower accelerations than observed and in our LCDM simulations. The assumptions that MOND is based on (e.g., asymptotically flat rotation curves, zero intrinsic scatter in the RAR), are approximately, but not exactly, true in LCDM. Thus if one wishes to go beyond Newtonian dynamics there is more freedom in the RAR than assumed by MOND.

We have used the NIHAO simulations to explore how supernovae (SNe) affect star formation in galaxies. We find that SN feedback operates on all scales from the interstellar medium (ISM) to several virial radii. SNe regulate star formation by preventing condensation of HI into H$_2$ and by moving cold neutral gas to the hot HII phase. The first effect explains why the cold neutral gas in dwarf galaxies forms stars inefficiently. The second maintains the hot ISM of massive galaxies (HII vents out at lower masses). At $v_{\rm vir}\simeq 67{\rm\,km\,s}^{-1}$, the outflow rate follows the relation: $\dot{M}_{\rm out}=23\,(v_{\rm vir}/67{\rm\,km\,s}^{-1})^{-4.6}\,{\rm SFR}$. $20\%$ to $70\%$ of the gas expelled from galaxies escapes from the halo (ejective feedback) but outflows are dominated by cold swept-up gas, most of which falls back onto the galaxy on a $\sim 1\,$Gyr timescale. This `fountain feedback' reduces the masses of galaxies by a factor of two to four, since gas spends half to three quarter of its time in the fountain. Less than $10\%$ of the ejected gas mixes with the hot circumgalactic medium and this gas is usually not reaccreted. On scales as large as $6r_{\rm vir}$, galactic winds divert the incoming gas from cosmic filaments and prevent if from accreting onto galaxies (pre-emptive feedback). This process is the main reason for the low baryon content of ultradwarves.

In this third paper of the series, we investigate the effects of warm dark matter with a particle mass of $m_\mathrm{WDM}=3\,\mathrm{keV}$ on the smallest galaxies in our Universe. We present a sample of 21 hydrodynamical cosmological simulations of dwarf galaxies and 20 simulations of satellite-host galaxy interaction that we performed both in a Cold Dark Matter (CDM) and Warm Dark Matter (WDM) scenario. In the WDM simulations, we observe a higher critical mass for the onset of star formation. Structure growth is delayed in WDM, as a result WDM haloes have a stellar population on average two Gyrs younger than their CDM counterparts. Nevertheless, despite this delayed star formation, CDM and WDM galaxies are both able to reproduce the observed scaling relations for velocity dispersion, stellar mass, size, and metallicity at $z=0$. WDM satellite haloes in a Milky Way mass host are more susceptible to tidal stripping due to their lower concentrations, but their galaxies can even survive longer than the CDM counterparts if they live in a dark matter halo with a steeper central slope. In agreement with our previous CDM satellite study we observe a steepening of the WDM satellites' central dark matter density slope due to stripping. The difference in the average stellar age for satellite galaxies, between CDM and WDM, could be used in the future for disentangling these two models.

We study the components of cool and warm/hot gas in the circumgalactic medium (CGM) of simulated galaxies and address the relative production of OVI by photoionization versus collisional ionization, as a function of halo mass, redshift, and distance from the galaxy halo center. This is done utilizing two different suites of zoom-in hydro-cosmological simulations, VELA (6 halos; $z>1$) and NIHAO (18 halos; to $z=0$), which provide a broad theoretical basis because they use different codes and physical recipes for star formation and feedback. In all halos studied in this work, we find that collisional ionization by thermal electrons dominates at high redshift, while photoionization of cool or warm gas by the metagalactic radiation takes over near $z\sim2$. In halos of $\sim 10^{12}M_{\odot}$ and above, collisions become important again at $z<0.5$, while photoionization remains significant down to $z=0$ for less massive halos. In halos with $M_{\textrm v}>3\times10^{11}~M_{\odot}$, at $z\sim 0$ most of the photoionized OVI is in a warm, not cool, gas phase ($T\lesssim 3\times 10^5$~K). We also find that collisions are dominant in the central regions of halos, while photoionization is more significant at the outskirts, around $R_{\textrm v}$, even in massive halos. This too may be explained by the presence of warm gas or, in lower mass halos, by cool gas inflows.

We study ultra-diffuse galaxies (UDGs) in zoom in cosmological simulations, seeking the origin of UDGs in the field versus galaxy groups. We find that while field UDGs arise from dwarfs in a characteristic mass range by multiple episodes of supernova feedback (Di Cintio et al. 2017), group UDGs may also form by tidal puffing up and they become quiescent by ram-pressure stripping. The field and group UDGs share similar properties, independent of distance from the group centre. Their dark-matter haloes have ordinary spin parameters and centrally dominant dark-matter cores. Their stellar components tend to have a prolate shape with a Sersic index n~1 but no significant rotation. Ram pressure removes the gas from the group UDGs when they are at pericentre, quenching star formation in them and making them redder. This generates a colour/star-formation-rate gradient with distance from the centre, as observed in clusters. We find that ~20 per cent of the field UDGs that fall into a massive halo survive as satellite UDGs. In addition, normal field dwarfs on highly eccentric orbits can become UDGs near pericentre due to tidal puffing up, contributing about half of the group-UDG population. We interpret our findings using simple toy models, showing that gas stripping is mostly due to ram pressure rather than tides. We estimate that the energy deposited by tides in the bound component of a satellite over one orbit can cause significant puffing up provided that the orbit is sufficiently eccentric.

We use 85 pairs of high resolution LCDM cosmological simulations from the NIHAO project to investigate why in dwarf galaxies neutral hydrogen (HI) linewidths measured at 50% of the peak flux W_50 /2 (from the hydrodynamical simulations) tend to underpredict the maximum circular velocity VmaxDMO (from the corresponding dark matter only simulations). There are two main contributing processes. 1) Lower mass galaxies are less rotationally supported. This is confirmed observationally from the skewness of linewidths in bins of HI mass in both ALFALFA and HIPASS observations. 2) The HI distributions are less extended (relative to the dark matter halo) in dwarf galaxies. Coupled to the lower baryon-to-halo ratio results in rotation curves that are still rising at the last measured data point, in agreement with observations from SPARC. Combining these two effects, in both simulations and observations lower mass galaxies have on average lower W_50 / W_20. Additionally, mass loss driven by supernovae and projection effects (dwarf galaxies are in general not thin disks) further reduce the linewidths. The implied HI linewidth velocity function from NIHAO is in good agreement with observations in the nearby Universe of dwarf galaxies: 10 < W_50 /2 < 80 km/s. The dark matter only slope of -2.9 is reduced to -1.0 in the hydro simulations. Future radio observations of unbiased samples with higher spatial resolution will enable stricter tests of the simulations, and thus of the LCDM model.

The NIHAO cosmological simulations form a collection of a hundred high-resolution galaxies. We used these simulations to test the impact of stellar feedback on the morphology of the HI distribution in galaxies. We ran a subsample of twenty of the galaxies with different parameterizations of stellar feedback, looking for differences in the HI spatial distribution and morphology. We found that different feedback models do leave a signature in HI, and can potentially be compared with current and future observations. These findings can help inform future modeling efforts in the parameterization of stellar feedback in cosmological simulations of galaxy formation and evolution.

We used high resolution dark matter only cosmological simulations to investigate the structural properties of Lambda Cold Dark Matter ($\Lambda$CDM) haloes over cosmic time. The haloes in our study range in mass from $\sim 10^{10}$ to $\sim 10^{12} \mathrm{M}_\odot$, and are resolved with $10^5$ to $10^7$ particles. We fit the spherically averaged density profiles of DM haloes with the three parameter Einasto function. For our sample of haloes, the Einasto shape parameter, $\alpha$, is uncorrelated with the concentration, $c$, at fixed halo mass, and at all redshifts. Previous reports of an anti-correlation are traced to fitting degeneracies, which our fits are less sensitive to due to our higher spatial resolution. However, for individual haloes the evolution in $\alpha$ and $c$ is anti-correlated: at redshift $z=7$, $\alpha \simeq 0.4$ and decreases with time, while $c\simeq 3$ and increases with time. The evolution in structure is primarily due to accretion of mass at larger radii. We suggest that $\alpha$ traces the evolutionary state of the halo, with dynamically young haloes having high $\alpha$ (closer to a top-hat: $\alpha^{-1}=0$), and dynamically relaxed haloes having low $\alpha$ (closer to isothermal: $\alpha=0$). Such an evolutionary dependence reconciles the increase of $\alpha$ vs peak height, $\nu$, with the dependence on the slope of the power spectrum of initial density fluctuations found by previous studies.

We study the origin of bulge stars and their angular momentum (AM) evolution in 10 spiral galaxies with baryonic masses above $10^{10}$M$_\odot$ in the NIHAO galaxy formation simulations. The simulated galaxies are in good agreement with observations of the relation between specific AM and mass of the baryonic component and the stellar bulge-to-total ratio ($B/T$). We divide the star particles at $z=0$ into disc and bulge components using a hybrid photometric/kinematic decomposition method that identifies all central mass above an exponential disc profile as the `bulge'. By tracking the bulge star particles back in time, we find that on average 95\% of the bulge stars formed {\it in situ}, 3\% formed {\it ex situ} in satellites of the same halo, and only 2\% formed {\it ex situ} in external galaxies. The evolution of the AM distribution of the bulge stars paints an interesting picture: the higher the final $B/T$ ratio, the more the specific AM remains preserved during the bulge formation. In all cases, bulge stars migrate significantly towards the central region, reducing their average galactocentric radius by roughly a factor 2, independently of the final $B/T$ value. However, in the higher $B/T$ ($\gtrsim0.2$) objects, the velocity of the bulge stars increases and the AM of the bulge is almost conserved, whereas at lower $B/T$ values, the velocity of the bulge stars decreases and the AM of bulge reduces. The correlation between the evolution of the AM and $B/T$ suggests that bulge and disc formation are closely linked and cannot be treated as independent processes.

We analyze 88 independent high-resolution cosmological zoom-in simulations of disk galaxies in the NIHAO simulations suite to explore the connection between the atomic gas fraction and angular momentum of baryons throughout cosmic time. The study is motivated by the analytic model of \citet{obreschkow16}, which predicts a relation between the atomic gas fraction $f_{\rm atm}$ and the global atomic stability parameter $q \equiv j\sigma / (GM)$, where $M$ and $j$ are the mass and specific angular momentum of the galaxy (stars+cold gas) and $\sigma$ is the velocity dispersion of the atomic gas. We show that the simulated galaxies follow this relation from their formation ($z\simeq4$) to present within $\sim 0.5$ dex. To explain this behavior, we explore the evolution of the local Toomre stability and find that $90\%$--$100\%$ of the atomic gas in all simulated galaxies is stable at any time. In other words, throughout the entire epoch of peak star formation until today, the timescale for accretion is longer than the timescale to reach equilibrium, thus resulting in a quasi-static equilibrium of atomic gas at any time. Hence, the evolution of $f_{\rm atm}$ depends on the complex hierarchical growth history primarily via the evolution of $q$. An exception are galaxies subject to strong environmental effects.

We study the impact of the host on dwarf galaxy properties using four new Milky Way-like, ultra high-resolution simulations, ($N_{\rm part} >10^7$) from the NIHAO project. We split our sample in satellite ($R<R_{200}$), nearby ($1<R/R_{200}<2.5$), and field ($R>2.5 R_{\rm 200}$) galaxies. Simulated galaxies from all three groups are in excellent agreement with Local Group dwarf galaxies in terms of: stellar mass-velocity dispersion, stellar mass-metallicity relation, star formation histories, and stellar mass functions. Satellites and nearby galaxies show lower velocity dispersions and gas fractions compared to field galaxies. While field galaxies follow global abundance matching relations, satellites and nearby galaxies deviate from them, showing lower dark matter masses for given stellar mass. The reason for this deficit in dark matter mass is substantial mass loss experienced by satellites and $\sim80$\% of the nearby galaxies, while orbiting inside $R_{200}$ at earlier times. However, both satellites and nearby objects fall back onto the relation for field galaxies if we use the maximum of their virial mass instead of the present-day value. This allows us to provide estimates for the peak masses of observed Local Group galaxies. Finally, using radial velocities, distances, and the velocity dispersion-stellar mass relation from our simulations, we derive a metric to distinguish between galaxies harassed by the central object and unaffected ones. Applying this metric to observed objects we find that even far away dwarf galaxies like Eri II ($D\approx$ 370 kpc) have a strong probability ($\approx 83$\%) of having been affected by the Milky Way in the past. This naturally explains the lack of gas and recent star formation seen in Eri II.

We show for the first time, that a fully cosmological hydrodynamical simulation can reproduce key properties of the innermost region of the Milky Way. Our high resolution simulation matches the profile and kinematics of the Milky Way's boxy/peanut-shaped bulge, and hence we can use it to reconstruct and understand the bulge assembly. In particular, the age dependence of the X-shape morphology of the simulated bulge parallels the observed metallicity dependent split in the red clump stars of the inner Galaxy. We use this feature to derive an observational metric that allows us to quantify when the bulge formed from the disk. The metric we propose can be employed with upcoming survey data to constrain the age of the Milky Way bar. From the split in stellar counts we estimate the formation of the 4~kpc scale bar in the simulation to have happened $t^{\rm bar}_{\rm form}\sim8^{+2}_{-2}$ Gyr ago, in good agreement with conventional methods to measure bar formation in simulations. We test the prospects for observationally differentiating the stars that belong to the bulge/bar compared to the surrounding disk, and find that the inner disk and bulge are practically indistinguishable in both chemistry and ages.

We present the first results of applying Gaussian Mixture Models in the stellar kinematic space of normalized angular momentum and binding energy on NIHAO high resolution galaxies to separate the stars into multiple components. We exemplify this method using a simulated Milky Way analogue, whose stellar component hosts: thin and thick discs, classical and pseudo bulges, and a stellar halo. The properties of these stellar structures are in good agreement with observational expectations in terms of sizes, shapes and rotational support. Interestingly, the two kinematic discs show surface mass density profiles more centrally concentrated than exponentials, while the bulges and the stellar halo are purely exponential. We trace back in time the Lagrangian mass of each component separately to study their formation history. Between z~3 and the end of halo virialization, z~1.3, all components lose a fraction of their angular momentum. The classical bulge loses the most (~95%) and the thin disc the least (~60%). Both bulges formed their stars in-situ at high redshift, while the thin disc formed ~98% in-situ, but with a constant SFR~1.5M$_{\rm\odot}$/yr$^{\rm-1}$ over the last ~11 Gyr. Accreted stars (6% of total stellar mass) are mainly incorporated to the thick disc or the stellar halo, which formed ex-situ 8% and 45% of their respective masses. Our analysis pipeline is freely available at https://github.com/aobr/gsf.

The similarity between the distributions of spins for galaxies ($\lambda_{\rm g}$) and for dark-matter haloes ($\lambda_{\rm h}$), indicated both by simulations and observations, is naively interpreted as a one-to-one correlation between the spins of a galaxy and its host halo. This is used to predict galaxy sizes in semi-analytic models via $R_{\rm e}\simeq\lambda_{\rm h} R_{\rm v}$, with $R_{\rm e}$ the half-mass radius of the galaxy and $R_{\rm v}$ the halo radius. Utilizing two different suites of zoom-in cosmological simulations, we find that $\lambda_{\rm g}$ and $\lambda_{\rm h}$ are in fact only barely correlated, especially at $z\geq 1$. A general smearing of this correlation is expected based on the different spin histories, where the more recently accreted baryons through streams gain and then lose significant angular momentum compared to the gradually accumulated dark matter. Expecting the spins of baryons and dark matter to be correlated at accretion into $R_{\rm v}$, the null correlation at the end reflects an anti-correlation between $\lambda_{\rm g}/\lambda_{\rm h}$ and $\lambda_{\rm h}$, which can partly arise from mergers and a compact star-forming phase that many galaxies undergo. On the other hand, the halo and galaxy spin vectors tend to be aligned, with a median $\cos\theta=0.6$-0.7 between galaxy and halo, consistent with instreaming within a preferred plane. The galaxy spin is better correlated with the spin of the inner halo, but this largely reflects the effect of the baryons on the halo. Following the null spin correlation, $\lambda_{\rm h}$ is not a useful proxy for $R_{\rm e}$. While our simulations reproduce a general relation of the sort $R_{\rm e}=AR_{\rm vir}$, in agreement with observational estimates, the relation becomes tighter with $A=0.02(c/10)^{-0.7}$, where $c$ is the halo concentration, which in turn introduces a dependence on mass and redshift.

Galaxy formation entails the hierarchical assembly of mass, along with the condensation of baryons and the ensuing, self-regulating star formation. The stars form a collisionless system whose orbit distribution retains dynamical memory that can constrain a galaxy's formation history. The ordered-rotation dominated orbits with near maximum circularity $\lambda_z \simeq1$ and the random-motion dominated orbits with low circularity $\lambda_z \simeq0$ are called kinematically cold and kinematically hot, respectively. The fraction of stars on `cold' orbits, compared to the fraction of stars on `hot' orbits, speaks directly to the quiescence or violence of the galaxies' formation histories. Here we present such orbit distributions, derived from stellar kinematic maps via orbit-based modelling for a well defined, large sample of 300 nearby galaxies. The sample, drawn from the CALIFA survey, includes the main morphological galaxy types and spans the total stellar mass range from $10^{8.7}$ to $10^{11.9}$ solar masses. Our analysis derives the orbit-circularity distribution as a function of galaxy mass, $p(\lambda_z~|~M_\star)$, and its volume-averaged total distribution, $p(\lambda_z)$. We find that across most of the considered mass range and across morphological types, there are more stars on `warm' orbits defined as $0.25\le \lambda_z \le 0.8$ than on either `cold' or `hot' orbits. This orbit-based "Hubble diagram" provides a benchmark for galaxy formation simulations in a cosmological context.

The significant diversity of rotation curve (RC) shapes in dwarf galaxies has recently emerged as a challenge to LCDM: in dark matter (DM) only simulations, DM halos have a universal cuspy density profile that results in self-similar RC shapes. We compare RC shapes of simulated galaxies from the NIHAO project with observed galaxies from the homogeneous SPARC dataset. The DM halos of the NIHAO galaxies can expand to form cores, with the degree of expansion depending on their stellar-to-halo mass ratio. By means of the $V_{\rm 2kpc}-V_{\rm Rlast}$ relation (where $V_{\rm Rlast}$ is the outermost measured rotation velocity), we show that both the average trend and the scatter in RC shapes of NIHAO galaxies are in reasonable agreement with SPARC: this represents a significant improvement compared to simulations that do not result in DM core formation, suggesting that halo expansion is a key process in matching the diversity of dwarf galaxy RCs. Note that NIHAO galaxies can reproduce even the extremely slowly rising RCs of IC 2574 and UGC 5750. Revealingly, the range where observed galaxies show the highest diversity corresponds to the range where core formation is most efficient in NIHAO simulations, 50$<V_{\rm Rlast}$\km s$^{-1}<$100. A few observed galaxies in this range cannot be matched by any NIHAO RC nor by simulations that predict a universal halo profile. Interestingly, the majority of these are starbursts or emission-line galaxies, with steep RCs and small effective radii. Such galaxies represent an interesting observational target providing new clues to the process/viability of cusp-core transformation, the relationship between starburst and inner potential well, and the nature of DM.

In an effort to search for Ly$\alpha$ emission from circum- and intergalactic gas on scales of hundreds of kpc around $z\sim3$ quasars, and thus characterise the physical properties of the gas in emission, we have initiated an extensive fast-survey with the Multi Unit Spectroscopic Explorer (MUSE): Quasar Snapshot Observations with MUse: Search for Extended Ultraviolet eMission (QSO MUSEUM). In this work, we report the discovery of an enormous Ly$\alpha$ nebula (ELAN) around the quasar SDSS~J102009.99+104002.7 at $z=3.164$, which we followed-up with deeper MUSE observations. This ELAN spans $\sim297$ projected kpc, has an average Ly$\alpha$ surface brightness ${\rm SB}_{\rm Ly\alpha}\sim 6.04\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ (within the $2\sigma$ isophote), and is associated with an additional four, previously unknown embedded sources: two Ly$\alpha$ emitters and two faint active galactic nuclei (one Type-1 and one Type-2 quasar). By mapping at high significance the line-of-sight velocity in the entirety of the observed structure, we unveiled a large-scale coherent rotation-like pattern spanning $\sim300$ km s$^{-1}$ with a velocity dispersion of $<270$ km s$^{-1}$, which we interpret as a signature of the inspiraling accretion of substructures within the quasar's host halo. Future multiwavelength data will complement our MUSE observations, and are definitely needed to fully characterise such a complex system. None the less, our observations reveal the potential of new sensitive integral-field spectrographs to characterise the dynamical state of diffuse gas on large scales in the young Universe, and thereby witness the assembly of galaxies.

In the first paper we presented 27 hydrodynamical cosmological simulations of galaxies with total masses between $5 \times 10^8$ and $10^{10}\,\mathrm{M}_\odot$. In this second paper we use a subset of these cosmological simulations as initial conditions (ICs) for more than forty hydrodynamical simulations of satellite and host galaxy interaction. Our cosmological ICs seem to suggest that galaxies on these mass scales have very little rotational support and are velocity dispersion ($\sigma$) dominated. Accretion and environmental effects increase the scatter in the galaxy scaling relations (e.g. size - velocity dispersion) in very good agreement with observations. Star formation is substantially quenched after accretion. Mass removal due to tidal forces has several effects: it creates a very flat stellar velocity dispersion profiles, and it reduces the dark matter content at all scales (even in the centre), which in turn lowers the stellar velocity on scales around 0.5 kpc even when the galaxy does not lose stellar mass. Satellites that start with a cored dark matter profile are more prone to either be destroyed or to end up in a very dark matter poor galaxy. Finally, we found that tidal effects always increase the "cuspyness" of the dark matter profile, even for haloes that infall with a core.

The satellites of the Milky Way and Andromeda represent the smallest galaxies we can observe in our Universe. In this series of papers we aim to shed light on their formation and evolution using cosmological hydrodynamical simulations. In this first paper we focus on the galaxy properties before accretion, by simulating twenty seven haloes with masses between $5\times 10^8$ and $10^{10} M_\odot$. Out of this set nineteen haloes successfully form stars, while eight remain dark. The simulated galaxies match quite well present day observed scaling relations between stellar mass, size and metallicity, showing that such relations are in place before accretion. Our galaxies show a large variety of star formation histories, from extended star formation periods to single bursts. As in more massive galaxies, large star formation bursts are connected with major mergers events, which greatly contribute to the overall stellar mass build up. The intrinsic stochasticity of mergers induces a large scatter in the stellar mass halo mass relation, up to two orders of magnitude. Despite the bursty star formation history, on these mass scales baryons are very ineffective in modifying the dark matter profiles, and galaxies with a stellar mass below $\approx 10^6 M_\odot$ retain their cuspy central dark matter distribution, very similar to results from pure N-body simulations.

The radial profiles of stars in disc galaxies are observed to be either purely exponential (Type-I), truncated (Type-II) or anti-truncated (Type-III) exponentials. Controlled formation simulations of isolated galaxies can reproduce all of these profile types by varying a single parameter, the initial halo spin. In this paper we examine these simulations in more detail in an effort to identify the physical mechanism that leads to the formation of Type-III profiles. The stars in the anti-truncated outskirts of such discs are now on eccentric orbits, but were born on near-circular orbits at much smaller radii. We show that, and explain how, they were driven to the outskirts via non-linear interactions with a strong and long-lived central bar, which greatly boosted their semi-major axis but also their eccentricity. While bars have been known to cause radial heating and outward migration to stellar orbits, we link this effect to the formation of Type-III profiles. This predicts that the anti-truncated parts of galaxies have unusual kinematics for disc-like stellar configurations: high radial velocity dispersions and slow net rotation. Whether such discs exist in nature, can be tested by future observations.

We present here parameter distributions and fundamental scaling relations for 190 galaxies as part of the Spectroscopy and H-bang Imaging of Virgo cluster galaxies (SHIVir) survey. We find the distribution of galaxy velocities to be bimodal about $V_{\rm circ} \sim 125$ km ${\rm s^{-1}}$, hinting at the existence of dynamically unstable modes in the inner regions of galaxies. An analysis of the Tully-Fisher relation (TFR) of late-type galaxies (LTGs) and fundamental plane (FP) of early-type galaxies (ETGs) is also presented, yielding a compendium of galaxy scaling relations. The slope and zero-point of the Virgo TFR match those of field galaxies, while scatter differences likely reflect distinct evolutionary histories. The velocities minimizing scatter for the TFR and FP are measured at large apertures where the baryonic fraction becomes subdominant. While TFR residuals remain independent of any galaxy parameters, FP residuals (i.e. the FP "tilt") correlate strongly with the dynamical-to-stellar mass ratio, yielding stringent galaxy formation constraints. Furthermore, we construct a stellar-to-total mass relation (STMR) for ETGs and LTGs and find linear but distinct trends over the range $M_{*} = 10^{8-11} M_{\odot}$. Stellar-to-halo mass relations (SHMRs), which probe the extended dark matter halo, can be scaled down to masses estimated within the optical radius, showing a tight match with the Virgo STMR at low masses; however, possibly inadequate halo abundance matching prescriptions and broad radial scalings complicate this comparison at all masses. While ETGs appear to be more compact than LTGs of the same stellar mass in projected space, their mass-size relations in physical space are identical. The trends reported here call for validation through well-resolved numerical simulations.

Many massive star forming disc galaxies in the redshift range 3 to 0.5 are observed to have a clumpy morphology showing giant clumps of size $\sim$1 kpc and masses of about $10^7M_{\odot}$ to $10^{10} M_{\odot}$. The nature and fate of these giant clumps is still under debate. In this work we use 19 high-resolution simulations of disc galaxies from the NIHAO sample to study the formation and the evolution of clumps in the discs of high redshift galaxies. We use mock HST - CANDELS observations created with the radiative transfer code GRASIL-3D to carry out, for the first time, a quantitative comparison of the observed fraction of clumpy galaxies and its evolution with redshift with simulations. We find a good agreement between the observed clumpy fraction and the one of the NIHAO galaxies. We find that dust attenuation can suppress intrinsically bright clumps and enhance less luminous ones. In our galaxy sample we only find clumps in light (u-band) from young stars but not in stellar mass surface density maps. This means that the NIHAO sample does not show clumpy stellar discs but rather a clumpy light distribution originating from clumpy star formation events. The clumps found in the NIHAO sample match observed age/color gradients as a function of distance from the galaxy center but they show no sign of inward migration. Clumps in our simulations disperse on timescales of a about a hundred Myr and their contribution to bulge growth is negligible.

We perform Jeans anisotropic modeling (JAM) on elliptical and spiral galaxies from the MaNGA DR13 sample. By comparing the stellar mass-to-light ratios estimated from stellar population synthesis (SPS) and from JAM, we find a similar systematic variation of the initial mass function (IMF) as in the earlier $\rm ATLAS^{3D}$ results. Early type galaxies (elliptical and lenticular) with lower velocity dispersions within one effective radius are consistent with a Chabrier-like IMF while galaxies with higher velocity dispersions are consistent with a more bottom heavy IMF such as the Salpeter IMF. Spiral galaxies have similar systematic IMF variations, but with slightly different slopes and larger scatters, due to the uncertainties caused by higher gas fractions and extinctions for these galaxies. Furthermore, we examine the effects of stellar mass-to-light ratio gradients on our JAM modeling, and find that the trends from our results becomes stronger after considering the gradients.

We present useful functions for the profiles of dark-matter (DM) haloes with a free inner slope, from cusps to cores, where the profiles of density, mass-velocity and potential are simple analytic expressions. Analytic velocity is obtained by expressing the mean density as a simple functional form, and deriving the local density by differentiation. The function involves four shape parameters, with only two or three free: a concentration parameter $c$, inner and outer asymptotic slopes $\alpha$ and $\bar{\gamma}$, and a middle shape parameter $\beta$. Analytic expressions for the potential and velocity dispersion exist for $\bar{\gamma}=3$ and for $\beta$ a natural number. We match the models to the DM haloes in cosmological simulations, with and without baryons, ranging from steep cusps to flat cores. Excellent fits are obtained with three free parameters ($c$, $\alpha$, $\bar{\gamma}$) and $\beta=2$. For an analytic potential, similar fits are obtained for $\bar{\gamma}=3$ and $\beta=2$ with only two free parameters ($c$, $\alpha$); this is our favorite model. A linear combination of two such profiles, with an additional free concentration parameter, provides excellent fits also for $\beta=1$, where the expressions are simpler. The fit quality is comparable to non-analytic popular models. An analytic potential is useful for modeling the inner-halo evolution due to gas inflows and outflows, studying environmental effects on the outer halo, and generating halo potentials or initial conditions for simulations. The analytic velocity can quantify simulated and observed rotation curves without numerical integrations.

We use a sample of 83 high-resolution cosmological zoom-in simulations and a semi-analytic model (SAM) to study the stochasticity of galaxy formation in haloes ranging from dwarf to Milky Way masses. Our simulated galaxies reproduce the observed inefficiency of galaxy formation as expressed through the stellar, gas and baryonic Tully-Fisher relations. For HI velocities in the range (70 < V < 220 km/s), the scatter is just 0.08 to 0.14 dex, consistent with the observed intrinsic scatter at these scales. At low velocities (20 < V < 70 km/s), the simulated scatter is 0.2-0.25 dex, which could be tested with future observations. The scatter in the stellar mass versus dark halo velocity relation is constant for 30 < V < 180 km/s, and smaller (~0.17 dex) when using the maximum circular velocity of the dark matter only simulation, Vmax, compared to the virial velocity ($V_{200}$ or $V_{200}^{DMO}$). The scatter in stellar mass is correlated with halo concentration, and is minimized when using a circular velocity at a fixed fraction of the virial radius $\simeq 0.4 R_{200}$ or with $V_{\alpha}=V_{200}^{DMO} (V_{max}^{DMO} / V_{200}^{DMO})^\alpha$ with \alpha ~ 0.7, consistent with constraints from halo clustering. Using the SAM we show the correlation between halo formation time and concentration is essential in order to reproduce this result. This uniformity in galaxy formation efficiency we see in our hydrodynamical simulations and a semi-analytic model proves the simplicity and self-regulating nature of galaxy formation in a \Lambda Cold Dark Matter universe.

We use the NIHAO galaxy formation simulations to make predictions for the baryonic budget in present day galaxies ranging from dwarf to Milky Way masses. The sample is made of 88 independent high resolution cosmological zoom-in simulations. NIHAO galaxies reproduce key properties of observed galaxies, such as the stellar mass vs halo mass and cold gas vs stellar mass relations. Thus they make plausible predictions for the baryon budget. We present the mass fractions of stars, cold gas ($T<10^4$K), cool gas ($10^4 < T < 10^5$K), warm-hot gas ($10^5 < T < 5\times10^6$K), and hot gas (T$> 5\times10^6$K) inside the virial radius, $R_{200}$. Compared to the predicted baryon mass, using the dark halo mass and the universal baryon fraction, $f_{\rm b}\equiv \Omega_{\rm b}/\Omega_{\rm m}=0.15$, we find that all of our haloes are missing baryons. The missing mass has been relocated past 2 virial radii, and cool gas dominates the corona at low mass (M$_{200} < 3 \times 10^{11} M_\odot$) while the warm-hot gas dominates at high mass (M$_{200} > 3 \times 10^{11} M_\odot$). Haloes of mass $\sim 10^{10} M_\odot$ are missing $\sim 90\%$ of their baryons. More massive haloes ($\sim 10^{12} M_\odot$) retain a higher fraction of their baryons, with $\sim 30\%$ missing, consistent with recent observational estimates. Moreover, these more massive haloes reproduce the observed fraction of cold, warm-hot and hot gas. The fraction of cool gas we predict ($0.11\pm0.06$) is significantly lower than the observation from COS-HALOs (0.3-0.47), but agrees with the alternative analysis of Stern et al. 2016.

The decrease in star formation (SF) and the morphological change necessary to produce the $z=0$ elliptical galaxy population are commonly ascribed to a sudden quenching event, which is able to rid the central galaxy of its cold gas reservoir in a short time. Following this event, the galaxy is able to prevent further SF and stay quiescent via a maintenance mode. We test whether such a quenching event is truly necessary using a simple model of quiescence. In this model, hot gas (all gas above a temperature threshold) in a $\sim10^{12} M_{\odot}$ halo mass galaxy at redshift $z\sim3$ is prevented from cooling. The cool gas continues to form stars at a decreasing rate and the galaxy stellar mass, morphology, velocity dispersion and position on the color magnitude diagram (CMD) proceed to evolve. By $z=0$, the halo mass has grown to $10^{13} M_{\odot}$ and the galaxy has attained characteristics typical of an observed $z=0$ elliptical galaxy. Our model is run in the framework of a cosmological, smooth particle hydrodynamic code which includes SF, early stellar feedback, supernova feedback, metal cooling and metal diffusion. Additionally, we post-process our simulations with a radiative transfer code to create a mock CMD. In contrast to previous assumptions that a pure "fade away" model evolves too slowly to account for the sparsity of galaxies in the "green valley", we demonstrate crossing times of $<$1 Gyr. We conclude that no sudden quenching event is necessary to produce such rapid colour transitions.

We address the origin of Ultra-Diffuse Galaxies (UDGs), which have stellar masses typical of dwarf galaxies but effective radii of Milky Way-sized objects. Their formation mechanism, and whether they are failed $\rm L_{\star}$ galaxies or diffuse dwarfs, are challenging issues. Using zoom-in cosmological simulations from the NIHAO project, we show that UDG analogues form naturally in medium-mass haloes due to episodes of gas outflows associated with star formation. The simulated UDGs live in isolated haloes of masses $10^{10-11}\rm M_{\odot}$, have stellar masses of $10^{7-8.5}\rm M_{\odot}$, effective radii larger than 1 kpc and dark matter cores. They show a broad range of colors, an average S\'ersic index of 0.83, a typical distribution of halo spin and concentration, and a non-negligible HI gas mass of $10^{7-9}\rm M_{\odot}$, which correlates with the extent of the galaxy. Gas availability is crucial to the internal processes that form UDGs: feedback driven gas outflows, and subsequent dark matter and stellar expansion, are the key to reproduce faint, yet unusually extended, galaxies. This scenario implies that UDGs represent a dwarf population of low surface brightness galaxies and should exist in the field. The largest isolated UDGs should contain more HI gas than less extended dwarfs of similar $\rm M_{\star}$.

We study the hot and cold circum-galactic medium (CGM) of 86 galaxies of the cosmological, hydrodynamical simulation suite NIHAO. NIHAO allows a study of how the $z=0$ CGM varies across 5 orders of magnitude of stellar mass using OVI and HI as proxies for hot and cold gas. The cool HI covering fraction and column density profiles match observations well, particularly in the inner CGM. OVI shows increasing column densities with mass, a trend seemingly echoed in the observations. As in multiple previous simulations, the OVI column densities in simulations are lower than observed and optically thick HI does not extend as far out as in observations. We take a look at the collisional ionisation fraction of OVI as a function of halo mass. We make observable predictions of the bipolarity of outflows and their effect on the general shape of the CGM. Bipolar outflows can be seen out to around 40 kpc in intermediate and low mass halos ($M_{\mathrm{halo}}<10^{11}M_{\mathrm{sun}}$), but outside that radius, the CGM is too well mixed to detect an elongated shape. Larger halos have extended gas discs beyond the stellar disc that dominate the shape of the inner CGM. The simulated CGM is remarkably spherical even in low mass simulations. The chemical enrichment of both halo and disc gas follow expected increasing trends as a function of halo mass that are well fit with power laws. These relations can be used in non-hydrodynamic models, such as semi-analytic models.

We present a useful function for describing the profiles of dark-matter haloes with a varying asymptotic inner slope -\alpha, ranging from a cusp to a core, where the profiles of density, mass-velocity and potential are simple analytic expressions for any \alpha. The idea is to express the mean-density profile as a simple functional form, and obtain the local density by derivative. The model involves a concentration parameter c analogous to the NFW profile. More flexibility is provided by a sum of two such functions, with concentrations c_1 and c_2. An optional additional parameter is the asymptotic outer slope -\gamma, which allows more flexibility in the outskirts of haloes. For \gamma different than 3, there are analytic expressions for the profiles of density and mass-velocity but not for the potential. We match the proposed function, in different variants, to the dark-matter profiles of haloes in cosmological simulations, with and without baryons, spanning inner profiles that range from steep cusps to flat cores, and find excellent fits. The analytic potential profile is useful in modeling the evolution of the inner halo due to episodes of gas inflow and outflow. The analytic profile, especially with a free \gamma, is useful for a study of environmental effects in the outer halo. In general, these analytic profiles can serve for straightforwardly quantifying the shapes of simulated and observed rotation curves, without the need for numerical integrations.

We investigate the matter distribution of a spiral galaxy with a counter-rotating stellar core, SDSS J1331+3628 (J1331), independently with gravitational lensing and stellar dynamical modelling. By fitting a gravitational potential model to a quadruplet of lensing images around J1331's bulge, we tightly constrain the mass inside the Einstein radius R_ein = (0.91 +/- 0.02)'' (~= 1.83 +/- 0.04 kpc) to within 4%: M_ein = (7.8 +/- 0.3) x 10^10 M_Sun. We model observed long-slit major axis stellar kinematics in J1331's central regions by finding Multi-Gaussian Expansion (MGE) models for the stellar and dark matter distribution that solve the axisymmetric Jeans equations. The lens and dynamical model are independently derived, but in very good agreement with each other around ~R_ein. We find that J1331's center requires a steep total mass-to-light ratio gradient. A dynamical model including an NFW halo (with virial velocity v_200 ~= 240 +/- 40 km/s and concentration c_200 ~= 8 +/- 2) and moderate tangential velocity anisotropy (beta_z ~= -0.4 +/- 0.1) can reproduce the signatures of J1331's counter-rotating core and predict the stellar and gas rotation curve at larger radii. However, our models do not agree with the observed velocity dispersion at large radii. We speculate that the reason could be a non-trivial change in structure and kinematics due to a possible merger event in J1331's recent past.

We used 87 high resolution hydrodynamical cosmological simulations from the NIHAO suite to investigate the relation between the maximum circular velocity (Vmax) of a dark matter halo in a collisionless simulation and the velocity width of the HI gas in the same halo in the hydrodynamical simulation. These two quantities are normally used to compare theoretical and observational velocity functions and have led to a possible discrepancy between observations and predictions based on the Cold Dark Matter (CDM) model. We show that below 100 km/s, there is clear bias between HI based velocities and Vmax, that leads to an underestimation of the actual circular velocity of the halo. When this bias is taken into account the CDM model has no trouble in reproducing the observed velocity function and no lack of low velocity galaxies is actually present. Our simulations also reproduce the linewidth - stellar mass (Tully-Fisher) relation and HI sizes, indicating that the HI gas in our simulations is as extended as observed. The physical reason for the lower than expected linewidths is that, in contrast to high mass galaxies, low mass galaxies no longer have extended thin HI rotating disks, as is commonly assumed.

We use the NIHAO (Numerical Investigation of Hundred Astrophysical Objects) cosmological simulations to study the effects of galaxy formation on key properties of dark matter (DM) haloes. NIHAO consists of $\simeq 90$ high-resolution SPH simulations that include (metal-line) cooling, star formation, and feedback from massive stars and SuperNovae, and cover a wide stellar and halo mass range: $10^6 < M_* / M_{\odot} < 10^{11}$ ( $10^{9.5} < M_{\rm halo} / M_{\odot} < 10^{12.5}$). When compared to DM-only simulations, the NIHAO haloes have similar shapes at the virial radius, R_{\rm vir}, but are substantially rounder inside $\simeq 0.1R_{\rm vir}$. In NIHAO simulations $c/a$ increases with halo mass and integrated star formation efficiency, reaching $\sim 0.8$ at the Milky Way mass (compared to 0.5 in DM-only), providing a plausible solution to the long-standing conflict between observations and DM-only simulations. The radial profile of the phase-space $Q$ parameter ($\rho/\sigma^3$) is best fit with a single power law in DM-only simulations, but shows a flattening within $\simeq 0.1R_{\rm vir}$ for NIHAO for total masses $M>10^{11} M_{\odot}$. Finally, the global velocity distribution of DM is similar in both DM-only and NIHAO simulations, but in the solar neighborhood, NIHAO galaxies deviate substantially from Maxwellian. The distribution is more symmetric, roughly Gaussian, with a peak that shifts to higher velocities for Milky Way mass haloes. We provide the distribution parameters which can be used for predictions for direct DM detection experiments. Our results underline the ability of the galaxy formation processes to modify the properties of dark matter haloes.

We use ~100 cosmological galaxy formation zoom-in simulations using the smoothed particle hydrodynamics code {\sc gasoline} to study the effect of baryonic processes on the mass profiles of cold dark matter haloes. The haloes in our study range from dwarf (M_{200}~10^{10}Msun) to Milky Way (M_{200}~10^{12}Msun) masses. Our simulations exhibit a wide range of halo responses, primarily varying with mass, from expansion to contraction, with up to factor ~10 changes in the enclosed dark matter mass at one per cent of the virial radius. Confirming previous studies, the halo response is correlated with the integrated efficiency of star formation: e_SF=(M_{star}/M_{200})/(\Omega_b/\Omega_m). In addition we report a new correlation with the compactness of the stellar system: e_R=r_{1/2}/R_{200}. We provide an analytic formula depending on e_SF and e_R for the response of cold dark matter haloes to baryonic processes. An observationally testable prediction is that, at fixed mass, larger galaxies experience more halo expansion, while the smaller galaxies more halo contraction. This diversity of dark halo response is captured by a toy model consisting of cycles of adiabatic inflow (causing contraction) and impulsive gas outflow (causing expansion). For net outflow, or equal inflow and outflow fractions, f, the overall effect is expansion, with more expansion with larger f. For net inflow, contraction occurs for small f (large radii), while expansion occurs for large f (small radii), recovering the phenomenology seen in our simulations. These regularities in the galaxy formation process provide a step towards a fully predictive model for the structure of cold dark matter haloes.

A large fraction of the dwarf satellites orbiting the Andromeda galaxy are surprisingly aligned in a thin, extended and apparently kinematically coherent planar structure. Such a structure is not easily found in simulations based on the Cold Dark Matter model ({\Lambda}CDM). Using 21 high resolution cosmological simulations we analyse the kinematics of planes of satellites similar to the one around Andromeda. We find good agreement when co-rotation is characterized by the line-of-sight velocity. At the same time, when co-rotation is inferred by the angular momenta of the satellites, the planes are in agreement with the plane around our Galaxy. We find such planes to be common in our high concentration haloes. The number of co-rotating satellites obtained from the sign of the line-of-sight velocity shows large variations depending on the viewing angle and is consistent with that obtained from a sample with random velocities. We find that the clustering of angular momentum vectors of the satellites in the plane is a better measure of the kinematic coherence. Thus we conclude that the line-of- sight velocity is not well suited as a proxy for the kinematical coherence of the plane. Analysis of the kinematics of our planes shows a fraction of $\sim$30% chance aligned satellites. Tracking the satellites in the plane back in time reveals that these planes are a transient feature and not kinematically coherent as would appear at first sight. Thus we expect some of the satellites in the plane around Andromeda to have high velocities perpendicular to the plane.

Detailed studies of galaxy formation require clear definitions of the structural components of galaxies. Precisely defined components also enable better comparisons between observations and simulations. We use a subsample of eighteen cosmological zoom-in simulations from the NIHAO project to derive a robust method for defining stellar kinematic discs in galaxies. Our method uses Gaussian Mixture Models in a 3D space of dynamical variables. The NIHAO galaxies have the right stellar mass for their halo mass, and their angular momenta and S\'ersic indices match observations. While the photometric disc-to-total ratios are close to 1 for all the simulated galaxies, the kinematic ratios are around ~0.5. Thus, exponential structure does not imply a cold kinematic disc. Above log(M*)~9.5, the decomposition leads to thin discs and spheroids that have clearly different properties, in terms of angular momentum, rotational support, ellipticity, [Fe/H] and [O/Fe]. At log(M*)<9.5, the decomposition selects discs and spheroids with less distinct properties. At these low masses, both the discs and spheroids have exponential profiles with high minor-to-major axes ratios, i.e. thickened discs.

We compare the half-light circular velocities, V_{1/2}, of dwarf galaxies in the Local Group to the predicted circular velocity curves of galaxies in the NIHAO suite of LCDM simulations. We use a subset of 34 simulations in which the central galaxy has a stellar luminosity in the range 0.5 x 10^5 < L_V < 2 x 10^8 L_{sun}. The NIHAO galaxy simulations reproduce the relation between stellar mass and halo mass from abundance matching, as well as the observed half-light size vs luminosity relation. The corresponding dissipationless simulations over-predict the V_{1/2}, recovering the problem known as too big to fail (TBTF). By contrast, the NIHAO simulations have expanded dark matter haloes, and provide an excellent match to the distribution of V_{1/2} for galaxies with L_V > 2 x 10^6 L_{sun}. For lower luminosities our simulations predict very little halo response, and tend to over predict the observed circular velocities. In the context of LCDM, this could signal the increased stochasticity of star formation in haloes below M_{halo} \sim 10^{10} M_{sun}, or the role of environmental effects. Thus, haloes that are "too big to fail", do not fail LCDM, but haloes that are "too small to pass" (the galaxy formation threshold) provide a future test of LCDM.

We investigate the compact, early-type galaxy NGC 1281 with integral field unit observations to map the stellar LOSVD out to 5 effective radii and construct orbit-based dynamical models to constrain its dark and luminous matter content. Under the assumption of mass-follows-light, the H-band stellar mass-to-light ratio (M/L) is {\Upsilon} = 2.7(+-0.1) {\Upsilon}_{sun}, higher than expected from our stellar population synthesis fits with either a canonical Kroupa ({\Upsilon} = 1.3 {\Upsilon}_{sun}) or Salpeter ({\Upsilon} = 1.7 {\Upsilon}_{sun}) stellar initial mass function. Such models also cannot reproduce the details of the LOSVD. Models with a dark halo recover the kinematics well and indicate that NGC 1281 is dark matter dominated, making up ~ 90 per cent of the total enclosed mass within the kinematic bounds. Parameterised as a spherical NFW profile, the dark halo mass is 11.5 < log(M_{DM}/M_{sun}) < 11.8 and the stellar M/L is 0.6 < {\Upsilon} < 1.1. However, this stellar M/L is lower than predicted by its old stellar population. Moreover, the halo mass within the kinematic extent is ten times larger than expected based on {\Lambda}CDM predictions, and an extrapolation yields cluster sized dark halo masses. Adopting {\Upsilon} = 1.7 {\Upsilon}_{sun} yields more moderate dark halo virial masses, but these models fit the kinematics worse. A non-NFW model might solve the discrepancy between the unphysical consequences of the best-fitting dynamical models and models based on more reasonable assumptions for the dark halo and stellar mass-to-light ratio, which are disfavoured according to our parameter estimation.

We introduce project NIHAO (Numerical Investigation of a Hundred Astrophysical Objects), a set of 100 cosmological zoom-in hydrodynamical simulations performed using the gasoline code, with an improved implementation of the SPH algorithm. The haloes in our study range from dwarf to Milky Way masses, and represent an unbiased sampling of merger histories, concentrations and spin parameters. The particle masses and force softenings are chosen to resolve the mass profile to below 1% of the virial radius at all masses, ensuring that galaxy half-light radii are well resolved. Using the same treatment of star formation and stellar feedback for every object, the simulated galaxies reproduce the observed inefficiency of galaxy formation across cosmic time as expressed through the stellar mass vs halo mass relation, and the star formation rate vs stellar mass relation. We thus conclude that stellar feedback is the chief piece of physics required to limit the efficiency of star formation in galaxies less massive than the Milky Way.

We use cosmological hydrodynamical zoom-in simulations with the SPH code gasoline of four haloes of mass M_{200} \sim 10^{13}\Msun to study the response of the dark matter to elliptical galaxy formation. Our simulations include metallicity dependent gas cooling, star formation, and feedback from massive stars and supernovae, but not active galactic nuclei (AGN). At z=2 the progenitor galaxies have stellar to halo mass ratios consistent with halo abundance matching, assuming a Salpeter initial mass function. However by z=0 the standard runs suffer from the well known overcooling problem, overpredicting the stellar masses by a factor of > 4. To mimic a suppressive halo quenching scenario, in our forced quenching (FQ) simulations, cooling and star formation are switched off at z=2. The resulting z=0 galaxies have stellar masses, sizes and circular velocities close to what is observed. Relative to the control simulations, the dark matter haloes in the FQ simulations have contracted, with central dark matter density slopes d\log\rho/d\log r \sim -1.5, showing that dry merging alone is unable to fully reverse the contraction that occurs at z>2. Simulations in the literature with AGN feedback however, have found expansion or no net change in the dark matter halo. Thus the response of the dark matter halo to galaxy formation may provide a new test to distinguish between ejective and suppressive quenching mechanisms.

Recently it has been shown that a large fraction of the dwarf satellite galaxies orbiting the Andromeda galaxy are surprisingly aligned in a thin, extended and kinematically coherent planar structure. The presence of such a structure seems to challenge the current Cold Dark Matter paradigm of structure formation, which predicts a more uniform distribution of satellites around central objects. We show that it is possible to obtain a thin, extended, rotating plane of satellites resembling the one in Andromeda in cosmological collisionless simulations based on the Cold Dark Matter model. Our new high resolution simulations show a correlation between the formation time of the dark matter halo and the thickness of the plane of satellites. Our simulations have a high incidence of satellite planes as thin, extended, and as rich as the one in Andromeda and with a very coherent kinematic structure when we select high concentration/early forming halos. By tracking the formation of the satellites in the plane we show that they have been mainly accreted onto the main object along thin dark matter filaments at high redshift. Our results show that the presence of a thin, extended, rotating plane of satellites is not a challenge for the Cold Dark Matter paradigm, but actually supports one of the predictions of this paradigm related to the presence of filaments of dark matter around galaxies at high redshift.

We show that the cool gas masses of galactic discs reach a steady state that lasts many Gyr after their last major merger in cosmological hydrodynamic simulations. The mass of disc gas, M$_{\rm gas}$, depends upon a galaxy halo's spin and virial mass, but not upon stellar feedback. Halos with low spin have high star formation efficiency and lower disc gas mass. Similarly, lower stellar feedback leads to more star formation so the gas mass ends up nearly the same irregardless of stellar feedback strength. Even considering spin, the M$_{\rm gas}$ relation with halo mass, M$_{200}$ only shows a factor of 3 scatter. The M$_{\rm gas}$--M$_{200}$ relation show a break at M$_{200}$=$2\times10^{11}$ M$_\odot$ that corresponds to an observed break in the M$_{\rm gas}$--M$_\star$ relation. The constant disc mass stems from a shared halo gas density profile in all the simulated galaxies. In their outer regions, the profiles are isothermal. Where the profile rises above $n=10^{-3}$ cm$^{-3}$, the gas readily cools and the profile steepens. Inside the disc, rotation supports gas with a flatter density profile except where supernova explosions disrupt the disc. Energy injection from stellar feedback also provides pressure support to the halo gas to prevent runaway cooling flows. The resulting constant gas mass makes simpler models for galaxy formation possible, either using a "bathtub" model for star formation rates or when modeling chemical evolution.

The radial density profiles of stellar galaxy discs can be well approximated as an exponential. Compared to this canonical form, however, the profiles in the majority of disc galaxies show downward or upward breaks at large radii. Currently, there is no coherent explanation in a galaxy formation context of the radial profile per se, along with the two types of profile breaks. Using a set of controlled hydrodynamic simulations of disc galaxy formation, we find a correlation between the host halo's initial angular momentum and the resulting radial profile of the stellar disc: galaxies that live in haloes with a low spin parameter {\lambda} <~ 0.03 show an up-bending break in their disc density profiles, while galaxies in haloes of higher angular momentum show a down-bending break. We find that the case of pure exponential profiles ({\lambda} ~ 0.035) coincides with the peak of the spin parameter distribution from cosmological simulations. Our simulations not only imply an explanation of the observed behaviours, but also suggest that the physical origin of this effect is related to the amount of radial redistribution of stellar mass, which is anti-correlated with {\lambda}.

We present a summary of recent assessments of the mass distribution in disk galaxies. Of issue in order to characterize galaxy formation models is to determine the relative fraction of baryons and dark matter at all radii in galaxies. For disk galaxies, various measurements of the mass distribution in galaxies based on vertical kinematics, strong lensing, residuals of scaling relations, fluid dynamical modeling, bar strength and pattern speed, warps, and others, have called for either a maximal or sub-maximal contribution of the baryons in the inner parts of the disk. We propose a global picture whereby all galaxies are typically baryon-dominated (maximal) at the center and dark-matter dominated (sub-maximal) in their outskirts. Using as a fiducial radius the peak of the rotation curve of a pure baryonic exponential disk at R_2.2=2.2 Rd, where Rd is the disk scale length, the transition from maximal to sub-maximal baryons occurs within or near R_2.2 for low-mass disk galaxies (with V_tot < 200 km/s) and beyond R_2.2 for more massive systems. The mass fraction's inverse dependence on circular velocity in disk galaxies is largely insensitive to the presence of a bar at R_2.2. The mean mass fractions of late- and early-type galaxies are shown to be different at the same fiducial radius and circular velocity, pointing to different galaxy formation mechanisms. Feedback, dynamical friction, size, and IMF variations are key ingredients for understanding these differences.

We present an overview of a new integral field spectroscopic survey called MaNGA (Mapping Nearby Galaxies at Apache Point Observatory), one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) that began on 2014 July 1. MaNGA will investigate the internal kinematic structure and composition of gas and stars in an unprecedented sample of 10,000 nearby galaxies. We summarize essential characteristics of the instrument and survey design in the context of MaNGA's key science goals and present prototype observations to demonstrate MaNGA's scientific potential. MaNGA employs dithered observations with 17 fiber-bundle integral field units that vary in diameter from 12" (19 fibers) to 32" (127 fibers). Two dual-channel spectrographs provide simultaneous wavelength coverage over 3600-10300 A at R~2000. With a typical integration time of 3 hr, MaNGA reaches a target r-band signal-to-noise ratio of 4-8 (per A, per 2" fiber) at 23 AB mag per sq. arcsec, which is typical for the outskirts of MaNGA galaxies. Targets are selected with stellar mass greater than 1e9 Msun using SDSS-I redshifts and i-band luminosity to achieve uniform radial coverage in terms of the effective radius, an approximately flat distribution in stellar mass, and a sample spanning a wide range of environments. Analysis of our prototype observations demonstrates MaNGA's ability to probe gas ionization, shed light on recent star formation and quenching, enable dynamical modeling, decompose constituent components, and map the composition of stellar populations. MaNGA's spatially resolved spectra will enable an unprecedented study of the astrophysics of nearby galaxies in the coming 6 yr.

Galaxies are not randomly distributed across the universe but showing different kinds of alignment on different scales. On small scales satellite galaxies have a tendency to distribute along the major axis of the central galaxy, with dependence on galaxy properties that both red satellites and centrals have stronger alignment than their blue counterparts. On large scales, it is found that the major axes of Luminous Red Galaxies (LRGs) have correlation up to 30Mpc/h. Using hydro-dynamical simulation with star formation, we investigate the origin of galaxy alignment on different scales. It is found that most red satellite galaxies stay in the inner region of dark matter halo inside which the shape of central galaxy is well aligned with the dark matter distribution. Red centrals have stronger alignment than blue ones as they live in massive haloes and the central galaxy-halo alignment increases with halo mass. On large scales, the alignment of LRGs is also from the galaxy-halo shape correlation, but with some extent of mis-alignment. The massive haloes have stronger alignment than haloes in filament which connect massive haloes. This is contrary to the naive expectation that cosmic filament is the cause of halo alignment.

Observations have shown that the spatial distribution of satellite galaxies is not random, but rather is aligned with the major axes of central galaxies (CGs). The strength of the alignment is dependent on the properties of both the satellites and centrals. Theoretical studies using dissipationless N-body simulations are limited by their inability to directly predict the shape of CGs. Using hydrodynamical simulations including gas cooling, star formation, and feedback, we carry out a study of galaxy alignment and its dependence on the galaxy properties predicted directly from the simulations.We found that the observed alignment signal is well produced, as is the color dependence: red satellites and red centrals both show stronger alignments than their blue counterparts. The reason for the stronger alignment of red satellites is that most of them stay in the inner region of the dark matter halo where the shape of the CG better traces the dark matter distribution. The dependence of alignment on the color of CGs arises from the halo mass dependence, since the alignment between the shape of the central stellar component and the inner halo increases with halo mass. We also find that the alignment of satellites is most strongly dependent on their metallicity, suggesting that the metallicity of satellites, rather than color, is a better tracer of galaxy alignment on small scales. This could be tested in future observational studies.

We present the evolution of the structure of relaxed cold dark matter haloes in the cosmology from the Planck satellite. Our simulations cover 5 decades in halo mass, from dwarf galaxies to galaxy clusters. Due to the increased matter density and power spectrum normalization the concentration mass relation in the Planck cosmology has a 20 percent higher normalization at redshift z=0 compared to WMAP cosmology. We confirm that CDM haloes are better described by the Einasto profile; for example, at scales near galaxy half-light radii CDM haloes have significantly steeper density profiles than implied by NFW fits. There is a scatter of 0.2 dex in the Einasto shape parameter at fixed halo mass, adding further to the diversity of CDM halo profiles. The evolution of the concentration mass relation in our simulations is not reproduced by any of the analytic models in the literature. We thus provide a simple fitting formula that accurately describes the evolution between redshifts z=5 to z=0 for both NFW and Einasto fits. Finally, the observed concentrations and halo masses of spiral galaxies, groups and clusters of galaxies at low redshifts are in good agreement with our simulations, suggesting only mild halo response to galaxy formation on these scales.

Galaxy masses play a fundamental role in our understanding of structure formation models. This review addresses the variety and reliability of mass estimators that pertain to stars, gas, and dark matter. The different sections on masses from stellar populations, dynamical masses of gas-rich and gas-poor galaxies, with some attention paid to our Milky Way, and masses from weak and strong lensing methods, all provide review material on galaxy masses in a self-consistent manner.

We introduce a mass dependent density profile to describe the distribution of dark matter within galaxies, which takes into account the stellar-to-halo mass dependence of the response of dark matter to baryonic processes. The study is based on the analysis of hydrodynamically simulated galaxies from dwarf to Milky Way mass, drawn from the MaGICC project, which have been shown to match a wide range of disk scaling relationships. We find that the best fit parameters of a generic double power-law density profile vary in a systematic manner that depends on the stellar-to-halo mass ratio of each galaxy. Thus, the quantity Mstar/Mhalo constrains the inner ($\gamma$) and outer ($\beta$) slopes of dark matter density, and the sharpness of transition between the slopes($\alpha$), reducing the number of free parameters of the model to two. Due to the tight relation between stellar mass and halo mass, either of these quantities is sufficient to describe the dark matter halo profile including the effects of baryons. The concentration of the haloes in the hydrodynamical simulations is consistent with N-body expectations up to Milky Way mass galaxies, at which mass the haloes become twice as concentrated as compared with pure dark matter runs. This mass dependent density profile can be directly applied to rotation curve data of observed galaxies and to semi analytic galaxy formation models as a significant improvement over the commonly used NFW profile.

Observations have shown that the spatial distribution of satellite galaxies is not random, but aligned with the major axes of central galaxies. This alignment is dependent on galaxy properties, such that red satellites are more strongly aligned than blue satellites. Theoretical work done to interpret this phenomena has found that it is due to the non-spherical nature of dark matter halos. However, most studies over-predict the alignment signal under the assumption that the central galaxy shape follows the shape of the host halo. It is also not clear whether the color dependence of alignment is due to an assembly bias or an evolution effect. In this paper we study these problems using a cosmological N-body simulation. Subhalos are used to trace the positions of satellite galaxies. It is found that the shape of dark matter halos are mis-aligned at different radii. If the central galaxy shares the same shape as the inner host halo, then the alignment effect is weaker and agrees with observational data. However, it predicts almost no dependence of alignment on the color of satellite galaxies, though the late accreted subhalos show stronger alignment with the outer layer of the host halo than their early accreted counterparts. We find that this is due to the limitation of pure N-body simulations that satellites galaxies without associated subhalos ('orphan galaxies') are not resolved. These orphan (mostly red) satellites often reside in the inner region of host halos and should follow the shape of the host halo in the inner region.

Recent studies have shown that massive elliptical galaxies have total mass density profiles within an effective radius that can be approximated as \rho_{tot}\propto r^{-\gamma'}, with mean slope <\gamma'>=2.08 \pm 0.03 and scatter \sigma_\gamma'=0.16 \pm 0.02. The small scatter of the slope (known as the bulge-halo conspiracy) is not generic in LCDM based models and therefore contains information about the galaxy formation process. We compute the distribution of \gamma' for LCDM-based models that reproduce the observed correlations between stellar mass, velocity dispersion, and effective radius of early-type galaxies in the SDSS. The models have a range of stellar initial mass functions (IMFs) and dark halo responses to galaxy formation. The observed distribution of \gamma' is well reproduced by a model with cosmologically motivated but uncontracted dark matter haloes, and a Salpeter-type IMF. Other models are on average ruled out by the data, even though they may happen in individual cases. Models with adiabatic halo contraction (and lighter IMFs) predict too small values of \gamma'. Models with halo expansion, or mass-follows-light predict too high values of \gamma'. Our study shows that the non-homologous structure of massive early-type galaxies can be precisely reproduced by LCDM models if the IMF is not universal and if mechanisms such as feedback from active galactic nuclei, or dynamical friction, effectively on average counterbalance the contraction of the halo expected as a result of baryonic cooling.

The long-standing assumption that the stellar initial mass function (IMF) is universal has recently been challenged by a number of observations. Several studies have shown that a "heavy" IMF (e.g., with a Salpeter-like abundance of low mass stars and thus normalisation) is preferred for massive early-type galaxies, while this IMF is inconsistent with the properties of less massive, later-type galaxies. These discoveries motivate the hypothesis that the IMF may vary (possibly very slightly) across galaxies and across components of individual galaxies (e.g. bulges vs discs). In this paper we use a sample of 19 late-type strong gravitational lenses from the SWELLS survey to investigate the IMFs of the bulges and discs in late-type galaxies. We perform a joint analysis of the galaxies' total masses (constrained by strong gravitational lensing) and stellar masses (constrained by optical and near-infrared colours in the context of a stellar population synthesis [SPS] model, up to an IMF normalisation parameter). Using minimal assumptions apart from the physical constraint that the total stellar mass within any aperture must be less than the total mass within the aperture, we find that the bulges of the galaxies cannot have IMFs heavier (i.e. implying high mass per unit luminosity) than Salpeter, while the disc IMFs are not well constrained by this data set. We also discuss the necessity for hierarchical modelling when combining incomplete information about multiple astronomical objects. This modelling approach allows us to place upper limits on the size of any departures from universality. More data, including spatially resolved kinematics (as in paper V) and stellar population diagnostics over a range of bulge and disc masses, are needed to robustly quantify how the IMF varies within galaxies.

We use 31 simulated galaxies from the MaGICC project to investigate the effects of baryonic feedback on the density profiles of dark matter (DM) haloes. The sample covers a wide mass range: 9.4e9<Mhalo/Msun<7.8e11, hosting galaxies with stellar masses: 5.0e5<M*/Msun<8.3e10, i.e. from dwarf to L*. The galaxies are simulated with several baryonic prescriptions, including a range of stellar feedbacks. The main result is a clear dependence of the inner slope of the DM density profile, \alpha\ in \rho r^\alpha, on the ratio between stellar-to-halo mass (M*/Mhalo). This relation is independent of the stellar feedback scheme, allowing a prediction for cusp vs core formation. When M*/Mhalo is low, ~0.01%, energy from stellar feedback is insufficient to significantly alter the inner DM density and the galaxy retains a cuspy profile. At higher M*/Mhalo, feedback drives the expansion of the DM and generates cored profiles. The flattest profiles form where M*/Mhalo~0.5%. Above this ratio, stars formed in the central regions deepen the gravitational potential enough to oppose this supernova-driven expansion process, resulting in smaller cores and cuspier profiles. Combining the dependence of \alpha\ on M*/Mhalo with the abundance matching relation between M* and Mhalo provides a prediction for how \alpha\ varies with M*. Further, using the Tully-Fisher relation allows a prediction for the dependence of the DM inner slope on the observed rotation velocity of galaxies. The most cored galaxies are expected to have Vrot~50km/s, with \alpha\ decreasing for more massive disc galaxies: spirals with Vrot~150km/s have central slopes \alpha<-0.8, approaching the NFW profile. This novel prediction for the dependence of \alpha\ on disc galaxy mass can be tested using current observational data sets, and can be applied to theoretical modeling of mass profiles and populations of disc galaxies.

Several independent lines of evidence suggest that the stellar initial mass function (IMF) in early-type galaxies becomes increasingly `bottom-heavy' with increasing galaxy mass and/or velocity dispersion, sigma. Here we consider evidence for IMF variation in a sample of relatively compact early-type galaxies drawn from the Sloan Digital Sky Survey. These galaxies are of sufficiently high stellar density that a dark halo likely makes a minor contribution to the total dynamical mass, Mdyn, within one effective radius. We fit our detailed stellar population synthesis models to the stacked absorption line spectra of these galaxies in bins of sigma and find evidence from IMF-sensitive spectral features for a bottom-heavy IMF at high sigma. We also apply simple `mass-follows-light' dynamical models to the same data and find that Mdyn is significantly higher than what would be expected if these galaxies were stellar dominated and had a universal Milky Way IMF. Adopting Mdyn M_* therefore implies that the IMF is `heavier' at high sigma. Most importantly, the quantitative amount of inferred IMF variation is very similar between the two techniques, agreeing to within 0.1 dex in mass. The agreement between two independent techniques, when applied to the same data, provides compelling evidence for systematic variation in the IMF as a function of early-type galaxy velocity dispersion. Any alternative explanations must reproduce both the results from dynamical and stellar population-based techniques.

We use the relations between aperture stellar velocity dispersion (\sigma_ap), stellar mass (M_sps), and galaxy size (R_e) for a sample of \sim 150,000 early-type galaxies from SDSS/DR7 to place constraints on the stellar initial mass function (IMF) and dark halo response to galaxy formation. We build LCDM based mass models that reproduce, by construction, the relations between galaxy size, light concentration and stellar mass, and use the spherical Jeans equations to predict \sigma_ap. Given our model assumptions (including those in the stellar population synthesis models), we find that reproducing the median \sigma_ap vs M_sps relation is not possible with {\it both} a universal IMF and a universal dark halo response. Significant departures from a universal IMF and/or dark halo response are required, but there is a degeneracy between these two solutions. We show that this degeneracy can be broken using the strength of the correlation between residuals of the velocity-mass (\Delta log \sigma_ap) and size-mass (\Delta log R_e) relations. The slope of this correlation, d_vr \equiv \Delta log \sigma_ap/\Delta log R_e, varies systematically with galaxy mass from d_vr \simeq -0.45 at M_sps \sim 10^{10}M_sun, to d_vr \simeq -0.15 at M_sps \sim 10^{11.6} M_sun. The virial fundamental plane (FP) has d_vr=-1/2, and thus we find the tilt of the observed FP is mass dependent. Reproducing this tilt requires {\it both} a non-universal IMF and a non-universal halo response. Our best model has mass-follows-light at low masses (Msps < 10^{11.2}M_sun) and unmodified NFW haloes at M_sps \sim 10^{11.5} M_sun. The stellar masses imply a mass dependent IMF which is "lighter" than Salpeter at low masses and "heavier" than Salpeter at high masses.

In this paper we combine high resolution N-body simulations with a semi analytical model of galaxy formation to study the effects of a possible Warm Dark Matter (WDM) component on the observable properties of galaxies. We compare three WDM models with a dark matter mass of 0.5, 0.75 and 2.0 keV, with the standard Cold Dark Matter case. For a fixed set of parameters describing the baryonic physics the WDM models predict less galaxies at low (stellar) masses, as expected due to the suppression of power on small scales, while no substantial difference is found at the high mass end. However these differences in the stellar mass function, vanish when different set of parameters are used to describe the (largely unknown) galaxy formation processes. We show that is possible to break this degeneracy between DM properties and the parameterization of baryonic physics by combining observations on the stellar mass function with the Tully-Fisher relation (the relation between stellar mass and the rotation velocity at large galactic radii as probed by resolved HI rotation curves). WDM models with a too warm candidate (m<0.75 keV) cannot simultaneously reproduce the stellar mass function and the Tully-Fisher relation. We conclude that accurate measurements of the galaxy stellar mass function and the link between galaxies and dark matter haloes down to the very low-mass end can give very tight constraints on the nature of DM candidates.

The shutdown of star formation in galaxies is generally termed `quenching'. Although quenching may occur through a variety of processes, the exact mechanism(s) that is in fact responsible for quenching is still in question. This paper addresses quenching by searching for traces of possible quenching processes through their effects on galaxy structural parameters such as surface stellar mass density and Sersic index (n). We analyze the rest-frame U-B color correlations versus these structural parameters using a sample of galaxies in the redshift range 0.5< z<0.8 from the DEEP2/AEGIS survey. We find that Sersic index (n) has the smallest overlap region among all tested parameters and resembles a step-function with a threshold value of n=2.3. There exists, however, a significant population of outliers with blue colors yet high n values that seem to contradict this behavior. We hypothesize that their Sersic values may be distorted by bursts of star formation, AGNs, and/or poor fits, leading us to consider central surface stellar mass density as an alternative to Sersic index. Not only does it correct the outliers, it also forms a tight relationship with color, suggesting that the innermost structure of galaxies is most physically linked with quenching. Furthermore, at z~0.65, the majority of the blue cloud galaxies cannot simply fade onto the red sequence since their GIM2D bulge masses are only half as large on average as the bulge masses of similar red sequence galaxies, thus demonstrating that stellar mass must absolutely increase at the centers of galaxies as they quench. We discuss a two-stage model for quenching in which galaxy star formation rates are controlled by their dark halos while they are still in the blue cloud and a second quenching process sets in later, associated with the central stellar mass build-up.

We study the dependence of star-formation quenching on galaxy mass and environment, in the SDSS (z~0.1) and the AEGIS (z~1). It is crucial that we define quenching by low star-formation rate rather than by red colour, given that one third of the red galaxies are star forming. We address stellar mass M*, halo mass Mh, density over the nearest N neighbours deltaN, and distance to the halo centre D. The fraction of quenched galaxies appears more strongly correlated with Mh at fixed M* than with M* at fixed Mh, while for satellites quenching also depends on D. We present the M*-Mh relation for centrals at z~1. At z~1, the dependence of quenching on M* at fixed Mh is somewhat more pronounced than at z~0, but the quenched fraction is low (10%) and the haloes are less massive. For satellites, M*-dependent quenching is noticeable at high D, suggesting a quenching dependence on sub-halo mass for recently captured satellites. At small D, where satellites likely fell in more than a few Gyr ago, quenching strongly depends on Mh, and not on M*. The Mh-dependence of quenching is consistent with theoretical wisdom where virial shock heating in massive haloes shuts down accretion and triggers ram-pressure stripping, causing quenching. The interpretation of deltaN is complicated by the fact that it depends on the number of observed group members compared to N, motivating the use of D as a better measure of local environment.

Recent work has suggested that the stellar initial mass function (IMF) is not universal, but rather is correlated with galaxy stellar mass, stellar velocity dispersion, or morphological type. In this paper, we investigate variations of the IMF within individual galaxies. For this purpose, we use strong lensing and gas kinematics to measure independently the normalisation of the IMF of the bulge and disk components of a sample of 5 massive spiral galaxies with substantial bulge components taken from the SWELLS survey. We find that the stellar mass of the bulges are tightly constrained by the lensing and kinematic data. A comparison with masses based on stellar population synthesis models fitted to optical and near infrared photometry favors a Salpeter-like normalisation of the IMF. Conversely, the disk masses are less well constrained due to degeneracies with the dark matter halo, but are consistent with Milky Way type IMFs in agreement with previous studies. The disks are submaximal at 2.2 disk scale lengths, but due to the contribution of the bulges, the galaxies are baryon dominated at 2.2 disk scale lengths. Globally, our inferred IMF normalisation is consistent with that found for early-type galaxies of comparable stellar mass (> 10^11 M_sun). Our results suggest a non-universal IMF within the different components of spiral galaxies, adding to the well-known differences in stellar populations between disks and bulges.

We present evidence from a sample of 544 galaxies from the DEEP2 Survey for evolution of the internal kinematics of blue galaxies with stellar masses ranging 8.0 < log M* (M_Sun) < 10.7 over 0.2<z<1.2. DEEP2 provides galaxy spectra and Hubble imaging from which we measure emission-line kinematics and galaxy inclinations, respectively. Our large sample allows us to overcome scatter intrinsic to galaxy properties in order to examine trends in kinematics. We find that at a fixed stellar mass galaxies systematically decrease in disordered motions and increase in rotation velocity and potential well depth with time. Massive galaxies are the most well-ordered at all times examined, with higher rotation velocities and less disordered motions than less massive galaxies. We quantify disordered motions with an integrated gas velocity dispersion corrected for beam smearing (sigma_g). It is unlike the typical pressure-supported velocity dispersion measured for early type galaxies and galaxy bulges. Because both seeing and the width of our spectral slits comprise a significant fraction of the galaxy sizes, sigma_g integrates over velocity gradients on large scales which can correspond to non-ordered gas kinematics. We compile measurements of galaxy kinematics from the literature over 1.2<z<3.8 and do not find any trends with redshift, likely for the most part because these datasets are biased toward the most highly star-forming systems. In summary, over the last ~8 billion years since z=1.2, blue galaxies evolve from disordered to ordered systems as they settle to become the rotation-dominated disk galaxies observed in the Universe today, with the most massive galaxies being the most evolved at any time.

Most of the baryons in the Universe are not in the form of stars and cold gas in galaxies. Galactic outflows driven by supernovae/stellar winds are the leading mechanism for explaining this fact. The scaling relation between galaxy mass and outer rotation velocity (also known as the baryonic Tully-Fisher relation, BTF) has recently been used as evidence against this viewpoint. We use a LCDM based semi-analytic disk galaxy formation model to investigate these claims. In our model, galaxies with less efficient star formation and higher gas fractions are more efficient at ejecting gas from galaxies. This is due to the fact that galaxies with less efficient star formation and higher gas fractions tend to live in dark matter haloes with lower circular velocities, from which less energy is required to escape the potential well. In our model the intrinsic scatter in the BTF is 0.15 dex, and mostly reflects scatter in dark halo concentration. The observed scatter, equal to 0.24 dex, is dominated by measurement errors. The best estimate for the intrinsic scatter is that it is less than 0.15 dex, and thus our LCDM based model (which does not include all possible sources of scatter) is only just consistent with this. In our model, gas rich galaxies, at fixed virial velocity (V_vir), with lower stellar masses have lower baryonic masses. This is consistent with the expectation that galaxies with lower stellar masses have had less energy available to drive an outflow. However, when the outer rotation velocity (V_flat) is used the correlation has the opposite sign, with a slope in agreement with observations. This is due to scatter in the relation between V_flat and V_vir. In summary, contrary to some previous claims, we show that basic features of the BTF are consistent with a LCDM based model in which the low efficiency of galaxy formation is determined by galactic outflows.

We construct a fully self-consistent mass model for the lens galaxy J2141 at z=0.14, and use it to improve on previous studies by modelling its gravitational lensing effect, gas rotation curve and stellar kinematics simultaneously. We adopt a very flexible axisymmetric mass model constituted by a generalized NFW dark matter halo and a stellar mass distribution obtained by deprojecting the MGE fit to the high-resolution K'-band LGSAO imaging data of the galaxy, with the (spatially constant) M/L ratio as a free parameter. We model the stellar kinematics by solving the anisotropic Jeans equations. We find that the inner logarithmic slope of the dark halo is weakly constrained (gamma = 0.82^{+0.65}_{-0.54}), and consistent with an unmodified NFW profile. We infer the galaxy to have (i) a dark matter fraction within 2.2 disk radii of 0.28^{+0.15}_{-0.10}, independent of the galaxy stellar population, implying a maximal disk for J2141; (ii) an apparently uncontracted dark matter halo, with concentration c_{-2} = 7.7_{-2.5}^{+4.2} and virial velocity v_{vir} = 242_{-39}^{+44} km/s, consistent with LCDM predictions; (iii) a slightly oblate halo (q_h = 0.75^{+0.27}_{-0.16}), consistent with predictions from baryon-affected models. Comparing the stellar mass inferred from the combined analysis (log_{10} Mstar/Msun = 11.12_{-0.09}^{+0.05}) with that inferred from SPS modelling of the galaxies colours, and accounting for a cold gas fraction of 20+/-10%, we determine a preference for a Chabrier IMF over Salpeter IMF by a Bayes factor of 5.7 (substantial evidence). We infer a value beta_{z} = 1 - sigma^2_{z}/sigma^2_{R} = 0.43_{-0.11}^{+0.08} for the orbital anisotropy parameter in the meridional plane, in agreement with most studies of local disk galaxies, and ruling out at 99% CL that the dynamics of this system can be described by a two-integral distribution function. [Abridged]

We describe the design and data sample from the DEEP2 Galaxy Redshift Survey, the densest and largest precision-redshift survey of galaxies at z ~ 1 completed to date. The survey has conducted a comprehensive census of massive galaxies, their properties, environments, and large-scale structure down to absolute magnitude M_B = -20 at z ~ 1 via ~90 nights of observation on the DEIMOS spectrograph at Keck Observatory. DEEP2 covers an area of 2.8 deg^2 divided into four separate fields, observed to a limiting apparent magnitude of R_AB=24.1. Objects with z < 0.7 are rejected based on BRI photometry in three of the four DEEP2 fields, allowing galaxies with z > 0.7 to be targeted ~2.5 times more efficiently than in a purely magnitude-limited sample. Approximately sixty percent of eligible targets are chosen for spectroscopy, yielding nearly 53,000 spectra and more than 38,000 reliable redshift measurements. Most of the targets which fail to yield secure redshifts are blue objects that lie beyond z ~ 1.45. The DEIMOS 1200-line/mm grating used for the survey delivers high spectral resolution (R~6000), accurate and secure redshifts, and unique internal kinematic information. Extensive ancillary data are available in the DEEP2 fields, particularly in the Extended Groth Strip, which has evolved into one of the richest multiwavelength regions on the sky. DEEP2 surpasses other deep precision-redshift surveys at z ~ 1 in terms of galaxy numbers, redshift accuracy, sample number density, and amount of spectral information. We also provide an overview of the scientific highlights of the DEEP2 survey thus far. This paper is intended as a handbook for users of the DEEP2 Data Release 4, which includes all DEEP2 spectra and redshifts, as well as for the publicly-available DEEP2 DEIMOS data reduction pipelines. [Abridged]

We present a public catalog of galaxy groups constructed from the spectroscopic sample of galaxies in the fourth data release from the DEEP2 Galaxy Redshift Survey, including the Extended Groth Strip (EGS). The catalog contains 1165 groups with two or more members in the EGS over the redshift range 0<z<1.5 and 1295 groups at z>0.6 in the rest of DEEP2. 25% of EGS galaxies and 14% of high-z DEEP2 galaxies are assigned to galaxy groups. The groups were detected using the Voronoi-Delaunay Method, after it has been optimized on mock DEEP2 catalogs following similar methods to those employed in Gerke et al. (2005). In the optimization effort, we have taken particular care to ensure that the mock catalogs resemble the data as closely as possible, and we have fine-tuned our methods separately on mocks constructed for the EGS and the rest of DEEP2. We have also probed the effect of the assumed cosmology on our inferred group-finding efficiency by performing our optimization on three different mock catalogs with different background cosmologies, finding large differences in the group-finding success we can achieve for these different mocks. Using the mock catalog whose background cosmology is most consistent with current data, we estimate that the DEEP2 group catalog is 72% complete and 61% pure (74% and 67% for the EGS) and that the group-finder correctly classifies 70% of galaxies that truly belong to groups, with an additional 46% of interloper galaxies contaminating the catalog (66% and 43% for the EGS). (Abridged)

We present gravitational lens models for 20 strong gravitational lens systems observed as part of the Sloan WFC Edge-on Late-type Lens Survey (SWELLS) project. Fifteen of the lenses are taken from paper I while five are newly discovered systems. The systems are galaxy-galaxy lenses where the foreground deflector has an inclined disc, with a wide range of morphological types, from late-type spiral to lenticular. For each system, we compare the total mass inside the critical curve inferred from gravitational lens modelling to the stellar mass inferred from stellar population synthesis (SPS) models, computing the stellar mass fraction f* = M(SPS)/M(lens). We find that, for the lower mass SWELLS systems, adoption of a Salpeter stellar initial mass function (IMF) leads to estimates of f* that exceed 1. This is unphysical, and provides strong evidence against the Salpeter IMF being valid for these systems. Taking the lower mass end of the SWELLS sample sigma(SIE) < 230 km/s, we find that the IMF is lighter (in terms of stellar mass-to-light ratio) than Salpeter with 98% probability, and consistent with the Chabrier IMF and IMFs between the two. This result is consistent with previous studies of spiral galaxies based on independent techniques. In combination with recent studies of massive early-type galaxies that have favoured a heavier Salpeter-like IMF, this result strengthens the evidence against a universal stellar IMF.

We determine an absolute calibration of stellar mass-to-light ratios for the densest \simeq 3% of early-type galaxies in the local universe (redshift z\simeq 0.08) from SDSS DR7. This sample of \sim 4000 galaxies has, assuming a Chabrier IMF, effective stellar surface densities, Sigma_e > 2500 M_sun/pc^2, stellar population synthesis (SPS) stellar masses log_10(M_sps/M_sun)<10.8, and aperture velocity dispersions of sigma_ap=168^{+37}_{-34} km/s (68% range). In contrast to typical early-type galaxies, we show that these dense early-type galaxies follow the virial fundamental plane, which suggests that mass-follows-light. With the additional assumption that any dark matter does not follow the light, the dynamical masses of dense galaxies provide a direct measurement of stellar masses. Our dynamical masses (M_dyn), obtained from the spherical Jeans equations, are only weakly sensitive to the choice of anisotropy (\beta) due to the relatively large aperture of the SDSS fiber for these galaxies: R_ap \simeq 1.5 R_e. Assuming isotropic orbits (\beta=0) we find a median log_{10} (M_dyn/M_sps) = 0.233 \pm 0.003, consistent with a Salpeter IMF, while more bottom heavy IMFs and standard Milky-Way IMFs are strongly disfavored. Our results are consistent with, but do not require, a dependence of the IMF on dynamical mass or velocity dispersion. We find evidence for a color dependence to the IMF such that redder galaxies have heavier IMFs with M_dyn/M_sps \propto (g-r)^{1.13\pm0.09}. This may reflect a more fundamental dependence of the IMF on the age or metallicity of a stellar population, or the density at which the stars formed.

We combine constraints on the galaxy-dark matter connection with structural and dynamical scaling relations to investigate the angular momentum content of disc galaxies. For haloes with masses in the interval 10^{11.3} < M_vir/M_sun < 10^{12.7} we find that the galaxy spin parameters are independent of halo mass with <\lambda'_gal> = (J_gal/M_gal) / (\sqrt{2} R_vir V_vir) = 0.019^{+0.004}_{-0.003} (1sigma). This is significantly lower than for relaxed LCDM haloes, which have an average spin parameter <\lambda'_halo> = 0.031. The average ratio between the specific angular momentum of disk galaxies and their host dark matter haloes is therefore R_j = \lambda'_gal/\lambda'_halo = 0.61^{+0.13}_{-0.10}. This calls into question a standard assumption made in the majority of all (semi-analytical) models for (disc) galaxy formation, namely that R_j=1. Using simple disc formation models we show that it is particularly challenging to understand why R_j is independent of halo mass, while the galaxy formation efficiency (\epsilon_GF, proportional to the ratio of galaxy mass to halo mass) reveals a strong halo mass dependence. We argue that the empirical scaling relations between \epsilon_GF, R_j and halo mass require both feedback (i.e., galactic outflows) and angular momentum transfer from the baryons to the dark matter (i.e., dynamical friction). The efficiency of angular momentum loss need to decrease with increasing halo mass. Such a mass dependence may reflect a bias against forming stable discs in high mass, low spin haloes or a transition from cold-mode accretion in low mass haloes to hot-mode accretion at the massive end. However, current hydrodynamical simulations of galaxy formation, which should include these processes, seem unable to reproduce the empirical relation between \epsilon_GF and R_j. We conclude that the angular momentum build-up of galactic discs remains poorly understood.

We have compiled a sample of 3041 spiral galaxies with multi-band gri imaging from the Sloan Digital Sky Survey (SDSS) Data Release 7 and available galaxy rotational velocities derived from HI line widths. We compare the data products provided through the SDSS imaging pipeline with our own photometry of the SDSS images, and use the velocities (V) as an independent metric to determine ideal galaxy sizes (R) and luminosities (L). Our radial and luminosity parameters improve upon the SDSS DR7 Petrosian radii and luminosities through the use of isophotal fits to the galaxy images. This improvement is gauged via VL and RV relations whose respective scatters are reduced by ~8% and ~30% compared to similar relations built with SDSS parameters. The tightest VRL relations are obtained with the i-band radius, R235i, measured at 23.5 mag/arcsec^-2, and the luminosity L235i, measured within R235i. Our VRL scaling relations compare well, both in scatter and slope, with similar studies (such comparisons however depend sensitively on the nature and size of the compared samples). The typical slopes, b, and observed scatters, sigma, of the i-band VL, RL and RV relations are bVL=0.27+/-0.01, bRL=0.41+/-0.01, bRV=1.52+/-0.07, and sigmaVL=0.074, sigmaRL=0.071, sigmaRV=0.154 dex. Similar results for the SDSS g and r bands are also provided. Smaller scatters may be achieved for more pruned samples. We also compute scaling relations in terms of the baryonic mass (stars + gas), Mbar, ranging from 10^8.7 Msol to 10^11.6 Msol. Our baryonic velocity-mass (VM) relation has slope 0.29+/-0.01 and a measured scatter sigma_meas = 0.076 dex. While the observed VL and VM relations have comparable scatter, the stellar and baryonic VM relations may be intrinsically tighter, and thus potentially more fundamental, than other VL relations of spiral galaxies.

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5\sigma\ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5\sigma\ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

Using data drawn from the DEEP2 and DEEP3 Galaxy Redshift Surveys, we investigate the relationship between the environment and the structure of galaxies residing on the red sequence at intermediate redshift. Within the massive (10 < log(M*/Msun) < 11) early-type population at 0.4 < z <1.2, we find a significant correlation between local galaxy overdensity (or environment) and galaxy size, such that early-type systems in higher-density regions tend to have larger effective radii (by ~0.5 kpc or 25% larger) than their counterparts of equal stellar mass and Sersic index in lower-density environments. This observed size-density relation is consistent with a model of galaxy formation in which the evolution of early-type systems at z < 2 is accelerated in high-density environments such as groups and clusters and in which dry, minor mergers (versus mechanisms such as quasar feedback) play a central role in the structural evolution of the massive, early-type galaxy population.

We present uncontaminated rest-frame u - R colors of 78 X-ray-selected AGN hosts at 0.5 < z < 1.5 in the Chandra Deep Fields measured with HST/ACS/NICMOS and VLT/ISAAC imaging. We also present spatially-resolved NUV - R color gradients for a subsample of AGN hosts imaged by HST/WFC3. Integrated, uncorrected photometry is not reliable for comparing the mean properties of soft and hard AGN host galaxies at z ~ 1 due to color contamination from point-source AGN emission. We use a cloning simulation to develop a calibration between concentration and this color contamination and use this to correct host galaxy colors. The mean u - R color of the unobscured/soft hosts beyond ~6 kpc is statistically equivalent to that of the obscured/hard hosts (the soft sources are 0.09 +/- 0.16 magnitudes bluer). Furthermore, the rest-frame V - J colors of the obscured and unobscured hosts beyond ~6 kpc are statistically equivalent, suggesting that the two populations have similar distributions of dust extinction. For the WFC3/IR sample, the mean NUV - R color gradients of unobscured and obscured sources differ by less than ~0.5 magnitudes for r > 1.1 kpc. These three observations imply that AGN obscuration is uncorrelated with the star formation rate beyond ~1 kpc. These observations favor a unification scenario for intermediate-luminosity AGNs in which obscuration is determined geometrically. Scenarios in which the majority of intermediate-luminosity AGN at z ~ 1 are undergoing rapid, galaxy-wide quenching due to AGN-driven feedback processes are disfavored.

The relative contribution of baryons and dark matter to the inner regions of spiral galaxies provides critical clues to their formation and evolution, but it is generally difficult to determine. For spiral galaxies that are strong gravitational lenses, however, the combination of lensing and kinematic observations can be used to break the disk-halo degeneracy. In turn, such data constrain fundamental parameters such as i) the mass density profile slope and axis ratio of the dark matter halo, and by comparison with dark matter-only numerical simulations the modifications imposed by baryons; ii) the mass in stars and therefore the overall star formation efficiency, and the amount of feedback; iii) by comparison with stellar population synthesis models, the normalization of the stellar initial mass function. In this first paper of a series, we present a sample of 16 secure, 1 probable, and 6 possible strong lensing spiral galaxies, for which multi-band high-resolution images and rotation curves were obtained using the Hubble Space Telescope and Keck-II Telescope as part of the Sloan WFC Edge-on Late-type Lens Survey (SWELLS). The sample includes 8 newly discovered secure systems. [abridged] We find that the SWELLS sample of secure lenses spans a broad range of morphologies (from lenticular to late-type spiral), spectral types (quantified by Halpha emission), and bulge to total stellar mass ratio (0.22-0.85), while being limited to M_*>10^{10.5} M_sun. The SWELLS sample is thus well-suited for exploring the relationship between dark and luminous matter in a broad range of galaxies. We find that the deflector galaxies obey the same size-mass relation as that of a comparison sample of elongated non-lens galaxies selected from the SDSS survey. We conclude that the SWELLS sample is consistent with being representative of the overall population of high-mass high-inclination disky galaxies.

We investigate the origin of the relations between stellar mass and optical circular velocity for early-type (ETG) and late-type (LTG) galaxies --- the Faber-Jackson (FJ) and Tully-Fisher (TF) relations. We combine measurements of dark halo masses (from satellite kinematics and weak lensing), and the distribution of baryons in galaxies (from a new compilation of galaxy scaling relations), with constraints on dark halo structure from cosmological simulations. The principle unknowns are the halo response to galaxy formation and the stellar initial mass function (IMF). The slopes of the TF and FJ relations are naturally reproduced for a wide range of halo response and IMFs. However, models with a universal IMF and universal halo response cannot simultaneously reproduce the zero points of both the TF and FJ relations. For a model with a universal Chabrier IMF, LTGs require halo expansion, while ETGs require halo contraction. A Salpeter IMF is permitted for high mass (sigma > 180 km/s) ETGs, but is inconsistent for intermediate masses, unless V_circ(R_e)/sigma_e > 1.6. If the IMF is universal and close to Chabrier, we speculate that the presence of a major merger may be responsible for the contraction in ETGs while clumpy accreting streams and/or feedback leads to expansion in LTGs. Alternatively, a recently proposed variation in the IMF disfavors halo contraction in both types of galaxies. Finally we show that our models naturally reproduce flat and featureless circular velocity profiles within the optical regions of galaxies without fine-tuning.

The degeneracy among the disk, bulge and halo contributions to galaxy rotation curves prevents an understanding of the distribution of baryons and dark matter in disk galaxies. In an attempt to break this degeneracy, we present an analysis of the spiral galaxy strong gravitational lens SDSS J2141-0001, discovered as part of the SLACS survey. We present new Hubble Space Telescope multicolor imaging, gas and stellar kinematics data derived from long-slit spectroscopy, and K-band LGS adaptive optics imaging, both from the Keck telescopes. We model the galaxy as a sum of concentric axisymmetric bulge, disk and halo components and infer the contribution of each component, using information from gravitational lensing and gas kinematics. This analysis yields a best-fitting total (disk plus bulge) stellar mass of log_{10}(Mstar/Msun) = 10.99(+0.11,-0.25). The photometric data combined with stellar population synthesis models yield log_{10}(Mstar/Msun) = 10.97\pm0.07, and 11.21\pm0.07 for the Chabrier and Salpeter IMFs, respectively. Accounting for the expected gas fraction of \simeq 20% reduces the lensing plus kinematics stellar mass by 0.10\pm0.05 dex, resulting in a Bayes factor of 11.9 in favor of a Chabrier IMF. The dark matter halo is roughly spherical, with minor to major axis ratio q_{halo}=0.91(+0.15,-0.13). The dark matter halo has a maximum circular velocity of V_{max}=276(+17,-18) km/s, and a central density parameter of log_{10}\Delta_{V/2}=5.9(+0.9,-0.5). This is higher than predicted for uncontracted dark matter haloes in LCDM cosmologies, log_{10}\Delta_{V/2}=5.2, suggesting that either the halo has contracted in response to galaxy formation, or that the halo has a higher than average concentration. At 2.2 disk scale lengths the dark matter fraction is f_{DM}=0.55(+0.20,-0.15), suggesting that SDSS J2141-0001 is sub-maximal.

We study the mass and redshift dependence of the concentration parameter in Nbody simulations spanning masses from $10^{10} \hMsun$ to $10^{15} \hMsun$ and redshifts from 0 to 2. We present a series of fitting formulas that accurately describe the time evolution of the concentration-mass relation since z=2. Using arguments based on the spherical collapse model we study the behaviour of the scale length of the density profile during the assembly history of haloes, obtaining physical insights on the origin of the observed time evolution of the concentration mass relation. We present preliminary results of the implementation of this model in the prediction of the values of the concentration parameter for different masses and redshifts.

We study the evolution of the scaling relations between maximum circular velocity, stellar mass and optical half-light radius of star-forming disk-dominated galaxies in the context of LCDM-based galaxy formation models. Using data from the literature combined with new data from the DEEP2 and AEGIS surveys we show that there is a consistent observational and theoretical picture for the evolution of these scaling relations from z\sim 2 to z=0. The evolution of the observed stellar scaling relations is weaker than that of the virial scaling relations of dark matter haloes, which can be reproduced, both qualitatively and quantitatively, with a simple, cosmologically-motivated model for disk evolution inside growing NFW dark matter haloes. In this model optical half-light radii are smaller, both at fixed stellar mass and maximum circular velocity, at higher redshifts. This model also predicts that the scaling relations between baryonic quantities evolve even more weakly than the corresponding stellar relations. We emphasize, though, that this weak evolution does not imply that individual galaxies evolve weakly. On the contrary, individual galaxies grow strongly in mass, size and velocity, but in such a way that they move largely along the scaling relations. Finally, recent observations have claimed surprisingly large sizes for a number of star-forming disk galaxies at z \sim 2, which has caused some authors to suggest that high redshift disk galaxies have abnormally high spin parameters. However, we argue that the disk scale lengths in question have been systematically overestimated by a factor \sim 2, and that there is an offset of a factor \sim 1.4 between H\alpha sizes and optical sizes. Taking these effects into account, there is no indication that star forming galaxies at high redshifts (z\sim 2) have abnormally high spin parameters.

We present a detailed study of the redshift evolution of dark matter halo structural parameters in a LambdaCDM cosmology. We study the mass and redshift dependence of the concentration, shape and spin parameter in Nbody simulations spanning masses from 10^{10} Msun/h to 10^{15} Msun/h and redshifts from 0 to 2. We present a series of fitting formulas that accurately describe the time evolution of the concentration-mass relation since z=2. Using arguments based on the spherical collapse model we study the behaviour of the scale length of the density profile during the assembly history of haloes, obtaining physical insights on the origin of the observed time evolution of the concentration mass relation. We also investigate the evolution with redshift of dark matter halo shape and its dependence on mass. Within the studied redshift range the relation between halo shape and mass can be well fitted by a power law. Finally we show that although for z=0 the spin parameter is practically mass independent, at increasing redshift it shows a increasing correlation with mass.

Using estimates of dark halo masses from satellite kinematics, weak gravitational lensing, and halo abundance matching, combined with the Tully-Fisher and Faber-Jackson relations, we derive the mean relation between the optical, V_opt, and virial, V_200, circular velocities of early- and late-type galaxies at redshift z~0. For late-type galaxies V_opt ~ V_200 over the velocity range V_opt=90-260 km/s, and is consistent with V_opt = V_maxh (the maximum circular velocity of NFW dark matter haloes in the concordance LCDM cosmology). However, for early-type galaxies V_opt \ne V_200, with the exception of early-type galaxies with V_opt simeq 350 km/s. This is inconsistent with early-type galaxies being, in general, globally isothermal. For low mass (V_opt < 250 km/s) early-types V_opt > V_maxh, indicating that baryons have modified the potential well, while high mass (V_opt > 400 km/s) early-types have V_opt < V_maxh. Folding in measurements of the black hole mass - velocity dispersion relation, our results imply that the supermassive black hole - halo mass relation has a logarithmic slope which varies from ~1.4 at halo masses of ~10^{12} Msun/h to ~0.65 at halo masses of 10^{13.5} Msun/h. The values of V_opt/V_200 we infer for the Milky Way and M31 are lower than the values currently favored by direct observations and dynamical models. This offset is due to the fact that the Milky Way and M31 have higher V_opt and lower V_200 compared to typical late-type galaxies of the same stellar masses. We show that current high resolution cosmological hydrodynamical simulations are unable to form galaxies which simultaneously reproduce both the V_opt/V_200 ratio and the V_opt-M_star (Tully-Fisher/Faber-Jackson) relation.

We use a semi-analytic model for disk galaxies to explore the origin of the time evolution and small scatter of the galaxy SFR sequence -- the tight correlation between star-formation rate (SFR) and stellar mass (M_star). The steep decline of SFR from z~2 to the present, at fixed M_star, is a consequence of the following: First, disk galaxies are in a steady state with the SFR following the net (i.e., inflow minus outflow) gas accretion rate. The evolution of the SFR sequence is determined by evolution in the cosmological specific accretion rates, \propto (1+z)^{2.25}, but is found to be independent of feedback. Although feedback determines the outflow rates, it shifts galaxies along the SFR sequence, leaving its zero point invariant. Second, the conversion of accretion rate to SFR is materialized through gas density, not gas mass. Although the model SFR is an increasing function of both gas mass fraction and gas density, only the gas densities are predicted to evolve significantly with redshift. Third, star formation is fueled by molecular gas. Since the molecular gas fraction increases monotonically with increasing gas density, the model predicts strong evolution in the molecular gas fractions, increasing by an order of magnitude from z=0 to z~2. On the other hand, the model predicts that the effective surface density of atomic gas is ~10 M_sun pc^{-2}, independent of redshift, stellar mass or feedback. Our model suggests that the scatter in the SFR sequence reflects variations in the gas accretion history, and thus is insensitive to stellar mass, redshift or feedback. The large scatter in halo spin contributes negligibly, because it scatters galaxies along the SFR sequence. An observational consequence of this is that the scatter in the SFR sequence is independent of the size (both stellar and gaseous) of galaxy disks.

We use a disk formation model to study the effects of galactic outflows (a.k.a. feedback) on the rotation velocity - stellar mass - disk size, gas fraction - stellar mass, and gas phase metalicity - stellar mass scaling relations of disk galaxies. We show that models without outflows are unable to explain these scaling relations, having both the wrong slopes and normalization. The problem can be traced to the model galaxies having too many baryons. Models with outflows can solve this "over-cooling" problem by removing gas before it has time to turn into stars. Models with both momentum and energy driven winds can reproduce the observed scaling relations. However, these models predict different slopes which, with better observations, may be used to discriminate between these models.

There is a large observational scatter toward low velocities in the stellar mass Tully-Fisher relation if disturbed and compact objects are included. However, this scatter can be eliminated if one replaces rotation velocity with $\rm S_{\rm 0.5}$, a quantity that includes a velocity dispersion term added in quadrature with the rotation velocity. In this work we use a large suite of hydrodynamic N-body galaxy merger simulations to explore a possible mechanism for creating the observed relations. Using mock observations of the simulations, we test for the presence of observational effects and explore the relationship between $\rm S_{\rm 0.5}$ and intrinsic properties of the galaxies. We find that galaxy mergers can explain the scatter in the TF as well as the tight $\rm S_{\rm 0.5}$-stellar mass relation. Furthermore, $\rm S_{\rm 0.5}$ is correlated with the total central mass of a galaxy, including contributions due to dark matter.

The quasar SDSS J105041.35+345631.3 (z = 0.272) has broad emission lines blueshifted by 3500 km/s relative to the narrow lines and the host galaxy. Such an object may be a candidate for a recoiling supermassive black hole, binary black hole, a superposition of two objects, or an unusual geometry for the broad emission-line region. The absence of narrow lines at the broad line redshift argues against superposition. New Keck spectra of J1050+3546 place tight constraints on the binary model. The combination of large velocity shift and symmetrical H-beta profile, as well as aspects of the narrow line spectrum, make J1050+3546 an interesting candidate for black hole recoil. Other aspects of the spectrum, however, suggest that the object is most likely an extreme case of a ``double-peaked emitter.'' We discuss possible observational tests to determine the true nature of this exceptional object.

We present a new technique for observationally identifying galaxy mergers spectroscopically rather than through host galaxy imaging. Our technique exploits the dynamics of supermassive black holes (SMBHs) powering active galactic nuclei (AGNs) in merger-remnant galaxies. Because structure in the universe is built up through galaxy mergers and nearly all galaxies host a central SMBH, some galaxies should possess two SMBHs near their centers as the result of a recent merger. These SMBHs spiral to the center of the resultant merger-remnant galaxy, and one or both of the SMBHs may power AGNs. Using the DEEP2 Galaxy Redshift Survey, we have examined 1881 red galaxies, of which 91 exhibit [O III] and Hbeta emission lines indicative of Seyfert 2 activity. Of these, 32 AGNs have [O III] emission-line redshifts significantly different from the redshifts of the host galaxies' stars, corresponding to velocity offsets of ~50 km/s to ~300 km/s. Two of these AGNs exhibit double-peaked [O III] emission lines, while the remaining 30 AGNs each exhibit a single set of velocity-offset [O III] emission lines. After exploring a variety of physical models for these velocity offsets, we argue that the most likely explanation is inspiralling SMBHs in merger-remnant galaxies. Based on this interpretation, we find that roughly half of the red galaxies hosting AGNs are also merger remnants, which implies that mergers may trigger AGN activity in red galaxies. The AGN velocity offsets we find imply a merger fraction of ~30% and a merger rate of ~3 mergers/Gyr for red galaxies at redshifts 0.34 < z < 0.82.

We use a disk galaxy evolution model to investigate the impact of mass outflows (a.k.a. feedback) on disk galaxy scaling relations. Our model follows the accretion, cooling, star formation and ejection of baryonic mass inside growing dark matter haloes, with cosmologically motivated specific angular momentum distributions. Models without feedback produce disks that are too small and rotate too fast. Feedback reduces the baryonic masses of galaxies, resulting in larger disks with lower rotation velocities. Models with feedback can reproduce the zero points of the scaling relations between rotation velocity, stellar mass and disk size, but only in the absence of adiabatic contraction. Our feedback mechanism is maximally efficient in expelling mass, but our successful models require 25% of the SN energy, or 100% of the SN momentum, to drive the outflows. It remains to be seen whether such high efficiencies are realistic or not. Our energy and momentum driven wind models result in different slopes of various scaling relations, such as size - stellar mass, stellar mass - halo mass, and metallicity - stellar mass. Observations favor the energy driven wind at stellar masses below Mstar = 10^{10.5} Msun, but the momentum driven wind model at high masses. The ratio between the specific angular momentum of the baryons to that of the halo, (j_gal/m_gal), is not unity in our models. Yet this is the standard assumption in models of disk galaxy formation. Feedback preferentially ejects low angular momentum material because star formation is more efficient at smaller galactic radii. This results in (j_gal/m_gal) increasing with decreasing halo mass. This effect helps to resolve the discrepancy between the high spin parameters observed for dwarf galaxies with the low spin parameters predicted from LCDM. [Abridged]

We use a disk galaxy evolution model to investigate whether galaxies with exponential surface brightness profiles can be produced in a cosmologically motivated framework for disk galaxy formation. Our model follows the accretion, cooling, and ejection of baryonic mass, as a function of radius, inside growing dark matter haloes. The surface density profile of the disk is determined by detailed angular momentum conservation, starting from the distribution of specific angular momentum as found in cosmological simulations. Exponential and quasi-exponential disks can be produced by our model through a combination of supernova driven galactic outflows (which preferentially remove low angular momentum material), intrinsic variation in the angular momentum distribution of the halo gas, and the inefficiency of star formation at large radii. We use observations from the SDSS NYU-VAGC to show that the median Sersic index of blue galaxies is a strong function of stellar mass. For blue galaxies, low mass galaxies have a median n=1.3, while high mass galaxies have a median n=4. Our model with energy driven outflows correctly reproduces this trend, whereas our models with momentum driven outflows and no outflows over predict the Sersic indices in low mass galaxies. We show that the observed fraction of "bulge-less" exponential galaxies is a strong function of stellar mass. For Milky-Way mass galaxies less than 0.1% of blue galaxies are bulge-less, whereas for M33 mass galaxies bulge-less and quasi-bulgeless galaxies are typical. These results suggest that the difficulty of hierarchical formation models to produce bulge-less Milky-Way mass galaxies is in fact not a problem. However, the problem of producing M33 like galaxies remains, and will provide a key test for hierarchical galaxy formation models. [Abridged]

Whereas considerable effort has been afforded in understanding the properties of galaxies, a full physical picture, connecting their baryonic and dark-matter content, super-massive black holes, and (metric) theories of gravity, is still ill-defined. Strong gravitational lensing furnishes a powerful method to probe gravity in the central regions of galaxies. It can (1) provide a unique detection-channel of dark-matter substructure beyond the local galaxy group, (2) constrain dark-matter physics, complementary to direct-detection experiments, as well as metric theories of gravity, (3) probe central super-massive black holes, and (4) provide crucial insight into galaxy formation processes from the dark matter point of view, independently of the nature and state of dark matter. To seriously address the above questions, a considerable increase in the number of strong gravitational-lens systems is required. In the timeframe 2010-2020, a staged approach with radio (e.g. EVLA, e-MERLIN, LOFAR, SKA phase-I) and optical (e.g. LSST and JDEM) instruments can provide 10^(2-4) new lenses, and up to 10^(4-6) new lens systems from SKA/LSST/JDEM all-sky surveys around ~2020. Follow-up imaging of (radio) lenses is necessary with moderate ground/space-based optical-IR telescopes and with 30-50m telescopes for spectroscopy (e.g. TMT, GMT, ELT). To answer these fundamental questions through strong gravitational lensing, a strong investment in large radio and optical-IR facilities is therefore critical in the coming decade. In particular, only large-scale radio lens surveys (e.g. with SKA) provide the large numbers of high-resolution and high-fidelity images of lenses needed for SMBH and flux-ratio anomaly studies.

There is a vast menagerie of plausible candidates for the constituents of dark matter, both within and beyond extensions of the Standard Model of particle physics. Each of these candidates may have scattering (and other) cross section properties that are consistent with the dark matter abundance, BBN, and the most scales in the matter power spectrum; but which may have vastly different behavior at sub-galactic "cutoff" scales, below which dark matter density fluctuations are smoothed out. The only way to quantitatively measure the power spectrum behavior at sub-galactic scales at distances beyond the local universe, and indeed over cosmic time, is through probes available in multiply imaged strong gravitational lenses. Gravitational potential perturbations by dark matter substructure encode information in the observed relative magnifications, positions, and time delays in a strong lens. Each of these is sensitive to a different moment of the substructure mass function and to different effective mass ranges of the substructure. The time delay perturbations, in particular, are proving to be largely immune to the degeneracies and systematic uncertainties that have impacted exploitation of strong lenses for such studies. There is great potential for a coordinated theoretical and observational effort to enable a sophisticated exploitation of strong gravitational lenses as direct probes of dark matter properties. This opportunity motivates this white paper, and drives the need for: a) strong support of the theoretical work necessary to understand all astrophysical consequences for different dark matter candidates; and b) tailored observational campaigns, and even a fully dedicated mission, to obtain the requisite data.

Galactic winds are a prime suspect for the metal enrichment of the intergalactic medium and may have a strong influence on the chemical evolution of galaxies and the nature of QSO absorption line systems. We use a sample of 1406 galaxy spectra at z~1.4 from the DEEP2 redshift survey to show that blueshifted Mg II 2796, 2803 A absorption is ubiquitous in starforming galaxies at this epoch. This is the first detection of frequent outflowing galactic winds at z~1. The presence and depth of absorption are independent of AGN spectral signatures or galaxy morphology; major mergers are not a prerequisite for driving a galactic wind from massive galaxies. Outflows are found in coadded spectra of galaxies spanning a range of 30x in stellar mass and 10x in star formation rate (SFR), calibrated from K-band and from MIPS IR fluxes. The outflows have column densities of order N_H ~ 10^20 cm^-2 and characteristic velocities of ~ 300-500 km/sec, with absorption seen out to 1000 km/sec in the most massive, highest SFR galaxies. The velocities suggest that the outflowing gas can escape into the IGM and that massive galaxies can produce cosmologically and chemically significant outflows. Both the Mg II equivalent width and the outflow velocity are larger for galaxies of higher stellar mass and SFR, with V_wind ~ SFR^0.3, similar to the scaling in low redshift IR-luminous galaxies. The high frequency of outflows in the star-forming galaxy population at z~1 indicates that galactic winds occur in the progenitors of massive spirals as well as those of ellipticals. The increase of outflow velocity with mass and SFR constrains theoretical models of galaxy evolution that include feedback from galactic winds, and may favor momentum-driven models for the wind physics.

We investigate the effects of changes in the cosmological parameters between the WMAP 1st, 3rd, and 5th year results on the structure of dark matter haloes. We use a set of simulations that cover 5 decades in halo mass ranging from the scales of dwarf galaxies (V_c ~30 km/s) to clusters of galaxies (V_c ~ 1000 km/s). We find that the concentration mass relation is a power law in all three cosmologies. However the slope is shallower and the zero point is lower moving from WMAP1 to WMAP5 to WMAP3. For haloes of mass log(M_200/Msun) = 10, 12, and 14 the differences in the concentration parameter between WMAP1 and WMAP3 are a factor of 1.55, 1.41, and 1.29, respectively. As we show, this brings the central densities of dark matter haloes in good agreement with the central densities of dwarf and low surface brightness galaxies inferred from their rotation curves, for both the WMAP3 and WMAP5 cosmologies. We also show that none of the existing toy models for the concentration-mass relation can reproduce our simulation results over the entire range of masses probed. In particular, the model of Bullock et al (B01) fails at the higher mass end (M > 1e13 Msun), while the NFW model of Navarro, Frenk & White (1997) fails dramatically at the low mass end (M < 1e12 Msun). We present a new model, based on a simple modification of that of B01, which reproduces the concentration-mass relations in our simulations over the entire range of masses probed (1e10 Msun < M < 1e15 Msun). Haloes in the WMAP3 cosmology (at a fixed mass) are more flatted compared to the WMAP1 cosmology, with a medium to long axis ration reduced by ~10 %. Finally, we show that the distribution of halo spin parameters is the same for all three cosmologies.

The scaling relations between rotation velocity, size and luminosity form a benchmark test for any theory of disk galaxy formation. We confront recent theoretical models of disk formation to a recent large compilation of such scaling relations. We stress the importance of achieving a fair comparison between models and observations.

A long standing problem for hierarchical disk galaxy formation models has been the simultaneous matching of the zero point of the Tully-Fisher relation and the galaxy luminosity function (LF). We illustrate this problem for a typical disk galaxy and discuss three solutions: low stellar mass-to-light ratios, low initial dark halo concentrations, and no halo contraction. We speculate that halo contraction may be reversed through a combination of mass ejection through feedback and angular momentum exchange brought about by dynamical friction between baryons and dark matter during the disk formation process.

We construct a large data set of global structural parameters for 1300 field and cluster spiral galaxies and explore the joint distribution of luminosity L, optical rotation velocity V, and disk size R at I- and 2MASS K-bands. The I- and K-band velocity-luminosity (VL) relations have log-slopes of 0.29 and 0.27, respectively with sigma_ln(VL)~0.13, and show a small dependence on color and morphological type in the sense that redder, early-type disk galaxies rotate faster than bluer, later-type disk galaxies for most luminosities. The VL relation at I- and K-bands is independent of surface brightness, size and light concentration. The log-slope of the I- and K-band RL relations is a strong function of morphology and varies from 0.25 to 0.5. The average dispersion sigma_ln(RL) decreases from 0.33 at I-band to 0.29 at K, likely due to the 2MASS selection bias against lower surface brightness galaxies. Measurement uncertainties are sigma_ln(V)~0.09, sigma_ln(L)~0.14 and somewhat larger and harder to estimate for ln(R). The color dependence of the VL relation is consistent with expectations from stellar population synthesis models. The VL and RL residuals are largely uncorrelated with each other; the RV-RL residuals show only a weak positive correlation. These correlations suggest that scatter in luminosity is not a significant source of the scatter in the VL and RL relations. The observed scaling relations can be understood in the context of a model of disk galaxies embedded in dark matter halos that invokes low mean spin parameters and dark halo expansion, as we describe in our companion paper (Dutton et al. 2007). We discuss in two appendices various pitfalls of standard analytical derivations of galaxy scaling relations, including the Tully-Fisher relation with different slopes. (Abridged).

We use a series of cosmological N-body simulations for a flat LCDM cosmology to investigate the properties of dark matter haloes in the mass range 3.0e9-3.0e13 Msun. These properties include the concentration parameter (c), the spin parameter (lambda) and the mean axis ratio (q). For the concentration-mass relation we find c~M^(-0.11) in agreement with the model proposed by Bullock et al. even if we find a lower normalization (15%). The results for lambda and q are in good agreement with previous studies, while c and lambda are anti-correlated. In an attempt to remove unrelaxed haloes, we use the offset parameter (xoff), defined as the distance between the most bound particle and the center of mass. Removing haloes with large xoff increases the c by ~10%, lowers the lambda by ~15%, and removes the most prolate haloes. In addition, it largely removes the anti-correlation between c and lambda though not entirely. We also investigate the effects of the large-scale environment. We find that more concentrated haloes live in denser environments. Note, however, that the trend is weak compared to the scatter. For the spin parameters we find no environment dependence, while there is a weak indication that the most spherical haloes reside in denser region. Finally, using a simple model for disk galaxy formation we show that haloes that host low surface brightness galaxies are expected to be hosted by a biased sub-set of haloes. Not only do these haloes have spin parameters that are larger than average, they also have c that are 15% percent lower than the average at a given halo mass. We discuss the implications of all these findings for the claimed disagreement between halo concentrations inferred from LSB rotation curves, and those expected for a LCDM cosmology. (abridged)

We use observed rotation velocity-luminosity (VL) and size-luminosity (RL) relations to single out a specific scenario for disk galaxy formation in the LCDM cosmology. Our model involves four independent log-normal random variables: dark-halo concentration c, disk spin lam_gal, disk mass fraction m_gal, and stellar mass-to-light ratio M/L_I. A simultaneous match of the VL and RL zero points with adiabatic contraction requires low-c halos, but this model has V_2.2~1.8 V_vir (where V_2.2 and V_vir are the circular velocity at 2.2 disk scale lengths and the virial radius, respectively) which will be unable to match the luminosity function (LF). Similarly models without adiabatic contraction but standard c also predict high values of V_2.2/V_vir. Models in which disk formation induces an expansion rather than the commonly assumed contraction of the dark-matter halos have V_2.2~1.2 V_vir which allows a simultaneous fit of the LF. This may result from non-spherical, clumpy gas accretion, where dynamical friction transfers energy from the gas to the dark matter. This model requires low lam_gal and m_gal values, contrary to naive expectations. However, the low lam_gal is consistent with the notion that disk galaxies predominantly survive in halos with a quiet merger history, while a low m_gal is also indicated by galaxy-galaxy lensing. The smaller than expected scatter in the RL relation, and the lack of correlation between the residuals of the VL and RL relations, respectively, imply that the scatter in lam_gal and in c need to be smaller than predicted for LCDM halos, again consistent with the idea that disk galaxies preferentially reside in halos with a quiet merger history.

This paper addresses available constraints on mass models fitted to rotation curves. Mass models of disk galaxies have well-known degeneracies, that prevent a unique mass decomposition. The most notable is due to the unknown value of the stellar mass-to-light ratio (the disk-halo degeneracy); even with this known, degeneracies between the halo parameters themselves may prevent an unambiguous determination of the shape of the dark halo profile, which includes the inner density slope of the dark matter halo. The latter is often referred to as the ``cusp-core degeneracy''. We explore constraints on the disk and halo parameters and apply these to four mock and six observed disk galaxies with high resolution and extended rotation curves. Our full set of constraints consists of mass-to-light (M/L) ratios from stellar population synthesis models based on B-R colors, constraints on halo parameters from N-body simulations, and constraining the halo virial velocity to be less than the maximum observed velocity. These constraints are only partially successful in lifting the cusp-core degeneracy. The effect of adiabatic contraction of the halo by the disk is to steepen cores into cusps and reduce the best-fit halo concentration and M/L values (often significantly). We also discuss the effect of disk thickness, halo flattening, distance errors, and rotation curve error values on mass modeling. Increasing the imposed minimum rotation curve error from typically low, underestimated values to more realistic estimates decreases the chi-square substantially and makes distinguishing between a cuspy or cored halo profile even more difficult. In spite of the degeneracies and uncertainties present, our constrained mass modeling favors sub-maximal disks (i.e., a dominant halo) at 2.2 disk scale lengths, with Vdisk/Vtot < 0.6.

We use a simple model of disc-galaxy formation in the LCDM cosmology to simultaneously reproduce the slopes, zero-points, scatter and uncorrelated residuals in the observed velocity-luminosity (VL) and radius-luminosity (RL) relations. Observed I-band luminosities are converted to stellar masses using the IMF-dependent relation between stellar mass-to-light ratio, M/L, and color. The model treats halo concentration, spin parameter, and disc mass fraction as independent log-normal random variables. Our main conclusion is that the VL and RL zero-points and the uncorrelated residuals can only be reproduced simultaneously if adiabatic contraction is avoided. One or the other could be fixed by appealing to unrealistic values for M/L, halo concentration c, and spin parameter lambda, but not all together. The small VL scatter is naturally determined by the predicted scatter in c and in M/L, quite independent of the large scatter in lambda. However, the RL scatter, driven by the scatter in lambda, can be as low as observed only if sigma_lambda simeq 0.25, about half the value predicted for CDM haloes. This may indicate that discs form in a special subset of haloes. The VL slope is reproduced once star formation occurs only above a threshold surface density, and the RL slope implies that the disc mass fraction is increasing with halo mass, consistent with feedback effects. A model that incorporates the above ingredients provides a simultaneous fit to all the observed features. In particular, the elimination of adiabatic contraction allows 80% disc contribution to the observed rotation velocity of bright discs at 2.2 disc scale lengths. The lack of halo contraction may indicate that disc formation is not as smooth as typically envisioned, but instead may involve clumpy, cold streams.