Very Long Baseline Interferometry Array Research Articles

Sgr A*, M87* and OJ287 as Extreme Central Binaries of Supermassive Black Holes- -Confirmation and extension of decameter radio wave observation results through VLBI data comparison

Based on a deduced 2153±5 s orbiting period of the extreme central binary of the supermassive black hole (ECB-SMBH), at Sgr A*,that is confirmed being based on the decameter radio wave pulse observations compared with 1.3 mm wavelength very long baseline interferometry (VLBI) data from the event horizon telescope (EHT-Data), we depict moving images of Sgr A* from EHT-Data, confirming further the existence of ECB-SMBH. The results show that radio wave images of the two members of ECB-SMBH are radio-bright objects due to emissions caused by source agents accreting from a thick disk toward ECB-SMBH. We extend the confirmed concept of ECB-SMBH to the interpretation of M87*, based on two images released by Miyoshi et al. in 2022 using EHT-Data and by Lu et al. in 2023 using the 3.5 mm wavelength data from the Global Millimeter VLBI Arrays as core telescopes. For the double brilliant spots displayed in two different images of M87*, we consider them the same two objects that shift locations during the 369-day interval of the observation periods of the two images. It is concluded that M87* A and M87* B (temporarily named) are orbiting with a period of 168.8 d at velocities of (6.28±0.54)% and (16.7±1.5)% of the speed of light, respectively. For the existence of ECB-SMBH, the condition of “no gravitational wave radiation” from the SMBH is essential. Regarding the secondary BH of blazar OJ 287, we find around 58~60-year period oscillation of the orbiting period, which is recently in a shortening phase. By applying the concept of ECB-SMBH to the OJ287 primary BH, whose orbiting period is approximately 1.6 yr, we interpret that the orbital period variation of the secondary BH is not due to the radiation of gravitational waves but is caused by the orbital motions of ECB-SMBH

SIMULATIONS OF LINEAR POLARIZATION OF PRECESSING AGN JETS AT PARSEC SCALES

The latest results of the most detailed analysis of multi-epoch polarization-sensitive observations of active galactic nuclei (AGN) jets at parsecs scales by very long baseline interferometry (VLBI) reveal several characteristic patterns of linear polarization distribution and its variability [1, 2]. Some of the observed profiles can be reproduced by a simple model of a jet threaded by a helical magnetic field. However, none of the models presented to date can explain the observed polarization profiles with an increase in its degree towards the edges of the jet, and accompanied by a “fountain” type electrical vector pattern and its high temporal variability in the center. Based on simulations of the VLBI observations of relativistic jets, we show here that the observed transverse linear polarization profiles, atypical for the simple magnetic field models can be naturally produced assuming the finite resolution of VLBI arrays and precession of a jet on ten-years scales, observational indications of which are found in an increasing number of AGN. In our simulations, we qualitatively reproduce the distribution of the electric vector and its variability, though the polarization images are characterized by a bright spine due to weak smearing, which is poorly consistent with observations. More effective depolarization can be obtained in models with the suppressed emission of the jet spine.

Dynamic and Polarimetric VLBI imaging with a multiscalar approach

Context. Due to the limited number of antennas and the limited observation time, an array of antennas in very long baseline interfer-ometry (VLBI) often samples the Fourier domain only very sparsely. Powerful deconvolution algorithms are needed to compute a final image. Multiscale imaging approaches such as DoG-HiT have recently been developed to solve the VLBI imaging problem and show promising performance: they are fast, accurate, unbiased, and automatic. Aims. We extend the multiscalar imaging approach to polarimetric imaging, to reconstructions of dynamically evolving sources, and finally to dynamic polarimetric reconstructions. Methods. These extensions (mr-support imaging) utilize a multiscalar approach. The time-averaged Stokes I image was decomposed by a wavelet transform into single subbands. We used the set of statistically significant wavelet coefficients, the multiresolution support (mr-support), computed by DoG-HiT as a prior in a constrained minimization manner; we fitted the single-frame (polarimetric) observables by only varying the coefficients in the multiresolution support. Results. The Event Horizon Telescope (EHT) is a VLBI array imaging supermassive black holes. We demonstrate on synthetic data that mr-support imaging offers ample regularization and is able to recover simple geometric dynamics at the horizon scale in a typical EHT setup. The approach is relatively lightweight, fast, and largely automatic and data driven. The ngEHT is a planned extension of the EHT designed to recover movies at the event horizon scales of a supermassive black hole. We benchmark the performance of mr-support imaging for the denser ngEHT configuration demonstrating the major improvements the additional ngEHT antennas will bring to dynamic polarimetric reconstructions. Conclusions. Current and upcoming instruments offer the observational possibility to do polarimetric imaging of dynamically evolving structural patterns with the highest spatial and temporal resolution. State-of-the-art dynamic reconstruction methods can capture this motion with a range of temporal regularizers and priors. With this work, we add an additional simpler regularizer to the list: constraining the reconstruction to the multiresolution support.

Astrometric VLBI observations of H2O masers in an extreme OH/IR star candidate NSV 17351

Abstract The results of astrometric very long baseline interferometry (VLBI) observations toward an extreme OH/IR star candidate NSV 17351 are presented. Using the VERA (VLBI Exploration of Radio Astrometry) VLBI array , we observed 22 GHz H2O masers of NSV 17351 and derived an annual parallax of 0.247 ± 0.035 mas, which corresponds to a distance of 4.05 ± 0.59 kpc, from the observation. By averaging proper motions of 15 maser spots further, we determined the systemic proper motion of NSV 17351 to be (μαcos δ, μδ)avg. = (−1.19 ± 0.11, 1.30 ± 0.19) mas yr−1. The maser spots spread out over a region of 20 × 30 mas, which can be converted to a spatial distribution of ∼80 × 120 au at the source distance. Internal motions of the maser spots suggest an outward-moving maser region with respect to the estimated position of the central star. From single-dish monitoring of the H2O maser emission, we estimate the pulsation period of NSV 17351 to be 1122 ± 24 d. This is the first report of the periodic activity of NSV 17351, indicating that NSV 17351 could have a mass of ∼4 M⊙. We confirmed that the time variation of H2O masers can be used as a period estimator of variable OH/IR stars. Furthermore, by inspecting dozens of double-peaked H2O maser spectra for the last 40 years, we discovered the long-term acceleration in the radial velocity of the circumstellar matter to be 0.17 ± 0.03 km s−1 yr−1. We finally determined the position and kinematics in the Milky Way Galaxy and found that NSV 17351 is located in an interarm region between the Outer and Perseus arms. We note that the astrometric VLBI observation toward extreme OH/IR stars shows us a useful sample of the Galactic dynamics.

Millimeter/Submillimeter VLBI with a Next Generation Large Radio Telescope in the Atacama Desert

The proposed next generation Event Horizon Telescope (ngEHT) concept envisions the imaging of various astronomical sources on scales of microarcseconds in unprecedented detail with at least two orders of magnitude improvement in the image dynamic ranges by extending the Event Horizon Telescope (EHT). A key technical component of ngEHT is the utilization of large aperture telescopes to anchor the entire array, allowing the connection of less sensitive stations through highly sensitive fringe detections to form a dense network across the planet. Here, we introduce two projects for planned next generation large radio telescopes in the 2030s on the Chajnantor Plateau in the Atacama desert in northern Chile, the Large Submillimeter Telescope (LST) and the Atacama Large Aperture Submillimeter Telescope (AtLAST). Both are designed to have a 50-meter diameter and operate at the planned ngEHT frequency bands of 86, 230 and 345 GHz. A large aperture of 50 m that is co-located with two existing EHT stations, the Atacama Large Millimeter/Submillimeter Array (ALMA) and the Atacama Pathfinder Experiment (APEX) Telescope in the excellent observing site of the Chajnantor Plateau, will offer excellent capabilities for highly sensitive, multi-frequency, and time-agile millimeter very long baseline interferometry (VLBI) observations with accurate data calibration relevant to key science cases of ngEHT. In addition to ngEHT, its unique location in Chile will substantially improve angular resolutions of the planned Next Generation Very Large Array in North America or any future global millimeter VLBI arrays if combined. LST and AtLAST will be a key element enabling transformative science cases with next-generation millimeter/submillimeter VLBI arrays.

Collimation of the Relativistic Jet in the Quasar 3C 273

The collimation of relativistic jets launched from the vicinity of supermassive black holes (SMBHs) at the centers of active galactic nuclei (AGNs) is one of the key questions to understand the nature of AGN jets. However, little is known about the detailed jet structure for AGN like quasars since very high angular resolutions are required to resolve these objects. We present very long baseline interferometry (VLBI) observations of the archetypical quasar 3C 273 at 86 GHz, performed with the Global Millimeter VLBI Array, for the first time including the Atacama Large Millimeter/submillimeter Array. Our observations achieve a high angular resolution down to ∼60 μas, resolving the innermost part of the jet ever on scales of ∼105 Schwarzschild radii. Our observations, including close-in-time High Sensitivity Array observations of 3C 273 at 15, 22, and 43 GHz, suggest that the inner jet collimates parabolically, while the outer jet expands conically, similar to jets from other nearby low-luminosity AGNs. We discovered the jet collimation break around 107 Schwarzschild radii, providing the first compelling evidence for structural transition in a quasar jet. The location of the collimation break for 3C 273 is farther downstream from the sphere of gravitational influence (SGI) from the central SMBH. With the results for other AGN jets, our results show that the end of the collimation zone in AGN jets is governed not only by the SGI of the SMBH but also by the more diverse properties of the central nuclei.

An Updated View of the Milky Way from Maser Astrometry

Astrometric observations of maser sources in the Milky Way, using the Very Long Baseline Interferometry (VLBI) technique, have been exploited to determine the spiral structure of our Galaxy. Several major spiral arms have now been pinpointed in the first and second Galactic quadrants. Fundamental Galactic parameters such as the distance to the Galactic Centre and the rotation curve and speed have been determined. In this review, we discuss the latest results from the Bar and Spiral Structure Legacy survey, the VLBI Exploration of Radio Astrometry survey and other VLBI arrays and compare them with astrometric measurements of stars from the Gaia mission. In particular, we present the peculiarities of the individual spiral arms and a thorough discussion of the methods to determine different Galactic parameters as well as the obtained values.

Inverse Multiview. I. Multicalibrator Inverse Phase Referencing for Microarcsecond Very Long Baseline Interferometry Astrometry

Very Long Baseline Interferometry (VLBI) astrometry is a well established technique for achieving ±10 μas parallax accuracies at frequencies well above 10 GHz. At lower frequencies, uncompensated interferometer delays associated with the ionosphere play the dominant role in limiting the astrometric accuracy. Multiview is a novel VLBI calibration method, which uses observations of multiple quasars to accurately model and remove time-variable, directional-dependent changes to the interferometer delay. Here we extend the Multiview technique by phase-referencing data to the target source (“inverse Multiview”) and test its performance. Multiple observations with a four-antenna VLBI array operating at 8.3 GHz show single-epoch astrometric accuracies near 20 μas for target–reference quasar separations up to about 7°. This represents an improvement in astrometric accuracy by up to an order of magnitude compared to standard phase-referencing.

Nonthermal Radio Continuum Emission from Young Nearby Stars

Young and magnetically active low-mass stars often exhibit nonthermal coronal radio emission owing to the gyration of electrons in their magnetized chromospheres. This emission is easily detectable at centimeter wavelengths with the current sensitivity of large radio interferometers like the Very Large Array (VLA). With the aim of identifying nearby stars adequate for future accurate radio astrometric monitoring using very long baseline interferometry (VLBI), we have used the VLA in its B configuration to search for radio emission at ν ≃ 6 GHz (λ ≃ 5 cm) toward a sample of 170 nearby (<130 pc), mostly young (5–500 Myr) stars of spectral types between F4 and M2. At our mean 3σ detection limit of ≃50 μJy, we identify 31 young stars with coronal radio emission (an 18% system detection rate) and more than 600 background (most likely extragalactic) sources. Among the targeted stars, we find a significant decline of the detection rate with age from 56% ± 20% for stars with ages ≤10 Myr to 10% ± 3% for stars with ages 100–200 Myr. No star older than 200 Myr was detected. The detection rate also declines with T eff from 36% ± 10% for stars with T eff < 4000 K to 13% ± 3% for earlier spectral types with T eff > 5000 K. The binarity fraction among the radio-bright stars is at least twice as high as among the radio-quiet stars. The radio-bright nearby young stars identified here provide an interesting sample for future astrometric studies using VLBI arrays aimed at searching for hitherto-unknown tight binary components or even exoplanets.

Superradiance evolution of black hole shadows revisited

Black hole (BH) shadows can be used to probe new physics in the form of ultra-light particles via the phenomenon of superradiant instability. By directly affecting the BH mass and spin, superradiance can lead to a time evolution of the BH shadow, which nonetheless has been argued to be unobservable through Very Long Baseline Interferometry (VLBI) over realistic observation timescales. We revisit the superradiance-induced BH shadow evolution including the competing effects of gas accretion and gravitational wave (GW) emission and, as a first step towards modelling realistic new physics scenarios which predict the existence of multiple ultra-light species, we study the system in the presence of two ultra-light bosons, whose combined effect could help reducing the shadow evolution timescale. We find that accretion and GW emission play a negligible role in our results (justifying previous simplified analyses), and that contrary to our intuition the inclusion of an additional ultra-light boson does not shorten the BH shadow evolution timescale and hence improve detection prospects. However, we point out an important subtlety concerning the observationally meaningful definition of the superradiance-induced BH shadow evolution timescale, which reduces the latter by about an order of magnitude, opening up the possibility of observing the superradiance-induced BH shadow evolution with upcoming VLBI arrays, provided angular resolutions just below the $\mu{\rm as}$ level can be reached. As a concrete example, we show that the angular size of the shadow of SgrA$^*$ can change by up to $0.6\,\mu{\rm as}$ over a period as short as $16$ years, which further strengthens the scientific case for targeting the shadow of SgrA$^*$ with next-generation VLBI arrays.

Radio emission in a nearby, ultra-cool dwarf binary: A multifrequency study

Context.The substellar triple system VHS J125601.92−125723.9 (hereafter VHS 1256−1257) is composed of an equal-mass M7.5 brown dwarf binary and an L7 low-mass substellar object. In Guirado et al. (2018, A&amp;A, 610, A23) we published the detection of radio emission at 8.4 GHz coming from the central binary and making it an excellent target for further observations.Aims.We aim to identify the origin of the radio emission occurring in the central binary of VHS 1256−1257 while discussing the expected mechanisms involved in the radio emission of ultra-cool dwarfs.Methods.We observed this system with theKarl G. JanskyVery Large Array, the European very-long-baseline interferometry (VLBI) Network, the enhanced Multi-Element Remotely Linked Interferometer Network, the NOrthern Extended Millimeter Array, and the Atacama Large Millimetre Array at frequencies ranging from 5 GHz up to 345 GHz in several epochs during 2017, 2018, and 2019.Results.We found radio emission at 6 GHz and 33 GHz coincident with the expected position of the central binary of VHS 1256−1257. The StokesIdensity fluxes detected were 73 ± 4 μJy and 83 ± 13 μJy, respectively, with no detectable circular polarisation or pulses. No emission is detected at higher frequencies (230 GHz and 345 GHz), nor at 5 GHz with VLBI arrays. The emission appears to be stable over almost three years at 6 GHz. To explain the constraints obtained both from the detections and non-detections, we considered multiple scenarios including thermal and nonthermal emission, and different contributions from each component of the binary.Conclusions.Our results can be well explained by nonthermal gyrosynchrotron emission originating at radiation belts with a low plasma density (ne = 300−700 cm−3), a moderate magnetic field strength (B ≈ 140 G), and an energy distribution of electrons following a power-law (dN/dE ∝ E−δ) withδfixed at 1.36. These radiation belts would need to be present in both components and also be viewed equatorially.

MeqSilhouette v2: spectrally resolved polarimetric synthetic data generation for the event horizon telescope

ABSTRACTWe 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.

East Asian VLBI Network observations of active galactic nuclei jets: imaging with KaVA+Tianma+Nanshan

The East Asian Very Long Baseline Interferometry (VLBI) Network (EAVN) is a rapidly evolving international VLBI array that is currently promoted under joint efforts among China, Japan and Korea. EAVN aims at forming a joint VLBI Network by combining a large number of radio telescopes distributed over East Asian regions. After the combination of the Korean VLBI Network (KVN) and the VLBI Exploration of Radio Astrometry (VERA) into KaVA, further expansion with the joint array in East Asia is actively promoted. Here we report the first imaging results (at 22 and 43 GHz) of bright radio sources obtained with KaVA connected to Tianma 65-m and Nanshan 26-m Radio Telescopes in China. To test the EAVN imaging performance for different sources, we observed four active galactic nuclei (AGN) having different brightness and morphology. As a result, we confirmed that the Tianma 65-m Radio Telescope (TMRT) significantly enhances the overall array sensitivity, a factor of 4 improvement in baseline sensitivity and 2 in image dynamic range compared to the case of KaVA only. The addition of the Nanshan 26-m Radio Telescope (NSRT) further doubled the east-west angular resolution. With the resulting high-dynamic-range, high-resolution images with EAVN (KaVA+TMRT+NSRT), various fine-scale structures in our targets, such as the counter-jet in M87, a kink-like morphology of the 3C 273 jet and the weak emission in other sources are successfully detected. This demonstrates the powerful capability of EAVN to study AGN jets and to achieve other science goals in general. Ongoing expansion of EAVN will further enhance the angular resolution, detection sensitivity and frequency coverage of the network.

Persistent Non-Gaussian Structure in the Image of Sagittarius A* at 86 GHz

Abstract Observations of the Galactic Center supermassive black hole Sagittarius A* (Sgr A*) with very long-baseline interferometry (VLBI) are affected by interstellar scattering along our line of sight. At long radio observing wavelengths (≲1 cm), the scattering heavily dominates image morphology. At 3.5 mm (86 GHz), the intrinsic source structure is no longer sub-dominant to scattering, and thus the intrinsic emission from Sgr A* is resolvable with the Global Millimeter VLBI Array (GMVA). Long-baseline detections to the phased Atacama Large Millimeter/submillimeter Array (ALMA) in 2017 provided new constraints on the intrinsic and scattering properties of Sgr A*, but the stochastic nature of the scattering Requires multiple observing epochs to reliably estimate its statistical properties. We present new observations with the GMVA+ALMA, taken in 2018, which confirm non-Gaussian structure in the scattered image seen in 2017. In particular, the ALMA–GBT baseline shows more flux density than expected for an anistropic Gaussian model, providing a tight constraint on the source size and an upper limit on the dissipation scale of interstellar turbulence. We find an intrinsic source extent along the minor axis of ∼100 μas both via extrapolation of longer wavelength scattering constraints and direct modeling of the 3.5 mm observations. Simultaneously fitting for the scattering parameters, we find an at-most modestly asymmetrical (major-to-minor axis ratio of 1.5 ± 0.2) intrinsic source morphology for Sgr A*.

Polarization calibration techniques for the new-generation VLBI

The calibration and analysis of polarization observations in very long baseline interferometry (VLBI) requires the use of specific algorithms that suffer from several limitations, closely related to assumptions in the data properties that may not hold in observations taken with new-generation VLBI equipment. Currently, the instantaneous bandwidth achievable with VLBI backends can be as high as several gigahertz, covering several radio bands simultaneously. In addition, the sensitivity of VLBI observations with the most updated equipment may reach dynamic ranges of tens of thousands, both in total intensity and polarization. In this paper, we discuss the impact of the limitations of common VLBI polarimetry algorithms on narrow-field observations taken with modern VLBI arrays, from the VLBI Global Observing System to the Event Horizon Telescope, and present new software that overcomes these limitations. In particular, our software is able to perform a simultaneous fit of multiple calibrator sources, include nonlinear terms in the model of the instrumental polarization and use a self-calibration approach for the estimate of the polarization leakage in the antenna receivers.

Yebes 40 m radio telescope and the broad band NANOCOSMOS receivers at 7 mm and 3 mm for line surveys.

Yebes 40m radio telescope is the main and largest observing instrument at Yebes Observatory and it is devoted to Very Long Baseline Interferometry (VLBI) and single dish observations since 2010. It has been covering frequency bands between 2 GHz and 90 GHz in discontinuous and narrow windows in most of the cases, to match the current needs of the European VLBI Network (EVN) and the Global Millimeter VLBI Array (GMVA). Nanocosmos project, a European Union funded synergy grant, opened the possibility to increase the instantaneous frequency coverage to observe many molecular transitions with single tunnings in single dish mode. This reduces the observing time and maximises the output from the telescope. We present the technical specifications of the recently installed 31.5 - 50GHz (Q band) and 72 - 90.5 GHz (W band) receivers along with the main characteristics of the telescope at these frequency ranges. We have observed IRC+10216, CRL 2688 and CRL 618, which harbour a rich molecular chemistry, to demonstrate the capabilities of the new instrumentation for spectral observations in single dish mode. The results show the high sensitivity of the telescope in the Q band. The spectrum of IRC+10126 offers a signal to noise ratio never seen before for this source in this band. On the other hand, the spectrum normalised by the continuum flux towards CRL 618 in the W band demonstrates that the 40 m radio telescope produces comparable results to those from the IRAM 30 m radio telescope, although with a smaller sensitivity. The new receivers fulfil one of the main goals of Nanocosmos and open the possibility to study the spectrum of different astrophysical media with unprecedented sensitivity.

GPCAL: A Generalized Calibration Pipeline for Instrumental Polarization in VLBI Data

Abstract We present the Generalized Polarization CALibration pipeline (GPCAL), an automated pipeline for instrumental polarization calibration of very long baseline interferometry (VLBI) data. The pipeline is designed to achieve a high calibration accuracy by means of fitting the instrumental polarization model, including the second-order terms, to data from multiple calibrator sources simultaneously. It also allows for using more accurate linear polarization models of calibrators for D-term estimation compared to the conventional way that assumes similar linear polarization and total intensity structures. This assumption has been widely used in the existing packages for instrumental polarization calibration but could be a source of significant uncertainties when there is no suitable calibrator satisfying the assumption. We demonstrate the capabilities of GPCAL by using simulated data, archival Very Long Baseline Array (VLBA) data of many active galactic nucleus (AGN) jets at 15 and 43 GHz, and our Korean VLBI Network (KVN) observations of many AGN jets at 86, 95, 130, and 142 GHz. The pipeline could reproduce the complex linear polarization structures of several sources shown in the previous studies using the same VLBA data. In addition, GPCAL reveals a complex linear polarization structure in the flat-spectrum radio quasar 3C 273 from the KVN data at all four frequencies. These results demonstrate that GPCAL can achieve a high calibration accuracy for various VLBI arrays.

Multiband RadioAstron space VLBI imaging of the jet in quasar S5 0836+710

Context. Detailed studies of relativistic jets in active galactic nuclei (AGN) require high-fidelity imaging at the highest possible resolution. This can be achieved using very long baseline interferometry (VLBI) at radio frequencies, combining worldwide (global) VLBI arrays of radio telescopes with a space-borne antenna on board a satellite. Aims. We present multiwavelength images made of the radio emission in the powerful quasar S5 0836+710, obtained using a global VLBI array and the antenna Spektr-R of the RadioAstron mission of the Russian Space Agency, with the goal of studying the internal structure and physics of the relativistic jet in this object. Methods. The RadioAstron observations at wavelengths of 18 cm, 6 cm, and 1.3 cm are part of the Key Science Program for imaging radio emission in strong AGN. The internal structure of the jet is studied by analyzing transverse intensity profiles and modeling the structural patterns developing in the flow. Results. The RadioAstron images reveal a wealth of structural detail in the jet of S5 0836+710 on angular scales ranging from 0.02 mas to 200 mas. Brightness temperatures in excess of 1013 K are measured in the jet, requiring Doppler factors of ≥100 for reconciling them with the inverse Compton limit. Several oscillatory patterns are identified in the ridge line of the jet and can be explained in terms of the Kelvin–Helmholtz (KH) instability. The oscillatory patterns are interpreted as the surface and body wavelengths of the helical mode of the KH instability. The interpretation provides estimates of the jet Mach number and of the ratio of the jet to the ambient density, which are found to be Mj ≈ 12 and η ≈ 0.33. The ratio of the jet to the ambient density should be conservatively considered an upper limit because its estimate relies on approximations.

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

Context. 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 the most 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 direct comparison with observational data. Aims. 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. Methods. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a millimetre 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, compared to the addition of thermal noise only. Using synthetic data based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M 87, we performed case studies to assess the image quality that can be obtained with the current and future EHT array for different weather conditions. Results. Our synthetic observations 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 our GRMHD source models can be recovered robustly with the EHT2017 array after performing calibration steps, which include fringe fitting, a priori amplitude and network calibration, and self-calibration. With the planned addition of new stations to the EHT array in the coming years, 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.

A heatwave of accretion energy traced by masers in the G358-MM1 high-mass protostar

High-mass stars are thought to accumulate much of their mass via short, infrequent bursts of disk-aided accretion. Such accretion events are rare and difficult to observe directly but are known to drive enhanced maser emission. In this Letter we report high-resolution, multi-epoch methanol maser observations toward G358.93-0.03, which reveal an interesting phenomenon: the subluminal propagation of a thermal radiation ‘heatwave’ emanating from an accreting high-mass protostar. The extreme transformation of the maser emission implies a sudden intensification of thermal infrared radiation from within the inner (40-mas, 270-au) region. Subsequently, methanol masers trace the radial passage of thermal radiation through the environment at ≥4% of the speed of light. Such a high translocation rate contrasts with the ≤10 km s-1 physical gas motions of methanol masers typically observed using very-long-baseline interferometry (VLBI). The observed scenario can readily be attributed to an accretion event in the high-mass protostar G358.93-0.03-MM1. While being the third case in its class, G358.93-0.03-MM1 exhibits unique attributes hinting at a possible ‘zoo’ of accretion burst types. These results promote the advantages of maser observations in understanding high-mass-star formation, both through single-dish maser monitoring campaigns and via their international cooperation as VLBI arrays.