Source author record

Michael D. Johnson

Michael D. Johnson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

35works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

35 published item(s)

preprint2026arXiv

The Role of Interferometric Phase in Measuring Black Hole Photon Rings

The Event Horizon Telescope (EHT) captured the first images of a black hole using Very Long Baseline Interferometry (VLBI). In the near future, extensions of the EHT such as the Black Hole Explorer (BHEX) will allow access to finer-scale features, such as a black hole's ''photon ring.'' In the Kerr spacetime, this image structure arises from strong gravitational lensing near the black hole that results in a series of increasingly demagnified images of each emitting region that exponentially converge to a limiting critical curve. Exotic black hole alternatives, such as wormholes, can introduce additional photon rings. Hence, precisely characterizing multi-ring images is a promising pathway for measuring black hole parameters, such as spin, as well as exploring non-Kerr spacetimes. Here, we examine the interferometric response of multi-ring systems using a series of 1) simple geometric toy models, 2) synthetic BHEX and EHT observations of geometric models, and 3) semi-analytic accretion models with ray-tracing in the Kerr spacetime. We find that interferometric amplitude is more sensitive to the shape of the photon ring, while interferometric phase is more sensitive to its displacement, which is most sensitive to black hole spin. We find that for models similar to Messier 87* (M87*), the relative displacement of the first strongly lensed image from the weakly lensed direct image is approximately $1\,μ{\rm as}$ per unit dimensionless spin, yielding an expected phase signature on a 25 G$λ$ baseline of $\sim44^\circ$ per unit spin.

preprint2022arXiv

Millimeter light curves of Sagittarius A* observed during the 2017 Event Horizon Telescope campaign

The Event Horizon Telescope (EHT) observed the compact radio source, Sagittarius A* (Sgr A*), in the Galactic Center on 2017 April 5-11 in the 1.3 millimeter wavelength band. At the same time, interferometric array data from the Atacama Large Millimeter/submillimeter Array and the Submillimeter Array were collected, providing Sgr A* light curves simultaneous with the EHT observations. These data sets, complementing the EHT very-long-baseline interferometry, are characterized by a cadence and signal-to-noise ratio previously unattainable for Sgr A* at millimeter wavelengths, and they allow for the investigation of source variability on timescales as short as a minute. While most of the light curves correspond to a low variability state of Sgr A*, the April 11 observations follow an X-ray flare, and exhibit strongly enhanced variability. All of the light curves are consistent with a red noise process, with a power spectral density (PSD) slope measured to be between -2 and -3 on timescales between 1 min and several hours. Our results indicate a steepening of the PSD slope for timescales shorter than 0.3 h. The spectral energy distribution is flat at 220 GHz and there are no time-lags between the 213 and 229 GHz frequency bands, suggesting low optical depth for the event horizon scale source. We characterize Sgr A*'s variability, highlighting the different behavior observed just after the X-ray flare, and use Gaussian process modeling to extract a decorrelation timescale and a PSD slope. We also investigate the systematic calibration uncertainties by analyzing data from independent data reduction pipelines.

preprint2022arXiv

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

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

preprint2022arXiv

Toward determining the number of observable supermassive black hole shadows

We present estimates for the number of shadow-resolved supermassive black hole (SMBH) systems that can be detected using radio interferometers, as a function of angular resolution, flux density sensitivity, and observing frequency. Accounting for the distribution of SMBHs across mass, redshift, and accretion rate, we use a new semi-analytic spectral energy distribution model to derive the number of SMBHs with detectable and optically thin horizon-scale emission. We demonstrate that (sub)millimeter interferometric observations with ${\sim}0.1$ $μ$as resolution and ${\sim}1$ $μ$Jy sensitivity could access ${>}10^6$ SMBH shadows. We then further decompose the shadow source counts into the number of black holes for which we could expect to observe the first- and second-order lensed photon rings. Accessing the bulk population of first-order photon rings requires ${\lesssim}2$ $μ$as resolution and ${\lesssim}0.5$ mJy sensitivity, while doing the same for second-order photon rings requires ${\lesssim}0.1$ $μ$as resolution and ${\lesssim}5$ $μ$Jy sensitivity. Our model predicts that with modest improvements to sensitivity, as many as $\sim$5 additional horizon-resolved sources should become accessible to the current Event Horizon Telescope (EHT), while a next-generation EHT observing at 345 GHz should have access to ${\sim}$3 times as many sources. More generally, our results can help guide enhancements of current arrays and specifications for future interferometric experiments that aim to spatially resolve a large population of SMBH shadows or higher-order photon rings.

preprint2021arXiv

Evaluation of New Submillimeter VLBI Sites for the Event Horizon Telescope

The Event Horizon Telescope (EHT) is a very long baseline interferometer built to image supermassive black holes on event-horizon scales. In this paper, we investigate candidate sites for an expanded EHT array with improved imaging capabilities. We use historical meteorology and radiative transfer analysis to evaluate site performance. Most of the existing sites in the EHT array have median zenith opacity less than 0.2 at 230 GHz during the March/April observing season. Seven of the existing EHT sites have 345 GHz opacity less than 0.5 during observing months. Out of more than forty candidate new locations analyzed, approximately half have 230 GHz opacity comparable to the existing EHT sites, and at least seventeen of the candidate sites would be comparably good for 345 GHz observing. A group of new sites with favorable transmittance and geographic placement leads to greatly enhanced imaging and science on horizon scales.

preprint2021arXiv

The intrinsic structure of Sagittarius A* at 1.3 cm and 7 mm

Sagittarius A* (Sgr A*), the Galactic Center supermassive black hole (SMBH), is one of the best targets to resolve the innermost region of SMBH with very long baseline interferometry (VLBI). In this study, we have carried out observations toward Sgr A* at 1.349 cm (22.223 GHz) and 6.950 mm (43.135 GHz) with the East Asian VLBI Network, as a part of the multi-wavelength campaign of the Event Horizon Telescope (EHT) in 2017 April. To mitigate scattering effects, the physically motivated scattering kernel model from Psaltis et al. (2018) and the scattering parameters from Johnson et al. (2018) have been applied. As a result, a single, symmetric Gaussian model well describes the intrinsic structure of Sgr A* at both wavelengths. From closure amplitudes, the major-axis sizes are ~704$\pm$102 $μ$as (axial ratio $\sim$1.19$^{+0.24}_{-0.19}$) and $\sim$300$\pm$25 $μ$as (axial ratio $\sim$1.28$\pm$0.2) at 1.349 cm and 6.95 mm respectively. Together with a quasi-simultaneous observation at 3.5 mm (86 GHz) by Issaoun et al. (2019), we show that the intrinsic size scales with observing wavelength as a power-law, with an index $\sim$1.2$\pm$0.2. Our results also provide estimates of the size and compact flux density at 1.3 mm, which can be incorporated into the analysis of the EHT observations. In terms of the origin of radio emission, we have compared the intrinsic structures with the accretion flow scenario, especially the radiatively inefficient accretion flow based on the Keplerian shell model. With this, we show that a nonthermal electron population is necessary to reproduce the source sizes.

preprint2020arXiv

Closure statistics in interferometric data

Interferometric visibilities, reflecting the complex correlations between signals recorded at antennas in an interferometric array, carry information about the angular structure of a distant source. While unknown antenna gains in both amplitude and phase can prevent direct interpretation of these measurements, certain combinations of visibilities called closure phases and closure amplitudes are independent of antenna gains and provide a convenient set of robust observables. However, these closure quantities have subtle noise properties and are generally both linearly and statistically dependent. These complications have obstructed the proper use of closure quantities in interferometric analysis, and they have obscured the relationship between analysis with closure quantities and other analysis techniques such as self calibration. We review the statistics of closure quantities, noting common pitfalls that arise when approaching low signal-to-noise due to the nonlinear propagation of statistical errors. We then develop a strategy for isolating and fitting to the independent degrees of freedom captured by the closure quantities through explicit construction of linearly independent sets of quantities along with their noise covariance in the Gaussian limit, valid for moderate signal-to-noise, and we demonstrate that model fits have biased posteriors when this covariance is ignored. Finally, we introduce a unified procedure for fitting to both closure information and partially calibrated visibilities, and we demonstrate both analytically and numerically the direct equivalence of inference based on closure quantities to that based on self calibration of complex visibilities with unconstrained antenna gains.

preprint2020arXiv

Detection of pulses from the Vela pulsar at millimeter wavelengths with phased ALMA

We report on the first detection of pulsed radio emission from a radio pulsar with the ALMA telescope. The detection was made in the Band-3 frequency range (85-101 GHz) using ALMA in the phased-array mode developed for VLBI observations. A software pipeline has been implemented to enable a regular pulsar observing mode in the future. We describe the pipeline and demonstrate the capability of ALMA to perform pulsar timing and searching. We also measure the flux density and polarization properties of the Vela pulsar (PSR J0835$-$4510) at mm-wavelengths, providing the first polarimetric study of any ordinary pulsar at frequencies above 32 GHz. Finally, we discuss the lessons learned from the Vela observations for future pulsar studies with ALMA, particularly for searches near the supermassive black hole in the Galactic Center, and the potential of using pulsars for polarization calibration of ALMA.

preprint2020arXiv

On the approximation of the black hole shadow with a simple polar curve

A black hole embedded within a bright, optically thin emitting region imprints a nearly circular "shadow" on its image, corresponding to the observer's line-of-sight into the black hole. The shadow boundary depends on the black hole's mass and spin, providing an observable signature of both properties via high resolution images. However, standard expressions for the shadow boundary are most naturally parametrized by Boyer-Lindquist radii rather than by image coordinates. We explore simple, approximate parameterizations for the shadow boundary using ellipses and a family of curves known as limacons. We demonstrate that these curves provide excellent and efficient approximations for all black hole spins and inclinations. In particular, we show that the two parameters of the limacon naturally account for the three primary shadow deformations resulting from mass and spin: size, displacement, and asymmetry. These curves are convenient for parametric model fitting directly to interferometric data, they reveal the degeneracies expected when estimating black hole properties from images with practical measurement limitations, and they provide a natural framework for parametric tests of the Kerr metric using black hole images.

preprint2020arXiv

Universal Interferometric Signatures of a Black Hole's Photon Ring

The Event Horizon Telescope image of the supermassive black hole in the galaxy M87 is dominated by a bright, unresolved ring. General relativity predicts that embedded within this image lies a thin "photon ring," which is composed of an infinite sequence of self-similar subrings that are indexed by the number of photon orbits around the black hole. The subrings approach the edge of the black hole "shadow," becoming exponentially narrower but weaker with increasing orbit number, with seemingly negligible contributions from high order subrings. Here, we show that these subrings produce strong and universal signatures on long interferometric baselines. These signatures offer the possibility of precise measurements of black hole mass and spin, as well as tests of general relativity, using only a sparse interferometric array.

preprint2020arXiv

Universal Polarimetric Signatures of the Black Hole Photon Ring

Black hole images present an annular region of enhanced brightness. In the absence of propagation effects, this "photon ring" has universal features that are completely governed by general relativity and independent of the details of the emission. Here, we show that the polarimetric image of a black hole also displays universal properties. In particular, the photon ring exhibits a self-similar pattern of polarization that encodes the black hole spin. We explore the corresponding universal polarimetric signatures of the photon ring on long interferometric baselines, and propose a method for measuring the black hole spin using a sparse interferometric array. These signatures could enable spin measurements of the supermassive black hole in M87, as well as precision tests of general relativity in the strong field regime, via a future extension of the Event Horizon Telescope to space.

preprint2019arXiv

Detailed Abundances in the Ultra-faint Magellanic Satellites Carina II and III

We present the first detailed elemental abundances in the ultra-faint Magellanic satellite galaxies Carina II (Car II) and Carina III (Car III). With high-resolution Magellan/MIKE spectroscopy, we determined abundances of nine stars in Car II including the first abundances of an RR Lyrae star in an ultra-faint dwarf galaxy; and two stars in Car III. The chemical abundances demonstrate that both systems are clearly galaxies and not globular clusters. The stars in these galaxies mostly display abundance trends matching those of other similarly faint dwarf galaxies: enhanced but declining [alpha/Fe] ratios, iron-peak elements matching the stellar halo, and unusually low neutron-capture element abundances. One star displays a low outlying [Sc/Fe] = -1.0. We detect a large Ba scatter in Car II, likely due to inhomogeneous enrichment by low-mass AGB star winds. The most striking abundance trend is for [Mg/Ca] in Car II, which decreases from +0.4 to -0.4 and indicates clear variation in the initial progenitor masses of enriching core-collapse supernovae. So far, the only ultra-faint dwarf galaxies displaying a similar [Mg/Ca] trend are likely satellites of the Large Magellanic Cloud. We find two stars with [Fe/H] < -3.5, whose abundances likely trace the first generation of metal-free Population III stars and are well-fit by Population III core-collapse supernova yields. An appendix describes our new abundance uncertainty analysis that propagates line-by-line stellar parameter uncertainties.

preprint2016arXiv

Computational Imaging for VLBI Image Reconstruction

Very long baseline interferometry (VLBI) is a technique for imaging celestial radio emissions by simultaneously observing a source from telescopes distributed across Earth. The challenges in reconstructing images from fine angular resolution VLBI data are immense. The data is extremely sparse and noisy, thus requiring statistical image models such as those designed in the computer vision community. In this paper we present a novel Bayesian approach for VLBI image reconstruction. While other methods often require careful tuning and parameter selection for different types of data, our method (CHIRP) produces good results under different settings such as low SNR or extended emission. The success of our method is demonstrated on realistic synthetic experiments as well as publicly available real data. We present this problem in a way that is accessible to members of the community, and provide a dataset website (vlbiimaging.csail.mit.edu) that facilitates controlled comparisons across algorithms.

preprint2016arXiv

Extreme Brightness Temperatures and Refractive Substructure in 3C273 with RadioAstron

Earth-space interferometry with RadioAstron provides the highest direct angular resolution ever achieved in astronomy at any wavelength. RadioAstron detections of the classic quasar 3C273 on interferometric baselines up to 171,000 km suggest brightness temperatures exceeding expected limits from the "inverse-Compton catastrophe" by two orders of magnitude. We show that at 18 cm, these estimates most probably arise from refractive substructure introduced by scattering in the interstellar medium. We use the scattering properties to estimate an intrinsic brightness temperature of 7*10^12 K, which is consistent with expected theoretical limits, but which is ~15 times lower than estimates that neglect substructure. At 6.2 cm, the substructure influences the measured values appreciably but gives an estimated brightness temperature that is comparable to models that do not account for the substructure. At 1.35 cm, the substructure does not affect the extremely high inferred brightness temperatures, in excess of 10^13 K. We also demonstrate that for a source having a Gaussian surface brightness profile, a single long-baseline estimate of refractive substructure determines an absolute minimum brightness temperature, if the scattering properties along a given line of sight are known, and that this minimum accurately approximates the apparent brightness temperature over a wide range of total flux densities.

preprint2016arXiv

High Resolution Linear Polarimetric Imaging for the Event Horizon Telescope

Images of the linear polarization of synchrotron radiation around Active Galactic Nuclei (AGN) identify their projected magnetic field lines and provide key data for understanding the physics of accretion and outflow from supermassive black holes. The highest resolution polarimetric images of AGN are produced with Very Long Baseline Interferometry (VLBI). Because VLBI incompletely samples the Fourier transform of the source image, any image reconstruction that fills in unmeasured spatial frequencies will not be unique and reconstruction algorithms are required. In this paper, we explore extensions of the Maximum Entropy Method (MEM) to linear polarimetric VLBI imaging. In contrast to previous work, our polarimetric MEM algorithm combines a Stokes I imager that uses only bispectrum measurements that are immune to atmospheric phase corruption with a joint Stokes Q and U imager that operates on robust polarimetric ratios. We demonstrate the effectiveness of our technique on 7- and 3-mm wavelength quasar observations from the VLBA and simulated 1.3-mm Event Horizon Telescope observations of Sgr A* and M87. Consistent with past studies, we find that polarimetric MEM can produce superior resolution compared to the standard CLEAN algorithm when imaging smooth and compact source distributions. As an imaging framework, MEM is highly adaptable, allowing a range of constraints on polarization structure. Polarimetric MEM is thus an attractive choice for image reconstruction with the EHT.

preprint2016arXiv

Interstellar Scintillation and the Radio Counterpart of the Fast Radio Burst FRB150418

Keane et al. (2016) have recently reported the discovery of a new fast radio burst, FRB150418, with a promising radio counterpart at 5.5 and 7.5 GHz -- a rapidly decaying source, falling from 200-300 $μ$Jy to 100 $μ$Jy on timescales of $\sim$6 d. This transient source may be associated with an elliptical galaxy at redshift $z=0.492$, providing the first firm spectroscopic redshift for a FRB and the ability to estimate the density of baryons in the intergalactic medium via the combination of known redshift and radio dispersion of the FRB. An alternative explanation, first suggested by Williams & Berger (2016b), is that the identified counterpart may instead be a compact AGN. The putative counterpart's variation may then instead be extrinsic, caused by refractive scintillation in the ionized interstellar medium of the Milky Way, which would invalidate the association with FRB150418. We examine this latter explanation in detail and show that the reported observations are consistent with scintillating radio emission from the core of a radio-loud active galactic nucleus (AGN) having a brightness temperature $T_{\rm b} \gtrsim 10^9 {\rm K}$. Using numerical simulations of the expected scattering for the line of sight to FRB150418, we provide example images and light curves of such an AGN at 5.5 and 7.5 GHz. These results can be compared with continued radio monitoring to conclusively determine the importance of scintillation for the observed radio variability, and they show that scintillation is a critical consideration for continued searches for FRB counterparts at radio wavelengths.

preprint2016arXiv

Modeling Seven Years of Event Horizon Telescope Observations with Radiatively Inefficient Accretion Flow Models

An initial three-station version of the Event Horizon Telescope, a millimeter-wavelength very-long baseline interferometer, has observed Sagittarius A* (Sgr A*) repeatedly from 2007 to 2013, resulting in the measurement of a variety of interferometric quantities. Of particular importance, there is now a large set of closure phases, measured over a number of independent observing epochs. We analyze these observations within the context of a realization of semi-analytic radiatively inefficient disk models, implicated by the low luminosity of Sgr A*. We find a broad consistency among the various observing epochs and between different interferometric data types, with the latter providing significant support for this class of models of Sgr A*. The new data significantly tighten existing constraints on the spin magnitude and its orientation within this model context, finding a spin magnitude of $a=0.10^{+0.30+0.56}_{-0.10-0.10}$, an inclination with respect to the line of sight of $θ={60^\circ}^{+5^\circ+10^\circ}_{-8^\circ-13^\circ}$, and a position angle of $ξ={156^\circ}^{+10^\circ+14^\circ}_{-17^\circ-27^\circ}$ east of north. These are in good agreement with previous analyses. Notably, the previous $180^\circ$ degeneracy in the position angle has now been conclusively broken by the inclusion of the closure phase measurements. A reflection degeneracy in the inclination remains, permitting two localizations of the spin vector orientation, one of which is in agreement with the orbital angular momentum of the infrared gas cloud G2 and the clockwise disk of young stars. This possibly supports a relationship between Sgr A*'s accretion flow and these larger-scale features.

preprint2016arXiv

Observing---and Imaging---Active Galactic Nuclei with the Event Horizon Telescope

Originally developed to image the shadow region of the central black hole in Sagittarius A* and in the nearby galaxy M87, the Event Horizon Telescope (EHT) provides deep, very high angular resolution data on other AGN sources too. The challenges of working with EHT data have spurred the development of new image reconstruction algorithms. This work briefly reviews the status of the EHT and its utility for observing AGN sources, with emphasis on novel imaging techniques that offer the promise of better reconstructions at 1.3 mm and other wavelengths.

preprint2016arXiv

Persistent Asymmetric Structure of Sagittarius A* on Event Horizon Scales

The Galactic Center black hole Sagittarius A* (Sgr A*) is a prime observing target for the Event Horizon Telescope (EHT), which can resolve the 1.3 mm emission from this source on angular scales comparable to that of the general relativistic shadow. Previous EHT observations have used visibility amplitudes to infer the morphology of the millimeter-wavelength emission. Potentially much richer source information is contained in the phases. We report on 1.3 mm phase information on Sgr A* obtained with the EHT on a total of 13 observing nights over 4 years. Closure phases, the sum of visibility phases along a closed triangle of interferometer baselines, are used because they are robust against phase corruptions introduced by instrumentation and the rapidly variable atmosphere. The median closure phase on a triangle including telescopes in California, Hawaii, and Arizona is nonzero. This result conclusively demonstrates that the millimeter emission is asymmetric on scales of a few Schwarzschild radii and can be used to break 180-degree rotational ambiguities inherent from amplitude data alone. The stability of the sign of the closure phase over most observing nights indicates persistent asymmetry in the image of Sgr A* that is not obscured by refraction due to interstellar electrons along the line of sight.

preprint2016arXiv

Stochastic Optics: A Scattering Mitigation Framework for Radio Interferometric Imaging

Just as turbulence in the Earth's atmosphere can severely limit the angular resolution of optical telescopes, turbulence in the ionized interstellar medium fundamentally limits the resolution of radio telescopes. We present a scattering mitigation framework for radio imaging with very long baseline interferometry (VLBI) that partially overcomes this limitation. Our framework, "stochastic optics," derives from a simplification of strong interstellar scattering to separate small-scale ("diffractive") effects from large-scale ("refractive") effects, thereby separating deterministic and random contributions to the scattering. Stochastic optics extends traditional synthesis imaging by simultaneously reconstructing an unscattered image and its refractive perturbations. Its advantages over direct imaging come from utilizing the many deterministic properties of the scattering -- such as the time-averaged "blurring," polarization independence, and the deterministic evolution in frequency and time -- while still accounting for the stochastic image distortions on large scales. These distortions are identified in the image reconstructions through regularization by their time-averaged power spectrum. Using synthetic data, we show that this framework effectively removes the blurring from diffractive scattering while reducing the spurious image features from refractive scattering. Stochastic optics can provide significant improvements over existing scattering mitigation strategies and is especially promising for imaging the Galactic Center supermassive black hole, Sagittarius A*, with the Global mm-VLBI Array and with the Event Horizon Telescope.

preprint2016arXiv

Testing General Relativity with the Shadow Size of Sgr A*

In general relativity, the angular radius of the shadow of a black hole is primarily determined by its mass-to-distance ratio and depends only weakly on its spin and inclination. If general relativity is violated, however, the shadow size may also depend strongly on parametric deviations from the Kerr metric. Based on a reconstructed image of Sagittarius A* (Sgr A*) from a simulated one-day observing run of a seven-station Event Horizon Telescope (EHT) array, we employ a Markov chain Monte Carlo algorithm to demonstrate that such an observation can measure the angular radius of the shadow of Sgr A* with an uncertainty of ~1.5 uas (6%). We show that existing mass and distance measurements can be improved significantly when combined with upcoming EHT measurements of the shadow size and that tight constraints on potential deviations from the Kerr metric can be obtained.

preprint2016arXiv

The Intrinsic Shape of Sagittarius A* at 3.5-mm Wavelength

The radio emission from Sgr A$^\ast$ is thought to be powered by accretion onto a supermassive black hole of $\sim\! 4\times10^6~ \rm{M}_\odot$ at the Galactic Center. At millimeter wavelengths, Very Long Baseline Interferometry (VLBI) observations can directly resolve the bright innermost accretion region of Sgr A$^\ast$. Motivated by the addition of many sensitive, long baselines in the north-south direction, we developed a full VLBI capability at the Large Millimeter Telescope Alfonso Serrano (LMT). We successfully detected Sgr A$^\ast$ at 3.5~mm with an array consisting of 6 Very Long Baseline Array telescopes and the LMT. We model the source as an elliptical Gaussian brightness distribution and estimate the scattered size and orientation of the source from closure amplitude and self-calibration analysis, obtaining consistent results between methods and epochs. We then use the known scattering kernel to determine the intrinsic two dimensional source size at 3.5 mm: $(147\pm7~μ\rm{as}) \times (120\pm12~μ\rm{as})$, at position angle $88^\circ\pm7^\circ$ east of north. Finally, we detect non-zero closure phases on some baseline triangles, but we show that these are consistent with being introduced by refractive scattering in the interstellar medium and do not require intrinsic source asymmetry to explain.

preprint2016arXiv

The Optics of Refractive Substructure

Newly recognized effects of refractive scattering in the ionized interstellar medium have broad implications for very long baseline interferometry (VLBI) at extreme angular resolutions. Building upon work by Blandford & Narayan (1985), we present a simplified, geometrical optics framework, which enables rapid, semi-analytic estimates of refractive scattering effects. We show that these estimates exactly reproduce previous results based on a more rigorous statistical formulation. We then derive new expressions for the scattering-induced fluctuations of VLBI observables such as closure phase, and we demonstrate how to calculate the fluctuations for arbitrary quantities of interest using a Monte Carlo technique.

preprint2016arXiv

Theory and Simulations of Refractive Substructure in Resolved Scatter-Broadened Images

At radio wavelengths, scattering in the interstellar medium distorts the appearance of astronomical sources. Averaged over a scattering ensemble, the result is a blurred image of the source. However, Narayan & Goodman (1989) and Goodman & Narayan (1989) showed that for an incomplete average, scattering introduces refractive substructure in the image of a point source that is both persistent and wideband. We show that this substructure is quenched but not smoothed by an extended source. As a result, when the scatter-broadening is comparable to or exceeds the unscattered source size, the scattering can introduce spurious compact features into images. In addition, we derive efficient strategies to numerically compute realistic scattered images, and we present characteristic examples from simulations. Our results show that refractive substructure is an important consideration for ongoing missions at the highest angular resolutions, and we discuss specific implications for RadioAstron and the Event Horizon Telescope.

preprint2015arXiv

230 GHz VLBI observations of M87: event-horizon-scale structure at the enhanced very-high-energy $\rm γ$-ray state in 2012

We report on 230 GHz (1.3 mm) VLBI observations of M87 with the Event Horizon Telescope using antennas on Mauna Kea in Hawaii, Mt. Graham in Arizona and Cedar Flat in California. For the first time, we have acquired 230 GHz VLBI interferometric phase information on M87 through measurement of closure phase on the triangle of long baselines. Most of the measured closure phases are consistent with 0$^{\circ}$ as expected by physically-motivated models for 230 GHz structure such as jet models and accretion disk models. The brightness temperature of the event-horizon-scale structure is $\sim 1 \times 10^{10}$ K derived from the compact flux density of $\sim 1$ Jy and the angular size of $\sim 40 $ $\rm μ$as $\sim$ 5.5 $R_{\rm s}$, which is broadly consistent with the peak brightness of the radio cores at 1-86 GHz located within $\sim 10^2$ $R_{\rm s}$. Our observations occurred in the middle of an enhancement in very-high-energy (VHE) $\rm γ$-ray flux, presumably originating in the vicinity of the central black hole. Our measurements, combined with results of multi-wavelength observations, favor a scenario in which the VHE region has an extended size of $\sim$20-60 $R_{\rm s}$.

preprint2015arXiv

Measuring the Direction and Angular Velocity of a Black Hole Accretion Disk via Lagged Interferometric Covariance

We show that interferometry can be applied to study irregular, rapidly rotating structures, as are expected in the turbulent accretion flow near a black hole. Specifically, we analyze the lagged covariance between interferometric baselines of similar lengths but slightly different orientations. For a flow viewed close to face-on, we demonstrate that the peak in the lagged covariance indicates the direction and angular velocity of the emission pattern from the flow. Even for moderately inclined flows, the covariance robustly estimates the flow direction, although the estimated angular velocity can be significantly biased. Importantly, measuring the direction of the flow as clockwise or counterclockwise on the sky breaks a degeneracy in accretion disk inclinations when analyzing time-averaged images alone. We explore the potential efficacy of our technique using three-dimensional, general relativistic magnetohydrodynamic (GRMHD) simulations, and we highlight several baseline pairs for the Event Horizon Telescope (EHT) that are well-suited to this application. These results indicate that the EHT may be capable of estimating the direction and angular velocity of the emitting material near Sagittarius A*, and they suggest that a rotating flow may even be utilized to improve imaging capabilities.

preprint2015arXiv

Resolved Magnetic-Field Structure and Variability Near the Event Horizon of Sagittarius A*

Near a black hole, differential rotation of a magnetized accretion disk is thought to produce an instability that amplifies weak magnetic fields, driving accretion and outflow. These magnetic fields would naturally give rise to the observed synchrotron emission in galaxy cores and to the formation of relativistic jets, but no observations to date have been able to resolve the expected horizon-scale magnetic-field structure. We report interferometric observations at 1.3-millimeter wavelength that spatially resolve the linearly polarized emission from the Galactic Center supermassive black hole, Sagittarius A*. We have found evidence for partially ordered fields near the event horizon, on scales of ~6 Schwarzschild radii, and we have detected and localized the intra-hour variability associated with these fields.

preprint2015arXiv

Unbound Debris Streams and Remnants Resulting From the Tidal Disruptions of Stars by Supermassive Black Holes

The kinetic energy of a star in orbit about a supermassive black hole is a significant fraction of its rest mass energy when its periapse is comparable to its tidal radius. Upon its destruction, a fraction of this energy is extracted and injected into the stellar debris, half of which becomes unbound from the black hole, with the fastest material moving at $\sim 0.03 c$. In this paper, we present a formalism for determining the fate of these unbound debris streams (UDSs) as they depart from the black hole and interact with the surrounding gas. As the density and velocity varies along the length of a UDS, we find that hydrodynamical drag quickly shapes UDSs into loop-like structures, with the densest portions of the streams leading portions of lower density. As UDSs travel outwards, their drag against the ISM increases quadratically with distance, which causes UDSs to deposit their momentum and energy into the ambient medium before the surrounding shocked ISM has a chance to cool. This sudden injection of $\sim 10^{50}$ erg into the ambient medium generates a Sedov-like unbound debris remnant (UDR) that mimics supernova remnants (SNRs) in energetics and appearance, accelerates particles which will produce cosmic rays and synchrotron emission, and provides momentum feedback into the molecular clouds surrounding a black hole. We estimate that a few of these UDRs might be present within a couple degrees of the Galactic Center masquerading as SNRs, and that the UDR scenario is a plausible explanation for Sgr A East.

preprint2014arXiv

Imaging an Event Horizon: Mitigation of Scattering Toward Sagittarius A*

The image of the emission surrounding the black hole in the center of the Milky Way is predicted to exhibit the imprint of general relativistic (GR) effects, including the existence of a shadow feature and a photon ring of diameter ~50 microarcseconds. Structure on these scales can be resolved by millimeter-wavelength very long baseline interferometry (VLBI). However, strong-field GR features of interest will be blurred at lambda >= 1.3 mm due to scattering by interstellar electrons. The scattering properties are well understood over most of the relevant range of baseline lengths, suggesting that the scattering may be (mostly) invertible. We simulate observations of a model image of Sgr A* and demonstrate that the effects of scattering can indeed be mitigated by correcting the visibilities before reconstructing the image. This technique is also applicable to Sgr A* at longer wavelengths.

preprint2014arXiv

Relative Astrometry of Compact Flaring Structures in Sgr A* with Polarimetric VLBI

We demonstrate that polarimetric interferometry can be used to extract precise spatial information about compact polarized flares of Sgr A*. We show that, for a faint dynamical component, a single interferometric baseline suffices to determine both its polarization and projected displacement from the quiescent intensity centroid. A second baseline enables two-dimensional reconstruction of the displacement, and additional baselines can self-calibrate using the flare, enhancing synthesis imaging of the quiescent emission. We apply this technique to simulated 1.3-mm wavelength observations of a "hot spot" embedded in a radiatively inefficient accretion disk around Sgr A*. Our results indicate that, even with current sensitivities, polarimetric interferometry with the Event Horizon Telescope can achieve ~5 microarcsecond relative astrometry of compact flaring structures near Sgr A* on timescales of minutes.

preprint2013arXiv

Interferometric Visibility of a Scintillating Source: Statistics at the Nyquist Limit

We derive the distribution of interferometric visibility for a source exhibiting strong diffractive scintillation, with particular attention to spectral resolution at or near the Nyquist limit. We also account for arbitrary temporal averaging, intrinsic variability within the averaging time, and the possibility of spatially-extended source emission. We demonstrate that the interplay between scintillation and self-noise induces several remarkable features, such as a broad "skirt" in the visibility distribution. Our results facilitate the interpretation of interferometric observations of pulsars at meter and decimeter wavelengths.

preprint2012arXiv

Constraining the Vela Pulsar's Radio Emission Region Using Nyquist-Limited Scintillation Statistics

Using a novel technique, we achieve ~100 picoarcsecond resolution and set an upper bound of less than 4 km for the characteristic size of the Vela pulsar's emission region. Specifically, we analyze flux-density statistics of the Vela pulsar at 760 MHz. Because the pulsar exhibits strong diffractive scintillation, these statistics convey information about the spatial extent of the radio emission region. We measure both a characteristic size of the emission region and the emission sizes for individual pulses. Our results imply that the radio emission altitude for the Vela pulsar at this frequency is less than 340 km.

preprint2012arXiv

Optimal Correlation Estimators for Quantized Signals

Using a maximum-likelihood criterion, we derive optimal correlation strategies for signals with and without digitization. We assume that the signals are drawn from zero-mean Gaussian distributions, as is expected in radio-astronomical applications, and we present correlation estimators both with and without a priori knowledge of the signal variances. We demonstrate that traditional estimators of correlation, which rely on averaging products, exhibit large and paradoxical noise when the correlation is strong. However, we also show that these estimators are fully optimal in the limit of vanishing correlation. We calculate the bias and noise in each of these estimators and discuss their suitability for implementation in modern digital correlators.

preprint2012arXiv

Ultra-High Resolution Intensity Statistics of a Scintillating Source

We derive the distribution of flux density of a compact source exhibiting strong diffractive scintillation. Our treatment accounts for arbitrary spectral averaging, spatially-extended source emission, and the possibility of intrinsic variability within the averaging time, as is typical for pulsars. We also derive the modulation index and present a technique for estimating the self-noise of the distribution, which can be used to identify amplitude variations on timescales shorter than the spectral accumulation time. Our results enable a for direct comparison with ultra-high resolution observations of pulsars, particularly single-pulse studies with Nyquist-limited resolution, and can be used to identify the spatial emission structure of individual pulses at a small fraction of the diffractive scale.

preprint2011arXiv

Noise and Signal for Spectra of Intermittent Noiselike Emission

We show that intermittency of noiselike emission, after propagation through a scattering medium, affects the distribution of noise in the observed correlation function. Intermittency also affects correlation of noise among channels of the spectrum, but leaves the average spectrum, average correlation function, and distribution of noise among channels of the spectrum unchanged. Pulsars are examples of such sources: intermittent and affected by interstellar propagation. We assume that the source emits Gaussian white noise, modulated by a time-envelope. Propagation convolves the resulting time series with an impulse-response function that represents effects of dispersion, scattering, and absorption. We assume that this propagation kernel is shorter than the time for an observer to accumulate a single spectrum. We show that rapidly-varying intermittent emission tends to concentrate noise near the central lag of the correlation function. We derive mathematical expressions for this effect, in terms of the time envelope and the propagation kernel. We present examples, discuss effects of background noise, and compare our results with observations.