-
Atmospheric limitations for high-frequency ground-based VLBI
Authors:
Dominic W. Pesce,
Lindy Blackburn,
Ryan Chaves,
Sheperd S. Doeleman,
Mark Freeman,
Sara Issaoun,
Michael D. Johnson,
Greg Lindahl,
Iniyan Natarajan,
Scott N. Paine,
Daniel C. M. Palumbo,
Freek Roelofs,
Paul Tiede
Abstract:
Very long baseline interferometry (VLBI) provides the highest-resolution images in astronomy. The sharpest resolution is nominally achieved at the highest frequencies, but as the observing frequency increases so too does the atmospheric contribution to the system noise, degrading the sensitivity of the array and hampering detection. In this paper, we explore the limits of high-frequency VLBI obser…
▽ More
Very long baseline interferometry (VLBI) provides the highest-resolution images in astronomy. The sharpest resolution is nominally achieved at the highest frequencies, but as the observing frequency increases so too does the atmospheric contribution to the system noise, degrading the sensitivity of the array and hampering detection. In this paper, we explore the limits of high-frequency VLBI observations using ngehtsim, a new tool for generating realistic synthetic data. ngehtsim uses detailed historical atmospheric models to simulate observing conditions, and it employs heuristic visibility detection criteria that emulate single- and multi-frequency VLBI calibration strategies. We demonstrate the fidelity of ngehtsim's predictions using a comparison with existing 230 GHz data taken by the Event Horizon Telescope (EHT), and we simulate the expected performance of EHT observations at 345 GHz. Though the EHT achieves a nearly 100% detection rate at 230 GHz, our simulations indicate that it should expect substantially poorer performance at 345 GHz; in particular, observations of M87 at 345 GHz are predicted to achieve detection rates of $\lesssim$20% that may preclude imaging. Increasing the array sensitivity through wider bandwidths and/or longer integration times -- as enabled through, e.g., the simultaneous multi-frequency upgrades envisioned for the next-generation EHT -- can improve the 345 GHz prospects and yield detection levels that are comparable to those at 230 GHz. M87 and Sgr A* observations carried out in the atmospheric window around 460 GHz could expect to regularly achieve multiple detections on long baselines, but analogous observations at 690 and 875 GHz consistently obtain almost no detections at all.
△ Less
Submitted 1 April, 2024;
originally announced April 2024.
-
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.
-
Control of a pendulum system: From simulation to reality
Authors:
Iyer Venkataraman Natarajan
Abstract:
Control theory deals with the study of controlling dynamical systems. Robots today are growing increasingly complex and moving out of factory floors to real world environment. These robots have to interact with real world environment factors such as disturbances and this requires the robot to have a control system that is robust. Testing control algorithms on robots in real world environment can p…
▽ More
Control theory deals with the study of controlling dynamical systems. Robots today are growing increasingly complex and moving out of factory floors to real world environment. These robots have to interact with real world environment factors such as disturbances and this requires the robot to have a control system that is robust. Testing control algorithms on robots in real world environment can pose critical safety issues and can be financially expensive. This has resulted in a heavy emphasis on using simulation to test control algorithms before deploying them in real world environments. Designing control algorithms is an iterative process that starts with modelling the target system in simulation, designing a controller, testing the controller in simulation and then changing the controller parameters to design a better controller. This report explores how an approximated system model of a target hardware system can be developed, which can then be used to design a LQR controller for the target system. The controller is then tested under a disturbance, on hardware and in simulation, and the system response is recorded. The system response from hardware and simulation are then compared to validate the use of approximated system models in simulation for designing and testing control algorithms.
△ Less
Submitted 7 February, 2024; v1 submitted 8 December, 2023;
originally announced December 2023.
-
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.
-
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.
-
MeqSilhouette v2: Spectrally-resolved polarimetric synthetic data generation for the Event Horizon Telescope
Authors:
Iniyan Natarajan,
Roger Deane,
Iván Martí-Vidal,
Freek Roelofs,
Michael Janssen,
Maciek Wielgus,
Lindy Blackburn,
Tariq Blecher,
Simon Perkins,
Oleg Smirnov,
Jordy Davelaar,
Monika Moscibrodzka,
Andrew Chael,
Katherine L. Bouman,
Jae-Young Kim,
Gianni Bernardi,
Ilse van Bemmel,
Heino Falcke,
Feryal Özel,
Dimitrios Psaltis
Abstract:
We present MeqSilhouette v2.0 (MeqSv2), a fully polarimetric, time-and frequency-resolved synthetic data generation software for simulating millimetre (mm) wavelength very long baseline interferometry (VLBI) observations with heterogeneous arrays. Synthetic data are a critical component in understanding real observations, testing calibration and imaging algorithms, and predicting performance metri…
▽ More
We present MeqSilhouette v2.0 (MeqSv2), a fully polarimetric, time-and frequency-resolved synthetic data generation software for simulating millimetre (mm) wavelength very long baseline interferometry (VLBI) observations with heterogeneous arrays. Synthetic data are a critical component in understanding real observations, testing calibration and imaging algorithms, and predicting performance metrics of existing or proposed sites. MeqSv2 applies physics-based instrumental and atmospheric signal corruptions constrained by empirically-derived site and station parameters to the data. The new version is capable of applying instrumental polarization effects and various other spectrally-resolved effects using the Radio Interferometry Measurement Equation (RIME) formalism and produces synthetic data compatible with calibration pipelines designed to process real data. We demonstrate the various corruption capabilities of MeqSv2 using different arrays, with a focus on the effect of complex bandpass gains on closure quantities for the EHT at 230 GHz. We validate the frequency-dependent polarization leakage implementation by performing polarization self-calibration of synthetic EHT data using PolSolve. We also note the potential applications for cm-wavelength VLBI array analysis and design and future directions.
△ Less
Submitted 23 February, 2022;
originally announced February 2022.
-
Topological Data Analysis of Black Hole Images
Authors:
Pierre Christian,
Chi-kwan Chan,
Anthony Hsu,
Feryal Ozel,
Dimitrios Psaltis,
Iniyan Natarajan
Abstract:
Features such as photon rings, jets, or hot. spots can leave particular topological signatures in a black hole image. As such, topological data analysis can be used to characterize images resulting from high resolution observations (synthetic or real) of black holes in the electromagnetic sector. We demonstrate that persistent homology allows for this characterization to be made automatically by c…
▽ More
Features such as photon rings, jets, or hot. spots can leave particular topological signatures in a black hole image. As such, topological data analysis can be used to characterize images resulting from high resolution observations (synthetic or real) of black holes in the electromagnetic sector. We demonstrate that persistent homology allows for this characterization to be made automatically by counting the number of connected components and one-dimensional holes. Further, persistent homology also allows for the distance between connected components or diameter of holes to be extracted from the image. In order to apply persistent homology on synthetic black hole images, we also introduce metronization, a new algorithm to prepare black hole images into a form that is suitable for topological analysis.
△ Less
Submitted 8 October, 2022; v1 submitted 1 February, 2022;
originally announced February 2022.
-
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.
-
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.
-
Constraints on black-hole charges with the 2017 EHT observations of M87*
Authors:
Prashant Kocherlakota,
Luciano Rezzolla,
Heino Falcke,
Christian M. Fromm,
Michael Kramer,
Yosuke Mizuno,
Antonios Nathanail,
Hector Olivares,
Ziri Younsi,
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,
Lindy Blackburn,
Raymond Blundell,
Wilfred Boland
, et al. (212 additional authors not shown)
Abstract:
Our understanding of strong gravity near supermassive compact objects has recently improved thanks to the measurements made by the Event Horizon Telescope (EHT). We use here the M87* shadow size to infer constraints on the physical charges of a large variety of nonrotating or rotating black holes. For example, we show that the quality of the measurements is already sufficient to rule out that M87*…
▽ More
Our understanding of strong gravity near supermassive compact objects has recently improved thanks to the measurements made by the Event Horizon Telescope (EHT). We use here the M87* shadow size to infer constraints on the physical charges of a large variety of nonrotating or rotating black holes. For example, we show that the quality of the measurements is already sufficient to rule out that M87* is a highly charged dilaton black hole. Similarly, when considering black holes with two physical and independent charges, we are able to exclude considerable regions of the space of parameters for the doubly-charged dilaton and the Sen black holes.
△ Less
Submitted 19 May, 2021;
originally announced May 2021.
-
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.
-
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.
-
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.
-
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.
-
First M87 Event Horizon Telescope Results and the Role of ALMA
Authors:
Ciriaco Goddi,
Geoff Crew,
Violette Impellizzeri,
Ivan Marti-Vidal,
Lynn D. Matthews,
Hugo Messias,
Helge Rottmann,
Walter Alef,
Lindy Blackburn,
Thomas Bronzwaer,
Chi-Kwan Chan,
Jordy Davelaar,
Roger Deane,
Jason Dexter,
Shep Doeleman,
Heino Falcke,
Vincent L. Fish,
Raquel Fraga-Encinas,
Christian M. Fromm,
Ruben Herrero-Illana,
Sara Issaoun,
David James,
Michael Janssen,
Michael Kramer,
Thomas P. Krichbaum
, et al. (19 additional authors not shown)
Abstract:
In April 2019, the Event Horizon Telescope (EHT) collaboration revealed the first image of the candidate super-massive black hole (SMBH) at the centre of the giant elliptical galaxy Messier 87 (M87). This event-horizon-scale image shows a ring of glowing plasma with a dark patch at the centre, which is interpreted as the shadow of the black hole. This breakthrough result, which represents a powerf…
▽ More
In April 2019, the Event Horizon Telescope (EHT) collaboration revealed the first image of the candidate super-massive black hole (SMBH) at the centre of the giant elliptical galaxy Messier 87 (M87). This event-horizon-scale image shows a ring of glowing plasma with a dark patch at the centre, which is interpreted as the shadow of the black hole. This breakthrough result, which represents a powerful confirmation of Einstein's theory of gravity, or general relativity, was made possible by assembling a global network of radio telescopes operating at millimetre wavelengths that for the first time included the Atacama Large Millimeter/ submillimeter Array (ALMA). The addition of ALMA as an anchor station has enabled a giant leap forward by increasing the sensitivity limits of the EHT by an order of magnitude, effectively turning it into an imaging array. The published image demonstrates that it is now possible to directly study the event horizon shadows of SMBHs via electromagnetic radiation, thereby transforming this elusive frontier from a mathematical concept into an astrophysical reality. The expansion of the array over the next few years will include new stations on different continents - and eventually satellites in space. This will provide progressively sharper and higher-fidelity images of SMBH candidates, and potentially even movies of the hot plasma orbiting around SMBHs. These improvements will shed light on the processes of black hole accretion and jet formation on event-horizon scales, thereby enabling more precise tests of general relativity in the truly strong field regime.
△ Less
Submitted 22 October, 2019;
originally announced October 2019.
-
The Synergy between VLBI and Gaia astrometry
Authors:
Huib van Langevelde,
Luis Henry Quiroga-Nuñez,
Wouter Vlemmings,
Laurent Loinard,
Mareki Honma,
Akiharu Nakagawa,
Katharina Immer,
Ross Burns,
Ylva Pihlström,
Lorant Sjouwerman,
R. Michael Rich,
Iniyan Natarajan,
Roger Deane
Abstract:
With the publication of Gaia DR2, 1.3 billion stars now have public parallax and proper motion measurements. In this contribution, we compare the results for sources that have both optical and radio measurements, focusing on circumstellar masers. For these large, variable and bright AGB stars, the VLBI astrometry results can be more robust. Moreover, there are a number of applications where VLBI a…
▽ More
With the publication of Gaia DR2, 1.3 billion stars now have public parallax and proper motion measurements. In this contribution, we compare the results for sources that have both optical and radio measurements, focusing on circumstellar masers. For these large, variable and bright AGB stars, the VLBI astrometry results can be more robust. Moreover, there are a number of applications where VLBI astrometry provides unique data for studying stellar populations and Galactic structure. The BeSSel project not only provides parallax and proper motions at much larger distances than Gaia can reach, but it also uniquely samples the spiral arms of the Galaxy. The evolved stars in the BAaDE sample can potentially constrain the dynamics and stellar content of the inner bulge and bar of the Milky Way, not reachable in the optical.
△ Less
Submitted 31 January, 2019; v1 submitted 23 January, 2019;
originally announced January 2019.
-
Resolving the blazar CGRaBS J0809+5341 in the presence of telescope systematics
Authors:
Iniyan Natarajan,
Zsolt Paragi,
Jonathan Zwart,
Simon Perkins,
Oleg Smirnov,
Kurt van der Heyden
Abstract:
We analyse Very Long Baseline Interferometry (VLBI) observations of the blazar CGRaBS J0809+5341 using Bayesian inference methods. The observation was carried out at 5 GHz using 8 telescopes that form part of the European VLBI Network. Imaging and deconvolution using traditional methods imply that the blazar is unresolved. To search for source structure beyond the diffraction limit, we perform Bay…
▽ More
We analyse Very Long Baseline Interferometry (VLBI) observations of the blazar CGRaBS J0809+5341 using Bayesian inference methods. The observation was carried out at 5 GHz using 8 telescopes that form part of the European VLBI Network. Imaging and deconvolution using traditional methods imply that the blazar is unresolved. To search for source structure beyond the diffraction limit, we perform Bayesian model selection between three source models (point, elliptical Gaussian, and circular Gaussian). Our modelling jointly accounts for antenna-dependent gains and system equivalent flux densities. We obtain posterior distributions for the various source and instrumental parameters along with the corresponding uncertainties and correlations between them. We find that there is very strong evidence (>1e9 :1) in favour of elliptical Gaussian structure and using this model derive the apparent brightness temperature distribution of the blazar, accounting for uncertainties in the shape estimates. To test the integrity of our method, we also perform model selection on synthetic observations and use this to develop a Bayesian criterion for the minimum resolvable source size and consequently the maximum measurable brightness temperature for a given interferometer, dependent on the signal-to-noise ratio (SNR) of the data incorporating the aforementioned systematics. We find that calibration errors play an increasingly important role in determining the over-resolution limit for SNR>>100. We show that it is possible to exploit the resolving power of future VLBI arrays down to about 5 per cent of the size of the naturally-weighted restoring beam, if the gain calibration is precise to 1 per cent or less.
△ Less
Submitted 12 October, 2016;
originally announced October 2016.
-
Bayesian Inference for Radio Observations - Going beyond deconvolution
Authors:
Michelle Lochner,
Bruce A. Bassett,
Martin Kunz,
Iniyan Natarajan,
Nadeem Oozeer,
Oleg Smirnov,
Jon Zwart
Abstract:
Radio interferometers suffer from the problem of missing information in their data, due to the gaps between the antennas. This results in artifacts, such as bright rings around sources, in the images obtained. Multiple deconvolution algorithms have been proposed to solve this problem and produce cleaner radio images. However, these algorithms are unable to correctly estimate uncertainties in deriv…
▽ More
Radio interferometers suffer from the problem of missing information in their data, due to the gaps between the antennas. This results in artifacts, such as bright rings around sources, in the images obtained. Multiple deconvolution algorithms have been proposed to solve this problem and produce cleaner radio images. However, these algorithms are unable to correctly estimate uncertainties in derived scientific parameters or to always include the effects of instrumental errors. We propose an alternative technique called Bayesian Inference for Radio Observations (BIRO) which uses a Bayesian statistical framework to determine the scientific parameters and instrumental errors simultaneously directly from the raw data, without making an image. We use a simple simulation of Westerbork Synthesis Radio Telescope data including pointing errors and beam parameters as instrumental effects, to demonstrate the use of BIRO.
△ Less
Submitted 14 September, 2015;
originally announced September 2015.
-
Montblanc: GPU accelerated Radio Interferometer Measurement Equations in support of Bayesian Inference for Radio Observations
Authors:
Simon Perkins,
Patrick Marais,
Jonathan Zwart,
Iniyan Natarajan,
Cyril Tasse,
Oleg Smirnov
Abstract:
We present Montblanc, a GPU implementation of the Radio interferometer measurement equation (RIME) in support of the Bayesian inference for radio observations (BIRO) technique. BIRO uses Bayesian inference to select sky models that best match the visibilities observed by a radio interferometer. To accomplish this, BIRO evaluates the RIME multiple times, varying sky model parameters to produce mult…
▽ More
We present Montblanc, a GPU implementation of the Radio interferometer measurement equation (RIME) in support of the Bayesian inference for radio observations (BIRO) technique. BIRO uses Bayesian inference to select sky models that best match the visibilities observed by a radio interferometer. To accomplish this, BIRO evaluates the RIME multiple times, varying sky model parameters to produce multiple model visibilities. Chi-squared values computed from the model and observed visibilities are used as likelihood values to drive the Bayesian sampling process and select the best sky model.
As most of the elements of the RIME and chi-squared calculation are independent of one another, they are highly amenable to parallel computation. Additionally, Montblanc caters for iterative RIME evaluation to produce multiple chi-squared values. Modified model parameters are transferred to the GPU between each iteration.
We implemented Montblanc as a Python package based upon NVIDIA's CUDA architecture. As such, it is easy to extend and implement different pipelines. At present, Montblanc supports point and Gaussian morphologies, but is designed for easy addition of new source profiles.
Montblanc's RIME implementation is performant: On an NVIDIA K40, it is approximately 250 times faster than MeqTrees on a dual hexacore Intel E5-2620v2 CPU. Compared to the OSKAR simulator's GPU-implemented RIME components it is 7.7 and 12 times faster on the same K40 for single and double-precision floating point respectively. However, OSKAR's RIME implementation is more general than Montblanc's BIRO-tailored RIME.
Theoretical analysis of Montblanc's dominant CUDA kernel suggests that it is memory bound. In practice, profiling shows that is balanced between compute and memory, as much of the data required by the problem is retained in L1 and L2 cache.
△ Less
Submitted 19 June, 2015; v1 submitted 30 January, 2015;
originally announced January 2015.
-
Bayesian Inference for Radio Observations
Authors:
Michelle Lochner,
Iniyan Natarajan,
Jonathan T. L. Zwart,
Oleg Smirnov,
Bruce A. Bassett,
Nadeem Oozeer,
Martin Kunz
Abstract:
New telescopes like the Square Kilometre Array (SKA) will push into a new sensitivity regime and expose systematics, such as direction-dependent effects, that could previously be ignored. Current methods for handling such systematics rely on alternating best estimates of instrumental calibration and models of the underlying sky, which can lead to inadequate uncertainty estimates and biased results…
▽ More
New telescopes like the Square Kilometre Array (SKA) will push into a new sensitivity regime and expose systematics, such as direction-dependent effects, that could previously be ignored. Current methods for handling such systematics rely on alternating best estimates of instrumental calibration and models of the underlying sky, which can lead to inadequate uncertainty estimates and biased results because any correlations between parameters are ignored. These deconvolution algorithms produce a single image that is assumed to be a true representation of the sky, when in fact it is just one realization of an infinite ensemble of images compatible with the noise in the data. In contrast, here we report a Bayesian formalism that simultaneously infers both systematics and science. Our technique, Bayesian Inference for Radio Observations (BIRO), determines all parameters directly from the raw data, bypassing image-making entirely, by sampling from the joint posterior probability distribution. This enables it to derive both correlations and accurate uncertainties, making use of the flexible software MEQTREES to model the sky and telescope simultaneously. We demonstrate BIRO with two simulated sets of Westerbork Synthesis Radio Telescope data sets. In the first, we perform joint estimates of 103 scientific (flux densities of sources) and instrumental (pointing errors, beamwidth and noise) parameters. In the second example, we perform source separation with BIRO. Using the Bayesian evidence, we can accurately select between a single point source, two point sources and an extended Gaussian source, allowing for 'super-resolution' on scales much smaller than the synthesized beam.
△ Less
Submitted 21 May, 2015; v1 submitted 21 January, 2015;
originally announced January 2015.
