Models Of Galaxy Evolution Research Articles

The ALMA Spectroscopic Survey in the HUDF: Deep 1.2 mm Continuum Number Counts

Abstract We present the results from the 1.2 mm continuum image obtained as part of the Atacama Large Millimeter/submillimeter Array Spectroscopic Survey in the Hubble Ultra Deep Field. The 1.2 mm continuum image has a size of 2.9 (4.2) arcmin2 within a primary beam response of 50% (10%) and an rms value of . We detect 35 sources at high significance (Fidelity ≥0.5); 32 have well-characterized near-infrared Hubble Space Telescope counterparts. We estimate the 1.2 mm number counts to flux levels of in two different ways: we first use the detected sources to constrain the number counts and find a significant flattening of the counts below S ν ∼ 0.1 mJy. In a second approach, we constrain the number counts using a probability of deflection statistics (P(D)) analysis. For this latter approach, we describe new methods to accurately measure the noise in interferometric imaging (employing jackknifing in the cube and in the visibility plane). This independent measurement confirms the flattening of the number counts. Our analysis of the differential number counts shows that we are detecting ∼93% (∼100% if we include the lower fidelity detections) of the total continuum dust emission associated with galaxies in the Hubble Ultra Deep Field. The ancillary data allow us to study the dependence of the 1.2 mm number counts on redshift (z = 0−4), galaxy dust mass ( ), stellar mass ( ), and star formation rate ( ). In an accompanying paper we show that the number counts are crucial to constrain galaxy evolution models and the understanding of star-forming galaxies at high redshift.

The rise of active galactic nuclei in the galaxy evolution and assembly semi-analytic model

ABSTRACT We present a new implementation of the GAlaxy Evolution and Assembly (gaea) semi-analytic model, that features an improved modelling of the process of cold gas accretion on to supermassive black hole (SMBHs), derived from both analytic arguments and high-resolution simulations. We consider different scenarios for the loss of angular momentum required for the available cold gas to be accreted on to the central SMBHs, and we compare different combinations of triggering mechanisms, including galaxy mergers and disc instabilities in star-forming discs. We compare our predictions with the luminosity function (LF) observed for active galactic nuclei (AGNs) and we confirm that a non-instantaneous accretion time-scale (either in the form of a low-angular momentum reservoir or as an assumed light-curve evolution) is needed in order to reproduce the measured evolution of the AGN-LF and the so-called AGN-downsizing trend. Moreover, we also study the impact of AGN feedback, in the form of AGN-driven outflows, on the SF properties of model galaxies, using prescriptions derived both from empirical studies and from numerical experiments. We show that AGN-driven outflows are effective in suppressing the residual star formation rate in massive galaxies (>1011 M⊙) without changing their overall assembly history. These winds also affect the SFR of lower mass galaxies, resulting in a too large fraction of passive galaxies at <1010 M⊙. Finally, we study the Eddington ratio distribution as a function of SMBH mass, showing that only objects more massive than 108 M⊙ are already in a self-regulated state as inferred from observations.

The MUSEHubbleUltra Deep Field Survey

Context.The Lyαemitter (LAE) fraction,XLAE, is a potentially powerful probe of the evolution of the intergalactic neutral hydrogen gas fraction. However, uncertainties in the measurement ofXLAEare still under debate.Aims.Thanks to deep data obtained with the integral field spectrograph Multi Unit Spectroscopic Explorer (MUSE), we can measure the evolution of the LAE fraction homogeneously over a wide redshift range ofz ≈ 3–6 for UV-faint galaxies (down to UV magnitudes ofM1500 ≈ −17.75). This is a significantly fainter range than in former studies (M1500 ≤ −18.75) and it allows us to probe the bulk of the population of high-redshift star-forming galaxies.Methods.We constructed a UV-complete photometric-redshift sample following UV luminosity functions and measured the Lyαemission with MUSE using the latest (second) data release from the MUSEHubbleUltra Deep Field Survey.Results.We derived the redshift evolution ofXLAEforM1500 ∈ [ − 21.75; −17.75] for the first time with a equivalent width rangeEW(Lyα) ≥ 65 Å and found low values ofXLAE ≲ 30% atz ≲ 6. The best-fit linear relation isXLAE= 0.07+0.06−0.03z− 0.22+0.12−0.24. ForM1500 ∈ [ − 20.25; −18.75] andEW(Lyα) ≥ 25 Å, ourXLAEvalues are consistent with those in the literature within 1σatz ≲ 5, but our median values are systematically lower than reported values over the whole redshift range. In addition, we do not find a significant dependence ofXLAEonM1500forEW(Lyα) ≥ 50 Å atz ≈ 3–4, in contrast with previous work. The differences inXLAEmainly arise from selection biases for Lyman Break Galaxies (LBGs) in the literature: UV-faint LBGs are more easily selected if they have strong Lyαemission, henceXLAEis biased towards higher values when those samples are used.Conclusions.Our results suggest either a lower increase ofXLAEtowardsz ≈ 6 than previously suggested, or even a turnover ofXLAEatz ≈ 5.5, which may be the signature of a late or patchy reionization process. We compared our results with predictions from a cosmological galaxy evolution model. We find that a model with a bursty star formation (SF) can reproduce our observed LAE fractions much better than models where SF is a smooth function of time.

The ALMA view of the high-redshift relation between supermassive black holes and their host galaxies

Context. The existence of tight correlations between supermassive black holes (BHs) and their host galaxies’ properties in the local Universe suggests a closely linked evolution. Investigating these relations up to the high redshifts (z ≳ 6) is crucial in order to understand the interplay between star formation and BH growth across the cosmic time and to set constraints on galaxy formation and evolution models. In this work, we focus on the relation between BH mass (MBH) and the dynamical mass (Mdyn) of the host galaxy. Aims. Previous works suggest an evolution of the MBH−Mdyn relation with redshift indicating that BH growth precedes the galaxy mass assembly during their co-evolution at z > 3. However, dynamical galaxy masses at high redshift are often estimated through the virial theorem, thus introducing significant uncertainties. Within the scope of this work, our aim is to study the MBH−Mdyn relation of a sample of 2 < z < 7 quasars by constraining their galaxy masses through a full kinematical modelling of the cold gas kinematics, thus avoiding all possible biases and effects introduced by the rough estimates usually adopted so far. Methods. For this purpose, we retrieved public observations of 72 quasar host galaxies observed in [CII]158 μm or CO transitions with the Atacama Large Millimeter Array (ALMA). We then selected those quasars whose line emission is spatially resolved, and performed a kinematic analysis on ALMA observations. We estimated the dynamical mass of the systems by modelling the gas kinematics with a rotating disc, taking into account geometrical and instrumental effects. Our dynamical mass estimates, combined with MBH obtained from literature and our own new CIVλ1550 observations allowed us to investigate the MBH/Mdyn in the early Universe. Results. Overall, we obtained a sample of ten quasars at z ∼ 2−7, in which line emission is detected with high S/N (≳5−10) and the gas kinematics are spatially resolved and dominated by ordered rotation. The estimated dynamical masses place six out of ten quasars above the local relation yielding to MBH/Mdyn ratios ∼10× higher than those estimated in low-z galaxies. On the other hand, we found that four quasars at z ∼ 4−6 have dynamical-to-BH-mass ratios consistent with what is observed in early-type galaxies in the local Universe.

Absorption-line Abundances in the SMC-like Galaxy UGC 5282: Evidence of ISM Dilution from Inflows on Kiloparsec Scales*

Abstract We present a Hubble Space Telescope Cosmic Origins Spectrograph spectrum of the QSO SDSS J095109.12+330745.8 ( ) whose sightline passes through the SMC-like dwarf galaxy UGC 5282 ( , cz = 1577 km s−1), 1.2 kpc in projection from the central H ii region of the galaxy. Damped Lyα (DLA) absorption is detected at the redshift of UGC 5282 with log [N(H i) cm . Analysis of the accompanying S ii, P ii, and O i metal lines yields a neutral gas metallicity, , of [S/H] ≃ [P/H] = −0.80 ± 0.24. The metallicity of ionized gas from the central H ii region measured from its emission lines is [O/H] = −0.37 ± 0.10, a difference of +0.43 ± 0.26 from . This difference δ is consistent with that seen toward H ii regions in other star-forming galaxies and supports the idea that ionized gas near star-forming regions shows systematically higher metallicities than exist in the rest of a galaxy’s neutral interstellar medium (ISM). The positive values of δ found in UGC 5282 (and the other star-forming galaxies) is likely due to infalling low-metallicity gas from the intergalactic medium that mixes with the galaxy’s ISM on kiloparsec scales. This model is also consistent with broad Lyα emission detected at the bottom of the DLA absorption, offset by ∼125 km s−1 from the absorption velocity. Models of galaxy evolution that attempt to replicate population characteristics, such as the mass–metallicity relation, may need to start with a galaxy metallicity represented by rather than that measured traditionally from .

Milky Way Satellite Census. I. The Observational Selection Function for Milky Way Satellites in DES Y3 and Pan-STARRS DR1

Abstract We report the results of a systematic search for ultra-faint Milky Way satellite galaxies using data from the Dark Energy Survey (DES) and Pan-STARRS1 (PS1). Together, DES and PS1 provide multi-band photometry in optical/near-infrared wavelengths over ∼80% of the sky. Our search for satellite galaxies targets ∼25,000 deg2 of the high-Galactic-latitude sky reaching a 10σ point-source depth of ≳22.5 mag in the g and r bands. While satellite galaxy searches have been performed independently on DES and PS1 before, this is the first time that a self-consistent search is performed across both data sets. We do not detect any new high-significance satellite galaxy candidates, recovering the majority of satellites previously detected in surveys of comparable depth. We characterize the sensitivity of our search using a large set of simulated satellites injected into the survey data. We use these simulations to derive both analytic and machine-learning models that accurately predict the detectability of Milky Way satellites as a function of their distance, size, luminosity, and location on the sky. To demonstrate the utility of this observational selection function, we calculate the luminosity function of Milky Way satellite galaxies, assuming that the known population of satellite galaxies is representative of the underlying distribution. We provide access to our observational selection function to facilitate comparisons with cosmological models of galaxy formation and evolution.

UV and U-band luminosity functions from CLAUDS and HSC-SSP – I. Using four million galaxies to simultaneously constrain the very faint and bright regimes to z ∼ 3

ABSTRACT We constrain the rest-frame FUV (1546 Å), NUV (2345 Å), and U-band (3690 Å) luminosity functions (LFs) and luminosity densities (LDs) with unprecedented precision from z ∼ 0.2 to z ∼ 3 (FUV, NUV) and z ∼ 2 (U band). Our sample of over 4.3 million galaxies, selected from the CFHT Large Area U-band Deep Survey (CLAUDS) and HyperSuprime-Cam Subaru Strategic Program (HSC-SSP) data lets us probe the very faint regime (down to MFUV, MNUV, MU ≃ −15 at low redshift), while simultaneously detecting very rare galaxies at the bright end down to comoving densities <10−5 Mpc−3. Our FUV and NUV LFs are well fitted by single-Schechter functions, with faint-end slopes that are very stable up to z ∼ 2. We confirm, but self-consistently and with much better precision than previous studies, that the LDs at all three wavelengths increase rapidly with lookback time to z ∼ 1, and then much more slowly at 1 < z < 2–3. Evolution of the FUV and NUV LFs and LDs at z < 1 is driven almost entirely by the fading of the characteristic magnitude, $M^\star _{\rm UV}$, while at z > 1 it is due to the evolution of both $M^\star _{\rm UV}$ and the characteristic number density $\phi ^\star _{\rm UV}$. In contrast, the U-band LF has an excess of faint galaxies and is fitted with a double-Schechter form; $M^\star _{U}$, both $\phi ^\star _{U}$ components, and the bright-end slope evolve throughout 0.2 < z < 2, while the faint-end slope is constant over at least the measurable 0.05 < z < 0.6. We present tables of our Schechter parameters and LD measurements that can be used for testing theoretical galaxy evolution models and forecasting future observations.

Spatially Resolved UV Diagnostics of AGN Feedback: Radiation Pressure Dominates in a Prototypical Quasar-driven Superwind

Abstract Galactic-scale winds driven by active galactic nuclei (AGN) are often invoked to suppress star formation in galaxy evolution models, but the mechanisms driving these outflows are hotly debated. Two key AGN feedback models are (1) radiation pressure accelerating cool gas and (2) a hot outflowing wind entraining the interstellar medium (ISM). Highly ionized emission-line diagnostics represent a powerful means of differentiating these scenarios because of their sensitivity to the expected compression of the ISM clouds by the hot wind. Here, we report the first spatially resolved UV emission spectroscopy of a prototypical (radio-quiet) quasar-driven superwind around the obscured quasar SDSS J1356+1026 at z = 0.123. We observe ratios of O vi/C iv, N v/C iv, and C iv/He ii that are remarkably similar for outflowing gas clouds ≲100 pc and ≈10 kpc from the nucleus. Such similarity is expected for clouds with AGN radiation-pressure-dominated dynamics. Comparing the observed line emission to models of clouds in balance with radiation pressure and/or a hot wind, we rule out the presence of a dynamically important hot wind and constrain the ratio of hot gas pressure to radiation pressure to P hot/P rad ≲ 0.25 both at ≲100 pc and ≈10 kpc from the nucleus. Moreover, the predictions of the radiation pressure confined cloud models that best fit observed UV line ratios are consistent with the observed diffuse X-ray spectrum. These results indicate that this AGN superwind is driven by radiation pressure or was driven by a hot wind that has since dissipated despite ongoing AGN activity.

A search for the lenses in the Herschel Bright Sources (HerBS) sample

ABSTRACT Verifying that sub-mm galaxies are gravitationally lensed requires time-expensive observations with oversubscribed high-resolution observatories. Here, we aim to strengthen the evidence of gravitational lensing within the Herschel Bright Sources (HerBS) by cross-comparing their positions to optical (SDSS) and near-infrared (VIKING) surveys, in order to search for the foreground lensing galaxy candidates. Resolved observations of the brightest HerBS sources have already shown that most are lensed, and a galaxy evolution model predicts that ∼76 per cent of the total HerBS sources are lensed, although with the SDSS survey we are only able to identify the likely foreground lenses for 25 per cent of the sources. With the near-infrared VIKING survey, however, we are able to identify the likely foreground lenses for 57 per cent of the sources, and we estimate that 82 per cent of the HerBS sources have lenses on the VIKING images even if we cannot identify the lens in every case. We find that the angular offsets between lens and Herschel source are larger than that expected if the lensing is done by individual galaxies. We also find that the fraction of HerBS sources that are lensed falls with decreasing 500-micron flux density, which is expected from the galaxy evolution model. Finally, we apply our statistical VIKING cross-identification to the entire Herschel-ATLAS catalogue, where we also find that the number of lensed sources falls with decreasing 500-micron flux density.

The dust effects on galaxy scaling relations

ABSTRACT Accurate galaxy scaling relations are essential for a successful model of galaxy formation and evolution as they provide direct information about the physical mechanisms of galaxy assembly over cosmic time. We present here a detailed analysis of a sample of nearby spiral galaxies taken from the KINGFISH survey. The photometric parameters of the morphological components are obtained from bulge–disc decompositions using galfit data analysis algorithm, with surface photometry of the sample done beforehand. Dust opacities are determined using a previously discovered correlation between the central face-on dust opacity of the disc and the stellar mass surface density. The method and the library of numerical results previously obtained in Pastrav et al. (2013a,b) are used to correct the measured photometric and structural parameters for projection (inclination), dust, and decomposition effects in order to derive their intrinsic values. Galaxy disc scaling relations are then presented, both the measured (observed) and the intrinsic (corrected) ones, in the optical regime, to show the scale of the biases introduced by the aforementioned effects. The slopes of the size–luminosity relations and the dust versus stellar mass are in agreement with values found in other works. We derive mean dust optical depth and dust/stellar mass ratios of the sample, which we find to be consistent with previous studies of nearby spiral galaxies. While our sample is rather small, it is sufficient to quantify the influence of galaxy environment (dust, in this case) when deriving scaling relations.

On the origin of dust in galaxy clusters at low-to-intermediate redshift

ABSTRACT Stacked analyses of galaxy clusters at low-to-intermediate redshift show signatures attributable to dust, but the origin of this dust is uncertain. We test the hypothesis that the bulk of cluster dust derives from galaxy ejecta. To do so, we employ dust abundances obtained from detailed chemical evolution models of galaxies. We integrate the dust abundances over cluster luminosity functions (one-slope and two-slope Schechter functions). We consider both a hierarchical scenario of galaxy formation and an independent evolution of the three main galactic morphologies: elliptical/S0, spiral and irregular. We separate the dust residing within galaxies from the dust ejected in the intracluster medium. To the latter, we apply thermal sputtering. The model results are compared to low-to-intermediate redshift observations of dust masses. We find that in any of the considered scenarios, elliptical/S0 galaxies contribute negligibly to the present-time intracluster dust, despite producing the majority of gas-phase metals in galaxy clusters. Spiral galaxies, instead, provide both the bulk of the spatially unresolved dust and of the dust ejected into the intracluster medium. The total dust-to-gas mass ratio in galaxy clusters amounts to 10−4, while the intracluster medium dust-to-gas mass ratio amounts to 10−6 at most. These dust abundances are consistent with the estimates of cluster observations at 0.2 < z < 1. We propose that galactic sources, spiral galaxies in particular, are the major contributors to the cluster dust budget.

The Properties of the Interstellar Medium of Galaxies across Time as Traced by the Neutral Atomic Carbon [C i

Abstract We report Atacama Large Millimeter Array observations of the neutral atomic carbon transitions [C i] and multiple CO lines in a sample of ∼30 main-sequence galaxies at , including novel information on [C i] and CO for 7 of such normal objects. We complement our observations with a collection of >200 galaxies with coverage of similar transitions, spanning the z = 0–4 redshift interval and a variety of ambient conditions from local to high-redshift starbursts. We find systematic variations in the [C i]/IR and [C i]/high-J upper (J upper = 7) CO luminosity ratios among the various samples. We interpret these differences as increased dense molecular gas fractions and star formation efficiencies in the strongest high-redshift starbursts with respect to normal main-sequence galaxies. We further report constant / ratios across the galaxy populations and redshifts, suggesting that gas temperatures T exc traced by [C i] do not strongly vary. We find only a mild correlation with T dust and that, generally, T exc ≲ T dust. We fit the line ratios with classical photodissociation region models, retrieving consistently larger densities and intensities of the UV radiation fields in submillimeter galaxies than in main-sequence and local objects. However, these simple models fall short in representing the complexity of a multiphase interstellar medium and should be treated with caution. Finally, we compare our observations with the Santa Cruz semi-analytical model of galaxy evolution, recently extended to simulate submillimeter emission. While we confirm the success in reproducing the CO lines, we find systematically larger [C i] luminosities at fixed IR luminosity than predicted theoretically. This highlights the necessity of improving our understanding of the mechanisms regulating the [C i] emission on galactic scales. We release our data compilation to the community.

Thermal Regulation and the Star-forming Main Sequence

Abstract We argue that the interplay between cosmic rays, the initial mass function (IMF), and star formation plays a crucial role in regulating the star-forming “main sequence.” To explore these phenomena we develop a toy model for galaxy evolution in which star formation is regulated by a combination of a temperature-dependent IMF and heating due to starlight, cosmic rays, and (at very high redshift) the cosmic microwave background. This produces an attractor, near-equilibrium solution which is consistent with observations of the star-forming main sequence over a broad redshift range. Additional solutions to the same equations may correspond to other observed phases of galaxy evolution, including quiescent galaxies. This model makes several falsifiable predictions, including higher metallicities and dust masses than anticipated at high redshift and isotopic abundances in the Milky Way. It also predicts that stellar mass-to-light ratios are lower than produced using a Milky Way–derived IMF, such that inferences of stellar masses and star formation rates for high redshift galaxies are overestimated. In some cases, this may also transform inferred dark matter profiles from core-like to cusp-like.

SDSS-IV MaNGA: Variations in the N/O–O/H Relation Bias Metallicity Gradient Measurements

Abstract In this Letter we use strong line calibrations of the N/O and O/H in Mapping Galaxies at Apache Point Observatory (MaNGA) galaxy survey spaxel data to explore the systematics introduced by variations in N/O on various strong-line metallicity diagnostics. We find radial variations in N/O at fixed O/H that correlate with total galaxy stellar mass and that can induce ∼40% systematic uncertainties in oxygen abundance gradients when nitrogen-dependent abundance calibrations are used. Empirically, we find that departures from the expected N/O are associated with variation in the local star formation efficiency, as predicted by recent chemical evolution models for galaxies, but we cannot rule out other processes such as radial migration also playing a role.

Measuring galaxy-galaxy-galaxy-lensing with higher precision and accuracy

Context. Galaxy-galaxy-galaxy lensing (G3L) is a powerful tool for constraining the three-point correlation between the galaxy and matter distribution and thereby models of galaxy evolution. Aims. We propose three improvements to current measurements of G3L: (i) a weighting of lens galaxies according to their redshift difference, (ii) adaptive binning of the three-point correlation function, and (iii) accounting for the effect of lens magnification by the cosmic large-scale structure. Improvement (i) is designed to improve the precision of the G3L measurement, whereas improvements (ii) and (iii) remove biases of the estimator. We further show how the G3L signal can be converted from angular into physical scales. Methods. The improvements were tested on simple mock data and simulated data based on the Millennium Run with an implemented semi-analytic galaxy model. Results. Our improvements increase the signal-to-noise ratio by 35% on average at angular scales between 0.′1 and 10′ and physical scales between 0.02 and 2 h−1 Mpc. They also remove the bias of the G3L estimator at angular scales below 1′, which was originally up to 40%. The signal due to lens magnification is approximately 10% of the total signal.

The ISM scaling relations in DustPedia late-type galaxies: A benchmark study for the Local Universe

Aims.The purpose of this work is the characterization of the main scaling relations between all of the interstellar medium (ISM) components, namely dust, atomic, molecular, and total gas, and gas-phase metallicity, as well as other galaxy properties, such as stellar mass (Mstar) and galaxy morphology, for late-type galaxies in the Local Universe.Methods.This study was performed by extracting late-type galaxies from the entire DustPedia sample and by exploiting the large and homogeneous dataset available thanks to the DustPedia project. The sample consists of 436 galaxies with morphological stage spanning fromT= 1−10,Mstarfrom 6 × 107to 3 × 1011M⊙, star formation rate from 6 × 10−4to 60M⊙yr−1, and oxygen abundance from 12 + log(O/H) = 8−9.5. Molecular and atomic gas data were collected from the literature and properly homogenized. All the masses involved in our analysis refer to the values within the optical disks of galaxies. The scaling relations involving the molecular gas are studied by assuming both a constant and a metallicity-dependent CO-to-H2conversion factor (XCO). The analysis was performed by means of the survival analysis technique, in order to properly take into account the presence of both detection and nondetection in the data.Results.We confirm that the dust mass correlates very well with the total gas mass, and find –for the first time– that the dust mass correlates better with the atomic gas mass than with the molecular one. We characterize important mass ratios such as the gas fraction, the molecular-to-atomic gas mass ratio, the dust-to-total gas mass ratio (DGR), and the dust-to-stellar mass ratio, and study how they relate to each other, to galaxy morphology, and to gas-phase metallicity. Only the assumption of a metallicity-dependentXCOreproduces the expected decrease of the DGR with increasing morphological stage and decreasing gas-phase metallicity, with a slope of about 1. The DGR, the gas-phase metallicity, and the dust-to-stellar mass ratio are, for our galaxy sample, directly linked to galaxy morphology. The molecular-to-atomic gas mass ratio and the DGR show a positive correlation for low molecular gas fractions, but for galaxies rich in molecular gas this trend breaks down. To our knowledge, this trend has never been found before, and provides new constraints for theoretical models of galaxy evolution and a reference for high-redshift studies. We discuss several scenarios related to this finding.Conclusions.The DustPedia database of late-type galaxies is an extraordinary tool for the study of the ISM scaling relations, thanks to its homogeneous collection of data for the different ISM components. The database is made publicly available to the whole community.

Massive and old quiescent galaxies at high redshift

Aims. Questions of how massive quiescent galaxies rapidly assembled and how abundant they are at high redshift are increasingly important in the study of galaxy formation. Looking at these systems can shed light on the processes of galaxy mass assembly and quenching of the star formation at early epochs. In order to address these questions, we aim to identify and characterize massive quiescent galaxies from z ∼ 2.5 out to the highest redshifts at which these systems can be found. The final purpose is to compare the results with the predictions of state-of-the-art semi-analytical models of galaxy formation and evolution. Methods. We defined observer-frame color–color diagrams to optimally select quiescent galaxies at z > 2.5 and applied them to the COSMOS2015 catalog. We refined the spectral energy distribution (SED) fitting analysis for the selected candidates to confirm their quiescent nature, then derived their number density, mass density, and stellar mass functions. Finally, we compared the results with previous observations and some current semi-analytic models. Results. We selected candidates for quiescent galaxies in the redshift range 2.5 ≲ z ≲ 4.5 from the COSMOS2015 catalog by means of two color–color diagrams. The additional SED fitting analysis allowed us to select 128 galaxies, consistent with being massive (log(M*/M⊙)≥10.6), old (ages ≳0.5 Gyr), and quiescent (log(sSFR [yr−1]) ≤ −10.5) objects at high redshift (2.5 < z < 4.5). Their number and mass densities are in fair agreement with previous observations and, if confirmed, show a discrepancy with current semi-analytical models of galaxy formation and evolution, that underpredict the number of massive quiescent systems up to a factor of ∼12 at 2.5 ≤ z < 3.0 and ∼10 at z ∼ 4.0. The evolution of the stellar mass functions (SMFs) of these systems is similar to previous estimates and indicates a disagreement with models, particularly with regard to the shape of the SMF. Conclusions. The present results add further evidence to the possibility that massive and quiescent galaxies can exist out to at least z ∼ 4. If future spectroscopic observations carried out with, for example, the James Webb Space Telecope (JWST), confirm the substantial presence of such a population, further work on modeling the stellar mass assembly, as well as supermassive black hole accretion and feedback processes at early cosmic epochs, is needed to understand how these systems formed, evolved, and quenched their star formation.

The 21 cm bispectrum during reionization: a tracer of the ionization topology

ABSTRACT We compute the bispectra of the 21cm signal during the epoch of reionization for three different reionization scenarios that are based on a dark matter N-body simulation combined with a self-consistent, semi-numerical model of galaxy evolution and reionization. Our reionization scenarios differ in their trends of ionizing escape fractions (fesc) with the underlying galaxy properties and cover the physically plausible range, i.e. fesc effectively decreasing, being constant, or increasing with halo mass. We find the 21 cm bispectrum to be sensitive to the resulting ionization topologies that significantly differ in their size distribution of ionized and neutral regions throughout reionization. From squeezed to stretched triangles, the 21 cm bispectra features a change of sign from negative to positive values, with ionized and neutral regions representing below-average and above-average concentrations contributing negatively and positively, respectively. The position of the change of sign provides a tracer of the size distribution of the ionized and neutral regions, and allows us to identify three major regimes that the 21 cm bispectrum undergoes during reionization. In particular the regime during the early stages of reionization, where the 21 cm bispectrum tracks the peak of the size distribution of the ionized regions, provides exciting prospects for pinning down reionization with the forthcoming Square Kilometre Array.

The Physical Origins of the Identified and Still Missing Components of the Warm–Hot Intergalactic Medium: Insights from Deep Surveys in the Field of Blazar 1ES1553+113

Abstract The relationship between galaxies and the state/chemical enrichment of the warm–hot intergalactic medium (WHIM) expected to dominate the baryon budget at low-z provides sensitive constraints on structure formation and galaxy evolution models. We present a deep redshift survey in the field of 1ES1553+113, a blazar with a unique combination of ultraviolet (UV)+X-ray spectra for surveys of the circumgalactic/intergalactic medium (CGM/IGM). Nicastro et al. reported the detection of two O vii WHIM absorbers at z = 0.4339 and 0.3551 in its spectrum, suggesting that the WHIM is metal rich and sufficient to close the missing baryons problem. Our survey indicates that the blazar is a member of a z = 0.433 group and that the higher-z O vii candidate arises from its intragroup medium. The resulting bias precludes its use in baryon censuses. The z = 0.3551 candidate occurs in an isolated environment 630 kpc from the nearest galaxy (with stellar mass ), which we show is unexpected for the WHIM. Finally, we characterize the galactic environments of broad H i Lyα absorbers (Doppler widths of b = 40–80 km s−1; T ≲ 4 × 105 K) that provide metallicity-independent WHIM probes. On average, broad Lyα absorbers are ≈2× closer to the nearest luminous (L > 0.25L *) galaxy (700 kpc) than narrow (b < 30 km s−1; T ≲ 4 × 105 K) ones (1300 kpc) but ≈2× further than O vi absorbers (350 kpc). These observations suggest that gravitational collapse heats portions of the IGM to form the WHIM, but with feedback that does not enrich the IGM far beyond galaxy/group halos to levels currently observable in UV/X-ray metal lines.

PRIMUS: Clustering of Star-forming and Quiescent Central Galaxies at 0.2 < z < 0.9

Abstract Previous work has demonstrated that at a given stellar mass, quiescent galaxies are more strongly clustered than star-forming galaxies. The contribution to this signal from central, as opposed to satellite, galaxies is not known, which has strong implications for galaxy evolution models. To investigate the contribution from central galaxies, here we present measurements of the clustering of isolated primary (IP) galaxies, used as a proxy for central galaxies, at 0.2 < z < 0.9 with data from the PRIsm MUlti-Object Survey (PRIMUS) galaxy redshift survey. Using a sample of spectroscopic redshifts for ∼60,000 galaxies with covering 5 deg2 on the sky, we define IP galaxies using isolation cuts in spatial proximity and stellar mass of nearby galaxies. We find that at a fixed stellar mass, quiescent IP galaxies are more strongly clustered than star-forming IP galaxies at z ∼ 0.35 (10σ). Using mock galaxy catalogs based on the recent halo occupation models of Behroozi et al. and designed to replicate the parameters of the PRIMUS survey data set, we find that these clustering differences are due in part to quiescent central galaxies being more strongly clustered than star-forming central galaxies. This is consistent with either distinct stellar-to-halo mass relations for quiescent and star-forming central galaxies, and/or central galaxy assembly bias. We additionally use mock catalogs to assess the dependence of both incompleteness and satellite galaxy contamination in the IP galaxy samples on redshift, galaxy type, and stellar mass, and demonstrate how isolation criteria yield biased subsamples of central galaxies via environmental incompleteness, or the preferential exclusion of central galaxies in overdense environments.