Source author record

Greg L. Bryan

Greg L. Bryan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

57works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

57 published item(s)

preprint2025arXiv

RIGEL: Feedback-regulated cloud-scale star formation efficiency in a simulated dwarf galaxy merger

Major mergers of galaxies are likely to trigger bursty star formation activities. The accumulation of dense gas and the boost of star formation efficiency (SFE) are considered to be the two main drivers of starbursts. However, it remains unclear how each process operates on the scale of individual star-forming clouds. Here, we present a high-resolution (2 Msun) RHD simulation of a gas-rich dwarf galaxy merger using the RIGEL model to investigate how mergers affect the properties of the structure of dense star-forming gas and the cloud-scale SFE. We tracked the evolution of sub-virial dense clouds in the simulation by mapping them across successive snapshots taken at intervals of 0.2 Myr. We find that the merger triggers a 130 fold increase in the SFR and shortens the galaxy-wide gas-depletion time by two orders of magnitude compared to those in two isolated galaxies. However, the depletion time of individual clouds and their lifetime distribution remained unchanged over the simulation period. The cloud life cycles and cloud-scale SFE are determined by the local stellar feedback rather than such environmental factors as tidal fields regardless of the merger process, and the integrated SFE ($ε_{\rm int}$) of clouds in complex environments remains well-described by an $ε_{\rm int}-Σ_{\rm tot}$ relation found in idealized isolated-cloud experiments. During the peak of the starburst, the media SFE was lower by only 0.17-0.33 dex compared to the value when the galaxies were not interacting. The merger boosts the SFR through the accumulation and compression of dense gas fueling star formation. Strong tidal torques assemble $>10^5$ Msun clouds, which seed massive star clusters. The average separation between star-forming clouds decreases during the merger, which in turn decreases the cloud--cluster spatial de-correlation from >1 kpc to 0.1 kpc depicted in tuning fork diagrams.

preprint2022arXiv

Formation and evolution of young massive clusters in galaxy mergers: the SMUGGLE view

Galaxy mergers are known to host abundant young massive cluster (YMC) populations, whose formation mechanism is still not well-understood. Here, we present a high-resolution galaxy merger simulation with explicit star formation and stellar feedback prescriptions to investigate how mergers affect the properties of the interstellar medium and YMCs. Compared with a controlled simulation of an isolated galaxy, the mass fraction of dense and high-pressure gas is much higher in mergers. Consequently, the mass function of both molecular clouds and YMCs becomes shallower and extends to higher masses. Moreover, cluster formation efficiency is significantly enhanced and correlates positively with the star formation rate surface density and gas pressure. We track the orbits of YMCs and investigate the time evolution of tidal fields during the course of the merger. At an early stage of the merger, the tidal field strength correlates positively with YMC mass, $λ_{\rm tid}\propto M^{0.71}$, which systematically affects the shape of the mass function and age distribution of the YMCs. At later times, most YMCs closely follow the orbits of their host galaxies, gradually sinking into the center of the merger remnant due to dynamical friction, and are quickly dissolved via efficient tidal disruption. Interestingly, YMCs formed during the first passage, mostly in tidal tails and bridges, are distributed over a wide range of galactocentric radii, greatly increasing their survivability because of the much weaker tidal field in the outskirts of the merger system. These YMCs are promising candidates for globular clusters that survive to the present day.

preprint2022arXiv

On the impact of runaway stars on dwarf galaxies with resolved interstellar medium

About ten to 20 percent of massive stars may be kicked out of their natal clusters before exploding as supernovae. These "runaway stars" might play a crucial role in driving galactic outflows and enriching the circumgalactic medium with metals. To study this effect, we carry out high resolution dwarf galaxy simulations that include velocity kicks to massive O/B stars above 8 M$_{\odot}$. We consider two scenarios, one that adopts a power law velocity distribution for kick velocities, resulting in more stars with high velocity kicks, and a more moderate scenario with a Maxwellian velocity distribution. We explicitly resolve the multi-phase interstellar medium (ISM), and include non-equilibrium cooling and chemistry channels. We adopt a resolved feedback scheme (\textsc{Griffin}) where we sample individual massive stars from an IMF. We follow the lifetime of these stars and add their photoionising radiation, their UV radiation field, and their photoelectric heating rate to the surrounding gas. At the end of their lifetime we explode the massive population as core collapse supernovae (CCSN). In the simulations with runaway massive stars, we add additional (natal) velocity kicks that mimic two and three body interactions that cannot be fully resolved in our simulations. We find that the inclusion of runaway or walkaway star scenarios has an impact on mass, metal, momentum and energy outflows as well as the respective loading factors. We find an increase in mass, metal and momentum loading by a factor of 2-3, whereas we find an increase in the mean energy loading by a factor of 5 in the runaway case and a factor of 3 in the walkaway case. However, we find that the peak values are increased by a factor of up to 10, independent of the adopted velocity kick model. We conclude that the inclusion of runaway stars could have a significant impact on the global outflow properties of dwarf galaxies.

preprint2022arXiv

Robust field-level inference with dark matter halos

We train graph neural networks on halo catalogues from Gadget N-body simulations to perform field-level likelihood-free inference of cosmological parameters. The catalogues contain $\lesssim$5,000 halos with masses $\gtrsim 10^{10}~h^{-1}M_\odot$ in a periodic volume of $(25~h^{-1}{\rm Mpc})^3$; every halo in the catalogue is characterized by several properties such as position, mass, velocity, concentration, and maximum circular velocity. Our models, built to be permutationally, translationally, and rotationally invariant, do not impose a minimum scale on which to extract information and are able to infer the values of $Ω_{\rm m}$ and $σ_8$ with a mean relative error of $\sim6\%$, when using positions plus velocities and positions plus masses, respectively. More importantly, we find that our models are very robust: they can infer the value of $Ω_{\rm m}$ and $σ_8$ when tested using halo catalogues from thousands of N-body simulations run with five different N-body codes: Abacus, CUBEP$^3$M, Enzo, PKDGrav3, and Ramses. Surprisingly, the model trained to infer $Ω_{\rm m}$ also works when tested on thousands of state-of-the-art CAMELS hydrodynamic simulations run with four different codes and subgrid physics implementations. Using halo properties such as concentration and maximum circular velocity allow our models to extract more information, at the expense of breaking the robustness of the models. This may happen because the different N-body codes are not converged on the relevant scales corresponding to these parameters.

preprint2022arXiv

The Circumgalactic Medium from the CAMELS Simulations: Forecasting Constraints on Feedback Processes from Future Sunyaev-Zeldovich Observations

The cycle of baryons through the circumgalactic medium (CGM) is important to understand in the context of galaxy formation and evolution. In this study we forecast constraints on the feedback processes heating the CGM with current and future Sunyaev-Zeldovich (SZ) observations. To constrain these processes, we use a suite of cosmological simulations, the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS), that varies four different feedback parameters of two previously existing hydrodynamical simulations, IllustrisTNG and SIMBA. We capture the dependencies of SZ radial profiles on these feedback parameters with an emulator, calculate their derivatives, and forecast future constraints on these feedback parameters from upcoming experiments. We find that for a DESI-like (Dark Energy Spectroscopic Instrument) galaxy sample observed by the Simons Observatory all four feedback parameters are able to be constrained (some within the $10\%$ level), indicating that future observations will be able to further restrict the parameter space for these sub-grid models. Given the modeled galaxy sample and forecasted errors in this work, we find that the inner SZ profiles contribute more to the constraining power than the outer profiles. Finally, we find that, despite the wide range of AGN feedback parameter variation in the CAMELS simulation suite, we cannot reproduce the tSZ signal of galaxies selected by the Baryon Oscillation Spectroscopic Survey as measured by the Atacama Cosmology Telescope.

preprint2022arXiv

The low redshift Lyman-$α$ Forest as a constraint for models of AGN feedback

We study the sensitivity of the $z=0.1$ Lyman-$α$ Forest observables, such as the column density distribution function (CDD), flux PDF, flux power spectrum, and line width distribution, to sub-grid models of active galactic nuclei (AGN) feedback using the Illustris and IllustrisTNG (TNG) cosmological simulations. The two simulations share an identical Ultraviolet Background (UVB) prescription and similar cosmological parameters, but TNG features an entirely reworked AGN feedback model. Due to changes in the AGN radio mode model, the original Illustris simulations have a factor of 2-3 fewer Lyman-$α$ absorbers than TNG at column densities $N_{\rm HI}< 10^{15.5}$ cm$^{-2}$. We compare the simulated forest statistics to UV data from the Cosmic Origins Spectrograph (COS) and find that neither simulation can reproduce the slope of the absorber distribution. Both Illustris and TNG also produce significantly smaller line width distributions than observed in the COS data. We show that TNG is in much better agreement with the observed $z=0.1$ flux power spectrum than Illustris. We explore which statistics can disentangle the effects of AGN feedback from alternative UVB models by rescaling the UVB of Illustris to produce a CDD match to TNG. While this UVB rescaling is degenerate with the effect of AGN feedback on the CDD, the amplitude and shape of the flux PDF and 1D flux power spectrum change in a way distinct from a scaling of the UVB. Our study suggests that the $z=0.1$ Lyman-$α$ forest observables can be used as a diagnostic of AGN feedback models.

preprint2021arXiv

A simple model for mixing and cooling in cloud-wind interactions

We introduce a simple entropy-based formalism to characterize the role of mixing in pressure-balanced multiphase clouds, and demonstrate example applications using Enzo-E (magneto)hydrodynamic simulations. Under this formalism, the high-dimensional description of the system's state at a given time is simplified to the joint distribution of mass over pressure ($P$) and entropy ($K=P/ρ^γ$). As a result, this approach provides a way for (empirically and analytically) quantifying the impact of different initial conditions and sets of physics on the system evolution. We find that mixing predominantly alters the distribution along the $K$ direction and illustrate how the formalism can be used to model mixing and cooling for fluid elements originating in the cloud. We further confirm and generalize a previously suggested criterion for cloud growth in the presence of radiative cooling, and demonstrate that the shape of the cooling curve, particularly at the low temperature end, can play an important role in controlling condensation. Moreover, we discuss the capacity of our approach to generalize such a criterion to apply to additional sets of physics, and to build intuition for the impact of subtle higher order effects not directly addressed by the criterion.

preprint2021arXiv

The Structure of Multiphase Galactic Winds

We present a novel analytic framework to model the steady-state structure of multiphase galactic winds comprised of a hot, volume-filling component and a cold, clumpy component. We first derive general expressions for the structure of the hot phase for arbitrary mass, momentum, and energy sources terms. Next, informed by recent simulations, we parameterize the cloud-wind mass transfer rates, which are set by the competition between turbulent mixing and radiative cooling. This enables us to cast the cloud-wind interaction as a source term for the hot phase and thereby simultaneously solve for the evolution of both phases fully accounting for their bidirectional influence. With this model, we explore the nature of galactic winds over a broad range of conditions. We find that: (i) with realistic parameter choices, we naturally produce a hot, low-density wind that transports energy while entraining a significant flux of cold clouds, (ii) mixing dominates the cold cloud acceleration and decelerates the hot wind, (iii) during mixing thermalization of relative kinetic energy provides significant heating, (iv) systems with low hot-phase mass loading factors and/or star formation rates can sustain higher initial cold phase mass loading factors, but the clouds are quickly shredded, and (v) systems with large hot-phase mass loading factors and/or star formation rates cannot sustain large initial cold-phase mass loading factors, but the clouds tend to grow with radius. Our results highlight the necessity of accounting for the multiphase structure of galactic winds, both physically and observationally, and have important implications for feedback in galactic systems.

preprint2020arXiv

First results from SMAUG: Characterization of Multiphase Galactic Outflows from a Suite of Local Star-Forming Galactic Disk Simulations

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 (Σ_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 Σ_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 Σ_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.

preprint2020arXiv

First results from SMAUG: Insights into star formation conditions from spatially-resolved ISM properties in TNG50

Physical and chemical 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 to explore the physical parameter space of 8 resolved ISM properties in star-forming regions to constrain the areas of this hyperspace over which most star-forming environments exist. We deconstruct our simulated galaxies spanning a wide range of mass (M$_\star = 10^{7-11}$ M$_\odot$) and redshift ($0 \leq z \leq 3$) into kpc-sized regions, and statistically analyze the gas/stellar surface densities, gas metallicity, vertical stellar velocity dispersion, epicyclic frequency and dark-matter volumetric density representative of each region in the context of their star formation activity and galactic environment (radial galactocentric location). By examining the star formation rate (SFR) weighted distributions of these properties, we show that stars primarily form in two spatially distinct environmental regimes, which are brought about by an underlying bi-component radial SFR surface density profile in galaxies. We examine how the relative prominence of these two regimes depends on host galaxy mass and cosmic time. We also compare our findings with those from integral field spectroscopy observations and achieve a good overall agreement. Further, using dimensionality reduction, we characterise the aforementioned hyperspace to reveal a high-degree of multicollinearity in relationships amongst ISM properties that drive the distribution of star formation at kpc-scales. Based on this, we show that a reduced 3D representation underpinned by a multi-variate radius relationship is sufficient to capture most of the variance in the original 8D space.

preprint2020arXiv

Multiphase Gas and the Fractal Nature of Radiative Turbulent Mixing Layers

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 $ρ_{\rm cold}/ρ_{\rm hot}$ of the system. Entrainment is expected to be enhanced in environments with short $t_{\rm cool}$, large $v_{\rm rel}$, and large $ρ_{\rm cold}/ρ_{\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.

preprint2020arXiv

Properties of the Simulated Circumgalactic Medium

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.

preprint2020arXiv

Self-consistent Semi-analytic Modeling of Feedback During Primordial Star Formation and Reionization

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.

preprint2020arXiv

Simple Yet Powerful: Hot Galactic Outflows Driven by Supernovae

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 $η_E$, mass $η_m$, and metals $η_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 η_{E,h}/η_{m,h}$ is almost invariant over four orders of magnitude of star formation surface density. Also, we find tentatively that $η_{E,h}/η_{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.

preprint2020arXiv

Simulating Metal Mixing of Both Common and Rare Enrichment Sources in a Low Mass Dwarf Galaxy

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.

preprint2020arXiv

The Angular Momentum of the Circumgalactic Medium in the TNG100 Simulation

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.

preprint2020arXiv

The Impact of Type Ia Supernovae in Quiescent Galaxies: I. Formation of the Multiphase Interstellar medium

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.

preprint2020arXiv

The Impact of Type Ia Supernovae in Quiescent Galaxies: II. Energetics and Turbulence

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.

preprint2020arXiv

The Role of Stellar Feedback in the Chemical Evolution of a Low Mass Dwarf Galaxy

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.

preprint2019arXiv

Direct Detection of Black Hole-Driven Turbulence in the Centers of Galaxy Clusters

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.

preprint2019arXiv

The Fate of AGB Wind in Massive Galaxies and the ICM

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.

preprint2016arXiv

Cosmological Simulations of Dwarf Galaxies with Cosmic Ray Feedback

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.

preprint2016arXiv

Formation of Massive Population III Galaxies through Photoionization Feedback: A Possible Explanation for CR7

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~Å~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$α$ emitters as bright as CR7.

preprint2016arXiv

GMC Evolution in a Barred Spiral Galaxy with Star Formation and Thermal Feedback

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.

preprint2015arXiv

Cooling, AGN Feedback and Star Formation in Simulated Cool-Core Galaxy Clusters

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.

preprint2015arXiv

Direct collapse black hole formation from synchronized pairs of atomic cooling halos

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.

preprint2015arXiv

High-Redshift Star Formation in a Time-Dependent Lyman-Werner Background

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.

preprint2015arXiv

Kinetic energy from supernova feedback in high-resolution galaxy simulations

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.

preprint2015arXiv

Limits on Population III star formation in minihaloes implied by Planck

Recently, Planck measured a value of the cosmic microwave background (CMB) optical depth due to electron scattering of $τ=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.

preprint2015arXiv

Looking for Population III stars with He II line intensity mapping

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 Å intensity mapping signal, along with the signals from other lines such as Ly$α$, H$α$, 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 Å 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 Å with CO(1-0) line emission from galaxies and the other with 21 cm emission from the intergalactic medium (IGM).

preprint2015arXiv

Precipitation-Regulated Star Formation in Galaxies

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.

preprint2015arXiv

Supernova Feedback and the Hot Gas Filling Fraction of the Interstellar Medium

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.

preprint2014arXiv

A No-Go Theorem for Direct Collapse Black Holes Without a Strong Ultraviolet Background

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.

preprint2014arXiv

An adaptive particle-mesh gravity solver for ENZO

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.

preprint2014arXiv

Cosmological simulations of galaxy formation with cosmic rays

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.

preprint2014arXiv

H_2 Suppression with Shocking Inflows: Testing a Pathway for Supermassive Black Hole Formation

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.

preprint2014arXiv

Modeling AGN Feedback in Cool-Core Clusters: The Balance between Heating and Cooling

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α$ 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.

preprint2013arXiv

Cosmic-ray Driven Outflows in Global Galaxy Disk Models

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.

preprint2013arXiv

Enzo: An Adaptive Mesh Refinement Code for Astrophysics

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.

preprint2013arXiv

Modeling AGN Feedback in Cool-Core Clusters: The Formation of Cold Clumps

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α 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.

preprint2013arXiv

The effect of feedback and reionization on star formation in low-mass dwarf galaxy haloes

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.

preprint2012arXiv

Gas Accretion is Dominated by Warm Ionized Gas in Milky Way-Mass Galaxies at z ~ 0

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.

preprint2012arXiv

Star Formation in Ram Pressure Stripped Tails

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.

preprint2011arXiv

Gas Condensation in the Galactic Halo

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.

preprint2011arXiv

How To Light It Up: Simulating Ram-Pressure Stripped X-ray Bright Tails

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)

preprint2011arXiv

Photodissociation of H2 in Protogalaxies: Modeling Self-Shielding in 3D Simulations

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.

preprint2011arXiv

Simulating the Common Envelope Phase of a Red Giant Using SPH and Uniform Grid Codes

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.

preprint2011arXiv

Simulating the Cooling Flow of Cool-Core Clusters

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.

preprint2010arXiv

Tidal Torquing of Elliptical Galaxies in Cluster Environments

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.

preprint2009arXiv

Gas Stripping in Simulated Galaxies with a Multiphase ISM

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.

preprint2009arXiv

Radiative Feedback in Relic HII Regions at High-Redshift

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]

preprint2009arXiv

The Tail of the Stripped Gas that Cooled: HI, Halpha and X-ray Observational Signatures of Ram Pressure Stripping

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).

preprint2001arXiv

Regulation of the X-ray luminosity of clusters of galaxies by cooling and supernova feedback

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.

preprint2001arXiv

The Formation of the First Star in the Universe

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.

preprint2000arXiv

Baryons in the Warm-Hot Intergalactic Medium

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.

preprint1998arXiv

Cosmological Adaptive Mesh Refinement

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}/Δ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.

preprint1995arXiv

A Parallel Self--Consistent Field Code

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.