We analyze the cold dark matter density profiles of 54 galaxy halos simulated with FIRE-2 galaxy formation physics, each resolved within $0.5\%$ of the halo virial radius. These halos contain galaxies with masses that range from ultra-faint dwarfs ($M_\star \simeq 10^{4.5} M_{\odot}$) to the largest spirals ($M_\star \simeq 10^{11} M_{\odot}$) and have density profiles that are both cored and cuspy. We characterize our results using a new analytic density profile that extends the standard Einasto form to allow for a pronounced constant-density core in the resolved innermost radius. With one additional core-radius parameter, $r_{c}$, this "core-Einasto" profile is able to characterize the shape and normalization of our feedback-impacted dark matter halos. In order to enable comparisons with observations, we provide fitting functions for $r_{c}$ and other profile parameters as a function of both $M_\star$ and $M_{\star}/M_{\rm halo}$. In agreement with similar studies done in the literature, we find that dark matter core formation is most efficient at the characteristic stellar-mass to halo-mass ratio $M_\star/M_{\rm halo} \simeq 5 \times 10^{-3}$, or $M_{\star} \sim 10^9 \, M_{\odot}$, with cores that are roughly the size of the galaxy half-light radius, $r_{c} \simeq 1-5$ kpc. Furthermore, we find no evidence for core formation at radii $\gtrsim 100\ \rm pc$ in galaxies with $M_{\star}/M_{\rm halo} < 5\times 10^{-4}$ or $M_\star \lesssim 10^6 \, M_{\odot}$. For Milky Way-size galaxies, baryonic contraction often makes halos significantly more concentrated and dense at the stellar half-light radius than dark matter only runs. However, even at the Milky Way scale, FIRE-2 galaxy formation still produces small dark matter cores of $\simeq 0.5-2$ kpc in size. Recent evidence for a ${\sim} 2$ kpc core in the Milky Way's dark matter halo is consistent with this expectation.

The past two decades have seen a tremendous investment in observational facilities that promise to reveal new and unprecedented discoveries about the universe. In comparison, the investment in theoretical work is completely dwarfed, even though theory plays a crucial role in the interpretation of these observations, predicting new types of phenomena, and informing observing strategies. In this white paper, we argue that in order to reach the promised critical breakthroughs in astrophysics over the next decade and well beyond, the national agencies must take a serious approach to investment in theoretical astrophysics research. We discuss the role of theory in shaping our understanding of the universe, and then we provide a multi-level strategy, from the grassroots to the national, to address the current underinvestment in theory relative to observational work.

Observations of low-mass satellite galaxies in the nearby Universe point towards a strong dichotomy in their star-forming properties relative to systems with similar mass in the field. Specifically, satellite galaxies are preferentially gas poor and no longer forming stars, while their field counterparts are largely gas rich and actively forming stars. Much of the recent work to understand this dichotomy has been statistical in nature, determining not just that environmental processes are most likely responsible for quenching these low-mass systems but also that they must operate very quickly after infall onto the host system, with quenching timescales $\lesssim 2~ {\rm Gyr}$ at ${M}_{\star} \lesssim 10^{8}~{\rm M}_{\odot}$. This work utilizes the newly-available $Gaia$ DR2 proper motion measurements along with the Phat ELVIS suite of high-resolution, cosmological, zoom-in simulations to study low-mass satellite quenching around the Milky Way on an object-by-object basis. We derive constraints on the infall times for $37$ of the known low-mass satellite galaxies of the Milky Way, finding that $\gtrsim~70\%$ of the `classical' satellites of the Milky Way are consistent with the very short quenching timescales inferred from the total population in previous works. The remaining classical Milky Way satellites have quenching timescales noticeably longer, with $\tau_{\rm quench} \sim 6 - 8~{\rm Gyr}$, highlighting how detailed orbital modeling is likely necessary to understand the specifics of environmental quenching for individual satellite galaxies. Additionally, we find that the $6$ ultra-faint dwarf galaxies with publicly available $HST$-based star-formation histories are all consistent with having their star formation shut down prior to infall onto the Milky Way -- which, combined with their very early quenching times, strongly favors quenching driven by reionization.

We study star formation histories (SFHs) of $\simeq500$ dwarf galaxies (stellar mass $M_\ast = 10^5 - 10^9\,M_\odot$) from FIRE-2 cosmological zoom-in simulations. We compare dwarfs around individual Milky Way (MW)-mass galaxies, dwarfs in Local Group (LG)-like environments, and true field (i.e. isolated) dwarf galaxies. We reproduce observed trends wherein higher-mass dwarfs quench later (if at all), regardless of environment. We also identify differences between the environments, both in terms of "satellite vs. central" and "LG vs. individual MWvs. isolated dwarf central." Around the individual MW-mass hosts, we recover the result expected from environmental quenching: central galaxies in the "near field" have more extended SFHs than their satellite counterparts, with the former more closely resemble isolated ("true field") dwarfs (though near-field centrals are still somewhat earlier forming). However, this difference is muted in the LG-like environments, where both near-field centrals and satellites have similar SFHs, which resemble satellites of single MW-mass hosts. This distinction is strongest for $M_\ast = 10^6 - 10^7\,M_\odot$ but exists at other masses. Our results suggest that the paired halo nature of the LG may regulate star formation in dwarf galaxies even beyond the virial radii of the MW and Andromeda. Caution is needed when comparing zoom-in simulations targeting isolated dwarf galaxies against observed dwarf galaxies in the LG.

We study a suite of extremely high-resolution cosmological FIRE simulations of dwarf galaxies ($M_{\rm halo} \lesssim 10^{10}$$M_{\odot}$), run to $z=0$ with $30 M_{\odot}$ resolution, sufficient (for the first time) to resolve the internal structure of individual supernovae remnants within the cooling radius. Every halo with $M_{\rm halo} \gtrsim 10^{8.6} M_{\odot}$ is populated by a resolved {\em stellar} galaxy, suggesting very low-mass dwarfs may be ubiquitous in the field. Our ultra-faint dwarfs (UFDs; $M_{\ast}<10^{5}\,M_{\odot}$) have their star formation truncated early ($z\gtrsim2$), likely by reionization, while classical dwarfs ($M_{\ast}>10^{5} M_{\odot}$) continue forming stars to $z<0.5$. The systems have bursty star formation (SF) histories, forming most of their stars in periods of elevated SF strongly clustered in both space and time. This allows our dwarf with $M_{\ast}/M_{\rm halo} > 10^{-4}$ to form a dark matter core $>200$pc, while lower-mass UFDs exhibit cusps down to $\lesssim100$pc, as expected from energetic arguments. Our dwarfs with $M_{\ast}>10^{4}\,M_{\odot}$ have half-mass radii ($R_{\rm 1/2}$) in agreement with Local Group (LG) dwarfs; dynamical mass vs. $R_{1/2}$ and the degree of rotational support also resemble observations. The lowest-mass UFDs are below surface brightness limits of current surveys but are potentially visible in next-generation surveys (e.g. LSST). The stellar metallicities are lower than in LG dwarfs; this may reflect pre-enrichment of the LG by the massive hosts or Pop-III stars. Consistency with lower resolution studies implies that our simulations are numerically robust (for a given physical model).

The Feedback In Realistic Environments (FIRE) project explores feedback in cosmological galaxy formation simulations. Previous FIRE simulations used an identical source code (FIRE-1) for consistency. Motivated by the development of more accurate numerics - including hydrodynamic solvers, gravitational softening, and supernova coupling algorithms - and exploration of new physics (e.g. magnetic fields), we introduce FIRE-2, an updated numerical implementation of FIRE physics for the GIZMO code. We run a suite of simulations and compare against FIRE-1: overall, FIRE-2 improvements do not qualitatively change galaxy-scale properties. We pursue an extensive study of numerics versus physics. Details of the star-formation algorithm, cooling physics, and chemistry have weak effects, provided that we include metal-line cooling and star formation occurs at higher-than-mean densities. We present new resolution criteria for high-resolution galaxy simulations. Most galaxy-scale properties are robust to numerics we test, provided: (1) Toomre masses are resolved; (2) feedback coupling ensures conservation, and (3) individual supernovae are time-resolved. Stellar masses and profiles are most robust to resolution, followed by metal abundances and morphologies, followed by properties of winds and circum-galactic media (CGM). Central (~kpc) mass concentrations in massive (L*) galaxies are sensitive to numerics (via trapping/recycling of winds in hot halos). Multiple feedback mechanisms play key roles: supernovae regulate stellar masses/winds; stellar mass-loss fuels late star formation; radiative feedback suppresses accretion onto dwarfs and instantaneous star formation in disks. We provide all initial conditions and numerical algorithms used.

We investigate stellar metallicity distribution functions (MDFs), including Fe and ${\alpha}$-element abundances, in dwarf galaxies from the Feedback in Realistic Environments (FIRE) project. We examine both isolated dwarf galaxies and those that are satellites of a Milky Way-mass galaxy. In particular, we study the effects of including a sub-grid turbulent model for the diffusion of metals in gas. Simulations that include diffusion have narrower MDFs and abundance ratio distributions, because diffusion drives individual gas and star particles toward the average metallicity. This effect provides significantly better agreement with observed abundance distributions of dwarf galaxies in the Local Group, including the small intrinsic scatter in [${\alpha}$/Fe] vs. [Fe/H] (less than 0.1 dex). This small intrinsic scatter arises in our simulations because the interstellar medium (ISM) in dwarf galaxies is well-mixed at nearly all cosmic times, such that stars that form at a given time have similar abundances to within 0.1 dex. Thus, most of the scatter in abundances at z = 0 arises from redshift evolution and not from instantaneous scatter in the ISM. We find similar MDF widths and intrinsic scatter for satellite and isolated dwarf galaxies, which suggests that environmental effects play a minor role compared with internal chemical evolution in our simulations. Overall, with the inclusion of metal diffusion, our simulations reproduce abundance distribution widths of observed low-mass galaxies, enabling detailed studies of chemical evolution in galaxy formation.

We present nebular phase optical and near-infrared spectroscopy of the Type Ia supernova (SN) 2017cbv. The early light curves of SN~2017cbv showed a prominent blue bump in the $U$, $B$ and $g$ bands lasting for $\sim$5 d. One interpretation of the early light curve was that the excess blue light was due to shocking of the SN ejecta against a nondegenerate companion star -- a signature of the single degenerate scenario. If this is the correct interpretation, the interaction between the SN ejecta and the companion star could result in significant H$\alpha$ (or helium) emission at late times, possibly along with other species, depending on the companion star and its orbital separation. A search for H$\alpha$ emission in our +302 d spectrum yields a nondetection, with a $L_{H\alpha}$$<$8.0$\times$10$^{35}$ erg/s (given an assumed distance of $D$=12.3 Mpc), which we have verified by implanting simulated H$\alpha$ emission into our data. We make a quantitative comparison to models of swept-up material stripped from a nondegenerate companion star, and limit the mass of hydrogen that might remain undetected to $M_{\rm H} < 1 \times 10^{-4}$ $\rm M_{\odot}$. A similar analysis of helium star related lines yields a $M_{\rm He} < 5 \times 10^{-4}$ $\rm M_{\odot}$. Taken at face value, these results argue against a nondegenerate H or He-rich companion in Roche lobe overflow as the progenitor of SN 2017cbv. Alternatively, there could be weaknesses in the envelope-stripping and radiative transfer models necessary to interpret the strong H and He flux limits.

In the local Universe, there is a strong division in the star-forming properties of low-mass galaxies, with star formation largely ubiquitous amongst the field population while satellite systems are predominantly quenched. This dichotomy implies that environmental processes play the dominant role in suppressing star formation within this low-mass regime (${M}_{\star} \sim 10^{5.5-8}~{\rm M}_{\odot}$). As shown by observations of the Local Volume, however, there is a non-negligible population of passive systems in the field, which challenges our understanding of quenching at low masses. By applying the satellite quenching models of Fillingham et al. (2015) to subhalo populations in the Exploring the Local Volume In Simulations (ELVIS) suite, we investigate the role of environmental processes in quenching star formation within the nearby field. Using model parameters that reproduce the satellite quenched fraction in the Local Group, we predict a quenched fraction -- due solely to environmental effects -- of $\sim 0.52 \pm 0.26$ within $1< R/R_{\rm vir} < 2$ of the Milky Way and M31. This is in good agreement with current observations of the Local Volume and suggests that the majority of the passive field systems observed at these distances are quenched via environmental mechanisms. Beyond $2~R_{\rm vir}$, however, dwarf galaxy quenching becomes difficult to explain through an interaction with either the Milky Way or M31, such that more isolated, field dwarfs may be self-quenched as a result of star-formation feedback.

Galaxies covering several orders of magnitude in stellar mass and a variety of Hubble types have been shown to follow the "Radial Acceleration Relation" (RAR), a relationship between $g_{\rm obs}$, the observed circular acceleration of the galaxy, and $g_{\rm bar}$, the acceleration due to the total baryonic mass of the galaxy. For accelerations above $10^{10}~{\rm m \, s}^{-2}$, $g_{\rm obs}$ traces $g_{\rm bar}$, asymptoting to the 1:1 line. Below this scale, there is a break in the relation such that $\rm g_{\rm obs} \sim g_{\rm bar}^{1/2}$. We show that the RAR slope, scatter and the acceleration scale are all natural consequences of the well-known baryonic Tully-Fisher relation (BTFR). We further demonstrate that galaxies with a variety of baryonic and dark matter (DM) profiles and a wide range of dark halo and galaxy properties (well beyond those expected in CDM) lie on the RAR if we simply require that their rotation curves satisfy the BTFR. We explore conditions needed to break this degeneracy: sub-kpc resolved rotation curves inside of "cored" DM-dominated profiles and/or outside $\gg 100\,$kpc could lie on BTFR but deviate in the RAR, providing new constraints on DM.

We investigate the merger histories of isolated dwarf galaxies based on a suite of 15 high-resolution cosmological zoom-in simulations, all with masses of $M_{\rm halo} \approx 10^{10}\,{\rm M}_{\odot}$ (and M$_\star\sim10^5-10^7\,{\rm M}_{\odot}$) at $z=0$, from the Feedback in Realistic Environments (FIRE) project. The stellar populations of these dwarf galaxies at $z=0$ are formed essentially entirely "in situ": over 90$\%$ of the stellar mass is formed in the main progenitor in all but two cases, and all 15 of the galaxies have >70$\%$ of their stellar mass formed in situ. Virtually all galaxy mergers occur prior to $z\sim3$, meaning that accreted stellar populations are ancient. On average, our simulated dwarfs undergo 5 galaxy mergers in their lifetimes, with typical pre-merger galaxy mass ratios that are less than 1:10. This merger frequency is generally comparable to what has been found in dissipationless simulations when coupled with abundance matching. Two of the simulated dwarfs have a luminous satellite companion at $z=0$. These ultra-faint dwarfs lie at or below current detectability thresholds but are intriguing targets for next-generation facilities. The small contribution of accreted stars make it extremely difficult to discern the effects of mergers in the vast majority of dwarfs either photometrically or using resolved-star color-magnitude diagrams (CMDs). The important implication for near-field cosmology is that star formation histories of comparably massive galaxies derived from resolved CMDs should trace the build-up of stellar mass in one main system across cosmic time as opposed to reflecting the contributions of many individual star formation histories of merged dwarfs.

We present a suite of 15 cosmological zoom-in simulations of isolated dark matter halos, all with masses of $M_{\rm halo} \approx 10^{10}\,{\rm M}_\odot$ at $z=0$, in order to understand the relationship between halo assembly, galaxy formation, and feedback's effects on the central density structure in dwarf galaxies. These simulations are part of the Feedback in Realistic Environments (FIRE) project and are performed at extremely high resolution. The resultant galaxies have stellar masses that are consistent with rough abundance matching estimates, coinciding with the faintest galaxies that can be seen beyond the virial radius of the Milky Way ($M_\star/{\rm M}_\odot\approx 10^5-10^7$). This non-negligible spread in stellar mass at $z=0$ in halos within a narrow range of virial masses is strongly correlated with central halo density or maximum circular velocity $V_{\rm max}$. Much of this dependence of $M_\star$ on a second parameter (beyond $M_{\rm halo}$) is a direct consequence of the $M_{\rm halo}\sim10^{10}\,{\rm M}_\odot$ mass scale coinciding with the threshold for strong reionization suppression: the densest, earliest-forming halos remain above the UV-suppression scale throughout their histories while late-forming systems fall below the UV-suppression scale over longer periods and form fewer stars as a result. In fact, the latest-forming, lowest-concentration halo in our suite fails to form any stars. Halos that form galaxies with $M_\star\gtrsim2\times10^{6}\,{\rm M}_\odot$ have reduced central densities relative to dark-matter-only simulations, and the radial extent of the density modifications is well-approximated by the galaxy half-mass radius $r_{1/2}$. This apparent stellar mass threshold of $M_\star \approx 2\times 10^{6} \approx 2\times 10^{-4} \,M_{\rm halo}$ is broadly consistent with previous work and provides a testable prediction of FIRE feedback models in LCDM.

We perform a systematic Bayesian analysis of rotation vs. dispersion support ($v_{\rm rot} / \sigma$) in $40$ dwarf galaxies throughout the Local Volume (LV) over a stellar mass range $10^{3.5} M_{\rm \odot} < M_{\star} < 10^8 M_{\rm \odot}$. We find that the stars in $\sim 80\%$ of the LV dwarf galaxies studied -- both satellites and isolated systems -- are dispersion-supported. In particular, we show that $6/10$ *isolated* dwarfs in our sample have $v_{\rm rot} / \sigma < 1.0$. All have $v_{\rm rot} / \sigma \lesssim 2.0$. These results challenge the traditional view that the stars in gas-rich dwarf irregulars (dIrrs) are distributed in cold, rotationally-supported stellar disks, while gas-poor dwarf spheroidals (dSphs) are kinematically distinct in having dispersion-supported stars. We see no clear trend between $v_{\rm rot} / \sigma$ and distance to the closest $\rm L_{\star}$ galaxy, nor between $v_{\rm rot} / \sigma$ and $M_{\star}$ within our mass range. We apply the same Bayesian analysis to four FIRE hydrodynamic zoom-in simulations of isolated dwarf galaxies ($10^9 M_{\odot} < M_{\rm vir} < 10^{10} M_{\rm \odot}$) and show that the simulated *isolated* dIrr galaxies have stellar ellipticities and stellar $v_{\rm rot} / \sigma$ ratios that are consistent with the observed population of dIrrs *and* dSphs without the need to subject these dwarfs to any external perturbations or tidal forces. We posit that most dwarf galaxies form as puffy, dispersion-dominated systems, rather than cold, angular momentum-supported disks. If this is the case, then transforming a dIrr into a dSph may require little more than removing its gas.

Recent studies of galaxies in the local Universe, including those in the Local Group, find that the efficiency of environmental (or satellite) quenching increases dramatically at satellite stellar masses below ~ $10^8\ {\rm M}_{\odot}$. This suggests a physical scale where quenching transitions from a slow "starvation" mode to a rapid "stripping" mode at low masses. We investigate the plausibility of this scenario using observed HI surface density profiles for a sample of 66 nearby galaxies as inputs to analytic calculations of ram-pressure and viscous stripping. Across a broad range of host properties, we find that stripping becomes increasingly effective at $M_{*} < 10^{8-9}\ {\rm M}_{\odot}$, reproducing the critical mass scale observed. However, for canonical values of the circumgalactic medium density ($n_{\rm halo} < 10^{-3.5}$ ${\rm cm}^{-3}$), we find that stripping is not fully effective; infalling satellites are, on average, stripped of < 40 - 70% of their cold gas reservoir, which is insufficient to match observations. By including a host halo gas distribution that is clumpy and therefore contains regions of higher density, we are able to reproduce the observed HI gas fractions (and thus the high quenched fraction and short quenching timescale) of Local Group satellites, suggesting that a host halo with clumpy gas may be crucial for quenching low-mass systems in Local Group-like (and more massive) host halos.

The vast majority of dwarf satellites orbiting the Milky Way and M31 are quenched, while comparable galaxies in the field are gas-rich and star-forming. Assuming that this dichotomy is driven by environmental quenching, we use the ELVIS suite of N-body simulations to constrain the characteristic timescale upon which satellites must quench following infall into the virial volumes of their hosts. The high satellite quenched fraction observed in the Local Group demands an extremely short quenching timescale (~ 2 Gyr) for dwarf satellites in the mass range Mstar ~ 10^6-10^8 Msun. This quenching timescale is significantly shorter than that required to explain the quenched fraction of more massive satellites (~ 8 Gyr), both in the Local Group and in more massive host halos, suggesting a dramatic change in the dominant satellite quenching mechanism at Mstar < 10^8 Msun. Combining our work with the results of complementary analyses in the literature, we conclude that the suppression of star formation in massive satellites (Mstar ~ 10^8 - 10^11 Msun) is broadly consistent with being driven by starvation, such that the satellite quenching timescale corresponds to the cold gas depletion time. Below a critical stellar mass scale of ~ 10^8 Msun, however, the required quenching times are much shorter than the expected cold gas depletion times. Instead, quenching must act on a timescale comparable to the dynamical time of the host halo. We posit that ram-pressure stripping can naturally explain this behavior, with the critical mass (of Mstar ~ 10^8 Msun) corresponding to halos with gravitational restoring forces that are too weak to overcome the drag force encountered when moving through an extended, hot circumgalactic medium.

We present FIRE/Gizmo hydrodynamic zoom-in simulations of isolated dark matter halos, two each at the mass of classical dwarf galaxies ($M_{\rm vir} \simeq 10^{10} M_{\odot}$) and ultra-faint galaxies ($M_{\rm vir} \simeq 10^9 M_{\odot}$), and with two feedback implementations. The resultant central galaxies lie on an extrapolated abundance matching relation from $M_{\star} \simeq 10^6$ to $10^4 M_{\odot}$ without a break. Every host is filled with subhalos, many of which form stars. Our dwarfs with $M_{\star} \simeq 10^6 M_{\odot}$ each have 1-2 well-resolved satellites with $M_{\star} = 3-200 \times 10^3 M_{\odot}$. Even our isolated ultra-faint galaxies have star-forming subhalos. If this is representative, dwarf galaxies throughout the universe should commonly host tiny satellite galaxies of their own. We combine our results with the ELVIS simulations to show that targeting $\sim 50~ \rm kpc$ regions around nearby isolated dwarfs could increase the chances of discovering ultra-faint galaxies by $\sim 35\%$ compared to random halo pointings, and specifically identify the region around the Phoenix dwarf galaxy as a good potential target. The well-resolved ultra-faint galaxies in our simulations ($M_{\star} \simeq 3 - 30 \times 10^3 M_{\odot}$) form within $M_{\rm peak} \simeq 0.5 - 3 \times 10^9 M_{\odot}$ halos. Each has a uniformly ancient stellar population ($ > 10~ \rm Gyr$) owing to reionization-related quenching. More massive systems, in contrast, all have late-time star formation. Our results suggest that $M_{\rm halo} \simeq 5 \times 10^9 M_{\odot}$ is a probable dividing line between halos hosting reionization "fossils" and those hosting dwarfs that can continue to form stars in isolation after reionization.

Using the Sloan Digital Sky Survey, we examine the quenching of satellite galaxies around isolated Milky Way-like hosts in the local Universe. We find that the efficiency of satellite quenching around isolated galaxies is low and roughly constant over two orders of magnitude in satellite stellar mass ($M_{*}$ = $10^{8.5}-10^{10.5} \, M_{\odot}$), with only $\sim~20\%$ of systems quenched as a result of environmental processes. While largely independent of satellite stellar mass, satellite quenching does exhibit clear dependence on the properties of the host. We show that satellites of passive hosts are substantially more likely to be quenched than those of star-forming hosts, and we present evidence that more massive halos quench their satellites more efficiently. These results extend trends seen previously in more massive host halos and for higher satellite masses. Taken together, it appears that galaxies with stellar masses larger than about $10^{8}~M_{\odot}$ are uniformly resistant to environmental quenching, with the relative harshness of the host environment likely serving as the primary driver of satellite quenching. At lower stellar masses ($< 10^{8}~M_{\odot}$), however, observations of the Local Group suggest that the vast majority of satellite galaxies are quenched, potentially pointing towards a characteristic satellite mass scale below which quenching efficiency increases dramatically.

We present cosmological constraints from a joint analysis of type Ia supernova (SN Ia) observations obtained by the SDSS-II and SNLS collaborations. The data set includes several low-redshift samples (z<0.1), all 3 seasons from the SDSS-II (0.05 < z < 0.4), and 3 years from SNLS (0.2 <z < 1) and totals \ntotc spectroscopically confirmed type Ia supernovae with high quality light curves. We have followed the methods and assumptions of the SNLS 3-year data analysis except for the following important improvements: 1) the addition of the full SDSS-II spectroscopically-confirmed SN Ia sample in both the training of the SALT2 light curve model and in the Hubble diagram analysis (\nsdssc SNe), 2) inter-calibration of the SNLS and SDSS surveys and reduced systematic uncertainties in the photometric calibration, performed blindly with respect to the cosmology analysis, and 3) a thorough investigation of systematic errors associated with the SALT2 modeling of SN Ia light-curves. We produce recalibrated SN Ia light-curves and associated distances for the SDSS-II and SNLS samples. The large SDSS-II sample provides an effective, independent, low-z anchor for the Hubble diagram and reduces the systematic error from calibration systematics in the low-z SN sample. For a flat LCDM cosmology we find Omega_m=0.295+-0.034 (stat+sys), a value consistent with the most recent CMB measurement from the Planck and WMAP experiments. Our result is 1.8sigma (stat+sys) different than the previously published result of SNLS 3-year data. The change is due primarily to improvements in the SNLS photometric calibration. When combined with CMB constraints, we measure a constant dark-energy equation of state parameter w=-1.018+-0.057 (stat+sys) for a flat universe. Adding BAO distance measurements gives similar constraints: w=-1.027+-0.055.

We study dwarf satellite galaxy quenching using observations from the Geha et al. (2012) NSA/SDSS catalog together with LCDM cosmological simulations to facilitate selection and interpretation. We show that fewer than 30% of dwarfs (M* ~ 10^8.5-10^9.5 Msun) identified as satellites within massive host halos (Mhost ~ 10^12.5-10^14 Msun) are quenched, in spite of the expectation from simulations that half of them should have been accreted more than 6 Gyr ago. We conclude that whatever the action triggering environmental quenching of dwarf satellites, the process must be highly inefficient. We investigate a series of simple, one-parameter quenching models in order to understand what is required to explain the low quenched fraction and conclude that either the quenching timescale is very long (> 9.5 Gyr, a "slow starvation" scenario) or that the environmental trigger is not well matched to accretion within the virial volume. We discuss these results in light of the fact that most of the low mass dwarf satellites in the Local Group are quenched, a seeming contradiction that could point to a characteristic mass scale for satellite quenching.

We examine the star formation properties of bright (~0.1 L*) satellites around isolated ~L* hosts in the local Universe using spectroscopically confirmed systems in the Sloan Digital Sky Survey DR7. Our selection method is carefully designed with the aid of N-body simulations to avoid groups and clusters. We find that satellites are significantly more likely to be quenched than a stellar mass-matched sample of isolated galaxies. Remarkably, this quenching occurs only for satellites of hosts that are themselves quenched: while star formation is unaffected in the satellites of star-forming hosts, satellites around quiescent hosts are more than twice as likely to be quenched than stellar-mass matched field samples. One implication of this is that whatever shuts down star formation in isolated, passive L* galaxies also plays at least an indirect role in quenching star formation in their bright satellites. The previously-reported tendency for "galactic conformity" in color/morphology may be a by-product of this host-specific quenching dichotomy. The S\'ersic indices of quenched satellites are statistically identical to those of field galaxies with the same specific star formation rates, suggesting that environmental and secular quenching give rise to the same morphological structure. By studying the distribution of pairwise velocities between the hosts and satellites, we find dynamical evidence that passive host galaxies reside in dark matter halos that are ~45% more massive than those of star-forming host galaxies of the same stellar mass. We emphasize that even around passive hosts, the mere fact that galaxies become satellites does not typically result in star formation quenching: we find that only ~30% of ~0.1 L* galaxies that fall in from the field are quenched around passive hosts, compared with ~0% around star forming hosts.

Observations with the Hubble Space Telescope (HST), conducted since 1990, now offer an unprecedented glimpse into fast astrophysical shocks in the young remnant of supernova 1987A. Comparing observations taken in 2010 using the refurbished instruments on HST with data taken in 2004, just before the Space Telescope Imaging Spectrograph failed, we find that the Ly-a and H-a lines from shock emission continue to brighten, while their maximum velocities continue to decrease. We observe broad blueshifted Ly-a, which we attribute to resonant scattering of photons emitted from hotspots on the equatorial ring. We also detect NV~\lambda\lambda 1239,1243 A line emission, but only to the red of Ly-A. The profiles of the NV lines differ markedly from that of H-a, suggesting that the N^{4+} ions are scattered and accelerated by turbulent electromagnetic fields that isotropize the ions in the collisionless shock.

Infrared carbon stars without visible counterparts are generally known as extreme carbon stars. We have selected a subset of these stars with absorption features in the 10-13 $\mu$m range, which has been tentatively attributed to silicon carbide (SiC). We add three new objects meeting these criterion to the seven previously known, bringing our total sample to ten sources. We also present the result of radiative transfer modeling for these stars, comparing these results to those of previous studies. In order to constrain model parameters, we use published mass-loss rates, expansion velocities and theoretical dust condensation models to determine the dust condensation temperature. These show that the inner dust temperatures of the dust shells for these sources are significantly higher than previously assumed. This also implies that the dominant dust species should be graphite instead of amorphous carbon. In combination with the higher condensation temperature we show that this results in a much higher acceleration of the dust grains than would be expected from previous work. Our model results suggest that the very optically thick stage of evolution does not coincide with the timescales for the superwind, but rather, that this is a very short-lived phase. Additionally, we compare model and observational parameters in an attempt to find any correlations. Finally, we show that the spectrum of one source, IRAS 17534$-$3030, strongly implies that the 10-13 $\mu$m feature is due to a solid state rather than a molecular species.

We use a series of N-body/smoothed particle hydrodynamics simulations and analytic arguments to show that the presence of an effective temperature floor in the interstellar medium at T_F ~ 10^4 K naturally explains the tendency for low-mass galaxies to be more spheroidal, more gas rich, and less efficient in converting baryons into stars than larger galaxies. The trend arises because gas pressure support becomes important compared to angular momentum support in small dark matter haloes. We suggest that dwarf galaxies with rotational velocities ~ 40 km/s do not originate as thin discs, but rather are born as thick, puffy systems. If accreted on to larger haloes, tenuous dwarfs of this kind will be more susceptible to gas loss or tidal transformation than scaled-down versions of larger spirals. For a constant temperature floor, pressure support becomes less important in large haloes, and this produces a tendency for massive isolated galaxies to have thinner discs and more efficient star formation than their less massive counterparts, as observed.

JEDI will probe dark energy in three independent ways by measuring the expansion history of the universe: (1) using type Ia supernovae as cosmological standard candles over a range of distances, (2) using baryon oscillations as a cosmological standard ruler over a range of cosmic epochs, (3) mapping the weak gravitational lensing distortion by foreground galaxies of the images of background galaxies at different distances. JEDI will unravel the nature of dark energy with accuracy and precision. JEDI is a 2m-class space telescope with the unique ability of simultaneous wide-field imaging (0.8-4.2 micron in five bands) and multi-slit spectroscopy (0.8-3.2 micron) with a field of view of 1 square degree. What makes JEDI efficient is its ability to simultaneously obtain high signal-to-noise ratio, moderate resolution slit spectra for all supernovae and ~ 5000 galaxies in its wide field of view, and to combine imaging and spectroscopy so that the appropriate balance of time is devoted to each. JEDI will measure the cosmic expansion history H(z) as a free function to < 2% accuracy in redshift bins of 0.2-0.3. Assuming a flat universe and \sigma(\Omega_m)=0.01 (0.03), JEDI could measure a departure from a vanilla \Lambda CDM model (w_0=-1, w'=0) with \sigma(w_0)=0.013 (0.031) and \sigma(w')=0.046 (0.063). JEDI will obtain the well-sampled lightcurves in Z, J, H, K, L bands and spectra of ~ 14,000 type Ia supernovae with redshifts ranging from 0 to 1.7; the redshifts of ~ 10-100 million galaxies to H ~ 23 and z ~ 4 over 1000-10,000 square degrees; and measurements of the shapes of galaxies over 1000-10,000 square degrees in Z,J,H,K,L for > 10^9 galaxies to H \~ 25.

GRB 000926 has one of the best-studied afterglows to-date, with multiple X-ray observations, as well as extensive multi-frequency optical and radio coverage. Broadband afterglow observations, spanning from X-ray to radio frequencies, provide a probe of the density structure of the circumburst medium, as well as of the ejecta energetics, geometry, and the physical parameters of the relativistic blastwave resulting from the explosion. We present an analysis of {\em Chandra X-ray Observatory} observations of this event, along with {\em Hubble Space Telescope} and radio monitoring. We combine these data with ground-based optical and IR observations and fit the synthesized afterglow lightcurve using models where collimated ejecta expand into a surrounding medium. We find that we can explain the broadband lightcurve with reasonable physical parameters only if the cooling is dominated by inverse Compton scattering. Excess X-ray emission in the broadband spectrum indicates that we are directly observing a contribution from inverse Compton scattering. It is the first time this has been observed in a GRB afterglow, and it implies that the GRB exploded in a reasonably dense (n~30 cm^{-3}) medium, consistent with a diffuse interstellar cloud environment.

Due to their high luminosity at maximum and degree of homogeneity, Type Ia supernovae have been extensively used for cosmological purpouses, in particular to estimate extragalactic distances and the Hubble constant. Recently the number of Type Ia supernovae detected at high redshift has increased, opening the possibility of determining the mass density parameter, the cosmological constant and the deceleration parameter. The observed supernovae appear to be further than expected -even for an empty Universe-, implying a low density Universe and moreover an accelerating Universe. Among the various uncertainties, we address the possibility that old supernovae are not equal to current supernovae. From first principles, an evolution of progenitors with time is expected. Additionally, some observations show a dependence of the observed properties on galaxy type and colour. Our aim in this work is to study the outcome of exploding CO white dwarfs following the evolution of the progenitor intermediate mass stars with different masses and metallicities. Once this influence of the progentitor has been determined, the observations may be corrected. At the present stage of this project we are not able to quantify this effect properly. One result is clear, that the differences at maximum are expected to be small ($\sim$0.2 mag) but this is of the same order as all the evidence for a positive cosmological constant ($\sim$0.25 mag).

Models of infrared spectra of Type Ia supernovae around maximum light are presented. The underlying dynamic models are delayed detonation explosions in Chandrasekhar mass carbon/oxygen white dwarfs. In combination with the radiative transport codes employed here, these models provide plausible fits to the optical spectra of "normal" Type Ia supernova. Two independent radiative transport codes are used, one that assumes LTE and one that computes non-LTE excitations and ionization. The models are compared with infrared data available in the literature. The independent codes give a reasonable representation of the data and provide physical explanations for their origin independent of the detailed assumptions of the radiative transfer. The infrared gives an especially powerful diagnostic of the dynamic model because with strongly variable line blanketing opacity it probes different depths within the exploded white dwarf at the same epoch. The velocity of the transition zone between explosive oxygen and carbon burning can be directly determined. The velocity at which the burning to nickel stops can also be probed. These velocities are very sensitive to the explosion physics.