Reconstructing Orbits of Galaxies in Extreme Regions (ROGER) II: reliability of projected phase-space in our understanding of galaxy populations

We connect galaxy properties with their orbital classification by analysing a sample of galaxies with stellar mass M⋆ ≥ 108.5h−1M⊙ residing in and around massive and isolated galaxy clusters with mass M200 > 1015h−1M⊙ at redshift z = 0. The galaxy population is generated by applying the semi-analytic model of galaxy formation sag on the cosmological simulation MultiDark Planck 2. We classify galaxies considering their real orbits (3D) and their projected phase-space position using the roger code (2D). We define five categories: cluster galaxies, galaxies that have recently fallen into a cluster, backsplash galaxies, infalling galaxies, and interloper galaxies. For each class, we analyse the 0.1(g − r) colour, the specific star formation rate (sSFR), and the stellar age, as a function of the stellar mass. For the 3D classes, we find that cluster galaxies have the lowest sSFR, and are the reddest and the oldest, as expected from environmental effects. Backsplash galaxies have properties intermediate between the cluster and recent infaller galaxies. For each 2D class, we find an important contamination by other classes. We find it necessary to separate the galaxy populations in red and blue to perform a more realistic analysis of the 2D data. For the red population, the 2D results are in good agreement with the 3D predictions. Nevertheless, when the blue population is considered, the 2D analysis only provides reliable results for recent infallers, infalling galaxies and interloper galaxies.

The PAU survey: measurement of narrow-band galaxy properties with approximate bayesian computation

Narrow-band imaging surveys allow the study of the spectral characteristics of galaxies without the need of performing their spectroscopic follow-up. In this work, we forward-model the Physics of the Accelerating Universe Survey (PAUS) narrow-band data. The aim is to improve the constraints on the spectral coefficients used to create the galaxy spectral energy distributions (SED) of the galaxy population model in Tortorelli et al. 2020. In that work, the model parameters were inferred from the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) data using Approximate Bayesian Computation (ABC). This led to stringent constraints on the B-band galaxy luminosity function parameters, but left the spectral coefficients only broadly constrained. To address that, we perform an ABC inference using CFHTLS and PAUS data. This is the first time our approach combining forward-modelling and ABC is applied simultaneously to multiple datasets. We test the results of the ABC inference by comparing the narrow-band magnitudes of the observed and simulated galaxies using Principal Component Analysis, finding a very good agreement. Furthermore, we prove the scientific potential of the constrained galaxy population model to provide realistic stellar population properties by measuring them with the SED fitting code CIGALE. We use CFHTLS broad-band and PAUS narrow-band photometry for a flux-limited (i < 22.5) sample of galaxies up to redshift z ∼ 0.8. We find that properties like stellar masses, star-formation rates, mass-weighted stellar ages and metallicities are in agreement within errors between observations and simulations. Overall, this work shows the ability of our galaxy population model to correctly forward-model a complex dataset such as PAUS and the ability to reproduce the diversity of galaxy properties at the redshift range spanned by CFHTLS and PAUS.

Extensive Lensing Survey of Optical and Near-infrared Dark Objects (El Sonido): HST H-faint Galaxies behind 101 Lensing Clusters

We present a Spitzer/IRAC survey of H-faint (H 160 ≳ 26.4, < 5σ) sources in 101 lensing cluster fields. Across a CANDELS/Wide-like survey area of ∼648 arcmin2 (effectively ∼221 arcmin2 in the source plane), we have securely discovered 53 sources in the IRAC Channel-2 band (CH2, 4.5 μm; median CH2 = 22.46 ± 0.11 AB mag) that lack robust HST/WFC3-IR F160W counterparts. The most remarkable source in our sample, namely ES-009 in the field of Abell 2813, is the brightest H-faint galaxy at 4.5 μm known so far (CH2 = 20.48 ± 0.03 AB mag). We show that the H-faint sources in our sample are massive (median M star = 1010.3±0.3 M ⊙), star-forming (median star formation rate = 100 − 40 + 60 M ⊙ yr−1), and dust-obscured (A V = 2.6 ± 0.3) galaxies around a median photometric redshift of z = 3.9 ± 0.4. The stellar continua of 14 H-faint galaxies can be resolved in the CH2 band, suggesting a median circularized effective radius (R e,circ; lensing corrected) of 1.9 ± 0.2 kpc and <1.5 kpc for the resolved and whole samples, respectively. This is consistent with the sizes of massive unobscured galaxies at z ∼ 4, indicating that H-faint galaxies represent the dusty tail of the distribution of a wider galaxy population. Comparing with the ALMA dust continuum sizes of similar galaxies reported previously, we conclude that the heavy dust obscuration in H-faint galaxies is related to the compactness of both stellar and dust continua (R e,circ ∼ 1 kpc). These H-faint galaxies make up 16 − 7 + 13 % of the galaxies in the stellar-mass range of 1010 − 1011.2 M ⊙ at z = 3 ∼ 5, contributing to 8 − 4 + 8 % of the cosmic star formation rate density in this epoch and likely tracing the early phase of massive galaxy formation.

The eMERLIN and EVN View of FR 0 Radio Galaxies

We present the results from high-resolution observations carried out with the eMERLIN UK-array and the European VLBI network (EVN) for a sample of 15 FR 0s, i.e., compact core-dominated radio sources associated with nearby early-type galaxies (ETGs), which represent the bulk of the local radio galaxy population. The 5 GHz eMERLIN observations available for five objects exhibit sub-mJy core components and reveal pc-scale twin jets for four out of five FR 0s once the eMERLIN and JVLA archival visibilities data are combined. The 1.66 GHz EVN observations available for 10 FR 0s display one- and two-sided jetted morphologies and compact cores. The pc-scale core emission contributes, on average, to about one tenth of the total extended radio emission, although we noted an increasing core contribution for flat-/inverted-spectrum sources. We found an unprecedented linear correlation between the pc-scale core luminosity (∼1021.3–1023.6 W Hz−1) and [O III] line luminosity, generally considered as proxy of the accretion power, for a large sample of LINER-type radio-loud low-luminosity active nuclei, all hosted in massive ETGs, which include FR 0s and FR Is. This result represents further evidence of a common jet–disc coupling in FR 0s and FR Is, despite then differing in kpc-scale radio structure. For our objects and for other FR 0 samples reported in the literature, we estimated the jet brightness sidedness ratios, which typically range between one and three. This parameter roughly gauges the jet bulk Lorentz factor Γ, which turns out to range from 1 to 2.5 for most of the sample. This corroborates the scenario that FR 0s are characterized by mildly relativistic jets, possibly as a result of lower-spinning black holes (BHs) than the highly spinning BHs of relativistic-jetted radio galaxies, FR Is.

Assessing the sources of reionisation: a spectroscopic case study of a 30× lensed galaxy at z ∼ 5 with Lyα, C iv, Mg ii, and [Ne iii]

We present a detailed spectroscopic analysis of a galaxy at z ≃ 4.88 that is, by chance, magnified ∼30 × by gravitational lensing. Only three sources at z ≳ 5 are known with such high magnification. This particular source has been shown to exhibit widespread, high equivalent width ${\rm C\, {\small IV}}\, \lambda \, 1549$ emission, implying it is a unique example of a metal-poor galaxy with a hard radiation field, likely representing the galaxy population responsible for cosmic reionisation. Using UV nebular line ratio diagnostics, VLT/X-shooter observations rule out strong AGN activity, indicating a stellar origin of the hard radiation field instead. We present a new detection of ${[\rm Ne\, {\small III}]}\, \lambda \, 3870$ and use the [Ne iii]/[O ii] line ratio to constrain the ionisation parameter and gas-phase metallicity. Closely related to the commonly used [O iii]/[O ii] ratio, our [Ne iii]/[O ii] measurement shows this source is similar to local “Green Pea” galaxies and Lyman-continuum leakers. It furthermore suggests this galaxy is more metal poor than expected from the Fundamental Metallicity Relation, possibly as a consequence of excess gas accretion diluting the metallicity. Finally, we present the highest redshift detection of ${\rm Mg\, {\small II}}\, \lambda \, 2796$, observed at high equivalent width in emission, in contrast to more evolved systems predominantly exhibiting Mg ii absorption. Strong Mg ii emission has been observed in most z ∼ 0 Lyman-continuum leakers known and has recently been proposed as an indirect tracer of escaping ionising radiation. In conclusion, this strongly lensed galaxy, observed just 300 Myr after reionisation ends, enables testing of observational diagnostics proposed to constrain the physical properties of distant galaxies in the JWST/ELT era.

Determining the full satellite population of a Milky Way-mass halo in a highly resolved cosmological hydrodynamic simulation

We investigate the formation of the satellite galaxy population of a Milky Way-mass halo in a very highly resolved magneto-hydrodynamic cosmological zoom-in simulation (baryonic mass resolution mb = 800 $\rm M_{\odot }$). We show that the properties of the central star-forming galaxy, such as the radial stellar surface density profile and star formation history, are: i) robust to stochastic variations associated with the so-called ‘Butterfly Effect’; and ii) well converged over 3.5 orders of magnitude in mass resolution. We find that there are approximately five times as many satellite galaxies at this high resolution compared to a standard ($m_b\sim 10^{4-5}\, \rm M_{\odot }$) resolution simulation of the same system. This is primarily because 2/3rds of the high resolution satellites do not form at standard resolution. A smaller fraction (1/6th) of the satellites present at high resolution form and disrupt at standard resolution; these objects are preferentially low-mass satellites on intermediate- to low-eccentricity orbits with impact parameters ≲ 30 kpc. As a result, the radial distribution of satellites becomes substantially more centrally concentrated at higher resolution, in better agreement with recent observations of satellites around Milky Way-mass haloes. Finally, we show that our galaxy formation model successfully forms ultra-faint galaxies and reproduces the stellar velocity dispersion, half-light radii, and V-band luminosities of observed Milky Way and Local Group dwarf galaxies across 6 orders of magnitude in luminosity (103-$10^{9}\, \rm L_{\odot }$).

Quiescent galaxies in a virialized cluster at redshift 2: Evidence for accelerated size-growth

We present an analysis of the galaxy population in XLSSC 122, an X-ray selected, virialized cluster at redshift z = 1.98. We utilize HST WFC3 photometry to characterize the activity and morphology of spectroscopically confirmed cluster members. The quiescent fraction is found to be $88^{+4}_{-20}$ per cent within 0.5r500, significantly enhanced over the field value of $20^{+2}_{-2}$ per cent at z ∼ 2. We find an excess of “bulge-like” quiescent cluster members with Sersic index n &gt; 2 relative to the field. These galaxies are found to be larger than their field counterparts at 99.6 per cent confidence, being on average $63^{+31}_{-24}$ per cent larger at a fixed mass of M⋆ = 5 × 1010 M⊙. This suggests that these cluster member galaxies have experienced an accelerated size evolution relative to the field at z &gt; 2. We discuss minor mergers as a possible mechanism underlying this disproportionate size growth.

The dwarf galaxy population in nearby clusters from the KIWICS survey

ABSTRACT We analyse a sample of 12 galaxy clusters, from the Kapteyn IAC WEAVE INT Cluster Survey (KIWICS) looking for dwarf galaxy candidates. By using photometric data in the r and g bands from the Wide Field Camera (WFC) at the 2.5-m Isaac Newton Telescope (INT), we select a sample of bright dwarf galaxies (M$_r\, \le$ −15.5 mag) in each cluster and analyse their spatial distribution, stellar colour, and as well as their Sérsic index and effective radius. We quantify the dwarf fraction inside the R200 radius of each cluster, which ranges from ∼0.7 to ∼0.9. Additionally, when comparing the fraction in the inner region with the outermost region of the clusters, we find that the fraction of dwarfs tends to increase going to the outer regions. We also study the clustercentric distance distribution of dwarf and giant galaxies (M$_r\, \lt $ −19.0 mag), and in half of the clusters of our sample, the dwarfs are distributed in a statistically different way as the giants, with the giant galaxies being closer to the cluster centre. We analyse the stellar colour of the dwarf candidates and quantify the fraction of blue dwarfs inside the R200 radius, which is found to be less than ∼0.4, but increases with distance from the cluster centre. Regarding the structural parameters, the Sérsic index for the dwarfs we visually classify as early-type dwarfs tends to be higher in the inner region of the cluster. These results indicate the role that the cluster environment plays in shaping the observational properties of low-mass haloes.

I- A hydrodynamical Clone of the Virgo cluster of galaxies to confirm observationally-driven formation scenarios

At ∼16-17 Mpc from us, the Virgo cluster is a formidable source of information to study cluster formation and galaxy evolution in rich environments. Several observationally-driven formation scenarios arose within the past decade to explain the properties of galaxies that entered the cluster recently and the nature of the last significant merger that the cluster underwent. Confirming these scenarios requires extremely faithful numerical counterparts of the cluster. This paper presents the first Clone, Constrained LOcal and Nesting Environment, simulation of the Virgo cluster within a ∼15 Mpc radius sphere. This cosmological hydrodynamical simulation, with feedback from supernovae and active galactic nuclei, with a ∼3 × 107 M⊙ dark matter particle mass and a minimum cell size of 350 pc in the zoom region, reproduces Virgo within its large scale environment unlike a random cluster simulation. Overall the distribution of the simulated galaxy population matches the observed one including M87. The simulated cluster formation reveals exquisite agreements with observationally-driven scenarios: within the last Gigayear, about 300 small galaxies (M*&gt;107 M⊙) entered the cluster, most of them within the last 500 Myr. The last significant merger event occurred about 2 Gigayears ago: a group with a tenth of the mass of today’s cluster entered from the far side as viewed from the Milky Way. This excellent numerical replica of Virgo will permit studying different galaxy type evolution (jellyfish, backsplash, etc.) as well as feedback phenomena in the cluster core via unbiased comparisons between simulated and observed galaxies and hot gas phase profiles to understand this great physics laboratory.

Extremely deep 150 MHz source counts from the LoTSS Deep Fields

With the advent of new generation low-frequency telescopes, such as the LOw Frequency ARray (LOFAR), and improved calibration techniques, we have now started to unveil the subgigahertz radio sky with unprecedented depth and sensitivity. The LOFAR Two Meter Sky Survey (LoTSS) is an ongoing project in which the whole northern radio sky will be observed at 150 MHz with a sensitivity better than 100 μJy beam−1 at a resolution of 6′′. Additionally, deeper observations are planned to cover smaller areas with higher sensitivity. The Lockman Hole, the Boötes, and the Elais-N1 regions are among the most well known northern extra-galactic fields and the deepest of the LoTSS Deep Fields so far. We exploited these deep observations to derive the deepest radio source counts at 150 MHz to date. Our counts are in broad agreement with those from the literature and show the well known upturn at ≤1 mJy, mainly associated with the emergence of the star-forming galaxy population. More interestingly, our counts show, for the first time a very pronounced drop around S ~ 2 mJy, which results in a prominent “bump” at sub-mJy flux densities. Such a feature was not observed in previous counts’ determinations (neither at 150 MHz nor at a higher frequency). While sample variance can play a role in explaining the observed discrepancies, we believe this is mostly the result of a careful analysis aimed at deblending confused sources and removing spurious sources and artifacts from the radio catalogs. This “drop and bump” feature cannot be reproduced by any of the existing state-of-the-art evolutionary models, and it appears to be associated with a deficiency of active galactic nuclei (AGN) at an intermediate redshift (1 < z < 2) and an excess of low-redshift (z < 1) galaxies and/or AGN.