The circumgalactic medium (CGM) of the Milky Way is mostly obscured by nearby gas in position-velocity space because we reside inside the Galaxy. Substantial biases exist in most studies on the Milky Way's CGM that focus on easier-to-detect high-velocity gas. With mock observations on a Milky-Way analog from the FOGGIE simulation, we investigate four observational biases related to the Milky Way's CGM. First, QSO absorption-line studies probe a limited amount of the CGM mass: only 35% of the mass is at high Galactic latitudes $|b|>20$ degrees, of which only half is moving at $|v_{\rm LSR}|\gtrsim100$ km s$^{-1}$. Second, the inflow rate ($\dot{M}$) of the cold gas observable in HI 21cm is reduced by a factor of $\sim10$ as we switch from the local standard of rest to the galaxy's rest frame; meanwhile $\dot{M}$ of the cool and warm gas does not change significantly. Third, OVI and NV are promising ions to probe the Milky Way's outer CGM ($r\gtrsim$15 kpc), but CIV may be less sensitive. Lastly, the scatter in ion column density is a factor of 2 higher if the CGM is observed from inside-out than from external views because of the gas radial density profile. Our work highlights that observations of the Milky Way's CGM, especially those using HI 21cm and QSO absorption lines, are highly biased. We demonstrate that these biases can be quantified and calibrated through synthetic observations with simulated Milky-Way analogs.

We study ram pressure stripping in simulated Milky Way-like halos at z>=2 from the Figuring Out Gas & Galaxies In Enzo (FOGGIE) project. These simulations reach exquisite resolution in their circumgalactic medium (CGM) gas owing to FOGGIE's novel refinement scheme. The CGM of each halo spans a wide dynamic range in density and velocity over its volume---roughly 6 dex and 1000 km/s, respectively---translating into a 5 dex range in ram pressure imparted to interacting satellites. The ram pressure profiles of the simulated CGM are highly stochastic, owing to kpc-scale variations of the density and velocity fields of the CGM gas. As a result, the efficacy of ram pressure stripping depends strongly on the specific path a satellite takes through the CGM. The ram-pressure history of a single satellite is generally unpredictable and not well correlated with its approach vector with respect to the host galaxy. The cumulative impact of ram pressure on the simulated satellites is dominated by only a few short strong impulses---on average, 90% of the total surface momentum gained through ram pressure is imparted in 20% or less of the total orbital time. These results reveal an erratic mode of ram pressure stripping in Milky-Way like halos at high redshift---one that is not captured by a smooth spherically-averaged model of the circumgalactic medium.

Maximizing the public impact of astronomy projects in the next decade requires NSF-funded centers to support the development of online, mobile-friendly outreach and education activities. EPO teams with astronomy, education, and web development expertise should be in place to build accessible programs at scale and support astronomers doing outreach.

A major goal for the astronomy and astrophysics communities is the pursuit of diversity, equity, and inclusion (DEI) in all ranks, from students through professional scientific researchers. Large scientific collaborations - increasingly a primary place for both professional interactions and research opportunities - can play an important role in the DEI effort. Multimessenger astronomy, a new and growing field, is based on the principle that working collaboratively produces synergies, enabling advances that would not be possible without cooperation. The nascent Multimessenger Diversity Network (MDN) is extending this collaborative approach to include DEI initiatives. After we review of the current state of DEI in astronomy and astrophysics, we describe the strategies the MDN is developing and disseminating to support and increase DEI in the fields over the coming decade: provide opportunities (real and virtual) to share DEI knowledge and resources, include DEI in collaboration-level activities, including external reviews, and develop and implement ways to recognize the DEI work of collaboration members.

We present simulations from the new "Figuring Out Gas & Galaxies in Enzo" (FOGGIE) project. In contrast to most extant simulations of galaxy formation, which concentrate computational resources on galactic disks and spheroids with fluid and particle elements of fixed mass, the FOGGIE simulations focus on extreme spatial and mass resolution in the circumgalactic medium (CGM) surrounding galaxies. Using the Enzo code and a new refinement scheme, FOGGIE reaches spatial resolutions of 381 comoving $h^{-1}$ pc and resolves extremely low masses ($\lesssim 1$--$100$ Msun out to 100 comoving $h^{-1}$ kpc from the central halo. At these resolutions, cloud and filament-like structures giving rise to simulated absorption are smaller, and better resolved, than the same structures simulated with standard density-dependent refinement. Most of the simulated absorption arises in identifiable and well-resolved structures with masses $\lesssim 10^4$ Msun, well below the mass resolution of typical zoom simulations. However, integrated quantities such as mass surface density and ionic covering fractions change at only the $\lesssim 30$% level as resolution is varied. This relatively small changes in projected quantities---even when the sizes and distribution of absorbing clouds change dramatically---indicate that commonly used observables provide only weak constraints on the physical structure of the underlying gas. Comparing the simulated absorption features to the KODIAQ (Keck Observatory Database of Ionized Absorption toward Quasars) survey of $z \sim2$--$3.5$ Lyman limit systems, we show that high-resolution FOGGIE runs better resolve the internal kinematic structure of detected absorption, and better match the observed distribution of absorber properties. These results indicate that CGM resolution is key in properly testing simulations of galaxy evolution with circumgalactic observations.

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.

Observing the circumgalactic medium (CGM) in emission provides 3D maps of the spatial and kinematic extent of the gas that fuels galaxies and receives their feedback. We present mock emission-line maps of highly resolved CGM gas from the FOGGIE project (Figuring Out Gas & Galaxies in Enzo) and link these maps back to physical and spatial properties of the gas. By increasing the spatial resolution alone, the total luminosity of the line emission increases by an order of magnitude. This increase arises in the abundance of dense small-scale structure resolved when the CGM gas is simulated to < 100 pc scales. Current integral field unit instruments like KCWI and MUSE should be able to detect the brightest knots and filaments of such emission, and from this to infer the bulk kinematics of the CGM gas with respect to the galaxy. We conclude that accounting for small-scale structure well below the level of instrument spatial resolution is necessary to properly interpret such observations in terms of the underlying gas structure driving observable emission.

Ultrafaint dwarf galaxies (UFDs) are typically assumed to have simple, stellar populations with star formation ending at reionization. Yet as the observations of these galaxies continue to improve, their star formation histories (SFHs) are revealed to be more complicated than previously thought. In this paper, we study how star formation, chemical enrichment, and mixing proceed in small, dark matter halos at early times using a high-resolution, cosmological, hydrodynamical simulation. The goals are to inform the future use of analytic models and to explore observable properties of the simulated halos in the context of UFD data. Specifically, we look at analytic approaches that might inform metal enrichment within and beyond small galaxies in the early Universe. We find that simple assumptions for modeling the extent of supernova-driven winds agree with the simulation on average whereas inhomogeneous mixing and gas flows have a large effect on the spread in simulated stellar metallicities. In the context of the UFDs, this work demonstrates that simulations can form halos with a complex SFH and a large spread in the metallicity distribution function within a few hundred Myr in the early Universe. In particular, bursty and continuous star formation are seen in the simulation and both scenarios have been argued from the data. Spreads in the simulated metallicities, however remain too narrow and too metal-rich when compared to the UFDs. Future work is needed to help reduce these discrepancies and advance our interpretation of the data.

We present a deep search for HI 21-cm emission from the gaseous halo of Messier 31 as part of Project AMIGA, a large program Hubble Space Telescope program to study the circumgalactic medium of the Andromeda galaxy. Our observations with the Robert C. Byrd Green Bank Telesope target sight lines to 48 background AGNs, more than half of which have been observed in the ultraviolet with the Cosmic Origins Spectrograph, with impact parameters $25 \lesssim \rho \lesssim 330$ kpc ($0.1 \lesssim \rho / R_{\rm vir} \lesssim 1.1$). We do not detect any 21-cm emission toward these AGNs to limits of $N({\rm HI}) \approx 4 \times10^{17}$ cm$^{-2}$ ($5\sigma$, per 2 kpc diameter beam). This column density corresponds to an optical depth of $\sim2.5$ at the Lyman limit, thus our observations overlap with absorption line studies of Lyman limit systems at higher redshift. Our non-detections place a limit on the covering factor of such optically-thick gas around M31 to $f_c < 0.051$ (at 90\% confidence) for $\rho \leq R_{\rm vir}$. While individual clouds have previously been found in the region between M31 and M33, the covering factor of strongly optically-thick gas is quite small. Our upper limits on the covering factor are consistent with expectations from recent cosmological "zoom" simulations. Recent COS-Halos ultraviolet measurements of \HI\ absorption about an ensemble of galaxies at $z \approx 0.2$ show significantly higher covering factors within $\rho \lesssim 0.5 R_{\rm vir}$ at the same $N({\rm H I})$, although the metal ion-to-H I ratios appear to be consistent with those seen in M31.

The circumgalactic medium (CGM) remains one of the least constrained components of galaxies and as such has significant potential for advancing galaxy formation theories. In this work, we vary the extragalactic ultraviolet background for a high-resolution cosmological simulation of a Milky Way-like galaxy and examine the effect on the absorption and emission properties of metals in the CGM. We find that a reduced quasar background brings the column density predictions into better agreement with recent data. Similarly, when the observationally derived physical properties of the gas are compared to the simulation, we find that the simulation gas is always at temperatures approximately 0.5 dex higher. Thus, similar column densities can be produced from fundamentally different gas. However, emission maps can provide complementary information to the line-of-sight column densities to better derive gas properties. From the simulations, we find that the brightest emission is less sensitive to the extragalactic background and that it closely follows the fundamental filamentary structure of the halo. This becomes increasingly true as the galaxy evolves from z=1 to z=0 and the majority of the gas transitions to a hotter, more diffuse phase. For the brightest ions (CIII, CIV, OVI), detectable emission can extend as far as 120 kpc at z=0. Finally, resolution is a limiting factor for the conclusions we can draw from emission observations but with moderate resolution and reasonable detection limits, upcoming instrumentation should place constraints on the physical properties of the CGM.

We explore the impact of cosmic rays (CRs) on cosmological adaptive-mesh refinement simulations of a forming 10^12 Msolar halo, focusing on the circumgalactic medium (CGM), and its resulting low-redshift structure and composition. In contrast to a run with star formation and energetic feedback but no CRs, the CR-inclusive runs feature a CGM substantially enriched with CRs and with metals to roughly 0.1 Zsolar, thanks to robust, persistent outflows from the disk. The CR-inclusive CGMs also feature more diffuse gas at lower temperatures, down to 10^4 K, than the non-CR run, with diffuse material often receiving a majority of its pressure support from the CR proton fluid. We compare to recent observations of the CGM of L ~ L* galaxies at low redshift, including UV absorption lines within background quasar spectra. The combination of metal-enriched, CR-driven winds and large swaths of CR pressure-supported, cooler diffuse gas leads to a CGM that provides a better match to data from COS-Halos (for HI, SiIV, CIII and OVI) than the non-CR run. We also compare our models to recent, preliminary observations of diffuse gamma-ray emission in local group halos. For our lower CR-diffusion runs with kappa_CR in the range 0.3 to 1 x 10^28 cm^2/s, the CR enriched CGM produces an inconsistently high level of gamma emission. But the model with a relatively high kappa_CR = 3 x 10^28 cm^2/s provided a gamma-ray luminosity consistent with the extra-galactic gamma-ray background observed by FERMI and roughly consistent with preliminary measures of the emission from M31's CGM.

Recent observations suggest that abundance pattern differences exist between low metallicity stars in the Milky Way stellar halo and those in the dwarf satellite galaxies. This paper takes a first look at what role the early environment for pre-galactic star formation might have played in shaping these stellar populations. In particular, we consider whether differences in cross-pollution between the progenitors of the stellar halo and the satellites could help to explain the differences in abundance patterns. Using an N-body simulation, we find that the progenitor halos of the main halo are primarily clustered together at z=10 while the progenitors of the satellite galaxies remain on the outskirts of this cluster. Next, analytically modeled supernova-driven winds show that main halo progenitors cross-pollute each other more effectively while satellite galaxy progenitors remain more isolated. Thus, inhomogeneous cross-pollution as a result of different high-z spatial locations of each system's progenitors can help to explain observed differences in abundance patterns today. Conversely, these differences are a signature of the inhomogeneity of metal enrichment at early times.

Consistency between Carnegie Supernova Project (CSP) and SDSS-II supernova (SN) survey ugri measurements has been evaluated by comparing SDSS and CSP photometry for nine spectroscopically confirmed Type Ia supernova observed contemporaneously by both programs. The CSP data were transformed into the SDSS photometric system. Sources of systematic uncertainty have been identified, quantified, and shown to be at or below the 0.023 magnitude level in all bands. When all photometry for a given band is combined, we find average magnitude differences of equal to or less than 0.011 magnitudes in ugri, with rms scatter ranging from 0.043 to 0.077 magnitudes. The u band agreement is promising, with the caveat that only four of the nine supernovae are well-observed in u and these four exhibit an 0.038 magnitude supernova-to-supernova scatter in this filter.