We investigate how each aspect of a multi-channel stellar feedback model drives the chemodynamical evolution of a low-mass, isolated dwarf galaxy using a suite of high-resolution simulations. Our model follows individual star particles sampled randomly from an adopted initial mass function, considering independently feedback from: supernovae; stellar radiation causing photoelectric heating of dust grains, ionization and associated heating, Lyman-Werner (LW) dissociation of H$_2$, and radiation pressure; and winds from massive main sequence (neglecting their energy input) and asymptotic giant branch (AGB) stars. Radiative transfer is done by ray tracing. We consider the effects each of these processes have on regulating the star formation rate, global properties, multi-phase interstellar medium (ISM), and driving of galactic winds. We follow individual metal species from distinct nucleosynthetic enrichment channels (AGB winds, massive star stellar winds, core collapse and Type Ia supernovae) and pay particular attention to how these feedback processes regulate metal mixing in the ISM, the metal content of outflows, and the stellar abundance patterns in our galaxy. We find that---for a low-metallicity, low-mass dwarf galaxy ---stellar radiation, particularly ionizing radiation and LW radiation, are important sources of stellar feedback whose effects dominate over photoelectric heating and HI radiation pressure. However, feedback is coupled non-linearly, and the inclusion or exclusion of each process produces non-negligible effects. We find strong variations with: the star formation history; the ejection fractions of metals, mass, and energy; and the distribution of elements from different nucleosynthetic sources in both the gas and stars.

Physical properties of the interstellar medium (ISM) at sub-galactic ($\sim$kpc) scales play an indispensable role in controlling the ability of gas to form stars. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, in this paper, we use the TNG50 cosmological simulation of the IllustrisTNG project to explore the physical parameter space comprised of 8 resolved ISM properties in star-forming regions across a wide range of host galaxy mass (M$_\star = 10^{7-11}$ M$_\odot$) and redshift ($0 \leq z \leq 3$). We deconstruct our simulated galaxies into regions of 1 kpc in size, and measure the gas/stellar surface densities, gas metallicity, vertical stellar velocity dispersion, epicyclic frequency and dark-matter volumetric density representative of each region, and study them in the context of their star formation activity and galactic environment (radial galactocentric location). By examining the star formation rate weighted distributions of these properties, we show that stars primarily form in two distinct regimes, which are brought about by an underlying bi-component radial star formation rate profile in galaxies. We examine how the relative prominence of these two regimes depends on host galaxy mass and cosmic time, and compare our findings with observations from the MaNGA IFU survey. Further, using principal component analysis, we characterise the aforementioned parameter space to reveal a high-degree of multicollinearity in relationships between ISM properties that drive the distribution of star formation at kiloparsec scales. Based on this, we find that a reduced three dimensional representation underpinned in essence by a multi-variate radius relationship is sufficient to capture most of the variance in the original 8D space.

A cool phase of the interstellar medium has been observed in many giant elliptical galaxies, but its origin remains unclear. We propose that uneven heating from Type Ia supernovae (SNe Ia), together with radiative cooling, can lead to the formation of the cool phase. The basic idea is that since SNe Ia explode randomly, gas parcels which are not directly heated by SN shocks will cool, forming multiphase gas. We run a series of idealized high-resolution numerical simulations, and find that cool gas develops even when the overall SNe heating rate $H$ exceeds the cooling rate $C$ by a factor as large as 1.4. We also find that the time for multiphase gas development depends on the gas temperature. When the medium has a temperature $T = 3\times 10^6$ K, the cool phase forms within one cooling time \tc; however, the cool phase formation is delayed to a few times \tc\ for higher temperatures. The main reason for the delay is turbulent mixing. Cool gas formed this way would naturally have a metallicity lower than that of the hot medium. For constant $H/C$, there is more turbulent mixing for higher temperature gas. We note that this mechanism of producing cool gas cannot be captured in cosmological simulations, which usually fail to resolve individual SN remnants.

Type Ia supernovae (SNe Ia) provide unique and important feedback in quiescent galaxies, but their impact has been underappreciated. In this paper, we analyze a series of high-resolution simulations to examine the energetics and turbulence of the medium under SNe Ia. We find that when SN remnants are resolved, their effects differ distinctly from a volumetric heating term, as is commonly assumed in unresolved simulations. First, the net heating is significantly higher than expected, by 30$\pm$10\% per cooling time. This is because a large fraction of the medium is pushed into lower densities which cool inefficiently. Second, the medium is turbulent; the root-mean-squared (RMS) velocity of the gas to 20-50 km s$^{-1}$ on a driving scale of tens of parsec. The velocity field of the medium is dominated by compressional modes, which are larger than the solenoidal components by a factor of 3-7. Third, the hot gas has a very broad density distribution. The ratio between the density fluctuations and the RMS Mach number, parameterized as $b$, is 2-20. This is in contrast to previous simulations of turbulent media, which have found $b\lesssim$ 1. The reason for the difference is mainly caused by the \textit{localized} heating of SNe Ia, which creates a large density contrast. Last, the typical length scale of a density fluctuation grows with time, forming increasingly larger bubbles and filamentary ridges. These underlying density fluctuations need to be included when X-ray observations are interpreted.

Semi-analytic models (SAMs) are a promising means of tracking the physical processes associated with galaxy formation, but many of their approximations have not been rigorously tested. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, we compare predictions from the FIRE-2 hydrodynamical "zoom-in" simulations to those from the Santa Cruz SAM run on the same halo merger trees, with an emphasis on the global mass flow cycle. Our study includes 13 halos spanning low-mass dwarfs (M_vir~10^10 M_sun at z=0), intermediate-mass dwarfs (M_vir~10^11 M_sun) and Milky Way-mass galaxies (M_vir~10^12 M_sun). The SAM and FIRE-2 predictions agree relatively well with each other in terms of stellar and interstellar mass, but differ dramatically on circumgalactic mass (the SAM is lower than FIRE-2 by ~3 orders of magnitude for dwarfs). Strikingly, the SAM predicts higher gas accretion rates for dwarfs compared to FIRE-2 by factors of ~10-100, and this is compensated for with higher mass outflow rates in the SAM. We argue that the most severe model discrepancies are caused by the lack of preventative stellar feedback and the assumptions for halo gas cooling and recycling in the SAM. As a first step towards resolving these model tensions, we present a simple yet promising new preventative stellar feedback model in which the energy carried by supernova-driven winds is allowed to heat some fraction of gas outside of halos to at least the virial temperature such that accretion is suppressed.

We examine the properties of the circumgalactic medium (CGM) at low redshift in a range of simulated Milky Way mass halos. The sample is comprised of seven idealized simulations, an adaptive mesh refinement cosmological zoom-in simulation, and two groups of 50 halos with star forming or quiescent galaxies taken from the IllustrisTNG100 simulation. The simulations have very different setups, resolution, and feedback models, but are analyzed in a uniform manner. By comparing median radial profiles and mass distributions of CGM properties, we isolate key similarities and differences. In doing so, we advance the efforts of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project that aims to understand the inherently multiscale galaxy formation process. In the cosmological simulations, the CGM exhibits nearly flat temperature distributions, and broad pressure and radial velocity distributions. In the idealized simulations, similar distributions are found in the inner CGM ($\lesssim 0.5 \, r_{\rm 200c}$) when strong galactic feedback models are employed, but the outer CGM ($\gtrsim 0.5 \, r_{\rm 200c}$) has a much less prominent cold phase, and narrower pressure and velocity distributions even in models with strong feedback. This comparative analysis demonstrates the dominant role feedback plays in shaping the inner CGM and the increased importance of cosmological effects, such as nonspherical accretion and satellite galaxies, in the outer CGM. Furthermore, our findings highlight that while cosmological simulations are required to capture the multiphase structure of the CGM at large radii, idealized simulations provide a robust framework to study how galactic feedback interacts with the inner CGM and thereby provide a reliable avenue to constrain feedback prescriptions.

Large scale outflows in star-forming galaxies are observed to be ubiquitous, and are a key aspect of theoretical modeling of galactic evolution in a cosmological context, the focus of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project. Gas blown out from galactic disks, similar to gas within galaxies, consists of multiple phases with large contrasts of density, temperature, and other properties. To study multiphase outflows as emergent phenomena, we run a suite of ~pc-resolution local galactic disk simulations using the TIGRESS framework. Explicit modeling of the interstellar medium (ISM), including star formation and self-consistent radiative heating plus supernova feedback, regulates ISM properties and drives the outflow. We investigate the scaling of outflow mass, momentum, energy, and metal loading factors with galactic disk properties, including star formation rate (SFR) surface density (\Sigma_SFR~10^{-4}-1 M_sun/kpc^2/yr), gas surface density (~1-100 M_sun/pc^2), and total midplane pressure (or weight) (~10^3-10^6 k_B cm^{-3} K). The main components of outflowing gas are mass-delivering cool gas (T~10^4 K) and energy/metal-delivering hot gas (T~10^6 K). Cool mass outflow rates measured at outflow launch points (one or two scale heights) are 1-100 times the SFR (decreasing with \Sigma_SFR), although in massive galaxies most mass falls back due to insufficient outflow velocity. The hot galactic outflow carries mass comparable to 10% of the SFR, together with 10-20% of the energy and 30-60% of the metal mass injected by SN feedback. The characteristic outflow velocities of both phases scale very weakly with SFR, as v_out \propto \Sigma_SFR^{0.1~0.2}, consistent with observations. Importantly, our analysis demonstrates that in any physically-motivated cosmological wind model, it is crucial to include at least two distinct thermal wind components.

We present a new semi-analytic model of the formation of the first stars. Our method takes dark matter halo merger trees (including 3-dimensional spatial information) from cosmological N-body simulations as input and applies analytic prescriptions to compute both the Population III (Pop III) and metal-enriched star formation histories. We have developed a novel method to accurately compute the major feedback processes affecting Pop III star formation: H$_2$ photodissociation from Lyman-Werner (LW) radiation, suppression of star formation due to inhomogeneous reionization, and metal enrichment via supernovae winds. Our method utilizes a grid-based approach relying on fast Fourier transforms (FFTs) to rapidly track the LW intensity, ionization fraction, and metallicity in 3-dimensions throughout the simulation box. We present simulations for a wide range of astrophysical model parameters from $z\approx 30-6$. Initially long-range LW feedback and local metal enrichment and reionization feedback dominate. However, for $z \lesssim 15$ we find that the star formation rate density (SFRD) of Pop III stars is impacted by the combination of external metal enrichment (metals from one halo polluting other pristine halos) and inhomogeneous reionization. We find that the interplay of these processes is particularly important for the Pop III SFRD at $z \lesssim 10$. Reionization feedback delays star formation long enough for metal bubbles to reach halos that would otherwise form Pop III stars. Including these effects can lead to more than an order of magnitude decrease in the Pop III SFRD at $z=6$ compared to LW feedback alone.

Star formation in the universe's most massive galaxies proceeds furiously early in time but then nearly ceases. Plenty of hot gas remains available but does not cool and condense into star-forming clouds. Active galactic nuclei (AGN) release enough energy to inhibit cooling of the hot gas, but energetic arguments alone do not explain why quenching of star formation is most effective in high-mass galaxies. In fact, optical observations show that quenching is more closely related to a galaxy's central stellar velocity dispersion ($\sigma_v$) than to any other characteristic. Here, we show that high $\sigma_v$ is critical to quenching because a deep central potential well maximizes the efficacy of AGN feedback. In order to remain quenched, a galaxy must continually sweep out the gas ejected from its aging stars. Supernova heating can accomplish this task as long as the AGN sufficiently reduces the gas pressure of the surrounding circumgalactic medium (CGM). We find that CGM pressure acts as the control knob on a valve that regulates AGN feedback and suggest that feedback power self-adjusts so that it suffices to lift the CGM out of the galaxy's potential well. Supernova heating then drives a galactic outflow that remains homogeneous if $\sigma_v \gtrsim 240 \, {\rm km \, s^{-1}}$. AGN feedback can effectively quench galaxies with a comparable velocity dispersion, but feedback in galaxies with a much lower velocity dispersion tends to result in convective circulation and accumulation of multiphase gas within the galaxy.

We present an analysis of the angular momentum content of the circumgalactic medium (CGM) using TNG100, one of the flagship runs of the IllustrisTNG project. We focus on Milky Way-mass halos ($\sim 10^{12} \; M_{\odot}$) at $z=0$ but also analyze other masses and redshifts up to $z=5$. We find that the CGM angular momentum properties are strongly correlated with the stellar angular momentum of the corresponding galaxy: the CGM surrounding high-angular momentum galaxies has a systematically higher angular momentum and is better aligned to the rotational axis of the galaxy itself than the CGM surrounding low-angular momentum galaxies. Both the hot and cold phases of the CGM show this dichotomy, though it is stronger for colder gas. The CGM of high-angular momentum galaxies is characterized by a large wedge of cold gas with rotational velocities at least $\sim1/2$ of the halo's virial velocity, extending out to $\sim 1/2$ of the virial radius, and by biconical polar regions dominated by radial velocities suggestive of galactic fountains; both of these features are absent from the CGM of low-angular momentum galaxies. These conclusions are general to halo masses $\lesssim 10^{12} \; M_{\odot}$ and for $z \lesssim 2$, but they do not apply for more massive halos or at the highest redshift studied. By comparing simulations run with alterations to the fiducial feedback model, we identify the better alignment of the CGM to high-angular momentum galaxies as a feedback-independent effect and the galactic winds as a dominant influence on the CGM's angular momentum.

A common situation in galactic and intergalactic gas involves cold dense gas in motion relative to hot diffuse gas. Kelvin-Helmholtz instability creates a turbulent mixing layer and populates the intermediate-temperature phase, which often cools rapidly. The energy lost to cooling is balanced by the advection of hot high enthalpy gas into the mixing layer, resulting in growth and acceleration of the cold phase. This process may play a major role in determining the interstellar medium and circumgalactic medium phase structure, and accelerating cold gas in galactic winds and cosmic filaments. Cooling in these mixing layers occurs in a thin corrugated sheet, which we argue has an area with fractal dimension $D=5/2$ and a thickness that adjusts to match the hot phase mixing time to the cooling time. These cooling sheet properties form the basis of a new model for how the cooling rate and hot gas inflow velocity depend on the size $L$, cooling time $t_{\rm cool}$, relative velocity $v_{\rm rel}$, and density contrast $\rho_{\rm cold}/\rho_{\rm hot}$ of the system. Entrainment is expected to be enhanced in environments with short $t_{\rm cool}$, large $v_{\rm rel}$, and large $\rho_{\rm cold}/\rho_{\rm hot}$. Using a large suite of three dimensional hydrodynamic simulations, we demonstrate that this fractal cooling layer model accurately captures the energetics and evolution of turbulent interfaces and can therefore be used as a foundation for understanding multiphase mixing with strong radiative cooling.

The interaction between a YSO stellar magnetic field and its protostellar disc can result in stellar accretional flows and outflows from the inner disc rim. Gas flows with a velocity component perpendicular to disc midplane subject particles to centrifugal acceleration away from the protostar, resulting in particles being catapulted across the face of the disc. The ejected material can produce a "dust fan", which may be dense enough to mimic the appearance of a "puffed-up" inner disc rim. We derive analytic equations for the time dependent disc toroidal field, the disc magnetic twist, the size of the stable toroidal disc region, the jet speed and the disc region of maximal jet flow speed. We show how the observed infrared variability of the pre-transition disc system LRLL~31 can be modelled by a dust ejecta fan from the inner-most regions of the disc whose height is partially dependent on the jet flow speed. The greater the jet flow speed, the higher is the potential dust fan scale height. An increase in mass accretion onto the star tends to increase the height and optical depth of the dust ejection fan, increasing the amount of 1--8~$\mu$m radiation. The subsequent shadow reduces the amount of light falling on the outer disc and decreases the 8-- 40~$\mu$m radiation. A decrease in the accretion rate reverses this scenario, thereby producing the observed "see-saw" infrared variability.

Supernovae (SNe) drive multiphase galactic outflows, impacting galaxy formation; however, cosmological simulations mostly use \textit{ad hoc} feedback models for outflows, making outflow-related predictions from first principles problematic. Recent small-box simulations resolve individual SNe remnants in the interstellar medium (ISM), naturally driving outflows and permitting a determination of the wind loading factors of energy $\eta_E$, mass $\eta_m$, and metals $\eta_Z$. In this Letter, we compile small-box results, and find consensus that the hot outflows are much more powerful than the cool outflows: (i) hot outflows generally dominate the energy flux, and (ii) their specific energy $e_{s,h}$ is 10-1000 times higher than cool outflows. Moreover, the properties of hot outflows are remarkably simple: $e_{s,h} \propto \eta_{E,h}/\eta_{m,h}$ is almost invariant over four orders of magnitude of star formation surface density. Also, we find tentatively that $\eta_{E,h}/\eta_{Z,h} \sim$ 0.5. If corroborated by more simulation data, these correlations reduce the three hot phase loading factors into one. Finally, this one parameter is closely related to whether the ISM has a "breakout" condition. The narrow range of \esh\ indicates that hot outflows cannot escape dark matter halos with log $M_{\rm{halo}}\ [M_\odot] \gtrsim 12$. This mass is also where the galaxy mass-metallicity relation reaches its plateau, implying a deep connection between \textit{hot} outflows and galaxy formation. We argue that hot outflows should be included explicitly in cosmological simulations and (semi-)analytic modeling of galaxy formation.

The circumgalactic medium (CGM) is closely linked to galaxy formation and evolution, but difficult to characterize observationally and typically poorly resolved in cosmological simulations. We use spherically-symmetric, idealized, high-resolution simulations of the CGM in $10^{12}M_\odot$ and $10^{11}M_\odot$ dark matter halos to characterize the gas pressure, turbulent and radial velocities, and degree of thermal and effective dynamic pressure support in the overall CGM as well as in its high- and low-temperature phases. We find that the $10^{12}M_\odot$ halo contains a CGM mostly formed of a hot gas halo in hydrostatic equilibrium out of which cold gas condenses and falls onto the central galaxy, while the $10^{11}M_\odot$ halo's CGM is not in hydrostatic equilibrium, has a wider spread of properties at a given galactocentric radius, does not have a clear separation of hot and cold phases, and is dominated by bulk motions. We also find that the degree of pressure support in the $10^{11}M_\odot$ halo is strongly dependent on the parameters of the galactic winds of the central galaxy. These results promote the idea that there is no "average" CGM and care must be taken when setting the initial conditions for a small-box simulation of a patch of the CGM.

We study the thermal evolution of UV-irradiated atomic cooling haloes using high-resolution three-dimensional hydrodynamic simulations. We consider the effect of H^- photodetachment by Ly{\alpha} cooling radiation trapped in the optically-thick cores of three such haloes, a process which has not been included in previous simulations. H^- is a precursor of molecular hydrogen, and therefore, its destruction can diminish the H2 abundance and cooling. Using a simple high-end estimate for the trapped Ly{\alpha} energy density, we find that H^- photodetachment by Ly{\alpha} decreases the critical UV flux for suppressing H2-cooling by up to a factor of \approx 5. With a more conservative estimate of the Ly{\alpha} energy density, we find the critical flux is decreased only by ~15-50 percent. Our results suggest that Ly{\alpha} radiation may have an important effect on the thermal evolution of UV-irradiated haloes, and therefore on the potential for massive black hole formation.

One-zone models constructed to match observed stellar abundance patterns have been used extensively to constrain the sites of nucleosynthesis with sophisticated libraries of stellar evolution and stellar yields. The metal mixing included in these models is usually highly simplified, although it is likely to be a significant driver of abundance evolution. In this work we use high-resolution hydrodynamics simulations to investigate how metals from individual enrichment events with varying source energies $E_{\rm ej}$ mix throughout the multi-phase interstellar medium (ISM) of a low-mass ($M_{\rm gas}=2\times 10^{6}$~M$_{\odot}$), low-metallicity, isolated dwarf galaxy. These events correspond to the characteristic energies of both common and exotic astrophysical sites of nucleosynthesis, including: asymptotic giant branch winds ($E_{\rm ej}\sim$10$^{46}$~erg), neutron star-neutron star mergers ($E_{\rm ej}\sim$10$^{49}$~erg), supernovae ($E_{\rm ej}\sim$10$^{51}$~erg), and hypernovae ($E_{\rm ej}\sim$10$^{52}$~erg). We find the mixing timescales for individual enrichment sources in our dwarf galaxy to be long (100~Myr--1~Gyr), with a clear trend of increasing homogeneity for the more energetic events. Given these timescales, we conclude that the spatial distribution and frequency of events are important drivers of abundance homogeneity on large scales; rare, low $E_{\rm ej}$ events should be characterized by particularly broad abundance distributions. The source energy $E_{\rm ej}$ also correlates with the fraction of metals ejected in galactic winds, ranging anywhere from 60\% at the lowest energy to 95\% for hypernovae. We conclude by examining how the radial position, local ISM density, and global star formation rate influence these results.

Supermassive black holes (SMBHs) are thought to provide energy that prevents catastrophic cooling in the centers of massive galaxies and galaxy clusters. However, it remains unclear how this "feedback" process operates. We use high-resolution optical data to study the kinematics of multi-phase filamentary structures by measuring the velocity structure function (VSF) of the filaments over a wide range of scales in the centers of three nearby galaxy clusters: Perseus, Abell 2597 and Virgo. We find that the motions of the filaments are turbulent in all three clusters studied. There is a clear correlation between features of the VSFs and the sizes of bubbles inflated by SMBH driven jets. Our study demonstrates that SMBHs are the main driver of turbulent gas motions in the centers of galaxy clusters and suggests that this turbulence is an important channel for coupling feedback to the environment. Our measured amplitude of turbulence is in good agreement with Hitomi Doppler line broadening measurement and X-ray surface brightness fluctuation analysis, suggesting that the motion of the cold filaments is well-coupled to that of the hot gas. The smallest scales we probe are comparable to the mean free path in the intracluster medium (ICM). Our direct detection of turbulence on these scales provides the clearest evidence to date that isotropic viscosity is suppressed in the weakly-collisional, magnetized intracluster plasma.

A small subset of young stellar objects (YSOs) exhibit "see-saw" temporal variations in their mid-infrared SED; as the flux short-ward of a fulcrum wavelength ($\lambda_{f}$) increases the flux long-wards of this wavelength decreases (and vice-versa) over timescales of weeks to years. While previous studies have shown that an opaque, axisymmetric occulter of variable height can cause this behaviour in the SED of these objects, the conditions under which a single $\lambda_{f}$ occurs have not previously been determined, nor the factors determining its value. Using radiative transfer modelling, we conduct a parametric study of the exemplar of this class, LRLL 31 to explore this phenomenon, and confirm that the cause of this flux variation is likely due to the change in height of the optically thick inner rim of the accretion disc at the dust sublimation radius, or some other phenomenon which results in a similar appearance. We also determine that a fulcrum wavelength only occurs for high inclinations, where the line of sight intersects the accretion disc. Accepting that the disc of LRLL 31 is highly inclined, the inner rim radius, radial and vertical density profiles are independently varied to gauge what effect this had on $\lambda_{f}$ and its position relative to the silicate feature near $10 \mu$m. While $\lambda_{f}$ is a function of each of these parameters, it is found to be most strongly dependent on the vertical density exponent $\beta$. All other factors being held constant, only for flatter discs ($\beta < 1.2$) did we find a $\lambda_{f}$ beyond the silicate feature.

Population III stars forming in minihalos tend to be relatively inefficient, with each minihalo hosting one or a small number of stars which are more massive than local stars but still challenging to observe directly at high redshift. Here we explore a possible mechanisms for the generation of larger clusters of such stars: a nearby ionizing source which ionizes a late forming halo, delaying its collapse until the halo is sufficiently large that the core can self-shield and suffer runaway collapse. We use simulations with a simple but accurate model for the radiative ionizing flux and confirm the basic predictions of previous work: higher ionizing fluxes can delay the collapse to lower redshifts and higher masses, up to an order of magnitude above the atomic cooling limit. In a limited number of runs we also examine the fragmentation of the cores at even higher resolution, using both simple estimates and sink particles to show that the number of fragments is generally small, at most a handful, and that the mass accretion rate on the fragments is of order $10^{-3}$ M$_\odot$ yr$^{-1}$. This rate is sufficiently high that the descent on the main sequence (and hence the suppression of accretion) is delayed until the stellar masses are of order $100-1000$ M$_\odot$, but not high enough to produce direct collapse black holes of mass $\sim 10^5$ M$_{\odot}$. The resulting clusters are larger than those produced in minihalos, but are still likely to fall short of being easily detectable in JWST blind fields.

Cosmic gas cycles in and out of galaxies, but outside of galaxies it is difficult to observe except for the absorption lines that circumgalactic clouds leave in the spectra of background quasars. Using photoionization modeling of those lines to determine cloud pressures, we find that galaxies are surrounded by extended atmospheres that confine the clouds and have a radial pressure profile that depends on galaxy mass. Motivated by observations of the universe's most massive galaxies, we compare those pressure measurements with models predicting the critical pressure at which cooler clouds start to precipitate out of the hot atmosphere and rain toward the center. We find excellent agreement, implying that the precipitation limit applies to galaxies over a wide mass range.

During the last decade, numerous and varied observations, along with increasingly sophisticated numerical simulations, have awakened astronomers to the central role the circumgalactic medium (CGM) plays in regulating galaxy evolution. It contains the majority of the baryonic matter associated with a galaxy, along with most of the metals, and must continually replenish the star forming gas in galaxies that continue to sustain star formation. And while the CGM is complex, containing gas ranging over orders of magnitude in temperature and density, a simple emergent property may be governing its structure and role. Observations increasingly suggest that the ambient CGM pressure cannot exceed the limit at which cold clouds start to condense out and precipitate toward the center of the potential well. If feedback fueled by those clouds then heats the CGM and causes it to expand, the pressure will drop and the "rain" will diminish. Such a feedback loop tends to suspend the CGM at the threshold pressure for precipitation. The coming decade will offer many opportunities to test this potentially fundamental principle of galaxy evolution.

Massive 10^6-10^10 Msun black holes (BHs) are ubiquitous in local galactic nuclei. They were common by the time the Universe is several Gyr old, and many of them were in place within the first 1~Gyr after the Big Bang. Their quick assembly has been attributed to mechanisms such as the rapid collapse of gas into the nuclei of early protogalaxies, accretion and mergers of stellar-mass BHs accompanying structure formation at early times, and the runaway collapse of early, ultra-dense stellar clusters. The origin of the early massive BHs remains an intriguing and long-standing unsolved puzzle in astrophysics. Here we discuss strategies for discerning between BH seeding models using electromagnetic observations. We argue that the most direct answers will be obtained through detection of BHs with masses M<10^5 Msun at redshifts z>10, where we expect them to first form. Reaching out to these redshifts and down to these masses is crucial, because BHs are expected to lose the memory of their initial assembly by the time they grow well above 10^5 Msun and are incorporated into higher-mass galaxies. The best way to detect 10^4-10^5 Msun BHs at high redshifts is by a sensitive X-ray survey. Critical constraining power is augmented by establishing the properties and the environments of their host galaxies in deep optical/IR imaging surveys. Required OIR data can be obtained with the JWST and WFIRST missions. The required X-ray flux limits (down to 10^{-19} erg/s/cm^2) are accessible only with a next-generation X-ray observatory which has both high (sub-1") angular resolution and high throughput. A combination of deep X-ray and OIR surveys will be capable of probing several generic markers of the BH seed scenarios, and resolving the long-stanging puzzle of their origin. These electromagnetic observations are also highly synergistic with the information from LISA on high-z BH mergers.

The past decade has seen an explosion of discoveries and new insights into the diffuse gas within galaxies, galaxy clusters, and the filaments composing the Cosmic Web. A new decade will bring fresh opportunities to further this progress towards developing a comprehensive view of the composition, thermal state, and physical processes of diffuse gas in the Universe. Ultraviolet (UV) spectroscopy, probing diffuse 10^4-10^6 K gas at high spectral resolution, is uniquely poised to (1) witness environmental galaxy quenching processes in action, such as strangulation and tidal- and ram-pressure stripping, (2) directly account for the baryon content of galaxy clusters in the cold-warm (T<10^6 K) gas, (3) determine the phase structure and kinematics of gas participating in the equilibrium-regulating exchange of energy at the cores of galaxy clusters, and (4) map cold streams and filaments of the Cosmic Web that feed galaxies and clusters. With a substantial UV undertaking beyond the Hubble Space Telescope, all of the above would be achievable over the entire epoch of galaxy cluster formation. Such capabilities, coupled with already-planned advancements at other wavelengths, will transform extragalactic astronomy by revealing the dominant formation and growth mechanisms of gaseous halos over the mass spectrum, settling the debate between early- and late-time metal enrichment scenarios, and revealing how the ecosystems in which galaxies reside ultimately facilitate their demise.

Asymptotic Giant Branch (AGB) winds from evolved stars not only provide a non-trivial amount of mass and energy return, but also produce dust grains in massive elliptical galaxies. Due to the fast stellar velocity and the high ambient temperature, the wind is thought to form a comet-like tail, similar to Mira in the Local Bubble. Many massive elliptical galaxies and cluster central galaxies host extended dusty cold filaments. The fate of the cold dusty stellar wind and its relation to cold filaments are not well understood. In this work, we carry out both analytical and numerical studies of the interaction between an AGB wind and the surrounding hot gas. We find that the cooling time of the tail is inversely proportional to the ambient pressure. In the absence of cooling, or in low pressure environments (e.g., the outskirts of elliptical galaxies), AGB winds are quickly mixed into the hot gas, and all the AGB winds have similar appearance and head-to-tail ratio. In high pressure environments, such as the Local Bubble and the central regions of massive elliptical galaxies, some of the gas in the mixing layer between the stellar wind and the surrounding hot gas can cool efficiently and cause the tail to become longer. Our simulated tail of Mira itself has similar length and velocity to that observed, and appears similar to the simulated AGB tail in the central regions of massive galaxies. We speculate that instead of thermal instability, the induced condensation at the mixing layer of AGB winds may be the origin of cold filaments in massive galaxies and galaxy clusters. This naturally explains the existence of dust and PAH in the filaments.

The interaction between inflowing gas clouds and galactic outflows at the interface where the galactic disk transitions into the circumgalactic medium is an important process in galaxy fueling, yet remains poorly understood. Using a series of tall-box hydrodynamic ENZO simulations, we have studied the interaction between smooth gas inflow and supernovae-driven outflow at the disk-halo interface with pc-scale resolution. A realistic wind of outflowing material is generated by supernovae explosions in the disk, while inflowing gas is injected at the top boundary of the simulation box with an injection velocity ranging from $10-100 \rm \ km \ s^{-1}$. We find that cooling and hydrodynamic instabilities drive the injected gas to fragment into cold ($\sim 10^{3}$ K) cloud clumps with typical densities of $\sim 1 \rm \ cm^{-3}$. These clumps initially accelerate before interacting and partially mixing with the outflow and decelerating to velocities in the 50-100 $\rm km \ s^{-1}$ range. When the gas clumps hit the disk, $10\%-50 \%$ of the injected material is able to accrete (depending on the injection velocity). Clumps originating from gas injected with a higher initial velocity approach the disk with greater ram pressure, allowing them to penetrate through the disk in low density regions. We use (equilibrium) CLOUDY photoionization models to generate absorption and emission signatures of gas accretion, finding that our mock HI and H$\alpha$ observables are prominent and generally consistent with measurements in the Milky Way. We do not predict enhanced emission/absorption for higher ionization states such as OVI.

We study the chaotic-like behavior of cosmological simulations by quantifying how minute perturbations grow over time and manifest as macroscopic differences in galaxy properties. When we run pairs of 'shadow' simulations that are identical except for random minute initial displacements to particle positions (e.g. of order 1e-7pc), the results diverge from each other at the individual galaxy level (while the statistical properties of the ensemble of galaxies are unchanged). After cosmological times, the global properties of pairs of 'shadow' galaxies that are matched between the simulations differ from each other generally at a level of ~2-25%, depending on the considered physical quantity. We perform these experiments using cosmological volumes of (25-50Mpc/h)^3 evolved either purely with dark matter, or with baryons and star-formation but no feedback, or using the full feedback model of the IllustrisTNG project. The runs cover four resolution levels spanning a factor of 512 in mass. We find that without feedback the differences between shadow galaxies generally become smaller as the resolution increases, but with the IllustrisTNG model the results are mostly converging towards a 'floor'. This hints at the role of feedback in setting the chaotic properties of galaxy formation. Importantly, we compare the macroscopic differences between shadow galaxies to the overall scatter in various galaxy scaling relations, and conclude that for the star formation-mass and the Tully-Fisher relations the butterfly effect in our simulations contributes significantly to the overall scatter. We find that our results are robust to whether random numbers are used in the sub-grid models or not. We discuss the implications for galaxy formation theory in general and for cosmological simulations in particular.

Using a high resolution simulation of an isolated dwarf galaxy, accounting for multi-channel stellar feedback and chemical evolution on a star-by-star basis, we investigate how each of 15 metal species are distributed within our multi-phase interstellar medium (ISM) and ejected from our galaxy by galactic winds. For the first time, we demonstrate that the mass fraction probability distribution functions (PDFs) of individual metal species in the ISM are well described by a piecewise log-normal and power-law distribution. The PDF properties vary within each ISM phase. Hot gas is dominated by recent enrichment, with a significant power-law tail to high metal fractions, while cold gas is predominately log-normal. In addition, elements dominated by asymptotic giant branch (AGB) wind enrichment (e.g. N and Ba) mix less efficiently than elements dominated by supernova enrichment (e.g. $\alpha$ elements and Fe). This result is driven by the differences in source energetics and source locations, particularly the higher chance compared to massive stars for AGB stars to eject material into cold gas. Nearly all of the produced metals are ejected from the galaxy (only 4% are retained), but over 20% of metals dominated by AGB enrichment are retained. In dwarf galaxies, therefore, elements synthesized predominately through AGB winds should be both overabundant and have a larger spread compared to elements synthesized in either core collapse or Type Ia supernovae. We discuss the observational implications of these results, their potential use in developing improved models of galactic chemical evolution, and their generalization to more massive galaxies.

We investigate nonlinear structure formation in the fuzzy dark matter (FDM) model using both numerical and perturbative techniques. On the numerical side, we examine the virtues and limitations of a Schrodinger-Poisson solver (wave formulation) versus a fluid dynamics solver (Madelung formulation). We also carry out a perturbative computation of the one-loop mass power spectrum. We find that (1) in many cases, the fluid dynamics solver is capable of producing the expected interference patterns, but it fails where destructive interference causes the density to vanish which generally occurs in the nonlinear regime. (2) The Schrodinger-Poisson solver works well in all test cases, but it is demanding in resolution: one must resolve the small de Broglie scale to obtain the correct dynamics on large scales. (3) We compare the mass power spectrum from perturbation theory against that from the Schrodinger-Poisson solver, and find good agreement in the mildly nonlinear regime. Compared with fluid perturbation theory, wave perturbation theory has a more limited range of validity. (4) As an application, we compare the Lyman-alpha forest flux power spectrum obtained from the Schrodinger-Poisson solver versus one from an N-body simulation (which is often used as an approximate method to make predictions for FDM). At redshift 5, the two, starting from the same initial condition, agree to better than 10 % on observationally relevant scales as long as the FDM mass exceeds $2 \times 10^{-23}$ eV.

In order to better understand the relationship between feedback and galactic chemical evolution, we have developed a new model for stellar feedback at grid resolutions of only a few parsecs in global disk simulations, using the adaptive mesh refinement hydrodynamics code Enzo. For the first time in galaxy scale simulations, we simulate detailed stellar feedback from individual stars including asymptotic giant branch winds, photoelectric heating, Lyman-Werner radiation, ionizing radiation tracked through an adaptive ray-tracing radiative transfer method, and core collapse and Type Ia supernovae. We furthermore follow the star-by-star chemical yields using tracer fields for 15 metal species: C, N, O, Na, Mg, Si, S, Ca, Mn, Fe, Ni, As, Sr, Y, and Ba. We include the yields ejected in massive stellar winds, but greatly reduce the winds' velocities due to computational constraints. We describe these methods in detail in this work and present the first results from 500~Myr of evolution of an isolated dwarf galaxy with properties similar to a Local Group, low-mass dwarf galaxy. We demonstrate that our physics and feedback model is capable of producing a dwarf galaxy whose evolution is consistent with observations in both the Kennicutt-Schmidt relationship and extended Schmidt relationship. Effective feedback drives outflows with a greater metallicity than the ISM, leading to low metal retention fractions consistent with observations. Finally, we demonstrate that these simulations yield valuable information on the variation in mixing behavior of individual metal species within the multi-phase interstellar medium.

Effective stellar feedback is used in models of galaxy formation to drive realistic galaxy evolution. Models typically include energy injection from supernovae as the dominant form of stellar feedback, often in some form of sub-grid recipe. However, it has been recently suggested that pre-SN feedback (stellar winds or radiation) is necessary in high-resolution simulations of galaxy evolution to properly regulate star formation and properties of the interstellar medium (ISM). Following these processes is computationally challenging, so many prescriptions model this feedback approximately, accounting for the local destruction of dense gas clouds around newly formed stars in lieu of a full radiative transfer calculation. In this work we examine high resolution simulations (1.8~pc) of an isolated dwarf galaxy with detailed stellar feedback tracked on a star-by-star basis. By following stellar ionizing radiation with an adaptive ray-tracing radiative transfer method, we test its importance in regulating star formation and driving outflows in this galaxy. We find that including ionizing radiation reduces the star formation rate (SFR) by over a factor of 5, and is necessary to produce the ISM conditions needed for supernovae to drive significant outflows. We find that a localized approximation for radiation feedback is sufficient to regulate the SFR on short timescales, but does not allow significant outflows. Short and long range radiation effects are both important in driving the evolution of our low-metallicity, low-mass dwarf galaxy. Generalizing these results to more massive galaxies would be a valuable avenue of future research.

We explore the relative role of small-scale fragmentation and global collapse in low-metallicity clouds, pointing out that in such clouds the cooling time may be longer than the dynamical time, allowing the cloud to collapse globally before it can fragment. This, we suggest, may help to explain the formation of the low-metallicity globular cluster population, since such dense stellar systems need a large amount of gas to be collected in a small region (without significant feedback during the collapse). To explore this further, we carry out numerical simulations of low-metallicity Bonner-Ebert stable gas clouds, demonstrating that there exists a critical metallicity (between 0.001 and 0.01 $Z_\odot$) below which the cloud collapses globally without fragmentation. We also run simulations including a background radiative heating source, showing that this can also produce clouds that do not fragment, and that the critical metallicity -- which can exceed the no-radiation case -- increases with the heating rate.

We have developed a semi-analytic framework to model the large-scale evolution of the first Population III (Pop III) stars and the transition to metal-enriched star formation. Our model follows dark matter halos from cosmological N-body simulations, utilizing their individual merger histories and three-dimensional positions, and applies physically motivated prescriptions for star formation and feedback from Lyman-Werner (LW) radiation, hydrogen ionizing radiation, and external metal enrichment due to supernovae winds. This method is intended to complement analytic studies, which do not include clustering or individual merger histories, and hydrodynamical cosmological simulations, which include detailed physics, but are computationally expensive and have limited dynamic range. Utilizing this technique, we compute the cumulative Pop III and metal-enriched star formation rate density (SFRD) as a function of redshift at $z \geq 20$. We find that varying the model parameters leads to significant qualitative changes in the global star formation history. The Pop III star formation efficiency and the delay time between Pop III and subsequent metal-enriched star formation are found to have the largest impact. The effect of clustering (i.e. including the three-dimensional positions of individual halos) on various feedback mechanisms is also investigated. The impact of clustering on LW and ionization feedback is found to be relatively mild in our fiducial model, but can be larger if external metal enrichment can promote metal-enriched star formation over large distances.

We present the science case and observations plan of the MeerKAT Fornax Survey, an HI and radio continuum survey of the Fornax galaxy cluster to be carried out with the SKA precursor MeerKAT. Fornax is the second most massive cluster within 20 Mpc and the largest nearby cluster in the southern hemisphere. Its low X-ray luminosity makes it representative of the environment where most galaxies live and where substantial galaxy evolution takes place. Fornax's ongoing growth makes it an excellent laboratory for studying the assembly of clusters, the physics of gas accretion and stripping in galaxies falling in the cluster, and the connection between these processes and the neutral medium in the cosmic web. We will observe a region of 12 deg$^2$ reaching a projected distance of 1.5 Mpc from the cluster centre. This will cover a wide range of environment density out to the outskirts of the cluster, where gas-rich in-falling groups are found. We will: study the HI morphology of resolved galaxies down to a column density of a few times 1e+19 cm$^{-2}$ at a resolution of 1 kpc; measure the slope of the HI mass function down to M(HI) 5e+5 M(sun); and attempt to detect HI in the cosmic web reaching a column density of 1e+18 cm$^{-2}$ at a resolution of 10 kpc.

We provide an analytic framework for interpreting observations of multiphase circumgalactic gas that is heavily informed by recent numerical simulations of thermal instability and precipitation in cool-core galaxy clusters. We start by considering the local conditions required for the formation of multiphase gas via two different modes: (1) uplift of ambient gas by galactic outflows, and (2) condensation in a stratified stationary medium in which thermal balance is explicitly maintained. Analytic exploration of these two modes provides insights into the relationships between the local ratio of the cooling and freefall time scales (i.e., t_cool / t_ff), the large-scale gradient of specific entropy, and development of precipitation and multiphase media in circumgalactic gas. We then use these analytic findings to interpret recent simulations of circumgalactic gas in which global thermal balance is maintained. We show that long-lasting configurations of gas with 5 < t_cool / t_ff < 20 and radial entropy profiles similar to observations of local cool-core galaxy cluster cores are a natural outcome of precipitation-regulated feedback. We conclude with some observational predictions that follow from these models. This work focuses primarily on precipitation and AGN feedback in galaxy cluster cores, because that is where the observations of multiphase gas around galaxies are most complete. However, many of the physical principles that govern condensation in those environments apply to circumgalactic gas around galaxies of all masses.

Recent observational evidence for initial mass function (IMF) variations in massive quiescent galaxies at $z = 0$ challenges the long-established paradigm of a universal IMF. While a few theoretical models relate the IMF to birth cloud conditions, the physical driver underlying these putative IMF variations is still largely unclear. Here we use post-processing analysis of the Illustris cosmological hydrodynamical simulation to investigate possible physical origins of IMF variability with galactic properties. We do so by tagging stellar particles in the simulation (each representing a stellar population of $\approx10^{6}~\mathrm{M}_{\odot}$) with individual IMFs that depend on various physical conditions, such as velocity dispersion, metallicity, or SFR, at the time and place the stars are formed. We then follow the assembly of these populations throughout cosmic time, and reconstruct the overall IMF of each $z=0$ galaxy from the many distinct IMFs it is comprised of. Our main result is that applying the observed relations between IMF and galactic properties to the conditions at the star-formation sites does not result in strong enough IMF variations between $z = 0$ galaxies. Steeper physical IMF relations are required for reproducing the observed IMF trends, and some stellar populations must form with more extreme IMFs than those observed. The origin of this result is the hierarchical nature of massive galaxy assembly, and it has implications for the reliability of the strong observed trends, for the ability of cosmological simulations to capture certain physical conditions in galaxies, and for theories of star-formation aiming to explain the physical origin of a variable IMF.

Galactic outflows are ubiquitously observed in star-forming disk galaxies and are critical for galaxy formation. Supernovae (SNe) play the key role in driving the outflows, but there is no consensus as to how much energy, mass and metal they can launch out of the disk. We perform 3D, high-resolution hydrodynamic simulations to study SNe-driven outflows from stratified media. Assuming SN rate scales with gas surface density $\Sigma_{\rm{gas}}$ as in the Kennicutt-Schmidt (KS) relation, we find the mass loading factor, $\eta_m$, defined as the mass outflow flux divided by the star formation surface density, decreases with increasing $\Sigma_{\rm{gas}}$ as $\eta_m \propto \Sigma^{-0.61}_{\rm{gas}}$. Approximately $\Sigma_{\rm{gas}}$ $\lesssim$ 50 $M_\odot/\rm{pc}^2$ marks when $\eta_m \gtrsim$1. About 10-50\% of the energy and 40-80\% of the metals produced by SNe end up in the outflows. The tenuous hot phase ($T>3\times 10^5$ K), which fills 60-80\% of the volume at mid-plane, carries the majority of the energy and metals in outflows. We discuss how various physical processes, including vertical distribution of SNe, photoelectric heating, external gravitational field and SN rate, affect the loading efficiencies. The relative scale height of gas and SNe is a very important factor in determining the loading efficiencies.

We utilize cosmological hydrodynamic simulations to study the formation of Population III (Pop III) stars in dark matter halos exposed to strong ionizing radiation. We simulate the formation of three halos subjected to a wide range of ionizing fluxes, and find that for high flux, ionization and photoheating can delay gas collapse and star formation up to halo masses significantly larger than the atomic cooling threshold. The threshold halo mass at which gas first collapses and cools increases with ionizing flux for intermediate values, and saturates at a value approximately an order of magnitude above the atomic cooling threshold for extremely high flux (e.g. $\approx 5 \times 10^8 ~ M_\odot$ at $z\approx6$). This behavior can be understood in terms of photoheating, ionization/recombination, and Ly$\alpha$ cooling in the pressure-supported, self-shielded gas core at the center of the growing dark matter halo. We examine the spherically-averaged radial velocity profiles of collapsing gas and find that a gas mass of up to $\approx 10^{6}~ M_\odot$ can reach the central regions within $3~{\rm Myr}$, providing an upper limit on the amount of massive Pop III stars that can form. The ionizing radiation increases this limit by a factor of a few compared to strong Lyman-Werner (LW) radiation alone. We conclude that the bright HeII 1640 \AA\ emission recently observed from the high-redshift galaxy CR7 cannot be explained by Pop III stars alone. However, in some halos, a sufficient number of Pop III stars may form to be detectable with future telescopes such as the James Webb Space Telescope (JWST).

Observed galactic disks have specific angular momenta similar to expectations for typical dark matter halos in $\Lambda$CDM. Cosmological hydrodynamical simulations have recently reproduced this similarity in large galaxy samples by including strong galactic winds, but the exact mechanism that achieves this is not yet clear. Here we present an analysis of key aspects contributing to this relation: angular momentum selection and evolution of Lagrangian mass elements as they accrete onto dark matter halos, condense into Milky Way-scale galaxies, and join the $z=0$ stellar phase. We contrast this evolution in the Illustris simulation with that in a simulation without galactic winds, where the $z=0$ angular momentum is $\approx0.6$ dex lower. We find that winds induce differences between these simulations in several ways: increasing angular momentum, preventing angular momentum loss, and causing $z=0$ stars to sample the accretion-time angular momentum distribution of baryons in a biased way. In both simulations, gas loses on average $\approx0.4$ dex between accreting onto halos and first accreting onto central galaxies. In Illustris, this is followed by $\approx0.2$ dex gains in the `galactic wind fountain' and no further net evolution past the final accretion onto the galaxy. Without feedback, further losses of $\approx0.2$ dex occur in the gas phase inside the galaxies. An additional $\approx0.15$ dex difference arises from feedback preferentially selecting higher angular momentum gas at accretion by expelling gas that is poorly aligned. These and additional effects of similar magnitude are discussed, suggesting a complex origin of the similarity between the specific angular momenta of galactic disks and typical halos.

The Direct Collapse Black Hole (DCBH) scenario provides a solution for forming the massive black holes powering bright quasars observed in the early Universe. A prerequisite for forming a DCBH is that the formation of (much less massive) Population III stars be avoided - this can be achieved by destroying H$_2$ via Lyman-Werner (LW) radiation (E$_{\rm{LW}}$ = 12.6 eV). We find that two conditions must be met in the proto-galaxy that will host the DCBH. First, prior star formation must be delayed; this can be achieved with a background LW flux of J$_{\rm BG} \gtrsim 100\ J_{21}$. Second, an intense burst of LW radiation from a neighbouring star-bursting proto-galaxy is required, just before the gas cloud undergoes gravitational collapse, to finally suppress star formation completely. We show here for the first time using high-resolution hydrodynamical simulations, including full radiative transfer, that this low-level background, combined with tight synchronisation and irradiation of a secondary proto-galaxy by a primary proto-galaxy, inevitably moves the secondary proto-galaxy onto the isothermal atomic cooling track, without the deleterious effects of either photo-evaporating the gas or polluting it by heavy elements. These, atomically cooled, massive proto-galaxies are expected to ultimately form a DCBH of mass $10^4 - 10^5 M_{\odot}$.

Suppression of H2-cooling in early protogalaxies has important implications for the formation of supermassive black holes seeds, the first generation of stars, and the epoch of reionization. This suppression can occur via photodissociation of H2 (by ultraviolet Lyman-Werner [LW] photons) or by photodetachment of H-, a precursor in H2 formation (by infrared [IR] photons). Previous studies have typically adopted idealised spectra, with a blackbody or a power-law shape, in modeling the chemistry of metal-free protogalaxies, and utilised a single parameter, the critical UV flux, or Jcrit, to determine whether H2-cooling is prevented. This can be misleading, and that independent of the spectral shape, there is a critical curve in the (kLW , kH^- ) plane, where kLW and kH^- are the H2-dissocation rates by LW and IR photons, which determines whether a protogalaxy can cool below ~1000 Kelvin. We use a one-zone model to follow the chemical and thermal evolution of gravitationally collapsing protogalactic gas, to compute this critical curve, and provide an accurate analytical fit for it. We improve on previous works by considering a variety of more realistic Pop III or Pop II-type spectra from population synthesis models and perform fully frequency-dependent calculations of the H2-photodissociation rates for each spectrum. We compute the ratio kLW/kH^- for each spectrum, as well as the minimum stellar mass M_star, for various IMFs and metallicities, required to prevent cooling in a neighboring halo a distance d away. We provide critical M_star/d2 values for suppression of H2-cooling, with analytic fits, which can be used in future studies.

The OVI $\lambda\lambda$1032, 1038\AA\ doublet emission traces collisionally ionized gas with $T\approx 10^{5.5}$ K, where the cooling curve peaks for metal-enriched plasma. This warm-hot phase is usually not well-resolved in numerical simulations of the multiphase interstellar medium (ISM), but can be responsible for a significant fraction of the emitted energy. Comparing simulated OVI emission to observations is therefore a valuable test of whether simulations predict reasonable cooling rates from this phase. We calculate OVI $\lambda$1032\AA\ emission, assuming collisional ionization equilibrium, for our small-box simulations of the stratified ISM regulated by supernovae. We find that the agreement is very good for our solar neighborhood model, both in terms of emission flux and mean OVI density seen in absorption. We explore runs with higher surface densities and find that, in our simulations, the OVI emission from the disk scales roughly linearly with the star formation rate. Observations of OVI emission are rare for external galaxies, but our results do not show obvious inconsistency with the existing data. Assuming the solar metallicity, OVI emission from the galaxy disk in our simulations accounts for roughly 0.5\% of supernovae heating.

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.

We analyze heating and cooling processes in an idealized simulation of a cool-core cluster, where momentum-driven AGN feedback balances radiative cooling in a time-averaged sense. We find that, on average, energy dissipation via shock waves is almost an order of magnitude higher than via turbulence. Most of the shock waves in the simulation are very weak shocks with Mach numbers smaller than 1.5, but the stronger shocks, although rare, dissipate energy more effectively. We find that shock dissipation is a steep function of radius, with most of the energy dissipated within 30 kpc, while radiative cooling loses area less concentrated. However, adiabatic processes and mixing (of post-shock materials and the surrounding gas) are able to redistribute the heat throughout the core. A considerable fraction of the AGN energy also escapes the core region. The cluster goes through cycles of AGN outbursts accompanied by periods of enhanced precipitation and star formation, over Gyr timescales. The cluster core is under-heated at the end of each cycle, but over-heated at the peak of the AGN outburst. During the heating-dominant phase, turbulent dissipation alone is often able to balance radiative cooling at every radius but, when this is occurs, shock waves inevitably dissipate even more energy. Our simulation explains why some clusters, such as Abell 2029, are cooling dominated, while in some other clusters, such as Perseus, various heating mechanisms including shock heating, turbulent dissipation and bubble mixing can all individually balance cooling, and together, overheat the core.

Over the last decade block-structured adaptive mesh refinement (SAMR) has found increasing use in large, publicly available codes and frameworks. SAMR frameworks have evolved along different paths. Some have stayed focused on specific domain areas, others have pursued a more general functionality, providing the building blocks for a larger variety of applications. In this survey paper we examine a representative set of SAMR packages and SAMR-based codes that have been in existence for half a decade or more, have a reasonably sized and active user base outside of their home institutions, and are publicly available. The set consists of a mix of SAMR packages and application codes that cover a broad range of scientific domains. We look at their high-level frameworks, and their approach to dealing with the advent of radical changes in hardware architecture. The codes included in this survey are BoxLib, Cactus, Chombo, Enzo, FLASH, and Uintah.

We explore the impact of star formation and thermal stellar feedback on the giant molecular cloud (GMC) population forming in a M83-type barred spiral galaxy. We compare three high-resolution simulations (1.5 pc cell size) with different star formation/feedback models: one with no star formation, one with star formation but no feedback, and one with star formation and thermal energy injection. We analyze the resulting population of clouds, finding that we can identify the same population of massive, virialized clouds and transient, low-surface density clouds found in our previous work (that did not include star formation or feedback). Star formation and feedback can affect the mix of clouds we identify. In particular, star formation alone simply converts dense cloud gas into stars with only a small change to the cloud populations, principally resulting in a slight decrease in the transient population. Feedback, however, has a stronger impact: while it is not generally sufficient to entirely destroy the clouds, it does eject gas out of them, increasing the gas density in the inter-cloud region. This decreases the number of massive clouds, but substantially increases the transient cloud population. We also find that feedback tends to drive a net radial inflow of massive clouds, leading to an increase in the star formation rate in the bar region. We examine a number of possible reasons for this and conclude that it is possible that the drag force from the enhanced intercloud density could be responsible.

We perform zoom-in cosmological simulations of a suite of dwarf galaxies, examining the impact of cosmic-rays generated by supernovae, including the effect of diffusion. We first look at the effect of varying the uncertain cosmic ray parameters by repeatedly simulating a single galaxy. Then we fix the comic ray model and simulate five dwarf systems with virial masses range from 8-30 $\times 10^{10}$ Msun. We find that including cosmic ray feedback (with diffusion) consistently leads to disk dominated systems with relatively flat rotation curves and constant star formation rates. In contrast, our purely thermal feedback case results in a hot stellar system and bursty star formation. The CR simulations very well match the observed baryonic Tully-Fisher relation, but have a lower gas fraction than in real systems. We also find that the dark matter cores of the CR feedback galaxies are cuspy, while the purely thermal feedback case results in a substantial core.

We explore the formation of massive high-redshift Population III (Pop III) galaxies through photoionization feedback. We consider dark matter halos formed from progenitors that have undergone no star formation as a result of early reionization and photoevaporation caused by a nearby galaxy. Once such a halo reaches $\approx 10^9~M_\odot$, corresponding to the Jeans mass of the photoheated intergalactic medium (IGM) at $z\approx 7$, pristine gas is able to collapse into the halo, potentially producing a massive Pop III starburst. We suggest that this scenario may explain the recent observation of strong He II 1640~\AA~line emission in CR7, which is consistent with $\sim 10^7~M_\odot$ of young Pop III stars. Such a large mass of Pop III stars is unlikely without the photoionization feedback scenario, because star formation is expected to inject metals into halos above the atomic cooling threshold ($\sim 10^8~M_\odot$ at $z \approx 7$). We use merger trees to analytically estimate the abundance of observable Pop III galaxies formed through this channel, and find a number density of $\approx 10^{-7}~{\rm Mpc^{-3}}$ at $z=6.6$ (the redshift of CR7). This is approximately a factor of ten lower than the density of Ly$\alpha$ emitters as bright as CR7.

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.

Supernovae (SN), the most energetic stellar feedback mechanism, are crucial for regulating the interstellar medium (ISM) and launching galactic winds. We explore how supernova remnants (SNRs) create a multiphase medium by performing 3D hydrodynamical simulations at various SN rates, $S$, and ISM average densities, $\bar{n}$. The evolution of a SNR in a self-consistently generated three-phase ISM is qualitatively different from that in a uniform or a two-phase warm/cold medium. By travelling faster and further in the low-density hot phase, the domain of a SNR increases by $>10^{2.5}$. Varying $\bar{n}$ and $S$, we find that a steady state can only be achieved when the hot gas volume fraction $f_{\rm{V,hot}}\lesssim 0.6 \pm 0.1 $. Above that level, overlapping SNRs render connecting topology of the hot gas, and the ISM is subjected to thermal runaway. Photoelectric heating (PEH) has a surprisingly strong impact on $f_{\rm{V,hot}}$. For $\bar{n}\gtrsim 3 \cm-3 $, a reasonable PEH rate is able to suppress the thermal runaway. Overall, we determine the critical SN rate for the onset of thermal runaway to be $S_{\rm{crit}} = 200 (\bar{n}/1\cm-3)^k (E_{\rm{SN}}/10^{51}\erg)^{-1} \kpc^{-3} \myr-1$, where $k = (1.2,2.7)$ for $\bar{n} \leq 1$ and $> 1\cm-3 $, respectively. We present a fitting formula of the ISM pressure $P(\bar{n}$, $S$), which can be used as an effective equation of state in cosmological simulations. Despite the 5 orders of magnitude span of $(\bar{n},S)$, the average Mach number varies little: $\mathcal{M} \approx \ 0.5\pm 0.2, \ 1.2\pm 0.3,\ 2.3\pm 0.9$ for the hot, warm and cold phases, respectively.

Recently, Planck measured a value of the cosmic microwave background (CMB) optical depth due to electron scattering of $\tau=0.066 \pm 0.016$. Here we show that this low value leaves essentially no room for an early partial reionisation of the intergalactic medium (IGM) by high-redshift Population III (Pop III) stars, expected to have formed in low-mass minihaloes. We perform semi-analytic calculations of reionisation which include the contribution from Pop II stars in atomic cooling haloes, calibrated with high-redshift galaxy observations, and Pop III stars in minihaloes with feedback due to Lyman-Werner (LW) radiation and metal enrichment. We find that without LW feedback or prompt metal enrichment (and assuming a minihalo escape fraction of 0.5) the Pop III star formation efficiency cannot exceed $\sim{\rm a~few}\times 10^{-4}$, without violating the constraints set by Planck data. This excludes massive Pop III star formation in typical $10^6 M_\odot$ minihaloes. Including LW feedback and metal enrichment alleviates this tension, allowing large Pop III stars to form early on before they are quenched by feedback. We find that the total density of Pop III stars formed across cosmic time is $\lesssim 10^{4-5}~M_\odot~{\rm Mpc^{-3}}$ and does not depend strongly on the feedback prescription adopted. Additionally, we perform a simple estimate of the possible impact on reionisation of X-rays produced by accretion onto black hole remnants of Pop III stars. We find that unless the accretion duty cycle is very low ($\lesssim 0.01$), this could lead to an optical depth inconsistent with Planck.

Understanding gas flows into and out of the most massive dark matter structures in our Universe, galaxy clusters, is fundamental to understanding their evolution. Gas in clusters is well studied in the hot ($>$ 10$^{6}$ K) and cold ($<$ 10$^{4}$ K) regimes, but the warm gas component (10$^{4}$ - 10$^{6}$ K) is poorly constrained. It is challenging to observe directly, but can be probed through Ly$\alpha$ absorption studies. We produce the first systematic study of the warm gas content of galaxy clusters through synthetic Ly$\alpha$ absorption studies using cosmological simulations of two galaxy clusters produced with Enzo. We explore the spatial and kinematic properties of our cluster absorbers, and show that the majority of the identified absorbers are due to fast moving gas associated with cluster infall from IGM filaments. Towards the center of the clusters, however, the warm IGM filaments are no longer dominant and the absorbers tend to have higher column densities and metallicities, representing stripped galaxy material. We predict that the absorber velocity distribution should generally be bi-modal and discuss the effects of cluster size, mass, and morphology on the properties of the identified absorbers, and the overall cluster warm gas content. We find tentative evidence for a change in the well known increasing N$_{HI}$ with decreasing impact parameter for the most massive dark matter halos. Our results are compared directly to observations of Ly$\alpha$ absorbers in the Virgo cluster, and provide predictions for future Ly$\alpha$ absorption studies.

Recent observations have constrained the orbit and structure of the Large Magellanic Cloud (LMC), implying a well-constrained pericentric passage about the Milky Way (MW) ~ 50 Myr ago. In this scenario, the LMC's gaseous disk has recently experienced maximal ram pressure stripping, suggesting the current extent of its HI disk directly probes the medium in which it is moving. From the observed stellar and HI distributions of the system we find evidence of a truncated gas profile along the windward ``leading edge' of the LMC disk, despite a far more extended stellar component. We explore the implications of this ram pressure stripping signature, using both analytic prescriptions and full 3-dimensional hydrodynamic simulations of the LMC. Our simulations subject the system to a headwind whose velocity components correspond directly to the recent orbital history of the LMC. We vary the density of this headwind, using a variety of sampled parameters for a Beta-profile for a theoretical MW circumgalactic medium (CGM), comparing the resulting HI morphology directly to observations of the LMC HI and stellar components. This model can match the radial extent of the LMC's leading (windward) edge only in scenarios where the MW CGM density at pericentric passage is n(R = 48.2 +/- 5 kpc) = 1.1 (+.44/-.45) x 1e-4 cm^-3. The implied pericentric density proves insensitive to both the broader CGM structure and temperature profile, thus providing a model-independent constraint on the local gas density. This result imposes an important constraint on the density profile of the MW's CGM, and thus the total baryon content of the MW. From our work, assuming a Beta-profile valid to ~ Rvir, we infer a total diffuse CGM mass M(300 kpc) = 2.6 +/- 1.4 x 1e10 Msun or approximately 15% of a 1e12 Msun MW's baryonic mass budget.

Constraining the properties of Population III (Pop III) stars will be very challenging because they reside in small galaxies at high redshift which will be difficult to directly detect. In this paper, we suggest that intensity mapping may be a promising method to study Pop III stars. Intensity mapping is a technique proposed to measure large-scale fluctuations of galaxy line emission in three dimensions without resolving individual sources. This technique is well suited for observing many faint galaxies because it can measure their cumulative emission even if they cannot be directly detected. We focus on intensity mapping of He II recombination lines. These lines are much stronger in Pop III stars than Pop II stars because the harder spectra of Pop III stars are expected to produce many He II ionizing photons. Measuring the He II 1640 \AA{} intensity mapping signal, along with the signals from other lines such as Ly$\alpha$, H$\alpha$, and metal lines, could give constraints on the initial mass function (IMF) and star formation rate density (SFRD) of Pop III stars as a function of redshift. To demonstrate the feasibility of these observations, we estimate the strength of the Pop III He II 1640 \AA{} intensity mapping signal from $z=10-20$. We show that at $z\approx10$, the signal could be measured accurately by two different hypothetical future instruments, one which cross-correlates He II 1640 \AA{} with CO(1-0) line emission from galaxies and the other with 21 cm emission from the intergalactic medium (IGM).

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.

Galaxy growth depends critically on the interplay between radiative cooling of cosmic gas and the resulting energetic feedback that cooling triggers. This interplay has proven exceedingly difficult to model, even with large supercomputer simulations, because of its complexity. Nevertheless, real galaxies are observed to obey simple scaling relations among their primary observable characteristics. Here we show that a generic emergent property of the interplay between cooling and feedback can explain the observed scaling relationships between a galaxy's stellar mass, its total mass, and its chemical enrichment level, as well as the relationship between the average orbital velocity of its stars and the mass of its central black hole. These relationships naturally result from any feedback mechanism that strongly heats a galaxy's circumgalactic gas in response to precipitation of colder clouds out of that gas, because feedback then suspends the gas in a marginally precipitating state.

Most of the massive elliptical galaxies in the universe stopped forming stars billions of years ago, even though plenty of hot gas remains available for star formation. Here we present compelling evidence indicating that quenching of star formation depends on both black-hole feedback and Type Ia supernova heating. We analyze Chandra X-ray observations of ten massive ellipticals, five with extended, potentially star-forming multiphase gas and five single-phase ellipticals with no star formation. The ratio of cooling time to freefall time at 1--10 kpc in the multiphase galaxies is tc/tff ~10, indicating that precipitation-driven feedback limits cooling but does not eliminate condensation. In the same region of the single-phase galaxies, the radial profiles of gas entropy are consistent with a thermally stable (tc/tff > 20) supernova-driven outflow that sweeps stellar ejecta out of the galaxy. However, in one of those single-phase ellipticals (NGC 4261) we find tc/tff < 10 at < 300 pc. Notably, its jets are ~50 times more powerful than in the other nine ellipticals, in agreement with models indicating that precipitation near the black hole should switch its fueling mode from Bondi-like accretion to cold chaotic accretion. We conclude by hypothesizing that particularly strong black-hole outbursts can shut off star formation in massive elliptical galaxies by boosting the entropy of the hot gas and flipping the system into the supernova-sweeping state.

Numerical simulations of active galactic nuclei (AGN) feedback in cool-core galaxy clusters have successfully avoided classical cooling flows, but often produce too much cold gas. We perform adaptive mesh simulations that include momentum-driven AGN feedback, self-gravity, star formation and stellar feedback, focusing on the interplay between cooling, AGN heating and star formation in an isolated cool-core cluster. Cold clumps triggered by AGN jets and turbulence form filamentary structures tens of kpc long. This cold gas feeds both star formation and the supermassive black hole (SMBH), triggering an AGN outburst that increases the entropy of the ICM and reduces its cooling rate. Within 1-2 Gyr, star formation completely consumes the cold gas, leading to a brief shutoff of the AGN. The ICM quickly cools and redevelops multiphase gas, followed by another cycle of star formation/AGN outburst. Within 6.5 Gyr, we observe three such cycles. There is good agreement between our simulated cluster and the observations of cool-core clusters. ICM cooling is dynamically balanced by AGN heating, and a cool-core appearance is preserved. The minimum cooling time to free-fall time ratio typically varies between a few and $\gtrsim 20$. The star formation rate (SFR) covers a wide range, from 0 to a few hundred $\rm M_{\odot}\, yr^{-1}$, with an average of $\sim 40 \,\rm M_{\odot}\, yr^{-1}$. The instantaneous SMBH accretion rate shows large variations on short timescales, but the average value correlates well with the SFR. Simulations without stellar feedback or self-gravity produce qualitatively similar results, but a lower SMBH feedback efficiency (0.1% compared to 1%) results in too many stars.

The universe's largest galaxies reside at the centers of galaxy clusters and are embedded in hot gas that, if left unchecked, would cool prodigiously and create many more new stars than are actually observed. Cooling can be regulated by feedback from accretion of cooling gas onto the central black hole, but requires an accretion rate finely tuned to the thermodynamic state of the hot gas. Theoretical models in which cold clouds precipitate out of the hot gas via thermal instability and accrete onto the black hole exhibit the necessary tuning. We have recently presented observational evidence showing that the abundance of cold gas in the central galaxy increases rapidly near the predicted threshold for instability. Here we present observations showing that this threshold extends over a large range in cluster radius, cluster mass, and cosmic time, and incorporate the precipitation threshold into a comprehensive framework of theoretical models for the thermodynamic state of hot gas in galaxy clusters. According to that framework, precipitation regulates star formation in some giant galaxies, while thermal conduction prevents star formation in others, if it can compensate for radiative cooling and shut off precipitation.

The first generation of stars produces a background of Lyman-Werner (LW) radiation which can photo-dissociate molecular hydrogen, increasing the mass of dark matter halos required to host star formation. Previous studies have determined the critical mass required for efficient molecular cooling with a constant LW background. However, the true background is expected to increase rapidly at early times. Neglecting this evolution could underestimate star formation in small halos that may have started to cool in the past when the LW intensity was much lower. Background evolution is a large source of uncertainty in pre-reionization predictions of the cosmological 21cm signal, which can be observed with future radio telescopes. To address this, we perform zero-dimentional one-zone calculations that follow the density, chemical abundances, and temperature of gas in the central regions of dark matter halos, including hierarchical growth and an evolving LW background. We begin by studying the physics of halos subjected to a background that increases exponentially with redshift. We find that when the intensity increases more slowly than $J_{\rm LW}(z) \propto 10^{-z/5}$, cooling in the past is a relatively small effect. We then self-consistently compute the cosmological LW background over $z=15-50$ and find that cooling in the past due to an evolving background has a modest impact. Finally, we compare these results to three-dimensional hydrodynamical cosmological simulations with varying LW histories. While only a small number of halos were simulated, the results are consistent with our one-zone calculations.

High-redshift quasar observations imply that supermassive black holes (SMBHs) larger than $\sim 10^9 ~ M_\odot$ formed before $z=6$. That such large SMBHs formed so early in the Universe remains an open theoretical problem. One possibility is that gas in atomic cooling halos exposed to strong Lyman-Werner (LW) radiation forms $10^4-10^6 ~ M_\odot$ supermassive stars which quickly collapse into black holes. We propose a scenario for direct collapse black hole (DCBH) formation based on synchronized pairs of pristine atomic cooling halos. We consider halos at very small separation with one halo being a subhalo of the other. The first halo to surpass the atomic cooling threshold forms stars. Soon after these stars are formed, the other halo reaches the cooling threshold and due to its small distance from the newly formed galaxy, is exposed to the critical LW intensity required to form a DCBH. The main advantage of this scenario is that synchronization can potentially prevent photoevaporation and metal pollution in DCBH-forming halos. Since the halos reach the atomic cooling threshold at nearly the same time, the DCBH-forming halo is only exposed to ionizing radiation for a brief period. Tight synchronization could allow the DCBH to form before stars in the nearby galaxy reach the end of their lives and generate supernovae winds. We use N-body simulations to estimate the abundance of DCBHs formed in this way. The largest source of uncertainty in our estimate is the initial mass function (IMF) of metal free stars formed in atomic cooling halos. We find that even for tight synchronization, the density of DCBHs formed in this scenario could explain the SMBHs implied by $z=6$ quasar observations. Metal pollution and photoevaporation could potentially reduce the abundance of DCBHs below that required to explain the observations in other models that rely on a high LW flux.

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 describe and implement an adaptive particle-mesh algorithm to solve the Poisson equation for grid-based hydrodynamics codes with nested grids. The algorithm is implemented and extensively tested within the astrophysical code Enzo against the multigrid solver available by default. We find that while both algorithms show similar accuracy for smooth mass distributions, the adaptive particle-mesh algorithm is more accurate for the case of point masses, and is generally less noisy. We also demonstrate that the two-body problem can be solved accurately in a configuration with nested grids. In addition, we discuss the effect of subcycling, and demonstrate that evolving all the levels with the same timestep yields even greater precision.

Explaining the existence of supermassive black holes (SMBHs) larger than $\sim 10^9 M_\odot$ at redshifts $z >\sim 6$ remains an open theoretical question. One possibility is that gas collapsing rapidly in pristine atomic cooling halos ($T_{\rm vir} >\sim 10^4 \rm{K}$) produces $10^4-10^6 M_\odot$ black holes. Previous studies have shown that the formation of such a black hole requires a strong UV background to prevent molecular hydrogen cooling and gas fragmentation. Recently it has been proposed that a high UV background may not be required for halos that accrete material extremely rapidly or for halos where gas cooling is delayed due to a high baryon-dark matter streaming velocity. In this work, we point out that building up a halo with $T_{\rm vir} >\sim 10^4 \rm{K}$ before molecular cooling becomes efficient is not sufficient for forming a direct collapse black hole (DCBH). Though molecular hydrogen formation may be delayed, it will eventually form at high densities leading to efficient cooling and fragmentation. The only obvious way that molecular cooling could be avoided in the absence of strong UV radiation, is for gas to reach high enough density to cause collisional dissociation of molecular hydrogen ($\sim 10^4 ~ {\rm cm}^{-3}$) before cooling occurs. However, we argue that the minimum core entropy, set by the entropy of the intergalactic medium (IGM) when it decouples from the CMB, prevents this from occurring for realistic halo masses. This is confirmed by hydrodynamical cosmological simulations without radiative cooling. We explain the maximum density versus halo mass in these simulations with simple entropy arguments. The low densities found suggest that DCBH formation indeed requires a strong UV background.

We study the long-term evolution of an idealized cool-core galaxy cluster under the influence of momentum-driven AGN feedback using three-dimensional high-resolution (60 pc) adaptive mesh refinement (AMR) simulations. The momentum-driven AGN feedback is modeled with a pair of (small-angle) precessing jets, and the jet power is calculated based on the accretion rate of the cold gas in the vicinity of the Supermassive Black Hole (SMBH). The ICM first cools into clumps along the propagation direction of the AGN jets. As the jet power increases, gas condensation occurs isotropically, forming spatially extended (up to a few tens kpc) structures that resemble the observed $\rm H\alpha$ filaments in Perseus and many other cool-core cluster. Jet heating elevates the gas entropy and cooling time, halting clump formation. The cold gas that is not accreted onto the SMBH settles into a rotating disk of $\sim 10^{11}$ M$_{\odot}$. The hot gas cools directly onto the cold disk while the SMBH accretes from the innermost region of the disk, powering the AGN that maintains a thermally balanced steady state for a few Gyr. The mass cooling rate averaged over 7 Gyr is $\sim 30$ M$_{\odot}$/yr, an order of magnitude lower than the classic cooling flow value (which we obtain in runs without the AGN). Medium resolution simulations produce similar results, but when the resolution is lower than 0.5 kpc, the cluster experiences cycles of gas condensation and AGN outbursts. Owing to its self-regulating mechanism, AGN feedback can successfully balance cooling with a wide range of model parameters. Besides suppressing cooling, our model produces cold structures in early stages (up to $\sim 2$ Gyr) that are in good agreement with the observations. However, the long-lived massive cold disk is unrealistic, suggesting that additional physical processes are still needed.

The presence of quasars at redshifts z > 6 indicates the existence of supermassive black holes (SMBHs) as massive as a few times 10^9 Msun, challenging models for SMBH formation. One pathway is through the direct collapse of gas in T_{vir} > 10^4 K halos; however, this requires the suppression of H_2 cooling to prevent fragmentation. In this paper, we examine a proposed new mechanism for this suppression which relies on cold-mode accretion flows leading to shocks at high densities (n > 10^4 cm^{-3}) and temperatures (T > 10^4 K). In such gas, H_2 is efficiently collisionally dissociated. We use high-resolution numerical simulations to test this idea, demonstrating that such halos typically have lower temperature progenitors, in which cooling is efficient. Those halos do show filamentary flows; however, the gas shocks at or near the virial radius (at low densities), thus preventing the proposed collisional mechanism from operating. We do find that, if we artificially suppress H_2 formation with a high UV background, so as to allow gas in the halo center to enter the high-temperature, high-density "zone of no return", it will remain there even if the UV flux is turned off, collapsing to high density at high temperature. Due to computational limitations, we simulated only three halos. However, we demonstrate, using Monte Carlo calculations of 10^6 halo merger histories, that a few rare halos could assemble rapidly enough to avoid efficient H_2 cooling in all of their progenitor halos, provided that the UV background exceeds J_{21} ~ few at redshifts as high as z ~ 20.

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 perform high-resolution (15-30 pc) adaptive mesh simulations to study the impact of momentum-driven AGN feedback in cool-core clusters, focusing in this paper on the formation of cold clumps. The feedback is jet-driven with an energy determined by the amount of cold gas within 500 pc of the SMBH. When the intra-cluster medium (ICM) in the core of the cluster becomes marginally stable to radiative cooling, with the thermal instability to the free-fall timescale ratio t_{TI}/t_{ff} < 3-10, cold clumps of gas start to form along the propagation direction of the AGN jets. By tracing the particles in the simulations, we find that these cold clumps originate from low entropy (but still hot) gas that is accelerated by the jet to outward radial velocities of a few hundred km/s. This gas is out of hydrostatic equilibrium and so can cool. The clumps then grow larger as they decelerate and fall towards the center of the cluster, eventually being accreted onto the super-massive black hole. The general morphology, spatial distribution and estimated H{\alpha} morphology of the clumps are in reasonable agreement with observations, although we do not fully replicate the filamentary morphology of the clumps seen in the observations, probably due to missing physics.

Galactic-scale winds are a generic feature of massive galaxies with high star formation rates across a broad range of redshifts. Despite their importance, a detailed physical understanding of what drives these mass-loaded global flows has remained elusive. In this paper, we explore the dynamical impact of cosmic rays by performing the first three-dimensional, adaptive mesh refinement simulations of an isolated starbursting galaxy that includes a basic model for the production, dynamics and diffusion of galactic cosmic rays. We find that including cosmic rays naturally leads to robust, massive, bipolar outflows from our 10^12 Msun halo, with a mass-loading factor Mout/SFR = 0.3 for our fiducial run. Other reasonable parameter choices led to mass-loading factors above unity. The wind is multiphase and is accelerated to velocities well in excess of the escape velocity. We employ a two-fluid model for the thermal gas and relativistic CR plasma and model a range of physics relevant to galaxy formation, including radiative cooling, shocks, self-gravity, star formation, supernovae feedback into both the thermal and CR gas, and isotropic CR diffusion. Injecting cosmic rays into star-forming regions can provide significant pressure support for the interstellar medium, suppressing star formation and thickening the disk. We find that CR diffusion plays a central role in driving superwinds, rapidly transferring long-lived CRs from the highest density regions of the disk to the ISM at large, where their pressure gradient can smoothly accelerate the gas out of the disk.

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.

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.

We simulate the evolution of a 10^9 Msun dark matter halo in a cosmological setting with an adaptive-mesh refinement code as an analogue to local low luminosity dwarf irregular and dwarf spheroidal galaxies. The primary goal of our study is to investigate the roles of reionization and supernova feedback in determining the star formation histories of low mass dwarf galaxies. We include a wide range of physical effects, including metal cooling, molecular hydrogen formation and cooling, photoionization and photodissociation from a metagalactic background, a simple prescription for self-shielding, star formation, and a simple model for supernova driven energetic feedback. We carry out simulations excluding each major effect in turn. We find that reionization is primarily responsible for expelling most of the gas in our simulations, but that supernova feedback is required to disperse the dense, cold gas in the core of the halo. Moreover, we show that the timing of reionization can produce an order of magnitude difference in the final stellar mass of the system. For our full physics run with reionization at z=9, we find a stellar mass of about 10^5 Msun at z=0, and a mass-to-light ratio within the half-light radius of approximately 130 Msun/Lsun, consistent with observed low-luminosity dwarfs. However, the resulting median stellar metallicity is 0.06 Zsun, considerably larger than observed systems. In addition, we find star formation is truncated between redshifts 4 and 7, at odds with the observed late time star formation in isolated dwarf systems but in agreement with Milky Way ultrafaint dwarf spheroidals. We investigate the efficacy of energetic feedback in our simple thermal-energy driven feedback scheme, and suggest that it may still suffer from excessive radiative losses, despite reaching stellar particle masses of about 100 Msun, and a comoving spatial resolution of 11 pc.

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.

The first systematic study of the warm gas (T=10^4-5 K) distribution across a galaxy cluster is presented using multiple background QSOs to the Virgo Cluster. We detect 25 Lya absorbers (N_HI = 10^13.1-15.4 cm^-2) in the Virgo velocity range toward 9 of 12 QSO sightlines observed with Cosmic Origin Spectrograph, with a cluster impact parameter range of 0.36-1.65 Mpc (0.23-1.05R_vir). Including 18 Lya absorbers previously detected by STIS or GHRS toward 7 of 11 background QSOs in and around the Virgo Cluster, we establish a sample of 43 absorbers towards a total of 23 background probes for studying the incidence of Lya absorbers in and around the Virgo Cluster. With these absorbers, we find: (1) warm gas is predominantly in the outskirts of the cluster and avoids the X-ray detected hot ICM. Also, Lya absorption strength increases with cluster impact parameter. (2) Lya absorbing warm gas traces cold HI emitting gas in the substructures of the Virgo Cluster. (3) Including the absorbers associated with the surrounding substructures, the warm gas covering fraction (100% for N_HI > 10^13.1 cm^-2) is in agreement with cosmological simulations. We speculate that the observed warm gas is part of large-scale gas flows feeding the cluster both the ICM and galaxies.

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 investigate the impact of star formation and feedback on ram pressure stripping using high-resolution adaptive mesh simulations, building on a previous series of papers that systematically investigated stripping using a realistic model for the interstellar medium, but without star formation. We find that star formation does not significantly affect the rate at which stripping occurs, and only has a slight impact on the density and temperature distribution of the stripped gas, indicating that our previous (gas-only) results are unaffected. For our chosen (moderate) ram pressure strength, stripping acts to truncate star formation in the disk over a few hundred million years, and does not lead to a burst of star formation. Star formation in the bulge is slightly enhanced, but the resulting change in the bulge-to-disk ratio is insignificant. We find that stars do form in the tail, primarily from gas that is ablated from the disk and the cools and condenses in the turbulent wake. The star formation rate in the tail is low, and any contribution to the intracluster light is likely to be very small. We argue that star formation in the tail depends primarily on the pressure in the intracluster medium, rather than the ram pressure strength. Finally, we compare to observations of star formation in stripped tails, finding that many of the discrepancies between our simulation and observed wakes can be accounted for by different intracluster medium pressures.

We study the buildup of magnetic fields during the formation of Population III star-forming regions, by conducting cosmological simulations from realistic initial conditions and varying the Jeans resolution. To investigate this in detail, we start simulations from identical initial conditions, mandating 16, 32 and 64 zones per Jeans length, and studied the variation in their magnetic field amplification. We find that, while compression results in some amplification, turbulent velocity fluctuations driven by the collapse can further amplify an initially weak seed field via dynamo action, provided there is sufficient numerical resolution to capture vortical motions (we find this requirement to be 64 zones per Jeans length, slightly larger than, but consistent with previous work run with more idealized collapse scenarios). We explore saturation of amplification of the magnetic field, which could potentially become dynamically important in subsequent, fully-resolved calculations. We have also identified a relatively surprising phenomena that is purely hydrodynamic: the higher-resolved simulations possess substantially different characteristics, including higher infall-velocity, increased temperatures inside 1000 AU, and decreased molecular hydrogen content in the innermost region. Furthermore, we find that disk formation is suppressed in higher-resolution calculations, at least at the times that we can follow the calculation. We discuss the effect this may have on the buildup of disks over the accretion history of the first clump to form as well as the potential for gravitational instabilities to develop and induce fragmentation.

We carry out high-resolution adaptive mesh refinement simulations of a cool core cluster, resolving the flow from Mpc scales down to pc scales. We do not (yet) include any AGN heating, focusing instead on cooling in order to understand how gas gets to the supermassive black hole (SMBH) at the center of the cluster. We find that, as the gas cools, the cluster develops a very flat temperature profile, undergoing a cooling catastrophe only in the central 10-100 pc of the cluster. Outside of this region, the flow is smooth, with no local cooling instabilities, and naturally produces very little low-temperature gas (below a few keV), in agreement with observations. The gas cooling in the center of the cluster rapidly forms a thin accretion disk. The amount of cold gas produced at the very center grows rapidly until a reasonable estimate of the resulting AGN heating rate (assuming even a moderate accretion efficiency) would overwhelm cooling. We argue that this naturally produces a thermostat which links the cooling of gas out to 100 kpc with the cold gas accretion in the central 100 pc, potentially closing the loop between cooling and heating. Isotropic heat conduction does not affect the result significantly, but we show that including the potential well of the brightest cluster galaxy is necessary to obtain the correct result. Also, we found that the outcome is sensitive to resolution, requiring very high mass resolution to correctly reproduce the small transition radius.

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 use three-dimensional hydrodynamical simulations to study the rapid infall phase of the common envelope interaction of a red giant branch star of mass equal to 0.88 \msun and a companion star of mass ranging from 0.9 down to 0.1 \msun. We first compare the results obtained using two different numerical techniques with different resolutions, and find overall very good agreement. We then compare the outcomes of those simulations with observed systems thought to have gone through a common envelope. The simulations fail to reproduce those systems in the sense that most of the envelope of the donor remains bound at the end of the simulations and the final orbital separations between the donor's remnant and the companion, ranging from 26.8 down to 5.9 \rsun, are larger than the ones observed. We suggest that this discrepancy vouches for recombination playing an essential role in the ejection of the envelope and/or significant shrinkage of the orbit happening in the subsequent phase.

The ability of primordial gas to cool in proto-galactic haloes exposed to Lyman-Werner (LW) radiation is critically dependent on the self-shielding of H_2. We perform radiative transfer calculations of LW line photons, post-processing outputs from three-dimensional adaptive mesh refinement (AMR) simulations of haloes with T_vir > 10^4 K at redshifts around z=10. We calculate the optically thick photodissociation rate numerically, including the effects of density, temperature, and velocity gradients in the gas, as well as line overlap and shielding of H_2 by HI, over a large number of sight-lines. In low-density regions (n<10^4 cm^-3) the dissociation rates exceed those obtained using most previous approximations by more than an order of magnitude; the correction is smaller at higher densities. We trace the origin of the deviations primarily to inaccuracies of (i) the most common fitting formula (Draine & Bertoldi 1996) for the suppression of the dissociation rate and (ii) estimates for the effective shielding column density from local properties of the gas. The combined effects of gas temperature and velocity gradients are comparatively less important, typically altering the spherically averaged rate only by a factor of less than two. We present a simple modification to the DB96 fitting formula for the optically thick rate which improves agreement with our numerical results to within approx. 15 per cent, and can be adopted in future simulations. We find that estimates for the effective shielding column can be improved by using the local Sobolev length. Our correction to the H_2 self-shielding reduces the critical LW flux to suppress H_2-cooling in T_vir>10^4 K haloes by an order of magnitude; this increases the number of such haloes in which supermassive (approx. M=10^5 M_sun) black holes may have formed.

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.

Some tails of ram-pressure stripped galaxies are detected in HI, some in Halpha, and some in X-ray (but never all three so far). We use numerical simulations to probe the conditions for the production of X-ray bright tails, demonstrating that the primary requirement is a high pressure intracluster medium (ICM). This is because the stripped tail is mostly in pressure equilibrium with the ICM, but mixing leaves it with densities and temperatures intermediate between the cold gas in the disk and the hot ICM. Given a high enough ICM pressure, this mixed gas lies in the X-ray bright region of the phase diagram. We compare the simulations to observations of the ram pressure stripped tail of ESO 137-001, showing excellent agreement in the total measured X-ray and Halpha emission and non-flaring morphology of the tail, and consistent HI measurements. Using these comparisons we constrain the level of mixing and efficiency of heat conduction in the intracluster medium (ICM)

Clusters of galaxies have been used extensively to determine cosmological parameters. A major difficulty in making best use of Sunyaev-Zel'dovich (SZ) and X-ray observations of clusters for cosmology is that using X-ray observations it is difficult to measure the temperature distribution and therefore determine the density distribution in individual clusters of galaxies out to the virial radius. Observations with the new generation of SZ instruments are a promising alternative approach. We use clusters of galaxies drawn from high-resolution adaptive mesh refinement (AMR) cosmological simulations to study how well we should be able to constrain the large-scale distribution of the intra-cluster gas (ICG) in individual massive relaxed clusters using AMiBA in its configuration with 13 1.2-m diameter dishes (AMiBA13) along with X-ray observations. We show that non-isothermal beta models provide a good description of the ICG in our simulated relaxed clusters. We use simulated X-ray observations to estimate the quality of constraints on the distribution of gas density, and simulated SZ visibilities (AMiBA13 observations) for constraints on the large-scale temperature distribution of the ICG. We find that AMiBA13 visibilities should constrain the scale radius of the temperature distribution to about 50% accuracy. We conclude that the upgraded AMiBA, AMiBA13, should be a powerful instrument to constrain the large-scale distribution of the ICG.

Accurate mass determination of clusters of galaxies is crucial if they are to be used as cosmological probes. However, there are some discrepancies between cluster masses determined based on gravitational lensing, and X-ray observations assuming strict hydrostatic equilibirium (i.e., the equilibrium gas pressure is provided entirely by thermal pressure). Cosmological simulations suggest that turbulent gas motions remaining from hierarchical structure formation may provide a significant contribution to the equilibrium pressure in clusters. We analyze a sample of massive clusters of galaxies drawn from high resolution cosmological simulations, and find a significant contribution (20%-45%) from non-thermal pressure near the center of relaxed clusters, and, in accord with previous studies, a minimum contribution at about 0.1 Rvir, growing to about 30%-45% at the virial radius, Rvir. Our results strongly suggest that relaxed clusters should have significant non-thermal support in their core region. As an example, we test the validity of strict hydrostatic equilibirium in the well-studied massive galaxy cluster Abell 1689 using the latest high resolution gravitational lensing and X-ray observations. We find a contribution of about 40% from non-thermal pressure within the core region of A1689, suggesting an alternate explanation for the mass discrepancy: the strict hydrostatic equilibirium is not valid in this region.

Observational studies of galaxy isophotal shapes have shown that galaxy orientations are anisotropic: a galaxy's long axis tends to be oriented toward the center of its host. This radial alignment is seen across a wide range of scales, from galaxies in massive clusters to small Milky Way type satellite systems. Recently, this effect has also been detected in dark matter simulations of cosmological structure, but the degree of alignment of dark matter substructures in these studies is significantly stronger than seen in observations. In this paper we attempt to reconcile these two results by performing high-resolution numerical experiments on N-body multi-component models of triaxial galaxies orbiting in an external analytical potential. The large number of particles employed allows us to probe deep into the inner structure of the galaxy: we show that the discrepancy between observed galaxies and simulated dark matter halos is a natural consequence of induced radial shape twisting in the galaxy by the external potential. The degree of twisting depends strongly on the orbital phase and eccentricity of the satellite, and it can, under certain conditions, be significant at radii smaller than the dark matter scale radius. Such internal misalignments will have important consequences, both for the dynamical evolution of the galaxy itself, and for mass modeling of galaxies in clustered environments.

Galaxies moving through the intracluster medium (ICM) of a cluster of galaxies can lose gas via ram pressure stripping. This stripped gas forms a tail behind the galaxy which is potentially observable. In this paper, we carry out hydrodynamical simulations of a galaxy undergoing stripping with a focus on the gas properties in the wake and their observational signatures. We include radiative cooling in an adaptive hydrocode in order to investigate the impact of a clumpy, multi-phase interstellar medium. We find that including cooling results in a very different morphology for the gas in the tail, with a much wider range of temperatures and densities. The tail is significantly narrower in runs with radiative cooling, in agreement with observed wakes. In addition, we make detailed predictions of H I, Halpha and X-ray emission for the wake, showing that we generally expect detectable H I and Halpha signatures, but no observable X-ray emission (at least for our chosen ram-pressure strength and ICM conditions). We find that the relative strength of the Halpha diagnostic depends somewhat on our adopted minimum temperature floor (below which we set cooling to zero to mimic physics processes not included in the simulation).

UV radiation from early astrophysical sources could have a large impact on subsequent star formation in nearby protogalaxies. Here we study the radiative feedback from the first, short-lived stars using hydrodynamical simulations with transient UV backgrounds (UVBs) and persistent Lyman-Werner backgrounds (LWBs) of varying intensity. We extend our prior work in Mesinger et al. (2006), by studying a more typical region whose proto-galaxies form at lower redshifts, z~13-20, in the epoch likely preceding the bulk of reionization. We confirm our previous results that feedback in the relic HII regions resulting from such transient radiation, is itself transient. Feedback effects dwindle away after ~30% of the Hubble time, and the same critical specific intensity of J_UV~0.1 x 10^{-21} ergs/s/cm^2/Hz/sr separates positive and negative feedback regimes. Additionally, we discover a second episode of eventual positive feedback in halos which have not yet collapsed when their progenitor regions were exposed to the transient UVB. This eventual positive feedback appears in all runs, regardless of the strength of the UVB. However, this feedback regime is very sensitive to the presence of Lyman-Werner radiation, and notable effects disappear under fairly modest background intensities of J_LW>10^{-3} x 10^{-21} ergs/s/cm^2/Hz/sr. We conclude that UV radiative feedback in relic HII regions, although a complicated process, seems unlikely to have a major impact on the progress of cosmological reionization, provided that present estimates of the lifetime and luminosity of a PopIII star are accurate. More likely is that the build-up of the LWB ultimately governs the feedback strength until a persistent UV background can be established. [abridged]

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.

In the absence of H_2 molecules, the primordial gas in early dark matter halos with virial temperatures just above T_vir >~ 10^4 K cools by collisional excitation of atomic H. Although it cools efficiently, this gas remains relatively hot, at a temperature near T ~ 8000 K, and consequently might be able to avoid fragmentation and collapse directly into a supermassive black hole (SMBH). In order for H_2--formation and cooling to be strongly suppressed, the gas must be irradiated by a sufficiently intense ultraviolet (UV) flux. We performed a suite of three--dimensional hydrodynamical adaptive mesh refinement (AMR) simulations of gas collapse in three different protogalactic halos with T_vir >~ 10^4 K, irradiated by a UV flux with various intensities and spectra. We determined the critical specific intensity, Jcrit, required to suppress H_2 cooling in each of the three halos. For a hard spectrum representative of metal--free stars, we find (in units of 10^{-21} erg s^{-1} Hz^{-1} sr^{-1} cm^{-2}) 10^4<Jcrit<10^5, while for a softer spectrum, which is characteristic of a normal stellar population, and for which H^{-} --dissociation is important, we find 30<Jcrit<300. These values are a factor of 3--10 lower than previous estimates. We attribute the difference to the higher, more accurate H_2 collisional dissociation rate we adopted. The reduction in Jcrit exponentially increases the number of rare halos exposed to super--critical radiation. When H_2 cooling is suppressed, gas collapse starts with a delay, but it ultimately proceeds more rapidly. The infall velocity is near the increased sound speed, and an object as massive as M ~ 10^5 solar mass may form at the center of these halos, compared to the M ~ 10^2 solar mass stars forming when H_2--cooling is efficient.

Great strides have been made in the last two decades in determining how galaxies evolve from their initial dark matter seeds to the complex structures we observe at z=0. The role of mergers has been documented through both observations and simulations, numerous satellites that may represent these initial dark matter seeds have been discovered in the Local Group, high redshift galaxies have been revealed with monstrous star formation rates, and the gaseous cosmic web has been mapped through absorption line experiments. Despite these efforts, the dark matter simulations that include baryons are still unable to accurately reproduce galaxies. One of the major problems is our incomplete understanding of how a galaxy accretes its baryons and subsequently forms stars. Galaxy formation simulations have been unable to accurately represent the required gas physics on cosmological timescales, and observations have only just begun to detect the star formation fuel over a range of redshifts and environments. How galaxies obtain gas and subsequently form stars is a major unsolved, yet tractable problem in contemporary extragalactic astrophysics. In this paper we outline how progress can be made in this area in the next decade.

Cold dark matter (CDM) hierarchical structure formation models predict the existence of large-scale accretion shocks between the virial and turnaround radii of clusters of galaxies. Kocsis et al. (2005) suggest that the Sunyaev-Zel'dovich (SZ) signal associated with such shocks might be observable with the next generation radio interferometer, ALMA. We study the three--dimensional distribution of accretion shocks around individual clusters of galaxies drawn from adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) simulations of LCDM (dark energy dominated CDM) models. In relaxed clusters, we find two distinct sets of shocks. One set ("virial shocks"), with Mach numbers of 2.5-4, is located at radii 0.9-1.3 Rvir, where Rvir is the spherical infall estimate of the virial radius, covering about 40-50% of the total surface area around clusters at these radii. Another set of stronger shocks ("external shocks") is located farther out, at about 3 Rvir, with large Mach numbers (~100), covering about 40-60% of the surface area. We simulate SZ surface brightness maps of relaxed massive galaxy clusters drawn from high resolution AMR runs, and conclude that ALMA should be capable of detecting the virial shocks in massive clusters of galaxies. More simulations are needed to improve estimates of astrophysical noise and to determine optimal observational strategies.

Cluster galaxies moving through the intracluster medium (ICM) are expected to lose some of their interstellar medium (ISM) through ISM-ICM interactions. We perform high resolution (40 pc) three-dimensional hydrodynamical simulations of a galaxy undergoing ram pressure stripping including radiative cooling in order to investigate stripping of a multiphase medium. The clumpy, multiphase ISM is self-consistently produced by the inclusion of radiative cooling, and spans six orders of magnitude in gas density. We find no large variations in the amount of gas lost whether or not cooling is involved, although the gas in the multiphase galaxy is stripped more quickly and to a smaller radius. We also see significant differences in the morphology of the stripped disks. This occurs because the multiphase medium naturally includes high density clouds set inside regions of lower density. We find that the lower density gas is stripped quickly from any radius of the galaxy, and the higher density gas can then be ablated. If high density clouds survive, through interaction with the ICM they lose enough angular momentum to drift towards the center of the galaxy where they are no longer stripped. Finally, we find that low ram pressure values compress gas into high density clouds that could lead to enhanced star formation, while high ram pressure leads to a smaller amount of high-density gas.

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.

We test four commonly used astrophysical simulation codes; Enzo, Flash, Gadget and Hydra, using a suite of numerical problems with analytic initial and final states. Situations similar to the conditions of these tests, a Sod shock, a Sedov blast and both a static and translating King sphere occur commonly in astrophysics, where the accurate treatment of shocks, sound waves, supernovae explosions and collapsed haloes is a key condition for obtaining reliable validated simulations. We demonstrate that comparable results can be obtained for Lagrangian and Eulerian codes by requiring that approximately one particle exists per grid cell in the region of interest. We conclude that adaptive Eulerian codes, with their ability to place refinements in regions of rapidly changing density, are well suited to problems where physical processes are related to such changes. Lagrangian methods, on the other hand, are well suited to problems where large density contrasts occur and the physics is related to the local density itself rather than the local density gradient.

We use numerical simulations to investigate the importance of HD formation and cooling on the first generation of metal-free stars in a LCDM cosmology. We have implemented and tested non-equilibrium HD chemistry in an adaptive mesh refinement simulation code and applied it to two situations. (1) It is first applied to the formation of 10^5 - 10^6 Msun halos which form in the absence of any ionizing source. We show, in agreement with previous work, that HD cooling is of only marginal importance for most halos; however, we find that for the lowest mass halos, HD cooling can equal or surpass the H2 cooling rate. This leads to a population of stars formed in halos with effective HD cooling that are less massive by a factor of ~6 compared to halos dominated by H2 cooling. (2) In the second part of the paper, we ionize the halos in order to explore the impact of HD cooling in the presence of an ample population of free electrons. This leads to cooler temperatures (due to the electron- catalyzed production of H2) and somewhat lower resulting proto-stellar mass. Adding HD chemistry lowers the temperature further, to the level of the CMB. We find that HD cooling dominates over H2 cooling in the density range 10^2 cm^-3 to 10^6 cm^-3, but above this density, the temperature rises and H2 cooling dominates again. Because of this, the accretion rate on to the protostar is almost the same as in the H2 case (at least for accreted masses below 50-100 Msun), therefore we argue that HD cooling in ionized halos will probably not result in a population of significantly lower mass stars.

Cluster galaxies moving through the intracluster medium (ICM) are expected to lose some of their interstellar medium (ISM) through ram pressure stripping and related ISM-ICM interactions. Using high-resolution cosmological simulations of a large galaxy cluster including star formation, we show that the ram pressure a galaxy experiences at a fixed distance from the cluster center can vary by well over an order of magnitude We find that this variation in ram pressure is due in almost equal parts to variation in the ICM density and in the relative velocity between the galaxy and the ICM. We also find that the ICM and galaxy velocities are weakly correlated for in-falling galaxies.

Until recently, simulations that modeled entire galaxies were restricted to an isothermal or fixed 2- or 3-phase interstellar medium (ISM). This obscured the full role of the ISM in shaping the observed galactic-scale star formation relations. In particular the Kennicutt relation suggests that star formation rates depend in a simple way on global galactic quantities, such as mean gas mass surface density and dynamical time. Contrary to this, observations of nearby star-forming regions, including images from the Spitzer telescope, show that all the way down to ``local-scales'' ($\sim$ parsecs) star formation is a highly clustered process with the gas existing at a wide range of densities and multiple phases, and seemingly decoupled from the larger-scale Galaxy. Many different physical processes appear to be influencing the star formation rate, including heating and cooling of the gas, ``turbulent'' energy injection from outflows, winds and supernova blastwaves, ionization, and magnetic field support. Here we present results from global simulations of disk galaxy star formation that include a fully multiphase ISM and compare to simpler isothermal models. We also investigate the effect of background heating. We discuss the process of giant molecular cloud formation as one of the possible links between the galactic and local scales of star formation.

Observational evidence for the radial alignment of satellites with their dark matter host has been accumulating steadily in the past few years. The effect is seen over a wide range of scales, from massive clusters of galaxies down to galaxy-sized systems, yet the underlying physical mechanism has still not been established. To this end, we have carried out a detailed analysis of the shapes and orientations of dark matter substructures in high-resolution N-body cosmological simulations. We find a strong tendency for radial alignment of the substructure with its host halo: the distribution of halo major axes is very anisotropic, with the majority pointing towards the center of mass of the host. The alignment peaks once the sub-halo has passed the virial radius of the host for the first time, but is not subsequently diluted, even after the halos have gone through as many as four pericentric passages. This evidence points to the existence of a very rapid dynamical mechanism acting on these systems and we argue that tidal torquing throughout their orbits is the most likely candidate.

Modelling global disc galaxies is a difficult task which has previously resulted in the small scale physics of the interstellar medium being greatly simplified. In this talk, I compare simulations of galaxies with different ISM properties to determine the importance of the ISM structure in the star formation properties of the disc.

We studied the effect of interstellar gas conditions on global galaxy simulations by considering three different models for the ISM. Our first model included only radiative cooling down to 300 K, our second model added an additional background heating term due to photoelectric heating, and our third model uses an isothermal equation of state with a temperature of 10^4 K and no explicit heating or cooling. Two common prescriptions for star formation are implemented in each case. The first is based on cosmological simulations with a low threshold for star formation but also a low efficiency. The second assumes stars form only in high density regions but with a higher efficiency. We also explore the effects of including feedback from type II supernovae. We find that the different ISM types produce marked differences in the structure of the disk and temperature phases present in the gas, although inclusion of feedback largely dominates these effects. In particular, size of the star-forming clumps was increased both by background heating and by enforcing an isothermal ISM. We also looked at the one dimensional profiles and found that a lognormal PDF provides a good fit for all our simulations over several orders of magnitude in density. Overall, despite noticeable structural differences, the star formation properties in the disk are largely insensitive to ISM type and agree reasonably well with observations.

Galaxies in clusters are gas-deficient and a number of possible explanations for this observation have been advanced, including galaxy-cluster tidal interactions, galaxy harassment, and ISM-ICM gas stripping. In this paper, we use a cosmological simulation of cluster formation and evolution in order to examine this issue from a theoretical standpoint. We follow a large number of galaxies over time and track each galaxy's gas and stellar mass changes to discover what mechanism(s) dominate the evolution of the cluster galaxies. We find that while gas is lost due to a wide variety of mechanisms, the most common way is via a gas-only stripping event, and the amount of gas lost correlates with the ram-pressure the galaxy is experiencing. Although this gas-stripping occurs primarily in the central region (r < 1 Mpc), it is an important mechanism out to the virial radius of the cluster. This is due to the wide scatter in ram-pressure strength that a galaxy experiences at fixed radius. We find that the timescale for complete gas removal is > 1 Gyr. In addition, we find that galaxies in the field and in the cluster periphery (r > 2.4 Mpc) often accrete cool gas; the accretion stops between 1-2.4 Mpc, possibly indicating the onset of galaxy starvation.

There is a strong connection between the formation of a disk galaxy and the properties of the interstellar medium (ISM). Theoretical work has typically either focused on the cosmological buildup of a galaxy with a relatively crude model for the gas physics, or examined local processes in the ISM and ignored the global evolution of the galaxy itself. Here, I briefly review what has been learned from both of these approaches, and what can be done to bridge the gap between them. I argue that cosmological simulations need to learn from observational and theoretical work on local ISM properties and adopt more sophisticated models for the processes that they cannot resolve. Since the ISM is still incompletely understood, there are a number of reasonable approaches for these "subgrid" models, and I will discuss the strengths and limitations of each.

We perform a set of non--radiative cosmological simulations of a preheated intracluster medium in which the entropy of the gas was uniformly boosted at high redshift. The results of these simulations are used first to test the current analytic techniques of preheating via entropy input in the smooth accretion limit. When the unmodified profile is taken directly from simulations, we find that this model is in excellent agreement with the results of our simulations. This suggests that preheated efficiently smoothes the accreted gas, and therefore a shift in the unmodified profile is a good approximation even with a realistic accretion history. When we examine the simulation results in detail, we do not find strong evidence for entropy amplification, at least for the high-redshift preheating model adopted here. In the second section of the paper, we compare the results of the preheating simulations to recent observations. We show -- in agreement with previous work -- that for a reasonable amount of preheating, a satisfactory match can be found to the mass-temperature and luminosity-temperature relations. However -- as noted by previous authors -- we find that the entropy profiles of the simulated groups are much too flat compared to observations. In particular, while rich clusters converge on the adiabatic self--similar scaling at large radius, no single value of the entropy input during preheating can simultaneously reproduce both the core and outer entropy levels. As a result, we confirm that the simple preheating scenario for galaxy cluster formation, in which entropy is injected universally at high redshift, is inconsistent with observations.

In this paper we describe our massively parallel version of Enzo, a multiphysics, parallel, AMR application for simulating cosmological structure formation developed at UCSD and Columbia. We describe its physics, numerical algorithms, implementation, and performance on current terascale platforms. We also discuss our future plans and some of the challenges we face as we move to the petascale.

Forthcoming large galaxy cluster surveys will yield tight constraints on cosmological models. It has been shown that in an idealized survey, containing > 10,000 clusters, statistical errors on dark energy and other cosmological parameters will be at the percent level. It has also been shown that through "self-calibration", parameters describing the mass-observable relation and cosmology can be simultaneously determined, though at a loss in accuracy by about an order of magnitude. Here we examine the utility of an alternative approach of self-calibration, in which a parametrized ab-initio physical model is used to compute cluster structure and the resulting mass-observable relations. As an example, we use a modified-entropy ("pre-heating") model of the intracluster medium, with the history and magnitude of entropy injection as unknown input parameters. Using a Fisher matrix approach, we evaluate the expected simultaneous statistical errors on cosmological and cluster model parameters. We study two types of surveys, in which a comparable number of clusters are identified either through their X-ray emission or through their integrated Sunyaev-Zel'dovich (SZ) effect. We find that compared to a phenomenological parametrization of the mass-observable relation, using our physical model yields significantly tighter constraints in both surveys, and offers substantially improved synergy when the two surveys are combined. These results suggest that parametrized physical models of cluster structure will be useful when extracting cosmological constraints from SZ and X-ray cluster surveys. (abridged)

This paper presents MoLUSC, a new method for generating mock galaxy catalogs from a large scale ($\approx 1000^3$ Mpc$^3$) dark matter simulation, that requires only modest CPU time and memory allocation. The method uses a small-scale ($\approx 256^3$ Mpc$^3$) dark matter simulation on which the \galics semi-analytic code has been run in order to define the transformation from dark matter density to galaxy density transformation using a probabilistic treatment. MoLUSC is then applied to a large-scale dark matter simulation in order to produce a realistic distribution of galaxies and their associated spectra. This permits the fast generation of large-scale mock surveys using relatively low-resolution simulations. We describe various tests which have been conducted to validate the method, and demonstrate a first application to generate a mock Sloan Digital Sky Survey redshift survey.

The concordance model of cosmology and structure formation predicts the formation of isolated very massive stars at high redshifts in dark matter dominated halos of 10^5 to 10^6 Msun. These stars photo-ionize their host primordial molecular clouds, expelling all the baryons from their halos. When the stars die, a relic HII region is formed within which large amounts of molecular hydrogen form which will allow the gas to cool efficiently when gravity assembles it into larger dark matter halos. The filaments surrounding the first star hosting halo are largely shielded and provide the pathway for gas to stream into the halo when the star has died. We present the first fully three dimensional cosmological radiation hydrodynamical simulations that follow all these effects. A novel adaptive ray casting technique incorporates the time dependent radiative transfer around point sources. This approach is fast enough so that radiation transport, kinetic rate equations, and hydrodynamics are solved self-consistently. It retains the time derivative of the transfer equation and is explicitly photon conserving. This method is integrated with the cosmological adaptive mesh refinement code enzo, and runs on distributed and shared memory parallel architectures. Where applicable the three dimensional calculation not only confirm expectations from earlier one dimensional results but also illustrate the multi-fold hydrodynamic complexities of HII regions. In the absence of stellar winds the circumstellar environments of the first supernovae and putative early gamma-ray bursts will be of low density ~1 cm^-3. Albeit marginally resolved, ionization front instabilities lead to cometary and elephant trunk like small scale structures reminiscent of nearby star forming regions.

We perform cosmological simulations of galaxies forming at z=3 using the hydrodynamics grid code, Enzo. By selecting the largest galaxies in the volume to correspond to Lyman-break galaxies, we construct observational spectra of the HI flux distribution around these objects, as well as column densities of CIV and OVI throughout a refined region. We successfully reproduce the most recent observations of the mean HI flux in the close vicinity of Lyman-break galaxies but see no evidence for the proximity effect in earlier observations. While our galaxies do return metals to the IGM, their quantity and volume appears to be somewhat less than observed. We conclude that either we do not adequately resolve galactic winds, or that at least some of the intergalactic metal enrichment is by early epoch objects whose mass is smaller than our minimum resolved halo mass.

We use three-dimensional hydrodynamic simulations to investigate the effects of a transient photoionizing ultraviolet (UV) flux on the collapse and cooling of pregalactic clouds. These clouds have masses in the range 10^5 -10^7 M_sun, form at high redshifts (z>18), are assumed to lie within the short-lived cosmological HII regions around the first generation of stars. In addition, we study the combined effects of this transient UV flux and a persistent Lyman-Werner (LW) background from distant sources. In the absence of a LW background, we find that a critical specific intensity of J_UV ~ 0.1 x 10^-21 ergs s^-1 cm^-2 Hz^-1 sr^-1 demarcates a transition from net negative to positive feedback for the halo population. A weaker UV flux stimulates subsequent star formation inside the fossil HII regions, by enhancing the H_2 molecule abundance. A stronger UV flux significantly delays star-formation by reducing the gas density, and increasing the cooling time, at the centers of collapsing halos. At a fixed J_UV, the sign of the feedback also depends strongly on the density of the gas at the time of UV illumination. Regardless of the whether the feedback is positive or negative, we find that once the UV flux is turned off, its impact stars to diminish after ~30% of the Hubble time. In the more realistic case when a LW background is present, with J_LW > 0.01 x 10^-21 ergs s^-1 cm^-2 Hz^-1 sr^-1, strong suppression persists down to the lowest redshift (z=18) in our simulations. Finally, we find evidence that heating and photoevaporation by the transient UV flux renders the ~10^6 M_sun halos inside fossil HII regions more vulnerable to subsequent H_2 photo-dissociation by a LW background.

The primordial gas in the earliest dark matter halos, collapsing at redshifts around z=20, with masses M_halo=10^6 M_sun, and virial temperatures T_vir<10^4K, relied on the presence of molecules for cooling. Several theoretical studies have suggested that gas contraction and star-formation in these minihalos was suppressed by radiative, chemical, thermal, and dynamical feedback processes. The recent measurement by the Wilkinson Microwave Anisotropy Probe (WMAP) of the optical depth to electron scattering, tau=0.09+/-0.03, provides the first empirical evidence for this suppression. The new WMAP result is consistent with vanilla models of reionization, in which ionizing sources populate cold dark matter (CDM) halos down to a virial temperature of T_vir=10^4K. On the other hand, we show that in order to avoid overproducing the optical depth, the efficiency for the production of ionizing photons in minihalos must have been about an order of magnitude lower than expected and lower than the efficiency in large halos that can cool via atomic hydrogen (T_vir > 10^4K). This conclusion is insensitive to assumptions about the efficiency of ionizing photon production in the large halos, as long as reionization ends by z=6, as required by the spectra of bright quasars at z<6. Our conclusion is strengthened if the clumping of the ionized gas evolves with redshift, as suggested by semi-analytical predictions and three-dimensional numerical simulations.

Recent far-ultraviolet (FUV) absorption line measurements of low-redshift quasars have unveiled a population of intervening broad HI Lya absorbers (BLAs) with large Doppler parameters (b> 40 km/s). If the large width of these lines is dominated by thermal line broadening, the BLAs may trace highly-ionized gas in the warm-hot intergalactic medium (WHIM) in the temperature range T ~ 10^5-10^6 K, a gas phase that is expected to contain a large fraction of the baryons at low redshift. In this paper we use a hydrodynamical simulation to study frequency, distribution, physical conditions, and baryon content of the BLAs at z=0. From our simulated spectra we derive a number of BLAs per unit redshift of (dN/dz)_BLA ~ 38 for HI absorbers with log (N(cm^-2)/b(km/s))>10.7, b>40 km/s, and log N(HII)<20.5. The baryon content of these systems is Omega_b(BLA)=0.0121/h_65, which represents ~25 percent of the total baryon budget in our simulation. Our results thus support the idea that BLAs represent a significant baryon reservoir at low redshift. BLAs predominantly trace shock-heated collisionally ionized WHIM gas at temperatures log T~4.4-6.2. About 27 percent of the BLAs in our simulation originate in the photoionized Lya forest (log T<4.3) and their large line widths are determined by non-thermal broadening effects such as unresolved velocity structure and macroscopic turbulence. Our simulation implies that for a large-enough sample of BLAs in FUV spectra it is possible to obtain a reasonable approximation of the baryon content of these systems solely from the measured HI column densities and b values.

We use a high-resolution grid-based hydrodynamics method to simulate the multi-phase interstellar medium in a Milky Way-size quiescent disk galaxy. The models are global and three-dimensional, and include a treatment of star formation and feedback. We examine the formation of gravitational instabilities and show that a form of the Toomre instability criterion can successfully predict where star formation will occur. Two common prescriptions for star formation are investigated. The first is based on cosmological simulations and has a relatively low threshold for star formation, but also enforces a comparatively low efficiency. The second only permits star formation above a number density of 1000 cm^-3 but adopts a high efficiency. We show that both methods can reproduce the observed slope of the relationship between star formation and gas surface density (although at too high a rate for our adopted parameters). A run which includes feedback from type II supernovae is successful at driving gas out of the plane, most of which falls back onto the disk. This feedback also substantially reduces the star formation rate. Finally, we examine the density and pressure distribution of the ISM, and show that there is a rough pressure equilibrium in the disk, but with a wide range of pressures at a given location (and even wider for the case including feedback

The radial entropy profile of the hot gas in clusters of galaxies tends to follow a power law in radius outside of the cluster core. Here we present a simple formula giving both the normalization and slope for the power-law entropy profiles of clusters that form in the absence of non-gravitational processes such as radiative cooling and subsequent feedback. It is based on seventy-one clusters drawn from four separate cosmological simulations, two using smoothed-particle hydrodynamics (SPH) and two using adaptive-mesh refinement (AMR), and can be used as a baseline for assessing the impact of non-gravitational processes on the intracluster medium outside of cluster cores. All the simulations produce clusters with self-similar structure in which the normalization of the entropy profile scales linearly with cluster temperature, and these profiles are in excellent agreement outside of 0.2 r_200. Because the observed entropy profiles of clusters do not scale linearly with temperature, our models confirm that non-gravitational processes are necessary to break the self-similarity seen in the simulations. However, the core entropy levels found by the two codes used here significantly differ, with the AMR code producing nearly twice as much entropy at the centre of a cluster.

The amount and nature of dark energy (DE) can be tightly constrained by measuring the spatial correlation features and evolution of a sample of ~ 100,000 galaxy clusters over the redshift range 0<z < 1.5. Such an X-ray survey will discover all collapsed structures with mass above 3.5e14 solar masss at redshifts z<2 (i.e. the full range where such objects are expected) in the high Galactic latitude sky. Above this mass threshold the tight correlations between X-ray observables and mass allow direct interpretation of the data. We describe the constraints on Dark Energy that can be inferred from such a survey, using powerful self-calibration techniques to relate the X-ray observables (luminosity and temperature) to the underlying mass.

There are (at least) two unsolved problems concerning the current state of the thermal gas in clusters of galaxies. The first is identifying the source of the heating which offsets cooling in the centers of clusters with short cooling times (the ``cooling flow'' problem). The second is understanding the mechanism which boosts the entropy in cluster and group gas. Since both of these problems involve an unknown source of heating it is tempting to identify them with the same process, particular since AGN heating is observed to be operating at some level in a sample of well-observed ``cooling flow'' clusters. Here we show, using numerical simulations of cluster formation, that much of the gas ending up in clusters cools at high redshift and so the heating is also needed at high-redshift, well before the cluster forms. This indicates that the same process operating to solve the cooling flow problem may not also resolve the cluster entropy problem.

As a first step to a more complete understanding of the local physical processes which determine star formation rates (SFRs) in the interstellar medium (ISM), we have performed controlled numerical experiments consisting of hydrodynamical simulations of a kilo-parsec scale, periodic, highly supersonic and "turbulent" three-dimensional flow. Using simple but physically motivated recipes for identifying star forming regions, we convert gas into stars which we follow self-consistently as they impact their surroundings through supernovae and stellar winds. We investigate how various processes (turbulence, radiative cooling, self-gravity, and supernovae feedback) structure the ISM, determine its energetics, and consequently affect its SFR. We find that the one-point statistical measurement captured by the probability density function (PDF) is sensitive to the simulated physics. Aside from its effect on the density structure of the medium, including self-gravity and/or supernovae feedback changes the dynamics of the medium by halting the decay of the kinetic energy. Since we find that the SFR depends most strongly on the underlying velocity field, the SFR declines in the runs lacking a means to sustain the kinetic energy, and the subsequent high density constrasts. This strong dependence on the gas velocity dispersion is in agreement with Silk's formula for the SFR (Silk 2001) which also takes the hot gas porosity, and the average gas density as important parameters. Measuring the porosity of the hot gas for the runs with supernovae feedback, we compare Silk's model for the SFR to our measured SFR and find agreement to better than a factor two.

In this paper we introduce Enzo, a 3D MPI-parallel Eulerian block-structured adaptive mesh refinement cosmology code. Enzo is designed to simulate cosmological structure formation, but can also be used to simulate a wide range of astrophysical situations. Enzo solves dark matter N-body dynamics using the particle-mesh technique. The Poisson equation is solved using a combination of fast fourier transform (on a periodic root grid) and multigrid techniques (on non-periodic subgrids). Euler's equations of hydrodynamics are solved using a modified version of the piecewise parabolic method. Several additional physics packages are implemented in the code, including several varieties of radiative cooling, a metagalactic ultraviolet background, and prescriptions for star formation and feedback. We also show results illustrating properties of the adaptive mesh portion of the code. Enzo is publicly available and can be downloaded at http://cosmos.ucsd.edu/enzo/ .

Active galactic nuclei are clearly heating gas in `cooling flows'. The effectiveness and spatial distribution of the heating are controversial. We use three-dimensional simulations on adaptive grids to study the impact on a cooling flow of weak, subrelativistic jets. The simulations show cavities and vortex rings as in the observations. The cavities are fast-expanding dynamical objects rather than buoyant bubbles as previously modelled, but shocks still remain extremely hard to detect with X-rays. At late times the cavities turn into overdensities that strongly excite the cluster's g-modes. These modes damp on a long timescale. Radial mixing is shown to be an important phenomenon, but the jets weaken the metallicity gradient only very near the centre. The central entropy density is modestly increased by the jets. We use a novel algorithm to impose the jets on the simulations.

Cooling flows, cluster mergers and motions of galaxies through the cluster gas with supersonic and sonic velocities must lead to large scale motions of the intracluster medium (ICM). A high-resolution numerical simulation of X-ray cluster formation by Norman and Bryan predicts cluster-wide turbulence with v_turb ~ 300-600 ~km/s and eddy scales l_outer ~ 100-500 kpc, the larger numbers being characteristic of turbulence near the virial radius, while the smaller numbers pertain to the core. The simulation also predicts the existence of ordered bulk flows in the core with v ~ 400 km/s on scales of several hundred kpc. In this paper we consider the observability of such fluid motions via the distortions it induces in the CMB via the kinematic SZ effect, as well as via Doppler broadening and shifting of metal lines in the X-ray spectrum. We estimate |\Delta T/T|_kinematic < few X 10^{-6} -- at or below current limits of detectability. However, we find that an energy resolution of a few eV is sufficient to detect several Doppler shifted components in the 6.7 keV Fe line in the core of the cluster.

We present results from recent simulations of the formation and evolution of clusters of galaxies in a LambdaCDM cosmology. These simulations contain our most physically complete input physics to date including radiative cooling, star formation that transforms rapidly cooling material into aggregate star particles and we also model the thermal feedback from resulting supernovae in the star particles. We use an adaptive mesh refinement (AMR) Eulerian hydrodynamics scheme to obtain very high spatial resolution (~ 2 kpc) in a computational volume 256 Mpc on a side with mass resolution for dark matter and star particles of ~ 10^8 M_solar. We examine in detail the appearance and evolution of the core region of our simulated clusters.

We present a new model for the creation of cool cores in rich galaxy clusters within a LambdaCDM cosmological framework using the results from high spatial dynamic range, adaptive mesh hydro/N-body simulations. It is proposed that cores of cool gas first form in subclusters and these subclusters merge to create rich clusters with cool, central X-Ray excesses. The rich cool clusters do not possess ``cooling flows'' due to the presence of bulk velocities in the intracluster medium in excess of 1000 km/sec produced by on-going accretion of gas from supercluster filaments. This new model has several attractive features including the presence of substantial core substructure within the cool cores, and it predicts the appearance of cool bullets, cool fronts, and cool filaments all of which have been recently observed with X-Ray satellites. This hierarchical formation model is also consistent with the observation that cool cores in Abell clusters occur preferentially in dense supercluster environments. On the other hand, our simulations overproduce cool cores in virtually all of our numerical clusters, the central densities are high, and physical core temperatures are often below 1 keV (in contrast to recent observations). We will discuss additional preliminary simulations to ``soften'' the cool cores involving star formation and supernova feedback.

The entropy distribution of the intracluster medium and the shape of its confining potential well completely determine the X-ray properties of a relaxed cluster of galaxies, motivating us to explore the origin of intracluster entropy and to describe how it develops in terms of some simple models. We present an analytical model for smooth accretion, including both preheating and radiative cooling, that links a cluster's entropy distribution to its mass accretion history and shows that smooth accretion overproduces the entropy observed in massive clusters by a factor ~2-3, depending on the mass accretion rate. Any inhomogeneity in the accreting gas reduces entropy production at the accretion shock; thus, smoothing of the gas accreting onto a cluster raises its entropy level. Because smooth accretion produces more entropy than hierarchical accretion, we suggest that some of the observed differences between clusters and groups may arise because preheating smooths the smaller-scale lumps of gas accreting onto groups more effectively than it smooths the larger-scale lumps accreting onto clusters. This effect may explain why entropy levels at the outskirts of groups are ~2-3 times larger than expected from self-similar scaling arguments. The details of how the density distribution of accreting gas affects the entropy distribution of a cluster are complex, and we suggest how to explore the relevant physics with numerical simulations.

Strong cluster eccentricity evolution for $z \le 0.13$ has appeared in a variety of observational data sets. We examine the evolution of eccentricity in simulated galaxy clusters using a variety of simulation methodologies, amplitude normalizations, and background cosmologies. We do not find find such evolution for $z < 0.1$ in any of our simulation ensembles. We suggest a systematic error in the form of a redshift-dependent selection effect in cluster catalogs or missing physics in cluster simulations important enough to modify the cluster morphology.

Measurements of the Hubble constant to distant galaxy clusters using the Sunyaev-Zeldovich effect are systematically low in comparison to values obtained by other means. These measurements usually assume a spherical isothermal $\beta$ model for the intracluster medium (ICM). We present the results of a statistical analysis of {\em temperature bias} in $H_0$ determinations in a sample of 27 numerically simulated X-ray clusters drawn from a $\Lambda$CDM model at z=0.5. We employ adaptive mesh refinement to provide high resolution (15.6 h$^{-1} kpc$) in cluster cores which dominate the X-ray and radio signals. Fitting synthetic X-ray and y-parameter maps to the standard isothermal beta model, we find a broad, skewed distribution in $f\equiv H_0(SZ)/H_0(actual)$ with a mean, median, and standard deviation of 0.89, 0.83 and 0.32 respectively, where $H_0(SZ)$ is the value of $H_0$ derived by using Sunyaev-Zeldovich effect method and $H_0(actual)$ is the value used in the cosmological simulation. We find that the clusters' declining temperature profiles systematically lower estimates of $H_0$ by 10% to 20%. The declining temperature profile of our adiabatic system is consistent with the result including radiative cooling and supernovae feedback (Loken et al. 2002). We thereby introduce a non-isothermal $\beta$ model as an improvement. Applying the non-isothermal $\beta$ model to the refined sample with well-fitted temperature profile, the value of f improves 9% relative to the actual value. The study of the morphology and the clumping effects conclude that these two factors combine to induce scatter in f of $\pm$ 30 %.

We present a new scenario for the formation of cool cores in rich galaxy clusters based on results from recent high spatial dynamic range, adaptive mesh Eulerian hydrodynamic simulations of large-scale structure formation. We find that cores of cool gas, material that would be identified as a classical cooling flow based on its X-ray luminosity excess and temperature profile, are built from the accretion of discrete, stable subclusters. Any ``cooling flow'' present is overwhelmed by the velocity field within the cluster - the bulk flow of gas through the cluster typically has speeds up to about 2,000 km s^-1 and significant rotation is frequently present in the cluster core. The inclusion of consistent initial cosmological conditions for the cluster within its surrounding supercluster environment is crucial when simulating the evolution of cool cores in rich galaxy clusters. This new model for the hierarchical assembly of cool gas naturally explains the high frequency of cool cores in rich galaxy clusters despite the fact that a majority of these clusters show evidence of substructure which is believed to arise from recent merger activity. Furthermore, our simulations generate complex cluster cores in concordance with recent X-ray observations of cool fronts, cool ``bullets'', and filaments in a number of galaxy clusters. Our simulations were computed with a coupled N-body, Eulerian, adaptive mesh refinement, hydrodynamics cosmology code that properly treats the effects of shocks and radiative cooling by the gas. We employ up to seven levels of refinement to attain a peak resolution of 15.6 h^-1 kpc within a volume 256 h-1 Mpc on a side and assume a standard LambdaCDM cosmology.

We develop a new method to estimate the redshift of galaxy clusters through resolved images of the Sunyaev-Zel'dovich effect (SZE). Our method is based on morphological observables which can be measured by actual and future SZE experiments. We test the method with a set of high resolution hydrodynamical simulations of galaxy clusters at different redshifts. Our method combines the observables in a principal component analysis. After calibrating the method with an independent redshift estimation for some of the clusters, we show - using a Bayesian approach - how the method can give an estimate of the redshift of the galaxy clusters. Although the error bars given by the morphological redshift estimation are large, it should be useful for future SZE surveys where thousands of clusters are expected to be detected; a first preselection of the high redshift candidates could be done using our proposed morphological redshift estimator. Although not considered in this work, our method should also be useful to give an estimate of the redshift of clusters in X-ray and optical surveys.

We present first results from three-dimensional hydrodynamic simulations of the high redshift formation of dwarf galaxies. The simulations use an Eulerian adaptive mesh refinement technique to follow the non-equilibrium chemistry of hydrogen and helium with cosmological initial conditions drawn from a popular Lambda-dominated CDM model. We include the effects of reionization using a uniform radiation field, a phenomenological description of the effect of star formation and, in a separate simulation, the effects of stellar feedback. The results highlight the effects of stellar feedback and photoionization on the baryon content and star formation of galaxies with virial temperatures of approximately 10^4K. Dwarf sized dark matter halos that assemble prior to reionization are able to form stars. Most halos of similar mass that assemble after reionization do not form stars by redshift of three. Dwarf galaxies that form stars show large variations in their gas content because of stellar feedback and photoionization effects. Baryon-to-dark matter mass ratios are found to lie below the cosmic mean as a result of stellar feedback. The supposed substructure problem of CDM is critically assessed on the basis of these results. The star formation histories modulated by radiative and stellar feedbacks are discussed. In addition, metallicities of individual objects are shown to be naturally correlated with their mass-to-light ratios as is also evident in the properties of local dwarf galaxies.

We use three dimensional hydrodynamic simulations to investigate the effects of a soft X-ray background, that could have been produced by an early generation of mini-quasars, on the subsequent cooling and collapse of high redshift pregalactic clouds. The simulations use an Eulerian adaptive mesh refinement technique with initial conditions drawn from a flat Lambda-dominated cold dark matter model cosmology to follow the nonequilibrium chemistry of nine chemical species in the presence of both a soft ultraviolet Lyman-Werner H_2 photodissociating flux and a soft X-ray background extending to 7.2 keV. Although we vary the normalization of the X-ray background by two orders of magnitude, the positive feedback effect of the X-rays on cooling and collapse of the pregalactic cloud expected due to the increased electron fraction is quite mild, only weakly affecting the mass threshold for collapse and the fraction of gas within the cloud that is able to cool. Inside most of the cloud we find that H_2 is in photodissociation equilibrium with the soft UV flux. The net buildup of the electron density needed to enhance H_2 formation occurs too slowly compared to the H_2 photodissociation and dynamical timescales within the cloud to overcome the negative impact of the soft UV photodissociating flux on cloud collapse. However, we find that even in the most extreme cases the first objects to form do rely on molecular hydrogen as coolant and stress that our results do not justify the neglect of these objects in models of galaxy formation.

We develop a new method for estimating the redshift of galaxy clusters through resolved images of the Sunyaev-Zel'dovich effect (SZE). Our method is based on morphological observables which can be measured by actual and future SZE experiments. The method is tested using a set of high resolution hydrodynamical simulations of galaxy clusters at different redshifts. The method combines the observables in a principal component analysis. We show how this can give an estimate of the redshift of the galaxy clusters. Although the uncertainty in the redshift estimation is large, the method should be useful for future SZE surveys where hundreds of clusters are expected to be detected. A first preselection of the high redshift candidates could be done using our proposed morphological redshift estimator.

We investigate the predicted present-day temperature profiles of the hot, X-ray emitting gas in galaxy clusters for two cosmological models - a current best-guess LCDM model and standard cold dark matter (SCDM). Our numerically-simulated "catalogs" of clusters are derived from high-resolution (15/h kpc) simulations which make use of a sophisticated, Eulerian-based, Adaptive Mesh-Refinement (AMR) code that faithfully captures the shocks which are essential for correctly modelling cluster temperatures. We show that the temperature structure on Mpc-scales is highly complex and non-isothermal. However, the temperature profiles of the simulated LCDM and SCDM clusters are remarkably similar and drop-off as $T +AFw-propto (1+-r/a_x)^{-+AFw-delta}$ where $a_x +AFw-sim r_{vir}/1.5$ and $+AFw-delta +AFw-sim 1.6$. This decrease is in good agreement with the observational results of Markevitch et al.(1998) but diverges, primarily in the innermost regions, from their fit which assumes a polytropic equation of state. Our result is also in good agreement with a recent sample of clusters observed by BeppoSAX though there is some indication of missing physics at small radii ($r<0.2 r_{vir}$). We discuss the interpretation of our results and make predictions for new x-ray observations that will extend to larger radii than previously possible. Finally, we show that, for $r>0.2 r_{vir}$, our universal temperature profile is consistent with our most recent simulations which include both radiative cooling and supernovae feedback.

We present a set of cluster models that link the present-day properties of clusters to the processes that govern galaxy formation. These models treat the entropy distribution of the intracluster medium as its most fundamental property. Because convection strives to establish an entropy gradient that rises with radius, the observable properties of a relaxed cluster depend entirely on its dark-matter potential and the entropy distribution of its uncondensed gas. Guided by simulations, we compute the intracluster entropy distribution that arises in the absence of radiative cooling and supernova heating by assuming that the gas-density distribution would be identical to that of the dark matter. The lowest-entropy gas would then fall below a critical entropy threshold at which the cooling time equals a Hubble time. Radiative cooling and whatever feedback is associated with it must modify the entropy of that low-entropy gas, changing the overall entropy distribution function and thereby altering the observable properties of the cluster. Using some phenomenological prescriptions for entropy modification based on the existence of this cooling threshold, we construct a remarkably realistic set of cluster models. The surface-brightness profiles, mass-temperature relation, and luminosity-temperature relation of observed clusters all naturally emerge from these models. By introducing a single adjustable parameter related to the amount of intracluster gas that can cool within a Hubble time, we can also reproduce the observed temperature gradients of clusters and the deviations of cooling-flow clusters from the standard luminosity-temperature relation.

As an entry for the 2001 Gordon Bell Award in the "special" category, we describe our 3-d, hybrid, adaptive mesh refinement (AMR) code, Enzo, designed for high-resolution, multiphysics, cosmological structure formation simulations. Our parallel implementation places no limit on the depth or complexity of the adaptive grid hierarchy, allowing us to achieve unprecedented spatial and temporal dynamic range. We report on a simulation of primordial star formation which develops over 8000 subgrids at 34 levels of refinement to achieve a local refinement of a factor of 10^12 in space and time. This allows us to resolve the properties of the first stars which form in the universe assuming standard physics and a standard cosmological model. Achieving extreme resolution requires the use of 128-bit extended precision arithmetic (EPA) to accurately specify the subgrid positions. We describe our EPA AMR implementation on the IBM SP2 Blue Horizon system at the San Diego Supercomputer Center.

We describe results from a fully self-consistent three dimensional hydrodynamical simulation of the formation of one of the first stars in the Universe. Dark matter dominated pre-galactic objects form because of gravitational instability from small initidal density perturbations. As they assemble via hierarchical merging, primordial gas cools through ro-vibrational lines of hydrogen molecules and sinks to the center of the dark matter potential well. The high redshift analog of a molecular cloud is formed. When the dense, central parts of the cold gas cloud become self-gravitating, a dense core of approximately 100 solar mass undergoes rapid contraction. At densities n>10^9 cm^-3 a one solar mass proto-stellar core becomes fully molecular due to three-body H_2 formation. Contrary to analytical expectations this process does not lead to renewed fragmentation and only one star is formed. The calculation is stopped when optical depth effects become important, leaving the final mass of the fully formed star somewhat uncertain. At this stage the protostar is acreting material very rapidly (~0.01 solar masses per year). Radiative feedback from the star will not only halt its growth but also inhibit the formation of other stars in the same pre-galactic object (at least until the first star ends its life, presumably as a supernova). We conclude that at most one massive (M >> 1 solar mass) metal free star forms per pre-galactic halo, consistent with recent abundance measurements of metal poor galactic halo stars.

Clusters of galaxies are thought to contain about ten times as much dark matter as baryonic matter. The dark component therefore dominates the gravitational potential of the cluster, and the baryons confined by this potential radiate X-rays with a luminosity that depends mainly on the gas density in the cluster's core. Predictions of the X-rays' properties based on models of cluster formation do not, however, agree with the observations. If the models ignore the condensation of cooling gas into stars and feedback from the associated supernovae, they overestimate the X-ray luminosity because the density of the core gas is too high. An early episode of uniformly distributed supernova feedback could rectify this by heating the uncondensed gas and therefore making it harder to compress into the core, but such a process seems to require an implausibly large number of supernovae. Here we show how radiative cooling of intergalactic gas and subsequent supernova heating conspire to eliminate highly compressible low-entropy gas from the intracluster medium. This brings the core entropy and X-ray luminosities of clusters into agreement with the observations, in a way that depends little on the efficiency of supernova heating in the early Universe.

Recent FUSE and STIS observations of OVI absorption at low redshift along the sightline to distant quasars have been interpreted as the signature of a warm-hot diffuse component of the intergalactic medium (IGM). In these proceedings we show that the predicted numbers of such absorbers in numerical simulations agree with the observed characteristics, lending support to this idea. We find that collisionally ionized lines tend to be stronger and wider, while photo-ionized absorbers are weaker and narrower. We also find that a comparison of the predicted distribution of line widths to that observed is marginally too low.

We present the high-resolution X-ray spectrum of the low-redshift quasar H 1821+643 and its surrounding hot cluster observed with the Chandra High Energy Transmission Grating Spectrometer (HETGS). An iron emission line attributed to the quasar at ~6.43 keV (rest frame) is clearly resolved, with an equivalent width of ~100 eV. Although we cannot rule out contributions to the line from a putative torus, the diskline model provides an acceptable fit to this iron line. We also detect a weak emission feature at ~6.9 keV (rest frame). We suggest that both lines could originate in an accretion disk comprised of a highly ionized optically thin atmosphere sitting atop a mostly neutral disk. We search for absorption features from a warm/hot component of the intergalactic medium along the ~1.5Gpc/h line of sight to the quasar. No absorption features are detected at or above the 3 sigma level while a total of six OVI intervening absorption systems have been detected with HST and FUSE. Based on the lack of OVII and OVIII absorption lines and by assuming collisionally ionization, we constrain the gas temperature of a typical OVI absorber to 10^5 < T < 10^6 K, which is consistent with the results from hydrodynamic simulations of the intergalactic medium. The zeroth order image reveals the extended emission from the surrounding cluster. We have been able to separate the moderate CCD X-ray spectrum of the surrounding cluster from the central quasar and find that this is a hot cluster with a temperature of ~10 keV and a metal abundance of ~0.3 Zo. We also independently obtain the redshift of the cluster, which is consistent with the optical results. We estimate that the cluster makes negligible contributions to the 6.9 keV iron K line flux.

Recent STIS and FUSE observations have detected OVI absorption lines at low redshift that are not clearly associated with any galactic system. In this paper, we argue that these lines are due to metal enriched hot gas in the intergalactic medium. Using numerical simulations of a cosmological-constant dominated cosmology, combined with reasonable assumptions about the metallicity distribution, we show that the number density and internal characteristics of these lines are correctly predicted. We find that the OVI is primarily produced by collisional ionization from gas at a few times 10^5 K for lines with equivalent widths larger than 40 mA, while weaker lines can also be produced by photo-ionization. The absorption occurs in diffuse gas in filaments at moderate overdensity (delta ~ 5-100).

Numerical simulations predict that a large number of baryons reside in intergalactic space at temperatures between 10^5-10^7 K. Highly-ionized metals, such as O VII and O VIII, are good tracers of this ``warm-hot intergalactic medium'', or WHIM. For collisionally-ionized gas, the ionization fraction of each ion peaks at some particular temperature (``peak temperatures''), so different ions can therefore trace the IGM at different temperatures. We performed a hydrodynamic simulation to study the metal distributions in the IGM. We then draw random lines-of-sight across the simulated region and synthesize resonance absorption line spectra in a similar way to simulating the Ly-alpha forest. By studying the distribution functions of H- and He-like O, Si and Fe in a collisionally-ionized IGM and comparing with semi-analytic results based on the Press-Schechter formalism, we find: (1) ions with higher peak temperatures (for instance, Fe XXVI) tend to concentrate around virialized halos, which can be well described by the Press-Schechter distribution, ions with lower peak temperatures are found both in small halos (such as groups of galaxies) and in filaments; (2) lower peak temperature ions are more abundant and should be easily observed; (3) peculiar velocities contribute a significant part to the broadening of the resonant absorption lines.

Numerical simulations predict that a large number of baryons reside in intergalactic space at temperatures between 10^5-10^7 K. Highly-ionized metals, such as O VII and O VIII, are good tracers of this ``warm-hot intergalactic medium'', or WHIM. For collisionally-ionized gas, the ionization fraction of each ion peaks at some particular temperature (``peak temperatures''), so different ions can therefore trace the IGM at different temperatures. We performed a hydrodynamic simulation to study the metal distributions in the IGM. By studying the distribution functions of H- and He-like O, Si and Fe in a collisionally-ionized IGM and comparing with semi-analytic results based on the Press-Schechter formalism, we find: (1) ions with higher peak temperatures (for instance, Fe XXVI) tend to concentrate around virialized halos, which can be well described by the Press-Schechter distribution, ions with lower peak temperatures are found both in small halos (such as groups of galaxies) and in filaments; (2) lower peak temperature ions are more abundant and should be easily observed; (3) peculiar velocities contribute a significant part to the broadening of the resonant absorption lines.

Numerical hydrodynamical simulations are used to predict the expected absorption properties of the Lyman-Alpha forest for a variety of Cold Dark Matter dominated cosmological scenarios: CHDM, OCDM, LCDM, SCDM, and tCDM. Synthetic spectra are constructed duplicating the resolution, signal-to-noise ratio, and wavelength coverage of several published high resolution spectra, and their statistical properties compared on the basis of the flux distribution of the spectra, the distribution of coefficients in a wavelet decomposition of the spectra, and the distributions of absorption line profile parameters. Agreement between the measured and predicted cumulative distributions is found at the few to several percent level. The best-fitting models to the flux distribution correspond to normalizations on the scale of the cosmological Jeans length of 1.3 < sigma_J < 1.7 at z=3. No single model provides an acceptable match to all the distributions. Larger median Doppler parameters are found in the measured spectra than predicted by all but the lowest normalization models (CHDM and tCDM), which provide poor fits to the flux distributions. The discrepancy in Doppler parameters is particularly large for absorption systems optically thin at the Lyman-Alpha line-centre. This may indicate the need to introduce additional energy injection throughout the IGM, as may be provided by late HeII reionization (z_HeII~3.5) or supernovae-driven winds from young galaxies. The models require a hydrogen ionization rate at redshifts 1.7<z<3.5 within a factor of 2 of that predicted from QSOs alone as the sources of the UV photoionization background, although with a slower rate of decline with redshift at z>3.5 than predicted from current QSO counts.

We observed MS 1054-0321, the highest redshift cluster of galaxies in the Einstein Medium Sensitivity Survey (EMSS), with the Chandra ACIS-S detector. We find the X-ray temperature of the cluster to be 10.4^{+1.7}_{-1.5} keV, slightly lower than that inferred previously and in agreement with the observed velocity dispersion and that found from weak lensing. We confirm significant substructure in the form of two distinct clumps in the X-ray distribution. We find the virial mass of the cluster to be at least 4.5 x 10^{14} h^{-1} M_{\odot} and the existence of this cluster constrains Omega_m to be less than one. If Omega_m = 1 and assuming Gaussian initial perturbations, we find the probability of observing MS 1054 in the EMSS is ~ 7 x 10^{-4}. We also present the first high-resolution X-ray spectra of two high-redshift quasars, S5 0836+710 and PKS 2149-306, obtained with the Chandra HETGS. The primary goal of this observation is to use the high spectral resolving power of the HETGS to detect X-ray absorption produced by a hot intergalactic medium. No significant emission or absorption feature is detected in either source at +/-3sigma level. Based on the detection limits we constrain the properties of possible emitters and absorbers. We propose a method to constrain the cosmological parameters via the X-ray forest theory, but current data do not give significant constraints. We also find that combined with the constraints from the distortion of the CMB spectrum, the X-ray Gunn-Peterson test can marginally constrain a uniform, enriched IGM.

We report the first high-resolution X-ray spectra of two high-redshift quasars, S5 0836+710 and PKS 2149-306, obtained with the Chandra High Energy Transmission Grating Spectrometer (HETGS). The primary goal of this observation is to use the high spectral resolving power of the HETGS to detect X-ray absorption produced by a hot intergalactic medium. The continuum of both quasars can be fitted by absorbed power laws. The power law photon index (Gamma) of S5 0836+710 is consistent with that found in previous observations with ROSAT and ASCA, while the power law of PKS 2149-306 is harder than found in a previous observation by ASCA. Excess continuum absorption above the Galactic value is found in S5 0836+710, as evidenced in ASCA and ROSAT observations. No significant emission or absorption feature is detected in either source at 3 sigma level. Based on the detection limits we constrain the properties of possible emitters and absorbers. The upper limit of the equivalent width of Fe K emission lines could be as low as ~ 10 eV. Absorbers with a column density higher than 8e16 cm^{-2} for O VIII or 5e16 cm^{-2} for Si XIV would have been detected. We propose a method to constrain the cosmological parameters (namely Omega_0 and Omega_b) via the X-ray forest theory, but current data do not give significant constraints. However, we estimate that observing ~ 7 bright quasars should give at least one O VIII and more O VII absorption lines at 95% possibility. We also find that combined with the constraints from the distortion of the CMB spectrum, the X-ray Gunn-Peterson test can marginally constrain a uniform, enriched IGM.

Hot intergalactic gas in clusters, groups, and filaments emanates a continuous background of 0.5-2.0 keV X-rays that ought to be detectable with the new generation of X-ray observatories. Here we present selected results from a program to simulate the surface-brightness distribution of this background with an adaptive-mesh cosmological hydrodynamics code. We show that the bright end of this distribution is well approximated by combining the cluster temperature function with a beta-model for surface brightness and appropriate luminosity-temperature and core radius-luminosity relations. Our simulations verify that the X-ray background from hot gas vastly exceeds observational limits if non-gravitational processes do not modify the intergalactic entropy distribution. An entropy floor of ~100 keV cm^2, which could be established by either heating or cooling, appears necessary to reconcile the simulated background with observations. Because the X-ray background distribution is so sensitive to the effects of non-gravitational processes, it offers a way to constrain the thermal history of the intergalactic medium that is independent of the uncertainties associated with surveys of clusters and groups.

We use simulations to predict the X-ray surface brightness distribution arising from hot, cosmologically distributed diffuse gas. The distribution is computed for two bands: 0.5-2 keV and 0.1-0.4 keV, using a cosmological-constant dominated cosmology that fits many other observations. We examine a number of numerical issues such as resolution, simulation volume and pixel size and show that the predicted mean background is sensitive to resolution such that higher resolution systematically increases the mean predicted background. Although this means that we can compute only lower bounds to the predicted level, these bounds are already quite restrictive. Since the observed extra-galactic X-ray background is mostly accounted for by compact sources, the amount of the observed background attributable to diffuse gas is tightly constrained. We show that without physical processes in addition to those included in the simulations (such as radiative cooling or non-gravitational heating), both bands exceed observational limits. In order to examine the effect of non-gravitational heating we explore a simple modeling of energy injection and show that substantial amounts of heating are required (i.e. 5 keV per particle when averaged over all baryons). Finally, we also compute the distribution of surface brightness on the sky and show that it has a well-resolved characteristic shape. This shape is substantially modified by non-gravitational heating and can be used as a probe of such energy injection.

Most of the baryons in the present-day universe are thought to reside in intergalactic space at temperatures of 10^5-10^7 K. X-ray emission from these baryons contributes a modest (~10%) fraction of the ~ 1 keV background whose prominence within the large-scale cosmic web depends on the amount of non-gravitational energy injected into intergalactic space by supernovae and AGNs. Here we show that the virialized regions of groups and clusters cover over a third of the sky, creating a source-confusion problem that may hinder X-ray searches for individual intercluster filaments and contaminate observations of distant groups.

Hydrodynamic simulations of Lyman alpha clouds based on ab inito cosmological models have produced results that are in broad agreement with observations. However, further analyses have revealed that, with progressively higher numerical resolution, the median or cutoff line width of the simulated Lyman alpha clouds (i.e. the Doppler parameter) appears to converge to a value significantly below what is observed at z~3 (by about a factor of 1.5). These convergence test simulations do not include feedback from star formation. Given the observed metallicity in the Lyman alpha clouds we suggest that supernovae, which presumably polluted the IGM with metals, may have deposited a sufficient amount of energy in the IGM to reconcile the theory with observations. Simple arguments immediately narrow the redshift range of pollution down to 4<z_{dep}<9. It seems quite certain that dwarf and sub-dwarf galaxies with total masses in the range 10^{6.5-9.0}Msun have to be largely responsible for the pollution. Furthermore, it is implied that either star formation is very efficient or metal yield is very high for these early dwarf galaxies, if the mean metallicity in the universe at z=3 is as high as <Z>=0.01Zsun. Finally, assuming the specific supernova heating energy is proportional to the metallicity of a gas, we note that the picture proposed here would be consistent with supernovae being the apparently needed heating source for the intra-cluster gas, if the required heating of the intra-cluster gas is no greater than 1 keV per particle.

The X-ray luminosity and temperature of clusters and groups of galaxies do not scale in a self-similar manner. This has often been interpreted as a sign that the intracluster medium has been substantially heated by non-gravitational sources. In this paper, we propose a simple model which, instead, uses the properties of galaxy formation to explain the observations. Drawing on available observations, we show that there is evidence that the efficiency of galaxy formation was higher in groups than in clusters. If confirmed, this would deplete the low-entropy gas in groups, increase their central entropy and decrease their X-ray luminosity. A simple, empirical, hydrostatic model appears to match both the luminosity-temperature relation of clusters and properties of their internal structure as well.

Approximately 30-40% of all baryons in the present day universe reside in a warm-hot intergalactic medium (WHIM), with temperatures between 10^5<T<10^7 K. This is a generic prediction from six hydrodynamic simulations of currently favored structure formation models having a wide variety of numerical methods, input physics, volumes, and spatial resolutions. Most of these warm-hot baryons reside in diffuse large-scale structures with a median overdensity around 10-30, not in virialized objects such as galaxy groups or galactic halos. The evolution of the WHIM is primarily driven by shock heating from gravitational perturbations breaking on mildly nonlinear, non-equilibrium structures such as filaments. Supernova feedback energy and radiative cooling play lesser roles in its evolution. WHIM gas is consistent with observations of the 0.25 keV X-ray background without being significantly heated by non-gravitational processes because the emitting gas is very diffuse. Our results confirm and extend previous work by Cen & Ostriker and Dave' et al.

We present results from three-dimensional hydrodynamic simulations of the high redshift collapse of pregalactic clouds including feedback effects from a soft H2 photodissociating UV radiation field. The simulations use an Eulerian adaptive mesh refinement technique to follow the nonequilibrium chemistry of nine chemical species with cosmological initial conditions drawn from a popular Lambda-dominated cold dark matter model. The results confirm that the soft UV background can delay the cooling and collapse of small halos (~10^6 Msun). For reasonable values of the photo-dissociating flux, the H2 fraction is in equilibrium throughout most of the objects we simulate. We determine the mass threshold for collapse for a range of soft-UV fluxes and also derive a simple analytic expression. Continuing the simulations beyond the point of initial collapse demonstrates that the fraction of gas which can cool depends mostly on the virial mass of the halo and the amount of soft-UV flux, with remarkably little scatter. We parameterize this relation, for use in semi-analytic models.

We discuss the results of recent 3D simulations of first structure formation in relationship to the formation of the first stars. On the basis of a new, high-resolution AMR simulation (spatial dynamic range = 30,000,000), we conclude that the first stars are likely to be massive.

We have simulated the evolution of a large sample of X-ray clusters in a Lambda-CDM universe at high spatial resolution using adaptive mesh refinement and placed the results in an online archive for public access. The Simulated Cluster Archive website http://sca.ncsa.uiuc.edu provides tools for interactive 2D and 3D analysis of gas and dark matter fields, X-ray and SZ imaging, and data export. We encourage community use and solicit their feedback.

Many questions in physical cosmology regarding the thermal history of the intergalactic medium, chemical enrichment, reionization, etc. are thought to be intimately related to the nature and evolution of pregalactic structure. In particular the efficiency of primordial star formation and the primordial IMF are of special interest. We present results from high resolution three--dimensional adaptive mesh refinement simulations that follow the collapse of primordial molecular clouds and their subsequent fragmentation within a cosmologically representative volume. Comoving scales from 128 kpc down to 1 pc are followed accurately. Dark matter dynamics, hydrodynamics and all relevant chemical and radiative processes (cooling) are followed self-consistently for a cluster normalized CDM structure formation model. Primordial molecular clouds with ~10^5 solar masses are assembled by mergers of multiple objects that have formed hydrogen molecules in the gas phase with a fractional abundance of ~10^-4. As the subclumps merge cooling lowers the temperature to ~200 Kelvin in a `cold pocket' at the center of the halo. Within this cold pocket, a quasi-hydrostatically contracting core with ~200 solar mass and number densities > 10^5 cm^-3 is found. We find that less than 1% of the primordial gas in such small scale structures cools and collapses to sufficiently high densities to be available for primordial star formation. Furthermore, it is worthwhile to note that this study achieved the highest dynamic range covered by structured adaptive mesh techniques in cosmological hydrodynamics to date.

A large fraction of the baryons today are predicted to be in hot, filamentary gas, which has yet to be detected. In this paper, we use numerical simulations of dark matter and gas to determine if these filaments and groups of galaxy will be observable be XMM. The simulated maps include free-free and line emission, galactic absorption, the XMM response function, photon noise, the extragalactic source distribution, and vignetting. A number of cosmological models are examined as well as a range of very simple prescriptions to account for the effect of supernovae feedback (preheating). We show that XMM has a good chance of observing emission from strong filaments at z~0.1. This becomes much more difficult by z~0.5. The primary difficulties lie in detecting such a large, diffuse object and in selecting an appropriate field. We also describe the range of group sizes that XMM should be sensitive to (both for detection and spectral analysis), although this is more dependent on the unknown nature of the feedback. Such observations will greatly improve our understanding of feedback and should also provide stronger cosmological constraints.

New adaptive mesh refinement N-body + hydrodynamics numerical simulations are used to illustrate the complex and changing cluster environments in which many radio galaxies live and evolve. Groups and clusters of galaxies form at the intersections of filaments where they continue to accrete gas and dark matter to the present day. The accretion process produces shocks, turbulence, and transonic bulk flows forming a kind of stormy weather within the intracluster medium (ICM). Radio sources embedded within the stormy ICM form distorted, complex morphologies as observed in recent VLA cluster surveys. We show that the bending of wide-angle tailed radio sources can be understood as the result of recent cluster-subcluster mergers. We use new MHD simulations to illustrate how cluster radio halos can be formed by the shocks and turbulence produced during cluster mergers. Finally, we discuss new observations of distant Abell clusters that reveal a class of weak radio sources, probably starbursts, likely produced during the formation of the clusters as they accrete material from the supercluster environment.

We investigate a method to determine the temperature-density relation of the intergalactic medium (IGM) at z~2-4 using quasar absorption line systems. Using a simple model combined with numerical simulations we show that there is a lower cutoff in the distribution of column density (NHI) and line width (b parameter). The location of this cutoff can be used to determine the temperature-density relation (under certain conditions). We describe and test an algorithm to do this. The method works as long as the amplitude of fluctuations on these scales (~100 kpc) is sufficiently large. Models with less power can mimic higher temperatures. A preliminary application is made to data from two quasar lines-of-sight, and we determine an upper limit to the temperature of the IGM. Finally, we examine the full distribution of b-parameters and show that this is completely specified by just two parameters: the temperature of the gas and the amplitude of the power spectrum. Using the temperature upper limit measured with the NHI-b cutoff method, we derive an upper limit to the amplitude of the power spectrum.

We investigate the properties of the Lyman alpha forest as predicted by numerical simulations for a range of currently viable cosmological models. This is done in order to understand the dependencies of the forest on cosmological parameters. Focusing on the redshift range from two to four, we show that: (1) most of the evolution in the distributions of optical depth, flux and column density can be understood by simple scaling relations, (2) the shape of optical depth distribution is a sensitive probe of the amplitude of density fluctuations on scales of a few hundred kpc, (3) the mean of the b distribution (a measure of the width of the absorption lines) is also very sensitive to fluctuations on these scales, and decreases as they increase. We perform a preliminary comparison to observations, where available. A number of other properties are also examined, including the evolution in the number of lines, the two-point flux distribution and the HeII opacity.

We have simulated the formation of an X-ray cluster in a cold dark matter universe using 12 different codes. The codes span the range of numerical techniques and implementations currently in use, including SPH and grid methods with fixed, deformable or multilevel meshes. The goal of this comparison is to assess the reliability of cosmological gas dynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at different epochs, global properties such as mass, temperature and X-ray luminosity, and radial profiles of various dynamical and thermodynamical quantities. On the whole, the agreement among the various simulations is gratifying although a number of discrepancies exist. Agreement is best for properties of the dark matter and worst for the total X-ray luminosity. Even in this case, simulations that adequately resolve the core radius of the gas distribution predict total X-ray luminosities that agree to within a factor of two. Other quantities are reproduced to much higher accuracy. For example, the temperature and gas mass fraction within the virial radius agree to about 10%, and the ratio of specific kinetic to thermal energies of the gas agree to about 5%. Various factors contribute to the spread in calculated cluster properties, including differences in the internal timing of the simulations. Based on the overall consistency of results, we discuss a number of general properties of the cluster we have modelled.

The comparison of quasar absorption spectra with numerically simulated spectra from hierarchical cosmological models of structure formation promises to be a valuable tool to discriminate among these models. We present simulation results for the column density, Doppler b parameter, and optical depth probability distributions for five popular cosmological models.

In this paper we critically examine predictions of the Lya forest within the standard cold dark matter (SCDM) model, paying particular attention to the low end of the column-density distribution. We show in particular that the width of these lines, typically measured by the b-parameter of a Voigt profile, is sensitive to spatial resolution in numerical simulations and has previously been overestimated. The new result, which predicts a distribution with a median b of around 20-22 km/s at z=3, is substantially below that observed. We examine a number of possible causes of this discrepancy and argue that it is unlikely to be rectified by an increase in the thermal broadening of the absorbing gas, but is instead telling us something about the distribution of matter on these scales. Although the median differs, the shape of the b-parameter distribution agrees quite well with that observed, and the high-end tail is naturally produced by the filamentary nature of gravitational collapse in these models. In particular, we demonstrate that lines of sight which obliquely intersect a filament or sheet tend to produce absorption lines with larger b parameters. We also examine the physical nature of the gas which is responsible for the forest, showing that for lines with neutral column densities below N_HI ~ 10^14 cm^-2 (for this model at z=3), the peculiar infall velocity is actually slower than the Hubble flow, while larger lines have, on average, turned around and are collapsing.

Many questions in physical cosmology regarding the thermal history of the intergalactic medium, chemical enrichment, reionization, etc. are thought to be intimately related to the nature and evolution of pregalactic structure. In particular the efficiency of primordial star formation and the primordial IMF are of special interest. We present results from high resolution three-dimensional adaptive mesh refinement simulations that follow the collapse of primordial molecular clouds and their subsequent fragmentation within a cosmologically representative volume. Comoving scales from 128 kpc down to 0.5 pc are followed accurately. Dark matter dynamics, hydrodynamics and all relevant chemical and radiative processes (cooling) are followed self-consistently for a cluster normalized CDM structure formation model. Primordial molecular clouds with ~100,000 solar masses are assembled by mergers of multiple objects that have formed hydrogen molecules in the gas phase with a fractional abundance of <~ 10^-4. As the subclumps merge cooling decreases the temperature to ~200 Kelvin in multiple "cold pockets" of the merger product. Within these cold pockets, quasi-hydrostatically contracting cores with masses \~ 100 M_sun and number densities > 100,000/cm^3 are found. We find that less than 1% of the primordial gas in such small scale structures cools and collapses to sufficiently high densities to be available for primordial star formation. Furthermore, our results indicate that the formation of very massive objects, massive black holes, fragmentation of a large fraction of baryons into brown dwars or Jupiter size fragments seems, in contrast to various claims in the literature, very unlikely. The expected escape fraction of UV photons with (h nu) > 11eV is very small.

In this paper we review the importance of the Lyman alpha forest as a probe of structure formation in the universe. We first discuss the statistics used to describe the Lyman alpha forest and the numerical techniques used to produce simulated spectra of the forest from a given cosmological model. We then discuss the physical picture of the absorbing structures that emerges from these numerical simulations. Finally, we comment on how two of the statistics, the slope of the column density distribution and the b parameter distribution, may be used to constrain competing cosmologies.

We describe a grid-based numerical method for 3D hydrodynamic cosmological simulations which is adaptive in space and time and combines the best features of higher order--accurate Godunov schemes for Eulerian hydrodynamics with adaptive particle--mesh methods for collisionless particles. The basis for our method is the structured adaptive mesh refinement (AMR) algorithm of Berger & Collela (1989), which we have extended to cosmological hydro + N-body simulations. The resulting multiscale hybrid method is a powerful alternative to particle-based methods in current use. The choices we have made in constructing this algorithm are discussed, and its performance on the Zeldovich pancake test problem is given. We present a sample application of our method to the problem of first structure formation. We have achieved a spatial dynamic range $L_{box}/\Delta x > 250,000$ in a 3D multispecies gas + dark matter calculation, which is sufficient to resolve the formation of primordial protostellar cloud cores starting from linear matter fluctuations in an expanding FRW universe.

We report on results of recent, high resolution hydrodynamic simulations of the formation and evolution of X-ray clusters of galaxies carried out within a cosmological framework. We employ the highly accurate piecewise parabolic method (PPM) on fixed and adaptive meshes which allow us to resolve the flow field in the intracluster gas. The excellent shock capturing and low numerical viscosity of PPM represent a substantial advance over previous studies using SPH. We find that in flat, hierarchical cosmological models, the ICM is in a turbulent state long after turbulence generated by the last major merger should have decayed away. Turbulent velocites are found to vary slowly with cluster radius, being $\sim 25%$ of $\sigma_{vir}$ in the core, increasing to $\sim 60%$ at the virial radius. We argue that more frequent minor mergers maintain the high level of turbulence found in the core where dynamical times are short. Turbulent pressure support is thus significant throughout the cluster, and results in a somewhat cooler cluster ($T/T_{vir} \sim .8$) for its mass. Some implications of cluster turbulence are discussed.

We describe an application of Berger and Colella's structured adaptive mesh refinement algorithm to the field of cosmological astrophysics. Simulations in this area must include not just a gaseous component which follows the hyperbolic equations of compressible gas dynamics, but also a collisionless component (such as dark matter or stars) described by the Newtonian dynamical equations. The two fluids interact via gravity which requires an elliptic solver. The challenge for AMR is twofold. First, the collisionless material is most easily modeled by following trajectories of individual objects, a method often referred to as an N-body scheme. The introduction of particles poses a number of difficulties, both physical (how do they interact accurately with the fluid in the mesh) and algorithmic (how to efficiently add a new data structure). The second challenge is to incorporate the gravitational interaction between the two components. We discuss our solutions to both of these issues and briefly present very encouraging results.

Gravitational instabilities naturally give rise to multi-scale structure, which is difficult for traditional Eulerian hydrodynamic methods to accurately evolve. This can be circumvented by adaptively adding resolution (in the form of multiple levels of finer meshes) to relatively small volumes as required. We describe an application of this adaptive mesh refinement (AMR) technique to cosmology, focusing on the formation and evolution of X-ray clusters. A set of simulations are performed on a single cluster, varying the initial resolution and refinement criteria. We find that although new, small scale structure continues to appear as the resolution is increased, bulk properties and radial profiles appear to converge at an effective resolution of $8192^3$.

We compare the results of Eulerian hydrodynamic simulations of cluster formation against virial scaling relations between four bulk quantities: the cluster mass, the dark matter velocity dispersion, the gas temperature and the cluster luminosity. The comparison is made for a large number of clusters at a range of redshifts in three different cosmological models (CHDM, CDM and OCDM). We find that the analytic formulae provide a good description of the relations between three of the four numerical quantities. The fourth (luminosity) also agrees once we introduce a procedure to correct for the fixed numerical resolution. We also compute the normalizations for the virial relations and compare extensively to the existing literature, finding remarkably good agreement. The Press-Schechter prescription is calibrated with the simulations, again finding results consistent with other authors. We also examine related issues such as the size of the scatter in the virial relations, the effect of metallicity with a fixed pass-band, and the structure of the halos. All of this is done in order to establish a firm groundwork for the use of clusters as cosmological probes. Implications for the models are briefly discussed.

We present the first X-ray Luminosity Function (XLF) for an optically-selected sample of 49 nearby poor clusters of galaxies and a sample of 67 Abell clusters with z < 0.15. We have extended the measured cluster XLF by more than a factor of 10 in X-ray luminosity. Our poor cluster sample was drawn from an optical catalog of groups with 0.01<z<0.03 composed of Zwicky galaxies. The X-ray emission was measured from the ROSAT all-sky survey. About 45% of the poor clusters were detected with 0.5-2.0 keV luminosities from (1.7-65)x 10^{41} h^{-2} ergs/sec. These are among the X-ray brightest, optically-selected poor clusters in the northern hemisphere. For this sample, the poor cluster XLF was found to be smooth extrapolation of the rich cluster XLF. A new Hydro/N-body simulation of a Hot +Cold dark matter model with \Omega_{total}=1, \Omega_{\nu}=0.2, and a baryon fraction of 7.5% was used to model and understand our observational selection effects. We found that the observed cluster Gas Mass Function was consistent with our model.

We describe a version of an algorithm for evolving self-gravitating collections of particles that should be nearly ideal for parallel architectures. Our method is derived from the ``self-consistent field'' (SCF) approach suggested previously by Clutton-Brock and others. Owing to the use of a global description of the gravitational field, the particles in an SCF simulation do not interact with one another directly, minimizing communications overhead between nodes in a parallel implementation. Ideal load-balancing is achieved since precisely the same number of operations are needed to compute the acceleration for each particle. Consequently, the SCF technique is perfectly scalable and the size of feasible applications will grow in simple proportion to advances in computational hardware. We describe an SCF code developed for and tested on a Connection Machine 5. Empirical tests demonstrate the efficient and scalable nature of the algorithm. Depending on the application, simulations with particle numbers in the range $N\sim 10^7 - 10^{8.5}$ are now possible. Larger platforms should make simulations with billions of particles feasible in the near future. Specific astrophysical applications are discussed in the context of collisionless dynamics.

We present the first self-consistently computed results on the evolution of X-ray properties of galaxy clusters in a Cold + Hot Dark Matter (CHDM) model. We have performed a hydrodynamic plus N-body simulation for the COBE-compatible CHDM model with standard mass components: Omega(hot) = 0.3, Omega(cold) = 0.6 and Omega(baryon) = 0.1 (h = 0.5). In contrast with the CDM model, which fails to reproduce the observed temperature distribution function dN/dT (Bryan et al. 1994b), the CHDM model fits the observational dN/dT quite well. Our results on X-ray luminosity are less firm but even more intriguing. We find that the resulting X-ray luminosity functions at redshifts z = 0.0, 0.2, 0.4, 0.7 are well fit by observations, where they overlap. The fact that both temperatures and luminosities provide a reasonable fit to the available observational data indicates that, unless we are missing some essential physics, there is neither room nor need for a large fraction of gas in rich clusters: 10% (or less) in baryons is sufficient to explain their X-ray properties. We also see a tight correlation between X-ray luminosity and gas temperature.

We present a detailed comparison of the simulation results of various cosmological hydrodynamic codes. Starting with identical initial conditions based on the Cold Dark Matter scenario for the growth of structure, we integrate from redshift $z=20$ to $z=0$ to determine the physical state within a representative volume of size $L^3$ where $L=64 h^{-1} {\rm Mpc}$. Five independent codes are compared: three of them Eulerian mesh based and two variants of the Smooth Particle Hydrodynamics "SPH" Lagrangian approach. The Eulerian codes were run at $N^3=(32^3,~64^3,~128^3,~{\rm and},~256^3)$ cells, the SPH codes at $N^3= 32^3$ and $64^3$ particles. Results were then rebinned to a $16^3$ grid with the expectation that the rebinned data should converge, by all techniques, to a common and correct result as $N \rightarrow \infty$. We find that global averages of various physical quantities do, as expected, tend to converge in the rebinned model, but that uncertainties in even primitive quantities such as $\langle T \rangle$, $\langle \rho^2\rangle^{1/2}$ persists at the 3\%-17\% level after completion of very large simulations. The two SPH codes and the two shock capturing Eulerian codes achieve comparable and satisfactory accuracy for comparable computer time in their treatment of the high density, high temperature regions as measured in the rebinned data; the variance among the five codes (at highest resolution) for the mean temperature (as weighted by $\rho^2$) is only 4.5\%. Overall the comparison allows us to better estimate errors, it points to ways of improving this current generation of hydrodynamic codes and of suiting their use to problems which exploit their individually best features.