-
Localized Deviations from the CO-PAH Relation in PHANGS-JWST Galaxies: Faint PAH Emission or Elevated CO Emissivity?
Authors:
Jaeyeon Kim,
Adam K. Leroy,
Karin Sandstrom,
Sharon E. Meidt,
Yu-Hsuan Teng,
Miguel Querejeta,
Eva Schinnerer,
Susan E. Clark,
Ryan Chown,
Simon C. O. Glover,
Daniel A. Dale,
Dalya Baron,
Jessica Sutter,
Ashley T. Barnes,
Jakob den Brok,
Rupali Chandar,
I-Da Chiang,
Oleg V. Egorov,
Kathryn Grasha,
Ralf S. Klessen,
Kathryn Kreckel,
Eric W. Koch,
Hannah Koziol,
Lukas Neumann,
Hsi-An Pan
, et al. (3 additional authors not shown)
Abstract:
Polycyclic aromatic hydrocarbon (PAH) emission is widely used to trace the distribution of molecular gas in the interstellar medium (ISM), exhibiting a tight correlation with CO(2-1) emission across nearby galaxies. Using PHANGS-JWST and PHANGS-ALMA data, we identify localized regions where this correlation fails, with CO flux exceeding that predicted from 7.7$μ$m PAH emission by more than an orde…
▽ More
Polycyclic aromatic hydrocarbon (PAH) emission is widely used to trace the distribution of molecular gas in the interstellar medium (ISM), exhibiting a tight correlation with CO(2-1) emission across nearby galaxies. Using PHANGS-JWST and PHANGS-ALMA data, we identify localized regions where this correlation fails, with CO flux exceeding that predicted from 7.7$μ$m PAH emission by more than an order of magnitude. These outlier regions are found in 20 out of 70 galaxies and are mostly located in galaxy centers and bars, without signs of massive star formation. We explore two scenarios to explain the elevated CO-to-PAH ratios, which can either be due to suppressed PAH emission or enhanced CO emissivity. We examine PAH emission in other bands (3.3$μ$m and 11.3$μ$m) and the dust continuum dominated bands (10$μ$m and 21$μ$m), finding consistently high CO-to-PAH (or CO-to-dust continuum) emission ratios, suggesting that 7.7$μ$m PAH emission is not particularly suppressed. In some outlier regions, PAH size distributions and spectral energy distribution of the illuminating radiation differ slightly compared to nearby control regions with normal CO-to-PAH ratios, though without a consistent trend. We find that the outlier regions generally show higher CO velocity dispersions ($Δv_{\mathrm{CO}}$). This increase in $Δv_{\mathrm{CO}}$ may lower CO optical depth and raise its emissivity for a given gas mass. Our results favor a scenario where shear along the bar lanes and shocks at the tips of the bar elevate CO emissivity, leading to the breakdown of the CO-PAH correlation.
△ Less
Submitted 18 November, 2025;
originally announced November 2025.
-
ReplicationBench: Can AI Agents Replicate Astrophysics Research Papers?
Authors:
Christine Ye,
Sihan Yuan,
Suchetha Cooray,
Steven Dillmann,
Ian L. V. Roque,
Dalya Baron,
Philipp Frank,
Sergio Martin-Alvarez,
Nolan Koblischke,
Frank J Qu,
Diyi Yang,
Risa Wechsler,
Ioana Ciuca
Abstract:
Frontier AI agents show increasing promise as scientific research assistants, and may eventually be useful for extended, open-ended research workflows. However, in order to use agents for novel research, we must first assess the underlying faithfulness and correctness of their work. To evaluate agents as research assistants, we introduce ReplicationBench, an evaluation framework that tests whether…
▽ More
Frontier AI agents show increasing promise as scientific research assistants, and may eventually be useful for extended, open-ended research workflows. However, in order to use agents for novel research, we must first assess the underlying faithfulness and correctness of their work. To evaluate agents as research assistants, we introduce ReplicationBench, an evaluation framework that tests whether agents can replicate entire research papers drawn from the astrophysics literature. Astrophysics, where research relies heavily on archival data and computational study while requiring little real-world experimentation, is a particularly useful testbed for AI agents in scientific research. We split each paper into tasks which require agents to replicate the paper's core contributions, including the experimental setup, derivations, data analysis, and codebase. Each task is co-developed with the original paper authors and targets a key scientific result, enabling objective evaluation of both faithfulness (adherence to original methods) and correctness (technical accuracy of results). ReplicationBench is extremely challenging for current frontier language models: even the best-performing language models score under 20%. We analyze ReplicationBench trajectories in collaboration with domain experts and find a rich, diverse set of failure modes for agents in scientific research. ReplicationBench establishes the first benchmark of paper-scale, expert-validated astrophysics research tasks, reveals insights about agent performance generalizable to other domains of data-driven science, and provides a scalable framework for measuring AI agents' reliability in scientific research.
△ Less
Submitted 23 November, 2025; v1 submitted 28 October, 2025;
originally announced October 2025.
-
Polycyclic aromatic hydrocarbons destruction in star-forming regions across 42 nearby galaxies
Authors:
Oleg V. Egorov,
Adam K. Leroy,
Karin Sandstrom,
Kathryn Kreckel,
Dalya Baron,
Francesco Belfiore,
Ryan Chown,
Jessica Sutter,
Médéric Boquien,
Mar Canal i Saguer,
Enrico Congiu,
Daniel A. Dale,
Evgeniya Egorova,
Michael Huber,
Jing Li,
Thomas G. Williams,
Jérémy Chastenet,
I-Da Chiang,
Ivan Gerasimov,
Hamid Hassani,
Hwihyun Kim,
Hannah Koziol,
Janice C. Lee,
Rebecca L. McClain,
José Eduardo Méndez Delgado
, et al. (8 additional authors not shown)
Abstract:
Polycyclic aromatic hydrocarbons (PAHs) are widespread in the interstellar medium (ISM) of Solar metallicity galaxies, where they play a critical role in ISM heating, cooling, and reprocessing stellar radiation. The PAH fraction, the abundance of PAHs relative to total dust mass, is a key parameter in ISM physics. Using JWST and MUSE observations of 42 galaxies from the PHANGS survey, we analyze t…
▽ More
Polycyclic aromatic hydrocarbons (PAHs) are widespread in the interstellar medium (ISM) of Solar metallicity galaxies, where they play a critical role in ISM heating, cooling, and reprocessing stellar radiation. The PAH fraction, the abundance of PAHs relative to total dust mass, is a key parameter in ISM physics. Using JWST and MUSE observations of 42 galaxies from the PHANGS survey, we analyze the PAH fraction in over 17 000 H II regions spanning a gas-phase oxygen abundance of 12+log(O/H) = 8.0-8.8 (Z ~ 0.2-1.3 Zsun), and ~400 isolated supernova remnants (SNRs). We find a significantly lower PAH fraction toward H II regions compared to a reference sample of diffuse ISM areas at matched metallicity. At 12+log(O/H) > 8.2, the PAH fraction toward H II regions is strongly anti-correlated with the local ionization parameter, suggesting that PAH destruction is correlated with ionized gas and/or hydrogen-ionizing UV radiation. At lower metallicities, the PAH fraction declines steeply in both H II regions and the diffuse ISM, likely reflecting less efficient PAH formation in metal-poor environments. Carefully isolating dust emission from the vicinity of optically-identified supernova remnants, we see evidence for selective PAH destruction from measurements of lower PAH fractions, which is, however, indistinguishable at ~50 pc scales. Overall, our results point to ionizing radiation as the dominant agent of PAH destruction within H II regions, with metallicity playing a key role in their global abundance in galaxies.
△ Less
Submitted 17 September, 2025;
originally announced September 2025.
-
Oxyster: A Circumgalactic Low-ionization Oxygen Nebula next to a Starburst Galaxy at $z\sim1$
Authors:
Pengjun Lu,
Mingyu Li,
Dalya Baron,
Minghao Yue,
Song Huang,
Zheng Cai
Abstract:
Extended emission line nebulae around galaxies or active galactic nuclei (AGNs) provide a unique window to investigate the galactic ecosystem through the circumgalactic medium (CGM). Using Subaru Hyper-Suprime Cam narrow-band imaging and spectroscopic follow-up, we serendipitously discover "Oxyster" - a large ionized nebula next to an interacting starburst galaxy at $z=0.924$. The nebula is traced…
▽ More
Extended emission line nebulae around galaxies or active galactic nuclei (AGNs) provide a unique window to investigate the galactic ecosystem through the circumgalactic medium (CGM). Using Subaru Hyper-Suprime Cam narrow-band imaging and spectroscopic follow-up, we serendipitously discover "Oxyster" - a large ionized nebula next to an interacting starburst galaxy at $z=0.924$. The nebula is traced by extended [OII]3726,3729 ($\sim 30$ kpc) and [OIII]5007 ($\sim 20$ kpc) emission lines. On the nebula luminosity-size plane, Oxyster surpasses the extended narrow-line regions around low-$z$ AGNs, resembling a higher-$z$ analog of "Hanny's Voorwerp". However, its uniformly low [OIII]/[OII] ratio (O32) sets it apart from typical AGN light echoes. For the host galaxy, HST and JWST images reveal a disturbed red disk galaxy with a single blue spiral "arm". Spectral energy distribution (SED) fitting suggests the $2-6\times 10^{10} ~\rm M_{\odot}$ host galaxy sits above the star-forming main sequence with an ongoing starburst, especially in the "arm", and have $<5\%$ luminosity contribution from AGN, consistent with X-ray non-detection and radio continuum. Standard photoionization and shock models struggle to explain simultaneously Oxyster's emission line luminosities, low O32 ratio, and the non-detection of H$β$ line. A plausible explanation could involve the combination of a recent ($<10^8$ yrs) starburst and a low-luminosity AGN ($L_{\rm{bol}} \sim 1\times10^{42}$ erg/s). While Oxyster's nature awaits future investigation, its discovery highlights the potential of ground-based narrow-band imaging to uncover extended emission line nebulae around non-AGN systems, opening new avenues for studying the CGM of normal galaxies in the early Universe.
△ Less
Submitted 3 November, 2025; v1 submitted 15 April, 2025;
originally announced April 2025.
-
Empirical SED Templates for Star Clusters Observed with HST and JWST: No Strong PAH or IR Dust Emission after Five Myr
Authors:
Bradley C. Whitmore,
Rupali Chandar,
Janice C. Lee,
Kiana F. Henny,
M. Jimena Rodriguez,
Dalya Baron,
F. Bigiel,
Mederic Boquien,
Melanie Chevance,
Ryan Chown,
Daniel A. Dale,
Matthew Floyd,
Kathryn Grasha,
Simon C. O. Glover,
Oleg Gnedin,
Hamid Hassani,
Remy Indebetouw,
Anand Utsav Kapoor,
Kirsten L. Larson,
Adam K. Leroy,
Daniel Maschmann,
Fabian Scheuermann,
Jessica Sutter,
Eva Schinnerer,
Sumit K. Sarbadhicary
, et al. (3 additional authors not shown)
Abstract:
JWST observations, when combined with HST data, promise to improve age estimates of star clusters in nearby spiral galaxies. However, feedback from young cluster stars pushes out the natal gas and dust, making cluster formation and evolution a challenge to model. Here, we use JWST + HST observations of the nearby spiral galaxy NGC 628 to produce spectral energy distribution (SED) templates of comp…
▽ More
JWST observations, when combined with HST data, promise to improve age estimates of star clusters in nearby spiral galaxies. However, feedback from young cluster stars pushes out the natal gas and dust, making cluster formation and evolution a challenge to model. Here, we use JWST + HST observations of the nearby spiral galaxy NGC 628 to produce spectral energy distribution (SED) templates of compact star clusters spanning 275 nm through 21 μm. These preliminary SEDs capture the cluster stars and associated gas and dust within radii of 0.12" to 0.67" (corresponding to 6 to 33 pc at the distance of NGC 628). One important finding is that the SEDs of 1, 2, 3, and 4 Myr clusters can be differentiated in the infrared. Another is that in 80-90% of the cases we study, the PAH and H_alpha emission track one another, with the dust responsible for the 3.3 μm PAH emission largely removed by 4 Myr, consistent with pre-supernova stellar feedback acting quickly on the surrounding gas and dust. Nearly-embedded cluster candidates have infrared SEDs which are quite similar to optically visible 1 to 3 Myr clusters. In nearly all cases we find there is a young star cluster within a few tenths of an arcsec (10 - 30 pc) of the nearly embedded cluster, suggesting the formation of the cluster was triggered by its presence. The resulting age estimates from the empirical templates are compatible both with dynamical estimates based on CO superbubble expansion velocities, and the TODDLERS models which track spherical evolution of homogeneous gas clouds around young stellar clusters.
△ Less
Submitted 26 March, 2025; v1 submitted 22 March, 2025;
originally announced March 2025.
-
Linking stellar populations to HII regions across nearby galaxies. II. Infrared Reprocessed and UV Direct Radiation Pressure in HII Regions
Authors:
Debosmita Pathak,
Adam Leroy,
Todd Thompson,
Laura Lopez,
Ashley Barnes,
Daniel Dale,
Ian Blackstone,
Simon C. O. Glover,
Shyam Menon,
Jessica Sutter,
Thomas Williams,
Dalya Baron,
Francesco Belfiore,
Frank Bigiel,
Alberto Bolatto,
Mederic Boquien,
Rupali Chandar,
Mélanie Chevance,
Ryan Chown,
Kathryn Grasha,
Brent Groves,
Ralf Klessen,
Kathryn Kreckel,
Jing Li,
José Méndez-Delgado
, et al. (5 additional authors not shown)
Abstract:
Radiation pressure is a key mechanism by which stellar feedback disrupts molecular clouds and drives HII region expansion. This includes direct radiation pressure exerted by UV photons on dust grains, pressure associated with photoionization, and infrared (IR) radiation pressure on grains due to dust-reprocessed IR photons. We present a new method that combines high resolution mid-IR luminosities…
▽ More
Radiation pressure is a key mechanism by which stellar feedback disrupts molecular clouds and drives HII region expansion. This includes direct radiation pressure exerted by UV photons on dust grains, pressure associated with photoionization, and infrared (IR) radiation pressure on grains due to dust-reprocessed IR photons. We present a new method that combines high resolution mid-IR luminosities from JWST-MIRI, optical attenuation and nebular line measurements from VLT-MUSE, and HST H$α$-based region sizes to estimate the strength of radiation pressure in $\approx 18,000$ HII regions across 19 nearby star-forming galaxies. This is the most extensive and direct estimate of these terms beyond the Local Group to date. In the disks of galaxies, we find that the total reprocessed IR pressure is on average 5% of the direct UV radiation pressure. This fraction rises to 10% in galaxy centers. We expect reprocessed IR radiation pressure to dominate over UV radiation pressure in regions where $L_{\rm F2100W}/L_{\rm Hα}^{\rm corr} \gtrsim 75$. Radiation pressure due to H ionizations is lower than pressure on dust in our sample, but appears likely to dominate the radiation pressure budget in dwarf galaxies similar to the Small Magellanic Cloud. The contribution from all radiation pressure terms appears to be subdominant compared to thermal pressure from ionized gas, reinforcing the view that radiation pressure is most important in compact, heavily embedded, and young regions.
△ Less
Submitted 25 March, 2025; v1 submitted 31 January, 2025;
originally announced February 2025.
-
PAH Feature Ratios Around Stellar Clusters and Associations in 19 Nearby Galaxies
Authors:
Daniel A. Dale,
Gabrielle B. Graham,
Ashley T. Barnes,
Dalya Baron,
Frank Bigiel,
Médéric Boquien,
Rupali Chandar,
Jérémy Chastenet,
Ryan Chown,
Oleg V. Egorov,
Simon C. O. Glover,
Lindsey Hands,
Kiana F. Henny,
Remy Indebetouw,
Ralf S. Klessen,
Kirsten L. Larson,
Janice C. Lee,
Adam K. Leroy,
Daniel Maschmann,
Debosmita Pathak,
M. Jimena Rodríguez,
Erik Rosolowsky,
Karin Sandstrom,
Eva Schinnerer,
Jessica Sutter
, et al. (5 additional authors not shown)
Abstract:
We present a comparison of observed polycyclic aromatic hydrocarbon (PAH) feature ratios in 19 nearby galaxies with a grid of theoretical expectations for near- and mid-infrared dust emission. The PAH feature ratios are drawn from Cycle 1 JWST observations and are measured for 7224 stellar clusters and 29176 stellar associations for which we have robust ages and mass estimates from HST five-band p…
▽ More
We present a comparison of observed polycyclic aromatic hydrocarbon (PAH) feature ratios in 19 nearby galaxies with a grid of theoretical expectations for near- and mid-infrared dust emission. The PAH feature ratios are drawn from Cycle 1 JWST observations and are measured for 7224 stellar clusters and 29176 stellar associations for which we have robust ages and mass estimates from HST five-band photometry. Though there are galaxy-to-galaxy variations, the observed PAH feature ratios largely agree with the theoretical models, particularly those that are skewed toward more ionized and larger PAH size distributions. For each galaxy we also extract PAH feature ratios for 200 pc-wide circular regions in the diffuse interstellar medium, which serve as a non-cluster/association control sample. Compared to what we find for stellar clusters and associations, the 3.3um/7.7um and 3.3um/11.3um ratios from the diffuse interstellar medium are $\sim 0.10-0.15$ dex smaller. When the observed PAH feature ratios are compared to the radiation field hardness as probed by the [OIII]/H$β$ ratio, we find anti-correlations for nearly all galaxies in the sample. These results together suggest that the PAH feature ratios are driven by the shape intensity of the radiation field, and that the smallest PAHs -- observed via JWST F335M imaging -- are increasingly 'processed' or destroyed in regions with the most intense and hard radiation fields.
△ Less
Submitted 17 January, 2025;
originally announced January 2025.
-
Polycyclic Aromatic Hydrocarbon and CO(2-1) Emission at 50-150 pc Scales in 70 Nearby Galaxies
Authors:
Ryan Chown,
Adam K. Leroy,
Karin Sandstrom,
Jeremy Chastenet,
Jessica Sutter,
Eric W. Koch,
Hannah B. Koziol,
Lukas Neumann,
Jiayi Sun,
Thomas G. Williams,
Dalya Baron,
Gagandeep S. Anand,
Ashley T. Barnes,
Zein Bazzi,
Francesco Belfiore,
Alberto Bolatto,
Mederic Boquien,
Frank Bigiel,
Yixian Cao,
Melanie Chevance,
Dario Colombo,
Daniel A. Dale,
Jakob den Brok,
Oleg V. Egorov,
Cosima Eibensteiner
, et al. (22 additional authors not shown)
Abstract:
Combining Atacama Large Millimeter/sub-millimeter Array CO(2-1) mapping and JWST near- and mid-infrared imaging, we characterize the relationship between CO(2-1) and polycyclic aromatic hydrocarbon (PAH) emission at ~100 pc resolution in 70 nearby star-forming galaxies. Leveraging a new Cycle 2 JWST treasury program targeting nearby galaxies, we expand the sample size by more than an order of magn…
▽ More
Combining Atacama Large Millimeter/sub-millimeter Array CO(2-1) mapping and JWST near- and mid-infrared imaging, we characterize the relationship between CO(2-1) and polycyclic aromatic hydrocarbon (PAH) emission at ~100 pc resolution in 70 nearby star-forming galaxies. Leveraging a new Cycle 2 JWST treasury program targeting nearby galaxies, we expand the sample size by more than an order of magnitude compared to previous ~100 pc resolution CO-PAH comparisons. Focusing on regions of galaxies where most of the gas is likely to be molecular, we find strong correlations between CO(2-1) and 3.3 um, 7.7 um, and 11.3 um PAH emission, estimated from JWST's F335M, F770W, and F1130W filters. We derive power law relations between CO(2-1) and PAH emission, which have indices in the range 0.8-1.3, implying relatively weak variations in the observed CO-to-PAH ratios across the regions that we study. We find that CO-to-PAH ratios and scaling relationships near HII regions are similar to those in diffuse sight lines. The main difference between the two types of regions is that sight lines near HII regions show higher intensities in all tracers. Galaxy centers, on the other hand, show higher overall intensities and enhanced CO-to-PAH ratios compared to galaxy disks. Individual galaxies show 0.19 dex scatter in the normalization of CO at fixed I_PAH, and this normalization anti-correlates with specific star formation rate (sSFR) and correlates with stellar mass. We provide a prescription that accounts for these galaxy-to-galaxy variations and represents our best current empirical predictor to estimate CO(2-1) intensity from PAH emission, which allows one to take advantage of JWST's excellent sensitivity and resolution to trace cold gas.
△ Less
Submitted 19 March, 2025; v1 submitted 7 October, 2024;
originally announced October 2024.
-
PHANGS-ML: the universal relation between PAH band and optical line ratios across nearby star-forming galaxies
Authors:
Dalya Baron,
Karin Sandstrom,
Jessica Sutter,
Hamid Hassani,
Brent Groves,
Adam Leroy,
Eva Schinnerer,
Médéric Boquien,
Matilde Brazzini,
Jérémy Chastenet,
Daniel Dale,
Oleg Egorov,
Simon Glover,
Ralf Klessen,
Debosmita Pathak,
Erik Rosolowsky,
Frank Bigiel,
Mélanie Chevance,
Kathryn Grasha,
Annie Hughes,
J. Eduardo Méndez-Delgado,
Jérôme Pety,
Thomas Williams,
Stephen Hannon,
Sumit Sarbadhicary
Abstract:
The structure and chemistry of the dusty interstellar medium (ISM) are shaped by complex processes that depend on the local radiation field, gas composition, and dust grain properties. Of particular importance are Polycyclic Aromatic Hydrocarbons (PAHs), which emit strong vibrational bands in the mid-infrared, and play a key role in the ISM energy balance. We recently identified global correlation…
▽ More
The structure and chemistry of the dusty interstellar medium (ISM) are shaped by complex processes that depend on the local radiation field, gas composition, and dust grain properties. Of particular importance are Polycyclic Aromatic Hydrocarbons (PAHs), which emit strong vibrational bands in the mid-infrared, and play a key role in the ISM energy balance. We recently identified global correlations between PAH band and optical line ratios across three nearby galaxies, suggesting a connection between PAH heating and gas ionization throughout the ISM. In this work, we perform a census of the PAH heating -- gas ionization connection using $\sim$700,000 independent pixels that probe scales of 40--150 pc in nineteen nearby star-forming galaxies from the PHANGS survey. We find a universal relation between $\log$PAH(11.3 \mic/7.7 \mic) and $\log$([SII]/H$α$) with a slope of $\sim$0.2 and a scatter of $\sim$0.025 dex. The only exception is a group of anomalous pixels that show unusually high (11.3 \mic/7.7 \mic) PAH ratios in regions with old stellar populations and high starlight-to-dust emission ratios. Their mid-infrared spectra resemble those of elliptical galaxies. AGN hosts show modestly steeper slopes, with a $\sim$10\% increase in PAH(11.3 \mic/7.7 \mic) in the diffuse gas on kpc scales. This universal relation implies an emerging simplicity in the complex ISM, with a sequence that is driven by a single varying property: the spectral shape of the interstellar radiation field. This suggests that other properties, such as gas-phase abundances, gas ionization parameter, and grain charge distribution, are relatively uniform in all but specific cases.
△ Less
Submitted 3 October, 2024;
originally announced October 2024.
-
Metallicity calibrations based on auroral lines from PHANGS-MUSE data
Authors:
Matilde Brazzini,
Francesco Belfiore,
Michele Ginolfi,
Brent Groves,
Kathryn Kreckel,
Ryan J. Rickards Vaught,
Dalya Baron,
Frank Bigiel,
Guillermo A. Blanc,
Daniel A. Dale,
Kathryn Grasha,
Eric Habjan,
Ralf S. Klessen,
J. Eduardo Méndez-Delgado,
Karin Sandstrom,
Thomas G. Williams
Abstract:
We present a chemical analysis of selected HII regions from the PHANGS-MUSE nebular catalogue. Our intent is to empirically re-calibrate strong-line diagnostics of gas-phase metallicity, applicable across a wide range of metallicities within nearby star-forming galaxies. To ensure reliable measurements of auroral line fluxes, we carried out a new spectral fitting procedure whereby only restricted…
▽ More
We present a chemical analysis of selected HII regions from the PHANGS-MUSE nebular catalogue. Our intent is to empirically re-calibrate strong-line diagnostics of gas-phase metallicity, applicable across a wide range of metallicities within nearby star-forming galaxies. To ensure reliable measurements of auroral line fluxes, we carried out a new spectral fitting procedure whereby only restricted wavelength regions around the emission lines of interest are taken into account: this assures a better fit for the stellar continuum. No prior cuts to nebulae luminosity were applied to limit biases in auroral line detections. Ionic abundances of O+, O++, N+, S+, and S++ were estimated by applying the direct method. We integrated the selected PHANGS-MUSE sample with other existing auroral line catalogues, appropriately re-analysed to obtain a homogeneous dataset. This was used to derive strong-line diagnostic calibrations that span from 12+log(O/H) = 7.5 to 8.8. We investigate their dependence on the ionisation parameter and conclude that it is likely the primary cause of the significant scatter observed in these diagnostics. We apply our newly calibrated strong-line diagnostics to the total sample of HII regions from the PHANGS-MUSE nebular catalogue, and we exploit these indirect metallicity estimates to study the radial metallicity gradient within each of the 19 galaxies of the sample. We compare our results with the literature and find good agreement, validating our procedure and findings. With this paper, we release the full catalogue of auroral and nebular line fluxes for the selected HII regions from the PHANGS-MUSE nebular catalogue. This is the first catalogue of direct chemical abundance measurements carried out with PHANGS-MUSE data.
△ Less
Submitted 30 September, 2024;
originally announced October 2024.
-
PHANGS-ML: dissecting multiphase gas and dust in nearby galaxies using machine learning
Authors:
Dalya Baron,
Karin M. Sandstrom,
Erik Rosolowsky,
Oleg V. Egorov,
Ralf S. Klessen,
Adam K. Leroy,
Médéric Boquien,
Eva Schinnerer,
Francesco Belfiore,
Brent Groves,
Jérémy Chastenet,
Daniel A. Dale,
Guillermo A. Blanc,
José E. Méndez-Delgado,
Eric W. Koch,
Kathryn Grasha,
Mélanie Chevance,
David A. Thilker,
Dario Colombo,
Thomas G. Williams,
Debosmita Pathak,
Jessica Sutter,
Toby Brown,
John F. Wu,
J. E. G. Peek
, et al. (3 additional authors not shown)
Abstract:
The PHANGS survey uses ALMA, HST, VLT, and JWST to obtain an unprecedented high-resolution view of nearby galaxies, covering millions of spatially independent regions. The high dimensionality of such a diverse multi-wavelength dataset makes it challenging to identify new trends, particularly when they connect observables from different wavelengths. Here we use unsupervised machine learning algorit…
▽ More
The PHANGS survey uses ALMA, HST, VLT, and JWST to obtain an unprecedented high-resolution view of nearby galaxies, covering millions of spatially independent regions. The high dimensionality of such a diverse multi-wavelength dataset makes it challenging to identify new trends, particularly when they connect observables from different wavelengths. Here we use unsupervised machine learning algorithms to mine this information-rich dataset to identify novel patterns. We focus on three of the PHANGS-JWST galaxies, for which we extract properties pertaining to their stellar populations; warm ionized and cold molecular gas; and Polycyclic Aromatic Hydrocarbons (PAHs), as measured over 150 pc-scale regions. We show that we can divide the regions into groups with distinct multiphase gas and PAH properties. In the process, we identify previously-unknown galaxy-wide correlations between PAH band and optical line ratios and use our identified groups to interpret them. The correlations we measure can be naturally explained in a scenario where the PAHs and the ionized gas are exposed to different parts of the same radiation field that varies spatially across the galaxies. This scenario has several implications for nearby galaxies: (i) The uniform PAH ionized fraction on 150 pc scales suggests significant self-regulation in the ISM, (ii) the PAH 11.3/7.7 \mic~ band ratio may be used to constrain the shape of the non-ionizing far-ultraviolet to optical part of the radiation field, and (iii) the varying radiation field affects line ratios that are commonly used as PAH size diagnostics. Neglecting this effect leads to incorrect or biased PAH sizes.
△ Less
Submitted 6 February, 2024;
originally announced February 2024.
-
Not so windy after all: MUSE disentangles AGN-driven winds from merger-induced flows in rapidly-transitioning galaxies
Authors:
Dalya Baron,
Hagai Netzer,
Dieter Lutz,
Ric I. Davies,
J. Xavier Prochaska
Abstract:
Post-starburst galaxies are believed to be in a rapid transition between major merger starbursts and quiescent ellipticals, where AGN feedback is suggested as one of the processes responsible for the quenching. To study the role of AGN feedback, we constructed a sample of post-starburst candidates with AGN and indications of ionized outflows. We use MUSE/VLT observations to resolve the properties…
▽ More
Post-starburst galaxies are believed to be in a rapid transition between major merger starbursts and quiescent ellipticals, where AGN feedback is suggested as one of the processes responsible for the quenching. To study the role of AGN feedback, we constructed a sample of post-starburst candidates with AGN and indications of ionized outflows. We use MUSE/VLT observations to resolve the properties of the stars and multi-phased gas in five of them. All the galaxies show signatures of interaction/merger in their stellar or gas properties, with some galaxies at an early stage of interaction with companions at distances $\sim$50 kpc, suggesting that optical post-starburst signatures may be present well before the final starburst and coalescence. We detect narrow and broad kinematic components in multiple transitions in all the galaxies. Our detailed analysis of their kinematics and morphology suggests that, contrary to our expectation, the properties of the broad kinematic components are inconsistent with AGN-driven winds in 3 out of 5 galaxies. The two exceptions are also the only galaxies in which spatially-resolved NaID P-Cygni profiles are detected. In some cases, the observations are more consistent with interaction-induced galactic-scale flows, an often overlooked process. These observations raise the question of how to interpret broad kinematic components in interacting and perhaps also in active galaxies, in particular when spatially-resolved observations are not available or cannot rule out merger-induced galactic-scale motions. We suggest that NaID P-Cygni profiles are more effective outflow tracers, and use them to estimate the energy that is carried by the outflow.
△ Less
Submitted 17 January, 2024;
originally announced January 2024.
-
Star formation and molecular gas properties of post-starburst galaxies
Authors:
Dalya Baron,
Hagai Netzer,
K. Decker French,
Dieter Lutz,
Ric I. Davies,
J. Xavier Prochaska
Abstract:
Post-starburst galaxies are believed to be in a rapid transition between major merger starbursts and quiescent ellipticals. Their optical spectrum is dominated by A-type stars, suggesting a starburst that was quenched recently. While optical observations suggest little ongoing star formation, some have been shown to host significant molecular gas reservoirs. This led to the suggestion that gas dep…
▽ More
Post-starburst galaxies are believed to be in a rapid transition between major merger starbursts and quiescent ellipticals. Their optical spectrum is dominated by A-type stars, suggesting a starburst that was quenched recently. While optical observations suggest little ongoing star formation, some have been shown to host significant molecular gas reservoirs. This led to the suggestion that gas depletion is not required to end the starburst, and that star formation is suppressed by other processes. We present NOEMA CO(1-0) observations of 15 post-starburst galaxies with emission lines consistent with Active Galactic Nucleus (AGN) photoionization. We collect post-starburst candidates with molecular gas measurements from the literature, with some classified as classical E+A, while others with line ratios consistent with AGN and/or shock ionization. Using far-infrared observations, we show that systems that were reported to host exceptionally-large molecular gas reservoirs host in fact obscured star formation, with some systems showing star formation rates comparable to ULIRGs. Among E+A galaxies with molecular gas measurements, 7 out of 26 (26\%) host obscured starbursts. Using far-infrared observations, post-starburst candidates show similar SFR-$M_{\mathrm{H_2}}$ and Kennicutt-Schmidt relations to those observed in star-forming and starburst galaxies. In particular, there is no need to hypothesize star formation quenching by processes other than the consumption of molecular gas by star formation. The combination of optical, far-infrared, and CO observations indicates that some regions within these galaxies have been recently quenched, while others are still forming stars in highly obscured regions. All this calls into question the traditional interpretation of such galaxies.
△ Less
Submitted 28 May, 2023; v1 submitted 25 April, 2022;
originally announced April 2022.
-
Multi-phase outflows in post starburst E+A galaxies -- I. General sample properties and the prevalence of obscured starbursts
Authors:
Dalya Baron,
Hagai Netzer,
Dieter Lutz,
J. Xavier Prochaska,
Ric I. Davies
Abstract:
E+A galaxies are believed to be a short phase connecting major merger ULIRGs with red and dead elliptical galaxies. Their optical spectrum suggests a massive starburst that was quenched abruptly, and their bulge-dominated morphologies with tidal tails suggest that they are merger remnants. AGN-driven winds are believed to be one of the processes responsible for the sudden quenching of star formati…
▽ More
E+A galaxies are believed to be a short phase connecting major merger ULIRGs with red and dead elliptical galaxies. Their optical spectrum suggests a massive starburst that was quenched abruptly, and their bulge-dominated morphologies with tidal tails suggest that they are merger remnants. AGN-driven winds are believed to be one of the processes responsible for the sudden quenching of star formation and for the expulsion and/or destruction of the remaining molecular gas. Little is known about AGN-driven winds in this short-lived phase. In this paper we present the first and unique sample of post starburst galaxy candidates with AGN that show indications of ionized outflows in their optical emission lines. Using IRAS-FIR observations, we study the star formation in these systems and find that many systems selected to have post starburst signatures in their optical spectrum are in fact obscured starbursts. Using SDSS spectroscopy, we study the stationary and outflowing ionized gas. We also detect neutral gas outflows in 40\% of the sources with mass outflow rates 10--100 times more massive than in the ionized phase. The mean mass outflow rate and kinetic power of the ionized outflows in our sample ($\dot{M}\sim 1\, \mathrm{M_{\odot}/yr}$, $\dot{E}\sim 10^{41}\, \mathrm{erg/sec}$) are larger than those derived for active galaxies of similar AGN luminosity and stellar mass. For the neutral outflow ($\dot{M}\sim 10\, \mathrm{M_{\odot}/yr}$, $\dot{E}\sim 10^{42}\, \mathrm{erg/sec}$), their mean is smaller than that observed in (U)LIRGs with and without AGN.
△ Less
Submitted 3 October, 2021; v1 submitted 17 May, 2021;
originally announced May 2021.
-
AGN-driven outflows and the AGN feedback efficiency in young radio galaxies
Authors:
F. Santoro,
C. Tadhunter,
D. Baron,
R. Morganti,
J. Holt
Abstract:
Active galactic nuclei (AGN) feedback operated by the expansion of radio jets can play a crucial role in driving gaseous outflows on galaxy scales. Galaxies hosting young radio AGN, whose jets are in the first phases of expansion through the surrounding interstellar medium (ISM), are the ideal targets to probe the energetic significance of this mechanism. In this paper, we characterise the warm io…
▽ More
Active galactic nuclei (AGN) feedback operated by the expansion of radio jets can play a crucial role in driving gaseous outflows on galaxy scales. Galaxies hosting young radio AGN, whose jets are in the first phases of expansion through the surrounding interstellar medium (ISM), are the ideal targets to probe the energetic significance of this mechanism. In this paper, we characterise the warm ionised gas outflows in a sample of nine young radio sources from the 2Jy sample, combining X-shooter spectroscopy and Hubble Space Telescope (HST) imaging data. We find that the warm outflows have similar radial extents (~0.06-2 kpc) as radio sources, consistent with the idea that `jet mode' AGN feedback is the dominant driver of the outflows detected in young radio galaxies. Exploiting the broad spectral coverage of the X-shooter data, we have used the ratios of trans-auroral emission lines of [SII] and [OII] to estimate the electron densities, finding that most of the outflows have gas densities ($\log( n_e~cm^{-3})~3-4.8 $), which we speculate could be the result of compression by jet-induced shocks. Combining our estimates of the emission-line luminosities, radii, and densities, we find that the kinetic powers of the warm outflows are a relatively small fraction of the energies available from the accretion of material onto the central supermassive black hole (SMBH), reflecting AGN feedback efficiencies below 1% in most cases. Overall, the warm outflows detected in our sample are strikingly similar to those found in nearby ultraluminous infrared galaxies (ULIRGs), but more energetic and with a high feedback efficiencies on average than the general population of nearby AGN of similar bolometric luminosity; this is likely to reflect a high degree of coupling between the jets and the near-nuclear ISM in the early stages of radio source evolution.
△ Less
Submitted 23 September, 2020;
originally announced September 2020.
-
Extracting the main trend in a dataset: the Sequencer algorithm
Authors:
Dalya Baron,
Brice Ménard
Abstract:
Scientists aim to extract simplicity from observations of the complex world. An important component of this process is the exploration of data in search of trends. In practice, however, this tends to be more of an art than a science. Among all trends existing in the natural world, one-dimensional trends, often called sequences, are of particular interest as they provide insights into simple phenom…
▽ More
Scientists aim to extract simplicity from observations of the complex world. An important component of this process is the exploration of data in search of trends. In practice, however, this tends to be more of an art than a science. Among all trends existing in the natural world, one-dimensional trends, often called sequences, are of particular interest as they provide insights into simple phenomena. However, some are challenging to detect as they may be expressed in complex manners. We present the Sequencer, an algorithm designed to generically identify the main trend in a dataset. It does so by constructing graphs describing the similarities between pairs of observations, computed with a set of metrics and scales. Using the fact that continuous trends lead to more elongated graphs, the algorithm can identify which aspects of the data are relevant in establishing a global sequence. Such an approach can be used beyond the proposed algorithm and can optimize the parameters of any dimensionality reduction technique. We demonstrate the power of the Sequencer using real-world data from astronomy, geology as well as images from the natural world. We show that, in a number of cases, it outperforms the popular t-SNE and UMAP dimensionality reduction techniques. This approach to exploratory data analysis, which does not rely on training nor tuning of any parameter, has the potential to enable discoveries in a wide range of scientific domains. The source code is available on github and we provide an online interface at \url{http://sequencer.org}.
△ Less
Submitted 24 June, 2020;
originally announced June 2020.
-
Multi-phase outflows in post starburst E+A galaxies -- II. A direct connection between the neutral and ionized outflow phases
Authors:
Dalya Baron,
Hagai Netzer,
Ric I. Davies,
J. Xavier Prochaska
Abstract:
Post starburst E+A galaxies are systems that hosted a powerful starburst that was quenched abruptly. Simulations suggest that these systems provide the missing link between major merger ULIRGs and red and dead ellipticals, where AGN feedback is responsible for the expulsion or destruction of the molecular gas. However, many details remain unresolved and little is known about AGN-driven winds in th…
▽ More
Post starburst E+A galaxies are systems that hosted a powerful starburst that was quenched abruptly. Simulations suggest that these systems provide the missing link between major merger ULIRGs and red and dead ellipticals, where AGN feedback is responsible for the expulsion or destruction of the molecular gas. However, many details remain unresolved and little is known about AGN-driven winds in this short-lived phase. We present spatially-resolved IFU spectroscopy with MUSE/VLT of SDSS J124754.95-033738.6, a post starburst E+A galaxy with a recent starburst that started 70 Myrs ago and ended 30 Myrs ago, with a peak SFR of $\sim 150\,\mathrm{M_{\odot}/yr}$. We detect disturbed gas throughout the entire field of view, suggesting triggering by a minor merger. We detect fast-moving multi-phased gas clouds, embedded in a double-cone face-on outflow, which are traced by ionized emission lines and neutral NaID \emph{emission} and absorption lines. We find remarkable similarities between the kinematics, spatial extents, and line luminosities of the ionized and neutral gas phases, and propose a model in which they are part of the same outflowing clouds, which are exposed to both stellar and AGN radiation. Our photoionization model provides consistent ionized line ratios, NaID absorption optical depths and EWs, and dust reddening. Using the model, we estimate, for the first time, the neutral-to-ionized gas mass ratio (about 20), the sodium neutral fraction, and the size of the outflowing clouds. This is one of the best ever observed direct connections between the neutral and ionized outflow phases in AGN.
△ Less
Submitted 19 April, 2020; v1 submitted 9 April, 2020;
originally announced April 2020.
-
Ionized outflows in local luminous AGN: what are the real densities and outflow rates?
Authors:
R. Davies,
D. Baron,
T. Shimizu,
H. Netzer,
L. Burtscher,
P. T. de Zeeuw,
R. Genzel,
E. K. S. Hicks,
M. Koss,
M. -Y. Lin,
D. Lutz,
W. Maciejewski,
F. Müller-Sánchez,
G. Orban de Xivry,
C. Ricci,
R. Riffel,
R. A. Riffel,
D. Rosario,
M. Schartmann,
A. Schnorr-Müller,
J. Shangguan,
A. Sternberg,
E. Sturm,
T. Storchi-Bergmann,
L. Tacconi
, et al. (1 additional authors not shown)
Abstract:
We report on the determination of electron densities, and their impact on the outflow masses and rates, measured in the central few hundred parsecs of 11 local luminous active galaxies. We show that the peak of the integrated line emission in the AGN is significantly offset from the systemic velocity as traced by the stellar absorption features, indicating that the profiles are dominated by outflo…
▽ More
We report on the determination of electron densities, and their impact on the outflow masses and rates, measured in the central few hundred parsecs of 11 local luminous active galaxies. We show that the peak of the integrated line emission in the AGN is significantly offset from the systemic velocity as traced by the stellar absorption features, indicating that the profiles are dominated by outflow. In contrast, matched inactive galaxies are characterised by a systemic peak and weaker outflow wing. We present three independent estimates of the electron density in these AGN, discussing the merits of the different methods. The electron density derived from the [SII] doublet is significantly lower than than that found with a method developed in the last decade using auroral and transauroral lines, as well as a recently introduced method based on the ionization parameter. The reason is that, for gas photoionized by an AGN, much of the [SII] emission arises in an extended partially ionized zone where the implicit assumption that the electron density traces the hydrogen density is invalid. We propose ways to deal with this situation and we derive the associated outflow rates for ionized gas, which are in the range 0.001--0.5 M$_{\odot}$yr$^{-1}$ for our AGN sample. We compare these outflow rates to the relation between $\dot{M}_{out}$ and $L_{AGN}$ in the literature, and argue that it may need to be modified and rescaled towards lower mass outflow rates.
△ Less
Submitted 11 August, 2020; v1 submitted 13 March, 2020;
originally announced March 2020.
-
Exploring the Diversity of Type 1 Active Galactic Nuclei Identified in SDSS-IV/SPIDERS
Authors:
Julien Wolf,
Mara Salvato,
Damien Coffey,
Andrea Merloni,
Johannes Buchner,
Riccardo Arcodia,
Dalya Baron,
Francisco J. Carrera,
Johan Comparat,
Donald P. Schneider,
Kirpal Nandra
Abstract:
We present a statistical analysis of the optical properties of an X-ray selected Type 1 AGN sample, using high signal to-noise ratio ($S/N>20$) spectra of the counterparts of the ROSAT/2RXS sources in the footprint of the SDSS-IV/SPIDERS (Spectroscopic IDentification of eROSITA Sources) programme. The final sample contains 2100 sources. It significantly extends the redshift and luminosity ranges (…
▽ More
We present a statistical analysis of the optical properties of an X-ray selected Type 1 AGN sample, using high signal to-noise ratio ($S/N>20$) spectra of the counterparts of the ROSAT/2RXS sources in the footprint of the SDSS-IV/SPIDERS (Spectroscopic IDentification of eROSITA Sources) programme. The final sample contains 2100 sources. It significantly extends the redshift and luminosity ranges ($\rm z \sim 0.01 - 0.80$ and $\rm L_{0.1-2.4 \,keV}$ ~ $\rm 2.0 \times 10^{41}$ - $\rm 1.0 \times 10^{46} \, erg s^{-1}$ ) used so far in this kind of analysis. By means of a Principal Component Analysis, we derive Eigenvector (EV) 1 and 2 in an eleven dimensional optical and X-ray parameter space, which are consistent with previous results. The validity of the correlations of the Eddington ratio $\rm L/L_{Edd}$ with EV1 and the black hole mass with EV2 are strongly confirmed These results imply that $\rm L/L_{Edd}$ and black hole mass are related to the diversity of the optical properties of Type 1 AGN. Investigating the relation of the width and asymmetry of $\rm Hβ$ and the relative strength of the iron emission $\rm r_{FeII}$ , we show that our analysis supports the presence of a distinct kinematic region: the Very Broad Line Region. Furthermore, comparing sources with a red-asymmetric broad $\rm Hβ$ emission line to sources for which it is blue-asymmetric, we find an intriguing difference in the correlation of the FeII and the continuum emission strengths. We show that this contrasting behaviour is consistent with a flattened, stratified model of the Broad Line Region, in which the FeII emitting region is shielded from the central source.
△ Less
Submitted 5 November, 2019;
originally announced November 2019.
-
The multiphase gas structure and kinematics in the circumnuclear region of NGC 5728
Authors:
T. Taro Shimizu,
R. I. Davies,
D. Lutz,
L. Burtscher,
M. Lin,
D. Baron,
R. L. Davies,
R. Genzel,
E. K. S. Hicks,
M. Koss,
W. Maciejewski,
F. Mu ller-Sanchez,
G. Orban de Xivry,
S. H. Price,
C. Ricci,
R. Riffel,
R. A. Riffel,
D. Rosario,
M. Schartmann,
A. Schnorr-Muller,
A. Sternberg,
E. Sturm,
T. Storchi-Bergmann,
L. Tacconi,
S. Veilleux
Abstract:
We report on our combined analysis of HST, VLT/MUSE, VLT/SINFONI, and ALMA observations of the local Seyfert 2 galaxy, NGC 5728 to investigate in detail the feeding and feedback of the AGN. The datasets simultaneously probe the morphology, excitation, and kinematics of the stars, ionized gas, and molecular gas over a large range of spatial scales (10 pc--10 kpc). NGC 5728 contains a large stellar…
▽ More
We report on our combined analysis of HST, VLT/MUSE, VLT/SINFONI, and ALMA observations of the local Seyfert 2 galaxy, NGC 5728 to investigate in detail the feeding and feedback of the AGN. The datasets simultaneously probe the morphology, excitation, and kinematics of the stars, ionized gas, and molecular gas over a large range of spatial scales (10 pc--10 kpc). NGC 5728 contains a large stellar bar which is driving gas along prominent dust lanes to the inner 1 kpc where the gas settles into a circumnuclear ring. The ring is strongly star forming and contains a substantial population of young stars as indicated by the lowered stellar velocity dispersion and gas excitation consistent with HII regions. We model the kinematics of the ring using the velocity field of the CO (2--1) emission and stars and find it is consistent with a rotating disk. The outer regions of the disk, where the dust lanes meet the ring, show signatures of inflow at a rate of 1 M$_{\sun}$ yr$^{-1}$. Inside the ring, we observe three molecular gas components corresponding to the circular rotation of the outer ring, a warped disk, and the nuclear stellar bar. The AGN is driving an ionized gas outflow that reaches a radius of 250 pc with a mass outflow rate of 0.08 M$_{\sun}$ yr$^{-1}$ consistent with its luminosity and scaling relations from previous studies. While we observe distinct holes in CO emission which could be signs of molecular gas removal, we find that largely the AGN is not disrupting the structure of the circumnuclear region.
△ Less
Submitted 2 October, 2019; v1 submitted 8 July, 2019;
originally announced July 2019.
-
Machine Learning in Astronomy: a practical overview
Authors:
Dalya Baron
Abstract:
Astronomy is experiencing a rapid growth in data size and complexity. This change fosters the development of data-driven science as a useful companion to the common model-driven data analysis paradigm, where astronomers develop automatic tools to mine datasets and extract novel information from them. In recent years, machine learning algorithms have become increasingly popular among astronomers, a…
▽ More
Astronomy is experiencing a rapid growth in data size and complexity. This change fosters the development of data-driven science as a useful companion to the common model-driven data analysis paradigm, where astronomers develop automatic tools to mine datasets and extract novel information from them. In recent years, machine learning algorithms have become increasingly popular among astronomers, and are now used for a wide variety of tasks. In light of these developments, and the promise and challenges associated with them, the IAC Winter School 2018 focused on big data in Astronomy, with a particular emphasis on machine learning and deep learning techniques. This document summarizes the topics of supervised and unsupervised learning algorithms presented during the school, and provides practical information on the application of such tools to astronomical datasets. In this document I cover basic topics in supervised machine learning, including selection and preprocessing of the input dataset, evaluation methods, and three popular supervised learning algorithms, Support Vector Machines, Random Forests, and shallow Artificial Neural Networks. My main focus is on unsupervised machine learning algorithms, that are used to perform cluster analysis, dimensionality reduction, visualization, and outlier detection. Unsupervised learning algorithms are of particular importance to scientific research, since they can be used to extract new knowledge from existing datasets, and can facilitate new discoveries.
△ Less
Submitted 15 April, 2019;
originally announced April 2019.
-
Discovering AGN-driven winds through their infrared emission: II. Mass outflow rate and energetics
Authors:
Dalya Baron,
Hagai Netzer
Abstract:
The global influence of AGN-driven outflows remains uncertain, due to a lack of large samples with accurately-determined outflow properties. In the second paper of this series, we determine the mass and energetics of ionized outflows is 234 type II AGN, the largest such sample to date, by combining the infrared emission of the dust in the wind (paper I) with the emission line properties. We provid…
▽ More
The global influence of AGN-driven outflows remains uncertain, due to a lack of large samples with accurately-determined outflow properties. In the second paper of this series, we determine the mass and energetics of ionized outflows is 234 type II AGN, the largest such sample to date, by combining the infrared emission of the dust in the wind (paper I) with the emission line properties. We provide new general expressions for the properties of the outflowing gas, which depend on the ionization state of the gas. We also present a novel method to estimate the electron density in the outflow, based on optical line ratios and on the known location of the wind. The inferred electron densities, $n_{\mathrm{e}} \sim 10^{4.5}\,\mathrm{(cm^{-3})}$, are two orders of magnitude larger than typically found in most other cases of ionized outflows. We argue that the discrepancy is due to the fact that the commonly-used [SII]-based method underestimates the true density by a large factor. As a result, the inferred mass outflow rates and kinetic coupling efficiencies are $\dot{M}_{\mathrm{out}} \sim 10^{-2}\, \mathrm{(M_{\odot}/yr)}$ and $ε= \dot{E}_{\mathrm{kin}}/L_{\mathrm{bol}} \sim 10^{-5}$ respectively, 1--2 orders of magnitude lower than previous estimates. Our analysis suggests the existence of a significant amount of neutral atomic gas at the back of the outflowing ionized gas clouds, with mass that is a factor of a few larger than the observed ionized gas mass. This has significant implications for the estimated mass and energetics of such flows.
△ Less
Submitted 26 March, 2019;
originally announced March 2019.
-
Black hole mass estimation for Active Galactic Nuclei from a new angle
Authors:
Dalya Baron,
Brice Ménard
Abstract:
The scaling relations between supermassive black holes and their host galaxy properties are of fundamental importance in the context black hole-host galaxy co-evolution throughout cosmic time. In this work, we use a novel algorithm that identifies smooth trends in complex datasets and apply it to a sample of 2,000 type I active galactic nuclei (AGN) spectra. We detect a sequence in emission line s…
▽ More
The scaling relations between supermassive black holes and their host galaxy properties are of fundamental importance in the context black hole-host galaxy co-evolution throughout cosmic time. In this work, we use a novel algorithm that identifies smooth trends in complex datasets and apply it to a sample of 2,000 type I active galactic nuclei (AGN) spectra. We detect a sequence in emission line shapes and strengths which reveals a correlation between the narrow L([OIII])/L(H$β$) line ratio and the width of the broad H$α$. This scaling relation ties the kinematics of the gas clouds in the broad line region to the ionisation state of the narrow line region, connecting the properties of gas clouds kiloparsecs away from the black hole to material gravitationally bound to it on sub-parsec scales. This relation can be used to estimate black hole masses from narrow emission lines only. It therefore enables black hole mass estimation for obscured type II AGN and allows us to explore the connection between black holes and host galaxy properties for thousands of objects, well beyond the local Universe. Using this technique, we present the $M_{\mathrm{BH}} - σ$ and $M_{\mathrm{BH}} - M_{*}$ scaling relations for a sample of about 10,000 type II AGN from SDSS. These relations are remarkably consistent with those observed for type I AGN, suggesting that this new method may perform as reliably as the classical estimate used in non-obscured type I AGN. These findings open a new window for studies of black hole-host galaxy co-evolution throughout cosmic time.
△ Less
Submitted 5 March, 2019;
originally announced March 2019.
-
Probabilistic Random Forest: A machine learning algorithm for noisy datasets
Authors:
Itamar Reis,
Dalya Baron,
Sahar Shahaf
Abstract:
Machine learning (ML) algorithms become increasingly important in the analysis of astronomical data. However, since most ML algorithms are not designed to take data uncertainties into account, ML based studies are mostly restricted to data with high signal-to-noise ratio. Astronomical datasets of such high-quality are uncommon. In this work we modify the long-established Random Forest (RF) algorit…
▽ More
Machine learning (ML) algorithms become increasingly important in the analysis of astronomical data. However, since most ML algorithms are not designed to take data uncertainties into account, ML based studies are mostly restricted to data with high signal-to-noise ratio. Astronomical datasets of such high-quality are uncommon. In this work we modify the long-established Random Forest (RF) algorithm to take into account uncertainties in the measurements (i.e., features) as well as in the assigned classes (i.e., labels). To do so, the Probabilistic Random Forest (PRF) algorithm treats the features and labels as probability distribution functions, rather than deterministic quantities. We perform a variety of experiments where we inject different types of noise to a dataset, and compare the accuracy of the PRF to that of RF. The PRF outperforms RF in all cases, with a moderate increase in running time. We find an improvement in classification accuracy of up to 10% in the case of noisy features, and up to 30% in the case of noisy labels. The PRF accuracy decreased by less then 5% for a dataset with as many as 45% misclassified objects, compared to a clean dataset. Apart from improving the prediction accuracy in noisy datasets, the PRF naturally copes with missing values in the data, and outperforms RF when applied to a dataset with different noise characteristics in the training and test sets, suggesting that it can be used for Transfer Learning.
△ Less
Submitted 14 November, 2018;
originally announced November 2018.
-
Discovering AGN-driven winds through their infrared emission: I. General method and wind location
Authors:
Dalya Baron,
Hagai Netzer
Abstract:
Large scale outflows of different gas phases are ubiquitous in the host galaxies of active galactic nuclei (AGN). Despite their many differences, they share a common property - they all contain dust. The dust is carried with the outflow, heated by the AGN, and emits at infrared wavelengths. This paper shows that the infrared emission of this dust can be used to detect AGN outflows and derive their…
▽ More
Large scale outflows of different gas phases are ubiquitous in the host galaxies of active galactic nuclei (AGN). Despite their many differences, they share a common property - they all contain dust. The dust is carried with the outflow, heated by the AGN, and emits at infrared wavelengths. This paper shows that the infrared emission of this dust can be used to detect AGN outflows and derive their properties. We use a sample of about 4000 type II AGN and compare the infrared properties of systems that show spectroscopic signature of ionized gas outflows to systems that do not. We detect an additional mid-infrared emission component in galaxies with spectroscopically discovered winds, and attribute it to the dust in the outflow. This new component offers novel constraints on the outflow properties, such as its mean location and covering factor. We measure the location of the outflow for roughly 1700 systems, with the distribution showing a prominent peak around 500 pc, a tail that extends to large distances (10 kpc), and no objects with a location smaller than 50 pc. The covering factor of the wind shows a wide distribution which is centered around 0.1, with 24% (8%) of the winds showing covering factors larger than 0.2 (0.5). The dust emission is not sensitive to various systematics affecting optically-selected outflows, and can be used to estimate the mass outflow rate in thousands of galaxies with only 1D spectra.
△ Less
Submitted 15 October, 2018;
originally announced October 2018.
-
Direct evidence of AGN-feedback: a post starburst galaxy stripped of its gas by AGN-driven winds
Authors:
Dalya Baron,
Hagai Netzer,
J. Xavier Prochaska,
Zheng Cai,
Sebastiano Cantalupo,
D. Christopher Martin,
Mateusz Matuszewski,
Anna M. Moore,
Patrick Morrissey,
James D. Neill
Abstract:
Post starburst E+A galaxies show indications of a powerful starburst that was quenched abruptly. Their disturbed, bulge-dominated morphologies suggest that they are merger remnants. The more massive E+A galaxies are suggested to be quenched by AGN feedback, yet little is known about AGN-driven winds in this short-lived phase. We present spatially-resolved IFU spectroscopy by the Keck Cosmic Web Im…
▽ More
Post starburst E+A galaxies show indications of a powerful starburst that was quenched abruptly. Their disturbed, bulge-dominated morphologies suggest that they are merger remnants. The more massive E+A galaxies are suggested to be quenched by AGN feedback, yet little is known about AGN-driven winds in this short-lived phase. We present spatially-resolved IFU spectroscopy by the Keck Cosmic Web Imager of SDSS J003443.68+251020.9, at z=0.118. The system consists of two galaxies, the larger of which is a post starburst E+A galaxy hosting an AGN. Our modelling suggests a 400 Myrs starburst, with a peak star formation rate of 120 Msun/yr. The observations reveal stationary and outflowing gas, photoionized by the central AGN. We detect gas outflows to a distance of 17 kpc from the central galaxy, far beyond the region of the stars (about 3 kpc), inside a conic structure with an opening angle of 70 degrees. We construct self-consistent photoionization and dynamical models for the different gas components and show that the gas outside the galaxy forms a continuous flow, with a mass outflow rate of about 24 Msun/yr. The gas mass in the flow, roughly $10^9$ Msun, is larger than the total gas mass within the galaxy, some of which is outflowing too. The continuity of the flow puts a lower limit of 60 Myrs on the duration of the AGN feedback. Such AGN are capable of removing, in a single episode, most of the gas from their host galaxies and expelling enriched material into the surrounding CGM.
△ Less
Submitted 9 April, 2018;
originally announced April 2018.
-
On the limitations of statistical absorption studies with the Sloan Digital Sky Surveys I--III
Authors:
Ting-Wen Lan,
Brice Ménard,
Dalya Baron,
Sean Johnson,
Dovi Poznanski,
J. Xavier Prochaska,
John O'Meara
Abstract:
We investigate the limitations of statistical absorption measurements with the SDSS optical spectroscopic surveys. We show that changes in the data reduction strategy throughout different data releases have led to a better accuracy at long wavelengths, in particular for sky line subtraction, but a degradation at short wavelengths with the emergence of systematic spectral features with an amplitude…
▽ More
We investigate the limitations of statistical absorption measurements with the SDSS optical spectroscopic surveys. We show that changes in the data reduction strategy throughout different data releases have led to a better accuracy at long wavelengths, in particular for sky line subtraction, but a degradation at short wavelengths with the emergence of systematic spectral features with an amplitude of about one percent. We show that these features originate from inaccuracy in the fitting of modeled F-star spectra used for flux calibration. The best-fit models for those stars are found to systematically over-estimate the strength of metal lines and under-estimate that of Lithium. We also identify the existence of artifacts due to masking and interpolation procedures at the wavelengths of the hydrogen Balmer series leading to the existence of artificial Balmer $α$ absorption in all SDSS optical spectra. All these effects occur in the rest-frame of the standard stars and therefore present Galactic longitude variations due to the rotation of the Galaxy. We demonstrate that the detection of certain weak absorption lines reported in the literature are solely due to calibration effects. Finally, we discuss new strategies to mitigate these issues.
△ Less
Submitted 4 April, 2018; v1 submitted 6 February, 2018;
originally announced February 2018.
-
Detecting outliers and learning complex structures with large spectroscopic surveys - a case study with APOGEE stars
Authors:
Itamar Reis,
Dovi Poznanski,
Dalya Baron,
Gail Zasowski,
Sahar Shahaf
Abstract:
In this work we apply and expand on a recently introduced outlier detection algorithm that is based on an unsupervised random forest. We use the algorithm to calculate a similarity measure for stellar spectra from the Apache Point Observatory Galactic Evolution Experiment (APOGEE). We show that the similarity measure traces non-trivial physical properties and contains information about complex str…
▽ More
In this work we apply and expand on a recently introduced outlier detection algorithm that is based on an unsupervised random forest. We use the algorithm to calculate a similarity measure for stellar spectra from the Apache Point Observatory Galactic Evolution Experiment (APOGEE). We show that the similarity measure traces non-trivial physical properties and contains information about complex structures in the data. We use it for visualization and clustering of the dataset, and discuss its ability to find groups of highly similar objects, including spectroscopic twins. Using the similarity matrix to search the dataset for objects allows us to find objects that are impossible to find using their best fitting model parameters. This includes extreme objects for which the models fail, and rare objects that are outside the scope of the model. We use the similarity measure to detect outliers in the dataset, and find a number of previously unknown Be-type stars, spectroscopic binaries, carbon rich stars, young stars, and a few that we cannot interpret. Our work further demonstrates the potential for scientific discovery when combining machine learning methods with modern survey data.
△ Less
Submitted 28 May, 2018; v1 submitted 31 October, 2017;
originally announced November 2017.
-
Evidence of ongoing AGN-driven feedback in a quiescent post starburst E+A galaxy
Authors:
Dalya Baron,
Hagai Netzer,
Dovi Poznanski,
J. Xavier Prochaska,
Natascha M. Forster Schreiber
Abstract:
Post starburst E+A galaxies are thought to have experienced a significant starburst that was quenched abruptly. Their disturbed, bulge-dominated morphologies suggest that they are merger remnants. We present ESI/Keck observations of SDSS J132401.63+454620.6, a post starburst galaxy at redshift z = 0.125, with a starburst that started 400 Myr ago, and other properties, like star formation rate (SFR…
▽ More
Post starburst E+A galaxies are thought to have experienced a significant starburst that was quenched abruptly. Their disturbed, bulge-dominated morphologies suggest that they are merger remnants. We present ESI/Keck observations of SDSS J132401.63+454620.6, a post starburst galaxy at redshift z = 0.125, with a starburst that started 400 Myr ago, and other properties, like star formation rate (SFR) consistent with what is measured in ultra luminous infrared galaxies (ULRIGs). The galaxy shows both zero velocity narrow lines, and blueshifted broader Balmer and forbidden emission lines (FWHM=1350 +- 240 km/s). The narrow component is consistent with LINER-like emission, and the broader component with Seyfert-like emission, both photoionized by an active galactic nucleus (AGN) whose properties we measure and model. The velocity dispersion of the broad component exceeds the escape velocity, and we estimate the mass outflow rate to be in the range 4-120 Mo/yr. This is the first reported case of AGN-driven outflows, traced by ionized gas, in post starburst E+A galaxies. We show, by ways of a simple model, that the observed AGN-driven winds can consistently evolve a ULIRG into the observed galaxy. Our findings reinforce the evolutionary scenario where the more massive ULIRGs are quenched by negative AGN feedback, evolve first to post starburst galaxies, and later become typical red and dead ellipticals.
△ Less
Submitted 10 May, 2017;
originally announced May 2017.
-
The Effect of Interstellar Absorption on Measurements of the Baryon Acoustic Peak in the Lyman-α Forest
Authors:
Yishay Vadai,
Dovi Poznanski,
Dalya Baron,
Peter E. Nugent,
David Schlegel
Abstract:
In recent years, the autocorrelation of the hydrogen Lyman-α forest has been used to observe the baryon acoustic peak at redshift 2 < z < 3.5 using tens of thousands of QSO spectra from the BOSS survey. However, the interstellar medium of the Milky-Way introduces absorption lines into the spectrum of any extragalactic source. These lines, while weak and undetectable in a single BOSS spectrum, coul…
▽ More
In recent years, the autocorrelation of the hydrogen Lyman-α forest has been used to observe the baryon acoustic peak at redshift 2 < z < 3.5 using tens of thousands of QSO spectra from the BOSS survey. However, the interstellar medium of the Milky-Way introduces absorption lines into the spectrum of any extragalactic source. These lines, while weak and undetectable in a single BOSS spectrum, could potentially bias the cosmological signal. In order to examine this, we generate absorption line maps by stacking over a million spectra of galaxies and QSOs. We find that the systematics introduced are too small to affect the current accuracy of the baryon acoustic peak, but might be relevant to future surveys such as the Dark Energy Spectroscopic Instrument (DESI). We outline a method to account for this with future datasets.
△ Less
Submitted 17 August, 2017; v1 submitted 9 May, 2017;
originally announced May 2017.
-
The weirdest SDSS galaxies: results from an outlier detection algorithm
Authors:
Dalya Baron,
Dovi Poznanski
Abstract:
How can we discover objects we did not know existed within the large datasets that now abound in astronomy? We present an outlier detection algorithm that we developed, based on an unsupervised Random Forest. We test the algorithm on more than two million galaxy spectra from the Sloan Digital Sky Survey and examine the 400 galaxies with the highest outlier score. We find objects which have extreme…
▽ More
How can we discover objects we did not know existed within the large datasets that now abound in astronomy? We present an outlier detection algorithm that we developed, based on an unsupervised Random Forest. We test the algorithm on more than two million galaxy spectra from the Sloan Digital Sky Survey and examine the 400 galaxies with the highest outlier score. We find objects which have extreme emission line ratios and abnormally strong absorption lines, objects with unusual continua, including extremely reddened galaxies. We find galaxy-galaxy gravitational lenses, double-peaked emission line galaxies, and close galaxy pairs. We find galaxies with high ionisation lines, galaxies which host supernovae, and galaxies with unusual gas kinematics. Only a fraction of the outliers we find were reported by previous studies that used specific and tailored algorithms to find a single class of unusual objects. Our algorithm is general and detects all of these classes, and many more, regardless of what makes them peculiar. It can be executed on imaging, time-series, and other spectroscopic data, operates well with thousands of features, is not sensitive to missing values, and is easily parallelisable.
△ Less
Submitted 22 November, 2016;
originally announced November 2016.
-
Evidence that most type 1 AGN are reddened by dust in the host ISM
Authors:
Dalya Baron,
Jonathan Stern,
Dovi Poznanski,
Hagai Netzer
Abstract:
The typical optical-UV continuum slopes observed in many type 1 AGN are redder than expected from thin accretion disk models. A possible resolution to this conundrum is that many AGN are reddened by dust along the line of sight. To explore this possibility, we stack 5000 SDSS AGN with luminosity L~10^45erg/s and redshift z~0.4 in bins of optical continuum slope alpha_opt and width of the broad H…
▽ More
The typical optical-UV continuum slopes observed in many type 1 AGN are redder than expected from thin accretion disk models. A possible resolution to this conundrum is that many AGN are reddened by dust along the line of sight. To explore this possibility, we stack 5000 SDSS AGN with luminosity L~10^45erg/s and redshift z~0.4 in bins of optical continuum slope alpha_opt and width of the broad H$β$ emission line. We measure the EW of the NaID absorption feature in each stacked spectrum. We find a linear relation between alpha_opt and EW(NaID), such that EW(NaID) increases as alpha_opt becomes redder. In the bin with the smallest H$β$ width, objects with the bluest slopes that are similar to accretion disk predictions are found to have EW(NaID)=0, supporting the line-of-sight dust hypothesis. This conclusion is also supported by the dependence of the $Hα/Hβ$ line ratio on alpha_opt. The implied relationship between alpha_opt and dust reddening is given by E(B-V)~0.2(-0.1-alpha_opt), and the implied reddening of a typical type 1 AGN with alpha_opt=-0.5 is E(B-V)~0.08mag. Photoionization calculations show that the dusty gas responsible for reddening is too ionized to produce the observed features. Therefore, we argue that the sodium absorption arises in regions of the host ISM which are shielded from the AGN radiation, and the correlation with alpha_opt arises since ISM columns along shielded and non-shielded sightlines are correlated. This scenario is supported by the similarity of the E(B-V)-NaID relation implied by our results with the relation in the Milky-Way found by previous studies.
△ Less
Submitted 31 August, 2016; v1 submitted 22 March, 2016;
originally announced March 2016.
-
Using machine learning to classify the diffuse interstellar bands
Authors:
Dalya Baron,
Dovi Poznanski,
Darach Watson,
Yushu Yao,
Nick L. J. Cox,
J. Xavier Prochaska
Abstract:
Using over a million and a half extragalactic spectra we study the correlations of the Diffuse Interstellar Bands (DIBs) in the Milky Way. We measure the correlation between DIB strength and dust extinction for 142 DIBs using 24 stacked spectra in the reddening range E(B-V) < 0.2, many more lines than ever studied before. Most of the DIBs do not correlate with dust extinction. However, we find 10…
▽ More
Using over a million and a half extragalactic spectra we study the correlations of the Diffuse Interstellar Bands (DIBs) in the Milky Way. We measure the correlation between DIB strength and dust extinction for 142 DIBs using 24 stacked spectra in the reddening range E(B-V) < 0.2, many more lines than ever studied before. Most of the DIBs do not correlate with dust extinction. However, we find 10 weak and barely studied DIBs with correlations that are higher than 0.7 with dust extinction and confirm the high correlation of additional 5 strong DIBs. Furthermore, we find a pair of DIBs, 5925.9A and 5927.5A which exhibits significant negative correlation with dust extinction, indicating that their carrier may be depleted on dust. We use Machine Learning algorithms to divide the DIBs to spectroscopic families based on 250 stacked spectra. By removing the dust dependency we study how DIBs follow their local environment. We thus obtain 6 groups of weak DIBs, 4 of which are tightly associated with C2 or CN absorption lines.
△ Less
Submitted 2 May, 2015; v1 submitted 19 January, 2015;
originally announced January 2015.
-
Dusting off the diffuse interstellar bands: DIBs and dust in extragalactic SDSS spectra
Authors:
Dalya Baron,
Dovi Poznanski,
Darach Watson,
Yushu Yao,
J. Xavier Prochaska
Abstract:
Using over a million and a half extragalactic spectra we study the properties of the mysterious Diffuse Interstellar Bands (DIBs) in the Milky Way. These data provide us with an unprecedented sampling of the skies at high Galactic-latitude and low dust-column-density. We present our method, study the correlation of the equivalent width of 8 DIBs with dust extinction and with a few atomic species,…
▽ More
Using over a million and a half extragalactic spectra we study the properties of the mysterious Diffuse Interstellar Bands (DIBs) in the Milky Way. These data provide us with an unprecedented sampling of the skies at high Galactic-latitude and low dust-column-density. We present our method, study the correlation of the equivalent width of 8 DIBs with dust extinction and with a few atomic species, and the distribution of four DIBs - 5780.6A, 5797.1A, 6204.3A, and 6613.6A - over nearly 15000 squared degrees. As previously found, DIBs strengths correlate with extinction and therefore inevitably with each other. However, we show that DIBs can exist even in dust free areas. Furthermore, we find that the DIBs correlation with dust varies significantly over the sky. DIB under- or over-densities, relative to the expectation from dust, are often spread over hundreds of square degrees. These patches are different for the four DIBs, showing that they are unlikely to originate from the same carrier, as previously suggested.
△ Less
Submitted 1 December, 2014; v1 submitted 26 June, 2014;
originally announced June 2014.