Showing 1–50 of 75 results for author: Lopez-Rodriguez, E

AtLAST Science Overview Report

Authors: Mark Booth, Pamela Klaassen, Claudia Cicone, Tony Mroczkowski, Martin A. Cordiner, Luca Di Mascolo, Doug Johnstone, Eelco van Kampen, Minju M. Lee, Daizhong Liu, John Orlowski-Scherer, Amélie Saintonge, Matthew W. L. Smith, Alexander Thelen, Sven Wedemeyer, Kazunori Akiyama, Stefano Andreon, Doris Arzoumanian, Tom J. L. C. Bakx, Caroline Bot, Geoffrey Bower, Roman Brajša, Chian-Chou Chen, Elisabete da Cunha, David Eden , et al. (59 additional authors not shown)

Abstract: Submillimeter and millimeter wavelengths provide a unique view of the Universe, from the gas and dust that fills and surrounds galaxies to the chromosphere of our own Sun. Current single-dish facilities have presented a tantalising view of the brightest (sub-)mm sources, and interferometers have provided the exquisite resolution necessary to analyse the details in small fields, but there are still… ▽ More Submillimeter and millimeter wavelengths provide a unique view of the Universe, from the gas and dust that fills and surrounds galaxies to the chromosphere of our own Sun. Current single-dish facilities have presented a tantalising view of the brightest (sub-)mm sources, and interferometers have provided the exquisite resolution necessary to analyse the details in small fields, but there are still many open questions that cannot be answered with current facilities. In this report we summarise the science that is guiding the design of the Atacama Large Aperture Submillimeter Telescope (AtLAST). We demonstrate how tranformational advances in topics including star formation in high redshift galaxies, the diffuse circumgalactic medium, Galactic ecology, cometary compositions and solar flares motivate the need for a 50m, single-dish telescope with a 1-2 degree field of view and a new generation of highly multiplexed continuum and spectral cameras. AtLAST will have the resolution to drastically lower the confusion limit compared to current single-dish facilities, whilst also being able to rapidly map large areas of the sky and detect extended, diffuse structures. Its high sensitivity and large field of view will open up the field of submillimeter transient science by increasing the probability of serendipitous detections. Finally, the science cases listed here motivate the need for a highly flexible operations model capable of short observations of individual targets, large surveys, monitoring programmes, target of opportunity observations and coordinated observations with other observatories. AtLAST aims to be a sustainable, upgradeable, multipurpose facility that will deliver orders of magnitude increases in sensitivity and mapping speeds over current and planned submillimeter observatories. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 47 pages, 12 figures. For further details on AtLAST see https://atlast.uio.no

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

Forgotten treasures in the HST/FOC UV imaging polarimetric archives of active galactic nuclei. II. Mrk 463E

Authors: Thibault Barnouin, Frédéric Marin, Enrique Lopez-Rodriguez

Abstract: The Mrk 463 system is known to host two powerful sources separated by about 4 kpc, both identified as active galactic nuclei (AGN). This makes the Mrk 463 system a unique laboratory to study the geometry and dynamics of galaxy merging and its relation to AGN duty cycles. The eastern nuclei, Mrk 463E, is the brightest of the two and thus a prime target for a polarimetric study. It is classified as… ▽ More The Mrk 463 system is known to host two powerful sources separated by about 4 kpc, both identified as active galactic nuclei (AGN). This makes the Mrk 463 system a unique laboratory to study the geometry and dynamics of galaxy merging and its relation to AGN duty cycles. The eastern nuclei, Mrk 463E, is the brightest of the two and thus a prime target for a polarimetric study. It is classified as a Seyfert 2 galaxy, meaning that one could expect large polarization degrees from scattering off electrons and dust in the polar winds. In the continuity of our series of papers, we reduced archived and previously unpublished polarization observations obtained with the Faint Object Camera (FOC) onboard the Hubble Space Telescope (HST), to obtain a high resolution near ultraviolet (near-UV) polarization map of the Mrk 463E nuclei. We coupled this map to near infrared (NIR) and X-ray observations to get a clear picture of the geometric arrangement of matter around the core of Mrk 463E. We found that the nucleus location is further South from the optical peak flux than previously estimated. The strongly polarized conical wind has a half-opening angle of ~15° and display three main periods of mass ejection. Its polarization allowed us to estimate the AGN inclination towards the observer (~55°) Finally, our maps revealed a gas streamer connecting Mrk 463E and Mrk 463W, with a tentative detection of a large kpc-scale ordered magnetic field connecting both galaxies. This unpublished observation turned out to offer more than the original proposal asked for and allowed to derive tight geometric and dynamical constraints for Mrk 463E. High resolution radio maps and IR polarimetry are now necessary to further study the jet and the newly discovered gas streamer. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: 9 pages, 7 figures, submitted Astronomy & Astrophysics, comments are welcome

MSC Class: 85-04 (Primary); 85-08; 85A25 (Secondary) ACM Class: J.2; I.4.1; I.4.7; I.6.4

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

Towards an observationally motivated AGN dusty torus model. I. Dust chemical composition from the modeling of Spitzer spectra

Authors: Omar Ulises Reyes-Amador, Jacopo Fritz, Omaira González-Martín, Sundar Srinivasan, Maarten Baes, Enrique Lopez-Rodriguez, Natalia Osorio-Clavijo, Cesar Iván Victoria-Ceballos, Marko Stalevski, C. Ramos Almeida

Abstract: Spectral energy distribution (SED) fitting is one of most commonly used techniques to study the dust properties in Active Galactic Nuclei (AGN). Works implementing this technique commonly use radiative transfer models that assume a variety of dust properties. Despite the key role of this aspect, limited effort has been put forward to explore the chemical composition, the role of different optical… ▽ More Spectral energy distribution (SED) fitting is one of most commonly used techniques to study the dust properties in Active Galactic Nuclei (AGN). Works implementing this technique commonly use radiative transfer models that assume a variety of dust properties. Despite the key role of this aspect, limited effort has been put forward to explore the chemical composition, the role of different optical properties, and the grain size distribution of dust, all of which can have a substantial impact on the theoretical radiative transfer calculations. In this work, we explore the role of the dust chemical composition in the AGN dusty torus through SED fitting to \emph{Spitzer}/IRS spectra of a sample of 49 nearby AGN with silicate features in emission. We implement a mineralogy model including the popular astronomical silicates and a set of oxides and amorphous silicates with different grain sizes. We find that best fits use principally porous alumina, periclase, and olivine. In terms of mass fractions, $\sim99\%$ of the dust is composed of dust grains of size $\rm{0.1 μm}$, with a $<1\%$ contribution from $\rm{3 μm}$ grains. Moreover, the astronomical silicates have a very low occurrence in the best fits, suggesting that they are not the most suited dust species to reproduce the silicate features in our sample. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 16 pages, 13 figures, Accepted by MNRAS

SAUNAS I: Searching for Low Surface Brightness X-ray Emission with Chandra/ACIS

Authors: Alejandro S. Borlaff, Pamela M. Marcum, Mehmet Alpaslan, Pasquale Temi, Nushkia Chamba, Drew S. Chojnowski, Michael N. Fanelli, Anton M. Koekemoer, Seppo Laine, Enrique Lopez-Rodriguez, Aneta Siemiginowska

Abstract: We present SAUNAS (Selective Amplification of Ultra Noisy Astronomical Signal), a pipeline designed for detecting diffuse X-ray emission in the data obtained with the Advanced CCD Imaging Spectrometer (ACIS) of the Chandra X-ray Observatory. SAUNAS queries the available observations in the Chandra archive, performs photometric calibration, PSF (point spread function) modeling, and deconvolution, p… ▽ More We present SAUNAS (Selective Amplification of Ultra Noisy Astronomical Signal), a pipeline designed for detecting diffuse X-ray emission in the data obtained with the Advanced CCD Imaging Spectrometer (ACIS) of the Chandra X-ray Observatory. SAUNAS queries the available observations in the Chandra archive, performs photometric calibration, PSF (point spread function) modeling, and deconvolution, point-source removal, adaptive smoothing, and background correction. This pipeline builds on existing and well-tested software including CIAO, VorBin, and LIRA. We characterize the performance of SAUNAS through several quality performance tests, and demonstrate the broad applications and capabilities of SAUNAS using two galaxies already known to show X-ray emitting structures. SAUNAS successfully detects the 30 kpc X-ray super-wind of NGC 3079 using Chandra/ACIS datasets, matching the spatial distribution detected with more sensitive XMM-Newton observations. The analysis performed by SAUNAS reveals an extended low surface brightness source in the field of UGC 5101 in the 0.3-1.0 keV and 1.0-2.0 keV bands. This source is potentially a background galaxy cluster or a hot gas plume associated with UGC 5101. SAUNAS demonstrates its ability to recover previously undetected structures in archival data, expanding exploration into the low surface brightness X-ray universe with Chandra/ACIS. △ Less

Submitted 2 May, 2024; originally announced May 2024.

Comments: Accepted for publication in The Astrophysical Journal. 27 pages, 19 figures. The SAUNAS X-ray surface brightness and signal-to-noise maps of NGC 3079 and UGC 5101 are publicly available in Zenodo in FITS format: https://zenodo.org/records/10892485. DOI: 10.5281/zenodo.10892485. The list of Chandra datasets, obtained by the Chandra X-ray Observatory, are contained in DOI: 10.25574/cdc.225

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!

Confusion of extragalactic sources in the far infrared: a baseline assessment of the performance of PRIMAger in intensity and polarization

Authors: Matthieu Béthermin, Alberto D. Bolatto, François Boulanger, Charles M. Bradford, Denis Burgarella, Laure Ciesla, James Donnellan, Brandon S. Hensley, Jason Glenn, Guilaine Lagache, Enrique Lopez-Rodriguez, Seb Oliver, Alexandra Pope, Marc Sauvage

Abstract: Because of their limited angular resolution, far-infrared telescopes are usually affected by confusion phenomenon. Since several galaxies can be located in the same instrumental beam, only the brightest objects emerge from the fluctuations caused by fainter sources. The probe far-infrared mission for astrophysics imager (PRIMAger) will observe the mid- and far-infrared (25-235 $μ$m) sky both in in… ▽ More Because of their limited angular resolution, far-infrared telescopes are usually affected by confusion phenomenon. Since several galaxies can be located in the same instrumental beam, only the brightest objects emerge from the fluctuations caused by fainter sources. The probe far-infrared mission for astrophysics imager (PRIMAger) will observe the mid- and far-infrared (25-235 $μ$m) sky both in intensity and polarization. We aim to provide predictions of the confusion level and its consequences for future surveys. We produced simulated PRIMAger maps affected only by the confusion noise using the simulated infrared extragalactic sky (SIDES) semi-empirical simulation. We then estimated the confusion limit in these maps and extracted the sources using a basic blind extractor. By comparing the input galaxy catalog and the extracted source catalog, we derived various performance metrics as completeness, purity, and the accuracy of various measurements. In intensity, we predict that the confusion limit increases rapidly with increasing wavelength. The confusion limit in polarization is more than 100x lower. The measured flux density is dominated by the brightest galaxy in the beam, but other objects also contribute at longer wavelength (~30% at 235 $μ$m). We also show that galaxy clustering has a mild impact on confusion in intensity (up to 25%), while it is negligible in polarization. In intensity, a basic blind extraction will be sufficient to detect galaxies at the knee of the luminosity function up to z~3 and 10$^{11}$ M$_\odot$ main-sequence galaxies up to z~5. In polarization for a conservative sensitivity, we expect ~8 000 detections up to z=2.5 opening a totally new window on the high-z dust polarization. Finally, we show that intensity surveys at short wavelength and polarization surveys at long wavelength tend to reach confusion at similar depth. There is thus a strong synergy. △ Less

Submitted 28 April, 2024; v1 submitted 5 April, 2024; originally announced April 2024.

Comments: 21 pages, 15 figures, 4 tables, submitted to A&A, minor update on the names of the sensitivity scenarii (for consistency with other papers)

Atacama Large Aperture Submillimeter Telescope (AtLAST) science: Gas and dust in nearby galaxies

Authors: Daizhong Liu, Amelie Saintonge, Caroline Bot, Francisca Kemper, Enrique Lopez-Rodriguez, Matthew W. L. Smith, Thomas Stanke, Paola Andreani, Alessandro Boselli, Claudia Cicone, Timothy A. Davis, Bendix Hagedorn, Akhil Lasrado, Ann Mao, Serena Viti, Mark Booth, Pamela Klaassen, Tony Mroczkowski, Frank Bigiel, Melanie Chevance, Martin A. Cordiner, Luca Di Mascolo, Doug Johnstone, Minju M. Lee, Thomas Maccarone , et al. (3 additional authors not shown)

Abstract: Understanding the physical processes that regulate star formation and galaxy evolution are major areas of activity in modern astrophysics. Nearby galaxies offer unique opportunities to inspect interstellar medium (ISM), star formation (SF), radiative, dynamic and magnetic physics in great detail from sub-galactic (kpc) scales to sub-cloud (sub-pc) scales, from quiescent galaxies to starbursts, and… ▽ More Understanding the physical processes that regulate star formation and galaxy evolution are major areas of activity in modern astrophysics. Nearby galaxies offer unique opportunities to inspect interstellar medium (ISM), star formation (SF), radiative, dynamic and magnetic physics in great detail from sub-galactic (kpc) scales to sub-cloud (sub-pc) scales, from quiescent galaxies to starbursts, and from field galaxies to overdensities. In this case study, we discuss the major breakthroughs in this area of research that will be enabled by the Atacama Large Aperture Submillimeter Telescope (AtLAST), a proposed 50-m single-dish submillimeter telescope. The new discovery space of AtLAST comes from its exceptional sensitivity, in particular to extended low surface brightness emission, a very large 2 degree field of view, and correspondingly high mapping efficiency. This paper focuses on four themes which will particularly benefit from AtLAST: 1) the LMC and SMC, 2) extragalactic magnetic fields, 3) the physics and chemistry of the interstellar medium, and 4) star formation and galaxy evolution. With ~1000-2000h surveys each, AtLAST could deliver deep dust continuum maps of the entire LMC and SMC fields at parsec-scale resolution, high-resolution maps of the magnetic field structure, gas density, temperature and composition of the dense and diffuse ISM in ~100 nearby galaxies, as well as the first large-scale blind CO survey in the nearby Universe, delivering molecular gas masses for up to 10^6 galaxies (3 orders of magnitude more than current samples). Through such observing campaigns, AtLAST will have a profound impact on our understanding of the baryon cycle and star formation across a wide range of environments. △ Less

Submitted 2 March, 2024; originally announced March 2024.

Comments: 29 pages, 11 figues, submitted to Open Research Europe as part of the AtLAST collection: https://open-research-europe.ec.europa.eu/collections/atlast/about

Probing the dust grain alignment mechanisms in spiral galaxies with M51 as the case study

Authors: Enrique Lopez-Rodriguez, Le Ngoc Tram

Abstract: Magnetic fields (B-fields) in galaxies have recently been traced using far-infrared and sub-mm polarimetric observations with SOFIA, JCMT, and ALMA. The main assumption is that dust grains are magnetically aligned with the local B-field in the interstellar medium (ISM). However, the range of physical conditions of the ISM, dust grain sizes, and B-field strengths in galaxies where this assumption i… ▽ More Magnetic fields (B-fields) in galaxies have recently been traced using far-infrared and sub-mm polarimetric observations with SOFIA, JCMT, and ALMA. The main assumption is that dust grains are magnetically aligned with the local B-field in the interstellar medium (ISM). However, the range of physical conditions of the ISM, dust grain sizes, and B-field strengths in galaxies where this assumption is valid has not been characterized yet. Here, we use the well-studied spiral galaxy M51 as a case study. We find that the timescale for the alignment mechanism arising from magnetically aligned dust grains (B-RAT) dominates over other alignment mechanisms, including radiative precession (k-RAT) and mechanical alignment (v-MAT), as well as the randomization effect (gas damping). We estimate the sizes of the aligned dust grain to be in the range of 0.009-0.182 $μ$m and 0.019-0.452 $μ$m for arms and inter-arms, respectively. We show that the difference in the polarization fraction between arms and interarms can arise from the enhancement of small dust grain sizes in the arms as an effect of the grain alignment disruption (RAT-D). We argue that the RAT-D mechanism needs to have additional effects, e.g., intrinsic variations of the B-field structure and turbulence, in the galaxy's components to fully explain the polarization fraction variations within the arms and inter-arms. △ Less

Submitted 28 February, 2024; originally announced February 2024.

Comments: 12 pages, 3 figures. Submitted to ApJ (comments are welcome)

Extragalactic Magnetism with SOFIA (SALSA Legacy Program). VII. A Tomographic View of Far-infrared and Radio Polarimetric Observations through MHD Simulations of Galaxies

Authors: Sergio Martin-Alvarez, Enrique Lopez-Rodriguez, Tara Dacunha, Susan E. Clark, Alejandro S. Borlaff, Rainer Beck, Francisco Rodríguez Montero, S. Lyla Jung, Julien Devriendt, Adrianne Slyz, Julia Roman-Duval, Evangelia Ntormousi, Mehrnoosh Tahani, Kandaswamy Subramanian, Daniel A. Dale, Pamela M. Marcum, Konstantinos Tassis, Ignacio del Moral-Castro, Le Ngoc Tram, Matt J. Jarvis

Abstract: The structure of magnetic fields in galaxies remains poorly constrained, despite the importance of magnetism in the evolution of galaxies. Radio synchrotron and far-infrared (FIR) polarization and polarimetric observations are the best methods to measure galactic scale properties of magnetic fields in galaxies beyond the Milky Way. We use synthetic polarimetric observations of a simulated galaxy t… ▽ More The structure of magnetic fields in galaxies remains poorly constrained, despite the importance of magnetism in the evolution of galaxies. Radio synchrotron and far-infrared (FIR) polarization and polarimetric observations are the best methods to measure galactic scale properties of magnetic fields in galaxies beyond the Milky Way. We use synthetic polarimetric observations of a simulated galaxy to identify and quantify the regions, scales, and interstellar medium (ISM) phases probed at FIR and radio wavelengths. Our studied suite of magnetohydrodynamical cosmological zoom-in simulations features high-resolutions (10 pc full-cell size) and multiple magnetization models. Our synthetic observations have a striking resemblance to those of observed galaxies. We find that the total and polarized radio emission extends to approximately double the altitude above the galactic disk (half-intensity disk thickness of $h_\text{I radio} \sim h_\text{PI radio} = 0.23 \pm 0.03$ kpc) relative to the total FIR and polarized emission that are concentrated in the disk midplane ($h_\text{I FIR} \sim h_\text{PI FIR} = 0.11 \pm 0.01$ kpc). Radio emission traces magnetic fields at scales of $\gtrsim 300$ pc, whereas FIR emission probes magnetic fields at the smallest scales of our simulations. These scales are comparable to our spatial resolution and well below the spatial resolution ($<300$ pc) of existing FIR polarimetric measurements. Finally, we confirm that synchrotron emission traces a combination of the warm neutral and cold neutral gas phases, whereas FIR emission follows the densest gas in the cold neutral phase in the simulation. These results are independent of the ISM magnetic field strength. The complementarity we measure between radio and FIR wavelengths motivates future multiwavelength polarimetric observations to advance our knowledge of extragalactic magnetism. △ Less

Submitted 24 March, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Comments: Accepted for publication in the ApJ. Replaced to match final version. 35 pages, 17 figures, 4 tables

Cygnus A Obscuring Torus: Ionized, Atomic or Molecular?

Authors: Atharva Bagul, Patrick Ogle, Robert Antonucci, Philip Maloney, Enrique Lopez-Rodriguez

Abstract: The prototypical powerful FR \Romannum{2} radio galaxy Cygnus A fits extremely well into the quasar/radio galaxy unified model: high polarization with an angle almost perpendicular to the radio jet and polarized flux showing broad permitted lines. It has been claimed that ionized gas in the torus reveals a very clear torus shape via Bremmstrahlung emission. We rule out the later with an energetic… ▽ More The prototypical powerful FR \Romannum{2} radio galaxy Cygnus A fits extremely well into the quasar/radio galaxy unified model: high polarization with an angle almost perpendicular to the radio jet and polarized flux showing broad permitted lines. It has been claimed that ionized gas in the torus reveals a very clear torus shape via Bremmstrahlung emission. We rule out the later with an energetic argument, and we constrain the molecular and atomic gas properties with existing observations. The atomic absorption against the core has been shown to match the X-ray column only if the spin temperature is an implausible $T_{\rm s} = 1\times 10^6$ K. This points to a molecular medium for the X-ray column $\log(N_{\rm H} ~[\rm{cm^{-2}}]) \sim 23.5$. Yet not low-J CO absorption is detected to sensitive limits. The non-detection is surprising given that this powerful radio galaxy hosts a luminous, dust-obscured active nucleus and copious warm molecular hydrogen. These conditions suggest a detectable level of emission. Furthermore, the torus X-ray column density suggests detectable absorption. We explore various possibilities to explain the lack of a signature from warm CO (200-250K). Specifically, that the radiative excitation by the radio core renders low-J CO absorption below current sensitivities, and that high-J levels are well populated and conducive to producing absorption. We test this hypothesis using archival \textit{Hershel}/SPIRE FTS observations of Cygnus A of high-J CO lines ($14 \geq J \geq 4$ transitions). Still high-J CO lines are not detected. We suggest that ALMA observations near its high frequency limit can be critical to obtain the signature of molecular line of the torus of Cygnus A. △ Less

Submitted 30 October, 2023; originally announced October 2023.

The Galaxy Activity, Torus, and Outflow Survey (GATOS). III: Revealing the inner icy structure in local AGN

Authors: I. García-Bernete, A. Alonso-Herrero, D. Rigopoulou, M. Pereira-Santaella, T. Shimizu, R. Davies, F. R. Donnan, P. F. Roche, O. González-Martín, C. Ramos Almeida, E. Bellocchi, P. Boorman, F. Combes, A. Efstathiou, D. Esparza-Arredondo, S. García-Burillo, E. González-Alfonso, E. K. S. Hicks, S. Hönig, A. Labiano, N. A. Levenson, E. López-Rodríguez, C. Ricci, C. Packham, D. Rouan , et al. (2 additional authors not shown)

Abstract: We use JWST/MIRI MRS spectroscopy of a sample of six local obscured type 1.9/2 active galactic nuclei (AGN) to compare their nuclear mid-IR absorption bands with the level of nuclear obscuration traced by X-rays. This study is the first to use sub-arcsecond angular resolution data of local obscured AGN to investigate the nuclear mid-IR absorption bands with a wide wavelength coverage (4.9-28.1… ▽ More We use JWST/MIRI MRS spectroscopy of a sample of six local obscured type 1.9/2 active galactic nuclei (AGN) to compare their nuclear mid-IR absorption bands with the level of nuclear obscuration traced by X-rays. This study is the first to use sub-arcsecond angular resolution data of local obscured AGN to investigate the nuclear mid-IR absorption bands with a wide wavelength coverage (4.9-28.1 $μ$m). All the nuclei show the 9.7 $μ$m silicate band in absorption. We compare the strength of the 9.7 and 18 $μ$m silicate features with torus model predictions. The observed silicate features are generally well explained by clumpy and smooth torus models. We report the detection of the 6 $μ$m dirty water ice band (i.e., a mix of water and other molecules such as CO and CO$_2$) at sub-arcsecond scales ($\sim$0.26 arcsec at 6 $μ$m; inner $\sim$50 pc) in a sample of local AGN with different levels of nuclear obscuration in the range log N$_{\rm H}^{\rm X-Ray}$(cm$^{-2}$)$\sim22-25$. We find a good correlation between the 6 $μ$m water ice optical depths and N$_{\rm H}^{\rm X-Ray}$. This result indicates that the water ice absorption might be a reliable tracer of the nuclear intrinsic obscuration in AGN. The weak water ice absorption in less obscured AGN (log N$_H^{X-ray}$ (cm$^{-2}$)$\lesssim$23.0 cm$^{-2}$) might be related to the hotter dust temperature ($>$T$_{sub}^{H_2O}\sim$110 K) expected to be reached in the outer layers of the torus due to their more inhomogeneous medium. Our results suggest it might be necessary to include the molecular content, such as, H$_2$O, aliphatic hydrocarbons (CH-) and more complex PAH molecules in torus models to better constrain key parameters such as the torus covering factor (i.e. nuclear obscuration). △ Less

Submitted 4 December, 2023; v1 submitted 13 October, 2023; originally announced October 2023.

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

An X-rays-to-radio investigation of the nuclear polarization from the radio-galaxy Centaurus A

Authors: Frédéric Marin, Thibault Barnouin, Steven R. Ehlert, Abel Lawrence Peirson, Enrique Lopez-Rodriguez, Maria Petropoulou, Kinwah Wu, Iván Martí-Vidal

Abstract: Centaurus A is one of the closest radio-galaxies to Earth. Its proximity allowed us to extensively study its active galactic nucleus but the core emission mechanism remains elusive because of local strong dust and gas obscuration. The capability of polarimetry to shave-off contaminating emission has been exploited without success in the near-infrared by previous studies but the very recent measure… ▽ More Centaurus A is one of the closest radio-galaxies to Earth. Its proximity allowed us to extensively study its active galactic nucleus but the core emission mechanism remains elusive because of local strong dust and gas obscuration. The capability of polarimetry to shave-off contaminating emission has been exploited without success in the near-infrared by previous studies but the very recent measurement of the 2 - 8 keV polarization by the Imaging X-ray Polarimetry Explorer (IXPE) brought the question back to the fore. To determine what is the prevalent photon generation mechanism to the multi-wavelength emission from the core of Centaurus A, we retrieved from the archives the panchromatic polarization measurements of the central compact component. We built the total and polarized flux spectral energy distributions of the core and demonstrated that synchrotron self-Compton models nicely fit the polarized flux from the radio to the X-ray band. The linear polarization of the synchrotron continuum is perpendicular to the jet radio axis from the optical to the radio band, and parallel to it at higher energies. The observed smooth rotation of the polarization angle in the ultraviolet band is attributed to synchrotron emission from regions that are getting closer to the particle acceleration site, where the orientation of the jet's magnetic fields become perpendicular to the jet axis. This phenomenon support the shock acceleration mechanism for particle acceleration in Centaurus A, in line with IXPE observations of several high-synchrotron peak blazars. △ Less

Submitted 6 October, 2023; originally announced October 2023.

Comments: 10 pages, 6 figures, accepted for publication in MNRAS

MSC Class: 85-06 ACM Class: J.2.3; J.2.9

Forgotten treasures in the HST/FOC UV imaging polarimetric archives of active galactic nuclei. I. Pipeline and benchmarking against NGC~1068 and exploring IC~5063

Authors: Thibault Barnouin, Frédéric Marin, Enrique Lopez-Rodriguez, Léo Huber, Makoto Kishimoto

Abstract: Over its 13 years of operation (1990 -- 2002), the Faint Object Camera (FOC) on board the Hubble Space Telescope (HST) observed 26 individual active galactic nuclei (AGNs) in ultraviolet (UV) imaging polarimetry. However, not all of the observations have been reduced and analyzed or set within a standardized framework. We plan to reduce and analyze the AGN observations that have been neglected in… ▽ More Over its 13 years of operation (1990 -- 2002), the Faint Object Camera (FOC) on board the Hubble Space Telescope (HST) observed 26 individual active galactic nuclei (AGNs) in ultraviolet (UV) imaging polarimetry. However, not all of the observations have been reduced and analyzed or set within a standardized framework. We plan to reduce and analyze the AGN observations that have been neglected in the FOC archives using a consistent, novel, and open-access reduction pipeline of our own. We then extend the method to the full AGN sample, thus leading to potential discoveries in the near future. We developed a new pipeline in Python that will be able to reduce all the FOC observations in imaging polarimetry in a homogeneous way. Most of the previously published reduced observations are dispersed throughout the literature, with the range of different analyses and approaches making it difficult to fully interpret the FOC AGN sample. By standardizing the method, we have enabled a coherent comparison among the different observational sets. In this first paper of a series exploring the full HST/FOC AGN sample, we present an exhaustively detailed account of how to properly reduce the observational data. Current progress in data-analysis is implemented in and has provided state-of-the-art UV polarimetric maps. We compare our new maps to the benchmark AGN case of NGC~1068 and successfully reproduce the main results previously published, while pushing the polarimetric exploration of this AGN futher, thanks to a finer resolution and a higher signal-to-noise ratio (S/N) than previously reported. We also present, for the first time, an optical polarimetric map of the radio-loud AGN IC~5063 and we examine the complex interactions between the AGN outflows and the surrounding interstellar medium (ISM). △ Less

Submitted 5 September, 2023; originally announced September 2023.

Comments: 17 pages, 19 figures, to be published in Astronomy & Astrophysics

MSC Class: 85-04 (Primary); 85-08; 85A25 (Secondary) ACM Class: J.2; I.4.1; I.4.7; I.6.4

Journal ref: A&A 678, A143 (2023)

Polarized thermal emission from dust in a galaxy at redshift 2.6

Authors: J. E. Geach, E. Lopez-Rodriguez, M. J. Doherty, Jianhang Chen, R. J. Ivison, G. J. Bendo, S. Dye, K. E. K. Coppin

Abstract: Magnetic fields are fundamental to the evolution of galaxies, playing a key role in the astrophysics of the interstellar medium and star formation. Large-scale ordered magnetic fields have been mapped in the Milky Way and nearby galaxies, but it is not known how early in the Universe such structures form. Here we report the detection of linearly polarized thermal emission from dust grains in a str… ▽ More Magnetic fields are fundamental to the evolution of galaxies, playing a key role in the astrophysics of the interstellar medium and star formation. Large-scale ordered magnetic fields have been mapped in the Milky Way and nearby galaxies, but it is not known how early in the Universe such structures form. Here we report the detection of linearly polarized thermal emission from dust grains in a strongly lensed, intrinsically luminous galaxy that is forming stars at a rate more than a thousand times that of the Milky Way at redshift 2.6, within 2.5 Gyr of the Big Bang. The polarized emission arises from the alignment of dust grains with the local magnetic field. The median polarization fraction is of order one per cent, similar to nearby spiral galaxies. Our observations support the presence of a 5 kiloparsec-scale ordered magnetic field with a strength of around 500uG or lower, orientated parallel to the molecular gas disk. This confirms that such structures can be rapidly formed in galaxies, early in cosmic history. △ Less

Submitted 5 September, 2023; originally announced September 2023.

Comments: Published in Nature. Online version available at https://www.nature.com/articles/s41586-023-06346-4

The magnetic fields of starburst galaxies. I. Identification and characterization of the thermal polarization in the galactic disk and outflow

Authors: Enrique Lopez-Rodriguez

Abstract: Far-infrared polarized emission by means of magnetically aligned dust grains is an excellent tracer of the magnetic fields (B-fields) in the cold phase of the galactic outflows in starburst galaxies. We present a comprehensive study of the B-fields in three nearby ($3.5$-$17.2$ Mpc) starbursts (M82, NGC 253, and NGC 2146) at $5$ pc-$1.5$ kpc resolutions using publicly available $53$-$890$ $μ$m ima… ▽ More Far-infrared polarized emission by means of magnetically aligned dust grains is an excellent tracer of the magnetic fields (B-fields) in the cold phase of the galactic outflows in starburst galaxies. We present a comprehensive study of the B-fields in three nearby ($3.5$-$17.2$ Mpc) starbursts (M82, NGC 253, and NGC 2146) at $5$ pc-$1.5$ kpc resolutions using publicly available $53$-$890$ $μ$m imaging polarimetric observations with SOFIA/HAWC+, JCMT/POL-2, and ALMA. We find that the polarized spectral energy distributions (SEDs) of the full galaxies are dominated by the polarized SEDs of the outflows with dust temperatures of $T_{\rm{d,outflow}}^{\rm{PI}}\sim45$ K and emissive index of $β_{\rm{outflow}}^{\rm{PI}}\sim2.3$. The disks are characterized by low $T_{\rm{d,disk}}^{\rm{PI}}=[24,31]$ K and $β_{\rm{disk}}^{\rm{PI}}\sim1$. We show that disk- and outflow-dominated galaxies can be better distinguished by using polarized SEDs instead of total SEDs. We compute the $53$-$850$ $μ$m polarization spectrum of the disk and outflow and find that dust models of the diffuse ISM can reproduce the fairly constant polarization spectrum of the disk, $\langle P_{\rm{disk}} \rangle=1.2\pm0.5$%. The dust models of heterogenous clouds and two temperature components are required to explain the polarization spectrum of the outflow ($2$-$4$% at $53$ $μ$m, $\sim1$% at $850$ $μ$m, and a minimum within $89$-$154$ $μ$m). We conclude that the polarized dust grains in the outflow arise from a dust population with higher dust temperature and emissivities than those from the total flux. The B-fields of the outflows have maximum extensions within $89$-$214$ $μ$m reaching heights of $\sim4$ kpc, and flatter polarized fluxes than total fluxes. The extension of the B-field permeating the circumgalactic medium increases with increasing the star formation rate. △ Less

Submitted 16 June, 2023; originally announced June 2023.

Comments: 21 pages, 14 figures, 4 tables. Submitted to ApJ with positive review by the anonymous referee. Comments are welcome

The origin of dust polarization in the Orion Bar

Authors: Valentin J. M. Le Gouellec, B-G Andersson, Archana Soam, Thiébaut Schirmer, Joseph M. Michail, Enrique Lopez-Rodriguez, Sophia Flores, David T. Chuss, John E. Vaillancourt, Thiem Hoang, Alex Lazarian

Abstract: The linear polarization of thermal dust emission provides a powerful tool to probe interstellar and circumstellar magnetic fields, because aspherical grains tend to align themselves with magnetic field lines. While the Radiative Alignment Torque (RAT) mechanism provides a theoretical framework to this phenomenon, some aspects of this alignment mechanism still need to be quantitatively tested. One… ▽ More The linear polarization of thermal dust emission provides a powerful tool to probe interstellar and circumstellar magnetic fields, because aspherical grains tend to align themselves with magnetic field lines. While the Radiative Alignment Torque (RAT) mechanism provides a theoretical framework to this phenomenon, some aspects of this alignment mechanism still need to be quantitatively tested. One such aspect is the possibility that the reference alignment direction changes from the magnetic field ("B-RAT") to the radiation field k-vector ("k-RAT") in areas of strong radiation fields. We investigate this transition toward the Orion Bar PDR, using multi-wavelength SOFIA HAWC+ dust polarization observations. The polarization angle maps show that the radiation field direction is on average not the preferred grain alignment axis. We constrain the grain sizes for which the transition from B-RAT to k-RAT occur in the Orion Bar (grains > 0.1 μm toward the most irradiated locations), and explore the radiatively driven rotational disruption that may take place in the high-radiation environment of the Bar for large grains. While the grains susceptible to rotational disruption should be in supra-thermal rotation and aligned with the magnetic field, k-RAT aligned grains would rotate at thermal velocities. We find that the grain size at which the alignment shifts from B-RAT to k-RAT corresponds to grains too large to survive the rotational disruption. Therefore, we expect a large fraction of grains to be aligned at supra-thermal rotation with the magnetic field, and potentially be subject to rotational disruption depending on their tensile strength. △ Less

Submitted 3 May, 2023; originally announced May 2023.

The Strength of the Sheared Magnetic Field in the Galactic's Circum-Nuclear Disk

Authors: J. A. Guerra, E. Lopez-Rodriguez, D. T. Chuss, N. O. Butterfield, J. T. Schmelz

Abstract: Recent high-resolution 53-$μ$m polarimetric observations from SOFIA/HAWC+ have revealed the inferred plane-of-the-sky magnetic field (B-field) orientation in the Galactic center's Circum-Nuclear Disk (CND). The B-field is mostly aligned with the steamers of ionized material falling onto Sgr A* at large, differential velocities (shear). In such conditions, estimating the B-field strength with the `… ▽ More Recent high-resolution 53-$μ$m polarimetric observations from SOFIA/HAWC+ have revealed the inferred plane-of-the-sky magnetic field (B-field) orientation in the Galactic center's Circum-Nuclear Disk (CND). The B-field is mostly aligned with the steamers of ionized material falling onto Sgr A* at large, differential velocities (shear). In such conditions, estimating the B-field strength with the ``classical" Davis-Chandrasekhar-Fermi (DCF) method does not provide accurate results. We derive a ``modified'' DCF method by solving the ideal MHD equations from first principles considering the effects of a large-scale, shear flow on the propagation of a fast magnetosonic wave. In the context of the DCF approximation, both the value of the shear and its Laplacian affect the inferred B-field strength. Using synthetic polarization data from MHD simulations for a medium dominated by shear flows, we find that the ``classical'' DCF determines B-field strengths only within $>50$\% of the true value where the ``modified" DCF results are improved significantly ($\sim$3-22\%). Applying our ``modified'' DCF method to the CND revealed B-field strengths of 1 - 16 mG in the northern arm, 1 - 13 mG in the eastern arm, and 3 - 27 mG in the western arm at spatial scales $\lesssim1$ pc, with median values of $5.1\pm0.8$, $4.0\pm1.2$, and $8.5\pm2.3$ mG, respectively. The balance between turbulent gas energy (kinetic plus hydrostatic) and turbulent magnetic energy densities suggest that, along the magnetic-field-flow direction, magnetic effects become less dominant as the shear flow increases and weakens the B-field via magnetic convection. Our results indicate that the transition from magnetically to gravitationally dominated accretion of material onto Sgr A* starts at distances $\sim$ 1 pc. △ Less

Submitted 10 June, 2023; v1 submitted 13 April, 2023; originally announced April 2023.

Comments: Accepted for publication in the Astronomical Journal

Extragalactic magnetism with SOFIA (SALSA Legacy Program) -- V: First results on the magnetic field orientation of galaxies

Authors: Alejandro S. Borlaff, Enrique Lopez-Rodriguez, Rainer Beck, Susan E. Clark, Evangelia Ntormousi, Konstantinos Tassis, Sergio Martin-Alvarez, Mehrnoosh Tahani, Daniel A. Dale, Ignacio del Moral Castro, Julia Roman-Duval, Pamela M. Marcum, John E. Beckman, Kandaswamy Subramanian, Sarah Eftekharzadeh, Leslie Proudfit

Abstract: We present the analysis of the magnetic field ($B$-field) structure of galaxies measured with far-infrared (FIR) and radio (3 and 6 cm) polarimetric observations. We use the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA) of 14 nearby ($<20$ Mpc) galaxies with resolved (5 arcsec-18 arcsec; $90$ pc--$1$ kpc) imaging polarimetric observations using HAWC+/SOFIA from… ▽ More We present the analysis of the magnetic field ($B$-field) structure of galaxies measured with far-infrared (FIR) and radio (3 and 6 cm) polarimetric observations. We use the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA) of 14 nearby ($<20$ Mpc) galaxies with resolved (5 arcsec-18 arcsec; $90$ pc--$1$ kpc) imaging polarimetric observations using HAWC+/SOFIA from $53$ to $214$ \um. We compute the magnetic pitch angle ($Ψ_{B}$) profiles as a function of the galactrocentric radius. We introduce a new magnetic alignment parameter ($ζ$) to estimate the disordered-to-ordered $B$-field ratio in spiral $B$-fields. We find FIR and radio wavelengths to not generally trace the same $B$-field morphology in galaxies. The $Ψ_{B}$ profiles tend to be more ordered with galactocentric radius in radio ($ζ_{\rm{6cm}} = 0.93\pm0.03$) than in FIR ($ζ_{\rm{154μm}} = 0.84\pm0.14$). For spiral galaxies, FIR $B$-fields are $2-75$\% more turbulent than the radio $B$-fields. For starburst galaxies, we find that FIR polarization is a better tracer of the $B$-fields along the galactic outflows than radio polarization. Our results suggest that the $B$-fields associated with dense, dusty, turbulent star-forming regions, those traced at FIR, are less ordered than warmer, less-dense regions, those traced at radio, of the interstellar medium. The FIR $B$-fields seem to be more sensitive to the activity of the star-forming regions and the morphology of the molecular clouds within a vertical height of few hundred pc in the disk of spiral galaxies than the radio $B$-fields. △ Less

Submitted 13 June, 2023; v1 submitted 23 March, 2023; originally announced March 2023.

Comments: 28 pages, 14 figures. Accepted for publication in The Astrophysical Journal

The structure of magnetic fields in spiral galaxies: a radio and far-infrared polarimetric analysis

Authors: William Jeffrey Surgent, Enrique Lopez-Rodriguez, Susan E. Clark

Abstract: We propose and apply a method to quantify the morphology of the large-scale ordered magnetic fields (B-fields) in galaxies. This method is adapted from the analysis of Event Horizon Telescope polarization data. We compute a linear decomposition of the azimuthal modes of the polarization field in radial galactocentric bins. We apply this approach to five low-inclination spiral galaxies with both fa… ▽ More We propose and apply a method to quantify the morphology of the large-scale ordered magnetic fields (B-fields) in galaxies. This method is adapted from the analysis of Event Horizon Telescope polarization data. We compute a linear decomposition of the azimuthal modes of the polarization field in radial galactocentric bins. We apply this approach to five low-inclination spiral galaxies with both far-infrared (FIR: 154 $μ$m) dust polarimetric observations taken from the Survey of ExtragALactic magnetiSm with SOFIA (SALSA) and radio (6 cm) synchrotron polarization observations. We find that the main contribution to the B-field structure of these spiral galaxies comes from the $m=2$ and $m=0$ modes at FIR wavelengths and the $m=2$ mode at radio wavelengths. The $m=2$ mode has a spiral structure and is directly related to the magnetic pitch angle, while $m=0$ has a constant B-field orientation. The FIR data tend to have a higher relative contribution from other modes than the radio data. The extreme case is NGC 6946: all modes contribute similarly in the FIR, while $m=2$ still dominates in the radio. The average magnetic pitch angle in the FIR data is smaller and has greater angular dispersion than in the radio, indicating that the B-fields in the disk midplane traced by FIR dust polarization are more tightly wound and more chaotic than the B-field structure in the radio, which probes a larger volume. We argue that our approach is more flexible and model-independent than standard techniques, while still producing consistent results where directly comparable. △ Less

Submitted 7 June, 2023; v1 submitted 14 February, 2023; originally announced February 2023.

Comments: 14 pages, 8 figures, 3 tables. Manuscript submitted to ApJ (Comments are welcome)

Dissecting the active galactic nucleus in Circinus -- III. VLT/FORS2 polarimetry confirms dusty cone illuminated by a tilted accretion disc

Authors: Marko Stalevski, Santiago González-Gaitán, Đorđe Savić, Makoto Kishimoto, Ana Mourão, Enrique Lopez-Rodriguez, Daniel Asmus

Abstract: We present polarimetric maps of the Circinus galaxy nucleus in the $BVRI$ bands, obtained with VLT/FORS2. Circinus is the closest Seyfert 2 galaxy and harbours an archetypal obscured active galactic nucleus (AGN). Recent high angular resolution imaging revealed that a major fraction of its mid-infrared (MIR) emission is coming from the polar region. Previously, we demonstrated that these observati… ▽ More We present polarimetric maps of the Circinus galaxy nucleus in the $BVRI$ bands, obtained with VLT/FORS2. Circinus is the closest Seyfert 2 galaxy and harbours an archetypal obscured active galactic nucleus (AGN). Recent high angular resolution imaging revealed that a major fraction of its mid-infrared (MIR) emission is coming from the polar region. Previously, we demonstrated that these observations are consistent with a model of a compact dusty disc and a hyperboloid shell, resembling a hollow cone on larger scales. Here we focus on the AGN core, up to 40 pc from the central engine, and compare the observations to the radiative transfer models. Polarization maps reveal a conical structure, coinciding with the ionization cone. The wavelength-dependence of the polarization degree indicates that scattering on dust grains is producing polarization. The observed polarization degree ($\approx1-3\%$) is lower than predicted by the models; however, this is only a lower limit, since stellar emission dominates the total flux in the optical. The observed polarization angle ($\approx30$ degrees) is reproduced by the model of a dusty disc with a hollow cone that is illuminated by a tilted anisotropic central source. An accretion disc aligned with the ionization cone axis, and alternative dust geometries, such as a paraboloid shell, or a torus enveloped by ambient dust, are inconsistent with the data. We conclude that the optical polarimetric imaging supports earlier evidence for the presence of dust in the polar region, tentatively associated with dusty outflows. △ Less

Submitted 20 December, 2022; originally announced December 2022.

Comments: MNRAS (accepted)

The Pandora project. I: the impact of radiation and cosmic rays on baryonic and dark matter properties of dwarf galaxies

Authors: Sergio Martin-Alvarez, Debora Sijacki, Martin G. Haehnelt, Marion Farcy, Yohan Dubois, Vasily Belokurov, Joakim Rosdahl, Enrique Lopez-Rodriguez

Abstract: Enshrouded in several well-known controversies, dwarf galaxies have been extensively studied to learn about the underlying cosmology, notwithstanding that physical processes regulating their properties are poorly understood. To shed light on these processes, we introduce the Pandora suite of 17 high-resolution (3.5 parsec half-cell side) dwarf galaxy formation cosmological simulations. Commencing… ▽ More Enshrouded in several well-known controversies, dwarf galaxies have been extensively studied to learn about the underlying cosmology, notwithstanding that physical processes regulating their properties are poorly understood. To shed light on these processes, we introduce the Pandora suite of 17 high-resolution (3.5 parsec half-cell side) dwarf galaxy formation cosmological simulations. Commencing with thermo-turbulent star formation and mechanical supernova feedback, we gradually increase the complexity of physics incorporated leading to full-physics models combining magnetism, on-the-fly radiative transfer and the corresponding stellar photoheating, and SN-accelerated cosmic rays. We investigate combinations of these processes, comparing them with observations to constrain what are the main mechanisms determining dwarf galaxy properties. We find hydrodynamical `SN feedback-only' simulations struggle to produce realistic dwarf galaxies, leading either to overquenched or too centrally concentrated, dispersion dominated systems when compared to observed field dwarfs. Accounting for radiation with cosmic rays results in extended and rotationally-supported systems. Spatially `distributed' feedback leads to realistic stellar and HI masses as well as kinematics. Furthermore, resolved kinematic maps of our full-physics models predict kinematically distinct clumps and kinematic misalignments of stars, HI and HII after star formation events. Episodic star formation combined with its associated feedback induces more core-like dark matter central profiles, which our `SN feedback-only' models struggle to achieve. Our results demonstrate the complexity of physical processes required to capture realistic dwarf galaxy properties, making tangible predictions for integral field unit surveys, radio synchrotron emission, and for galaxy and multi-phase interstellar medium properties that JWST will probe. △ Less

Submitted 25 September, 2023; v1 submitted 16 November, 2022; originally announced November 2022.

Comments: Replaced to match the MNRAS published version

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 525, Issue 3, November 2023, Pages 3806-3830

Extragalactic magnetism with SOFIA (SALSA Legacy Program). VI. The magnetic fields in the multi-phase interstellar medium of the Antennae galaxies

Authors: Enrique Lopez-Rodriguez, Alejandro S. Borlaff, Rainer Beck, William T. Reach, Sui Ann Mao, Evangelia Ntormousi, Konstantinos Tassis, Sergio Martin-Alvarez, Susan E. Clark, Daniel A. Dale, Ignacio del Moral-Castro

Abstract: Mergers are thought to be a fundamental channel for galaxy growth, perturbing the gas dynamics and the magnetic fields (B-fields) in the interstellar medium (ISM). However, the mechanisms that amplify and dissipate B-fields during a merger remain unclear. We characterize the morphology of the ordered B-fields in the multi-phase ISM of the closest merger of two spiral galaxies, the Antennae galaxie… ▽ More Mergers are thought to be a fundamental channel for galaxy growth, perturbing the gas dynamics and the magnetic fields (B-fields) in the interstellar medium (ISM). However, the mechanisms that amplify and dissipate B-fields during a merger remain unclear. We characterize the morphology of the ordered B-fields in the multi-phase ISM of the closest merger of two spiral galaxies, the Antennae galaxies. We compare the inferred B-fields using $154~μ$m thermal dust and $11$ cm radio synchrotron emission polarimetric observations. We find that the $154~μ$m B-fields are more ordered across the Antennae galaxies than the $11$ cm B-fields. The turbulent-to-ordered $154~μ$m B-field increases at the galaxy cores and star-forming regions. The relic spiral arm has an ordered spiral $154~μ$m B-field, while the $11$ cm B-field is radial. The $154~μ$m B-field may be dominated by turbulent dynamos with high $^{12}$CO(1-0) velocity dispersion driven by star-forming regions, while the $11$ cm B-field is cospatial with high HI velocity dispersion driven by galaxy interaction. This result shows the dissociation between the warm gas mainly disturbed by the merger, and the dense gas still following the dynamics of the relic spiral arm. We find a $\sim8.9$ kpc scale ordered B-field connecting the two galaxies. The base of the tidal tail is cospatial with the HI and $^{12}$CO(1-0) emission and has compressed and/or sheared $154~μ$m and $11$ cm B-fields driven by the merger. We suggest that amplify B-fields, with respect to the rest of the system and other spiral galaxies, may be supporting the gas flow between both galaxies and the tidal tail. △ Less

Submitted 9 December, 2022; v1 submitted 31 October, 2022; originally announced November 2022.

Comments: 11 pages, 5 figures, Accepted for publication in ApJ Letters

Torus and polar dust dependence on AGN properties

Authors: I. García-Bernete, O. González-Martín, C. Ramos Almeida, A. Alonso-Herrero, M. Martínez-Paredes, M. J. Ward, P. F. Roche, J. A. Acosta-Pulido, E. López-Rodríguez, D. Rigopoulou, D. Esparza-Arredondo

Abstract: We present a statistical analysis of the properties of the obscuring material around active galactic nuclei (AGN). This study represents the first of its kind for an ultra-hard X-ray (14-195keV; Swift/BAT) volume-limited (DL<40 Mpc) sample of 24 Seyfert (Sy) galaxies (BCS40 sample) using high angular resolution infrared data and various torus models: smooth, clumpy and two-phase torus models and c… ▽ More We present a statistical analysis of the properties of the obscuring material around active galactic nuclei (AGN). This study represents the first of its kind for an ultra-hard X-ray (14-195keV; Swift/BAT) volume-limited (DL<40 Mpc) sample of 24 Seyfert (Sy) galaxies (BCS40 sample) using high angular resolution infrared data and various torus models: smooth, clumpy and two-phase torus models and clumpy disc+wind models. We find that the smooth, clumpy and two-phase torus models (i.e. without including the polar dusty wind component) and disc+wind models provide best fits for a comparable number of galaxies, 8/24 (33.3%) and 9/24 (37.5%), respectively. We find that the best-fitted models depend on the hydrogen column density (NH), which is related to the X-ray (unobscured/obscured) and/or optical (Sy1/Sy2) classification. In particular, smooth, clumpy and two-phase torus models best reproduce the infrared (IR) emission of AGN with relatively high hydrogen column density (median value of log (NH)=23.5+-0.8; i.e. Sy2). However, clumpy disc+wind models provide the best fits to the nuclear IR spectral energy distributions (SEDs) of Sy1/1.8/1.9 (median value of log (NH)=21.0+-1.0), specifically in the near-IR (NIR) range. The success of the disc+wind models in fitting the NIR emission of Sy1 galaxies is due to the combination of adding large graphite grains to the dust composition and self-obscuration effects caused by the wind at intermediate inclinations. In general, we find that the Seyfert galaxies having unfavourable (favourable) conditions, i.e. nuclear hydrogen column density and Eddington ratio, for launching IR dusty polar outflows are best-fitted with smooth, clumpy and two-phase torus (disk+wind) models confirming the predictions from simulations. Therefore, our results indicate that the nature of the inner dusty structure in AGN depend on the intrinsic AGN properties. △ Less

Submitted 7 October, 2022; originally announced October 2022.

Comments: Accepted for publication in A&A, 43 pages, 43 figures, 16 tables

Journal ref: A&A 667, A140 (2022)

On the origin of radio-loudness in active galactic nuclei using far-infrared polarimetric observations

Authors: Enrique Lopez-Rodriguez, Makoto Kishimoto, Robert Antonucci, Mitchell C. Begelman, Noemie Globus, Roger Blandford

Abstract: The dichotomy between radio-loud (RL) and radio-quiet (RQ) active galactic nuclei (AGN) is thought to be intrinsically related to radio jet production. This difference may be explained by the presence of a strong magnetic field (B-field) that enhances, or is the cause of, the accretion activity and the jet power. Here, we report the first evidence of an intrinsic difference in the dust polarized e… ▽ More The dichotomy between radio-loud (RL) and radio-quiet (RQ) active galactic nuclei (AGN) is thought to be intrinsically related to radio jet production. This difference may be explained by the presence of a strong magnetic field (B-field) that enhances, or is the cause of, the accretion activity and the jet power. Here, we report the first evidence of an intrinsic difference in the dust polarized emission cores of four RL and five RQ obscured AGN using 89 $μ$m polarization with HAWC+/SOFIA. We find that the thermal polarized emission increases with the nuclear radio-loudness, $R_{20} = L_{\rm 5GHz}/ L_{\rm 20μm}$. The dust emission cores of RL AGN are measured to be polarized, $\sim5-11$%, while RQ AGN are unpolarized, $<1$%. For RQ AGN, our results are consistent with the observed region being filled with an unmagnetized or highly turbulent, disk and/or expanding outflow at scales of $5-130$ pc from the AGN. For RL AGN, the measured $89$ $μ$m polarization arises primarily from magnetically aligned dust grains associated with a $5-130$ pc-scale dusty obscuring structure with a toroidal B-field orientation highly offset, $65\pm22^{\circ}$, with respect to the jet axis. Our results indicate that the size and strength of the B-fields surrounding the AGN are intrinsically related to the strength of the jet power -- the stronger the jet power is, the larger and stronger the toroidal B-field is. The detection of a $\le130$ pc-scale ordered toroidal B-field suggests that a) the infalling gas that fuels RL AGN is magnetized, b) there is a magnetohydrodynamic wind that collimates the jet, and/or c) the jet is able to magnetize its surroundings. △ Less

Submitted 7 April, 2023; v1 submitted 19 July, 2022; originally announced July 2022.

Comments: 23 pages, 6 figures, 7 tables. Accepted for publication; ApJ

The Role of Magnetic Fields in the Stability and Fragmentation of Filamentary Molecular Clouds: Two Case Studies at OMC-3 and OMC-4

Authors: Pak Shing Li, Enrique Lopez-Rodriguez, Archana Soam, Richard I. Klein

Abstract: We present the stability analysis of two regions, OMC-3 and OMC-4, in the massive and long molecular cloud complex of Orion A. We obtained $214~μ$m HAWC+/SOFIA polarization data, and we make use of archival data for the column density and C$^{18}$O (1-0) emission line. We find clear depolarization in both observed regions and that the polarization fraction is anti-correlated with the column densit… ▽ More We present the stability analysis of two regions, OMC-3 and OMC-4, in the massive and long molecular cloud complex of Orion A. We obtained $214~μ$m HAWC+/SOFIA polarization data, and we make use of archival data for the column density and C$^{18}$O (1-0) emission line. We find clear depolarization in both observed regions and that the polarization fraction is anti-correlated with the column density and the polarization-angle dispersion function. We find that the filamentary cloud and dense clumps in OMC-3 are magnetically supercritical and strongly subvirial. This region should be in the gravitational collapse phase and is consistent with many young stellar objects (YSOs) forming in the region. Our histogram of relative orientations (HROs) analysis shows that the magnetic field is dynamically sub-dominant in the dense gas structures of OMC-3. We present the first polarization map of OMC-4. We find that the observed region is generally magnetically subcritical except for an elongated dense core, which could be a result of projection effect of a filamentary structure aligned close to the line-of-sight. The relative large velocity dispersion and the unusual positive shape parameters at high column densities in the HROs analysis suggest that our viewing angle may be close to axes of filamentary substructures in OMC-4. The dominating strong magnetic field in OMC-4 is unfavorable for star formation and is consistent with much fewer YSOs than in OMC-3. △ Less

Submitted 31 May, 2022; originally announced June 2022.

Comments: 19 pages, 11 figures, accepted by MNRAS

SOFIA observations of 30 Doradus: II -- Magnetic fields and large scale gas kinematics

Authors: Le Ngoc Tram, Lars Bonne, Yue Hu, Enrique Lopez-Rodriguez, Jordan A. Guerra, Pierre Lesaffre, Antoine Gusdorf, Thiem Hoang, Min-Young Lee, Alex Lazarian, B-G Andersson, Simon Coude, Archana Soam, William D. Vacca, Hyeseung Lee, Michael Gordon

Abstract: The heart of the Large Magellanic Cloud, 30 Doradus, is a complex region with a clear core-halo structure. Feedback from the stellar cluster R$\,$136 has been shown to be the main source of energy creating multiple pc-scale expanding-shells in the outer region, and carving a nebula core in the proximity of the ionization source. We present the morphology and strength of the magnetic fields (B-fiel… ▽ More The heart of the Large Magellanic Cloud, 30 Doradus, is a complex region with a clear core-halo structure. Feedback from the stellar cluster R$\,$136 has been shown to be the main source of energy creating multiple pc-scale expanding-shells in the outer region, and carving a nebula core in the proximity of the ionization source. We present the morphology and strength of the magnetic fields (B-fields) of 30 Doradus inferred from the far-infrared polarimetric observations by SOFIA/HAWC+ at 89, 154, and 214$\,μ$m. The B-field morphology is complex, showing bending structures around R$\,$136. In addition, we use high spectral and angular resolution [\textsc{CII}] observations from SOFIA/GREAT and CO(2-1) from APEX. The kinematic structure of the region correlates with the B-field morphology and shows evidence of multiple expanding shells. Our B-field strength maps, estimated using the Davis-Chandrasekhar-Fermi method and structure-function, show variations across the cloud within a maximum of 600, 450, and 350$\,μ$G at 89, 154, and 214$\,μ$m, respectively. We estimated that the majority of the 30 Doradus clouds are sub-critical and sub-Alfvénic. The probability distribution function of the gas density shows that the turbulence is mainly compressively driven, while the plasma beta parameter indicates supersonic turbulence. We show that the B-field is sufficient to hold the cloud structure integrity under feedback from R$\,$136. We suggest that supersonic compressive turbulence enables the local gravitational collapse and triggers a new generation of stars to form. The velocity gradient technique (VGT) using [\textsc{CII}] and CO(2-1) is likely to confirm these results. △ Less

Submitted 12 December, 2022; v1 submitted 24 May, 2022; originally announced May 2022.

Comments: 26 pages, 17 figures, 3 tables, accepted to ApJ

Extragalactic magnetism with SOFIA (SALSA Legacy Program) -- IV: Program overview and first results on the polarization fraction

Authors: Enrique Lopez-Rodriguez, Sui Ann Mao, Rainer Beck, Alejandro S. Borlaff, Evangelia Ntormousi, Konstantinos Tassis, Daniel A. Dale, Julia Roman-Duval, Kandaswamy Subramanian, Sergio Martin-Alvarez, Pamela M. Marcum, Susan E. Clark, William T. Reach, Doyal A. Harper, Ellen G. Zweibel

Abstract: We present the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program) with a set of 14 nearby ($<20$ Mpc) galaxies with resolved imaging polarimetric observations using HAWC+ from $53$ to $214$ $μ$m at a resolution of $5-18$" ($90$ pc $-$ $1$ kpc). We introduce the definitions and background on extragalactic magnetism, and present the scientific motivation an… ▽ More We present the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program) with a set of 14 nearby ($<20$ Mpc) galaxies with resolved imaging polarimetric observations using HAWC+ from $53$ to $214$ $μ$m at a resolution of $5-18$" ($90$ pc $-$ $1$ kpc). We introduce the definitions and background on extragalactic magnetism, and present the scientific motivation and sample selection of the program. Here, we focus on the general trends in the emissive polarization fraction. Far-infrared polarimetric observations trace the thermal polarized emission of magnetically aligned dust grains across the galaxy disks with polarization fractions of $P=0-15$% in the cold, $T_{\rm d} = [19,48]$ K, and dense, $\log_{10}(N_{\rm HI+H_{2}}) = [19.96,22.91]$, interstellar medium. The spiral galaxies show a median $\langle P_{154μm} \rangle = 3.3\pm0.9 $% across the disks. We report the first polarized spectrum of starburst galaxies showing a minimum within $89-154$ $μ$m. The falling $53-154$ $μ$m polarized spectrum may be due to a decrease in the dust grain alignment efficiency produced by variations in dust temperatures along the line-of-sight in the galactic outflow. We find that the starburst galaxies and the star-forming regions within normal galaxies have the lowest polarization fractions. We find that 50% (7 out of 14) of the galaxies require a broken power-law in the $P-N_{HI+H_{2}}$ and $P-T_{d}$ relations with three different trends. Group 1 has a relative increase of anisotropic random B-fields produced by compression or shear of B-fields in the galactic outflows, starburst rings, and inner-bar of galaxies; and Groups 2 and 3 have a relative increase of isotropic random B-fields driven by star-forming regions in the spiral arms, and/or an increase of dust grain alignment efficiency caused by shock-driven regions or evolutionary stages of a galaxy. △ Less

Submitted 7 July, 2022; v1 submitted 2 May, 2022; originally announced May 2022.

Comments: 43 pages, 28 figures. Accepted for publication to ApJ

Extragalactic magnetism with SOFIA (SALSA Legacy Program) -- III: First data release and on-the-fly polarization mapping characterization

Authors: Enrique Lopez-Rodriguez, Melanie Clarke, Sachin Shenoy, William Vacca, Simon Coude, Ryan Arneson, Peter Ashton, Sarah Eftekharzadeh, Rainer Beck, John E. Beckman, Alejandro S. Borlaff, Susan E. Clark, Daniel A. Dale, Sergio Martin-Alvarez, Evangelia Ntormousi, William T. Reach, Julia Roman-Duval, Konstantinos Tassis, Doyal A. Harper, Pamela M. Marcum

Abstract: We describe the data processing of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program). This first data release presents 33% (51.34h out of 155.7h, including overheads) of the total awarded time taken from January 2020 to December 2021. Our observations were performed using the newly implemented on-the-fly mapping (OTFMAP) technique in the polarimetric mode. We present the pipe… ▽ More We describe the data processing of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program). This first data release presents 33% (51.34h out of 155.7h, including overheads) of the total awarded time taken from January 2020 to December 2021. Our observations were performed using the newly implemented on-the-fly mapping (OTFMAP) technique in the polarimetric mode. We present the pipeline steps to obtain homogeneously reduced high-level data products of polarimetric maps of galaxies for use in scientific analysis. Our approach has a general design and can be applied to sources smaller than the field-of-view of the HAWC+ array in any given band. We estimate that the OTFMAP polarimetric mode offers a reduction of observing overheads by a factor 2.34, and an improvement in sensitivity by a factor 1.80 when compared to previously obtained polarimetric observations using the chopping and nodding mode. The OTFMAP is a significant optimization of the polarimetric mode of HAWC+ as it ultimately reduces the cost of operations of SOFIA/HAWC+ by increasing the science collected per hour of observation up to an overall factor of 2.49. The OTFMAP polarimetric mode is the standard observing strategy of SALSA. The results and quantitative analysis of this first data release are presented in Papers IV and V of the series. △ Less

Submitted 28 April, 2022; originally announced April 2022.

Comments: 26 pages, 16 figures, 5 tables. Submitted to ApJ (Comments welcome)

The complex infrared dust continuum emission of NGC1068: ground-based N- and Q-band spectroscopy and new radiative transfer models

Authors: César Ivan Victoria-Ceballos, Omaira González-Martín, Jacopo Fritz, Cristina Ramos Almeida, Enrique López-Rodríguez, Santiago García-Burillo, Almudena Alonso-Herrero, Mariela Martínez-Paredes, Donaji Esparza-Arredondo, Natalia Osorio-Clavijo

Abstract: Thanks to ground-based infrared and sub-mm observations the study of the dusty torus of nearby AGN has greatly advanced in the last years. With the aim of further investigating the nuclear mid-infrared emission of the archetypal Seyfert 2 galaxy NGC1068, here we present a fitting to the N- and Q-band Michelle/Gemini spectra. We initially test several available SED libraries, including a smooth, cl… ▽ More Thanks to ground-based infrared and sub-mm observations the study of the dusty torus of nearby AGN has greatly advanced in the last years. With the aim of further investigating the nuclear mid-infrared emission of the archetypal Seyfert 2 galaxy NGC1068, here we present a fitting to the N- and Q-band Michelle/Gemini spectra. We initially test several available SED libraries, including a smooth, clumpy and two phase torus models, and a clumpy disk plus wind model. Although, the smooth torus model describe the spectra of NGC1068 if we allow to vary some model parameters among the two spectral bands. Motivated by this result, we produced new SEDs using the radiative transfer code SKIRT. We use two concentric tori that allow us to test a more complex geometry. We test different values for the inner and outer radii, half opening angle, radial and polar exponent of the power-law density profile, opacity, and viewing angle. Furthermore, we also test the dust grains size and different optical and calorimetric properties of silicate grains. The best fitting model consists of two concentric components with outer radii of 1.8 and 28 pc, respectively. We find that the size and the optical and calorimetric properties of graphite and silicate grains in the dust structure are key to reproduce the spectra of NGC1068. We conclude that the dust in NGC1068 reaches different scales, where the highest contribution to the mid-infrared is given by a central and compact component. A less dense and extended component is present, which can be either part of the same torus (conforming a flared disk) or can represent the emission of a polar dust component, as already suggested from interferometric observations. △ Less

Submitted 27 January, 2022; originally announced January 2022.

The best broths are cooked in the oldest pans: revisiting the archival HST/FOC observations of quasars

Authors: F. Marin, T. Barnouin, E. Lopez-Rodriguez

Abstract: The Faint Object Camera (FOC) aboard the Hubble Space Telescope (HST) observed 26 individual active galactic nuclei (AGNs) in ultraviolet imaging polarimetry between 1990 and 2002. Tremendous progresses have been made thanks to those high spatial resolution, high signal-to-noise ratio observations, such as the identification of the location of hidden active nuclei and the three dimensional arrange… ▽ More The Faint Object Camera (FOC) aboard the Hubble Space Telescope (HST) observed 26 individual active galactic nuclei (AGNs) in ultraviolet imaging polarimetry between 1990 and 2002. Tremendous progresses have been made thanks to those high spatial resolution, high signal-to-noise ratio observations, such as the identification of the location of hidden active nuclei and the three dimensional arrangement of polar material within the first hundred of parsecs around the central core. However, not all AGN observations have been reduced and analyzed, and none in a standardized framework. In this lecture note, we present our project of downloading, reducing and analyzing all the AGN HST/FOC observations that were achieved using a consistent, novel and open-access reduction pipeline. We briefly present our methodology and show the first, preliminary result from our reduction pipeline: NGC 1068. △ Less

Submitted 17 January, 2022; originally announced January 2022.

Comments: 4 pages, 2 figures, to be published in the proceedings of the French SF2A

Mapping the magnetic field in the Taurus/B211 filamentary cloud with SOFIA HAWC+ and comparing with simulation

Authors: Pak Shing Li, Enrique Lopez-Rodriguez, Hamza Ajeddig, Philippe André, Christopher F. McKee, Jeonghee Rho, Richard I. Klein

Abstract: Optical and infrared polarization mapping and recent Planck observations of the filamentary cloud L1495 in Taurus show that the large-scale magnetic field is approximately perpendicular to the long axis of the cloud. We use the HAWC+ polarimeter on SOFIA to probe the complex magnetic field in the B211 part of the cloud. Our results reveal a dispersion of polarization angles of $36^\circ$, about fi… ▽ More Optical and infrared polarization mapping and recent Planck observations of the filamentary cloud L1495 in Taurus show that the large-scale magnetic field is approximately perpendicular to the long axis of the cloud. We use the HAWC+ polarimeter on SOFIA to probe the complex magnetic field in the B211 part of the cloud. Our results reveal a dispersion of polarization angles of $36^\circ$, about five times that measured on a larger scale by Planck. Applying the Davis-Chandrasekhar-Fermi (DCF) method with velocity information obtained from IRAM 30m C$^{18}$O(1-0) observations, we find two distinct sub-regions with magnetic field strengths differing by more than a factor 3. The quieter sub-region is magnetically critical and sub-Alfvénic; the field is comparable to the average field measured in molecular clumps based on Zeeman observations. The more chaotic, super-Alfvénic sub-region shows at least three velocity components, indicating interaction among multiple substructures. Its field is much less than the average Zeeman field in molecular clumps, suggesting that the DCF value of the field there may be an underestimate. Numerical simulation of filamentary cloud formation shows that filamentary substructures can strongly perturb the magnetic field. DCF and true field values in the simulation are compared. Pre-stellar cores are observed in B211 and are seen in our simulation. The appendices give a derivation of the standard DCF method that allows for a dispersion in polarization angles that is not small, present an alternate derivation of the structure function version of the DCF method, and treat fragmentation of filaments. △ Less

Submitted 24 November, 2021; originally announced November 2021.

Comments: 27 pages, 12 figures

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

Extragalactic magnetism with SOFIA (Legacy Program) -- II: A magnetically-driven flow in the starburst ring of NGC 1097

Authors: Enrique Lopez-Rodriguez, Rainer Beck, Susan E. Clark, Annie Hughes, Alejandro S. Borlaff, Evangelia Ntormousi, Lucas Grosset, Konstantinos Tassis, John E. Beckman, Kandaswamy Subramanian, Daniel Dale, Tanio Diaz-Santos, Legacy Team

Abstract: Galactic bars are frequent in disk galaxies and they may support the transfer of matter towards the central engine of active nuclei. The barred galaxy NGC 1097 has magnetic forces controlling the gas flow at several kpc scales, which suggest that magnetic fields (B-fields) are dynamically important along the bar and nuclear ring. However, the effect of the B-field on the gas flows in the central k… ▽ More Galactic bars are frequent in disk galaxies and they may support the transfer of matter towards the central engine of active nuclei. The barred galaxy NGC 1097 has magnetic forces controlling the gas flow at several kpc scales, which suggest that magnetic fields (B-fields) are dynamically important along the bar and nuclear ring. However, the effect of the B-field on the gas flows in the central kpc scale has not been characterized. Using thermal polarized emission at $89$ $μ$m with HAWC+/SOFIA, here, we measure that the polarized flux is spatially located at the contact regions of the outer-bar with the starburst ring. The linear polarization decomposition analysis shows that the $89$ $μ$m and radio ($3.5$ and $6.2$ cm) polarization traces two different modes, $m$, of the B-field: a constant B-field orientation and dominated by $m=0$ at $89$ $μ$m, and a spiral B-field dominated by $m=2$ at radio. We show that the B-field at 89 $μ$m is concentrated in the warmest region of a shock driven by the galactic-bar dynamics in the contact regions between the outer-bar with the starburst ring. Radio polarization traces a superposition of the spiral B-field outside and within the starburst ring. According to Faraday rotation measures between $3.5$ and $6.2$ cm, the radial component of the B-field along the contact regions points toward the galaxy's center on both sides. We conclude that gas streams outside and within the starburst ring follow the B-field, which feeds the black hole with matter from the host galaxy. △ Less

Submitted 7 October, 2021; v1 submitted 19 July, 2021; originally announced July 2021.

Comments: 17 pages, 10 figures, Accepted for publication to 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

SOFIA observations of 30 Doradus: I -- Far-Infrared dust polarization and implications for grain alignment and disruption by radiative torques

Authors: Le Ngoc Tram, Thiem Hoang, Enrique Lopez-Rodriguez, Simon Coudé, Archana Soam, B-G Andersson, Min-Young Lee, Lars Bonne, William D. Vacca, Hyeseung Lee

Abstract: Located in the Large Magellanic cloud and mostly irradiated by a massive-star cluster R$\,$136, 30 Doradus is an ideal target to test the leading theory of the grain alignment and rotational disruption by RAdiative Torques (RATs). Here, we use publicly available polarized thermal dust emission observations of 30 Doradus at 89, 154, and 214$\,μ$m using SOFIA/HAWC+. We analyse the variation of the d… ▽ More Located in the Large Magellanic cloud and mostly irradiated by a massive-star cluster R$\,$136, 30 Doradus is an ideal target to test the leading theory of the grain alignment and rotational disruption by RAdiative Torques (RATs). Here, we use publicly available polarized thermal dust emission observations of 30 Doradus at 89, 154, and 214$\,μ$m using SOFIA/HAWC+. We analyse the variation of the dust polarization degree ($p$) with the total emission intensity ($I$), the dust temperature ($T_{\rm d}$), and the gas column density ($N_{\rm H}$) constructed from ${\it Herschel}$ data. The 30 Doradus complex is divided into two main regions relative to R$\,$136, namely North and South. In the North, we find that the polarization degree first decreases and then increases before decreasing again when the dust temperature increases toward the irradiating cluster R$\,$136. The first depolarization likely arises from the decrease of grain alignment efficiency toward the dense medium due to the attenuation of the interstellar radiation field and the increase of the gas density. The second trend (the increase of $p$ with $T_{\rm d}$) is consistent with the RAT alignment theory. The final trend (the decrease of $p$ with $T_{\rm d}$) is consistent with the RAT alignment theory only when the grain rotational disruption by RATs is taken into account. In the South, we find that the polarization degree is nearly independent of the dust temperature, while the grain alignment efficiency is higher around the peak of the gas column density and decreases toward the radiation source. The latter feature is also consistent with the prediction of the rotational disruption by RATs. △ Less

Submitted 31 July, 2021; v1 submitted 20 May, 2021; originally announced May 2021.

Comments: 19 pages, 13 figures, 1 table, accepted to ApJ

Extragalactic Magnetism with SOFIA (Legacy Program) -- I: The magnetic field in the multi-phase interstellar medium of M51

Authors: Alejandro S. Borlaff, Enrique Lopez-Rodriguez, Rainer Beck, Rodion Stepanov, Eva Ntormousi, Annie Hughes, Konstantinos Tassis, Pamela M. Marcum, Lucas Grosset, John E. Beckman, Leslie Proudfit, Susan E. Clark, Tanio Díaz-Santos, Sui Ann Mao, William T. Reach, Julia Roman-Duval, Kandaswamy Subramanian, Le Ngoc Tram, Ellen G. Zweibel, SOFIA Legacy Team

Abstract: The recent availability of high-resolution far-infrared (FIR) polarization observations of galaxies using HAWC+/SOFIA has facilitated studies of extragalactic magnetic fields in the cold and dense molecular disks.We investigate if any significant structural differences are detectable in the kpc-scale magnetic field of the grand design face-on spiral galaxy M51 when traced within the diffuse (radio… ▽ More The recent availability of high-resolution far-infrared (FIR) polarization observations of galaxies using HAWC+/SOFIA has facilitated studies of extragalactic magnetic fields in the cold and dense molecular disks.We investigate if any significant structural differences are detectable in the kpc-scale magnetic field of the grand design face-on spiral galaxy M51 when traced within the diffuse (radio) and the dense and cold (FIR) interstellar medium (ISM). Our analysis reveals a complex scenario where radio and FIR polarization observations do not necessarily trace the same magnetic field structure. We find that the magnetic field in the arms is wrapped tighter at 154um than at 3 and 6 cm; statistically significant lower values for the magnetic pitch angle are measured at FIR in the outskirts (R > 7 kpc) of the galaxy. This difference is not detected in the interarm region. We find strong correlations of the polarization fraction and total intensity at FIR and radio with the gas column density and 12CO(1-0) velocity dispersion. We conclude that the arms show a relative increase of small-scale turbulent B-fields at regions with increasing column density and dispersion velocities of the molecular gas. No correlations are found with HI neutral gas. The star formation rate shows a clear correlation with the radio polarized intensity, which is not found in FIR, pointing to a small-scale dynamo-driven B-field amplification scenario. This work shows that multi-wavelength polarization observations are key to disentangling the interlocked relation between star formation, magnetic fields, and gas kinematics in the multi-phase ISM. △ Less

Submitted 8 September, 2021; v1 submitted 19 May, 2021; originally announced May 2021.

Comments: 35 pages, 25 figures, 5 tables. Accepted for publication in ApJ

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)

The Magnetic Field Across the Molecular Warped Disk of Centaurus A

Authors: Enrique Lopez-Rodriguez

Abstract: Magnetic fields are amplified as a consequence of galaxy formation and turbulence-driven dynamos. Galaxy mergers can potentially amplify the magnetic fields from their progenitors, making the magnetic fields dynamically important. However, the effect of mergers on magnetic fields is still poorly understood. We use thermal polarized emission observations to trace the magnetic fields in the molecula… ▽ More Magnetic fields are amplified as a consequence of galaxy formation and turbulence-driven dynamos. Galaxy mergers can potentially amplify the magnetic fields from their progenitors, making the magnetic fields dynamically important. However, the effect of mergers on magnetic fields is still poorly understood. We use thermal polarized emission observations to trace the magnetic fields in the molecular disk of the nearest radio active galaxy, Centaurus A, which is thought to be the remnant of a merger. Here, we detect that the magnetic field orientations in the plane of the sky are tightly following the $\sim3.0$ kpc-scale molecular warped disk. Our simple regular large-scale axisymmetric spiral magnetic field model can explain, to some extent, the averaged magnetic field orientations across the disk projected on the sky. Our observations also suggest the presence of small-scale turbulent fields, whose relative strength increases with velocity dispersion and column density. These results have strong implications for understanding the generation and role of magnetic fields in the formation of galaxies across cosmic time. △ Less

Submitted 1 April, 2021; originally announced April 2021.

Comments: 54 pages, 11 figures, 5 tables. Accepted for publication in Nature Astronomy

The strength and structure of the magnetic field in the galactic outflow of M82

Authors: Enrique Lopez-Rodriguez, Jordan Guerra, Mahboubeh Asgari-Targhi, Joan T. Schmelz

Abstract: Galactic outflows driven by starbursts can modify the galactic magnetic fields and drive them away from the galactic planes. Here, we quantify how these fields may magnetize the intergalactic medium. We estimate the strength and structure of the fields in the starburst galaxy M82 using thermal polarized emission observations from SOFIA/HAWC+ and a potential field extrapolation commonly used in sol… ▽ More Galactic outflows driven by starbursts can modify the galactic magnetic fields and drive them away from the galactic planes. Here, we quantify how these fields may magnetize the intergalactic medium. We estimate the strength and structure of the fields in the starburst galaxy M82 using thermal polarized emission observations from SOFIA/HAWC+ and a potential field extrapolation commonly used in solar physics. We modified the Davis-Chandrasekhar-Fermi method to account for the large-scale flow and the turbulent field. Results show that the observed magnetic fields arise from the combination of a large-scale ordered potential field associated with the outflow and a small-scale turbulent field associated with bow-shock-like features. Within the central $900$ pc radius, the large-scale field accounts for $53\pm4$% of the observed turbulent magnetic energy with a median field strength of $305\pm15$ $μ$G, while small-scale turbulent magnetic fields account for the remaining $40\pm5$% with a median field strength of $222\pm19$ $μ$G. We estimate that the turbulent kinetic and turbulent magnetic energies are in close equipartition up to $\sim2$ kpc (measured), while the turbulent kinetic energy dominates at $\sim7$ kpc (extrapolated). We conclude that the fields are frozen into the ionized outflowing medium and driven away kinetically. The magnetic field lines in the galactic wind of M82 are `open,' providing a direct channel between the starburst core and the intergalactic medium. Our novel approach offers the tools needed to quantify the effects of outflows on galactic magnetic fields as well as their influence on the intergalactic medium and evolution of energetic particles. △ Less

Submitted 15 April, 2021; v1 submitted 5 February, 2021; originally announced February 2021.

Comments: 17 pages, 9 figures, Accepted for publication in ApJ

Magnetized filamentary gas flows feeding the young embedded cluster in Serpens South

Authors: Thushara G. S. Pillai, Dan P. Clemens, Stefan Reissl, Philip C. Myers, Jens Kauffmann, Enrique Lopez-Rodriguez, Felipe O. Alves, Gabriel P. Franco, Jonathan D. Henshaw, Karl M. Menten, Fumitaka Nakamura, Daniel Seifried, Koji Sugitani, Helmut Wiesemeyer

Abstract: Observations indicate that molecular clouds are strongly magnetized, and that magnetic fields influence the formation of stars. A key observation supporting the conclusion that molecular clouds are significantly magnetized is that the orientation of their internal structure is closely related to that of the magnetic field. At low column densities the structure aligns parallel with the field, where… ▽ More Observations indicate that molecular clouds are strongly magnetized, and that magnetic fields influence the formation of stars. A key observation supporting the conclusion that molecular clouds are significantly magnetized is that the orientation of their internal structure is closely related to that of the magnetic field. At low column densities the structure aligns parallel with the field, whereas at higher column densities, the gas structure is typically oriented perpendicular to magnetic fields, with a transition at visual extinctions $A_V\gtrsim{}3~\rm{}mag$. Here we use far-infrared polarimetric observations from the HAWC+ polarimeter on SOFIA to report the discovery of a further transition in relative orientation, i.e., a return to parallel alignment at $A_V\gtrsim{}21~\rm{}mag$ in parts of the Serpens South cloud. This transition appears to be caused by gas flow and indicates that magnetic supercriticality sets in near $A_V\gtrsim{}21~\rm{}mag$, allowing gravitational collapse and star cluster formation to occur even in the presence of relatively strong magnetic fields. △ Less

Submitted 29 September, 2020; originally announced September 2020.

Comments: 16 pages, 4 figures, Published in Nature Astronomy (August 2020). This is the authors' version before final edits. Link to the NA publication: https://www.nature.com/articles/s41550-020-1172-6

Journal ref: Nature Astronomy, August 2020

The polarized spectral energy distribution of NGC 4151

Authors: F. Marin, J. Le Cam, E. Lopez-Rodriguez, M. Kolehmainen, B. L. Babler, M. R. Meade

Abstract: NGC 4151 is among the most well-studied Seyfert galaxies that does not suffer from strong obscuration along the observer's line-of-sight. This allows to probe the central active galactic nucleus (AGN) engine with photometry, spectroscopy, reverberation mapping or interferometry. Yet, the broadband polarization from NGC 4151 has been poorly examined in the past despite the fact that polarimetry giv… ▽ More NGC 4151 is among the most well-studied Seyfert galaxies that does not suffer from strong obscuration along the observer's line-of-sight. This allows to probe the central active galactic nucleus (AGN) engine with photometry, spectroscopy, reverberation mapping or interferometry. Yet, the broadband polarization from NGC 4151 has been poorly examined in the past despite the fact that polarimetry gives us a much cleaner view of the AGN physics than photometry or spectroscopy alone. In this paper, we compile the 0.15 -- 89.0 $μ$m total and polarized fluxes of NGC 4151 from archival and new data in order to examine the physical processes at work in the heart of this AGN. We demonstrate that, from the optical to the near-infrared (IR) band, the polarized spectrum of NGC 4151 shows a much bluer power-law spectral index than that of the total flux, corroborating the presence of an optically thick, locally heated accretion flow, at least in its near-IR emitting radii. Specific signatures from the atmosphere of the accretion structure are tentatively found at the shortest ultraviolet (UV) wavelengths, before the onset of absorption opacity. Otherwise, dust scattering appears to be the dominant contributor from the near-UV to near-IR polarized spectrum, superimposed onto a weaker electron component. We also identify a change in the polarization processes from the near-IR to the mid-IR, most likely associated with the transition from Mie scattering to dichroic absorption from aligned dust grains in the dusty torus or narrow-line region. Finally, we present and dicuss the very first far-infrared polarization measurement of NGC 4151 at 89 $μ$m. △ Less

Submitted 28 May, 2020; originally announced May 2020.

Comments: 8 pages, 4 figures, accepted for publication in MNRAS

MSC Class: 85-06 ACM Class: J.2.3; J.2.9

SOFIA/HAWC+ traces the magnetic fields in NGC 1068

Authors: E. Lopez-Rodriguez, C. D. Dowell, T. J. Jones, D. A. Harper, M. Berthoud, D. Chuss, D. A. Dale, J. A. Guerra, R. T. Hamilton, L. W. Looney, J. M. Michail, R. Nikutta, G. Novak, F. P. Santos, K. Sheth, J. Siah, J. Staguhn, I. W. Stephens, K. Tassis, C. Q. Trinh, D. Ward-Thompson, M. Werner, E. J. Wollack, E. Zweibel

Abstract: We report the first detection of galactic spiral structure by means of thermal emission from magnetically aligned dust grains. Our 89 $μ$m polarimetric imaging of NGC 1068 with the High-resolution Airborne Wideband Camera/Polarimeter (HAWC+) on NASA's Stratospheric Observatory for Infrared Astronomy (SOFIA) also sheds light on magnetic field structure in the vicinity of the galaxy's inner-bar and… ▽ More We report the first detection of galactic spiral structure by means of thermal emission from magnetically aligned dust grains. Our 89 $μ$m polarimetric imaging of NGC 1068 with the High-resolution Airborne Wideband Camera/Polarimeter (HAWC+) on NASA's Stratospheric Observatory for Infrared Astronomy (SOFIA) also sheds light on magnetic field structure in the vicinity of the galaxy's inner-bar and active galactic nucleus (AGN). We find correlations between the 89 $μ$m magnetic field vectors and other tracers of spiral arms, and a symmetric polarization pattern as a function of the azimuthal angle arising from the projection and inclination of the disk field component in the plane of the sky. The observations can be fit with a logarithmic spiral model with pitch angle of $16.9^{+2.7}_{-2.8}$$^{\circ}$ and a disk inclination of $48\pm2^{\circ}$. We infer that the bulk of the interstellar medium from which the polarized dust emission originates is threaded by a magnetic field that closely follows the spiral arms. Inside the central starburst disk ($<1.6$ kpc), the degree of polarization is found to be lower than for far-infrared sources in the Milky Way, and has minima at the locations of most intense star formation near the outer ends of the inner-bar. Inside the starburst ring, the field direction deviates from the model, becoming more radial along the leading edges of the inner-bar. The polarized flux and dust temperature peak $\sim 3-6$" NE of the AGN at the location of a bow shock between the AGN outflow and the surrounding interstellar medium, but the AGN itself is weakly polarized ($< 1$%) at both 53 and 89 \um. △ Less

Submitted 15 November, 2019; v1 submitted 15 July, 2019; originally announced July 2019.

Comments: 13 pages, 7 figures, Accepted for publication by ApJ

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)

A Compact Jet at the Infrared Heart of the Prototypical Low-Luminosity AGN in NGC 1052

Authors: J. A. Fernández-Ontiveros, N. López-Gonzaga, M. A. Prieto, J. A. Acosta-Pulido, E. Lopez-Rodriguez, D. Asmus, K. R. W. Tristram

Abstract: The feeble radiative efficiency characteristic of Low-Luminosity Active Galactic Nuclei (LLAGN) is ascribed to a sub-Eddington accretion rate, typically at $\log(L_{\rm bol}/L_{\rm edd}) \lesssim -3$. At the finest angular resolutions that are attainable nowadays using mid-infrared (mid-IR) interferometry, the prototypical LLAGN in NGC 1052 remains unresolved down to $< 5\, \rm{mas}$ (… ▽ More The feeble radiative efficiency characteristic of Low-Luminosity Active Galactic Nuclei (LLAGN) is ascribed to a sub-Eddington accretion rate, typically at $\log(L_{\rm bol}/L_{\rm edd}) \lesssim -3$. At the finest angular resolutions that are attainable nowadays using mid-infrared (mid-IR) interferometry, the prototypical LLAGN in NGC 1052 remains unresolved down to $< 5\, \rm{mas}$ ($0.5\, \rm{pc}$). This is in line with non-thermal emission from a compact jet, a scenario further supported by a number of evidences: the broken power-law shape of the continuum distribution in the radio-to-UV range; the $\sim 4\%$ degree of polarisation measured in the nuclear mid-IR continuum, together with the mild optical extinction ($A_V \sim 1\, \rm{mag}$); and the "harder when brighter" behaviour of the X-ray spectrum, indicative of self-Compton synchrotron radiation. A remarkable feature is the steepness of the IR-to-UV core continuum, characterised by a power-law index of $\sim 2.6$, as compared to the canonical value of $0.7$. Alternatively, to explain the interferometric data by thermal emission would require an exceptionally compact dust distribution when compared to those observed in nearby AGN, with $A_V \gtrsim 2.8\, \rm{mag}$ to account for the IR polarisation. This is in contrast with several observational evidences against a high extinction along the line of sight, including the detection of the nucleus in the UV range and the well defined shape of the power-law continuum. The case of NGC 1052 shows that compact jets can dominate the nuclear emission in LLAGN across the whole electromagnetic spectrum, a scenario that might be common among this class of active nuclei. △ Less

Submitted 25 March, 2019; v1 submitted 12 March, 2019; originally announced March 2019.

Comments: Accepted for publication by MNRAS

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

SOFIA Community Science I: HAWC+ Polarimetry of 30 Doradus

Authors: M. S. Gordon, E. Lopez-Rodriguez, B. -G. Andersson, M. Clarke, S. Coude, A. Moullet, S. N. Richards, R. Y. Shuping, W. Vacca, H. Yorke

Abstract: The Stratospheric Observatory for Infrared Astronomy (SOFIA) is a Boeing 747SP aircraft modified to accommodate a 2.7 meter gyro-stabilized telescope, which is mainly focused to studying the Universe at infrared wavelengths. As part of the Strategic Director's Discretionary Time (S-DDT) program, SOFIA performs observations of relevant science cases and immediately offers science-ready data product… ▽ More The Stratospheric Observatory for Infrared Astronomy (SOFIA) is a Boeing 747SP aircraft modified to accommodate a 2.7 meter gyro-stabilized telescope, which is mainly focused to studying the Universe at infrared wavelengths. As part of the Strategic Director's Discretionary Time (S-DDT) program, SOFIA performs observations of relevant science cases and immediately offers science-ready data products to the astronomical community. We present the first data release of the S-DDT program on far-infrared imaging polarimetric observations of 30 Doradus using the High-resolution Airborne Wideband Camera-Plus (HAWC+) at 53, 89, 154, and 214 micron. We present the status and quality of the observations, an overview of the SOFIA data products, and examples of working with HAWC+ polarimetric data that will enhance the scientific analysis of this, and future, data sets. These observations illustrate the potential influence of magnetic fields and turbulence in a star-forming region within the Tarantula Nebula. △ Less

Submitted 7 November, 2018; originally announced November 2018.

Comments: SOFIA Science Center White Paper, 20 pages, 14 figures, 3 tables