Showing 1–48 of 48 results for author: Herrero-Illana, R

An ALCHEMI inspection of sulphur-bearing species towards the central molecular zone of NGC 253

Authors: M. Bouvier, S. Viti, E. Behrens, J. Butterworth, K. -Y. Huang, J. G. Mangum, N. Harada, S. Martín, V. M. Rivilla, S. Muller, K. Sakamoto, Y. Yoshimura, K. Tanaka, K. Nakanishi, R. Herrero-Illana, L. Colzi, M. D. Gorski, C. Henkel, P. K. Humire, D. S. Meier, P. P. van der Werf, Y. T. Yan

Abstract: Sulphur-bearing species are detected in various environments within Galactic star-forming regions and are particularly abundant in the gas phase of outflows and shocks, and photo-dissociation regions. In this work, we aim to investigate the nature of the emission from the most common sulphur-bearing species observable at millimetre wavelengths towards the nuclear starburst of the galaxy NGC 253. W… ▽ More Sulphur-bearing species are detected in various environments within Galactic star-forming regions and are particularly abundant in the gas phase of outflows and shocks, and photo-dissociation regions. In this work, we aim to investigate the nature of the emission from the most common sulphur-bearing species observable at millimetre wavelengths towards the nuclear starburst of the galaxy NGC 253. We intend to understand which type of regions are probed by sulphur-bearing species and which process(es) dominate(s) the release of sulphur into the gas phase. We used the high-angular resolution (1.6" or 27 pc) observations from the ALCHEMI ALMA Large Program to image several sulphur-bearing species towards the central molecular zone (CMZ) of NGC 253. We performed local thermodynamic equilibrium (LTE) and non-LTE large velocity gradient (LVG) analyses to derive the physical conditions of the gas in which S-bearing species are emitted, and their abundance ratios across the CMZ. Finally, we compared our results with previous ALCHEMI studies and a few selected Galactic environments. We found that not all sulphur-bearing species trace the same type of gas: strong evidence indicates that H2S and part of the emission of OCS, H2CS, and SO, are tracing shocks whilst part of SO and CS emission rather trace the dense molecular gas. For some species, such as CCS and SO2, we could not firmly conclude on their origin of emission. The present analysis indicates that the emission from most sulphur-bearing species throughout the CMZ is likely dominated by shocks associated with ongoing star formation. In the inner part of the CMZ where the presence of super star clusters was previously indicated, we could not distinguish between shocks or thermal evaporation as the main process releasing the S-bearing species. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 44 pages, 20 figures, Accepted for publication in A&A

Molecular isotopologue measurements toward super star clusters and the relation to their ages in NGC253 with ALCHEMI

Authors: J. Butterworth, S. Viti, P. P. Van der Werf, J. G. Mangum, S. Martín, N. Harada, K. L. Emig, S. Muller, K. Sakamoto, Y. Yoshimura, K. Tanaka, R. Herrero-Illana, L. Colzi, V. M. Rivilla, K. Y. Huang, M. Bouvier, E. Behrens, C. Henkel, Y. T. Yan, D. S. Meier, D. Zhou

Abstract: Determining the evolution of the CNO isotopes in the interstellar medium (ISM) of starburst galaxies can yield important constraints on the ages of superstar clusters (SSCs), or on other aspects and contributing factors of their evolution. Due to the time-dependent nature of the abundances of isotopes within the ISM as they are supplied from processes such as nucleosynthesis or chemical fractionat… ▽ More Determining the evolution of the CNO isotopes in the interstellar medium (ISM) of starburst galaxies can yield important constraints on the ages of superstar clusters (SSCs), or on other aspects and contributing factors of their evolution. Due to the time-dependent nature of the abundances of isotopes within the ISM as they are supplied from processes such as nucleosynthesis or chemical fractionation, this provides the possible opportunity to probe the ability of isotopes ratios to trace the ages of high star forming regions, such as SSCs. The goal of this study is to investigate whether the isotopic variations in SSC regions within NGC253 are correlated with their different ages as derived from stellar population modelling. We have measured abundance ratios of CO, HCN and HCO$^+$ isotopologues in six regions containing SSCs within NGC253 using high spatial resolution (1.6",$\sim 28$pc) data from the ALCHEMI (ALma Comprehensive High-resolution Extragalactic Molecular Inventory) ALMA Large program. We have then analysed these ratios using RADEX radiative transfer modelling, with the parameter space sampled using the nested sampling Monte Carlo algorithm MLFriends. These abundance ratios were then compared to ages predicted in each region via the fitting of observed star formation tracers (such as Br$γ$) to starburst stellar population evolution models. We do not find any significant trend with age for the CO and HCN isotopologue ratios on the timescales for the ages of the SSC* regions observed. The driving factors of these ratios within SSCs could be the Initial Mass Function as well as possibly fractionation effects. To further probe these effects in SSCs over time a larger sample of SSCs must be observed spanning a larger age range. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: 44 pages, 43 Figures, Accepted for Publication to A&A

Journal ref: A&A 686, A31 (2024)

The ALCHEMI atlas: principal component analysis reveals starburst evolution in NGC 253

Authors: Nanase Harada, David S. Meier, Sergio Martín, Sebastien Muller, Kazushi Sakamoto, Toshiki Saito, Mark D. Gorski, Christian Henkel, Kunihiko Tanaka, Jeffrey G. Mangum, Susanne Aalto, Rebeca Aladro, Mathilde Bouvier, Laura Colzi, Kimberly L. Emig, Rubén Herrero-Illana, Ko-Yun Huang, Kotaro Kohno, Sabine König, Kouichiro Nakanishi, Yuri Nishimura, Shuro Takano, Víctor M. Rivilla, Serena Viti, Yoshimasa Watanabe , et al. (2 additional authors not shown)

Abstract: Molecular lines are powerful diagnostics of the physical and chemical properties of the interstellar medium (ISM). These ISM properties, which affect future star formation, are expected to differ in starburst galaxies from those of more quiescent galaxies. We investigate the ISM properties in the central molecular zone of the nearby starburst galaxy NGC 253 using the ultra-wide millimeter spectral… ▽ More Molecular lines are powerful diagnostics of the physical and chemical properties of the interstellar medium (ISM). These ISM properties, which affect future star formation, are expected to differ in starburst galaxies from those of more quiescent galaxies. We investigate the ISM properties in the central molecular zone of the nearby starburst galaxy NGC 253 using the ultra-wide millimeter spectral scan survey from the ALMA Large Program ALCHEMI. We present an atlas of velocity-integrated images at a 1".6 resolution of 148 unblended transitions from 44 species, including the first extragalactic detection of HCNH$^+$ and the first interferometric images of C$_3$H$^+$, NO, HCS$^+$. We conduct a principal component analysis (PCA) on these images to extract correlated chemical species and to identify key groups of diagnostic transitions. To the best of our knowledge, our dataset is currently the largest astronomical set of molecular lines to which PCA has been applied. The PCA can categorize transitions coming from different physical components in NGC 253 such as i) young starburst tracers characterized by high-excitation transitions of HC$_3$N and complex organic molecules (COMs) versus tracers of on-going star formation (radio recombination lines) and high-excitation transitions of CCH and CN tracing PDRs, ii) tracers of cloud-collision-induced shocks (low-excitation transitions of CH$_3$OH, HNCO, HOCO$^+$, and OCS) versus shocks from star-formation-induced outflows (high-excitation transitions of SiO), as well as iii) outflows showing emission from HOC$^+$, CCH, H$_3$O$^+$, CO isotopologues, HCN, HCO$^+$, CS, and CN. Our findings show these intensities vary with galactic dynamics, star formation activities, and stellar feedback. △ Less

Submitted 4 January, 2024; originally announced January 2024.

Comments: 65 pages, 39 figures. Accepted for publication in ApJS

A ring-like accretion structure in M87 connecting its black hole and jet

Authors: Ru-Sen Lu, Keiichi Asada, Thomas P. Krichbaum, Jongho Park, Fumie Tazaki, Hung-Yi Pu, Masanori Nakamura, Andrei Lobanov, Kazuhiro Hada, Kazunori Akiyama, Jae-Young Kim, Ivan Marti-Vidal, José L. Gómez, Tomohisa Kawashima, Feng Yuan, Eduardo Ros, Walter Alef, Silke Britzen, Michael Bremer, Avery E. Broderick, Akihiro Doi, Gabriele Giovannini, Marcello Giroletti, Paul T. P. Ho, Mareki Honma , et al. (96 additional authors not shown)

Abstract: The nearby radio galaxy M87 is a prime target for studying black hole accretion and jet formation^{1,2}. Event Horizon Telescope observations of M87 in 2017, at a wavelength of 1.3 mm, revealed a ring-like structure, which was interpreted as gravitationally lensed emission around a central black hole^3. Here we report images of M87 obtained in 2018, at a wavelength of 3.5 mm, showing that the comp… ▽ More The nearby radio galaxy M87 is a prime target for studying black hole accretion and jet formation^{1,2}. Event Horizon Telescope observations of M87 in 2017, at a wavelength of 1.3 mm, revealed a ring-like structure, which was interpreted as gravitationally lensed emission around a central black hole^3. Here we report images of M87 obtained in 2018, at a wavelength of 3.5 mm, showing that the compact radio core is spatially resolved. High-resolution imaging shows a ring-like structure of 8.4_{-1.1}^{+0.5} Schwarzschild radii in diameter, approximately 50% larger than that seen at 1.3 mm. The outer edge at 3.5 mm is also larger than that at 1.3 mm. This larger and thicker ring indicates a substantial contribution from the accretion flow with absorption effects in addition to the gravitationally lensed ring-like emission. The images show that the edge-brightened jet connects to the accretion flow of the black hole. Close to the black hole, the emission profile of the jet-launching region is wider than the expected profile of a black-hole-driven jet, suggesting the possible presence of a wind associated with the accretion flow. △ Less

Submitted 25 April, 2023; originally announced April 2023.

Comments: 50 pages, 18 figures, 3 tables, author's version of the paper published in Nature

Reconstructing the shock history in the CMZ of NGC 253 with ALCHEMI

Authors: K. -Y. Huang, S. Viti, J. Holdship, J. G. Mangum, S. Martín, N. Harada, S. Muller, K. Sakamoto, K. Tanaka, Y. Yoshimura, R. Herrero-Illana, D. S. Meier, E. Behrens, P. P. van der Werf, C. Henkel, S. García-Burillo, V. M. Rivilla, K. L. Emig, L. Colzi, P. K. Humire, R. Aladro, M. Bouvier

Abstract: HNCO and SiO are well known shock tracers and have been observed in nearby galaxies, including the nearby (D=3.5 Mpc) starburst galaxy NGC 253. The simultaneous detection of these two species in regions where the star formation rate is high may be used to study the shock history of the gas. We perform a multi-line molecular study using these two shock tracers (SiO and HNCO) with the aim of charact… ▽ More HNCO and SiO are well known shock tracers and have been observed in nearby galaxies, including the nearby (D=3.5 Mpc) starburst galaxy NGC 253. The simultaneous detection of these two species in regions where the star formation rate is high may be used to study the shock history of the gas. We perform a multi-line molecular study using these two shock tracers (SiO and HNCO) with the aim of characterizing the gas properties. We also explore the possibility of reconstructing the shock history in NGC 253's Central Molecular Zone (CMZ). Six SiO transitions and eleven HNCO transitions were imaged at high resolution $1''.6$ (28 pc) with the Atacama Large Millimeter/submillimeter Array (ALMA) as part of the ALCHEMI Large Programme. Both non-LTE radiative transfer analysis and chemical modelling were performed in order to characterize the gas properties, and to investigate the chemical origin of the emission. The non-LTE radiative transfer analysis coupled with Bayesian inference shows clear evidence that the gas traced by SiO has different densities and temperatures than that traced by HNCO, with an indication that shocks are needed to produce both species. Chemical modelling further confirms such a scenario and suggests that fast and slow shocks are responsible for SiO and HNCO production, respectively, in most GMCs. We are also able to infer the physical characteristics of the shocks traced by SiO and HNCO for each GMC. Radiative transfer and chemical analysis of the SiO and HNCO in the CMZ of NGC 253 reveal a complex picture whereby most of the GMCs are subjected to shocks. We speculate on the possible shock scenarios responsible for the observed emission and provide potential history and timescales for each shock scenario. Higher spatial resolution observations of these two species are required in order to quantitatively differentiate between scenarios. △ Less

Submitted 22 March, 2023; originally announced March 2023.

Journal ref: A&A 675, A151 (2023)

Tracing Interstellar Heating: An ALCHEMI Measurement of the HCN Isomers in NGC 253

Authors: Erica Behrens, Jeffrey G. Mangum, Jonathan Holdship, Serena Viti, Nanase Harada, Sergio Martin, Kazushi Sakamoto, Sebastien Muller, Kunihiko Tanaka, Kouichiro Nakanishi, Ruben Herrero-Illana, Yuki Yoshimura, Rebeca Aladro, Laura Colzi, Kimberly L. Emig, Christian Henkel, Ko-Yun Huang, P. K. Humire, David S. Meier, Victor M. Rivilla

Abstract: We analyze HCN and HNC emission in the nearby starburst galaxy NGC 253 to investigate its effectiveness in tracing heating processes associated with star formation. This study uses multiple HCN and HNC rotational transitions observed using ALMA via the ALCHEMI Large Program. To understand the conditions and associated heating mechanisms within NGC 253's dense gas, we employ Bayesian nested samplin… ▽ More We analyze HCN and HNC emission in the nearby starburst galaxy NGC 253 to investigate its effectiveness in tracing heating processes associated with star formation. This study uses multiple HCN and HNC rotational transitions observed using ALMA via the ALCHEMI Large Program. To understand the conditions and associated heating mechanisms within NGC 253's dense gas, we employ Bayesian nested sampling techniques applied to chemical and radiative transfer models which are constrained using our HCN and HNC measurements. We find that the volume density $n_{\text{H}_{2}}$ and cosmic ray ionization rate (CRIR) $ζ$ are enhanced by about an order of magnitude in the galaxy's central regions as compared to those further from the nucleus. In NGC 253's central GMCs, where observed HCN/HNC abundance ratios are lowest, $n \sim 10^{5.5}$ cm$^{-3}$ and $ζ\sim 10^{-12}$ s$^{-1}$ (greater than $10^4$ times the average Galactic rate). We find a positive correlation in the association of both density and CRIR with the number of star formation-related heating sources (supernova remnants, HII regions, and super hot cores) located in each GMC, as well as a correlation between CRIRs and supernova rates. Additionally, we see an anticorrelation between the HCN/HNC ratio and CRIR, indicating that this ratio will be lower in regions where $ζ$ is higher. Though previous studies suggested HCN and HNC may reveal strong mechanical heating processes in NGC 253's CMZ, we find cosmic ray heating dominates the heating budget, and mechanical heating does not play a significant role in the HCN and HNC chemistry. △ Less

Submitted 8 November, 2022; v1 submitted 13 September, 2022; originally announced September 2022.

Comments: 33 pages, 23 figures, accepted for publication by the Astrophysical Journal

ALCHEMI finds a "shocking" carbon footprint in the starburst galaxy NGC~253

Authors: Nanase Harada, Sergio Martin, Jeff Mangum, Kazushi Sakamoto, Sebastian Muller, Victor Rivilla, Christian Henkel, David Meier, Laura Colzi, Mitsuyoshi Yamagishi, Kunihiko Tanaka, Kouichiro Nakanishi, Ruben Herrero-Illana, Yuki Yoshimura, Pedro Humire, Rebeca Aladro, Paul van der Werf, Kim Emig

Abstract: Centers of starburst galaxies may be characterized by a specific gas and ice chemistry due to their gas dynamics and the presence of various ice desorption mechanisms. This may result in a peculiar observable composition. We analyze abundances of $CO_2$, a reliable tracer of ice chemistry, from data collected as part of the ALMA large program ALCHEMI, a wide-frequency spectral scan toward the star… ▽ More Centers of starburst galaxies may be characterized by a specific gas and ice chemistry due to their gas dynamics and the presence of various ice desorption mechanisms. This may result in a peculiar observable composition. We analyze abundances of $CO_2$, a reliable tracer of ice chemistry, from data collected as part of the ALMA large program ALCHEMI, a wide-frequency spectral scan toward the starburst galaxy NGC~253 with an angular resolution of 1.6$''$. We constrain the $CO_2$ abundances in the gas phase using its protonated form $HOCO^+$. The distribution of $HOCO^+$ is similar to that of methanol, which suggests that $HOCO^+$ is indeed produced from the protonation of $CO_2$ sublimated from ice. The $HOCO^+$ fractional abundances are found to be $(1-2)\times10^{-9}$ at the outer part of the central molecular zone (CMZ), while they are lower ($\sim10^{-10}$) near the kinematic center. This peak fractional abundance at the outer CMZ is comparable to that in the Milky Way CMZ, and orders of magnitude higher than that in Galactic disk star-forming regions. From the range of $HOCO^+/CO_2$ ratios suggested from chemical models, the gas-phase $CO_2$ fractional abundance is estimated to be $(1-20)\times10^{-7}$ at the outer CMZ, and orders of magnitude lower near the center. We estimate the $CO_2$ ice fractional abundances at the outer CMZ to be $(2-5)\times10^{-6}$ from the literature. A comparison between the ice and gas $CO_2$ abundances suggests an efficient sublimation mechanism. This sublimation is attributed to large-scale shocks at the orbital intersections of the bar and CMZ. △ Less

Submitted 30 August, 2022; originally announced August 2022.

Comments: 22 pages, 9 figures. Accepted for publication in the Astrophysical Journal

Methanol masers in NGC 253 with ALCHEMI

Authors: P. K. Humire, C. Henkel, A. Hernández-Gómez, S. Martín, J. Mangum, N. Harada, S. Muller, K. Sakamoto, K. Tanaka, Y. Yoshimura, K. Nakanishi, S. Mühle, R. Herrero-Illana, D. S. Meier, E. Caux, R. Aladro, R. Mauersberger, S. Viti, L. Colzi, V. M. Rivilla, M. Gorski, K. M. Menten, K. -Y. Huang, S. Aalto, P. P. van der Werf , et al. (1 additional authors not shown)

Abstract: Context: Methanol masers of Class I (collisionally-pumped) and Class II (radiatively-pumped) have been studied in great detail in our Galaxy in a variety of astrophysical environments such as shocks and star-forming regions and are helpful to analyze the properties of the dense interstellar medium. However, the study of methanol masers in external galaxies is still in its infancy. Aims: Our main g… ▽ More Context: Methanol masers of Class I (collisionally-pumped) and Class II (radiatively-pumped) have been studied in great detail in our Galaxy in a variety of astrophysical environments such as shocks and star-forming regions and are helpful to analyze the properties of the dense interstellar medium. However, the study of methanol masers in external galaxies is still in its infancy. Aims: Our main goal is to search for methanol masers in the central molecular zone (CMZ; inner 500 pc) of the nearby starburst galaxy NGC 253. Methods: Covering a frequency range between 84 and 373 GHz ($λ$ = 3.6 to 0.8 mm) at high angular (1.6"$\sim$27 pc) and spectral ($\sim$8--9 km s$^{-1}$) resolution with the ALMA large program ALCHEMI, we have probed different regions across the CMZ of NGC 253. In order to look for methanol maser candidates, we employed the rotation diagram method and a set of radiative transfer models. Results: We detect for the first time masers above 84 GHz in NGC 253, covering an ample portion of the $J_{-1}\rightarrow(J-$ 1)$_{0}-E$ line series (at 84, 132, 229, and 278 GHz) and the $J_{0}\rightarrow(J-$ 1)$_{1}-A$ series (at 95, 146, and 198 GHz). This confirms the presence of the Class I maser line at 84 GHz, already reported but now being detected in more than one location. For the $J_{-1}\rightarrow(J-$ 1)$_{0}-E$ line series, we observe a lack of Class I maser candidates in the central star-forming disk. Conclusions: The physical conditions for maser excitation in the $J_{-1}\rightarrow(J-$ 1)$_{0}-E$ line series can be weak shocks and cloud-cloud collisions as suggested by shock tracers (SiO and HNCO) in bi-symmetric shock/active regions located in the outskirts of the CMZ. On the other hand, the presence of photodissociation regions due to a high star-formation rate would be needed to explain the lack of Class I masers in the very central regions. △ Less

Submitted 6 May, 2022; originally announced May 2022.

Comments: Accepted for publication in A&A. 29 pages, 17 figures (4 in Appendix)

Energizing Star Formation: The Cosmic Ray Ionization Rate in NGC 253 Derived From ALCHEMI Measurements of H$_3$O$^+$ and SO

Authors: Jonathan Holdship, Jeffrey G. Mangum, Serena Viti, Erica Behrens, Nanase Harada, Sergio Martín, Kazushi Sakamoto, Sebastien Muller, Kunihiko Tanaka, Kouichiro Nakanishi, Rubén Herrero-Illana, Yuki Yoshimura, Rebeca Aladro, Laura Colzi, Kimberly L. Emig, Christian Henkel, Yuri Nishimura, Víctor M. Rivilla, Paul P. van der Werf

Abstract: The cosmic ray ionization rate (CRIR) is a key parameter in understanding the physical and chemical processes in the interstellar medium. Cosmic rays are a significant source of energy in star formation regions, which impacts the physical and chemical processes which drive the formation of stars. Previous studies of the circum-molecular zone (CMZ) of the starburst galaxy NGC 253 have found evidenc… ▽ More The cosmic ray ionization rate (CRIR) is a key parameter in understanding the physical and chemical processes in the interstellar medium. Cosmic rays are a significant source of energy in star formation regions, which impacts the physical and chemical processes which drive the formation of stars. Previous studies of the circum-molecular zone (CMZ) of the starburst galaxy NGC 253 have found evidence for a high CRIR value; $10^3-10^6$ times the average cosmic ray ionization rate within the Milky Way. This is a broad constraint and one goal of this study is to determine this value with much higher precision. We exploit ALMA observations towards the central molecular zone of NGC 253 to measure the CRIR. We first demonstrate that the abundance ratio of H$_3$O$^+$ and SO is strongly sensitive to the CRIR. We then combine chemical and radiative transfer models with nested sampling to infer the gas properties and CRIR of several star-forming regions in NGC 253 due to emission from their transitions. We find that each of the four regions modelled has a CRIR in the range $(1-80)\times10^{-14}$ s$^{-1}$ and that this result adequately fits the abundances of other species that are believed to be sensitive to cosmic rays including C$_2$H, HCO$^+$, HOC$^+$, and CO. From shock and PDR/XDR models, we further find that neither UV/X-ray driven nor shock dominated chemistry are a viable single alternative as none of these processes can adequately fit the abundances of all of these species. △ Less

Submitted 11 May, 2022; v1 submitted 7 April, 2022; originally announced April 2022.

Comments: 24 pages, 15 figures, accepted for publication in ApJ, updated to include changes made for proofs

ALCHEMI: an ALMA Comprehensive High-resolution Extragalactic Molecular Inventory. Survey presentation and first results from the ACA array

Authors: S. Martín, J. G. Mangum, N. Harada, F. Costagliola, K. Sakamoto, S. Muller, R. Aladro, K. Tanaka, Y. Yoshimura, K. Nakanishi, R. Herrero-Illana, S. Mühle, S. Aalto, E. Behrens, L. Colzi, K. L. Emig, G. A. Fuller, S. García-Burillo, T. R. Greve, C. Henkel, J. Holdship, P. Humire, L. Hunt, T. Izumi, K. Kohno , et al. (10 additional authors not shown)

Abstract: We used the Atacama Large Millimeter/submillimeter Array (ALMA), covering a nearly contiguous 289 GHz frequency range between 84.2 and 373.2 GHz, to image the continuum and spectral line emission at 1.6\arcsec ($\sim 28$ pc) resolution down to a sensitivity of $30-50$ mK. This article describes the ALMA Comprehensive High-resolution Extragalactic Molecular Inventory (ALCHEMI) Large Program. We foc… ▽ More We used the Atacama Large Millimeter/submillimeter Array (ALMA), covering a nearly contiguous 289 GHz frequency range between 84.2 and 373.2 GHz, to image the continuum and spectral line emission at 1.6\arcsec ($\sim 28$ pc) resolution down to a sensitivity of $30-50$ mK. This article describes the ALMA Comprehensive High-resolution Extragalactic Molecular Inventory (ALCHEMI) Large Program. We focus on the analysis of the spectra extracted from the $15''$ ($\sim255$ pc) resolution ALMA Compact Array data. We model the molecular emission assuming local thermodynamic equilibrium with 78 species detected. Additionally, multiple hydrogen and helium recombination lines are identified. Spectral lines contribute 5 to 36\% of the total emission in frequency bins of 50 GHz. We report the first extragalactic detections of C$_2$H$_5$OH, HOCN, HC$_3$HO, and several rare isotopologues. Isotopic ratios of carbon, oxygen, sulfur, nitrogen and silicon were measure with multiple species. Infrared pumped vibrationaly excited HCN, HNC, and HC$_3$N emission, originating in massive star formation locations, is clearly detected at low resolution, while we do not detect it for HCO$^+$. We suggest high temperature conditions in these regions driving a seemingly "carbon-rich" chemistry which may also explain the observed high abundance of organic species close to those in Galactic hot cores. The $L_{vib}/L_{IR}$ ratio is used as a proxy to estimate a $3\%$ contribution from proto super star cluster to the global infrared emission. Measured isotopic ratios with high dipole moment species agree with those within the central kiloparsec of the Galaxy, while those derived from $\rm^{13}C^{18}O$ are a factor of 5 larger, confirming the existence of multiple ISM components within NGC 253 with different degrees of nucleosynthesis enrichment. ALCHEMI provides a template for early Universe galaxies. △ Less

Submitted 17 September, 2021; originally announced September 2021.

Comments: 60 pages, 32 figures, 15 tables, Accepted for publication at A&A

Journal ref: A&A 656, A46 (2021)

Starburst Energy Feedback Seen Through HCO$^+$/HOC$^+$ Emission in NGC 253 from ALCHEMI

Authors: Nanase Harada, Sergio Martín, Jeffrey G. Mangum, Kazushi Sakamoto, Sebastien Muller, Kunihiko Tanaka, Kouichiro Nakanishi, Rubén Herrero-Illana, Yuki Yoshimura, Stefanie Mühle, Rebeca Aladro, Laura Colzi, Víctor M. Rivilla, Susanne Aalto, Erica Behrens, Christian Henkel, Jonathan Holdship, P. K. Humire, David S. Meier, Yuri Nishimura, Paul P. van der Werf, Serena Viti

Abstract: Molecular abundances are sensitive to UV-photon flux and cosmic-ray ionization rate. In starburst environments, the effects of high-energy photons and particles are expected to be stronger. We examine these astrochemical signatures through multiple transitions of HCO$^+$ and its metastable isomer HOC$^+$ in the center of the starburst galaxy NGC 253 using data from the ALMA large program ALCHEMI.… ▽ More Molecular abundances are sensitive to UV-photon flux and cosmic-ray ionization rate. In starburst environments, the effects of high-energy photons and particles are expected to be stronger. We examine these astrochemical signatures through multiple transitions of HCO$^+$ and its metastable isomer HOC$^+$ in the center of the starburst galaxy NGC 253 using data from the ALMA large program ALCHEMI. The distribution of the HOC$^+$(1-0) integrated intensity shows its association with "superbubbles", cavities created either by supernovae or expanding HII regions. The observed HCO$^+$/HOC$^+$ abundance ratios are $\sim 10-150$, and the fractional abundance of HOC$^+$ relative to H$_2$ is $\sim 1.5\times 10^{-11} - 6\times 10^{-10}$, which implies that the HOC$^+$ abundance in the center of NGC 253 is significantly higher than in quiescent spiral-arm dark clouds in the Galaxy and the Galactic center clouds. Comparison with chemical models implies either an interstellar radiation field of $G_0\gtrsim 10^3$ if the maximum visual extinction is $\gtrsim 5$, or a cosmic-ray ionization rate of $ζ\gtrsim 10^{-14}$ s$^{-1}$ (3-4 orders of magnitude higher than that within clouds in the Galactic spiral-arms) to reproduce the observed results. From the difference in formation routes of HOC$^+$, we propose that a low-excitation line of HOC$^+$ traces cosmic-ray dominated regions, while high-excitation lines trace photodissociation regions. Our results suggest that the interstellar medium in the center of NGC 253 is significantly affected by energy input from UV-photons and cosmic rays, sources of energy feedback. △ Less

Submitted 14 September, 2021; originally announced September 2021.

Comments: 57 pages, 32 figures, accepted for publication in ApJ

Subarcsecond LOFAR imaging of Arp299 at 150 MHz. Tracing the nuclear and diffuse extended emission of a bright LIRG

Authors: Naím Ramírez-Olivencia, Eskil Varenius, Miguel Pérez-Torres, Antonio Alberdi, John Conway, Almudena Alonso-Herrero, Miguel Pereira-Santaella, Rubén Herrero-Illana

Abstract: We study for the first time the low-frequency ($\sim$150 MHz) radio brightness distribution of Arp~299 at subarcsecond resolution, tracing in both compact and extended emission regions the local spectral energy distribution (SED) in order to characterize the dominant emission and absorption processes. We analysed the spatially resolved emission of Arp 299 revealed by 150 MHz international baseline… ▽ More We study for the first time the low-frequency ($\sim$150 MHz) radio brightness distribution of Arp~299 at subarcsecond resolution, tracing in both compact and extended emission regions the local spectral energy distribution (SED) in order to characterize the dominant emission and absorption processes. We analysed the spatially resolved emission of Arp 299 revealed by 150 MHz international baseline Low-Frequency Array (LOFAR) and 1.4, 5.0, and 8.4 GHz Very Large Array (VLA) observations. We present the first subarcsecond (0.4"$\sim$100~pc) image of the whole Arp~299 system at 150~MHz. The high surface brightness sensitivity of our LOFAR observations ($\sim$100 $μ$Jy/beam) allowed us to detect all of the nuclear components detected at higher frequencies, as well as the extended steep-spectrum emission surrounding the nuclei. We obtained spatially resolved, two-point spectral index maps for the whole galaxy: the compact nuclei show relatively flat spectra, while the extended, diffuse component shows a steep spectrum. We fitted the radio SED of the nuclear regions using two different models: a continuous free-free medium model and a clumpy model. The continuous model can explain the SED of the nuclei assuming a population of relativistic electrons subjected to synchrotron, bremsstrahlung, and ionization losses. The clumpy model fits assuming relativistic electrons with negligible energy losses, and thermal fractions that are more typical of star-forming galaxies than those required for the continuous model. Our results confirm the usefulness of combining spatially resolved radio imaging at both MHz and GHz frequencies to characterize in detail the radio emission properties of LIRGs from the central 100 pc out to the kiloparsec galaxy-wide scales. △ Less

Submitted 16 August, 2021; originally announced August 2021.

Comments: Accepted / Submitted to a special issue of A&A on sub-arcsecond imaging with LOFAR

Journal ref: A&A 658, A4 (2022)

The Distribution and Origin of C$_2$H in NGC 253 from ALCHEMI

Authors: J. Holdship, S. Viti, S. Martín, N. Harada, J. Mangum, K. Sakamoto, S. Muller, K. Tanaka, Y. Yoshimura, K. Nakanishi, R. Herrero-Illana, S. Mühle, R. Aladro, L. Colzi, K. L. Emig, S. García-Burillo, C. Henkel, P. Humire, D. S. Meier, V. M. Rivilla, P. van der Werf

Abstract: Observations of chemical species can provide an insight into the physical conditions of the emitting gas but it is important to understand how their abundances and excitation vary within different heating environments. C$_2$H is a molecule typically found in PDR regions of our own Galaxy but there is evidence to suggest it also traces other regions undergoing energetic processing in extragalactic… ▽ More Observations of chemical species can provide an insight into the physical conditions of the emitting gas but it is important to understand how their abundances and excitation vary within different heating environments. C$_2$H is a molecule typically found in PDR regions of our own Galaxy but there is evidence to suggest it also traces other regions undergoing energetic processing in extragalactic environments. As part of the ALCHEMI ALMA large program, the emission of C$_2$H in the central molecular zone of the nearby starburst galaxy NGC 253 was mapped at 1.6 " (28 pc) resolution and characterized to understand its chemical origins. Spectral modelling of the N=1-0 through N=4-3 rotational transitions of C$_2$H was used to derive the C$_2$H column densities towards the dense clouds in NGC 253. Chemical modelling, including PDR, dense cloud, and shock models were then used to investigate the chemical processes and physical conditions that are producing the molecular emission. We find high C$_2$H column densities of $\sim 10^{15} cm^{-3}$ detected towards the dense regions of NGC 253. We further find that these column densities cannot be reproduced by assuming that the emission arises from the PDR regions at the edge of the clouds. Instead, we find that the C$_2$H abundance remains high even in the high visual extinction interior of these clouds and that this is most likely caused by a high cosmic-ray ionization rate. △ Less

Submitted 16 July, 2021; v1 submitted 9 July, 2021; originally announced July 2021.

Comments: 18 pages, 13 figures, accepted for publication in A&A

Journal ref: A&A 654, A55 (2021)

Polarimetric properties of Event Horizon Telescope targets from ALMA

Authors: Ciriaco Goddi, Ivan Marti-Vidal, Hugo Messias, Geoffrey C. Bower, Avery E. Broderick, Jason Dexter, Daniel P. Marrone, Monika Moscibrodzka, Hiroshi Nagai, Juan Carlos Algaba, Keiichi Asada, Geoffrey B. Crew, Jose L. Gomez, C. M. Violette Impellizzeri, Michael Janssen, Matthias Kadler, Thomas P. Krichbaum, Rocco Lico, Lynn D. Matthews, Antonios Nathanail, Angelo Ricarte, Eduardo Ros, Ziri Younsi, The Event Horizon Telescope Collaboration, Gabriele Bruni , et al. (9 additional authors not shown)

Abstract: We present the results from a full polarization study carried out with ALMA during the first VLBI campaign, which was conducted in Apr 2017 in the $λ$3mm and $λ$1.3mm bands, in concert with the Global mm-VLBI Array (GMVA) and the Event Horizon Telescope (EHT), respectively. We determine the polarization and Faraday properties of all VLBI targets, including Sgr A*, M87, and a dozen radio-loud AGN.… ▽ More We present the results from a full polarization study carried out with ALMA during the first VLBI campaign, which was conducted in Apr 2017 in the $λ$3mm and $λ$1.3mm bands, in concert with the Global mm-VLBI Array (GMVA) and the Event Horizon Telescope (EHT), respectively. We determine the polarization and Faraday properties of all VLBI targets, including Sgr A*, M87, and a dozen radio-loud AGN. We detect high linear polarization fractions (2-15%) and large rotation measures (RM $>10^{3.3}-10^{5.5}$ rad m$^{-2}$). For Sgr A* we report a mean RM of $(-4.2\pm0.3) \times10^5$ rad m$^{-2}$ at 1.3 mm, consistent with measurements over the past decade, and, for the first time, an RM of $(-2.1\pm0.1) \times10^5$ rad m$^{-2}$ at 3 mm, suggesting that about half of the Faraday rotation at 1.3 mm may occur between the 3 mm photosphere and the 1.3 mm source. We also report the first unambiguous measurement of RM toward the M87 nucleus at mm wavelengths, which undergoes significant changes in magnitude and sign reversals on a one year time-scale, spanning the range from -1.2 to 0.3 $\times\,10^5$ rad m$^{-2}$ at 3 mm and -4.1 to 1.5 $\times\,10^5$ rad m$^{-2}$ at 1.3 mm. Given this time variability, we argue that, unlike the case of Sgr A*, the RM in M87 does not provide an accurate estimate of the mass accretion rate onto the black hole. We put forward a two-component model, comprised of a variable compact region and a static extended region, that can simultaneously explain the polarimetric properties observed by both the EHT and ALMA. These measurements provide critical constraints for the calibration, analysis, and interpretation of simultaneously obtained VLBI data with the EHT and GMVA. △ Less

Submitted 5 May, 2021; originally announced May 2021.

Comments: 62 pages, 26 figures, published in The Astrophysical Journal Letters on March 24, 2021

Journal ref: The Astrophysical Journal Letters, Volume 910, Issue 1, id.L14, 54 pp., 2021 March 24

Radio observations of massive stars in the Galactic centre: The Arches Cluster

Authors: A. T. Gallego-Calvente, R. Schoedel, A. Alberdi, R. Herrero-Illana, F. Najarro, F. Yusef-Zadeh, H. Dong, J. Sanchez-Bermudez, B. Shahzamanian, F. Nogueras-Lara, E. Gallego-Cano

Abstract: We present high-angular-resolution radio observations of the Arches cluster in the Galactic centre, one of the most massive young clusters in the Milky Way. The data were acquired in two epochs and at 6 and 10 GHz with the Karl G. Jansky Very Large Array (JVLA). The rms noise reached is three to four times better than during previous observations and we have almost doubled the number of known radi… ▽ More We present high-angular-resolution radio observations of the Arches cluster in the Galactic centre, one of the most massive young clusters in the Milky Way. The data were acquired in two epochs and at 6 and 10 GHz with the Karl G. Jansky Very Large Array (JVLA). The rms noise reached is three to four times better than during previous observations and we have almost doubled the number of known radio stars in the cluster. Nine of them have spectral indices consistent with thermal emission from ionised stellar winds, one is a confirmed colliding wind binary (CWB), and two sources are ambiguous cases. Regarding variability, the radio emission appears to be stable on timescales of a few to ten years. Finally, we show that the number of radio stars can be used as a tool for constraining the age and/or mass of a cluster and also its mass function. △ Less

Submitted 13 January, 2021; originally announced January 2021.

Comments: Accepted for publication by A&A

Journal ref: A&A 647, A110 (2021)

Rapid Variability of Sgr A* across the Electromagnetic Spectrum

Authors: G. Witzel, G. Martinez, S. P. Willner, E. E. Becklin, 4 H. Boyce, T. Do, A. Eckart, G. G. Fazio, A. Ghez, M. A. Gurwell, D. Haggard, R. Herrero-Illana, J. L. Hora, Z. Li, J. Liu, N. Marchili, Mark R. Morris, Howard A. Smith, M. Subroweit, J. A. Zensus

Abstract: Sagittarius A* (Sgr A*) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We have analyzed a comprehensive submillimeter (including new observations simultaneous with NIR monitoring), NIR, and 2-8 keV dataset. Submillimeter variations tend to lag those in the NIR by $\sim$30 minutes. An approximate Bayesian computation (ABC)… ▽ More Sagittarius A* (Sgr A*) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We have analyzed a comprehensive submillimeter (including new observations simultaneous with NIR monitoring), NIR, and 2-8 keV dataset. Submillimeter variations tend to lag those in the NIR by $\sim$30 minutes. An approximate Bayesian computation (ABC) fit to the X-ray first-order structure function shows significantly less power at short timescales in the X-rays than in the NIR. Less X-ray variability at short timescales combined with the observed NIR-X-ray correlations means the variability can be described as the result of two strictly correlated stochastic processes, the X-ray process being the low-pass-filtered version of the NIR process. The NIR--X-ray linkage suggests a simple radiative model: a compact, self-absorbed synchrotron sphere with high-frequency cutoff close to NIR frequencies plus a synchrotron self-Compton scattering component at higher frequencies. This model, with parameters fit to the submillimeter, NIR, and X-ray structure functions, reproduces the observed flux densities at all wavelengths, the statistical properties of all light curves, and the time lags between bands. The fit also gives reasonable values for physical parameters such as magnetic flux density $B\approx13$ G, source size $L \approx2.2R_{S}$, and high-energy electron density $n_{e}\approx4\times10^{7}$ cm$^{-3}$. An animation illustrates typical light curves, and we make public the parameter chain of our Bayesian analysis, the model implementation, and the visualization code. △ Less

Submitted 7 June, 2021; v1 submitted 18 November, 2020; originally announced November 2020.

Comments: Accepted by ApJS, June 3, 2021. The journal version of figure 21 is animated. The animation can also be found here: https://doi.org/10.17617/1.kctx3s25. This version (version 2) has been revised according to the referee's suggestions and includes optimized figures and some changes and corrections of the text; no scientific conclusions have changed

A Hard X-ray Test of HCN Enhancements as a Tracer of Embedded Black Hole Growth

Authors: G. C. Privon, C. Ricci, S. Aalto, S. Viti, L. Armus, T. Díaz-Santos, E. González-Alfonso, K. Iwasawa, D. L. Jeff, E. Treister, F. Bauer, A. S. Evans, P. Garg, R. Herrero-Illana, J. M. Mazzarella, K. Larson, L. Blecha, L. Barcos-Muñoz, V. Charmandaris, S. Stierwalt, M. A. Pérez-Torres

Abstract: Enhanced emission from the dense gas tracer HCN (relative to HCO$^+ $) has been proposed as a signature of active galactic nuclei (AGN). In a previous single-dish millimeter line survey we identified galaxies with HCN/HCO$ ^+ $ (1-0) intensity ratios consistent with those of many AGN but whose mid-infrared spectral diagnostics are consistent with little to no ( $\lesssim15\% $) contribution of an… ▽ More Enhanced emission from the dense gas tracer HCN (relative to HCO$^+ $) has been proposed as a signature of active galactic nuclei (AGN). In a previous single-dish millimeter line survey we identified galaxies with HCN/HCO$ ^+ $ (1-0) intensity ratios consistent with those of many AGN but whose mid-infrared spectral diagnostics are consistent with little to no ( $\lesssim15\% $) contribution of an AGN to the bolometric luminosity. To search for putative heavily obscured AGN, we present and analyze \nustar hard X-ray (3-79 keV) observations of four such galaxies from the Great Observatories All-sky LIRG Survey. We find no X-ray evidence for AGN in three of the systems and place strong upper limits on the energetic contribution of any heavily obscured ($N_{\rm H}>10^{24}$ cm$^{-2}$) AGN to their bolometric luminosity. The X-ray flux upper limits are presently an order of magnitude below what XDR-driven chemistry model predict are necessary to drive HCN enhancements. In a fourth system we find a hard X-ray excess consistent with the presence of an AGN, but contributing only $\sim3\%$ of the bolometric luminosity. It is also unclear if the AGN is spatially associated with the HCN enhancement. We further explore the relationship between HCN/HCO$^+$ (for several $\mathrm{J}_\mathrm{upper}$ levels) and $L_\mathrm{AGN}/L_\mathrm{IR}$ for a larger sample of systems in the literature. We find no evidence for correlations between the line ratios and the AGN fraction derived from X-rays, indicating that HCN/HCO$^+$ intensity ratios are not driven by the energetic dominance of AGN, nor are they reliable indicators of whether SMBH accretion is ongoing. △ Less

Submitted 24 April, 2020; v1 submitted 6 April, 2020; originally announced April 2020.

Comments: 14 pages, 4 figures. accepted for publication in The Astrophysical Journal

Journal ref: ApJ 893 149 (2020)

The Quest for Dual and Binary Supermassive Black Holes: A Multi-Messenger View

Authors: Alessandra De Rosa, Cristian Vignali, Tamara Bogdanović, Pedro R. Capelo, Maria Charisi, Massimo Dotti, Bernd Husemann, Elisabeta Lusso, Lucio Mayer, Zsolt Paragi, Jessie Runnoe, Alberto Sesana, Lisa Steinborn, Stefano Bianchi, Monica Colpi, Luciano Del Valle, Sándor Frey, Krisztina É. Gabányi, Margherita Giustini, Matteo Guainazzi, Zoltan Haiman, Noelia Herrera Ruiz, Rubén Herrero-Illana, Kazushi Iwasawa, S. Komossa , et al. (5 additional authors not shown)

Abstract: The quest for binary and dual supermassive black holes (SMBHs) at the dawn of the multi-messenger era is compelling. Detecting dual active galactic nuclei (AGN) -- active SMBHs at projected separations larger than several parsecs -- and binary AGN -- probing the scale where SMBHs are bound in a Keplerian binary -- is an observational challenge. The study of AGN pairs (either dual or binary) also r… ▽ More The quest for binary and dual supermassive black holes (SMBHs) at the dawn of the multi-messenger era is compelling. Detecting dual active galactic nuclei (AGN) -- active SMBHs at projected separations larger than several parsecs -- and binary AGN -- probing the scale where SMBHs are bound in a Keplerian binary -- is an observational challenge. The study of AGN pairs (either dual or binary) also represents an overarching theoretical problem in cosmology and astrophysics. The AGN triggering calls for detailed knowledge of the hydrodynamical conditions of gas in the imminent surroundings of the SMBHs and, at the same time, their duality calls for detailed knowledge on how galaxies assemble through major and minor mergers and grow fed by matter along the filaments of the cosmic web. This review describes the techniques used across the electromagnetic spectrum to detect dual and binary AGN candidates and proposes new avenues for their search. The current observational status is compared with the state-of-the-art numerical simulations and models for formation of dual and binary AGN. Binary SMBHs are among the loudest sources of gravitational waves (GWs) in the Universe. The search for a background of GWs at nHz frequencies from inspiralling SMBHs at low redshifts, and the direct detection of signals from their coalescence by the Laser Interferometer Space Antenna in the next decade, make this a theme of major interest for multi-messenger astrophysics. This review discusses the future facilities and observational strategies that are likely to significantly advance this fascinating field. △ Less

Submitted 17 January, 2020; originally announced January 2020.

Comments: 69 pages, 21 figures. Accepted for publication in New Astronomy Reviews

Journal ref: New Astronomy Reviews, Vol. 86, id. 101525 (2019)

First M87 Event Horizon Telescope Results and the Role of ALMA

Authors: Ciriaco Goddi, Geoff Crew, Violette Impellizzeri, Ivan Marti-Vidal, Lynn D. Matthews, Hugo Messias, Helge Rottmann, Walter Alef, Lindy Blackburn, Thomas Bronzwaer, Chi-Kwan Chan, Jordy Davelaar, Roger Deane, Jason Dexter, Shep Doeleman, Heino Falcke, Vincent L. Fish, Raquel Fraga-Encinas, Christian M. Fromm, Ruben Herrero-Illana, Sara Issaoun, David James, Michael Janssen, Michael Kramer, Thomas P. Krichbaum , et al. (19 additional authors not shown)

Abstract: In April 2019, the Event Horizon Telescope (EHT) collaboration revealed the first image of the candidate super-massive black hole (SMBH) at the centre of the giant elliptical galaxy Messier 87 (M87). This event-horizon-scale image shows a ring of glowing plasma with a dark patch at the centre, which is interpreted as the shadow of the black hole. This breakthrough result, which represents a powerf… ▽ More In April 2019, the Event Horizon Telescope (EHT) collaboration revealed the first image of the candidate super-massive black hole (SMBH) at the centre of the giant elliptical galaxy Messier 87 (M87). This event-horizon-scale image shows a ring of glowing plasma with a dark patch at the centre, which is interpreted as the shadow of the black hole. This breakthrough result, which represents a powerful confirmation of Einstein's theory of gravity, or general relativity, was made possible by assembling a global network of radio telescopes operating at millimetre wavelengths that for the first time included the Atacama Large Millimeter/ submillimeter Array (ALMA). The addition of ALMA as an anchor station has enabled a giant leap forward by increasing the sensitivity limits of the EHT by an order of magnitude, effectively turning it into an imaging array. The published image demonstrates that it is now possible to directly study the event horizon shadows of SMBHs via electromagnetic radiation, thereby transforming this elusive frontier from a mathematical concept into an astrophysical reality. The expansion of the array over the next few years will include new stations on different continents - and eventually satellites in space. This will provide progressively sharper and higher-fidelity images of SMBH candidates, and potentially even movies of the hot plasma orbiting around SMBHs. These improvements will shed light on the processes of black hole accretion and jet formation on event-horizon scales, thereby enabling more precise tests of general relativity in the truly strong field regime. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: 11 pages + cover page, 6 figures

Journal ref: The Messenger, 177, 25 (2019)

Molecular gas and dust properties of galaxies from the Great Observatories All-sky LIRG Survey

Authors: R. Herrero-Illana, G. C. Privon, A. S. Evans, T. Díaz-Santos, M. Á. Pérez-Torres, V. U, A. Alberdi, K. Iwasawa, L. Armus, S. Aalto, J. Mazzarella, J. Chu, D. B. Sanders, L. Barcos-Muñoz, V. Charmandaris, S. T. Linden, I. Yoon, D. T. Frayer, H. Inami, D. -C. Kim, H. J. Borish, J. Conway, E. J Murphy, Y. Song, S. Stierwalt , et al. (1 additional authors not shown)

Abstract: We present IRAM-30m Telescope $^{12}$CO and $^{13}$CO observations of a sample of 55 luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) in the local universe. This sample is a subset of the Great Observatory All-Sky LIRG Survey (GOALS), for which we use ancillary multi-wavelength data to better understand their interstellar medium and star formation properties. Fifty-three (96%) of th… ▽ More We present IRAM-30m Telescope $^{12}$CO and $^{13}$CO observations of a sample of 55 luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) in the local universe. This sample is a subset of the Great Observatory All-Sky LIRG Survey (GOALS), for which we use ancillary multi-wavelength data to better understand their interstellar medium and star formation properties. Fifty-three (96%) of the galaxies are detected in $^{12}$CO, and 29 (52%) are also detected in $^{13}$CO above a 3$σ$ level. The median full width at zero intensity (FWZI) velocity of the CO line emission is 661km s$^{-1}$, and $\sim$54% of the galaxies show a multi-peak CO profile. Herschel photometric data is used to construct the far-IR spectral energy distribution of each galaxy, which are fit with a modified blackbody model that allows us to derive dust temperatures and masses, and infrared luminosities. We make the assumption that the gas-to-dust mass ratio of (U)LIRGs is comparable to local spiral galaxies with a similar stellar mass (i.e., gas/dust of mergers is comparable to their progenitors) to derive a CO-to-H$_2$ conversion factor of $\langleα\rangle=1.8^{+1.3}_{-0.8}M_\odot$(K km s$^{-1}$pc$^{2}$)$^{-1}$; such a value is comparable to that derived for (U)LIRGs based on dynamical mass arguments. We derive gas depletion times of $400-600$Myr for the (U)LIRGs, compared to the 1.3Gyr for local spiral galaxies. Finally, we re-examine the relationship between the $^{12}$CO/$^{13}$CO ratio and dust temperature, confirming a transition to elevated ratios in warmer systems. △ Less

Submitted 8 July, 2019; originally announced July 2019.

Comments: 32 pages, 10 figures. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 628, A71 (2019)

A VLBI study of the wind-wind collision region in the massive multiple HD 167971

Authors: J. Sanchez-Bermudez, A. Alberdi, R. Schödel, W. Brandner, R. Galván-Madrid, J. C. Guirado, R. Herrero-Illana, C. A. Hummel, J. M. Marcaide, M. A. Pérez-Torres

Abstract: Context. Colliding winds in massive binaries are able to accelerate particles up to relativistic speeds as the result of the interaction between the winds of the different stellar components. HD 167971 exhibits this phenomenology which makes it a strong radio source. Aims. We aim at characterizing the morphology of the radio emission and its dependence on the orbital motion, traced independently b… ▽ More Context. Colliding winds in massive binaries are able to accelerate particles up to relativistic speeds as the result of the interaction between the winds of the different stellar components. HD 167971 exhibits this phenomenology which makes it a strong radio source. Aims. We aim at characterizing the morphology of the radio emission and its dependence on the orbital motion, traced independently by NIR-interferometry, of the spectroscopic binary and the tertiary component that conforms HD 167971. Methods. We analyze 2006 and 2016 very long baseline interferometric data at C and X bands. We complement our analysis with a geometrical model of the wind-wind collision region, and with an astrometric description of the system. Results. We confirm that the detected non-thermal radio emission is associated with the wind-wind collision region of the spectroscopic binary and the tertiary component in HD 167971. The wind-wind collision region changes orientation in agreement with the orbital motion of the tertiary around the spectroscopic binary. The total intensity also changes between the two observing epochs in a way inversely proportional to the separation between the SB and T, with a negative-steep spectral index typical of an optically thin synchrotron emission possibly steepened by an inverse Compton cooling effect. The wind-wind collision bow-shock shape and its position with respect to the stars indicates that the wind momentum from the spectroscopic binary is stronger than that of the tertiary. Finally, the astrometric solution derived for the stellar system and the wind-wind collision region is consistent with independent Gaia data. △ Less

Submitted 29 January, 2019; originally announced January 2019.

Comments: Accepted for publication in Astronomy and Astrophysics, 7 pages, 6 figures

Journal ref: A&A 624, A55 (2019)

Calibration of ALMA as a phased array: ALMA observations during the 2017 VLBI campaign

Authors: Ciriaco Goddi, Ivan Marti-Vidal, Hugo Messias, Geoff Crew, Ruben Herrero-Illana, Violette Impellizzeri, Helge Rottmann, Jan Wagner, Ed Fomalont, Lynn D. Matthews, Dirk Petry, Neil Phillips, Remo Tilanus, Eric Villard, Lindy Blackburn, Michael Janßen, Maciek Wielgus

Abstract: We present a detailed description of the special procedures for calibration and quality assurance of Atacama Large Millimeter/submillimeter Array (ALMA) observations in Very Long Baseline Interferometry (VLBI) mode. These procedures are required to turn the phased ALMA array into a fully calibrated VLBI station. As an illustration of these methodologies, we present full-polarization observations c… ▽ More We present a detailed description of the special procedures for calibration and quality assurance of Atacama Large Millimeter/submillimeter Array (ALMA) observations in Very Long Baseline Interferometry (VLBI) mode. These procedures are required to turn the phased ALMA array into a fully calibrated VLBI station. As an illustration of these methodologies, we present full-polarization observations carried out with ALMA as a phased array at 3mm (Band 3) and 1.3mm (Band 6) as part of Cycle-4. These are the first VLBI science observations conducted with ALMA and were obtained during a 2017 VLBI campaign in concert with other telescopes worldwide as part of the Global mm-VLBI Array (GMVA, April 1-3) and the Event Horizon Telescope (EHT, April 5-11) in ALMA Bands 3 and 6, respectively. △ Less

Submitted 28 March, 2019; v1 submitted 28 January, 2019; originally announced January 2019.

Comments: 39 pages, 20 figures, 10 tables, accepted by PASP

A Serendipitous Hard X-ray Detection of the Blazar LBQS 1319+0039

Authors: G. C. Privon, C. Ricci, F. E. Bauer, M. A. Perez-Torres, R. Herrero-Illana, E. Treister, S. Aalto

Abstract: We report a serendipitous hard X-ray (3-24 keV) detection of the blazar LBQS 1319+0039. The X-ray spectrum is consistent with powerlaw emission with a photon index of $Γ= 1.78$. The observed emission and published redshift imply an observed 2-10 keV luminosity of $3.4\times10^{45}$ erg s$^{-1}$. This is consistent with extrapolations from a previous soft X-ray detection, thus there is no evidence… ▽ More We report a serendipitous hard X-ray (3-24 keV) detection of the blazar LBQS 1319+0039. The X-ray spectrum is consistent with powerlaw emission with a photon index of $Γ= 1.78$. The observed emission and published redshift imply an observed 2-10 keV luminosity of $3.4\times10^{45}$ erg s$^{-1}$. This is consistent with extrapolations from a previous soft X-ray detection, thus there is no evidence for strong variability. △ Less

Submitted 26 September, 2018; originally announced September 2018.

Comments: published in Research Notes of the American Astronomical Society

Journal ref: G. C. Privon et al. 2018 Res. Notes AAS 2 177

Disclosing the properties of low-redshift dual AGN through XMM-Newton and SDSS spectroscopy

Authors: A. De Rosa, C. Vignali, B. Husemann, S. Bianchi, T. Bogdanovic, M. Guainazzi, R. Herrero-Illana, S. Komossa, E. Kun, N. Loiseau, Z. Paragi, M. Perez-Torres, E. Piconcelli

Abstract: We report on an optical (SDSS) and X-ray (XMM) study of an optically selected sample of four dual AGN systems at projected separations of 30--60~kpc. All sources are detected in the X-ray band (0.3-10 keV); seven objects are optically identified as Seyfert, while one source, optically classified as a LINER, is likely powered by accretion in virtue of its relatively high X-ray luminosity (1.2… ▽ More We report on an optical (SDSS) and X-ray (XMM) study of an optically selected sample of four dual AGN systems at projected separations of 30--60~kpc. All sources are detected in the X-ray band (0.3-10 keV); seven objects are optically identified as Seyfert, while one source, optically classified as a LINER, is likely powered by accretion in virtue of its relatively high X-ray luminosity (1.2$\times10^{41}$ erg/s). Six of the eight objects are obscured in X-rays with N$_{\rm H} \geq$ 10$^{23}$ cm$^{-2}$; three of these, whose X-ray spectrum is dominated by a reflection component, are likely Compton-thick (N$_{\rm H} \geq$ 10$^{24}$ cm$^{-2}$). This finding is in agreement with the hypothesis that galaxy encounters are effective in driving gas inflow toward the nuclear region, thus increasing the obscuration. We compare the absorption properties in our dual AGN with those in larger samples observed in X-rays but selected in different ways (optical, IR and hard X-rays). We find that the obscured (N$_{\rm H} \geq$ 10$^{22}$ cm$^{-2}$) AGN fraction within the larger sample is 84$\pm$4 per cent (taking into account the 90 per cent error on the N$_{\rm H}$ measure) up to large pair separations ($\sim$100~kpc). This is statistically higher than the fraction of obscured AGN in isolated galaxies found in X-ray surveys. △ Less

Submitted 9 July, 2018; originally announced July 2018.

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

A dust-enshrouded tidal disruption event with a resolved radio jet in a galaxy merger

Authors: S. Mattila, M. Pérez-Torres, A. Efstathiou, P. Mimica, M. Fraser, E. Kankare, A. Alberdi, M. Á. Aloy, T. Heikkilä, P. G. Jonker, P. Lundqvist, I. Martí-Vidal, W. P. S. Meikle, C. Romero-Cañizales, S. J. Smartt, S. Tsygankov, E. Varenius, A. Alonso-Herrero, M. Bondi, C. Fransson, R. Herrero-Illana, T. Kangas, R. Kotak, N. Ramírez-Olivencia, P. Väisänen , et al. (11 additional authors not shown)

Abstract: Tidal disruption events (TDEs) are transient flares produced when a star is ripped apart by the gravitational field of a supermassive black hole (SMBH). We have observed a transient source in the western nucleus of the merging galaxy pair Arp 299 that radiated >1.5x10^52 erg in the infrared and radio, but was not luminous at optical or X-ray wavelengths. We interpret this as a TDE with much of its… ▽ More Tidal disruption events (TDEs) are transient flares produced when a star is ripped apart by the gravitational field of a supermassive black hole (SMBH). We have observed a transient source in the western nucleus of the merging galaxy pair Arp 299 that radiated >1.5x10^52 erg in the infrared and radio, but was not luminous at optical or X-ray wavelengths. We interpret this as a TDE with much of its emission re-radiated at infrared wavelengths by dust. Efficient reprocessing by dense gas and dust may explain the difference between theoretical predictions and observed luminosities of TDEs. The radio observations resolve an expanding and decelerating jet, probing the jet formation and evolution around a SMBH. △ Less

Submitted 14 June, 2018; originally announced June 2018.

Comments: Published online by Science on 14 Jun 2018

Variability Timescale and Spectral Index of Sgr A* in the Near Infrared: Approximate Bayesian Computation Analysis of the Variability of the Closest Supermassive Black Hole

Authors: G. Witzel, G. Martinez, J. Hora, S. P. Willner, M. R. Morris, C. Gammie, E. E. Becklin, M. L. N. Ashby, F. Baganoff, S. Carey, T. Do, G. G. Fazio, A. Ghez, W. J. Glaccum, D. Haggard, R. Herrero-Illana, J. Ingalls, R. Narayan, H. A. Smith

Abstract: Sagittarius A* (Sgr A*) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We present an analysis of the most comprehensive NIR variability dataset of Sgr A* to date: eight 24-hour epochs of continuous monitoring of Sgr A* at 4.5 $μ$m with the IRAC instrument on the Spitzer Space Telescope, 93 epochs of 2.18 $μ$m data from Na… ▽ More Sagittarius A* (Sgr A*) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We present an analysis of the most comprehensive NIR variability dataset of Sgr A* to date: eight 24-hour epochs of continuous monitoring of Sgr A* at 4.5 $μ$m with the IRAC instrument on the Spitzer Space Telescope, 93 epochs of 2.18 $μ$m data from Naos Conica at the Very Large Telescope, and 30 epochs of 2.12 $μ$m data from the NIRC2 camera at the Keck Observatory, in total 94,929 measurements. A new approximate Bayesian computation method for fitting the first-order structure function extracts information beyond current Fast Fourier Transformation (FFT) methods of power spectral density (PSD) estimation. With a combined fit of the data of all three observatories, the characteristic coherence timescale of Sgr A* is $τ_{b} = 243^{+82}_{-57}$ minutes ($90\%$ credible interval). The PSD has no detectable features on timescales down to 8.5 minutes ($95\%$ credible level), which is the ISCO orbital frequency for a dimensionless spin parameter $a = 0.92$. One light curve measured simultaneously at 2.12 and 4.5 $μ$m during a low flux-density phase gave a spectral index $α_s = 1.6 \pm 0.1$ ($F_ν\propto ν^{-α_s}$). This value implies that the Sgr A* NIR color becomes bluer during higher flux-density phases. The probability densities of flux densities of the combined datasets are best fit by log-normal distributions. Based on these distributions, the Sgr A* spectral energy distribution is consistent with synchrotron radiation from a non-thermal electron population from below 20 GHz through the NIR. △ Less

Submitted 5 September, 2018; v1 submitted 1 June, 2018; originally announced June 2018.

Comments: Accepted for publication in ApJ on May 30, 2018. A machine readable version of the light curve data is included in the journal's online publication. Version 2 includes proof corrections

Sub-arcsecond imaging of Arp\,299-A at 150 MHz with LOFAR: Evidence for a starburst-driven outflow

Authors: Naím Ramírez-Olivencia, Eskil Varenius, Miguel Pérez-Torres, Antonio Alberdi, Enrique Pérez, Almudena Alonso-Herrero, Adam Deller, Rubén Herrero-Illana, Javier Moldón, Loreto Barcos-Muñoz, Iván Martí-Vidal

Abstract: We report on the first sub-arcsecond (0.44 $\times$ 0.41 arcsec$\rm ^2$) angular resolution image at 150 MHz of the A-nucleus in the Luminous Infrared Galaxy Arp$\,$299, from International Low Frequency Array (LOFAR) Telescope observations. The most remarkable finding is that of an intriguing two-sided, filamentary structure emanating from A-nucleus, which we interpret as an outflow that extends u… ▽ More We report on the first sub-arcsecond (0.44 $\times$ 0.41 arcsec$\rm ^2$) angular resolution image at 150 MHz of the A-nucleus in the Luminous Infrared Galaxy Arp$\,$299, from International Low Frequency Array (LOFAR) Telescope observations. The most remarkable finding is that of an intriguing two-sided, filamentary structure emanating from A-nucleus, which we interpret as an outflow that extends up to at least 14 arcseconds from the A-nucleus in the N-S direction ($\approx$ 5 kpc deprojected size) and accounts for almost 40% of the extended emission of the entire galaxy system. We also discuss HST/NICMOS [FeII] 1.64 $\rm μm$ and H$\rm_2$ 2.12 $\rm μm$ images of Arp$\,$299-A, which show similar features to those unveiled by our 150 MHz LOFAR observations, thus giving string morphological support for the outflow scenario. Finally, we discuss unpublished NaI D spectra that confirm the outflow nature of this structure. From energetic arguments, we rule out the low-luminosity active galactic nucleus in Arp$\,$299-A as a driver for the outflow. On the contrary, the powerful, compact starburst in the central regions of Arp$\,$299-A provides plenty of mechanical energy to sustain an outflow, and we conclude that the intense supernova (SN) activity in the nuclear region of Arp299-A is driving the observed outflow. We estimate that the starburst wind can support a mass-outflow rate in the range (11-63) $\rm M_{\odot} yr^{-1}$ at speeds of up to (370 - 890) $\rm km \, s^{-1}$, and is relatively young, with an estimated kinematic age of (3 - 7) Myr. Those results open an avenue to the use of low-frequency (150 MHz), sub-arcsecond imaging with LOFAR to detect outflows in the central regions of local luminous infrared galaxies. △ Less

Submitted 20 February, 2018; v1 submitted 9 February, 2018; originally announced February 2018.

First results from GeMS/GSAOI for project SUNBIRD: Supernovae UNmasked By Infra-Red Detection

Authors: E. C. Kool, S. Ryder, E. Kankare, S. Mattila, T. Reynolds, R. M. McDermid, M. A. Pérez-Torres, R. Herrero-Illana, M. Schirmer, A. Efstathiou, F. E. Bauer, J. Kotilainen, P. Väisänen, C. Baldwin, C. Romero-Cañizales, A. Alberdi

Abstract: Core collapse supernova (CCSN) rates suffer from large uncertainties as many CCSNe exploding in regions of bright background emission and significant dust extinction remain unobserved. Such a shortfall is particularly prominent in luminous infrared galaxies (LIRGs), which have high star formation (and thus CCSN) rates and host bright and crowded nuclear regions, where large extinctions and reduced… ▽ More Core collapse supernova (CCSN) rates suffer from large uncertainties as many CCSNe exploding in regions of bright background emission and significant dust extinction remain unobserved. Such a shortfall is particularly prominent in luminous infrared galaxies (LIRGs), which have high star formation (and thus CCSN) rates and host bright and crowded nuclear regions, where large extinctions and reduced search detection efficiency likely lead to a significant fraction of CCSNe remaining undiscovered. We present the first results of project SUNBIRD (Supernovae UNmasked By Infra-Red Detection), where we aim to uncover CCSNe that otherwise would remain hidden in the complex nuclear regions of LIRGs, and in this way improve the constraints on the fraction that is missed by optical seeing-limited surveys. We observe in the near-infrared 2.15 μm $K_s$-band, which is less affected by dust extinction compared to the optical, using the multi-conjugate adaptive optics imager GeMS/GSAOI on Gemini South, allowing us to achieve a spatial resolution that lets us probe close in to the nuclear regions. During our pilot program and subsequent first full year we have discovered three CCSNe and one candidate with projected nuclear offsets as small as 200 pc. When compared to the total sample of LIRG CCSNe discovered in the near-IR and optical, we show that our method is singularly effective in uncovering CCSNe in nuclear regions and we conclude that the majority of CCSNe exploding in LIRGs are not detected as a result of dust obscuration and poor spatial resolution. △ Less

Submitted 24 September, 2017; originally announced September 2017.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society, 19 pages, 11 figures

No AGN evidence in NGC1614 from deep radio VLBI observations

Authors: Rubén Herrero-Illana, Antxon Alberdi, Miguel Ángel Pérez-Torres, Almudena Alonso-Herrero, Daniel González-Millán, Miguel Pereira-Santaella

Abstract: We present deep dual-band 5.0 and 8.4GHz European VLBI Network (EVN) observations of NGC1614, a local luminous infrared galaxy with a powerful circumnuclear starburst ring, and whose nuclear engine origin is still controversial. We aim at detecting and characterizing compact radio structures both in the nuclear region and in the circumnuclear ring. We do not find any compact source in the central… ▽ More We present deep dual-band 5.0 and 8.4GHz European VLBI Network (EVN) observations of NGC1614, a local luminous infrared galaxy with a powerful circumnuclear starburst ring, and whose nuclear engine origin is still controversial. We aim at detecting and characterizing compact radio structures both in the nuclear region and in the circumnuclear ring. We do not find any compact source in the central 200pc region, setting a very tight 5 sigma upper limit of $3.7\times10^{36}$erg s$^{-1}$ and $5.8\times10^{36}$erg s$^{-1}$, at 5.0 and 8.4GHz, respectively. However, we report a clear detection at both frequencies of a compact structure in the circumnuclear ring, 190pc to the north of the nucleus, whose luminosity and spectral index are compatible with a core-collapse supernova, giving support to the high star formation rate in the ring. Our result favors the pure starburst scenario, even for the nucleus of NGC1614, and shows the importance of radio VLBI observations when dealing with the obscured environments of dusty galaxies. △ Less

Submitted 30 May, 2017; originally announced May 2017.

Comments: Letter accepted for publication in MNRAS. 5 pages, 2 figures

Star formation and AGN activity in a sample of local Luminous Infrared Galaxies through multi-wavelength characterization

Authors: Rubén Herrero-Illana, Miguel Á. Pérez-Torres, Zara Randriamanakoto, Antxon Alberdi, Andreas Efstathiou, Petri Väisänen, Erkki Kankare, Erik Kool, Seppo Mattila, Rajin Ramphul, Stuart Ryder

Abstract: Nuclear starbursts and AGN activity are the main heating processes in luminous infrared galaxies (LIRGs) and their relationship is fundamental to understand galaxy evolution. In this paper, we study the star-formation and AGN activity of a sample of 11 local LIRGs imaged with subarcsecond angular resolution at radio (8.4GHz) and near-infrared ($2.2μ$m) wavelengths. This allows us to characterize t… ▽ More Nuclear starbursts and AGN activity are the main heating processes in luminous infrared galaxies (LIRGs) and their relationship is fundamental to understand galaxy evolution. In this paper, we study the star-formation and AGN activity of a sample of 11 local LIRGs imaged with subarcsecond angular resolution at radio (8.4GHz) and near-infrared ($2.2μ$m) wavelengths. This allows us to characterize the central kpc of these galaxies with a spatial resolution of $\simeq100$pc. In general, we find a good spatial correlation between the radio and the near-IR emission, although radio emission tends to be more concentrated in the nuclear regions. Additionally, we use an MCMC code to model their multi-wavelength spectral energy distribution (SED) using template libraries of starburst, AGN and spheroidal/cirrus models, determining the luminosity contribution of each component, and finding that all sources in our sample are starburst-dominated, except for NGC6926 with an AGN contribution of $\simeq64$\%. Our sources show high star formation rates ($40$ to $167M_\odot\mathrm{yr}^{-1}$), supernova rates (0.4 to $2.0\mathrm{SN}\mathrm{yr}^{-1}$), and similar starburst ages (13 to $29\mathrm{Myr}$), except for the young starburst (9Myr) in NGC6926. A comparison of our derived star-forming parameters with estimates obtained from different IR and radio tracers shows an overall consistency among the different star formation tracers. AGN tracers based on mid-IR, high-ionization line ratios also show an overall agreement with our SED model fit estimates for the AGN. Finally, we use our wide-band VLA observations to determine pixel-by-pixel radio spectral indices for all galaxies in our sample, finding a typical median value ($α\simeq-0.8$) for synchrotron-powered LIRGs. △ Less

Submitted 3 July, 2017; v1 submitted 26 May, 2017; originally announced May 2017.

Comments: Accepted for publication in MNRAS. 20 pages, 12 figures

Constraining magnetic field amplification in SN shocks using radio observations of SNe 2011fe and 2014J

Authors: Esha Kundu, Peter Lundqvist, Miguel A. Pérez-Torres, Rubén Herrero-Illana, Antxon Alberdi

Abstract: We modeled the radio non-detection of two Type Ia supernovae (SNe) 2011fe and 2014J considering synchrotron emission from the interaction between SN ejecta and the circumstellar medium. For an ejecta with the outer part having a power law density structure we compare synchrotron emission with radio observations. Assuming that 20$\%$ of the bulk shock energy is being shared equally between electron… ▽ More We modeled the radio non-detection of two Type Ia supernovae (SNe) 2011fe and 2014J considering synchrotron emission from the interaction between SN ejecta and the circumstellar medium. For an ejecta with the outer part having a power law density structure we compare synchrotron emission with radio observations. Assuming that 20$\%$ of the bulk shock energy is being shared equally between electrons and magnetic fields we found a very low density medium around both the SNe. A less tenuous medium with particle density $\sim$ 1 $\rm cm^{-3}$, which could be expected around both SNe, can be estimated when the magnetic field amplification is less than that presumed for energy equipartition. This conclusion also holds if the progenitor of SN 2014J was a rigidly rotating white dwarf (WD) with a main sequence (MS) or red giant companion. For a He star companion, or a MS for SN 2014J, with 10$\%$ and 1$\%$ of bulk kinetic energy in magnetic fields, we obtain a mass loss rate $< 10^{-9}$ and $< \sim 4\times 10^{-9}$ M$_{\odot}$yr$^{-1}$ for a wind velocity of 100 km/s. The former requires a mass accretion efficiency $>$ 99$\%$ onto the WD, but is less restricted for the latter case. However, if the tenuous medium is due to a recurrent nova it is difficult from our model to predict synchrotron luminosities. Although the formation channels of SNe 2011fe and 2014J are not clear, the null detection in radio wavelengths could point toward a low amplification efficiency for magnetic fields in SN shocks. △ Less

Submitted 11 May, 2017; originally announced May 2017.

Comments: 10 pages, 4 figures; Accepted for publication in ApJ

Journal ref: ApJ 842, 17 (2017)

Unveiling the AGN in IC 883: discovery of a parsec-scale radio jet

Authors: Cristina Romero-Cañizales, Antxon Alberdi, Claudio Ricci, Patricia Arévalo, Miguel A. Pérez-Torres, John E. Conway, Rob J. Beswick, Marco Bondi, Tom W. B. Muxlow, Megan K. Argo, Franz E. Bauer, Andreas Efstathiou, Rubén Herrero-Illana, Seppo Mattila, Stuart D. Ryder

Abstract: IC883 is a luminous infrared galaxy (LIRG) classified as a starburst-active galactic nucleus (AGN) composite. In a previous study we detected a low-luminosity AGN (LLAGN) radio candidate. Here we report on our radio follow-up at three frequencies which provides direct and unequivocal evidence of the AGN activity in IC883. Our analysis of archival X-ray data, together with the detection of a transi… ▽ More IC883 is a luminous infrared galaxy (LIRG) classified as a starburst-active galactic nucleus (AGN) composite. In a previous study we detected a low-luminosity AGN (LLAGN) radio candidate. Here we report on our radio follow-up at three frequencies which provides direct and unequivocal evidence of the AGN activity in IC883. Our analysis of archival X-ray data, together with the detection of a transient radio source with luminosity typical of bright supernovae, give further evidence of the ongoing star formation activity, which dominates the energetics of the system. At sub-parsec scales, the radio nucleus has a core-jet morphology with the jet being a newly ejected component showing a subluminal proper motion of 0.6c-1c. The AGN contributes less than two per cent of the total IR luminosity of the system. The corresponding Eddington factor is ~1E-3, suggesting this is a low-accretion rate engine, as often found in LLAGNs. However, its high bolometric luminosity (~10E44erg/s) agrees better with a normal AGN. This apparent discrepancy may just be an indication of the transition nature of the nucleus from a system dominated by star-formation, to an AGN-dominated system. The nucleus has a strongly inverted spectrum and a turnover at ~4.4GHz, thus qualifying as a candidate for the least luminous (L_5.0GHz ~ 6.3E28erg/s/Hz) and one of the youngest (~3000yr) gigahertz-peaked spectrum (GPS) sources. If the GPS origin for the IC883 nucleus is confirmed, then advanced mergers in the LIRG category are potentially key environments to unveil the evolution of GPS sources into more powerful radio galaxies. △ Less

Submitted 24 January, 2017; originally announced January 2017.

Comments: 13 pages, 8 figures. Accepted for publication in MNRAS

Cosmic evolution and metal aversion in super-luminous supernova host galaxies

Authors: S. Schulze, T. Krühler, G. Leloudas, J. Gorosabel, A. Mehner, J. Buchner, S. Kim, E. Ibar, R. Amorín, R. Herrero-Illana, J. P. Anderson, F. E. Bauer, L. Christensen, M. de Pasquale, A. de Ugarte Postigo, A. Gallazzi, J. Hjorth, N. Morrell, D. Malesani, M. Sparre, B. Stalder, A. A. Stark, C. C. Thöne, J. C. Wheeler

Abstract: The SUperluminous Supernova Host galaxIES (SUSHIES) survey aims to provide strong new constraints on the progenitors of superluminous supernovae (SLSNe) by understanding the relationship to their host galaxies. We present the photometric properties of 53 H-poor and 16 H-rich SLSN host galaxies out to $z\sim4$. We model their spectral energy distributions to derive physical properties, which we com… ▽ More The SUperluminous Supernova Host galaxIES (SUSHIES) survey aims to provide strong new constraints on the progenitors of superluminous supernovae (SLSNe) by understanding the relationship to their host galaxies. We present the photometric properties of 53 H-poor and 16 H-rich SLSN host galaxies out to $z\sim4$. We model their spectral energy distributions to derive physical properties, which we compare with other galaxy populations. At low redshift, H-poor SLSNe are preferentially found in very blue, low-mass galaxies with high average specific star-formation rates. As redshift increases, the host population follows the general evolution of star-forming galaxies towards more luminous galaxies. After accounting for secular evolution, we find evidence for differential evolution in galaxy mass, but not in the $B$-band and the far UV luminosity ($3σ$ confidence). Most remarkable is the scarcity of hosts with stellar masses above $10^{10}~M_\odot$ for both classes of SLSNe. In the case of H-poor SLSNe, we attribute this to a stifled production efficiency above $\sim0.4$ solar metallicity. However, we argue that, in addition to low metallicity, a short-lived stellar population is also required to regulate the SLSN production. H-rich SLSNe are found in a very diverse population of star-forming galaxies. Still, the scarcity of massive hosts suggests a stifled production efficiency above $\sim0.8$ solar metallicity. The large dispersion of the H-rich SLSNe host properties is in stark contrast to those of gamma-ray burst, regular core-collapse SN, and H-poor SLSNe host galaxies. We propose that multiple progenitor channels give rise to this sub-class. △ Less

Submitted 3 November, 2017; v1 submitted 18 December, 2016; originally announced December 2016.

Comments: Accepted on 8 September 2017. Currently in press at MNRAS. Substantially revised version. 61 pages, 31 pages main text, 30 pages appendices, 22 figures and 9 tables

Unveiling multiple AGN activity in galaxy mergers

Authors: A. De Rosa, S. Bianchi, T. Bogdanovic, R. Decarli, J. Heidt, R. Herrero-Illana, B. Husemann, S. Komossa, E. Kun, N. Loiseau, M. Guainazzi, Z. Paragi, M. Perez-Torres, E. Piconcelli, K. Schawinski, C. Vignali

Abstract: In this paper we present an overview of the MAGNA (Multiple AGN Activity) project aiming at a comprehensive study of multiple supemassive black hole systems. With the main goal to characterize the sources in merging systems at different stages of evolution, we selected a sample of objects optically classified as multiple systems on the basis of emission line diagnostics and started a massive multi… ▽ More In this paper we present an overview of the MAGNA (Multiple AGN Activity) project aiming at a comprehensive study of multiple supemassive black hole systems. With the main goal to characterize the sources in merging systems at different stages of evolution, we selected a sample of objects optically classified as multiple systems on the basis of emission line diagnostics and started a massive multiband observational campaign. Here we report on the discovery of the exceptionally high AGN density compact group SDSS~J0959+1259. A multiband study suggests that strong interactions are taking place among its galaxies through tidal forces, therefore this system represents a case study for physical mechanisms that trigger nuclear activity and star formation. We also present a preliminary analysis of the multiple AGN system SDSS~J1038+3921.} △ Less

Submitted 2 November, 2016; originally announced November 2016.

Comments: Proceedings paper for "XMM-Newton: The Next Decade" (held in 2016 May); accepted to be published in Astronomische Nachrichten / Astronomical Notes (AN)

Radio continuum and X-ray emission from the most extreme FIR-excess galaxy NGC 1377: An extremely obscured AGN revealed

Authors: F. Costagliola, R. Herrero-Illana, A. Lohfink, M. Pérez-Torres, S. Aalto, S. Muller, A. Alberdi

Abstract: Galaxies which strongly deviate from the radio-far IR correlation are of great importance for studies of galaxy evolution as they may be tracing early, short-lived stages of starbursts and active galactic nuclei (AGNs). The most extreme FIR-excess galaxy NGC1377 has long been interpreted as a young dusty starburst, but millimeter observations of CO lines revealed a powerful collimated molecular ou… ▽ More Galaxies which strongly deviate from the radio-far IR correlation are of great importance for studies of galaxy evolution as they may be tracing early, short-lived stages of starbursts and active galactic nuclei (AGNs). The most extreme FIR-excess galaxy NGC1377 has long been interpreted as a young dusty starburst, but millimeter observations of CO lines revealed a powerful collimated molecular outflow which cannot be explained by star formation alone. We present new radio observations at 1.5 and 10 GHz obtained with the Jansky Very Large Array (JVLA) and Chandra X-ray observations towards NGC1377. The observations are compared to synthetic starburst models to constrain the properties of the central energy source. We obtained the first detection of the cm radio continuum and X-ray emission in NGC1377. We find that the radio emission is distributed in two components, one on the nucleus and another offset by 4$"$.5 to the South-West. We confirm the extreme FIR-excess of the galaxy, with a $q_\mathrm{FIR}\simeq$4.2, which deviates by more than 7-$σ$ from the radio-FIR correlation. Soft X-ray emission is detected on the off-nucleus component. From the radio emission we estimate for a young ($<10$ Myr) starburst a star formation rate SFR$<$0.1 M$_\odot$ yr$^{-1}$. Such a SFR is not sufficient to power the observed IR luminosity and to drive the CO outflow. We find that a young starburst cannot reproduce all the observed properties of the nucleus of NGC1377. We suggest that the galaxy may be harboring a radio-quiet, obscured AGN of 10$^6$M$_\odot$, accreting at near-Eddington rates. We speculate that the off-nucleus component may be tracing an hot-spot in the AGN jet. △ Less

Submitted 14 July, 2016; originally announced July 2016.

Comments: 14 pages, accepted for publication on Astronomy and Astrophysics on 08/07/2016

Journal ref: A&A 594, A114 (2016)

Probing dust-obscured star formation in the most massive Gamma-Ray Burst host galaxies

Authors: Jochen Greiner, Michal J. Michalowski, Sylvio Klose, Leslie K. Hunt, Gianfranco Gentile, Peter Kamphuis, Ruben Herrero-Illana, Mark Wieringa, Thomas Krühler, Patricia Schady, Jonathan Elliott, John F. Graham, Eduardo Ibar, Fabian Knust, Ana Nicuesa Guelbenzu, Eliana Palazzi, Andrea Rossi, Sandra Savaglio

Abstract: Due to their relation to massive stars, long-duration gamma-ray bursts (GRBs) allow pinpointing star formation in galaxies independently of redshift, dust obscuration, or galaxy mass/size, thus providing a unique tool to investigate the star-formation history over cosmic time. About half of the optical afterglows of long-duration GRBs are missed due to dust extinction, and are primarily located in… ▽ More Due to their relation to massive stars, long-duration gamma-ray bursts (GRBs) allow pinpointing star formation in galaxies independently of redshift, dust obscuration, or galaxy mass/size, thus providing a unique tool to investigate the star-formation history over cosmic time. About half of the optical afterglows of long-duration GRBs are missed due to dust extinction, and are primarily located in the most massive GRB hosts. In order to understand this bias it is important to investigate the amount of obscured star-formation in these GRB host galaxies. Radio emission of galaxies correlates with star-formation, but does not suffer extinction as do the optical star-formation estimators. We selected 11 GRB host galaxies with either large stellar mass or large UV-/optical-based star-formation rates (SFRs) and obtained radio observations of these with the Australia Telescope Compact Array and the Karl Jansky Very Large Array. Despite intentionally selecting GRB hosts with expected high SFRs, we do not find any star-formation-related radio emission in any of our targets. Our upper limit for GRB 100621A implies that the earlier reported radio detection was due to afterglow emission. We do detect radio emission from the position of GRB 020819B, but argue that it is in large parts, if not all, due to afterglow contamination. Half of our sample has radio-derived SFR limits which are only a factor 2--3 above the optically measured SFRs. This supports other recent studies that the majority of star formation in GRB hosts is not obscured by dust. △ Less

Submitted 27 June, 2016; originally announced June 2016.

Comments: 12 pages, 5 figs., A&A accepted

Journal ref: A&A 593, A17 (2016)

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

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

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

Submitted 25 September, 2015; originally announced September 2015.

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

Excitation Mechanisms for HCN (1-0) and HCO+ (1-0) in Galaxies from the Great Observatories All-sky LIRG Survey

Authors: G. C. Privon, R. Herrero-Illana, A. S. Evans, K. Iwasawa, M. A. Perez-Torres, L. Armus, T. Diaz-Santos, E. J. Murphy, S. Stierwalt, S. Aalto, J. M. Mazzarella, L. Barcos-Munoz, H. J. Borish, H. Inami, D. -C. Kim, E. Treister, J. A. Surace, S. Lord, J. Conway, D. T. Frayer, A. Alberdi

Abstract: We present new IRAM 30m spectroscopic observations of the $\sim88$ GHz band, including emission from the CCH (n=1-0) multiplet, HCN (1-0), HCO+ (1-0), and HNC (1-0), for a sample of 58 local luminous and ultraluminous infrared galaxies from the Great Observatories All-sky LIRG Survey (GOALS). By combining our new IRAM data with literature data and Spitzer/IRS spectroscopy, we study the corresponde… ▽ More We present new IRAM 30m spectroscopic observations of the $\sim88$ GHz band, including emission from the CCH (n=1-0) multiplet, HCN (1-0), HCO+ (1-0), and HNC (1-0), for a sample of 58 local luminous and ultraluminous infrared galaxies from the Great Observatories All-sky LIRG Survey (GOALS). By combining our new IRAM data with literature data and Spitzer/IRS spectroscopy, we study the correspondence between these putative tracers of dense gas and the relative contribution of active galactic nuclei (AGN) and star formation to the mid-infrared luminosity of each system. We find the HCN (1-0) emission to be enhanced in AGN-dominated systems ($\langle$L'$_{HCN (1-0)}$/L'$_{HCO^+ (1-0)}\rangle=1.84$), compared to composite and starburst-dominated systems ($\langle$L'$_{HCN (1-0)}$/L'$_{HCO^+ (1-0)}\rangle=1.14$, and 0.88, respectively). However, some composite and starburst systems have L'$_{HCN (1-0)}$/L'$_{HCO^+ (1-0)}$ ratios comparable to those of AGN, indicating that enhanced HCN emission is not uniquely associated with energetically dominant AGN. After removing AGN-dominated systems from the sample, we find a linear relationship (within the uncertainties) between $\log_{10}$(L'$_{HCN (1-0)}$) and $\log_{10}$(L$_{IR}$), consistent with most previous findings. L'$_{HCN (1-0)}$/L$_{IR}$, typically interpreted as the dense gas depletion time, appears to have no systematic trend with L$_{IR}$ for our sample of luminous and ultraluminous infrared galaxies, and has significant scatter. The galaxy-integrated HCN (1-0) and HCO+ (1-0) emission do not appear to have a simple interpretation, in terms of the AGN dominance or the star formation rate, and are likely determined by multiple processes, including density and radiative effects. △ Less

Submitted 8 December, 2016; v1 submitted 24 September, 2015; originally announced September 2015.

Comments: fix presentation of numerical value of y-intercept in eqs 1-3 to match forthcoming ApJ erratum. figure 6 and the conclusions are unaffected

Journal ref: ApJ, 814, 39 (2015)

Multiple AGN in the crowded field of the compact group SDSSJ0959+1259

Authors: Alessandra De Rosa, Stefano Bianchi, Tamara Bogdanovic, Roberto Decarli, Ruben Herrero-Illana, Bernd Husemann, Stefanie Komossa, Emma Kun, Nora Loiseau, Zsolt Paragi, Miguel Perez-Torres, Enrico Piconcelli, Kevin Schawinski, Cristian Vignali

Abstract: We present a multi-wavelength study of a newly discovered compact group (CG), SDSS J0959+1259, based data from XMM-Newton, SDSS and the Calar Alto optical imager BUSCA. With a maximum velocity offset of 500 km s$^{-1}$, a mean redshift of 0.035, and a mean spatial extension of 480 kpc, this CG is exceptional in having the highest concentration of nuclear activity in the local Universe, established… ▽ More We present a multi-wavelength study of a newly discovered compact group (CG), SDSS J0959+1259, based data from XMM-Newton, SDSS and the Calar Alto optical imager BUSCA. With a maximum velocity offset of 500 km s$^{-1}$, a mean redshift of 0.035, and a mean spatial extension of 480 kpc, this CG is exceptional in having the highest concentration of nuclear activity in the local Universe, established with a sensitivity limit L$_{X}>4\times $10$^{40}$ erg s$^{-1}$ in 2--10 keV band and R-band magnitude $M_R < -19$. The group is composed of two type-2 Seyferts, one type-1 Seyfert, two LINERs and three star forming galaxies. Given the high X-ray luminosity of LINERs which reaches $\sim 10^{41}$ erg s$^{-1}$, it is likely that they are also accretion driven, bringing the number of active nuclei in this group to to 5 out of 8 (AGN fraction of 60\%). The distorted shape of one member of the CG suggests that strong interactions are taking place among its galaxies through tidal forces. Therefore, this system represents a case study for physical mechanisms that trigger nuclear activity and star formation in CGs. △ Less

Submitted 19 July, 2015; v1 submitted 16 July, 2015; originally announced July 2015.

Comments: 8 pages, 4 figures and 4 tables. Accepted for publication in MNRAS

The Spanish Square Kilometre Array White Book

Authors: J. A. Acosta-Pulido, I. Agudo, A. Alberdi, J. Alcolea, E. J. Alfaro, A. Alonso-Herrero, G. Anglada, P. Arnalte-Mur, Y. Ascasibar, B. Ascaso, R. Azulay, R. Bachiller, A. Baez-Rubio, E. Battaner, J. Blasco, C. B. Brook, V. Bujarrabal, G. Busquet, M. D. Caballero-Garcia, C. Carrasco-Gonzalez, J. Casares, A. J. Castro-Tirado, L. Colina, F. Colomer, I. de Gregorio-Monsalvo , et al. (94 additional authors not shown)

Abstract: The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on… ▽ More The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project. △ Less

Submitted 17 June, 2015; v1 submitted 10 June, 2015; originally announced June 2015.

Comments: Editors: Miguel A. Perez-Torres (Editor-in-chief), L. Verdes-Montenegro, J. C. Guirado, A. Alberdi, J. Martin-Pintado, R. Bachiller, D. Herranz, J. M. Girart, J. Gorgas, C. Hernandez-Monteagudo, S. Migliari and J. M. Rodriguez Espinosa; 289 pages, published in June 2015 by Sociedad Española de Astronomía

Diversity in extinction laws of Type Ia supernovae measured between $0.2$ and $2\,μ\mathrm{m}$

Authors: R. Amanullah, J. Johansson, A. Goobar, R. Ferretti, S. Papadogiannakis, T. Petrushevska, P. J. Brown, Y. Cao, C. Contreras, H. Dahle, N. Elias-Rosa, J. P. U. Fynbo, J. Gorosabel, L. Guaita, L. Hangard, D. A. Howell, E. Y. Hsiao, E. Kankare, M. Kasliwal, G. Leloudas, P. Lundqvist, S. Mattila, P. Nugent, M. M. Phillips, A. Sandberg , et al. (14 additional authors not shown)

Abstract: We present ultraviolet (UV) observations of six nearby Type~Ia supernovae (SNe~Ia) obtained with the Hubble Space Telescope, three of which were also observed in the near-IR (NIR) with Wide-Field Camera~3. UV observations with the Swift satellite, as well as ground-based optical and near-infrared data provide complementary information. The combined data-set covers the wavelength range $0.2$--… ▽ More We present ultraviolet (UV) observations of six nearby Type~Ia supernovae (SNe~Ia) obtained with the Hubble Space Telescope, three of which were also observed in the near-IR (NIR) with Wide-Field Camera~3. UV observations with the Swift satellite, as well as ground-based optical and near-infrared data provide complementary information. The combined data-set covers the wavelength range $0.2$--$2~μ$m. By also including archival data of SN 2014J, we analyse a sample spanning observed colour excesses up to $E(B-V)=1.4~$mag. We study the wavelength dependent extinction of each individual SN and find a diversity of reddening laws when characterised by the total-to-selective extinction $R_V$. In particular, we note that for the two SNe with $E(B-V)\gtrsim1~$mag, for which the colour excess is dominated by dust extinction, we find $R_V=1.4\pm0.1$ and $R_V=2.8\pm0.1$. Adding UV photometry reduces the uncertainty of fitted $R_V$ by $\sim50\,$% allowing us to also measure $R_V$ of individual low-extinction objects which point to a similar diversity, currently not accounted for in the analyses when SNe~Ia are used for studying the expansion history of the universe. △ Less

Submitted 8 July, 2015; v1 submitted 8 April, 2015; originally announced April 2015.

Comments: 31 pages, 28 figures, accepted for publication in MNRAS

ALMA Observations of Warm Dense Gas in NGC 1614 --- Breaking of Star Formation Law in the Central kpc

Authors: C. K. Xu, C. Cao, N. Lu, Y. Gao, T. Diaz-Santos, R. Herrero-Illana, R. Meijerink, G. Privon, Y. -H. Zhao, A. S. Evans, S. König, J. M. Mazzarella, S. Aalto, P. Appleton, L. Armus, V. Charmandaris, J. Chu, S. Haan, H. Inami, E. J. Murphy, D. B. Sanders, B. Schulz, P. van der Werf

Abstract: We present ALMA Cycle-0 observations of the CO (6-5) line emission and of the 435um dust continuum emission in the central kpc of NGC 1614, a local luminous infrared galaxy (LIRG) at a distance of 67.8 Mpc (1 arcsec = 329 pc). The CO emission is well resolved by the ALMA beam (0".26 x 0".20) into a circum-nuclear ring, with an integrated flux of f_{CO(6-5)} = 898 (+-153) Jy km/s, which is 63(+-12)… ▽ More We present ALMA Cycle-0 observations of the CO (6-5) line emission and of the 435um dust continuum emission in the central kpc of NGC 1614, a local luminous infrared galaxy (LIRG) at a distance of 67.8 Mpc (1 arcsec = 329 pc). The CO emission is well resolved by the ALMA beam (0".26 x 0".20) into a circum-nuclear ring, with an integrated flux of f_{CO(6-5)} = 898 (+-153) Jy km/s, which is 63(+-12)% of the total CO(6-5) flux measured by Herschel. The molecular ring, located between 100pc < r < 350pc from the nucleus, looks clumpy and includes seven unresolved (or marginally resolved) knots with median velocity dispersion of 40 km/s. These knots are associated with strong star formation regions with Σ_{SFR} 100 M_\sun/yr/kpc^{2} and Σ_{Gas} 1.0E4 M_\sun/pc^{2}. The non-detections of the nucleus in both the CO (6-5) line emission and the 435um continuum rule out, with relatively high confidence, a Compton-thick AGN in NGC 1614. Comparisons with radio continuum emission show a strong deviation from an expected local correlation between Σ_{Gas} and Σ_{SFR}, indicating a breakdown of the Kennicutt-Schmidt law on the linear scale of 100 pc. △ Less

Submitted 4 November, 2014; originally announced November 2014.

Comments: 13 pages, 10 figures; accepted by ApJ

Core-collapse and Type Ia supernovae with the SKA

Authors: M. A. Pérez-Torres, A. Alberdi, R. J. Beswick, P. Lundqvist, R. Herrero-Illana, C. Romero-Cañizales, S. Ryder, M. della Valle, J. Conway, J. M. Marcaide, S. Mattila, T. Murphy, E. Ros

Abstract: Core-collapse SNe (CCSNe): Systematic searches of radio emission from CCSNe are still lacking, and only targeted searches of radio emission from just some of the optically discovered CCSNe in the local universe have been carried out. Optical searches miss a significant fraction of CCSNe due to dust obscuration; CCSN radio searches are thus more promising for yielding the complete, unobscured star-… ▽ More Core-collapse SNe (CCSNe): Systematic searches of radio emission from CCSNe are still lacking, and only targeted searches of radio emission from just some of the optically discovered CCSNe in the local universe have been carried out. Optical searches miss a significant fraction of CCSNe due to dust obscuration; CCSN radio searches are thus more promising for yielding the complete, unobscured star-formation rates in the local universe. The SKA yields the possibility to piggyback for free in this area of research by carrying out commensal, wide-field, blind transient survey observations. SKA1-SUR should be able to discover several hundreds of CCSNe in just one year, compared to about a dozen CCSNe that the VLASS would be able to detect in one year, at most. SKA, with an expected sensitivity ten times that of SKA1, is expected to detect CCSNe in the local Universe by the thousands. Therefore, commensal SKA observations could easily result in an essentially complete census of all CCSNe in the local universe, thus yielding an accurate determination of the volumetric CCSN rate. Type Ia SNe: We advocate for the use of the SKA to search for the putative prompt (~first few days after the explosion) radio emission of any nearby type Ia SN, via target-of-opportunity observations. The huge improvement in sensitivity of the SKA with respect to its predecessors will allow to unambiguously discern which progenitor scenario (single-degenerate vs. double-degenerate) applies to them. △ Less

Submitted 4 February, 2015; v1 submitted 5 September, 2014; originally announced September 2014.

Comments: Proceedings of Science: Advancing Astrophysics with the Square Kilometre Array. Accepted version, including comments by the referee

A multi-wavelength view of the central kiloparsec region in the Luminous Infrared Galaxy NGC1614

Authors: Rubén Herrero-Illana, Miguel Á. Pérez-Torres, Almudena Alonso-Herrero, Antxon Alberdi, Luis Colina, Andreas Efstathiou, Lorena Hernández-García, Daniel Miralles-Caballero, Petri Väisänen, Christopher C. Packham, Vinesh Rajpaul, Albert A. Zijlstra

Abstract: The Luminous Infrared Galaxy NGC1614 hosts a prominent circumnuclear ring of star formation. However, the nature of the dominant emitting mechanism in its central ~100 pc is still under debate. We present sub-arcsecond angular resolution radio, mid-infrared, Pa-alpha, optical, and X-ray observations of NGC1614, aimed at studying in detail both the circumnuclear ring and the nuclear region. The 8.4… ▽ More The Luminous Infrared Galaxy NGC1614 hosts a prominent circumnuclear ring of star formation. However, the nature of the dominant emitting mechanism in its central ~100 pc is still under debate. We present sub-arcsecond angular resolution radio, mid-infrared, Pa-alpha, optical, and X-ray observations of NGC1614, aimed at studying in detail both the circumnuclear ring and the nuclear region. The 8.4 GHz continuum emission traced by the Very Large Array (VLA) and the Gemini/T-ReCS 8.7 micron emission, as well as the Pa-alpha line emission, show remarkable morphological similarities within the star-forming ring, suggesting that the underlying emission mechanisms are tightly related. We used an HST/NICMOS Pa-alpha map of similar resolution to our radio maps to disentangle the thermal free-free and non-thermal synchrotron radio emission, from which we obtained the intrinsic synchrotron power-law for each individual region within the central kpc of NGC1614. The radio ring surrounds a relatively faint, steep-spectrum source at the very center of the galaxy, suggesting that the central source is not powered by an AGN, but rather by a compact (r < 90 pc) starburst. Chandra X-ray data also show that the central kpc region is dominated by starburst activity, without requiring the existence of an AGN. We also used publicly available infrared data to model-fit the spectral energy distribution of both the starburst ring and a putative AGN in NGC1614. In summary, we conclude that there is no need to invoke an AGN to explain the observed bolometric properties of the galaxy. △ Less

Submitted 27 March, 2014; originally announced March 2014.

Comments: 13 pages, 7 figures, 5 tables. Accepted for publication in ApJ

The nature of supernovae 2010O and 2010P in Arp 299 - II. Radio emission

Authors: C. Romero-Cañizales, R. Herrero-Illana, M. A. Pérez-Torres, A. Alberdi, E. Kankare, F. E. Bauer, S. D. Ryder, S. Mattila, J. E. Conway, R. J. Beswick, T. W. B. Muxlow

Abstract: We report radio observations of two stripped-envelope supernovae (SNe), 2010O and 2010P, which exploded within a few days of each other in the luminous infrared galaxy Arp 299. Whilst SN 2010O remains undetected at radio frequencies, SN 2010P was detected (with an astrometric accuracy better than 1 milli arcsec in position) in its optically thin phase in epochs ranging from ~1 to ~3yr after its ex… ▽ More We report radio observations of two stripped-envelope supernovae (SNe), 2010O and 2010P, which exploded within a few days of each other in the luminous infrared galaxy Arp 299. Whilst SN 2010O remains undetected at radio frequencies, SN 2010P was detected (with an astrometric accuracy better than 1 milli arcsec in position) in its optically thin phase in epochs ranging from ~1 to ~3yr after its explosion date, indicating a very slow radio evolution and a strong interaction of the SN ejecta with the circumstellar medium. Our late-time radio observations toward SN 2010P probe the dense circumstellar envelope of this SN, and imply a mass-loss rate (Msun/yr) to wind velocity (in units of 10 km/s) ratio of (3.0-5.1)E-05, with a 5 GHz peak luminosity of ~1.2E+27 erg/s/Hz on day ~464 after explosion. This is consistent with a Type IIb classification for SN 2010P, making it the most distant and most slowly evolving Type IIb radio SN detected to date. △ Less

Submitted 5 March, 2014; originally announced March 2014.

Comments: 14 pages, 8 tables and 7 figures. Accepted for publication in MNRAS

The nature of supernovae 2010O and 2010P in Arp 299 - I. Near-infrared and optical evolution

Authors: E. Kankare, S. Mattila, S. Ryder, M. Fraser, A. Pastorello, N. Elias-Rosa, C. Romero-Canizales, A. Alberdi, V. -P. Hentunen, R. Herrero-Illana, J. Kotilainen, M. -A. Perez-Torres, P. Vaisanen

Abstract: We present near-infrared and optical photometry, plus optical spectroscopy of two stripped-envelope supernovae (SNe) 2010O and 2010P that exploded in two different components of an interacting luminous infrared galaxy Arp 299 within only a few days of one another. SN 2010O is found to be photometrically and spectroscopically similar to many normal Type Ib SNe and our multiwavelength observations o… ▽ More We present near-infrared and optical photometry, plus optical spectroscopy of two stripped-envelope supernovae (SNe) 2010O and 2010P that exploded in two different components of an interacting luminous infrared galaxy Arp 299 within only a few days of one another. SN 2010O is found to be photometrically and spectroscopically similar to many normal Type Ib SNe and our multiwavelength observations of SN 2010P suggest it to be a Type IIb SN. No signs of clear hydrogen features or interaction with the circumstellar medium are evident in the optical spectrum of SN 2010P. We derive estimates for the host galaxy line-of-sight extinctions for both SNe, based on both light-curve and spectroscopic comparison finding consistent results. These methods are also found to provide much more robust estimates of the SN host galaxy reddening than the commonly used empirical relations between extinction and equivalent width of Na I D absorption features. The SN observations also suggest that different extinction laws are present in different components of Arp 299. For completeness, we study high-resolution pre-explosion images of Arp 299 and find both SNe to be close to, but not coincident with, extended sources that are likely massive clusters. A very simple model applied to the bolometric light curve of SN 2010O implies a rough estimate for the explosion parameters of $E_{\mathrm{k}} \approx 3 \times 10^{51}$ erg, $M_{\mathrm{ej}} \approx 2.9$ M$_{\odot}$ and $M_{\mathrm{Ni}} \approx 0.16$ M$_{\odot}$. △ Less

Submitted 30 April, 2014; v1 submitted 25 November, 2013; originally announced November 2013.

Comments: 16 pages, 11 figures. Minor changes in text from v1. Published in MNRAS

Journal ref: 2014, MNRAS, 440, 1052

Evidence of nuclear disks in starburst galaxies from their radial distribution of supernovae

Authors: R. Herrero-Illana, M. A. Pérez-Torres, A. Alberdi

Abstract: Galaxy-galaxy interactions are expected to be responsible for triggering massive star formation and possibly accretion onto a supermassive black hole, by providing large amounts of dense molecular gas down to the central kiloparsec region. Several scenarios to drive the gas further down to the central ~100 pc, have been proposed, including the formation of a nuclear disk around the black hole, whe… ▽ More Galaxy-galaxy interactions are expected to be responsible for triggering massive star formation and possibly accretion onto a supermassive black hole, by providing large amounts of dense molecular gas down to the central kiloparsec region. Several scenarios to drive the gas further down to the central ~100 pc, have been proposed, including the formation of a nuclear disk around the black hole, where massive stars would produce supernovae. Here, we probe the radial distribution of supernovae and supernova remnants in the nuclear regions of the starburst galaxies M82, Arp 299-A, and Arp 220, by using high-angular resolution (< 0."1) radio observations published in the literature (for M82 and Arp 220), or obtained by ourselves from the European VLBI Network (Arp 299-A). Our main goal was to characterize the nuclear starbursts in those galaxies and thus test scenarios that propose that nuclear disks of sizes ~100 pc form in the central regions of starburst galaxies. We obtained the radial distribution of supernovae (SNe) in the nuclear starbursts of M82, Arp 299-A, and Arp 220, and derived scale-length values for the putative nuclear disks powering the bursts in those central regions. The scale lengths for the (exponential) disks range from ~20-30 pc for Arp 299-A and Arp 220, up to ~140 pc for M82. The radial distribution of SNe for the nuclear disks in Arp 299-A and Arp 220 is also consistent with a power-law surface density profile of exponent gamma=1, as expected from detailed hydrodynamical simulations of nuclear disks. Our results support scenarios where a nuclear disk of size ~100 pc is formed in (U)LIRGs, and sustained by gas pressure, in which case the accretion onto the black hole could be lowered by supernova feedback. △ Less

Submitted 13 March, 2012; originally announced March 2012.

Comments: Accepted for publication in Astronomy and Astrophysics

The nuclear starburst in Arp 299-A: From the 5.0 GHz VLBI radio light-curves to its core-collapse supernova rate

Authors: M. Bondi, M. A. Perez-Torres, R. Herrero-Illana, A. Alberdi

Abstract: The nuclear region of the Luminous Infra-red Galaxy Arp 299-A hosts a recent ($\simeq 10$ Myr), intense burst of massive star formation which is expected to lead to numerous core-collapse supernovae (CCSNe). Previous VLBI observations, carried out with the EVN at 5.0 GHz and with the VLBA at 2.3 and 8.4 GHz, resulted in the detection of a large number of compact, bright, non-thermal sources in a r… ▽ More The nuclear region of the Luminous Infra-red Galaxy Arp 299-A hosts a recent ($\simeq 10$ Myr), intense burst of massive star formation which is expected to lead to numerous core-collapse supernovae (CCSNe). Previous VLBI observations, carried out with the EVN at 5.0 GHz and with the VLBA at 2.3 and 8.4 GHz, resulted in the detection of a large number of compact, bright, non-thermal sources in a region $\lsim$150 pc in size. We aim at establishing the nature of all non-thermal, compact components in Arp 299-A, as well as estimating its core-collapse supernova rate. We use multi-epoch European VLBI Network (EVN) observations taken at 5.0 GHz to image with milliarcsecond resolution the compact radio sources in the nuclear region of Arp 299-A. We also use one single-epoch 5.0 GHz Multi-Element Radio Linked Interferometer Network (MERLIN) observation to image the extended emission in which the compact radio sources --traced by our EVN observations-- are embedded. Twenty-six compact sources are detected, 8 of them are new objects not previously detected. The properties of all detected objects are consistent with them being a mixed population of CCSNe and SNRs. We find clear evidence for at least two new CCSNe, implying a lower limit to the CCSN rate of $ν_{\rm SN}\gsim$0.80 SN/yr indicating that the bulk of the current star formation in Arp 299-A is taking place in the innermost $\sim 150$ pc. Our MERLIN observations trace a region of diffuse, extended emission which is co-spatial to the region where all compact sources are found. From this diffuse, non-thermal radio emission we obtain an independent estimate for the core-collapse supernova rate, which is in the range $ν_{\rm SN}=0.40$ - 0.65 SN/yr, roughly in agreement with previous estimates and our direct estimate of the CCSN rate from the compact radio emission. △ Less

Submitted 16 January, 2012; originally announced January 2012.

Comments: 13 pages, 5 figures, accepted for publication on Astronomy & Astrophysics