Dr Harley Katz

Abstract:

To understand the impact of radiation feedback during the formation of a globular cluster (GC), we simulate a head-on collision of two turbulent giant molecular clouds (GMCs). A series of idealized radiation-hydrodynamic simulations is performed, with and without stellar radiation or Type II supernovae. We find that a gravitationally bound, compact star cluster of mass M_GC ∼ 10⁵ M_⊙ forms within ≈3 Myr when two GMCs with mass M_GMC = 3.6 × 10⁵ M⊙ collide. The GC candidate does not form during a single collapsing event but emerges due to the mergers of local dense gas clumps and gas accretion. The momentum transfer due to the absorption of the ionizing radiation is the dominant feedback process that suppresses the gas collapse, and photoionization becomes efficient once a sufficient number of stars form. The cluster mass is larger by a factor of ∼2 when the radiation feedback is neglected, and the difference is slightly more pronounced (16%) when extreme Lyα feedback is considered in the fiducial run. In the simulations with radiation feedback, supernovae explode after the star-forming clouds are dispersed, and their metal ejecta are not instantaneously recycled to form stars.

Abstract:

Our understanding of the process through which magnetic fields reached their observed strengths in present-day galaxies remains incomplete. One of the advocated solutions is a turbulent dynamo mechanism that rapidly amplifies weak magnetic field seeds to the order of ∼μG. However, simulating the turbulent dynamo is a very challenging computational task due to the demanding span of spatial scales and the complexity of the required numerical methods. In particular, turbulent velocity and magnetic fields are extremely sensitive to the spatial discretization of simulated domains. To explore how refinement schemes affect galactic turbulence and amplification of magnetic fields in cosmological simulations, we compare two refinement strategies. A traditional quasi-Lagrangian adaptive mesh refinement approach focusing spatial resolution on dense regions, and a new refinement method that resolves the entire galaxy with a high resolution quasi-uniform grid. Our new refinement strategy yields much faster magnetic energy amplification than the quasi-Lagrangian method, which is also significantly greater than the adiabatic compressional estimate indicating that the extra amplification is produced through stretching of magnetic field lines. Furthermore, with our new refinement the magnetic energy growth factor scales with resolution following ∝Δx−1/2max⁠, in much better agreement with small-scale turbulent box simulations. Finally, we find evidence suggesting most magnetic amplification in our simulated galaxies occurs in the warm phase of their interstellar medium, which has a better developed turbulent field with our new refinement strategy.