-
Polycyclic Aromatic Hydrocarbon Emission in Galaxies as seen with JWST
Authors:
Dimitra Rigopoulou,
Fergus R. Donnan,
Ismael García-Bernete,
Miguel Pereira-Santaella,
Almudena Alonso-Herrero,
Ric Davies,
Leslie K. Hunt,
Patrick F. Roche,
Taro Shimizu
Abstract:
We present a systematic study of mid-infrared spectra of galaxies including star-forming galaxies and Active Galactic Nuclei observed with JWST MIRI-MRS and NIRSpec-IFU. We focus on the relative variations of the 3.3, 6.2, 7.7, 11.3, 12.7 and 17 micron Polycyclic Aromatic Hydrocarbon (PAH) features within spatially resolved regions of galaxies including NGC 3256, NGC 7469, VV 114, II Zw96 and NGC…
▽ More
We present a systematic study of mid-infrared spectra of galaxies including star-forming galaxies and Active Galactic Nuclei observed with JWST MIRI-MRS and NIRSpec-IFU. We focus on the relative variations of the 3.3, 6.2, 7.7, 11.3, 12.7 and 17 micron Polycyclic Aromatic Hydrocarbon (PAH) features within spatially resolved regions of galaxies including NGC 3256, NGC 7469, VV 114, II Zw96 and NGC 5728. Using theoretical PAH models and extending our earlier work, we introduce a new PAH diagnostic involving the 17 micron PAH feature. To determine the drivers of PAH band variations in galaxies, we compare observed PAH spectral bands to predictions from theoretical PAH models. We consider extinction, dehydrogenation and PAH size and charge as possible drivers of PAH band variations. We find a surprising uniformity in PAH size distribution among the spatially resolved regions of the galaxies studied here, with no evidence for preferential destruction of the smallest grains, contrary to earlier findings. Neither extinction nor dehydrogenation play a crucial role in setting the observed PAH bands. Instead, we find that PAH charge plays a significant role in PAH inter-band variations. We find a tight relation between PAH charge and the intensity of the radiation field as traced by the [NeIII]/[NeII] maps. In agreement with recent JWST results, we find a predominance of neutral PAH molecules in the nuclei of Active Galaxies and their outflows. Ionised PAHs are the dominant population in star-forming galaxies. We discuss the implications of our findings for the use of PAHs as ISM tracers in high redshift galaxies.
△ Less
Submitted 17 June, 2024;
originally announced June 2024.
-
The Kinematics of Polycyclic Aromatic Hydrocarbons (PAHs) in Galaxies revealed by Principal Component Analysis (PCA) tomography with JWST/NIRSpec
Authors:
Fergus R. Donnan,
Dimitra Rigopoulou,
Ismael García-Bernete
Abstract:
Polycyclic Aromatic Hydrocarbons (PAHs) are organic molecules which comprise the smallest particles of dust in the interstellar medium (ISM). Due to their broad/complex emission profiles, obtaining kinematics is a challenge with traditional methods, especially before the advent of the JWST. In this work we employ Principal Component Analysis (PCA) tomography to analyse JWST/NIRSpec IFU data of 3 n…
▽ More
Polycyclic Aromatic Hydrocarbons (PAHs) are organic molecules which comprise the smallest particles of dust in the interstellar medium (ISM). Due to their broad/complex emission profiles, obtaining kinematics is a challenge with traditional methods, especially before the advent of the JWST. In this work we employ Principal Component Analysis (PCA) tomography to analyse JWST/NIRSpec IFU data of 3 nearby luminous infrared galaxies (LIRGs), namely, NGC 3256 N, NGC 3256 S and NGC 7469. We detect the signature of rotation in the second principal component of the 3.3 $μ$m PAH feature in all three targets. We construct velocity maps from the principal components for the 3.3 $μ$m PAH feature, Br$β$ (2.625$μ$m) and molecular hydrogen, H$_2$ 1-0 S(1) (2.12$μ$m). We find that in each target, the PAHs qualitatively follow the rotation of the galaxy, consistent with the rotational signature derived from Br$β$ and H$_2$. There are however some differences, with the PAH rotation in NGC 3256 N appearing at a different position angle, which suggest differences in the motion of the dust as compared to the gas. This kind of analysis opens a new window into this key component of the ISM.
△ Less
Submitted 22 May, 2024;
originally announced May 2024.
-
Dust beyond the torus: Revealing the mid-infrared heart of local Seyfert ESO 428-G14 with JWST/MIRI
Authors:
Houda Haidar,
David J. Rosario,
Almudena Alonso-Herrero,
Miguel Pereira-Santaella,
Ismael García-Bernete,
Stephanie Campbell,
Sebastian F. Hönig,
Cristina Ramos Almeida,
Erin Hicks,
Daniel Delaney,
Richard Davies,
Claudio Ricci,
Chris M. Harrison,
Mason Leist,
Enrique Lopez-Rodriguez,
Santiago Garcia-Burillo,
Lulu Zhang,
Chris Packham,
Poshak Gandhi,
Anelise Audibert,
Enrica Bellocchi,
Peter Boorman,
Andrew Bunker,
Françoise Combes,
Tanio Diaz Santos
, et al. (12 additional authors not shown)
Abstract:
Polar dust has been discovered in a number of local Active Galactic Nuclei (AGN), with radiation-driven torus models predicting a wind to be its main driver. However, little is known about its characteristics, spatial extent, or connection to the larger scale outflows. We present the first JWST/MIRI study aimed at imaging polar dust by zooming onto the heart of ESO 428-G14, part of the GATOS surve…
▽ More
Polar dust has been discovered in a number of local Active Galactic Nuclei (AGN), with radiation-driven torus models predicting a wind to be its main driver. However, little is known about its characteristics, spatial extent, or connection to the larger scale outflows. We present the first JWST/MIRI study aimed at imaging polar dust by zooming onto the heart of ESO 428-G14, part of the GATOS survey of local AGN. We clearly detect extended mid-infrared (MIR) emission within 200 pc from the nucleus. This polar structure is co-linear with a radio jet and lies perpendicular to a molecular gas lane that feeds and obscures the nucleus. The morphology of the MIR structure bears a striking resemblance to that of gas ionised by the AGN in the narrow-line region (NLR). We demonstrate that part of this spatial correspondence is due to contamination within the JWST filter bands from strong emission lines. Correcting for the contamination using ancillary spectroscopy, we find the morphology of the dust continuum to be asymmetric around the nucleus and more compact, though still clearly extended out to r ~ 100 pc. We estimate a temperature of the emitting dust of ~ 120 K. Using simple models, we find that the heating of small dust grains (~ 0.01 microns) by the radiation from the central AGN and/or radiative jet-induced shocks is responsible for the extended MIR emission. Large-grained dust, predicted by models of radiation-driven dusty winds from the torus, is unlikely to be important. This has important implications for scales to which AGN winds can carry dust and dense gas out into their host galaxies.
△ Less
Submitted 24 April, 2024;
originally announced April 2024.
-
Peeling Back the Layers of Extinction of Dusty Galaxies in the Era of JWST: Modelling Joint NIRSpec + MIRI Spectra at rest-frame 1.5-28 $μ$m
Authors:
F. R. Donnan,
I. García-Bernete,
D. Rigopoulou,
M. Pereira-Santaella,
P. F. Roche,
A. Alonso-Herrero
Abstract:
We present an analysis of the combined NIRSpec and MIRI spectra of dusty galaxies between 1.5 - 28 $μ$m restframe by implementing a differential extinction model, where the strength of extinction varies across the spectrum as different layers of the obscuring dust are probed. Our model is able to recover a 2D distribution of dust temperature and extinction which allows inference of the physical na…
▽ More
We present an analysis of the combined NIRSpec and MIRI spectra of dusty galaxies between 1.5 - 28 $μ$m restframe by implementing a differential extinction model, where the strength of extinction varies across the spectrum as different layers of the obscuring dust are probed. Our model is able to recover a 2D distribution of dust temperature and extinction which allows inference of the physical nature of the dust in these environments. We show that differential extinction is necessary to reproduce the spectra of 4 highly obscured Luminous Infrared Galaxies observed with NIRSpec IFU and MIRI MRS, where simple screen or uniformly mixed dust distributions fail to fit the data. We additionally compare the extinction of HII regions in these galaxies via hydrogen recombination lines, the extinction of molecular gas via the H$_2$ lines, Polycyclic Aromatic Hydrocarbons via the 12.7/11.3 PAH ratio and the stellar continuum. We find that the molecular gas is deeply buried with the HII regions in star-forming regions, with a similar extinction to the hottest dust components. However we find the cooler dust to be less obscured, at a similar extinction to the stellar continuum and PAHs. The nuclei show a complex dust distribution with VV114 NE, NGC 3256 S, IIZw96 SW showing a deeply buried continuum source relative to the molecular gas/HII regions. Additionally, NGC 3256 S, NGC 7469 and VV114 SW show an isolated hot dust component, indicative of AGN heating, where NGC 3256 S and NGC 7469 are previously known AGN.
△ Less
Submitted 27 February, 2024;
originally announced February 2024.
-
Structures Of Dust and gAs (SODA): Constraining the innermost dust properties of II Zw96 with JWST observations of H$_2$O and CO
Authors:
I. García-Bernete,
M. Pereira-Santaella,
E. González-Alfonso,
D. Rigopoulou,
A. Efstathiou,
F. R. Donnan,
N. Thatte
Abstract:
We analyze JWST NIRSpec$+$MIRI/MRS observations of the infrared (IR) gas-phase molecular bands of the most enshrouded source (D1) within the interacting system and luminous IR galaxy II Zw 096. We report the detection of rovibrational lines of H$_2$O $ν_2$=1-0 ($\sim$5.3-7.2 $μ$m) and $^{12}$CO $ν$=1-0 ($\sim$4.45-4.95 $μ$m) in D1. The CO band shows the R- and P-branches in emission and the spectr…
▽ More
We analyze JWST NIRSpec$+$MIRI/MRS observations of the infrared (IR) gas-phase molecular bands of the most enshrouded source (D1) within the interacting system and luminous IR galaxy II Zw 096. We report the detection of rovibrational lines of H$_2$O $ν_2$=1-0 ($\sim$5.3-7.2 $μ$m) and $^{12}$CO $ν$=1-0 ($\sim$4.45-4.95 $μ$m) in D1. The CO band shows the R- and P-branches in emission and the spectrum of the H$_2$O band shows the P-branch in emission and the R-branch in absorption. The H$_2$O R-branch in absorption unveils an IR-bright embedded compact source in D1 and the CO broad component features a highly turbulent environment. From both bands, we also identified extended intense star-forming (SF) activity associated with circumnuclear photodissociation regions (PDRs), consistent with the strong emission of the ionised 7.7 $μ$m polycyclic aromatic hydrocarbon band in this source. By including the 4.5-7.0 $μ$m continuum information derived from the H$_2$O and CO analysis, we modelled the IR emission of D1 with a dusty torus and SF component. The torus is very compact (diameter of $\sim$3 pc at 5 $μ$m) and characterised by warm dust ($\sim$ 370 K), giving an IR surface brightness of $\sim$3.6$\times$10$^{8}$ L$_{\rm sun}$/pc$^2$. This result suggests the presence of a dust-obscured active galactic nucleus (AGN) in D1, which has an exceptionally high covering factor that prevents the direct detection of AGN emission. Our results open a new way to investigate the physical conditions of inner dusty tori via modelling the observed IR molecular bands.
△ Less
Submitted 5 January, 2024;
originally announced January 2024.
-
JWST detection of extremely excited outflowing CO and H2O in VV 114 E SW: a possible rapidly accreting IMBH
Authors:
Eduardo González-Alfonso,
Ismael García-Bernete,
Miguel Pereira-Santaella,
David A. Neufeld,
Jacqueline Fischer,
Fergus R. Donnan
Abstract:
Mid-infrared (mid-IR) gas-phase molecular bands are powerful diagnostics of the warm interstellar medium. We report the James Webb Space Telescope detection of the CO v=1-0 (4.4-5.0 um) and H2O nu2=1-0 (5.0-7.8um) ro-vibrational bands, both in absorption, toward the ``s2'' core in the southwest nucleus of the merging galaxy VV 114 E. All ro-vibrational CO lines up to J_low=33 (E_low~3000 K) are de…
▽ More
Mid-infrared (mid-IR) gas-phase molecular bands are powerful diagnostics of the warm interstellar medium. We report the James Webb Space Telescope detection of the CO v=1-0 (4.4-5.0 um) and H2O nu2=1-0 (5.0-7.8um) ro-vibrational bands, both in absorption, toward the ``s2'' core in the southwest nucleus of the merging galaxy VV 114 E. All ro-vibrational CO lines up to J_low=33 (E_low~3000 K) are detected, as well as a forest of H2O lines up to 13_{0,13} (E_low~2600 K). The highest-excitation lines are blueshifted by ~180 km s^{-1} relative to the extended molecular cloud, which is traced by the rotational CO J=3-2 346 GHz line observed with the Atacama Large Millimeter/submillimeter Array. The bands also show absorption in a low-velocity component (blueshifted by ~30 km s^{-1}) with lower excitation. The analysis shows that the bands are observed against a continuum with effective temperature of T_bck~550 K extinguished with tau_6um^ext~ 2.5-3 (A_k~6.9-8.3 mag). The high-excitation CO and H2O lines are consistent with v=0 thermalization with T_rot~450 K and column densities of N_CO~(1.7-3.5)x10^{19} cm^{-2} and N_H2O~(1.5-3.0)x10^{19} cm$^{-2}$. Thermalization of the v=0 levels of H2O requires either an extreme density of n_H2>~10^9 cm^{-3}, or radiative excitation by the mid-IR field in a very compact (<1 pc) optically thick source emitting ~10^{10} L_sun. The latter alternative is favored, implying that the observed absorption probes the very early stages of a fully enshrouded active black hole (BH). On the basis of a simple model for BH growth and applying a lifetime constraint to the s2 core, an intermediate-mass BH (IMBH, M_BH~4.5x10^4 M_sun) accreting at super-Eddington rates is suggested, where the observed feedback has not yet been able to break through the natal cocoon.
△ Less
Submitted 8 December, 2023;
originally announced December 2023.
-
The Galaxy Activity, Torus, and Outflow Survey (GATOS). III: Revealing the inner icy structure in local AGN
Authors:
I. García-Bernete,
A. Alonso-Herrero,
D. Rigopoulou,
M. Pereira-Santaella,
T. Shimizu,
R. Davies,
F. R. Donnan,
P. F. Roche,
O. González-Martín,
C. Ramos Almeida,
E. Bellocchi,
P. Boorman,
F. Combes,
A. Efstathiou,
D. Esparza-Arredondo,
S. García-Burillo,
E. González-Alfonso,
E. K. S. Hicks,
S. Hönig,
A. Labiano,
N. A. Levenson,
E. López-Rodríguez,
C. Ricci,
C. Packham,
D. Rouan
, et al. (2 additional authors not shown)
Abstract:
We use JWST/MIRI MRS spectroscopy of a sample of six local obscured type 1.9/2 active galactic nuclei (AGN) to compare their nuclear mid-IR absorption bands with the level of nuclear obscuration traced by X-rays. This study is the first to use sub-arcsecond angular resolution data of local obscured AGN to investigate the nuclear mid-IR absorption bands with a wide wavelength coverage (4.9-28.1…
▽ More
We use JWST/MIRI MRS spectroscopy of a sample of six local obscured type 1.9/2 active galactic nuclei (AGN) to compare their nuclear mid-IR absorption bands with the level of nuclear obscuration traced by X-rays. This study is the first to use sub-arcsecond angular resolution data of local obscured AGN to investigate the nuclear mid-IR absorption bands with a wide wavelength coverage (4.9-28.1 $μ$m). All the nuclei show the 9.7 $μ$m silicate band in absorption. We compare the strength of the 9.7 and 18 $μ$m silicate features with torus model predictions. The observed silicate features are generally well explained by clumpy and smooth torus models. We report the detection of the 6 $μ$m dirty water ice band (i.e., a mix of water and other molecules such as CO and CO$_2$) at sub-arcsecond scales ($\sim$0.26 arcsec at 6 $μ$m; inner $\sim$50 pc) in a sample of local AGN with different levels of nuclear obscuration in the range log N$_{\rm H}^{\rm X-Ray}$(cm$^{-2}$)$\sim22-25$. We find a good correlation between the 6 $μ$m water ice optical depths and N$_{\rm H}^{\rm X-Ray}$. This result indicates that the water ice absorption might be a reliable tracer of the nuclear intrinsic obscuration in AGN. The weak water ice absorption in less obscured AGN (log N$_H^{X-ray}$ (cm$^{-2}$)$\lesssim$23.0 cm$^{-2}$) might be related to the hotter dust temperature ($>$T$_{sub}^{H_2O}\sim$110 K) expected to be reached in the outer layers of the torus due to their more inhomogeneous medium. Our results suggest it might be necessary to include the molecular content, such as, H$_2$O, aliphatic hydrocarbons (CH-) and more complex PAH molecules in torus models to better constrain key parameters such as the torus covering factor (i.e. nuclear obscuration).
△ Less
Submitted 4 December, 2023; v1 submitted 13 October, 2023;
originally announced October 2023.
-
AGN STORM 2. IV. Swift X-ray and ultraviolet/optical monitoring of Mrk 817
Authors:
Edward M. Cackett,
Jonathan Gelbord,
Aaron J. Barth,
Gisella De Rosa,
Rick Edelson,
Michael R. Goad,
Yasaman Homayouni,
Keith Horne,
Erin A. Kara,
Gerard A. Kriss,
Kirk T. Korista,
Hermine Landt,
Rachel Plesha,
Nahum Arav,
Misty C. Bentz,
Benjamin D. Boizelle,
Elena Dalla Bonta,
Maryam Dehghanian,
Fergus Donnan,
Pu Du,
Gary J. Ferland,
Carina Fian,
Alexei V. Filippenko,
Diego H. Gonzalez Buitrago,
Catherine J. Grier
, et al. (26 additional authors not shown)
Abstract:
The AGN STORM 2 campaign is a large, multiwavelength reverberation mapping project designed to trace out the structure of Mrk 817 from the inner accretion disk to the broad emission line region and out to the dusty torus. As part of this campaign, Swift performed daily monitoring of Mrk 817 for approximately 15 months, obtaining observations in X-rays and six UV/optical filters. The X-ray monitori…
▽ More
The AGN STORM 2 campaign is a large, multiwavelength reverberation mapping project designed to trace out the structure of Mrk 817 from the inner accretion disk to the broad emission line region and out to the dusty torus. As part of this campaign, Swift performed daily monitoring of Mrk 817 for approximately 15 months, obtaining observations in X-rays and six UV/optical filters. The X-ray monitoring shows that Mrk 817 was in a significantly fainter state than in previous observations, with only a brief flare where it reached prior flux levels. The X-ray spectrum is heavily obscured. The UV/optical light curves show significant variability throughout the campaign and are well correlated with one another, but uncorrelated with the X-rays. Combining the Swift UV/optical light curves with Hubble UV continuum light curves, we measure interband continuum lags, $τ(λ)$, that increase with increasing wavelength roughly following $τ(λ) \propto λ^{4/3}$, the dependence expected for a geometrically thin, optically thick, centrally illuminated disk. Modeling of the light curves reveals a period at the beginning of the campaign where the response of the continuum is suppressed compared to later in the light curve - the light curves are not simple shifted and scaled versions of each other. The interval of suppressed response corresponds to a period of high UV line and X-ray absorption, and reduced emission line variability amplitudes. We suggest that this indicates a significant contribution to the continuum from the broad line region gas that sees an absorbed ionizing continuum.
△ Less
Submitted 26 September, 2023; v1 submitted 30 June, 2023;
originally announced June 2023.
-
The Sloan Digital Sky Survey Reverberation Mapping Project: Key Results
Authors:
Yue Shen,
Catherine J. Grier,
Keith Horne,
Zachary Stone,
Jennifer I. Li,
Qian Yang,
Yasaman Homayouni,
Jonathan R. Trump,
Scott F. Anderson,
W. N. Brandt,
Patrick B. Hall,
Luis C. Ho,
Linhua Jiang,
Patrick Petitjean,
Donald P. Schneider,
Charling Tao,
Fergus. R. Donnan,
Yusra AlSayyad,
Matthew A. Bershady,
Michael R. Blanton,
Dmitry Bizyaev,
Kevin Bundy,
Yuguang Chen,
Megan C. Davis,
Kyle Dawson
, et al. (22 additional authors not shown)
Abstract:
We present the final data from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, a precursor to the SDSS-V Black Hole Mapper Reverberation Mapping program. This data set includes 11-year photometric and 7-year spectroscopic light curves for 849 broad-line quasars over a redshift range of 0.1<z<4.5 and a luminosity range of Lbol=1E44-47.5 erg/s, along with spectral and variabili…
▽ More
We present the final data from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, a precursor to the SDSS-V Black Hole Mapper Reverberation Mapping program. This data set includes 11-year photometric and 7-year spectroscopic light curves for 849 broad-line quasars over a redshift range of 0.1<z<4.5 and a luminosity range of Lbol=1E44-47.5 erg/s, along with spectral and variability measurements. We report 23, 81, 125, and 110 reverberation mapping lags (relative to optical continuum variability) for broad Halpha, Hbeta, MgII and CIV using the SDSS-RM sample, spanning much of the luminosity and redshift ranges of the sample. Using 30 low-redshift RM AGNs with dynamical-modeling black hole masses, we derive a new estimate of the average virial factor of <log f>=0.62+-0.07 for the line dispersion measured from the RMS spectrum. The intrinsic scatter of individual virial factors is 0.31+-0.07 dex, indicating a factor of two systematic uncertainty in RM black hole masses. Our lag measurements reveal significant R-L relations for Hbeta and MgII at high redshift, consistent with the latest measurements based on heterogeneous samples. While we are unable to robustly constrain the slope of the R-L relation for CIV given the limited dynamical range in luminosity, we found substantially larger scatter in CIV lags at fixed L1350. Using the SDSS-RM lag sample, we derive improved single-epoch (SE) mass recipes for Hbeta, MgII and CIV, which are consistent with their respective RM masses as well as between the SE recipes from two different lines, over the luminosity range probed by our sample. The new Hbeta and MgII recipes are approximately unbiased estimators at given RM masses, but there are systematic biases in the CIV recipe. The intrinsic scatter of SE masses around RM masses is ~0.45 dex for Hbeta and MgII, increasing to ~0.58 dex for CIV.
△ Less
Submitted 1 April, 2024; v1 submitted 1 May, 2023;
originally announced May 2023.
-
Testing Super-Eddington Accretion onto a Supermassive Black Hole: Reverberation Mapping of PG 1119+120
Authors:
Fergus R. Donnan,
Juan V. Hernández Santisteban,
Keith Horne,
Chen Hu,
Pu Du,
Yan-Rong Li,
Ming Xiao,
Luis C. Ho,
Jesús Aceituno,
Jian-Min Wang,
Wei-Jian Guo,
Sen Yang,
Bo-Wei Jiang,
Zhu-Heng Yao
Abstract:
We measure the black hole mass and investigate the accretion flow around the local ($z=0.0502$) quasar PG 1119+120. Spectroscopic monitoring with Calar Alto provides H$β$ lags and linewidths from which we estimate a black hole mass of $\log \left(M_{\bullet}/\mathrm{M}_{\odot} \right) = 7.0$, uncertain by $\sim0.4$ dex. High cadence photometric monitoring over two years with the Las Cumbres Observ…
▽ More
We measure the black hole mass and investigate the accretion flow around the local ($z=0.0502$) quasar PG 1119+120. Spectroscopic monitoring with Calar Alto provides H$β$ lags and linewidths from which we estimate a black hole mass of $\log \left(M_{\bullet}/\mathrm{M}_{\odot} \right) = 7.0$, uncertain by $\sim0.4$ dex. High cadence photometric monitoring over two years with the Las Cumbres Observatory provides lightcurves in 7 optical bands suitable for intensive continuum reverberation mapping. We identify variability on two timescales. Slower variations on a 100-day timescale exhibit excess flux and increased lag in the $u'$ band and are thus attributable to diffuse bound-free continuum emission from the broad line region. Faster variations that we attribute to accretion disc reprocessing lack a $u'$-band excess and have flux and delay spectra consistent with either $τ\propto λ^{4/3}$, as expected for a temperature structure of $T(R) \propto R^{-3/4}$ for a thin accretion disc, or $τ\propto λ^{2}$ expected for a slim disc. Decomposing the flux into variable (disc) and constant (host galaxy) components, we find the disc SED to be flatter than expected with $f_ν \sim \rm{const}$. Modelling the SED predicts an Eddington ratio of $λ_{\rm Edd} > 1$, where the flat spectrum can be reproduced by a slim disc with little dust extinction or a thin disc which requires more dust extinction. While this accretion is super-Eddington, the geometry is still unclear, however a slim disc is expected due to the high radiation pressure at these accretion rates, and is entirely consistent with our observations.
△ Less
Submitted 16 May, 2023; v1 submitted 18 February, 2023;
originally announced February 2023.
-
A Detailed Look at the Most Obscured Galactic Nuclei in the Mid-Infrared
Authors:
F. R. Donnan,
D. Rigopoulou,
I. García-Bernete,
M. Pereira-Santaella,
A. Alonso-Herrero,
P. F. Roche,
S. Aalto,
A. Hernán-Caballero,
H. W. W. Spoon
Abstract:
Context. Compact Obscured Nuclei (CONs) are an extreme phase of galaxy evolution where rapid supermassive black hole growth and$/$or compact star-forming activity is completely obscured by gas and dust. Aims. We investigate the properties of CONs in the mid-infrared and explore techniques aimed at identifying these objects such as through the equivalent width (EW) ratios of their Polycyclic Aromat…
▽ More
Context. Compact Obscured Nuclei (CONs) are an extreme phase of galaxy evolution where rapid supermassive black hole growth and$/$or compact star-forming activity is completely obscured by gas and dust. Aims. We investigate the properties of CONs in the mid-infrared and explore techniques aimed at identifying these objects such as through the equivalent width (EW) ratios of their Polycyclic Aromatic Hydrocarbon (PAH) features. Methods. We model Spitzer spectra by decomposing the continua into nuclear and star-forming components from which we then measure the nuclear optical depth, $τ_N$, of the $9.8 μ$m silicate absorption feature. We also use Spitzer spectral maps to investigate how PAH EW ratios vary with aperture size for objects hosting CONs. Results. We find that the nuclear optical depth, $τ_N$, strongly correlates with the HCN-vib emission line in the millimetre for CONs with a Pearson correlation coefficient of 0.91. We find the PAH EW ratios technique to be effective at selecting CONs and robust against highly inclined galaxies where strong dust lanes may mimic a CON like spectrum by producing a high $τ_N$. Our analysis of the Spitzer spectral maps showed that the efficacy of the PAH EW ratios to isolate CONs is reduced when there is a strong star-forming component from the host galaxy. In addition, we find that the use of the inferred nuclear optical depth is a reliable method to identify CONs in $36^{+8}_{-7}\%$ of ULIRGs and $17^{+3}_{-3}\%$ of LIRGs, consistent with previous work. Conclusions. We confirm mid-IR spectra to be a powerful diagnostic of CONs where the increased sensitivity of JWST will allow identification of CONs at cosmic noon revealing this extreme but hidden phase of galaxy evolution.
△ Less
Submitted 17 November, 2022;
originally announced November 2022.
-
The Obscured Nucleus and Shocked Environment of VV 114E Revealed by JWST/MIRI Spectroscopy
Authors:
F. R. Donnan,
I. García-Bernete,
D. Rigopoulou,
M. Pereira-Santaella,
A. Alonso-Herrero,
P. F. Roche,
A. Hernán-Caballero,
H. W. W. Spoon
Abstract:
Compact Obscured Nuclei (CONs) potentially hide extreme supermassive black hole (SMBH) growth behind large column densities of gas/dust. We present a spectroscopic analysis of the heavily obscured nucleus and the surrounding environment of the eastern region of the nearby ($z = 0.02007$) interacting galaxy VV 114 with the JWST Mid-InfraRed Instrument (MIRI). We model the spectrum from 4.9 - 28…
▽ More
Compact Obscured Nuclei (CONs) potentially hide extreme supermassive black hole (SMBH) growth behind large column densities of gas/dust. We present a spectroscopic analysis of the heavily obscured nucleus and the surrounding environment of the eastern region of the nearby ($z = 0.02007$) interacting galaxy VV 114 with the JWST Mid-InfraRed Instrument (MIRI). We model the spectrum from 4.9 - 28 $μ$m to extract Polycyclic Aromatic Hydrocarbon (PAH) emission and the underlying obscured continuum. We find that the NE nucleus (A) is highly obscured where the low PAH equivalent width (EW) ratio, EW(12.7)/EW(11.3), reveals a dust enshrouded continuum source. This is confirmed by decomposing the continuum into nuclear and star-forming where the nuclear component is found to be typical of CONs. The 11.3/6.2 PAH flux ratio is consistent with originating in star-forming regions rather than typical AGN. The second nucleus (B) is much less obscured, with PAH flux ratios also typical of star-forming regions. We do not detect any high ionisation lines such as [Ne V] or [Ne VI] which suggests that if an AGN is present it must be highly obscured. Additionally, we detect a shock front south of the secondary nucleus (B) in the [Fe II] (5.34 $μ$m) line and in warm molecular hydrogen. The 6.2 PAH emission does not spatially coincide with the low-J transitions of H$_2$ but rather appears strong at the shock front which may suggest destruction of the ionised PAHs in the post-shock gas behind the shock front.
△ Less
Submitted 15 December, 2022; v1 submitted 10 October, 2022;
originally announced October 2022.
-
A high angular resolution view of the PAH emission in Seyfert galaxies using JWST/MRS data
Authors:
I. García-Bernete,
D. Rigopoulou,
A. Alonso-Herrero,
F. R. Donnan,
P. F. Roche,
M. Pereira-Santaella,
A. Labiano,
L. Peralta de Arriba,
T. Izumi,
C. Ramos Almeida,
T. Shimizu,
S. Hönig,
S. García-Burillo,
D. J. Rosario,
M. J. Ward,
E. Bellocchi,
E. K. S. Hicks,
L. Fuller,
C. Packham
Abstract:
Polycyclic Aromatic Hydrocarbons (PAHs) are carbon-based molecules that are ubiquitous in a variety of astrophysical objects and environments. In this work, we use JWST/MIRI MRS spectroscopy of three Seyferts to compare their nuclear PAH emission with that of star-forming regions. This study represents the first of its kind using sub-arcsecond angular resolution data of local luminous Seyferts (Lb…
▽ More
Polycyclic Aromatic Hydrocarbons (PAHs) are carbon-based molecules that are ubiquitous in a variety of astrophysical objects and environments. In this work, we use JWST/MIRI MRS spectroscopy of three Seyferts to compare their nuclear PAH emission with that of star-forming regions. This study represents the first of its kind using sub-arcsecond angular resolution data of local luminous Seyferts (Lbol>10^44.46 erg/s) on a wide wavelength coverage (4.9-28.1 micron). We present an analysis of their nuclear PAH properties by comparing the observed ratios with PAH diagnostic model grids, derived from theoretical spectra. Our results show that a suite of PAH features is present in the innermost parts (~0.45 arcsec at 12 micron; in the inner ~142-245 pc) of luminous Seyfert galaxies. We find that the nuclear regions of AGN lie at different positions of the PAH diagnostic diagrams, whereas the SF regions are concentrated around the average values of SF galaxies. In particular, we find that the nuclear PAH emission mainly originates in neutral PAHs. In contrast, PAH emission originating in the SF regions favours ionised PAH grains. The observed PAH ratios in the nuclear region of AGN-dominated galaxy NGC 6552 indicate the presence of larger-sized PAH molecules compared with those of the SF regions. Therefore, our results provide evidence that the AGN have a significant impact on the ionization state (and probably the size) of the PAH grains on scales of ~142-245 pc.
△ Less
Submitted 28 September, 2022; v1 submitted 24 August, 2022;
originally announced August 2022.
-
Bayesian Analysis of Quasar Lightcurves with a Running Optimal Average: New Time Delay Measurements of COSMOGRAIL Gravitationally Lensed Quasars
Authors:
Fergus R. Donnan,
Keith Horne,
Juan V. Hernández Santisteban
Abstract:
We present a new method of modelling time-series data based on the running optimal average (ROA). By identifying the effective number of parameters for the ROA model, in terms of the shape and width of its window function and the times and accuracies of the data, we enable a Bayesian analysis, optimising the ROA width, along with other model parameters, by minimising the Bayesian Information Crite…
▽ More
We present a new method of modelling time-series data based on the running optimal average (ROA). By identifying the effective number of parameters for the ROA model, in terms of the shape and width of its window function and the times and accuracies of the data, we enable a Bayesian analysis, optimising the ROA width, along with other model parameters, by minimising the Bayesian Information Criterion (BIC) and sampling joint posterior parameter distributions using MCMC methods. For analysis of quasar lightcurves, our implementation of ROA modelling can measure time delays among lightcurves at different wavelengths or from different images of a lensed quasar and, in future work, be used to inter-calibrate lightcurve data from different telescopes and estimate the shape and thus the power-density spectrum of the lightcurve. Our noise model implements a robust treatment of outliers and error-bar adjustments to account for additional variance or poorly-quantified uncertainties. Tests with simulated data validate the parameter uncertainty estimates. We compare ROA delay measurements with results from cross-correlation and from JAVELIN, which models lightcurves with a prior on the power-density spectrum. We analyse published COSMOGRAIL lightcurves of multi-lensed quasar lightcurves and present the resulting measurements of the inter-image time delays and detection of microlensing effects.
△ Less
Submitted 7 October, 2021; v1 submitted 26 July, 2021;
originally announced July 2021.
-
The asymmetric inner disk of the Herbig Ae star HD 163296 in the eyes of VLTI/MATISSE: evidence for a vortex?
Authors:
J. Varga,
M. Hogerheijde,
R. van Boekel,
L. Klarmann,
R. Petrov,
L. B. F. M. Waters,
S. Lagarde,
E. Pantin,
Ph. Berio,
G. Weigelt,
S. Robbe-Dubois,
B. Lopez,
F. Millour,
J. -C. Augereau,
H. Meheut,
A. Meilland,
Th. Henning,
W. Jaffe,
F. Bettonvil,
P. Bristow,
K. -H. Hofmann,
A. Matter,
G. Zins,
S. Wolf,
F. Allouche
, et al. (111 additional authors not shown)
Abstract:
Context. The inner few au region of planet-forming disks is a complex environment. High angular resolution observations have a key role in understanding the disk structure and the dynamical processes at work. Aims. In this study we aim to characterize the mid-infrared brightness distribution of the inner disk of the young intermediate-mass star HD 163296, from VLTI/MATISSE observations. Methods. W…
▽ More
Context. The inner few au region of planet-forming disks is a complex environment. High angular resolution observations have a key role in understanding the disk structure and the dynamical processes at work. Aims. In this study we aim to characterize the mid-infrared brightness distribution of the inner disk of the young intermediate-mass star HD 163296, from VLTI/MATISSE observations. Methods. We use geometric models to fit the data. Our models include a smoothed ring, a flat disk with inner cavity, and a 2D Gaussian. The models can account for disk inclination and for azimuthal asymmetries as well. We also perform numerical hydro-dynamical simulations of the inner edge of the disk. Results. Our modeling reveals a significant brightness asymmetry in the L-band disk emission. The brightness maximum of the asymmetry is located at the NW part of the disk image, nearly at the position angle of the semimajor axis. The surface brightness ratio in the azimuthal variation is $3.5 \pm 0.2$. Comparing our result on the location of the asymmetry with other interferometric measurements, we confirm that the morphology of the $r<0.3$ au disk region is time-variable. We propose that this asymmetric structure, located in or near the inner rim of the dusty disk, orbits the star. For the physical origin of the asymmetry, we tested a hypothesis where a vortex is created by Rossby wave instability, and we find that a unique large scale vortex may be compatible with our data. The half-light radius of the L-band emitting region is $0.33\pm 0.01$ au, the inclination is ${52^\circ}^{+5^\circ}_{-7^\circ}$, and the position angle is $143^\circ \pm 3^\circ$. Our models predict that a non-negligible fraction of the L-band disk emission originates inside the dust sublimation radius for $μ$m-sized grains. Refractory grains or large ($\gtrsim 10\ μ$m-sized) grains could be the origin for this emission.
△ Less
Submitted 10 December, 2020;
originally announced December 2020.