Showing 1–50 of 50 results for author: Emig, K

JWST Observations of Starbursts: Cold Clouds and Plumes Launching in the M82 Outflow

Authors: Deanne B. Fisher, Alberto D. Bolatto, John Chisholm, Drummond Fielding, Rebecca C. Levy, Elizabeth Tarantino, Martha L. Boyer, Serena A. Cronin, Laura A. Lopez, J. D. Smith, Danielle A. Berg, Sebastian Lopez, Sylvain Veilleux, Paul P. van der Werf, Torsten Böker, Leindert A. Boogaard, Laura Lenkić, Simon C. O. Glover, Vicente Villanueva, Divakara Mayya, Thomas S. -Y. Lai, Daniel A. Dale, Kimberly L. Emig, Fabian Walter, Monica Relaño , et al. (6 additional authors not shown)

Abstract: In this paper we study the filamentary substructure of 3.3 $μ$m PAH emission from JWST/NIRCam observations in the base of the M82 star-burst driven wind. We identify plume-like substructure within the PAH emission with widths of $\sim$50 pc. Several of the plumes extend to the edge of the field-of-view, and thus are at least 200-300 pc in length. In this region of the outflow, the vast majority (… ▽ More In this paper we study the filamentary substructure of 3.3 $μ$m PAH emission from JWST/NIRCam observations in the base of the M82 star-burst driven wind. We identify plume-like substructure within the PAH emission with widths of $\sim$50 pc. Several of the plumes extend to the edge of the field-of-view, and thus are at least 200-300 pc in length. In this region of the outflow, the vast majority ($\sim$70\%) of PAH emission is associated with the plumes. We show that those structures contain smaller scale "clouds" with widths that are $\sim$5-15 pc, and they are morphologically similar to the results of "cloud-crushing" simulations. We estimate the cloud-crushing time-scales of $\sim$0.5-3 Myr, depending on assumptions. We show this time scale is consistent with a picture in which these observed PAH clouds survived break-out from the disk rather than being destroyed by the hot wind. The PAH emission in both the midplane and the outflow is shown to tightly correlate with that of Pa$α$ emission (from HST/NICMOS data), at the scale of both plumes and clouds, though the ratio of PAH-to-Pa$α$ increases at further distances from the midplane. Finally, we show that the outflow PAH emission is suppressed in regions of the M82 wind that are bright in X-ray emission. Overall, our results are broadly consistent with a picture in which cold gas in galactic outflows is launched via hierarchically structured plumes, and those small scale clouds are more likely to survive the wind environment when collected into the larger plume structure. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: Submitted to MNRAS

Physical Properties of the Southwest Outflow Streamer in the Starburst Galaxy NGC 253 with ALCHEMI

Authors: Min Bao, Nanase Harada, Kotaro Kohno, Yuki Yoshimura, Fumi Egusa, Yuri Nishimura, Kunihiko Tanaka, Kouichiro Nakanishi, Sergio Martín, Jeffrey G. Mangum, Kazushi Sakamoto, Sébastien Muller, Mathilde Bouvier, Laura Colzi, Kimberly L. Emig, David S. Meier, Christian Henkel, Pedro Humire, Ko-Yun Huang, Víctor M. Rivilla, Paul van der Werf, Serena Viti

Abstract: The physical properties of galactic molecular outflows are important as they could constrain outflow formation mechanisms. We study the properties of the southwest (SW) outflow streamer including gas kinematics, optical depth, dense gas fraction, and shock strength in the central molecular zone of the starburst galaxy NGC 253. We image the molecular emission at a spatial resolution of $\sim$27 pc… ▽ More The physical properties of galactic molecular outflows are important as they could constrain outflow formation mechanisms. We study the properties of the southwest (SW) outflow streamer including gas kinematics, optical depth, dense gas fraction, and shock strength in the central molecular zone of the starburst galaxy NGC 253. We image the molecular emission at a spatial resolution of $\sim$27 pc based on data from the ALCHEMI program. We trace the kinematics of molecular gas with CO(1-0) line. We constrain the optical depth of CO emission with CO/$^{13}$CO(1-0) ratio, the dense gas fraction with HCN/CO(1-0) ratio, as well as the shock strength with SiO(2-1)/$^{13}$CO(1-0) ratio. The CO/$^{13}$CO(1-0) integrated intensity ratio is $\sim$21 in the SW streamer region, which approximates the C/$^{13}$C isotopic abundance ratio. The higher integrated intensity ratio compared to the disk can be attributed to the optically thinner environment for CO(1-0) emission inside the SW streamer. The HCN/CO(1-0) and SiO(2-1)/$^{13}$CO(1-0) integrated intensity ratios both approach $\sim$0.2 in three giant molecular clouds (GMCs) at the base of the outflow streamers, which implies the higher dense gas fraction and enhanced strength of fast shocks in those GMCs than in the disk. The contours of those two integrated intensity ratios are extended towards the directions of outflow streamers, which connects the enhanced dense gas fraction and shock strength with molecular outflow. Moreover, the molecular gas with enhanced dense gas fraction and shock strength located at the base of the SW streamer shares the same velocity with the outflow. These phenomena suggest that the star formation inside the GMCs can trigger the shocks and further drive the molecular outflow. △ Less

Submitted 6 April, 2024; originally announced April 2024.

Comments: Accepted for publication in A&A, 14 pages, 11 figures

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)

JWST Observations of Starbursts: Polycyclic Aromatic Hydrocarbon Emission at the Base of the M 82 Galactic Wind

Authors: Alberto D. Bolatto, Rebecca C. Levy, Elizabeth Tarantino, Martha L. Boyer, Deanne B. Fisher, Adam K. Leroy, Serena A. Cronin, Ralf S. Klessen, J. D. Smith, Dannielle A. Berg, Torsten Boeker, Leindert A. Boogaard, Eve C. Ostriker, Todd A. Thompson, Juergen Ott, Laura Lenkic, Laura A. Lopez, Daniel A. Dale, Sylvain Veilleux, Paul P. van der Werf, Simon C. O. Glover, Karin M. Sandstrom, Evan D. Skillman, John Chisholm, Vicente Villanueva , et al. (15 additional authors not shown)

Abstract: We present new observations of the central 1 kpc of the M 82 starburst obtained with the James Webb Space Telescope (JWST) near-infrared camera (NIRCam) instrument at a resolution ~0.05"-0.1" (~1-2 pc). The data comprises images in three mostly continuum filters (F140M, F250M, and F360M), and filters that contain [FeII] (F164N), H2 v=1-0 (F212N), and the 3.3 um PAH feature (F335M). We find promine… ▽ More We present new observations of the central 1 kpc of the M 82 starburst obtained with the James Webb Space Telescope (JWST) near-infrared camera (NIRCam) instrument at a resolution ~0.05"-0.1" (~1-2 pc). The data comprises images in three mostly continuum filters (F140M, F250M, and F360M), and filters that contain [FeII] (F164N), H2 v=1-0 (F212N), and the 3.3 um PAH feature (F335M). We find prominent plumes of PAH emission extending outward from the central starburst region, together with a network of complex filamentary substructure and edge-brightened bubble-like features. The structure of the PAH emission closely resembles that of the ionized gas, as revealed in Paschen alpha and free-free radio emission. We discuss the origin of the structure, and suggest the PAHs are embedded in a combination of neutral, molecular, and photoionized gas. △ Less

Submitted 21 April, 2024; v1 submitted 29 January, 2024; originally announced January 2024.

Comments: Submitted to The Astrophysical Journal

Hidden Gems on a Ring: Infant Massive Clusters and Their Formation Timeline Unveiled by ALMA, HST, and JWST in NGC 3351

Authors: Jiayi Sun, Hao He, Kyle Batschkun, Rebecca C. Levy, Kimberly Emig, M. Jimena Rodriguez, Hamid Hassani, Adam K. Leroy, Eva Schinnerer, Eve C. Ostriker, Christine D. Wilson, Alberto D. Bolatto, Elisabeth A. C. Mills, Erik Rosolowsky, Janice C. Lee, Daniel A. Dale, Kirsten L. Larson, David A. Thilker, Leonardo Ubeda, Bradley C. Whitmore, Thomas G. Williams, Ashley. T. Barnes, Frank Bigiel, Melanie Chevance, Simon C. O. Glover , et al. (16 additional authors not shown)

Abstract: We study young massive clusters (YMCs) in their embedded "infant" phase with $\sim0.\!^{\prime\prime}1$ ALMA, HST, and JWST observations targeting the central starburst ring in NGC 3351, a nearby Milky Way analog galaxy. Our new ALMA data reveal 18 bright and compact (sub-)millimeter continuum sources, of which 8 have counterparts in JWST images and only 6 have counterparts in HST images. Based on… ▽ More We study young massive clusters (YMCs) in their embedded "infant" phase with $\sim0.\!^{\prime\prime}1$ ALMA, HST, and JWST observations targeting the central starburst ring in NGC 3351, a nearby Milky Way analog galaxy. Our new ALMA data reveal 18 bright and compact (sub-)millimeter continuum sources, of which 8 have counterparts in JWST images and only 6 have counterparts in HST images. Based on the ALMA continuum and molecular line data, as well as ancillary measurements for the HST and JWST counterparts, we identify 14 sources as infant star clusters with high stellar and/or gas masses (${\sim}10^5\;\mathrm{M_\odot}$), small radii (${\lesssim}\,5\;\mathrm{pc}$), large escape velocities ($6{-}10\;\mathrm{km/s}$), and short free-fall times ($0.5{-}1\;\mathrm{Myr}$). Their multiwavelength properties motivate us to divide them into four categories, likely corresponding to four evolutionary stages from starless clumps to exposed HII region-cluster complexes. Leveraging age estimates for HST-identified clusters in the same region, we infer an evolutionary timeline going from $\sim$1-2 Myr before cluster formation as starless clumps, to $\sim$4-6 Myr after as exposed HII region-cluster complexes. Finally, we show that the YMCs make up a substantial fraction of recent star formation across the ring, exhibit an non-uniform azimuthal distribution without a very coherent evolutionary trend along the ring, and are capable of driving large-scale gas outflows. △ Less

Submitted 10 April, 2024; v1 submitted 25 January, 2024; originally announced January 2024.

Comments: 27 pages, 12 figures; ApJ accepted

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

Volume density structure of the NGC253 CMZ through ALCHEMI excitation analysis

Authors: Kunihiko Tanaka, Jeffrey G. Mangum, Serena Viti, Sergio Martin, Nanase Harada, Kazushi Sakamoto, Sebastien Muller, Yuki Yoshimura, Kouichiro Nakanishi, Ruben Herrero Illana, Kimberly L. Emig, S. Muhle, Hiroyuki Kaneko, Tomoka Tosaki, Erica Behrens, Victor M. Rivilla, Laura Colzi, Yuri Nishimura, P. K. Humire, Mathilde Bouvier, Ko-Yun Huang, Joshua Butterworth, David S. Meier, Paul P. van der Werf

Abstract: We present a spatially-resolved excitation analysis for the central molecular zone (CMZ) of the starburst galaxy NGC 253 using the data from the ALMA Large program ALCHEMI, whereby we explore parameters distinguishing NGC 253 from the quiescent Milky Way's Galactic Center (GC). Non-LTE analyses employing a hierarchical Bayesian framework are applied to Band 3-7 transitions from nine molecular spec… ▽ More We present a spatially-resolved excitation analysis for the central molecular zone (CMZ) of the starburst galaxy NGC 253 using the data from the ALMA Large program ALCHEMI, whereby we explore parameters distinguishing NGC 253 from the quiescent Milky Way's Galactic Center (GC). Non-LTE analyses employing a hierarchical Bayesian framework are applied to Band 3-7 transitions from nine molecular species to delineate the position-position-velocity distributions of column density ($N_\mathrm{H_2}$), volume density ($n_\mathrm{H_2}$), and temperature ($T_\mathrm{kin}$) at 27 pc resolution. Two distinct components are detected: a low-density component with $(n_\mathrm{H_2},\ T_\mathrm{kin})\sim(10^{3.3}\ \mathrm{cm}^{-3}, 85 K)$ and a high-density component with $(n_\mathrm{H_2},\ T_\mathrm{kin})\sim (10^{4.4}\ \mathrm{cm}^{-3}, 110\ \mathrm{K})$, separated at $n_\mathrm{H_2}\sim10^{3.8}\ \mathrm{cm}^{-3}$. NGC 253 has $\sim10$ times the high-density gas mass and $\sim3$ times the dense-gas mass fraction of the GC. These properties are consistent with their HCN/CO ratio but cannot alone explain the factor of $\sim30$ difference in their star formation efficiencies (SFEs), contradicting the dense-gas mass to star formation rate scaling law. The $n_\mathrm{H_2}$ histogram toward NGC 253 exhibits a shallow declining slope up to $n_\mathrm{H_2}\sim10^6\ \mathrm{cm}^{-3}$, while that of the GC steeply drops in $n_\mathrm{H_2}\gtrsim10^{4.5}\ \mathrm{cm}^{-3}$ and vanishes at $10^5\ \mathrm{cm}^{-3}$. Their dense-gas mass fraction ratio becomes consistent with their SFEs when the threshold $n_\mathrm{H_2}$ for the dense gas is taken at $\sim 10^{4.2\mbox{-}4.6}\ \mathrm{cm}^{-3}$. The rich abundance of gas above this density range in the NGC 253 CMZ, or its scarcity in the GC, is likely to be the critical difference characterizing the contrasting star formation in the centers of the two galaxies. △ Less

Submitted 20 November, 2023; originally announced November 2023.

Comments: 49 pages, 27 figures, 7 tables, accepted for publication in the Astrophysical Journal

MALS discovery of a rare HI 21-cm absorber at $z\sim1.35$: origin of the absorbing gas in powerful AGN

Authors: P. P. Deka, N. Gupta, H. W. Chen, S. D. Johnson, P. Noterdaeme, F. Combes, E. Boettcher, S. A. Balashev, K. L. Emig, G. I. G. Józsa, H. -R. Klöckner, J-. K. Krogager, E. Momjian, P. Petitjean, G. C. Rudie, J. Wagenveld, F. S. Zahedy

Abstract: We report a new, rare detection of HI 21-cm absorption associated with a quasar (only six known at $1<z<2$) here towards J2339-5523 at $z_{em}$ = 1.3531, discovered through the MeerKAT Absorption Line Survey (MALS). The absorption profile is broad ($\sim 400$ km/s), and the peak is redshifted by $\sim 200$ km/s, from $z_{em}$. Interestingly, optical/FUV spectra of the quasar from Magellan-MIKE/HST… ▽ More We report a new, rare detection of HI 21-cm absorption associated with a quasar (only six known at $1<z<2$) here towards J2339-5523 at $z_{em}$ = 1.3531, discovered through the MeerKAT Absorption Line Survey (MALS). The absorption profile is broad ($\sim 400$ km/s), and the peak is redshifted by $\sim 200$ km/s, from $z_{em}$. Interestingly, optical/FUV spectra of the quasar from Magellan-MIKE/HST-COS spectrographs do not show any absorption features associated with the 21-cm absorption. This is despite the coincident presence of the optical quasar and the radio `core' inferred from a flat spectrum component of flux density $\sim 65$ mJy at high frequencies ($>5$ GHz). The simplest explanation would be that no large HI column (N(HI)$>10^{17}$ cm$^{-2}$) is present towards the radio `core' and the optical AGN. Based on the joint optical and radio analysis of a heterogeneous sample of 16 quasars ($z_{median}$ = 0.7) and 15 radio galaxies ($z_{median}$ = 0.3) with HI 21-cm absorption detection and matched in 1.4 GHz luminosity (L$_{\rm 1.4\,GHz}$), a consistent picture emerges where quasars are primarily tracing the gas in the inner circumnuclear disk and cocoon created by the jet-ISM interaction. These exhibit L$_{1.4\,\rm GHz}$ - $ΔV_{\rm null}$ correlation, and frequent mismatch between the radio and optical spectral lines. The radio galaxies show no such correlation and likely trace the gas from the cocoon and the galaxy-wide ISM outside the photoionization cone. The analysis presented here demonstrates the potential of radio spectroscopic observations to reveal the origin of the absorbing gas associated with AGN that may be missed in optical observations. △ Less

Submitted 19 February, 2024; v1 submitted 1 November, 2023; originally announced November 2023.

Comments: 10 pages, 8 figures, accepted for publication in A&A

The FLASH pilot survey: an HI absorption search against MRC 1-Jy radio sources

Authors: J. N. H. S. Aditya, Hyein Yoon, James R. Allison, Tao An, Rajan Chhetri, Stephen J. Curran, Jeremy Darling, Kimberly L. Emig, Marcin Glowacki, Emily Kerrison, Bärbel S. Koribalski, Elizabeth K. Mahony, Vanessa A. Moss, John Morgan, Elaine M. Sadler, Roberto Soria, Renzhi Su, Simon Weng, Matthew Whiting

Abstract: We report an ASKAP search for associated HI 21-cm absorption against bright radio sources from the Molonglo Reference Catalogue (MRC) 1-Jy sample. The search uses pilot survey data from the ASKAP First Large Absorption Survey in \hi (FLASH) covering the redshift range $0.42 < z < 1.00$. From a sample of 62 MRC 1-Jy radio galaxies and quasars in this redshift range we report three new detections of… ▽ More We report an ASKAP search for associated HI 21-cm absorption against bright radio sources from the Molonglo Reference Catalogue (MRC) 1-Jy sample. The search uses pilot survey data from the ASKAP First Large Absorption Survey in \hi (FLASH) covering the redshift range $0.42 < z < 1.00$. From a sample of 62 MRC 1-Jy radio galaxies and quasars in this redshift range we report three new detections of associated HI 21-cm absorption, yielding an overall detection fraction of $1.8\%^{+4.0\%}_{-1.5\%}$. The detected systems comprise two radio galaxies (MRC 2216$-$281 at $z=0.657$ and MRC 0531$-$237 at $z=0.851$) and one quasar (MRC 2156$-$245 at $z=0.862$). The MRC 0531$-$237 absorption system is the strongest found to date, with a velocity integrated optical depth of $\rm 143.8 \pm 0.4 \ km \ s^{-1}$. All three objects with detected HI 21-cm absorption are peaked-spectrum or compact steep-spectrum (CSS) radio sources, classified based on our SED fits to the spectra. Two of them show strong interplanetary scintillation at 162 MHz, implying that the radio continuum source is smaller than 1 arcsec in size even at low frequencies. Among the class of peaked-spectrum and compact steep-spectrum radio sources, the HI detection fraction is $23\%^{+22\%}_{-13\%}$. This is consistent within $1σ$ with a detection fraction of $\approx 42\%^{+21\%}_{-15\%}$ in earlier reported GPS and CSS samples at intermediate redshifts ($0.4 < z < 1.0$). All three detections have a high 1.4 GHz radio luminosity, with MRC 0531$-$237 and MRC 2216$-$281 having the highest values in the sample, $\rm > 27.5 \ W \ Hz^{-1}$. The preponderance of extended radio sources in our sample could partially explain the overall low detection fraction, while the effects of a redshift evolution in gas properties and AGN UV luminosity on the neutral gas absorption still need to be investigated. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 28 pages, 9 figures and 7 Tables. Submitted to MNRAS

[CII] Spectral Mapping of the Galactic Wind and Starbursting Disk of M82 with SOFIA

Authors: Rebecca C. Levy, Alberto D. Bolatto, Elizabeth Tarantino, Adam K. Leroy, Lee Armus, Kimberly L. Emig, Rodrigo Herrera-Camus, Daniel P. Marrone, Elisabeth Mills, Oliver Ricken, Juergen Stutzki, Sylvain Veilleux, Fabian Walter

Abstract: M82 is an archetypal starburst galaxy in the local Universe. The central burst of star formation, thought to be triggered by M82's interaction with other members in the M81 group, is driving a multiphase galaxy-scale wind away from the plane of the disk that has been studied across the electromagnetic spectrum. Here, we present new velocity-resolved observations of the [CII] 158$μ$m line in the ce… ▽ More M82 is an archetypal starburst galaxy in the local Universe. The central burst of star formation, thought to be triggered by M82's interaction with other members in the M81 group, is driving a multiphase galaxy-scale wind away from the plane of the disk that has been studied across the electromagnetic spectrum. Here, we present new velocity-resolved observations of the [CII] 158$μ$m line in the central disk and the southern outflow of M82 using the upGREAT instrument onboard SOFIA. We also report the first detections of velocity-resolved ($ΔV = 10$ km s$^{-1}$) [CII] emission in the outflow of M82 at projected distances of $\approx1-2$ kpc south of the galaxy center. We compare the [CII] line profiles to observations of CO and HI and find that likely the majority ($>55$%) of the [CII] emission in the outflow is associated with the neutral atomic medium. We find that the fraction of [CII] actually outflowing from M82 is small compared to the bulk gas outside the midplane (which may be in a halo or tidal streamers), which has important implications for observations of [CII] outflows at higher redshift. Finally, by comparing the observed ratio of the [CII] and CO intensities to models of photodissociation regions, we estimate that the far-ultraviolet (FUV) radiation field in the disk is $\sim10^{3.5}~G_0$, in agreement with previous estimates. In the outflow, however, the FUV radiation field is 2-3 orders of magnitudes lower, which may explain the high fraction of [CII] arising from the neutral medium in the wind. △ Less

Submitted 27 September, 2023; originally announced September 2023.

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

The MeerKAT Absorption Line Survey (MALS) data release I: Stokes I image catalogs at 1-1.4 GHz

Authors: P. P. Deka, N. Gupta, P. Jagannathan, S. Sekhar, E. Momjian, S. Bhatnagar, J. Wagenveld, H. -R. Klöckner, J. Jose, S. A. Balashev, F. Combes, M. Hilton, D. Borgaonkar, A. Chatterjee, K. L. Emig, A. N. Gaunekar, G. I. G. Józsa, D. Y. Klutse, K. Knowles, J-. K. Krogager, A. Mohapatra, K. Moodley, Sébastien Muller, P. Noterdaeme, P. Petitjean , et al. (2 additional authors not shown)

Abstract: The MeerKAT Absorption Line Survey (MALS) has observed 391 telescope pointings at L-band (900 - 1670 MHz) at $δ\lesssim$ $+20°$. We present radio continuum images and a catalog of 495,325 (240,321) radio sources detected at a signal-to-noise ratio (SNR) $>$5 over an area of 2289 deg$^2$ (1132 deg$^2$) at 1006 MHz (1381 MHz). Every MALS pointing contains a central bright radio source (… ▽ More The MeerKAT Absorption Line Survey (MALS) has observed 391 telescope pointings at L-band (900 - 1670 MHz) at $δ\lesssim$ $+20°$. We present radio continuum images and a catalog of 495,325 (240,321) radio sources detected at a signal-to-noise ratio (SNR) $>$5 over an area of 2289 deg$^2$ (1132 deg$^2$) at 1006 MHz (1381 MHz). Every MALS pointing contains a central bright radio source ($S_{1\,\mathrm{GHz}} \gtrsim 0.2$ Jy). The median spatial resolution is $12^{\prime\prime}$ ($8^{\prime\prime}$). The median rms noise away from the pointing center is 25 $μ$Jy beam$^{-1}$ (22 $μ$Jy beam$^{-1}$) and is within $\sim$ 15% of the achievable theoretical sensitivity. The flux density scale ratio and astrometric accuracy deduced from multiply observed sources in MALS are less than 1% (8% scatter) and $1^{\prime\prime}$, respectively. Through comparisons with NVSS and FIRST at 1.4 GHz, we establish the catalog's accuracy in the flux density scale and astrometry to be better than 6% (15% scatter) and $0.8^{\prime\prime}$, respectively. The median flux density offset is higher (9%) for an alternate beam model based on holographic measurements. The MALS radio source counts at 1.4 GHz are in agreement with literature. We estimate spectral indices ($α$) of a subset of 125,621 sources (SNR$>$8), confirm the flattening of spectral indices with decreasing flux density and identify 140 ultra steep-spectrum ($α<-1.3$) sources as prospective high-$z$ radio galaxies ($z>2$). We have identified 1308 variable and 122 transient radio sources comprising primarily of AGN that demonstrate long-term (26 years) variability in their observed flux densities. The MALS catalogs and images are publicly available at https://mals.iucaa.in. △ Less

Submitted 23 August, 2023; originally announced August 2023.

Comments: 64 pages, 25 figures, accepted for publication in the ApJS (full version of the paper with complete tables is available at DR1 release notes)

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)

The MeerKAT Absorption Line Survey: Homogeneous continuum catalogues towards a measurement of the cosmic radio dipole

Authors: J. D. Wagenveld, H. -R. Klöckner, N. Gupta, P. P. Deka, P. Jagannathan, S. Sekhar, S. A. Balashev, E. Boettcher, F. Combes, K. L. Emig, M. Hilton, G. I. G. Józsa, P. Kamphuis, D. Y. Klutse, K. Knowles, J. -K. Krogager, A. Mohapatra, E. Momjian, K. Moodley, S. Muller, P. Petitjean, P. Salas, S. Sikhosana, R. Srianand

Abstract: The number counts of homogeneous samples of radio sources are a tried and true method of probing the large scale structure of the Universe, as most radio sources outside the galactic plane are at cosmological distances. As such they are expected to trace the cosmic radio dipole, an anisotropy analogous to the dipole seen in the cosmic microwave background (CMB). Results have shown that although th… ▽ More The number counts of homogeneous samples of radio sources are a tried and true method of probing the large scale structure of the Universe, as most radio sources outside the galactic plane are at cosmological distances. As such they are expected to trace the cosmic radio dipole, an anisotropy analogous to the dipole seen in the cosmic microwave background (CMB). Results have shown that although the cosmic radio dipole matches the direction of the CMB dipole, it has a significantly larger amplitude. This result challenges our assumption of the Universe being isotropic, which can have large repercussions for the current cosmological paradigm. Though significant measurements have been made, sensitivity to the radio dipole is generally hampered by systematic effects that can cause large biases in the measurement. Here we assess these systematics with data from the MeerKAT Absorption Line Survey (MALS). We present the analysis of ten MALS pointings, focusing on systematic effects that could lead to an inhomogeneous catalogue. We describe the calibration and creation of full band continuum images and catalogues, producing a combined catalogue containing 16,313 sources and covering 37.5 square degrees of sky down to a sensitivity of 10 $μ$Jy/beam. We measure the completeness, purity, and flux recovery statistics for these catalogues using simulated data. We investigate different source populations in the catalogues by looking at flux densities and spectral indices, and how they might influence source counts. Using the noise characteristics of the pointings, we find global measures that can be used to correct for the incompleteness of the catalogue, producing corrected number counts down to 100 - 200 $μ$Jy. We show that we can homogenise the catalogues and properly account for systematic effects. We determine that we can measure the dipole to $3σ$ significance with 100 MALS pointings. △ Less

Submitted 21 February, 2023; originally announced February 2023.

Comments: 30 pages, 25 figures. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 673, A113 (2023)

The LOFAR LBA Sky Survey II. First data release

Authors: F. de Gasperin, H. W. Edler, W. L. Williams, J. R. Callingham, B. Asabere, M. Bruggen, G. Brunetti, T. J. Dijkema, M. J. Hardcastle, M. Iacobelli, A. Offringa, M. J. Norden, H. J. A. Rottgering, T. Shimwell, R. J. van Weeren, C. Tasse, D. J. Bomans, A. Bonafede, A. Botteon, R. Cassano, K. T. Chyzy, V. Cuciti, K. L. Emig, M. Kadler, G. Miley , et al. (5 additional authors not shown)

Abstract: The Low Frequency Array (LOFAR) is the only existing radio interferometer able to observe at ultra-low frequencies (<100 MHz) with high resolution (<15") and high sensitivity (<1 mJy/beam). To exploit these capabilities, the LOFAR Surveys Key Science Project is using the LOFAR Low Band Antenna (LBA) to carry out a sensitive wide-area survey at 41-66 MHz named the LOFAR LBA Sky Survey (LoLSS). LoLS… ▽ More The Low Frequency Array (LOFAR) is the only existing radio interferometer able to observe at ultra-low frequencies (<100 MHz) with high resolution (<15") and high sensitivity (<1 mJy/beam). To exploit these capabilities, the LOFAR Surveys Key Science Project is using the LOFAR Low Band Antenna (LBA) to carry out a sensitive wide-area survey at 41-66 MHz named the LOFAR LBA Sky Survey (LoLSS). LoLSS is covering the whole northern sky above declination 24 deg with a resolution of 15" and a sensitivity of 1-2 mJy/beam (1 sigma) depending on declination, field properties, and observing conditions. Here we present the first data release. An automated pipeline was used to reduce the 95 fields included in this data release. The data reduction procedures developed for this project have general application and are currently being used to process LOFAR LBA interferometric observations. Compared to the preliminary release, direction-dependent errors have been corrected for during the calibration process. This results in a typical sensitivity of 1.55 mJy/beam at the target resolution of 15". The first data release of the LOFAR LBA Sky Survey covers 650 sqdeg in the HETDEX spring field. The resultant data products released to the community include mosaic images (I and V Stokes) of the region, and a catalogue of 42463 detected sources and related Gaussian components used to describe sources' morphologies. Separate catalogues for 6 in-band frequencies are also released. The first data release of LoLSS shows that, despite the influences of the ionosphere, LOFAR can conduct large-scale surveys in the frequency window 42-66 MHz with unprecedentedly high sensitivity and resolution. The data can be used to derive unique information on the low-frequency spectral properties of many thousands of sources with a wide range of applications in extragalactic and galactic astronomy. △ Less

Submitted 30 January, 2023; originally announced January 2023.

Comments: 20 pages, 22 figures, images and catalogues available at https://www.lofar-surveys.org/lolss.html

Journal ref: A&A 673, A165 (2023)

Discovery of Hydrogen Radio Recombination Lines at z=0.89 towards PKS 1830-211

Authors: Kimberly L. Emig, Neeraj Gupta, Pedro Salas, Sebastien Muller, Sergei A. Balashev, Francoise Combes, Emmanuel Momjian, Yiqing Song, Preshanth Jagannathan, Partha P. Deka, Gyula I. G. Jozsa, Hans-Rainer Klockner, Abhisek Mohapatra, Pasquier Noterdaeme, Patrick Petitjean, Raghunathan Srianand, Jonah D. Wagenveld

Abstract: We report the detection of stimulated hydrogen radio recombination line (RRL) emission from ionized gas in a $z=0.89$ galaxy using 580--1670 MHz observations from the MeerKAT Absorption Line Survey (MALS). The RRL emission originates in a galaxy that intercepts and strongly lenses the radio blazar PKS 1830-211 ($z=2.5$). This is the second detection of RRLs outside of the local universe and the fi… ▽ More We report the detection of stimulated hydrogen radio recombination line (RRL) emission from ionized gas in a $z=0.89$ galaxy using 580--1670 MHz observations from the MeerKAT Absorption Line Survey (MALS). The RRL emission originates in a galaxy that intercepts and strongly lenses the radio blazar PKS 1830-211 ($z=2.5$). This is the second detection of RRLs outside of the local universe and the first clearly associated with hydrogen. We detect effective H144$α$ (and H163$α$) transitions at observed frequencies of 1156 (798) MHz by stacking 17 (27) RRLs with 21$σ$ (14$σ$) significance. The RRL emission contains two main velocity components and is coincident in velocity with HI 21 cm and OH 18 cm absorption. We use the RRL spectral line energy distribution and a Bayesian analysis to constrain the density ($n_e$) and the volume-averaged pathlength ($\ell$) of the ionized gas. We determine $\log( n_e ) = 2.0_{-0.7}^{+1.0}$ cm$^{-3}$ and $\log( \ell ) = -0.7_{-1.1}^{+1.1}$ pc towards the north east (NE) lensed image, likely tracing the diffuse thermal phase of the ionized ISM in a thin disk. Towards the south west (SW) lensed image, we determine $\log( n_e ) = 3.2_{-1.0}^{+0.4}$ cm$^{-3}$ and $\log( \ell ) = -2.7_{-0.2}^{+1.8}$ pc, tracing gas that is more reminiscent of H II regions. We estimate a star formation (surface density) rate of $Σ_{\mathrm{SFR}} \sim 0.6$ M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$ or SFR $\sim 50$ M$_{\odot}$ yr$^{-1}$, consistent with a star-forming main sequence galaxy of $M_{\star} \sim 10^{11}$ M$_{\odot}$. The discovery presented here opens up the possibility of studying ionized gas at high redshifts using RRL observations from current and future (e.g., SKA and ngVLA) radio facilities. △ Less

Submitted 8 January, 2023; originally announced January 2023.

Comments: 17 pages, 10 figures, 2 tables, accepted in ApJ

V-LoTSS: The Circularly-Polarised LOFAR Two-metre Sky Survey

Authors: J. R. Callingham, T. W. Shimwell, H. K. Vedantham, C. G. Bassa, S. P. O'Sullivan, T. W. H. Yiu, S. Bloot, P. N. Best, M. J. Hardcastle, M. Haverkorn, R. D. Kavanagh, L. Lamy, B. J. S. Pope, H. J. A. Röttgering, D. J. Schwarz, C. Tasse, R. J. van Weeren, G. J. White, P. Zarka, D. J. Bomans, A. Bonafede, M. Bonato, A. Botteon, M. Bruggen, K. T. Chyży , et al. (22 additional authors not shown)

Abstract: We present the detection of 68 sources from the most sensitive radio survey in circular polarisation conducted to date. We use the second data release of the 144 MHz LOFAR Two-metre Sky Survey to produce circularly-polarised maps with median 140 $μ$Jy beam$^{-1}$ noise and resolution of 20$''$ for $\approx$27% of the northern sky (5634 deg$^{2}$). The leakage of total intensity into circular polar… ▽ More We present the detection of 68 sources from the most sensitive radio survey in circular polarisation conducted to date. We use the second data release of the 144 MHz LOFAR Two-metre Sky Survey to produce circularly-polarised maps with median 140 $μ$Jy beam$^{-1}$ noise and resolution of 20$''$ for $\approx$27% of the northern sky (5634 deg$^{2}$). The leakage of total intensity into circular polarisation is measured to be $\approx$0.06%, and our survey is complete at flux densities $\geq1$ mJy. A detection is considered reliable when the circularly-polarised fraction exceeds 1%. We find the population of circularly-polarised sources is composed of four distinct classes: stellar systems, pulsars, active galactic nuclei, and sources unidentified in the literature. The stellar systems can be further separated into chromospherically-active stars, M dwarfs, and brown dwarfs. Based on the circularly-polarised fraction and lack of an optical counterpart, we show it is possible to infer whether the unidentified sources are likely unknown pulsars or brown dwarfs. By the completion of this survey of the northern sky, we expect to detect 300$\pm$100 circularly-polarised sources. △ Less

Submitted 19 December, 2022; originally announced December 2022.

Comments: Accepted for publication in A&A, 15 pages, 8 figures, 1 table. The catalogue will be publicly available at http://lofar-surveys.org/ and via Vizier shortly

Journal ref: A&A 670, A124 (2023)

PKS1413+135: OH and HI at z = 0.247 with MeerKAT

Authors: F. Combes, N. Gupta, S. Muller, S. Balashev, P. Deka, K. Emig, H. -R. Kloeckner, D. Klutse, K. Knowles, A. Mohapatra, E. Momjian, P. Noterdaeme, P. Petitjean, P. Salas, R. Srianand, J. Wagenveld

Abstract: The BL Lac PKS 1413+135 was observed by the Large Survey Project "MeerKAT Absorption Line Survey" (MALS) in the L-band, at 1139 MHz and 1293-1379 MHz, targeting the HI and OH lines in absorption at z = 0.24671. The radio continuum is thought to come from a background object at redshift lower than 0.5, as suggested by the absence of gravitational images. The HI absorption line is detected at high s… ▽ More The BL Lac PKS 1413+135 was observed by the Large Survey Project "MeerKAT Absorption Line Survey" (MALS) in the L-band, at 1139 MHz and 1293-1379 MHz, targeting the HI and OH lines in absorption at z = 0.24671. The radio continuum is thought to come from a background object at redshift lower than 0.5, as suggested by the absence of gravitational images. The HI absorption line is detected at high signal-to-noise, with a narrow central component, and a red wing, confirming previous results. The OH 1720 MHz line is clearly detected in (maser) emission, peaking at a velocity shifted by -10 to -15 km/s with respect to the HI peak. The 1612 MHz line is lost due to radio interferences. The OH 1667 MHz main line is tentatively detected in absorption, but not the 1665 MHz one. Over 30 years, a high variability is observed in optical depths, due to the rapid changes of the line of sight, caused by the superluminal motions of the radio knots. The HI line has varied by 20 per cent in depth, while the OH-1720 MHz depth has varied by a factor 4. The position of the central velocity and the widths also varied. The absorbing galaxy is an early-type spiral (maybe S0) seen edge-on, with a prominent dust lane, covering the whole disk. Given the measured mass concentration, and the radio continuum size at centimeter wavelengths (100 mas corresponding to 400 pc at z = 0.25), the width of absorption lines from the nuclear regions are expected up to 250 km/S. The narrowness of the observed lines (< 15 km/s) suggest that the absorption comes from an outer gas ring, as frequently observed in S0 galaxies. The millimetric lines are even narrower (< 1 km/s), which corresponds to the continuum size restricted to the core. The core source is covered by individual 1 pc molecular clouds, of column density a few 10^22 cm-2, which is compatible with the gas screen detected in X-rays. △ Less

Submitted 2 January, 2023; v1 submitted 17 November, 2022; originally announced November 2022.

Comments: 7 pages, 6 figures, accepted in A&A

Journal ref: A&A 671, A43 (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

Characterizing Compact 15-33 GHz Radio Continuum Sources in Local U/LIRGs

Authors: Y. Song, S. T. Linden, A. S. Evans, L. Barcos-Munoz, E. J. Murphy, E. Momjian, T. Diaz-Santos, K. L. Larson, G. C. Privon, X. Huang, L. Armus, J. M. Mazzarella, V. U, H. Inami, V. Charmandaris, C. Ricci, K. L. Emig, J. McKinney, I. Yoon, D. Kunneriath, T. S. -Y. Lai, E. E. Rodas-Quito, A. Saravia, T. Gao, W. Meynardie , et al. (1 additional authors not shown)

Abstract: We present the analysis of $\sim 100$pc-scale compact radio continuum sources detected in 63 local (Ultra) Luminous Infrared Galaxies (U/LIRGs; $L_{\rm IR} \ge 10^{11} L_\odot$), using FWHM $\lesssim 0''.1 - 0''.2$ resolution 15 and 33 GHz observations with the Karl G. Jansky Very Large Array. We identify a total of 133 compact radio sources with effective radii of 8 - 170pc, which are classified… ▽ More We present the analysis of $\sim 100$pc-scale compact radio continuum sources detected in 63 local (Ultra) Luminous Infrared Galaxies (U/LIRGs; $L_{\rm IR} \ge 10^{11} L_\odot$), using FWHM $\lesssim 0''.1 - 0''.2$ resolution 15 and 33 GHz observations with the Karl G. Jansky Very Large Array. We identify a total of 133 compact radio sources with effective radii of 8 - 170pc, which are classified into four main categories -- "AGN" (AGN), "AGN/SBnuc" (AGN-starburst composite nucleus), "SBnuc" (starburst nucleus) and "SF" (star-forming clumps) -- based on ancillary datasets and the literature. We find that "AGN" and "AGN/SBnuc" more frequently occur in late-stage mergers and have up to 3 dex higher 33 GHz luminosities and surface densities compared with "SBnuc" and "SF", which may be attributed to extreme nuclear starburst and/or AGN activity in the former. Star formation rates (SFRs) and surface densities ($Σ_{\rm SFR}$) are measured for "SF" and "SBnuc" using both the total 33 GHz continuum emission (SFR $\sim 0.14 - 13$ M$_\odot$ yr$^{-1}$, $Σ_{\rm SFR} \sim 13 - 1600$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$) and the thermal free-free emission from HII regions (median SFR$_{\rm th} \sim 0.4$ M$_\odot$ yr$^{-1}$, $Σ_{\rm SFR_{th}} \sim 44$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$). These values are 1 - 2 dex higher than those measured for similar-sized clumps in nearby normal (non-U/LIRGs). The latter also have much flatter median 15 - 33 GHz spectral index ($\sim -0.08$) compared with "SBnuc" and "SF" ($\sim -0.46$), which may reflect higher non-thermal contribution from supernovae and/or ISM densities in local U/LIRGs that directly result from and/or lead to their extreme star-forming activities on 100\,pc scales. △ Less

Submitted 8 September, 2022; originally announced September 2022.

Comments: Accepted for publication in ApJ

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

Emergence of a new HI 21-cm absorption component at z~1.1726 towards the gamma-ray blazar PKS~2355-106

Authors: Raghunathan Srianand, Neeraj Gupta, Patrick Petitjean, Emmanuel Momjian, Sergei A. Balashev, Francoise Combes, Hsiao-Wen Chen, Jens-Kristian Krogager, Pasquier Noterdaeme, Hadi Rahmani, Andrew J. Baker, Kimberly L. Emig, Gyula I. G. Jozsa, Hans-Rainer Kloeckner, Kavilan Moodley

Abstract: We report the emergence of a new HI 21-cm absorption at z_abs = 1.172635 in the damped Lyman-alpha absorber (DLA) towards the gamma-ray blazar PKS 2355-106 (z_em~1.639) using science verification observations (June 2020) from the MeerKAT Absorption Line Survey (MALS). Since 2006, this DLA is known to show a narrow HI 21-cm absorption at z_abs = 1.173019 coinciding with a distinct metal absorption… ▽ More We report the emergence of a new HI 21-cm absorption at z_abs = 1.172635 in the damped Lyman-alpha absorber (DLA) towards the gamma-ray blazar PKS 2355-106 (z_em~1.639) using science verification observations (June 2020) from the MeerKAT Absorption Line Survey (MALS). Since 2006, this DLA is known to show a narrow HI 21-cm absorption at z_abs = 1.173019 coinciding with a distinct metal absorption line component. We do not detect significant HI 21-cm optical depth variations from this known HI component. A high resolution optical spectrum (August 2010) shows a distinct Mg I absorption at the redshift of the new HI 21-cm absorber. However, this component is not evident in the profiles of singly ionized species. We measure the metallicity ([Zn/H] = -(0.77\pm0.11) and [Si/H]= -(0.96\pm0.11)) and depletion ([Fe/Zn] = -(0.63\pm0.16)) for the full system. Using the apparent column density profiles of Si II, Fe II and Mg I we show that the depletion and the N(Mg I)/N(Si II) column density ratio systematically vary across the velocity range. The region with high depletion tends to have slightly larger N(Mg I)/N(Si II) ratio. The two HI 21-cm absorbers belong to this velocity range. The emergence of z_abs = 1.172635 can be understood if there is a large optical depth gradient over a length scale of ~0.35 pc. However, the gas producing the z_abs = 1.173019 component must be nearly uniform over the same scale. Systematic uncertainties introduced by the absorption line variability has to be accounted for in experiments measuring the variations of fundamental constants and cosmic acceleration even when the radio emission is apparently compact as in PKS 2355-106. △ Less

Submitted 5 July, 2022; originally announced July 2022.

Comments: 8 pages, 7 figures and accepted for publication in MNRAS

The Morpho-Kinematic Architecture of Super Star Clusters in the Center of NGC253

Authors: Rebecca C. Levy, Alberto D. Bolatto, Adam K. Leroy, Mattia C. Sormani, Kimberly L. Emig, Mark Gorski, Laura Lenkić, Elisabeth A. C. Mills, Elizabeth Tarantino, Peter Teuben, Sylvain Veilleux, Fabian Walter

Abstract: The center of the nearby galaxy NGC\,253 hosts a population of more than a dozen super star clusters (SSCs) which are still in the process of forming. The majority of the star formation of the burst is concentrated in these SSCs, and the starburst is powering a multiphase outflow from the galaxy. In this work, we measure the 350~GHz dust continuum emission towards the center of NGC\,253 at 47~mill… ▽ More The center of the nearby galaxy NGC\,253 hosts a population of more than a dozen super star clusters (SSCs) which are still in the process of forming. The majority of the star formation of the burst is concentrated in these SSCs, and the starburst is powering a multiphase outflow from the galaxy. In this work, we measure the 350~GHz dust continuum emission towards the center of NGC\,253 at 47~milliarcsecond (0.8 pc) resolution using data from the Atacama Large Millimeter/submillimeter Array (ALMA). We report the detection of 350~GHz (dust) continuum emission in the outflow for the first time, associated with the prominent South-West streamer. In this feature, the dust emission has a width of $\approx$~8~pc, is located at the outer edge of the CO emission, and corresponds to a molecular gas mass of $\sim~(8-17)\times10^6$~M$_\odot$. In the starburst nucleus, we measure the resolved radial profiles, sizes, and molecular gas masses of the SSCs. Compared to previous work at somewhat lower spatial resolution, the SSCs here break apart into smaller substructures with radii $0.4-0.7$~pc. In projection, the SSCs, dust, and dense molecular gas appear to be arranged as a thin, almost linear, structure roughly 155~pc in length. The morphology and kinematics of this structure can be well explained as gas following $x_2$ orbits at the center of a barred potential. We constrain the morpho-kinematic arrangement of the SSCs themselves, finding that an elliptical, angular momentum-conserving ring is a good description of the both morphology and kinematics of the SSCs. △ Less

Submitted 22 June, 2022; v1 submitted 9 June, 2022; originally announced June 2022.

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

Filamentary structures of ionized gas in Cygnus X

Authors: K. L. Emig, G. J. White, P. Salas, R. L. Karim, R. J. van Weeren, P. J. Teuben, A. Zavagno, P. Chiu, M. Haverkorn, J. B. R. Oonk, E. Orrú, I. M. Polderman, W. Reich, H. J. A. Röttgering, A. G. G. M. Tielens

Abstract: Ionized gas probes the influence of massive stars on their environment. The Cygnus X region (d~1.5 kpc) is one of the most massive star forming complexes in our Galaxy, in which the Cyg OB2 association (age of 3-5 Myr and stellar mass $2 \times 10^{4}$ M$_{\odot}$) has a dominant influence. We observe the Cygnus X region at 148 MHz using the Low Frequency Array (LOFAR) and take into account short-… ▽ More Ionized gas probes the influence of massive stars on their environment. The Cygnus X region (d~1.5 kpc) is one of the most massive star forming complexes in our Galaxy, in which the Cyg OB2 association (age of 3-5 Myr and stellar mass $2 \times 10^{4}$ M$_{\odot}$) has a dominant influence. We observe the Cygnus X region at 148 MHz using the Low Frequency Array (LOFAR) and take into account short-spacing information during image deconvolution. Together with data from the Canadian Galactic Plane Survey, we investigate the morphology, distribution, and physical conditions of low-density ionized gas in a $4^{\circ} \times 4^{\circ}$ (100 pc $\times$ 100 pc) region at a resolution of 2' (0.9 pc). The Galactic radio emission in the region analyzed is almost entirely thermal (free-free) at 148 MHz, with emission measures of $10^3 < EM~{\rm[pc~cm^{-6}]} < 10^6$. As filamentary structure is a prominent feature of the emission, we use DisPerSE and FilChap to identify filamentary ridges and characterize their radial ($EM$) profiles. The distribution of radial profiles has a characteristic width of 4.3 pc and a power-law distribution ($β= -1.8 \pm 0.1$) in peak $EM$ down to our completeness limit of 4200 pc cm$^{-6}$. The electron densities of the filamentary structure range from $10 < n_e~{\rm[cm^{-3}]} < 400$ with a median value of 35 cm$^{-3}$, remarkably similar to [N II] surveys of ionized gas. Cyg OB2 may ionize at most two-thirds of the total ionized gas and the ionized gas in filaments. More than half of the filamentary structures are likely photoevaporating surfaces flowing into a surrounding diffuse (~5 cm$^{-3}$) medium. However, this is likely not the case for all ionized gas ridges. A characteristic width in the distribution of ionized gas points to the stellar winds of Cyg OB2 creating a fraction of the ionized filaments through swept-up ionized gas or dissipated turbulence. △ Less

Submitted 18 May, 2022; originally announced May 2022.

Comments: 19 pages, 14 figures, 1 table; accepted for publication in A&A

Journal ref: A&A 664, A88 (2022)

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

The LOFAR Two-metre Sky Survey -- V. Second data release

Authors: T. W. Shimwell, M. J. Hardcastle, C. Tasse, P. N. Best, H. J. A. Röttgering, W. L. Williams, A. Botteon, A. Drabent, A. Mechev, A. Shulevski, R. J. van Weeren, L. Bester, M. Brüggen, G. Brunetti, J. R. Callingham, K. T. Chyży, J. E. Conway, T. J. Dijkema, K. Duncan, F. de Gasperin, C. L. Hale, M. Haverkorn, B. Hugo, N. Jackson, M. Mevius , et al. (81 additional authors not shown)

Abstract: In this data release from the LOFAR Two-metre Sky Survey (LoTSS) we present 120-168MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44$^\circ$30' and 1h00m +28$^\circ$00' and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451hrs (7.6PB) of LOFAR High Band Antenna data which were corrected for th… ▽ More In this data release from the LOFAR Two-metre Sky Survey (LoTSS) we present 120-168MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44$^\circ$30' and 1h00m +28$^\circ$00' and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451hrs (7.6PB) of LOFAR High Band Antenna data which were corrected for the direction-independent instrumental properties as well as direction-dependent ionospheric distortions during extensive, but fully automated, data processing. A catalogue of 4,396,228 radio sources is derived from our total intensity (Stokes I) maps, where the majority of these have never been detected at radio wavelengths before. At 6" resolution, our full bandwidth Stokes I continuum maps with a central frequency of 144MHz have: a median rms sensitivity of 83$μ$Jy/beam; a flux density scale accuracy of approximately 10%; an astrometric accuracy of 0.2"; and we estimate the point-source completeness to be 90% at a peak brightness of 0.8mJy/beam. By creating three 16MHz bandwidth images across the band we are able to measure the in-band spectral index of many sources, albeit with an error on the derived spectral index of +/-0.2 which is a consequence of our flux-density scale accuracy and small fractional bandwidth. Our circular polarisation (Stokes V) 20" resolution 120-168MHz continuum images have a median rms sensitivity of 95$μ$Jy/beam, and we estimate a Stokes I to Stokes V leakage of 0.056%. Our linear polarisation (Stokes Q and Stokes U) image cubes consist of 480 x 97.6 kHz wide planes and have a median rms sensitivity per plane of 10.8mJy/beam at 4' and 2.2mJy/beam at 20"; we estimate the Stokes I to Stokes Q/U leakage to be approximately 0.2%. Here we characterise and publicly release our Stokes I, Q, U and V images in addition to the calibrated uv-data. △ Less

Submitted 23 February, 2022; originally announced February 2022.

Comments: 23 figures, 1 table and 29 pages. The catalogues, images and uv-data associated with this data release are publicly available via https://lofar-surveys.org/

First extragalactic detection of a phosphorus-bearing molecule with ALCHEMI: phosphorus nitride (PN)

Authors: D. Haasler, V. M. Rivilla, S. Martín, J. Holdship, S. Viti, N. Harada, J. Mangum, K. Sakamoto, S. Muller, K. Tanaka, Y. Yoshimura, K. Nakanishi, L. Colzi, L. Hunt, K. L. Emig, R. Aladro, P. Humire, C. Henkel, P. van der Werf

Abstract: Phosphorus (P) is a crucial element for life given its central role in several biomolecules. P-bearing molecules have been discovered in different regions of the Milky Way, but not yet towards an extragalactic environment. We have searched for P-bearing molecules towards the nearby starburst Galaxy NGC 253. Using observations from the ALMA Comprehensive High-resolution Extragalactic Molecular Inve… ▽ More Phosphorus (P) is a crucial element for life given its central role in several biomolecules. P-bearing molecules have been discovered in different regions of the Milky Way, but not yet towards an extragalactic environment. We have searched for P-bearing molecules towards the nearby starburst Galaxy NGC 253. Using observations from the ALMA Comprehensive High-resolution Extragalactic Molecular Inventory (ALCHEMI) project, we used the MADCUBA package to model the emission of P-bearing molecules assuming Local Thermodynamic Equilibrium (LTE) conditions. We have also performed a non-LTE analysis using SpectralRadex. We report the detection of a P-bearing molecule, phosphorus nitride (PN), for the first time in an extragalactic environment, towards two giant molecular clouds (GMCs) of NGC 253. The LTE analysis yields total PN beam-averaged column densities $N$=(1.20$\pm$0.09)$\times$10$^{13}$ cm$^{-2}$ and $N$=(6.5$\pm$1.6)$\times$10$^{12}$ cm$^{-2}$, which translate into abundances with respect to H$_2$ of $χ$=(8.0$\pm$1.0)$\times$10$^{-12}$ and $χ$=(4.4$\pm$1.2)$\times$10$^{-12}$. We derived a low excitation temperature of $T_{\rm ex}$=(4.4$\pm$1.3) K towards the GMC with the brightest PN emission, which indicates that PN is sub-thermally excited. The non-LTE analysis results in column densities consistent with the LTE values. We have also searched for other P-bearing molecules (PO, PH$_{3}$, CP and CCP), and upper limits were derived. The derived PO/PN ratios are $<$1.3 and $<$1.7. The abundance ratio between PN and the shock-tracer SiO derived towards NGC 253 follows the same trend previously found towards Galactic sources. Comparison of the observations with chemical models indicates that the derived molecular abundances of PN in NGC 253 can be explained by shock-driven chemistry followed by cosmic-ray-driven photochemistry. △ Less

Submitted 15 January, 2022; v1 submitted 9 December, 2021; originally announced December 2021.

Comments: 14 pages, 14 figures, 4 tables. Accepted for publication in A&A

Journal ref: A&A 659, A158 (2022)

Arecibo-Green Bank-LOFAR Carbon Radio Recombination Line observations toward cold HI Clouds

Authors: D. Anish Roshi, W. M. Peters, K. L. Emig, P. Salas, J. B. R. Oonk, M. E. Lebrón, J. M. Dickey

Abstract: We present results from a search for radio recombination lines in three HI self-absorbing (HISA) clouds at 750 MHz and 321 MHz with the Robert C. Byrd Green Bank Telescope (GBT), and in three Galactic Plane positions at 327 MHz with the Arecibo Telescope. We detect Carbon Recombination Lines (CRRLs) in the direction of DR4 and DR21, as well as in the galactic plane position G34.94+0.0. We addition… ▽ More We present results from a search for radio recombination lines in three HI self-absorbing (HISA) clouds at 750 MHz and 321 MHz with the Robert C. Byrd Green Bank Telescope (GBT), and in three Galactic Plane positions at 327 MHz with the Arecibo Telescope. We detect Carbon Recombination Lines (CRRLs) in the direction of DR4 and DR21, as well as in the galactic plane position G34.94+0.0. We additionally detect Hydrogen Recombination Lines (HRRLs) in emission in five of the six sightlines, and a Helium line at 750 MHz towards DR21. Combining our new data with 150 MHz LOFAR detections of CRRL absorption towards DR4 and DR21, we estimate the electron densities of the line forming regions by modeling the line width as a function of frequency. The estimated densities are in the range 1.4 $\rightarrow$ 6.5 cm$^{-3}$ towards DR4, for electron temperatures 200 $\rightarrow$ 20 K. A dual line forming region with densities between 3.5 $\rightarrow$ 24 cm$^{-3}$ and 0.008 $\rightarrow$ 0.3 cm$^{-3}$ could plausibly explain the observed line width as a function of frequency on the DR21 sightline. The central velocities of the CRRLs compare well with CO emission and HISA lines in these directions. The cloud densities estimated from the CO lines are smaller (at least a factor of 5) than those of the CRRL forming regions. It is likely that the CRRL forming and HI self-absorbing gas is located in a denser, shocked region either at the boundary of or within the CO emitting cloud. △ Less

Submitted 2 November, 2021; originally announced November 2021.

Comments: 24 pages, 10 figures, to appear in ApJ

The First Large Absorption Survey in HI (FLASH): I. Science Goals and Survey Design

Authors: J. R. Allison, E. M. Sadler, A. D. Amaral, T. An, S. J. Curran, J. Darling, A. C. Edge, S. L. Ellison, K. L. Emig, B. M. Gaensler, L. Garratt-Smithson, M. Glowacki, K. Grasha, B. S. Koribalski, C. del P. Lagos, P. Lah, E. K. Mahony, S. A. Mao, R. Morganti, V. A. Moss, M. Pettini, K. A. Pimbblet, C. Power, P. Salas, L. Staveley-Smith , et al. (5 additional authors not shown)

Abstract: We describe the scientific goals and survey design of the First Large Absorption Survey in HI (FLASH), a wide field survey for 21-cm line absorption in neutral atomic hydrogen (HI) at intermediate cosmological redshifts. FLASH will be carried out with the Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope and is planned to cover the sky south of $δ\approx +40$deg at frequencies b… ▽ More We describe the scientific goals and survey design of the First Large Absorption Survey in HI (FLASH), a wide field survey for 21-cm line absorption in neutral atomic hydrogen (HI) at intermediate cosmological redshifts. FLASH will be carried out with the Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope and is planned to cover the sky south of $δ\approx +40$deg at frequencies between 711.5 and 999.5MHz. At redshifts between $z = 0.4$ and $1.0$ (look back times of 4 - 8Gyr), the HI content of the Universe has been poorly explored due to the difficulty of carrying out radio surveys for faint 21-cm line emission and, at ultra-violet wavelengths, space-borne searches for Damped Lyman-$α$ absorption in quasar spectra. The ASKAP wide field of view and large spectral bandwidth, in combination with a radio-quiet site, will enable a search for absorption lines in the radio spectra of bright continuum sources over 80% of the sky. This survey is expected to detect at least several hundred intervening 21-cm absorbers, and will produce an HI-absorption-selected catalogue of galaxies rich in cool, star-forming gas, some of which may be concealed from optical surveys. Likewise, at least several hundred associated 21-cm absorbers are expected to be detected within the host galaxies of radio sources at $0.4 < z < 1.0$, providing valuable kinematical information for models of gas accretion and jet-driven feedback in radio-loud active galactic nuclei. FLASH will also detect OH 18-cm absorbers in diffuse molecular gas, megamaser OH emission, radio recombination lines, and stacked HI emission. △ Less

Submitted 23 January, 2022; v1 submitted 1 October, 2021; originally announced October 2021.

Comments: 31 pages, 10 figures, 6 tables, accepted for publication in PASA

ALMA Imaging of a Galactic Molecular Outflow in NGC4945

Authors: Alberto D. Bolatto, Adam K. Leroy, Rebecca C. Levy, David S. Meier, Elisabeth A. C. Mills, Todd A. Thompson, Kimberly L. Emig, Sylvain Veilleux, Juergen Ott, Mark Gorski, Fabian Walter, Laura A. Lopez, Laura Lenkic

Abstract: We present the ALMA detection of molecular outflowing gas in the central regions of NGC4945, one of the nearest starbursts and also one of the nearest hosts of an active galactic nucleus (AGN). We detect four outflow plumes in CO (3-2) at ~0.3" resolution that appear to correspond to molecular gas located near the edges of the known ionized outflow cone and its (unobserved) counterpart behind the… ▽ More We present the ALMA detection of molecular outflowing gas in the central regions of NGC4945, one of the nearest starbursts and also one of the nearest hosts of an active galactic nucleus (AGN). We detect four outflow plumes in CO (3-2) at ~0.3" resolution that appear to correspond to molecular gas located near the edges of the known ionized outflow cone and its (unobserved) counterpart behind the disk. The fastest and brightest of these plumes has emission reaching observed line-of-sight projected velocities of over 450 km/s beyond systemic, equivalent to an estimated physical outflow velocity v>600 km/s for the fastest emission. Most of these plumes have corresponding emission in HCN or HCO+ (4-3). We discuss a kinematic model for the outflow emission where the molecular gas has the geometry of the ionized gas cone and shares the rotation velocity of the galaxy when ejected. We use this model to explain the velocities we observe, constrain the physical speed of the ejected material, and account for the fraction of outflowing gas that is not detected due to confusion with the galaxy disk. We estimate a total molecular mass outflow rate dMmol/dt~20 Msun/yr flowing through a surface within 100 pc of the disk midplane, likely driven by a combination of the central starburst and AGN. △ Less

Submitted 1 October, 2021; v1 submitted 21 September, 2021; originally announced September 2021.

Comments: 20 pages, 12 figures, accepted by The Astrophysical Journal on September 28, 2021

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)

Searching for a Galactic component in the IceCube track-like neutrino events

Authors: Gregory S. Vance, Kimberly L. Emig, Cecilia Lunardini, Rogier A. Windhorst

Abstract: Searches for spatial associations between high-energy neutrinos observed at the IceCube Neutrino Observatory and known astronomical objects may hold the key to establishing the neutrinos' origins and the origins of hadronic cosmic rays. While extragalactic sources like the blazar TXS 0506+056 merit significant attention, Galactic sources may also represent part of the puzzle. Here, we explore whet… ▽ More Searches for spatial associations between high-energy neutrinos observed at the IceCube Neutrino Observatory and known astronomical objects may hold the key to establishing the neutrinos' origins and the origins of hadronic cosmic rays. While extragalactic sources like the blazar TXS 0506+056 merit significant attention, Galactic sources may also represent part of the puzzle. Here, we explore whether open clusters and supernova remnants in the Milky Way contribute measurably to the IceCube track-like neutrino events above 200 TeV. By searching for positional coincidences with catalogs of known astronomical objects, we can identify and investigate neutrino events whose origins are potentially Galactic. We use Monte Carlo randomization together with models of the Galactic plane in order to determine whether these coincidences are more likely to be causal associations or random chance. In all analyses presented, the number of coincidences detected was found to be consistent with the null hypothesis of chance coincidence. Our results imply that the combined contribution of Galactic open clusters and supernova remnants to the track-like neutrino events detected at IceCube is well under 30%. This upper limit is compatible with the results presented in other Galactic neutrino studies. △ Less

Submitted 3 August, 2021; originally announced August 2021.

Comments: 26 pages, 4 figures, 7 tables; to be submitted to JCAP

The ionization fraction in OMC-2 and OMC-3

Authors: P. Salas, M. R. Rugel, K. L. Emig, J. Kauffmann, K. M. Menten, F. Wyrowski, A. G. G. M. Tielens

Abstract: The electron density ($n_{e^{-}}$) plays an important role in setting the chemistry and physics of the interstellar medium. However, measurements of $n_{e^{-}}$ in neutral clouds have been directly obtained only toward a few lines of sight or they rely on indirect determinations. We use carbon radio recombination lines and the far-infrared lines of C$^{+}$ to directly measure $n_{e^{-}}$ and the g… ▽ More The electron density ($n_{e^{-}}$) plays an important role in setting the chemistry and physics of the interstellar medium. However, measurements of $n_{e^{-}}$ in neutral clouds have been directly obtained only toward a few lines of sight or they rely on indirect determinations. We use carbon radio recombination lines and the far-infrared lines of C$^{+}$ to directly measure $n_{e^{-}}$ and the gas temperature in the envelope of the integral shaped filament (ISF) in the Orion A molecular cloud. We observed the C$102α$ and C$109α$ carbon radio recombination lines (CRRLs) using the Effelsberg 100m telescope at ~2' resolution toward five positions in OMC-2 and OMC-3. Since the CRRLs have similar line properties, we averaged them to increase the signal-to-noise ratio of the spectra. We compared the intensities of the averaged CRRLs, and the 158 μm-[CII] and [$^{13}$CII] lines to the predictions of a homogeneous model for the C$^{+}$/C interface in the envelope of a molecular cloud and from this comparison we determined the electron density, temperature and C$^{+}$ column density of the gas. We detect the CRRLs toward four positions, where their velocity and widths (FWHM 2.3 km s$^{-1}$) confirms that they trace the envelope of the ISF. Toward two positions we detect the CRRLs, and the [CII] and [$^{13}$CII] lines with a signal-to-noise ratio >5, and we find $n_{e^{-}}=0.65\pm0.12$ cm$^{-3}$ and $0.95\pm0.02$ cm$^{-3}$, which corresponds to a gas density $n_{H}\approx5\times10^{3}$ cm$^{-3}$ and a thermal pressure of $p_{th}\approx4\times10^{5}$ K cm$^{-3}$. We also constrained the ionization fraction in the denser portions of the molecular cloud using the HCN(1-0) and C$_{2}$H(1-0) lines to $x(e^{-})<3\times10^{-6}$. The derived electron densities and ionization fraction imply that $x(e^{-})$ drops by a factor >100 between the C$^{+}$ layer and the regions probed by HCN(1-0). △ Less

Submitted 14 July, 2021; originally announced July 2021.

Comments: 14 pages, 6 figures, 8 tables, accepted for publication in A&A

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)

Clustered Star Formation in the center of NGC 253 Contributes to Driving the Ionized Nuclear Wind

Authors: E. A. C. Mills, M. Gorski, K. L. Emig, A. D. Bolatto, R. C. Levy, A. K. Leroy, A. Ginsburg, J. D. Henshaw, L. K. Zschaechner, S. Veilleux, K. Tanaka, D. S. Meier, F. Walter, N. Krieger, J. Ott

Abstract: We present new 3 mm observations of the ionized gas toward the nuclear starburst in the nearby (D ~ 3.5 Mpc) galaxy NGC 253. With ALMA, we detect emission from the H40-alpha and He40-alpha lines in the central 200 pc of this galaxy on spatial scales of ~4 pc. The recombination line emission primarily originates from a population of approximately a dozen embedded super star clusters in the early st… ▽ More We present new 3 mm observations of the ionized gas toward the nuclear starburst in the nearby (D ~ 3.5 Mpc) galaxy NGC 253. With ALMA, we detect emission from the H40-alpha and He40-alpha lines in the central 200 pc of this galaxy on spatial scales of ~4 pc. The recombination line emission primarily originates from a population of approximately a dozen embedded super star clusters in the early stages of formation. We find that emission from these clusters is characterized by electron temperatures ranging from 7000-10000 K and measure an average singly-ionized helium abundance <Y+> = 0.25 +/- 0.06, both of which are consistent with values measured for HII regions in the center of the Milky Way. We also report the discovery of unusually broad-linewidth recombination line emission originating from seven of the embedded clusters. We suggest that these clusters contribute to the launching of the large-scale hot wind observed to emanate from the central starburst. Finally, we use the measured recombination line fluxes to improve the characterization of overall embedded cluster properties, including the distribution of cluster masses and the fractional contribution of the clustered star formation to the total starburst, which we estimate is at least 50%. △ Less

Submitted 28 June, 2021; originally announced June 2021.

Comments: ApJ Accepted. 20 pages, 8 figures

NOEMA High Fidelity Imaging of the Molecular Gas in and around M82

Authors: Nico Krieger, Fabian Walter, Alberto D. Bolatto, Pierre Guillard, Matthew Lehnert, Adam Leroy, Jérôme Pety, Kimberly L. Emig, Rebecca C. Levy, Melanie Krips, Hans-Walter Rix, Dragan Salak, Axel Weiss, Sylvain Veilleux

Abstract: We present a 154 pointing IRAM NOEMA mosaic of the CO(1-0) line emission in and around the nearby starburst galaxy M82. The observations, complemented by zero--spacing observations, reach a spatial resolution of $\sim$30 pc ($\sim 1.9^{\prime\prime}$) at 5.0 km s$^{-1}$ spectral resolution, sufficient to resolve the molecular gas in the central starburst disk, the outflow, as well as the tidal str… ▽ More We present a 154 pointing IRAM NOEMA mosaic of the CO(1-0) line emission in and around the nearby starburst galaxy M82. The observations, complemented by zero--spacing observations, reach a spatial resolution of $\sim$30 pc ($\sim 1.9^{\prime\prime}$) at 5.0 km s$^{-1}$ spectral resolution, sufficient to resolve the molecular gas in the central starburst disk, the outflow, as well as the tidal streamers. The resulting moment and peak brightness maps show a striking amount of structure. Using a clump decomposition algorithm, we analyse the physical properties (e.g., radii $R$, line widths $σ$, and masses $M$) of $\sim2000$ molecular clouds. To first order, the clouds' properties are very similar, irrespective of their environment. This also holds for the size-line width relations of the clouds. The distribution of clouds in the $σ^2/R$ vs. column density $Σ$ space suggests that external pressure does not play a significant role in setting their physical parameters in the outflow and the streamers. We find that the clouds in the streamers stay approximately constant in size ($R \sim 50$ pc) and mass ($M \sim 10^5$ M$_\odot$) and do not vary with their projected distance from M82's center. The clouds in the outflow, on the other hand, appear to decrease in size and mass with distance towards the Southern outflow. The reduction in the molecular gas luminosity could be indicative of cloud evaporation of embedded clouds in the hot outflow. △ Less

Submitted 16 May, 2021; originally announced May 2021.

Comments: 16 pages, 10 figures, 2 tables; accepted for publication in ApJ

The LOFAR LBA Sky Survey I. survey description and preliminary data release

Authors: F. de Gasperin, W. L. Williams, P. Best, M. Bruggen, G. Brunetti, V. Cuciti, T. J. Dijkema, M. J. Hardcastle, M. J. Norden, A. Offringa, T. Shimwell, R. van Weeren, D. Bomans, A. Bonafede, A. Botteon, J. R. Callingham, R. Cassano, K. T. Chyzy, K. L. Emig, H. Edler, M. Haverkorn, G. Heald, V. Heesen, M. Iacobelli, H. T. Intema , et al. (16 additional authors not shown)

Abstract: LOFAR is the only radio telescope that is presently capable of high-sensitivity, high-resolution (<1 mJy/b and <15") observations at ultra-low frequencies (<100 MHz). To utilise these capabilities, the LOFAR Surveys Key Science Project is undertaking a large survey to cover the entire northern sky with Low Band Antenna (LBA) observations. The LOFAR LBA Sky Survey (LoLSS) aims to cover the entire n… ▽ More LOFAR is the only radio telescope that is presently capable of high-sensitivity, high-resolution (<1 mJy/b and <15") observations at ultra-low frequencies (<100 MHz). To utilise these capabilities, the LOFAR Surveys Key Science Project is undertaking a large survey to cover the entire northern sky with Low Band Antenna (LBA) observations. The LOFAR LBA Sky Survey (LoLSS) aims to cover the entire northern sky with 3170 pointings in the frequency range 42-66 MHz, at a resolution of 15 arcsec and at a sensitivity of 1 mJy/b. Here we outline the survey strategy, the observational status, the current calibration techniques, and briefly describe several scientific motivations. We also describe the preliminary public data release. The preliminary images were produced using a fully automated pipeline that aims to correct all direction-independent effects in the data. Whilst the direction-dependent effects, such as those from the ionosphere, are not yet corrected, the images presented in this work are still 10 times more sensitive than previous surveys available at these low frequencies. The preliminary data release covers 740 sqdeg around the HETDEX spring field region at a resolution of 47" with a median noise level of 5 mJy/b. The images and the catalogue with 25,247 sources are publicly released. We demonstrate that the system is capable of reaching an rms noise of 1 mJy/b and the resolution of 15" once direction-dependent effects are corrected for. LoLSS will provide the ultra-low-frequency information for hundreds of thousands of radio sources, providing critical spectral information and producing a unique dataset that can be used for a wide range of science topics such as: the search for high redshift galaxies and quasars, the study of the magnetosphere of exoplanets, and the detection of the oldest populations of cosmic-rays in galaxies, clusters of galaxies, and from AGN activity. △ Less

Submitted 18 February, 2021; originally announced February 2021.

Comments: 19 pages, 14 figures, Accepted A&A, catalogue and images on www.lofar-surveys.org

Extremely deep 150 MHz source counts from the LoTSS Deep Fields

Authors: S. Mandal, I. Prandoni, M. J. Hardcastle, T. W. Shimwell, H. T. Intema, C. Tasse, R. J. van Weeren, H. Algera, K. L. Emig, H. J. A. Röttgering, D. J. Schwarz, T. M. Siewert, P. N. Best, M. Bonato, M. Bondi, M. J. Jarvis, R. Kondapally, S. K. Leslie, V. H. Mahatma, J. Sabater, E. Retana-Montenegro, W. L. Williams

Abstract: With the advent of new generation low-frequency telescopes, such as the LOw Frequency ARray (LOFAR), and improved calibration techniques, we have now started to unveil the sub GHz radio sky with unprecedented depth and sensitivity. The LOFAR Two Meter Sky Survey (LoTSS) is an ongoing project in which the whole northern radio sky will be observed at 150 MHz with a sensitivity better than 100 $μ$Jy… ▽ More With the advent of new generation low-frequency telescopes, such as the LOw Frequency ARray (LOFAR), and improved calibration techniques, we have now started to unveil the sub GHz radio sky with unprecedented depth and sensitivity. The LOFAR Two Meter Sky Survey (LoTSS) is an ongoing project in which the whole northern radio sky will be observed at 150 MHz with a sensitivity better than 100 $μ$Jy beam$^{-1}$ at a resolution of \asec{6}. Additionally, deeper observations are planned to cover smaller areas with higher sensitivity. The Lockman Hole, the Boötes and the Elais-N1 regions are among the most well known northern extra-galactic fields, and the deepest of the LoTSS Deep Fields so far. We exploit these deep observations to derive the deepest radio source counts at 150~MHz to date. Our counts are in broad agreement with those from the literature, and show the well known upturn at $\leq$ few mJy, mainly associated with the emergence of the star-forming galaxy population. More interestingly, our counts show for the first time a very pronounced drop around S$\sim$2 mJy, which results in a prominent `bump' at sub-mJy flux densities. Such a feature was not observed in previous counts' determinations (neither at 150 MHz nor at higher frequency). While sample variance can play a role in explaining the observed discrepancies, we believe this is mostly the result of a careful analysis aimed at deblending confused sources and removing spurious sources and artifacts from the radio catalogues. This `drop and bump' feature cannot be reproduced by any of the existing state-of-the-art evolutionary models, and appears to be associated with a deficiency of AGN at intermediate redshift ($1<z<2$) and an excess of low-redshift ($z<1$) galaxies and/or AGN. △ Less

Submitted 5 February, 2021; v1 submitted 17 November, 2020; originally announced November 2020.

Comments: Accepted for publication in A&A special issue. This paper is part of the 1st data release of the LoTSS Deep Fields

Journal ref: A&A 648, A5 (2021)

Outflows from Super Star Clusters in the Central Starburst of NGC253

Authors: Rebecca C. Levy, Alberto D. Bolatto, Adam K. Leroy, Kimberly L. Emig, Mark Gorski, Nico Krieger, Laura Lenkic, David S. Meier, Elisabeth A. C. Mills, Juergen Ott, Erik Rosolowsky, Elizabeth Tarantino, Sylvain Veilleux, Fabian Walter, Axel Weiss, Martin A. Zwaan

Abstract: Young massive clusters play an important role in the evolution of their host galaxies, and feedback from the high-mass stars in these clusters can have profound effects on the surrounding interstellar medium. The nuclear starburst in the nearby galaxy NGC253 at a distance of 3.5 Mpc is a key laboratory in which to study star formation in an extreme environment. Previous high resolution (1.9 pc) du… ▽ More Young massive clusters play an important role in the evolution of their host galaxies, and feedback from the high-mass stars in these clusters can have profound effects on the surrounding interstellar medium. The nuclear starburst in the nearby galaxy NGC253 at a distance of 3.5 Mpc is a key laboratory in which to study star formation in an extreme environment. Previous high resolution (1.9 pc) dust continuum observations from ALMA discovered 14 compact, massive super star clusters (SSCs) still in formation. We present here ALMA data at 350 GHz with 28 milliarcsecond (0.5 pc) resolution. We detect blueshifted absorption and redshifted emission (P-Cygni profiles) towards three of these SSCs in multiple lines, including CS 7$-$6 and H$^{13}$CN 4$-$3, which represents direct evidence for previously unobserved outflows. The mass contained in these outflows is a significant fraction of the cluster gas masses, which suggests we are witnessing a short but important phase. Further evidence of this is the finding of a molecular shell around the only SSC visible at near-IR wavelengths. We model the P-Cygni line profiles to constrain the outflow geometry, finding that the outflows must be nearly spherical. Through a comparison of the outflow properties with predictions from simulations, we find that none of the available mechanisms completely explains the observations, although dust-reprocessed radiation pressure and O star stellar winds are the most likely candidates. The observed outflows will have a very substantial effect on the clusters' evolution and star formation efficiency. △ Less

Submitted 6 March, 2021; v1 submitted 10 November, 2020; originally announced November 2020.

Comments: Accepted to ApJ

Super Star Clusters in the Central Starburst of NGC 4945

Authors: Kimberly L. Emig, Alberto D. Bolatto, Adam K. Leroy, Elisabeth A. C. Mills, Maria J. Jimenez Donaire, Alexander G. G. M. Tielens, Adam Ginsburg, Mark Gorski, Nico Krieger, Rebecca C. Levy, David S. Meier, Jurgen Ott, Erik Rosolowsky, Todd A. Thompson, Sylvain Veilleux

Abstract: NGC 4945 is a nearby (3.8 Mpc) galaxy hosting a nuclear starburst and Seyfert Type 2 AGN. We use the Atacama Large Millimeter/submillimeter Array (ALMA) to image the 93 GHz (3.2 mm) free-free continuum and hydrogen recombination line emission (H40$α$ and H42$α$) at 2.2 pc (0.12'') resolution. Our observations reveal 27 bright, compact sources with FWHM sizes of 1.4 - 4.0 pc, which we identify as c… ▽ More NGC 4945 is a nearby (3.8 Mpc) galaxy hosting a nuclear starburst and Seyfert Type 2 AGN. We use the Atacama Large Millimeter/submillimeter Array (ALMA) to image the 93 GHz (3.2 mm) free-free continuum and hydrogen recombination line emission (H40$α$ and H42$α$) at 2.2 pc (0.12'') resolution. Our observations reveal 27 bright, compact sources with FWHM sizes of 1.4 - 4.0 pc, which we identify as candidate super star clusters. Recombination line emission, tracing the ionizing photon rate of the candidate clusters, is detected in 15 sources, 6 of which have a significant synchrotron component to the 93 GHz continuum. Adopting an age of ~5 Myr, the stellar masses implied by the ionizing photon luminosities are $\log_{10}$($M_{\star}$/M$_{\odot}$) $\approx$ 4.7 - 6.1. We fit a slope to the cluster mass distribution and find $β= -1.8 \pm 0.4$. The gas masses associated with these clusters, derived from the dust continuum at 350 GHz, are typically an order of magnitude lower than the stellar mass. These candidate clusters appear to have already converted a large fraction of their dense natal material into stars and, given their small free-fall times of ~0.05 Myr, are surviving an early volatile phase. We identify a point-like source in 93 GHz continuum emission which is presumed to be the AGN. We do not detect recombination line emission from the AGN and place an upper limit on the ionizing photons which leak into the starburst region of $Q_0 < 10^{52}$ s$^{-1}$. △ Less

Submitted 10 September, 2020; originally announced September 2020.

Comments: accepted for publication in ApJ. 27 pages, 11 figures, 6 tables

Searching for the largest bound atoms in space

Authors: K. L. Emig, P. Salas, F. de Gasperin, J. B. R. Oonk, M. C. Toribio, A. P. Mechev, H. J. A. Rottgering, A. G. G. M. Tielens

Abstract: (abridged) Radio recombination lines (RRLs) at frequencies $ν$ < 250 MHz trace the cold, diffuse phase of the ISM. Next generation low frequency interferometers, such as LOFAR, MWA and the future SKA, with unprecedented sensitivity, resolution, and large fractional bandwidths, are enabling the exploration of the extragalactic RRL universe. We observed the radio quasar 3C 190 (z~1.2) with the LOFAR… ▽ More (abridged) Radio recombination lines (RRLs) at frequencies $ν$ < 250 MHz trace the cold, diffuse phase of the ISM. Next generation low frequency interferometers, such as LOFAR, MWA and the future SKA, with unprecedented sensitivity, resolution, and large fractional bandwidths, are enabling the exploration of the extragalactic RRL universe. We observed the radio quasar 3C 190 (z~1.2) with the LOFAR HBA. In reducing this data for spectroscopic analysis, we have placed special emphasis on bandpass calibration. We devised cross-correlation techniques to significantly identify the presence of RRLs in a low frequency spectrum. We demonstrate the utility of this method by applying it to existing low-frequency spectra of Cassiopeia A and M 82, and to the new observations of 3C 190. RRLs have been detected in the foreground of 3C 190 at z = 1.12355 (assuming a carbon origin), owing to the first detection of RRLs outside of the local universe (first reported in Emig et al. 2019). Towards the Galactic supernova remnant Cas A, we uncover three new detections: (1) C$ε$-transitions ($Δ$n = 5) for the first time at low radio frequencies, (2) H$α$-transitions at 64 MHz with a FWHM of 3.1 km/s, the most narrow and one of the lowest frequency detections of hydrogen to date, and (3) C$α$ at v$_{LSR}$ = 0 km/s in the frequency range 55-78 MHz for the first time. Additionally we recover C$α$, C$β$, C$γ$, and C$δ$ from the -47 km/s and -38 km/s components. In the nearby starburst galaxy, M 82, we do not find a significant feature. Our current searches for RRLs in LOFAR observations are limited to narrow (< 100 km/s) features, owing to the relatively small number of channels available for continuum estimation. Future strategies making use of larger contiguous frequency coverage would aid calibration to deeper sensitivities and broader features. △ Less

Submitted 14 October, 2019; originally announced October 2019.

Comments: 21 pages, 21 figures, accepted in A&A

Journal ref: A&A 634, A138 (2020)

Carbon radio recombination lines from gigahertz to megahertz frequencies towards Orion A

Authors: P. Salas, J. B. R. Oonk, K. L. Emig, C. Pabst, M. C. Toribio, H. J. A. Röttgering, A. G. G. M. Tielens

Abstract: Context. The combined use of carbon radio recombination lines (CRRLs) and the 158 $μ$m-[CII] line is a powerful tool for the study of the energetics and physical conditions (e.g., temperature and density) of photodissociation regions (PDRs). However, there are few observational studies that exploit this synergy. Aims. Here we explore the relation between CRRLs and the 158 $μ$m-[CII] line in light… ▽ More Context. The combined use of carbon radio recombination lines (CRRLs) and the 158 $μ$m-[CII] line is a powerful tool for the study of the energetics and physical conditions (e.g., temperature and density) of photodissociation regions (PDRs). However, there are few observational studies that exploit this synergy. Aims. Here we explore the relation between CRRLs and the 158 $μ$m-[CII] line in light of new observations and models. Methods. We present new and existing observations of CRRLs in the frequency range 0.15--230 GHz with ALMA, VLA, the GBT, Effelsberg 100m, and LOFAR towards Orion~A (M42). We complement these observations with SOFIA observations of the 158 $μ$m-[CII] line. We studied two PDRs: the foreground atomic gas, known as the Veil, and the dense PDR between the HII region and the background molecular cloud. Results. In the Veil we are able to determine the gas temperature and electron density, which we use to measure the ionization parameter and the photoelectric heating efficiency. In the dense PDR, we are able to identify a layered PDR structure at the surface of the molecular cloud to the south of the Trapezium cluster. There we find that the radio lines trace the colder portion of the ionized carbon layer, the C$^{+}$/C/CO interface. By modeling the emission of the $158$~$μ$m-[CII] line and CRRLs as arising from a PDR we derive a thermal pressure $>5\times10^{7}$ K cm$^{-3}$ and a radiation field $G_{0}\approx10^{5}$ close to the Trapezium. Conclusions. This work provides additional observational support for the use of CRRLs and the 158 $μ$m-[CII] line as complementary tools to study dense and diffuse PDRs, and highlights the usefulness of CRRLs as probes of the C$^{+}$/C/CO interface. △ Less

Submitted 7 May, 2019; originally announced May 2019.

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

Are starburst galaxies a common source of high energy neutrinos and cosmic rays?

Authors: Cecilia Lunardini, Gregory S. Vance, Kimberly L. Emig, Rogier A. Windhorst

Abstract: A recent analysis of cosmic ray air showers observed at the Pierre Auger Observatory indicates that nearby starburst galaxies (SBGs) might be the cause of ~10% of the Ultra-High-Energy Cosmic Ray flux at energies E > 39 EeV. Since high energy neutrinos are a direct product of cosmic ray interactions, we investigate SBGs as a possible source of some of the ~ 10^-2 - 1 PeV neutrinos observed at IceC… ▽ More A recent analysis of cosmic ray air showers observed at the Pierre Auger Observatory indicates that nearby starburst galaxies (SBGs) might be the cause of ~10% of the Ultra-High-Energy Cosmic Ray flux at energies E > 39 EeV. Since high energy neutrinos are a direct product of cosmic ray interactions, we investigate SBGs as a possible source of some of the ~ 10^-2 - 1 PeV neutrinos observed at IceCube. A statistical analysis is performed to establish the degree of positional correlation between the observed neutrinos and a set of 45 nearby radio- and infrared-bright SBGs. Our results are consistent with no causal correlation. However, a scenario where ~10% of the High Energy Starting Events (HESE) in the detector are coming from the candidate SBGs is not excluded. The same conclusion is reached for different data subsets, as well as two different subsets of SBGs motivated by the Pierre Auger Observatory analysis. △ Less

Submitted 3 November, 2019; v1 submitted 25 February, 2019; originally announced February 2019.

Comments: 13 pages, including 1 figure and 1 table. Expanded statistical analysis, expanded text. Conclusions unchanged

The first detection of radio recombination lines at cosmological distances

Authors: K. L. Emig, P. Salas, F. de Gasperin, J. B. R. Oonk, M. C. Toribio, H. J. A. Rottgering, A. G. G. M. Tielens

Abstract: Recombination lines involving high principal quantum numbers populate the radio spectrum in large numbers. Low-frequency (<1 GHz) observations of radio recombination lines (RRLs) primarily from carbon and hydrogen offer a new, if not unique, way to probe cold, largely atomic gas and warm, ionised gas in other galaxies. Furthermore, RRLs can be used to determine the physical state of the emitting r… ▽ More Recombination lines involving high principal quantum numbers populate the radio spectrum in large numbers. Low-frequency (<1 GHz) observations of radio recombination lines (RRLs) primarily from carbon and hydrogen offer a new, if not unique, way to probe cold, largely atomic gas and warm, ionised gas in other galaxies. Furthermore, RRLs can be used to determine the physical state of the emitting regions, such as temperature and density. These properties make RRLs, potentially, a powerful tool of extragalactic ISM physics. At low radio frequencies, its conceivable to detect RRLs out to cosmological distances when illuminated by a strong radio continuum. However, they are extremely faint (tau ~ 1e-3 -- 1e-4) and have so far eluded detection outside of the local universe. With LOFAR observations of the radio quasar 3C 190 (z=1.1946), we aim to demonstrate that the ISM can be explored out to great distances through low-frequency RRLs. We report the detection of RRLs in the frequency range 112--163 MHz in the spectrum of 3C 190. Stacking 13 a-transitions with principal quantum numbers n=266-301, a peak 6sigma feature of optical depth, tau(peak) = (1.0 +\- 0.2) x 1e-3 and FWHM = 31.2 +/- 8.3 km/s was found at z=1.124. This corresponds to a velocity offset of -9965 km/s with respect to the systemic redshift of 3C 190. We consider three interpretations of the origin of the RRL emission: an intervening dwarf-like galaxy, an AGN-driven outflow, and the inter-galactic medium. We argue that the RRLs most likely originate in a dwarf-like galaxy (M ~ 1e9 Msun) along the line of sight, although we cannot rule out an AGN-driven outflow. We do find the RRLs to be inconsistent with an inter-galactic medium origin. With this detection, we have opened up a new way to study the physical properties of cool, diffuse gas out to cosmological distances. △ Less

Submitted 20 November, 2018; originally announced November 2018.

Comments: 5 figures, 2 tables and 9 pages. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available on lofar-surveys.org upon publication of the journal

Journal ref: A&A 622, A7 (2019)

Systematic effects in LOFAR data: A unified calibration strategy

Authors: F. de Gasperin, T. J. Dijkema, A. Drabent, M. Mevius, D. Rafferty, R. van Weeren, M. Brüggen, J. R. Callingham, K. L. Emig, G. Heald, H. T. Intema, L. K. Morabito, A. R. Offringa, R. Oonk, E. Orrù, H. Röttgering, J. Sabater, T. Shimwell, A. Shulevski, W. Williams

Abstract: Context: New generation low-frequency telescopes are exploring a new parameter space in terms of depth and resolution. The data taken with these interferometers, for example with the LOw Frequency ARray (LOFAR), are often calibrated in a low signal-to-noise ratio regime and the removal of critical systematic effects is challenging. The process requires an understanding of their origin and properti… ▽ More Context: New generation low-frequency telescopes are exploring a new parameter space in terms of depth and resolution. The data taken with these interferometers, for example with the LOw Frequency ARray (LOFAR), are often calibrated in a low signal-to-noise ratio regime and the removal of critical systematic effects is challenging. The process requires an understanding of their origin and properties. Aim: In this paper we describe the major systematic effects inherent to next generation low-frequency telescopes, such as LOFAR. With this knowledge, we introduce a data processing pipeline that is able to isolate and correct these systematic effects. The pipeline will be used to calibrate calibrator observations as the first step of a full data reduction process. Methods: We processed two LOFAR observations of the calibrator 3C196: the first using the Low Band Antenna (LBA) system at 42-66 MHz and the second using the High Band Antenna (HBA) system at 115-189 MHz. Results: We were able to isolate and correct for the effects of clock drift, polarisation misalignment, ionospheric delay, Faraday rotation, ionospheric scintillation, beam shape, and bandpass. The designed calibration strategy produced the deepest image to date at 54 MHz. The image has been used to confirm that the spectral energy distribution of the average radio source population tends to flatten at low frequencies. Conclusions: We prove that LOFAR systematic effects can be described by a relatively small number of parameters. Furthermore, the identification of these parameters is fundamental to reducing the degrees of freedom when the calibration is carried out on fields that are not dominated by a strong calibrator. △ Less

Submitted 19 November, 2018; originally announced November 2018.

Comments: 20 pages, 15 figures. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available on lofar-surveys.org upon publication of the journal

The LOFAR Two-metre Sky Survey - II. First data release

Authors: T. W. Shimwell, C. Tasse, M. J. Hardcastle, A. P. Mechev, W. L. Williams, P. N. Best, H. J. A. Röttgering, J. R. Callingham, T. J. Dijkema, F. de Gasperin, D. N. Hoang, B. Hugo, M. Mirmont, J. B. R. Oonk, I. Prandoni, D. Rafferty, J. Sabater, O. Smirnov, R. J. van Weeren, G. J. White, M. Atemkeng, L. Bester, E. Bonnassieux, M. Brüggen, G. Brunetti , et al. (82 additional authors not shown)

Abstract: The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45… ▽ More The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45$^\circ$00$'$00$''$ to 57$^\circ$00$'$00$''$) were mapped using a fully automated direction-dependent calibration and imaging pipeline that we developed. A total of 325,694 sources are detected with a signal of at least five times the noise, and the source density is a factor of $\sim 10$ higher than the most sensitive existing very wide-area radio-continuum surveys. The median sensitivity is S$_{\rm 144 MHz} = 71\,μ$Jy beam$^{-1}$ and the point-source completeness is 90% at an integrated flux density of 0.45mJy. The resolution of the images is 6$''$ and the positional accuracy is within 0.2$''$. This data release consists of a catalogue containing location, flux, and shape estimates together with 58 mosaic images that cover the catalogued area. In this paper we provide an overview of the data release with a focus on the processing of the LOFAR data and the characteristics of the resulting images. In two accompanying papers we provide the radio source associations and deblending and, where possible, the optical identifications of the radio sources together with the photometric redshifts and properties of the host galaxies. These data release papers are published together with a further $\sim$20 articles that highlight the scientific potential of LoTSS. △ Less

Submitted 19 November, 2018; originally announced November 2018.

Comments: 16 figures, 1 table and 22 pages. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available on lofar-surveys.org upon publication of the journal

Mapping low frequency carbon radio recombination lines towards Cassiopeia A at 340, 148, 54 and 43 MHz

Authors: Pedro Salas, J. B. Raymond Oonk, Reinout J. van Weeren, Mark G. Wolfire, Kimberly L. Emig, M. Carmen Toribio, Huub J. A. Röttgering, Alexander G. G. M. Tielens

Abstract: Quantitative understanding of the interstellar medium requires knowledge of its physical conditions. Low frequency carbon radio recombination lines (CRRLs) trace cold interstellar gas, and can be used to determine its physical conditions (e.g., electron temperature and density). In this work we present spatially resolved observations of the low frequency ($\leq390$ MHz) CRRLs centered around C… ▽ More Quantitative understanding of the interstellar medium requires knowledge of its physical conditions. Low frequency carbon radio recombination lines (CRRLs) trace cold interstellar gas, and can be used to determine its physical conditions (e.g., electron temperature and density). In this work we present spatially resolved observations of the low frequency ($\leq390$ MHz) CRRLs centered around C$268α$, C$357α$, C$494α$ and C$539α$ towards Cassiopeia A on scales of $\leq1.2$ pc. We compare the spatial distribution of CRRLs with other ISM tracers. This comparison reveals a spatial offset between the peak of the CRRLs and other tracers, which is very characteristic for photodissociation regions and that we take as evidence for CRRLs being preferentially detected from the surfaces of molecular clouds. Using the CRRLs we constrain the gas electron temperature and density. These constraints on the gas conditions suggest variations of less than a factor of two in pressure over $\sim1$ pc scales, and an average hydrogen density of $200$-$470$ cm$^{-3}$. From the electron temperature and density maps we also constrain the ionized carbon emission measure, column density and path length. Based on these, the hydrogen column density is larger than $10^{22}$ cm$^{-2}$, with a peak of $\sim4\times10^{22}$ cm$^{-2}$ towards the South of Cassiopeia A. Towards the southern peak the line of sight length is $\sim40$ pc over a $\sim2$ pc wide structure, which implies that the gas is a thin surface layer on a large (molecular) cloud that is only partially intersected by Cassiopeia A. These observations highlight the utility of CRRLs as tracers of low density extended HI and CO-dark gas halo's around molecular clouds. △ Less

Submitted 16 January, 2018; originally announced January 2018.

Comments: 18 pages, 15 figures, 4 tables, accepted by MNRAS

LOFAR observations of decameter carbon radio recombination lines towards Cassiopeia A

Authors: P. Salas, J. B. R. Oonk, R. J. van Weeren, F. Salgado, L. K. Morabito, M. C. Toribio, K. Emig, H. J. A. Röttgering, A. G. G. M. Tielens

Abstract: We present a study of carbon radio recombination lines towards Cassiopeia A using LOFAR observations in the frequency range 10-33 MHz. Individual carbon $α$ lines are detected in absorption against the continuum at frequencies as low as 16 MHz. Stacking several C$α$ lines we obtain detections in the 11-16 MHz range. These are the highest signal-to-noise measurements at these frequencies. The peak… ▽ More We present a study of carbon radio recombination lines towards Cassiopeia A using LOFAR observations in the frequency range 10-33 MHz. Individual carbon $α$ lines are detected in absorption against the continuum at frequencies as low as 16 MHz. Stacking several C$α$ lines we obtain detections in the 11-16 MHz range. These are the highest signal-to-noise measurements at these frequencies. The peak optical depth of the C$α$ lines changes considerably over the 11-33 MHz range with the peak optical depth decreasing from 4$\times10^{-3}$ at 33 MHz to 2$\times10^{-3}$ at 11 MHz, while the line width increases from 20 km s$^{-1}$ to 150 km s$^{-1}$. The combined change in peak optical depth and line width results in a roughly constant integrated optical depth. We interpret this as carbon atoms close to local thermodynamic equilibrium. In this work we focus on how the 11-33 MHz carbon radio recombination lines can be used to determine the gas physical conditions. We find that the ratio of the carbon radio recombination lines to that of the 158 $μ$m [CII] fine-structure line is a good thermometer, while the ratio between low frequency carbon radio recombination lines provides a good barometer. By combining the temperature and pressure constraints with those derived from the line width we are able to constrain the gas properties (electron temperature and density) and radiation field intensity. Given the 1$σ$ uncertainties in our measurements these are; $T_{e}\approx68$-$98$ K, $n_{e}\approx0.02$-$0.035$ cm$^{-3}$ and $T_{r,100}\approx1500$-$1650$ K. Despite challenging RFI and ionospheric conditions, our work demonstrates that observations of carbon radio recombination lines in the 10-33 MHz range can provide insight into the gas conditions. △ Less

Submitted 30 January, 2017; originally announced January 2017.

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

Extra-galactic background light measurements from the far-UV to the far-IR from deep ground and space-based galaxy counts

Authors: Simon P. Driver, Stephen K. Andrews, Luke J. Davies, Aaron S. G. Robotham, Angus H. Wright, Rogier A. Windhorst, Seth Cohen, Kim Emig, Rolf A. Jansen, Loretta Dunne

Abstract: We combine wide and deep galaxy number-count data from GAMA, COSMOS/G10, HST ERS, HST UVUDF and various near-, mid- and far- IR datasets from ESO, Spitzer and Herschel. The combined data range from the far-UV (0.15microns) to far-IR (500microns), and in all cases the contribution to the integrated galaxy light (IGL) of successively fainter galaxies converges. Using a simple spline fit, we derive t… ▽ More We combine wide and deep galaxy number-count data from GAMA, COSMOS/G10, HST ERS, HST UVUDF and various near-, mid- and far- IR datasets from ESO, Spitzer and Herschel. The combined data range from the far-UV (0.15microns) to far-IR (500microns), and in all cases the contribution to the integrated galaxy light (IGL) of successively fainter galaxies converges. Using a simple spline fit, we derive the IGL and the extrapolated-IGL in all bands. We argue undetected low surface brightness galaxies and intra-cluster/group light is modest, and that our extrapolated-IGL measurements are an accurate representation of the extra-galactic background light. Our data agree with most earlier IGL estimates and with direct measurements in the far-IR, but disagree strongly with direct estimates in the optical. Close agreement between our results and recent very high-energy experiments (H.E.S.S. and MAGIC), suggest that there may be an additional foreground affecting the direct estimates. The most likely culprit could be the adopted Zodiacal light model. Finally we use a modified version of the two-component model to integrate the EBL and obtain measurements of the Cosmic Optical Background (COB) and Cosmic Infrared Background (CIB) of: $24^{+4}_{-4}$nWm$^{-2}$sr$^{-1}$ and $26^{+5}_{-5}$nWm$^{-2}$sr$^{-1}$ respectively (48:52%). Over the next decade, upcoming space missions such as Euclid and WFIRST, have the capacity to reduce the COB error to $<1\%$, at which point comparisons to the very high energy data could have the potential to provide a direct detection and measurement of the reionisation field. △ Less

Submitted 5 May, 2016; originally announced May 2016.

Comments: Accepted for publication in The Astrophysical Journal

Do high energy astrophysical neutrinos trace star formation?

Authors: Kimberly Emig, Cecilia Lunardini, Rogier Windhorst

Abstract: The IceCube Neutrino Observatory has provided the first map of the high energy (~ 0.01 -- 1 PeV) sky in neutrinos. Since neutrinos propagate undeflected, their arrival direction is an important identifier for sources of high energy particle acceleration. Reconstructed arrival directions are consistent with an extragalactic origin, with possibly a galactic component, of the neutrino flux. We presen… ▽ More The IceCube Neutrino Observatory has provided the first map of the high energy (~ 0.01 -- 1 PeV) sky in neutrinos. Since neutrinos propagate undeflected, their arrival direction is an important identifier for sources of high energy particle acceleration. Reconstructed arrival directions are consistent with an extragalactic origin, with possibly a galactic component, of the neutrino flux. We present a statistical analysis of positional coincidences of the IceCube neutrinos with known astrophysical objects from several catalogs. When considering starburst galaxies with the highest flux in gamma-rays and infrared radiation, up to n=8 coincidences are found, representing an excess over the ~4 predicted for the randomized, or "null" distribution. The probability that this excess is realized in the null case, the p-value, is p=0.042. This value falls to p=0.003 for a partial subset of gamma-ray-detected starburst galaxies and superbubble regions in the galactic neighborhood. Therefore, it is possible that starburst galaxies, and the typically hundreds of superbubble regions within them, might account for a portion of IceCube neutrinos. The physical plausibility of such correlation is discussed briefly. △ Less

Submitted 15 December, 2015; v1 submitted 21 July, 2015; originally announced July 2015.

Comments: 28 pages, 5 figures, 2 tables, accepted in JCAP

Report number: INT-PUB-15-030