Cosmological Zoom-in Simulations Research Articles

Stellar angular momentum can be controlled from cosmological initial conditions

ABSTRACT The angular momentum of galaxies controls the kinematics of their stars, which in turn drives observable quantities such as the apparent radius, the bulge fraction, and the alignment with other nearby structures. To show how angular momentum of galaxies is determined, we build high (35 pc) resolution numerical experiments in which we increase or decrease the angular momentum of the Lagrangian patches in the early universe. We perform cosmological zoom-in simulations of three galaxies over their histories from z = 200 to z = 2, each with five different choices for the angular momentum (15 simulations in total). Our results show that altering early universe angular momentum changes the timing and orbital parameters of mergers, which in turn changes the total stellar angular momentum within a galaxy’s virial radius in a predictable manner. Of our three galaxies, one has no large satellite at z = 2; in this case, the specific angular momentum is concentrated in the central galaxy. Our changes to the initial conditions result in its stellar angular momentum changing over 0.7 dex (from 61 to ${320}\, {\rm kpc\, km\, s}^{-1}$) at z = 2. This causes its effective radius to grow by 40 per cent, its v/σ parameter to grow by a factor of 2.6, and its bulge fraction to decrease from 0.72 to 0.57. This proof of concept illustrates how causal studies can contribute to a better understanding of the origin of galaxy scaling relations and intrinsic alignments.

The merger–starburst connection across cosmic times

ABSTRACT The correspondence between galaxy major mergers and starburst activity is well established observationally and in simulations of low-redshift galaxies. However, the evolution of the properties of interactions and of the galaxies involved suggests that the starburst response of galaxies to merger events could vary across cosmic time. Using the vintergatan cosmological zoom-in simulation of a Milky Way-like galaxy, we show here that starbursts, i.e. episodes of fast star formation, are connected with the onset of tidal compression, itself induced by mergers. However, this compression becomes strong enough to trigger starbursts only after the formation of the galactic disc. As a consequence, starburst episodes are only found during a precise phase of galaxy evolution, after the formation of the disc and until the last major merger. As the depletion time quantifies the instantaneous star formation activity, while the specific star formation rate involves the integrated result of the past activity (via the stellar mass), starburst episodes do not necessarily coincide with elevated specific star formation rate. This suggests that not all starburst galaxies are outliers above the main sequence of galaxy formation.ok

Sub-parsec resolution cosmological simulations of star-forming clumps at high redshift with feedback of individual stars

ABSTRACT We introduce a new set of zoom-in cosmological simulations with sub-pc resolution, intended to model extremely faint, highly magnified star-forming stellar clumps, detected at z = 6.14 thanks to gravitational lensing. The simulations include feedback from individual massive stars (in both the pre-supernova and supernova phases), generated via stochastic, direct sampling of the stellar initial mass function. We adopt a modified ‘delayed cooling’ feedback scheme, specifically created to prevent artificial radiative loss of the energy injected by individual stars in very dense gas (n ∼ 103–105 cm−3). The sites where star formation ignites are characterized by maximum densities of the order of 105 cm−3 and gravitational pressures Pgrav/k >107 K cm−3, corresponding to the values of the local, turbulent regions where the densest stellar aggregates form. The total stellar mass at z = 6.14 is 3.4$\times 10^7~\rm M_{\odot }$, in satisfactory agreement with the observed stellar mass of the observed systems. The most massive clumps have masses of $\sim 10^6~\rm M_{\odot }$ and half-mass sizes of ∼100 pc. These sizes are larger than the observed ones, including also other samples of lensed high-redshift clumps, and imply an average density one orders of magnitude lower than the observed one. In the size–mass plane, our clumps populate a sequence that is intermediate between the ones of observed high-redshift clumps and local dSph galaxies.

KODIAQ-Z: Metals and Baryons in the Cool Intergalactic and Circumgalactic Gas at 2.2 ≲ z ≲ 3.6

We present the KODIAQ-Z survey aimed to characterize the cool, photoionized gas at 2.2 ≲ z ≲ 3.6 in 202 H i-selected absorbers with 14.6 ≤ < 20 that probe the interface between galaxies and the intergalactic medium (IGM). We find that gas with at 2.2 ≲ z ≲ 3.6 can be metal-rich (−1.6 ≲ [X/H] ≲ − 0.2) as seen in damped Lyα absorbers (DLAs); it can also be very metal-poor ([X/H] < − 2.4) or even pristine ([X/H] < − 3.8), which is not observed in DLAs but is common in the IGM. For absorbers, the frequency of pristine absorbers is about 1%–10%, while for absorbers it is 10%–20%, similar to the diffuse IGM. Supersolar gas is extremely rare (<1%) at these redshifts. The factor of several thousand spread from the lowest to highest metallicities and large metallicity variations (a factor of a few to >100) between absorbers separated by less than Δv < 500 km s−1 imply that the metals are poorly mixed in gas. We show that these photoionized absorbers contribute to about 14% of the cosmic baryons and 45% of the cosmic metals at 2.2 ≲ z ≲ 3.6. We find that the mean metallicity increases with N H i , consistent with what is found in z < 1 gas. The metallicity of gas in this column density regime has increased by a factor ∼8 from 2.2 ≲ z ≲ 3.6 to z < 1, but the contribution of the absorbers to the total metal budget of the universe at z < 1 is a quarter of that at 2.2 ≲ z ≲ 3.6. We show that FOGGIE cosmological zoom-in simulations have a similar evolution of [X/H] with N H i , which is not observed in lower-resolution simulations. In these simulations, very metal-poor absorbers with [X/H] < − 2.4 at z ∼ 2–3 are tracers of inflows, while higher-metallicity absorbers are a mixture of inflows and outflows.

X-ray morphology of cluster-mass haloes in self-interacting dark matter

ABSTRACT We perform cosmological zoom-in simulations of 19 relaxed cluster-mass haloes with the inclusion of adiabatic gas in the cold dark matter (CDM) and self-interacting dark matter (SIDM) models. These clusters are selected as dynamically relaxed clusters from a parent simulation with $M_{\rm 200} \simeq (1\!-\!3)\times 10^{15}{\, \rm M_\odot }$. Both the dark matter and the intracluster gas distributions in SIDM appear more spherical than their CDM counterparts. Mock X-ray images are generated based on the simulations and are compared to the real X-ray images of 84 relaxed clusters selected from the Chandra and ROSAT archives. We perform ellipse fitting for the isophotes of mock and real X-ray images and obtain the ellipticities at cluster-centric radii of $r\simeq 0.1\!-\!0.2R_{\rm 200}$. The X-ray isophotes in SIDM models with increasing cross-sections are rounder than their CDM counterparts, which manifests as a systematic shift in the distribution function of ellipticities. Unexpectedly, the X-ray morphology of the observed non-cool-core clusters agrees better with SIDM models with cross-section $(\sigma /m)= 0.5\!-\!1\, {\rm cm}^2\, {\rm g}^{-1}$ than CDM and SIDM with $(\sigma /m)=0.1\, {\rm cm}^2\, {\rm g}^{-1}$. Our statistical analysis indicates that the latter two models are disfavoured at the $68{{\ \rm per\ cent}}$ confidence level (as conservative estimates). This conclusion is not altered by shifting the radial range of measurements or applying a temperature selection criterion. However, the primary uncertainty originates from the lack of baryonic physics in the adiabatic model, such as cooling, star formation and feedback effects, which still have the potential to reconcile CDM simulations with observations.

VINTERGATAN IV: Cosmic phases of star formation in Milky Way-like galaxies

ABSTRACT The star formation history of a galaxy is modulated by a plethora of internal processes and environmental conditions. The details of how these evolve and couple together are not fully understood yet. In this work, we study the effects that galaxy mergers and morphological transformations have on setting different modes of star formation at galactic scales and across cosmic time. We monitor the global properties of vintergatan, a 20 pc resolution cosmological zoom-in simulation of a Milky Way-type galaxy. Between redshifts 1 and 5, we find that major mergers trigger multiple starburst episodes, corresponding to a tenfold drop of the gas depletion time down to 100 Myr. Bursty star formation is enabled by the emergence of a galactic disc, when the rotational velocity of gas starts to dominate over its velocity dispersion. Coherent motions of gas then outweigh disordered ones, such that the galaxy responds to merger-induced forcings by redistributing large amounts of gas towards high densities. As a result, the overall star formation rate (SFR) is enhanced with an associated decrease in the depletion time. Before redshift 5, mergers were expected to be even more frequent. However, a more turbulent interstellar medium is incapable of reacting in such a collective manner so as to spark rapid star formation. Thus, a constant long depletion time of 1 Gyr is kept, along with a low, but gradually increasing SFR. After the last major merger at redshift 1, vintergatan spends the next 8 Gyr evolving secularly. It has a settled and adiabatically growing disc, and a constant SFR with gas depletion times of 1–2 Gyr. Our results are compatible with the observed rapid transition between different modes of star formation when galaxies leave the main sequence.

Two can play at that game: constraining the role of supernova and AGN feedback in dwarf galaxies with cosmological zoom-in simulations

ABSTRACT There is growing observational evidence for dwarf galaxies hosting active galactic nuclei (AGN), including hints of AGN-driven outflows in dwarfs. However, in the common theoretical model of galaxy formation, efficient supernova (SN) feedback is the tool of choice for regulating star formation in the low-mass regime. In this paper, we present a suite of high-resolution cosmological dwarf zoom-in simulations relaxing the assumption of strong SN feedback, with the goal to determine whether more moderate SN feedback in combination with an efficient AGN could be a suitable alternative. Importantly, we find that there are sufficient amounts of gas to power brief Eddington-limited accretion episodes in dwarfs. This leads to a variety of outcomes depending on the AGN accretion model: from no additional suppression to moderate regulation of star formation to catastrophic quenching. Efficient AGN can drive powerful outflows, depleting the gas reservoir of their hosts via ejective feedback and then maintaining a quiescent state through heating the circumgalactic medium. Moderate AGN outflows can be as efficient as the strong SN feedback commonly employed, leading to star formation regulation and H i gas masses in agreement with observations of field dwarfs. All efficient AGN set-ups are associated with overmassive black holes (BHs) compared to the (heavily extrapolated) observed BH mass–stellar mass scaling relations, with future direct observational constraints in this mass regime being crucially needed. Efficient AGN activity is mostly restricted to high redshifts, with hot, accelerated outflows and high X-ray luminosities being the clearest tell-tale signs for future observational campaigns.

Exploring metallicity-dependent rates of Type Ia supernovae and their impact on galaxy formation

ABSTRACT Type Ia supernovae are critical for feedback and elemental enrichment in galaxies. Recent surveys like the All-Sky Automated Survey for Supernova (ASAS-SN) and the Dark Energy Survey (DES) find that the specific supernova Ia rate at z ∼ 0 may be ≲20–50× higher in lower mass galaxies than at Milky Way-mass. Independently, observations show that the close-binary fraction of solar-type Milky Way stars is higher at lower metallicity. Motivated by these observations, we use the FIRE-2 cosmological zoom-in simulations to explore the impact of metallicity-dependent rate models on galaxies of $M_* \sim 10^7\!-\!10^{11}\, \rm {M}_{\odot }$. First, we benchmark our simulated star formation histories against observations, and show that the assumed stellar mass functions play a major role in determining the degree of tension between observations and metallicity-independent rate models, potentially causing ASAS-SN and DES observations to agree more than might appear. Models in which the supernova Ia rate increases with decreasing metallicity ($\propto Z^{-0.5 \,\, \rm {to} \,\, -1}$) provide significantly better agreement with observations. Encouragingly, these rate increases (≳10× in low-mass galaxies) do not significantly impact galaxy masses and morphologies, which remain largely unaffected except for our most extreme models. We explore implications for both [Fe/H] and [$\alpha /\rm {Fe}$] enrichment; metallicity-dependent rate models can improve agreement with the observed stellar mass–metallicity relations in low-mass galaxies. Our results demonstrate that a range of metallicity-dependent rate models are viable for galaxy formation and motivate future work.

Cold-mode and hot-mode accretion in galaxy formation: an entropy approach

ABSTRACT We have analysed two cosmological zoom simulations with $M_{\rm vir}\sim 10^{12}{\rm \, M}_\odot$ from the Numerical Investigation of a Hundred Astrophysical Objects (NIHAO) series, both with and without feedback. We show that an entropy criterion based on the equation of state of the intergalactic medium can successfully separate cold- and hot-mode accretion. The shock-heated gas has non-negligible turbulent support and cools inefficiently. In the simulations without feedback, only a small fraction (≲20 per cent) of the stellar mass comes from baryons that have been in the hot circumgalactic medium, although quantitative conclusions should be taken with caution due to our small-number statistics. With feedback, the fraction is larger because of the reaccretion of gas heated by supernovae, which has lower entropies and shorter cooling times than the gas heated by accretion shocks. We have compared the results of NIHAO to predictions of the GalICS 2.1 semi-analytic model of galaxy formation. The shock-stability criterion implemented in GalICS 2.1 successfully reproduces the transition from cold- to hot-mode accretion.

Modeling Evolution of Galactic Bars at Cosmic Dawn

We study the evolution of galactic bars using a suite of very high-resolution zoom-in cosmological simulations of galaxies at z ∼ 9–2. Our models were chosen to lie within similar-mass dark matter (DM) halos, log(M vir/M ⊙) ∼ 11.65 ± 0.05, at z = 6, 4, and 2, in high- and low-overdensity environments. We apply two galactic wind feedback mechanisms for each model. All galaxies develop subkiloparsec stellar bars differing in their properties. We find that (1) the high-z bars form in response to various perturbations: mergers, close flybys, cold accretion inflows along the cosmological filaments, etc.; (2) these bars account for the large mass fraction of galaxies; (3) bars display large corotation-to-bar size ratios, are weaker compared to their low-redshift counterparts by measuring their Fourier amplitudes, and are very gas-rich; (4) their pattern speed does not exhibit monotonic decline with time owing to braking against DM, as at low z; (5) bar properties, including their stellar population (star formation rates and metal enrichment), depend sensitively on prevailing feedback; and (6) bars can weaken substantially during cosmological evolution, becoming weak oval distortions—hence bars are destroyed and reformed multiple times, unlike their low-z counterparts. In all cases, bars in our simulations have been triggered by interactions. In summary, not only do stellar bars appear to be a contemporary phenomenon, but based on increased frequency of mergers, flybys, and the strength of cold accretion flows at high z, we also expect them to be ubiquitous at redshifts ≳2—the epoch of rapid galaxy growth and larger stellar dispersion velocities.

On the Small-scale Turbulent Dynamo in the Intracluster Medium: A Comparison to Dynamo Theory* * Released on August, 19th, 2021.

We present non-radiative, cosmological zoom-in simulations of galaxy-cluster formation with magnetic fields and (anisotropic) thermal conduction of one massive galaxy cluster with M vir ∼ 2 × 1015 M ⊙ at z ∼ 0. We run the cluster on three resolution levels (1×, 10×, 25×), starting with an effective mass resolution of 2 × 108 M ⊙, subsequently increasing the particle number to reach 4 × 106 M ⊙. The maximum spatial resolution obtained in the simulations is limited by the gravitational softening reaching ϵ = 1.0 kpc at the highest resolution level, allowing one to resolve the hierarchical assembly of the structures in fine detail. All simulations presented are carried out with the SPMHD code gadget3 with an updated SPMHD prescription. The primary focus of this paper is to investigate magnetic field amplification in the intracluster medium. We show that the main amplification mechanism is the small-scale turbulent dynamo in the limit of reconnection diffusion. In our two highest resolution models we start to resolve the magnetic field amplification driven by the dynamo and we explicitly quantify this with the magnetic power spectra and the curvature of the magnetic field lines, consistent with dynamo theory. Furthermore, we investigate the ∇ · B = 0 constraint within our simulations and show that we achieve comparable results to state-of-the-art AMR or moving-mesh techniques, used in codes such as enzo and arepo. Our results show for the first time in a cosmological simulation of a galaxy cluster that dynamo action can be resolved with modern numerical Lagrangian magnetohydrodynamic methods, a study that is currently missing in the literature.

Reliability and limitations of inferring birth radii in the Milky Way disc

ABSTRACT Recovering the birth radii of observed stars in the Milky Way is one of the ultimate goals of Galactic Archaeology. One method to infer the birth radius and the evolution of the interstellar medium (ISM) metallicity assumes a linear relation between the ISM metallicity with radius at any given look-back time. Here, we test the reliability of this assumption by using four zoom-in cosmological hydrodynamic simulations from the NIHAO-UHD project. We find that one can infer precise birth radii only when the stellar disc starts to form, which for our modelled galaxies happens ∼10 Gyr ago, in agreement with recent estimates for the Milky Way. With a current day measurement of ISM metallicity gradient of −0.05 dex and a dispersion of 0.03 dex, the intrinsic uncertainty in inferring Rbirth is ∼0.6 kpc. At later times, the linear correlation between the ISM metallicity and radius increases, as stellar motions become more ordered and the azimuthal variations of the ISM metallicity start to drop. The formation of a central bar and perturbations from mergers can increase this uncertainty in the inner and outer disc, respectively.

Stellar migration in the Auriga simulations

ABSTRACT We study the presence and importance of stellar migration in the evolution of 17 Milky-Way like disc galaxies with stellar mass 10 &amp;lt; log(M*/M⊙) &amp;lt; 11 from the Auriga suite of zoom-in cosmological hydrodynamical simulations. We compare the birth radii of the stars to their radii at z = 0 for each system and present mean values of the strength of stellar migration as a function of radius and stellar age which vary between 1–4 kpc. We also investigate the effect of migration on age and metallicity radial profiles in the discs. We find several cases of age gradient flattening due to migration, but significant changes to metallicity profiles only for older stellar populations and discs that develop a strong bar. Furthermore, we study stellar migration from the perspective of the change of the galactocentric radius (ΔR) and orbital guiding centre radius (ΔRg) of stellar particles between given time intervals. We find that stars migrate approximately as a diffusion process only in the outer parts of the discs and for particular galaxies that have a weak bar. Strongly barred galaxies in our sample show larger stellar migration but its timestep evolution is slower-than-diffusion. Finally, we give parametrizations that encapsulate the dependence of the strength of the radial migration as a function of time and radius, for incorporation into (semi-)analytic models of galaxy evolution.

Machine-guided exploration and calibration of astrophysical simulations

ABSTRACT We apply a novel method with machine learning to calibrate sub-grid models within numerical simulation codes to achieve convergence with observations and between different codes. It utilizes active learning and neural density estimators. The hyper parameters of the machine are calibrated with a well-defined projectile motion problem. Then, using a set of 22 cosmological zoom simulations, we tune the parameters of a popular star formation and feedback model within Enzo to match observations. The parameters that are adjusted include the star formation efficiency, coupling of thermal energy from stellar feedback, and volume into which the energy is deposited. This number translates to a factor of more than three improvements over manual calibration. Despite using fewer simulations, we obtain a better agreement to the observed baryon makeup of a Milky Way (MW)-sized halo. Switching to a different strategy, we improve the consistency of the recommended parameters from the machine. Given the success of the calibration, we then apply the technique to reconcile metal transport between grid-based and particle-based simulation codes using an isolated galaxy. It is an improvement over manual exploration while hinting at a less-known relation between the diffusion coefficient and the metal mass in the halo region. The exploration and calibration of the parameters of the sub-grid models with a machine learning approach is concluded to be versatile and directly applicable to different problems.

The In Situ Origins of Dwarf Stellar Outskirts in FIRE-2

Extended, old, and round stellar halos appear to be ubiquitous around high-mass dwarf galaxies (108.5 < M ⋆/M ⊙ < 109.6) in the observed universe. However, it is unlikely that these dwarfs have undergone a sufficient number of minor mergers to form stellar halos that are composed of predominantly accreted stars. Here, we demonstrate that FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulations are capable of producing dwarf galaxies with realistic structures, including both a thick disk and round stellar halo. Crucially, these stellar halos are formed in situ, largely via the outward migration of disk stars. However, there also exists a large population of “nondisky” dwarfs in FIRE-2 that lack a well-defined disk/halo and do not resemble the observed dwarf population. These nondisky dwarfs tend to be either more gas-poor or to have burstier recent star formation histories than the disky dwarfs, suggesting that star formation feedback may be preventing disk formation. Both classes of dwarfs underscore the power of a galaxy’s intrinsic shape—which is a direct quantification of the distribution of the galaxy’s stellar content—to interrogate the feedback implementation in simulated galaxies.

The Signatures of Self-interacting Dark Matter and Subhalo Disruption on Cluster Substructure

The abundance, distribution, and inner structure of satellites of galaxy clusters can be sensitive probes of the properties of dark matter. We run 30 cosmological zoom-in simulations with self-interacting dark matter (SIDM), with a velocity-dependent cross section, to study the properties of subhalos within cluster-mass hosts. We find that the abundance of subhalos that survive in the SIDM simulations are suppressed relative to their cold dark matter (CDM) counterparts. Once the population of disrupted subhalos—which may host orphan galaxies—is taken into account, satellite galaxy populations in CDM and SIDM models can be reconciled. However, even in this case, the inner structures of subhalos are significantly different in the two dark matter models. We study the feasibility of using the weak-lensing signal from the subhalo density profiles to distinguish between the cold and self-interacting dark matter while accounting for the potential contribution of orphan galaxies. We find that the effects of self-interactions on the density profile of subhalos can appear degenerate with subhalo disruption in CDM, when orphans are accounted for. With current error bars from the Subaru Hyper Suprime-Cam Strategic Program, we find that subhalos in the outskirts of clusters (where disruption is less prevalent) can be used to constrain dark matter physics. In the future, the Vera C. Rubin Observatory Legacy Survey of Space and Time will give precise measurements of the weak-lensing profile and can be used to constrain σ T /m at the ∼1 cm2 g−1 level at v ∼ 2000 km s−1.

Effects of stellar-mass primordial black holes on first star formation

ABSTRACT We use cosmological hydrodynamic zoom-in simulations and semi-analytical models to study the effects of primordial black holes (PBHs) on first star formation. Our models self-consistently combine two competing effects: initial (isocurvature) perturbations induced by PBHs and BH accretion feedback. Focusing on PBHs with masses $\sim 30\, \rm M_{\odot }$, we find that the standard picture of first star formation in molecular-cooling minihaloes is not changed by PBHs, as the simulated star-forming gas clouds in the central parsec are very similar to those in the $\rm \Lambda CDM$ case when PBHs make up fPBH ∼ 10−4−0.1 of dark matter. With a dynamical friction time-scale of $\sim 2{\!-\!}10\ \rm Myr$ when the central gas density reaches $10^{5}\ \rm cm^{-3}$, it is also unlikely that PBHs can sink into star-forming discs and affect the evolution of protostars, although they may interact with the stars during the main-sequence stage. At larger scales, PBHs tend to shift star formation to more massive haloes and accelerate structure formation. The latter effect is stronger in regions with higher initial overdensities. For fPBH ∼ 10−4−0.01 (allowed by observational constraints), the collapsed mass fraction of haloes hosting Population III stars is similar (within a factor of ∼2 at z ≲ 30) to that in $\rm \Lambda CDM$, implying that the impact of stellar-mass PBHs on the cosmic star formation history at z ≳ 10 is small. We also find that the Lyman–Werner photons from PBH accretion in atomic-cooling haloes may facilitate the formation of direct-collapse BHs.

Feedback effect on the observable properties of z &amp;gt; 6 AGN

ABSTRACT Active galactic nuclei (AGNs) feedback has a major impact on to the supermassive black hole (SMBH) growth, the properties of the host galaxies, and their cosmic evolution. We investigate the effects of different kinetic feedback prescriptions on the observable properties of AGN and their host galaxies at z &amp;gt; 6 in a suite of zoom-in cosmological simulations. We find that kinetic feedback decreases the column density of the interstellar medium (ISM) in the host galaxy by up to a factor of ≈10, especially when the SMBHs reach high accretion rates ($\approx 10{-}30\, \mathrm{M_\odot \, yr^{-1}}$). In particular, kinetic feedback is required to extend the ISM size to &amp;gt;1 kpc and match the observed sizes of the gas reservoirs in z &amp;gt; 6 AGN host galaxies. Moreover, it produces unobscured lines of sight along which the AGN can be detected in the rest-frame UV band with magnitudes consistent with observed values of z &amp;gt; 6 AGN. The assumed geometry of the outflow plays an important role in shaping the observed properties of high-redshift AGN. We find that a biconical geometry is favoured over a spherical one to reproduce the observed properties, but it overestimates the number of multiple AGN systems detectable in X-ray observations. This result suggests that simplistic BH seeding recipes widely employed in cosmological simulations produce too many X-ray detectable multiple AGN at z = 6−7, thus soliciting the adoption of more physically motivated seeding prescriptions.

Enhanced star formation in z ∼ 6 quasar companions

ABSTRACT Quasars powered by supermassive black holes (&amp;gt;108 M⊙) at z ∼ 6 are predicted to reside in cosmic overdense regions. However, observations so far could not confirm this expectation due to limited statistics. The picture is further complicated by the possible effects of quasar outflows (i.e. feedback) that could either suppress or stimulate the star formation rate (SFR) of companion galaxies, thus modifying the expected bias. Here, we quantify feedback effects on the properties and detectability of companions by comparing cosmological zoom-in simulations of a quasar in which feedback is either included or turned-off. With respect to the no-feedback case, companions (a) directly impacted by the outflow have their SFR increased by a factor of 2−3, and (b) tend to be more massive. Both effects shift the [C ii] 158 μm and UV luminosity functions towards brighter magnitudes. This leads us to conclude that quasar feedback slightly increases the effective quasar bias, boosting the number density of observable quasar companions, in agreement with what has been found around the brightest quasars of recent Atacama Large Millimeter Array (ALMA) [C ii] surveys. Deeper observations performed with James Webb Space Telescope and/or ALMA will improve the statistical significance of this result by detecting a larger number of fainter quasar companions.

Hαemission in local galaxies: star formation, time variability, and the diffuse ionized gas

ABSTRACTThe nebular recombination line H α is widely used as a star formation rate (SFR) indicator in the local and high-redshift Universe. We present a detailed H α radiative transfer study of high-resolution isolated Milky-Way and Large Magellanic Cloud simulations that include radiative transfer, non-equilibrium thermochemistry, and dust evolution. We focus on the spatial morphology and temporal variability of the H α emission, and its connection to the underlying gas and star formation properties. The H α and H β radial and vertical surface brightness profiles are in excellent agreement with observations of nearby galaxies. We find that the fraction of H α emission from collisional excitation amounts to fcol ∼ 5–$10{{\ \rm per\ cent}}$, only weakly dependent on radius and vertical height, and that scattering boosts the H α luminosity by $\sim 40{{\ \rm per\ cent}}$. The dust correction via the Balmer decrement works well (intrinsic H α emission recoverable within 25 per cent), though the dust attenuation law depends on the amount of attenuation itself both on spatially resolved and integrated scales. Important for the understanding of the H α–SFR connection is the dust and helium absorption of ionizing radiation (Lyman continuum [LyC] photons), which are about $f_{\rm abs}\approx 28{{\ \rm per\ cent}}$ and $f_{\rm He}\approx 9{{\ \rm per\ cent}}$, respectively. Together with an escape fraction of $f_{\rm esc}\approx 6{{\ \rm per\ cent}}$, this reduces the available budget for hydrogen line emission by nearly half ($f_{\rm H}\approx 57{{\ \rm per\ cent}}$). We discuss the impact of the diffuse ionized gas, showing – among other things – that the extraplanar H α emission is powered by LyC photons escaping the disc. Future applications of this framework to cosmological (zoom-in) simulations will assist in the interpretation of spectroscopy of high-redshift galaxies with the upcoming James Webb Space Telescope.