Beecroft Institute for Particle Astrophysics and Cosmology

Abstract:

Hydrodynamical cosmological simulations are increasing their level of realism by considering more physical processes and having greater resolution or larger statistics. However, usually either the statistical power of such simulations or the resolution reached within galaxies are sacrificed. Here, we introduce the NEWHORIZON project in which we simulate at high resolution a zoom-in region of ∼(16 Mpc)3 that is larger than a standard zoom-in region around a single halo and is embedded in a larger box. A resolution of up to 34 pc, which is typical of individual zoom-in, up-to-date resimulated halos, is reached within galaxies; this allows the simulation to capture the multi-phase nature of the interstellar medium and the clumpy nature of the star formation process in galaxies. In this introductory paper, we present several key fundamental properties of galaxies and their black holes, including the galaxy mass function, cosmic star formation rate, galactic metallicities, the Kennicutt–Schmidt relation, the stellar-to-halo mass relation, galaxy sizes, stellar kinematics and morphology, gas content within galaxies and its kinematics, and the black hole mass and spin properties over time. The various scaling relations are broadly reproduced by NEWHORIZON with some differences with the standard observables. Owing to its exquisite spatial resolution, NEWHORIZON captures the inefficient process of star formation in galaxies, which evolve over time from being more turbulent, gas rich, and star bursting at high redshift. These high-redshift galaxies are also more compact, and they are more elliptical and clumpier until the level of internal gas turbulence decays enough to allow for the formation of discs. The NEWHORIZON simulation gives access to a broad range of galaxy formation and evolution physics at low-to-intermediate stellar masses, which is a regime that will become accessible in the near future through surveys such as the LSST.

Abstract:

We present the first results from SPHINX-MHD, a suite of cosmological radiation-magnetohydrodynamics simulations designed to study the impact of primordial magnetic fields (PMFs) on galaxy formation and the evolution of the intergalactic medium during the epoch of reionization. The simulations are among the first to employ on-the-fly radiation transfer and constrained transport ideal MHD in a cosmological context to simultaneously model the inhomogeneous process of reionization and the growth of PMFs. We run a series of $(5{\rm Mpc})^3$ cosmological volumes, varying both the strength of the seed magnetic field and its spectral index. We find that PMFs with a spectral index ($n_B$) and a comoving amplitude ($B_0$) that have $n_B>-0.562\log_{10}(B_0/1{\rm n}G) - 3.35$ produce electron optical depths ($\tau_e$) that are inconsistent with CMB constraints due to the unrealistically early collapse of low-mass dwarf galaxies. For $n_B\geq-2.9$, our constraints are considerably tighter than the $\sim{\rm n}G$ constraints from Planck. PMFs that do not satisfy our constraints have little impact on the reionization history or the shape of the UV luminosity function. Likewise, detecting changes in the Ly$\alpha$ forest due to PMFs will be challenging because photoionisation and photoheating efficiently smooth the density field. However, we find that the first absorption feature in the global 21cm signal is a particularly sensitive indicator of the properties of the PMFs, even for those that satisfy our $\tau_e$ constraint. Furthermore, strong PMFs can increase the escape of LyC photons by up to 25% and shrink the effective radii of galaxies by 44% which could increase the completeness fraction of galaxy surveys. Finally, our simulations show that surveys with a magnitude limit of ${\rm M_{UV,1500{\rm A}}=-13}$ can probe the sources that provide the 50% of photons for reionization out to $z=12$.

Abstract:

Due to the absence of symmetries under time and spatial translations in a general curved spacetime, the energy and momentum of matter is not conserved as it is in flat space. This means, for example, that the flux of matter energy through a surface is in general not balanced by an equal increase in the energy of the matter contained within the enclosed volume—there is an additional 'source' resulting from the curvature of spacetime acting on the matter (and vice versa). One can calculate this source term and reconcile the flux and energy accumulation over time in an arbitrary volume, although a foliation of the spacetime must be chosen, making the quantities inherently coordinate dependent. Despite this dependence, these quantities are practically useful in numerical relativity simulations for a number of reasons. We provide expressions for general matter sources in a form appropriate for implementation in the Arnowitt Deser Misner decomposition, and discuss several applications in simulations of compact object dynamics and cosmology.

Abstract:

We present the serendipitous discovery of a galaxy group in the XMM-LSS field with MIGHTEE Early Science observations. 20 galaxies are detected in H I in this z ∼ 0.044 group, with a 3σ column density sensitivity of NHI=1.6×1020cm−2⁠. This group has not been previously identified, despite residing in a well-studied extragalactic legacy field. We present spatially resolved H I total intensity and velocity maps for each of the objects which reveal environmental influence through disturbed morphologies. The group has a dynamical mass of log10(Mdyn/M⊙)=12.32⁠, and is unusually gas-rich, with an H I-to-stellar mass ratio of log10(f∗HI)=−0.2⁠, which is 0.7 dex greater than expected. The group’s high H I content, spatial, velocity, and identified galaxy type distributions strongly suggest that it is in the early stages of its assembly. The discovery of this galaxy group is an example of the importance of mapping spatially resolved H I in a wide range of environments, including galaxy groups. This scientific goal has been dramatically enhanced by the high sensitivity, large field-of-view, and wide instantaneous bandwidth of the MeerKAT telescope.

Abstract:

Utilizing optical and near-infrared broad-band photometry covering >5 deg2 in two of the most well-studied extragalactic legacy fields (COSMOS and XMM-LSS), we measure the galaxy stellar mass function (GSMF) between 0.1 < z < 2.0. We explore in detail the effect of two source extraction methods (SExtractor and ProFound) in addition to the inclusion/exclusion of Spitzer IRAC 3.6 and 4.5 μm photometry when measuring the GSMF. We find that including IRAC data reduces the number of massive (log10(M/M⊙) > 11.25) galaxies found due to improved photometric redshift accuracy, but has little effect on the more numerous lower-mass galaxies. We fit the resultant GSMFs with double Schechter functions down to log10(M/M⊙) = 7.75 (9.75) at z = 0.1 (2.0) and find that the choice of source extraction software has no significant effect on the derived best-fitting parameters. However, the choice of methodology used to correct for the Eddington bias has a larger impact on the high-mass end of the GSMF, which can partly explain the spread in derived M* values from previous studies. Using an empirical correction to model the intrinsic GSMF, we find evidence for an evolving characteristic stellar mass with δlog10(M*/M⊙)/δz = −0.16±0.05(−0.11±0.05)⁠, when using SExtractor (ProFound). We argue that with widely quenched star formation rates in massive galaxies at low redshift (z < 0.5), additional growth via mergers is required in order to sustain such an evolution to a higher characteristic mass.