Population Synthesis Models Research Articles

UVCANDELS: Catalogs of Photometric Redshifts and Galaxy Physical Properties

The UltraViolet imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey Fields (UVCANDELS) program provides deep Hubble Space Telescope (HST) F275W and F435W imaging over four CANDELS fields (GOODS-N, GOODS-S, COSMOS, and Extended Groth Strip). We combine this newly acquired UV imaging with existing HST imaging from CANDELS as well as existing ancillary data to obtain robust photometric redshifts and reliable estimates for galaxy physical properties for over 150,000 galaxies in the ∼430 arcmin2 UVCANDELS area. Here, we leverage the power of the new UV photometry to not only improve the photometric redshift measurements in these fields, but also constrain the full redshift probability distribution combining multiple redshift-fitting tools. Furthermore, using the full UV-to-IR photometric data set, we measure the galaxy physical properties by fitting templates from population synthesis models with two different parameterizations (flexible and fixed form) of the star formation histories (SFHs). Compared to the flexible SFH parameterization, we find that the fixed-form SFHs systematically underestimate the galaxy stellar masses, both at the low-mass (≲109 M ⊙) and high-mass (≳1010 M ⊙) end, by as much as ∼0.5 dex. This underestimation is primarily due the limited ability of fixed-form SFH parameterization to simultaneously capture the chaotic nature of star formation in these galaxies.

An upper limit on the spins of merging binary black holes formed through isolated binary evolution

As the sensitivity of ground-based gravitational wave detectors progressively increases, observations of black hole mergers will provide us with the joint distribution of their masses and spins. This will be a critical benchmark to validate different formation scenarios. Merging binary black holes formed through the evolution of isolated binary systems require both components to be stripped of their hydrogen envelopes before core-collapse. The rotation rates of such stripped stars are constrained by the critical rotation limit at their surface, including its deviation from the Keplerian value owing to the outward force provided by radiation. This sets a restriction on their angular momentum content at core-collapse. We aim to determine if this restriction plays a role in the spins of binary black hole mergers. We used detailed calculations of stripped stars with the MESA code at low metallicities ($Z=Z_ $Z_ and $Z_ to determine the dimensionless spins of black holes produced by critically rotating stellar progenitors. To study how such progenitors can arise, we considered their formation through chemically homogeneous evolution (CHE) in binary stars. We used a semi-analytical model to study the physical processes that determine the final angular momentum of CHE binaries, and compared our results against available population synthesis models that rely on detailed binary evolution calculations. We find that above black hole masses of $ 25M_ the dimensionless spin parameter of critically rotating stripped stars ($a=Jc/(GM^2)$) is below unity. This results in an exclusion region at high chirp masses and effective spins that cannot be populated by isolated binary evolution. CHE can produce binaries where both black holes hit this limit, producing a pileup at the boundary of the excluded region. High-spin black holes arise from very low-metallicity CHE systems with short delay times, which merge at higher redshifts. On the other hand, the contribution of CHE to merging binary black holes detected in the third observing run of the LVK collaboration is expected to be dominated by systems with low spins ($ eff <0.5$) that merge near redshift zero. Owing to its higher projected sensitivity and runtime, the fourth observing run of the LVK collaboration can potentially place constraints on the high-spin population and the existence of a limit set by critical rotation.

Counting stars from the integrated spectra of galaxies

Over the last few decades, evolutionary population synthesis models have powered an unmatched leap forward in our understanding of galaxies. From dating the age of the first galaxies in the Universe to providing detailed measurements of the chemical composition of nearby galaxies, the success of this approach built upon simple stellar population (SSP) spectro-photometric models is unquestionable. However, the internal constraints inherent to the construction of SSP models can hinder our ability to analyze the integrated spectra of galaxies in situations where the SSP assumption does not sufficiently hold. Thus, here we revisit the possibilities of fitting galaxy spectra as a linear combination of stellar templates without assuming any a priori knowledge on stellar evolution. We showcase the sensitivity of this alternative approach to changes in the stellar population properties, in particular the direct connection to variations in the stellar initial mass function, as well as its advantages when dealing with noncanonical integrated populations and semi-resolved observations. Furthermore, our analysis demonstrates that the absorption spectra of galaxies can be used to independently constrain stellar evolution theory beyond the limited conditions of the solar neighborhood.

The outcome of collisions between gaseous clumps formed by disk instability

The disk instability model is a promising pathway for giant planet formation in various conditions. At the moment, population synthesis models are used to investigate the outcomes of this theory, where a key ingredient of the disk population evolution are collisions of self-gravitating clumps formed by the disk instabilities. In this study, we explored the wide range of dynamics between the colliding clumps by performing state-of-the-art smoothed particle hydrodynamics simulations with a hydrogen-helium mixture equation of state and investigated the parameter space of collisions between clumps of different ages, masses (1--10 Jupiter mass), various impact conditions (head-on to oblique collisions) and a range of relative velocities. We find that the perfect merger assumption used in population synthesis models is rarely satisfied and that the outcomes of most of the collisions lead to erosion, disruption or a hit-and-run. We also show that in some cases collisions can initiate the dynamical collapse of the clump. We conclude that population synthesis models should abandon the simplifying assumption of perfect merging. Relaxing this assumption will significantly affect the inferred population of planets resulting from the disk instability model.

Constraining the overcontact phase in massive binary evolution. III. Period stability of known B+B and O+B overcontact systems

Binary systems play a crucial role in massive star evolution. Systems composed of B-type and O-type stars are of particular interest due to their potential to lead to very energetic phenomena or the merging of exotic compact objects. We aim to determine the orbital period variations of a sample of B+B and O+B massive overcontact binaries, with the primary objectives of characterizing the evolutionary timescales of these systems and addressing the existing discrepancy between observational data and theoretical predictions derived from population synthesis models. We used Period04 to analyze archival photometric data going back a century for a sample of seven binary systems to measure their orbital periods. We then determine the period variations using a linear fit. We find that the period variation timescales of five truly overcontact binary systems align with the nuclear timescale, in agreement with previous findings for more massive overcontact binaries. Additionally, we noticed a clear distinction between the five systems that had been unambiguously classified as overcontact systems and both SV Cen and VFTS 066, which seem to be evolving on thermal timescales and might be misclassified as overcontact systems. In the case of the five overcontact binaries, our results indicate a noticeable mismatch between the observational data and the theoretical predictions derived from population synthesis models. Furthermore, our results suggest that additional physical mechanisms must be investigated to compare the observed variations more thoroughly with theoretical predictions.

A trifecta of modelling tools: a Bayesian binary black hole model selection combining population synthesis and galaxy formation models

ABSTRACT Gravitational waves (GWs) have revealed surprising properties of binary black hole (BBH) populations, but there is still mystery surrounding how these compact objects evolve. We apply Bayesian inference and an efficient method to calculate the BBH merger rates in the Shark host galaxies, to determine the combination of COMPAS parameters that outputs a population most like the GW sources from the LIGO, Virgo, and KAGRA (LVK) transient catalogue. For our COMPAS models, we calculate the likelihood with and without the dependence on the predicted number of BBH merger events. We find strong correlations between hyper-parameters governing the specific angular momentum (AM) of mass lost during mass transfer, the mass-loss rates of Wolf–Rayet stars via winds and the chemically homogeneous evolution (CHE) formation channel. We conclude that analysing the marginalized and unmarginalized likelihood is a good indicator of whether the population parameters distribution and number of observed events reflect the LVK data. In doing so, we see that the majority of the models preferred in terms of the population-level parameters of the BBHs greatly overpredict the number of events we should have observed to date. Looking at the smaller number of models that perform well with both likelihoods, we find that those with no CHE, AM loss occurring closer to the donor during the first mass-transfer event, and/or higher rates of mass-loss from Wolf–Rayet winds are generally preferred by current data. We find these conclusions to be robust to our choice of selection criteria.

Constraining the Origin of the Nanohertz Gravitational-wave Background by Pulsar Timing Array Observations of Both the Background and Individual Supermassive Binary Black Holes

The gravitational waves (GWs) from supermassive binary black holes (BBHs) have long been sought by pulsar timing array (PTA) experiments, in the forms of both a stochastic GW background (GWB) and individual sources. Evidence for a GWB was reported recently by several PTAs with origins to be determined. Here we use a BBH population synthesis model to investigate the detection probability of individual BBHs by the Chinese PTA (CPTA) and the constraint on the GWB origin that may be obtained by PTA observations of both GWB and individual BBHs. If the detected GWB signal is entirely due to BBHs, a significantly positive redshift evolution (∝ (1 + z)2.07) of the mass scaling relation between supermassive black holes and their host galaxies is required. In this case, we find that the detection probability of individual BBHs is ∼85% or 64% if using a period of 3.4 yr of CPTA observation data, with an expectation of ∼1.9 or 1.0 BBHs detectable with a signal-to-noise ratio ≥3 or 5, and it is expected to increase to >95% if the observation period is extended to 5 yr or longer. Even if the contribution from BBHs to the GWB power signal is as small as ∼10%, a positive detection of individual BBHs can still be expected within an observation period of ∼10 yr. A nondetection of individual BBHs within several years from now jointly with the detected GWB signal can put a strong constraint on the upper limit of the BBH contribution to the GWB signal and help identify/falsify a cosmological origin.

Predicting gravitational wave signals from BPASS White Dwarf Binary and Black Hole Binary populations of a Milky Way-like galaxy model for LISA

Abstract Galactic white dwarf binaries (WDBs) and black hole binaries (BHBs) will be gravitational wave (GW) sources for LISA. Their detection will provide insights into binary evolution and the evolution of our Galaxy through cosmic history. Here, we make predictions of the expected WDB and BHB population within our Galaxy. We combine predictions of the compact remnant binary populations expected by stellar evolution from the detailed Binary Population and Spectral Synthesis (BPASS) code, with a Milky Way analogue galaxy model from the Feedback In Realistic Environment (FIRE) simulations. We use PhenomA and legwork to simulate LISA observations. Both packages make similar predictions that on average four Galactic BHBs and 673 Galactic WDBs are above the signal-to-noise ratio (SNR) threshold of 7 after a four-year mission. We compare these predictions to earlier results using the Binary Star Evolution (BSE) code with the same FIRE model galaxy. We find that BPASS predicts a few more LISA observable Galactic BHBs and a twentieth of the Galactic WDBs. The differences are due to the different physical assumptions that have gone into the binary evolution calculations. These results indicate that the expected population of compact binaries that LISA will detect depends very sensitively on the binary population synthesis models used and thus observations of the LISA population will provide tight constraints on our modelling of binary stars. Finally, from our synthetic populations we have created mock LISA signals that can be used to test and refine data processing methods of the eventual LISA observations.

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.

Pulsar Population Synthesis with Magnetorotational Evolution: Constraining the Decay of the Magnetic Field

We present a population synthesis model for normal radio pulsars in the Galaxy incorporating the latest developments in the field and the magnetorotational evolution processes. Our model considers spin-down with a force-free magnetosphere and the decay of the magnetic field strength and its inclination angle. The simulated pulsar population is fit to a large observation sample that covers the majority of radio surveys using the Markov Chain Monte Carlo technique. We compare the distributions of four major observables—spin period (P), spin-down rate ( Ṗ ), dispersion measure, and radio flux density—using accurate high-dimensional Kolmogorov–Smirnov statistics. We test two B-field decay scenarios, an exponential model motivated by ohmic dissipation and a power-law model motivated by the Hall effect. The former clearly provides a better fit, and it can successfully reproduce the observed pulsar distributions with a decay timescale of 8.3−3.0+3.9 Myr. The result suggests that significant B-field decay in aged pulsars and ohmic dissipation could be the dominant process.

Population synthesis models indicate a need for early and ubiquitous disk substructures

Context. Large millimeter surveys of star-forming regions enable the study of entire populations of planet-forming disks and reveal correlations between their observable properties. The ever-increasing number of these surveys has led to a flourishing of population study, a valuable tool and approach that is spreading in ever more fields. Population studies of disks have shown that the correlation between disk size and millimeter flux could be explained either through disks with strong substructure, or alternatively by the effects of radial inward drift of growing dust particles. Aims. This study aims to constrain the parameters and initial conditions of planet-forming disks and address the question of the need for the presence of substructures in disks and, if needed, their predicted characteristics, based on the large samples of disk sizes, millimeter fluxes, and spectral indices available. Methods. We performed a population synthesis of the continuum emission of disks, exploiting a two-population model (two-pop-py), considering the influence of viscous evolution, dust growth, fragmentation, and transport, varying the initial conditions of the disk and substructure to find the best match with the observed distributions. Disks both with and without substructure have been examined. We obtained the simulated population distribution for the disk sizes, millimeter fluxes, and spectral indices by post-processing the resulting disk profiles (surface density, maximum grain size, and disk temperature). Results. We show that the observed distributions of spectral indices, sizes, and luminosities together can be best reproduced by disks with significant substructure; namely, a perturbation that is strong enough to be able to trap particles, that is formed early in the evolution of the disk, and that is within 0.4 Myr. Agreement is reached by relatively high initial disk masses (10−2.3 M⋆ ⩽ Mdisk ⩽ 10−0.5 M⋆) and moderate levels of turbulence (10−3.5 ⩽ α ⩽ 10−2.5). Other disk parameters play a weaker role. Only opacities with a high absorption efficiency can reproduce the observed spectral indices. Conclusions. Disk population synthesis is a precious tool for investigating and constraining the parameters and initial conditions of planet-forming disks. The generally low observed spectral indices call for significant substructure, like that which planets in the mass range of Saturn to a few Jupiters would induce, to already be present before 0.4 Myr. Our results indicate that substructure, which so far has only been assessed in individual disks, is likely ubiquitous and extends to the whole population, and imply that most “smooth” disks hide unresolved substructure.

Wide Post-common Envelope Binaries from Gaia: Orbit Validation and Formation Models

Astrometry from Gaia DR3 has enabled the discovery of a sample of 3000+ binaries containing white dwarfs (WD) and main-sequence (MS) stars in relatively wide orbits, with orbital periods P orb = (100–1000) days. This population was not predicted by binary population synthesis models before Gaia and—if the Gaia orbits are robust—likely requires very efficient envelope ejection during common envelope evolution (CEE). To assess the reliability of the Gaia solutions, we measured multi-epoch radial velocities (RVs) of 31 WD+MS binary candidates with P orb = (40–300) days and AstroSpectroSB1 orbital solutions. We jointly fit the RVs and astrometry, allowing us to validate the Gaia solutions and tighten constraints on component masses. We find a high success rate for the Gaia solutions, with only 2 out of the 31 systems showing significant discrepancies between their Gaia orbital solutions and our RVs. Joint fitting of RVs and astrometry allows us to directly constrain the secondary-to-primary flux ratio S , and we find S≲0.02 for most objects, confirming the companions are indeed WDs. We tighten constraints on the binaries’ eccentricities, finding a median e ≈ 0.1. These eccentricities are much lower than those of normal MS+MS binaries at similar periods, but much higher than predicted for binaries formed via stable mass transfer. We present MESA single and binary evolution models to explore how the binaries may have formed. The orbits of most binaries in the sample can be produced through CEE that begins when the WD progenitor is an AGB star, corresponding to initial separations of 2–5 au. Roughly 50% of all post-common envelope binaries are predicted to have first interacted on the AGB, ending up in wide orbits like these systems.

Big Galaxies and Big Black Holes: The Massive Ends of the Local Stellar and Black Hole Mass Functions and the Implications for Nanohertz Gravitational Waves

We construct the z = 0 galaxy stellar mass function (GSMF) by combining the GSMF at stellar masses M * ≲ 1011.3 M ⊙ from the census study of Leja et al. and the GSMF of massive galaxies at M * ≳ 1011.5 M ⊙ from the volume-limited MASSIVE galaxy survey. To obtain a robust estimate of M * for local massive galaxies, we use MASSIVE galaxies with M * measured from detailed dynamical modeling or stellar population synthesis modeling (incorporating a bottom-heavy initial mass function) with high-quality spatially resolved spectroscopy. These two independent sets of M * agree to within ∼7%. Our new z = 0 GSMF has a higher amplitude at M * ≳ 1011.5 M ⊙ than previous studies, alleviating prior concerns of a lack of mass growth in massive galaxies between z ∼ 1 and 0. We derive a local black hole mass function (BHMF) from this GSMF and the scaling relation of supermassive black holes (SMBHs) and galaxy masses. The inferred abundance of local SMBHs above ∼1010 M ⊙ is consistent with the number of currently known systems. The predicted amplitude of the nanohertz stochastic gravitational-wave background is also consistent with the levels reported by Pulsar Timing Array teams. Our z = 0 GSMF therefore leads to concordant results in the high-mass regime of the local galaxy and SMBH populations and the gravitational-wave amplitude from merging SMBHs. An exception is that our BHMF yields a z = 0 SMBH mass density that is notably higher than the value estimated from quasars at higher redshifts.

CHAOS. VIII. Far-ultraviolet Spectra of M101 and the Impact of Wolf–Rayet Stars* ∗Based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the Data Archive at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555.

We investigate the stellar and nebular properties of nine H ii regions in the spiral galaxy M101 with far-ultraviolet (FUV; ∼900–2000 Å) and optical (∼3200–10,000 Å) spectra. We detect significant C iii] λλ1907,1909 nebular emission in seven regions, but O iii] λ1666 only in the lowest-metallicity region. We produce new analytic functions of the carbon ionization correction factors as a function of metallicity in order to perform a preliminary C/O abundance analysis. The FUV spectra also contain numerous stellar emission and P-Cygni features that we fit with luminosity-weighted combinations of single-burst Starburst99 and BPASS models. We find that the best-fit Starburst99 models closely match the observed very-high-ionization P-Cygni features, requiring very hot, young (≲3 Myr), metal-enriched massive stars. The youngest stellar populations are strongly correlated with broad He ii emission, nitrogen Wolf–Rayet (WR) FUV and optical spectral features, and enhanced N/O gas abundances. Thus, the short-lived WR phase may be driving excess emission in several N P-Cygni wind features (λ955, λ991, λ1720) that bias the stellar continuum fits to higher metallicities relative to the gas-phase metallicities. Accurate characterization of these H ii regions requires additional inclusion of WR stars in the stellar population synthesis models. Our FUV spectra demonstrate that the ∼900–1200 Å FUV can provide a strong test bed for future WR atmosphere and evolution models.

SMA 200–400 GHz Survey for Dust Properties in the Icy Class II Disks in the Taurus Molecular Cloud

We present a new Submillimeter Array survey of 47 Class II sources in the Taurus–Auriga region. Our observations made 12 independent samples of flux densities over the 200–400 GHz frequency range. We tightly constrained the spectral indices of most sources to a narrow range of 2.0 ± 0.2; only a handful of spatially resolved (e.g., diameter >250 au) disks present larger spectral indices. The simplest interpretation for this result is that the (sub)millimeter luminosities of all of the observed target sources are dominated by very optically thick (e.g., τ ≳ 5) dust thermal emission. Some previous works that were based on the optically thin assumption thus might have underestimated optical depths by at least 1 order of magnitude. Assuming DSHARP dust opacities, this corresponds to underestimates of dust masses by a similar factor. For our specific selected sample, the lower limits of dust masses implied by the optically thick interpretation are 1–3 times higher than those previous estimates that were made based on the optically thin assumption. Moreover, some population synthesis models show that, to explain the observed, narrowly distributed spectral indices, the disks in our selected sample need to have very similar dust temperatures (T dust). Given a specific assumption of median T dust, the maximum grain sizes ( amax ) can also be constrained, which is a few times smaller than 0.1 mm for T dust ∼ 100 K and a few millimeters for T dust ∼ 24 K. The results may indicate that dust grain growth outside the water snow line is limited by the bouncing/fragmentation barriers. This is consistent with the recent laboratory experiments, which indicated that the coagulation of water-ice-coated dust is not efficient, and the water-ice-free dust is stickier and thus can coagulate more efficiently. In the Class II disks, the dust mass budget outside of the water snow line may be largely retained instead of being mostly consumed by planet formation. While Class II disks still possess sufficient dust masses to feed planet formation at a later time, it is unknown whether or not dust coagulation and planet formation can be efficient or natural outside of the water snow line.

Cluster population demographics in NGC 628 derived from stochastic population synthesis models

ABSTRACT The physical properties of star cluster populations offer valuable insights into their birth, evolution, and disruption. However, individual stars in clusters beyond the nearest neighbours of the Milky Way are unresolved, forcing analyses of star cluster demographics to rely on integrated light, a process fraught with uncertainty. Here, we infer the demographics of the cluster population in the benchmark galaxy NGC 628 using data from the Legacy Extra-galactic UV Survey (LEGUS) coupled to a novel Bayesian forward-modelling technique. Our method analyses all 1178 clusters in the LEGUS catalogue, $\sim 4$ times more than prior studies severely affected by completeness cuts. Our results indicate that the cluster mass function is either significantly steeper than the commonly observed slope of $-2$ or is truncated at $\approx 10^{4.5}$ M$_\odot$; the latter possibility is consistent with proposed relations between truncation mass and star formation surface density. We find that cluster disruption is relatively mild for the first $\approx 200$ Myr of cluster evolution; no evidence for mass-dependent disruption is found. We find suggestive but not incontrovertible evidence that inner galaxy clusters may be more prone to disruption and outer galaxy clusters have a more truncated mass function, but confirming or refuting these findings will require larger samples from future observations of outer galaxy clusters. Finally, we find that current stellar track and atmosphere models, along with common forms for cluster mass and age distributions, cannot fully capture all features in the multidimensional photometric distribution of star clusters. While our forward-modelling approach outperforms earlier backward-modelling approaches, some systematic differences persist between observed and modelled photometric distributions.

Fate of supernova progenitors in massive binary systems

ABSTRACT How massive stars end their lives depends on the core mass, core angular momentum, and hydrogen envelopes at death. However, these key physical facets of stellar evolution can be severely affected by binary interactions. In turn, the effectiveness of binary interactions itself varies greatly depending on the initial conditions of the binaries, making the situation much more complex. We investigate systematically how binary interactions influence core–collapse progenitors and their fates. Binary evolution simulations are performed to survey the parameter space of supernova progenitors in solar metallicity binary systems and to delineate major evolutionary paths. We first study fixed binary mass ratios ($q=M_2/M_1$ = 0.5, 0.7, and 0.9) to elucidate the impacts of initial mass and initial separation on the outcomes, treating separately Type Ibc supernova, Type II supernova, accretion-induced collapse (AIC), rapidly rotating supernova (Ibc-R), black hole formation, and long gamma ray burst (long GRB). We then conduct 12 binary population synthesis model calculations, varying the initial condition distributions and binary evolution parameters, to estimate various supernova fractions. We obtain a Milky Way supernova rate $R_{\rm SN} = (1.78$–$2.47) \times 10^{-2} \, {\rm yr}^{-1}$ which is consistent with observations. We find the rates of AIC, Ibc-R, and long GRB to be $\sim 1/100$ the rate of regular supernovae. Our estimated long GRB rates are higher than the observed long GRB rate and close to the low luminosity GRB rate, although care must be taken considering our models are computed with solar metallicity. Furthering binary modelling and improving the inputs one by one will enable more detailed studies of these and other transients associated with massive stars.

Heavily obscured AGN detection: A radio versus X-ray challenge

Context. In the supermassive black hole (SMBH)-galaxy coevolution scenario, heavily obscured active galactic nuclei (AGN) represent a fundamental phase of SMBH growth during which most of the BH mass is accreted and the scaling relations with the host galaxy are set. Obscured nuclei are thought to constitute a major fraction of the whole AGN population, but their statistics and evolution across cosmic time are still highly uncertain. Therefore, it is pivotal to identify new ways to detect this vast and hidden population of growing SMBHs. A promising way to select heavily obscured AGN is through radio emission, which is largely unaffected by obscuration and can be used as a proxy for nuclear activity. Aims. In this work, we study the AGN radio detection effectiveness in the major deep extragalactic surveys, considering different AGN obscuration levels, redshift, and AGN bolometric luminosities. We particularly focus on comparing their radio and X-ray detectability, making predictions for present and future radio surveys. Methods. We extrapolated the predictions of the AGN population synthesis model of the cosmic X-ray background (CXB) to the radio band, by deriving the 1.4 GHz luminosity functions of unobscured (i.e., with hydrogen column densities log NH &lt; 22), obscured (22 &lt; log NH &lt; 24), and Compton-thick (CTK, log NH &gt; 24) AGN. We then used these functions to forecast the number of detectable AGN based on the area, flux limit, and completeness of a given radio survey and compare it with the AGN number resulting from X-ray predictions. Results. When applied to deep extragalactic fields covered both by radio and X-ray observations, we show that, while X-ray selection is generally more effective in detecting unobscured AGN, the surface density of CTK AGN radio detected is on average ten times larger than the X-ray one, and even greater at high redshifts, considering the current surveys and facilities. Our results suggest that thousands of CTK AGN are already present in current radio catalogs, but most of them escaped any detection in the corresponding X-ray observations. We also present expectations for the number of AGN to be detected by the Square Kilometer Array Observatory (SKAO) in its future deep and wide radio continuum surveys, finding that it will be able to detect more than 2000 AGN at z &gt; 6 and tens of them at z &gt; 10, more than half of which are expected to be CTK.

Structure and Kinematics of Star-forming Elliptical Galaxies in SDSS-MaNGA

Using “spatially” resolved spectroscopy, we investigated the characteristics and different modes of formation of stars in elliptical galaxies. We identified an unusual population of 59 star-forming elliptical (SF-E) galaxies in SDSS-MaNGA, our primary sample. To identify these rare star-forming ellipticals, we combined GSWLC-A2 containing outputs of stellar population synthesis models with morphological results from the deep-learning catalog and resolved and integrated properties from the MaNGA Pipe3D value-added catalog. We have also constructed two control samples of star-forming spirals (SF-Sps; 2419 galaxies) and quenched ellipticals (Q-Es; 684 galaxies) to compare with our primary sample of SF-Es. Hα emission line flux of SF-Es is similar to spiral galaxies. The D4000 spectral index indicates that SF-Es have a mixture of old and young stellar populations. Mass-weighted stellar age and metallicity for the SF-Es are lower than the Q-Es and 67% of stellar- and gas-velocity maps of the primary sample show signs of kinematic disturbance. All of these indicate that SF-Es have acquired metal-poor gas through recent mergers or interactions with other galaxies and are forming a new generation of stars. Further, we subdivide our primary sample of SF-Es into four classes based on their bulge to total luminosity ratio (B/T) and spin parameter λre . These four classes have their distinct evolutionary history and modes of formation. Based on these results, we suggest that the Hubble diagram does not accurately capture galaxy evolution processes, and we need a revised morphology diagram like the comb morphology diagram to get a complete picture of the galaxy evolution processes.

Surveying the Onset and Evolution of Supermassive Black Holes at High-z with AXIS

The nature and origin of supermassive black holes (SMBHs) remain an open matter of debate within the scientific community. While various theoretical scenarios have been proposed, each with specific observational signatures, the lack of sufficiently sensitive X-ray observations hinders the progress of observational tests. In this white paper, we present how AXIS will contribute to solving this issue. With an angular resolution of 1.5″ on-axis and minimal off-axis degradation, we designed a deep survey capable of reaching flux limits in the [0.5–2] keV range of approximately 2 × 10−18 erg s−1 cm−2 over an area of 0.13 deg2 in approximately 7 million seconds (7 Ms). Furthermore, we planned an intermediate depth survey covering approximately 2 deg2 and reaching flux limits of about 2 × 10−17 erg s−1 cm−2 in order to detect a significant number of SMBHs with X-ray luminosities (LX) of approximately 1042 erg s−1 up to z∼10. These observations will enable AXIS to detect SMBHs with masses smaller than 105 M⊙, assuming Eddington-limited accretion and a typical bolometric correction for Type II AGN. AXIS will provide valuable information on the seeding and population synthesis models of SMBHs, allowing for more accurate constraints on their initial mass function (IMF) and accretion history from z∼0–10. To accomplish this, AXIS will leverage the unique synergy of survey telescopes such as the JWST, Roman, Euclid, Vera Rubin Telescope, and the new generation of 30 m class telescopes. These instruments will provide optical identification and redshift measurements, while AXIS will discover the smoking gun of nuclear activity, particularly in the case of highly obscured AGN or peculiar UV spectra as predicted and recently observed by the JWST in the early Universe.