Early Cosmic Times Research Articles

Crimson Behemoth: A massive clumpy structure hosting a dusty AGN at z=4.91

Abstract The current paradigm for the co-evolution of galaxies and their supermassive black holes postulates that dust-obscured active galactic nuclei (AGNs) represent a transitional phase towards a more luminous and unobscured state. However, our understanding of dusty AGNs and their host galaxies at early cosmic times is inadequate due to observational limitations. Here, we present JWST observations of CID-931, an X-ray-detected AGN at a spectroscopic redshift of $z_{\rm spec}=4.91$. Multiband NIRCam imaging from the COSMOS-Web program reveals an unresolved red core, similar to JWST-discovered dusty AGNs. Strikingly, the red core is surrounded by at least eight massive star-forming clumps spread over ${1{^{\prime \prime}_{.}}6} \approx 10\,\,{\rm kpc}$, each of which has a stellar mass of $10^9$–$10^{10}\, M_{\odot }$ and a radius of $\sim$0.1–1 kpc. The whole system amounts to $10^{11}\, M_{\odot }$ in stellar mass, higher than typical star-forming galaxies at the same epoch. In this system, gas inflows and/or complex merger events may trigger clump formation and AGN activity, thus leading to the rapid formation of a massive galaxy hosting a supermassive black hole. Future follow-up observations will provide new insights into the evolution of the galaxy–black hole relationship during such transitional phases in the early universe.

Discovery and Characterization of Two Ultrafaint Dwarfs outside the Halo of the Milky Way: Leo M and Leo K

We report the discovery of two ultrafaint dwarf galaxies, Leo M and Leo K, that lie outside the halo of the Milky Way (MW). Using Hubble Space Telescope imaging of the resolved stars, we create color–magnitude diagrams reaching the oldest main-sequence turnoff of each system and (i) fit for structural parameters of the galaxies; (ii) measure their distances using the luminosity of the horizontal branch stars; (iii) estimate integrated magnitudes and stellar masses; and (iv) reconstruct the star formation histories. Based on their location in the Local Group, neither galaxy is currently within the halo of the MW although Leo K is located ∼26 kpc from the low-mass galaxy Leo T and these two systems may have had a past interaction. Leo M and Leo K have stellar masses of 1.8−0.2+0.3×104 M ⊙ and 1.2 ± 0.2 × 104 M ⊙, and were quenched 10.6−1.1+2.2 Gyr and 12.8−4.2+0.1 Gyr ago, respectively. Given that the galaxies are at farther distances from the MW, it is unlikely that they were quenched by environmental processing. Instead, given their low stellar masses, their early quenching timescales are consistent with the scenario that a combination of reionization and stellar feedback shut down star formation at early cosmic times.

Strong damped Lyman-α absorption in young star-forming galaxies at redshifts 9 to 11.

Primordial neutral atomic gas, mostly composed of hydrogen, is the raw material for star formation in galaxies. However, there are few direct constraints on the amount of neutral atomic hydrogen (Hi) in galaxies at early cosmic times. We analyzed James Webb Space Telescope (JWST) near-infrared spectroscopy of distant galaxies, at redshifts ≳8. From a sample of 12 galaxies, we identified three that show strong damped Lyman-α absorption due to Hi in their local surroundings. The galaxies are located at spectroscopic redshifts of 8.8, 10.2, and 11.4, corresponding to 400 to 600 million years after the Big Bang. They have Hi column densities ≳1022 cm-2, which is an order of magnitude higher than expected for a fully neutral intergalactic medium, and constitute a gas-rich population of young star-forming galaxies.

A JWST investigation into the bar fraction at redshifts 1 ≤ z ≤ 3

ABSTRACT The presence of a stellar bar in a disc galaxy indicates that the galaxy hosts in its main part a dynamically settled disc and that bar-driven processes are taking place in shaping its evolution. Studying the cosmic evolution of the bar fraction in disc galaxies is therefore essential to understand galaxy evolution in general. Previous studies have found, using the Hubble Space Telescope (HST), that the bar fraction significantly declines from the local Universe to redshifts near one. Using the first four pointings from the JWST Cosmic Evolution Early Release Science Survey and the initial public observations for the Public Release Imaging for Extragalactic Research, we extend the studies of the bar fraction in disc galaxies to redshifts 1 ≤ $z$ ≤ 3, that is, for the first time beyond redshift two. We only use galaxies that are also present in the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey on the Extended Groth Strip and Ultra Deep Survey HST observations. An optimized sample of 368 close-to-face-on galaxies is visually classified to find the fraction of bars in disc galaxies in two redshift bins: 1 ≤ $z$ ≤ 2 and 2 < $z$ ≤ 3. The bar fraction decreases from $\approx 17.8^{+ 5.1}_{- 4.8}$ per cent to $\approx 13.8^{+ 6.5}_{- 5.8}$ per cent (from the lower to the higher redshift bin), but is about twice the bar fraction found using bluer HST filters. Our results show that bar-driven evolution might commence at early cosmic times and that dynamically settled discs are already present at a lookback time of ∼11 Gyr.

Lyman continuum leaker candidates at z ∼ 3–4 in the HDUV based on a spectroscopic sample of MUSE LAEs

Context. In recent years, a number of Lyman continuum (LyC) leaker candidates have been found at intermediate redshifts, providing insight into how the Universe was reionised at early cosmic times. There are now around 100 known LyC leakers at all redshifts, which enables us to analyse their properties statistically. Aims. Here, we identify new LyC leaker candidates at z ≈ 3 − 4.5 and compare them to objects from the literature to get an overview of the different observed escape fractions and their relation to the properties of the Lyman α (Lyα) emission line. The aim of this work is to test the indicators (or proxies) for LyC leakage suggested in the literature and to improve our understanding of the kinds of galaxies from which LyC radiation can escape. Methods. We used data from the Hubble Deep Ultraviolet (HDUV) legacy survey to search for LyC emission based on a sample of ≈2000 Lyα emitters (LAEs) detected previously in two surveys with the Multi-Unit Spectroscopic Explorer (MUSE), namely MUSE-Deep and MUSE-Wide. Based on the redshifts and positions of the LAEs, we look for potential LyC leakage in the WFC3/UVIS F336W band of the HDUV survey. The escape fractions are measured and compared in different ways, including spectral energy distribution (SED) fitting performed using the CIGALE software. Results. We add 12 objects to the sample of known LyC leaker candidates (5 highly likely leakers and 7 potential ones), 1 of which was previously known, and compare their Lyα properties to their escape fractions. We find escape fractions of between ∼20% and ∼90%, assuming a high transmission in the intergalactic medium (IGM). We present a method whereby the number of LyC leaker candidates we find is used to infer the underlying average escape fraction of galaxies, which is ≈12%. Conclusion. Based on their Lyα properties, we conclude that LyC leakers are not very different from other high-z LAEs and suggest that most LAEs could be leaking LyC even if this cannot always be detected because of the direction of emission and the transmission properties of the IGM.

Cosmological constraints on f(Q)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$f(Q)$$\\end{document} gravity with redshift space distortion data

In this paper, we explore one of Einstein’s alternative formulations which involves the non-metricity scalar, Q, within the framework of f(Q) theory. Our study focuses on solving the modified Friedmann equations for the case of dust matter, ρ=ρm\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\rho =\\rho _{m}$$\\end{document}, and a form of f(Q)=α+βQn\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$f(Q) = \\alpha + \\beta Q^n$$\\end{document}. We investigate the behavior of our model in both linear (n=1)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$(n=1)$$\\end{document} and nonlinear (n≠1)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$(n\ e 1)$$\\end{document} scenarios at the background and perturbation levels. By employing the Markov chain Monte Carlo (MCMC) method, we constrain our model using observational datasets including redshift space distortion, cosmic chronometers, and Pantheon+\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^+$$\\end{document}. Without using any parameterization of the growth rate index which quantifies the deviation from the Λ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Lambda $$\\end{document}CDM model, both models exhibit good accuracy in describing the redshift space distortion. We further analyze the dynamics of the Universe using cosmography parameters, where our model exhibits a phase transition between deceleration and acceleration phases at z=0.789\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$z=0.789$$\\end{document}. Our findings reveal that our model exhibits a phantom-like behavior based on statefinder diagnostic analysis. Interestingly, the model demonstrates a rich variety of behaviors, resembling either a quintessence-like scenario for (n<1)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$(n<1)$$\\end{document} or phantom-like scenario for (n≥1)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$(n\\ge 1)$$\\end{document}. Using the MCMC best fit and parameterizing the growth index, the evolution of the growth index also depends on the parameter n, either remaining constant (in the linear case) or showing a decreasing trend (in the nonlinear case), indicating a weaker growth rate of density perturbations during earlier cosmic times. Finally, we compare our findings of the growth index with the values obtained in the literature.

A recently quenched galaxy 700 million years after the Big Bang

Local and low-redshift (z < 3) galaxies are known to broadly follow a bimodal distribution: actively star-forming galaxies with relatively stable star-formation rates and passive systems. These two populations are connected by galaxies in relatively slow transition. By contrast, theory predicts that star formation was stochastic at early cosmic times and in low-mass systems1–4. These galaxies transitioned rapidly between starburst episodes and phases of suppressed star formation, potentially even causing temporary quiescence—so-called mini-quenching events5,6. However, the regime of star-formation burstiness is observationally highly unconstrained. Directly observing mini-quenched galaxies in the primordial Universe is therefore of utmost importance to constrain models of galaxy formation and transformation7,8. Early quenched galaxies have been identified out to redshift z < 5 (refs. 9–12) and these are all found to be massive (M⋆ > 1010 M⊙) and relatively old. Here we report a (mini-)quenched galaxy at z = 7.3, when the Universe was only 700 Myr old. The JWST/NIRSpec spectrum is very blue (U–V = 0.16 ± 0.03 mag) but exhibits a Balmer break and no nebular emission lines. The galaxy experienced a short starburst followed by rapid quenching; its stellar mass (4–6 × 108 M⊙) falls in a range that is sensitive to various feedback mechanisms, which can result in perhaps only temporary quenching.

The JWST Advanced Deep Extragalactic Survey: Discovery of an Extreme Galaxy Overdensity at z = 5.4 with JWST/NIRCam in GOODS-S

We report the discovery of an extreme galaxy overdensity at z = 5.4 in the GOODS-S field using James Webb Space Telescope (JWST)/NIRCam imaging from JADES and JEMS alongside JWST/NIRCam wide-field slitless spectroscopy from FRESCO. We identified potential members of the overdensity using Hubble Space Telescope+JWST photometry spanning λ = 0.4–5.0 μm. These data provide accurate and well-constrained photometric redshifts down to m ≈ 29–30 mag. We subsequently confirmed N = 81 galaxies at 5.2 < z < 5.5 using JWST slitless spectroscopy over λ = 3.9–5.0 μm through a targeted line search for Hα around the best-fit photometric redshift. We verified that N = 42 of these galaxies reside in the field, while N = 39 galaxies reside in a density around ∼10 times that of a random volume. Stellar populations for these galaxies were inferred from the photometry and used to construct the star-forming main sequence, where protocluster members appeared more massive and exhibited earlier star formation (and thus older stellar populations) when compared to their field galaxy counterparts. We estimate the total halo mass of this large-scale structure to be 12.6≲log10Mhalo/M⊙≲12.8 using an empirical stellar mass to halo mass relation, which is likely an underestimate as a result of incompleteness. Our discovery demonstrates the power of JWST at constraining dark matter halo assembly and galaxy formation at very early cosmic times.

Emulation of the cosmic dawn 21-cm power spectrum and classification of excess radio models using an artificial neural network

ABSTRACT The cosmic 21-cm line of hydrogen is expected to be measured in detail by the next generation of radio telescopes. The enormous data set from future 21-cm surveys will revolutionize our understanding of early cosmic times. We present a machine learning approach based on an artificial neural network that uses emulation in order to uncover the astrophysics in the epoch of reionization and cosmic dawn. Using a seven-parameter astrophysical model that covers a very wide range of possible 21-cm signals, over the redshift range 6 to 30 and wavenumber range 0.05 to $1 \ \rm {Mpc}^{-1}$ we emulate the 21-cm power spectrum with a typical accuracy of $10 - 20~{{\ \rm per\ cent}}$. As a realistic example, we train an emulator using the power spectrum with an optimistic noise model of the square kilometre array (SKA). Fitting to mock SKA data results in a typical measurement accuracy of 2.8 per cent in the optical depth to the cosmic microwave background, 34 per cent in the star-formation efficiency of galactic haloes, and a factor of 9.6 in the X-ray efficiency of galactic haloes. Also, with our modelling we reconstruct the true 21-cm power spectrum from the mock SKA data with a typical accuracy of $15 - 30~{{\ \rm per\ cent}}$. In addition to standard astrophysical models, we consider two exotic possibilities of strong excess radio backgrounds at high redshifts. We use a neural network to identify the type of radio background present in the 21-cm power spectrum, with an accuracy of 87 per cent for mock SKA data.

Ultra-deep Keck/MOSFIRE Spectroscopic Observations of z ∼ 2 Galaxies: Direct Oxygen Abundances and Nebular Excitation Properties

Using deep near-infrared Keck/MOSFIRE observations, we analyze the rest-optical spectra of eight star-forming galaxies in the COSMOS and GOODS-N fields. We reach integration times of ∼10 hr in the deepest bands, pushing the limits on current ground-based observational capabilities. The targets fall into two redshift bins, of five galaxies at z ∼ 1.7 and three galaxies at z ∼ 2.5, and were selected as likely to yield significant auroral-line detections. Even with long integration times, detection of the auroral lines remains challenging. We stack the spectra together into subsets based on redshift, improving the signal-to-noise ratio on the [O iii]λ4364 auroral emission line and, in turn, enabling a direct measurement of the oxygen abundance for each stack. We compare these measurements to commonly employed strong-line ratios alongside measurements from the literature. We find that the stacks fall within the distribution of z > 1 literature measurements, but a larger sample size is needed to robustly constrain the relationships between strong-line ratios and oxygen abundance at high redshift. We additionally report detections of [O i]λ6302 for nine individual galaxies and composite spectra of 21 targets in the MOSFIRE pointings. We plot their line ratios on the [O iii]λ5008/Hβ versus [O i]λ6302/Hα diagnostic diagram, comparing our targets to local galaxies and H ii regions. We find that the [O i]/Hα ratios in our sample of galaxies are consistent with being produced in gas ionized by α-enhanced massive stars, as has been previously inferred for rapidly forming galaxies at early cosmic times.

The most luminous blue quasars at 3.0 &lt; z &lt; 3.3

We present the analysis of the rest frame ultraviolet and optical spectra of 30 bright blue quasars at z ∼ 3, selected to examine the suitability of active galactic nuclei as cosmological probes. In our previous works, based on pointed XMM-Newton observations, we found an unexpectedly high fraction (≈25%) of X-ray weak quasars in the sample. The latter sources also display a flatter UV continuum and a broader and fainter C IV profile in the archival UV data with respect to their X-ray normal counterparts. Here we present new observations with the Large Binocular Telescope in both the zJ (covering the rest frame ≃2300–3100 Å) and the KS (≃4750–5350 Å) bands. We estimated black hole masses (MBH) and Eddington ratios (λEdd) from the available rest frame optical and UV emission lines (Hβ, Mg II), finding that our z ∼ 3 quasars are on average highly accreting (⟨λEdd⟩≃1.2 and ⟨MBH⟩≃109.7 M⊙), with no difference in λEdd or MBH between X-ray weak and X-ray normal quasars. From the zJ spectra, we derived the properties (e.g. flux, equivalent width) of the main emission lines (Mg II, Fe II), finding that X-ray weak quasars display higher Fe II/Mg II ratios with respect to typical quasars. Fe II/Mg II ratios of X-ray normal quasars are instead consistent with other estimates up to z ≃ 6.5, corroborating the idea of already chemically mature broad line regions at early cosmic time. From the KS spectra, we find that all the X-ray weak quasars present generally weaker [O III] emission (EW &lt; 10 Å) than the normal ones. The sample as a whole, however, abides by the known X-ray-[O III] luminosity correlation, hence the different [O III] properties are likely due to an intrinsically weaker [O III] emission in X-ray weak objects, associated to the shape of the spectral energy distribution. We interpret these results in the framework of accretion-disc winds.

Anisotropic Satellite Galaxy Quenching: A Unique Signature of Energetic Feedback by Supermassive Black Holes?

The quenched fraction of satellite galaxies is aligned with the orientation of the halo’s central galaxy, such that on average, satellites form stars at a lower rate along the major axis of the central. This effect, called anisotropic satellite galaxy quenching (ASGQ), has been found in observational data and cosmological simulations. Analyzing the IllustrisTNG simulation, Martín-Navarro et al. recently argued that ASGQ is caused by anisotropic energetic feedback and constitutes “compelling observational evidence for the role of black holes in regulating galaxy evolution.” In this Letter, we study the causes of ASGQ in state-of-the-art galaxy formation simulations to evaluate this claim. We show that cosmological simulations predict that on average, satellite galaxies along the major axis of the dark matter halo tend to have been accreted at earlier cosmic times and are hosted by subhalos of larger peak halo masses. As a result, a modulation of the quenched fraction with respect to the major axis of the central galaxy is a natural prediction of hierarchical structure formation. We show that ASGQ is predicted by the UniverseMachine galaxy formation model, a model without anisotropic feedback. Furthermore, we demonstrate that even in the IllustrisTNG simulation, anisotropic satellite accretion properties are the main cause of ASGQ. Ultimately, we argue that ASGQ is not a reliable indicator of supermassive black hole feedback in galaxy formation simulations and, thus, should not be interpreted as such in observational data.

Star formation histories of UV-luminous galaxies at z ≃ 6.8: implications for stellar mass assembly at early cosmic times

ABSTRACT The variety of star formation histories (SFHs) of z ≳ 6 galaxies provides important insights into early star formation, but has been difficult to systematically quantify. Some observations suggest that many z ∼ 6–9 galaxies are dominated by ≳200 Myr stellar populations, implying significant star formation at z ≳ 9, while others find that most reionization era galaxies are ≲10 Myr, consistent with little z ≳ 9 star formation. Here, we quantify the distribution of ages of UV-bright ($-22.5\lesssim M_{\rm \small UV}\lesssim -21$) galaxies colour-selected to lie at z ≃ 6.6–6.9, an ideal redshift range to systematically study the SFHs of reionization era galaxies with ground-based observatories and Spitzer. We infer galaxy properties with two SED modelling codes and compare results, finding that stellar masses are largely insensitive to the model, but the inferred ages can vary by an order of magnitude. We infer a distribution of ages assuming a simple, parametric SFH model, finding a median age of ∼30–70 Myr depending on SED model. We quantify the fractions of ≤10 and ≥250 Myr galaxies, finding that these systems comprise ∼15–30 per cent and ∼20–25 per cent of the population, respectively. With a flexible SFH model, the shapes of the SFHs are consistent with those implied by the simple model (e.g. young galaxies have rapidly rising SFHs). However, stellar masses can differ significantly, with those of young systems sometimes being more than an order of magnitude larger with the flexible SFH. We quantify the implications of these results for z ≳ 9 stellar mass assembly and discuss improvements expected from JWST.

Contribution of flares from tidal disruptions of stars to high-redshift AGN

We explore the contribution of tidal disruption events (TDEs) to the flares in Active Galactic Nuclei (AGN) at high redshifts. Using the latest data available from X-ray and optical observations of high-redshift galaxies, in combination with the evolution of their central supermassive black holes, we calculate the contribution of TDE to AGNs as a function of their luminosities. We find that at low redshifts (z &lt; 1), a few percent of all AGN with bolometric luminosities LBol ≲ 1044 erg s−1 may be attributable to possible TDEs. However, this fraction can significantly increase at earlier cosmic times, including up to several tens of percent of the population of AGN at z ≳ 3. TDEs may comprise a significant fraction of the Compton-Thick AGN population at z ≳ 3. The above findings motivate further calibrations with upcoming X-ray missions and spectroscopic surveys targeting TDE-AGN.

An ACA 1 mm survey of HzRGs in the ELAIS-S1: survey description and first results

ABSTRACT Radio-emitting jets might be one of the main ingredients shaping the evolution of massive galaxies in the Universe since early cosmic times. However, identifying early radio active galactic nuclei (AGNs) and confirming this scenario have been hard to accomplish, with studies of samples of radio AGN hosts at z &amp;gt; 2 becoming routinely possible only recently. With the above in mind, we have carried out a survey with the Atacama Compact Array (ACA or Morita Array) at 1.3 mm (rms = 0.15 mJy) of 36 high-redshift radio AGN candidates found within 3.9 deg2 in the ELAIS-S1 field. The work presented here describes the survey and showcases a preliminary set of results. The selection of the sample was based on three criteria making use of infrared (IR) and radio fluxes only. The criterion providing the highest selection rate of high-redshift sources (86 per cent at z &amp;gt; 0.8) is one combining an IR colour cut and radio flux cut (S5.8μm/S3.6μm &amp;gt; 1.3 and $S_{\rm 1.4\, GHz}\gt 1\,$ mJy). Among the sample of 36 sources, 16 show a millimetre (mm) detection. In eight of these cases, the emission has a non-thermal origin. A zsp = 1.58 object, with a mm detection of non-thermal origin, shows a clear spatial offset between the jet-dominated mm continuum emission and that of the host’s molecular gas, as traced by serendipitously detected CO(5-4) emission. Among the objects with serendipitous line detections there is a source with a narrow jet-like region, as revealed by CS(6-5) emission stretching 20 kpc out of the host galaxy.

Dust, Gas, and Metal Content in Star-forming Galaxies at z ∼ 3.3 Revealed with ALMA and Near-IR Spectroscopy

Abstract We conducted submillimeter observations with the Atacama Large Millimeter/submillimeter Array (ALMA) of star-forming galaxies at z ∼ 3.3, whose gas-phase metallicities have been measured previously. We investigated the dust and gas contents of the galaxies at z ∼ 3.3 and studied the interaction of galaxies with their circumgalactic or intergalactic medium at this epoch by probing their gas mass fractions and gas-phase metallicities. Single-band dust continuum emission tracing dust mass and the relation between the gas-phase metallicity and gas-to-dust mass ratio were used to estimate the gas masses. The estimated gas mass fractions and depletion timescales are f gas= 0.20–0.75 and t dep= 0.09–1.55 Gyr. Although the galaxies appear to be tightly distributed around the star-forming main sequence at z ∼ 3.3, both quantities show a wider spread at a fixed stellar mass than expected from the scaling relation, suggesting a large diversity of fundamental gas properties in star-forming galaxies that apparently lie on the main sequence. When we compared gas mass fraction and gas-phase metallicity in star-forming galaxies at z ∼ 3.3 and at lower redshifts, star-forming galaxies at z ∼ 3.3 appear to be more metal poor than local galaxies with similar gas mass fractions. Using the gas regulator model to interpret this offset, we find that this can be explained by a higher mass-loading factor, suggesting that the mass-loading factor in outflows increases at earlier cosmic times.

Tidal disruption events in the first billion years of a galaxy

ABSTRACT Accretion of stars on massive black holes (MBHs) can feed MBHs and generate tidal disruption events (TDEs). We introduce a new physically motivated model to self-consistently treat TDEs in cosmological simulations, and apply it to the assembly of a galaxy with final mass $3\times 10^{10}\, \mathrm{M}_{\odot }$ at z = 6. This galaxy exhibits a TDE rate of $\sim 10^{-5}\, \mathrm{yr}^{-1}$, consistent with local observations but already in place when the Universe was one billion year old. A fraction of the disrupted stars participate in the growth of MBHs, dominating it until the MBH reaches mass $\sim 5 \times 10^5 \, \mathrm{M}_{\odot }$, but their contribution then becomes negligible compared to gas. TDEs could be a viable mechanism to grow light MBH seeds, but fewer TDEs are expected when the MBH becomes sufficiently massive to reach the luminosity of, and be detected as, an active galactic nucleus. Galaxy mergers bring multiple MBHs in the galaxy, resulting in an enhancement of the global TDE rate in the galaxy by ∼1 order of magnitude during $100\, \mathrm{Myr}$ around mergers. This enhancement is not on the central MBH, but caused by the presence of MBHs in the infalling galaxies. This is the first self-consistent study of TDEs in a cosmological environment and highlights that accretion of stars and TDEs are a natural process occurring in a Milky Way-mass galaxy at early cosmic times.

Inception of a first quasar at cosmic dawn

ABSTRACT Earliest quasars at the cosmic dawn are powered by mass accretion on to supermassive black holes of a billion solar masses. Massive black hole (MBH) seeds forming through the direct collapse mechanism are considered the most promising candidates but how do they grow and coevolve with their host galaxies at early cosmic times remains unknown. We here present results from a cosmological radiation hydrodynamical simulation including self-consistent modelling of both Population III (Pop III) and Population II (Pop II) star formation, their radiative and supernova feedback in the host galaxy along with X-ray feedback from an accreting MBH of $\rm 10^5 \, M_{\odot }$ in a halo of $\rm 2 \times 10^9 \, M_{\odot }$ from z = 26 down to z = 16. Our results show that energy deposition from X-rays in the proximity of MBH suppresses Pop III star formation for about 12 Myr while at the same time these X-rays catalyse $\rm H_2$ formation that leads to the formation of a Pop III star cluster of 500 $\rm M_{\odot }$ in the close vicinity of the MBH. We find that mode of star formation for Pop III is episodic and bursty due to the clumpy accretion, while for Pop II it is continuous. The stellar mass of the host galaxy at z ∼ 16 is $\rm 2 \times 10^7 \, M_{\odot }$ with a star formation rate of ${\sim} 0.1\!-\!1 \, \mathrm{ M}_{\odot }\, \mathrm{ yr}^{-1}$. In total, the MBH accretes $\rm 1.5 \times 10^6\, M_{\odot }$ during 120 Myr with the mean accretion rate of ${\sim} 0.01\, \mathrm{ M}_{\odot }\, \mathrm{ yr}^{-1}$ corresponding to an average Eddington fraction of 50 per cent.

Detections of Simultaneous Brightening of γ-Ray and Optical Emissions of a Distant Blazar GB 1508+5714 at Redshift 4.3

Abstract GB 1508+5714 is a high-redshift blazar (z = 4.3), and a spectrally soft γ-ray source has been detected in its direction. By analyzing 11.4 yr Fermi Large Area Telescope data, significant long-term variability of the γ-ray source is confirmed. More importantly, a γ-ray emission flare appeared in an epoch of several tens of days in 2018, when the flux was about four times the value from the global fit. Meanwhile, optical flares were displayed in both the r and i bands from the Zwicky Transient Facility light curves. Detections of the simultaneous γ-ray and optical brightening provide decisive evidence to pin down the association between the γ-ray source and GB 1508+5714, which makes it the first identified γ-ray blazar beyond redshift 4. A broadband spectral energy distribution in the high flux state is constructed, and the origin of the multiwavelength brightening is also briefly discussed. Upcoming wide–deep–fast optical telescopes together with the γ-ray sky surveyors will shed light on the role that the active galactic nucleus jets play in the early cosmic time.

The age–chemical abundance structure of the Galactic disc – II. α-dichotomy and thick disc formation

ABSTRACT We extend our previous work on the age–chemical abundance structure of the Galactic outer disc to the inner disc (4 &amp;lt; r &amp;lt; 8 kpc) based on the SDSS/APOGEE survey. Different from the outer disc, the inner disc stars exhibit a clear bimodal distribution in the [Mg/Fe]–[Fe/H] plane. While a number of scenarios have been proposed in the literature, it remains challenging to recover this bimodal distribution with theoretical models. To this end, we present a chemical evolution model embedding a complex multiphase inner disc formation scenario that matches the observed bimodal [Mg/Fe]–[Fe/H] distribution. In this scenario, the formation of the inner disc is dominated by two main starburst episodes $6\,$Gyr apart with secular, low-level star formation activity in between. In our model, the first starburst occurs at early cosmic times ($t\sim 1\,$ Gyr) and the second one $6\,$ Gyr later at a cosmic time of $t\sim 7\,$ Gyr. Both these starburst episodes are associated with gas accretion events in our model, and are quenched rapidly. The first starburst leads to the formation of the high-α sequence, and the second starburst leads to the formation of the metal-poor low-α sequence. The metal-rich low-α stars, instead, form during the secular evolution phase between the two bursts. Our model shows that the α-dichotomy originates from the rapid suppression of star formation after the first starburst. The two starburst episodes are likely to be responsible for the formation of the geometric thick disc (z &amp;gt;1 kpc), with the old inner thick disc and the young outer thick disc forming during the first and the second starbursts, respectively.