Search

Dark Matter (DM) remains poorly probed on critical, sub-galactic scales, where predictions from different models diverge in terms of abundance and density profiles of halos. Gravitational lens systems on milli-arcsecond scales (milli-lenses) are expected for a population of dense DM halos, or free-floating supermassive black holes (SMBHs), that might be comprised of primordial black holes (PBHs), in the mass range of $10^6$ to $10^9 M_\odot$. In this paper, we aim to look for milli-lens systems via a systematic search in a large sample of radio-loud AGN observed with very-long-baseline interferometry (VLBI). We present the observational strategy to discriminate milli-lenses from contaminant objects mimicking a milli-lens morphology. In a pilot project, we have investigated VLBI images from 13,828 sources from the Astrogeo VLBI image database and reduced the number of candidates to 40 in a first step. We present here the images and analysis of sensitive follow-up observations with the EVN at 5 and 22 GHz, and streamline our analysis to reject milli-lens candidates. Using constraints such as the surface brightness ratio, conservation of spectral shape, stability of flux ratios over time, and changes in morphology, we can confidently discriminate between milli-lenses and their mimickers. Using the above constraints, we rule out 32 out of our initial 40 candidates as milli-lenses, demonstrating the power of our approach. Also, we find many new candidates for compact symmetric objects, that are thought to be short-lived, jetted radio sources. This serves as a pathfinder for the final sample used for the SMILE (Search for MIlli-LEnses) project, which will allow us to constrain DM models by comparing the results to theoretical predictions. This SMILE sample will consist of $\sim$5,000 sources based on the VLA CLASS survey, including many observations obtained for this project specifically., Comment: 35 pages, 36 figures, submitted for publication in Astronomy & Astrophysics

The spatial scales of relativistic radio jets, probed by relativistic magneto-hydrodynamic jet launching simulations (RMHDs) and by most very-long-baseline interferometry (VLBI) observations differ by an order of magnitude. Bridging the gap between these RMHD simulations and VLBI observations requires selecting nearby active galactic nuclei (AGN), the parsec-scale region of which can be resolved. 3C 84 is a nearby bright AGN fulfilling the necessary requirements: it is launching a powerful, relativistic jet powered by a central supermassive black hole, while also being very bright. Using 22 GHz global VLBI measurements of 3C 84 we aim to study its sub-parsec region in both total intensity and linear polarisation, to explore the properties of this jet, with a linear resolution of $\sim0.1$ parsec. We test different simulation setups by altering the bulk Lorentz factor $\Gamma$ of the jet, as well as the magnetic field configuration (toroidal, poloidal, helical). We confirm the persistence of a limb brightened structure, which reaches deep into the sub-parsec region. The corresponding electric vector position angles (EVPAs) follow the bulk jet flow inside but tend to be orthogonal to it near the edges. Our state-of-the-art RMHD simulations show that this geometry is consistent with a spine-sheath model, associated with a mildly relativistic flow and a toroidal magnetic field configuration., Comment: 10 pages, 8 figures, accepted for publication in Astronomy and Astrophysics

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 $\nu_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., Comment: 15 pages, 6 figures, published in A&A

The Event Horizon Telescope (EHT) observations carried out in 2018 April at 1.3 mm wavelengths included 9 stations in the array, comprising 7 single-dish telescopes and 2 phased arrays. The metadata package for the 2018 EHT observing campaign contains calibration tables required for the a-priori amplitude calibration of the 2018 April visibility data. This memo is the official documentation accompanying the release of the 2018 EHT metadata package, providing an overview of the contents of the package. We describe how telescope sensitivities, gain curves and other relevant parameters for each station in the EHT array were collected, processed, and validated to produce the calibration tables., Comment: 26 pages, 7 figures, EHT Memo Series 2023-L1-01

We report four novel position angle measurements of the core region of M81* at 5GHz and 8GHz, which confirm the presence of sinusoidal jet precession of the M81 jet region as suggested by \cite{Marti-Vidal2011}. The model makes three testable predictions on the evolution of the jet precession, which we test in our data with observations in 2017, 2018, and 2019. Our data confirms a precession period of $\sim7~\mathrm{yr}$ on top of a small linear drift. We further show that two 8 GHz observation are consistent with a precession period of $\sim 7~\mathrm{yr}$, but show a different time-lag w.r.t. to the 5 GHz and 1.7 GHz observations. We do not find a periodic modulation of the light curve with the jet precession, and therefore rule out a Doppler nature of the historic 1998-2002 flare. Our observations are consistent with either a binary black hole origin of the precession or the Lense-Thirring effect.

The discovery of high-energy astrophysical neutrinos by IceCube has opened a new window to the Universe. However, the origin of these neutrinos is still a mystery, and some of them could be a result of dark matter interactions such as decay. Next generation gigaton water-Cherenkov neutrino telescope, KM3NeT, is expected to offer significantly improved energy resolution in the cascade channel, and advantageous viewing condition to the Galactic Center; both important for searches of dark matter decay signals. We study the sensitivity of KM3NeT on dark matter decays by performing a mock likelihood analysis for both cascade and track type events, taking into account both angular and energy information. We find that, combining both channels, KM3NeT is expected to produce world leading limits on dark matter decay lifetime in the PeV mass range, and could test some of the dark matter hints in the current IceCube data., Comment: 9 pages, 4 figures. Comments welcome!

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., Comment: 20 pages, 15 figures, accepted for publication in A&A

We present multi-wavelength measurements of the thermal, chemical, and cloud contrasts associated with the visibly dark formations (also known as 5-$\mu$m hot spots) and intervening bright plumes on the boundary between Jupiter's Equatorial Zone (EZ) and North Equatorial Belt (NEB). Observations made by the TEXES 5-20 $\mu$m spectrometer at the Gemini North Telescope in March 2017 reveal the upper-tropospheric properties of 12 hot spots, which are directly compared to measurements by Juno using the Microwave Radiometer (MWR), JIRAM at 5 $\mu$m, and JunoCam visible images. MWR and thermal-infrared spectroscopic results are consistent near 0.7 bar. Mid-infrared-derived aerosol opacity is consistent with that inferred from visible-albedo and 5-$\mu$m opacity maps. Aerosol contrasts, the defining characteristics of the cloudy plumes and aerosol-depleted hot spots, are not a good proxy for microwave brightness. The hot spots are neither uniformly warmer nor ammonia-depleted compared to their surroundings at $p<1$ bar. At 0.7 bar, the microwave brightness at the edges of hot spots is comparable to other features within the NEB. Conversely, hot spots are brighter at 1.5 bar, signifying either warm temperatures and/or depleted NH$_3$ at depth. Temperatures and ammonia are spatially variable within the hot spots, so the precise location of the observations matters to their interpretation. Reflective plumes sometimes have enhanced NH$_3$, cold temperatures, and elevated aerosol opacity, but each plume appears different. Neither plumes nor hot spots had microwave signatures in channels sensing $p>10$ bars, suggesting that the hot-spot/plume wave is a relatively shallow feature., Comment: 50 pages, 16 figures, accepted for publication in JGR-Planets

The galaxy M 87 harbors a kiloparsec-scale relativistic jet, whose origin coincides with a supermassive black hole. Observational mm-VLBI campaigns are capable of resolving the jet-launching region at the scale of the event horizon. In order to provide a context for interpreting these observations, realistic general-relativistic magnetohydrodynamical (GRMHD) models of the accretion flow are constructed. The characteristics of the observed spectral-energy distribution (SED) depend on the shape of the electrons' energy-distribution function (eDF). The dependency on the eDF is omitted in the modeling of the first Event Horizon Telescope results. In this work, we aim to model the M 87 SED from radio up to NIR/optical frequencies using a thermal-relativistic Maxwell- J\"uttner distribution, as well as a relativistic $\kappa$-distribution function. The electrons are injected based on sub-grid, particle-in-cell parametrizations for sub-relativistic reconnection. A GRMHD simulation in Cartesian-Kerr-Schild coordinates, using eight levels of adaptive mesh refinement (AMR), forms the basis of our model. To obtain spectra and images, the GRMHD data is post-processed with the ray-tracing code RAPTOR, which is capable of ray tracing through AMR GRMHD simulation data. We obtain radio spectra in both the thermal-jet and $\kappa$-jet models consistent with radio observations. Additionally, the $\kappa$-jet models also recover the NIR/optical emission. The models recover the observed source sizes and core shifts and obtain a jet power of $\approx 10^{43}$ ergs/s. In the $\kappa$-jet models, both the accretion rates and jet powers are approximately two times lower than the thermal-jet model. The frequency cut-off observed at $\nu \approx 10^{15}$ Hz is recovered when the accelerator size is $10^6$ - $10^8$ cm, this could potentially point to an upper limit for plasmoid sizes in the jet of M 87., Comment: 15 pages, 10 figures

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., Comment: 6 pages, 1 figure, EVN 2018 symposium proceedings

The Galactic Center supermassive black hole Sagittarius A* (Sgr A*) is one of the most promising targets to study the dynamics of black hole accretion and outflow via direct imaging with very long baseline interferometry (VLBI). At 3.5 mm (86 GHz), the emission from Sgr A* is resolvable with the Global Millimeter VLBI Array (GMVA). We present the first observations of Sgr A* with the phased Atacama Large Millimeter/submillimeter Array (ALMA) joining the GMVA. Our observations achieve an angular resolution of ~87{\mu}as, improving upon previous experiments by a factor of two. We reconstruct a first image of the unscattered source structure of Sgr A* at 3.5 mm, mitigating effects of interstellar scattering. The unscattered source has a major axis size of 120 $\pm$ 34{\mu}as (12 $\pm$ 3.4 Schwarzschild radii), and a symmetrical morphology (axial ratio of 1.2$^{+0.3}_{-0.2}$), which is further supported by closure phases consistent with zero within 3{\sigma}. We show that multiple disk-dominated models of Sgr A* match our observational constraints, while the two jet-dominated models considered are constrained to small viewing angles. Our long-baseline detections to ALMA also provide new constraints on the scattering of Sgr A*, and we show that refractive scattering effects are likely to be weak for images of Sgr A* at 1.3 mm with the Event Horizon Telescope. Our results provide the most stringent constraints to date for the intrinsic morphology and refractive scattering of Sgr A*, demonstrating the exceptional contribution of ALMA to millimeter VLBI., Comment: Accepted for publication in ApJ

The enhancement and detection of elongated structures in noisy image data is relevant for many biomedical imaging applications. To handle complex crossing structures in 2D images, 2D orientation scores $U: \mathbb{R} ^ 2\times S ^ 1 \rightarrow \mathbb{C}$ were introduced, which already showed their use in a variety of applications. Here we extend this work to 3D orientation scores $U: \mathbb{R} ^ 3 \times S ^ 2\rightarrow \mathbb{C}$. First, we construct the orientation score from a given dataset, which is achieved by an invertible coherent state type of transform. For this transformation we introduce 3D versions of the 2D cake-wavelets, which are complex wavelets that can simultaneously detect oriented structures and oriented edges. Here we introduce two types of cake-wavelets, the first uses a discrete Fourier transform, the second is designed in the 3D generalized Zernike basis, allowing us to calculate analytical expressions for the spatial filters. Finally, we show two applications of the orientation score transformation. In the first application we propose an extension of crossing-preserving coherence enhancing diffusion via our invertible orientation scores of 3D images which we apply to real medical image data. In the second one we develop a new tubularity measure using 3D orientation scores and apply the tubularity measure to both artificial and real medical data.

Einstein's General Theory of Relativity (GR) successfully describes gravity. The most fundamental predictions of GR are black holes (BHs), but in spite of many convincing BH candidates in the Universe, there is no conclusive experimental proof of their existence using astronomical observations in the electromagnetic spectrum. Are BHs real astrophysical objects? Does GR hold in its most extreme limit or are alternatives needed? The prime target to address these fundamental questions is in the center of our own Galaxy, which hosts the closest and best-constrained supermassive BH candidate in the Universe, Sagittarius A* (Sgr A*). Three different types of experiments hold the promise to test GR in a strong-field regime using observations of Sgr A* with new-generation instruments. The first experiment aims to image the relativistic plasma emission which surrounds the event horizon and forms a "shadow" cast against the background, whose predicted size (~50 microarcseconds) can now be resolved by upcoming VLBI experiments at mm-waves such as the Event Horizon Telescope (EHT). The second experiment aims to monitor stars orbiting Sgr A* with the upcoming near-infrared interferometer GRAVITY at the Very Large Telescope (VLT). The third experiment aims to time a radio pulsar in tight orbit about Sgr A* using radio telescopes (including the Atacama Large Millimeter Array or ALMA). The BlackHoleCam project exploits the synergy between these three different techniques and aims to measure the main BH parameters with sufficient precision to provide fundamental tests of GR and probe the spacetime around a BH in any metric theory of gravity. Here, we review our current knowledge of the physical properties of Sgr A* as well as the current status of such experimental efforts towards imaging the event horizon, measuring stellar orbits, and timing pulsars around Sgr A*., Comment: review paper, 36 pages, 12 figures (v2 is the published version)