We present new measurements of the spatial distribution and kinematics of neutral hydrogen in the circumgalactic and intergalactic medium surrounding star-forming galaxies at z ~ 2. Using the spectra of ~ 3000 galaxies with redshifts <z> +/- 0.4 from the Keck Baryonic Structure Survey (KBSS), we assemble a sample of more than 200,000 distinct foreground-background pairs with projected angular separations of 3 - 500 arcsec and spectroscopic redshifts, with <$z_{fg}$> = 2.23 and <$z_{bg}$> = 2.57. The ensemble of sightlines and foreground galaxies is used to construct a 2D map of the mean excess Ly$\alpha$ optical depth relative to the intergalactic mean as a function of projected galactocentric distance (20 < $D_{tran}$/pkpc < 4000) and line-of-sight velocity. We provide information on the line-of-sight kinematics of H I gas as a function of projected distance $D_{tran}$. We compare the map with cosmological zoom-in simulation, finding qualitative agreement between them. A simple two-component (accretion, outflow) analytical model generally reproduces the observed line-of-sight kinematics and projected spatial distribution of H I. The best-fitting model suggests that galaxy-scale outflows with initial velocity $v_{out}$ ~ 600 km/s dominate the kinematics of circumgalactic H I out to $D_{tran}$ ~ 50 kpc, while H I at $D_{tran}$ > 100 kpc is dominated by infall with characteristic $v_{in}$ < $v_c$, where $v_c$ is the circular velocity of the host halo ($M_h$ ~ $10^{12} M_\odot$). Over the impact parameter range 80 < $D_{tran}$/pkpc < 200, the H I line-of-sight velocity range reaches a minimum, with a corresponding flattening in the rest-frame Ly$\alpha$ equivalent width. These observations can be naturally explained as the transition between outflow-dominated and accretion-dominated flows. Beyond $D_{tran}$ ~ 300 kpc, the line of sight kinematics are dominated by Hubble expansion.

We present a suite of high-resolution cosmological simulations, using the FIRE-2 feedback physics together with explicit treatment of magnetic fields, anisotropic conduction and viscosity, and cosmic rays (CRs) injected by supernovae (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultra-faint dwarf ($M_{\ast}\sim 10^{4}\,M_{\odot}$, $M_{\rm halo}\sim 10^{9}\,M_{\odot}$) through Milky Way masses, systematically vary CR parameters (e.g. the diffusion coefficient $\kappa$ and streaming velocity), and study an ensemble of galaxy properties (masses, star formation histories, mass profiles, phase structure, morphologies). We confirm previous conclusions that magnetic fields, conduction, and viscosity on resolved ($\gtrsim 1\,$pc) scales have small effects on bulk galaxy properties. CRs have relatively weak effects on all galaxy properties studied in dwarfs ($M_{\ast} \ll 10^{10}\,M_{\odot}$, $M_{\rm halo} \lesssim 10^{11}\,M_{\odot}$), or at high redshifts ($z\gtrsim 1-2$), for any physically-reasonable parameters. However at higher masses ($M_{\rm halo} \gtrsim 10^{11}\,M_{\odot}$) and $z\lesssim 1-2$, CRs can suppress star formation by factors $\sim 2-4$, given relatively high effective diffusion coefficients $\kappa \gtrsim 3\times10^{29}\,{\rm cm^{2}\,s^{-1}}$. At lower $\kappa$, CRs take too long to escape dense star-forming gas and lose energy to hadronic collisions, producing negligible effects on galaxies and violating empirical constraints from $\gamma$-ray emission. But around $\kappa\sim 3\times10^{29}\,{\rm cm^{2}\,s^{-1}}$, CRs escape the galaxy and build up a CR-pressure-dominated halo which supports dense, cool ($T\ll 10^{6}$ K) gas that would otherwise rain onto the galaxy. CR heating (from collisional and streaming losses) is never dominant.

We study the effects of cosmic rays (CRs) on outflows from star-forming galaxies in the circum and inter-galactic medium (CGM/IGM), in high-resolution, fully-cosmological FIRE-2 simulations (accounting for mechanical and radiative stellar feedback, magnetic fields, anisotropic conduction/viscosity/CR diffusion and streaming, and CR losses). We showed previously that massive ($M_{\rm halo}\gtrsim 10^{11}\,M_{\odot}$), low-redshift ($z\lesssim 1-2$) halos can have CR pressure dominate over thermal CGM pressure and balance gravity, giving rise to a cooler CGM with an equilibrium density profile. This dramatically alters outflows. Absent CRs, high gas thermal pressure in massive halos "traps" galactic outflows near the disk, so they recycle. With CRs injected in supernovae as modeled here, the low-pressure halo allows "escape" and CR pressure gradients continuously accelerate this material well into the IGM in "fast" outflows, while lower-density gas at large radii is accelerated in-situ into "slow" outflows that extend to $>$Mpc scales. CGM/IGM outflow morphologies are radically altered: they become mostly volume-filling (with inflow in a thin mid-plane layer) and coherently biconical from the disk to $>$Mpc. The CR-driven outflows are primarily cool ($T\sim10^{5}\,$K) and low-velocity. All of these effects weaken and eventually vanish at lower halo masses ($\lesssim 10^{11}\,M_{\odot}$) or higher redshifts ($z\gtrsim 1-2$), reflecting the ratio of CR to thermal+gravitational pressure in the outer halo. We present a simple analytic model which explains all of the above phenomena.

We explore the survival of cool clouds in multi-phase circum-galactic media. We revisit the "cloud crushing problem" in a large survey of simulations including radiative cooling, self-shielding, self-gravity, magnetic fields, and anisotropic Braginskii conduction and viscosity (with saturation). We explore a wide range of parameters including cloud size, velocity, ambient temperature and density, as well as a variety of magnetic field configurations and cloud turbulence. We find that realistic magnetic fields and turbulence have weaker effects on cloud survival; the most important physics is radiative cooling and conduction. Self-gravity and self-shielding are important for clouds which are initially Jeans-unstable, but largely irrelevant otherwise. Non-self-gravitating, realistically magnetized clouds separate into four regimes: (1) At low column densities, clouds evaporate rapidly via conduction. (2) A "failed pressure confinement" regime, where the ambient hot gas cools too rapidly to provide pressure confinement for the cloud. (3) An "infinitely long-lived" regime, in which the cloud lifetime becomes longer than the cooling time of gas swept up in the leading bow shock, so the cloud begins to accrete and grow. (4) A "classical cloud destruction" regime, where clouds are eventually destroyed by instabilities. In the final regime, the cloud lifetime can exceed the naive cloud-crushing time owing to conduction-induced compression. However, small and/or slow-moving clouds can also evaporate more rapidly than the cloud-crushing time. We develop simple analytic models that explain the simulated cloud destruction times in this regime.

We address the spatial scale, ionization structure, mass and metal content of gas at the Milky Way disk-halo interface detected as absorption in the foreground of seven closely-spaced, high-latitude halo blue horizontal branch stars (BHBs) with heights z = 3 - 14 kpc. We detect transitions that trace multiple ionization states (e.g. CaII, FeII, SiIV, CIV) with column densities that remain constant with height from the disk, indicating that the gas most likely lies within z < 3.4 kpc. The intermediate ionization state gas traced by CIV and SiIV is strongly correlated over the full range of transverse separations probed by our sightlines, indicating large, coherent structures greater than 1 kpc in size. The low ionization state material traced by CaII and FeII does not exhibit a correlation with either N$_{\rm HI}$ or transverse separation, implying cloudlets or clumpiness on scales less than 10 pc. We find that the observed ratio log(N_SiIV/ N_CIV), with a median value of -0.69+/-0.04, is sensitive to the total carbon content of the ionized gas under the assumption of either photoionization or collisional ionization. The only self-consistent solution for photoionized gas requires that Si be depleted onto dust by 0.35 dex relative to the solar Si/C ratio, similar to the level of Si depletion in DLAs and in the Milky Way ISM. The allowed range of values for the areal mass infall rate of warm, ionized gas at the disk-halo interface is 0.0003 < dM_gas / dtdA [M_sun kpc^-2 yr^-] < 0.006. Our data support a physical scenario in which the Milky Way is fed by complex, multiphase processes at its disk-halo interface that involve kpc-scale ionized envelopes or streams containing pc-scale, cool clumps.

We investigate the impact of cosmic rays (CRs) on the circumgalactic medium (CGM) in FIRE-2 simulations, for ultra-faint dwarf through Milky Way (MW)-mass halos hosting star-forming (SF) galaxies. Our CR treatment includes injection by supernovae, anisotropic streaming and diffusion along magnetic field lines, collisional and streaming losses, with constant parallel diffusivity $\kappa\sim3\times10^{29}\,\mathrm{cm^2\ s^{-1}}$ chosen to match $\gamma$-ray observations. With this, CRs become more important at larger halo masses and lower redshifts, and dominate the pressure in the CGM in MW-mass halos at $z\lesssim 1-2$. The gas in these ``CR-dominated'' halos differs significantly from runs without CRs: the gas is primarily cool (a few $\sim10^{4}\,$K), and the cool phase is volume-filling and has a thermal pressure below that needed for virial or local thermal pressure balance. Ionization of the ``low'' and ``mid'' ions in this diffuse cool gas is dominated by photo-ionization, with O VI columns $\gtrsim 10^{14.5}\,\mathrm{cm^{-2}}$ at distances $\gtrsim 150\,\mathrm{kpc}$. CR and thermal gas pressure are locally anti-correlated, maintaining total pressure balance, and the CGM gas density profile is determined by the balance of CR pressure gradients and gravity. Neglecting CRs, the same halos are primarily warm/hot ($T\gtrsim 10^{5}\,$K) with thermal pressure balancing gravity, collisional ionization dominates, O VI columns are lower and Ne VIII higher, and the cool phase is confined to dense filaments in local thermal pressure equilibrium with the hot phase.

Software is a critical part of modern research, and yet there are insufficient mechanisms in the scholarly ecosystem to acknowledge, cite, and measure the impact of research software. The majority of academic fields rely on a one-dimensional credit model whereby academic articles (and their associated citations) are the dominant factor in the success of a researcher's career. In the petabyte era of astronomical science, citing software and measuring its impact enables academia to retain and reward researchers that make significant software contributions. These highly skilled researchers must be retained to maximize the scientific return from petabyte-scale datasets. Evolving beyond the one-dimensional credit model requires overcoming several key challenges, including the current scholarly ecosystem and scientific culture issues. This white paper will present these challenges and suggest practical solutions for elevating the role of software as a product of the research enterprise.

We use a particle tracking analysis to study the origins of the circumgalactic medium (CGM), separating it into (1) accretion from the intergalactic medium (IGM), (2) wind from the central galaxy, and (3) gas ejected from other galaxies. Our sample consists of 21 FIRE-2 simulations, spanning the halo mass range log(Mh/Msun) ~ 10-12 , and we focus on z=0.25 and z=2. Owing to strong stellar feedback, only ~L* halos retain a baryon mass >~50% of their cosmic budget. Metals are more efficiently retained by halos, with a retention fraction >~50%. Across all masses and redshifts analyzed >~60% of the CGM mass originates as IGM accretion (some of which is associated with infalling halos). Overall, the second most important contribution is wind from the central galaxy, though gas ejected or stripped from satellites can contribute a comparable mass in ~L* halos. Gas can persist in the CGM for billions of years, resulting in well-mixed halo gas. Sight lines through the CGM are therefore likely to intersect gas of multiple origins. For low-redshift ~L* halos, cool gas (T<10^4.7 K) is distributed on average preferentially along the galaxy plane, however with strong halo-to-halo variability. The metallicity of IGM accretion is systematically lower than the metallicity of winds (typically by >~1 dex), although CGM and IGM metallicities depend significantly on the treatment of subgrid metal diffusion. Our results highlight the multiple physical mechanisms that contribute to the CGM and will inform observational efforts to develop a cohesive picture.

The majority of baryons reside beyond the optical extent of a galaxy in the circumgalactic and intergalactic media (CGM/IGM). Gaseous halos are inextricably linked to the appearance of their host galaxies through a complex story of accretion, feedback, and continual recycling. The energetic processes, which define the state of gas in the CGM, are the same ones that 1) regulate stellar growth so that it is not over-efficient, and 2) create the diversity of today's galaxy colors, SFRs, and morphologies spanning Hubble's Tuning Fork Diagram. They work in concert to set the speed of growth on the star-forming Main Sequence, transform a galaxy across the Green Valley, and maintain a galaxy's quenched appearance on the Red Sequence. Most baryons in halos more massive than 10^12 Msolar along with their high-energy physics and dynamics remain invisible because that gas is heated above the UV ionization states. We argue that information on many of the essential drivers of galaxy evolution is primarily contained in this "missing" hot gas phase. Completing the picture of galaxy formation requires uncovering the physical mechanisms behind stellar and SMBH feedback driving mass, metals, and energy into the CGM. By opening galactic hot halos to new wavebands, we not only obtain fossil imprints of >13 Gyrs of evolution, but observe on-going hot-mode accretion, the deposition of superwind outflows into the CGM, and the re-arrangement of baryons by SMBH feedback. A description of the flows of mass, metals, and energy will only be complete by observing the thermodynamic states, chemical compositions, structure, and dynamics of T>=10^6 K halos. These measurements are uniquely possible with a next-generation X-ray observatory if it provides the sensitivity to detect faint CGM emission, spectroscopic power to measure absorption lines and gas motions, and high spatial resolution to resolve structures.

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.

We present results from the Keck Baryonic Structure Survey (KBSS) including the first detailed measurements of the column densities, kinematics, and internal energy of metal-bearing gas within the virial radius (35-100 physical kpc) of eight ~$L^*$ galaxies at $z\sim2$. From our full sample of 130 metal-bearing absorbers, we infer that halo gas is kinematically complex when viewed in singly, doubly, and triply ionized species. Broad OVI and CIV absorbers are detected at similar velocities to the lower-ionization gas but with very different kinematic structure indicating that the circumgalactic medium (CGM) is multi-phase. There is a high covering fraction of metal-bearing gas within 100 kpc including highly ionized gas such as OVI; however, observations of a single galaxy probed by a lensed background QSO suggest the size of metal-bearing clouds is small (<400 pc for all but the OVI-bearing gas). The mass in metals found within the halo is substantial, equivalent to $\gtrsim$ 25% of the metal mass within the interstellar medium. The gas kinematics unambiguously show that 70% of galaxies with detected metal absorption have some unbound metal-enriched gas, suggesting galactic winds may commonly eject gas from halos at $z\sim2$. Significant thermal broadening is detected in CGM absorbers which dominates the internal energy of the gas. 40% of the detected gas has temperatures in the range $10^{4.5-5.5}$ K where cooling times are short, suggesting the CGM is dynamic, with constant heating and/or cooling to produce this short-lived thermal phase.

Traditional cosmological hydrodynamics simulations fail to spatially resolve the circumgalatic medium (CGM), the reservoir of tenuous gas surrounding a galaxy and extending to its virial radius. We introduce the technique of Enhanced Halo Resolution (EHR), enabling more realistic physical modeling of the simulated CGM by consistently forcing gas refinement to smaller scales throughout the virial halo of a simulated galaxy. We investigate the effects of EHR in the Tempest simulations, a suite of Enzo-based cosmological zoom simulations following the evolution of an L* galaxy, resolving spatial scales of 500 comoving pc out to 100 comoving kpc in galactocentric radius. Among its many effects, EHR (1) changes the thermal balance of the CGM, increasing its cool gas content and decreasing its warm/hot gas content; (2) preserves cool gas structures for longer periods; and (3) enables these cool clouds to exist at progressively smaller size scales. Observationally, this results in a boost in "low ions" like H I and a drop in "high ions" like O VI throughout the CGM. These effects of EHR do not converge in the Tempest simulations, but extrapolating these trends suggests that the CGM in reality is a mist consisting of ubiquitous, small, long-lived, cool clouds suspended in a hot medium at the virial temperature of the halo. Additionally, we explore the physical mechanisms to explain why EHR produces the above effects, proposing that it works both by (1) better sampling the distribution of CGM phases enabling runaway cooling in the denser, cooler tail of the phase distribution; and (2) preventing cool gas clouds from artificially mixing with the ambient hot halo and evaporating. Evidence is found for both EHR mechanisms occurring in the Tempest simulations.

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.

Cosmological simulations predict that a significant fraction of the low-$z$ baryon budget resides in large-scale filaments in the form of a diffuse plasma at temperatures $T \sim 10^{5} - 10^{7}$ K. However, direct observation of this so-called warm-hot intergalactic medium (WHIM) has been elusive. In the $\Lambda$CDM paradigm, galaxy clusters correspond to the nodes of the cosmic web at the intersection of several large-scale filamentary threads. In previous work, we used HST/COS data to conduct the first survey of broad HI Ly$\alpha$ absorbers (BLAs) potentially produced by WHIM in inter-cluster filaments. We targeted a single QSO, namely Q1410, whose sight-line intersects $7$ independent inter-cluster axes at impact parameters $<3$ Mpc (co-moving), and found a tentative excess of a factor of ${\sim}4$ with respect to the field. Here, we further investigate the origin of these BLAs by performing a blind galaxy survey within the Q1410 field using VLT/MUSE. We identified $77$ sources and obtained the redshifts for $52$ of them. Out of the total sample of $7$ BLAs in inter-cluster axes, we found $3$ without any galaxy counterpart to stringent luminosity limits ($\sim 4 \times 10^{8}$ L$_{\odot}$ ${\sim} 0.01$ L$_{*}$), providing further evidence that these BLAs may represent genuine WHIM detections. We combined this sample with other suitable BLAs from the literature and inferred the corresponding baryon mean density for these filaments in the range $\Omega^{\rm fil}_{\rm bar}= 0.02-0.04$. Our rough estimates are consistent with the predictions from numerical simulations but still subject to large systematic uncertainties, mostly from the adopted geometry, ionization corrections and density profile.

Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

We present the extended GALEX Arecibo SDSS Survey (xGASS), a gas fraction-limited census of the atomic (HI) gas content of 1179 galaxies selected only by stellar mass ($M_\star =10^{9}-10^{11.5} M_\odot$) and redshift ($0.01<z<0.05$). This includes new Arecibo observations of 208 galaxies, for which we release catalogs and HI spectra. In addition to extending the GASS HI scaling relations by one decade in stellar mass, we quantify total (atomic+molecular) cold gas fractions and molecular-to-atomic gas mass ratios, $R_{mol}$, for the subset of 477 galaxies observed with the IRAM 30 m telescope. We find that atomic gas fractions keep increasing with decreasing stellar mass, with no sign of a plateau down to $\log M_\star/M_\odot = 9$. Total gas reservoirs remain HI-dominated across our full stellar mass range, hence total gas fraction scaling relations closely resemble atomic ones, but with a scatter that strongly correlates with $R_{mol}$, especially at fixed specific star formation rate. On average, $R_{mol}$ weakly increases with stellar mass and stellar surface density $\mu_\star$, but individual values vary by almost two orders of magnitude at fixed $M_\star$ or $\mu_\star$. We show that, for galaxies on the star-forming sequence, variations of $R_{mol}$ are mostly driven by changes of the HI reservoirs, with a clear dependence on $\mu_\star$. Establishing if galaxy mass or structure plays the most important role in regulating the cold gas content of galaxies requires an accurate separation of bulge and disk components for the study of gas scaling relations.

Hydrodynamical simulations are increasingly able to accurately model physical systems on stellar, galactic, and cosmological scales, however, the utility of these simulations is often limited by our ability to directly compare them with the datasets produced by observers: spectra, photometry, etc. To address this problem, we have created Trident}, a Python-based, open-source tool for post-processing hydrodynamical simulations to produce synthetic absorption spectra and related data. Trident} can (i) create absorption-line spectra for any trajectory through a simulated dataset mimicking both background quasar and down-the-barrel configurations, (ii) reproduce the spectral characteristics of common instruments like the Cosmic Origins Spectrograph, (iii) operate across the ultraviolet, optical and infrared using customizable absorption line lists, (iv) trace simulated physical structures directly to spectral features, (v) approximate the presence of ion species absent from the simulation outputs, (vi) generate column density maps for any ion, and (vii) provide support for all major astrophysical hydrodynamical codes. The focus of Trident's development is for using simulated datasets to better interpret observations of the circumgalactic medium (CGM) and intergalactic medium (IGM), but it remains a general tool applicable in other contexts.

We present the Grackle chemistry and cooling library for astrophysical simulations and models. Grackle provides a treatment of non-equilibrium primordial chemistry and cooling for H, D, and He species, including H2 formation on dust grains; tabulated primordial and metal cooling; multiple UV background models; and support for radiation transfer and arbitrary heat sources. The library has an easily implementable interface for simulation codes written in C, C++, and Fortran as well as a Python interface with added convenience functions for semi-analytical models. As an open-source project, Grackle provides a community resource for accessing and disseminating astrochemical data and numerical methods. We present the full details of the core functionality, the simulation and Python interfaces, testing infrastructure, performance, and range of applicability. Grackle is a fully open-source project and new contributions are welcome.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

We describe a new method for adding a prescribed amount of kinetic energy to simulated gas modeled on a cartesian grid by directly altering grid cells' mass and velocity in a distributed fashion. The method is explored in the context of supernova feedback in high-resolution ($\sim 10$ pc) hydrodynamic simulations of galaxy formation. Resolution-dependence is a primary consideration in our application of the method and simulations of isolated explosions (performed at different resolutions) motivate a resolution-dependent scaling for the injected fraction of kinetic energy that we apply in cosmological simulations of a $10^9$ Msun dwarf halo. We find that in high density media ($\gtrsim$ 50 cm$^{-3}$) with coarse resolution ($\gtrsim 4$ pc per cell), results are sensitive to the initial kinetic energy fraction due to early and rapid cooling. In our galaxy simulations, the deposition of small amounts of supernova energy in kinetic form (as little as 1%) has a dramatic impact on the evolution of the system, resulting in an order of magnitude suppression of stellar mass. The overall behavior of the galaxy in the two highest resolution simulations we perform appears to converge. We discuss the resulting distribution of stellar metallicities, an observable sensitive to galactic wind properties, and find that while the new method demonstrates increased agreement with observed systems, significant discrepancies remain, likely due to simplistic assumptions that neglect contributions from Type Ia supernovae and stellar winds.

We investigate the dynamical impact of cosmic rays in cosmological simulations of galaxy formation using adaptive-mesh refinement simulations of a $10^{12}$ solar mass halo. In agreement with previous work, a run with only our standard thermal energy feedback model results in a massive spheroid and unrealistically peaked rotation curves. However, the addition of a simple two-fluid model for cosmic rays drastically changes the morphology of the forming disk. We include an isotropic diffusive term and a source term tied to star formation due to (unresolved) supernova-driven shocks. Over a wide range of diffusion coefficients, the CRs generate thin, extended disks with a significantly more realistic (although still not flat) rotation curve. We find that the diffusion of CRs is key to this process, as they escape dense star forming clumps and drive outflows within the more diffuse ISM.

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

We present the final data release from the GALEX Arecibo SDSS Survey (GASS), a large Arecibo program that measured the HI properties for an unbiased sample of ~800 galaxies with stellar masses greater than 10^10 Msun and redshifts 0.025<z<0.05. This release includes new Arecibo observations for 250 galaxies. We use the full GASS sample to investigate environmental effects on the cold gas content of massive galaxies at fixed stellar mass. The environment is characterized in terms of dark matter halo mass, obtained by cross-matching our sample with the SDSS group catalog of Yang et al. Our analysis provides, for the first time, clear statistical evidence that massive galaxies located in halos with masses of 10^13-10^14 Msun have at least 0.4 dex less HI than objects in lower density environments. The process responsible for the suppression of gas in group galaxies most likely drives the observed quenching of the star formation in these systems. Our findings strongly support the importance of the group environment for galaxy evolution, and have profound implications for semi-analytic models of galaxy formation, which currently do not allow for stripping of the cold interstellar medium in galaxy groups.

This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.

Cosmological hydrodynamical simulations of galaxy evolution are increasingly able to produce realistic galaxies, but the largest hurdle remaining is in constructing subgrid models that accurately describe the behavior of stellar feedback. As an alternate way to test and calibrate such models, we propose to focus on the circumgalactic medium. To do so, we generate a suite of adaptive-mesh refinement (AMR) simulations for a Milky-Way-massed galaxy run to z=0, systematically varying the feedback implementation. We then post-process the simulation data to compute the absorbing column density for a wide range of common atomic absorbers throughout the galactic halo, including H I, Mg II, Si II, Si III, Si IV, C IV, N V, O VI, and O VII. The radial profiles of these atomic column densities are compared against several quasar absorption line studies, to determine if one feedback prescription is favored. We find that although our models match some of the observations (specifically those ions with lower ionization strengths), it is particularly difficult to match O VI observations. There is some indication that the models with increased feedback intensity are better matches. We demonstrate that sufficient metals exist in these halos to reproduce the observed column density distribution in principle, but the simulated circumgalactic medium lacks significant multiphase substructure and is generally too hot. Furthermore, we demonstrate the failings of inflow-only models (without energetic feedback) at populating the CGM with adequate metals to match observations even in the presence of multiphase structure. Additionally, we briefly investigate the evolution of the CGM from z=3 to present. Overall, we find that quasar absorption line observations of the gas around galaxies provide a new and important constraint on feedback models.

We present the second data release from the GALEX Arecibo SDSS Survey (GASS), an ongoing large Arecibo program to measure the HI properties for an unbiased sample of ~1000 galaxies with stellar masses greater than 10^10 Msun and redshifts 0.025<z<0.05. GASS targets are selected from the Sloan Digital Sky Survey (SDSS) spectroscopic and Galaxy Evolution Explorer (GALEX) imaging surveys, and are observed until detected or until a gas mass fraction limit of a few per cent is reached. This second data installment includes new Arecibo observations of 240 galaxies, and marks the 50% of the complete survey. We present catalogs of the HI, optical and ultraviolet parameters for these galaxies, and their HI-line profiles. Having more than doubled the size of the sample since the first data release, we also revisit the main scaling relations of the HI mass fraction with galaxy stellar mass, stellar mass surface density, concentration index, and NUV-r color, as well as the gas fraction plane introduced in our earlier work.

We have obtained long-slit spectra of 174 star-forming galaxies with stellar masses greater than 10^10 M_\odot from the GALEX Arecibo SDSS (GASS) survey. These galaxies have both HI and H_2 mass measurements. The average metallicity profile is strikingly flat out to R_90, the radius enclosing 90% of the r-band light. Metallicity profiles which decline steadily with radius are found primarily for galaxies in our sample with low stellar mass (Log(M_*)<10.2), concentration, and/or mean stellar mass density. Beyond ~R_90, however, around 10 percent of the galaxies in our sample exhibit a sharp downturn in metallicity. Remarkably, we find that the magnitude of the outer metallicity drop is well correlated with the total HI content of the galaxy (measured as f_HI=M_HI/M_*). We examine the radial profiles of stellar population ages and star formation rate densities, and conclude that the galaxies with largest outer metallicity drops are actively growing their stellar disks, with mass doubling times across the whole disk only one third as long as a typical GASS galaxy. We also describe a correlation between local stellar mass density and metallicity, which is valid across all galaxies in our sample. We argue that much of the recent stellar mass growth at the edges of these galaxies can be linked to the accretion or radial transport of relatively pristine gas from beyond the galaxies' stellar disks.

We carry out adaptive mesh refinement (AMR) cosmological simulations of Milky-Way mass halos in order to investigate the formation of disk-like galaxies in a {\Lambda}-dominated Cold Dark Matter model. We evolve a suite of five halos to z = 0 and find gaseous-disk formation in all; however, in agreement with previous SPH simulations (that did not include a subgrid feedback model), the rotation curves of all halos are centrally peaked due to a massive spheroidal component. Our standard model includes radiative cooling and star formation, but no feedback. We further investigate this angular momentum problem by systematically modifying various simulation parameters including: (i) spatial resolution, ranging from 1700 to 212 pc; (ii) an additional pressure component to ensure that the Jeans length is always resolved; (iii) low star formation efficiency, going down to 0.1%; (iv) fixed physical resolution as opposed to comoving resolution; (v) a supernova feedback model which injects thermal energy to the local cell; and (vi) a subgrid feedback model which suppresses cooling in the immediate vicinity of a star formation event. Of all of these, we find that only the last (cooling suppression) has any impact on the massive spheroidal component. In particular, a simulation with cooling suppression and feedback results in a rotation curve that, while still peaked, is considerably reduced from our standard runs.

We present dynamical scaling relations for a homogeneous and representative sample of ~500 massive galaxies, selected only by stellar mass (>10^10 Msun) and redshift (0.025<z<0.05) as part of the ongoing GALEX Arecibo SDSS Survey. We compare baryonic Tully-Fisher (BTF) and Faber-Jackson (BFJ) relations for this sample, and investigate how galaxies scatter around the best fits obtained for pruned subsets of disk-dominated and bulge-dominated systems. The BFJ relation is significantly less scattered than the BTF when the relations are applied to their maximum samples, and is not affected by the inclination problems that plague the BTF. Disk-dominated, gas-rich galaxies systematically deviate from the BFJ relation defined by the spheroids. We demonstrate that by applying a simple correction to the stellar velocity dispersions that depends only on the concentration index of the galaxy, we are able to bring disks and spheroids onto the same dynamical relation -- in other words, we obtain a generalized BFJ relation that holds for all the galaxies in our sample, regardless of morphology, inclination or gas content, and has a scatter smaller than 0.1 dex. We find that disks and spheroids are offset in the stellar dispersion-size relation, and that the offset is removed when corrected dispersions are used instead. The generalized BFJ relation represents a fundamental correlation between the global dark matter and baryonic content of galaxies, which is obeyed by all (massive) systems regardless of morphology. [abridged]

We use measurements of the HI content, stellar mass and star formation rates in ~190 massive galaxies with stellar masses greater than 10^10 Msun, obtained from the Galex Arecibo SDSS Survey (GASS) described in Paper I (Catinella et al. 2010) to explore the global scaling relations associated with the bin-averaged ratio of the star formation rate over the HI mass, which we call the HI-based star formation efficiency (SFE). Unlike the mean specific star formation rate, which decreases with stellar mass and stellar mass surface density, the star formation efficiency remains relatively constant across the sample with a value close to SFE = 10^-9.5 yr^-1 (or an equivalent gas consumption timescale of ~3 Gyr). Specifically, we find little variation in SFE with stellar mass, stellar mass surface density, NUV-r color and concentration. We interpret these results as an indication that external processes or feedback mechanisms that control the gas supply are important for regulating star formation in massive galaxies. An investigation into the detailed distribution of SFEs reveals that approximately 5% of the sample shows high efficiencies with SFE > 10^-9 yr^-1, and we suggest that this is very likely due to a deficiency of cold gas rather than an excess star formation rate. Conversely, we also find a similar fraction of galaxies that appear to be gas-rich for their given specific star-formation rate, although these galaxies show both a higher than average gas fraction and lower than average specific star formation rate. Both of these populations are plausible candidates for "transition" galaxies, showing potential for a change (either decrease or increase) in their specific star formation rate in the near future. We also find that 36+/-5% of the total HI mass density and 47+/-5% of the total SFR density is found in galaxies with stellar mass greater than 10^10 Msun. [abridged]

We introduce the GALEX Arecibo SDSS Survey (GASS), an on-going large program that is gathering high quality HI-line spectra using the Arecibo radio telescope for an unbiased sample of ~1000 galaxies with stellar masses greater than 10^10 Msun and redshifts 0.025<z<0.05, selected from the SDSS spectroscopic and GALEX imaging surveys. The galaxies are observed until detected or until a low gas mass fraction limit (1.5-5%) is reached. This paper presents the first Data Release, consisting of ~20% of the final GASS sample. We use this data set to explore the main scaling relations of HI gas fraction with galaxy structure and NUV-r colour. A large fraction (~60%) of the galaxies in our sample are detected in HI. We find that the atomic gas fraction decreases strongly with stellar mass, stellar surface mass density and NUV-r colour, but is only weakly correlated with galaxy bulge-to-disk ratio (as measured by the concentration index of the r-band light). We also find that the fraction of galaxies with significant (more than a few percent) HI decreases sharply above a characteristic stellar surface mass density of 10^8.5 Msun kpc^-2. The fraction of gas-rich galaxies decreases much more smoothly with stellar mass. One of the key goals of GASS is to identify and quantify the incidence of galaxies that are transitioning between the blue, star-forming cloud and the red sequence of passively-evolving galaxies. Likely transition candidates can be identified as outliers from the mean scaling relations between gas fraction and other galaxy properties. [abridged]

We follow Paper I with predictions of how gas leaking through the lunar surface could influence the regolith, as might be observed via optical Transient Lunar Phenomena (TLPs) and related effects. We touch on several processes, but concentrate on low and high flow rate extremes, perhaps the most likely. We model explosive outgassing for the smallest gas overpressure at the regolith base that releases the regolith plug above it. This disturbance's timescale and affected area are consistent with observed TLPs; we also discuss other effects. For slow flow, escape through the regolith is prolonged by low diffusivity. Water, found recently in deep magma samples, is unique among candidate volatiles, capable of freezing between the regolith base and surface, especially near the lunar poles. For major outgassing sites, we consider the possible accumulation of water ice. Over geological time ice accumulation can evolve downward through the regolith. Depending on gases additional to water, regolith diffusivity might be suppressed chemically, blocking seepage and forcing the ice zone to expand to larger areas, up to square km scales, again, particularly at high latitudes. We propose an empirical path forward, wherein current and forthcoming technologies provide controlled, sensitive probes of outgassing. The optical transient/outgassing connection, addressed via Earth-based remote sensing, suggests imaging and/or spectroscopy, but aspects of lunar outgassing might be more covert, as indicated above. TLPs betray some outgassing, but does outgassing necessarily produces TLPs? We also suggest more intrusive techniques from radar to in-situ probes. Many of these approaches should be practiced in a pristine lunar atmosphere, before significant confusing signals likely to be produced upon humans returning to the Moon.

We consider the implications from Paper I on how gas leaking through the lunar surface might interact with the regolith, and in what respects this might affect or cause the appearance of optical Transient Lunar Phenomena (TLPs). We consider briefly a range of phenomena, but concentrate at the extremes of high and low gas flow rate, which might represent the more likely behaviors. Extremely fast i.e., explosive, expulsion of gas from the surface is investigated by examining the minimal amount of gas needed to displace a plug of regolith above a site of gaseous overpressure at the regolith's base. The area and timescale of this disturbance, it is consistent with observed TLPs. Furthermore there are several ways in which such an explosion might be expected to change the lunar surface appearance in a way consistent with many TLPs, including production of obscuration, brightening and color changes. At the slow end of the volatile flow range, gas seeping from the interior is retained below the surface for extensive times due to the low diffusivity of regolith material. A special circumstance arises if the volatile flow contains water vapor, because water is uniquely capable of freezing as it passes from the base to the surface of the regolith. For a large TLP site, it is plausible to think of areas on the square-km scale accumulating significant bodies of water ice. Furthermore, as the system evolves over geological time, the ice accumulation zone will evolve downwards into the regolith. Since many reactions possible between the volatiles and regolith, depending on the additional gases besides water, can act to decrease diffusivity in the regolith, it is plausible that the volatiles produce a barrier between the seepage source and vacuum, forcing the ice zone to expand to larger areas.

We present moderate resolution (~6 km/s) spectroscopy of 284 M giant candidates selected from the Two Micron All Sky Survey photometry. Radial velocities (RVs) are presented for stars mainly in the south, with a number having positions consistent with association to the trailing tidal tail of the Sagittarius (Sgr) dwarf galaxy. The latter show a clear RV trend with orbital longitude, as expected from models of the orbit and destruction of Sgr. A minimum 8 kpc width of the trailing stream about the Sgr orbital midplane is implied by verified RV members. The coldness of this stream (dispersion ~10 km/s) provides upper limits on the combined contributions of stream heating by a lumpy Galactic halo and the intrinsic dispersion of released stars, which is a function of the Sgr core mass. The Sgr trailing arm is consistent with a Galactic halo containing one dominant, LMC-like lump, however some lumpier halos are not ruled out. An upper limit to the total M/L of the Sgr core is 21 in solar units. A second structure that roughly mimics expectations for wrapped, leading Sgr arm debris crosses the trailing arm in the Southern Hemisphere; however, this may also be an unrelated tidal feature. Among the <13 kpc M giants toward the South Galactic Pole are some with large RVs that identify them as halo stars, perhaps part of the Sgr leading arm near the Sun. The positions and RVs of Southern Hemisphere M giants are compared with those of southern globular clusters potentially stripped from the Sgr system and support for association of Pal 2 and Pal 12 with Sgr debris is found. Our discussion includes description of a masked-filtered cross-correlation methodology that achieves better than 1/20 of a resolution element RVs in moderate resolution spectra.

We describe first results of a spectroscopic probe of selected fields from the Grid Giant Star Survey. Multifiber spectroscopy of several hundred stars in a strip of eleven fields along delta approximately -17^{circ}, in the range 12 <~ alpha <~ 17 hours, reveals a group of 8 giants that have kinematical characteristics differing from the main field population, but that as a group maintain coherent, smoothly varying distances and radial velocities with position across the fields. Moreover, these stars have roughly the same abundance, according to their MgH+Mgb absorption line strengths. Photometric parallaxes place these stars in a semi-loop structure, arcing in a contiguous distribution between 5.7 and 7.9 kpc from the Galactic center. The spatial, kinematical, and abundance coherence of these stars suggests that they are part of a diffuse stream of tidal debris, and one roughly consistent with a wrapped, leading tidal arm of the Sagittarius dwarf spheroidal galaxy.