Galaxy Formation Research Articles

The host of GRB 171205A in 3D

Long GRB hosts at z < 1 are usually low-mass, low-metallicity star-forming galaxies. Here we present the most detailed, spatially resolved study of the host of GRB 171205A so far, a grand-design barred spiral galaxy at z = 0.036. Our analysis includes MUSE integral field spectroscopy complemented with high-spatial-resolution UV/VIS HST imaging and CO(1−0) and H I 21 cm data. The GRB is located in a small star-forming region in a spiral arm of the galaxy at a deprojected distance of ∼8 kpc from the center. The galaxy shows a smooth negative metallicity gradient and the metallicity at the GRB site is half solar, slightly below the mean metallicity at the corresponding distance from the center. Star formation in this galaxy is concentrated in a few H II regions between 5 and 7 kpc from the center and at the end of the bar, inwards from the GRB region; however the H II region hosting the GRB is in the top 10% of the regions with the highest specific star-formation rate. The stellar population at the GRB site has a very young component (< 5 Myr) that contributes a significant part of the light. Ionized and molecular gas show only minor deviations at the end of the bar. A parallel study found an asymmetric H I distribution and some additional gas near the position of the GRB, which might explain the star-forming region of the GRB site. Our study shows that long GRBs can occur in many types of star-forming galaxies; however the actual GRB sites have consistently low metallicity, high star formation rates, and a young population. Furthermore, gas inflow or interactions triggering the star formation producing the GRB progenitor might not be evident in ionized or even molecular gas but only in H I.

Exploring the Origin of Cold Gas and Star Formation in a Rare Population of Strongly Bulge-dominated Early-type Galaxies

We analyze the properties of a rare population, the strongly bulge-dominated early-type galaxies (sBDEs) with significant H i gas, using the databases from the FAST All Sky H i survey (FASHI) and the Arecibo Legacy Fast ALFA (ALFALFA) survey. We select the sBDEs from the Sloan Digital Sky Survey and crossmatch with the FASHI-ALFALFA combined H i sample, resulting in 104 H i-rich sBDEs. These sBDEs tend to have extremely high H i reservoirs, which is rare in previous studies such as ATLAS3D. A total of 70% of the selected sBDEs are classified as quiescent galaxies, even though they have a large H i reservoir. We study the properties of these sBDEs from four main aspects: stellar population, gas-phase metallicity, stacked H i spectra, and environment. The majority of H i-rich sBDEs appear to show lower gas-phase metallicity and are located in significantly lower-density environments, suggesting an external origin for their H i gas. We find that star-forming sBDEs exhibit statistically higher star formation efficiency and slightly older stellar populations compared to normal star-forming galaxies, suggesting a recent star formation on the Gyr timescale. They also show narrower and more concentrated H i profiles compared to control star-forming galaxies, which may explain their higher star formation efficiency.

The ALMA-CRISTAL survey

We present the morphological parameters and global properties of dust-obscured star formation in typical star-forming galaxies at z = 4–6. Among 26 galaxies composed of 20 galaxies observed by the Cycle-8 ALMA Large Program, CRISTAL, and 6 galaxies from archival data, we individually detect rest-frame 158 μm dust continuum emission from 19 galaxies, 9 of which are reported for the first time. The derived far-infrared luminosities are in the range log10LIR [L⊙] = 10.9 − 12.4, an order of magnitude lower than previously detected massive dusty star-forming galaxies (DSFGs). We find the average relationship between the fraction of dust-obscured star formation (fobs) and the stellar mass to be consistent with previous results at z = 4–6 in a mass range of log10M* [M⊙]∼9.5 − 11.0 and to show potential evolution from z = 6 − 9. The individual fobs exhibits significant diversity, and we find a potential correlation with the spatial offset between the dust and UV continuum, suggesting that inhomogeneous dust reddening may cause the source-to-source scatter in fobs. The effective radii of the dust emission are on average ∼1.5 kpc and are about two times more extended than those seen in rest-frame UV. The infrared surface densities of these galaxies (ΣIR ∼ 2.0 × 1010 L⊙ kpc−2) are one order of magnitude lower than those of DSFGs that host compact central starbursts. On the basis of the comparable contribution of dust-obscured and dust-unobscured star formation along with their similar spatial extent, we suggest that typical star-forming galaxies at z = 4 − 6 form stars throughout the entirety of their disks.

MICONIC: JWST/MIRI MRS observations of the nuclear and circumnuclear regions of Mrk 231

We present JWST/MIRI MRS spatially resolved ∼5 − 28 μm observations of the central ∼4 − 8 kpc of the ultraluminous infrared galaxy and broad absorption line quasar Mrk 231. These are part of the Mid-Infrared Characterization of Nearby Iconic galaxy Centers (MICONIC) program of the MIRI European Consortium guaranteed time observations. No high excitation lines (i.e., [Mg V] at 5.61 μm or [Ne V] at 14.32 μm) typically associated with the presence of an active galactic nucleus (AGN) are detected in the nuclear region of Mrk 231. This is likely due to the intrinsically X-ray weak nature of its quasar. Some intermediate ionization potential lines, for instance, [Ar III] at 8.99 μm and [S IV] at 10.51 μm, are not detected either, even though they are clearly observed in a star-forming region ∼920 pc south-east of the AGN. Thus, the strong nuclear mid-infrared (mid-IR) continuum is also in part hampering the detection of faint lines in the nuclear region. The nuclear [Ne III]/[Ne II] line ratio is consistent with values observed in star-forming galaxies. Moreover, we resolve for the first time the nuclear starburst in the mid-IR low-excitation line emission (size of ∼400 pc, FWHM). Several pieces of evidence also indicate that it is partly obscured even at these wavelengths. At the AGN position, the ionized and warm molecular gas emission lines have modest widths (W80 ∼ 300 km s−1). There are, however, weak blueshifted wings reaching velocities v02 ≃ − 400 km s−1 in [Ne II]. The nuclear starburst is at the center of a large (∼8 kpc), massive rotating disk with widely-spread, low velocity outflows. Given the high star formation rate of Mrk 231, we speculate that part of the nuclear outflows and the large-scale non-circular motions observed in the mid-IR are driven by its powerful nuclear starburst.

The ALMA-CRISTAL survey: Dust temperature and physical conditions of the interstellar medium in a typical galaxy at z=5.66

We present new $ rest =77$ mu m dust continuum observations from the Atacama Large Millimeter/submillimeter Array of HZ10 (CRISTAL-22). This dusty main sequence galaxy at $z$=5.66 was observed as part of the Resolved Ism in STar-forming Alma Large program (CRISTAL). The high angular resolution of the ALMA Band 7 and new Band 9 data ($ ''.4$) reveals the complex structure of HZ10, which comprises two main components (HZ10-C and HZ10-W), along with a bridge-like dusty emission between them (i.e., "the bridge"). Using a modified blackbody function to model the dust spectral energy distribution (SED), we constrained the physical conditions of the interstellar medium (ISM) and its variations among the different components identified in HZ10. We find that HZ10-W (the more UV-obscured component) has an SED dust temperature of $T_ SED $sim 51.2$ K; this was found to be sim 5 K higher (which is statistically insignificant; i.e., less than 1sigma ) than that of the central component and previous global estimations for HZ10. Our new ALMA data allow us to reduce the uncertainties of global SED $ measurements by a factor of sim 2.3, compared to previous studies. The HZ10 components have (FIR) luminosity ratios and FIR surface densities values that are consistent with local starburst galaxies. However, HZ10-W shows a lower CII /FIR ratio compared to the other two components (albeit still within the uncertainties), which may suggest a harder radiation field destroying polycyclic aromatic hydrocarbon associated with emission (e.g., active galactic nuclei or young stellar populations). While HZ10-C appears to follow the tight IRX-$ UV $ relation seen in local UV-selected starburst galaxies and high-$z$ star-forming galaxies, we find that both HZ10-W and the bridge depart from this relation and are well described by dust-screen models with holes in front of a hard UV radiation field. This suggests that the UV emission, which is likely coming from young stellar populations, is strongly attenuated in the "dustier" components of the HZ10 system.

Physical conditions in Centaurus A’s northern filaments

We present the first observations of HCO+(1–0) and HCN(1–0) emission in the northern filaments of Centaurus A with ALMA. HCO+(1–0) is detected in nine clumps of the Horseshoe complex, with similar velocities as the CO(1–0) emission. Conversely, HCN(1–0) is not detected, and we derive upper limits on the flux. At a resolution of ∼40 pc, the line ratio of the velocity-integrated intensities IHCO+/ICO varies between 0.03 and 0.08, while IHCO+/IHCN is higher than unity, with an average lower limit of 1.51. These ratios are significantly higher than what is observed in nearby star-forming galaxies. Moreover, the ratio IHCO+/ICO decreases with increasing CO-integrated intensity, contrary to what is observed in the star-forming galaxies. This indicates that the HCO+ emission is enhanced and may not arise from dense gas within the Horseshoe complex. This hypothesis is strengthened by the average line ratio IHCN/ICO < 0.03, which suggests that the gas density is rather low. Using non-local thermal equilibrium, large velocity gradient modelling with RADEX, we explored two possible phases of the gas, which we call ‘diffuse’ and ‘dense’ and are characterised by a significant difference in the HCO+ abundance relative to CO, respectively NHCO+/NCO = 10−3 and NHCO+/NCO = 3 × 10−5. The average CO(1–0) and HCO+(1–0) integrated intensities and the upper limit on HCN(1–0) are compatible with both diffuse (nH = 103 cm−3, Tkin = 15 − 165 K) and dense gas (nH = 104 cm−3, Tkin > 65 K). The spectral setup of the present observations also covers SiO(2–1). While undetected, the upper limit on SiO(2–1) is not compatible with the RADEX predictions for the dense gas. We conclude that the nine molecular clouds detected in HCO+(1–0) are likely dominated by diffuse molecular gas. While the exact origin of the HCO+(1–0) emission remains to be investigated, it is likely related to the energy injection within the molecular gas that prevents gravitational collapse and star formation.

Linking High-z and Low-z: Are We Observing the Progenitors of the Milky Way with JWST?

The recent JWST observation of the Firefly Sparkle at z = 8.3 offers a unique opportunity to link the high- and the low-z Universe. Indeed, the claim of it being a Milky Way (MW) type of assembly at the cosmic dawn opens the possibility of interpreting the observation with locally calibrated galaxy-formation models. Here, we use the a state-of-the-art MW-evolution model to perform forward modeling of our Galaxy's progenitors at high-z. We build a set of mock spectra for the MW building blocks to make predictions for JWST and to interpret the Firefly Sparkle observation. First, we find that the most massive MW progenitor becomes detectable in a deep survey like JADES from z ≈ 8.2, meaning that we could have already observed MW analogs that still need interpretation. Second, we provide predictions for the number of detectable MW progenitors in lensed surveys like the CAnadian NIRISS Unbiased Cluster Survey, and interpret the Firefly Sparkle as a group of MW building blocks. Both the number of detections and the observed NIRCam photometry are consistent with our predictions. By identifying the MW progenitors whose mock photometry best fits the data, we find bursty and extended star formation histories, lasting >150–300 Myr, and estimate their properties: M h ≈ 108−9 M ⊙, M ⋆ ≈ 106.2−7.5 M ⊙, SFR ≈ 0.04–0.20 M ⊙ yr−1, and Z gas ≈ 0.04–0.24 Z ⊙. Uncovering the properties of MW analogs at cosmic dawn by combining JWST observations and locally constrained models will allow us to understand our Galaxy's formation, linking the high- and low-z perspectives.

Two-process Model and Residual Abundance Analysis of the Milky Way Massive Satellites

The “two-process model” is a promising technique for interpreting stellar chemical abundance data from large-scale surveys (e.g., the Sloan Digital Sky Survey IV/V and the Galactic Archeology with HERMES survey), enabling more quantitative empirical studies of differences in chemical enrichment history between galaxies without relying on detailed yield and evolution models. In this work, we fit two-process model parameters to (1) a luminous giant Milky Way (MW) sample and (2) stars comprising the Sagittarius dwarf galaxy (Sgr). We then use these two sets of model parameters to predict the abundances of 14 elements of stars belonging to the MW and in five of its massive satellite galaxies, analyzing the residuals between the predicted and observed abundances. We find that the model fit to (1) results in large residuals (0.1–0.3 dex) for most metallicity-dependent elements in the metal-rich ([Mg/H] > −0.8) stars of the satellite galaxies. However, the model fit to (2) results in small or no residuals for all elements across all satellite galaxies. Therefore, despite the wide variation in [X/Mg]–[Mg/H] abundance patterns of the satellite galaxies, the two-process framework provides an accurate characterization of their abundance patterns across many elements, but these multielement patterns are systematically different between the dwarf galaxy satellites and the MW disks. We consider a variety of scenarios for the origin of this difference, highlighting the possibility that a large inflow of pristine gas to the MW disk diluted the metallicity of star-forming gas without changing abundance ratios.

MEASNet. I. A Model for Barium Star Identification and s-process Abundance Estimation from LAMOST DR10 Low-resolution Survey

Barium stars are peculiar stars with enhanced slow neutron capture process (s-process) elements. Abundance analysis of them aids in better understanding the chemical evolution of the Milky Way. In this paper, we introduce a data-driven method named the memory-enhanced adaptive spectral network (MEASNet) to search for barium candidates in the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) low-resolution survey (LRS) and estimate the abundance of five s-process elements: Sr, Y, Ba, Ce, and Nd. MEASNet, trained using spectra from common stars in both LAMOST and the Galactic Archaeology with HERMES survey, showcases notable performance: for the classification task, precision = 98.22% and recall = 94.12%; in prediction, the mean absolute error for the seven elements range between 0.07 and 0.15 dex. After training, we apply the model to 4,083,003 stellar spectra from LAMOST DR10 LRS, successfully identifying 1,803,670 spectra of barium candidates ([Ba/Fe] ≥ 0.25 dex) along with their five s-process elemental abundances. The catalog enlarges the sample size, providing a wealth of data for further statistical analysis of the formation and evolution of barium stars. Meanwhile, this work highlights the potential value of MEASNet in star classification and abundance estimation, offering a strong reference for future data-driven models.

Retrieval of the physical parameters of galaxies from WEAVE-StePS-like data using machine learning

Context. The William Herschel Telescope Enhanced Area Velocity Explorer (WEAVE) is a new, massively multiplexing spectrograph that allows us to collect about one thousand spectra over a 3 square degree field in one observation. The WEAVE Stellar Population Survey (WEAVE-StePS) in the next 5 years will exploit this new instrument to obtain high-S/N spectra for a magnitude-limited (IAB = 20.5) sample of ∼25 000 galaxies at moderate redshifts (z ≥ 0.3), providing insights into galaxy evolution in this as yet unexplored redshift range. Aims. We aim to test novel techniques for retrieving the key physical parameters of galaxies from WEAVE-StePS spectra using both photometric and spectroscopic (spectral indices) information for a range of noise levels and redshift values. Methods. We simulated ∼105 000 galaxy spectra assuming star formation histories with an exponentially declining star formation rate, covering a wide range of ages, stellar metallicities, specific star formation rates (sSFRs), and dust extinction values. We considered three redshifts (i.e. z = 0.3, 0.55, and 0.7), covering the redshift range that WEAVE-StePS will observe. We then evaluated the ability of the random forest and K-nearest neighbour algorithms to correctly predict the average age, metallicity, sSFR, dust attenuation, and time since the bulk of formation, assuming no measurement errors. We also checked how much the predictive ability deteriorates for different noise levels, with S/NI,obs = 10, 20, and 30, and at different redshifts. Finally, the retrieved sSFR was used to classify galaxies as part of the blue cloud, green valley, or red sequence. Results. We find that both the random forest and K-nearest neighbour algorithms accurately estimate the mass-weighted ages, u-band-weighted ages, and metallicities with low bias. The dispersion varies from 0.08–0.16 dex for age and 0.11–0.25 dex for metallicity, depending on the redshift and noise level. For dust attenuation, we find a similarly low bias and dispersion. For the sSFR, we find a very good constraining power for star-forming galaxies, log sSFR ≳ −11, where the bias is ∼0.01 dex and the dispersion is ∼0.10 dex. However, for more quiescent galaxies, with log sSFR ≲ −11, we find a higher bias, ranging from 0.61 to 0.86 dex, and a higher dispersion, ∼0.4 dex, depending on the noise level and redshift. In general, we find that the random forest algorithm outperforms the K-nearest neighbours. Finally, we find that the classification of galaxies as members of the green valley is successful across the different redshifts and S/Ns. Conclusions. We demonstrate that machine learning algorithms can accurately estimate the physical parameters of simulated galaxies for a WEAVE-StePS-like dataset, even at relatively low S/NI, obs = 10 per Å spectra with available ancillary photometric information. A more traditional approach, Bayesian inference, yields comparable results. The main advantage of using a machine learning algorithm is that, once trained, it requires considerably less time than other methods.

Exploring the connection between AGN radiative feedback and massive black hole spin

We present a novel implementation for active galactic nucleus (AGN) feedback through ultrafast winds in the code GIZMO. Our feedback recipe accounts for the angular dependence of radiative feedback on black hole spin. We self-consistently evolve in time (i) the gas-accretion process from resolved scales to a smaller scale unresolved (subgrid) AGN disk, (ii) the evolution of the spin of the massive black hole (MBH), (iii) the injection of AGN-driven winds into the resolved scales, and (iv) the spin-induced anisotropy of the overall feedback process. We tested our implementation by following the propagation of the wind-driven outflow into an homogeneous medium, and here we present a comparison of the results against simple analytical models. We also considered an isolated galaxy setup, where the galaxy is thought to be formed from the collapse of a spinning gaseous halo, and there we studied the impact of the AGN feedback on the evolution of the MBH and of the host galaxy. We find that: (i) AGN feedback limits the gas inflow that powers the MBH, with a consequent weak impact on the host galaxy characterized by a suppression of star formation by about a factor of two in the nuclear (≲kpc) region; (ii) the impact of AGN feedback on the host galaxy and on MBH growth is primarily determined by the AGN luminosity rather than by its angular pattern set by the MBH spin (i.e., more luminous AGNs more efficiently suppress central star formation (SF), clearing wider central cavities and driving outflows with larger semiopening angles); (iii) the imprint of the angular pattern of AGN radiation emission is detected more clearly at high (i.e., Eddington) accretion rates. At such high rates, the more isotropic angular patterns, as occur for high spin values, sweep away gas in the nuclear region more easily, therefore causing a slower MBH mass and spin growths and a higher quenching of SF. We argue that the influence of spin-dependent anisotropy of AGN feedback on MBH and galaxy evolution is likely to be relevant in those scenarios characterized by high and prolonged MBH accretion episodes and by high AGN wind–galaxy coupling. Such conditions are more frequently met in galaxy mergers and/or high-redshift galaxies.

Spatially–resolved gas-phase metallicity in Seyfert galaxies

Abstract We explore the relations between the gas-phase metallicity radial profiles (few hundred inner parsec) and multiple galaxy properties for 15 Seyfert galaxies from the AGNIFS (Active Galactic Nuclei Integral Field Spectroscopy) sample using optical Integral Field Unit (IFU) observations from Gemini Multi-Object Spectrographs (GMOS) and Multi Unit Spectroscopic Explorer (MUSE) processed archival data. The data were selected at z ≲ 0.013 within black hole mass range [6 < log (MBH/M⊙) < 9] with moderate 14–150 keV X-ray luminosities $\left[42\, \lesssim \, \log L_X (\rm erg\, s^{-1})\, \lesssim \, 44\right]$. We estimated the gas-phase metallicity using the strong-line methods and found mean values for the oxygen dependent (Z ∼ 0.75 Z⊙) and nitrogen dependent (Z ∼ 1.14 Z⊙) calibrations. These estimates show excellent agreement with ΔZ ≈ 0.19 dex and ΔZ ≈ 0.18 dex between the mean values from the two strong-line calibrations for GMOS and MUSE respectively, consistent with the order of metallicity uncertainty via the strong-line methods. We contend that our findings align with a scenario wherein local Seyferts have undergone seamless gas accretion histories, resulting in positive metallicity profile over an extended period of time, thereby providing insights into galaxy evolution and the chemical enrichment or depletion of the universe. Additionally, we argue that metal-poor gas inflow from the local interstellar medium (ISM) and accreted through the circumgalactic medium (CGM) onto the galaxy systems regulates the star formation processes by diluting their central metallicity and inverting their metallicity gradients, producing a more prominent anti-correlation between gas-phase metallicity and Eddington ratio.

The hot circumgalactic medium in the eROSITA All-Sky Survey

Context. The circumgalactic medium (CGM) provides the material needed for galaxy formation and influences galaxy evolution. The hot (T > 106K) CGM is poorly detected around galaxies with stellar masses (M*) lower than 3 × 1011 M⊙ due to the low surface brightness. Aims. We aim to detect the X-ray emission from the hot CGM around Milky Way-mass (MW-mass, log(M*/M⊙) = 10.5 − 11.0) and M31-mass (log(M*/M⊙) = 11.0 − 11.25) galaxies, in addition to measuring the X-ray surface brightness profile of the hot CGM. Methods. We applied a stacking technique to gain enough statistics to detect the hot CGM. We used the X-ray data from the first four SRG/eROSITA All-Sky Surveys (eRASS:4). We discussed how the satellite galaxies could bias the stacking and the method we used to carefully build the central galaxy samples. Based on the SDSS spectroscopic survey and halo-based group finder algorithm, we selected central galaxies with spectroscopic redshifts of zspec < 0.2 and stellar masses of 10.0 < log(M*/M⊙) < 11.5 (85 222 galaxies) – or halo masses of 11.5 < log(M200m/M⊙) < 14.0 (125,512 galaxies). By stacking the X-ray emission around galaxies, we obtained the mean X-ray surface brightness profiles. We masked the detected X-ray point sources and carefully modeled the X-ray emission from the unresolved active galactic nuclei (AGN) and X-ray binaries (XRB) to obtain the X-ray emission from the hot CGM. Results. We measured the X-ray surface brightness profiles for central galaxies of log(M*/M⊙) > 10.0 or log(M200m/M⊙) > 11.5. We detected the X-ray emission around MW-mass and more massive central galaxies extending up to the virial radius (Rvir). The signal-to-noise ratio (S/N) of the extended emission around MW-mass (M31-mass) galaxy is about 3.1σ (4.7σ) within Rvir. We used a β model to describe the X-ray surface brightness profile of the hot CGM (SX, CGM). We obtained a central surface brightness of log(SX,0[erg s−1 kpc−2]) = 36.7−0.4+1.4 (37.1−0.4+1.5) and β = 0.43−0.06+0.10 (0.37−0.02+0.04) for MW-mass (M31-mass) galaxies. For galaxies with log(M200m/M⊙) > 12.5, the extended X-ray emission is detected with S/N > 2.8σ and the SX, CGM can be described by a β model with β ≈ 0.4 and log(SX,0[erg s−1 kpc−2]) > 37.2. We estimated the baryon budget of the hot CGM and obtained a value that is lower than the prediction of ΛCDM cosmology, indicating significant gas depletion in these halos. We extrapolated the hot CGM profile measured within Rvir to larger radii and found that within ≈3Rvir, the baryon budget is close to the ΛCDM cosmology prediction. Conclusions. We measured the extended X-ray emission from representative populations of central galaxies around and above MW-mass out to Rvir. Our results set a firm footing for the presence of the hot CGM around such galaxies. These measurements constitute a new benchmark for galaxy evolution models and possible implementations of feedback processes therein.

Ubiquitous radio emission in quasars: predominant AGN origin and a connection to jets, dust, and winds

We present a comprehensive study of the physical origin of radio emission in optical quasars at redshifts $z<2.5$. We focus particularly on the associations between compact radio emission, dust reddening, and outflows identified in our earlier work. Leveraging the deepest low-frequency radio data available to date (LoTSS Deep DR1), we achieve radio detection fractions of up to 94<!PCT!>, demonstrating the virtual ubiquity of radio emission in quasars, and a continuous distribution in radio loudness. Through our analysis of radio properties, combined with spectral energy distribution modelling of deep multiwavelength photometry, we establish that the primary source of radio emission in quasars is the active galactic nucleus (AGN), rather than star formation. Modelling the dust reddening of the accretion disc emission shows a continuous increase in radio detection in quasars as a function of the reddening parameter E(B-V) suggesting a causal link between radio emission and dust reddening. Confirming previous findings, we observe that the radio excess in red quasars is most pronounced for sources with compact radio morphologies and intermediate radio loudness. We find a significant increase in Oiii and Civ outflow velocities for red quasars not seen in our control sample, with particularly powerful Oiii winds in those around the threshold from radio-quiet to radio-loud. Based on the combined characterisation of radio, reddening, and outflow properties in our sample, we favour a model in which the compact radio emission observed in quasars originates in compact radio jets and their interaction with a dusty, circumnuclear environment. In particular, our results align with the theory that jet-induced winds and shocks resulting from this interaction are the origin of the enhanced radio emission in red quasars. Further investigation of this model is crucial for advancing our understanding of quasar feedback mechanisms and their role in galaxy evolution.

Disentangling the Faraday rotation sky

Context. Magnetic fields permeate the diffuse interstellar medium (ISM) of the Milky Way, and are essential to explain the dynamical evolution and current shape of the Galaxy. Magnetic fields reveal themselves via their influence on the surrounding matter, and as such are notoriously hard to measure independently of other tracers. Aims. In this work, we attempt to disentangle an all-sky map of the line-of-sight (LoS)-parallel component of the Galactic magnetic field from the Faraday effect, utilizing several tracers of the Galactic electron density, ne. Additionally, we aim to produce a Galactic electron dispersion measure map and quantify several tracers of the structure of the ionized medium of the Milky Way. Methods. The method developed to reach these aims is based on information field theory, a Bayesian inference framework for fields, which performs well when handling noisy and incomplete data and constraining high-dimensional-parameter spaces. We rely on compiled catalogs of extragalactic Faraday rotation measures and Galactic pulsar dispersion measures, a well as data on bremsstrahlung and the hydrogen α spectral line to trace the ionized medium of the Milky Way. Results. We present the first full sky map of the LoS-averaged Galactic magnetic field. Within this map, we find LoS-parallel and LoS-averaged magnetic field strengths of up to 4 µG, with an all-sky root mean square of 1.1 µG, which is consistent with previous local measurements and global magnetic field models. Additionally, we produce a detailed electron dispersion measure map that agrees with existing parametric models at high latitudes but suffers from systematic effects in the disk. Further analysis of our results with regard to the 3D structure of ne reveals that it follows a Kolmogorov-type turbulence for most of the sky. From the reconstructed dispersion measure and emission measure maps, we construct several tracers of variability in ne along the LoS. Conclusions. This work demonstrates the power of consistent joint statistical analysis including multiple datasets and physical quantities and defines a road map toward a full three-dimensional joint reconstruction of the Galactic magnetic field and the ionized ISM.

Observations of Extremely Metal-poor O Stars: Weak Winds and Constraints for Evolution Models

Metal-poor massive stars drive the evolution of low-mass galaxies, both locally and at high redshift. However, quantifying the feedback they impart to their local surroundings remains uncertain because models of stellar evolution, mass loss, and ionizing spectra are unconstrained by observations below 20% solar metallicity (Z ⊙). We present new Keck Cosmic Web Imager optical spectroscopy of three O-type stars in the nearby dwarf galaxies Leo P, Sextans A, and WLM, which have gas-phase oxygen abundances of 3%–14% Z ⊙. To characterize their fundamental stellar properties and radiation-driven winds, we fit PoWR atmosphere models to the optical spectra simultaneously with Hubble Space Telescope far-ultraviolet (FUV) spectra and multiwavelength photometry. We find that all three stars have effective temperatures consistent with their spectral types and surface gravities typical of main-sequence dwarf stars. Yet, the combination of those inferred parameters and luminosity for the two lower-Z stars is not reproduced by stellar evolution models, even those that include rotation or binary interactions. The scenario of multiple-star systems is difficult to reconcile with all available data, suggesting that these observations pose a challenge to current evolution models. We highlight the importance of validating the relationship between stellar mass, temperature, and luminosity at very low Z for accurate estimates of ionizing photon production and spectral hardness. Finally, all three stars’ FUV wind profiles reveal low mass-loss rates and terminal wind velocities in tension with expectations from widely adopted radiation-driven wind models. These results provide empirical benchmarks for future development of mass-loss and evolution models for metal-poor stellar populations.

HR-GO

Context. Dwarf galaxy streams encode vast amounts of information essential to understanding early galaxy formation and nucleosynthesis channels. Due to the variation in the timescales of star formation history in their progenitors, stellar streams serve as ‘snapshots’ that record different stages of galactic chemical evolution. Aims. This study focusses on the Cetus stream, stripped from a low-mass dwarf galaxy. We aim to uncover its chemical evolution history as well as the different channels of its element production from detailed elemental abundances. Methods. We carried out a comprehensive analysis of the chemical composition of 22 member stars based on their high-resolution spectra. We derived abundances for up to 28 chemical species from C to Dy and, for 20 of them, we account for the departures from local thermodynamic equilibrium (NLTE effects). Results. We confirm that the Cetus stream has a mean metallicity of [Fe/H] = −2.11 ± 0.21. All observed Cetus stars are α enhanced with [α/Fe] ≃ 0.3. The absence of the α-‘knee’ implies that star formation stopped before iron production in type Ia supernovae (SNe Ia) became substantial. Neutron capture element abundances suggest that both the rapid (r-) and the main slow (s-) processes contributed to their origin. The decrease in [Eu/Ba] from a typical r-process value of [Eu/Ba] = 0.7–0.3 with increasing [Ba/H] indicates a distinct contribution of the r- and s-processes to the chemical composition of different Cetus stars. For barium, the r-process contribution varies from 100 to 20% in different sample stars, with an average value of 50%. Conclusions. Our abundance analysis indicates that the star formation in the Cetus progenitor ceased after the onset of the main s-process in low- to intermediate-mass asymptotic giant branch stars but before SNe Ia played an important role. A distinct evolution scenario is revealed by comparing the abundances in the Ursa Minor dwarf spheroidal galaxy, showing the diversity in – and uniqueness of – the chemical evolution of low-mass dwarf galaxies.

From 100 MHz to 10 GHz: Unveiling the spectral evolution of the X-shaped radio galaxy in Abell 3670

Context. X-shaped radio galaxies (XRGs) are characterised by two pairs of misaligned lobes: active lobes hosting radio jets and the wings. None of the formation mechanisms proposed thus far are able to exhaustively reproduce the diverse features observed among XRGs. Emerging evidence has proposed the existence of sub-populations of XRGs forming via different processes. Aims. The brightest cluster galaxy (BCG) in Abell 3670 (A3670) is a dumbbell system hosting the XRG MRC 2011-298. The morphological and spectral properties of this interesting XRG were first characterised based on Karl G. Jansky Very Large Array (JVLA) data at 1–10 GHz. In the present work, we follow up on MRC 2011-298 with the upgraded Giant Metrewave Radio Telescope (uGMRT) at 120–800 MHz to further constrain its properties and origin. Methods. We carried out a detailed spectral analysis sampling different spatial scales. Integrated radio spectra, spectral index maps, radio colour-colour diagrams, and radiative age maps of both the active lobes and prominent wings were employed to test the origin of the source. Results. We confirm a progressive spectral steepening from the lobes to the wings. The maximum radiative age of the source is ~80 Myr, with the wings being older than the lobes by ≳30 Myr in their outermost regions. Conclusions. The observed properties are in line with an abrupt reorientation of the jets by ~90 deg from the direction of the wings to their present position. This formation mechanism is further supported by the comparison with numerical simulations in the literature, which additionally highlight the role of hydrodynamic processes in the evolution of large wings such as those of MRC 2011-298. It is plausible that the coalescence of supermassive black holes could have triggered the spin-flip of the jets. Moreover, we show that the S-shape of the radio jets is likely driven by precession with a period of P ~ 10 Myr.

The Velocity Dispersion Function for Quiescent Galaxies in Massive Clusters from IllustrisTNG

We derive the central stellar velocity dispersion function (VDF) for quiescent galaxies in 280 massive clusters with log(M200/M⊙)>14 in IllustrisTNG300. The VDF is an independent tracer of the dark matter mass distribution of subhalos in galaxy clusters. Based on the IllustrisTNG cluster catalog, we select quiescent member subhalos with a specific star formation rate <2 × 10−11 yr−1 and stellar mass log(M*/M⊙)>9 . We then simulate fiber spectroscopy to measure the stellar velocity dispersion of the simulated galaxies; we compute the line-of-sight velocity dispersions of star particles within a cylindrical volume that penetrates the core of each subhalo. We construct the VDFs for quiescent subhalos within R 200. The simulated cluster VDF exceeds the simulated field VDF for logσ*>2.2 , indicating the preferential formation of large velocity dispersion galaxies in dense environments. The excess is similar in simulations and in the observations. We also compare the simulated VDF for the three most massive clusters with log(M200/M⊙)>15 with the observed VDF for the two most massive clusters in the local Universe, Coma and A2029. Intriguingly, the simulated VDFs are significantly lower for logσ*>2.0 . This discrepancy results from (1) a smaller number of subhalos with log(M*/M⊙)>10 in TNG300 compared to the observed clusters, and (2) a significant offset between the observed and simulated M *–σ * relations. The consistency in the overall shape of the observed and simulated VDFs offers a unique window into galaxy and structure formation in simulations.

Core Formation by Binary Scouring and Gravitational Wave Recoil in Massive Elliptical Galaxies

Scouring by supermassive black hole (SMBH) binaries is the most accepted mechanism for the formation of the cores seen in giant elliptical galaxies. However, an additional mechanism is required to explain the largest observed cores. Gravitational-wave (GW) recoil is expected to trigger further growth of the core, as subsequent heating from dynamical friction of the merged SMBH removes stars from the central regions. We model core formation in massive elliptical galaxies from both binary scouring and heating by GW recoil and examine their unique signatures. We aim to determine whether the nature of cores in 3D space density can be attributed uniquely to either process and whether the magnitude of the kick can be inferred. We perform N-body simulations of galactic mergers of multicomponent galaxies, based on the observed parameters of four massive elliptical galaxies with cores >0.5 kpc. After binary scouring and hardening, the merged SMBH remnant is given a range of GW recoil kicks with 0.5–0.9 of the escape speed of the galaxy. We find that binary scouring alone can form the cores of NGC 1600 and A2147-BCG, which are <1.3 kpc in size. However, the >2 kpc cores in NGC 6166 and A2261-BCG require heating from GW recoil kicks of <0.5 of the galaxy escape speed. A unique feature of GW recoil heating is flatter cores in surface brightness, corresponding to truly flat cores in 3D space density. It also preferentially removes stars on low angular momentum orbits from the galactic nucleus.