Search

Supermassive black holes (SMBHs) of $\sim 10^9\, \rm M_{\odot }$ are generally believed to be the central engines of the luminous quasars observed at z ≳ 6, but their astrophysical origin remains elusive. To investigate the formation of these distant quasars, we perform a suite of zoom-in simulations on a favourable halo, with a mass of $\sim 10^{13}\, \rm M_{\odot }$ at z = 6 and a history of multiple major mergers, ideal for BH growth. We test BH seeds of $10 {\!-\!} 10^6\, \rm M_{\odot }$, and various accretion and feedback models, including thin-disc and slim-disc accretion. We find that abundant gas supply and efficient angular momentum transport by gravitational torques are critical to BH accretion, and that the final BH mass depends strongly on the seed mass and radiative efficiency which in turn depends on feedback model. In our simulations, heavy seeds of $\gtrsim 10^4\, \rm M_{\odot }$ can grow rapidly to SMBHs of $10^8 {\!-\!} 10^9\, \rm M_{\odot }$ by z ∼6 in both thin and slim disc accretion models, provided that the duty cycle of near-Eddington accretion with an Eddington ratio λEdd ≳ 0.6 is maintained at $\gtrsim 40{{\ \rm per\ cent}}$. In particular, for a $10^5\, \rm M_{\odot }$ seed, both our fiducial model and the supercritical accretion model with no spin produce bright quasars with SMBHs of $\sim 10^9\, \rm M_{\odot }$ at z ∼ 7.5. On the other hand, the light seeds of $\lesssim 10^3\, \rm M_{\odot }$ failed to grow to $10^8\, \rm M_{\odot }$ by z ∼ 6 in our simulations. We caution that the results are inconclusive due to limitations of our numerical methods and physical models, and we stress the need for further investigations on the growth of BH seeds with more advanced methods.

Aims. The primary difficulty in understanding the sources and processes that powered cosmic reionization is that it is not possible to directly probe the ionizing Lyman-continuum (LyC) radiation at that epoch as those photons have been absorbed by the intervening neutral hydrogen. It is therefore imperative to build a model to accurately predict LyC emission using other properties of galaxies in the reionization era. Methods. In recent years, studies have shown that the LyC emission from galaxies may be correlated to their Lyman-alpha (Lyα) emission. In this paper we study this correlation by analyzing thousands of simulated galaxies at high redshift in the SPHINX cosmological simulation. We post-process these galaxies with the Lyα radiative transfer code RASCAS and analyze the Lyα – LyC connection. Results. We find that the Lyα and LyC luminosities are strongly correlated with each other, although with dispersion. There is a positive correlation between the escape fractions of Lyα and LyC radiations in the brightest Lyman-alpha emitters (LAEs; escaping Lyα luminosity LescLyα > 1041 erg s−1), similar to that reported by recent observational studies. However, when we also include fainter LAEs, the correlation disappears, which suggests that the observed relation may be driven by selection effects. We also find that the brighter LAEs are dominant contributors to reionization, with LescLyα > 1040 erg s−1 galaxies accounting for > 90% of the total amount of LyC radiation escaping into the intergalactic medium in the simulation. Finally, we build predictive models using multivariate linear regression, where we use the physical and Lyα properties of simulated reionization era galaxies to predict their LyC emission. We build a set of models using different sets of galaxy properties as input parameters and predict their intrinsic and escaping LyC luminosity with a high degree of accuracy (the adjusted R2 of these predictions in our fiducial model are 0.89 and 0.85, respectively, where R2 is a measure of how much of the response variance is explained by the model). We find that the most important galaxy properties for predicting the escaping LyC luminosity of a galaxy are its LescLyα, gas mass, gas metallicity, and star formation rate. Conclusions. These results and the predictive models can be useful for predicting the LyC emission from galaxies using their physical and Lyα properties and can thus help us identify the sources of reionization.

ART$^2$ is a 3D multi-wavelength Monte Carlo radiative transfer (RT) code that couples continuum and emission lines to track the propagation of photons and their interactions with the interstellar medium (ISM). The previous version of ART$^2$, which included continuum and Ly$\alpha$ line, has been extensively applied to hydrodynamics simulations to study multi-band properties of galaxies and ISM. Here, we describe new implementations of non-local thermodynamic equilibrium RT of molecular and atomic fine structure emission lines, and the parallelization of the code using a number of novel methods. The new ART$^2$ can efficiently and self-consistently produce a full spectrum that includes both continuum and lines such as [CII], [NII], [OIII], Ly$\alpha$, and CO. These essential features, together with the multi-phase ISM model and the adaptive grid, make ART$^2$ a multi-purpose code to study multi-wavelength properties of a wide range of astrophysical systems from planetary disks to large-scale structures. To demonstrate the capability of the new ART$^2$, we applied it to two hydrodynamics simulations: the zoom-in Milky Way Simulation to obtain panchromatic properties of individual galaxies, and the large-scale IllustrisTNG100 Simulation to obtain global properties such as the line intensity mappings. These products are essential for a broad array of studies such as the correlations between physical and panchromatic properties and their evolution. By enabling direct comparison between numerical simulations and multi-band observations, ART$^2$ provides a crucial theoretical framework for the interpretations of existing observations, the plan for future surveys, and the synergy between multi-band galaxy surveys and line intensity mappings. Therefore, ART$^2$ is a powerful and versatile tool to bridge the gap between theories and observations of cosmic structures., Comment: Accepted for publication in MNRAS; 20 pages, 26 figures

The "Disk of satellites" (DoS) around Milky Way is a highly debated topic with conflicting interpretations of observations and their theoretical models. We perform a comprehensive analysis of all dwarfs detected in the Milky Way and find that the DoS structure depends strongly on the plane identification method and the sample size. In particular, we demonstrate that a smaller sample size produces a higher anisotropy of the spatial distribution and a stronger clustering of the angular momentum of the satellites. Moreover, we calculate the evolution of the 11 classical satellites with proper motion measurements and find that the thin DoS they currently reside in is transient. Furthermore, we analyze two cosmological simulations using the same initial conditions of a Milky Way-sized galaxy, an N-body run with dark matter only and a hydrodynamic one with both baryonic and dark matter, and find that the hydrodynamic simulation produces more anisotropic distribution of satellites than the N-body one. Our results suggest that an anisotropic distribution of satellites in galaxies can originate from baryonic processes in the hierarchical structure formation model, but the claimed highly-flattened, coherently-rotating DoS of the Milky Way may be a small-number selection effect. These findings may help resolve the contradictory claims of DoS in other galaxies and the discrepancy among numerical simulations., Comment: Published in ApJ. 7 pages, 5 figures

Observations of globular clusters show that they have universal lognormal mass functions with a characteristic peak at $\sim 2\times 10^{5}\, {\rm{M_{\odot}}}$, but the origin of this peaked distribution is highly debated. Here we investigate the formation and evolution of star clusters in interacting galaxies using high-resolution hydrodynamical simulations performed with two different codes in order to mitigate numerical artifacts. We find that massive star clusters in the range of $\sim 10^{5.5} - 10^{7.5}\, {\rm{M_{\odot}}}$ form preferentially in the highly-shocked regions produced by galaxy interactions. The nascent cluster-forming clouds have high gas pressures in the range of $P/k \sim 10^8 - 10^{12}\, \rm{K}\,\rm{cm^{-3}}$, which is $\sim 10^4 - 10^8$ times higher than the typical pressure of the interstellar medium but consistent with recent observations of a pre-super star cluster cloud in the Antennae Galaxies. Furthermore, these massive star clusters have quasi-lognormal initial mass functions with a peak around $\sim 10^{6}\, {\rm{M_{\odot}}}$. The number of clusters declines with time due to destructive processes, but the shape and the peak of the mass functions do not change significantly during the course of galaxy collisions. Our results suggest that gas-rich galaxy mergers may provide a favorable environment for the formation of massive star clusters such as globular clusters, and that the lognormal mass functions and the unique peak may originate from the extreme high-pressure conditions of the birth clouds and may survive the dynamical evolution., Comment: Accepted in ApJ. 16 pages, 9 figures

We study the impact of baryons on the distribution of dark matter in a Milky Way-size halo by comparing a high-resolution, moving-mesh cosmological simulation with its dark matter-only counterpart. We identify three main processes related to baryons -- adiabatic contraction, tidal disruption and reionization -- which jointly shape the dark matter distribution in both the main halo and its subhalos. The relative effect of each baryonic process depends strongly on the subhalo mass. For massive subhalos with maximum circular velocity $v_{\rm max} > 35 km/s$, adiabatic contraction increases the dark matter concentration, making these halos less susceptible to tidal disruption. For low-mass subhalos with $v_{\rm max} < 20 km/s$, reionization effectively reduces their mass on average by $\approx$ 30% and $v_{\rm max}$ by $\approx$ 20%. For intermediate subhalos with $20 km/s < v_{\rm max} < 35 km/s$, which share a similar mass range as the classical dwarf spheroidals, strong tidal truncation induced by the main galaxy reduces their $v_{\rm max}$. Moreover, the stellar disk of the main galaxy effectively depletes subhalos near the central region. As a combined result of reionization and increased tidal disruption, the total number of low-mass subhalos in the hydrodynamic simulation is nearly halved compared to that of the $\textit{N-}$body simulation. We do not find dark matter cores in dwarf galaxies, unlike previous studies that employed bursty feedback-driven outflows. The substantial impact of baryons on the abundance and internal structure of subhalos suggests that galaxy formation and evolution models based on $\textit{N}$-body simulations should include these physical processes as major components., 22 pages, accepted for publication in MNRAS