We study the origins of 122 ultra-diffuse galaxies (UDGs) in the {\sc RomulusC} zoom-in cosmological simulation of a galaxy cluster (M$_{200} = 1.15\times10^{14}$ M$_{\odot}$), one of the only such simulations capable of resolving the evolution and structure of dwarf galaxies (M$_{\star} < 10^9$ M$_{\odot}$). We find broad agreement with observed cluster UDGs and predict that they are not separate from the overall cluster dwarf population. UDGs in cluster environments form primarily from dwarf galaxies that experienced early cluster in-fall and subsequent quenching due to ram pressure. The ensuing dimming of these dwarf galaxies due to passive stellar evolution results in a population of very low surface brightness galaxies that are otherwise typical dwarfs. UDGs and non-UDGs alike are affected by tidal interactions with the cluster potential. Tidal stripping of dark matter, as well as mass loss from stellar evolution, results in the adiabatic expansion of stars, particularly in the lowest mass dwarfs. High mass dwarf galaxies show signatures of tidal heating while low mass dwarfs that survive until $z=0$ typically have not experienced such impulsive interactions. There is little difference between UDGs and non-UDGs in terms of their dark matter halos, stellar morphology, colors, and location within the cluster. In most respects cluster UDG and non-UDGs alike are similar to isolated dwarf galaxies, except for the fact that they are typically quenched.

We use the Romulus25 cosmological simulation volume to identify the largest-ever simulated sample of field ultra-diffuse galaxies (UDGs). At z=0, we find that isolated UDGs have average star formation rates, colors, and virial masses for their stellar masses and environment. Although they tend to have slightly elevated HI masses, UDGs are consistent with the general isolated dwarf galaxy population and make up ~20% of all field galaxies with 10$^7$<M$_\star$/M$_\odot$<10$^{9}$. The properties of our UDGs are consistent with existing observations of field UDGs, but we predict that many isolated UDGs have been missed by current surveys. Despite their isolation at z=0, the UDGs in our sample are the products of major mergers at z>1. Mergers are no more common in UDG than non-UDG progenitors, but mergers that create UDGs tend to happen earlier, produce a temporary boost in spin, and cause star formation to migrate to the outskirts of galaxies, resulting in lower central star formation rates. The centers of the galaxies fade as their central stellar populations age, but their global star formation rates are maintained through bursts of star formation at larger radii. This formation channel is unique relative to other proposals for UDG formation in isolated galaxies, or for UDGs in groups and clusters, demonstrating that UDGs can potentially be formed through multiple mechanisms.

Using cosmological hydrodynamic simulations with physically motivated models of super-massive black hole (SMBH) formation and growth, we compare the assembly of Milky Way-mass ( $M_{\mathrm{vir}} \approx 7 \times 10^{11}$ $M_{\odot}$ at $z = 0$) galaxies in cold dark matter (CDM) and self-interacting dark matter (SIDM) models. Our SIDM model adopts a constant cross-section of 1 cm$^2$/g. We find that SMBH formation is suppressed in the early universe due to SIDM coring. SMBH-SMBH mergers are also suppressed in SIDM as a consequence of the lower number of SMBHs formed. Lack of initial merger-driven SMBH growth in turn delays SMBH growth by billions of years in SIDM compared to CDM. Further, we find that this delayed growth suppresses SMBH accretion in the largest progenitors of the main SIDM galaxies during the first 5 Gyrs of their evolution. Nonetheless, by $z = 0.8$ the CDM and SIDM SMBH masses differ only by around 0.2 dex, so that both remain compatible with the $M_{BH}-M_{*}$ relation. We show that the reduced accretion causes the SIDM SMBHs to less aggressively regulate star formation in their host galaxies than their CDM counterparts, resulting in a factor of 3 or more stars being produced over the lifetime of the SIDM galaxies compared to the CDM galaxies. Our results highlight a new way in which SIDM can affect the growth and merger history of SMBHs and ultimately give rise to very different galaxy evolution compared to the classic CDM model.

We construct a sample of 10,680 wall galaxies and 3,064 void galaxies with $M_R \lesssim -20$ by cross referencing a void catalog from literature with Baryon Oscillation Spectroscopic Survey (BOSS) CMASS and WiggleZ galaxies, where the CMASS survey targets redder galaxies and the WiggleZ survey targets bluer galaxies. Comparing the density profiles of the red and blue galaxies as a function of the void radius, we find that the number ratio of red-to-blue galaxies increases with distances from the void centers, suggesting a deficit of luminous and normal red galaxies in voids. We find a mean (g -- r) magnitude color of 1.298 and 1.210 for the wall and void galaxies, respectively, when considering the combined red and blue sample, which is found to be a significant difference. However, when considering the blue and red samples separately, we find no significant color difference. We conclude that the constituents galaxies of each population, rather than intrinsic color difference, is the main driver in the apparent average color difference of galaxies in voids and walls, indicating a deficit of luminous and normal red galaxies in voids. Our analysis suggests that the primary environmental-dependence effect on galaxy evolution for normal and luminous galaxies between void and wall regions is manifested in the number of red galaxies, which depends on the environmental-dependent merger history. Using a semi-analytic simulation model, we can successfully reproduce the apparent color difference between the void and wall galaxies.

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.

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.

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.

We predict the stellar mass-halo mass (SMHM) relationship for dwarf galaxies and their satellites residing in halos down to M$_{halo} =$ 10$^7$ M$_{\odot}$ with 10$^4$ M$_{\odot} <$ M$_{star}$($z=0$) $< 10^8$ M$_{\odot}$, and quantify the predicted scatter in the relation at the low mass end, using cosmological simulations. The galaxies were drawn from a cosmological simulation of dwarf galaxies, run with the N-body + SPH code, ChaNGA, at a high resolution of 60 pc. For M$_{halo} > 10^9$ M$_{\odot}$, the simulated SMHM relationship agrees with literature determinations, including exhibiting a small scatter. However, the scatter in the SMHM relation increases dramatically for lower-mass halos. We find that some of this scatter is due to {\em dark dwarfs}, halos devoid of stars. However, even when only considering well-resolved halos that contain a stellar population, the scatter in stellar mass reaches nearly 1 dex for M$_{halo}$($z=0$) 10$^7$ M$_{\odot}$. Much of this scatter is due to including satellites of the dwarf galaxies that have had their halo masses reduced through tidal stripping. The fraction of dark dwarfs (those that contain no stars) increases substantially with decreasing halo mass. When these dark halos are considered, the true scatter in the SMHM at low masses is even larger. At the faintest end of the SMHM relation probed by our simulations, a galaxy cannot be assigned a unique halo mass based solely on its luminosity. We provide a formula to stochastically populate low-mass halos following our results. Our predicted large scatter at low halo masses increases the slope of the resulting stellar mass function on the ultra-faint dwarf galaxy scales currently being probed by such surveys as the Dark Energy Survey or the Hyper-Suprime Cam Subaru Strategic Program, and in the future by the Large Synoptic Survey Telescope.

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 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 follow-up observations and analysis of the recently discovered short period low-mass eclipsing binary, SDSS J001641-000925. With an orbital period of 0.19856 days, this system has one of the shortest known periods for an M dwarf binary system. Medium-resolution spectroscopy and multi-band photometry for the system are presented. Markov chain Monte Carlo modeling of the light curves and radial velocities yields estimated masses for the stars of M1 = 0.54 +/- 0.07 Msun and M2 = 0.34 +/- 0.04 Msun, and radii of R1 = 0.68 +/- 0.03 Rsun and R2 = 0.58 +/- 0.03 Rsun respectively. This solution places both components above the critical Roche overfill limit, providing strong evidence that SDSS J001641-000925 is the first verified M-dwarf contact binary system. Within the follow-up spectroscopy we find signatures of non-solid body rotation velocities, which we interpret as evidence for mass transfer or loss within the system. In addition, our photometry samples the system over 9 years, and we find strong evidence for period decay at the rate of dP/dt ~8 s/yr. Both of these signatures raise the intriguing possibility that the system is in over-contact, and actively losing angular momentum, likely through mass loss. This places SDSS J001641-000925 as not just the first M-dwarf over-contact binary, but one of the few systems of any spectral type known to be actively undergoing coalescence. Further study SDSS J001641-000925 is on-going to verify the nature of the system, which may prove to be a unique astrophysical laboratory.

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 present ACS, NICMOS, and Keck AO-assisted photometry of 20 Type Ia supernovae SNe Ia from the HST Cluster Supernova Survey. The SNe Ia were discovered over the redshift interval 0.623 < z < 1.415. Fourteen of these SNe Ia pass our strict selection cuts and are used in combination with the world's sample of SNe Ia to derive the best current constraints on dark energy. Ten of our new SNe Ia are beyond redshift $z=1$, thereby nearly doubling the statistical weight of HST-discovered SNe Ia beyond this redshift. Our detailed analysis corrects for the recently identified correlation between SN Ia luminosity and host galaxy mass and corrects the NICMOS zeropoint at the count rates appropriate for very distant SNe Ia. Adding these supernovae improves the best combined constraint on the dark energy density \rho_{DE}(z) at redshifts 1.0 < z < 1.6 by 18% (including systematic errors). For a LambdaCDM universe, we find \Omega_\Lambda = 0.724 +0.015/-0.016 (68% CL including systematic errors). For a flat wCDM model, we measure a constant dark energy equation-of-state parameter w = -0.985 +0.071/-0.077 (68% CL). Curvature is constrained to ~0.7% in the owCDM model and to ~2% in a model in which dark energy is allowed to vary with parameters w_0 and w_a. Tightening further the constraints on the time evolution of dark energy will require several improvements, including high-quality multi-passband photometry of a sample of several dozen z>1 SNe Ia. We describe how such a sample could be efficiently obtained by targeting cluster fields with WFC3 on HST.

We present a new survey strategy to discover and study high redshift Type Ia supernovae (SNe Ia) using the Hubble Space Telescope (HST). By targeting massive galaxy clusters at 0.9<z<1.5, we obtain a twofold improvement in the efficiency of finding SNe compared to an HST field survey and a factor of three improvement in the total yield of SN detections in relatively dust-free red-sequence galaxies. In total, sixteen SNe were discovered at z>0.95, nine of which were in galaxy clusters. This strategy provides a SN sample that can be used to decouple the effects of host galaxy extinction and intrinsic color in high redshift SNe, thereby reducing one of the largest systematic uncertainties in SN cosmology.

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.

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.