We examine the impact of dark matter particle resolution on the formation of a baryonic core in high resolution adaptive mesh refinement simulations. We test the effect that both particle smoothing and particle splitting have on the hydrodynamic prop
erties of a collapsing halo at high redshift (z > 20). Furthermore, we vary the background field intensity, with energy below the Lyman limit (< 13.6 eV), as may be relevant for the case of metal-free star formation and super-massive black hole seed formation. We find that using particle splitting methods greatly increases our particle resolution without introducing any numerical noise and allows us to achieve converged results over a wide range of external background fields. Additionally, we find that for lower values of the background field a lower dark matter particle mass is required. We define the radius of the core as the point at which the enclosed baryonic mass dominates over the enclosed dark matter mass. For our simulations this results in $rm{R_{core} sim 5 pc}$. We find that in order to produce converged results which are not affected by dark matter particles requires that the relationship ${M_{rm{core}} / M_{rm{DM}}} > 100.0$ be satisfied, where ${M_{rm{core}}}$ is the enclosed baryon mass within the core and $M_{rm{DM}}$ is the minimum dark matter particle mass. This ratio should provide a very useful starting point for conducting convergence tests before any production run simulations. We find that dark matter particle smoothing is a useful adjunct to already highly resolved simulations.
We explore structure formation in the dark ages ($zsim 30-6$) using two well-known methods for initializing cosmological $N$-body simulations. Overall, both the Zeldovich approximation (za) and second order Lagrangian perturbation theory (lpt) are kn
own to produce accurate present-day dark matter halo mass functions. However, since the lpt method drives more rapid evolution of dense regions, it increases the occurrence of rare massive objects -- an effect that is most pronounced at high redshift. We find that lpt produces more halos that could harbor Population III stars and their black hole remnants, and they produce them earlier. Although the differences between the lpt and za mass functions are nearly erased by $z=6$, this small boost to the number and mass of black holes more than doubles the reionized volume of the early Universe. We discuss the implications for reionization and massive black hole growth.
There is mounting evidence for the coevolution of galaxies and their embedded massive black holes (MBHs) in a hierarchical structure formation paradigm. To tackle the nonlinear processes of galaxy-MBH interaction, we describe a self-consistent numeri
cal framework which incorporates both galaxies and MBHs. The high-resolution adaptive mesh refinement (AMR) code Enzo is modified to model the formation and feedback of molecular clouds at their characteristic scale of 15.2 pc and the accretion of gas onto a MBH. Two major channels of MBH feedback, radiative feedback (X-ray photons followed through full 3D adaptive ray tracing) and mechanical feedback (bipolar jets resolved in high-resolution AMR), are employed. We investigate the coevolution of a 9.2e11 Msun galactic halo and its 1e5 Msun embedded MBH at redshift 3 in a cosmological LCDM simulation. The MBH feedback heats the surrounding ISM up to 1e6 K through photoionization and Compton heating and locally suppresses star formation in the galactic inner core. The feedback considerably changes the stellar distribution there. This new channel of feedback from a slowly growing MBH is particularly interesting because it is only locally dominant, and does not require the heating of gas globally on the disk. The MBH also self-regulates its growth by keeping the surrounding ISM hot for an extended period of time.
In hierarchical structure formation, merging of galaxies is frequent and known to dramatically affect their properties. To comprehend these interactions high-resolution simulations are indispensable because of the nonlinear coupling between pc and Mp
c scales. To this end, we present the first adaptive mesh refinement (AMR) simulation of two merging, low mass, initially gas-rich galaxies (1.8e10 Ms each), including star formation and feedback. With galaxies resolved by ~2e7 total computational elements, we achieve unprecedented resolution of the multiphase interstellar medium, finding a widespread starburst in the merging galaxies via shock-induced star formation. The high dynamic range of AMR also allows us to follow the interplay between the galaxies and their embedding medium depicting how galactic outflows and a hot metal-rich halo form. These results demonstrate that AMR provides a powerful tool in understanding interacting galaxies.
It has been argued that low-luminosity dwarf galaxies are the dominant source of ionizing radiation during cosmological reionization. The fraction of ionizing radiation that escapes into the intergalactic medium from dwarf galaxies with masses less t
han ~10^9.5 solar masses plays a critical role during this epoch. Using an extensive suite of very high resolution (0.1 pc), adaptive mesh refinement, radiation hydrodynamical simulations of idealized and cosmological dwarf galaxies, we characterize the behavior of the escape fraction in galaxies between 3 x 10^6 and 3 x 10^9 solar masses with different spin parameters, amounts of turbulence, and baryon mass fractions. For a given halo mass, escape fractions can vary up to a factor of two, depending on the initial setup of the idealized halo. In a cosmological setting, we find that the time-averaged photon escape fraction always exceeds 25% and reaches up to 80% in halos with masses above 10^8 solar masses with a top-heavy IMF. The instantaneous escape fraction can vary up to an order of magnitude in a few million years and tend to be positively correlated with star formation rate. We find that the mean of the star formation efficiency times ionizing photon escape fraction, averaged over all atomic cooling (T_vir > 8000 K) galaxies, ranges from 0.02 for a normal IMF to 0.03 for a top-heavy IMF, whereas smaller, molecular cooling galaxies in minihalos do not make a significant contribution to reionizing the universe due to a much lower star formation efficiency. These results provide the physical basis for cosmological reionization by stellar sources, predominately atomic cooling dwarf galaxies.
The first stars form in dark matter halos of masses ~10^6 M_sun as suggested by an increasing number of numerical simulations. Radiation feedback from these stars expels most of the gas from their shallow potential well of their surrounding dark matt
er halos. We use cosmological adaptive mesh refinement simulations that include self-consistent Population III star formation and feedback to examine the properties of assembling early dwarf galaxies. Accurate radiative transport is modeled with adaptive ray tracing. We include supernova explosions and follow the metal enrichment of the intergalactic medium. The calculations focus on the formation of several dwarf galaxies and their progenitors. In these halos, baryon fractions in 10^8 solar mass halos decrease by a factor of 2 with stellar feedback and by a factor of 3 with supernova explosions. We find that radiation feedback and supernova explosions increase gaseous spin parameters up to a factor of 4 and vary with time. Stellar feedback, supernova explosions, and H_2 cooling create a complex, multi-phase interstellar medium whose densities and temperatures can span up to 6 orders of magnitude at a given radius. The pair-instability supernovae of Population III stars alone enrich the halos with virial temperatures of 10^4 K to approximately 10^{-3} of solar metallicity. We find that 40% of the heavy elements resides in the intergalactic medium (IGM) at the end of our calculations. The highest metallicity gas exists in supernova remnants and very dilute regions of the IGM.
(Abridged) Using ab initio cosmological Eulerian adaptive mesh refinement radiation hydrodynamical calculations, we discuss how very massive stars start the process of cosmological reionization. The models include non-equilibrium primordial gas chemi
stry and cooling processes and accurate radiation transport in the Case B approximation using adaptively ray traced photon packages, retaining the time derivative in the transport equation. Supernova feedback is modeled by thermal explosions triggered at parsec scales. All calculations resolve the local Jeans length by at least 16 grid cells at all times and as such cover a spatial dynamic range of ~10^6. These first sources of reionization are highly intermittent and anisotropic and first photoionize the small scales voids surrounding the halos they form in, rather than the dense filaments they are embedded in. As the merging objects form larger, dwarf sized galaxies, the escape fraction of UV radiation decreases and the HII regions only break out on some sides of the galaxies making them even more anisotropic. In three cases, SN blast waves induce star formation in overdense regions that were formed earlier from ionization front instabilities. These stars form tens of parsecs away from the center of their parent DM halo. Approximately 5 ionizing photons are needed per sustained ionization when star formation in 10^6 M_sun halos are dominant in the calculation. As the halos become larger than ~10^7 M_sun, the ionizing photon escape fraction decreases, which in turn increases the number of photons per ionization to 15-50, in calculations with stellar feedback only. Supernova feedback in these more massive halos creates a more diffuse medium, allowing the stellar radiation to escape more easily and maintaining the ratio of 5 ionizing photons per sustained ionization.
Numerous cosmological hydrodynamic studies have addressed the formation of galaxies. Here we choose to study the first stages of galaxy formation, including non-equilibrium atomic primordial gas cooling, gravity and hydrodynamics. Using initial condi
tions appropriate for the concordance cosmological model of structure formation, we perform two adaptive mesh refinement simulations of ~10^8 M_sun galaxies at high redshift. The calculations resolve the Jeans length at all times with more than 16 cells and capture over 14 orders of magnitude in length scales. In both cases, the dense, 10^5 solar mass, one parsec central regions are found to contract rapidly and have turbulent Mach numbers up to 4. Despite the ever decreasing Jeans length of the isothermal gas, we only find one site of fragmentation during the collapse. However, rotational secular bar instabilities transport angular momentum outwards in the central parsec as the gas continues to collapse and lead to multiple nested unstable fragments with decreasing masses down to sub-Jupiter mass scales. Although these numerical experiments neglect star formation and feedback, they clearly highlight the physics of turbulence in gravitationally collapsing gas. The angular momentum segregation seen in our calculations plays an important role in theories that form supermassive black holes from gaseous collapse.
