We investigate angular momentum acquisition in Milky Way-sized galaxies by comparing five high resolution zoom-in simulations, each implementing identical cosmological initial conditions, but utilizing different hydrodynamic codes: Enzo, Art, Ramses, Arepo, and Gizmo-PSPH. Each code implements a distinct set of feedback and star formation prescriptions. We find that while many galaxy and halo properties vary between the different codes (and feedback prescriptions), there is qualitative agreement on the process of angular momentum acquisition in the galaxy's halo. In all simulations, cold filamentary gas accretion to the halo results in ~4 times more specific angular momentum in cold halo gas ($\lambda_{cold} \gtrsim 0.1$) than in the dark matter halo. At z>1, this inflow takes the form of inspiraling cold streams that are co-directional in the halo of the galaxy and are fueled, aligned, and kinematically connected to filamentary gas infall along the cosmic web. Due to the qualitative agreement among disparate simulations, we conclude that the buildup of high angular momentum halo gas and the presence of these inspiraling cold streams are robust predictions of Lambda Cold Dark Matter galaxy formation, though the detailed morphology of these streams is significantly less certain. A growing body of observational evidence suggests that this process is borne out in the real universe.

We assess the fraction of the Milky Way's circumgalactic medium (CGM) eluding detection due to its velocity being similar to gas in the disk. This is achieved using synthetic observations of the CGM in a simulated MW-mass galaxy that shows similar CGM kinematics to the MW and external L$\sim$L$_*$ galaxies. As viewed by a mock observer at a location similar to the Sun, only 50$\%$ (by mass) of the gas moves at high velocity ($|v_{\rm LSR}|\geq$100 km s$^{-1}$ or $|v_{\rm DEV}|\geq$50 km s$^{-1}$) in the simulated CGM and would be observable. The low velocity gas is thermodynamically similar to the high velocity gas, indicating the 50$\%$ observable fraction is applicable to spectral lines from the radio to the ultraviolet. We apply the observable mass fraction (50$\%$) to current estimates of the MW's CGM, and find a corrected total mass of 2.8$\times$10$^{8} M_{\odot}$ for gas below 10$^6$K within $\sim15$ kpc (this excludes the Magellanic System). This is less than the total mass of the CGM extending out to $\sim$150 kpc in other L$\sim$L$_*$ galaxies. However, we find similar OVI column densities when the discrepancy in path length between the MW and external galaxies is considered. The coherent spatial and kinematic distribution of low velocity gas in the simulated CGM suggests that current HI observations of the MW's CGM may miss large low velocity HI complexes. In addition, current mass estimates of the MW's CGM based on high-velocity line observations with distance constraints may miss a non-negligible fraction of gas in the outer halo which can be obscured if it moves at a velocity similar to the gas in the lower halo.

Galactic halo gas traces inflowing star formation fuel and feedback from a galaxy's disk and is therefore crucial to our understanding of galaxy evolution. In this review, we summarize the multi-wavelength observational properties and origin models of Galactic and low redshift spiral galaxy halo gas. Galactic halos contain multiphase gas flows that are dominated in mass by the ionized component and extend to large radii. The densest, coldest halo gas observed in neutral hydrogen (HI) is generally closest to the disk (< 20 kpc), and absorption line results indicate warm and warm-hot diffuse halo gas is present throughout a galaxy's halo. The hot halo gas detected is not a significant fraction of a galaxy's baryons. The disk-halo interface is where the multiphase flows are integrated into the star forming disk, and there is evidence for both feedback and fueling at this interface from the temperature and kinematic gradient of the gas and HI structures. The origin and fate of halo gas is considered in the context of cosmological and idealized local simulations. Accretion along cosmic filaments occurs in both a hot (> 10^5.5 K) and cold mode in simulations, with the compressed material close to the disk the coldest and densest, in agreement with observations. There is evidence in halo gas observations for radiative and mechanical feedback mechanisms, including escaping photons from the disk, supernova-driven winds, and a galactic fountain. Satellite accretion also leaves behind abundant halo gas. This satellite gas interacts with the existing halo medium, and much of this gas will become part of the diffuse halo before it can reach the disk. The accretion rate from cold and warm halo gas is generally below a galaxy disk's star formation rate, but gas at the disk-halo interface and stellar feedback may be important additional fuel sources.

We examine galaxy formation in a cosmological AMR simulation, which includes two high resolution boxes, one centered on a 3 \times 10^14 M\odot cluster, and one centered on a void. We examine the evolution of 611 massive (M\ast > 10^10M\odot) galaxies. We find that the fraction of the final stellar mass which is accreted from other galaxies is between 15 and 40% and increases with stellar mass. The accreted fraction does not depend strongly on environment at a given stellar mass, but the galaxies in groups and cluster environments are older and underwent mergers earlier than galaxies in lower density environments. On average, the accreted stars are ~2.5 Gyrs older, and ~0.15 dex more metal poor than the stars formed in-situ. Accreted stellar material typically lies on the outskirts of galaxies; the average half-light radius of the accreted stars is 2.6 times larger than that of the in-situ stars. This leads to radial gradients in age and metallicity for massive galaxies, in qualitative agreement with observations. Massive galaxies grow by mergers at a rate of approximately 2.6% per Gyr. These mergers have a median (mass-weighted) mass ratio less than 0.26 \pm 0.21, with an absolute lower limit of 0.20, for galaxies with M\ast ~ 10^12 M\odot. This suggests that major mergers do not dominate in the accretion history of massive galaxies. All of these results agree qualitatively with results from SPH simulations by Oser et al. (2010, 2012).

We perform high-resolution hydrodynamic simulations of a Milky Way-mass galaxy in a fully cosmological setting using the adaptive mesh refinement code, Enzo, and study the kinematics of gas in the simulated galactic halo. We find that the gas inflow occurs mostly along filamentary structures in the halo. The warm-hot (10^5 K < T < 10^6 K) and hot (T > 10^6 K) ionized gases are found to dominate the overall mass accretion in the system (with dM/dt = 3-5 M_solar/yr) over a large range of distances, extending from the virial radius to the vicinity of the disk. Most of the inflowing gas (by mass) does not cool, and the small fraction that manages to cool does so primarily close to the galaxy (R <~ 20 kpc), perhaps comprising the neutral gas that may be detectable as, e.g., high-velocity clouds. The neutral clouds are embedded within larger, accreting filamentary flows, and represent only a small fraction of the total mass inflow rate. The inflowing gas has relatively low metallicity (Z/Z_solar < 0.2). The outer layers of the filamentary inflows are heated due to compression as they approach the disk. In addition to the inflow, we find high-velocity, metal-enriched outflows of hot gas driven by supernova feedback. Our results are consistent with observations of halo gas at low z.

We analyze an adaptive mesh refinement hydrodynamic cosmological simulation of a Milky Way-sized galaxy to study the cold gas in the halo. HI observations of the Milky Way and other nearby spirals have revealed the presence of such gas in the form of clouds and other extended structures, which indicates on-going accretion. We use a high-resolution simulation (136-272 pc throughout) to study the distribution of cold gas in the halo, compare it with observations, and examine its origin. The amount (10^8 Msun in HI), covering fraction, and spatial distribution of the cold halo gas around the simulated galaxy at z=0 are consistent with existing observations. At z=0 the HI mass accretion rate onto the disk is 0.2 Msun/yr. We track the histories of the 20 satellites that are detected in HI in the redshift interval 0.5>z>0 and find that most of them are losing gas, with a median mass loss rate per satellite of 3.1 x 10^{-3} Msun/yr. This stripped gas is a significant component of the HI gas seen in the simulation. In addition, we see filamentary material coming into the halo from the IGM at all redshifts. Most of this gas does not make it directly to the disk, but part of the gas in these structures is able to cool and form clouds. The metallicity of the gas allows us to distinguish between filamentary flows and satellite gas. We find that the former accounts for at least 25-75% of the cold gas in the halo seen at any redshift analyzed here. Placing constraints on cloud formation mechanisms allows us to better understand how galaxies accrete gas and fuel star formation at z=0.

Stellar feedback drives the circulation of matter from the disk to the halo of galaxies. We perform three-dimensional magnetohydrodynamic simulations of a vertical column of the interstellar medium with initial conditions typical of the solar circle in which supernovae drive turbulence and determine the vertical stratification of the medium. The simulations were run using a stable, positivity-preserving scheme for ideal MHD implemented in the FLASH code. We find that the majority (\approx 90 %) of the mass is contained in thermally-stable temperature regimes of cold molecular and atomic gas at T < 200 K or warm atomic and ionized gas at 5000 K < T < 10^{4.2} K, with strong peaks in probability distribution functions of temperature in both the cold and warm regimes. The 200 - 10^{4.2} K gas fills 50-60 % of the volume near the plane, with hotter gas associated with supernova remnants (30-40 %) and cold clouds (< 10 %) embedded within. At |z| ~ 1-2 kpc, transition-temperature (10^5 K) gas accounts for most of the mass and volume, while hot gas dominates at |z| > 3 kpc. The magnetic field in our models has no significant impact on the scale heights of gas in each temperature regime; the magnetic tension force is approximately equal to and opposite the magnetic pressure, so the addition of the field does not significantly affect the vertical support of the gas. The addition of a magnetic field does reduce the fraction of gas in the cold (< 200 K) regime with a corresponding increase in the fraction of warm (~ 10^4 K) gas. However, our models lack rotational shear and thus have no large-scale dynamo, which reduces the role of the field in the models compared to reality. The supernovae drive oscillations in the vertical distribution of halo gas, with the period of the oscillations ranging from ~ 30 Myr in the T < 200 K gas to ~ 100 Myr in the 10^6 K gas, in line with predictions by Walters & Cox.

We present new 3D magnetohydrodynamic (MHD) simulations of a supernova-driven, stratified interstellar medium. These simulations were run using the Waagan (2009) positivity preserving scheme for ideal MHD implemented in the Flash code. The scheme is stable even for the Mach numbers approaching 100 found in this problem. We have previously shown that the density distribution arising from hydrodynamical versions of these simulations creates low-density pathways through which Lyman continuum photons can travel to heights |z| > 1 kpc. This naturally produces the warm ionized medium through photoionization due primarily to O stars near the plane. However, our earlier models reproduce the peak but not the width of the observed emission measure distribution. Here, we examine whether inclusion of magnetic fields and a greater vertical extent to the simulation domain produce a gas distribution that better matches the observations. We further study the change of magnetic energy over time in our models, showing that it appears to reach a steady state after a few hundred megayears, presumably supported by a turbulent dynamo driven by the supernova explosions.

Using adaptive mesh refinement (AMR) hydrodynamic simulations of vertically stratified hot halo gas, we examine the conditions under which clouds can form and condense out of the hot halo medium to potentially fuel star formation in the gaseous disk. We find that halo clouds do not develop from linear isobaric perturbations. This is a regime where the cooling time is longer than the Brunt-Vaisala time, confirming previous linear analysis. We extend the analysis into the nonlinear regime by considering mildly or strongly nonlinear perturbations with overdensities up to 100, also varying the initial height, the cloud size, and the metallicity of the gas. Here, the result depends on the ratio of cooling time to the time required to accelerate the cloud to the sound speed (similar to the dynamical time). If the ratio exceeds a critical value near unity, the cloud is accelerated without further cooling and gets disrupted by Kelvin-Helmholtz and/or Rayleigh-Taylor instabilities. If it is less than the critical value, the cloud cools and condenses before disruption. Accreting gas with overdensities of 10-20 is expected to be marginally unstable; the cooling fraction will depend on the metallicity, the size of the incoming cloud, and the distance to the galaxy. Locally enhanced overdensities within cold streams have a higher likelihood of cooling out. Our results have implications on the evolution of clouds seeded by cold accretion that are barely resolved in current cosmological hydrodynamic simulations and absorption line systems detected in galaxy halos.

We analyze a (120 h^{-1} Mpc)^3 adaptive mesh refinement hydrodynamic simulation that contains a higher-resolution 31 x 31 x 35 h^{-3} Mpc subvolume centered on a ~30 Mpc diameter void. Our detailed ~1 kpc resolution allows us to identify 1300 galaxies within this void to a limiting halo mass of ~10^{10} M_sun. Nearly 1000 galaxies are found to be in underdense regions, with 300 galaxies residing in regions less than half the mean density of the simulation volume. We construct mock observations of the stellar and gas properties of these systems, and reproduce the range of colors and luminosities observed in the SDSS for nearby (z < 0.03) galaxies. We find no trends with density for the most luminous (M_r < -18) galaxies, however our dwarf void galaxies (M_r > -16), though they are less reliably resolved, typically appear bluer, with higher rates of star formation and specific star formation and lower mean stellar ages than galaxies in average density environments. We find a significant population of low luminosity (M_r ~ -14) dwarf galaxies that is preferentially located in low density regions and specifically in the void center. This population may help to reduce, but not remove, the discrepancy between the predicted and observed number of void galaxies.

We compare the predictions of three physical models for the origin of the hot halo gas with the observed halo X-ray emission, derived from 26 high-latitude XMM-Newton observations of the soft X-ray background between $l=120\degr$ and $l=240\degr$. These observations were chosen from a much larger set of observations as they are expected to be the least contaminated by solar wind charge exchange emission. We characterize the halo emission in the XMM-Newton band with a single-temperature plasma model. We find that the observed halo temperature is fairly constant across the sky (~1.8e6-2.3e6 K), whereas the halo emission measure varies by an order of magnitude (~0.0005-0.006 cm^-6 pc). When we compare our observations with the model predictions, we find that most of the hot gas observed with XMM-Newton does not reside in isolated extraplanar supernova remnants -- this model predicts emission an order of magnitude too faint. A model of a supernova-driven interstellar medium, including the flow of hot gas from the disk into the halo in a galactic fountain, gives good agreement with the observed 0.4-2.0 keV surface brightness. This model overpredicts the halo X-ray temperature by a factor of ~2, but there are a several possible explanations for this discrepancy. We therefore conclude that a major (possibly dominant) contributor to the halo X-ray emission observed with XMM-Newton is a fountain of hot gas driven into the halo by disk supernovae. However, we cannot rule out the possibility that the extended hot halo of accreted material predicted by disk galaxy formation models also contributes to the emission.

We investigate models for the photoionization of the widespread diffuse ionized gas in galaxies. In particular we address the long standing question of the penetration of Lyman continuum photons from sources close to the galactic midplane to large heights in the galactic halo. We find that recent hydrodynamical simulations of a supernova-driven interstellar medium have low density paths and voids that allow for ionizing photons from midplane OB stars to reach and ionize gas many kiloparsecs above the midplane. We find ionizing fluxes throughout our simulation grids are larger than predicted by one dimensional slab models, thus allowing for photoionization by O stars of low altitude neutral clouds in the Galaxy that are also detected in Halpha. In previous studies of such clouds the photoionization scenario had been rejected and the Halpha had been attributed to enhanced cosmic ray ionization or scattered light from midplane H II regions. We do find that the emission measure distributions in our simulations are wider than those derived from Halpha observations in the Milky Way. In addition, the horizontally averaged height dependence of the gas density in the hydrodynamical models is lower than inferred in the Galaxy. These discrepancies are likely due to the absence of magnetic fields in the hydrodynamic simulations and we discuss how magnetohydrodynamic effects may reconcile models and observations. Nevertheless, we anticipate that the inclusion of magnetic fields in the dynamical simulations will not alter our primary finding that midplane OB stars are capable of producing high altitude diffuse ionized gas in a realistic three-dimensional interstellar medium.

Feedback from massive stars is one of the least understood aspects of galaxy formation. We perform a suite of vertically stratified local interstellar medium (ISM) simulations in which supernova rates and vertical gas column densities are systematically varied based on the Schmidt-Kennicutt law. Our simulations have a sufficiently high spatial resolution (1.95 pc) to follow the hydrodynamic interactions among multiple supernovae that structure the ISM. At a given supernova rate, we find that the mean mass-weighted sound speed and velocity dispersion decrease as the inverse square root of gas density, indicating that both thermal and turbulent pressures are nearly constant in the midplane, so the effective equation of state is isobaric. In contrast, across our four models having supernova rates that range from one to 512 times the Galactic supernova rate, the mass-weighted velocity dispersion remains in the range 4-6 km/s. Hence, gas averaged over ~100 pc regions follows an effective equation of state that is close to isothermal. Simulated H I emission lines have widths of 10-18 km/s, comparable to observed values. In our highest supernova rate model, superbubble blow-outs occur, and the turbulent pressure on large scales is >~4 times higher than the thermal pressure. We find a tight correlation between the thermal and turbulent pressures averaged over ~100 pc regions in the midplane of each model, as well as across the four ISM models. We construct a subgrid model for turbulent pressure based on analytic arguments and explicitly calibrate it against our stratified ISM simulations. The subgrid model provides a simple yet physically motivated way to include supernova feedback in cosmological simulations.

Stellar feedback in galactic bulges plays an essential role in shaping the evolution of galaxies. To quantify this role and facilitate comparisons with X-ray observations, we conduct 3D hydrodynamical simulations with the adaptive mesh refinement code, FLASH, to investigate the physical properties of hot gas inside a galactic bulge, similar to that of our Galaxy or M31. We assume that the dynamical and thermal properties of the hot gas are dominated by mechanical energy input from SNe, primarily Type Ia, and mass injection from evolved stars as well as iron enrichment from SNe. We study the bulge-wide outflow as well as the SN heating on scales down to ~4 pc. An embedding scheme that is devised to plant individual SNR seeds, allows to examine, for the first time, the effect of sporadic SNe on the density, temperature, and iron ejecta distribution of the hot gas as well as the resultant X-ray morphology and spectrum. We find that the SNe produce a bulge wind with highly filamentary density structures and patchy ejecta. Compared with a 1D spherical wind model, the non-uniformity of simulated gas density, temperature, and metallicity substantially alters the spectral shape and increases the diffuse X-ray luminosity. The differential emission measure as a function of temperature of the simulated gas exhibits a log-normal distribution, with a peak value much lower than that of the corresponding 1D model. The bulk of the X-ray emission comes from the relatively low temperature and low abundance gas shells associated with SN blastwaves. SN ejecta are not well mixed with the ambient medium, at least in the bulge region. These results, at least partly, account for the apparent lack of evidence for iron enrichment in the soft X-ray-emitting gas in galactic bulges and intermediate-mass elliptical galaxies.[...]

Using state-of-the-art adaptive mesh refinement cosmological hydrodynamic simulations with a spatial resolution of proper 0.21kpc/h in refined subregions embedded within a comoving cosmological volume (27.4Mpc/h)^3, we investigate the sizes of galaxies at z=3 in the standard cold dark matter model where reionization is assumed to complete at zri~6. Our simulated galaxies are found to be significantly smaller than the observed ones: while more than one half of the galaxies observed by HST and VLT ranging from rest-frame UV to optical bands with stellar masses larger than 2E10 Msun have half-light radii larger than 2kpc/h, none of the simulated massive galaxies in the same mass range have half-light radii larger than 2kpc/h, after taking into account dust extinction. Corroborative evidence is provided by the rotation curves of the simulated galaxies with total masses of 1E11-1E12Msun, which display values 300-1000km/s at small radii (0.5kpc/h) due to high stellar concentration in the central regions, larger than those of any well observed galaxies. Possible physical mechanisms to resolve this serious problem include: (1) an early reionization at zri>>6 to suppress gas condensation hence star formation, (2) a strong, internal energetic feedback from stars or central black holes to reduce the overall star formation efficiency, or (3) a substantial small-scale cutoff in the matter power spectrum.

To study how supernova feedback structures the turbulent interstellar medium, we construct 3D models of vertically stratified gas stirred by discrete supernova explosions, including vertical gravitational field and parametrized heating and cooling. The models reproduce many observed characteristics of the Galaxy such as global circulation of gas (i.e., galactic fountain) and the existence of cold dense clouds in the galactic disk. Global quantities of the model such as warm and hot gas filling factors in the midplane, mass fraction of thermally unstable gas, and the averaged vertical density profile are compared directly with existing observations, and shown to be broadly consistent. We find that energy injection occurs over a broad range of scales. There is no single effective driving scale, unlike the usual assumption for idealized models of incompressible turbulence. However, >90% of the total kinetic energy is contained in wavelengths shortward of 200 pc. The shape of the kinetic energy spectrum differs substantially from that of the velocity power spectrum, which implies that the velocity structure varies with the gas density. Velocity structure functions demonstrate that the phenomenological theory proposed by Boldyrev is applicable to the medium. We show that it can be misleading to predict physical properties such as the stellar initial mass function based on numerical simulations that do not include self-gravity of the gas. Even if all the gas in turbulently Jeans unstable regions in our simulation is assumed to collapse and form stars in local freefall times, the resulting total collapse rate is significantly lower than the value consistent with the input supernova rate. Supernova-driven turbulence inhibits star formation globally rather than triggering it.

We study magnetic field steepening due to ambipolar diffusion (Brandenburg & Zweibel 1994) in protoplanetary disk environments and draw the following conclusions. Current sheets are generated in magnetically active regions of the disk where the ionization fraction is high enough for the magnetorotational instability to operate. In late stages of solar nebula evolution, the surface density is expected to have lowered and dust grains to have gravitationally settled to the midplane. If the local dust-to-gas mass ratio near the midplane is increased above cosmic abundances by factors > 10^3, current sheets reach high enough temperatures to melt millimeter-sized dust grains, and hence may provide the mechanism to form meteoritic chondrules. In addition, these current sheets possibly explain the near-infrared excesses observed in spectral energy distributions (SEDs) of young stellar objects. Direct imaging of protoplanetary disks via a nulling interferometer or, in the future, a multi-band, adaptive optics coronagraph can test this hypothesis.

We improve upon the radiative, hydrostatic equilibrium models of passive circumstellar disks constructed by Chiang & Goldreich (1997). New features include (1) account for a range of particle sizes, (2) employment of laboratory-based optical constants of representative grain materials, and (3) numerical solution of the equations of radiative and hydrostatic equilibrium within the original 2-layer (disk surface + disk interior) approximation. We explore how the spectral energy distribution (SED) of a face-on disk depends on grain size distributions, disk geometries and surface densities, and stellar photospheric temperatures. Observed SEDs of 3 Herbig Ae and 2 T Tauri stars, including spectra from the Long Wavelength Spectrometer (LWS) aboard the Infrared Space Observatory (ISO), are fitted with our models. Silicate emission bands from optically thin, superheated disk surface layers appear in nearly all systems. Water ice emission bands appear in LWS spectra of 2 of the coolest stars. Infrared excesses in several sources are consistent with vertical settling of photospheric grains. While this work furnishes further evidence that passive reprocessing of starlight by flared disks adequately explains the origin of infrared-to-millimeter wavelength excesses of young stars, we emphasize how the SED alone does not provide sufficient information to constrain particle sizes and disk masses uniquely.