No Arabic abstract
We post-process galaxies in the IllustrisTNG simulations with SKIRT radiative transfer calculations to make predictions for the rest-frame near-infrared (NIR) and far-infrared (FIR) properties of galaxies at $zgeq 4$. The rest-frame $K$- and $z$-band galaxy luminosity functions from TNG are overall consistent with observations, despite a $sim 0.4,mathrm{dex}$ underprediction at $z=4$ for $M_{rm z}lesssim -24$. Predictions for the JWST MIRI observed galaxy luminosity functions and number counts are given. We show that the next-generation survey conducted by JWST can detect 500 (30) galaxies in F1000W in a survey area of $500,{rm arcmin}^{2}$ at $z=6$ ($z=8$). As opposed to the consistency in the UV, optical and NIR, we find that TNG, combined with our dust modelling choices, significantly underpredicts the abundance of most dust-obscured and thus most luminous FIR galaxies. As a result, the obscured cosmic star formation rate density (SFRD) and the SFRD contributed by optical/NIR dark objects are underpredicted. The discrepancies discovered here could provide new constraints on the sub-grid feedback models, or the dust contents, of simulations. Meanwhile, although the TNG predicted dust temperature and its relations with IR luminosity and redshift are qualitatively consistent with observations, the peak dust temperature of $zgeq 6$ galaxies are overestimated by about $20,{rm K}$. This could be related to the limited mass resolution of our simulations to fully resolve the porosity of the interstellar medium (or specifically its dust content) at these redshifts.
The James Webb Space Telescop (JWST) promises to revolutionise our understanding of the early Universe, and contrasting its upcoming observations with predictions of the $Lambda$CDM model requires detailed theoretical forecasts. Here, we exploit the large dynamic range of the IllustrisTNG simulation suite, TNG50, TNG100, and TNG300, to derive multi-band galaxy luminosity functions from $z=2$ to $z=10$. We put particular emphasis on the exploration of different dust attenuation models to determine galaxy luminosity functions for the rest-frame ultraviolet (UV), and apparent wide NIRCam bands. Our most detailed dust model is based on continuum Monte Carlo radiative transfer calculations employing observationally calibrated dust properties. This calibration results in constraints on the redshift evolution of the dust attenuation normalisation and dust-to-metal ratios yielding a stronger redshift evolution of the attenuation normalisation compared to most previous theoretical studies. Overall we find good agreement between the rest-frame UV luminosity functions and observational data for all redshifts, also beyond the regimes used for the dust-model calibrations. Furthermore, we also recover the observed high redshift ($z=4-6$) UV luminosity versus stellar mass relation, the H$alpha$ versus star formation rate relation, and the H$alpha$ luminosity function at $z=2$. The bright end ($M_{rm UV}>-19.5$) cumulative galaxy number densities are consistent with observational data. For the F200W NIRCam band, we predict that JWST will detect $sim 80$ ($sim 200$) galaxies with a signal-to-noise ratio of $10$ ($sim 5$) within the NIRCam field of view, $2.2times2.2 ,{rm arcmin}^{2}$, for a total exposure time of $10^5{rm s}$ in the redshift range $z=8 pm 0.5$. These numbers drop to $sim 10$ ($sim 40$) for an exposure time of $10^4{rm s}$.
We present predictions for high redshift ($z=2-10$) galaxy populations based on the IllustrisTNG simulation suite and a full Monte Carlo dust radiative transfer post-processing. Specifically, we discuss the ${rm H}_{alpha}$ and ${rm H}_{beta}$ + $[rm O ,III]$ luminosity functions up to $z=8$. The predicted ${rm H}_{beta}$ + $[rm O ,III]$ luminosity functions are consistent with present observations at $zlesssim 3$ with $lesssim 0.1,{rm dex}$ differences in luminosities. However, the predicted ${rm H}_{alpha}$ luminosity function is $sim 0.3,{rm dex}$ dimmer than the observed one at $zsimeq 2$. Furthermore, we explore continuum spectral indices, the Balmer break at $4000$AA (D4000) and the UV continuum slope $beta$. The median D4000 versus sSFR relation predicted at $z=2$ is in agreement with the local calibration despite a different distribution pattern of galaxies in this plane. In addition, we reproduce the observed $A_{rm UV}$ versus $beta$ relation and explore its dependence on galaxy stellar mass, providing an explanation for the observed complexity of this relation. We also find a deficiency in heavily attenuated, UV red galaxies in the simulations. Finally, we provide predictions for the dust attenuation curves of galaxies at $z=2-6$ and investigate their dependence on galaxy colors and stellar masses. The attenuation curves are steeper in galaxies at higher redshifts, with bluer colors, or with lower stellar masses. We attribute these predicted trends to dust geometry. Overall, our results are consistent with present observations of high redshift galaxies. Future JWST observations will further test these predictions.
We present here a three-dimesional hydrodynamical simulation for star formation. Our aim is to explore the effect of the metal-line cooling on the thermodynamics of the star-formation process. We explore the effect of changing the metallicty of the gas from $Z/Z_{odot}=10^{-4}$ to $Z/Z_{odot}=10^{-2}$. Furthermore, we explore the implications of using the observational abundance pattern of a CEMP-no star, which have been considered to be the missing second-generation stars, the so-called Pop. III.2 stars. In order to pursue our aim, we modelled the microphysics by employing the public astrochemistry package KROME, using a chemical network which includes sixteen chemical species (H, H$^{+}$, H$^{-}$, He, He$^{+}$, He$^{++}$, e$^{-}$, H$_{2}$, H$_{2}^{+}$, C, C$^{+}$, O, O$^{+}$, Si, Si$^{+}$, and Si$^{++}$). We couple KROME with the fully three-dimensional Smoothed-particle hydrodynamics (SPH) code GRADSPH. With this framework we investigate the collapse of a metal-enhanced cloud, exploring the fragmentation process and the formation of stars. We found that the metallicity has a clear impact on the thermodynamics of the collapse, allowing the cloud to reach the CMB temperature floor for a metallicity $Z/Z_{odot}=10^{-2}$, which is in agreement with previous work. Moreover, we found that adopting the abundance pattern given by the star SMSS J031300.36-670839.3 the thermodynamics behavior is very similar to simulations with a metallicity of $Z/Z_{odot}=10^{-2}$, due to the high carbon abundance. As long as only metal line cooling is considered, our results support the metallicity threshold proposed by previous works, which will very likely regulate the first episode of fragmentation and potentially determine the masses of the resulting star clusters.
Hot, dust-obscured galaxies (Hot DOGs) are a population of hyper-luminous obscured quasars identified by WISE. We present ALMA observations of the [CII] fine-structure line and underlying dust continuum emission in a sample of seven of the most extremely luminous (EL; L$_{rm bol}$ $ge$ 10$^{14}$ L$_odot$) Hot DOGs, at redshifts z ~ 3.0-4.6. The [CII] line is robustly detected in four objects, tentatively in one, and likely red-shifted out of the spectral window in the remaining two based on additional data. On average, [CII] is red-shifted by ~ 780 km/s from rest-frame ultraviolet emission lines. EL Hot DOGs exhibit consistently very high ionized gas surface densities, with $Sigma_{rm [CII]}$ ~ 1-2 x 10$^{9}$ L$_odot$ kpc$^{-2}$; as high as the most extreme cases seen in other high-redshift quasars. As a population, EL Hot DOG hosts seem to be roughly centered on the main-sequence of star forming galaxies, but the uncertainties are substantial and individual sources can fall above and below. The average, intrinsic [CII] and dust continuum sizes (FWHMs) are ~ 2.1 kpc and ~ 1.6 kpc, respectively, with a very narrow range of line-to-continuum size ratios, 1.61 $pm$ 0.10, suggesting they could be linearly proportional. The [CII] velocity fields of EL Hot DOGs are diverse: from barely rotating structures, to resolved hosts with ordered, circular motions, to complex, disturbed systems that are likely the result of ongoing mergers. In contrast, all sources display large line-velocity dispersions, FWHM $gtrsim$ 500 km/s, which on average are larger than optically and IR-selected quasars at similar or higher redshifts. We argue that one possible hypothesis for the lack of a common velocity structure, the systematically large dispersion of the ionized gas, and the presence of nearby companion galaxies may be that, rather than a single event, the EL Hot DOG phase could be recurrent.
The bright emission from high-redshift quasars completely conceals their host galaxies in the rest-frame ultraviolet/optical, with detection of the hosts in these wavelengths eluding even the Hubble Space Telescope (HST) using detailed point spread function (PSF) modelling techniques. In this study we produce mock images of a sample of z=7 quasars extracted from the BlueTides simulation, and apply Markov Chain Monte Carlo-based PSF modelling to determine the detectability of their host galaxies with the James Webb Space Telescope (JWST). While no statistically significant detections are made with HST, we predict that at the same wavelengths and exposure times JWST NIRCam imaging will detect ~50% of quasar host galaxies. We investigate various observational strategies, and find that NIRCam wide-band imaging in the long-wavelength filters results in the highest fraction of successful quasar host detections, detecting >80% of the hosts of bright quasars in exposure times of 5 ks. Exposure times of ~5 ks are required to detect the majority of host galaxies in the NIRCam wide-band filters, however even 10 ks exposures with MIRI result in <30% successful host detections. We find no significant trends between galaxy properties and their detectability. The PSF modelling can accurately recover the host magnitudes, radii, and spatial distribution of the larger-scale emission, when accounting for the central core being contaminated by residual quasar flux. Care should be made when interpreting the host properties measured using PSF modelling.