Showing 1–22 of 22 results for author: Kettenis, M

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

Reference Array and Design Consideration for the next-generation Event Horizon Telescope

Authors: Sheperd S. Doeleman, John Barrett, Lindy Blackburn, Katherine Bouman, Avery E. Broderick, Ryan Chaves, Vincent L. Fish, Garret Fitzpatrick, Antonio Fuentes, Mark Freeman, José L. Gómez, Kari Haworth, Janice Houston, Sara Issaoun, Michael D. Johnson, Mark Kettenis, Laurent Loinard, Neil Nagar, Gopal Narayanan, Aaron Oppenheimer, Daniel C. M. Palumbo, Nimesh Patel, Dominic W. Pesce, Alexander W. Raymond, Freek Roelofs , et al. (4 additional authors not shown)

Abstract: We describe the process to design, architect, and implement a transformative enhancement of the Event Horizon Telescope (ngEHT). This program - the next-generation Event Horizon Telescope (ngEHT) - will form a networked global array of radio dishes capable of making high-fidelity real-time movies of supermassive black holes (SMBH) and their emanating jets. This builds upon the EHT principally by d… ▽ More We describe the process to design, architect, and implement a transformative enhancement of the Event Horizon Telescope (ngEHT). This program - the next-generation Event Horizon Telescope (ngEHT) - will form a networked global array of radio dishes capable of making high-fidelity real-time movies of supermassive black holes (SMBH) and their emanating jets. This builds upon the EHT principally by deploying additional modest-diameter dishes to optimized geographic locations to enhance the current global mm/submm wavelength Very Long Baseline Interferometric (VLBI) array, which has, to date, utilized mostly pre-existing radio telescopes. The ngEHT program further focuses on observing at three frequencies simultaneously for increased sensitivity and Fourier spatial frequency coverage. Here, the concept, science goals, design considerations, station siting and instrument prototyping are discussed, and a preliminary reference array to be implemented in phases is described. △ Less

Submitted 17 August, 2023; v1 submitted 14 June, 2023; originally announced June 2023.

Comments: Submitted to the journal Galaxies

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

CASA, the Common Astronomy Software Applications for Radio Astronomy

Authors: THE CASA TEAM, Ben Bean, Sanjay Bhatnagar, Sandra Castro, Jennifer Donovan Meyer, Bjorn Emonts, Enrique Garcia, Robert Garwood, Kumar Golap, Justo Gonzalez Villalba, Pamela Harris, Yohei Hayashi, Josh Hoskins, Mingyu Hsieh, Preshanth Jagannathan, Wataru Kawasaki, Aard Keimpema, Mark Kettenis, Jorge Lopez, Joshua Marvil, Joseph Masters, Andrew McNichols, David Mehringer, Renaud Miel, George Moellenbrock , et al. (24 additional authors not shown)

Abstract: CASA, the Common Astronomy Software Applications, is the primary data processing software for the Atacama Large Millimeter/submillimeter Array (ALMA) and the Karl G. Jansky Very Large Array (VLA), and is frequently used also for other radio telescopes. The CASA software can handle data from single-dish, aperture-synthesis, and Very Long Baseline Interferometery (VLBI) telescopes. One of its core f… ▽ More CASA, the Common Astronomy Software Applications, is the primary data processing software for the Atacama Large Millimeter/submillimeter Array (ALMA) and the Karl G. Jansky Very Large Array (VLA), and is frequently used also for other radio telescopes. The CASA software can handle data from single-dish, aperture-synthesis, and Very Long Baseline Interferometery (VLBI) telescopes. One of its core functionalities is to support the calibration and imaging pipelines for ALMA, VLA, VLA Sky Survey (VLASS), and the Nobeyama 45m telescope. This paper presents a high-level overview of the basic structure of the CASA software, as well as procedures for calibrating and imaging astronomical radio data in CASA. CASA is being developed by an international consortium of scientists and software engineers based at the National Radio Astronomical Observatory (NRAO), the European Southern Observatory (ESO), the National Astronomical Observatory of Japan (NAOJ), and the Joint Institute for VLBI European Research Infrastructure Consortium (JIV-ERIC), under the guidance of NRAO. △ Less

Submitted 5 October, 2022; originally announced October 2022.

Comments: Accepted for publication in PASP (20 pages, 4 figures). Joint publication with CASA-VLBI paper

CASA on the fringe -- Development of VLBI processing capabilities for CASA

Authors: Ilse M. van Bemmel, Mark Kettenis, Des Small, Michael Janssen, George A. Moellenbrock, Dirk Petry, Ciriaco Goddi, Justin D. Linford, Kazi L. J. Rygl, Elisabetta Liuzzo, Benito Marcote, Olga S. Bayandina, Neal Schweigart, Marjolein Verkouter, Aard Keimpema, Arpad Szomoru, Huib Jan van Langevelde

Abstract: New functionality to process Very Long Baseline Interferometry (VLBI) data has been implemented in the CASA package. This includes two new tasks to handle fringe fitting and VLBI-specific amplitude calibration steps. Existing tasks have been adjusted to handle VLBI visibility data and calibration meta-data properly. With these updates, it is now possible to process VLBI continuum and spectral line… ▽ More New functionality to process Very Long Baseline Interferometry (VLBI) data has been implemented in the CASA package. This includes two new tasks to handle fringe fitting and VLBI-specific amplitude calibration steps. Existing tasks have been adjusted to handle VLBI visibility data and calibration meta-data properly. With these updates, it is now possible to process VLBI continuum and spectral line observations in CASA. This article describes the development and implementation, and presents an outline for the workflow when calibrating European VLBI Network or Very Long Baseline Array data in CASA. Though the CASA VLBI functionality has already been vetted extensively as part of the Event Horizon Telescope data processing, in this paper we compare results for the same dataset processed in CASA and AIPS. We find identical results for the two packages and conclude that CASA in some cases performs better, though it cannot match AIPS for single-core processing time. The new functionality in CASA allows for easy development of pipelines or Jupyter notebooks, and thus contributes to raising VLBI data processing to present day standards for accessibility, reproducibility, and reusability. △ Less

Submitted 5 October, 2022; originally announced October 2022.

Comments: PASP in press, joint submission with the CASA reference paper; 12 pages, 5 figures

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

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)

Radio data archives round table

Authors: K. A. Lutz, J. Dempsey, Y. G. Grange, M. Kettenis, M. Lacy, C. Schollar

Abstract: With SKA precursor and pathfinder operations in full swing, radio and (sub-)mm astronomy is entering the era of super big data. The big questions is how to make (sub-)mm and radio data available to the astronomical community, preferably using FAIR (findable, accessible, interoperable and re-useable; Wilkinson et al. 2016) principles. There are already a lot of efforts going on around the globe: fa… ▽ More With SKA precursor and pathfinder operations in full swing, radio and (sub-)mm astronomy is entering the era of super big data. The big questions is how to make (sub-)mm and radio data available to the astronomical community, preferably using FAIR (findable, accessible, interoperable and re-useable; Wilkinson et al. 2016) principles. There are already a lot of efforts going on around the globe: facilities such as ALMA, LOFAR, MWA, NRAO and ASKAP are already publishing much of their imaging data in the form of "science ready" products, SKA regional centres are being formed and a radio astronomy interest group has been initiated within the IVOA. In addition, also non-imaging data, like timing, pulsar or beam forming data, needs to be available to the community as well. This BoF intended to bring everyone interested in this topic around one virtual table to hear about and discuss the following questions: What is the status of efforts to expose both visibility and science ready data? What is already there, maybe has been used for decades by traditional observatories? What is still missing? Where do we want to go next? △ Less

Submitted 22 December, 2020; originally announced December 2020.

Comments: 4 pages, to appear in the proceedings of Astronomical Data Analysis Software and Systems XXX published by ASP

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

A probabilistic approach to phase calibration: I. Effects of source structure on fringe-fitting

Authors: Iniyan Natarajan, Roger Deane, Ilse van Bemmel, Huib Jan van Langevelde, Des Small, Mark Kettenis, Zsolt Paragi, Oleg Smirnov, Arpad Szomoru

Abstract: We propose a probabilistic framework for performing simultaneous estimation of source structure and fringe-fitting parameters in Very Long Baseline Interferometry (VLBI) observations. As a first step, we demonstrate this technique through the analysis of synthetic short-duration Event Horizon Telescope (EHT) observations of various geometric source models at 230 GHz, in the presence of baseline-de… ▽ More We propose a probabilistic framework for performing simultaneous estimation of source structure and fringe-fitting parameters in Very Long Baseline Interferometry (VLBI) observations. As a first step, we demonstrate this technique through the analysis of synthetic short-duration Event Horizon Telescope (EHT) observations of various geometric source models at 230 GHz, in the presence of baseline-dependent thermal noise. We perform Bayesian parameter estimation and model selection between the different source models to obtain reliable uncertainty estimates and correlations between various source and fringe-fitting related model parameters. We also compare the Bayesian posteriors with those obtained using widely-used VLBI data reduction packages such as CASA and AIPS, by fringe-fitting 200 Monte Carlo simulations of each source model with different noise realisations, to obtain distributions of the Maximum A Posteriori (MAP) estimates. We find that, in the presence of resolved asymmetric source structure and a given array geometry, the traditional practice of fringe-fitting with a point source model yields appreciable offsets in the estimated phase residuals, potentially biasing or limiting the dynamic range of the starting model used for self-calibration. Simultaneously estimating the source structure earlier in the calibration process with formal uncertainties improves the precision and accuracy of fringe-fitting and establishes the potential of the available data especially when there is little prior information. We also note the potential applications of this method to astrometry and geodesy for specific science cases and the planned improvements to the computational performance and analyses of more complex source distributions. △ Less

Submitted 26 May, 2020; originally announced May 2020.

Comments: accepted for publication in MNRAS

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

rPICARD: A CASA-based Calibration Pipeline for VLBI Data

Authors: M. Janssen, C. Goddi, H. Falcke, D. van Rossum, I. van Bemmel, M. Kettenis, D. Small, I. Martí-Vidal

Abstract: Currently, HOPS and AIPS are the primary choices for the time-consuming process of (millimeter) Very Long Baseline Interferometry (VLBI) data calibration. However, for a full end-to-end pipeline, they either lack the ability to perform easily scriptable incremental calibration or do not provide full control over the workflow with the ability to manipulate and edit calibration solutions directly. T… ▽ More Currently, HOPS and AIPS are the primary choices for the time-consuming process of (millimeter) Very Long Baseline Interferometry (VLBI) data calibration. However, for a full end-to-end pipeline, they either lack the ability to perform easily scriptable incremental calibration or do not provide full control over the workflow with the ability to manipulate and edit calibration solutions directly. The Common Astronomy Software Application (CASA) offers all these abilities, together with a secure development future and an intuitive Python interface, which is very attractive for young radio astronomers. Inspired by the recent addition of a global fringe-fitter, the capability to convert FITS-IDI files to measurement sets, and amplitude calibration routines based on ANTAB metadata, we have developed the the CASA-based Radboud PIpeline for the Calibration of high Angular Resolution Data (rPICARD). The pipeline will be able to handle data from multiple arrays: EHT, GMVA, VLBA and the EVN in the first release. Polarization and phase-referencing calibration are supported and a spectral line mode will be added in the future. The large bandwidths of future radio observatories ask for a scalable reduction software. Within CASA, a message passing interface (MPI) implementation is used for parallelization, reducing the total time needed for processing. The most significant gain is obtained for the time-consuming fringe-fitting task where each scan be processed in parallel. △ Less

Submitted 6 May, 2019; originally announced May 2019.

Comments: 6 pages, 1 figure, EVN 2018 symposium proceedings

Journal ref: PoS(EVN2018)080

CASA on the fringe: VLBI data processing in the CASA software package

Authors: Ilse van Bemmel, Des Small, Mark Kettenis, Arpad Szomoru, George Moellenbrock, Michael Janssen

Abstract: In recent years new functionality for VLBI data processing has been added to the CASA package. This paper presents the new CASA tasks 'fringefit' and 'accor', which are closely matched to their AIPS counterparts FRING and ACCOR. Several CASA tasks received upgrades to handle VLBI specific metadata. With the current CASA release VLBI data processing is possible, and functionality will be expanded i… ▽ More In recent years new functionality for VLBI data processing has been added to the CASA package. This paper presents the new CASA tasks 'fringefit' and 'accor', which are closely matched to their AIPS counterparts FRING and ACCOR. Several CASA tasks received upgrades to handle VLBI specific metadata. With the current CASA release VLBI data processing is possible, and functionality will be expanded in the upcoming release. Longer term developments include fringe fitting of broad, non-continuous frequency bands and dispersive delays, which will ensure that the number of use cases for VLBI calibration will increase in future CASA releases. △ Less

Submitted 26 April, 2019; originally announced April 2019.

Comments: 5 pages, 2 figures, EVN 2018 symposium proceedings

Journal ref: PoS(EVN2018)079

rPICARD: A CASA-based Calibration Pipeline for VLBI Data. Calibration and imaging of 7 mm VLBA observations of the AGN jet in M87

Authors: Michael Janssen, Ciriaco Goddi, Ilse M. van Bemmel, Mark Kettenis, Des Small, Elisabetta Liuzzo, Kazi Rygl, Ivan Martí-Vidal, Lindy Blackburn, Maciek Wielgus, Heino Falcke

Abstract: (Abridged) The CASA software suite, can now reduce very long baseline interferometry (VLBI) data with the recent addition of a fringe fitter. Here, we present the Radboud PIpeline for the Calibration of high Angular Resolution Data (rPICARD), which is an open-source VLBI calibration and imaging pipeline built on top of the CASA framework. The pipeline is capable of reducing data from different VLB… ▽ More (Abridged) The CASA software suite, can now reduce very long baseline interferometry (VLBI) data with the recent addition of a fringe fitter. Here, we present the Radboud PIpeline for the Calibration of high Angular Resolution Data (rPICARD), which is an open-source VLBI calibration and imaging pipeline built on top of the CASA framework. The pipeline is capable of reducing data from different VLBI arrays. It can be run non-interactively after only a few non-default input parameters are set and delivers high-quality calibrated data. CPU scalability based on a message-passing interface (MPI) implementation ensures that large bandwidth data from future arrays can be processed within reasonable computing times. Phase calibration is done with a Schwab-Cotton fringe fit algorithm. For the calibration of residual atmospheric effects, optimal solution intervals are determined based on the signal-to-noise ratio (S/N) of the data for each scan. Different solution intervals can be set for different antennas in the same scan to increase the number of detections in the low S/N regime. These novel techniques allow rPICARD to calibrate data from different arrays, including high-frequency and low-sensitivity arrays. The amplitude calibration is based on standard telescope metadata, and a robust algorithm can solve for atmospheric opacity attenuation in the high-frequency regime. Standard CASA tasks are used for CLEAN imaging and self-calibration. In this work we demonstrate the capabilities of rPICARD by calibrating and imaging 7 mm VLBA data of the central radio source in the M87 galaxy. The reconstructed jet image reveals a complex collimation profile and edge-brightened structure. A potential counter-jet is detected that has 10 % of the brightness of the approaching jet. This constrains jet speeds close to the radio core to about half the speed of light for small inclination angles. △ Less

Submitted 15 May, 2019; v1 submitted 5 February, 2019; originally announced February 2019.

Comments: Accepted for publication Astronomy and Astrophysics. Code available at https://bitbucket.org/M_Janssen/picard

Journal ref: A&A 626, A75 (2019)

Minimal Re-computation for Exploratory Data Analysis in Astronomy

Authors: Bojan Nikolic, Des Small, Mark Kettenis

Abstract: We present a technique to automatically minimise the re-computation when a data analysis program is iteratively changed, or added to, as is often the case in exploratory data analysis in astronomy. A typical example is flagging and calibration of demanding or unusual observations where visual inspection suggests improvement to the processing strategy. The technique is based on memoization and refe… ▽ More We present a technique to automatically minimise the re-computation when a data analysis program is iteratively changed, or added to, as is often the case in exploratory data analysis in astronomy. A typical example is flagging and calibration of demanding or unusual observations where visual inspection suggests improvement to the processing strategy. The technique is based on memoization and referentially transparent tasks. We describe the implementation of this technique for the CASA radio astronomy data reduction package. We also propose a technique for optimising efficiency of storage of memoized intermediate data products using copy-on-write and block level de-duplication and measure their practical efficiency. We find the minimal recomputation technique improves the efficiency of data analysis while reducing the possibility for user error and improving the reproducibility of the final result. It also aids exploratory data analysis on batch-schedule cluster computer systems. △ Less

Submitted 6 September, 2018; originally announced September 2018.

Comments: Accepted for publication in Astronomy and Computing. Software described in the paper is available at: http://www.mrao.cam.ac.uk/~bn204/sw/recipe . This refereed paper is an expanded version of the conference proceedings arxiv:1711.06124

Minimal Re-computation for Exploratory Data Analysis in Astronomy

Authors: Bojan Nikolic, Des Small, Mark Kettenis

Abstract: We present a technique to automatically minimise the re-computation when a data processing program is iteratively changed, or added to, as is often the case in exploratory data analysis in radio astronomy. A typical example is flagging and calibration of demanding or unusual observations where visual inspection suggests improvement to the processing strategy. The technique is based on memoization… ▽ More We present a technique to automatically minimise the re-computation when a data processing program is iteratively changed, or added to, as is often the case in exploratory data analysis in radio astronomy. A typical example is flagging and calibration of demanding or unusual observations where visual inspection suggests improvement to the processing strategy. The technique is based on memoization and referentially transparent tasks. We describe a prototype implementation for the CASA data reduction package. This technique improves the efficiency of data analysis while reducing the possibility for user error and improving the reproducibility of the final result. △ Less

Submitted 16 November, 2017; originally announced November 2017.

Comments: Submitted to ASP Conference Series: Proceedings of ADASS XVII, Santiago, Chile 2017

Planetary Radio Interferometry and Doppler Experiment (PRIDE) technique: A test case of the Mars Express Phobos fly-by

Authors: Dmitry A. Duev, Sergei V. Pogrebenko, Giuseppe Cimò, Guifré Molera Calvés, Tatiana M. Bocanegra Bahamón, Leonid I. Gurvits, Mark M. Kettenis, Joseph Kania, Valeriu Tudose, Pascal Rosenblatt, Jean-Charles Marty, Valery Lainey, Pablo de Vicente, Jonathan Quick, Marisa Nickola, Alexander Neidhardt, Gerhard Kronschnabl, Christian Plötz, Rüdiger Haas, Michael Lindqvist, Andrea Orlati, Alexander V. Ipatov, Mikhail A. Kharinov, Andrey G. Mikhailov, Jim Lovell , et al. (11 additional authors not shown)

Abstract: The closest ever fly-by of the Martian moon Phobos, performed by the European Space Agency's Mars Express spacecraft, gives a unique opportunity to sharpen and test the Planetary Radio Interferometry and Doppler Experiments (PRIDE) technique in the interest of studying planet - satellite systems. The aim of this work is to demonstrate a technique of providing high precision positional and Doppler… ▽ More The closest ever fly-by of the Martian moon Phobos, performed by the European Space Agency's Mars Express spacecraft, gives a unique opportunity to sharpen and test the Planetary Radio Interferometry and Doppler Experiments (PRIDE) technique in the interest of studying planet - satellite systems. The aim of this work is to demonstrate a technique of providing high precision positional and Doppler measurements of planetary spacecraft using the Mars Express spacecraft. The technique will be used in the framework of Planetary Radio Interferometry and Doppler Experiments in various planetary missions, in particular in fly-by mode. We advanced a novel approach to spacecraft data processing using the techniques of Doppler and phase-referenced very long baseline interferometry spacecraft tracking. We achieved, on average, mHz precision (30 μm/s at a 10 seconds integration time) for radial three-way Doppler estimates and sub-nanoradian precision for lateral position measurements, which in a linear measure (at a distance of 1.4 AU) corresponds to ~50 m. △ Less

Submitted 19 June, 2016; originally announced June 2016.

Comments: 9 pages, 14 figures. Astronomy and Astrophysics, accepted on 2016/05/31

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

The SFXC software correlator for Very Long Baseline Interferometry: Algorithms and Implementation

Authors: A. Keimpema, M. M. Kettenis, S. V. Pogrebenko, R. M. Campbell, G. Cimó, D. A. Duev, B. Eldering, N. Kruithof, H. J. van Langevelde, D. Marchal, G. Molera Calvés, H. Ozdemir, Z. Paragi, Y. Pidopryhora, A. Szomoru, J. Yang

Abstract: In this paper a description is given of the SFXC software correlator, developed and maintained at the Joint Institute for VLBI in Europe (JIVE). The software is designed to run on generic Linux-based computing clusters. The correlation algorithm is explained in detail, as are some of the novel modes that software correlation has enabled, such as wide-field VLBI imaging through the use of multiple… ▽ More In this paper a description is given of the SFXC software correlator, developed and maintained at the Joint Institute for VLBI in Europe (JIVE). The software is designed to run on generic Linux-based computing clusters. The correlation algorithm is explained in detail, as are some of the novel modes that software correlation has enabled, such as wide-field VLBI imaging through the use of multiple phase centres and pulsar gating and binning. This is followed by an overview of the software architecture. Finally, the performance of the correlator as a function of number of CPU cores, telescopes and spectral channels is shown. △ Less

Submitted 2 February, 2015; originally announced February 2015.

Comments: Accepted by Experimental Astronomy

Using VO tools to investigate distant radio starbursts hosting obscured AGN in the HDF(N) region

Authors: A. M. S. Richards, T. W. B. Muxlow, R. Beswick, M. G. Allen, K. Benson, R. C. Dickson, M. A. Garrett, S. T. Garrington, E. Gonzalez-Solarez, P. A. Harrison, A. J. Holloway, M. M. Kettenis, R. A. Laing, E. A. Richards, H. Thrall, H. J. van Langevelde, N. A. Walton, P. N. Wilkinson, N. Winstanley, .

Abstract: A 10-arcmin field around the HDF(N) contains 92 radio sources >40 uJy, resolved by MERLIN+VLA at 0".2-2".0 resolution. 55 have Chandra X-ray counterparts including 18 with a hard X-ray photon index and high luminosity characteristic of a type-II (obscured) AGN. >70% of the radio sources have been classified as starbursts or AGN using radio morphologies, spectral indices and comparisons with opti… ▽ More A 10-arcmin field around the HDF(N) contains 92 radio sources >40 uJy, resolved by MERLIN+VLA at 0".2-2".0 resolution. 55 have Chandra X-ray counterparts including 18 with a hard X-ray photon index and high luminosity characteristic of a type-II (obscured) AGN. >70% of the radio sources have been classified as starbursts or AGN using radio morphologies, spectral indices and comparisons with optical appearance and MIR emission. Starbursts outnumber radio AGN 3:1. This study extends the VO methods previously used to identify X-ray-selected obscured type-II AGN to investigate whether very luminous radio and X-ray emission originates from different phenomena in the same galaxy. The high-redshift starbursts have typical sizes of 5--10 kpc and star formation rates of ~1000 Msun/yr. There is no correlation between radio and X-ray luminosities nor spectral indices at z>~1.3. ~70% of both the radio-selected AGN and the starburst samples were detected by Chandra. The X-ray luminosity indicates the presence of an AGN in at least half of the 45 cross-matched radio starbursts, of which 11 are type-II AGN including 7 at z>1.5. This distribution overlaps closely with the X-ray detected radio sources which were also detected by SCUBA. Stacked 1.4-GHz emission at the positions of radio-faint X-ray sources is correlated with X-ray hardness. Most extended radio starbursts at z>1.3 host X-ray selected obscured AGN. Radio emission from most of these ultra-luminous objects is dominated by star formation but it contributes less than 1/3 of their X-ray luminosity. Our results support the inferences from SCUBA and IR data, that at z>1.5, star formation is an order of magnitude more extended and more copious, it is closely linked to AGN activity and it is triggered differently, compared with star formation at lower redshifts. △ Less

Submitted 26 June, 2007; originally announced June 2007.

Comments: 24 pages, 12 figures, uses graphicx, rotating, natbib, supertabular packages and aa.cls. Accepted for publication in A&A

High-z radio starbursts host obscured X-ray AGN

Authors: A. M. S. Richards, R. J. Beswick, S. T. Garrington, T. W. B. Muxlow, H. Thrall, M. A. Garrett, M. Kettenis, H. J. van Langevelde, E. Gonzalez-Solarez, N. A. Walton, M. G. Allen

Abstract: We use Virtual Observatory methods to investigate the association between radio and X-ray emission at high redshifts. Fifty-five of the 92 HDF(N) sources resolved by combining MERLIN+VLA data were detected by Chandra, of which 18 are hard enough and bright enough to be obscured AGN. The high-z population of microJy radio sources is dominated by starbursts an order of magnitude more active and m… ▽ More We use Virtual Observatory methods to investigate the association between radio and X-ray emission at high redshifts. Fifty-five of the 92 HDF(N) sources resolved by combining MERLIN+VLA data were detected by Chandra, of which 18 are hard enough and bright enough to be obscured AGN. The high-z population of microJy radio sources is dominated by starbursts an order of magnitude more active and more extended than any found at z<1 and at least a quarter of these simultaneously host highly X-ray-luminous obscured AGN. △ Less

Submitted 4 December, 2006; originally announced December 2006.

Comments: 4 pages, 2 figures, To appear in the proceedings of 'At the Edge of the Universe' (9-13 October 2006, Sintra, Portugal)