Showing 1–50 of 99 results for author: Levenson, N

GATOS: missing molecular gas in the outflow of NGC5728 revealed by JWST

Authors: R. Davies, T. Shimizu, M. Pereira-Santaella, A. Alonso-Herrero, A. Audibert, E. Bellocchi, P. Boorman, S. Campbell, Y. Cao, F. Combes, D. Delaney, T. Diaz-Santos, F. Eisenhauer, D. Esparza Arredondo, H. Feuchtgruber, N. M. Forster Schreiber, L. Fuller, P. Gandhi, I. Garcia-Bernete, S. Garcia-Burillo, B. Garcia-Lorenzo, R. Genzel, S. Gillessen, O. Gonzalez Martin, H. Haidar , et al. (27 additional authors not shown)

Abstract: The ionisation cones of NGC5728 have a deficit of molecular gas based on millimetre observations of CO(2-1) emission. Although photoionisation from the active nucleus may lead to suppression of this transition, warm molecular gas can still be present. We report the detection of eight mid-infrared rotational H$_2$ lines throughout the central kiloparsec, including the ionisation cones, using integr… ▽ More The ionisation cones of NGC5728 have a deficit of molecular gas based on millimetre observations of CO(2-1) emission. Although photoionisation from the active nucleus may lead to suppression of this transition, warm molecular gas can still be present. We report the detection of eight mid-infrared rotational H$_2$ lines throughout the central kiloparsec, including the ionisation cones, using integral field spectroscopic observations with JWST/MIRI MRS. The H$_2$ line ratios, characteristic of a power-law temperature distribution, indicate that the gas is warmest where it enters the ionisation cone through disk rotation, suggestive of shock excitation. In the nucleus, where the data can be combined with an additional seven ro-vibrational H$_2$ transitions, we find that moderate velocity (30 km s$^{-1}$) shocks in dense ($10^5$ cm$^{-3}$) gas, irradiated by an external UV field ($G_0 = 10^3$), do provide a good match to the full set. The warm molecular gas in the ionisation cone that is traced by the H$_2$ rotational lines has been heated to temperatures $>200$ K. Outside of the ionisation cone the molecular gas kinematics are undisturbed. However, within the ionisation cone, the kinematics are substantially perturbed, indicative of a radial flow, but one that is quantitatively different from the ionised lines. We argue that this outflow is in the plane of the disk, implying a short 50 pc acceleration zone up to speeds of about 400 km s$^{-1}$ followed by an extended deceleration over $\sim$700 pc where it terminates. The deceleration is due to both the radially increasing galaxy mass, and mass-loading as ambient gas in the disk is swept up. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: A&A accepted; 16 pages

Deciphering the imprint of AGN feedback in Seyfert galaxies: Nuclear-scale molecular gas deficits

Authors: S. García-Burillo, E. K. S. Hicks, A. Alonso-Herrero, M. Pereira-Santaella, A. Usero, M. Querejeta, O. González-Martin, D. Delaney, C. Ramos Almeida, F. Combes, D. Anglés-Alcázar, A. Audibert, E. Bellocchi, R. I. Davies, T. A. Davis, J. S. Elford, I. García-Bernete, S. Hönig, A. Labiano, M. T. Leist, N. A. Levenson, E. López-Rodríguez, J. Mercedes-Feliz, C. Packham, C. Ricci , et al. (4 additional authors not shown)

Abstract: We use a sample of 64 nearby (D=7-45 Mpc) disk galaxies including 45 AGN and 19 non-AGN, that have high spatial resolution multiline CO observations obtained with the ALMA and/or PdBI arrays to study the distribution of cold molecular gas in their circumunuclear disks (CND). We analyze whether the concentration of cold molecular gas changes as a function of the X-ray luminosity in the 2-10 keV ran… ▽ More We use a sample of 64 nearby (D=7-45 Mpc) disk galaxies including 45 AGN and 19 non-AGN, that have high spatial resolution multiline CO observations obtained with the ALMA and/or PdBI arrays to study the distribution of cold molecular gas in their circumunuclear disks (CND). We analyze whether the concentration of cold molecular gas changes as a function of the X-ray luminosity in the 2-10 keV range ($L_{\rm X}$). We also study the concentration of the hot molecular gas using NIR data obtained for the H2 1-0S(1) line. We find a turnover in the distribution of the cold molecular gas concentration as a function of $L_{\rm X}$ with a breakpoint which divides the sample into two branches: the AGN build-up branch ($L_{\rm X}\leq10^{41.5\pm0.3}$erg/s) and the AGN feedback branch ($L_{\rm X}\geq10^{41.5\pm0.3}$erg/s) . Lower luminosity AGN and non-AGN of the AGN build-up branch show high cold molecular gas concentrations and centrally peaked radial profiles on nuclear ($r\leq50$~pc) scales. Higher luminosity AGN of the AGN feedback branch, show a sharp decrease in the concentration of molecular gas and flat or inverted radial profiles. The cold molecular gas concentration index ($CCI$), defined as the ratio of surface densities at $r\leq50$~pc and $r\leq200$~pc , namely $CCI \equiv$~log$_{\rm 10}(Σ^{\rm gas}_{\rm 50}/Σ^{\rm gas}_{\rm 200}$), spans a factor ~4-5 between the galaxies lying at the high end of the AGN build-up branch and the galaxies of the AGN feedback branch. The concentration and radial distributions of the hot molecular gas in our sample follow less extreme trends as a function of the X-ray luminosity. These observations confirm, on a three times larger sample, previous evidence found by the GATOS survey that the imprint of AGN feedback on the CND-scale distribution of molecular gas is more extreme in higher luminosity Seyfert galaxies of the local universe. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: Accepted by Astronomy & Astrophysics (A&A) (14/06/2024), 26 pages, 16 figures

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.

Comments: 17 pages, 8 figures, 1 table. Submitted to MNRAS, comments are welcome!

AGN feedback in the Local Universe: multiphase outflow of the Seyfert galaxy NGC 5506

Authors: Federico Esposito, Almudena Alonso-Herrero, Santiago García-Burillo, Viviana Casasola, Françoise Combes, Daniele Dallacasa, Richard Davies, Ismael García-Bernete, Begoña García-Lorenzo, Laura Hermosa Muñoz, Luis Peralta de Arriba, Miguel Pereira-Santaella, Francesca Pozzi, Cristina Ramos Almeida, Thomas Taro Shimizu, Livia Vallini, Enrica Bellocchi, Omaira González-Martín, Erin K. S. Hicks, Sebastian Hönig, Alvaro Labiano, Nancy A. Levenson, Claudio Ricci, David J. Rosario

Abstract: We present new optical GTC/MEGARA seeing-limited (0.9") integral-field observations of NGC 5506, together with ALMA observations of the CO(3-2) transition at a 0.2" (25 pc) resolution. NGC 5506 is a luminous (bolometric luminosity of $\sim 10^{44}$ erg/s) nearby (26 Mpc) Seyfert galaxy, part of the Galaxy Activity, Torus, and Outflow Survey (GATOS). We modelled the CO(3-2) kinematics with 3D-Barol… ▽ More We present new optical GTC/MEGARA seeing-limited (0.9") integral-field observations of NGC 5506, together with ALMA observations of the CO(3-2) transition at a 0.2" (25 pc) resolution. NGC 5506 is a luminous (bolometric luminosity of $\sim 10^{44}$ erg/s) nearby (26 Mpc) Seyfert galaxy, part of the Galaxy Activity, Torus, and Outflow Survey (GATOS). We modelled the CO(3-2) kinematics with 3D-Barolo, revealing a rotating and outflowing cold gas ring within the central 1.2 kpc. We derived an integrated cold molecular gas mass outflow rate for the ring of 8 M$_{\odot}$/yr. We fitted the optical emission lines with a maximum of two Gaussian components to separate rotation from non-circular motions. We detected high [OIII]$λ$5007 projected velocities (up to 1000 km/s) at the active galactic nucleus (AGN) position, decreasing with radius to an average 330 km/s around 350 pc. We also modelled the [OIII] gas kinematics with a non-parametric method, estimating the ionisation parameter and electron density in every spaxel, from which we derived an ionised mass outflow rate of 0.076 M$_{\odot}$/yr within the central 1.2 kpc. Regions of high CO(3-2) velocity dispersion, extending to projected distances of 350 pc from the AGN, appear to be the result from the interaction of the AGN wind with molecular gas in the galaxy's disc. Additionally, we find the ionised outflow to spatially correlate with radio and soft X-ray emission in the central kiloparsec. We conclude that the effects of AGN feedback in NGC 5506 manifest as a large-scale ionised wind interacting with the molecular disc, resulting in outflows extending to radial distances of 610 pc △ Less

Submitted 6 March, 2024; originally announced March 2024.

Comments: 24 pages, 30 figures, accepted for publication by Astronomy & Astrophysics

Deconvolution of JWST/MIRI Images: Applications to an AGN Model and GATOS Observations of NGC 5728

Authors: M. T. Leist, C. Packham, D. J. V. Rosario, D. A. Hope, A. Alonso-Herrero, E. K. S. Hicks, S. Hönig, L. Zhang, R. Davies, T. Díaz-Santos, O. Ganzález-Martín, E. Bellocchi, P. G. Boorman, F. Combes, I. García-Bernete, S. García-Burillo, B. García-Lorenzo, H. Haidar, K. Ichikawa, M. Imanishi, S. M. Jefferies, Á. Labiano, N. A. Levenson, R. Nikutta, M. Pereira-Santaella , et al. (10 additional authors not shown)

Abstract: The superb image quality, stability and sensitivity of the JWST permit deconvolution techniques to be pursued with a fidelity unavailable to ground-based observations. We present an assessment of several deconvolution approaches to improve image quality and mitigate effects of the complex JWST point spread function (PSF). The optimal deconvolution method is determined by using WebbPSF to simulate… ▽ More The superb image quality, stability and sensitivity of the JWST permit deconvolution techniques to be pursued with a fidelity unavailable to ground-based observations. We present an assessment of several deconvolution approaches to improve image quality and mitigate effects of the complex JWST point spread function (PSF). The optimal deconvolution method is determined by using WebbPSF to simulate JWST's complex PSF and MIRISim to simulate multi-band JWST/Mid-Infrared Imager Module (MIRIM) observations of a toy model of an active galactic nucleus (AGN). Five different deconvolution algorithms are tested: (1) Kraken deconvolution, (2) Richardson-Lucy, (3) Adaptive Imaging Deconvolution Algorithm, (4) Sparse regularization with the Condat-Vũ algorithm, and (5) Iterative Wiener Filtering and Thresholding. We find that Kraken affords the greatest FWHM reduction of the nuclear source of our MIRISim observations for the toy AGN model while retaining good photometric integrity across all simulated wavebands. Applying Kraken to Galactic Activity, Torus, and Outflow Survey (GATOS) multi-band JWST/MIRIM observations of the Seyfert 2 galaxy NGC 5728, we find that the algorithm reduces the FWHM of the nuclear source by a factor of 1.6-2.2 across all five filters. Kraken images facilitate detection of a SE to NW $\thicksim$2".5 ($\thicksim$470 pc, PA $\simeq$115°) extended nuclear emission, especially in the longest wavelengths. We demonstrate that Kraken is a powerful tool to enhance faint features otherwise hidden in the complex JWST PSF. △ Less

Submitted 14 February, 2024; v1 submitted 21 December, 2023; originally announced December 2023.

Comments: 32 pages, 23 figures, published in AJ 2024 February 7

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.

Comments: Accepted for publication in A&A. 13 pages, 12 Figures

AGN feedback in action in the molecular gas ring of the Seyfert galaxy NGC7172

Authors: A. Alonso-Herrero, S. Garcia-Burillo, M. Pereira-Santaella, T. Shimizu, F. Combes, E. K. S. Hicks, R. Davies, C. Ramos Almeida, I. Garcia-Bernete, S. F. Hoenig, N. A. Levenson, C. Packham, E. Bellocchi, L. K. Hunt, M. Imanishi, C. Ricci, P. Roche

Abstract: We present new ALMA observations of the CO(3-2) transition and 854micron continuum at 0.06-0.3" resolution, together with new VLT/SINFONI observations of NGC7172. This is a luminous (bolometric luminosity of ~10^44 erg/s) Seyfert galaxy that belongs to the Galaxy Activity, Torus, and Outflow Survey (GATOS). The CO(3-2) observations reveal the presence of a highly inclined cold molecular gas ring w… ▽ More We present new ALMA observations of the CO(3-2) transition and 854micron continuum at 0.06-0.3" resolution, together with new VLT/SINFONI observations of NGC7172. This is a luminous (bolometric luminosity of ~10^44 erg/s) Seyfert galaxy that belongs to the Galaxy Activity, Torus, and Outflow Survey (GATOS). The CO(3-2) observations reveal the presence of a highly inclined cold molecular gas ring with an approximate radius of 3-4"~540-720 pc, which is likely associated with an inner Lindblad resonance of a putative stellar bar. There are noncircular motions in the VLT/SINFONI [SiVI]1.96micron and H2 at 2.12micron, and ALMA CO(3-2) velocity fields. After subtracting the stellar velocity field, we detected [SiVI] blueshifted velocities of a few hundred km/s to the south of the AGN. They trace outflowing ionized gas outside the plane of the galaxy and out to projected distances of ~200 pc. The CO(3-2) position-velocity diagram along the kinematic minor axis displays noncircular motions with observed velocities of up to ~150 km/s. Assuming that these are taking place in the disk of the galaxy, the observed velocity signs imply that the molecular gas ring is not only rotating but also outflowing. We derived an integrated cold molecular gas mass outflow rate of ~40 Msun/yr for the ring. Using the 854micron map, we resolved a 32 pc radius torus with a gas mass of 8x10^5 Msun. These torus properties are similar to other Seyfert galaxies in the GATOS sample. We measured a decreased cold molecular gas concentration in the nuclear-torus region relative to the circumnuclear region when compared to other less luminous Seyfert galaxies. We conclude that the effects of AGN feedback in NGC7172, which are likely caused by the AGN wind and/or the moderate luminosity radio jet, are seen as a large-scale outflowing molecular gas ring and accompanying redistribution of molecular gas in the nuclear regions. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: Accepted for publication to A&A

A radio-jet driven outflow in the Seyfert 2 galaxy NGC 2110?

Authors: L. Peralta de Arriba, A. Alonso-Herrero, S. García-Burillo, I. García-Bernete, M. Villar-Martín, B. García-Lorenzo, R. Davies, D. J. Rosario, S. F. Hönig, N. A. Levenson, C. Packham, C. Ramos Almeida, M. Pereira-Santaella, A. Audibert, E. Bellocchi, E. K. S. Hicks, A. Labiano, C. Ricci, D. Rigopoulou

Abstract: We present a spatially-resolved study of the ionised gas in the central 2 kpc of the Seyfert 2 galaxy NGC 2110 and investigate the role of its moderate luminosity radio jet (kinetic radio power of $P_\mathrm{jet} = 2.3 \times 10^{43}\mathrm{erg\ s^{-1}}$). We use new optical integral-field observations taken with the MEGARA spectrograph at GTC. We fit the emission lines with a maximum of two Gauss… ▽ More We present a spatially-resolved study of the ionised gas in the central 2 kpc of the Seyfert 2 galaxy NGC 2110 and investigate the role of its moderate luminosity radio jet (kinetic radio power of $P_\mathrm{jet} = 2.3 \times 10^{43}\mathrm{erg\ s^{-1}}$). We use new optical integral-field observations taken with the MEGARA spectrograph at GTC. We fit the emission lines with a maximum of two Gaussian components, except at the AGN position where we used three. Aided by existing stellar kinematics, we use the observed velocity and velocity dispersion of the emission lines to classify the different kinematic components. The disc component is characterised by lines with $σ\sim 60-200\ \mathrm{km\ s^{-1}}$. The outflow component has typical values of $σ\sim 700\ \mathrm{km\ s^{-1}}$ and is confined to the central 400 pc, which is coincident with linear part of the radio jet detected in NGC 2110. At the AGN position, the [O III]$λ$5007 line shows high velocity components reaching at least $1000\ \mathrm{km\ s^{-1}}$. This and the high velocity dispersions indicate the presence of outflowing gas outside the galaxy plane. Spatially-resolved diagnostic diagrams reveal mostly LI(N)ER-like excitation in the outflow and some regions in the disc, which could be due to the presence of shocks. However, there is also Seyfert-like excitation beyond the bending of the radio jet, probably tracing the edge of the ionisation cone that intercepts with the disc of the galaxy. NGC 2110 follows well the observational trends between the outflow properties and the jet radio power found for a few nearby Seyfert galaxies. All these pieces of information suggest that part of observed ionised outflow in NGC 2110 might be driven by the radio jet. However, the radio jet was bent at radial distances of 200 pc (in projection) from the AGN, and beyond there, most of the gas in the galaxy disc is rotating. △ Less

Submitted 10 May, 2023; originally announced May 2023.

Comments: 25 pages, 22 figures, 3 tables. Accepted for publication in A&A

The James Webb Space Telescope Mission

Authors: Jonathan P. Gardner, John C. Mather, Randy Abbott, James S. Abell, Mark Abernathy, Faith E. Abney, John G. Abraham, Roberto Abraham, Yasin M. Abul-Huda, Scott Acton, Cynthia K. Adams, Evan Adams, David S. Adler, Maarten Adriaensen, Jonathan Albert Aguilar, Mansoor Ahmed, Nasif S. Ahmed, Tanjira Ahmed, Rüdeger Albat, Loïc Albert, Stacey Alberts, David Aldridge, Mary Marsha Allen, Shaune S. Allen, Martin Altenburg , et al. (983 additional authors not shown)

Abstract: Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astrono… ▽ More Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 years, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: Accepted by PASP for the special issue on The James Webb Space Telescope Overview, 29 pages, 4 figures

Absence of nuclear polycyclic aromatic hydrocarbon emission from a compact starburst: The case of the type-2 quasar Mrk 477

Authors: C. Ramos Almeida, D. Esparza-Arredondo, O. Gonzalez-Martin, I. Garcia-Bernete, M. Pereira-Santaella, A. Alonso-Herrero, J. A. Acosta-Pulido, P. S. Bessiere, N. A. Levenson, C. N. Tadhunter, D. Rigopoulou, M. Martinez-Paredes, S. Cazzoli, B. Garcia-Lorenzo

Abstract: Mrk 477 is the closest type-2 quasar (QSO2), at a distance of 163 Mpc. This makes it an ideal laboratory for studying the interplay between nuclear activity and star formation with a great level of detail and signal-to-noise. In this Letter we present new mid-infrared (mid-IR) imaging and spectroscopic data with an angular resolution of 0.4 arcsec (~300 pc) obtained with the Gran Telescopio Canari… ▽ More Mrk 477 is the closest type-2 quasar (QSO2), at a distance of 163 Mpc. This makes it an ideal laboratory for studying the interplay between nuclear activity and star formation with a great level of detail and signal-to-noise. In this Letter we present new mid-infrared (mid-IR) imaging and spectroscopic data with an angular resolution of 0.4 arcsec (~300 pc) obtained with the Gran Telescopio Canarias (GTC) instrument CanariCam. The N-band (8-13 micron) spectrum of the central ~400 pc of the galaxy reveals [S IV]10.51 micron emission, but no 8.6 or 11.3 micron PAH features, which are commonly used as tracers of recent star formation. This is in stark contrast with the presence of a nuclear starburst of ~300 pc in size, an age of 6 Myr, and a mass of 1.1x10^8 Msun, as constrained from ultraviolet Hubble Space Telescope observations. Considering this, we argue that even the more resilient, neutral molecules that mainly produce the 11.3 micron PAH band are most likely being destroyed in the vicinity of the active nucleus despite the relatively large X-ray column density, of log N_H=23.5 cm^-2, and modest X-ray luminosity, of 1.5x10^43 erg/s. This highlights the importance of being cautious when using PAH features as star formation tracers in the central region of galaxies to evaluate the impact of feedback from active galactic nuclei. △ Less

Submitted 3 January, 2023; v1 submitted 2 December, 2022; originally announced December 2022.

Comments: 7 pages, 3 figures. A&A Letters accepted version including proof corrections

Journal ref: A&A 669, L5 (2023)

The Science Performance of JWST as Characterized in Commissioning

Authors: Jane Rigby, Marshall Perrin, Michael McElwain, Randy Kimble, Scott Friedman, Matt Lallo, René Doyon, Lee Feinberg, Pierre Ferruit, Alistair Glasse, Marcia Rieke, George Rieke, Gillian Wright, Chris Willott, Knicole Colon, Stefanie Milam, Susan Neff, Christopher Stark, Jeff Valenti, Jim Abell, Faith Abney, Yasin Abul-Huda, D. Scott Acton, Evan Adams, David Adler , et al. (601 additional authors not shown)

Abstract: This paper characterizes the actual science performance of the James Webb Space Telescope (JWST), as determined from the six month commissioning period. We summarize the performance of the spacecraft, telescope, science instruments, and ground system, with an emphasis on differences from pre-launch expectations. Commissioning has made clear that JWST is fully capable of achieving the discoveries f… ▽ More This paper characterizes the actual science performance of the James Webb Space Telescope (JWST), as determined from the six month commissioning period. We summarize the performance of the spacecraft, telescope, science instruments, and ground system, with an emphasis on differences from pre-launch expectations. Commissioning has made clear that JWST is fully capable of achieving the discoveries for which it was built. Moreover, almost across the board, the science performance of JWST is better than expected; in most cases, JWST will go deeper faster than expected. The telescope and instrument suite have demonstrated the sensitivity, stability, image quality, and spectral range that are necessary to transform our understanding of the cosmos through observations spanning from near-earth asteroids to the most distant galaxies. △ Less

Submitted 10 April, 2023; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: 5th version as accepted to PASP; 31 pages, 18 figures; https://iopscience.iop.org/article/10.1088/1538-3873/acb293

Journal ref: PASP 135 048001 (2023)

Hypercubes of AGN Tori (HYPERCAT) -- II. Resolving the Torus with Extremely Large Telescopes

Authors: R. Nikutta, E. Lopez-Rodriguez, K. Ichikawa, N. A. Levenson, C. Packham, S. F. Hönig, A. Alonso-Herrero

Abstract: Recent infrared interferometric observations revealed sub-parsec scale dust distributions around active galactic nuclei (AGNs). Using images of CLUMPY torus models and NGC 1068 as an example, we demonstrate that the near- and mid-infrared nuclear emission of some nearby AGNs will be resolvable in direct imaging with the next generation of 30~m telescopes, potentially breaking degeneracies from pre… ▽ More Recent infrared interferometric observations revealed sub-parsec scale dust distributions around active galactic nuclei (AGNs). Using images of CLUMPY torus models and NGC 1068 as an example, we demonstrate that the near- and mid-infrared nuclear emission of some nearby AGNs will be resolvable in direct imaging with the next generation of 30~m telescopes, potentially breaking degeneracies from previous studies that used integrated spectral energy distributions of unresolved AGN tori. To that effect we model wavelength-dependent point spread functions from the pupil images of various telescopes: James Webb Space Telescope, Keck, Giant Magellan Telescope, Thirty Meter Telescope, and Extremely Large Telescope. We take into account detector pixel scales and noise, and apply deconvolution techniques for image recovery. We also model 2D maps of the 10-micron silicate feature strength, $S_{10}$, of NGC 1068 and compare with observations. When the torus is resolved, we find $S_{10}$ variations across the image. However, to reproduce the $S_{10}$ measurements of an unresolved torus a dusty screen of $A_V > 9$ mag is required. We also fit the first resolved image of the K-band emission in NGC 1068 recently published by the GRAVITY collaboration, deriving likely model parameters of the underlying dust distribution. We find that both 1) an elongated structure suggestive of a highly inclined emission ring, and 2) a geometrically thin but optically thick flared disk where the emission arises from a narrow strip of hot cloud surface layers on the far inner side of the torus funnel, can explain the observations. △ Less

Submitted 24 September, 2021; originally announced September 2021.

Comments: 21 pages, 11 figures, accepted for publication in ApJ

Hypercubes of AGN Tori (HYPERCAT) -- I. Models and Image Morphology

Authors: R. Nikutta, E. Lopez-Rodriguez, K. Ichikawa, N. A. Levenson, C. Packham, S. F. Hönig, A. Alonso-Herrero

Abstract: Near- and mid-infrared interferometers have resolved the dusty parsec-scale obscurer (torus) around nearby active galactic nuclei (AGNs). With the arrival of extremely large single-aperture telescopes, the emission morphology will soon be resolvable unambiguously, without modeling directly the underlying brightness distribution probed by interferometers today. Simulations must instead deliver the… ▽ More Near- and mid-infrared interferometers have resolved the dusty parsec-scale obscurer (torus) around nearby active galactic nuclei (AGNs). With the arrival of extremely large single-aperture telescopes, the emission morphology will soon be resolvable unambiguously, without modeling directly the underlying brightness distribution probed by interferometers today. Simulations must instead deliver the projected 2D brightness distribution as a result of radiative transfer through a 3D distribution of dusty matter around the AGN. We employ such physically motivated 3D dust distributions in tori around AGNs to compute 2D images of the emergent thermal emission using CLUMPY, a dust radiative transfer code for clumpy media. We demonstrate that CLUMPY models can exhibit morphologies with significant polar elongation in the mid-infrared (i.e. the emission extends perpendicular to the dust distribution) on scales of several parsecs, in line with observations in several nearby AGNs. We characterize the emission and cloud distribution morphologies. The observed emission from near- to mid-infrared wavelengths generally does not trace the bulk of the cloud distribution. The elongation of the emission is sensitive to the torus opening angle or scale height. For cloud distributions with a flat radial profile, polar extended emission is realized only at wavelengths shorter than ~18 micron, and shorter than ~5 micron for steep profiles. We make the full results available through HYPERCAT, a large hypercube of resolved AGN torus brightness maps computed with CLUMPY. HYPERCAT also comprises software to process and analyze such large data cubes and provides tools to simulate observations with various current and future telescopes. △ Less

Submitted 24 September, 2021; originally announced September 2021.

Comments: 28 pages, 18 figures, accepted for publication in ApJ

The Galaxy Activity, Torus and Outflow Survey (GATOS): II. Torus and polar dust emission in nearby Seyfert galaxies

Authors: A. Alonso-Herrero, S. García-Burillo, S. F. Hoenig, I. García-Bernete, C. Ramos Almeida, O. González-Martín, E. López-Rodríguez, P. G. Boorman, A. J. Bunker, L. Burtscher, F. Combes, R. Davies, T. Díaz-Santos, P. Gandhi, B. García-Lorenzo, E. K. S. Hicks, L. K. Hunt, K. Ichikawa, M. Imanishi, T. Izumi, A. Labiano, N. A. Levenson, C. Packham, M. Pereira-Santaella, C. Ricci , et al. (8 additional authors not shown)

Abstract: We compare mid-IR and ALMA far-IR images of 12 nearby Seyferts selected from GATOS. The mid-IR unresolved emission contributes more than 60% of the nuclear emission in most galaxies. By contrast, the ALMA 870micron continuum emission is mostly resolved and typically along the torus equatorial direction (Paper I, Garcia-Burillo et al. 2021). The Eddington ratios and nuclear hydrogen column densitie… ▽ More We compare mid-IR and ALMA far-IR images of 12 nearby Seyferts selected from GATOS. The mid-IR unresolved emission contributes more than 60% of the nuclear emission in most galaxies. By contrast, the ALMA 870micron continuum emission is mostly resolved and typically along the torus equatorial direction (Paper I, Garcia-Burillo et al. 2021). The Eddington ratios and nuclear hydrogen column densities NH of half the sample are favorable to launching polar and/or equatorial dusty winds, according to simulations. Six show mid-IR extended emission in the polar direction as traced by the NLR and perpendicular to the ALMA emission. In a few, the nuclear NH might be too high to uplift large quantities of dusty material along the polar direction. Five galaxies have low NH and/or Eddington ratios and thus polar dusty winds are not likely. We generate new CAT3D-WIND disk-wind model images. At low wind-to-disk cloud ratios the far-IR model images have disk- and ring-like morphologies. The X-shape associated with dusty winds is seen better in the far-IR at intermediate-high inclinations for the extended-wind configurations. In most models, the mid-IR emission comes from the inner part of the disk/cone. Extended bi-conical and one-sided polar mid-IR emission is seen in extended-wind configurations and high wind-to-disk cloud ratios. When convolved to our resolution, the model images reproduce qualitative aspects of the observed morphologies. Low-intermediate wind-to-disk ratios are required to account for the large fractions of unresolved mid-IR emission. This work and Paper I provide observational support for the torus+wind scenario. The wind component is more relevant at high Eddington ratios and/or AGN luminosities, and polar dust emission is predicted at NH of up to $10^{24}$cm$^{-2}$. The torus/disk component, on the other hand, prevails at low luminosities and/or Eddington ratios. (Abridged) △ Less

Submitted 2 July, 2021; v1 submitted 1 July, 2021; originally announced July 2021.

Comments: 27 pages, 27 figures. Accepted for publication in A&A

The Galaxy Activity, Torus and Outflow Survey (GATOS) I. ALMA images of dusty molecular tori in Seyfert galaxies

Authors: S. García-Burillo, A. Alonso-Herrero, C. Ramos Almeida, O. González-Martín, F. Combes, A. Usero, S. Hönig, M. Querejeta, E. K. S. Hicks, L. K. Hunt, D. Rosario, R. Davies, P. G. Boorman, A. J. Bunker, L. Burtscher, L. Colina, T. Díaz-Santos, P. Gandhi, I. García-Bernete, B. García-Lorenzo, K. Ichikawa, M. Imanishi, T. Izumi, A. Labiano, N. A. Levenson , et al. (10 additional authors not shown)

Abstract: We present the first results of the Galaxy Activity, Torus and Outflow Survey (GATOS), a project aimed at understanding the properties of the dusty molecular tori and their connection to the host galaxy in nearby Seyfert galaxies. Our project expands the range of AGN luminosities and Eddington ratios covered by previous surveys of Seyferts conducted by ALMA and allows us to study the gas feeding a… ▽ More We present the first results of the Galaxy Activity, Torus and Outflow Survey (GATOS), a project aimed at understanding the properties of the dusty molecular tori and their connection to the host galaxy in nearby Seyfert galaxies. Our project expands the range of AGN luminosities and Eddington ratios covered by previous surveys of Seyferts conducted by ALMA and allows us to study the gas feeding and feedback cycle in a combined sample of 19 Seyferts. We used ALMA to obtain new images of the emission of molecular gas and dust using the CO(3-2) and HCO+(4-3) lines as well as their underlying continuum emission at 870 microns with high spatial resolutions (0.1'' ~ 7 - 13 pc) in the CND of 10 nearby (D < 28 Mpc) Seyfert galaxies. Our new ALMA observations detect 870 micron continuum and CO line emission from spatially resolved disks located around the AGN in all the sources. The bulk of the continuum flux can be accounted for by thermal emission from dust in the majority of the targets. For most of the sources the disks show a preponderant orientation perpendicular to the AGN wind axes, as expected for dusty molecular tori. The median diameters and molecular gas masses of the tori are ~ 42 pc, and ~ 6 x 10**5 Msun, respectively. We find a positive correlation between the line-of-sight gas column densities responsible for the absorption of X-rays and the molecular gas column densities derived from CO towards the AGN in our sources. The radial distributions of molecular gas in the CND of our combined sample show signs of nuclear-scale molecular gas deficits. We also detect molecular outflows in the sources that show the most extreme nuclear-scale gas deficits in our sample. These observations find for the first time supporting evidence that the imprint of AGN feedback is more extreme in higher luminosity and/or higher Eddington ratio Seyfert galaxies. △ Less

Submitted 22 June, 2021; v1 submitted 20 April, 2021; originally announced April 2021.

Comments: 41 pages, new version revised and accepted by A&A on the 22nd of June 2021

Journal ref: A&A 652, A98 (2021)

Cold molecular gas and PAH emission in Seyfert galaxies

Authors: A. Alonso-Herrero, M. Pereira-Santaella, D. Rigopoulou, I. Garcia-Bernete, S. Garcia-Burillo, A. J. Dominguez-Fernandez, F. Combes, R. I. Davies, T. Diaz-Santos, D. Esparza-Arredondo, O. Gonzalez-Martin, A. Hernan-Caballero, E. K. S. Hicks, S. F. Hoenig, N. A. Levenson, C. Ramos Almeida, P. F. Roche, D. Rosario

Abstract: We investigate the relation between the detection of the $11.3\,μ$m PAH feature in the nuclear ($\sim 24-230\,$pc) regions of 22 nearby Seyfert galaxies and the properties of the cold molecular gas. For the former we use ground-based (0.3-0.6" resolution) mid-infrared (mid-IR) spectroscopy. The cold molecular gas is traced by ALMA and NOEMA high (0.2-1.1") angular resolution observations of the CO… ▽ More We investigate the relation between the detection of the $11.3\,μ$m PAH feature in the nuclear ($\sim 24-230\,$pc) regions of 22 nearby Seyfert galaxies and the properties of the cold molecular gas. For the former we use ground-based (0.3-0.6" resolution) mid-infrared (mid-IR) spectroscopy. The cold molecular gas is traced by ALMA and NOEMA high (0.2-1.1") angular resolution observations of the CO(2-1) transition. Galaxies with a nuclear detection of the $11.3\,μ$m PAH feature contain more cold molecular gas (median $1.6\times 10^7\,M_\odot$) and have higher column densities ($N({\rm H}_2) = 2 \times 10^{23}\,{\rm cm}^{-2}$) over the regions sampled by the mid-IR slits than those without a detection. This suggests that molecular gas plays a role in shielding the PAH molecules in the harsh environments of Seyfert nuclei. Choosing the PAH molecule naphthalene as an illustration, we compute its half-life in the nuclear regions of our sample when exposed to 2.5keV hard X-ray photons. We estimate shorter half-lives for naphthalene in nuclei without a $11.3\,μ$m PAH detection than in those with a detection. The Spitzer/IRS PAH ratios on circumnuclear scales ($\sim$ 4" $\sim$ 0.25-1.3kpc) are in between model predictions for neutral and partly ionized PAHs. However, Seyfert galaxies in our sample with the highest nuclear H$_2$ column densities are not generally closer to the neutral PAH tracks. This is because in the majority of our sample galaxies, the CO(2-1) emission in the inner $\sim$ 4" is not centrally peaked and in some galaxies traces circumnuclear sites of strong star formation activity. Spatially resolved observations with the MIRI medium-resolution spectrograph (MRS) on the James Webb Space Telescope will be able to distinguish the effects of an active galactic nucleus (AGN) and star formation on the PAH emission in nearby AGN. △ Less

Submitted 25 May, 2020; originally announced May 2020.

Comments: Accepted for publication in A&A

Searching for molecular gas in/outflows in the nuclear regions of five Seyfert galaxies

Authors: A. J. Domínguez-Fernández, A. Alonso-Herrero, S. García-Burillo, R. I. Davies, A. Usero, A. Labiano, N. A. Levenson, M. Pereira-Santaella, M. Imanishi, C. Ramos Almeida, D. Rigopoulou

Abstract: AGN-driven outflows are believed to play an important role in regulating the growth of galaxies mostly via negative feedback. However, their effects on their hosts are far from clear, especially for low and moderate luminosity Seyferts. To investigate this issue, we have obtained cold molecular gas observations, traced by the CO(2-1) transition, using the NOEMA interferometer of five nearby (dista… ▽ More AGN-driven outflows are believed to play an important role in regulating the growth of galaxies mostly via negative feedback. However, their effects on their hosts are far from clear, especially for low and moderate luminosity Seyferts. To investigate this issue, we have obtained cold molecular gas observations, traced by the CO(2-1) transition, using the NOEMA interferometer of five nearby (distances between 19 and 58 Mpc) Seyfert galaxies. The resolution of approx. 0.3-0.8 arcsec (approx. 30-100 pc) and field of view of NOEMA allowed us to study the CO(2-1) morphology and kinematics in the nuclear regions (approx. 100 pc) and up to radial distances of approx. 900 pc. We have detected CO(2-1) emission in all five galaxies with disky or circumnuclear ring like morphologies. We derived cold molecular gas masses on nuclear (approx. 100 pc) and circumnuclear (approx. 650 pc) scales in the range from $10^6$ to $10^7$M$_{\odot}$ and from $10^7$ to $10^8$ $M_{\odot}$, respectively. In all of our galaxies the bulk of this gas is rotating in the plane of the galaxy. However, non-circular motions are also present. In NGC 4253, NGC 4388 and NGC 7465, we can ascribe the streaming motions to the presence of a large-scale bar. In Mrk 1066 and NGC 4388, the non-circular motions in the nuclear regions are explained as outflowing material due to the interaction of the AGN wind with molecular gas in the galaxy disk. We conclude that for an unambiguous and precise interpretation of the kinematics of the cold molecular gas we need a detailed knowledge of the host galaxy (i.e., presence of bars, interactions, etc) as well as of the ionized gas kinematics and the ionization cone geometry. △ Less

Submitted 12 March, 2020; originally announced March 2020.

Comments: Accepted for publication in A&A, 23 pages, 17 figures

NuSTAR Uncovers an Extremely Local Compton-thick AGN in NGC 4968

Authors: Stephanie M. LaMassa, Tahir Yaqoob, Peter G. Boorman, Panayiotis Tzanavaris, N. A. Levenson, Poshak Gandhi, Andrew F. Ptak, Timothy M. Heckman

Abstract: We present the analysis of Chandra and NuSTAR spectra of NGC 4968, a local (D$\sim$44 Mpc) 12$μ$m-selected Seyfert 2 galaxy, enshrouded within Compton-thick layers of obscuring gas. We find no evidence of variability between the Chandra and NuSTAR observations (separated by 2 years), and between the two NuSTAR observations (separated by 10 months). Using self-consistent X-ray models, we rule out t… ▽ More We present the analysis of Chandra and NuSTAR spectra of NGC 4968, a local (D$\sim$44 Mpc) 12$μ$m-selected Seyfert 2 galaxy, enshrouded within Compton-thick layers of obscuring gas. We find no evidence of variability between the Chandra and NuSTAR observations (separated by 2 years), and between the two NuSTAR observations (separated by 10 months). Using self-consistent X-ray models, we rule out the scenario where the obscuring medium is nearly spherical and uniform, contradicting the results implied by the $<$10 keV Chandra spectrum. The line-of-sight column density, from intervening matter between the source and observer that intercepts the intrinsic AGN X-ray emission, is well within the Compton-thick regime, with a minimum column density of $2\times10^{24}$ cm$^{-2}$. The average global column density is high ($> 3\times10^{23}$ cm$^{-2}$), with both Compton-thick and Compton-thin solutions permitted depending on the X-ray spectral model. The spectral models provide a range of intrinsic AGN continuum parameters and implied 2-10 keV luminosities ($L_{\rm 2-10keV,intrinsic}$), where the higher end of $L_{\rm 2-10keV,intrinsic}$ is consistent with expectations from the 12$μ$m luminosity ($L_{\rm 2-10keV,intrinisc} \sim 7\times10^{42}$ erg s$^{-1}$). Compared with Compton-thick AGN previously observed by {\it NuSTAR}, NGC 4968 is among the most intrinsically X-ray luminous. However, despite its close proximity and relatively high intrinsic X-ray luminosity, it is undetected by the 105 month Swift-BAT survey, underscoring the importance of multi-wavelength selection for obtaining the most complete census of the most hidden black holes. △ Less

Submitted 13 November, 2019; originally announced November 2019.

Comments: accepted for publication in ApJ; 20 pages, 8 figures

A nuclear molecular outflow in the Seyfert galaxy NGC3227

Authors: A. Alonso-Herrero, S. Garcia-Burillo, M. Pereira-Santaella, R. I. Davies, F. Combes, M. Vestergaard, S. I. Raimundo, A. Bunker, T. Diaz-Santos, P. Gandhi, I. Garcia-Bernete, E. K. S. Hicks, S. F. Hönig, L. K. Hunt, M. Imanishi, T. Izumi, N. A. Levenson, W. Maciejewski, C. Packham, C. Ramos Almeida, C. Ricci, D. Rigopoulou, P. F. Roche, D. Rosario, M. Schartmann , et al. (2 additional authors not shown)

Abstract: We present ALMA observations of the CO(2-1) and CO(3-2) molecular gas transitions and associated (sub)-mm continua of the nearby Seyfert 1.5 galaxy NGC3227 with angular resolutions 0.085-0.21" (7-15pc). On large scales the cold molecular gas shows circular motions as well as streaming motions on scales of a few hundred parsecs associated with a large scale bar. We fitted the nuclear ALMA 1.3mm emi… ▽ More We present ALMA observations of the CO(2-1) and CO(3-2) molecular gas transitions and associated (sub)-mm continua of the nearby Seyfert 1.5 galaxy NGC3227 with angular resolutions 0.085-0.21" (7-15pc). On large scales the cold molecular gas shows circular motions as well as streaming motions on scales of a few hundred parsecs associated with a large scale bar. We fitted the nuclear ALMA 1.3mm emission with an unresolved component and an extended component. The 850$μ$m emission shows at least two extended components, one along the major axis of the nuclear disk and the other along the axis of the ionization cone. The molecular gas in the central region (1" ~73pc) shows several CO clumps with complex kinematics which appears to be dominated by non-circular motions. While we cannot demonstrate conclusively the presence of a warped nuclear disk, we also detected non-circular motions along the kinematic minor axis. They reach line-of-sight velocities of v-vsys =150-200km/s. Assuming that the radial motions are in the plane of the galaxy, then we interpret them as a nuclear molecular outflow due to molecular gas in the host galaxy being entrained by the AGN wind. We derive molecular outflow rates of $5\,M_\odot\,{\rm yr}^{-1}$ and $0.6\,M_\odot\,{\rm yr}^{-1}$ at projected distances of up to 30pc to the northeast and southwest of the AGN, respectively. At the AGN location we estimate a mass in molecular gas of $5\times 10^{5}\,M_\odot$ and an average column density $N({\rm H}_2) = 2-3\times 10^{23}\,{\rm cm}^{-2}$ in the inner 15pc. The nuclear molecular gas and sub-mm continuum emission of NGC3227 do not resemble the classical compact torus. Rather, these emissions extend for several tens of parsecs and appear connected with the circumnuclear ring in the host galaxy disk, as found in other local AGN. (Abridged) △ Less

Submitted 17 June, 2019; originally announced June 2019.

Comments: Accepted for publication in Astronomy and Astrophysics, 17 pages, 12 figures

ALMA polarimetry measures magnetically aligned dust grains in the torus of NGC 1068

Authors: Enrique Lopez-Rodriguez, Almudena Alonso-Herrero, Santiago García-Burillo, Michael S. Gordon, Kohei Ichikawa, Masatoshi Imanishi, Seiji Kameno, Nancy A. Levenson, Robert Nikutta, Chris Packham

Abstract: The obscuring structure surrounding active galactic nuclei (AGN) can be explained as a dust and gas flow cycle that fundamentally connects the AGN with their host galaxies. This structure is believed to be associated with dusty winds driven by radiation pressure. However, the role of magnetic fields, which are invoked in almost all models for accretion onto a supermassive black hole and outflows,… ▽ More The obscuring structure surrounding active galactic nuclei (AGN) can be explained as a dust and gas flow cycle that fundamentally connects the AGN with their host galaxies. This structure is believed to be associated with dusty winds driven by radiation pressure. However, the role of magnetic fields, which are invoked in almost all models for accretion onto a supermassive black hole and outflows, is not thoroughly studied. Here we report the first detection of polarized thermal emission by means of magnetically aligned dust grains in the dusty torus of NGC 1068 using ALMA Cycle 4 polarimetric dust continuum observations ($0.07"$, $4.2$ pc; 348.5 GHz, $860$ $μ$m). The polarized torus has an asymmetric variation across the equatorial axis with a peak polarization of $3.7\pm0.5$\% and position angle of $109\pm2^{\circ}$ (B-vector) at $\sim8$ pc east from the core. We compute synthetic polarimetric observations of magnetically aligned dust grains assuming a toroidal magnetic field and homogeneous grain alignment. We conclude that the measured 860 $μ$m continuum polarization arises from magnetically aligned dust grains in an optically thin region of the torus. The asymmetric polarization across the equatorial axis of the torus arises from 1) an inhomogeneous optical depth, and 2) a variation of the velocity dispersion, i.e. variation of the magnetic field turbulence at sub-pc scales, from the eastern to the western region of the torus. These observations and modeling constrain the torus properties beyond spectral energy distribution results. This study strongly supports that magnetic fields up to a few pc contribute to the accretion flow onto the active nuclei. △ Less

Submitted 16 March, 2020; v1 submitted 21 May, 2019; originally announced May 2019.

Comments: 19 pages, 11 figures (Accepted for Publication to ApJ)

Torus model properties of an ultra-hard X-ray selected sample of Seyfert galaxies

Authors: I. García-Bernete, C. Ramos Almeida, A. Alonso-Herrero, M. J. Ward, J. A. Acosta-Pulido, M. Pereira-Santaella, A. Hernán-Caballero, A. Asensio Ramos, O. González-Martín, N. A. Levenson, S. Mateos, F. J. Carrera, C. Ricci, P. Roche, I. Marquez, C. Packham, J. Masegosa, L. Fuller

Abstract: We characterize for the first time the torus properties of an ultra-hard X-ray (14-195 keV) volume-limited (DL<40 Mpc) sample of 24 Seyfert (Sy) galaxies (BCS40 sample). The sample was selected from the Swift/BAT nine month catalog. We use high angular resolution nuclear infrared (IR) photometry and N-band spectroscopy, the CLUMPY torus models and a Bayesian tool to characterize the properties of… ▽ More We characterize for the first time the torus properties of an ultra-hard X-ray (14-195 keV) volume-limited (DL<40 Mpc) sample of 24 Seyfert (Sy) galaxies (BCS40 sample). The sample was selected from the Swift/BAT nine month catalog. We use high angular resolution nuclear infrared (IR) photometry and N-band spectroscopy, the CLUMPY torus models and a Bayesian tool to characterize the properties of the nuclear dust. In the case of the Sy1s we estimate the accretion disk contribution to the subarcsecond resolution nuclear IR SEDs (~0.4'') which is, on average, 46+-28, 23+-13 and 11+-5% in the J-, H- and K-bands, respectively. This indicates that the accretion disk templates that assume a steep fall for longer wavelengths than 1 micron might underestimate its contribution to the near-IR emission. Using both optical (broad vs narrow lines) and X-ray (unabsorbed vs absorbed) classifications, we compare the global posterior distribution of the torus model parameters. We confirm that Sy2s have larger values of the torus covering factor (CT~0.95) than Sy1s (CT~0.65) in our volume-limited Seyfert sample. These findings are independent of whether we use an optical or X-ray classification. We find that the torus covering factor remains essentially constant within the errors in our luminosity range and there is no clear dependence with the Eddington ratio. Finally, we find tentative evidence that even an ultra hard X-ray selection is missing a significant fraction of highly absorbed type 2 sources with very high covering factor tori. △ Less

Submitted 7 April, 2019; originally announced April 2019.

Comments: 28 pages, 14 figures and 14 tables. Monthly Notices of the Royal Astronomical Society, Accepted 2019 April 4. Received 2019 March 14; in original form 2019 January 18, in press

Quantifying star formation activity in the inner 1kpc of local MIR bright QSOs

Authors: M. Martínez-Paredes, I. Aretxaga, O. Gonzalez-Martin, A. Alonso-Herrero, N. A. Levenson, C. Ramos Almeida, E. Lopez-Rodriguez

Abstract: We examine star formation activity in a distance- (z<0.1) and flux-limited sample of quasars (QSOs). Mid-infrared (MIR) spectral diagnostics at high spatial resolution (~0.4 arcsec) yield star formation rates (SFRs) in the inner regions (~300 pc to 1 kpc) for 13 of 20 of the sample members. We group these objects according to the size probed by the high angular resolution spectroscopy, with charac… ▽ More We examine star formation activity in a distance- (z<0.1) and flux-limited sample of quasars (QSOs). Mid-infrared (MIR) spectral diagnostics at high spatial resolution (~0.4 arcsec) yield star formation rates (SFRs) in the inner regions (~300 pc to 1 kpc) for 13 of 20 of the sample members. We group these objects according to the size probed by the high angular resolution spectroscopy, with characteristic scales of <0.7 and ~0.7-1 kpc. Using the polycyclic aromatic hydrocarbon (PAH) feature at 11.3 um, we measure SFRs around 0.2 and 1.6 M_{\odot}yr^{-1}. We also measure the larger aperture PAH-derived SFRs in the individual IRS/Spitzer spectra of the sample and obtain a clear detection in ~58 percent of them. We compare smaller and larger aperture measurements and find that they are similar, suggesting that star formation activity in these QSOs is more centrally concentrated, with the inner region (>~1 kpc) accounting for the majority of star formation measured on these scales, and that PAH molecules can be present in most local MIR-bright QSOs within a few hundred pc from the central engine. By comparison with merger simulations, we find that our estimation of the SFR and black hole (BH) accretion rates are consistent with a scenario in which the star formation activity is centrally peaked as predicted by simulations. △ Less

Submitted 8 January, 2019; originally announced January 2019.

Comments: Accepted for publication in ApJ

SOFIA/FORCAST resolves 30 - 40 um extended dust emission in nearby active galactic nuclei

Authors: Lindsay Fuller, Enrique Lopez-Rodriguez, Chris Packham, Kohei Ichikawa, Aditya Togi, Almudena Alonso-Herrero, Cristina Ramos-Almeida, Tanio Diaz-Santos, N. A. Levenson, James Radomski

Abstract: We present arcsecond-scale observations of the active galactic nuclei (AGNs) of seven nearby Seyfert galaxies observed from the Stratospheric Observatory For Infrared Astronomy (SOFIA) using the 31.5 and 37.1 um filters of the Faint Object infraRed CAmera for the SOFIA Telescope (FORCAST). We isolate unresolved emission from the torus and find extended diffuse emission in six 37.1 um residual imag… ▽ More We present arcsecond-scale observations of the active galactic nuclei (AGNs) of seven nearby Seyfert galaxies observed from the Stratospheric Observatory For Infrared Astronomy (SOFIA) using the 31.5 and 37.1 um filters of the Faint Object infraRed CAmera for the SOFIA Telescope (FORCAST). We isolate unresolved emission from the torus and find extended diffuse emission in six 37.1 um residual images in our sample. Using Spitzer/IRS spectra, we determine the dominant mid-infrared (MIR) extended emission source and attribute it to dust in the narrow line region (NLR) or star formation. We compare the optical NLR and radio jet axes to the extended 37.1 um emission and find coincident axes for three sources. For those AGNs with extended emission coincident with the optical axis, we find that spatial scales of the residual images are consistent with 0.1 - 1 kpc scale distances to which dust can be heated by the AGN. Using previously published subarcsecond 1 - 20 um imaging and spectroscopic data along with our new observations, we construct broadband spectral energy distributions (SEDs) of the AGNs at wavelengths 1 - 40 um. We find that three AGNs in our sample tentatively show a turnover in the SED between 30 - 40 um. Using results from Clumpy torus models and the Bayesian inference tool BayesClumpy, we find that the posterior outputs for AGNs with MIR turnover revealed by SOFIA/FORCAST have smaller uncertainties than AGNs that do not show a turnover. △ Less

Submitted 19 December, 2018; originally announced December 2018.

Comments: 16 pages, 9 figures

The origin of the mid-infrared nuclear polarization of active galactic nuclei

Authors: E. Lopez-Rodriguez, A. Alonso-Herrero, T. Diaz-Santos, O. Gonzalez-Martin, K. Ichikawa, N. A. Levenson, M. Martinez-Paredes, R. Nikutta, C. Packham, E. Perlman, C. Ramos Almeida, J. M. Rodriguez-Espinosa, C. M. Telesco

Abstract: We combine new (NGC 1275, NGC 4151, and NGC 5506) and previously published (Cygnus A, Mrk 231, and NGC 1068) sub-arcsecond resolution mid-infrared (MIR; 8-13 $μ$m) imaging- and spectro-polarimetric observations of six Seyfert galaxies using CanariCam on the 10.4-m Gran Telescopio CANARIAS. These observations reveal a diverse set of physical processes responsible for the nuclear polarization, and p… ▽ More We combine new (NGC 1275, NGC 4151, and NGC 5506) and previously published (Cygnus A, Mrk 231, and NGC 1068) sub-arcsecond resolution mid-infrared (MIR; 8-13 $μ$m) imaging- and spectro-polarimetric observations of six Seyfert galaxies using CanariCam on the 10.4-m Gran Telescopio CANARIAS. These observations reveal a diverse set of physical processes responsible for the nuclear polarization, and permit characterization of the origin of the MIR nuclear polarimetric signature of active galactic nuclei (AGN). For all radio quiet objects, we found that the nuclear polarization is low (<1 per cent), and the degree of polarization is often a few per cent over extended regions of the host galaxy where we have sensitivity to detect such extended emission (i.e., NGC 1068 and NGC 4151). We suggest that the higher degree of polarization previously found in lower resolution data arises only on the larger-than-nuclear scales. Only the radio-loud Cygnus A exhibits significant nuclear polarization ($\sim$11 per cent), attributable to synchrotron emission from the pc-scale jet close to the core. We present polarization models that suggest that the MIR nuclear polarization for highly obscured objects arises from a self-absorbed MIR polarized clumpy torus and/or dichroism from the host galaxy, while for unabsorbed cores, MIR polarization arises from dust scattering in the torus and/or surrounding nuclear dust. △ Less

Submitted 4 May, 2018; originally announced May 2018.

Comments: 9 pages, 4 figures. Accepted to MNRAS

Resolving the nuclear obscuring disk in the Compton-thick Seyfert galaxy NGC5643 with ALMA

Authors: A. Alonso-Herrero, M. Pereira-Santaella, S. García-Burillo, R. I. Davies, F. Combes, D. Asmus, A. Bunker, T. Díaz-Santos, P. Gandhi, O. González-Martín, A. Hernán-Caballero, E. Hicks, S. Hönig, A. Labiano, N. A. Levenson, C. Packham, C. Ramos Almeida, C. Ricci, D. Rigopoulou, D. Rosario, E. Sani, M. J. Ward

Abstract: We present ALMA Band 6 $^{12}$CO(2--1) line and rest-frame 232GHz continuum observations of the nearby Compton-thick Seyfert galaxy NGC5643 with angular resolutions 0.11-0.26arcsec (9-21pc). The CO(2--1) integrated line map reveals emission from the nuclear and circumnuclear region with a two-arm nuclear spiral extending 10arcsec on each side. The circumnuclear CO(2--1) kinematics can be fitted wi… ▽ More We present ALMA Band 6 $^{12}$CO(2--1) line and rest-frame 232GHz continuum observations of the nearby Compton-thick Seyfert galaxy NGC5643 with angular resolutions 0.11-0.26arcsec (9-21pc). The CO(2--1) integrated line map reveals emission from the nuclear and circumnuclear region with a two-arm nuclear spiral extending 10arcsec on each side. The circumnuclear CO(2--1) kinematics can be fitted with a rotating disk, although there are regions with large residual velocities and/or velocity dispersions. The CO(2--1) line profiles of these regions show two different velocity components. One is ascribed to the circular component and the other to the interaction of the AGN outflow, as traced by the [O III]5007AA emission, with molecular gas in the disk a few hundred parsecs from the AGN. On nuclear scales, we detected an inclined CO(2--1) disk (diameter 26 pc, FWHM) oriented almost in a north-south direction. The CO(2--1) nuclear kinematics can be fitted with a rotating disk which appears to be tilted with respect to the large scale disk. There are strong non-circular motions in the central 0.2-0.3 arcsec with velocities of up to 110km/s. In the absence of a nuclear bar, these motions could be explained as radial outflows in the nuclear disk. We estimate a total molecular gas mass for the nuclear disk of $M({\rm H}_2)=1.1\times 10^7\,M_\odot$ and an H$_2$ column density toward the location of the AGN of $N({\rm H}_2)\sim 5 \times 10^{23}\,{\rm cm}^{-2}$, for a standard CO-to-H$_2$ conversion factor. We interpret this nuclear molecular gas disk as the obscuring torus of NGC5643 as well as the collimating structure of the ionization cone. △ Less

Submitted 10 May, 2018; v1 submitted 13 April, 2018; originally announced April 2018.

Comments: 16 papes, 9 figures. Accepted to ApJ. Updated Fig.4

The emission and distribution of dust of the torus of NGC 1068

Authors: Enrique Lopez-Rodriguez, Lindsay Fuller, Almudena Alonso-Herrero, Andreas Efstathiou, Kohei Ichikawa, Nancy A. Levenson, Chris Packham, James Radomski, Cristina Ramos Almeida, Dominic J. Benford, Marc Berthoud, Ryan Hamilton, Doyal Harper, Attila Kovavcs, Fabio P. Santos, J. Staguhn, Terry Herter

Abstract: We present observations of NGC 1068 covering the $19.7-53.0$ $μ$m wavelength range using FORCAST and HAWC+ onboard SOFIA. Using these observations, high-angular resolution infrared (IR) and sub-mm observations, we find an observational turn-over of the torus emission in the $30-40$ $μ$m wavelength range with a characteristic temperature of $70-100$ K. This component is clearly different from the d… ▽ More We present observations of NGC 1068 covering the $19.7-53.0$ $μ$m wavelength range using FORCAST and HAWC+ onboard SOFIA. Using these observations, high-angular resolution infrared (IR) and sub-mm observations, we find an observational turn-over of the torus emission in the $30-40$ $μ$m wavelength range with a characteristic temperature of $70-100$ K. This component is clearly different from the diffuse extended emission in the narrow line and star formation regions at 10-100 $μ$m within the central 700 pc. We compute $2.2-432$ $μ$m 2D images using the best inferred \textsc{clumpy} torus model based on several nuclear spectral energy distribution (SED) coverages. We find that when $1-20$ $μ$m SED is used, the inferred result gives a small torus size ($<4$ pc radius) and a steep radial dust distribution. The computed torus using the $1-432$ $μ$m SED provides comparable torus sizes, $5.1^{+0.4}_{-0.4}$ pc radius, and morphology to the recently resolved 432 $μ$m ALMA observations. This result indicates that the $1-20$ $μ$m wavelength range is not able to probe the full extent of the torus. The characterization of the turn-over emission of the torus using the $30-60$ $μ$m wavelength range is sensitive to the detection of cold dust in the torus. The morphology of the dust emission in our 2D image at 432 $μ$m is spatially coincident with the cloud distribution, while the morphology of the emission in the $1-20$ $μ$m wavelength range shows an elongated morphology perpendicular to the cloud distribution. We find that our 2D \textsc{clumpy} torus image at 12 $μ$m can produce comparable results to those observed using IR interferometry. △ Less

Submitted 11 April, 2018; originally announced April 2018.

Comments: 12 pages, 7 Figures, Accepted for publication in ApJ

A mid-infrared statistical investigation of clumpy torus model predictions

Authors: J. García-González, A. Alonso-Herrero, S. F. Hönig, A. Hernán-Caballero, C. Ramos Almeida, N. A. Levenson, P. F. Roche, O. González-Martín, C. Packham, M. Kishimoto

Abstract: We present new calculations of the CAT3D clumpy torus models, which now include a more physical dust sublimation model as well as AGN anisotropic emission. These new models allow graphite grains to persist at temperatures higher than the silicate dust sublimation temperature. This produces stronger near-infrared emission and bluer mid-infrared (MIR) spectral slopes. We make a statistical compariso… ▽ More We present new calculations of the CAT3D clumpy torus models, which now include a more physical dust sublimation model as well as AGN anisotropic emission. These new models allow graphite grains to persist at temperatures higher than the silicate dust sublimation temperature. This produces stronger near-infrared emission and bluer mid-infrared (MIR) spectral slopes. We make a statistical comparison of the CAT3D model MIR predictions with a compilation of sub-arcsecond resolution ground-based MIR spectroscopy of 52 nearby Seyfert galaxies (median distance of 36 Mpc) and 10 quasars. We focus on the AGN MIR spectral index $α_{MIR}$ and the strength of the 9.7 $μ$m silicate feature $S_{Sil}$. As with other clumpy torus models, the new CAT3D models do not reproduce the Seyfert galaxies with deep silicate absorption ($S_{Sil}<-1$). Excluding those, we conclude that the new CAT3D models are in better agreement with the observed $α_{MIR}$ and $S_{Sil}$ of Seyfert galaxies and quasars. We find that Seyfert 2 are reproduced with models with low photon escape probabilities, while the quasars and the Seyfert 1-1.5 require generally models with higher photon escape probabilities. Quasars and Seyfert 1-1.5 tend to show steeper radial cloud distributions and fewer clouds along an equatorial line-of-sight than Seyfert 2. Introducing AGN anisotropic emission besides the more physical dust sublimation models alleviates the problem of requiring inverted radial cloud distributions (i.e., more clouds towards the outer parts of the torus) to explain the MIR spectral indices of type 2 Seyferts. △ Less

Submitted 23 June, 2017; v1 submitted 22 June, 2017; originally announced June 2017.

Comments: 23 pages, 17 figures, 5 tables, Accepted for publication in MNRAS

Chandra Reveals Heavy Obscuration and Circumnuclear Star Formation in Seyfert 2 Galaxy NGC 4968

Authors: Stephanie M. LaMassa, Tahir Yaqoob, N. A. Levenson, Peter Boorman, Timothy M. Heckman, Poshak Gandhi, Jane R. Rigby, C. Megan Urry, Andrew F. Ptak

Abstract: We present the Chandra imaging and spectral analysis of NGC 4968, a nearby (z = 0.00986) Seyfert 2 galaxy. We discover extended ($\sim$1 kpc) X-ray emission in the soft band (0.5 - 2 keV) that is neither coincident with the narrow line region nor the extended radio emission. Based on spectral modeling, it is linked to on-going star formation ($\sim$2.6-4 M$_{\sun}$ yr$^{-1}$). The soft emission at… ▽ More We present the Chandra imaging and spectral analysis of NGC 4968, a nearby (z = 0.00986) Seyfert 2 galaxy. We discover extended ($\sim$1 kpc) X-ray emission in the soft band (0.5 - 2 keV) that is neither coincident with the narrow line region nor the extended radio emission. Based on spectral modeling, it is linked to on-going star formation ($\sim$2.6-4 M$_{\sun}$ yr$^{-1}$). The soft emission at circumnuclear scales (inner $\sim$400 pc) originates from hot gas, with kT $\sim$ 0.7 keV, while the most extended thermal emission is cooler (kT $\sim$ 0.3 keV). We refine previous measurements of the extreme Fe K$α$ equivalent width in this source (EW = 2.5$^{+2.6}_{-1.0}$ keV), which suggests the central engine is completely embedded within Compton-thick levels of obscuration. Using physically motivated models fit to the Chandra spectrum, we derive a Compton-thick column density ($N_{\rm H} > 1.25\times10^{24}$ cm$^{-2}$) and an intrinsic hard (2-10 keV) X-ray luminosity of $\sim$3-8$\times 10^{42}$ erg s$^{-1}$ (depending on the presumed geometry of the obscurer), which is over two orders of magnitude larger than that observed. The large Fe K$α$ EW suggests a spherical covering geometry, which could be confirmed with X-ray measurements above 10 keV. NGC 4968 is similar to other active galaxies that exhibit extreme Fe K$α$ EWs (i.e., $>$2 keV) in that they also contain ongoing star formation. This work supports the idea that gas associated with nuclear star formation may increase the covering factor of the enshrouding gas and play a role in obscuring AGN. △ Less

Submitted 6 December, 2016; originally announced December 2016.

Comments: 11 pages, 8 figures, 4 tables. Accepted for publication in ApJ

Infrared polarimetry of Mrk 231: Scattering off hot dust grains in the central core

Authors: E. Lopez-Rodriguez, C. Packham, T. J. Jones, R. Siebenmorgen, P. F. Roche, N. A. Levenson, A. Alonso-Herrero, E. Perlman, K. Ichikawa, C. Ramos Almeida, O. Gonzalez-Martin, R. Nikutta, M. Martinez-Paredes, D. Shenoy, M. S. Gordon, C. M. Telesco

Abstract: We present high-angular (0.17$-$0.35 arcsec) resolution imaging polarimetric observations of Mrk 231 in the 3.1 $μ$m filter using MMT-Pol on the 6.5-m MMT, and in the 8.7 $μ$m, 10.3 $μ$m, and 11.6 $μ$m filters using CanariCam on the 10.4-m Gran Telescopio CANARIAS. In combination with already published observations, we compile the 1$-$12 $μ$m total and polarized nuclear spectral energy distributio… ▽ More We present high-angular (0.17$-$0.35 arcsec) resolution imaging polarimetric observations of Mrk 231 in the 3.1 $μ$m filter using MMT-Pol on the 6.5-m MMT, and in the 8.7 $μ$m, 10.3 $μ$m, and 11.6 $μ$m filters using CanariCam on the 10.4-m Gran Telescopio CANARIAS. In combination with already published observations, we compile the 1$-$12 $μ$m total and polarized nuclear spectral energy distribution (SED). The total flux SED in the central 400 pc is explained as the combination of 1) a hot (731 $\pm$ 4 K) dusty structure, directly irradiated by the central engine, which is at 1.6 $\pm$ 0.1 pc away and attributed to be in the pc-scale polar region, 2) an optically-thick, smooth and disk-like dusty structure (`torus') with an inclination of 48 $\pm$ 23$^{\circ}$ surrounding the central engine, and 3) an extinguished (A$_{\mbox{V}} =$ 36 $\pm$ 5 mag) starburst component. The polarized SED decreases from 0.77 $\pm$ 0.14 per cent at 1.2 $μ$m to 0.31 $\pm$ 0.15 per cent at 11.6 $μ$m and follows a power-law function, $λ^{\sim0.57}$. The polarization angle remains constant ($\sim$108$^{\circ}$) in the 1$-$12 $μ$m wavelength range. The dominant polarization mechanism is explained as scattering off hot dust grains in the pc-scale polar regions. △ Less

Submitted 28 September, 2016; originally announced September 2016.

Comments: 9 pages, 2 figures. Accepted for publication in MNRAS

The complex evolutionary paths of local infrared bright galaxies: a high angular resolution mid-infrared view

Authors: A. Alonso-Herrero, R. Poulton, P. F. Roche, A. Hernan-Caballero, I. Aretxaga, M. Martinez-Paredes, C. Ramos Almeida, M. Pereira-Santaella, T. Diaz-Santos, N. A. Levenson, C. Packham, L. Colina, P. Esquej, O. Gonzalez-Martin, K. Ichikawa, M. Imanishi, J. M. Rodriguez Espinosa, C. Telesco

Abstract: We investigate the evolutionary connection between local IR-bright galaxies ($\log L_{\rm IR}\ge 11.4\,L_\odot$) and quasars. We use high angular resolution ($\sim$ 0.3-0.4 arcsec $\sim$ few hundred parsecs) $8-13\,μ$m ground-based spectroscopy to disentangle the AGN mid-IR properties from those of star formation. The comparison between the nuclear $11.3\,μ$m PAH feature emission and that measured… ▽ More We investigate the evolutionary connection between local IR-bright galaxies ($\log L_{\rm IR}\ge 11.4\,L_\odot$) and quasars. We use high angular resolution ($\sim$ 0.3-0.4 arcsec $\sim$ few hundred parsecs) $8-13\,μ$m ground-based spectroscopy to disentangle the AGN mid-IR properties from those of star formation. The comparison between the nuclear $11.3\,μ$m PAH feature emission and that measured with Spitzer/IRS indicates that the star formation is extended over a few kpc in the IR-bright galaxies. The AGN contribution to the total IR luminosity of IR-bright galaxies is lower than in quasars. Although the dust distribution is predicted to change as IR-bright galaxies evolve to IR-bright quasars and then to optical quasars, we show that the AGN mid-IR emission of all the quasars in our sample is not significantly different. In contrast, the nuclear emission of IR-bright galaxies with low AGN contributions appears more heavily embedded in dust although there is no clear trend with the interaction stage or projected nuclear separation. This suggests that the changes in the distribution of the nuclear obscuring material may be taking place rapidly and at different interaction stages washing out the evidence of an evolutionary path. When compared to normal AGN, the nuclear star formation activity of quasars appears to be dimming whereas it is enhanced in some IR-bright nuclei, suggesting that the latter are in an earlier star-formation dominated phase. △ Less

Submitted 31 August, 2016; originally announced August 2016.

Comments: Accepted for publication in MNRAS

The nuclear and extended mid-infrared emission of Seyfert galaxies

Authors: I. García-Bernete, C. Ramos Almeida, J. A. Acosta-Pulido, A. Alonso-Herrero, O. González-Martín, A. Hernán-Caballero, M. Pereira-Santaella, N. A. Levenson, C. Packham, E. S. Perlman, K. Ichikawa, P. Esquej, T. Díaz-Santos

Abstract: We present subarcsecond resolution mid-infrared (MIR) images obtained with 8-10 m-class ground-based telescopes of a complete volume-limited (DL<40 Mpc) sample of 24 Seyfert galaxies selected from the Swift/BAT nine month catalog. We use those MIR images to study the nuclear and circumnuclear emission of the galaxies. Using different methods to classify the MIR morphologies on scales of ~400 pc, w… ▽ More We present subarcsecond resolution mid-infrared (MIR) images obtained with 8-10 m-class ground-based telescopes of a complete volume-limited (DL<40 Mpc) sample of 24 Seyfert galaxies selected from the Swift/BAT nine month catalog. We use those MIR images to study the nuclear and circumnuclear emission of the galaxies. Using different methods to classify the MIR morphologies on scales of ~400 pc, we find that the majority of the galaxies (75-83%) are extended or possibly extended and 17-25% are point-like. This extended emission is compact and it has low surface brightness compared with the nuclear emission, and it represents, on average, ~30% of the total MIR emission of the galaxies in the sample. We find that the galaxies whose circumnuclear MIR emission is dominated by star formation show more extended emission (650+-700 pc) than AGN-dominated systems (300+-100 pc). In general, the galaxies with point-like MIR morphologies are face-on or moderately inclined (b/a~0.4-1.0), and we do not find significant differences between the morphologies of Sy1 and Sy2. We used the nuclear and circumnuclear fluxes to investigate their correlation with different AGN and SF activity indicators. We find that the nuclear MIR emission (the inner ~70 pc) is strongly correlated with the X-ray emission (the harder the X-rays the better the correlation) and with the [O IV] lambda 25.89 micron emission line, indicating that it is AGN-dominated. We find the same results, although with more scatter, for the circumnuclear emission, which indicates that the AGN dominates the MIR emission in the inner ~400 pc of the galaxies, with some contribution from star formation. △ Less

Submitted 24 August, 2016; v1 submitted 23 August, 2016; originally announced August 2016.

Comments: 27 pages, 12 figures, accepted by MNRAS

Investigating the dusty torus of Seyfert galaxies using SOFIA/FORCAST photometry

Authors: Lindsay Fuller, Enrique Lopez-Rodriguez, Chris Packham, Cristina Ramos-Almeida, Almudena Alonso-Herrero, Nancy Levenson, James Radomski, Kohei Ichikawa, Ismael Garcia-Bernete, Omaira Gonzalez-Martin, Tanio Diaz-Santos, Mariela Martinez-Parades

Abstract: We present 31.5 micron imaging photometry of 11 nearby Seyfert galaxies observed from the Stratospheric Observatory For Infrared Astronomy (SOFIA) using the Faint Object infraRed CAmera for the SOFIA Telescope (FORCAST). We tentatively detect extended 31 micron emission for the first time in our sample. In combination with this new data set, subarcsecond resolution 1-18 micron imaging and 7.5-13 m… ▽ More We present 31.5 micron imaging photometry of 11 nearby Seyfert galaxies observed from the Stratospheric Observatory For Infrared Astronomy (SOFIA) using the Faint Object infraRed CAmera for the SOFIA Telescope (FORCAST). We tentatively detect extended 31 micron emission for the first time in our sample. In combination with this new data set, subarcsecond resolution 1-18 micron imaging and 7.5-13 micron spectroscopic observations were used to compute the nuclear spectral energy distribution (SED) of each galaxy. We found that the turnover of the torus emission does not occur at wavelengths <31.5 micron, which we interpret as a lower-limit for the wavelength of peak emission. We used CLUMPY torus models to fit the nuclear infrared (IR) SED and infer trends in the physical parameters of the AGN torus for the galaxies in the sample. Including the 31.5 micron nuclear flux in the SED 1) reduces the number of clumpy torus models compatible with the data, and 2) modifies the model output for the outer radial extent of the torus for 10 of the 11 objects. Specifically, six (60%) objects show a decrease in radial extent while four (40%) show an increase. We find torus outer radii ranging from <1pc to 8.4 pc △ Less

Submitted 26 July, 2016; originally announced July 2016.

Comments: 15 pages, 7 figures. Accepted for publication in MNRAS

About AGN ionization echoes, thermal echoes, and ionization deficits in low redshift Lyman-alpha blobs

Authors: Mischa Schirmer, Sangeeta Malhotra, Nancy A. Levenson, Hai Fu, Rebecca L. Davies, William C. Keel, Paul Torrey, Vardha N. Bennert, Anna Pancoast, James E. H. Turner

Abstract: We report the discovery of 14 Lyman-alpha blobs (LABs) at z~0.3, existing at least 4-7 billion years later in the Universe than all other LABs known. Their optical diameters are 20-70 kpc, and GALEX data imply Ly-alpha luminosities of (0.4-6.3)x10^43 erg/s. Contrary to high-z LABs, they live in low-density areas. They are ionized by AGN, suggesting that cold accretion streams as a power source mus… ▽ More We report the discovery of 14 Lyman-alpha blobs (LABs) at z~0.3, existing at least 4-7 billion years later in the Universe than all other LABs known. Their optical diameters are 20-70 kpc, and GALEX data imply Ly-alpha luminosities of (0.4-6.3)x10^43 erg/s. Contrary to high-z LABs, they live in low-density areas. They are ionized by AGN, suggesting that cold accretion streams as a power source must deplete between z=2 and z=0.3. We also show that transient AGN naturally explain the ionization deficits observed in many LABs: Their Ly-alpha and X-ray fluxes decorrelate below 10^6 years because of the delayed escape of resonantly scattering Ly-alpha photons. High Ly-alpha luminosities do not require currently powerful AGN, independent of obscuration. Chandra X-ray data reveal intrinsically weak AGN, confirming the luminous optical nebulae as impressive ionization echoes. For the first time, we also report mid-infrared thermal echoes from the dusty tori. We conclude that the AGN have faded by 3-4 orders of magnitude within the last 10^(4-5) years, leaving fossil UV, optical and thermal radiation behind. The host galaxies belong to the group of previously discovered Green Bean galaxies (GBs). Gemini optical imaging reveals smooth spheres, mergers, spectacular outflows and ionization cones. Because of their proximity and high flux densities, GBs are perfect targets to study AGN feedback, mode switching and the Ly-alpha escape. The fully calibrated, coadded optical FITS images are publicly available. △ Less

Submitted 21 July, 2016; originally announced July 2016.

Comments: 36 pages, 27 figures. Accepted for publication in MNRAS 2016 July 21. The coadded gri FITS images are publicly available at https://zenodo.org/record/56059

Journal ref: 2016 MNRAS, 463, 1554

A mid-infrared spectroscopic atlas of local active galactic nuclei on sub-arcsecond resolution using GTC/CanariCam

Authors: A. Alonso-Herrero, P. Esquej, P. F. Roche, C. Ramos Almeida, O. Gonzalez-Martin, C. Packham, N. A. Levenson, R. E. Mason, A. Hernan-Caballero, M. Pereira-Santaella, C. Alvarez, I. Aretxaga, E. Lopez-Rodriguez, L. Colina, T. Diaz-Santos, M. Imanishi, J. M. Rodriguez Espinosa, E. Perlman

Abstract: We present an atlas of mid-infrared (mid-IR) ~7.5-13micron spectra of 45 local active galactic nuclei (AGN) obtained with CanariCam on the 10.4m Gran Telescopio CANARIAS (GTC) as part of an ESO/GTC large program. The sample includes Seyferts and other low luminosity AGN (LLAGN) at a median distance of 35Mpc and luminous AGN, namely PG quasars, (U)LIRGs, and radio galaxies (RG) at a median distance… ▽ More We present an atlas of mid-infrared (mid-IR) ~7.5-13micron spectra of 45 local active galactic nuclei (AGN) obtained with CanariCam on the 10.4m Gran Telescopio CANARIAS (GTC) as part of an ESO/GTC large program. The sample includes Seyferts and other low luminosity AGN (LLAGN) at a median distance of 35Mpc and luminous AGN, namely PG quasars, (U)LIRGs, and radio galaxies (RG) at a median distance of 254Mpc. To date, this is the largest mid-IR spectroscopic catalog of local AGN at sub-arcsecond resolution (median 0.3arcsec). The goal of this work is to give an overview of the spectroscopic properties of the sample. The nuclear 12micron luminosities of the AGN span more than four orders of magnitude, nu*Lnu(12micron)~ 3e41-1e46erg/s. In a simple mid-IR spectral index vs. strength of the 9.7micron silicate feature diagram most LLAGN, Seyfert nuclei, PG quasars, and RGs lie in the region occupied by clumpy torus model tracks. However, the mid-IR spectra of some might include contributions from other mechanisms. Most (U)LIRG nuclei in our sample have deeper silicate features and flatter spectral indices than predicted by these models suggesting deeply embedded dust heating sources and/or contribution from star formation. The 11.3micron PAH feature is clearly detected in approximately half of the Seyfert nuclei, LLAGN, and (U)LIRGs. While the RG, PG quasars, and (U)LIRGs in our sample have similar nuclear 12micron luminosities, we do not detect nuclear PAH emission in the RGs and PG quasars. △ Less

Submitted 9 October, 2015; originally announced October 2015.

Comments: Accepted for publication in MNRAS, 22 pages, 8 figures

Sub-arcsec mid-IR observations of NGC 1614: Nuclear star-formation or an intrinsically X-ray weak AGN?

Authors: M. Pereira-Santaella, L. Colina, A. Alonso-Herrero, A. Usero, T. Díaz-Santos, S. García-Burillo, A. Alberdi, O. Gonzalez-Martin, R. Herrero-Illana, M. Imanishi, N. A. Levenson, M. A. Pérez-Torres, C. Ramos Almeida

Abstract: We present new mid-infrared N-band spectroscopy and Q-band photometry of the local luminous infrared galaxy NGC1614, one of the most extreme nearby starbursts. We analyze the mid-IR properties of the nucleus (central 150 pc) and four regions of the bright circumnuclear (diameter~600 pc) star-forming (SF) ring of this object. The nucleus differs from the circumnuclear SF ring by having a strong 8-1… ▽ More We present new mid-infrared N-band spectroscopy and Q-band photometry of the local luminous infrared galaxy NGC1614, one of the most extreme nearby starbursts. We analyze the mid-IR properties of the nucleus (central 150 pc) and four regions of the bright circumnuclear (diameter~600 pc) star-forming (SF) ring of this object. The nucleus differs from the circumnuclear SF ring by having a strong 8-12 micron continuum (low 11.3 micron PAH equivalent width). These characteristics, together with the nuclear X-ray and sub-mm properties, can be explained by an X-ray weak active galactic nucleus (AGN), or by peculiar SF with a short molecular gas depletion time and producing an enhanced radiation field density. In either case, the nuclear luminosity (L(IR) < 6e43 erg/s) is only <5% of the total bolometric luminosity of NGC1614. So this possible AGN does not dominate the energy output in this object. We also compare three star-formation rate (SFR) tracers (Pa$α$, 11.3 micron PAH, and 24 micron emissions) at 150 pc scales in the circumnuclear ring. In general, we find that the SFR is underestimated (overestimated) by a factor of 2-4 (2-3) using the 11.3 micron PAH (24 micron) emission with respect to the extinction corrected Pa$α$ SFR. The former can be explained because we do not include diffuse PAH emission in our measurements, while the latter might indicate that the dust temperature is particularly warmer in the central regions of NGC1614. △ Less

Submitted 25 September, 2015; originally announced September 2015.

Comments: Accepted for publication in MNRAS; 10 pages, 5 figures

A deep look at the nuclear region of UGC 5101 through high angular resolution mid-IR data with GTC/CanariCam

Authors: M. Martínez-Paredes, A. Alonso-Herrero, I. Aretxaga, C. Ramos Almeida, A. Hernán-Caballero, O. González-Martín, M. Pereira-Santaella, C. Packham, A. Asensio Ramos, T. Díaz-Santos, M. Elitzur, P. Esquej, I. García-Bernete, M. Imanishi, N. A. Levenson, J. M. Rodríguez Espinosa

Abstract: We present an analysis of the nuclear infrared (IR, 1.6 to 18 $μ$m) emission of the ultraluminous IR galaxy UGC 5101 to derive the properties of its active galactic nucleus (AGN) and its obscuring material. We use new mid-IR high angular resolution ($0.3-0.5$ arcsec) imaging using the Si-2 filter ($λ_{C}=8.7\, μ$m) and $7.5-13$ $μ$m spectroscopy taken with CanariCam (CC) on the 10.4m Gran Telescop… ▽ More We present an analysis of the nuclear infrared (IR, 1.6 to 18 $μ$m) emission of the ultraluminous IR galaxy UGC 5101 to derive the properties of its active galactic nucleus (AGN) and its obscuring material. We use new mid-IR high angular resolution ($0.3-0.5$ arcsec) imaging using the Si-2 filter ($λ_{C}=8.7\, μ$m) and $7.5-13$ $μ$m spectroscopy taken with CanariCam (CC) on the 10.4m Gran Telescopio CANARIAS. We also use archival HST/NICMOS and Subaru/COMICS imaging and Spitzer/IRS spectroscopy. We estimate the near- and mid-IR unresolved nuclear emission by modelling the imaging data with GALFIT. We decompose the Spitzer/IRS and CC spectra using a power-law component, which represents the emission due to dust heated by the AGN, and a starburst component, both affected by foreground extinction. We model the resulting unresolved near- and mid-IR, and the starburst subtracted CC spectrum with the CLUMPY torus models of Nenkova et al. The derived geometrical properties of the torus, including the large covering factor and the high foreground extinction needed to reproduce the deep $9.7\, μ$m silicate feature, are consistent with the lack of strong AGN signatures in the optical. We derive an AGN bolometric luminosity $L_{bol}\sim1.9\times10^{45}\,$erg s$^{-1}$ that is in good agreement with other estimates in the literature. △ Less

Submitted 15 September, 2015; originally announced September 2015.

Comments: Accepted for publication in MNRAS

Near-Infrared Polarimetric Adaptive Optics Observations of NGC 1068: A torus created by a hydromagnetic outflow wind

Authors: E. Lopez-Rodriguez, C. Packham, T. J. Jones, R. Nikutta, L. McMaster, R. E. Mason, M. Elvis, D. Shenoy, A. Alonso-Herrero, E. Ramirez, O. Gonzalez Martin, S. F. Hoenig, N. A. Levenson, C. Ramos Almeida, E. Perlman

Abstract: We present J' and K' imaging linear polarimetric adaptive optics observations of NGC 1068 using MMT-Pol on the 6.5-m MMT. These observations allow us to study the torus from a magnetohydrodynamical (MHD) framework. In a 0.5" (30 pc) aperture at K', we find that polarisation arising from the passage of radiation from the inner edge of the torus through magnetically aligned dust grains in the clumps… ▽ More We present J' and K' imaging linear polarimetric adaptive optics observations of NGC 1068 using MMT-Pol on the 6.5-m MMT. These observations allow us to study the torus from a magnetohydrodynamical (MHD) framework. In a 0.5" (30 pc) aperture at K', we find that polarisation arising from the passage of radiation from the inner edge of the torus through magnetically aligned dust grains in the clumps is the dominant polarisation mechanism, with an intrinsic polarisation of 7.0%$\pm$2.2%. This result yields a torus magnetic field strength in the range of 4$-$82 mG through paramagnetic alignment, and 139$^{+11}_{-20}$ mG through the Chandrasekhar-Fermi method. The measured position angle (P.A.) of polarisation at K$'$ is found to be similar to the P.A. of the obscuring dusty component at few parsec scales using infrared interferometric techniques. We show that the constant component of the magnetic field is responsible for the alignment of the dust grains, and aligned with the torus axis onto the plane of the sky. Adopting this magnetic field configuration and the physical conditions of the clumps in the MHD outflow wind model, we estimate a mass outflow rate $\le$0.17 M$_{\odot}$ yr$^{-1}$ at 0.4 pc from the central engine for those clumps showing near-infrared dichroism. The models used were able to create the torus in a timescale of $\geq$10$^{5}$ yr with a rotational velocity of $\leq$1228 km s$^{-1}$ at 0.4 pc. We conclude that the evolution, morphology and kinematics of the torus in NGC 1068 can be explained within a MHD framework. △ Less

Submitted 25 June, 2015; v1 submitted 23 June, 2015; originally announced June 2015.

Comments: 14 pages, 4 figures, Accepted by MNRAS

The nuclear and extended infrared emission of the Seyfert galaxy NGC 2992 and the interacting system Arp 245

Authors: I. García-Bernete, C. Ramos Almeida, J. A. Acosta-Pulido, A. Alonso-Herrero, M. Sánchez-Portal, M. Castillo, M. Pereira-Santaella, P. Esquej, O. González-Martín, T. Díaz-Santos, P. Roche, S. Fisher, M. Pović, A. M. Pérez García, I. Valtchanov, C. Packham, N. A. Levenson

Abstract: We present subarcsecond resolution infrared (IR) imaging and mid-IR spectroscopic observations of the Seyfert 1.9 galaxy NGC 2992, obtained with the Gemini North Telescope and the Gran Telescopio CANARIAS (GTC). The N-band image reveals faint extended emission out to ~3 kpc, and the PAH features detected in the GTC/CanariCam 7.5-13 micron spectrum indicate that the bulk of this extended emission i… ▽ More We present subarcsecond resolution infrared (IR) imaging and mid-IR spectroscopic observations of the Seyfert 1.9 galaxy NGC 2992, obtained with the Gemini North Telescope and the Gran Telescopio CANARIAS (GTC). The N-band image reveals faint extended emission out to ~3 kpc, and the PAH features detected in the GTC/CanariCam 7.5-13 micron spectrum indicate that the bulk of this extended emission is dust heated by star formation. We also report arcsecond resolution MIR and far-IR imaging of the interacting system Arp 245, taken with the Spitzer Space Telescope and the Herschel Space Observatory. Using these data, we obtain nuclear fluxes using different methods and find that we can only recover the nuclear fluxes obtained from the subarcsecond data at 20-25 micron, where the AGN emission dominates. We fitted the nuclear IR spectral energy distribution of NGC 2992, including the GTC/CanariCam nuclear spectrum (~50 pc), with clumpy torus models. We then used the best-fitting torus model to decompose the Spitzer/IRS 5-30 spectrum (~630 pc) in AGN and starburst components, using different starburst templates. We find that, whereas at shorter mid-IR wavelengths the starburst component dominates (64% at 6 micron), the AGN component reaches 90% at 20 micron. We finally obtained dust masses, temperatures and star formation rates for the different components of the Arp 245 system and find similar values for NGC 2992 and NGC 2993. These measurements are within those reported for other interacting systems in the first stages of the interaction. △ Less

Submitted 16 February, 2015; originally announced February 2015.

Comments: 20 pages, 12 figures, accepted by MNRAS

On the difference of torus geometry between hidden and non-hidden broad line active galactic nuclei

Authors: Kohei Ichikawa, Christopher Packham, Cristina Ramos Almeida, Andres Asensio Ramos, Almudena Alonso-Herrero, Omaira Gonzalez-Martin, Enrique Lopez-Rodriguez, Yoshihiro Ueda, Tanio Diaz-Santos, Moshe Elitzur, Sebastian F. Hoenig, Masatoshi Imanishi, Nancy A. Levenson, Rachel E. Mason, Eric S. Perlman, Crystal D. Alsip

Abstract: We present results from the fitting of infrared (IR) spectral energy distributions of 21 active galactic nuclei (AGN) with clumpy torus models. We compiled high spatial resolution ($\sim 0.3$--$0.7$ arcsec) mid-IR $N$-band spectroscopy, $Q$-band imaging and nuclear near- and mid-IR photometry from the literature. Combining these nuclear near- and mid-IR observations, far-IR photometry and clumpy t… ▽ More We present results from the fitting of infrared (IR) spectral energy distributions of 21 active galactic nuclei (AGN) with clumpy torus models. We compiled high spatial resolution ($\sim 0.3$--$0.7$ arcsec) mid-IR $N$-band spectroscopy, $Q$-band imaging and nuclear near- and mid-IR photometry from the literature. Combining these nuclear near- and mid-IR observations, far-IR photometry and clumpy torus models, enables us to put constraints on the torus properties and geometry. We divide the sample into three types according to the broad line region (BLR) properties; type-1s, type-2s with scattered or hidden broad line region (HBLR) previously observed, and type-2s without any published HBLR signature (NHBLR). Comparing the torus model parameters gives us the first quantitative torus geometrical view for each subgroup. We find that NHBLR AGN have smaller torus opening angles and larger covering factors than those of HBLR AGN. This suggests that the chance to observe scattered (polarized) flux from the BLR in NHBLR could be reduced by the dual effects of (a) less scattering medium due to the reduced scattering volume given the small torus opening angle and (b) the increased torus obscuration between the observer and the scattering region. These effects give a reasonable explanation for the lack of observed HBLR in some type-2 AGN. △ Less

Submitted 26 January, 2015; originally announced January 2015.

Comments: 13 pages, 5 figures, accepted for publication in ApJ

Mid-IR Spectra of Type Ia SN 2014J in M82 Spanning the First Four Months

Authors: Charles M. Telesco, Peter Höflich, Dan Li, Carlos Álvarez, Christopher M. Wright, Peter J. Barnes, Sergio Fernández, James H. Hough, N. A. Levenson, Naibí Mariñas, Christopher Packham, Eric Pantin, Rafael Rebolo, Patrick Roche, Han Zhang

Abstract: We present a time series of 8 - 13 $μ$m spectra and photometry for SN 2014J obtained 57, 81, 108, and 137 d after the explosion using CanariCam on the Gran Telescopio Canarias. This is the first mid-IR time series ever obtained for a Type Ia supernova. These observations can be understood within the framework of the delayed detonation model and the production of $\sim$0.6 $\rm M_\odot$ of $^{56}$N… ▽ More We present a time series of 8 - 13 $μ$m spectra and photometry for SN 2014J obtained 57, 81, 108, and 137 d after the explosion using CanariCam on the Gran Telescopio Canarias. This is the first mid-IR time series ever obtained for a Type Ia supernova. These observations can be understood within the framework of the delayed detonation model and the production of $\sim$0.6 $\rm M_\odot$ of $^{56}$Ni, consistent with the observed brightness, the brightness decline relation, and the $γ$-ray fluxes. The [Co III] line at 11.888 $μ$m is particularly useful for evaluating the time evolution of the photosphere and measuring the amount of $^{56}$Ni and thus the mass of the ejecta. Late-time line profiles of SN 2014J are rather symmetric and not shifted in the rest frame. We see Argon emission, which provides a unique probe of mixing in the transition layer between incomplete burning and nuclear statistical equilibrium. We may see [Fe III] and [Ni IV] emission, both of which are observed to be substantially stronger than indicated by our models. If the latter identification is correct, then we are likely observing stable Ni, which might imply central mixing. In addition, electron capture, also required for stable Ni, requires densities larger than $\sim$$1 \times 10^9$ g cm$^{-3}$, which are expected to be present only in white dwarfs close to the Chandrasekhar limit. This study demonstrates that mid-IR studies of Type Ia supernovae are feasible from the ground and provide unique information, but it also indicates the need for better atomic data. △ Less

Submitted 4 November, 2014; v1 submitted 7 September, 2014; originally announced September 2014.

Comments: 20 pages, 4 figures, accepted by ApJ

A mid-infrared view of the inner parsecs of the Seyfert galaxy Mrk 1066 using CanariCam/GTC

Authors: C. Ramos Almeida, A. Alonso-Herrero, P. Esquej, O. Gonzalez-Martin, R. A. Riffel, I. Garcia-Bernete, J. M. Rodriguez Espinosa, C. Packham, N. A. Levenson, P. Roche, T. Diaz-Santos, I. Aretxaga, C. Alvarez

Abstract: We present mid-infrared (MIR) imaging and spectroscopic data of the Seyfert 2 galaxy Mrk 1066 obtained with CanariCam (CC) on the 10.4 m Gran Telescopio CANARIAS (GTC). The galaxy was observed in imaging mode with an angular resolution of 0.24 arcsec (54 pc) in the Si-2 filter (8.7 micron). The image reveals a series of star-forming knots within the central ~400 pc, after subtracting the dominant… ▽ More We present mid-infrared (MIR) imaging and spectroscopic data of the Seyfert 2 galaxy Mrk 1066 obtained with CanariCam (CC) on the 10.4 m Gran Telescopio CANARIAS (GTC). The galaxy was observed in imaging mode with an angular resolution of 0.24 arcsec (54 pc) in the Si-2 filter (8.7 micron). The image reveals a series of star-forming knots within the central ~400 pc, after subtracting the dominant active galactic nucleus (AGN) component. We also subtracted this AGN unresolved component from the 8-13 micron spectra of the knots and the nucleus, and measured equivalent widths (EWs) of the 11.3 micron Polycyclic Aromatic Hydrocarbon (PAH) feature which are typical of pure starburst galaxies. This EW is larger in the nucleus than in the knots, confirming that, at least in the case of Mrk 1066, the AGN dilutes, rather than destroys, the molecules responsible for the 11.3 micron PAH emission. By comparing the nuclear GTC/CC spectrum with the Spitzer/IRS spectrum of the galaxy, we find that the AGN component that dominates the continuum emission at wavelengths below 15 micron on scales of ~60 pc (90-100%) decreases to 35-50% when the emission of the central ~830 pc is considered. On the other hand, the AGN contribution dominates the 15-25 micron emission (75%) on the scales probed by Spitzer/IRS. We find a good match between the MIR morphology of Mrk 1066 and the extended Pab, Brg and [O III] emission. This coincidence implies that the 8.7 micron emission is probing star formation, dust in the narrow-line region, and the oval structure previously detected in the near-infrared. On the other hand, the Chandra soft X-ray morphology does not match any of the previous, contrary to what it is generally assumed for Seyfert galaxies. △ Less

Submitted 2 September, 2014; originally announced September 2014.

Comments: 16 pages including appendix, accepted for publication in MNRAS

The Gemini Observatory Fast Turnaround Program

Authors: R. E. Mason, S. Cote, M. Kissler-Patig, N. A. Levenson, A. Adamson, C. Emmanuel, D. Crabtree

Abstract: Gemini's Fast Turnaround program is intended to greatly decrease the time from having an idea to acquiring the supporting data. The scheme will offer monthly proposal submission opportunities, and proposals will be reviewed by the principal investigators or co-investigators of other proposals submitted during the same round. Here, we set out the design of the system and outline the plan for its im… ▽ More Gemini's Fast Turnaround program is intended to greatly decrease the time from having an idea to acquiring the supporting data. The scheme will offer monthly proposal submission opportunities, and proposals will be reviewed by the principal investigators or co-investigators of other proposals submitted during the same round. Here, we set out the design of the system and outline the plan for its implementation, leading to the launch of a pilot program at Gemini North in January 2015. △ Less

Submitted 25 August, 2014; originally announced August 2014.

Comments: SPIE proceedings

Polarized Mid-Infrared Synchrotron Emission in the Core of Cygnus A

Authors: E. Lopez-Rodriguez, C. Packham, C. Tadhunter, R. Mason, E. Perlman, A. Alonso-Herrero, C. Ramos Almeida, K. Ichikawa, N. A. Levenson, J. M. Rodrıguez-Espinosa, C. A. Alvarez, E. A. Ramırez, C. M. Telesco

Abstract: We present high-angular (~0.4") resolution mid-infrared (MIR) polarimetric observations in the 8.7 $μ$m and 11.6 $μ$m filters of Cygnus A using CanariCam on the 10.4-m Gran Telescopio CANARIAS. A highly polarized nucleus is observed with a degree of polarization of 11${\pm}$3% and 12${\pm}$3% and position angle of polarization of 27${\pm}$8 degrees and 35${\pm}$8 degrees in a 0.38" (~380 pc) apert… ▽ More We present high-angular (~0.4") resolution mid-infrared (MIR) polarimetric observations in the 8.7 $μ$m and 11.6 $μ$m filters of Cygnus A using CanariCam on the 10.4-m Gran Telescopio CANARIAS. A highly polarized nucleus is observed with a degree of polarization of 11${\pm}$3% and 12${\pm}$3% and position angle of polarization of 27${\pm}$8 degrees and 35${\pm}$8 degrees in a 0.38" (~380 pc) aperture for each filter. The observed rising of the polarized flux density with increasing wavelength is consistent with synchrotron radiation from the pc-scale jet close to the core of Cygnus A. Based on our polarization model, the synchrotron emission from the pc-scale jet is estimated to be 14% and 17% of the total flux density in the 8.7 $μ$m and 11.6 $μ$m filters, respectively. A blackbody component with a characteristic temperature of 220 K accounts for >75% of the observed MIR total flux density. The blackbody emission arises from a combination of (1) dust emission in the torus; and (2) diffuse dust emission around the nuclear region, but the contributions of the two components cannot be well constrained in these observations. △ Less

Submitted 23 July, 2014; originally announced July 2014.

Comments: 8 pages, 2 figures (Accepted for publication in ApJ)

Nuclear 11.3$μ$m PAH emission in local active galactic nuclei

Authors: A. Alonso-Herrero, C. Ramos Almeida, P. Esquej, P. F. Roche, A. Hernan-Caballero, S. F. Hoenig, O. Gonzalez-Martin, I. Aretxaga, R. E. Mason, C. Packham, N. A. Levenson, J. M. Rodriguez Espinosa, R. Siebenmorgen, M. Pereira-Santaella, T. Diaz-Santos, L. Colina, C. Alvarez, C. M. Telesco

Abstract: We present Gran Telescopio CANARIAS CanariCam 8.7$μ$m imaging and 7.5-13$μ$m spectroscopy of six local systems known to host an active galactic nucleus (AGN) and have nuclear star formation. Our main goal is to investigate whether the molecules responsible for the 11.3$μ$m polyclyclic aromatic hydrocarbon (PAH) feature are destroyed in the close vicinity of an AGN. We detect 11.3$μ$m PAH feature e… ▽ More We present Gran Telescopio CANARIAS CanariCam 8.7$μ$m imaging and 7.5-13$μ$m spectroscopy of six local systems known to host an active galactic nucleus (AGN) and have nuclear star formation. Our main goal is to investigate whether the molecules responsible for the 11.3$μ$m polyclyclic aromatic hydrocarbon (PAH) feature are destroyed in the close vicinity of an AGN. We detect 11.3$μ$m PAH feature emission in the nuclear regions of the galaxies as well as extended PAH emission over a few hundred parsecs. The equivalent width (EW) of the feature shows a minimum at the nucleus but increases with increasing radial distances, reaching typical star-forming values a few hundred parsecs away from the nucleus. The reduced nuclear EW are interpreted as due to increased dilution from the AGN continuum rather than destruction of the PAH molecules. We conclude that at least those molecules responsible for the 11.3$μ$m PAH feature survive in the nuclear environments as close as 10pc from the AGN and for Seyfert-like AGN luminosities. We propose that material in the dusty tori, nuclear gas disks, and/or host galaxies of AGN is likely to provide the column densities necessary to protect the PAH molecules from the AGN radiation field. △ Less

Submitted 4 July, 2014; originally announced July 2014.

Comments: Accepted for publication in MNRAS, 19 pages, 13 figures

Near- to Mid- Infrared Imaging and Spectroscopy of Two Buried AGNs of the Nearby Merging Galaxy NGC 6240 with Subaru/IRCS+AO and GTC/CanariCam

Authors: Tamami I. Mori, Masatoshi Imanishi, Almudena Alonso-Herrero, Chris Packham, Cristina Ramos Almeida, Robert Nikutta, OmairaGonzález-Martín, Eric Perlman, Yuriko Saito, Nancy A. Levenson

Abstract: We report near-infrared K', L', and M' band imaging observations of the nearby merging galaxy NGC 6240 with the Infrared Camera and Spectrograph on the Subaru telescope. The observations were performed with the assistance of the Subaru Adaptive Optics System, and the achieved spatial resolutions were around 0.10--0.20$^{\prime\prime}$. We also obtained new mid-infrared imaging in the Si-2 filter b… ▽ More We report near-infrared K', L', and M' band imaging observations of the nearby merging galaxy NGC 6240 with the Infrared Camera and Spectrograph on the Subaru telescope. The observations were performed with the assistance of the Subaru Adaptive Optics System, and the achieved spatial resolutions were around 0.10--0.20$^{\prime\prime}$. We also obtained new mid-infrared imaging in the Si-2 filter band (8.7$μ$m) and N-band (7.5--13$μ$m) spectroscopy of this galaxy with the CanariCam on the Gran Telescopio Canarias with a spatial resolution of 0.4--0.5$^{\prime\prime}$. In the K' band image the two nuclei of the galaxy each show a double-peak suggesting the complex geometry of the source, while the L', M', and Si-2 band images show single compact structures in each of the two nuclei. Assuming that the center core observed at wavelengths longer than 3.8$μ$m is associated with dust heated by the buried AGN, we re-evaluated the spectral energy distributions (SEDs) of the southern nucleus from 2 to 30$μ$m with the additional literature values, and performed the SED+spectroscopy fitting using the clumpy torus models of Nenkova et al. (2008) and a Bayesian fitting approach. The model fit suggests that the high covering factor torus emission in the southern nucleus is also obscured by foreground dust in the host galaxy. The estimated AGN bolometric luminosity of the southern nucleus, $L_{\rm bol}({\rm AGN})\sim1\times10^{45}$ [erg$\cdot$s$^{-1}$], accounts for approximately 40% of the whole luminosity of the system. △ Less

Submitted 30 June, 2014; originally announced June 2014.

Comments: 15 pages, 3 figures, accepted for publication in PASJ

Subaru Spectroscopy and Spectral Modeling of Cygnus A

Authors: Matthew J. Merlo, Eric S. Perlman, Robert Nikutta, Christopher Packham, Moshe Elitzur, Masatoshi Imanishi, N. A. Levenson, James T. Radomski, Itziar Aretxaga

Abstract: We present high angular resolution ($\sim$0.5$^\prime$$^\prime$) MIR spectra of the powerful radio galaxy, Cygnus A, obtained with the Subaru telescope. The overall shape of the spectra agree with previous high angular resolution MIR observations, as well as previous Spitzer spectra. Our spectra, both on and off nucleus, show a deep silicate absorption feature. The absorption feature can be modele… ▽ More We present high angular resolution ($\sim$0.5$^\prime$$^\prime$) MIR spectra of the powerful radio galaxy, Cygnus A, obtained with the Subaru telescope. The overall shape of the spectra agree with previous high angular resolution MIR observations, as well as previous Spitzer spectra. Our spectra, both on and off nucleus, show a deep silicate absorption feature. The absorption feature can be modeled with a blackbody obscured by cold dust or a clumpy torus. The deep silicate feature is best fit by a simple model of a screened blackbody, suggesting foreground absorption plays a significant, if not dominant role, in shaping the spectrum of Cygnus A. This foreground absorption prevents a clear view of the central engine and surrounding torus, making it difficult to quantify the extent the torus attributes to the obscuration of the central engine, but does not eliminate the need for a torus in Cygnus A. △ Less

Submitted 1 April, 2014; originally announced April 2014.

Investigating the sensitivity of observed spectral energy distributions to clumpy torus properties in Seyfert galaxies

Authors: C. Ramos Almeida, A. Alonso-Herrero, N. A. Levenson, A. Asensio Ramos, J. M. Rodriguez Espinosa, O. Gonzalez-Martin, C. Packham, M. Martinez

Abstract: We present nuclear spectral energy distributions (SEDs) from 1 to 18 micron of a small sample of nearby, nearly face-on and undisturbed Seyfert galaxies without prominent nuclear dust lanes. These nuclear SEDs probe the central ~35 pc of the galaxies, on average, and include photometric and spectroscopic infrared (IR) data. We use these SEDs, the clumpy torus models of Nenkova et al. and a Bayesia… ▽ More We present nuclear spectral energy distributions (SEDs) from 1 to 18 micron of a small sample of nearby, nearly face-on and undisturbed Seyfert galaxies without prominent nuclear dust lanes. These nuclear SEDs probe the central ~35 pc of the galaxies, on average, and include photometric and spectroscopic infrared (IR) data. We use these SEDs, the clumpy torus models of Nenkova et al. and a Bayesian approach to study the sensitivity of different IR wavelengths to the torus parameters. We find that high angular resolution 8-13 micron spectroscopy alone reliably constrains the number of clumps and their optical depth (N0 and tau_V). On the other hand, we need a combination of mid- and near-IR subarcsecond resolution photometry to constrain torus width and inclination, as well as the radial distribution of the clouds (sigma, i and q). For flat radial profiles (q=0,1), it is possible to constrain the extent of the mid-IR-emitting dust within the torus (Y) when N-band spectroscopy is available, in addition to near-IR photometry. Finally, by fitting different combinations of average and individual Seyfert 1 and Seyfert 2 data, we find that, in general, for undisturbed, nearly face-on Seyferts without prominent nuclear dust lanes, the minimum combination of data necessary to reliably constrain all the torus parameters is J+K+M-band photometry + N-band spectroscopy. △ Less

Submitted 3 February, 2014; originally announced February 2014.

Comments: 15 pages, 11 figures, accepted in MNRAS

Uncovering the deeply embedded AGN activity in the nuclear regions of the interacting galaxy Arp299

Authors: A. Alonso-Herrero, P. F. Roche, P. Esquej, O. Gonzalez-Martin, M. Pereira-Santaella, C. Ramos Almeida, N. A. Levenson, C. Packham, A. Asensio Ramos, R. E. Mason, J. M. Rodriguez Espinosa, C. Alvarez, L. Colina, I. Aretxaga, T. Diaz-Santos, E. Perlman, C. M. Telesco

Abstract: We present mid-infrared (MIR) 8-13micron spectroscopy of the nuclear regions of the interacting galaxy Arp299 (IC694+NGC3690) obtained with CanariCam (CC) on the 10.4m Gran Telescopio Canarias (GTC). The high angular resolution (~0.3-0.6arcsec) of the data allows us to probe nuclear physical scales between 60 and 120pc, which is a factor of 10 improvement over previous MIR spectroscopic observatio… ▽ More We present mid-infrared (MIR) 8-13micron spectroscopy of the nuclear regions of the interacting galaxy Arp299 (IC694+NGC3690) obtained with CanariCam (CC) on the 10.4m Gran Telescopio Canarias (GTC). The high angular resolution (~0.3-0.6arcsec) of the data allows us to probe nuclear physical scales between 60 and 120pc, which is a factor of 10 improvement over previous MIR spectroscopic observations of this system. The GTC/CC spectroscopy displays evidence of deeply embedded Active Galactic Nucleus (AGN) activity in both nuclei. The GTC/CC nuclear spectrum of NGC3690/Arp299-B1 can be explained as emission from AGN-heated dust in a clumpy torus with both a high covering factor and high extinction along the line of sight. The estimated bolometric luminosity of the AGN in NGC3690 is 3.2(+/-0.6)x10^44 erg/s. The nuclear GTC/CC spectrum of IC694/Arp299-A shows 11.3micron polycyclic aromatic hydrocarbon (PAH) emission stemming from a deeply embedded (A_V~24mag) region of less than 120pc in size. There is also a continuum-emitting dust component. If associated with th putative AGN in IC694, we estimate that it would be approximately 5 times less luminous than the AGN in NGC3690. The presence of dual AGN activity makes Arp299 a good example to study such phenomenon in the early coalescence phase of interacting galaxies. △ Less

Submitted 14 November, 2013; originally announced November 2013.

Comments: Accepted for publication in ApJ Letters

Nuclear star formation activity and black hole accretion in nearby Seyfert galaxies

Authors: P. Esquej, A. Alonso-Herrero, O. González-Martín, S. F. Hönig, A. Hernán Caballero, P. F. Roche, C. Ramos Almeida, R. E. Mason, T. Díaz-Santos, N. A. Levenson, I. Aretxaga, J. M. Rodríguez Espinosa, C. Packham

Abstract: Recent theoretical and observational works indicate the presence of a correlation between the star formation rate (SFR) and the active galactic nuclei (AGN) luminosity (and, therefore, the black hole accretion rate) of Seyfert galaxies. This suggests a physical connection between the gas forming stars on kpc scales and the gas on sub-pc scales that is feeding the black hole. We compiled the larges… ▽ More Recent theoretical and observational works indicate the presence of a correlation between the star formation rate (SFR) and the active galactic nuclei (AGN) luminosity (and, therefore, the black hole accretion rate) of Seyfert galaxies. This suggests a physical connection between the gas forming stars on kpc scales and the gas on sub-pc scales that is feeding the black hole. We compiled the largest sample of Seyfert galaxies to date with high angular resolution (0.4-0.8 arcsec) mid-infrared (8-13 micron) spectroscopy. The sample includes 29 Seyfert galaxies drawn from the AGN Revised Shapley-Ames catalogue. At a median distance of 33 Mpc, our data allow us to probe nuclear regions on scales of 65 pc (median value). We found no general evidence of suppression of the 11.3 micron polycyclic aromatic hydrocarbon (PAH) emission in the vicinity of these AGN, and used this feature as a proxy for the SFR. We detected the 11.3 micron PAH feature in the nuclear spectra of 45% of our sample. The derived nuclear SFRs are, on average, five times lower than those measured in circumnuclear regions of 600 pc in size (median value). However, the projected nuclear SFR densities are a factor of 20 higher than those measured on circumnuclear scales. This indicates that the SF activity per unit area in the central 65 pc of Seyfert galaxies is much higher than at larger distances from their nuclei. We studied the connection between the nuclear SFR and the black hole accretion rate and showed that numerical simulations reproduce fairly well our observed relation. △ Less

Submitted 4 November, 2013; originally announced November 2013.

Comments: 16 pages, 8 figures. Accepted for publication in ApJ

The role of the accretion disk, dust, and jets in the IR emission of low-luminosity AGN

Authors: R. E. Mason, C. Ramos Almeida, N. A. Levenson, R. Nemmen, A. Alonso-Herrero

Abstract: We use recent high-resolution infrared (IR; 1 - 20 micron) photometry to examine the origin of the IR emission in low-luminosity active galactic nuclei (LLAGN). The data are compared with published model fits that describe the spectral energy distribution of LLAGN in terms of an advection-dominated accretion flow (ADAF), truncated thin accretion disk, and jet. The truncated disk in these models is… ▽ More We use recent high-resolution infrared (IR; 1 - 20 micron) photometry to examine the origin of the IR emission in low-luminosity active galactic nuclei (LLAGN). The data are compared with published model fits that describe the spectral energy distribution of LLAGN in terms of an advection-dominated accretion flow (ADAF), truncated thin accretion disk, and jet. The truncated disk in these models is usually not luminous enough to explain the observed IR emission, and in all cases its spectral shape is much narrower than the broad IR peaks in the data. Synchrotron radiation from the jet appears to be important in very radio-loud nuclei, but the detection of strong silicate emission features in many objects indicates that dust must also contribute. We investigate this point by fitting the IR SED of NGC 3998 using dusty torus and optically thin (tau_MIR ~ 1) dust shell models. While more detailed modeling is necessary, these initial results suggest that dust may account for the nuclear mid-IR emission of many LLAGN. △ Less

Submitted 7 October, 2013; originally announced October 2013.

Comments: ApJ, accepted; 6 pages, 4 figures