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.

We study stellar-halo formation using six Milky Way-mass galaxies in FIRE-2 cosmological zoom simulations. We find that $5-40\%$ of the outer ($50-300$ kpc) stellar halo in each system consists of $\textit{in-situ}$ stars that were born in outflows from the main galaxy. Outflow stars originate from gas accelerated by super-bubble winds, which can be compressed, cool, and form co-moving stars. The majority of these stars remain bound to the halo and fall back with orbital properties similar to the rest of the stellar halo at $z=0$.In the outer halo, outflow stars are more spatially homogeneous, metal rich, and alpha-element-enhanced than the accreted stellar halo. At the solar location, up to $\sim 10 \%$ of our kinematically-identified halo stars were born in outflows; the fraction rises to as high as $\sim 40\%$ for the most metal-rich local halo stars ([Fe/H] $> -0.5$). We conclude that the Milky Way stellar halo could contain local counterparts to stars that are observed to form in molecular outflows in distant galaxies. Searches for such a population may provide a new, near-field approach to constraining feedback and outflow physics. A stellar halo contribution from outflows is a phase-reversal of the classic halo formation scenario of Eggen, Lynden-Bell $\&$ Sandange, who suggested that halo stars formed in rapidly $\textit{infalling}$ gas clouds. Stellar outflows may be observable in direct imaging of external galaxies and could provide a source for metal-rich, extreme velocity stars in the Milky Way.

We present a suite of FIRE-2 cosmological zoom-in simulations of isolated field dwarf galaxies, all with masses of $M_\mathrm{halo} \approx 10^{10}\,$M$_\odot$ at $z=0$, across a range of dark matter models. For the first time, we compare how both self-interacting dark matter (SIDM) and/or warm dark matter (WDM) models affect the assembly histories as well as the central density structure in fully hydrodynamical simulations of dwarfs. Dwarfs with smaller stellar half-mass radii (r$_{1/2}<500$ pc) have lower $\sigma_\star/V_\mathrm{max}$ ratios, reinforcing the idea that smaller dwarfs may reside in halos that are more massive than is naively expected. The majority of dwarfs simulated with self-interactions actually experience contraction of their inner density profiles with the addition of baryons relative to the cores produced in dark-matter-only runs, though the simulated dwarfs are always less centrally dense than in $\Lambda$CDM. The V$_{1/2}-$r$_{1/2}$ relation across all simulations is generally consistent with observations of Local Field dwarfs, though compact objects such as Tucana provide a unique challenge. Spatially-resolved rotation curves in the central regions ($<400$ pc) of small dwarfs could provide a way to distinguish between CDM, WDM, and SIDM, however: at the masses probed in this simulation suite, cored density profiles in dwarfs with small r$_{1/2}$ values can only originate from dark matter self-interactions.

We introduce an extension of the ELVIS project to account for the effects of the Milky Way galaxy on its subhalo population. Our simulation suite, Phat ELVIS, consists of twelve high-resolution cosmological dark matter-only (DMO) zoom simulations of Milky Way-size $\Lambda$CDM~ haloes ($M_{\rm v} = 0.7-2 \times 10^{12} \,\mathrm{M}_\odot$) along with twelve re-runs with embedded galaxy potentials grown to match the observed Milky Way disk and bulge today. The central galaxy potential destroys subhalos on orbits with small pericenters in every halo, regardless of the ratio of galaxy mass to halo mass. This has several important implications. 1) Most of the $\mathtt{Disk}$ runs have no subhaloes larger than $V_{\rm max} = 4.5$ km s$^{-1}$ within $20$ kpc and a significant lack of substructure going back $\sim 8$ Gyr, suggesting that local stream-heating signals from dark substructure will be rare. 2) The pericenter distributions of Milky Way satellites derived from $\mathit{Gaia}$ data are remarkably similar to the pericenter distributions of subhaloes in the $\mathtt{Disk}$ runs, while the DMO runs drastically over-predict galaxies with pericenters smaller than 20 kpc. 3) The enhanced destruction produces a tension opposite to that of the classic `missing satellites' problem: in order to account for ultra-faint galaxies known within $30$ kpc of the Galaxy, we must populate haloes with $V_\mathrm{peak} \simeq 7$ km s$^{-1}$ ($M \simeq 3 \times 10^{7} \,\mathrm{M}_\odot$ at infall), well below the atomic cooling limit of $V_\mathrm{peak} \simeq 16$ km s$^{-1}$ ($M \simeq 5 \times 10^{8} \,\mathrm{M}_\odot$ at infall). 4) If such tiny haloes do host ultra-faint dwarfs, this implies the existence of $\sim 1000$ satellite galaxies within 300 kpc of the Milky Way.

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.

One of the most robust observations of the stellar initial mass function (IMF) is its near-universality in the Milky Way and neighboring galaxies. But recent observations of early-type galaxies can be interpreted to imply a bottom-heavy IMF, while others of ultra-faint dwarfs could imply a top-heavy IMF. This would impose powerful constraints on star formation models. We explore what sort of cloud-scale IMF models could possibly satisfy these constraints. We utilize simulated galaxies which reproduce (broadly) the observed galaxy properties, while they also provide the detailed star formation history and properties of each progenitor star-forming cloud. We then consider generic models where the characteristic mass of the IMF is some arbitrary power-law function of progenitor cloud properties, along with well-known literature IMF models, which scale with Jeans mass, turbulent Bonnor-Ebert mass, temperature, the opacity limit, metallicity, or the protostellar heating mass. We show that no IMF models currently in the literature - nor any model where the turnover mass is an arbitrary power-law function of a combination of cloud temperature/density/size/metallicity/velocity dispersion/magnetic field - can reproduce the claimed IMF variation in ellipticals or dwarfs without severely violating observational constraints in the Milky Way. Specifically, they predict too much variation in the extreme environments of the Galaxy, compared to that observed. Either the IMF varies in a more complicated manner, or alternative interpretations of the extragalactic observations must be explored.

We explore the radial variation of star formation histories in dwarf galaxies simulated with Feedback In Realistic Environments (FIRE) physics. The sample contains 9 low-mass field dwarfs with M_ star = 10^5 - 10^7 M_sun from previous FIRE results, and a new suite of 17 higher mass field dwarfs with M_star = 10^7 - 10^9 M_sun introduced here. We find that age gradients are common in our dwarfs, with older stars dominant at large radii. The strength of the gradient correlates with overall galaxy age such that earlier star formation produces a more pronounced gradient. The relation between formation time and strength of the gradient is driven by both mergers and star-formation feedback. Mergers can both steepen and flatten the age gradient depending on the timing of the merger and star formation history of the merging galaxy. In galaxies without significant mergers, early feedback pushes stars to the outskirts at early times. Interestingly, among galaxies without mergers, those with large dark matter cores have flatter age gradients because these galaxies have more late-time feedback. If real galaxies have age gradients as we predict, stellar population studies that rely on sampling a limited fraction of a galaxy can give a biased view of its global star formation history. We show that central fields can be biased young by a few Gyrs while outer fields are biased old. Fields positioned near the 2D half-light radius will provide the least biased measure of a dwarf galaxy's global star formation history.

A standard prediction of galaxy formation theory is that the ionizing background suppresses galaxy formation in haloes with peak circular velocities smaller than Vpeak ~ 20 km/s, rendering the majority of haloes below this scale completely dark. We use a suite of cosmological zoom simulations of Milky Way-like haloes that include central Milky Way disk galaxy potentials to investigate the relationship between subhaloes and ultrafaint galaxies. We find that there are far too few subhaloes within 50 kpc of the Milky Way that had Vpeak > 20 km/s to account for the number of ultrafaint galaxies already known within that volume today. In order to match the observed count, we must populate subhaloes down to Vpeak ~ 6 km/s with ultrafaint dwarfs. The required haloes have peak virial temperatures as low as 1,500 K, well below the atomic hydrogen cooling limit of 10^4 K. Allowing for the possibility that the Large Magellanic Cloud contributes several of the satellites within 50 kpc could potentially raise this threshold to 10 km/s (4,000 K), still below the atomic cooling limit and far below the nominal reionization threshold.

A fundamental prediction of the cold dark matter cosmology is the existence of a large number of dark subhalos around galaxies, most of which should be entirely devoid of stars. Confirming the existence of dark substructures stands among the most important empirical challenges in modern cosmology: if they are found and quantified with the mass spectrum expected, then this would close the door on a vast array of competing theories. But in order for observational programs of this kind to reach fruition, we need robust predictions. Here we explore substructure predictions for lensing using galaxy lens-like hosts at z=0.2 from the Illustris simulations both in full hydrodynamics and dark matter only. We quantify substructures more massive than ~ 10^9 M_sun, comparable to current lensing detections derived from HST, Keck, and ALMA. The addition of full hydrodynamics reduces the overall subhalo mass function by about a factor of two. Even for the dark matter only runs, most (~ 85%) lines of sight through projected cylinders of size close to an Einstein radius contain no substructures larger than 10^9 M_sun. The fraction of empty sight lines rises to ~ 95% in full physics simulations. This suggests we will likely need hundreds of strong lensing systems suitable for substructure studies, as well as predictions that include the effects of baryon physics on substructure, to properly constrain cosmological models. Fortunately, the field is poised to fulfill these requirements.

Among the most important goals in cosmology is detecting and quantifying small ($M_{\rm halo}\simeq10^{6-9}~\mathrm{M}_\odot$) dark matter (DM) subhalos. Current probes around the Milky Way (MW) are most sensitive to such substructure within $\sim20$ kpc of the halo center, where the galaxy contributes significantly to the potential. We explore the effects of baryons on subhalo populations in $\Lambda$CDM using cosmological zoom-in baryonic simulations of MW-mass halos from the Latte simulation suite, part of the Feedback In Realistic Environments (FIRE) project. Specifically, we compare simulations of the same two halos run using (1) DM-only (DMO), (2) full baryonic physics, and (3) DM with an embedded disk potential grown to match the FIRE simulation. Relative to baryonic simulations, DMO simulations contain $\sim2\times$ as many subhalos within 100 kpc of the halo center; this excess is $\gtrsim5\times$ within 25 kpc. At $z=0$, the baryonic simulations are completely devoid of subhalos down to $3\times10^6~\mathrm{M}_\odot$ within $15$ kpc of the MW-mass galaxy, and fewer than 20 surviving subhalos have orbital pericenters <20 kpc. Despite the complexities of baryonic physics, the simple addition of an embedded central disk potential to DMO simulations reproduces this subhalo depletion, including trends with radius, remarkably well. Thus, the additional tidal field from the central galaxy is the primary cause of subhalo depletion. Subhalos on radial orbits that pass close to the central galaxy are preferentially destroyed, causing the surviving subhalo population to have tangentially biased orbits compared to DMO predictions. Our method of embedding a disk potential in DMO simulations provides a fast and accurate alternative to full baryonic simulations, thus enabling suites of cosmological simulations that can provide accurate and statistical predictions of substructure populations.

We investigate the response of self-interacting dark matter (SIDM) halos to the growth of galaxy potentials using idealized simulations, each run in tandem with standard collisionless Cold Dark Matter (CDM). We find a greater diversity in the SIDM halo profiles compared to the CDM halo profiles. If the stellar gravitational potential strongly dominates in the central parts of a galaxy, then SIDM halos can be as dense as CDM halos on observable scales. For extreme cases with highly compact disks core collapse can occur, leading to SIDM halos that are denser and cuspier than their CDM counterparts. If the stellar potential is not dominant, then SIDM halos retain constant density cores with densities far below CDM predictions. When a disk potential is present, the inner SIDM halo becomes \em{more flattened} in the disk plane than the CDM halo. These results are in excellent quantitative agreement with the predictions of Kaplinghat et al. (2014). We also simulated a galaxy cluster halo with a central stellar distribution similar to the brightest central galaxy of the cluster A2667. A SIDM halo simulated with cross section over mass $\sigma/m = 0.1\ \mathrm{cm^2 g^{-1}}$ provides a good match to the measured dark matter density profile of A2667, while an adiabatically-contracted CDM halo is denser and cuspier. The cored profile of the same halo simulated with $\sigma/m = 0.5\ \mathrm{cm^2 g^{-1}}$ is not dense enough to match A2667. Our findings are in agreement with previous results that $\sigma/m \gtrsim 0.1\ \mathrm{cm^2 g^{-1}}$ is disfavored for dark matter collision velocities in excess of about 1500 km/s. More generally, the predictive cross-talk between baryonic potentials and SIDM density distributions offers new directions for constraining SIDM cross sections in massive galaxies where baryons are dynamically important.

We present $29\pm1$ classical Oe stars from RIOTS4, a spatially complete, spectroscopic survey of Small Magellanic Cloud (SMC) field OB stars. The two earliest are O6e stars, and four are earlier than any Milky Way (MW) Oe stars. We also find ten Ope stars, showing He~\textsc{i} infill and/or emission; five appear to be at least as hot as $\sim$O7.5e stars. The hottest, star 77616, shows He~\textsc{ii} disk emission, suggesting that even the hottest O stars can form decretion disks, and offers observational support for theoretical predictions that the hottest, fastest rotators can generate He$^+$-ionizing atmospheres. Our data also demonstrate that Ope stars correspond to Oe stars earlier than O7.5e with strong disk emission. We find that in the SMC, Oe stars extend to earlier spectral types than in the MW, and our SMC Oe/O frequency, $0.26\pm0.04$, is much greater than the MW value, $0.03\pm0.01$. These results are consistent with angular momentum transport by stronger winds suppressing decretion disk formation at higher metallicity. In addition, our SMC field Oe star frequency is indistinguishable from that for clusters, which is consistent with the similarity between rotation rates in these environments, and contrary to the pattern for MW rotation rates. Thus, our findings strongly support the viscous decretion disk model and confirm that Oe stars are the high-mass extension of the Be phenomenon. Additionally, we find that Fe~\textsc{ii} emission occurs among Oe stars later than O7.5e with massive disks, and we revise a photometric criterion for identifying Oe stars to $J-[3.6] \geq 0.1$.

We present the Runaways and Isolated O-Type Star Spectroscopic Survey of the SMC (RIOTS4), a spatially complete survey of uniformly selected field OB stars that covers the entire star-forming body of the SMC. Using the IMACS multislit spectrograph and MIKE echelle spectrograph on the Magellan telescopes, we obtained spectra of 374 early-type field stars that are at least 28 pc from any other OB candidates. We also obtained spectra of an additional 23 field stars in the SMC bar identified from slightly different photometric criteria. Here, we present the observational catalog of stars in the RIOTS4 survey, including spectral classifications and radial velocities. For three multi-slit fields covering 8% of our sample, we carried out monitoring observations over 9-16 epochs to study binarity, finding a spectroscopic, massive binary frequency of at least $\sim$60% in this subsample. Classical Oe/Be stars represent a large fraction of RIOTS4 (42%), occurring at much higher frequency than in the Galaxy, consistent with expectation at low metallicity. RIOTS4 confirmed a steep upper IMF in the field, apparently caused by the inability of the most massive stars to form in the smallest clusters. Our survey also yields evidence for in-situ field OB star formation, and properties of field emission-line star populations, including sgB[e] stars and classical Oe/Be stars. We also discuss the radial velocity distribution and its relation to SMC kinematics and runaway stars. RIOTS4 presents a first quantitative characterization of field OB stars in an external galaxy, including the contributions of sparse, but normal, star formation; runaway stars; and candidate isolated star formation.

We constrain the evolution of the galaxy stellar mass function from 2 < z < 5 for galaxies with stellar masses as low as 10^5 Msun by combining star formation histories of Milky Way satellite galaxies derived from deep Hubble Space Telescope observations with merger trees from the ELVIS suite of N-body simulations. This approach extends our understanding more than two orders of magnitude lower in stellar mass than is currently possible by direct imaging. We find the faint end slopes of the mass functions to be alpha= -1.42(+0.07/-0.05) at z = 2 and alpha = -1.57^(+0.06/-0.06) at z = 5, and show the slope only weakly evolves from z = 5 to z = 0. Our findings are in stark contrast to a number of direct detection studies that suggest slopes as steep as alpha = -1.9 at these epochs. Such a steep slope would result in an order of magnitude too many luminous Milky Way satellites in a mass regime that is observationally complete (Mstar > 2*10^5 Msun at z = 0). The most recent studies from ZFOURGE and CANDELS also suggest flatter faint end slopes that are consistent with our results, but with a lower degree of precision. This work illustrates the strong connections between low and high-z observations when viewed through the lens of LCDM numerical simulations.

No classical Oe/Be stars with spectral type earlier than O7.5e have been identified to date in the Milky Way (MW). This is consistent with the decretion disk model because strong stellar winds cause early-type O stars to lose angular momentum, thereby preventing them from rotating fast enough to spin out decretion disks. How- ever, metal-poor O stars have weaker stellar winds, allowing the stars to retain angular momentum. Therefore, low-metallicity environments should promote the formation of Oe stars, including those of earlier spectral types than observed in high-metallicity en- vironments. Using the RIOTS4 survey, a spatially complete spectroscopic survey of Small Magellanic Cloud (SMC) field OB stars taken with the IMACS multi-slit spec- trograph at the Magellan Baade Telescope, we identify 25-31 SMC field Oe stars, which account for 20-28% of SMC field O stars. This fraction is significantly higher than in the MW, where < 10-15% of O stars display the Be phenomenon. We also present 5-7 Oe stars of spectral type ranging from O5.5e to O7e, all earlier spectral types than the earliest MW Oe star. These early type Oe stars represent 20-23% of our SMC Oe stars, a dramatic increase compared to the MW, where no Oe stars have been identified with these early spectral types. Thus, the higher frequencies of Oe stars and their earlier spectral range in the metal-poor SMC are consistent with the decretion disk model.

We present a sample of 14 OB stars in the Small Magellanic Cloud that meet strong criteria for having formed under extremely sparse star-forming conditions in the field. These stars are a minimum of 28 pc in projection from other OB stars, and they are centered within symmetric, round HII regions. They show no evidence of bow shocks, implying that the targets are not transverse runaway stars. Their radial velocities relative to local HI also indicate that they are not line-of-sight runaway stars. A friends-of-friends analysis shows that 9 of the objects present a few low-mass companion stars, with typical mass ratios for the two highest-mass stars of around 0.1. This further substantiates that these OB stars formed in place, and that they can and do form in extremely sparse conditions. This poses strong constraints on theories of star formation and challenges proposed relations between cluster mass and maximum stellar mass.

Some theories of star formation suggest massive stars may only form in clustered environments, which would create a deficit of massive stars in low density environments. Observationally, Massey (2002) finds such a deficit in samples of the field population in the Small and Large Magellanic Clouds, with an IMF slope of {\Gamma} ~ 4. These IMF measurements represent some of the largest known deviations from the standard Salpeter IMF slope of {\Gamma}=1.35. Here, we carry out a comprehensive investigation of the mass function above 20 solar masses for the entire field population of the Small Magellanic Cloud, based on data from the Runaways and Isolated O Type Star Spectroscopic Survey of the SMC (RIOTS4). This is a spatially complete census of the entire field OB star population of the SMC obtained with the IMACS multi-object spectrograph and MIKE echelle spectrograph on the Magellan telescopes. Based on Monte-Carlo simulations of the evolved present-day mass function, we find the slope of the field IMF above 20 solar masses is {\Gamma}=2.3+/-0.4. We extend our IMF measurement to lower masses using BV photometry from the OGLE II survey. We use a statistical approach to generate a probability distribution for the mass of each star from the OGLE photometry, and we again find {\Gamma}=2.3+/-0.6 for stellar masses from 7-20 solar masses. The discovery and removal of ten runaways in our RIOTS4 sample steepens the field IMF slope to {\Gamma}=2.8+/-0.5. We discuss the possible effects of binarity and star-formation history on our results, and conclude that the steep field massive star IMF is most likely a real effect.

We present the discovery of three new B[e] supergiants (sgB[e] stars) in the Small Magellanic Cloud (SMC). All three stars (R15, R38, and R48) were identified in the course of our Runaways and Isolated O Type Star Spectroscopic Survey of the SMC (RIOTS4). The stars show optical spectra that closely resemble those of previously known B[e] stars, presenting numerous low-ionization forbidden and permitted emission lines such as [Fe II] and Fe II. Furthermore, our stars have luminosities of log(L/L_sun) > 4, demonstrating that they are supergiants. However, we find lower infrared excesses and weaker forbidden emission lines than for previously identified B[e] supergiants. Thus our stars appear to either have less material in their circumstellar disks than other sgB[e] stars, or the circumstellar material has lower dust content. We suggest that these may constitute a new subclass of dust-poor sgB[e] stars.

The population of field OB stars are an important component of a galaxy's stellar content, representing 20-30% of the massive stars. To study this population, we have undertaken the Runaways and Isolated O Type Star Spectroscopic Survey of the SMC (RIOTS4). RIOTS4 surveys a spatially complete sample of >350 field OB stars in the Small Magellanic Cloud and will serve as a key probe of runaways, binaries, and the stellar IMF in the field massive star population. Here, we focus on the field IMF, which provides an empirical probe of the star-forming process and is a fundamental property of a stellar population. Together with photometry from the OGLE survey, RIOTS4 will yield a definitive stellar IMF for the SMC field massive star population. We present preliminary results that suggests the field IMF is much steeper, Gamma = 2.9, than the canonical stellar IMF of Gamma = 1.35. Despite the steep slope, we see no evidence of a stellar upper mass limit, up to our highest mass star of 65 solar masses.