-
WISDOM project XX -- Strong shear tearing molecular clouds apart in NGC 524
Authors:
Anan Lu,
Daryl Haggard,
Martin Bureau,
Jindra Gensior,
Sarah Jeffreson,
Carmelle Robert,
Thomas G. Williams,
Fu-Heng Liang,
Woorak Choi,
Timothy A. Davis,
Sara Babic,
Hope Boyce,
Benjamin Cheung,
Laurent Drissen,
Jacob S. Elford,
Lijie Liu,
Thomas Martin,
Carter Rhea,
Laurie Rousseau-Nepton,
Ilaria Ruffa
Abstract:
Early-type galaxies (ETGs) are known to harbour dense spheroids of stars but scarce star formation (SF). Approximately a quarter of these galaxies have rich molecular gas reservoirs yet do not form stars efficiently. We study here the ETG NGC~524, with strong shear suspected to result in a smooth molecular gas disc and low star-formation efficiency (SFE). We present new spatially-resolved observat…
▽ More
Early-type galaxies (ETGs) are known to harbour dense spheroids of stars but scarce star formation (SF). Approximately a quarter of these galaxies have rich molecular gas reservoirs yet do not form stars efficiently. We study here the ETG NGC~524, with strong shear suspected to result in a smooth molecular gas disc and low star-formation efficiency (SFE). We present new spatially-resolved observations of the \textsuperscript{12}CO(2-1)-emitting cold molecular gas from the Atacama Large Millimeter/sub-millimeter Array (ALMA) and of the warm ionised-gas emission lines from SITELLE at the Canada-France-Hawaii Telescope. Although constrained by the resolution of the ALMA observations ($\approx37$~pc), we identify only $52$ GMCs with radii ranging from $30$ to $140$~pc, a low mean molecular gas mass surface density $\langleΣ_{\rm gas}\rangle\approx125$~M$_\odot$~pc$^{-2}$ and a high mean virial parameter $\langleα_{\rm obs,vir}\rangle\approx5.3$. We measure spatially-resolved molecular gas depletion times ($τ_{\rm dep}\equiv1/{\rm SFE}$) with a spatial resolution of $\approx100$~pc within a galactocentric distance of $1.5$~kpc. The global depletion time is $\approx2.0$~Gyr but $τ_{\rm dep}$ increases toward the galaxy centre, with a maximum $τ_{\rm dep,max}\approx5.2$~Gyr. However, no pure \ion{H}{II} region is identified in NGC~524 using ionised-gas emission-line ratio diagnostics, so the $τ_{\rm dep}$ inferred are in fact lower limits. Measuring the GMC properties and dynamical states, we conclude that shear is the dominant mechanism shaping the molecular gas properties and regulating SF in NGC~524. This is supported by analogous analyses of the GMCs in a simulated ETG similar to NGC~524.
△ Less
Submitted 3 June, 2024;
originally announced June 2024.
-
Reconstructing Robust Background IFU spectra using Machine Learning
Authors:
Carter Lee Rhea,
Julie Hlavacek-Larrondo,
Justine Giroux,
Auriane Thilloy,
Hyunseop Choi,
Laurie Rousseau-Nepton,
Marie-Lou Gendron-Marsolais,
Mario Pasquato,
Simon Prunet
Abstract:
In astronomy, spectroscopy consists of observing an astrophysical source and extracting its spectrum of electromagnetic radiation. Once extracted, a model is fit to the spectra to measure the observables, leading to an understanding of the underlying physics of the emission mechanism. One crucial, and often overlooked, aspect of this model is the background emission, which contains foreground and…
▽ More
In astronomy, spectroscopy consists of observing an astrophysical source and extracting its spectrum of electromagnetic radiation. Once extracted, a model is fit to the spectra to measure the observables, leading to an understanding of the underlying physics of the emission mechanism. One crucial, and often overlooked, aspect of this model is the background emission, which contains foreground and background astrophysical sources, intervening atmospheric emission, and artifacts related to the instrument such as noise. This paper proposes an algorithmic approach to constructing a background model for SITELLE observations using statistical tools and supervised machine learning algorithms. We apply a segmentation algorithm implemented in photutils to divide the data cube into background and source spaxels. After applying a principal component analysis (PCA) on the background spaxels, we train an artificial neural network to interpolate from the background to the source spaxels in the PCA coefficient space, which allows us to generate a local background model over the entire data cube. We highlight the performance of this methodology by applying it to SITELLE observations obtained of a SIGNALS galaxy, \NGC4449, and the Perseus galaxy cluster of galaxies, NGC 1275. We discuss the physical interpretation of the principal components and noise reduction in the resulting PCA-based reconstructions. Additionally, we compare the fit results using our new background modeling approach to standard methods used in the literature and find that our method better captures the emission from HII regions in NGC 4449 and the faint emission regions in NGC 1275. These methods also demonstrate that the background does change as a function of the position of the datacube.
△ Less
Submitted 1 April, 2024;
originally announced April 2024.
-
Unraveling the Mysteries of Galaxy Clusters: Recurrent Inference Deconvolution of X-ray Spectra
Authors:
Carter Rhea,
Julie Hlavacek-Larrondo,
Ralph Kraft,
Akos Bogdan,
Alexandre Adam,
Laurence Perreault-Levasseur
Abstract:
In the realm of X-ray spectral analysis, the true nature of spectra has remained elusive, as observed spectra have long been the outcome of convolution between instrumental response functions and intrinsic spectra. In this study, we employ a recurrent neural network framework, the Recurrent Inference Machine (RIM), to achieve the high-precision deconvolution of intrinsic spectra from instrumental…
▽ More
In the realm of X-ray spectral analysis, the true nature of spectra has remained elusive, as observed spectra have long been the outcome of convolution between instrumental response functions and intrinsic spectra. In this study, we employ a recurrent neural network framework, the Recurrent Inference Machine (RIM), to achieve the high-precision deconvolution of intrinsic spectra from instrumental response functions. Our RIM model is meticulously trained on cutting-edge thermodynamic models and authentic response matrices sourced from the Chandra X-ray Observatory archive. Demonstrating remarkable accuracy, our model successfully reconstructs intrinsic spectra well below the 1-sigma error level. We showcase the practical application of this novel approach through real Chandra observations of the galaxy cluster Abell 1550 - a vital calibration target for the recently launched X-ray telescope, XRISM. This work marks a significant stride in the domain of X-ray spectral analysis, offering a promising avenue for unlocking hitherto concealed insights into spectra.
△ Less
Submitted 29 November, 2023;
originally announced November 2023.
-
High-Spectral Resolution Observations of the Optical Filamentary Nebula in NGC 1275
Authors:
Benjamin Vigneron,
Julie Hlavacek-Larrondo,
Carter Lee Rhea,
Marie-Lou Gendron-Marsolais,
Jeremy Lim,
Jake Reinheimer,
Yuan Li,
Laurent Drissen,
Greg L. Bryan,
Megan Donahue,
Alastair Edge,
Andrew Fabian,
Stephen Hamer,
Thomas Martin,
Michael McDonald,
Brian McNamara,
Annabelle Richard-Lafferriere,
Laurie Rousseau-Nepton,
G. Mark Voit,
Tracy Webb,
Norbert Werner
Abstract:
We present new high-spectral resolution observations (R = $λ/Δλ$ = 7000) of the filamentary nebula surrounding NGC 1275, the central galaxy of the Perseus cluster. These observations have been obtained with SITELLE, an imaging Fourier transform spectrometer installed on the Canada-France-Hawai Telescope (CFHT) with a field of view of $11\text{ arcmin }\times 11 \text{ arcmin}$ encapsulating the en…
▽ More
We present new high-spectral resolution observations (R = $λ/Δλ$ = 7000) of the filamentary nebula surrounding NGC 1275, the central galaxy of the Perseus cluster. These observations have been obtained with SITELLE, an imaging Fourier transform spectrometer installed on the Canada-France-Hawai Telescope (CFHT) with a field of view of $11\text{ arcmin }\times 11 \text{ arcmin}$ encapsulating the entire filamentary structure of ionised gas despite its large size of $80 \text{ kpc}\times50 \text{ kpc}$. Here, we present renewed flux, velocity and velocity dispersion maps that show in great detail the kinematics of the optical nebula at \sii$\lambda6716$, \sii$\lambda6731$, \nii$\lambda6584$, H$α$(6563Å), and \nii$\lambda6548$. These maps reveal the existence of a bright flattened disk-shaped structure in the core extending to r $\sim 10$ kpc and dominated by a chaotic velocity field. This structure is located in the wake of X-ray cavities and characterised by a high mean velocity dispersion of $134$ km/s. The disk-shaped structure is surrounded by an extended array of filaments spread out to $r\sim 50$ kpc that are 10 times fainter in flux, remarkably quiescent and has a uniform mean velocity dispersion of $44$ km/s. This stability is puzzling given that the cluster core exhibits several energetic phenomena. Based on these results, we argue that there are two mechanisms to form multiphase gas in clusters of galaxies: a first triggered in the wake of X-ray cavities leading to more turbulent multiphase gas and a second, distinct mechanism, that is gentle and leads to large-scale multiphase gas spread throughout the core.
△ Less
Submitted 27 March, 2024; v1 submitted 27 November, 2023;
originally announced November 2023.
-
Rapid Spectral Parameter Prediction for Black Hole X-Ray Binaries using Physicalised Autoencoders
Authors:
Ethan Tregidga,
James Steiner,
Cecilia Garraffo,
Carter Rhea,
Mayeul Aubin
Abstract:
Black hole X-ray binaries (BHBs) offer insights into extreme gravitational environments and the testing of general relativity. The X-ray spectrum collected by NICER offers valuable information on the properties and behaviour of BHBs through spectral fitting. However, traditional spectral fitting methods are slow and scale poorly with model complexity. This paper presents a new semi-supervised auto…
▽ More
Black hole X-ray binaries (BHBs) offer insights into extreme gravitational environments and the testing of general relativity. The X-ray spectrum collected by NICER offers valuable information on the properties and behaviour of BHBs through spectral fitting. However, traditional spectral fitting methods are slow and scale poorly with model complexity. This paper presents a new semi-supervised autoencoder neural network for parameter prediction and spectral reconstruction of BHBs, showing an improvement of up to a factor of 2,700 in speed while maintaining comparable accuracy. The approach maps the spectral features from the numerous outbursts catalogued by NICER and generalises them to new systems for efficient and accurate spectral fitting. The effectiveness of this approach is demonstrated in the spectral fitting of BHBs and holds promise for use in other areas of astronomy and physics for categorising large datasets.
△ Less
Submitted 2 April, 2024; v1 submitted 26 October, 2023;
originally announced October 2023.
-
A Machine Learning Approach to Galactic Emission-Line Region Classification
Authors:
Carter Lee Rhea,
Laurie Rousseau-Nepton,
Ismael Moumen,
Simon Prunet,
Julie Hlavacek-Larrondo,
Kathryn Grasha,
Carmelle Roberts,
Christophe Morisset,
Grazyna Stasinska,
Natalia Vale-Asari,
Justine Giroux,
Anna McLeod,
Marie-Lou Gendron-Marsolais,
Junfeng Wang,
Joe Lyman,
Laurent Chemin
Abstract:
Diagnostic diagrams of emission-line ratios have been used extensively to categorize extragalactic emission regions; however, these diagnostics are occasionally at odds with each other due to differing definitions. In this work, we study the applicability of supervised machine-learning techniques to systematically classify emission-line regions from the ratios of certain emission lines. Using the…
▽ More
Diagnostic diagrams of emission-line ratios have been used extensively to categorize extragalactic emission regions; however, these diagnostics are occasionally at odds with each other due to differing definitions. In this work, we study the applicability of supervised machine-learning techniques to systematically classify emission-line regions from the ratios of certain emission lines. Using the Million Mexican Model database, which contains information from grids of photoionization models using \texttt{cloudy}, and from shock models, we develop training and test sets of emission line fluxes for three key diagnostic ratios. The sets are created for three classifications: classic \hii{} regions, planetary nebulae, and supernova remnants. We train a neural network to classify a region as one of the three classes defined above given three key line ratios that are present both in the SITELLE and MUSE instruments' band-passes: [{\sc O\,iii}]$\lambda5007$/H$β$, [{\sc N\,ii}]$\lambda6583$/H$α$, ([{\sc S\,ii}]$\lambda6717$+[{\sc S\,ii}]$\lambda6731$)/H$α$. We also tested the impact of the addition of the [{\sc O\,ii}]$\lambda3726,3729$/[{\sc O\,iii}]$\lambda5007$ line ratio when available for the classification. A maximum luminosity limit is introduced to improve the classification of the planetary nebulae. Furthermore, the network is applied to SITELLE observations of a prominent field of M33. We discuss where the network succeeds and why it fails in certain cases. Our results provide a framework for the use of machine learning as a tool for the classification of extragalactic emission regions. Further work is needed to build more comprehensive training sets and adapt the method to additional observational constraints.
△ Less
Submitted 20 June, 2023;
originally announced June 2023.
-
Constraining the LyC escape fraction from LEGUS star clusters with SIGNALS HII region observations: A pilot study of NGC 628
Authors:
J. W. Teh,
K. Grasha,
M. R. Krumholz,
A. Battisti,
D. Calzetti,
L. Rousseau-Nepton,
C. Rhea,
A. Adamo,
R. C. Kennicutt,
E. K. Grebel,
D. O. Cook,
F. Combes,
M. Messa,
S. Linden,
R. S. Klessen,
J. M. Vilchez,
M. Fumagalli,
A. F. McLeod,
L. J. Smith,
L. Chemin,
J. Wang,
E. Sabbi,
E. Sacchi,
A. Petric,
L. Della Bruna
, et al. (1 additional authors not shown)
Abstract:
The ionising radiation of young and massive stars is a crucial form of stellar feedback. Most ionising (Lyman-continuum; LyC, $λ< 912A$) photons are absorbed close to the stars that produce them, forming compact HII regions, but some escape into the wider galaxy. Quantifying the fraction of LyC photons that escape is an open problem. In this work, we present a semi-novel method to estimate the esc…
▽ More
The ionising radiation of young and massive stars is a crucial form of stellar feedback. Most ionising (Lyman-continuum; LyC, $λ< 912A$) photons are absorbed close to the stars that produce them, forming compact HII regions, but some escape into the wider galaxy. Quantifying the fraction of LyC photons that escape is an open problem. In this work, we present a semi-novel method to estimate the escape fraction by combining broadband photometry of star clusters from the Legacy ExtraGalactic UV Survey (LEGUS) with HII regions observed by the Star formation, Ionized gas, and Nebular Abundances Legacy Survey (SIGNALS) in the nearby spiral galaxy NGC 628. We first assess the completeness of the combined catalogue, and find that 49\% of HII regions lack corresponding star clusters as a result of a difference in the sensitivities of the LEGUS and SIGNALS surveys. For HII regions that do have matching clusters, we infer the escape fraction from the difference between the ionising power required to produce the observed HII luminosity and the predicted ionising photon output of their host star clusters; the latter is computed using a combination of LEGUS photometric observations and a stochastic stellar population synthesis code SLUG (Stochastically Lighting Up Galaxies). Overall, we find an escape fraction of $f_{esc} = 0.09^{+0.06}_{-0.06}$ across our sample of 42 HII regions; in particular, we find HII regions with high $f_{esc}$ are predominantly regions with low H$α$-luminosity. We also report possible correlation between $f_{esc}$ and the emission lines [O ii]/[N ii] and [O ii]/H$β$.
△ Less
Submitted 8 June, 2023;
originally announced June 2023.
-
Extended Line Emission in the BCG of Abell 2390
Authors:
Leo Y. Alcorn,
H. K. C Yee,
Laurent Drissen,
Carter Rhea,
Suresh Sivanandam,
Julie Hlavacek-Larrondo,
Bau-Ching Hsieh,
Lihwai Lin,
Yen-Ting Lin,
Qing Liu,
Adam Muzzin,
Allison Noble,
Irene Pintos-Castro
Abstract:
We report CFHT/SITELLE imaging Fourier Transform Spectrograph observations of the Brightest Cluster Galaxy (BCG) of galaxy cluster Abell 2390 at z=0.228. The BCG displays a prominent cone of emission in H$α$, H$β$, [NII], and [OII] to the North-West with PA = 42$^o$, 4.4 arcsec in length (15.9 kpc), which is associated with elongated and asymmetric Chandra soft X-ray emission. The H$α$ flux map al…
▽ More
We report CFHT/SITELLE imaging Fourier Transform Spectrograph observations of the Brightest Cluster Galaxy (BCG) of galaxy cluster Abell 2390 at z=0.228. The BCG displays a prominent cone of emission in H$α$, H$β$, [NII], and [OII] to the North-West with PA = 42$^o$, 4.4 arcsec in length (15.9 kpc), which is associated with elongated and asymmetric Chandra soft X-ray emission. The H$α$ flux map also contains a "hook" of H$α$ and [NII] emission resulting in a broadened northern edge to the cone. Using SITELLE/LUCI software we extract emission line flux, velocity, velocity dispersion, and continuum maps, and utilize them to derive flux ratio maps to determine ionization mechanisms and dynamical information in the BCG's emission line region. The Baldwin-Phillips-Terlevich diagnostics on the BCG cone indicate a composite ionization origin of photoionization due to star formation and shock. Strong LINER-like emission is seen in the nuclear region which hosts an AGN. As Abell 2390 is a cool-core cluster, we suggest that the cooling flow is falling onto the central BCG and interacting with the central AGN. The AGN produces jets that inflate "bubbles" of plasma in the ICM, as is often observed in local galaxy clusters. Furthermore, combining signs of AGN activities from radio, optical emission line and X-ray data over a large range of physical scale, we find evidence for three possible episodes of AGN activity in different epochs associated with the Abell 2390 BCG.
△ Less
Submitted 27 March, 2023;
originally announced March 2023.
-
The First High-Contrast Images of X-Ray Binaries: Detection of Candidate Companions in the $γ$ Cas Analog RX J1744.7$-$2713
Authors:
M. Prasow-Émond,
J. Hlavacek-Larrondo,
K. Fogarty,
J. Rameau,
D. Mawet,
L. -S. Guité,
P. Gandhi,
A. Rao,
J. Steiner,
É. Artigau,
D. Lafrenière,
A. Fabian,
D. Walton,
L. Weiss,
R. Doyon,
C. L. Rhea,
T. Bégin,
B. Vigneron,
M. -E. Naud
Abstract:
X-ray binaries provide exceptional laboratories for understanding the physics of matter under the most extreme conditions. Until recently, there were few, if any, observational constraints on the circumbinary environments of X-ray binaries at $\sim$ 100-5000 AU scales; it remains unclear how the accretion onto the compact objects or the explosions giving rise to the compact objects interact with t…
▽ More
X-ray binaries provide exceptional laboratories for understanding the physics of matter under the most extreme conditions. Until recently, there were few, if any, observational constraints on the circumbinary environments of X-ray binaries at $\sim$ 100-5000 AU scales; it remains unclear how the accretion onto the compact objects or the explosions giving rise to the compact objects interact with their immediate surroundings. Here, we present the first high-contrast adaptive optics images of X-ray binaries. These observations target all X-ray binaries within $\sim$ 3 kpc accessible with the Keck/NIRC2 vortex coronagraph. This paper focuses on one of the first key results from this campaign: our images reveal the presence of 21 sources potentially associated with the $γ$ Cassiopeiae analog high-mass X-ray binary RX J1744.7$-$2713. By conducting different analyses - a preliminary proper motion analysis, a color-magnitude diagram and a probability of chance alignment calculation - we found that three of these 21 sources have a high probability of being bound to the system. If confirmed, they would be in wide orbits ($\sim$ 450 AU to 2500 AU). While follow-up astrometric observations will be needed in $\sim$ 5-10 years to confirm further the bound nature of these detections, these discoveries emphasize that such observations may provide a major breakthrough in the field. In fact, they would be useful not only for our understanding of stellar multiplicity but also for our understanding of how planets, brown dwarfs and stars can form even in the most extreme environments.
△ Less
Submitted 10 May, 2022;
originally announced May 2022.
-
An Assessment of the In-Situ Growth of the Intracluster Light in the High Redshift Galaxy Cluster SpARCS1049+56
Authors:
Capucine Barfety,
Félix-Antoine Valin,
Tracy M. A. Webb,
Min Yun,
Heath Shipley,
Kyle Boone,
Brian Hayden,
Julie Hlavacek-Larrondo,
Adam Muzzin,
Allison G. Noble,
Saul Perlmutter,
Carter Rhea,
Gillian Wilson,
H. K. C Yee
Abstract:
The formation of the stellar mass within galaxy cluster cores is a poorly understood process. It features the complicated physics of cooling flows, AGN feedback, star formation and more. Here, we study the growth of the stellar mass in the vicinity of the Brightest Cluster Galaxy (BCG) in a z = 1.7 cluster, SpARCS1049+56. We synthesize a reanalysis of existing HST imaging, a previously published m…
▽ More
The formation of the stellar mass within galaxy cluster cores is a poorly understood process. It features the complicated physics of cooling flows, AGN feedback, star formation and more. Here, we study the growth of the stellar mass in the vicinity of the Brightest Cluster Galaxy (BCG) in a z = 1.7 cluster, SpARCS1049+56. We synthesize a reanalysis of existing HST imaging, a previously published measurement of the star formation rate, and the results of new radio molecular gas spectroscopy. These analyses represent the past, present and future star formation respectively within this system. We show that a large amount of stellar mass -- between $(2.2 \pm 0.5) \times 10^{10} \: M_\odot$ and $(6.6 \pm 1.2) \times 10^{10}\: M_\odot$ depending on the data processing -- exists in a long and clumpy tail-like structure that lies roughly 12 kpc off the BCG. Spatially coincident with this stellar mass is a similarly massive reservoir ($(1.0 \pm 0.7) \times 10^{11} \: M_\odot$) of molecular gas that we suggest is the fuel for the immense star formation rate of $860 \pm 130 \: M_\odot$/yr, as measured by infrared observations. Hlavacek-Larrondo et al. 2021 surmised that massive, runaway cooling of the hot intracluster X-ray gas was feeding this star formation, a process that had not been observed before at high-redshift. We conclude, based on the amount of fuel and current stars, that this event may be rare in the lifetime of a cluster, producing roughly 15 to 21% of the Intracluster Light (ICL) mass in one go, though perhaps a common event for all galaxy clusters.
△ Less
Submitted 25 March, 2022;
originally announced March 2022.
-
Extended radio emission in the galaxy cluster MS 0735.6+7421 detected with the Karl G. Jansky Very Large Array
Authors:
T. Bégin,
J. Hlavacek-Larrondo,
C. L. Rhea,
M. Gendron-Marsolais,
B. McNamara,
R. J. van Weeren,
A. Richard-Laferrière,
L. Guité,
M. Prasow-Émond,
D. Haggard
Abstract:
MS 0735.6+7421 ($z = 0.216$) is a massive cool core galaxy cluster hosting one of the most powerful active galactic nuclei (AGN) outbursts known. The radio jets of the AGN have carved out an unusually large pair of X-ray cavities, each reaching a diameter of $200$ kpc. This makes MS 0735.6+7421 a unique case to investigate active galactic nuclei feedback processes, as well as other cluster astroph…
▽ More
MS 0735.6+7421 ($z = 0.216$) is a massive cool core galaxy cluster hosting one of the most powerful active galactic nuclei (AGN) outbursts known. The radio jets of the AGN have carved out an unusually large pair of X-ray cavities, each reaching a diameter of $200$ kpc. This makes MS 0735.6+7421 a unique case to investigate active galactic nuclei feedback processes, as well as other cluster astrophysics at radio wavelengths. We present new low-radio-frequency observations of MS 0735.6+7421 taken with the Karl G. Jansky Very Large Array (VLA): 5 hours of P-band ($224-480$ MHz) and 5 hours of L-band ($1-2$ GHz) observations, both in C configuration. Our VLA P-band ($224-480$ MHz) observations reveal the presence of a new diffuse radio component reaching a scale of $\sim$ $900$ kpc in the direction of the jets and of $\sim$ $500$ kpc in the direction perpendicular to the jets. This component is centered on the cluster core and has a radio power scaled at $1.4$ GHz of $P_{1.4\text{ GHz}} = (4\pm2)\times 10^{24}$ WHz$^{-1}$. Its properties are consistent with those expected from a radio mini-halo as seen in other massive cool core clusters, although it may also be associated with radio plasma that has diffused out of the X-ray cavities. Observations at higher spatial resolution are needed to fully characterize the properties and nature of this component. We also suggest that if radio mini-halos originate from jetted activity, we may be witnessing the early stages of this process.
△ Less
Submitted 19 January, 2023; v1 submitted 2 February, 2022;
originally announced February 2022.
-
Updates to LUCI: A New Fitting Paradigm Using Mixture Density Networks
Authors:
Carter L. Rhea,
Julie Hlavacek-Larrondo,
Laurie Rousseau-Nepton,
Simon Prunet
Abstract:
LUCI is an general-purpose spectral line-fitting pipeline which natively integrates machine learning algorithms to initialize fit functions. LUCI currently uses point-estimates obtained from a convolutional neural network (CNN) to inform optimization algorithms; this methodology has shown great promise by reducing computation time and reducing the chance of falling into a local minimum using conve…
▽ More
LUCI is an general-purpose spectral line-fitting pipeline which natively integrates machine learning algorithms to initialize fit functions. LUCI currently uses point-estimates obtained from a convolutional neural network (CNN) to inform optimization algorithms; this methodology has shown great promise by reducing computation time and reducing the chance of falling into a local minimum using convex optimization methods. In this update to LUCI, we expand upon the CNN developed in Rhea et al. 2020 so that it outputs Gaussian posterior distributions of the fit parameters of interest (the velocity and broadening) rather than simple point-estimates. Moreover, these posteriors are then used to inform the priors in a Bayesian inference scheme, either emcee or dynesty. The code is publicly available at https://github.com/crhea93/LUCI.
△ Less
Submitted 24 November, 2021;
originally announced November 2021.
-
A Machine Learning Approach to Integral Field Unit Spectroscopy Observations: III. Disentangling Multiple Components in Hii regions
Authors:
Carter Lee Rhea,
Laurie Rousseau-Nepton,
Simon Prunet,
Julie Hlavacek-Larrondo,
R. Pierre Martin,
Kathryn Grasha,
Natalia Vale Asari,
Théophile Bégin,
Benjamin Vigneron,
Myriam Prasow-Émond
Abstract:
In the first two papers of this series (Rhea et al. 2020; Rhea et al. 2021), we demonstrated the dynamism of machine learning applied to optical spectral analysis by using neural networks to extract kinematic parameters and emission-line ratios directly from the spectra observed by the SITELLE instrument located at the Canada-France-Hawai'i Telescope. In this third installment, we develop a framew…
▽ More
In the first two papers of this series (Rhea et al. 2020; Rhea et al. 2021), we demonstrated the dynamism of machine learning applied to optical spectral analysis by using neural networks to extract kinematic parameters and emission-line ratios directly from the spectra observed by the SITELLE instrument located at the Canada-France-Hawai'i Telescope. In this third installment, we develop a framework using a convolutional neural network trained on synthetic spectra to determine the number of line-of-sight components present in the SN3 filter (656--683nm) spectral range of SITELLE. We compare this methodology to standard practice using Bayesian Inference. Our results demonstrate that a neural network approach returns more accurate results and uses less computational resources over a range of spectral resolutions. Furthermore, we apply the network to SITELLE observations of the merging galaxy system NGC2207/IC2163. We find that the closest interacting sector and the central regions of the galaxies are best characterized by two line-of-sight components while the outskirts and spiral arms are well-constrained by a single component. Determining the number of resolvable components is crucial in disentangling different galactic components in merging systems and properly extracting their respective kinematics.
△ Less
Submitted 1 October, 2021;
originally announced October 2021.
-
LUCI: A Python package for SITELLE spectral analysis
Authors:
Carter Lee Rhea,
Julie Hlavacek-Larrondo,
Laurie Rousseau-Nepton,
Benjamin Vigneron,
Louis-Simon Guité
Abstract:
High-resolution optical integral field units (IFUs) are rapidly expanding our knowledge of extragalactic emission nebulae in galaxies and galaxy clusters. By studying the spectra of these objects -- which include classic HII regions, supernova remnants, planetary nebulae, and cluster filaments -- we are able to constrain their kinematics (velocity and velocity dispersion). In conjunction with addi…
▽ More
High-resolution optical integral field units (IFUs) are rapidly expanding our knowledge of extragalactic emission nebulae in galaxies and galaxy clusters. By studying the spectra of these objects -- which include classic HII regions, supernova remnants, planetary nebulae, and cluster filaments -- we are able to constrain their kinematics (velocity and velocity dispersion). In conjunction with additional tools, such as the BPT diagram, we can further classify emission regions based on strong emission-line flux ratios. LUCI is a simple-to-use python module intended to facilitate the rapid analysis of IFU spectra. LUCI does this by integrating well-developed pre-existing python tools such as astropy and scipy with new machine learning tools for spectral analysis (Rhea et al. 2020). Furthermore, LUCI provides several easy-to-use tools to access and fit SITELLE data cubes.
△ Less
Submitted 27 August, 2021;
originally announced August 2021.
-
A Data-driven Approach to X-ray Spectral Fitting: Quasi-Deconvolution
Authors:
Carter Lee Rhea,
Julie Hlavacek-Larrondo,
Ralph Kraft,
Akos Bogdan,
Rudy Geelen
Abstract:
X-ray spectral fitting of astronomical sources requires convolving the intrinsic spectrum or model with the instrumental response. Standard forward modeling techniques have proven success in recovering the underlying physical parameters in moderate to high signal-to-noise regimes; however, they struggle to achieve the same level of accuracy in low signal-to-noise regimes. Additionally, the use of…
▽ More
X-ray spectral fitting of astronomical sources requires convolving the intrinsic spectrum or model with the instrumental response. Standard forward modeling techniques have proven success in recovering the underlying physical parameters in moderate to high signal-to-noise regimes; however, they struggle to achieve the same level of accuracy in low signal-to-noise regimes. Additionally, the use of machine learning techniques on X-ray spectra requires access to the intrinsic spectrum. Therefore, the measured spectrum must be effectively deconvolved from the instrumental response. In this note, we explore numerical methods for inverting the matrix equation describing X-ray spectral convolution. We demonstrate that traditional methods are insufficient to recover the intrinsic X-ray spectrum and argue that a novel approach is required.
△ Less
Submitted 13 May, 2021;
originally announced May 2021.
-
A Machine Learning Approach to Integral Field Unit Spectroscopy Observations: II. HII Region LineRatios
Authors:
Carter Rhea,
Laurie Rousseau-Nepton,
Simon Prunet,
Myriam Prasow-Emond,
Julie Hlavacek-Larrondo,
Natalia Vale Asari,
Kathryn Grasha,
Laurence Perreault-Levasseur
Abstract:
In the first paper of this series (Rhea et al. 2020), we demonstrated that neural networks can robustly and efficiently estimate kinematic parameters for optical emission-line spectra taken by SITELLE at the Canada-France-Hawaii Telescope. This paper expands upon this notion by developing an artificial neural network to estimate the line ratios of strong emission-lines present in the SN1, SN2, and…
▽ More
In the first paper of this series (Rhea et al. 2020), we demonstrated that neural networks can robustly and efficiently estimate kinematic parameters for optical emission-line spectra taken by SITELLE at the Canada-France-Hawaii Telescope. This paper expands upon this notion by developing an artificial neural network to estimate the line ratios of strong emission-lines present in the SN1, SN2, and SN3 filters of SITELLE. We construct a set of 50,000 synthetic spectra using line ratios taken from the Mexican Million Model database replicating Hii regions. Residual analysis of the network on the test set reveals the network's ability to apply tight constraints to the line ratios. We verified the network's efficacy by constructing an activation map, checking the [N ii] doublet fixed ratio, and applying a standard k-fold cross-correlation. Additionally, we apply the network to SITELLE observation of M33; the residuals between the algorithm's estimates and values calculated using standard fitting methods show general agreement. Moreover, the neural network reduces the computational costs by two orders of magnitude. Although standard fitting routines do consistently well depending on the signal-to-noise ratio of the spectral features, the neural network can also excel at predictions in the low signal-to-noise regime within the controlled environment of the training set as well as on observed data when the source spectral properties are well constrained by models. These results reinforce the power of machine learning in spectral analysis.
△ Less
Submitted 11 February, 2021;
originally announced February 2021.
-
The H$α$ star formation main sequence in cluster and field galaxies at $z\sim1.6$
Authors:
Julie Nantais,
Gillian Wilson,
Adam Muzzin,
Lyndsay J. Old,
Ricardo Demarco,
Pierluigi Cerulo,
Michael Balogh,
Gregory Rudnick,
Jeffrey Chan,
M. C. Cooper,
Ben Forrest,
Brian Hayden,
Chris Lidman,
Allison Noble,
Saul Perlmutter,
Carter Rhea,
Jason Surace,
Remco van der Burg,
Eelco van Kampen
Abstract:
We calculate H$α$-based star formation rates and determine the star formation rate-stellar mass relation for members of three SpARCS clusters at $z \sim 1.6$ and serendipitously identified field galaxies at similar redshifts to the clusters. We find similar star formation rates in cluster and field galaxies throughout our range of stellar masses. The results are comparable to those seen in other c…
▽ More
We calculate H$α$-based star formation rates and determine the star formation rate-stellar mass relation for members of three SpARCS clusters at $z \sim 1.6$ and serendipitously identified field galaxies at similar redshifts to the clusters. We find similar star formation rates in cluster and field galaxies throughout our range of stellar masses. The results are comparable to those seen in other clusters at similar redshifts, and consistent with our previous photometric evidence for little quenching activity in clusters. One possible explanation for our results is that galaxies in our $z \sim 1.6$ clusters have been accreted too recently to show signs of environmental quenching. It is also possible that the clusters are not yet dynamically mature enough to produce important environmental quenching effects shown to be important at low redshift, such as ram pressure stripping or harassment.
△ Less
Submitted 16 September, 2020;
originally announced September 2020.
-
A Novel Machine Learning Approach to Disentangle Multi-Temperature Regions in Galaxy Clusters
Authors:
Carter L. Rhea,
Julie Hlavacek-Larrondo,
Laurence Perreault-Levasseur,
Marie-Lou Gendron-Marsolais,
Ralph Kraft
Abstract:
The hot intra-cluster medium (ICM) surrounding the heart of galaxy clusters is a complex medium comprised of various emitting components. Although previous studies of nearby galaxy clusters, such as the Perseus, the Coma, or the Virgo cluster, have demonstrated the need for multiple thermal components when spectroscopically fitting the ICM's X-ray emission, no systematic methodology for calculatin…
▽ More
The hot intra-cluster medium (ICM) surrounding the heart of galaxy clusters is a complex medium comprised of various emitting components. Although previous studies of nearby galaxy clusters, such as the Perseus, the Coma, or the Virgo cluster, have demonstrated the need for multiple thermal components when spectroscopically fitting the ICM's X-ray emission, no systematic methodology for calculating the number of underlying components currently exists. In turn, underestimating or overestimating the number of components can cause systematic errors in the emission parameter estimations. In this paper, we present a novel approach to determining the number of components using an amalgam of machine learning techniques. Synthetic spectra containing a various number of underlying thermal components were created using well-established tools available from the \textit{Chandra} X-ray Observatory. The dimensions of the training set was initially reduced using the Principal Component Analysis and then categorized based on the number of underlying components using a Random Forest Classifier. Our trained and tested algorithm was subsequently applied to \textit{Chandra} X-ray observations of the Perseus cluster. Our results demonstrate that machine learning techniques can efficiently and reliably estimate the number of underlying thermal components in the spectra of galaxy clusters, regardless of the thermal model (MEKAL versus APEC). %and signal-to-noise ratio used. We also confirm that the core of the Perseus cluster contains a mix of differing underlying thermal components. We emphasize that although this methodology was trained and applied on \textit{Chandra} X-ray observations, it is readily portable to other current (e.g. XMM-Newton, eROSITA) and upcoming (e.g. Athena, Lynx, XRISM) X-ray telescopes. The code is publicly available at \url{https://github.com/XtraAstronomy/Pumpkin}.
△ Less
Submitted 1 September, 2020;
originally announced September 2020.
-
Machine Learning Approach to Integral Field Unit Spectroscopy Observations: I. HII Region Kinematics
Authors:
Carter L. Rhea,
Laurie Rousseau-Nepton,
Simon Prunet,
Julie Hlavacek-larrondo,
Sebastien Fabbro
Abstract:
SITELLE is a novel integral field unit spectroscopy instrument that has an impressive spatial (11 by 11 arcmin), spectral coverage, and spectral resolution (R=1-20000). SIGNALS is anticipated to obtain deep observations (down to 3.6x10-17ergs s-1cm-2) of 40 galaxies, each needing complex and substantial time to extract spectral information. We present a method that uses Convolution Neural Networks…
▽ More
SITELLE is a novel integral field unit spectroscopy instrument that has an impressive spatial (11 by 11 arcmin), spectral coverage, and spectral resolution (R=1-20000). SIGNALS is anticipated to obtain deep observations (down to 3.6x10-17ergs s-1cm-2) of 40 galaxies, each needing complex and substantial time to extract spectral information. We present a method that uses Convolution Neural Networks (CNN) for estimating emission line parameters in optical spectra obtained with SITELLE as part of the SIGNALS large program. Our algorithm is trained and tested on synthetic data representing typical emission spectra for HII regions based on Mexican Million Models database(3MdB) BOND simulations. The network's activation map demonstrates its ability to extract the dynamical (broadening and velocity) parameters from a set of 5 emission lines (e.g. Hα, N[II] doublet, and S[II] doublet) in the SN3 (651-685 nm) filter of SITELLE. Once trained, the algorithm was tested on real SITELLE observations in the SIGNALS program of one of the South West fields of M33. The CNN recovers the dynamical parameters with an accuracy better than 5 km s-1 in regions with a signal-to-noise ratio greater than 15 over the Hαline. More importantly, our CNN method reduces calculation time by over an order of magnitude on the spectral cube with native spatial resolution when compared with standard fitting procedures. These results clearly illustrate the power of machine learning algorithms for the use in future IFU-based missions. Subsequent work will explore the applicability of the methodology to other spectral parameters such as the flux of key emission lines.
△ Less
Submitted 25 August, 2020; v1 submitted 18 August, 2020;
originally announced August 2020.
-
Evidence of runaway gas cooling in the absence of supermassive black hole feedback at the epoch of cluster formation
Authors:
J. Hlavacek-Larrondo,
C. L. Rhea,
T. Webb,
M. McDonald,
A. Muzzin,
G. Wilson,
K. Finner,
F. Valin,
N. Bonaventura,
M. Cooper,
A. C. Fabian,
M. -L. Gendron-Marsolais,
M. J. Jee,
C. Lidman,
M. Mezcua,
A. Noble,
H. R. Russell,
J. Surace,
A. Trudeau,
H. K. C. Yee
Abstract:
Cosmological simulations, as well as mounting evidence from observations, have shown that supermassive black holes play a fundamental role in regulating the formation of stars throughout cosmic time. This has been clearly demonstrated in the case of galaxy clusters in which powerful feedback from the central black hole is preventing the hot intracluster gas from cooling catastrophically, thus redu…
▽ More
Cosmological simulations, as well as mounting evidence from observations, have shown that supermassive black holes play a fundamental role in regulating the formation of stars throughout cosmic time. This has been clearly demonstrated in the case of galaxy clusters in which powerful feedback from the central black hole is preventing the hot intracluster gas from cooling catastrophically, thus reducing the expected star formation rates by orders of magnitude. These conclusions have however been almost entirely based on nearby clusters. Based on new Chandra X-ray observations, we present the first observational evidence for massive, runaway cooling occurring in the absence of supermassive black hole feedback in the high-redshift galaxy cluster SpARCS104922.6+564032.5 ($z=1.709$). The hot intracluster gas appears to be fueling a massive burst of star formation ($\approx900$~M$_\odot$yr$^{-1}$) that is offset by dozens of kpc from the central galaxy. The burst is co-spatial with the coolest intracluster gas but not associated with any galaxy in the cluster. In less than 100 million years, such runaway cooling can form the same amount of stars as in the Milky Way. Intracluster stars are therefore not only produced by tidal stripping and the disruption of cluster galaxies, but can also be produced by runaway cooling of hot intracluster gas at early times. Overall, these observations show the dramatic impact when supermassive black hole feedback fails to operate in clusters. They indicate that in the highest overdensities such as clusters and proto-clusters, runaway cooling may be a new and important mechanism for fueling massive bursts of star formation in the early universe.
△ Less
Submitted 30 July, 2020;
originally announced July 2020.
-
On the relation between mini-halos and AGN feedback in clusters of galaxies
Authors:
A. Richard-Laferrière,
J. Hlavacek-Larrondo,
R. S. Nemmen,
C. L. Rhea,
G. B. Taylor,
M. Prasow-Émond,
M. Gendron-Marsolais,
M. Latulippe,
A. C. Edge,
A. C. Fabian,
J. S. Sanders,
M. T. Hogan,
G. Demontigny
Abstract:
A variety of large-scale diffuse radio structures have been identified in many clusters with the advent of new state-of-the-art facilities in radio astronomy. Among these diffuse radio structures, radio mini-halos are found in the central regions of cool core clusters. Their origin is still unknown and they are challenging to discover; less than thirty have been published to date. Based on new VLA…
▽ More
A variety of large-scale diffuse radio structures have been identified in many clusters with the advent of new state-of-the-art facilities in radio astronomy. Among these diffuse radio structures, radio mini-halos are found in the central regions of cool core clusters. Their origin is still unknown and they are challenging to discover; less than thirty have been published to date. Based on new VLA observations, we confirmed the mini-halo in the massive strong cool core cluster PKS 0745$-$191 ($z=0.1028$) and discovered one in the massive cool core cluster MACS J1447.4+0827 ($z=0.3755$). Furthermore, using a detailed analysis of all known mini-halos, we explore the relation between mini-halos and AGN feedback processes from the central galaxy. We find evidence of strong, previously unknown correlations between mini-halo radio power and X-ray cavity power, and between mini-halo and the central galaxy radio power related to the relativistic jets when spectrally decomposing the AGN radio emission into a component for past outbursts and one for on-going accretion. Overall, our study indicates that mini-halos are directly connected to the central AGN in clusters, following previous suppositions. We hypothesize that AGN feedback may be one of the dominant mechanisms giving rise to mini-halos by injecting energy into the intra-cluster medium and reaccelerating an old population of particles, while sloshing motion may drive the overall shape of mini-halos inside cold fronts. AGN feedback may therefore not only play a vital role in offsetting cooling in cool core clusters, but may also play a fundamental role in re-energizing non-thermal particles in clusters.
△ Less
Submitted 28 September, 2020; v1 submitted 2 July, 2020;
originally announced July 2020.
-
A Multiwavelength Study of the Massive Cool Core Cluster MACS J1447.4+0827
Authors:
M. Prasow-Émond,
J. Hlavacek-Larrondo,
C. L. Rhea,
M. Latulippe,
M. -L. Gendron-Marsolais,
A. Richard-Laferrière,
J. S. Sanders,
A. C. Edge,
S. W. Allen,
A. Mantz,
A. von der Linden
Abstract:
Clusters of galaxies are outstanding laboratories for understanding the physics of supermassive black hole feedback. Here, we present the first \textit{Chandra}, Karl G. Janksy Very Large Array and \textit{Hubble Space Telescope} analysis of MACS J1447.4+0827 ($z = 0.3755$), one of the strongest cool core clusters known, in which extreme feedback from its central supermassive black hole is needed…
▽ More
Clusters of galaxies are outstanding laboratories for understanding the physics of supermassive black hole feedback. Here, we present the first \textit{Chandra}, Karl G. Janksy Very Large Array and \textit{Hubble Space Telescope} analysis of MACS J1447.4+0827 ($z = 0.3755$), one of the strongest cool core clusters known, in which extreme feedback from its central supermassive black hole is needed to prevent the hot intracluster gas from cooling. Using this multiwavelength approach, including 70 ks of \textit{Chandra} X-ray observations, we detect the presence of collimated jetted-outflows that coincides with a southern and a northern X-ray cavity. The total mechanical power associated with these outflows ($P_{\mathrm{cav}} \approx 6 \times 10^{44}$ erg s$^{-1}$) is roughly consistent with the energy required to prevent catastrophic cooling of the hot intracluster gas ($L_{\mathrm{cool}} = 1.71 \pm 0.01 \times 10^{45}$ erg s$^{-1}$ for t$_\mathrm{cool}$ = 7.7 Gyrs); implying that powerful supermassive black hole feedback has been in place several Giga-years ago in MACS J1447.7+0827. In addition, we detect the presence of a radio mini-halo that extends over 300 kpc in diameter ($P_{1.4 \mathrm{GHz}} = 3.0 \pm 0.3 \times 10^{24}$ W Hz$^{-1}$). The X-ray observations also reveal a $\sim20$ kpc plume-like structure that coincides with optical dusty filaments that surround the central galaxy. Overall, this study demonstrates that the various physical phenomena occurring in the most nearby clusters of galaxies are also occurring in their more distant analogues.
△ Less
Submitted 8 June, 2020;
originally announced June 2020.
-
Gravitational Lensing Size Scales for Quasars
Authors:
G. Chartas,
C. Rhea,
C. Kochanek,
X. Dai,
C. Morgan,
J. Blackburne,
B. Chen,
A. Mosquera,
C. MacLeod
Abstract:
We review results from our monitoring observations of several lensed quasars performed in the optical, UV, and X-ray bands. Modeling of the multi-wavelength light curves provides constraints on the extent of the optical, UV, and X-ray emission regions. One of the important results of our analysis is that the optical sizes as inferred from the microlensing analysis are significantly larger than tho…
▽ More
We review results from our monitoring observations of several lensed quasars performed in the optical, UV, and X-ray bands. Modeling of the multi-wavelength light curves provides constraints on the extent of the optical, UV, and X-ray emission regions. One of the important results of our analysis is that the optical sizes as inferred from the microlensing analysis are significantly larger than those predicted by the theoretical-thin-disk estimate. In a few cases we also constrain the slope of the size-wavelength relation. Our size constraints of the soft and hard X-ray emission regions of quasars indicate that in some objects of our sample the hard X-ray emission region is more compact than the soft and in others the soft emission region is smaller. This difference may be the result of the relative strengths of the disk-reflected (harder and extended) versus corona-direct (softer and compact) components in the quasars of our sample. Finally, we present the analysis of several strong microlensing events where we detect an evolution of the relativistic Fe line profile as the magnification caustic traverses the accretion disk. These caustic crossings are used to provide constraints on the innermost stable circular orbit (ISCO) radius and the accretion disk inclination angle of the black hole in quasar RX J1131-1231.
△ Less
Submitted 17 September, 2015;
originally announced September 2015.