2025-04-03 | | Total: 36
Star clusters are crucial for understanding how stars evolve. Their colour-magnitude diagrams show the effects of stellar evolution of approximately coeval objects with the same chemical composition. Furthermore, the determination of their astrophysical parameters (age, distance, colour excess and metallicity) together with their spatial distribution provides information about the structure and the evolution of the Galaxy itself. Using data from the \textit{Gaia} DR3 and 2MASS catalogues, we develop methodologies for characterizing open clusters. Precise membership lists, mean astrometric parameters and radii are obtained. Using photometric data from both data sources, we carried out new age calibrations that rely on morphological indices based on colour (ΔBR) and magnitude (\Delta G) differences between the red clump and the turnoff for a sample of 34 open clusters with ages covering the interval 8.3 < \log[t({\rm yr})] < 9.9. A set of age calibration functions based on \textit{Gaia} morphological age indices are determined for the first time. We demonstrate their accuracy, obtaining a mean residual of 0.06 dex in \log[t(yr)]. Our results also show that stellar evolution models tend to predict the difference \Delta G. However, they typically overestimate the difference \Delta BR for objects younger than \log[t({\rm yr})] = 8.8.
Context: Using Gaia DR3 data, Jiménez-Arranz et al. (2024a) computed the LMC bar pattern speed using three different methods. One of them suggested that the LMC might be hosting a bar that barely rotates, and is slightly counter-rotating with respect to the LMC disc, with a pattern speed of \Omega_p = -1.0 \pm 0.5 km/s/kpc. Aims: To confirm that tidal interactions can trigger the LMC hosting a non-rotating bar due to its interaction with the SMC, which could cause the LMC bar to slow down significantly until it (momentarily) stops. Methods: We analyse a subset of models (K9 and K21) from the KRATOS suite (Jiménez-Arranz et al. 2024b) where we detected non-rotating bars. We make use of two different methods to track the evolution of the bar pattern speed: the program patternSpeed.py (Dehnen et al. 2023), and temporal finite-differences of the change rate the bar major axis' phase angle. Results: In the second LMC-SMC-like pericenter passage of K9, the bar of the LMC-like galaxy weakens to almost disappear and regenerates with a pattern speed that suffers a slowdown from \Omega_p \sim 20 km/s/kpc to \Omega_p \sim 0 km/s/kpc in less than 75 Myr. Then, the bar rotates at less than \Omega_p \sim 3-5 km/s/kpc for around 100 Myr, until it recovers the initial (before interaction) pattern speed of \Omega_p \sim 10 km/s/kpc. The results for the K21 simulation are comparable. Conclusions: This work is the first direct evidence that galactic bars can be slowed down or even stopped by tidal interaction, which strengthens the possibility of the LMC hosting a non-rotating bar, and can add an alternative formation scenario for observed slow-rotating bars.
Broadband photometric reverberation mapping (RM) provides a measure of the size of the continuum-emitting region in active galactic nuclei (AGN). Previous monitoring campaigns of PG 2130+099 disagree as to whether the continuum emitting region size is consistent with that predicted for a standard optically thick geometrically thin accretion disk. We present \sim6 months of observations from several robotic telescopes, providing the highest cadence and widest wavelength coverage photometric RM study of PG 2130+099 to date. Our results indicate that inferred size of the continuum-emitting region in PG 2130+099, like many recently observed AGN, is larger than the simplest predictions for an irradiated geometrically thin, optically thick accretion disk. We also perform a flux-flux analysis, finding a variable spectrum broadly consistent with a disk, and a constant component with enhanced \textit{i}-band emission, potentially due to H\alpha. We find some evidence of increasing lag with luminosity, but previous lag measurements are too uncertain to be definitive.
We present a comprehensive analysis of the MIRI Extremely Red Object Virgil, a Lyman-\alpha emitter at z_{spec} = 6.6379 \pm 0.0035 with the photometric properties of a Little Red Dot. Leveraging new JWST/MIRI imaging from the MIDIS and PAHSPECS programs, we confirm Virgil's extraordinary nature among galaxies in JADES/GOODS-South, exhibiting a strikingly red NIRCam-to-MIRI color (F444W - F1500W = 2.84\pm0.04~mag). Deep NIRSpec/PRISM spectroscopy from the OASIS program offers key insights into the host galaxy, revealing properties of an average star-forming galaxy during Cosmic Reionization, such as a subsolar metallicity, low-to-moderate dust content, and a relatively high ionization parameter and electron temperature. By estimating the star formation rate of Virgil from UV and H\alpha, we find evidence that the galaxy is either entering or fading out of a bursty episode. Although line-ratio diagnostics employed at high-z would classify Virgil as an Active Galactic Nucleus (AGN), this classification becomes ambiguous once redshift evolution is considered. Nonetheless, Virgil occupies the same parameter space as recently confirmed AGNs at similar redshifts. The new deep MIRI data at 15~\mum reinforce the AGN nature of Virgil, as inferred from multiple spectral energy distribution (SED) fitting codes. Virgil's rising infrared SED and UV excess resemble those of Dust-Obscured Galaxies (DOGs) studied with Spitzer at Cosmic Noon, particularly blue-excess HotDOGs. Our results highlight the need for a multi-wavelength approach incorporating MIRI to uncover such extreme sources at z\gtrsim6 and to shed light on the interplay between galaxy evolution and early black hole growth during Cosmic Reionization.
The nature and formation history of our Galaxy's largest and most enigmatic stellar cluster, known as Omega Centauri (ocen) remains debated. Here, we offer a novel approach to disentangling the complex stellar populations within ocen based on phylogenetics methodologies from evolutionary biology. These include the Gaussian Mixture Model and Neighbor-Joining clustering algorithms applied to a set of chemical abundances of ocen stellar members. Instead of using the classical approach in astronomy of grouping them into separate populations, we focused on how the stars are related to each other. In this way, we could identify stars that likely formed in globular clusters versus those originating from prolonged in-situ star formation and how these stars interconnect. Our analysis supports the hypothesis that ocen might be a nuclear star cluster of a galaxy accreted by the Milky Way with a mass of about 10^9M_sun. Furthermore, we revealed the existence of a previously unidentified in-situ stellar population with a distinct chemical pattern unlike any known population found in the Milky Way to date. Our analysis of ocen is an example of the success of cross-disciplinary research and shows the vast potential of applying evolutionary biology tools to astronomical datasets, opening new avenues for understanding the chemical evolution of complex stellar systems.
The vertical metallicity gradient of the Galactic disk offers valuable insights into the disk's formation and chemical evolution over time. We utilized the LAMOST-LRS young stellar sample to investigate this gradient and found that it approaches zero as stellar effective temperature (or age) increases (or decreases) across various Galactocentric distances. To validate this result, we analyzed 295 open clusters younger than 3 Gyr and 976 classical cepheids within the Galactic disk. The findings confirmed that, within a given narrow age range, the vertical metallicity gradient is effectively zero. This relationship between metallicity and age supports the ``upside-down'' disk formation theory, as it indicates that the youngest and most metal-rich stars dominate the midplane, while older and more metal-poor stars formed at larger vertical heights and currently tend to be at these heights. Overall, our results align well with theoretical predictions, offering further insight into the chemical evolution and structural properties of the Milky Way.
We use HDGAS hydrodynamic simulations to study the impact of active galactic nucleus (AGN) feedback on the conversion of atomic-gas to molecular-gas within the circumnuclear-disc (CND) of a typical AGN-dominated galaxy. The comparison of CI, CII, and CO line intensities and their ratios in the HDGAS post-processing radiative-transfer analysis reveals the complex interplay between AGN-activity, cold molecular gas properties, and the physical processes governing the evolution of star-formation in galaxies. Our results demonstrate that the CI/CO intensity ratio serves as a reliable indicator of the atomic-to-molecular gas transition. We present the probability distribution function (PDF) and abundance trends of various metal species related to molecular H2 gas, highlighting differences in clumpiness and intensity maps between AGN feedback and NoAGN models. The profile of the integrated intensity (moment-0) maps shows that the AGN-feedback model exhibits a lower CI/CO intensity ratio in the vicinity of the supermassive black hole (< 50 pc), indicating a smaller atomic-gas abundance and the presence of positive AGN-feedback. Our simulations have successfully predicted the presence of faint-CO emissions extending to larger radii from the galactic center. We also explore the relationships between CII/CO and CI/CII intensity ratios, as well as the ratios versus CO intensity, which provides insights into the "CO-dark" issues. One notable feature in the later time-scale of the AGN model is the presence of a "CO-dark" region, where the intensity of CO emission (\rm I_{CO}) is depleted relative to the H_2 column density (N_{\rm H_2}) compared to the NoAGN model.
Cas OB5 is an OB association located at a distance of 2.5-3 kpc that intercepts the Perseus spiral arm. It carries a moderate amount of reddening (A_V \sim 2-3 mag) and contains several well-known open clusters within its boundaries, such as King 12, NGC 7788, and NGC 7790. The availability of modern clustering algorithms, together with \textit{Gaia} DR3 kinematics and complementary spectroscopic data, makes it a suitable site for studies of Galactic structure. We seek to quantify the spatial scale of star formation in the spiral arms, using Cas OB5 as a pilot target before extending our study to more distant and extinguished regions of the Galaxy. We selected 129,695 candidate OBA stars in a 6x8 deg^2 region around Cas OB5. We applied a spectral energy distribution (SED) fitting process to this sample to derive the physical parameters. Through this process, we found 56 379 OBA stars, which we then clustered using HDBSCAN. We identified 17 open clusters inside this area, four of which appear to form a coherent structure that we identify as Cas OB5. Nevertheless, our findings suggest that these clusters belong to two different age groups despite sharing a similar position and kinematics. Spectroscopic observations confirm the youth of NGC 7788 (10-15 Myr) compared to NGC 7790 (110\pm15\:Myr). We have determined a spatial scale for star formation of a few tens of pc to a few hundreds of pc, comparing the clustered to the diffuse population of Cas~OB5 across this part of the Perseus arm. A spectroscopic analysis was required to complement the clustering algorithm, so that we could separate younger OCs (tracers of the spiral arm) from older ones. These results highlight the need to combine these techniques to fully disentangle the Milky Way structure.
Type-2 quasars (QSO2s) are active galactic nuclei (AGN) seen through a significant amount of dust and gas that obscures the central supermassive black hole and the broad line region. Here we present new mid-infrared spectra of the central kiloparsec of five optically-selected QSO2s at redshift z~0.1 obtained with JWST/MIRI/MRS. These QSO2s belong to the QSOFEED sample and they have log Lbol=45.5-46.0 erg/s, global SFRs that place them above the main sequence, and practically identical optical spectral shape and [OIII] luminosity, but their nuclear mid-infrared spectra exhibit an unexpected diversity of both continua and features. They show: 1) 9.7 micron silicate features going from emission (strength of S9.7=0.5) to relatively strong absorption (S9.7=-1.0) and 18 and 23 micron silicates either in emission or flat. In addition, two of the QSO2s show absorption bands of CO, H2O, and aliphatic grains, indicating different levels of nuclear obscuration across the sample. 2) [NeV]/[NeII] ratios ranging from 0.1 to 2.1 and [NeIII]/[NeII] from 1.0 to 3.5, indicating different coronal line and ionizing continuum strengths. 3) Warm molecular gas masses of 1-4x10^7 Msun and warm-to-cold gas mass ratios of 1-2%, with molecular gas excitation likely due to jet-induced shocks in J1430+1339, and to UV heating and/or turbulence in J1509+0434. 4) PAH emission features with equivalent widths ranging from <0.002 to 0.075 micron, from which we measure a larger contribution from neutral molecules (PAH 11.3/6.2=1.3-3.4) and SFRs<3-7 Msun/yr. This unprecedented dataset allowed us to start exploring the role of various AGN and galaxy properties including ionizing continuum, obscuration, electron density, and jet-ISM interactions on some of the spectral differences listed above, but larger samples are now required to fully understand the diversity of QSO2s' nuclear mid-infrared spectra.
We applied Stars' Galactic Origin (StarGO) algorithm in 7D space (i.e., [Fe/H], [Mg/Fe], [Al/Fe], L_{z}, J_{r}, J_{z}, E) to analyze stars in the inner halo with APOGEE DR17. We identified some known substructures, including Gaia-Sausage-Enceladus (GSE), Sagittarius Stream (Sgr), LMS-1 (Wukong), Thamnos, Metal-Weak Thick Disk (MWTD) and Aleph. Additionally, we identified an undefined metal-poor group (UDG, with [Fe/H] <-0.8 dex) probably linked to known substructure like Aleph, as well as a high \alpha-abundance substructure (HAS) associated with both the Nyx and Nyx-2 streams. Chemical abundance of the HAS supports the argument that Nyx and Nyx-2 share a common origin. We discovered three substructures, which we refer to as new-substructure candidate-1,2,3 (NSTC-1,2,3). Despite exhibiting disk-like dynamics, these NSTCs demonstrate notably low [Mg/Fe] (< 0.2 dex) and [Al/Fe] (< -0.15 dex), similar to the properties of dwarf galaxies. Their high orbital energy and low [\alpha/Fe] indicate the association with recent accretion events.
Using elemental abundances for 1.26 million K giants in the LAMOST DR8 value-added catalog, we analyze the chemical abundances of the Virgo Overdensity (VOD) and Hercules-Aquila Cloud (HAC). We find two distinct chemical populations in both overdensities, which is in disagreement with the mainstream hypothesis that both overdensities are composed of materials from a single merger event, namely Gaia-Sausage-Enceladus (GSE). The two populations show different chemical trends: one exhibits low metallicities and high \alpha abundances, and the other shows high metallicities and low \alpha abundances, which is associated with the recently discovered Nereus and Virgo Radial Merger (VRM) components in the local stellar halo, respectively. The Nereus component in these overdensities uniquely exhibits a decreasing trend in the [Fe/H]-[Mn/Fe] plane. Out of all observed Milky Way dwarf galaxies, this trend is only found in the Sculptor dwarf galaxy, which provides clues for the properties of Nereus progenitor. We also find that the velocity ellipse with high aniostropy parameters that is usually considered to be part of GSE are actually a mix of the two components. Both overdensities are well-mixed in kinematic spaces, confirming recent claims that the debris of merger pairs are kinematically indistinguishable in a recent simulation. We find that the velocity ellipses of the VRM stars in these overdensities have large inclination angles, which may be an indication of the merger time in simulated merger events.
In this work an analysis of the intracluster light (ICL) in the galaxy cluster SMACS J0723.3-7327 (hereafter, SMACS J0723) using JWST/NIRCam deep imaging in six filters (F090W to F444W) is presented. The images were processed for low surface brightness (LSB) science, with additional correction for instrumental scattering in the short-wavelength channels, and analysed using wavelet-based decomposition. The ICL, brightest cluster galaxy (BCG), and satellite galaxies were extracted and modelled, with 2D maps for each component. ICL and ICL+BCG fractions, computed across all filters within a 400 kpc radius, exhibit a flat trend with wavelength, averaging 28% and 34%, respectively. Flux ratios between the BCG and the next brightest members (M_{12}, M_{13} and M_{14}) also display minimal wavelength dependence. These results indicate that SMACS J0723 is a dynamically evolved cluster with a dominant BCG and well-developed ICL. Five prominent ICL substructures are analysed, contributing to 10-12% of the total ICL+BCG flux budget, slightly exceeding simulation predictions. Their short dynamical timescales suggest an instantaneous ICL injection rate of several 10^3 L_{\odot}\,{\rm yr}^{-1}, consistent with active dynamical assembly. These findings support a scenario where SMACS J0723's ICL growth is currently driven by galaxy mergers involving the BCG and other bright satellites, rather than by the accretion of pre-processed ICL from a recent cluster merger. However, extrapolating the current injection rate to the cluster's lifetime indicates that additional mechanisms are required to match the growth observed in other clusters over cosmic times.
The interstellar medium (ISM) is the material that fills the space between the stars in all galaxies; it is a multi-phase medium in pressure equilibrium, with densities and temperatures covering over 6 orders of magnitude. Although accounting for only a small fraction of the mass of any given galaxy, it is a vital component, since it holds the material responsible for galaxy growth through star formation. Studying the ISM requires careful observations at all wavelengths of the electromagnetic spectrum. This article describes the multi-phase nature of the ISM, and then puts it in the context of galaxy evolution models, emphasising the importance of the cycling of baryons in and out of galaxies. Within this framework, the ISM plays a central role: it connects the physical processes operating on very large physical- and time-scales which control the accretion of gas onto galaxies, and the small scale processes that regulate star formation.
We investigate the intrinsic and observational properties of z\gtrsim 6 galaxies hosting coalescing massive black holes (MBHs) that gives rise to gravitational waves (GWs) detectable with the Laser Interferometer Space Antenna (LISA). We adopt a zoom-in cosmological hydrodynamical simulation of galaxy formation and black hole (BH) co-evolution, zoomed-in on a M_h \sim 10^{12}~\rm M_{\odot} dark matter halo at z = 6, which hosts a fast accreting super-massive black hole (SMBH) and a star-forming galaxy. Following the SMBH formation backward in time, we identify the merging events that concurred to its formation and we pick up the ones that are detectable with LISA. Among these LISA detectable events (LDEs), we select those that, based on their intrinsic properties are expected to be bright in one or more electromagnetic (EM) bands. We post-process these events with dust radiative transfer calculations to make predictions about their spectral energy distributions and continuum maps in the JWST to ALMA wavelength range. We compare the spectra arising from galaxies hosting the merging MBHs with those arising from AGN powered by single accreting BHs. We find that it will be impossible to identify an LDE from the continuum SEDs because of the absence of specific imprints from the merging MBHs. We also compute the profile of the H_{\rm \alpha} line arising from LDEs, considering the contribution from their star-forming regions and the accreting MBHs. We find that the presence of two accreting MBHs would be difficult to infer even if both MBHs accrete at super-Eddington rates. We conclude that the combined detection of GW and EM signals from z\gtrsim 6 MBHs is challenging not only because of the poor sky-localization provided by LISA, but also because the loudest GW emitters are not massive enough to leave significant signatures in the emission lines arising from the broad line region.
The James Webb Space Telescope (JWST) is challenging our understanding of the nature of the very first galaxies in the Universe, having discovered a surprising abundance of very massive galaxies in early cosmic epochs. By applying a model of primeval dust, we estimate the far-infrared (FIR)-continuum luminosities for galaxies of different masses at redshifts z\gtrsim7. In particular, we predict observed fluxes for different available bands (3-10) of the Atacama Large Millimeter/sub-mm Array (ALMA), considering typical conservative values expected for the properties of first galaxies (e.g., gas-phase metallicities, dust-to-metal ratio, star formation efficiency). As expected, FIR fluxes increase with stellar mass for all ALMA bands, but with a steeper slope for bands 9 and 10. Encouragingly, first galaxies are affected by a strong negative-K correction, in such a way that sources with similar properties are brighter in bands 3-8 at higher redshifts. Such behavior is stronger for bands 6-7. Although the trends for bands 9-10 are not clear, the highest fluxes for such bands are reached towards extreme z\gtrsim 15. Counterintuitively, our results suggest that JWST sources with similar masses and dust properties would be more easily detectable with ALMA if they are located at higher z.
Studying the atomic-to-molecular transition is essential for understanding the evolution of interstellar medium. The linear edge of Taurus molecular cloud, clearly identified in the ^{13}CO(1-0) intensity map, serves as an ideal site for investigating this transition. Utilizing the Arizona Radio Observatory Sub-Millimeter Telescope, we obtained mapping observations of CO(2-1), ^{13}CO(2-1), and CO(3-2) across this linear edge. The intensity ratio between CO(2-1) and ^{13}CO(2-1) indicates a lower limit of the ^{12}C/^{13}C ratio of 54\pm 17. Based on multi-transition observations of CO and ^{13}CO, we performed Markov Chain Monte Carlo (MCMC) fit of the physical properties across this edge using non-Local Thermodynamic Equilibrium analysis with the RADEX code, based on the Large velocity Gradient (LVG) assumption. The number density profile exhibits a pronounced jump coinciding with the H_2 infrared emission peak. The cold HI gas within the molecular cloud, manifested as HI-Narrow Self-Absorption (HINSA) features, is detected along the cloud edge. Our quantitative comparison with numerical simulations provides tentative evidence that shocks induced by colliding gas flows may contribute to the atomic-to-molecular phase transition observed along the linear edge.
One of the crucial components in simulating the growth and evolution of galaxies within a cosmological framework is the modeling of star formation (SF) and its corresponding feedback. Traditionally, the implemented SF law follows the empirical Kennicutt-Schmidt relation, which links the SF rate (SFR) in a gas element to the total gas density. More recently, an even stronger correlation has been observed between the SFR and the amount of molecular hydrogen (\mathrm{H}_2). This opens up the question of whether molecular hydrogen is a necessary precursor for SF or is instead a tracer of the total amount of gas, both of which would explain the observed correlations. In this study, we examine the impact of using an \mathrm{H}_2-based SF law on the formation of Milky Way-mass galaxies using cosmological, hydrodynamical simulations. We find that the galaxy modeled with our approach exhibits a well-defined, disc-like morphology. Compared to the traditional recipe, our model delays the onset of SF by approximately 500 \, \mathrm{Myr} resulting in a lower SFR, a smaller disc size, and a higher proportion of neutral to ionized gas within the disc region. These findings highlight the importance of including sophisticated SF models which can be compared to several observations -- including those related to \mathrm{H}_2 -- to better understand the processes affecting galaxy formation and evolution.
We investigate the influence of the cosmic web on galaxy properties in the IllustrisTNG simulations. To disentangle the effects of galaxy groups and cosmic filaments, we classify the cosmic web environment into four categories: group, group-dominated, filament-dominated, and field. By controlling for stellar mass, we reveal evident differences in specific star formation rates (sSFR), quenched fraction, gas fractions, local density, and stellar ages among central galaxies in different cosmic web environments, particularly for lower-mass galaxies. However, these differences largely diminish when the effect of local overdensity is further accounted for, indicating its dominant role. Additionally, we observe distinct differences in these properties among satellite galaxies across environments, mainly driven by stellar mass, halo mass, and overdensity. Notably, residual differences between satellites in field and filament-dominated region persist even after controlling for these factors, suggesting a stronger susceptibility of satellite galaxies to filaments compared to centrals. Our findings highlight the importance of differentiating between central and satellite to accurately assess the environmental effects of the cosmic web. Our analysis suggests that the relationship between galaxy properties and their distance from filaments arises from a combination of factors, including stellar and halo mass, groups, overdensity, and the intrinsic influence of the cosmic web. Additionally, we find that the effect of the cosmic web on galaxy properties is reduced at z=0.5, compared to z=0. Furthermore, central galaxies near thick filaments tend to exhibit slightly to moderately lower sSFR and cold gas fractions compared to those near thin filaments.
Mass loss is a key aspect of stellar evolution, particularly in evolved massive stars, yet episodic mass loss remains poorly understood. To investigate this, we need evolved massive stellar populations across various galactic environments. However, spectral classifications are challenging to obtain in large numbers, especially for distant galaxies. We addressed this by leveraging machine-learning techniques. We combined \textit{Spitzer} photometry and Pan-STARRS1 optical data to classify point sources in 26 galaxies within 5 Mpc, and a metallicity range 0.07-1.36 Z_\odot. \textit{Gaia} DR3 astrometry was used to remove foreground sources. Classifications are derived using a machine-learning model developed by Maravelias et al. (2022). We report classifications for 1,147,650 sources, with 276,657 sources (\sim24\%) being robust. Among these are 120,479 Red Supergiants (RSGs; \sim11\%). The classifier performs well even at low metallicities (\sim0.1 Z_\odot) and distances under 1.5 Mpc, with a slight decrease in accuracy beyond \sim3 Mpc due to \textit{Spitzer}'s resolution limits. We also identified 21 luminous RSGs (\textrm{log}(L/L_\odot)\ge5.5), 159 dusty Yellow Hypergiants in M31 and M33, as well as 6 extreme RSGs (\textrm{log}(L/L_\odot)\ge6) in M31, challenging observed luminosity limits. Class trends with metallicity align with expectations, though biases exist. This catalog serves as a valuable resource for individual-object studies and \textit{James Webb} Space Telescope target selection. It enables follow-up on luminous RSGs and Yellow Hypergiants to refine our understanding of their evolutionary pathways. Additionally, we provide the largest spectroscopically confirmed catalog of massive stars and candidates to date, comprising 5,273 sources (including \sim330 other objects).
We investigate the capabilities of the mid-infrared instrument (MIRI) of James Webb Space Telescope (JWST) to advance our knowledge of AGN dust using the spectral fitting technique on an AGN collection of 21 nearby (z<0.05) AGN (7 type-1 and 14 type-2) observations obtained with the medium resolution spectroscopy (MRS) mode. This collection includes publicly available AGN and data from the collaboration of Galactic Activity, Torus, and Outflow Survey (GATOS). We developed a tool named MRSPSFisol that decomposes MRS cubes into point-like and extended contributions. We found statistically good fits for 12 targets with current AGN dust models. The model that provides good fits (chi2/dof<2) for {these 12 targets} assumes a combination of clumpy and smooth distribution of dust in a flare-disk geometry where the dust grain size is a free parameter. Still, two and one AGN statistically prefer the disk+wind and the classical clumpy torus model, respectively. However, the currently available models fail to reproduce 40% of the targets, likely due to the extreme silicate features not well reproduced by the models and signatures of water-ice and aliphatic hydrocarbon absorption features in most targets. New models exploring, for instance, new chemistry, are needed to explain the complexity of AGN dust continuum emission observed by JWST.
High resolution surveys reveal that the interstellar medium in the Milky Way and nearby galaxies consists of interlinked hierarchies of filamentary structure and superbubbles extending from galactic to subpc scales. The characterization of filament properties across this hierarchy is of fundamental importance for the origin of giant molecular clouds and their star clusters. In this paper we characterize the properties of filaments greater than 25 pc in length that are produced in the multi-scale galactic MHD simulations of Zhao et al. 2024. By adapting the FilFinder algorithm of Koch & Rosolowsky, 2015, we extract over 500 filaments ranging up to 10 kpc scales, to derive the probability distribution functions for filament masses and lengths, magnetic field orientations, and the gravitational stability and fragmentation patterns of filaments. We find power-law distributions for filament masses and lengths. The former has a power law index \alpha_m = 1.85 that is nearly identical to that of observed GMC mass functions in extragalactic and Galactic surveys, suggesting that GMC properties are inherited from their host filaments. The fragmentation of magnetized filaments on 200 pc scales or less occurs when they exceed an average critical line mass, as predicted by theory. On larger scales however, kpc filaments form out of the cold neutral medium (CNM) and fragmentation follows local variations in the critical line mass along spiral arms or at the boundaries of superbubbles.
A systematic search for HI 21-cm absorption in Quasar-Galaxy Pairs (QGPs) provides a powerful means to map the distribution of cold gas around high-redshift star-forming galaxies. Fiber spectroscopy of high-redshift quasars enables the serendipitous detection of foreground star-forming galaxies at extremely small impact parameters, forming a unique subset of QGPs known as Galaxies On Top Of Quasars (GOTOQs). In this study, we present results from a pilot upgraded Giant Metrewave Radio Telescope (uGMRT) survey of three GOTOQs, where we achieved a remarkable 100\% detection rate of HI 21-cm absorption. By combining our findings with existing literature, we establish that GOTOQs constitute a distinct population in terms of HI 21-cm absorption, with significantly higher detection rates than those observed in Damped Lyman-\alpha (DLA)-based or metal absorption-based searches. For the GOTOQs, we find a strong correlation between the line-of-sight reddening and the HI 21-cm optical depth, characterized by \int \tau\, dv\, (\rm{km\,s^{-1}}) = 13.58^{+2.75}_{-2.35} E(B-V) + 0.68^{+1.06}_{-1.27}, consistent with the Milky Way sightlines. We also show that the HI 21-cm detection rates and optical depth declines rapidly with the impact parameter. With upcoming wide-field spectroscopic surveys expected to substantially expand the catalog of known GOTOQs, the success of this pilot survey lays the foundation for constructing a statistically significant sample of intervening HI 21-cm absorbers.
Galaxies are not isolated systems; they continuously interact with their surroundings by ejecting gas via stellar feedback and accreting gas from the environment. Understanding the interplay between outflows from the disc and the surrounding circumgalactic medium (CGM) is key to learning how star-forming galaxies evolve. Our goal is to understand how gas in the CGM is accreted onto the inner regions of the disc, making it available for the formation of stars, exploring the connection between stellar feedback and gas accretion from the CGM in Milky Way-like galaxies. We focus on the distribution of vertical and radial gas flows to and from the disc as a function of galactocentric radius, and examine the implications of these processes for the evolution of such galaxies. We use the Arepo code coupled with the SMUGGLE sub-grid model to perform hydrodynamic N-body simulations of 9 different galaxies surrounded by a hot CGM. Each simulation features a gaseous disc with different mass and scale length, allowing us to examine how disc structure impacts gas dynamics. We find evidence of a crucial link between stellar feedback and gas accretion from the CGM, which together play an essential role in sustaining ongoing star formation in the disc. In particular, the ejection of gas from the disc plane by stellar feedback leads to the generation of a baryon cycle in which the CGM gas is mainly accreted onto the external regions of the disc ( \approx 3-10 M_{\odot} yr^{-1} of gas is accreted into the whole disc). From these regions it is then transported to the centre with radial mass rates \approx 1-4 M_{\odot} yr^{-1} on average, owing to angular momentum conservation, forming new stars and starting the whole cycle again. We find that both vertical accretion onto the inner regions of the disc and the radial transport of gas from the disc outskirts are necessary to sustain star formation.
Pulsar Timing Array (PTA) experiments have the potential to unveil continuous gravitational wave (CGW) signals from individual massive black hole binaries (MBHBs). Detecting them in both gravitational waves (GW) and the electromagnetic (EM) spectrum will open a new chapter in multimessenger astronomy. We investigate the feasibility of conducting multimessenger studies by combining the CGW detections from an idealized 30-year SKA PTA and the optical data from the forthcoming LSST survey. To this end, we employed the \texttt{L-Galaxies} semi-analytical model applied to the \texttt{Millennium} simulation. We generated 200 different all-sky lightcones that include galaxies, massive black holes, and MBHBs whose emission is modeled based on their star formation histories and gas accretion physics. We predict an average of \approx 33 CGW detections, with signal-to-noise ratios S/N > 5. The detected MBHBs are typically at z < 0.5, with masses of \sim 3 \times 10^{9} M_{\odot}, mass ratios > 0.6 and eccentricities \lesssim 0.2. In terms of EM counterparts, we find less than 15% of these systems to be connected with an AGN detectable by LSST, while their host galaxies are easily detectable ( < 23 mag) massive ( M_{\star} > 10^{11} M_{\odot}) ellipticals with typical star formation rates (10^{-15} yr^{-1} < sSRF < 10^{-10} yr^{-1}). Although the CGW-EM counterpart association is complicated by poor sky localization (only 35% of these CGWs are localized within \rm 100\, deg^2), the number of galaxy host candidates can be considerably reduced (thousands to tens) by applying priors based on the galaxy-MBH correlations. However, picking the actual host among these candidates is highly non-trivial, as they occupy a similar region in any optical color-color diagram. Our findings highlight the considerable challenges entailed in opening the low-frequency multimessenger GW sky.
We present a new dynamical measurement of the supermassive black hole mass and intrinsic shape of the stellar halo of the massive radio galaxy NGC 315 as part of the MASSIVE survey. High signal-to-noise ratio spectra from integral-field spectrographs at the Gemini and McDonald Observatories provide stellar kinematic measurements in 304 spatial bins from the central {\sim}0.3'' out to 30''. Using {\sim} 2300 kinematic constraints, we perform triaxial stellar orbit modeling with the TriOS code and search over {\sim}15,000 galaxy models with a Bayesian scheme to simultaneously measure six mass and intrinsic shape parameters. NGC 315 is triaxial and highly prolate, with middle-to-long and short-to-long axis ratios of p=0.854 and q=0.833 and a triaxiality parameter of T=0.89. The black hole mass inferred from our stellar kinematics is M_\mathrm{BH} = \left(3.0 {\pm} 0.3\right) {\times} 10^{9}\ M_\odot, which is higher than M_\mathrm{BH}=(1.96^{+0.30}_{-0.13}) {\times} 10^{9} M_\odot inferred from CO kinematics (scaled to our distance). When the seven galaxies with M_\mathrm{BH} measurements from both stellar and CO kinematics are compared, we find an intrinsic scatter of 0.28 dex in M_\mathrm{BH} from the two tracers and do not detect statistically significant biases between the two methods in the current data. The implied black hole shadow size ({\approx} 4.7\, \mu{\rm as}) and the relatively high millimeter flux of NGC 315 makes this galaxy a prime candidate for future horizon-size imaging studies.