Redshift Range Research Articles

A Census from JWST of Extreme Emission-line Galaxies Spanning the Epoch of Reionization in CEERS

We present a sample of 1165 extreme emission-line galaxies (EELGs) at 4 < z < 9 selected using James Webb Space Telescope (JWST) NIRCam photometry in the Cosmic Evolution Early Release Science (CEERS) program. We use a simple method to photometrically identify EELGs with Hβ + [O iii] (combined) or Hα emission of observed-frame equivalent width (EW) > 5000 Å. JWST/NIRSpec spectroscopic observations of a subset (34) of the photometrically selected EELGs validate our selection method: All spectroscopically observed EELGs confirm our photometric identification of extreme emission, including some cases where the spectral-energy-distribution-derived photometric redshifts are incorrect. We find that the medium-band F410M filter in CEERS is particularly efficient at identifying EELGs, both in terms of including emission lines in the filter and in correctly identifying the continuum between Hβ + [O iii] and Hα in the neighboring broadband filters. We present examples of EELGs that could be incorrectly classified as ultrahigh redshift (z > 12) as a result of extreme Hβ + [O iii] emission blended across the reddest photometric filters. We compare the EELGs to the broader (subextreme) galaxy population in the same redshift range and find that they are consistent with being the bluer, high-EW tail of a broader population of emission-line galaxies. The highest-EW EELGs tend to have more compact emission-line sizes than continuum sizes, suggesting that active galactic nuclei are responsible for at least some of the most extreme EELGs. The photometrically inferred emission-line ratios are consistent with interstellar medium conditions with high ionization and moderately low metallicity, consistent with previous spectroscopic studies.

Retrieval of the physical parameters of galaxies from WEAVE-StePS-like data using machine learning

Context. The William Herschel Telescope Enhanced Area Velocity Explorer (WEAVE) is a new, massively multiplexing spectrograph that allows us to collect about one thousand spectra over a 3 square degree field in one observation. The WEAVE Stellar Population Survey (WEAVE-StePS) in the next 5 years will exploit this new instrument to obtain high-S/N spectra for a magnitude-limited (IAB = 20.5) sample of ∼25 000 galaxies at moderate redshifts (z ≥ 0.3), providing insights into galaxy evolution in this as yet unexplored redshift range. Aims. We aim to test novel techniques for retrieving the key physical parameters of galaxies from WEAVE-StePS spectra using both photometric and spectroscopic (spectral indices) information for a range of noise levels and redshift values. Methods. We simulated ∼105 000 galaxy spectra assuming star formation histories with an exponentially declining star formation rate, covering a wide range of ages, stellar metallicities, specific star formation rates (sSFRs), and dust extinction values. We considered three redshifts (i.e. z = 0.3, 0.55, and 0.7), covering the redshift range that WEAVE-StePS will observe. We then evaluated the ability of the random forest and K-nearest neighbour algorithms to correctly predict the average age, metallicity, sSFR, dust attenuation, and time since the bulk of formation, assuming no measurement errors. We also checked how much the predictive ability deteriorates for different noise levels, with S/NI,obs = 10, 20, and 30, and at different redshifts. Finally, the retrieved sSFR was used to classify galaxies as part of the blue cloud, green valley, or red sequence. Results. We find that both the random forest and K-nearest neighbour algorithms accurately estimate the mass-weighted ages, u-band-weighted ages, and metallicities with low bias. The dispersion varies from 0.08–0.16 dex for age and 0.11–0.25 dex for metallicity, depending on the redshift and noise level. For dust attenuation, we find a similarly low bias and dispersion. For the sSFR, we find a very good constraining power for star-forming galaxies, log sSFR ≳ −11, where the bias is ∼0.01 dex and the dispersion is ∼0.10 dex. However, for more quiescent galaxies, with log sSFR ≲ −11, we find a higher bias, ranging from 0.61 to 0.86 dex, and a higher dispersion, ∼0.4 dex, depending on the noise level and redshift. In general, we find that the random forest algorithm outperforms the K-nearest neighbours. Finally, we find that the classification of galaxies as members of the green valley is successful across the different redshifts and S/Ns. Conclusions. We demonstrate that machine learning algorithms can accurately estimate the physical parameters of simulated galaxies for a WEAVE-StePS-like dataset, even at relatively low S/NI, obs = 10 per Å spectra with available ancillary photometric information. A more traditional approach, Bayesian inference, yields comparable results. The main advantage of using a machine learning algorithm is that, once trained, it requires considerably less time than other methods.

A new census of dust and polycyclic aromatic hydrocarbons at z = 0.7–2 with JWST MIRI

Aims. This paper utilises the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) to extend the observational studies of dust and polycyclic aromatic hydrocarbon (PAH) emission to a new mass and star formation rate (SFR) parameter space beyond our local Universe. The combination of fully sampled spectral energy distributions (SEDs) with multiple mid-infrared (mid-IR) bands and the unprecedented sensitivity of MIRI allows us to investigate dust obscuration and PAH behaviour from z = 0.7 up to z = 2 in typical main-sequence galaxies. Our focus is on constraining the evolution of PAH strength and the dust-obscured luminosity fraction before and during cosmic noon, the epoch of peak star formation activity in the Universe. Methods. In this study, we utilise MIRI multi-band imaging data from the SMILES survey (5 to 25 μm), complemented with NIRCam photometry from the JADES survey (1 to 5 μm), available HST photometry (0.4 to 0.9 μm), and spectroscopic redshifts from the FRESCO and JADES surveys in GOODS-S for 443 star-forming (without dominant active galactic nucleus (AGN)) galaxies at z = 0.7 − 2.0. This redshift range was chosen to ensure that the MIRI data cover mid-IR dust emission. Our methodology involved employing ultraviolet (UV) to IR energy balance SED fitting to robustly constrain the fraction of dust mass in PAHs and dust-obscured luminosity. Additionally, we inferred dust sizes from MIRI 15 μm imaging data, enhancing our understanding of the physical characteristics of dust within these galaxies. Results. We find a strong correlation between the fraction of dust in PAHs (PAH fraction, qPAH) with stellar mass. Moreover, the sub-sample with robust qPAH measurements (N = 216) shows a similar behaviour between qPAH and gas-phase metallicity to that at z ∼ 0, suggesting a universal relation: qPAH is constant (∼3.4%) above a metallicity of Z ∼ 0.5 Z⊙ and decreases to &lt; 1% at metallicities ≲0.3 Z⊙. This indicates that metallicity is a good indicator of the interstellar medium properties that affect the balance between the formation and destruction of PAHs. The lack of a redshift evolution from z ∼ 0 − 2 also implies that above Z ∼ 0.5 Z⊙ the PAH emission effectively traces obscured luminosity and the previous locally calibrated PAH-SFR calibrations remain applicable in this metallicity regime. We observe a strong correlation between the obscured UV luminosity fraction (ratio of obscured to total luminosity) and stellar mass. Above the stellar mass of M* &gt; 5 × 109 M⊙, on average, more than half of the emitted luminosity is obscured, while there exists a non-negligible population of lower-mass galaxies with &gt; 50% obscured fractions. At a fixed mass, the obscured fraction correlates with SFR surface density. This is a result of higher dust covering fractions in galaxies with more compact star-forming regions. Similarly, galaxies with high IRX (IR to UV luminosity) at a given mass or UV continuum slope (β) tend to have higher ΣSFR and shallower attenuation curves, owing to their higher effective dust optical depths and more compact star-forming regions.

An Estimate of the Impact of Reionization on Supermassive Black Hole Growth

The supermassive black holes (SMBHs) that power active galactic nuclei found at z ≥ 6 were formed during the Epoch of Reionization. Because reionization is an inhomogeneous process, the physical properties of SMBH host-galaxy environments will vary spatially during reionization. We construct a semi-analytic model to estimate the impact of reionization on SMBH growth. Using a series of merger trees, reionization models, and black hole (BH) growth models, we find that early reionization can reduce a SMBH’s mass by up to [50, 70, 90]% within dark matter halos of mass [1012, 1011, 1010] M ⊙ by z = 6. Our findings also suggest that the redshift range in which BH growth is impacted by reionization strongly depends on whether the Eddington accretion rate can be exceeded. If so, we find that BH masses are significantly suppressed principally during the early phases of reionization (z ≳ 10), while they are more readily suppressed across the full redshift range if super-Eddington growth is not allowed. We find that the global average impact of reionization may be to reduce the masses of BHs residing in ≲1011 M ⊙ halos by a factor of ≳2. The census of SMBHs being uncovered by the JWST may offer a means to test the basic prediction that more massive BHs reside in cosmological volumes that are reionized at later times.

Charting the Lyman-α escape fraction in the range 2.9 &lt; z &lt; 6.7 and consequences for the LAE reionisation contribution

Context. The escape of Lyman-α photons at redshifts greater than two is an ongoing subject of study and an important quantity to further understanding of Lyman-α emitters (LAEs), the transmission of Lyman-α photons through the interstellar medium and intergalactic medium, and the impact these LAEs have on cosmic reionisation. Aims. This study aims to assess the Lyman-α escape fraction, fesc, Lyα, over the redshift range 2.9 &lt; z &lt; 6.7, focusing on Very Large Telescope/Multi Unit Spectroscopic Explorer (VLT/MUSE) selected, gravitationally lensed, intrinsically faint LAEs. These galaxies are of particular interest as the potential drivers of cosmic reionisation. Methods. We assessed fesc, Lyα in two ways: through an individual study of 96 LAEs behind the A2744 lensing cluster, with James Webb Space Telescope/Near-Infrared Camera (JWST/NIRCam) and HST data, and through a study of the global evolution of fesc, Lyα using the state-of-the-art luminosity functions for LAEs and the UV-selected ‘parent’ population (dust-corrected). We compared these studies to those in the literature based on brighter samples. Results. We find a negligible redshift evolution of fesc, Lyα for our individual galaxies; it is likely that it was washed out by significant intrinsic scatter. We observed a more significant evolution towards higher escape fractions with decreasing UV magnitude and fit this relation. When comparing the two luminosity functions to derive fesc, Lyα in a global sense, we saw agreement with previous literature when integrating the luminosity functions to a bright limit. However, when integrating using a faint limit equivalent to the observational limits of our samples, we observed enhanced values of fesc, Lyα, particularly around z ∼ 6, where fesc, Lyα becomes consistent with 100% escape. This indicates for the faint regimes we sampled that galaxies towards reionisation tend to allow very large fractions of Lyman-α photons to escape. We interpret this as evidence of a lack of any significant dust in these populations; our sample is likely dominated by young, highly star-forming chemically unevolved galaxies. Finally, we assessed the contribution of the LAE population to reionisation using our latest values for fesc, Lyα and the LAE luminosity density. The dependence on the escape fraction of Lyman continuum photons is strong, but for values similar to those observed recently in z ∼ 3 LAEs and high-redshift analogues, LAEs could provide all the ionising emissivity necessary for reionisation.

Prior-informed Active Galactic Nucleus Host Spectral Decomposition Using PyQSOFit

We introduce an improved method for decomposing the emission of active galactic nuclei (AGN) and their host galaxies using templates from principal component analysis (PCA). This approach integrates prior information from PCA with a penalized pixel fitting mechanism that improves the precision and effectiveness of the decomposition process. Specifically, we have reduced the degeneracy and overfitting in AGN host decomposition, particularly for those with low signal-to-noise ratios (SNRs), where traditional methods tend to fail. By applying our method to 76,565 Sloan Digital Sky Survey Data Release 16 quasars with z < 0.8, we achieve a success rate of ≈94%, thus establishing the largest host-decomposed spectral catalog of quasars to date. Our fitting results consider the impact of the host galaxy on the overestimation of the AGN luminosity and black hole mass (M BH). Furthermore, we obtained stellar velocity dispersion (σ ⋆) measurements for 4137 quasars. The slope of the M BH−σ ⋆ relation in this subsample is generally consistent with previous quasar studies beyond the local Universe. Our method provides a robust and efficient approach to disentangle the AGN and host galaxy components across a wide range of SNRs and redshifts.

Type Ia Supernovae from First-generation Stars

Type Ia supernovae (SNe Ia) discovered at redshift z ≲ 2.5 are presumed to be produced from Population I/II stars. In this work, we investigate the production of SNe Ia from Population III binaries in the cosmological framework. We derive the SN Ia rate as a function of redshift in a theoretical context for the production of first-generation stars and evaluate the likelihood of their detection by the James Webb Space Telescope (JWST). Assuming the initial stellar mass function (IMF) favors low-mass stars, as from recent numerical simulations, we found that Population III stars may give rise to a considerable number of SNe Ia at high redshift, and Population III stars may even be the dominant SN Ia producer at z ≳ 6. In an optimistic scenario, we expect ∼1(2) SNe Ia from Population III stars at z ≈ 4(5) for a survey of area 300arcmin2 during a 3 yr period with JWST. The same survey may record more than ∼400 SNe Ia at lower redshift (z ≲ 2.5) but with only about one of them from Population III progenitors. There will be ∼six Population III SNe Ia in the same field of view at redshifts of 5–10. Observational constraints on SN Ia rates at the redshift range of 5–10 can place crucial constraints on the IMF of Population III stars.

Towards a new model-independent calibration of Gamma-Ray Bursts

Current data on baryon acoustic oscillations and Supernovae of Type Ia (SNIa) cover up to z∼2.5. These low-redshift observations play a very important role in the determination of cosmological parameters and have been widely used to constrain the ΛCDM and models beyond the standard, such as the ones with open curvature. To extend this investigation to higher redshifts, Gamma-Ray Bursts (GRBs) stand out as one of the most promising observables. In spite of being transient, they are extremely energetic and can be used to probe the universe up to z∼9.4. They exhibit characteristics that suggest they are potentially standardizable candles and this allows their use to extend the distance ladder beyond SNIa. The use of GRB correlations is still a challenge due to the spread in their intrinsic properties. One of the correlations that can be employed for the standardization is the fundamental plane relation between the peak prompt luminosity, the rest-frame end time of the plateau phase, and its corresponding luminosity, also known as the three-dimensional Dainotti correlation. In this work, we propose an innovative method of calibration of the Dainotti relation which is independent of cosmology. We employ state-of-the-art data on Cosmic Chronometers (CCH) at z≲2 and use the Gaussian Processes Bayesian reconstruction tool. To match the CCH redshift range, we select 20 long GRBs in the range 0.553≤z≤1.96 from the Platinum sample, which consists of well-defined GRB plateau properties that obey the fundamental plane relation. To ensure the generality of our method, we verify that the choice of priors on the parameters of the Dainotti relation and the modelling of CCH uncertainties and covariance have negligible impact on our results. Moreover, we consider the case in which the redshift evolution of the physical features of the plane is accounted for. We find that the use of CCH allows us to identify a sub-sample of GRBs that adhere even more closely to the fundamental plane relation, with an intrinsic scatter of σint=0.20−0.05+0.03 obtained in this analysis when evolutionary effects are considered. In an epoch in which we strive to reduce uncertainties on the variables of the GRB correlations in order to tighten constraints on cosmological parameters, we have found a novel model-independent approach to pinpoint a sub-sample that can thus represent a valuable set of standardizable candles. This allows us to extend the cosmic distance ladder presenting a new catalogue of calibrated luminosity distances up to z=5.

Evolution of HOD and galaxy properties in filaments and nodes of the cosmic web

ABSTRACT We study the evolution of the halo occupation distribution (HOD) and galaxy properties of nodes and filamentary structures obtained by disperse from the Illustris TNG300-1 hydrodynamical simulation, in the redshift range $0 \le z \le 2$. We compute the HOD in filaments and nodes and fit the HOD parameters to study their evolution for both faint and bright galaxies. In nodes, the number of faint galaxies increases with decreasing redshift in the low-mass haloes, while no significant differences are seen in high-mass haloes. Limiting the HOD to bright galaxies shows that haloes increase in mass more than the number of bright galaxies they accrete. For filaments, no large differences in HOD are found for faint galaxies, although for brighter galaxies the behaviour is the same as in nodes. The HOD parametrization suggests that filaments have no effect on the mass required to host a galaxy (central or satellite), whereas nodes do. The results of the study indicate that with this parametrization, filaments do not seem to affect the stellar mass content of galaxies. In contrast, nodes appear to affect haloes with masses below approximately $10^{12.5} h^{-1} {\rm M}_{\odot }$ at local redshift. The analysis of the galaxy colour evolution shows a reddening towards lower redshift, although these processes seem to be more efficient in massive haloes, with a strong effect on bright galaxies. The general evolution suggests that the building of galaxy population within haloes is influenced by both the accretion of faint galaxies and the mass growth of the bright ones.

The GLASS-JWST Early Release Science Program

We release fully reduced spectra obtained with NIRSpec onboard JWST as part of the GLASS-JWST Early Release Science Program and a follow-up Director’s Discretionary Time program 2756. From these 263 spectra of 245 unique sources, acquired with low (R = 30–300) and high dispersion (R ~ 2700) gratings, we derive redshifts for 200 unique sources in the redshift range ɀ = 0–10. We describe the sample selection and characterize its high completeness as a function of redshift and apparent magnitude. Comparison with independent estimates based on different methods and instruments shows that the redshifts are accurate, with 80% differing less than 0.005. We stack the GLASS-JWST spectra to produce the first high-resolution (R ~ 2700) JWST spectral template extending in the rest frame wavelength from 2000 Å to 20 000 Å. Catalogs, reduced spectra, and template are made publicly available to the community.

First measurement of the Mg ii forest correlation function in the Epoch of Reionization

Abstract Strong low-ionization transitions like the Mg ii λ2796, 2804Å doublet are believed to produce a detectable ‘metal-line forest’, if metals pollute the neutral IGM. We measure the auto-correlation of the Mg II forest transmission using ten ground based z ≥ 6.80 quasar spectra probing the redshift range $5.96 &amp;lt; z_{\rm Mg\, {\small II}} &amp;lt; 7.42$ ($z_{\rm Mg\, {\small II}, median} = 6.47$). The correlation function exhibits strong small-scale clustering and a pronounced peak at the doublet velocity (Δv = 768 km s−1) arising from discrete absorbers in the CGM of galaxies. After these strong absorbers are identified and masked the signal is consistent with noise. Our measurements are compared to a suite of models generated by combining a large hydrodynamical simulation with a semi-numerical reionization topology, assuming a simple uniform enrichment model. We obtain a 95% credibility upper limit of [Mg/H] &amp;lt; −3.73 at $z_{\rm Mg\, {\small II},median} = 6.47$, assuming uninformative priors on [Mg/H] and the IGM neutral fraction $x_{\rm {H\, {\small I}}}$. Splitting the data into low-z ($5.96 &amp;lt; z_{\rm Mg\, {\small II}} &amp;lt; 6.47$; $z_{\rm Mg\, {\small II},median} = 6.235$) and high-z ($6.47 &amp;lt; z_{\rm Mg\, {\small II}} &amp;lt; 7.42$; $z_{\rm Mg\, {\small II},median} = 6.72$) subsamples again yields null-detections and 95% upper limits of [Mg/H] &amp;lt; −3.75 and [Mg/H] &amp;lt; −3.45, respectively. These first measurements set the stage for making the Mg II forest an emerging tool to precisely constrain the Universe reionization and enrichment history.

Detection of gamma-ray burst Amati relation based on Hubble data set and Pantheon+ samples

Using gamma-ray bursts as standard candles for cosmological parameter constraints rely on their empirical luminosity relations and low-redshift calibration. In this paper, we examine the Amati relation and its potential corrections based on the A118 sample of higher-quality gamma-ray bursts, using both Hubble data set and Pantheon+ samples as calibration samples in the redshift range of z<1.965\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$z<1.965$$\\end{document}. In calibrating gamma-ray bursts using these two datasets, we employ Gaussian processes to obtain corresponding Hubble diagrams to avoid the dependence on cosmological models in the calibration process. We first divided the low-redshift sample of GRBs into two bins and examined the Amati relation and its potential modifications. We found that under both calibrations, the Amati relation did not show evidence of redshift evolution (68%\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\%$$\\end{document} confidence level). For the other two Amati relations that include redshift evolution terms, the central values of the redshift evolution coefficients deviated from 0, but due to the limitations of the sample size and the increase in the number of parameters, most of the redshift evolution coefficients were not able to be excluded from 0 at the 1σ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sigma $$\\end{document} level. Therefore, to assess their situation across the entire redshift range, we employed MCMC to globally fit three types of Amati relations. By computing AIC and BIC, we found that for the GRB A118 sample, the standard Amati relation remains the most fitting empirical luminosity formula, and no potential redshift evolution trend was observed for two different low-redshift calibrating sources.

Distance estimation of gamma-ray-emitting BL Lac objects from imaging observations

The direct redshift determination of BL Lac objects is highly challenging as the emission in the optical and near-infrared bands is largely dominated by the non-thermal emission from the relativistic jet, which points very close to our line of sight. Therefore, the optical spectra of BL Lac objects often show no spectral lines from the host galaxy. In this work, we aim to overcome this difficulty by attempting to detect the host galaxy and derive redshift constraints based on assumptions on the galaxy magnitude (`imaging redshifts'). Imaging redshifts were derived by obtaining deep optical images under good seeing conditions, making it possible to detect the host galaxy as a weak extension of the point-like source. We then derived the imaging redshift by using the host galaxy as a standard candle, employing two different methods. We determine the imaging redshift for 9 out of 17 blazars that we observed as part of this programme. The redshift range of these targets is 0.28-0.60, and the two methods used to derive the redshift give very consistent results within the uncertainties. We also performed a detailed comparison of the imaging redshifts with those obtained using other methods, such as direct spectroscopic constraints or looking for groups of galaxies close to the blazar. We show that the constraints from the different methods are consistent and that combining the three constraints narrows down the redshift. For example, in the case of J2156.0+1818, which is the most distant source for which we detect the host galaxy, the redshift range is narrowed to $0.63&lt;z&lt;0.71$. This makes the source interesting for future studies of extragalactic background light in the Cherenkov Telescope Array Observatory era.

JADES. The diverse population of infant black holes at 4&lt;z&lt;11: Merging, tiny, poor, but mighty

Spectroscopy with the James Webb Space Telescope has opened the possibility of identifying moderate-luminosity active galactic nuclei (AGNs) in the early Universe, at and beyond the epoch of re-ionisation, complementing previous surveys of much more luminous (and much rarer) quasars. We present 12 new AGNs at 4$&lt;$z$&lt;$7 in the JADES survey (in addition to the previously identified AGN in GN-z11 at z=10.6) revealed through the detection of a broad-line region (BLR) seen in the Balmer emission lines. The depth of JADES, together with the use of three different spectral resolutions, enables us to probe a lower-mass regime relative to previous studies. In a few cases, we find evidence for two broad components of Halpha , which suggests that these could be candidate merging black holes (BHs), although a complex BLR geometry cannot be excluded. The inferred BH masses range from $ 10^7 M_ down to $ 10^5 M_ interestingly probing the regime expected for direct collapse BHs. The inferred AGN bolometric luminosities ($ $ erg/s) imply accretion rates that are $&lt; 0.5$ times the Eddington rate in most cases. However, small BHs, with $ BH 10^6 M_ tend to accrete at Eddington or super-Eddington rates. These BHs at zsim 4-11 are over-massive relative to their host galaxies' stellar masses when compared to the local $ BH -M_ star $ relation, even approaching $ BH star $, as was expected from heavy BH seeds and/or super-Eddington accretion scenarios. However, we find that these early BHs tend to be more consistent with the local relation between $ BH $ and velocity dispersion, as well as between $ BH $ and dynamical mass, suggesting that these are more fundamental and universal relations. On the classical, optical narrow-line excitation-diagnostic diagrams, these AGNs are located in the region that is that is locally occupied by star-forming galaxies, implying that they would be missed by the standard classification techniques if they did not display broad lines. Their location on the diagram is consistent with what is expected for AGNs hosted in metal-poor galaxies ($ Z 0.1-0.2 Z_ The fraction of broad-line AGNs with $ AGN erg/s$ among galaxies in the redshift range of $4&lt;z&lt;6$ is about 10&lt;!PCT!&gt;, suggesting that the contribution of AGNs and their hosts to the re-ionisation of the Universe is $&gt;$10&lt;!PCT!&gt;.

The SRG/eROSITA All-Sky Survey: X-ray beacons at late cosmic dawn

The Spektrum-Roentgen-Gamma (SRG)/extended Roentgen Survey with an Imaging Telescope Array (eROSITA) All-Sky Survey (eRASS) is expected to contain $ quasars that emitted their light when the universe was less than a billion years old, that is, at $z&gt;5.6$. By selection, these quasars populate the bright end of the active galactic nuclei (AGN) X-ray luminosity function, and their space density offers a powerful demographic diagnostic of the parent super-massive black hole (SMBH) population. Of the $ 400$ quasars that have been discovered at $z&gt;5.6$ to date, less than 15&lt;!PCT!&gt; have been X-ray detected. We present a pilot survey to uncover the elusive X-ray luminous end of the distant quasar population. We have designed a quasar selection pipeline based on optical, infrared and X-ray imaging data from DES DR2, VHS DR5, CatWISE2020 and the eRASS (up to its four-pass cumulative version, eRASS:4). The core selection method relies on SED template fitting. We performed optical follow-up spectroscopy with the Magellan /LDSS3 instrument for the redshift confirmation of a subset of candidates. We have further obtained a deeper X-ray image of one of our candidates with Chandra ACIS-S. We report the discovery of five new quasars in the redshift range $5.6 &lt; z &lt; 6.1$. Two of these quasars are detected in eRASS and are, therefore, X-ray ultra-luminous by selection. We also report the detection of these quasars at radio frequencies. The first one is a broad absorption line quasar, which shows significant, order-of-magnitude X-ray dimming over 3.5 years, corresponding to six months in the quasar rest frame. The second X-ray detected quasar is a jetted source with compact morphology. We show that a blazar configuration is likely for this source, making it one of the most distant blazars known to date. With our pilot study, we demonstrate the power of eROSITA as a discovery machine for luminous quasars in the epoch of reionization. The X-ray emission of the two eROSITA detected quasars are likely to be driven by different high-energetic emission mechanisms, a diversity which we will further explore in a future systematic full-hemisphere survey.

The rate and contribution of mergers to mass assembly from NIRCam observations of galaxy candidates up to 13.3 billion years ago

ABSTRACT We present an analysis of the galaxy merger rate in the redshift range $4.0\lt z\lt 9.0$ (i.e. about 1.5 to 0.5 Gyr after the big bang) based on visually identified galaxy mergers from morphological parameter analysis. Our data set is based on high-resolution NIRCam JWST data (a combination of F150W and F2000W broad-band filters) in the low-to-moderate magnification ($\mu \lt 2$) regions of the Abell 2744 cluster field. From a parent set of 675 galaxies $(M_{U}\in [-26.6,-17.9])$, we identify 64 merger candidates from the Gini, $M_{20}$ and asymmetry morphological parameters, leading to a merger fraction $f_m=0.11\pm 0.04$. There is no evidence of redshift evolution of $f_m$ even at the highest redshift considered, thus extending well into the epoch of reionization the constant trend seen previously at $z\lesssim 6$. Furthermore, we investigate any potential redshift-dependent differences in the specific star formation rates between mergers and non-mergers. Our analysis reveals no significant correlation in this regard, with deviations in the studied redshift range typically falling within $(1-1.5)\sigma$ from the null hypothesis that can be attributed to sample variance and measurement errors. Finally, we also demonstrate that the classification of a merging system is robust with respect to the observed (and equivalently rest frame) wavelength of the high-quality JWST broad-band images used. This preliminary study highlights the potential for progress in quantifying galaxy assembly through mergers during the epoch of reionization, with significant sample size growth expected from upcoming large JWST infrared imaging data sets.

The Dark Energy Survey: Cosmology Results with ∼1500 New High-redshift Type Ia Supernovae Using the Full 5 yr Data Set

We present cosmological constraints from the sample of Type Ia supernovae (SNe Ia) discovered and measured during the full 5 yr of the Dark Energy Survey (DES) SN program. In contrast to most previous cosmological samples, in which SNe are classified based on their spectra, we classify the DES SNe using a machine learning algorithm applied to their light curves in four photometric bands. Spectroscopic redshifts are acquired from a dedicated follow-up survey of the host galaxies. After accounting for the likelihood of each SN being an SN Ia, we find 1635 DES SNe in the redshift range 0.10 < z < 1.13 that pass quality selection criteria sufficient to constrain cosmological parameters. This quintuples the number of high-quality z > 0.5 SNe compared to the previous leading compilation of Pantheon+ and results in the tightest cosmological constraints achieved by any SN data set to date. To derive cosmological constraints, we combine the DES SN data with a high-quality external low-redshift sample consisting of 194 SNe Ia spanning 0.025 < z < 0.10. Using SN data alone and including systematic uncertainties, we find ΩM = 0.352 ± 0.017 in flat ΛCDM. SN data alone now require acceleration (q 0 < 0 in ΛCDM) with over 5σ confidence. We find (ΩM,w)=(0.264−0.096+0.074,−0.80−0.16+0.14) in flat wCDM. For flat w 0 w a CDM, we find (ΩM,w0,wa)=(0.495−0.043+0.033,−0.36−0.30+0.36,−8.8−4.5+3.7) , consistent with a constant equation of state to within ∼2σ. Including Planck cosmic microwave background, Sloan Digital Sky Survey baryon acoustic oscillation, and DES 3 × 2pt data gives (ΩM, w) = (0.321 ± 0.007, −0.941 ± 0.026). In all cases, dark energy is consistent with a cosmological constant to within ∼2σ. Systematic errors on cosmological parameters are subdominant compared to statistical errors; these results thus pave the way for future photometrically classified SN analyses.

Chemical abundances along the quasar main sequence

Context. The main sequence of quasars has emerged as a powerful tool for organizing the observational and physical characteristics of type-1 active galactic nuclei (AGNs). Aims. In this study, we present a comprehensive analysis of the metallicity of the gas in the broad-line region, incorporating both new data and previously published findings, to assess the presence of any trend along the main sequence. Methods. We performed a multicomponent analysis on the strongest ultraviolet (UV) and optical emission lines for a sample of 13 radio quiet quasars in the 0.009 ≤ z ≤ 0.472 redshift range, selected based on the availability of multiwavelength data. We employed UV and optical data obtained from the Hubble Space Telescope (mainly from the Cosmic Origins Spectrograph and Faint Object Spectrograph) and several ground-based observatories, respectively. We then measured ten diagnostic ratios and compared them with the prediction of CLOUDY photoionization simulations, identifying the closest photoionization solution to the data. Results. Our investigation reveals a consistent pattern along the main sequence. We observe a systematic progression in metallicity, ranging from subsolar values to metallicity levels several times higher than solar values. Conclusions. These findings underscore the fundamental role of metallicity in correlating with the main sequence of quasars. Extreme metallicity values, at least several dozen times the solar metallicity, are confirmed in low-z AGNs radiating at a high Eddington ratio, although the origin of the extreme enrichment remains open to debate.

Pop-cosmos: A Comprehensive Picture of the Galaxy Population from COSMOS Data

We present pop-cosmos: a comprehensive model characterizing the galaxy population, calibrated to 140,938 (r < 25 selected) galaxies from the Cosmic Evolution Survey (COSMOS) with photometry in 26 bands from the ultraviolet to the infrared. We construct a detailed forward model for the COSMOS data, comprising: a population model describing the joint distribution of galaxy characteristics and its evolution (parameterized by a flexible score-based diffusion model); a state-of-the-art stellar population synthesis model connecting galaxies’ intrinsic properties to their photometry; and a data model for the observation, calibration, and selection processes. By minimizing the optimal transport distance between synthetic and real data, we are able to jointly fit the population and data models, leading to robustly calibrated population-level inferences that account for parameter degeneracies, photometric noise and calibration, and selection. We present a number of key predictions from our model of interest for cosmology and galaxy evolution, including the mass function and redshift distribution; the mass–metallicity-redshift and fundamental metallicity relations; the star-forming sequence; the relation between dust attenuation and stellar mass, star formation rate, and attenuation-law index; and the relation between gas-ionization and star formation. Our model encodes a comprehensive picture of galaxy evolution that faithfully predicts galaxy colors across a broad redshift (z < 4) and wavelength range.

Tully-Fisher relation of late-type galaxies at 0.6 ≤ z ≤ 2.5

We present a study of the stellar and baryonic Tully-Fisher relation within the redshift range of 0.6 ≤ z ≤ 2.5, utilizing observations of star-forming galaxies. This dataset comprises of disk-like galaxies spanning a stellar mass range of 8.89 ≤ log(Mstar [M⊙]) ≤ 11.5, a baryonic mass range of 9.0 ≤ log(Mbar [M⊙]) ≤ 11.5, and a circular velocity range of 1.65 ≤ log(Vc [km/s]) ≤ 2.85. We estimated the stellar masses of these objects using spectral energy distribution fitting techniques, while the gas masses were determined via scaling relations. Circular velocities were directly derived from the rotation curves (RCs), after meticulously correcting for beam smearing and pressure support. Our analysis confirms that our sample adheres to the fundamental mass-size relations of galaxies and reflects the evolution of velocity dispersion in galaxies, in line with previous findings. This reaffirms the reliability of our photometric and kinematic parameters (i.e., Mstar and Vc), thereby enabling a comprehensive examination of the Tully-Fisher relation. To attain robust results, we employed a novel orthogonal likelihood fitting technique designed to minimize intrinsic scatter around the best-fit line, as required at high redshifts. For the stellar Tully-Fisher relation, we obtained a slope of α = 3.03 ± 0.25, an offset of β = 3.34 ± 0.53, and an intrinsic scatter of ζint = 0.08 dex. Correspondingly, the baryonic Tully-Fisher relation yielded α = 3.21 ± 0.28, β = 3.16 ± 0.61, and ζint = 0.09 dex. Our findings indicate a subtle deviation in the stellar and baryonic Tully-Fisher relation with respect to local studies, which is most likely due to the evolutionary processes governing disk formation.