Showing 1–50 of 127 results for author: Conway, J

LOFAR Deep Fields: Probing the sub-mJy regime of polarized extragalactic sources in ELAIS-N1. I. The catalog

Authors: S. Piras, C. Horellou, J. E. Conway, M. Thomasson, S. del Palacio, T. W. Shimwell, S. P. O'Sullivan, E. Carretti, I. Šnidaric, V. Jelic, B. Adebahr, A. Berger, P. N. Best, M. Brüggen, N. Herrera Ruiz, R. Paladino, I. Prandoni, J. Sabater, V. Vacca

Abstract: The aim of this study is to probe the sub-mJy polarized source population with LOFAR. We present the method used to stack LOFAR polarization datasets, the resulting catalog of polarized sources, and the derived polarized source counts. The ELAIS-N1 field was selected for a polarimetric study at 114.9-177.4 MHz. A total area of 25 deg2 was imaged at 6"- resolution in the Stokes Q and U parameters.… ▽ More The aim of this study is to probe the sub-mJy polarized source population with LOFAR. We present the method used to stack LOFAR polarization datasets, the resulting catalog of polarized sources, and the derived polarized source counts. The ELAIS-N1 field was selected for a polarimetric study at 114.9-177.4 MHz. A total area of 25 deg2 was imaged at 6"- resolution in the Stokes Q and U parameters. Alignment of polarization angles was done both in frequency and in Faraday space before stacking datasets from 19 eight-hour-long epochs. A search for polarized sources was carried out in the final, stacked dataset, and the properties of the detected sources were examined. The depolarization level of sources known to be polarized at 1.4 GHz was quantified. A one-sigma noise level of 19 μJy/beam was reached in the central part of the field after stacking. Twenty-five polarized sources were detected above 8σ, five of which had not been detected in polarization at any other radio frequencies before. Seven additional polarized components were found by lowering the threshold to 6σat positions corresponding to sources known to be polarized at 1.4 GHz. In two radio galaxies, polarization was detected from both radio lobes, so the final number of associated radio continuum sources is 31. The detected sources are weakly polarized, with a median degree of polarization of 1.75% for the sample of sources detected in polarized emission. The sources previously detected in polarization at 1.4 GHz are significantly depolarized at 150 MHz. The catalog is used to derive the polarized source counts at 150 MHz. This is the deepest and highest-resolution polarization study at 150 MHz to date. △ Less

Submitted 12 June, 2024; originally announced June 2024.

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

The PARADIGM Project I: A Multiscale radio morphological analysis of local U/LIRGS

Authors: G. Lucatelli, R. Beswick, J. Moldon, M. Á. Pérez-Torres, J. E. Conway, A. Alberdi, C. Romero-Cañizales, E. Varenius, H. -R. Klöckner, L. Barcos-Muñoz, M. Bondi, S. T. Garrington, S. Aalto, W. A. Baan, Y. M. Pihlstrom

Abstract: Disentangling the radio flux contribution from star formation (SF) and active-galactic-nuclei (AGN) activity is a long-standing problem in extragalactic astronomy, since at frequencies of $\lesssim$ 10 GHz, both processes emit synchrotron radiation. We present in this work the general objectives of the PARADIGM Project, a multi-instrument concept to explore star-formation and mass assembly of gala… ▽ More Disentangling the radio flux contribution from star formation (SF) and active-galactic-nuclei (AGN) activity is a long-standing problem in extragalactic astronomy, since at frequencies of $\lesssim$ 10 GHz, both processes emit synchrotron radiation. We present in this work the general objectives of the PARADIGM Project, a multi-instrument concept to explore star-formation and mass assembly of galaxies. We introduce two novel general approaches for a detailed multiscale study of the radio emission in local U/LIRGs. In this work, we use archival interferometric data from the Very Large Array (VLA) centred at ~ 6 GHz (C band) and present new observations from the e-Multi-Element Radio-Linked Interferometer Network (e-MERLIN) for UGC5101, VV705, VV250 and UGC8696. Using our image decomposition methods, we robustly disentangle the radio emission into distinct components by combining information from the two interferometric arrays. We use e-MERLIN as a probe of the core-compact radio emission (AGN or starburst) at ~ 20 pc scales, and as a probe of nuclear diffuse emission, at scales ~ 100 - 200 pc. With VLA, we characterise the source morphology and the flux density on scales from 200 pc up to and above 1 kpc. As a result, we find deconvolved and convolved sizes for nuclear regions from ~ 10 pc to ~ 200 pc. At larger scales, we find sizes of 1.5 - 2 kpc for diffuse structures (with effective sizes of ~ 300 - 400 pc). We demonstrate that the radio emission from nuclear extended structures (~ 100 pc) can dominate over core-compact components, providing a significant fraction of the total multiscale SF output. We establish a multiscale radio tracer for star formation by combining information from different instruments. Consequently, this work sets a starting point to potentially correct for overestimations of AGN fractions and underestimates of SF activity. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: Accepted for publication in MNRAS. 31 pages (9 pages appendix), 13 figures, 9 tables. Online code repository: morphen - https://github.com/lucatelli/morphen

Ordered magnetic fields around the 3C 84 central black hole

Authors: G. F. Paraschos, J. -Y. Kim, M. Wielgus, J. Röder, T. P. Krichbaum, E. Ros, I. Agudo, I. Myserlis, M. Moscibrodzka, E. Traianou, J. A. Zensus, L. Blackburn, C. -K. Chan, S. Issaoun, M. Janssen, M. D. Johnson, V. L. Fish, K. Akiyama, A. Alberdi, W. Alef, J. C. Algaba, R. Anantua, K. Asada, R. Azulay, U. Bach , et al. (258 additional authors not shown)

Abstract: 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures a… ▽ More 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures and understand the physical conditions in the compact region of 3C84. We used EHT 228GHz observations and, given the limited (u,v)-coverage, applied geometric model fitting to the data. We also employed quasi-simultaneously observed, multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure. We report the detection of a highly ordered, strong magnetic field around the central, SMBH of 3C84. The brightness temperature analysis suggests that the system is in equipartition. We determined a turnover frequency of $ν_m=(113\pm4)$GHz, a corresponding synchrotron self-absorbed magnetic field of $B_{SSA}=(2.9\pm1.6)$G, and an equipartition magnetic field of $B_{eq}=(5.2\pm0.6)$G. Three components are resolved with the highest fractional polarisation detected for this object ($m_\textrm{net}=(17.0\pm3.9)$%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017-2018. We report a steeply negative slope of the spectrum at 228GHz. We used these findings to test models of jet formation, propagation, and Faraday rotation in 3C84. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C84. However, systematic uncertainties due to the limited (u,v)-coverage, however, cannot be ignored. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 15 pages, 6 figures, published in A&A

Journal ref: Issue: A&A Volume 682, February 2024; Article number: L3; Number of pages: 15

A search for pulsars around Sgr A* in the first Event Horizon Telescope dataset

Authors: Pablo Torne, Kuo Liu, Ralph P. Eatough, Jompoj Wongphechauxsorn, James M. Cordes, Gregory Desvignes, Mariafelicia De Laurentis, Michael Kramer, Scott M. Ransom, Shami Chatterjee, Robert Wharton, Ramesh Karuppusamy, Lindy Blackburn, Michael Janssen, Chi-kwan Chan, Geoffrey B. Crew, Lynn D. Matthews, Ciriaco Goddi, Helge Rottmann, Jan Wagner, Salvador Sanchez, Ignacio Ruiz, Federico Abbate, Geoffrey C. Bower, Juan J. Salamanca , et al. (261 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission… ▽ More The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission spectra - are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic Center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most-sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the Fast-Folding-Algorithm and single pulse search targeting both pulsars and burst-like transient emission; using the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction ($\lesssim$2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

Sub-arcsecond resolution imaging of M 51 with the International LOFAR Telescope

Authors: Deepika Venkattu, Peter Lundqvist, Miguel Pérez-Torres, Leah Morabito, Javier Moldón, John Conway, Poonam Chandra, Cyril Tasse

Abstract: We present an International LOFAR Telescope sub-arcsecond resolution image of the nearby galaxy M 51 with a beam size of 0.436" x 0.366" and rms of 46 $μ$Jy. We compare this image with an European VLBI Network study of M 51, and discuss the supernovae in this galaxy, which have not yet been probed at these low radio frequencies. We find a flux density of 0.97 mJy for SN 2011dh in the ILT image, wh… ▽ More We present an International LOFAR Telescope sub-arcsecond resolution image of the nearby galaxy M 51 with a beam size of 0.436" x 0.366" and rms of 46 $μ$Jy. We compare this image with an European VLBI Network study of M 51, and discuss the supernovae in this galaxy, which have not yet been probed at these low radio frequencies. We find a flux density of 0.97 mJy for SN 2011dh in the ILT image, which is about five times smaller than the flux density reported by the LOFAR Twometre Sky Survey at 6" resolution using the same dataset without the international stations. This difference makes evident the need for LOFAR international baselines to reliably obtain flux density measurements of compact objects in nearby galaxies. Our LOFAR flux density measurement of SN 2011dh directly translates into fitting the radio light curves for the supernova and constraining massloss rates of progenitor star. We do not detect two other supernovae in the same galaxy, SN 1994I and SN 2005cs, and our observations place limits on the evolution of both supernovae at radio wavelengths. We also discuss the radio emission from the centre of M 51, in which we detect the Active Galactic Nucleus and other parts of the nuclear emission in the galaxy, and a possible detection of Component N. We discuss a few other sources, including the detection of a High mass X-ray Binary not detected by LoTSS, but with a flux density in the ILT image that matches well with higher frequency catalogues. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: Accepted for publication in the Astrophysical Journal. 16 pages, 5 figures

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

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

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

Submitted 25 April, 2023; originally announced April 2023.

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

Comparison of Polarized Radiative Transfer Codes used by the EHT Collaboration

Authors: Ben S. Prather, Jason Dexter, Monika Moscibrodzka, Hung-Yi Pu, Thomas Bronzwaer, Jordy Davelaar, Ziri Younsi, Charles F. Gammie, Roman Gold, George N. Wong, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, Anne-Kathrin Baczko, David Ball, Mislav Baloković, John Barrett, Michi Bauböck, Bradford A. Benson, Dan Bintley , et al. (248 additional authors not shown)

Abstract: Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curve… ▽ More Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U , and V respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ

SKA Science Data Challenge 2: analysis and results

Authors: P. Hartley, A. Bonaldi, R. Braun, J. N. H. S. Aditya, S. Aicardi, L. Alegre, A. Chakraborty, X. Chen, S. Choudhuri, A. O. Clarke, J. Coles, J. S. Collinson, D. Cornu, L. Darriba, M. Delli Veneri, J. Forbrich, B. Fraga, A. Galan, J. Garrido, F. Gubanov, H. Håkansson, M. J. Hardcastle, C. Heneka, D. Herranz, K. M. Hess , et al. (83 additional authors not shown)

Abstract: The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed t… ▽ More The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed to familiarise the scientific community with SKAO data and to drive the development of new analysis techniques. We present the results from Science Data Challenge 2 (SDC2), which invited participants to find and characterise 233245 neutral hydrogen (Hi) sources in a simulated data product representing a 2000~h SKA MID spectral line observation from redshifts 0.25 to 0.5. Through the generous support of eight international supercomputing facilities, participants were able to undertake the Challenge using dedicated computational resources. Alongside the main challenge, `reproducibility awards' were made in recognition of those pipelines which demonstrated Open Science best practice. The Challenge saw over 100 participants develop a range of new and existing techniques, with results that highlight the strengths of multidisciplinary and collaborative effort. The winning strategy -- which combined predictions from two independent machine learning techniques to yield a 20 percent improvement in overall performance -- underscores one of the main Challenge outcomes: that of method complementarity. It is likely that the combination of methods in a so-called ensemble approach will be key to exploiting very large astronomical datasets. △ Less

Submitted 14 March, 2023; originally announced March 2023.

Comments: Under review by MNRAS; 28 pages, 16 figures

Intermediate-mass black holes: finding of episodic, large-scale and powerful jet activity in a dwarf galaxy

Authors: Jun Yang, Zsolt Paragi, Sándor Frey, Leonid I. Gurvits, Mai Liao, Xiang Liu, Lang Cui, Xiaolong Yang, Wen Chen, Eskil Varenius, John E. Conway, Rurong Chen, Ning Chang

Abstract: Dwarf galaxies are characterised by a very low luminosity and low mass. Because of significant accretion and ejection activity of massive black holes, some dwarf galaxies also host low-luminosity active galactic nuclei (AGNs). In a few dwarf AGNs, very long baseline interferometry (VLBI) observations have found faint non-thermal radio emission. SDSS J090613.77+561015.2 is a dwarf AGN owning an int… ▽ More Dwarf galaxies are characterised by a very low luminosity and low mass. Because of significant accretion and ejection activity of massive black holes, some dwarf galaxies also host low-luminosity active galactic nuclei (AGNs). In a few dwarf AGNs, very long baseline interferometry (VLBI) observations have found faint non-thermal radio emission. SDSS J090613.77+561015.2 is a dwarf AGN owning an intermediate-mass black hole (IMBH) with a mass of $M_{BH} = 3.6^{+5.9}_{-2.3} \times 10^5 M_{sun}$ and showing a rarely-seen two-component radio structure in its radio nucleus. To further probe their nature, i.e. the IMBH jet activity, we performed additional deep observations with the European VLBI Network (EVN) at 1.66 GHz and 4.99 GHz. We find the more diffuse emission regions and structure details. These new EVN imaging results allow us to reveal a two-sided jet morphology with a size up to about 150 mas (projected length $\sim$140 pc) and a radio luminosity of about $3\times10^{38}$ erg s$^{-1}$. The peak feature has an optically thin radio spectrum and thus more likely represents a relatively young ejecta instead of a jet base. The EVN study on SDSS J090613.77+561015.2 demonstrates the existence of episodic, relatively large-scale and powerful IMBH jet activity in dwarf AGNs. Moreover, we collected a small sample of VLBI-detected dwarf AGNs and investigated their connections with normal AGNs. We notice that these radio sources in the dwarf AGNs tend to have steep spectra and small linear sizes, and possibly represent ejecta from scaled-down episodic jet activity. △ Less

Submitted 13 February, 2023; originally announced February 2023.

Comments: Accepted for publication in MRNAS, 10 pages, 6 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Vol. 520, pp. 5964-5973 (2023)

The accuracy of mutual potential approximations in simulations of binary asteroids

Authors: Alex Ho, Margrethe Wold, Mohammad Poursina, John T. Conway

Abstract: Simulations of asteroid binaries commonly use mutual gravitational potentials approximated by series expansions, leading to truncation errors, and also preventing correct computations of the forces and torques when the bodies are close. We make of a recently developed method where the mutual potential is calculated with the use of surface integrals and is exact for bodies of ellipsoidal shapes. Th… ▽ More Simulations of asteroid binaries commonly use mutual gravitational potentials approximated by series expansions, leading to truncation errors, and also preventing correct computations of the forces and torques when the bodies are close. We make of a recently developed method where the mutual potential is calculated with the use of surface integrals and is exact for bodies of ellipsoidal shapes. The solutions produced by the surface integration method are compared with an approach that expands the mutual potential, truncated at second and fourth order. The approximate solutions are generated with the ``General Use Binary Asteroid Simulator'' (gubas). We find that the differences in the forces and torques are the largest when the bodies are nearly touching. These differences can exceed 1000% if the shape of the primary is highly elongated. Long term simulations show more than 100% difference in the dynamics if the bodies are initially close, while the differences are negligible if the bodies are initially far apart. For simulations with two triaxial ellipsoids, the computational efficiency of the surface integral method is comparable to fourth order approximations with gubas, and superior to potentials truncated to order eight or higher. △ Less

Submitted 24 January, 2023; originally announced January 2023.

Comments: 15 pages, 16 figures. Accepted for publication in A&A. The abstract has been abridged for arXiv submission

Journal ref: A&A 671, A38 (2023)

Identifying active galactic nuclei via brightness temperature with sub-arcsecond International LOFAR Telescope observations

Authors: Leah K. Morabito, F. Sweijen, J. F. Radcliffe, P. N. Best, Rohit Kondapally, Marco Bondi, Matteo Bonato, K. J. Duncan, Isabella Prandoni, T. W. Shimwell, W. L. Williams, R. J. van Weeren, J. E. Conway, G. Calistro Rivera

Abstract: Identifying active galactic nuclei (AGN) and isolating their contribution to a galaxy's energy budget is crucial for studying the co-evolution of AGN and their host galaxies. Brightness temperature ($T_b$) measurements from high-resolution radio observations at GHz frequencies are widely used to identify AGN. Here we investigate using new sub-arcsecond imaging at 144 MHz with the International LOF… ▽ More Identifying active galactic nuclei (AGN) and isolating their contribution to a galaxy's energy budget is crucial for studying the co-evolution of AGN and their host galaxies. Brightness temperature ($T_b$) measurements from high-resolution radio observations at GHz frequencies are widely used to identify AGN. Here we investigate using new sub-arcsecond imaging at 144 MHz with the International LOFAR Telescope to identify AGN using $T_b$ in the Lockman Hole field. We use ancillary data to validate the 940 AGN identifications, finding 83 percent of sources have AGN classifications from SED fitting and/or photometric identifications, yielding 160 new AGN identifications. Considering the multi-wavelength classifications, brightness temperature criteria select over half of radio-excess sources, 32 percent of sources classified as radio-quiet AGN, and 20 percent of sources classified as star-forming galaxies. Infrared colour-colour plots and comparison with what we would expect to detect based on peak brightness in 6 arcsec LOFAR maps, imply that the star-forming galaxies and sources at low flux densities have a mixture of star-formation and AGN activity. We separate the radio emission from star-formation and AGN in unresolved, $T_b$-identified AGN with no significant radio excess and find the AGN comprises $0.49\pm 0.16$ of the radio luminosity. Overall the non-radio excess AGN show evidence for having a variety of different radio emission mechanisms, which can provide different pathways for AGN and galaxy co-evolution. This validation of AGN identification using brightness temperature at low frequencies opens the possibility for securely selecting AGN samples where ancillary data is inadequate. △ Less

Submitted 26 July, 2022; originally announced July 2022.

Comments: Accepted for publication in MNRAS. 19 pages, 15 figures, 2 tables

Dynamics of asteroid systems post rotational fission

Authors: Alex Ho, Margrethe Wold, Mohammad Poursina, John T. Conway

Abstract: Asteroid binaries found amongst the Near-Earth objects are believed to have formed from rotational fission. In this paper, we aim to study the dynamical evolution of asteroid systems the moment after fission. The initial condition is modelled as a contact binary, similar to that of Boldrin et al. (2016). Both bodies are modelled as ellipsoids, and the secondary is given an initial rotation angle a… ▽ More Asteroid binaries found amongst the Near-Earth objects are believed to have formed from rotational fission. In this paper, we aim to study the dynamical evolution of asteroid systems the moment after fission. The initial condition is modelled as a contact binary, similar to that of Boldrin et al. (2016). Both bodies are modelled as ellipsoids, and the secondary is given an initial rotation angle about its body-fixed $y$-axis. Moreover, we consider six different cases, three where the density of the secondary varies, and three where we vary its shape. The simulations consider 45 different initial tilt angles of the secondary, each with 37 different mass ratios. We start the dynamical simulations at the moment the contact binary reaches a spin fission limit, and our model ensures that the closest distance between the surfaces of the two bodies is always kept at 1 cm. The forces, torques and gravitational potential between the two bodies are modelled using a newly developed surface integration scheme, giving exact results for two ellipsoids. We find that more than 80% of the simulations end with the two bodies impacting, and collisions between the bodies are more common when the density of the secondary is lower, or when it becomes more elongated. When comparing with data on asteroid pairs from Pravec et al. (2019) we find that variations in density and shape of the secondary can account for some of the spread seen in the rotation period for observed pairs. Furthermore, the secondary may also reach a spin limit for surface disruption, creating a ternary/multiple system. We find that secondary fission typically occurs within the first five hours after the contact binary separates, and is more common when the secondary is less dense or more elongated. △ Less

Submitted 23 June, 2022; originally announced June 2022.

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

Journal ref: A&A 665, A43 (2022)

The science case and challenges of space-borne sub-millimeter interferometry

Authors: Leonid I. Gurvits, Zsolt Paragi, Ricardo I. Amils, Ilse van Bemmel, Paul Boven, Viviana Casasola, John Conway, Jordy Davelaar, M. Carmen Díez-González, Heino Falcke, Rob Fender, Sándor Frey, Christian M. Fromm, Juan D. Gallego-Puyol, Cristina García-Miró, Michael A. Garrett, Marcello Giroletti, Ciriaco Goddi, José L. Gómez, Jeffrey van der Gucht, José Carlos Guirado, Zoltán Haiman, Frank Helmich, Ben Hudson, Elizabeth Humphreys , et al. (29 additional authors not shown)

Abstract: Ultra-high angular resolution in astronomy has always been an important vehicle for making fundamental discoveries. Recent results in direct imaging of the vicinity of the supermassive black hole in the nucleus of the radio galaxy M87 by the millimeter VLBI system Event Horizon Telescope and various pioneering results of the Space VLBI mission RadioAstron provided new momentum in high angular reso… ▽ More Ultra-high angular resolution in astronomy has always been an important vehicle for making fundamental discoveries. Recent results in direct imaging of the vicinity of the supermassive black hole in the nucleus of the radio galaxy M87 by the millimeter VLBI system Event Horizon Telescope and various pioneering results of the Space VLBI mission RadioAstron provided new momentum in high angular resolution astrophysics. In both mentioned cases, the angular resolution reached the values of about 10-20 microrcseconds. Further developments toward at least an order of magnitude "sharper" values are dictated by the needs of astrophysical studies and can only be achieved by placing millimeter and submillimeter wavelength interferometric systems in space. A concept of such the system, called Terahertz Exploration and Zooming-in for Astrophysics (THEZA), has been proposed in the framework of the ESA Call for White Papers for the Voayage 2050 long term plan in 2019. In the current paper we discuss several approaches for addressing technological challenges of the THEZA concept. In particular, we consider a novel configuration of a space-borne millimeter/sub-millimeter antenna which might resolve several bottlenecks in creating large precise mechanical structures. The paper also presents an overview of prospective space-qualified technologies of low-noise analogue front-end instrumentation for millimeter/sub-millimeter telescopes, data handling and processing. The paper briefly discusses approaches to the interferometric baseline state vector determination and synchronisation and heterodyning system. In combination with the original ESA Voyage 2050 White Paper, the current work sharpens the case for the next generation microarcsceond-level imaging instruments and provides starting points for further in-depth technology trade-off studies. △ Less

Submitted 27 April, 2022; v1 submitted 19 April, 2022; originally announced April 2022.

Comments: The paper on 29 pages contains 12 figures, accepted for publication in Acta Astronautica, based on the presentation at the 72nd IAC, Dubai, UAE, October 2021

Journal ref: Acta Astronautica 196 (2022) 314

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/

Event Horizon Telescope observations of the jet launching and collimation in Centaurus A

Authors: Michael Janssen, Heino Falcke, Matthias Kadler, Eduardo Ros, Maciek Wielgus, Kazunori Akiyama, Mislav Baloković, Lindy Blackburn, Katherine L. Bouman, Andrew Chael, Chi-kwan Chan, Koushik Chatterjee, Jordy Davelaar, Philip G. Edwards, Christian M. Fromm, José L. Gómez, Ciriaco Goddi, Sara Issaoun, Michael D. Johnson, Junhan Kim, Jun Yi Koay, Thomas P. Krichbaum, Jun Liu, Elisabetta Liuzzo, Sera Markoff , et al. (215 additional authors not shown)

Abstract: Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supe… ▽ More Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supermassive black holes (SMBHs) in Messier 87 and our galactic center. A large southern declination of $-43^{\circ}$ has however prevented VLBI imaging of Centaurus A below $λ1$cm thus far. Here, we show the millimeter VLBI image of the source, which we obtained with the Event Horizon Telescope at $228$GHz. Compared to previous observations, we image Centaurus A's jet at a tenfold higher frequency and sixteen times sharper resolution and thereby probe sub-lightday structures. We reveal a highly-collimated, asymmetrically edge-brightened jet as well as the fainter counterjet. We find that Centaurus A's source structure resembles the jet in Messier 87 on ${\sim}500r_g$ scales remarkably well. Furthermore, we identify the location of Centaurus A's SMBH with respect to its resolved jet core at $λ1.3$mm and conclude that the source's event horizon shadow should be visible at THz frequencies. This location further supports the universal scale invariance of black holes over a wide range of masses. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 27 pages, 9 figures. This is a post-peer-review, pre-copyedit version of an article published in Nature Astronomy. The final authenticated version is available online at: https://doi.org/10.1038/s41550-021-01417-w

Journal ref: Nature Astronomy, July 2021, Volume 5, p. 1017-1028

The Variability of the Black-Hole Image in M87 at the Dynamical Time Scale

Authors: Kaushik Satapathy, Dimitrios Psaltis, Feryal Ozel, Lia Medeiros, Sean T. Dougall, Chi-kwan Chan, Maciek Wielgus, Ben S. Prather, George N. Wong, Charles F. Gammie, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David R. Ball, Mislav Baloković, John Barrett, Bradford A. Benson, Dan Bintley, Lindy Blackburn, Raymond Blundell , et al. (213 additional authors not shown)

Abstract: The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expect… ▽ More The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure phase measurements on all six linearly independent non-trivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of $\sim3-5^\circ$. The only triangles that exhibit substantially higher variability ($\sim90-180^\circ$) are the ones with baselines that cross visibility amplitude minima on the $u-v$ plane, as expected from theoretical modeling. We used two sets of General Relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black-hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black-hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: Accepted for Publication in ApJ

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

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

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

Submitted 16 August, 2021; originally announced August 2021.

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

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

Sub-arcsecond imaging with the International LOFAR Telescope I. Foundational calibration strategy and pipeline

Authors: L. K. Morabito, N. J. Jackson, S. Mooney, F. Sweijen, S. Badole, P. Kukreti, D. Venkattu, C. Groeneveld, A. Kappes, E. Bonnassieux, A. Drabent, M. Iacobelli, J. H. Croston, P. N. Best, M. Bondi, J. R. Callingham, J. E. Conway, A. T. Deller, M. J. Hardcastle, J. P. McKean, G. K. Miley, J. Moldon, H. J. A. Röttgering, C. Tasse, T. W. Shimwell , et al. (49 additional authors not shown)

Abstract: [abridged] The International LOFAR Telescope is an interferometer with stations spread across Europe. With baselines of up to ~2,000 km, LOFAR has the unique capability of achieving sub-arcsecond resolution at frequencies below 200 MHz, although this is technically and logistically challenging. Here we present a calibration strategy that builds on previous high-resolution work with LOFAR. We give… ▽ More [abridged] The International LOFAR Telescope is an interferometer with stations spread across Europe. With baselines of up to ~2,000 km, LOFAR has the unique capability of achieving sub-arcsecond resolution at frequencies below 200 MHz, although this is technically and logistically challenging. Here we present a calibration strategy that builds on previous high-resolution work with LOFAR. We give an overview of the calibration strategy and discuss the special challenges inherent to enacting high-resolution imaging with LOFAR, and describe the pipeline, which is publicly available, in detail. We demonstrate the calibration strategy by using the pipeline on P205+55, a typical LOFAR Two-metre Sky Survey (LoTSS) pointing. We perform in-field delay calibration, solution referencing to other calibrators, self-calibration, and imaging of example directions of interest in the field. For this specific field and these ionospheric conditions, dispersive delay solutions can be transferred between calibrators up to ~1.5 degrees away, while phase solution transferral works well over 1 degree. We demonstrate a check of the astrometry and flux density scale. Imaging in 17 directions, the restoring beam is typically 0.3" x 0.2" although this varies slightly over the entire 5 square degree field of view. We achieve ~80 to 300 $μ$Jy/bm image rms noise, which is dependent on the distance from the phase centre; typical values are ~90 $μ$Jy/bm for the 8 hour observation with 48 MHz of bandwidth. Seventy percent of processed sources are detected, and from this we estimate that we should be able to image ~900 sources per LoTSS pointing. This equates to ~3 million sources in the northern sky, which LoTSS will entirely cover in the next several years. Future optimisation of the calibration strategy for efficient post-processing of LoTSS at high resolution (LoTSS-HR) makes this estimate a lower limit. △ Less

Submitted 16 August, 2021; originally announced August 2021.

Comments: Accepted to a special issue of A&A on sub-arcsecond imaging with LOFAR. 24 pages, 16 figures

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

Extended two-body problem for rotating rigid bodies

Authors: Alex Ho, Margrethe Wold, John T. Conway, Mohammad Poursina

Abstract: A new technique that utilizes surface integrals to find the force, torque and potential energy between two non-spherical, rigid bodies is presented. The method is relatively fast, and allows us to solve the full rigid two-body problem for pairs of spheroids and ellipsoids with 12 degrees of freedom. We demonstrate the method with two dimensionless test scenarios, one where tumbling motion develops… ▽ More A new technique that utilizes surface integrals to find the force, torque and potential energy between two non-spherical, rigid bodies is presented. The method is relatively fast, and allows us to solve the full rigid two-body problem for pairs of spheroids and ellipsoids with 12 degrees of freedom. We demonstrate the method with two dimensionless test scenarios, one where tumbling motion develops, and one where the motion of the bodies resemble spinning tops. We also test the method on the asteroid binary (66391) 1999 KW4, where both components are modelled either as spheroids or ellipsoids. The two different shape models have negligible effects on the eccentricity and semi-major axis, but have a larger impact on the angular velocity along the $z$-direction. In all cases, energy and total angular momentum is conserved, and the simulation accuracy is kept at the machine accuracy level. △ Less

Submitted 7 July, 2021; originally announced July 2021.

Comments: 24 pages, 9 figures, accepted for publication in Celestial Mechanics and Dynamical Astronomy

The planar two-body problem for spheroids and disks

Authors: M. Wold, J. T. Conway

Abstract: We outline a new method suggested by Conway (2016) for solving the two-body problem for solid bodies of spheroidal or ellipsoidal shape. The method is based on integrating the gravitational potential of one body over the surface of the other body. When the gravitational potential can be analytically expressed (as for spheroids or ellipsoids), the gravitational force and mutual gravitational potent… ▽ More We outline a new method suggested by Conway (2016) for solving the two-body problem for solid bodies of spheroidal or ellipsoidal shape. The method is based on integrating the gravitational potential of one body over the surface of the other body. When the gravitational potential can be analytically expressed (as for spheroids or ellipsoids), the gravitational force and mutual gravitational potential can be formulated as a surface integral instead of a volume integral, and solved numerically. If the two bodies are infinitely thin disks, the surface integral has an analytical solution. The method is exact as the force and mutual potential appear in closed-form expressions, and does not involve series expansions with subsequent truncation errors. In order to test the method, we solve the equations of motion in an inertial frame, and run simulations with two spheroids and two infinitely thin disks, restricted to torque-free planar motion. The resulting trajectories display precession patterns typical for non-Keplerian potentials. We follow the conservation of energy and orbital angular momentum, and also investigate how the spheroid model approaches the two cases where the surface integral can be solved analytically, i.e. for point masses and infinitely thin disks. △ Less

Submitted 25 June, 2021; originally announced June 2021.

Comments: 15 pages, 5 figures, accepted for publication in Celestial Mechanics and Dynamical Astronomy

Journal ref: Celest Mech Dyn Astr, 133, 27 (2021)

Burst timescales and luminosities link young pulsars and fast radio bursts

Authors: K. Nimmo, J. W. T. Hessels, F. Kirsten, A. Keimpema, J. M. Cordes, M. P. Snelders, D. M. Hewitt, R. Karuppusamy, A. M. Archibald, V. Bezukovs, M. Bhardwaj, R. Blaauw, S. T. Buttaccio, T. Cassanelli, J. E. Conway, A. Corongiu, R. Feiler, E. Fonseca, O. Forssen, M. Gawronski, M. Giroletti, M. A. Kharinov, C. Leung, M. Lindqvist, G. Maccaferri , et al. (23 additional authors not shown)

Abstract: Fast radio bursts (FRBs) are extragalactic radio flashes of unknown physical origin. Their high luminosities and short durations require extreme energy densities, like those found in the vicinity of neutron stars and black holes. Studying the burst intensities and polarimetric properties on a wide range of timescales, from milliseconds down to nanoseconds, is key to understanding the emission mech… ▽ More Fast radio bursts (FRBs) are extragalactic radio flashes of unknown physical origin. Their high luminosities and short durations require extreme energy densities, like those found in the vicinity of neutron stars and black holes. Studying the burst intensities and polarimetric properties on a wide range of timescales, from milliseconds down to nanoseconds, is key to understanding the emission mechanism. However, high-time-resolution studies of FRBs are limited by their unpredictable activity levels, available instrumentation and temporal broadening in the intervening ionised medium. Here we show that the repeating FRB 20200120E can produce isolated shots of emission as short as about 60 nanoseconds in duration, with brightness temperatures as high as $3\times 10^{41}$ K (excluding relativistic effects), comparable to `nano-shots' from the Crab pulsar. Comparing both the range of timescales and luminosities, we find that FRB 20200120E observationally bridges the gap between known Galactic young pulsars and magnetars, and the much more distant extragalactic FRBs. This suggests a common magnetically powered emission mechanism spanning many orders of magnitude in timescale and luminosity. In this work, we probe a relatively unexplored region of the short-duration transient phase space; we highlight that there likely exists a population of ultra-fast radio transients at nanosecond to microsecond timescales, which current FRB searches are insensitive to. △ Less

Submitted 29 September, 2021; v1 submitted 24 May, 2021; originally announced May 2021.

Comments: Submitted. Comments welcome

A repeating fast radio burst source in a globular cluster

Authors: F. Kirsten, B. Marcote, K. Nimmo, J. W. T. Hessels, M. Bhardwaj, S. P. Tendulkar, A. Keimpema, J. Yang, M. P. Snelders, P. Scholz, A. B. Pearlman, C. J. Law, W. M. Peters, M. Giroletti, Z. Paragi, C. Bassa, D. M. Hewitt, U. Bach, V. Bezrukovs, M. Burgay, S. T. Buttaccio, J. E. Conway, A. Corongiu, R. Feiler, O. Forssén , et al. (41 additional authors not shown)

Abstract: Fast radio bursts (FRBs) are exceptionally luminous flashes of unknown physical origin, reaching us from other galaxies (Petroff et al. 2019). Most FRBs have only ever been seen once, while others flash repeatedly, though sporadically (Spitler et al. 2016, CHIME/FRB Collaboration et al. 2021). Many models invoke magnetically powered neutron stars (magnetars) as the engines producing FRB emission (… ▽ More Fast radio bursts (FRBs) are exceptionally luminous flashes of unknown physical origin, reaching us from other galaxies (Petroff et al. 2019). Most FRBs have only ever been seen once, while others flash repeatedly, though sporadically (Spitler et al. 2016, CHIME/FRB Collaboration et al. 2021). Many models invoke magnetically powered neutron stars (magnetars) as the engines producing FRB emission (Margalit & Metzger 2018, CHIME/FRB Collaboration et al. 2020). Recently, CHIME/FRB announced the discovery (Bhardwaj et al. 2021) of the repeating FRB 20200120E, coming from the direction of the nearby grand design spiral galaxy M81. Four potential counterparts at other observing wavelengths were identified (Bhardwaj et al. 2021) but no definitive association with these sources, or M81, could be made. Here we report an extremely precise localisation of FRB 20200120E, which allows us to associate it with a globular cluster (GC) in the M81 galactic system and to place it ~2pc offset from the optical center of light of the GC. This confirms (Bhardwaj et al. 2021) that FRB 20200120E is 40 times closer than any other known extragalactic FRB. Because such GCs host old stellar populations, this association strongly challenges FRB models that invoke young magnetars formed in a core-collapse supernova as powering FRB emission. We propose, instead, that FRB 20200120E is a highly magnetised neutron star formed via either accretion-induced collapse of a white dwarf or via merger of compact stars in a binary system (Margalit et al. 2019). Alternative scenarios involving compact binary systems, efficiently formed inside globular clusters, could also be responsible for the observed bursts. △ Less

Submitted 29 September, 2021; v1 submitted 24 May, 2021; originally announced May 2021.

Comments: Submitted. Comments welcome

The Polarized Image of a Synchrotron Emitting Ring of Gas Orbiting a Black Hole

Authors: Ramesh Narayan, Daniel C. M. Palumbo, Michael D. Johnson, Zachary Gelles, Elizabeth Himwich, Dominic O. Chang, Angelo Ricarte, Jason Dexter, Charles F. Gammie, Andrew A. Chael, The Event Horizon Telescope Collaboration, :, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David Ball, Mislav Balokovic, John Barrett, Bradford A. Benson, Dan Bintley , et al. (215 additional authors not shown)

Abstract: Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equ… ▽ More Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov (2002) and conservation of the Walker-Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images. △ Less

Submitted 13 May, 2021; v1 submitted 4 May, 2021; originally announced May 2021.

Comments: 29 pages, 11 figures, published in ApJ on May 3

Journal ref: ApJ 912 35 (2021)

Low-frequency radio spectra of submillimetre galaxies in the Lockman Hole

Authors: J. Ramasawmy, J. E. Geach, M. J. Hardcastle, P. N. Best, M. Bonato, M. Bondi, G. Calistro Rivera, R. K. Cochrane, J. E. Conway, K. Coppin, K. J. Duncan, J. S. Dunlop, M. Franco, C. García-Vergara, M. J. Jarvis, R. Kondapally, I. McCheyne, I. Prandoni, H. J. A. Röttgering, D. J. B. Smith, C. Tasse, L. Wang

Abstract: We investigate the radio properties of a sample of 53 sources selected at 850 $μ$m from the SCUBA-2 Cosmology Legacy Survey using new deep, low-frequency radio imaging of the Lockman Hole field from the Low Frequency Array. Combining these data with additional radio observations from the GMRT and the JVLA, we find a variety of radio spectral shapes and luminosities within our sample despite their… ▽ More We investigate the radio properties of a sample of 53 sources selected at 850 $μ$m from the SCUBA-2 Cosmology Legacy Survey using new deep, low-frequency radio imaging of the Lockman Hole field from the Low Frequency Array. Combining these data with additional radio observations from the GMRT and the JVLA, we find a variety of radio spectral shapes and luminosities within our sample despite their similarly bright submillimetre flux densities. We characterise their spectral shapes in terms of multi-band radio spectral indices. Finding strong spectral flattening at low frequencies in ~20% of sources, we investigate the differences between sources with extremely flat low-frequency spectra and those with `normal' radio spectral indices. As there are no other statistically significant differences between the two subgroups of our sample as split by the radio spectral index, we suggest that any differences are undetectable in galaxy-averaged properties that we can observe with our unresolved images, and likely relate to galaxy properties that we cannot resolve, on scales $\lesssim$ 1 kpc. We attribute the observed spectral flattening in the radio to free-free absorption, proposing that those sources with significant low-frequency spectral flattening have a clumpy distribution of star-forming gas. We estimate an average spatial extent of absorbing material of at most several hundred parsecs to produce the levels of absorption observed in the radio spectra. This estimate is consistent with the highest-resolution observations of submillimetre galaxies in the literature, which find examples of non-uniform dust distributions on scales of ~100 pc, with evidence for clumps and knots in the interstellar medium. Additionally, we find two bright (> 6 mJy) submm sources undetected at all other wavelengths. We speculate that these objects may be very high redshift sources, likely residing at z > 4. △ Less

Submitted 17 March, 2021; originally announced March 2021.

Comments: 15 pages, 10 figures, accepted by A&A

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

On stable H-C-N-O compounds at high pressure

Authors: Lewis J. Conway, Chris J. Pickard, Andreas Hermann

Abstract: The make-up of the outer planets, and many of their moons, are dominated by matter from the H-C-N-O chemical space, commonly assumed to originate from mixtures of hydrogen and the planetary ices H$_2$O, CH$_4$, and NH$_3$. In their interiors, these ices experience extreme pressure conditions, around 5 Mbar at the Neptune mantle-core boundary, and it is expected that they undergo phase transitions,… ▽ More The make-up of the outer planets, and many of their moons, are dominated by matter from the H-C-N-O chemical space, commonly assumed to originate from mixtures of hydrogen and the planetary ices H$_2$O, CH$_4$, and NH$_3$. In their interiors, these ices experience extreme pressure conditions, around 5 Mbar at the Neptune mantle-core boundary, and it is expected that they undergo phase transitions, decompose, and form entirely new compounds. In turn, this determines planets' interior structure, thermal history, magnetic field generation, etc. Despite its importance, the H-C-N-O space has not been surveyed systematically. Asked simply: at high-pressure conditions, what compounds emerge within this space, and what governs their stability? Here, we report on results from an unbiased crystal structure search amongst H-C-N-O compounds at 5 Mbar to answer this question. △ Less

Submitted 26 November, 2020; originally announced November 2020.

Gravitational Test Beyond the First Post-Newtonian Order with the Shadow of the M87 Black Hole

Authors: Dimitrios Psaltis, Lia Medeiros, Pierre Christian, Feryal Ozel, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Keiichi Asada, Rebecca Azulay, David Ball, Mislav Balokovic, John Barrett, Dan Bintley, Lindy Blackburn, Wilfred Boland, Geoffrey C. Bower, Michael Bremer, Christiaan D. Brinkerink, Roger Brissenden, Silke Britzen, Dominique Broguiere, Thomas Bronzwaer, Do-Young Byun, John E. Carlstrom, Andrew Chael , et al. (163 additional authors not shown)

Abstract: The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the p… ▽ More The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the predicted black-hole shadows that are inconsistent with even the current EHT measurements. We use numerical calculations of regular, parametric, non-Kerr metrics to identify the common characteristic among these different parametrizations that control the predicted shadow size. We show that the shadow-size measurements place significant constraints on deviation parameters that control the second post-Newtonian and higher orders of each metric and are, therefore, inaccessible to weak-field tests. The new constraints are complementary to those imposed by observations of gravitational waves from stellar-mass sources. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: Physical Review Letters

VLBI20-30: a scientific roadmap for the next decade -- The future of the European VLBI Network

Authors: Tiziana Venturi, Zsolt Paragi, Michael Lindqvist, Anna Bartkiewicz, Rob Beswick, Tamara Bogdanović, Walter Brisken, Patrick Charlot, Francisco Colomer, John Conway, Sándor Frey, José Carlos Guirado, Leonid Gurvits, Huib van Langevelde, Andrei Lobanov, John McKean, Raffaella Morganti, Tom Muxlow, Miguel Pérez-Torres, Kazi Rygl, Robert Schulz, Arpad Szomoru, Pablo de Vicente, Tao An, Guillem Anglada , et al. (55 additional authors not shown)

Abstract: This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, dem… ▽ More This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, demonstrate the EVN's potential to generate new and exciting results that will transform our view of the cosmos. Together with e-MERLIN, the EVN provides a range of baseline lengths that permit unique studies of faint radio sources to be made over a wide range of spatial scales. The science cases are reviewed in six chapters that cover the following broad areas: cosmology, galaxy formation and evolution, innermost regions of active galactic nuclei, explosive phenomena and transients, stars and stellar masers in the Milky Way, celestial reference frames and space applications. The document concludes with identifying the synergies with other radio, as well as multi-band/multi-messenger instruments, and provide the recommendations for future improvements. The appendices briefly describe other radio VLBI arrays, the technological framework for EVN developments, and a selection of spectral lines of astrophysical interest below 100 GHz. The document includes a glossary for non-specialists, and a list of acronyms at the end. △ Less

Submitted 5 July, 2020; originally announced July 2020.

Comments: Editors: Tiziana Venturi, Zsolt Paragi, Michael Lindqvist. "EVN Vision Document (2020)", 7 chapters, 3 appendices, 191 pages, 64 figures. Chapter coordinators and all the contributors are listed in the document

Markov Chains for Horizons (MARCH). I. Identifying Biases in Fitting Theoretical Models to Event Horizon Telescope Observations

Authors: Dimitrios Psaltis, Feryal Ozel, Lia Medeiros, Pierre Christian, Junhan Kim, Chi-kwan Chan, Landen J. Conway, Carolyn A. Raithel, Dan Marrone, Tod R. Lauer

Abstract: We introduce a new Markov Chain Monte Carlo (MCMC) algorithm with parallel tempering for fitting theoretical models of horizon-scale images of black holes to the interferometric data from the Event Horizon Telescope (EHT). The algorithm implements forms of the noise distribution in the data that are accurate for all signal-to-noise ratios. In addition to being trivially parallelizable, the algorit… ▽ More We introduce a new Markov Chain Monte Carlo (MCMC) algorithm with parallel tempering for fitting theoretical models of horizon-scale images of black holes to the interferometric data from the Event Horizon Telescope (EHT). The algorithm implements forms of the noise distribution in the data that are accurate for all signal-to-noise ratios. In addition to being trivially parallelizable, the algorithm is optimized for high performance, achieving 1 million MCMC chain steps in under 20 seconds on a single processor. We use synthetic data for the 2017 EHT coverage of M87 that are generated based on analytic as well as General Relativistic Magnetohydrodynamic (GRMHD) model images to explore several potential sources of biases in fitting models to sparse interferometric data. We demonstrate that a very small number of data points that lie near salient features of the interferometric data exert disproportionate influence on the inferred model parameters. We also show that the preferred orientations of the EHT baselines introduce significant biases in the inference of the orientation of the model images. Finally, we discuss strategies that help identify the presence and severity of such biases in realistic applications. △ Less

Submitted 19 May, 2020; originally announced May 2020.

Comments: ApJ, submitted

SYMBA: An end-to-end VLBI synthetic data generation pipeline

Authors: F. Roelofs, M. Janssen, I. Natarajan, R. Deane, J. Davelaar, H. Olivares, O. Porth, S. N. Paine, K. L. Bouman, R. P. J. Tilanus, I. M. van Bemmel, H. Falcke, K. Akiyama, A. Alberdi, W. Alef, K. Asada, R. Azulay, A. Baczko, D. Ball, M. Baloković, J. Barrett, D. Bintley, L. Blackburn, W. Boland, G. C. Bower , et al. (183 additional authors not shown)

Abstract: Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabili… ▽ More Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images. △ Less

Submitted 2 April, 2020; originally announced April 2020.

Comments: 20 pages, 15 figures, accepted for publication in A&A

A parsec-scale radio jet launched by the central intermediate-mass black hole in the dwarf galaxy SDSS J090613.77+561015.2?

Authors: Jun Yang, Leonid I. Gurvits, Zsolt Paragi, Sandor Frey, John E. Conway, Xiang Liu, Lang Cui

Abstract: The population of intermediate-mass black holes (IMBHs) in nearby dwarf galaxies plays an important "ground truth" role in exploring black hole formation and growth in the early Universe. In the dwarf elliptical galaxy SDSS J090613.77+561015.2 (z=0.0465), an accreting IMBH has been revealed by optical and X-ray observations. Aiming to search for possible radio core and jet associated with the IMBH… ▽ More The population of intermediate-mass black holes (IMBHs) in nearby dwarf galaxies plays an important "ground truth" role in exploring black hole formation and growth in the early Universe. In the dwarf elliptical galaxy SDSS J090613.77+561015.2 (z=0.0465), an accreting IMBH has been revealed by optical and X-ray observations. Aiming to search for possible radio core and jet associated with the IMBH, we carried out very long baseline interferometry (VLBI) observations with the European VLBI Network (EVN) at 1.66 GHz. Our imaging results show that there are two 1-mJy components with a separation of about 52 mas (projected distance 47 pc) and the more compact component is located within the 1-sigma error circle of the optical centroid from available Gaia astrometry. Based on their positions, elongated structures and relatively high brightness temperatures, as well as the absence of star-forming activity in the host galaxy, we argue that the radio morphology originates from the jet activity powered by the central IMBH. The existence of the large-scale jet implies that violent jet activity might occur in the early epochs of black hole growth and thus help to regulate the co-evolution of black holes and galaxies. △ Less

Submitted 25 March, 2020; originally announced March 2020.

Comments: Accepted for publication as a Letter in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society Letters, Vol. 495, pp. L71-L75 (2020)

TeraHertz Exploration and Zooming-in for Astrophysics (THEZA): ESA Voyage 2050 White Paper

Authors: Leonid I. Gurvits, Zsolt Paragi, Viviana Casasola, John Conway, Jordy Davelaar, Heino Falcke, Rob Fender, Sándor Frey, Christian M. Fromm, Cristina García Miró, Michael A. Garrett, Marcello Giroletti, Ciriaco Goddi, José-Luis Gómez, Jeffrey van der Gucht, José Carlos Guirado, Zoltán Haiman, Frank Helmich, Elizabeth Humphreys, Violette Impellizzeri, Michael Kramer, Michael Lindqvist, Hendrik Linz, Elisabetta Liuzzo, Andrei P. Lobanov , et al. (10 additional authors not shown)

Abstract: This paper presents the ESA Voyage 2050 White Paper for a concept of TeraHertz Exploration and Zooming-in for Astrophysics (THEZA). It addresses the science case and some implementation issues of a space-borne radio interferometric system for ultra-sharp imaging of celestial radio sources at the level of angular resolution down to (sub-) microarcseconds. THEZA focuses at millimetre and sub-millime… ▽ More This paper presents the ESA Voyage 2050 White Paper for a concept of TeraHertz Exploration and Zooming-in for Astrophysics (THEZA). It addresses the science case and some implementation issues of a space-borne radio interferometric system for ultra-sharp imaging of celestial radio sources at the level of angular resolution down to (sub-) microarcseconds. THEZA focuses at millimetre and sub-millimetre wavelengths (frequencies above $\sim$300~GHz), but allows for science operations at longer wavelengths too. The THEZA concept science rationale is focused on the physics of spacetime in the vicinity of supermassive black holes as the leading science driver. The main aim of the concept is to facilitate a major leap by providing researchers with orders of magnitude improvements in the resolution and dynamic range in direct imaging studies of the most exotic objects in the Universe, black holes. The concept will open up a sizeable range of hitherto unreachable parameters of observational astrophysics. It unifies two major lines of development of space-borne radio astronomy of the past decades: Space VLBI (Very Long Baseline Interferometry) and mm- and sub-mm astrophysical studies with "single dish" instruments. It also builds upon the recent success of the Earth-based Event Horizon Telescope (EHT) -- the first-ever direct image of a shadow of the super-massive black hole in the centre of the galaxy M87. As an amalgam of these three major areas of modern observational astrophysics, THEZA aims at facilitating a breakthrough in high-resolution high image quality studies in the millimetre and sub-millimetre domain of the electromagnetic spectrum. △ Less

Submitted 28 May, 2021; v1 submitted 28 August, 2019; originally announced August 2019.

Comments: White Paper submitted in response to the ESA Call Voyage 2050

Journal ref: Experimental Astronomy, 2021

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

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

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

Submitted 8 July, 2019; originally announced July 2019.

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

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

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

The case for gravitational millilensing in the multiply--imaged quasar B1152+199

Authors: Saghar Asadi, Erik Zackrisson, Eskil Varenius, Emily Freeland, John Conway, Kaj Wiik

Abstract: Previous Very Long Baseline Interferometry (VLBI) observations of the quasar B1152+199 at 5GHz has revealed two images of a strongly lensed jet with seemingly discordant morphologies. Whereas the jet appears straight in one of the images, the other exhibits slight curvature on milliarcsecond scales. This is unexpected from the lensing solution and has been interpreted as possible evidence for seco… ▽ More Previous Very Long Baseline Interferometry (VLBI) observations of the quasar B1152+199 at 5GHz has revealed two images of a strongly lensed jet with seemingly discordant morphologies. Whereas the jet appears straight in one of the images, the other exhibits slight curvature on milliarcsecond scales. This is unexpected from the lensing solution and has been interpreted as possible evidence for secondary, small-scale lensing (millilensing) by a compact object with a mass of $~10^5$-$10^7\ M_\odot$ located close to the curved image. The probability for such a superposition is extremely low unless the millilens population has very high surface number density. Here, we revisit the case for millilensing in B1152+199 by combining new global-VLBI data at 8.4GHz with two datasets from the European VLBI Network (EVN) at 5GHz (archival) and at 22GHz (new dataset), and the previously published 5GHz Very Long Baseline Array (VLBA) data. We find that the new data with a more circular synthesized beam, exhibits no apparent milliarcsecond-scale curvature in image B. Various observations of the object spanning $\sim$15 years apart enable us to improve the constraints on lens system (thanks also to the improved astrometry resulting from 22GHz observations) to the point that the only plausible explanation left for the apparent curvature is the artifact due to the shape of the synthesized beam. △ Less

Submitted 14 November, 2018; originally announced November 2018.

Comments: 9 pages, 4 figures, submitted to MNRAS

LOFAR MSSS: Flattening low-frequency radio continuum spectra of nearby galaxies

Authors: K. T. Chyży, W. Jurusik, J. Piotrowska, B. Nikiel-Wroczyński, V. Heesen, V. Vacca, N. Nowak, R. Paladino, P. Surma, S. S. Sridhar, G. Heald, R. Beck, J. Conway, K. Sendlinger, M. Curyło, D. Mulcahy, J. W. Broderick, M. J. Hardcastle, J. R. Callingham, G. Gürkan, M. Iacobelli, H. J. A. Röttgering, B. Adebahr, A. Shulevski, R. -J. Dettmar , et al. (13 additional authors not shown)

Abstract: The shape of low-frequency radio continuum spectra of normal galaxies is not well understood, the key question being the role of physical processes such as thermal absorption in shaping them. In this work we take advantage of the LOFAR Multifrequency Snapshot Sky Survey (MSSS) to investigate such spectra for a large sample of nearby star-forming galaxies. Using the measured 150MHz flux densities f… ▽ More The shape of low-frequency radio continuum spectra of normal galaxies is not well understood, the key question being the role of physical processes such as thermal absorption in shaping them. In this work we take advantage of the LOFAR Multifrequency Snapshot Sky Survey (MSSS) to investigate such spectra for a large sample of nearby star-forming galaxies. Using the measured 150MHz flux densities from the LOFAR MSSS survey and literature flux densities at various frequencies we have obtained integrated radio spectra for 106 galaxies. The spectra are explained through the use of a three-dimensional model of galaxy radio emission, and radiation transfer dependent on the galaxy viewing angle and absorption processes. Spectra of our galaxies are generally flatter at lower compared to higher frequencies but as there is no tendency for the highly inclined galaxies to have more flattened low-frequency spectra, we argue that the observed flattening is not due to thermal absorption, contradicting the suggestion of Israel & Mahoney (1990). According to our modelled radio maps for M51-like galaxies, the free-free absorption effects can be seen only below 30MHz and in the global spectra just below 20MHz, while in the spectra of starburst galaxies, like M82, the flattening due to absorption is instead visible up to higher frequencies of about 150MHz. Locally, within galactic disks, the absorption effects are distinctly visible in M51-like galaxies as spectral flattening around 100-200MHz in the face-on objects, and as turnovers in the edge-on ones, while in M82-like galaxies there are strong turnovers at frequencies above 700MHz, regardless of viewing angle. Our modelling suggests that the weak spectral flattening observed in the nearby galaxies studied here results principally from synchrotron spectral curvature due to cosmic ray energy losses and propagation effects. △ Less

Submitted 31 August, 2018; v1 submitted 30 August, 2018; originally announced August 2018.

Comments: Accepted for publication in Astronomy and Astrophysics (corrected typo in the title)

Journal ref: A&A 619, A36 (2018)

Continuum and Spectral Line Radiation from a Random Clumpy Medium

Authors: John E. Conway, Moshe Elitzur, Rodrigo Parra

Abstract: We present a formalism for continuum and line emission from random clumpy media together with its application to problems of current interest, including CO spectral lines from ensembles of clouds and radio emission from HII regions, supernovae and star-forming regions. For line emission we find that the effects of clump opacity on observed line ratios can be indistinguishable from variations of in… ▽ More We present a formalism for continuum and line emission from random clumpy media together with its application to problems of current interest, including CO spectral lines from ensembles of clouds and radio emission from HII regions, supernovae and star-forming regions. For line emission we find that the effects of clump opacity on observed line ratios can be indistinguishable from variations of intrinsic line strengths, adding to the difficulties in determining abundances from line observations. Our formalism is applicable to arbitrary distributions of cloud properties, provided the cloud volume filling factor is small; numerical simulations show it to hold up to filling factors of about 10%. We show that irrespective of the complexity of the cloud ensemble, the radiative effect of clumpiness can be parametrized at each frequency by a single multiplicative correction to the overall optical depth; this multiplier is derived from appropriate averaging over individual cloud properties. Our main finding is that cloud shapes have only a negligible effect on radiation propagation in clumpy media; the results of calculations employing point-like clouds are practically indistinguishable from those for finite-size clouds with arbitrary geometrical shapes. △ Less

Submitted 22 August, 2018; originally announced August 2018.

Comments: ApJ, to be published

Reliable detection and characterization of low-frequency polarized sources in the LOFAR M51 field

Authors: A. Neld, C. Horellou, D. D. Mulcahy, R. Beck, S. Bourke, T. D. Carozzi, K. T. Chyży, J. E. Conway, J. S. Farnes, A. Fletcher, M. Haverkorn, G. Heald, A. Horneffer, B. Nikiel-Wroczyński, R. Paladino, S. S. Sridhar, C. L. Van Eck

Abstract: The new generation of broad-band radio continuum surveys will provide large data sets with polarization information. New algorithms need to be developed to extract reliable catalogs of linearly polarized sources that can be used to characterize those sources and produce a dense rotation measure (RM) grid to probe magneto-ionized structures along the line of sight via Faraday rotation. The aim of t… ▽ More The new generation of broad-band radio continuum surveys will provide large data sets with polarization information. New algorithms need to be developed to extract reliable catalogs of linearly polarized sources that can be used to characterize those sources and produce a dense rotation measure (RM) grid to probe magneto-ionized structures along the line of sight via Faraday rotation. The aim of the paper is to develop a computationally efficient and rigorously defined source-finding algorithm for linearly polarized sources. We used a calibrated data set from the LOw Frequency ARray (LOFAR) at 150 MHz centered on the nearby galaxy M51 to search for polarized background sources. With a new imaging software, we re-imaged the field at a resolution of 18''x15'' and cataloged a total of about 3000 continuum sources within 2.5 degrees of the center of M51. We made small Stokes Q and U images centered on each source brighter than 100 mJy in total intensity (201 sources) and used RM synthesis to create corresponding Faraday cubes that were analyzed individually. For each source, the noise distribution function was determined from a subset of measurements at high Faraday depths where no polarization is expected; the peaks in polarized intensity in the Faraday spectrum were identified and the p-value of each source was calculated. Finally, the false discovery rate method was applied to the list of p-values to produce a list of polarized sources and quantify the reliability of the detections. We also analyzed sources fainter than 100 mJy but that were reported as polarized in the literature at at least another radio frequency. Of the 201 sources that were searched for polarization, 6 polarized sources were detected confidently (with a false discovery rate of 5 percent). This corresponds to a number density of 1 polarized source per 3.3 square degrees, or 0.3 source/square degree. [Abridged] △ Less

Submitted 24 July, 2018; originally announced July 2018.

Comments: 17 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A 617, A136 (2018)

SEPIA - a new single pixel receiver at the APEX Telescope

Authors: V. Belitsky, I. Lapkin, M. Fredrixon, D. Meledin, E. Sundin, B. Billade, S. -E. Ferm, A. Pavolotsky, H. Rashid, M. Strandberg, V. Desmaris, A. Ermakov, S. Krause, M. Olberg, P. Aghdam, S. Shafiee, P. Bergman, E. De Beck, H. Olofsson, J. Conway, C. De Breuck, K. Immer, P. Yagoubov, F. M. Montenegro-Montes, K. Torstensson , et al. (9 additional authors not shown)

Abstract: Context: We describe the new SEPIA (Swedish-ESO PI Instrument for APEX) receiver, which was designed and built by the Group for Advanced Receiver Development (GARD), at Onsala Space Observatory (OSO) in collaboration with ESO. It was installed and commissioned at the APEX telescope during 2015 with an ALMA Band 5 receiver channel and updated with a new frequency channel (ALMA Band 9) in February 2… ▽ More Context: We describe the new SEPIA (Swedish-ESO PI Instrument for APEX) receiver, which was designed and built by the Group for Advanced Receiver Development (GARD), at Onsala Space Observatory (OSO) in collaboration with ESO. It was installed and commissioned at the APEX telescope during 2015 with an ALMA Band 5 receiver channel and updated with a new frequency channel (ALMA Band 9) in February 2016. Aims: This manuscript aims to provide, for observers who use the SEPIA receiver, a reference in terms of the hardware description, optics and performance as well as the commissioning results. Methods: Out of three available receiver cartridge positions in SEPIA, the two current frequency channels, corresponding to ALMA Band 5, the RF band 158--211 GHz, and Band 9, the RF band 600--722 GHz, provide state-of-the-art dual polarization receivers. The Band 5 frequency channel uses 2SB SIS mixers with an average SSB noise temperature around 45K with IF (intermediate frequency) band 4--8 GHz for each sideband providing total 4x4 GHz IF band. The Band 9 frequency channel uses DSB SIS mixers with a noise temperature of 75--125K with IF band 4--12 GHz for each polarization. Results: Both current SEPIA receiver channels are available to all APEX observers. △ Less

Submitted 20 December, 2017; originally announced December 2017.

Comments: Accepted for publication in Astronomy & Astrophysics

An Overview of the MHONGOOSE Survey: Observing Nearby Galaxies with MeerKAT

Authors: W. J. G. de Blok, E. A. K. Adams, P. Amram, E. Athanassoula, I. Bagetakos, C. Balkowski, M. A. Bershady, R. Beswick, F. Bigiel, S. -L. Blyth, A. Bosma, R. S. Booth, A. Bouchard, E. Brinks, C. Carignan, L. Chemin, F. Combes, J. Conway, E. C. Elson, J. English, B. Epinat, B. S. Frank, J. Fiege, F. Fraternali, J. S. Gallagher , et al. (43 additional authors not shown)

Abstract: MHONGOOSE is a deep survey of the neutral hydrogen distribution in a representative sample of 30 nearby disk and dwarf galaxies with HI masses from 10^6 to ~10^{11} M_sun, and luminosities from M_R ~ -12 to M_R ~ -22. The sample is selected to uniformly cover the available range in log(M_HI). Our extremely deep observations, down to HI column density limits of well below 10^{18} cm^{-2} - or a few… ▽ More MHONGOOSE is a deep survey of the neutral hydrogen distribution in a representative sample of 30 nearby disk and dwarf galaxies with HI masses from 10^6 to ~10^{11} M_sun, and luminosities from M_R ~ -12 to M_R ~ -22. The sample is selected to uniformly cover the available range in log(M_HI). Our extremely deep observations, down to HI column density limits of well below 10^{18} cm^{-2} - or a few hundred times fainter than the typical HI disks in galaxies - will directly detect the effects of cold accretion from the intergalactic medium and the links with the cosmic web. These observations will be the first ever to probe the very low-column density neutral gas in galaxies at these high resolutions. Combination with data at other wavelengths, most of it already available, will enable accurate modelling of the properties and evolution of the mass components in these galaxies and link these with the effects of environment, dark matter distribution, and other fundamental properties such as halo mass and angular momentum. MHONGOOSE can already start addressing some of the SKA-1 science goals and will provide a comprehensive inventory of the processes driving the transformation and evolution of galaxies in the nearby universe at high resolution and over 5 orders of magnitude in column density. It will be a Nearby Galaxies Legacy Survey that will be unsurpassed until the advent of the SKA, and can serve as a highly visible, lasting statement of MeerKAT's capabilities. △ Less

Submitted 25 September, 2017; originally announced September 2017.

Comments: To be published in "MeerKAT Science: On the Pathway to the SKA". Proceedings of Science. Workshop held 25-27 May, 2016 Stellenbosch, South Africa

An atomic hydrogen bridge fueling NGC 4418 with gas from VV 655

Authors: E. Varenius, F. Costagliola, H. -R. Klöckner, S. Aalto, H. Spoon, I. Martí-Vidal, J. E. Conway

Abstract: The galaxy NGC 4418 harbours a compact ($<20$ pc) core with a very high bolometric luminosity ($\sim10^{11}$L$_\odot$). As most of the galaxy's energy output comes from this small region, it is of interest to determine what fuels this intense activity. An interaction with VV 655 has been proposed, where gas aquired by NGC 4418 could trigger intense star formation and/or black hole accretion in the… ▽ More The galaxy NGC 4418 harbours a compact ($<20$ pc) core with a very high bolometric luminosity ($\sim10^{11}$L$_\odot$). As most of the galaxy's energy output comes from this small region, it is of interest to determine what fuels this intense activity. An interaction with VV 655 has been proposed, where gas aquired by NGC 4418 could trigger intense star formation and/or black hole accretion in the centre. We aim to constrain the interaction hypothesis by studying neutral hydrogen structures around the two galaxies. We present observations at 1.4 GHz with the Very Large Array of radio continuum as well as emission and absorption from atomic hydrogen. Gaussian distributions are fitted to observed HI emission and absorption spectra. An atomic HI bridge is seen in emission, connecting NGC 4418 to VV 655. While NGC 4418 is bright in continuum emission and seen in HI absorption, VV 655 is barely detected in the continuum but show bright HI emission (M$_\mathrm{HI}\sim10^9$ M$_\odot$). We estimate SFRs from 1.4 GHz of 3.2 M$_\odot$ yr$^{-1}$ and 0.13 M$_\odot$ yr$^{-1}$ for NGC 4418 and VV 655 respectively. Systemic HI velocities of 2202$\pm$20 km s$^{-1}$ (emission) and 2105.4$\pm$10 km s$^{-1}$ (absorption) are measured for VV 655 and NGC 4418 respectively. Redshifted HI absorption is seen towards NGC 4418, suggesting gas infall. Blueshifted HI-emission is seen north-west of NGC 4418, which we interpret as a continuation of the outflow previously discussed by Sakamoto et al. (2013). The morphology and velocity structure seen in HI is consistent with an interaction scenario, where gas was transferred from VV 655 to NGC 4418, and may fuel the activity in the centre. △ Less

Submitted 2 May, 2017; originally announced May 2017.

Comments: 7 pages, 6 figures

Journal ref: A&A 607, A43 (2017)

The population of SNe/SNRs in the starburst galaxy Arp 220. A self-consistent analysis of 20 years of VLBI monitoring

Authors: E. Varenius, J. E. Conway, F. Batejat, I. Martí-Vidal, M. A. Pérez-Torres, S. Aalto, A. Alberdi, C. J. Lonsdale, P. Diamond

Abstract: The nearby ultra-luminous infrared galaxy (ULIRG) Arp 220 is an excellent laboratory for studies of extreme astrophysical environments. For 20 years, Very Long Baseline Interferometry (VLBI) has been used to monitor a population of compact sources thought to be supernovae (SNe), supernova remnants (SNRs) and possibly active galactic nuclei (AGNs). Using new and archival VLBI data spanning 20 years… ▽ More The nearby ultra-luminous infrared galaxy (ULIRG) Arp 220 is an excellent laboratory for studies of extreme astrophysical environments. For 20 years, Very Long Baseline Interferometry (VLBI) has been used to monitor a population of compact sources thought to be supernovae (SNe), supernova remnants (SNRs) and possibly active galactic nuclei (AGNs). Using new and archival VLBI data spanning 20 years, we obtain 23 high-resolution radio images of Arp 220 at wavelengths from 18 cm to 2 cm. From model-fitting to the images we obtain estimates of flux densities and sizes of all detected sources. We detect radio continuum emission from 97 compact sources and present flux densities and sizes for all analysed observation epochs. We find evidence for a LD-relation within Arp 220, with larger sources being less luminous. We find a compact source LF $n(L)\propto L^β$ with $β=-2.19\pm0.15$, similar to SNRs in normal galaxies. Based on simulations we argue that there are many relatively large and weak sources below our detection threshold. The observations can be explained by a mixed population of SNe and SNRs, where the former expand in a dense circumstellar medium (CSM) and the latter interact with the surrounding interstellar medium (ISM). Nine sources are likely luminous, type IIn SNe. This number of luminous SNe correspond to few percent of the total number of SNe in Arp 220 which is consistent with a total SN-rate of 4 yr$^{-1}$ as inferred from the total radio emission given a normal stellar initial mass function (IMF). Based on the fitted luminosity function, we argue that emission from all compact sources, also below our detection threshold, make up at most 20\% of the total radio emission at GHz frequencies. △ Less

Submitted 15 February, 2019; v1 submitted 15 February, 2017; originally announced February 2017.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 623, A173 (2019)

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

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

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

Submitted 24 January, 2017; originally announced January 2017.

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

The LOFAR Two-metre Sky Survey - I. Survey Description and Preliminary Data Release

Authors: T. W. Shimwell, H. J. A. Röttgering, P. N. Best, W. L. Williams, T. J. Dijkema, F. de Gasperin, M. J. Hardcastle, G. H. Heald, D. N. Hoang, A. Horneffer, H. Intema, E. K. Mahony, S. Mandal, A. P. Mechev, L. Morabito, J. B. R. Oonk, D. Rafferty, E. Retana-Montenegro, J. Sabater, C. Tasse, R. J. van Weeren, M. Brüggen, G. Brunetti, K. T. Chyży, J. E. Conway , et al. (47 additional authors not shown)

Abstract: The LOFAR Two-metre Sky Survey (LoTSS) is a deep 120-168 MHz imaging survey that will eventually cover the entire Northern sky. Each of the 3170 pointings will be observed for 8 hrs, which, at most declinations, is sufficient to produce ~5arcsec resolution images with a sensitivity of ~0.1mJy/beam and accomplish the main scientific aims of the survey which are to explore the formation and evolutio… ▽ More The LOFAR Two-metre Sky Survey (LoTSS) is a deep 120-168 MHz imaging survey that will eventually cover the entire Northern sky. Each of the 3170 pointings will be observed for 8 hrs, which, at most declinations, is sufficient to produce ~5arcsec resolution images with a sensitivity of ~0.1mJy/beam and accomplish the main scientific aims of the survey which are to explore the formation and evolution of massive black holes, galaxies, clusters of galaxies and large-scale structure. Due to the compact core and long baselines of LOFAR, the images provide excellent sensitivity to both highly extended and compact emission. For legacy value, the data are archived at high spectral and time resolution to facilitate subarcsecond imaging and spectral line studies. In this paper we provide an overview of the LoTSS. We outline the survey strategy, the observational status, the current calibration techniques, a preliminary data release, and the anticipated scientific impact. The preliminary images that we have released were created using a fully-automated but direction-independent calibration strategy and are significantly more sensitive than those produced by any existing large-area low-frequency survey. In excess of 44,000 sources are detected in the images that have a resolution of 25arcsec, typical noise levels of less than 0.5 mJy/beam, and cover an area of over 350 square degrees in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45d00m00s to 57d00m00s). △ Less

Submitted 8 November, 2016; originally announced November 2016.

Comments: 23 pages, 15 figures, 3 tables, accepted for publication in A&A

Journal ref: A&A 598, A104 (2017)

Estimating sizes of faint, distant galaxies in the submillimetre regime

Authors: L. Lindroos, K. K. Knudsen, L. Fan, J. Conway, K. Coppin, R. Decarli, G. Drouart, J. A. Hodge, A. Karim, J. M. Simpson, J. Wardlow

Abstract: We measure the sizes of redshift ~2 star-forming galaxies by stacking data from the Atacama Large Millimeter/submillimeter Array (ALMA). We use a uv-stacking algorithm in combination with model fitting in the uv-domain and show that this allows for robust measures of the sizes of marginally resolved sources. The analysis is primarily based on the 344 GHz ALMA continuum observations centred on 88 s… ▽ More We measure the sizes of redshift ~2 star-forming galaxies by stacking data from the Atacama Large Millimeter/submillimeter Array (ALMA). We use a uv-stacking algorithm in combination with model fitting in the uv-domain and show that this allows for robust measures of the sizes of marginally resolved sources. The analysis is primarily based on the 344 GHz ALMA continuum observations centred on 88 sub-millimeter galaxies in the LABOCA ECDFS Submillimeter Survey (ALESS). We study several samples of galaxies at z~2 with $M_* \sim{} 5\times{}10^{10} $ M$_\odot$ , selected using near-infrared photometry (distant red galaxies, extremely red objects, sBzK-galaxies, and galaxies selected on photometric redshift). We find that the typical sizes of these galaxies are ~0.6 arcsec which corresponds to ~5 kpc at z~2, this agrees well with the median sizes measured in the near-infrared z-band (~0.6 arcsec). We find errors on our size estimates of ~0.1-0.2 arcsec, which agree well with the expected errors for model fitting at the given signal-to-noise ratio. With the uv-coverage of our observations (18-160 m), the size and flux density measurements are sensitive to scales out to 2 . We compare this to a simulated ALMA Cycle 3 dataset with intermediate length baseline coverage, and we find that, using only these baselines, the measured stacked flux density would be an order of magnitude fainter. This highlights the importance of short baselines to recover the full flux density of high-redshift galaxies. △ Less

Submitted 31 August, 2016; originally announced August 2016.

Comments: Accepted for publication in MNRAS

LBCS: the LOFAR Long-Baseline Calibrator Survey

Authors: N. Jackson, A. Tagore, A. Deller, J. Moldón, E. Varenius, L. Morabito, O. Wucknitz, T. Carozzi, J. Conway, A. Drabent, A. Kapinska, E. Orrù, M. Brentjens, R. Blaauw, G. Kuper, J. Sluman, J. Schaap, N. Vermaas, M. Iacobelli, L. Cerrigone, A. Shulevski, S. ter Veen, R. Fallows, R. Pizzo, M. Sipior , et al. (54 additional authors not shown)

Abstract: (abridged). We outline LBCS (the LOFAR Long-Baseline Calibrator Survey), whose aim is to identify sources suitable for calibrating the highest-resolution observations made with the International LOFAR Telescope, which include baselines >1000 km. Suitable sources must contain significant correlated flux density (50-100mJy) at frequencies around 110--190~MHz on scales of a few hundred mas. At least… ▽ More (abridged). We outline LBCS (the LOFAR Long-Baseline Calibrator Survey), whose aim is to identify sources suitable for calibrating the highest-resolution observations made with the International LOFAR Telescope, which include baselines >1000 km. Suitable sources must contain significant correlated flux density (50-100mJy) at frequencies around 110--190~MHz on scales of a few hundred mas. At least for the 200--300-km international baselines, we find around 1 suitable calibrator source per square degree over a large part of the northern sky, in agreement with previous work. This should allow a randomly selected target to be successfully phase calibrated on the international baselines in over 50% of cases. Products of the survey include calibrator source lists and fringe-rate and delay maps of wide areas -- typically a few degrees -- around each source. The density of sources with significant correlated flux declines noticeably with baseline length over the range 200--600~km, with good calibrators on the longest baselines appearing only at the rate of 0.5 per square degree. Coherence times decrease from 1--3 minutes on 200-km baselines to about 1 minute on 600-km baselines, suggesting that ionospheric phase variations contain components with scales of a few hundred kilometres. The longest median coherence time, at just over 3 minutes, is seen on the DE609 baseline, which at 227km is close to being the shortest. We see median coherence times of between 80 and 110 seconds on the four longest baselines (580--600~km), and about 2 minutes for the other baselines. The success of phase transfer from calibrator to target is shown to be influenced by distance, in a manner that suggests a coherence patch at 150-MHz of the order of 1 degree. △ Less

Submitted 9 September, 2016; v1 submitted 6 August, 2016; originally announced August 2016.

Comments: Accepted by Astronomy & Astrophysics. Error in figure 6 corrected

Journal ref: A&A 595, A86 (2016)

Detection of 183 GHz water megamaser emission towards NGC 4945

Authors: Elizabeth Humphreys, Wouter Vlemmings, Violette Impellizzeri, Maud Galametz, Michael Olberg, John Conway, Victor Belitsky, Carlos De Breuck

Abstract: Aim: The aim of this work is to search Seyfert 2 galaxy NGC 4945, a well-known 22 GHz water megamaser galaxy, for water (mega)maser emission at 183 GHz. Method: We used APEX SEPIA Band 5 to perform the observations. Results: We detected 183 GHz water maser emission towards NGC 4945 with a peak flux density of ~3 Jy near the galactic systemic velocity. The emission spans a velocity range of several… ▽ More Aim: The aim of this work is to search Seyfert 2 galaxy NGC 4945, a well-known 22 GHz water megamaser galaxy, for water (mega)maser emission at 183 GHz. Method: We used APEX SEPIA Band 5 to perform the observations. Results: We detected 183 GHz water maser emission towards NGC 4945 with a peak flux density of ~3 Jy near the galactic systemic velocity. The emission spans a velocity range of several hundred km/s. We estimate an isotropic luminosity of > 1000 Lsun, classifying the emission as a megamaser. A comparison of the 183 GHz spectrum with that observed at 22 GHz suggests that 183 GHz emission also arises from the active galactic nucleus (AGN) central engine. If the 183 GHz emission originates from the circumnuclear disk, then we estimate that a redshifted feature at 1084 km/s in the spectrum should arise from a distance of 0.022 pc from the supermassive black hole (1.6 x 10(5) Schwarzschild radii), i.e. closer than the water maser emission previously detected at 22 GHz. This is only the second time 183 GHz maser emission has been detected towards an AGN central engine (the other galaxy being NGC 3079). It is also the strongest extragalactic millimetre/submillimetre water maser detected to date. Conclusions: Strong millimetre 183 GHz water maser emission has now been shown to occur in an external galaxy. For NGC 4945, we believe that the maser emission arises, or is dominated by, emission from the AGN central engine. Emission at higher velocity, i.e. for a Keplerian disk closer to the black hole, has been detected at 183 GHz compared with that for the 22 GHz megamaser. This indicates that millimetre/submillimetre water masers can indeed be useful probes for tracing out more of AGN central engine structures and dynamics than previously probed. Future observations using ALMA Band 5 should unequivocally determine the origin of the emission in this and other galaxies. △ Less

Submitted 31 July, 2016; originally announced August 2016.

Comments: 4 pages, accepted by A&A Letters

Journal ref: A&A 592, L13 (2016)

Subarcsecond international LOFAR radio images of Arp 220 at 150 MHz: A kpc-scale star forming disk surrounding nuclei with shocked outflows

Authors: E. Varenius, J. E. Conway, I. Martí-Vidal, S. Aalto, L. Barcos-Muñoz, S. König, M. A. Pérez-Torres, A. T. Deller, J. Moldón, J. S. Gallagher, T. M. Yoast-Hull, C. Horellou, L. K. Morabito, A. Alberdi, N. Jackson, R. Beswick, T. D. Carozzi, O. Wucknitz, N. Ramírez-Olivencia

Abstract: We analyse new observations with the International Low Frequency Array (LOFAR) telescope, and archival data from the Multi-Element Radio Linked Interferometer Network (MERLIN) and the Karl G. Jansky Very Large Array (VLA). We model the spatially resolved radio spectrum of Arp 220 from 150 MHz to 33 GHz. We present an image of Arp 220 at 150 MHz with resolution $0.65''\times0.35''$, sensitivity 0.1… ▽ More We analyse new observations with the International Low Frequency Array (LOFAR) telescope, and archival data from the Multi-Element Radio Linked Interferometer Network (MERLIN) and the Karl G. Jansky Very Large Array (VLA). We model the spatially resolved radio spectrum of Arp 220 from 150 MHz to 33 GHz. We present an image of Arp 220 at 150 MHz with resolution $0.65''\times0.35''$, sensitivity 0.15 mJy beam$^{-1}$, and integrated flux density $394\pm59$ mJy. More than 80% of the detected flux comes from extended ($6''\approx$2.2 kpc) steep spectrum ($α=-0.7$) emission, likely from star formation in the molecular disk surrounding the two nuclei. We find elongated features extending $0.3''$ (110 pc) and $0.9''$ (330 pc) from the eastern and western nucleus respectively, which we interpret as evidence for outflows. The extent of radio emission requires acceleration of cosmic rays far outside the nuclei. We find that a simple three component model can explain most of the observed radio spectrum of the galaxy. When accounting for absorption at 1.4 GHz, Arp 220 follows the FIR/radio correlation with $q=2.36$, and we estimate a star formation rate of 220 M$_\odot\text{yr}^{-1}$. We derive thermal fractions at 1 GHz of less than 1% for the nuclei, which indicates that a major part of the UV-photons are absorbed by dust. International LOFAR observations shows great promise to detect steep spectrum outflows and probe regions of thermal absorption. However, in LIRGs the emission detected at 150 MHz does not necessarily come from the main regions of star formation. This implies that high spatial resolution is crucial for accurate estimates of star formation rates for such galaxies at 150 MHz. △ Less

Submitted 10 July, 2016; originally announced July 2016.

Comments: 17 pages, 9 figures

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

LOFAR/H-ATLAS: A deep low-frequency survey of the Herschel-ATLAS North Galactic Pole field

Authors: M. J. Hardcastle, G. Gürkan, R. J. van Weeren, W. L. Williams, P. N. Best, F. de Gasperin, D. A. Rafferty, S. C. Read, J. Sabater, T. W. Shimwell, D. J. B. Smith, C. Tasse, N. Bourne, M. Brienza, M. Brüggen, G. Brunetti, K. T. Chyży, J. Conway, L. Dunne, S. A. Eales, S. J. Maddox, M. J. Jarvis, E. K. Mahony, R. Morganti, I. Prandoni , et al. (3 additional authors not shown)

Abstract: We present LOFAR High-Band Array (HBA) observations of the Herschel-ATLAS North Galactic Pole survey area. The survey we have carried out, consisting of four pointings covering around 142 square degrees of sky in the frequency range 126--173 MHz, does not provide uniform noise coverage but otherwise is representative of the quality of data to be expected in the planned LOFAR wide-area surveys, and… ▽ More We present LOFAR High-Band Array (HBA) observations of the Herschel-ATLAS North Galactic Pole survey area. The survey we have carried out, consisting of four pointings covering around 142 square degrees of sky in the frequency range 126--173 MHz, does not provide uniform noise coverage but otherwise is representative of the quality of data to be expected in the planned LOFAR wide-area surveys, and has been reduced using recently developed `facet calibration' methods at a resolution approaching the full resolution of the datasets ($\sim 10 \times 6$ arcsec) and an rms off-source noise that ranges from 100 $μ$Jy beam$^{-1}$ in the centre of the best fields to around 2 mJy beam$^{-1}$ at the furthest extent of our imaging. We describe the imaging, cataloguing and source identification processes, and present some initial science results based on a 5-$σ$ source catalogue. These include (i) an initial look at the radio/far-infrared correlation at 150 MHz, showing that many Herschel sources are not yet detected by LOFAR; (ii) number counts at 150 MHz, including, for the first time, observational constraints on the numbers of star-forming galaxies; (iii) the 150-MHz luminosity functions for active and star-forming galaxies, which agree well with determinations at higher frequencies at low redshift, and show strong redshift evolution of the star-forming population; and (iv) some discussion of the implications of our observations for studies of radio galaxy life cycles. △ Less

Submitted 30 June, 2016; originally announced June 2016.

Comments: 26 pages, 22 colour figures. Resubmitted to MNRAS after referee comments

LOFAR 150-MHz observations of the Boötes field: Catalogue and Source Counts

Authors: W. L. Williams, R. J. van Weeren, H. J. A. Röttgering, P. Best, T. J. Dijkema, F. de Gasperin, M. J. Hardcastle, G. Heald, I. Prandoni, J. Sabater, T. W. Shimwell, C. Tasse, I. M. van Bemmel, M. Brüggen, G. Brunetti, J. E. Conway, T. Enßlin, D. Engels, H. Falcke, C. Ferrari, M. Haverkorn, N. Jackson, M. J. Jarvis, A. D. Kapinska, E. K. Mahony , et al. (10 additional authors not shown)

Abstract: We present the first wide area (19 deg$^2$), deep ($\approx120-150$ μJy beam$^{-1}$), high resolution ($5.6 \times 7.4$ arcsec) LOFAR High Band Antenna image of the Boötes field made at 130-169 MHz. This image is at least an order of magnitude deeper and 3-5 times higher in angular resolution than previously achieved for this field at low frequencies. The observations and data reduction, which inc… ▽ More We present the first wide area (19 deg$^2$), deep ($\approx120-150$ μJy beam$^{-1}$), high resolution ($5.6 \times 7.4$ arcsec) LOFAR High Band Antenna image of the Boötes field made at 130-169 MHz. This image is at least an order of magnitude deeper and 3-5 times higher in angular resolution than previously achieved for this field at low frequencies. The observations and data reduction, which includes full direction-dependent calibration, are described here. We present a radio source catalogue containing 6276 sources detected over an area of $19$\,deg$^2$, with a peak flux density threshold of $5σ$. As the first thorough test of the facet calibration strategy, introduced by van Weeren et al., we investigate the flux and positional accuracy of the catalogue. We present differential source counts that reach an order of magnitude deeper in flux density than previously achieved at these low frequencies, and show flattening at 150 MHz flux densities below 10 mJy associated with the rise of the low flux density star-forming galaxies and radio-quiet AGN. △ Less

Submitted 5 May, 2016; originally announced May 2016.

Comments: 29 pages, 27 figures, accepted for publication in MNRAS

A large light-mass component of cosmic rays at 10^{17} - 10^{17.5} eV from radio observations

Authors: S. Buitink, A. Corstanje, H. Falcke, J. R. Hörandel, T. Huege, A. Nelles, J. P. Rachen, L. Rossetto, P . Schellart, O. Scholten, S. ter Veen, S. Thoudam, T. N. G. Trinh, J. Anderson, A. Asgekar, I. M. Avruch, M. E. Bell, M. J. Bentum, G. Bernardi, P. Best, A. Bonafede, F. Breitling, J. W. Broderick, W. N. Brouw, M. Brüggen , et al. (79 additional authors not shown)

Abstract: Cosmic rays are the highest energy particles found in nature. Measurements of the mass composition of cosmic rays between 10^{17} eV and 10^{18} eV are essential to understand whether this energy range is dominated by Galactic or extragalactic sources. It has also been proposed that the astrophysical neutrino signal comes from accelerators capable of producing cosmic rays of these energies. Cosmic… ▽ More Cosmic rays are the highest energy particles found in nature. Measurements of the mass composition of cosmic rays between 10^{17} eV and 10^{18} eV are essential to understand whether this energy range is dominated by Galactic or extragalactic sources. It has also been proposed that the astrophysical neutrino signal comes from accelerators capable of producing cosmic rays of these energies. Cosmic rays initiate cascades of secondary particles (air showers) in the atmosphere and their masses are inferred from measurements of the atmospheric depth of the shower maximum, Xmax, or the composition of shower particles reaching the ground. Current measurements suffer from either low precision, or a low duty cycle and a high energy threshold. Radio detection of cosmic rays is a rapidly developing technique, suitable for determination of Xmax with a duty cycle of in principle nearly 100%. The radiation is generated by the separation of relativistic charged particles in the geomagnetic field and a negative charge excess in the shower front. Here we report radio measurements of Xmax with a mean precision of 16 g/cm^2 between 10^{17}-10^{17.5} eV. Because of the high resolution in $Xmax we can determine the mass spectrum and find a mixed composition, containing a light mass fraction of ~80%. Unless the extragalactic component becomes significant already below 10^{17.5} eV, our measurements indicate an additional Galactic component dominating at this energy range. △ Less

Submitted 1 May, 2016; v1 submitted 4 March, 2016; originally announced March 2016.

Comments: 35 pages, 11 figures, updated version: Pierre Auger Observatory data ICRC 2015 added to Fig 2

Journal ref: Nature 531, 70 (2016)