Low-mass End Research Articles

Tracing astrophysical black hole seeds and primordial black holes with LISA-Taiji network

ABSTRACT In this work, we discuss the improvement that the joint network of LISA and Taiji could provide on exploring two kinds of black hole formation mechanisms. For astrophysical origin, we consider light seed and heavy seed scenarios, and generate populations accordingly. We find that the joint network has the potential to observe growing light seeds in the range 15 < z < 20, while single detector can hardly see, which would shed light on the light seeding mechanism. For the heavy seeds, the joint network only improves the signal-to-noise ratio. For primordial origin, we calculate the event rate at z > 20 and detection rates of LISA and the joint network. The joint network expands LISA’s horizon towards lower mass end, where the event rate is high, so we have better chance observing primordial black holes with the joint network. We also estimate the parameters using Fisher matrices of LISA and the joint network, and find that the joint network significantly improves the estimation.

Probing the low-mass end of the companion mass function for O-type stars

Context. Past observations of O-type stars in the Galaxy have shown that almost all massive stars are part of a binary or higher-order multiple system. Given the wide range of separations at which these companions are found, several observational techniques have been adopted to characterize them. Despite the recent advancements in interferometric and adaptive optics observations, contrasts greater than 4 in the H band have never been reached between 100 and 1000 mas. Aims. Using new adaptive optics (AO) assisted coronagraphic observations, we aim to study the multiplicity properties of a sample of 18 dwarf (or sub-giant) O stars in the galactic field and in OB associations to probe the existence of stellar companions in the angular separation range from 0.″15 to 6″ down to very low mass ratios. Methods. We used VLT/SPHERE to observe simultaneously with the IRDIS and IFS sub-systems 18 O-type stars within 6 kpc and ages between 1 and 5 Myr. The IFS YJH band observations have allowed us to probe the presence of sub-solar companions in a 1.7″ × 1.7″ field-of-view down to magnitude limits of ΔH = 10 at 0.″4. In the wider 12″ × 12″ IRDIS field-of-view, we reached contrasts of ΔK = 12 at 1″, enabling us to look for even fainter companions at larger angular separations and to probe the source density of the surrounding portion of the sky. Results. This paper presents five newly discovered intermediate (< 1″) separation companions, three of which are smaller than 0.2 M⊙. If confirmed by future analyses of proper motions, these new detections represent the lowest-mass companions ever found around O-type stars. Additionally, 29 other sources are found in the IRDIS field-of-view with spurious association probabilities smaller than 5%. Assuming that all sources detected within 1″ are physically bound companions, the observed (uncorrected for bias) fraction of companions for O-type stars between 150 and 900 mas is 0.39 ± 0.15, whereas it increases to 1.6 ± 0.3 in the separation range from 0.″9 to 6″. Conclusions. These findings clearly support the notion that massive stars form almost exclusively in multiple systems, serving as proof of concept that supports the application of larger AO-assisted coronagraphic surveys as a crucial step in placing constraints on the multiplicity properties of massive star companions in regions of the parameter space that have previously gone unexplored. These results also demonstrate that the companion mass function is populated down to the lowest stellar masses.

Running late: testing delayed supermassive black hole growth models against the quasar luminosity function

ABSTRACT Observations of massive galaxies at low redshift have revealed approximately linear scaling relations between the mass of a supermassive black hole (SMBH) and properties of its host galaxy. How these scaling relations evolve with redshift and whether they extend to lower-mass galaxies, however, remain open questions. Recent galaxy formation simulations predict a delayed, or ‘two-phase,’ growth of SMBHs: slow, highly intermittent BH growth due to repeated gas ejection by stellar feedback in low-mass galaxies, followed by more sustained gas accretion that eventually brings BHs on to the local scaling relations. The predicted two-phase growth implies a steep increase, or ‘kink,’ in BH-galaxy scaling relations at a stellar mass $\rm {M}_{*}\sim 5\times 10^{10}$ M⊙. We develop a parametric, semi-analytic model to compare different SMBH growth models against observations of the quasar luminosity function (QLF) at z ∼ 0.5−4. We compare models in which the relation between SMBH mass and galaxy mass is purely linear versus two-phase models. The models are anchored to the observed galaxy stellar mass function, and the BH mass functions at different redshifts are consistently connected by the accretion rates contributing to the QLF. The best fits suggest that two-phase evolution is significantly preferred by the QLF data over a purely linear scaling relation. Moreover, when the model parameters are left free, the two-phase model fits imply a transition mass consistent with that predicted by simulations. Our analysis motivates further observational tests, including measurements of BH masses and active galactic nuclei activity at the low-mass end, which could more directly test two-phase SMBH growth.

Pair-instability Mass Loss for Top-down Compact Object Mass Calculations

Abstract Population synthesis relies on semi-analytic formulae to determine masses of compact objects from the (helium or carbon-oxygen) cores of collapsing stars. Such formulae are combined across mass ranges that span different explosion mechanisms, potentialy introducing artificial features in the compact object mass distribution. Such artifacts impair the interpretation of gravitational-wave observations. We propose a “top-down” remnant mass prescription where we remove mass from the star for each possible mass-loss mechanism, instead of relying on the fallback onto a “proto-compact-object” to get the final mass. For one of these mass-loss mechanisms, we fit the metallicity-dependent mass lost to pulsational-pair instability supernovae from numerical simulations. By imposing no mass loss in the absence of pulses, our approach recovers the existing compact object masses at the low mass end and ensures continuity across the core-collapse/pulsational-pair-instability regime.

Carina High-contrast Imaging Project for massive Stars (CHIPS)

Context.Most massive stars belong to multiple systems, yet the formation process leading to such high multiplicity remains insufficiently understood. To help constrain the different formation scenarios that exist, insights into the low-mass end of the companion mass function of such stars is crucial. However, this is a challenging endeavour as (sub-)solar mass companions at angular separations (ρ) below 1″ (corresponding to 1000–3000 au in nearby young open clusters and OB associations) are difficult to detect due to the large brightness contrast with the central star.Aims.With the Carina High-contrast Imaging Project of massive Stars (CHIPS), we aim to obtain statistically significant constraints on the presence and properties of low-mass companions around massive stars in a previously unreachable observing window (Δmag ≳ 10 atρ ≲ 1″). In the second paper of the series, we focus on the Trumpler 14 cluster, which harbours some of the youngest and most massive O-type stars in the Milky Way.Methods.We obtained VLT-SPHERE observations of seven O-type objects in Trumpler 14 using IRDIFS_EXT mode. These provide us with a 12″ × 12″ field of view (approximately ((3 × 3)×104au) centred on each O star and allow us to search for companions at separations larger than 0″​​.15 (approx. 360 au) and down to magnitude contrast > 10 mag in the near-infrared. We used angular and spectral differential imaging along with Point Spread Function (PSF) fitting to detect sources and measure their flux relative to that of the central object. We then used grids of ATLAS9 and PHOENIX Local Thermodynamic Equilibrium (LTE) atmosphere models combined with (pre-)main-sequence evolutionary tracks to estimate the mass of the detected candidate companions.Results.We detected 211 sources with near-infrared magnitude contrast in the range of 2–12. Given the large surface number density of stars in Trumpler 14, one cannot reliably distinguish between cluster members and genuine companions for most of the detected sources. The closest companion, at only 0″​​.26, is characterised as a 1.4M⊙star with an age of 0.6 Myr, in excellent agreement with previous age estimates for Tr 14. The mass function peaks at about 0.4M⊙and presents a dearth of stars in the 0.5–0.8M⊙mass range compared to previous estimates of the initial mass function in Tr 14. While statistically significant, part of these differences may result from contamination of theK-band fluxes by circumstellar material.Conclusions.SPHERE is clearly suitable to probe the low-mass end of the mass function in the vicinity of massive stars. Follow-up SPHERE observations to obtain the full Y to K spectral energy distribution would allow for better constraints on the masses of the detected sources, and to confirm (or invalidate) the curious mass function that we derived for low-mass stars in the vicinity of the O-type objects in Trumpler 14.

The low-end of the black hole mass function at cosmic dawn

ABSTRACT Understanding the formation and growth of supermassive black holes (SMBHs) at high redshift represents a major challenge for theoretical models. In this work, we investigate the early evolution of the first SMBHs by constraining their distribution in mass and luminosity at z > 4. In particular, we focus on the poorly explored low-mass end of the nuclear black hole (BH) distribution down to z ≃ 4, and explore its connection with the nature of the first BH seeds and the processes governing their mass growth. To this aim, we have developed CAT (Cosmic Archaeology Tool), a new semi-analytic model that describes the formation of the first stars and BHs in a self-consistent way and follows the co-evolution of nuclear BHs and their host galaxies for a representative population at z > 4. We find that current observational constraints favour models where the growth of BH seeds is Eddington limited and occurs at the Bondi–Hoyle–Lyttleton rate or where super-Eddington accretion occurs via a slim disc during gas-rich galaxy mergers. The main difference between these two model variants lies at the low end of the predicted mass and luminosity functions at 4 ≤ z ≤ 6, where a clear gap appears in the first model, reflecting the stunted growth of light BH seeds formed as remnants of the first stars. Detecting this signature will be extremely challenging even for the future generation of space observatories, such as JWST, Athena, and Lynx.

The dark matter haloes of HI selected galaxies

ABSTRACT We present the neutral hydrogen mass (MHI) function (HIMF) and velocity width (w50) function (HIWF) based on a sample of 7857 galaxies from the 40 per cent data release of the ALFALFA survey (α.40). The low mass (velocity width) end of the HIMF (HIWF) is dominated by the blue population of galaxies whereas the red population dominates the HIMF (HIWF) at the high mass (velocity width) end. We use a deconvolution method to estimate the HI rotational velocity (Vrot) functions (HIVF) from the HIWF for the total, red, and blue samples. The HIWF and HIVF for the red and blue samples are well separated at the knee of the function compared to their HIMFs. We then use recent stacking results from the ALFALFA survey to constrain the halo mass (Mh) function of HI-selected galaxies. This allows us to obtain various scaling relations between MHI−w50−Vrot−Mh, which we present. The MHI−Mh relation has a steep slope ∼2.10 at small masses and flattens to ∼0.34 at masses larger than a transition halo mass, $\log _{10}(M_{\rm{ht}}h_{70}^2/M_{\odot })=10.62$. Our scaling relation is robust and consistent with a volume-limited sample of α.40. The MHI−Mh relation is qualitatively similar to the Mstar−Mh relation but the transition halo mass is smaller by ∼1.4 dex compared to that of the Mstar−Mh relation. Our results suggest that baryonic processes like heating and feedback in larger mass haloes suppress HI gas on a shorter time-scale compared to star formation.

Deep Extragalactic VIsible Legacy Survey (DEVILS): evolution of the σSFR–M⋆ relation and implications for self-regulated star formation

ABSTRACT We present the evolution of the star formation dispersion–stellar mass relation (σSFR–M⋆) in the DEVILS D10 region using new measurements derived using the ProSpect spectral energy distribution fitting code. We find that σSFR–M⋆ shows the characteristic ‘U-shape’ at intermediate stellar masses from 0.1 < z < 0.7 for a number of metrics, including using the deconvolved intrinsic dispersion. A physical interpretation of this relation is the combination of stochastic star formation and stellar feedback causing large scatter at low stellar masses and AGN feedback causing asymmetric scatter at high stellar masses. As such, the shape of this distribution and its evolution encodes detailed information about the astrophysical processes affecting star formation, feedback and the lifecycle of galaxies. We find that the stellar mass that the minimum σSFR occurs evolves linearly with redshift, moving to higher stellar masses with increasing lookback time and traces the turnover in the star-forming sequence. This minimum σSFR point is also found to occur at a fixed specific star formation rate (sSFR) at all epochs (sSFR ∼ 10−9.6 Gyr−1). The physical interpretation of this is that there exists a maximum sSFR at which galaxies can internally self-regulate on the tight sequence of star formation. At higher sSFRs, stochastic stellar processes begin to cause galaxies to be pushed both above and below the star-forming sequence leading to increased SFR dispersion. As the Universe evolves, a higher fraction of galaxies will drop below this sSFR threshold, causing the dispersion of the low stellar mass end of the star-forming sequence to decrease with time.

Stellar feedback in a clumpy galaxy at z ∼ 3.4

ABSTRACT Giant star-forming regions (clumps) are widespread features of galaxies at z ≈ 1−4. Theory predicts that they can play a crucial role in galaxy evolution, if they survive to stellar feedback for >50 Myr. Numerical simulations show that clumps’ survival depends on the stellar feedback recipes that are adopted. Up to date, observational constraints on both clumps’ outflows strength and gas removal time-scale are still uncertain. In this context, we study a line-emitting galaxy at redshift z ≃ 3.4 lensed by the foreground galaxy cluster Abell 2895. Four compact clumps with sizes ≲280 pc and representative of the low-mass end of clumps’ mass distribution (stellar masses ≲2 × 108 M⊙) dominate the galaxy morphology. The clumps are likely forming stars in a starbursting mode and have a young stellar population (∼10 Myr). The properties of the Lyman-α (Lyα) emission and nebular far-ultraviolet absorption lines indicate the presence of ejected material with global outflowing velocities of ∼200–300 km s−1. Assuming that the detected outflows are the consequence of star formation feedback, we infer an average mass loading factor (η) for the clumps of ∼1.8–2.4 consistent with results obtained from hydrodynamical simulations of clumpy galaxies that assume relatively strong stellar feedback. Assuming no gas inflows (semiclosed box model), the estimates of η suggest that the time-scale over which the outflows expel the molecular gas reservoir (≃7 × 108 M⊙) of the four detected low-mass clumps is ≲50 Myr.

The CARMENES search for exoplanets around M dwarfs

Context. Stellar activity poses one of the main obstacles for the detection and characterisation of small exoplanets around cool stars, as it can induce radial velocity (RV) signals that can hide or mimic the presence of planetary companions. Several indicators of stellar activity are routinely used to identify activity-related signals in RVs, but not all indicators trace exactly the same activity effects, nor are any of them always effective in all stars. Aims. We evaluate the performance of a set of spectroscopic activity indicators for M dwarf stars with different masses and activity levels with the aim of finding a relation between the indicators and stellar properties. Methods. In a sample of 98 M dwarfs observed with CARMENES, we analyse the temporal behaviour of RVs and nine spectroscopic activity indicators: cross-correlation function (CCF) full-width-at-half-maximum (FWHM), CCF contrast, CCF bisector inverse slope (BIS), RV chromatic index (CRX), differential line width (dLW), and indices of the chromospheric lines Hα and calcium infrared triplet. Results. A total of 56 stars of the initial sample show periodic signals related to activity in at least one of these ten parameters. RV is the parameter for which most of the targets show an activity-related signal. CRX and BIS are effective activity tracers for the most active stars in the sample, especially stars with a relatively high mass, while for less active stars, chromospheric lines perform best. FWHM and dLW show a similar behaviour in all mass and activity regimes, with the highest number of activity detections in the low-mass, high-activity regime. Most of the targets for which we cannot identify any activity-related signals are stars at the low-mass end of the sample (i.e. with the latest spectral types). These low-mass stars also show the lowest RV scatter, which indicates that ultracool M dwarfs could be better candidates for planet searches than earlier types, which show larger RV jitter. Conclusions. Our results show that the spectroscopic activity indicators analysed behave differently, depending on the mass and activity level of the target star. This underlines the importance of considering different indicators of stellar activity when studying the variability of RV measurements. Therefore, when assessing the origin of an RV signal, it is critical to take into account a large set of indicators, or at least the most effective ones considering the characteristics of the star, as failing to do so may lead to false planet claims.

A search for active galactic nuclei in low-mass compact galaxies

ABSTRACT Low-mass compact galaxies (ultracompact dwarfs, UCDs, and compact ellipticals, cEs) populate the stellar size–mass plane between globular clusters and early-type galaxies. Known to be formed either in situ with an intrinsically low mass or resulting from the stripping of a more massive galaxy, the presence of a supermassive or an intermediate-mass black hole (BH) could help discriminate between these possible scenarios. With this aim, we have performed a multiwavelength search of active BH activity, i.e. active galactic nuclei (AGN), in a sample of 937 low-mass compact galaxies (580 UCDs and 357 cEs). This constitutes the largest study of AGN activity in these types of galaxies. Based on their X-ray luminosity, radio luminosity, and morphology, and/or optical emission line diagnostic diagrams, we find a total of 11 cEs that host an AGN. We also study for the first time the location of both low-mass compact galaxies (UCDs and cEs) and dwarf galaxies hosting AGN on the BH–galaxy scaling relations, finding that low-mass compact galaxies tend to be overmassive in the BH mass–stellar mass plane but not as much in the BH mass–stellar velocity dispersion correlation. This, together with available BH mass measurements for some of the low-mass compact galaxies, supports a stripping origin for the majority of these objects that would contribute to the scatter seen at the low-mass end of the BH–galaxy scaling relations. However, the differences are too large to be explained solely by this scatter, and thus our results suggest that a flattening at such low masses is also plausible, happening at a velocity dispersion of ∼20–40 km s−1.

UniverseMachine: Predicting Galaxy Star Formation over Seven Decades of Halo Mass with Zoom-in Simulations

Abstract We apply the empirical galaxy–halo connection model UniverseMachine to dark-matter-only zoom-in simulations of isolated Milky Way (MW)–mass halos, along with their parent cosmological simulations. This application extends UniverseMachine predictions into the ultrafaint dwarf galaxy regime (102 M ⊙ ≤ M * ≤ 105 M ⊙) and yields a well-resolved stellar mass–halo mass (SMHM) relation over the peak halo mass range of 108–1015 M ⊙. The extensive dynamic range provided by the zoom-in simulations allows us to assess specific aspects of dwarf galaxy evolution predicted by UniverseMachine. In particular, although UniverseMachine is not constrained for dwarf galaxies with M * ≲ 108 M ⊙, our predicted SMHM relation is consistent with that inferred for MW satellite galaxies at z = 0 using abundance matching. However, UniverseMachine predicts that nearly all galaxies are actively star-forming below M * ∼ 107 M ⊙ and that these systems typically form more than half of their stars at z ≲ 4, which is discrepant with the star formation histories of Local Group dwarf galaxies that favor early quenching. This indicates that the current UniverseMachine model does not fully capture galaxy quenching physics at the low-mass end. We highlight specific improvements necessary to incorporate environmental and reionization-driven quenching for dwarf galaxies, and we provide a new tool to connect dark matter accretion to star formation over the full dynamic range that hosts galaxies.

The EBLM project – VIII. First results for M-dwarf mass, radius, and effective temperature measurements using CHEOPS light curves

ABSTRACT The accuracy of theoretical mass, radius, and effective temperature values for M-dwarf stars is an active topic of debate. Differences between observed and theoretical values have raised the possibility that current theoretical stellar structure and evolution models are inaccurate towards the low-mass end of the main sequence. To explore this issue, we use the CHEOPS satellite to obtain high-precision light curves of eclipsing binaries with low-mass stellar companions. We use these light curves combined with the spectroscopic orbit for the solar-type companion to measure the mass, radius, and effective temperature of the M-dwarf star. Here, we present the analysis of three eclipsing binaries. We use the pycheops data analysis software to fit the observed transit and eclipse events of each system. Two of our systems were also observed by the TESS satellite – we similarly analyse these light curves for comparison. We find consistent results between CHEOPS and TESS, presenting three stellar radii and two stellar effective temperature values of low-mass stellar objects. These initial results from our on-going observing programme with CHEOPS show that we can expect to have ∼24 new mass, radius, and effective temperature measurements for very low-mass stars within the next few years.

Mild radial variations of the stellar IMF in the bulge of M31

ABSTRACT Using new, homogeneous, long-slit spectroscopy in the wavelength range from ∼0.35 to $\sim 1 \, \mu$m, we study radial gradients of optical and near-infrared (NIR) initial mass function (IMF)-sensitive features along the major axis of the bulge of M31, out to a galactocentric distance of ∼200 arcsec (∼800 pc). Based on state-of-the-art stellar population synthesis models with varying Na abundance ratio, we fit a number of spectral indices, from different chemical species (including TiO’s, Ca, and Na indices), to constrain the low-mass (≲0.5 M⊙) end slope (i.e. the fraction of low-mass stars) of the stellar IMF, as a function of galactocentric distance. Outside a radial distance of ∼10 arcsec, we infer an IMF similar to a Milky Way-like distribution, while at small galactocentric distances, an IMF radial gradient is detected, with a mildly bottom-heavy IMF in the few inner arcsec. We are able to fit Na features (both NaD and $\rm Na\,{\small I}8190$), without requiring extremely high Na abundance ratios. $\rm [Na/Fe]$ is ∼0.4 dex for most of the bulge, rising up to ∼0.6 dex in the innermost radial bins. Our results imply an overall, luminosity-weighted, IMF and mass-to-light ratio for the M31 bulge, consistent with those for a Milky Way-like distribution, in contrast to results obtained, in general, for most massive early-type galaxies.

Testing the role of environmental effects on the initial mass function of low-mass stars

ABSTRACT In the star-formation process, the vital impact of environmental factors such as feedback from massive stars and stellar density on the form of the initial mass function (IMF) at the low-mass end is yet to be understood. Hence a systematic highly sensitive observational analysis of a sample of regions under diverse environmental conditions is essential. We analyse the IMF of eight young clusters (<5 Myr), namely IC 1848-West, IC 1848-East, NGC 1893, NGC 2244, NGC 2362, NGC 6611, Stock 8, and Cygnus OB2, which are located at the Galactocentric distance (Rg) range ∼6–12 kpc along with the nearby cluster IC 348 using deep near-IR photometry and Gaia DR2. These clusters are embedded in massive stellar environments of radiation strength $\log(L_\mathrm{FUV}/\mathrm{L}_{\odot })\, \sim 2.6$–6.8, $\log(L_\mathrm{EUV})\, \sim$ 42.2–50.85 photon s−1, with stellar density in the range of ∼170–1220 star pc−2. After structural analysis and field decontamination we obtain an unbiased uniformly sensitive sample of pre-main-sequence members of the clusters down to the brown-dwarf regime. The lognormal fit to the IMF of nine clusters gives the mean characteristic mass (mc) and σ of 0.32 ± 0.02 M⊙ and 0.47 ± 0.02, respectively. We compare the IMF with that of low- and high-mass clusters across the Milky Way. We also check for any systematic variation with respect to the radiation field strength and the stellar density as well with Rg. We conclude that there is no strong evidence for an environmental effect in the underlying form of the IMF of these clusters.

A preserved high-z compact progenitor in the heart of NGC 3311 revealed with MUSE 2D stellar population analysis

Context. Massive early-type galaxies are believed to be the end result of an extended mass accretion history. The stars formed in situ very early on in the initial phase of the mass assembly might have originated from an extremely intense and rapid burst of star formation. These stars may still be found within the cores of such galaxies at z = 0, depending on their accretion and merger histories. Aims. We wish to investigate the presence of a surviving high-z compact progenitor component in the brightest galaxy of the Hydra I cluster, NGC 3311, by mapping its 2D kinematics and stellar population out to 2 effective radii. Our goal is to understand the formation of its several structural components and trace their mass assembly back in time. Methods. We combined MUSE observations, a customized and extended version of the state-of-the-art EMILES single stellar population models, and a newly developed parametric fully Bayesian framework to model the observed spectra using full-spectrum fitting. Results. We present 2D maps and radial profiles of the stellar velocity dispersion, age, total metallicity, α-element, sodium abundance ([Na/Fe]), and the initial mass function (IMF) slope. All properties have significant gradients, confirming the existence of multiple structural components, also including a “blue spot” characterized by younger and more metal-rich stars. We find that the component dominating the light budget of NGC 3311 within R ≲ 2.0 kpc is the surviving z = 0 analog of a high-z compact core. This concentrated structure has a relatively small velocity dispersion (σ* ≈ 180 km s−1), is very old (ages ≳ 11 Gyr), metal-rich ([Z/H] ∼ 0.2 and [Na/Fe] ∼ 0.4), and has a bottom-heavy IMF (with slope Γb ∼ 2.4). In the outer region, instead, the line-of-sight velocity distribution becomes increasingly broad, and the stars are younger. They are also more metal and sodium poor but are richer in α-elements. The low-mass end of the IMF slope becomes Chabrier-like with increasing galactocentric distance. Conclusions. The existence of multiple structural components in NGC 3311 from photometry, kinematics, and stellar populations confirms the predictions from the two-phase formation scenario for NGC 3311, according to which a first very short, high-z star-formation episode formed a compact stellar structure in its core, which then grew in size by the extended mass assembly of relatively massive satellites. Interestingly, the outer stellar population has an overabundant [α/Fe], most likely because NGC 3311, located at the center of the galaxy cluster, accreted stars from rapidly quenched satellites.

On the origin of the mass–metallicity gradient relation in the local Universe

ABSTRACT In addition to the well-known gas phase mass–metallicity relation (MZR), recent spatially resolved observations have shown that local galaxies also obey a mass–metallicity gradient relation (MZGR), whereby metallicity gradients can vary systematically with galaxy mass. In this work, we use our recently developed analytic model for metallicity distributions in galactic discs, which includes a wide range of physical processes – radial advection, metal diffusion, cosmological accretion, and metal-enriched outflows – to simultaneously analyse the MZR and MZGR. We show that the same physical principles govern the shape of both: centrally peaked metal production favours steeper gradients, and this steepening is diluted by the addition of metal-poor gas, which is supplied by inward advection for low-mass galaxies and by cosmological accretion for massive galaxies. The MZR and the MZGR both bend at galaxy stellar mass $\sim 10^{10}{-}10^{10.5}\, \rm {M_{\odot }}$, and we show that this feature corresponds to the transition of galaxies from the advection-dominated to the accretion-dominated regime. We also find that both the MZR and MZGR strongly suggest that low-mass galaxies preferentially lose metals entrained in their galactic winds. While this metal-enrichment of the galactic outflows is crucial for reproducing both the MZR and the MZGR at the low-mass end, we show that the flattening of gradients in massive galaxies is expected regardless of the nature of their winds.

Deep extragalactic visible legacy survey (DEVILS): stellar mass growth by morphological type since z = 1

ABSTRACT Using high-resolution Hubble Space Telescope imaging data, we perform a visual morphological classification of ∼36 000 galaxies at z < 1 in the deep extragalactic visible legacy survey/cosmological evolution survey region. As the main goal of this study, we derive the stellar mass function (SMF) and stellar mass density (SMD) sub-divided by morphological types. We find that visual morphological classification using optical imaging is increasingly difficult at z > 1 as the fraction of irregular galaxies and merger systems (when observed at rest-frame UV/blue wavelengths) dramatically increases. We determine that roughly two-thirds of the total stellar mass of the Universe today was in place by z ∼ 1. Double-component galaxies dominate the SMD at all epochs and increase in their contribution to the stellar mass budget to the present day. Elliptical galaxies are the second most dominant morphological type and increase their SMD by ∼2.5 times, while by contrast, the pure-disc population significantly decreases by $\sim 85{{\ \rm per\ cent}}$. According to the evolution of both high- and low-mass ends of the SMF, we find that mergers and in situ evolution in discs are both present at z < 1, and conclude that double-component galaxies are predominantly being built by the in situ evolution in discs (apparent as the growth of the low-mass end with time), while mergers are likely responsible for the growth of ellipticals (apparent as the increase of intermediate/high-mass end).

A little FABLE: exploring AGN feedback in dwarf galaxies with cosmological simulations

ABSTRACT Contrary to the standard lore, there is mounting observational evidence that feedback from active galactic nuclei (AGN) may also play a role at the low-mass end of the galaxy population. We investigate this using the cosmological simulation suite fable, with a particular focus on the dwarf regime (Mstellar < 109.5 M⊙). We find that overmassive black holes (BHs), with respect to the mean scaling relations with their host galaxies, drive hotter and faster outflows and lead to significantly reduced gas mass fractions. They are also more likely to display a kinematically misaligned ionized gas component in our mock MaNGA velocity maps, although we caution that cosmic inflows and mergers contribute to misalignments as well. While in the local Universe the majority of AGN in dwarfs are much dimmer than the stellar component, for z ≥ 2 there is a significant population that outshines their hosts. These high-redshift overmassive BHs contribute to the quenching of dwarfs, whereas at late cosmic times supernova (SN) feedback is more efficient. While our results are overall in good agreement with X-ray observations of AGN in dwarfs, the lack of high-luminosity X-ray AGN in fable at low redshifts highlights an interesting possibility that SN feedback could be too strong in fable’s dwarfs, curtailing AGN growth and feedback. We predict that future observations may uncover many more AGN in dwarfs with lower luminosities and at higher redshifts.

Discovery of soft and hard X-ray time lags in low-mass AGNs

ABSTRACT The scaling relations between the black hole (BH) mass and soft lag properties for both active galactic nuclei (AGNs) and BH X-ray binaries (BHXRBs) suggest the same underlying physical mechanism at work in accreting BH systems spanning a broad range of mass. However, the low-mass end of AGNs has never been explored in detail. In this work, we extend the existing scaling relations to lower mass AGNs, which serve as anchors between the normal-mass AGNs and BHXRBs. For this purpose, we construct a sample of low-mass AGNs ($M_{\rm BH}\lt 3\times 10^{6} \, \mathrm{M}_{\odot }$) from the XMM–Newton archive and measure frequency-resolved time-delays between the soft (0.3–1 keV) and hard (1–4 keV) X-ray emissions. We report that the soft band lags behind the hard band emission at high frequencies ∼[1.3−2.6] × 10−3 Hz, which is interpreted as a sign of reverberation from the inner accretion disc in response to the direct coronal emission. At low frequencies (∼[3−8] × 10−4 Hz), the hard-band lags behind the soft-band variations, which we explain in the context of the inward propagation of luminosity fluctuations through the corona. Assuming a lamppost geometry for the corona, we find that the X-ray source of the sample extends at an average height and radius of ∼10rg and ∼6rg, respectively. Our results confirm that the scaling relations between the BH mass and soft lag amplitude/frequency derived for higher mass AGNs can safely extrapolate to lower mass AGNs, and the accretion process is indeed independent of the BH mass.