Kinematic studies of disk galaxies, using individual stars in the Milky Way or statistical studies of global disk kinematics over time, provide insight into how disks form and evolve. We use a high-resolution, cosmological zoom-simulation of a Milky Way-mass disk galaxy h277 to tie together local disk kinematics and the evolution of the disk over time. The present-day stellar age-velocity relationship (AVR) of h277 is nearly identical to that of the analogous solar-neighborhood measurement in the Milky Way. A crucial element of this success is the simulation's dynamically cold multi-phase ISM, which allows young stars to form with a low velocity dispersion ($\sigma_{\mathrm{birth}}$$\sim 6 - 8\ \mathrm{km}\, \mathrm{s}^{-1}$) at late times. Older stars are born kinematically hotter (i.e., the disk settles over time in an "upside-down" formation scenario), and are subsequently heated after birth. The disk also grows "inside-out", and many of the older stars in the solar neighborhood at $z=0$ are present because of radial mixing. We demonstrate that the evolution of $\sigma_{\mathrm{birth}}$ in h277 can be explained by the same model used to describe the general decrease in velocity dispersion observed in disk galaxies from $z\sim 2-3$ to the present-day, in which the disk evolves in quasi-stable equilibrium and the ISM velocity dispersion decreases over time due to a decreasing gas fraction. Thus, our results tie together local observations of the Milky Way's AVR with observed kinematics of high $z$ disk galaxies.

Cosmological simulations are reaching the resolution necessary to study ultra-faint dwarf galaxies. Observations indicate that in small populations, the stellar initial mass function (IMF) is not fully populated; rather, stars are sampled in a way that can be approximated as coming from an underlying probability density function. To ensure the accuracy of cosmological simulations in the ultra-faint regime, we present an improved treatment of the IMF. We implement a self-consistent, stochastically populated IMF in cosmological hydrodynamic simulations. We test our method using high-resolution simulations of a Milky Way halo, run to $z=6$, yielding a sample of nearly 100 galaxies. We also use an isolated dwarf galaxy to investigate the resulting systematic differences in galaxy properties. We find that a stochastic IMF in simulations makes feedback burstier, strengthening feedback, and quenching star formation earlier in small dwarf galaxies. For galaxies in halos with mass $\lesssim10^{8.5}$ M$_\odot$, a stochastic IMF typically leads to lower stellar mass compared to a continuous IMF, sometimes by more than an order of magnitude. We show that existing methods of ensuring discrete supernovae incorrectly determine the mass of the star particle and its associated feedback. This leads to overcooling of surrounding gas, with at least ${\sim}$10 per cent higher star formation and ${\sim}$30 per cent higher cold gas content. Going forward, to accurately model dwarf galaxies and compare to observations, it will be necessary to incorporate a stochastically populated IMF that samples the full spectrum of stellar masses.

Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution.

The existence of ultra-faint dwarf (UFD) galaxies highlights the need to push our theoretical understanding of galaxies to extremely low mass. We examine the formation of UFDs by twice running a fully cosmological simulations of dwarf galaxies, but varying star formation. One run uses a temperature-density threshold for star formation, while the other uses an H$_{2}$-based sub-grid star formation model. The total number of dwarf galaxies that forms is different by a factor of 2 between the two runs, but most of these are satellites, leading to a factor of 5 difference in the number of luminous UFD companions around more massive, isolated dwarfs. The first run yields a 47\% chance of finding a satellite around a M$_{halo}$ $\sim 10^{10}$ M$_{\odot}$ host, while the H$_2$ run predicts only a 16\% chance. Metallicity is the primary physical parameter that creates this difference. As metallicity decreases, the formation of H$_2$ is slowed and relegated to higher-density material. Thus, our H$_2$ run is unable to form many (and often, any) stars before reionization removes gas. These results emphasize that predictions for UFD properties made using hydrodynamic simulations, in particular regarding the frequency of satellites around dwarf galaxies, the slope of the stellar mass function at low masses, as well as the properties of ultra-faint galaxies occupying the smallest halos, are extremely sensitive to the subgrid physics of star formation contained within the simulation. However, upcoming discoveries of ultra-faint dwarfs will provide invaluable constraining power on the physics of the first star formation.

Star formation is observed to be strongly correlated to dense regions of molecular gas. Although the exact nature of the link between star formation and molecular hydrogen is still unclear, some have suggested that shielding of dense gas by dust grains is the key factor enabling the presence of both. We present a sub-grid model for use in galaxy formation simulations in which star formation is linked explicitly to local dust shielding. We developed and tested our shielding and star formation models using smoothed particle hydrodynamic simulations of solar and sub-solar metallicity isolated Milky Way-mass disk galaxies. We compared our dust shielding-based star formation model to two other star formation recipes that used gas temperature and H$_2$ fraction as star formation criteria. We further followed the evolution of a dwarf galaxy within a cosmological context using both the shielding and H$_2$-based star formation models. We find that the shielding-based model allows for star formation at higher temperatures and lower densities than a model in which star formation is tied directly to H$_2$ abundance, as requiring H$_2$ formation leads the gas to undergo additional gravitational collapse before star formation. However, the resulting galaxies are very similar for both the shielding and H$_2$-based star formation models, and both models reproduce the resolved Kennicutt-Schmidt law. Therefore, both star formation models appear viable in the context of galaxy formation simulations.

We examine the effects of SMBH feedback on the CGM using a cosmological hydrodynamic simulation \citep[{\sc Romulus25};][]{Tremmel2017} and a set of four zoom-in `genetically modified' Milky Way-mass galaxies sampling different evolutionary paths. By tracing the distribution of metals in the circumgalactic medium (CGM), we show that \ion{O}{6} is a sensitive indicator of supermassive black hole (SMBH) feedback. First, we calculate the column densities of \ion{O}{6} in simulated Milky Way-mass galaxies and compare them with observations from the COS-Halos Survey. Our simulations show column densities of \ion{O}{6} in the CGM consistent with those of COS-Halos star forming and quenched galaxies. These results contrast with those from previous simulation studies which typically underproduce CGM column densities of \ion{O}{6}. We determine that a galaxy's star formation history and assembly record have little effect on the amount of \ion{O}{6} in its CGM. Instead, column densities of \ion{O}{6} are closely tied to galaxy halo mass and BH growth history. The set of zoom-in, genetically modified Milky Way-mass galaxies indicates that the SMBH drives highly metal-enriched material out into its host galaxy's halo which in turn elevates the column densities of \ion{O}{6} in the CGM.

Recent discoveries of massive black holes (MBHs) in dwarf galaxies suggest that they may have a more common presence than once thought. Systematic searches are revealing more candidates, but this process could be accelerated by predictions from simulations. We perform a study of several high-resolution, cosmological, zoom-in simulations focusing on dwarf galaxies that host massive black holes at z = 0, with the aim of determining when the black holes are most observable. Larger dwarf galaxies are more likely to host MBHs than those of lower mass. About 50% of the MBHs in dwarfs are not centrally located, but rather are wandering within a few kpc of the galaxy center. The accretion luminosities of MBHs in dwarfs are low throughout cosmic time, rendering them extremely difficult to detect. However, the merger history of these MBHs is optimal for gravitational wave detection by LISA.

The Local Group hosts a number of star-forming dwarf galaxies that show evidence of periods of little to no star formation. We use a suite of cosmological simulations to study how star formation is reignited in such galaxies. We focus on isolated galaxies at $z=0$ with halo masses between 9.2$\times$10$^8$ M$_\odot$ and 8.4$\times$10$^9$ M$_\odot$, where star formation is typically shut off by reionization or by supernova feedback. Nearly 20% of these simulated galaxies later restart star formation, due to interactions with streams of gas in the intergalactic medium, indicating that this mechanism is relatively common in this mass range and that many isolated dwarfs at $z=0$ may not have been isolated throughout their histories. The source of this gas is not necessarily cosmic filaments. Rather, the dwarfs interact with gas thrown off by nearby galaxy mergers or streams extending from neighboring galaxies. While high ram pressure interactions of this nature lead to stripping, the encounters that reignite star formation are low density and/or low velocity and thus low ram pressure, resulting in compression of the hot gas in the halos of our dwarfs. The gas mass bound up in hot halos can be substantial -- at least an order of magnitude greater than the mass contained in HI. Consequently, we find that dwarfs that have experienced reignition tend to be more HI-rich and have a higher M$_{HI}$/M$_{*}$ ratio at $z=0$ than galaxies with continuous star formation. Using this fact, we identify galaxies in the Local Volume that might have "gappy" star formation histories, and can be studied by the Hubble Space Telescope or the James Webb Space Telescope.

We analyze the metal accumulation in dwarf and spiral galaxies by following the history of metal enrichment and outflows in a suite of twenty high-resolution simulated galaxies. These simulations agree with the observed stellar and gas-phase mass-metallicity relation, an agreement that relies on large fractions of the produced metals escaping into the CGM. For instance, in galaxies with Mvir ~ 1e9.5 -- 1e10 solar masses, we find that about ~ 85% of the available metals are outside of the galactic disk at z = 0, although the fraction decreases to a little less than half in Milky Way-mass galaxies. In many cases, these metals are spread far beyond the virial radius. We analyze the metal deficit within the ISM and stars in the context of previous work tracking the inflow and outflow of baryons. Outflows are prevalent across the entire mass range, as is reaccretion. We find that between 40 and 80% of all metals removed from the galactic disk are later reaccreted. The outflows themselves are metal enriched relative to the ISM by a factor of 0.2 dex because of the correspondence between sites of metal enrichment and outflows. As a result, the metal mass loading factor scales as eta_metals \propto v_circ^-0.91, a somewhat shallower scaling than the total mass loading factor. We analyze the simulated galaxies within the context of analytic chemical evolution models by determining their net metal expulsion efficiencies, which encapsulate the rates of metal loss and reaccretion. We discuss these results in light of the inflow and outflow properties necessary for reproducing the mass-metallicity relation.

We study the birth of supermassive black holes from the direct collapse process and characterize the sites where these black hole seeds form. In the pre-reionization epoch, molecular hydrogen (H$_2$) is an efficient coolant, causing gas to fragment and form Population III stars, but Lyman-Werner radiation can suppress H$_2$ formation and allow gas to collapse directly into a massive black hole. The critical flux required to inhibit H$_2$ formation, $J_{\rm crit}$, is hotly debated, largely due to the uncertainties in the source radiation spectrum, H$_2$ self-shielding, and collisional dissociation rates. Here, we test the power of the direct collapse model in a self-consistent, time-dependant, non-uniform Lyman-Werner radiation field -- the first time such has been done in a cosmological volume -- using an updated version of the SPH+N-body tree code Gasoline with H$_2$ non-equilibrium abundance tracking, H$_2$ cooling, and a modern SPH implementation. We vary $J_{\rm crit} $ from $30$ to $10^3$ in units of $J_{21}$ to study how this parameter impacts the number of seed black holes and the type of galaxies which host them. We focus on black hole formation as a function of environment, halo mass, metallicity, and proximity of the Lyman-Werner source. Massive black hole seeds form more abundantly with lower $J_{\rm crit}$ thresholds, but regardless of $J_{\rm crit}$, these seeds typically form in halos that have recently begun star formation. Our results do not confirm the proposed atomic cooling halo pair scenario; rather black hole seeds predominantly form in low-metallicity pockets of halos which already host star formation.

We examine the scalings of galactic outflows with halo mass across a suite of twenty high-resolution cosmological zoom galaxy simulations covering halo masses from 10^9.5 - 10^12 M_sun. These simulations self-consistently generate outflows from the available supernova energy in a manner that successfully reproduces key galaxy observables including the stellar mass-halo mass, Tully-Fisher, and mass-metallicity relations. We quantify the importance of ejective feedback to setting the stellar mass relative to the efficiency of gas accretion and star formation. Ejective feedback is increasingly important as galaxy mass decreases; we find an effective mass loading factor that scales as v_circ^(-2.2), with an amplitude and shape that is invariant with redshift. These scalings are consistent with analytic models for energy-driven wind, based solely on the halo potential. Recycling is common: about half the outflow mass across all galaxy masses is later re-accreted. The recycling timescale is typically about 1 Gyr, virtually independent of halo mass. Recycled material is re-accreted farther out in the disk and with typically about 2-3 times more angular momentum. These results elucidate and quantify how the baryon cycle plausibly regulates star formation and alters the angular momentum distribution of disk material across the halo mass range where most of cosmic star formation occurs.

The velocity anisotropy parameter, beta, is a measure of the kinematic state of orbits in the stellar halo which holds promise for constraining the merger history of the Milky Way (MW). We determine global trends for beta as a function of radius from three suites of simulations, including accretion only and cosmological hydrodynamic simulations. We find that both types of simulations are consistent and predict strong radial anisotropy (<beta>~0.7) for Galactocentric radii greater than 10 kpc. Previous observations of beta for the MW's stellar halo claim a detection of an isotropic or tangential "dip" at r~20 kpc. Using the N-body+SPH simulations, we investigate the temporal persistence, population origin, and severity of "dips" in beta. We find dips in the in situ stellar halo are long-lived, while dips in the accreted stellar halo are short-lived and tied to the recent accretion of satellite material. We also find that a major merger as early as z~1 can result in a present day low (isotropic to tangential) value of beta over a wide range of radii and angular expanse. While all of these mechanisms are plausible drivers for the beta dip observed in the MW, in the simulations, each mechanism has a unique metallicity signature associated with it, implying that future spectroscopic surveys could distinguish between them. Since an accurate knowledge of beta(r) is required for measuring the mass of the MW halo, we note significant transient dips in beta could cause an overestimate of the halo's mass when using spherical Jeans equation modeling.

Within a Lambda Cold Dark Matter (LCDM) scenario, we use high resolution cosmological simulations spanning over four orders of magnitude in galaxy mass to understand the deficit of dwarf galaxies in observed velocity functions. We measure velocities in as similar a way as possible to observations, including generating mock HI data cubes for our simulated galaxies. We demonstrate that this apples-to-apples comparison yields an "observed" velocity function in agreement with observations, reconciling the large number of low-mass halos expected in a LCDM cosmological model with the low number of observed dwarfs at a given velocity. We then explore the source of the discrepancy between observations and theory, and conclude that the dearth of observed dwarf galaxies is primarily explained by two effects. The first effect is that galactic rotational velocities derived from the HI linewidth severely underestimate the maximum halo velocity. The second effect is that a large fraction of halos at the lowest masses are too faint to be detected by current galaxy surveys. We find that cored dark matter density profiles can contribute to the lower observed velocity of galaxies, but only for galaxies in which the velocity is measured interior to the size of the core (~3 kpc).

We analyze age-velocity dispersion relations (AVRs) from kinematics of individual stars in eight Local Group galaxies ranging in mass from Carina ($M_{*} \sim 10^{6}$) to M31 ($M_{*} \sim 10^{11}$). Observationally the $\sigma$ vs. stellar age trends can be interpreted as dynamical heating of the stars by GMCs, bars/spiral arms, or merging subhalos; alternatively the stars could have simply been born out of a more turbulent ISM at high redshift and retain that larger velocity dispersion till present day - consistent with recent IFU studies. To ascertain the dominant mechanism and better understand the impact of instabilities and feedback, we develop models based on observed SFHs of these Local Group galaxies in order to create an evolutionary formalism which describes the ISM velocity dispersion due to a galaxy's evolving gas fraction. These empirical models relax the common assumption that the stars are born from gas which has constant velocity dispersion at all redshifts. Using only the observed SFHs as input, the ISM velocity dispersion and a mid-plane scattering model fits the observed AVRs of low mass galaxies without fine tuning. Higher mass galaxies above $M_{vir} > 10^{11}$ need a larger contribution from latent dynamical heating processes (for example minor mergers), in excess of the ISM model. Using the SFHs we also find that supernovae feedback does not appear to be a dominant driver of the gas velocity dispersion compared to gravitational instabilities - at least for dispersions $\sigma \gtrsim 25$ km/s. Together our results point to stars being born with a velocity dispersion close to that of the gas at the time of their formation, with latent dynamical heating operating with a galaxy mass-dependent efficiency. These semi-empirical relations may help constrain the efficiency of feedback and its impact on the physics of disk settling in galaxy formation simulations.

We investigate whether the inclusion of baryonic physics influences the formation of thin, coherently rotating planes of satellites such as those seen around the Milky Way and Andromeda. For four Milky Way-mass simulations, each run both as dark matter-only and with baryons included, we are able to identify a planar configuration that significantly maximizes the number of plane satellite members. The maximum plane member satellites are consistently different between the dark matter-only and baryonic versions of the same run due to the fact that satellites are both more likely to be destroyed and to infall later in the baryonic runs. Hence, studying satellite planes in dark matter-only simulations is misleading, because they will be composed of different satellite members than those that would exist if baryons were included. Additionally, the destruction of satellites in the baryonic runs leads to less radially concentrated satellite distributions, a result that is critical to making planes that are statistically significant compared to a random distribution. Since all planes pass through the centre of the galaxy, it is much harder to create a plane from a random distribution if the satellites have a low radial concentration. We identify Andromeda's low radial satellite concentration as a key reason why the plane in Andromeda is highly significant. Despite this, when co-rotation is considered, none of the satellite planes identified for the simulated galaxies are as statistically significant as the observed planes around the Milky Way and Andromeda, even in the baryonic runs.

Cosmological simulations of galaxy formation can produce present-day galaxies with a large range of assembly and star formation histories. A detailed study of the metallicity evolution and star formation history of such simulations can assist in predicting LIGO-detectable compact object binary mergers. Recent simulations of compact binary evolution suggest the compact object merger rate depends sensitively on the progenitor's metallicity. Rare low-metallicity star formation during galaxy assembly can produce more detected compact binaries than typical star formation. Using detailed simulations of galaxy and chemical evolution, we determine how sensitively the compact binary populations of galaxies with similar present-day appearance depend on the details of their assembly. We also demonstrate by concrete example the extent to which dwarf galaxies overabundantly produce compact binary mergers, particularly binary black holes, relative to more massive galaxies. We discuss the implications for transient multimessenger astronomy with compact binary sources.

The latest generation of cosmological simulations are on the verge of being able to resolve the structure of bulges for the first time. Hence, we review the current state of bulge formation in cosmological simulations, and discuss open questions that can be addressed in the near future by simulators, with a particular focus on merger-driven bulge growth. Galaxy mergers have long been assumed to produce classical bulges in disk galaxies. Under this bulge-formation model, though, the high rates of mergers in Cold Dark Matter galaxy formation theory predict many more classical bulges than are observed. Furthermore, simulations of galaxy formation continue to generally produce too massive of bulges. Feedback offers a promising avenue for reducing merger-driven bulge growth by maintaining high gas fractions in galaxies and ejecting low-angular momentum gas driven to the centers of galaxies. After reviewing the results of relevant research that has been published to date, we use cosmological simulations to explore the ability of feedback to reduce or even prevent bulge growth during mergers. In dwarf galaxies, mergers actually reduce the central concentration of galaxies as the induced burst of star formation drives out low angular momentum material. This result shows the potential for feedback to reduce central mass growth. However, we also demonstrate that it is very difficult for current stellar feedback models to reproduce the small bulges observed in more massive disk galaxies like the Milky Way. We argue that feedback models need to be improved, or an additional source of feedback such as AGN is necessary to generate the required outflows.

The star formation histories (SFHs) of dwarf galaxies are thought to be \emph{bursty}, with large -- order of magnitude -- changes in the star formation rate on timescales similar to O-star lifetimes. As a result, the standard interpretations of many galaxy observables (which assume a slowly varying SFH) are often incorrect. Here, we use the SFHs from hydro-dynamical simulations to investigate the effects of bursty SFHs on sample selection and interpretation of observables and make predictions to confirm such SFHs in future surveys. First, because dwarf galaxies' star formation rates change rapidly, the mass-to-light ratio is also changing rapidly in both the ionizing continuum and, to a lesser extent, the non-ionizing UV continuum. Therefore, flux limited surveys are highly biased toward selecting galaxies in the \emph{burst} phase and very deep observations are required to detect all dwarf galaxies at a given stellar mass. Second, we show that a $\log_{10}[\nu L_{\nu}(1500{\rm \AA})/L_{{\rm H}\alpha}]>2.5$ implies a very recent quenching of star formation and can be used as evidence of stellar feedback regulating star formation. Third, we show that the ionizing continuum can be significantly higher than when assuming a constant SFH, which can affect the interpretation of nebular emission line equivalent widths and direct ionizing continuum detections. Finally, we show that a star formation rate estimate based on continuum measurements only (and not on nebular tracers such as the hydrogen Balmer lines) will not trace the rapid changes in star formation and will give the false impression of a star-forming main sequence with low dispersion.

We present robust constraints from the Sloan Digital Sky Survey (SDSS) on the shape and distribution of the dark matter halo within the Milky Way (MW). Using the number density distribution and kinematics of SDSS halo stars, we probe the dark matter distribution to heliocentric distances exceeding 10 kpc and galactocentric distances exceeding 20 kpc. Our analysis utilizes Jeans equations to generate two-dimensional acceleration maps throughout the volume; this approach is thoroughly tested on a cosmologically derived N-body+SPH simulation of a MW-like galaxy. We show that the known accelerations (gradients of the gravitational potential) can be successfully recovered in such a realistic system. Leveraging the baryonic gravitational potential derived by Bovy & Rix (2013), we show that the gravitational potential implied by the SDSS observations cannot be explained, assuming Newtonian gravity, by visible matter alone: the gravitational force experienced by stars at galactocentric distances of 20 kpc is as much as three times stronger than what can be attributed to purely visible matter. We also show that the SDSS data provide a strong constraint on the shape of the dark matter halo potential. Within galactocentric distances of 20 kpc, the dark matter halo potential is well described as an oblate halo with axis ratio qDM=0.7+/-0.1; this corresponds to an axis ratio qDM=0.4+/-0.1 for the dark matter density distribution. Because of our precise two-dimensional measurements of the acceleration of the halo stars, we can reject several MOND models as an explanation of the observed behavior.

We use high resolution Hydro$+$N-Body cosmological simulations to compare the assembly and evolution of a small field dwarf (stellar mass ~ 10$^{6-7}$ M$\odot$, total mass 10$^{10}$ M$\odot$ in $\Lambda$ dominated CDM and 2keV WDM cosmologies. We find that star formation (SF) in the WDM model is reduced and delayed by 1-2 Gyr relative to the CDM model, independently of the details of SF and feedback. Independent of the DM model, but proportionally to the SF efficiency, gas outflows lower the central mass density through `dynamical heating', such that all realizations have circular velocities $<$ 20kms at 500$~$pc, in agreement with local kinematic constraints. As a result of dynamical heating, older stars are less centrally concentrated than younger stars, similar to stellar population gradients observed in nearby dwarf galaxies. Introducing an important diagnostic of SF and feedback models, we translate our simulations into artificial color-magnitude diagrams and star formation histories in order to directly compare to available observations. The simulated galaxies formed most of their stars in many $\sim$10 Myr long bursts. The CDM galaxy has a global SFH, HI abundance and Fe/H and alpha-elements distribution well matched to current observations of dwarf galaxies. These results highlight the importance of directly including `baryon physics' in simulations when 1) comparing predictions of galaxy formation models with the kinematics and number density of local dwarf galaxies and 2) differentiating between CDM and non-standard models with different DM or power spectra.

The relation of central black hole mass and stellar spheroid velocity dispersion (the M-$\sigma$ relation) is one of the best-known and tightest correlations linking black holes and their host galaxies. There has been much scrutiny concerning the difficulty of obtaining accurate black hole measurements, and rightly so; however, it has been taken for granted that measurements of velocity dispersion are essentially straightforward. We examine five disk galaxies from cosmological SPH simulations and find that line-of-sight effects due to galaxy orientation can affect the measured $\sigma$ by 30%, and consequently black hole mass predictions by up to 1.0 dex. Face-on orientations correspond to systematically lower velocity dispersion measurements, while more edge-on orientations give higher velocity dispersions, due to contamination by disk stars when measuring line of sight quantities. We caution observers that the uncertainty of velocity dispersion measurements is at least 20 km/s, and can be much larger for moderate inclinations. This effect may account for some of the scatter in the locally measured M-$\sigma$ relation, particularly at the low-mass end. We provide a method for correcting observed $\sigma_{\rm los}$ values for inclination effects based on observable quantities.

We present a detailed comparison between the photometric properties of the bulges of two simulated galaxies and those of a uniform sample of observed galaxies. This analysis shows that the simulated galaxies have bulges with realistic surface brightnesses for their sizes and magnitude. These two field disc galaxies have rotational velocities ~ 100 km/s and were integrated to a redshift of zero in a fully cosmological Lambda cold dark matter context as part of high-resolution smoothed particle hydrodynamic simulations. We performed bulge-disc decompositions of the galaxies using artificial observations, in order to conduct a fair comparison to observations. We also dynamically decomposed the galaxies and compared the star formation histories of the bulges to those of the entire galaxies. These star formation histories showed that the bulges were primarily formed before z = 1 and during periods of rapid star formation. Both galaxies have large amounts of early star formation, which is likely related to the relatively high bulge-to-disc ratios also measured for them. Unlike almost all previous cosmological simulations, the realistically concentrated bulges of these galaxies do not lead to unphysically high rotational velocities, causing them to naturally lie along the observed Tully-Fisher relation.

We compare the central mass distribution of galaxies simulated with three different models of the interstellar medium (ISM) with increasing complexity: primordial (H+He) cooling down to 10^4K, additional cooling via metal lines and to lower temperatures, and molecular hydrogen (H_2) with shielding of atomic and molecular hydrogen, in addition to metal line cooling. In order to analyze the effect of these models, we follow the evolution of four field galaxies with V_peak < 120 km/s to a redshift of zero using high-resolution Smoothed Particle Hydrodynamic simulations in a fully cosmological LCDM context. The spiral galaxies produced in simulations with either primordial cooling or H_2 physics have realistic, rising rotation curves. In contrast, the simulations with metal line cooling and otherwise similar feedback and star formation produced galaxies with the peaked rotation curves typical of most previous LCDM simulations of spiral galaxies. The less-massive bulges and non-peaked rotation curves in the galaxies simulated with primordial cooling or H_2 are linked to changes in the angular momentum distribution of the baryons. These galaxies had smaller amounts of low-angular momentum baryons because of increased gas loss from stellar feedback. When there is only primordial cooling, the star forming gas is hotter and the feedback-heated gas cools more slowly than when metal line cooling is included and so requires less energy to be expelled. When H_2 is included, the accompanying shielding produces large amounts of clumpy, cold gas where H_2 forms. Star formation in clumpy gas results in more concentrated supernova feedback and greater efficiency of mass loss. The higher feedback efficiency causes a decrease of low-angular momentum material. (abridged)

We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.

We examine the pressure of the star-forming interstellar medium (ISM) of Milky-Way sized disk galaxies using fully cosmological SPH+N-body, high resolution simulations. These simulations include explicit treatment of metal-line cooling in addition to dust and self-shielding, $\mathrm{H_{2}}$ based star formation. The 4 simulated halos have masses ranging from a few times $10^{10}$ to nearly $10^{12}$ solar masses. Using a kinematic decomposition of these galaxies into present-day bulge and disk components, we find that the typical pressure of the star-forming ISM in the present-day bulge is higher than that in the present-day disk by an order of magnitude. We also find that pressure of the star-forming ISM at high redshift is on average, higher than ISM pressures at low redshift. This explains the why the bulge forms at higher pressures: the disk assembles at lower redshift, when the ISM is lower pressure and the bulge forms at high redshift, when the ISM is at higher pressure. If ISM pressure and IMF variation are tied together as suggested in studies like \cite{Conroy2012}, these results could indicate a time-dependent IMF in Milky-Way like systems, as well as a different IMF in the bulge and the disk.

We present a first study of the effect of local photoionising radiation on gas cooling in smoothed particle hydrodynamics simulations of galaxy formation. We explore the combined effect of ionising radiation from young and old stellar populations. The method computes the effect of multiple radiative sources using the same tree algorithm used for gravity, so it is computationally efficient and well resolved. The method foregoes calculating absorption and scattering in favour of a constant escape fraction for young stars to keep the calculation efficient enough to simulate the entire evolution of a galaxy in a cosmological context to the present day. This allows us to quantify the effect of the local photoionisation feedback through the whole history of a galaxy`s formation. The simulation of a Milky Way like galaxy using the local photoionisation model forms ~ 40 % less stars than a simulation that only includes a standard uniform background UV field. The local photoionisation model decreases star formation by increasing the cooling time of the gas in the halo and increasing the equilibrium temperature of dense gas in the disc. Coupling the local radiation field to gas cooling from the halo provides a preventive feedback mechanism which keeps the central disc light and produces slowly rising rotation curves without resorting to extreme feedback mechanisms. These preliminary results indicate that the effect of local photoionising sources is significant and should not be ignored in models of galaxy formation.

Using high resolution cosmological hydrodynamical simulations of Milky Way-massed disk galaxies, we demonstrate that supernovae feedback and tidal stripping lower the central masses of bright (-15 < M_V < -8) satellite galaxies. These simulations resolve high density regions, comparable to giant molecular clouds, where stars form. This resolution allows us to adopt a prescription for H_2 formation and destruction that ties star formation to the presence of shielded, molecular gas. Before infall, supernova feedback from the clumpy, bursty star formation captured by this physically motivated model leads to reduced dark matter (DM) densities and shallower inner density profiles in the massive satellite progenitors (Mvir > 10^9 Msun, Mstar > 10^7 Msun) compared to DM-only simulations. The progenitors of the lower mass satellites are unable to maintain bursty star formation histories, due to both heating at reionization and gas loss from initial star forming events, preserving the steep inner density profile predicted by DM-only simulations. After infall, tidal stripping acts to further reduce the central densities of the luminous satellites, particularly those that enter with cored dark matter halos, increasing the discrepancy in the central masses predicted by baryon+DM and DM-only simulations. We show that DM-only simulations, which neglect the baryonic effects described in this work, produce denser satellites with larger central velocities. We provide a simple correction to the central DM mass predicted for satellites by DM-only simulations. We conclude that DM-only simulations should be used with great caution when interpreting kinematic observations of the Milky Way's dwarf satellites.

We search for evidence of dark matter in the Milky Way by utilizing the stellar number density distribution and kinematics measured by the Sloan Digital Sky Survey (SDSS) to heliocentric distances exceeding ~10 kpc. We employ the cylindrically symmetric form of Jeans equations and focus on the morphology of the resulting acceleration maps, rather than the normalization of the total mass as done in previous, mostly local, studies. Jeans equations are first applied to a mock catalog based on a cosmologically derived N-body + SPH simulation, and the known acceleration (gradient of gravitational potential) is successfully recovered. The same simulation is also used to quantify the impact of dark matter on the total acceleration. We use Galfast, a code designed to quantitatively reproduce SDSS measurements and selection effects, to generate a synthetic stellar catalog. We apply Jeans equations to this catalog and produce two-dimensional maps of stellar acceleration. These maps reveal that in a Newtonian framework, the implied gravitational potential cannot be explained by visible matter alone. The acceleration experienced by stars at galactocentric distances of ~20 kpc is three times larger than what can be explained by purely visible matter. The application of an analytic method for estimating the dark matter halo axis ratio to SDSS data implies an oblate halo with q_DM = 0.47 +/- 0.14 within the same distance range. These techniques can be used to map the dark matter halo to much larger distances from the Galactic center using upcoming deep optical surveys, such as LSST.

We examine the present-day total stellar-to-halo mass (SHM) ratio as a function of halo mass for a new sample of simulated field galaxies using fully cosmological, LCDM, high resolution SPH + N-Body simulations.These simulations include an explicit treatment of metal line cooling, dust and self-shielding, H2 based star formation and supernova driven gas outflows. The 18 simulated halos have masses ranging from a few times 10^8 to nearly 10^12 solar masses. At z=0 our simulated galaxies have a baryon content and morphology typical of field galaxies. Over a stellar mass range of 2.2 x 10^3 to 4.5 x 10^10 solar masses, we find extremely good agreement between the SHM ratio in simulations and the present-day predictions from the statistical Abundance Matching Technique presented in Moster et al. (2012). This improvement over past simulations is due to a number systematic factors, each decreasing the SHM ratios: 1) gas outflows that reduce the overall SF efficiency but allow for the formation of a cold gas component 2) estimating the stellar masses of simulated galaxies using artificial observations and photometric techniques similar to those used in observations and 3) accounting for a systematic, up to 30 percent overestimate in total halo masses in DM-only simulations, due to the neglect of baryon loss over cosmic times. Our analysis suggests that stellar mass estimates based on photometric magnitudes can underestimate the contribution of old stellar populations to the total stellar mass, leading to stellar mass errors of up to 50 percent for individual galaxies. These results highlight the importance of using proper techniques to compare simulations with observations and reduce the perceived tension between the star formation efficiency in galaxy formation models and in real galaxies.

We discuss opportunities that may arise from subjecting high-multiplicity events in relativistic heavy ion collisions to an analysis similar to the one used in cosmology for the study of fluctuations of the Cosmic Microwave Background (CMB). To this end, we discuss examples of how pertinent features of heavy ion collisions including global characteristics, signatures of collective flow and event-wise fluctuations are visually represented in a Mollweide projection commonly used in CMB analysis, and how they are statistically analyzed in an expansion over spherical harmonic functions. If applied to the characterization of purely azimuthal dependent phenomena such as collective flow, the expansion coefficients of spherical harmonics are seen to contain redundancies compared to the set of harmonic flow coefficients commonly used in heavy ion collisions. Our exploratory study indicates, however, that these redundancies may offer novel opportunities for a detailed characterization of those event-wise fluctuations that remain after subtraction of the dominant collective flow signatures. By construction, the proposed approach allows also for the characterization of more complex collective phenomena like higher-order flow and other sources of fluctuations, and it may be extended to the characterization of phenomena of non-collective origin such as jets.

Although relatively common in the local Universe, only one grand-design spiral galaxy has been spectroscopically confirmed to lie at z>2 (HDFX 28; z=2.011), and may prove to be a major merger that simply resembles a spiral in projection. The rarity of spirals has been explained as a result of disks being dynamically 'hot' at z>2 which may instead favor the formation of commonly-observed clumpy structures. Alternatively, current instrumentation may simply not be sensitive enough to detect spiral structures comparable to those in the modern Universe. At redshifts <2, the velocity dispersion of disks decreases, and spiral galaxies are more numerous by z~1. Here we report observations of the grand design spiral galaxy Q2343-BX442 at z=2.18. Spectroscopy of ionized gas shows that the disk is dynamically hot, implying an uncertain origin for the spiral structure. The kinematics of the galaxy are consistent with a thick disk undergoing a minor merger, which can drive the formation of short-lived spiral structure. A duty cycle of < 100 Myr for such tidally-induced spiral structure in a hot massive disk is consistent with their rarity.

Motivated by the observed connection between molecular hydrogen (H2) and star formation, we present a method for tracking the non-equilibrium abundance and cooling processes of H2 and H2-based star formation in Smoothed Particle Hydrodynamic simulations. The local abundances of H2 are calculated by integrating over the hydrogen chemical network. This calculation includes the gas-phase and dust grain formation of H2, shielding of HI and H2, and photodissociation of H2 by Lyman-Werner radiation from nearby stellar populations. Because this model does not assume equilibrium abundances, it is particularly well suited for simulations that model low-metallicity environments, such as dwarf galaxies and the early Universe. We further introduce an explicit link between star formation and local H2 abundance. This link limits star formation to "star-forming regions," represented by areas with abundant H2. With this implementation, we determine the effect of H2 on star formation in a cosmological simulation of a dwarf galaxy. This simulation is the first cosmological simulation with non-equilibrium H2 abundances to be integrated to a redshift of zero or to include efficient SN feedback. We find that our simulations are consistent with the observed Tully-Fisher, global Kennicutt-Schmidt, and resolved Kennicutt-Schmidt relations. We find that the inclusion of shielding of both the atomic and molecular hydrogen and, to a lesser extent, the additional cooling from H2 at temperatures between 200 and 5000 K increases the amount of cold gas in the galaxies. The changes to the ISM result in an increased amount of cold, dense gas in the disk of the galaxy and the formation of a clumpier interstellar media (ISM). The explicit link between star formation and H2 and the clumpier ISM results in a bluer galaxy with a greater spatial distribution of star formation at a redshift of zero. (abridged)

We examine the evolution of the inner dark matter (DM) and baryonic density profile of a new sample of simulated field galaxies using fully cosmological, Lambda CDM, high resolution SPH + N-Body simulations. These simulations include explicit H2 and metal cooling, star formation (SF) and supernovae (SNe) driven gas outflows. Starting at high redshift, rapid, repeated gas outflows following bursty SF transfer energy to the DM component and significantly flatten the originally `cuspy' central DM mass profile of galaxies with present day stellar masses in the 10^4.5 -- 10^9.8 Msolar range. At z=0, the central slope of the DM density profile of our galaxies (measured between 0.3 and 0.7 kpc from their centre) is well fitted by rhoDM propto r^alpha with alpha \simeq -0.5 + 0.35 log_10(Mstar/10^8Msolar) where Mstar is the stellar mass of the galaxy and 4 < log_10 Mstar < 9.4. These values imply DM profiles flatter than those obtained in DM--only simulations and in close agreement with those inferred in galaxies from the THINGS and LITTLE THINGS survey. Only in very small halos, where by z=0 star formation has converted less than ~ 0.03% of the original baryon abundance into stars, outflows do not flatten the original cuspy DM profile out to radii resolved by our simulations. The mass (DM and baryonic) measured within the inner 500 pc of each simulated galaxy remains nearly constant over four orders of magnitudes in stellar mass for Mstar 10^9 Msolar. This finding is consistent with estimates for faint Local Group dwarfs and field galaxies. These results address one of the outstanding problems faced by the CDM model, namely the strong discrepancy between the original predictions of cuspy DM profiles and the shallower central DM distribution observed in galaxies.

We examine the effect of mass and force resolution on a specific star formation (SF) recipe using a set of N-body/Smooth Particle Hydrodynamic simulations of isolated galaxies. Our simulations span halo masses from 10^9 to 10^13 solar masses, more than four orders of magnitude in mass resolution, and two orders of magnitude in the gravitational softening length, epsilon, representing the force resolution. We examine the total global star formation rate, the star formation history, and the quantity of stellar feedback and compare the disk structure of the galaxies. Based on our analysis, we recommend using at least 10^4 particles each for the dark matter and gas component and a force resolution of epsilon approximately equal to 10^-3 R_vir when studying global SF and feedback. When the spatial distribution of stars is important, the number of gas and dark matter particles must be increased to at least 10^5 of each. Low mass resolution simulations with fixed softening lengths show particularly weak stellar disks due to two-body heating. While decreasing spatial resolution in low mass resolution simulations limits two-body effects, density and potential gradients cannot be sustained. Regardless of the softening, low-mass resolution simulations contain fewer high density regions where SF may occur. Galaxies of approximately 10^10 solar masses display unique sensitivity to both mass and force resolution. This mass of galaxy has a shallow potential and is on the verge of forming a disk. The combination of these factors give this galaxy the potential for strong gas outflows driven by supernova feedback and make it particularly sensitive to any changes to the simulation parameters.

Using high resolution SPH simulations in a fully cosmological Lambda CDM context we study the formation of a bright disk dominated galaxy that originates from a "wet" major merger at z=0.8. The progenitors of the disk galaxy are themselves disk galaxies that formed from early major mergers between galaxies with blue colors. A substantial thin stellar disk grows rapidly following the last major merger and the present day properties of the final remnant are typical of early type spiral galaxies, with an i band B/D ~0.65, a disk scale length of 7.2 kpc, g-r = 0.5 mag, an HI line width (W_{20}/2) of 238 km/sec and total magnitude i = -22.4. The key ingredients for the formation of a dominant stellar disk component after a major merger are: i) substantial and rapid accretion of gas through cold flows followed at late times by cooling of gas from the hot phase, ii) supernova feedback that is able to partially suppress star formation during mergers and iii) relative fading of the spheroidal component. The gas fraction of the progenitors' disks does not exceed 25% at z<3, emphasizing that the continuous supply of gas from the local environment plays a major role in the regrowth of disks and in keeping the galaxies blue. The results of this simulation alleviate the problem posed for the existence of disk galaxies by the high likelihood of interactions and mergers for galaxy sized halos at relatively low z.

The ACS Nearby Galaxy Survey Treasury (ANGST) is a systematic survey to establish a legacy of uniform multi-color photometry of resolved stars for a volume-limited sample of nearby galaxies (D<4 Mpc). The survey volume encompasses 69 galaxies in diverse environments, including close pairs, small & large groups, filaments, and truly isolated regions. The galaxies include a nearly complete range of morphological types spanning a factor of ~10^4 in luminosity and star formation rate. The survey data consists of images taken with ACS on HST, supplemented with archival data and new WFPC2 imaging taken after the failure of ACS. Survey images include wide field tilings covering the full radial extent of each galaxy, and single deep pointings in uncrowded regions of the most massive galaxies in the volume. The new wide field imaging in ANGST reaches median 50% completenesses of m_F475W=28.0 mag, m_F606W=27.3 mag, and m_F814W=27.3 mag, several magnitudes below the tip of the red giant branch (TRGB). The deep fields reach magnitudes sufficient to fully resolve the structure in the red clump. The resulting photometric catalogs are publicly accessible and contain over 34 million photometric measurements of >14 million stars. In this paper we present the details of the sample selection, imaging, data reduction, and the resulting photometric catalogs, along with an analysis of the photometric uncertainties (systematic and random), for both the ACS and WFPC2 imaging. We also present uniformly derived relative distances measured from the apparent magnitude of the TRGB.

In order to attract, recruit and retain underrepresented minority students to pursue Astronomy and related fields, we must ensure that there continues to be a well qualified pool of graduate and undergraduate students from which to recruit. This required pool of people are today's elementary, middle and high school students. The Astronomy community must be proactive in demonstrating the importance of pursing scientific study and careers to these students and their parents. Only by actively engaging these communities can U.S Astronomy hope to increase the numbers of minority PhDs and continue to be a leader in Astronomical discovery and knowledge.

If the ethnic makeup of the astronomy profession is to achieve parity with the general population within one generation (~30 years), the number of underrepresented minorities earning graduate degrees in astronomy and astrophysics must increase in the coming decade by a factor of 5 to 10. To accomplish this, the profession must develop and invest in mechanisms to more effectively move individuals across critical educational junctures to the PhD and beyond. Early and continuous research engagement starting in the undergraduate years is critical to this vision, in which the federally funded research internship programs (e.g. NSF REU, NASA GSRP) and national centers/observatories play a vital role. Regionally based partnerships with minority-serving institutions (MSIs) are crucial for tapping extant pools of minority talent, as are post-baccalaurate and/or masters degree "bridging" programs that provide critical stepping stones to the PhD. Because of the strong undergraduate physics, engineering, and computer science backgrounds of many students from MSIs, we suggest that instrument development and large scale computing/data-mining are particularly promising avenues for engagement in the coming decade.

Promoting racial and ethnic diversity is critically important to the future success and growth of the field of astronomy. The raw ability, drive and interest required to excel in the field is distributed without regard to race, gender, or socioeconomic background. By not actively promoting diversity in our field we risk losing talented people to other professions (or losing them entirely), which means that there will be astronomical discoveries that simply won't get made. There is demonstrated evidence that STEM fields benefit from diverse perspectives on problems that require more complex thought processes. This is especially relevant to a field like astronomy where more and more work is being done collaboratively. The lack of notable growth in African American, Hispanic, and Native American representation in astronomy indicates that the 'pipeline' for these individuals is systemically leaky at critical junctures. Substantially more effort must be directed toward improving the educational and career development of minorities to insure that these potential colleagues are supported through the process. However, simply recognizing that the pipeline is faulty is woefully inadequate. There must be very specific, targeted solutions to help improve the situation. With this in mind, we offer two position papers addressing specific areas of improvement that we identify as (a) essential for any foreseeable progress in the field, and (b) attainable in the 2010-2020 decade. These position papers focus primarily on African Americans, Hispanics, and Native Americans. Although we do not directly address issues of Asian Americans, Pacific Islanders, and other groups, many of the recommendations made here can be adapted to address issues faced by these groups as well.