Search

Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U , and V respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties., Comment: Accepted for publication in ApJ

In 1977, Blandford and Znajek showed that the electromagnetic field surrounding a rotating black hole can harvest its spin energy and use it to power a collimated astrophysical jet, such as the one launched from the center of the elliptical galaxy M87. Today, interferometric observations with the Event Horizon Telescope (EHT) are delivering high-resolution, event-horizon-scale, polarimetric images of the supermassive black hole M87* at the jet launching point. These polarimetric images offer an unprecedented window into the electromagnetic field structure around a black hole. In this paper, we show that a simple polarimetric observable that quantifies the magnetic field helicity -- the sign of $\angle\beta_2$ in a near-horizon image -- depends on the sign of the electromagnetic energy flux and therefore provides a direct probe of black hole energy extraction. In Boyer-Lindquist coordinates, the Poynting flux for axisymmetric electromagnetic fields is proportional to the product $B^\phi B^r$. The polarimetric observable $\angle\beta_2$ likewise depends on the ratio $B^\phi/B^r$, thereby enabling an observer to experimentally determine the direction of electromagnetic energy flow in the near-horizon environment. Data from the 2017 EHT observations of M87* are consistent with electromagnetic energy outflow. Currently envisioned multi-frequency observations of M87* will achieve higher dynamic range and angular resolution, and hence deliver measurements of $\angle\beta_2$ closer to the event horizon as well as better constraints on Faraday rotation. Such observations will enable a definitive test for energy extraction from the black hole M87*., Comment: 34 pages, 5 figures. Submitted to ApJL

The Event Horizon Telescope (EHT) has produced images of the plasma flow around the supermassive black holes in Sgr A* and M87* with a resolution comparable to the projected size of their event horizons. Observations with the next-generation Event Horizon Telescope (ngEHT) will have significantly improved Fourier plane coverage and will be conducted at multiple frequency bands (86, 230, and 345 GHz), each with a wide bandwidth. At these frequencies, both Sgr A* and M87* transition from optically thin to optically thick. Resolved spectral index maps in the near-horizon and jet-launching regions of these supermassive black hole sources can constrain properties of the emitting plasma that are degenerate in single-frequency images. In addition, combining information from data obtained at multiple frequencies is a powerful tool for interferometric image reconstruction, since gaps in spatial scales in single-frequency observations can be filled in with information from other frequencies. Here we present a new method of simultaneously reconstructing interferometric images at multiple frequencies along with their spectral index maps. The method is based on existing Regularized Maximum Likelihood (RML) methods commonly used for EHT imaging and is implemented in the eht-imaging Python software library. We show results of this method on simulated ngEHT data sets as well as on real data from the VLBA and ALMA. These examples demonstrate that simultaneous RML multi-frequency image reconstruction produces higher-quality and more scientifically useful results than is possible from combining independent image reconstructions at each frequency., Comment: 25 pages, 15 figures. Accepted to ApJ

We introduce a new library of 535,194 model images of the supermassive black holes and Event Horizon Telescope (EHT) targets Sgr A* and M87*, computed by performing general relativistic radiative transfer calculations on general relativistic magnetohydrodynamics simulations. Then, to infer underlying black hole and accretion flow parameters (spin, inclination, ion-to-electron temperature ratio, and magnetic field polarity), we train a random forest machine learning model on various hand-picked polarimetric observables computed from each image. Our random forest is capable of making meaningful predictions of spin, inclination, and the ion-to-electron temperature ratio, but has more difficulty inferring magnetic field polarity. To disentangle how physical parameters are encoded in different observables, we apply two different metrics to rank the importance of each observable at inferring each physical parameter. Details of the spatially resolved linear polarization morphology stand out as important discriminators between models. Bearing in mind the theoretical limitations and incompleteness of our image library, for the real M87* data, our machinery favours high-spin retrograde models with large ion-to-electron temperature ratios. Due to the time-variable nature of these targets, repeated polarimetric imaging will further improve model inference as the EHT and next-generation (EHT) continue to develop and monitor their targets., Comment: 24 pages, 27 figures

We investigate general relativistic magnetohydrodynamic simulations (GRMHD) to determine the physical origin of the twisty patterns of linear polarization seen in spatially resolved black hole images and explain their morphological dependence on black hole spin. By characterising the observed emission with a simple analytic ring model, we find that the twisty morphology is determined by the magnetic field structure in the emitting region. Moreover, the dependence of this twisty pattern on spin can be attributed to changes in the magnetic field geometry that occur due to the frame dragging. By studying an analytic ring model, we find that the roles of Doppler boosting and lensing are subdominant. Faraday rotation may cause a systematic shift in the linear polarization pattern, but we find that its impact is subdominant for models with strong magnetic fields and modest ion-to-electron temperature ratios. Models with weaker magnetic fields are much more strongly affected by Faraday rotation and have more complicated emission geometries than can be captured by a ring model. However, these models are currently disfavoured by the recent EHT observations of M87*. Our results suggest that linear polarization maps can provide a probe of the underlying magnetic field structure around a black hole, which may then be usable to indirectly infer black hole spins. The generality of these results should be tested with alternative codes, initial conditions, and plasma physics prescriptions., Comment: 26 pages, 20 figures, accepted for publication in ApJ

Full list of authors: Okino, Hiroki; Akiyama, Kazunori; Asada, Keiichi; Gomez, Jose L.; Hada, Kazuhiro; Honma, Mareki; Krichbaum, Thomas P.; Kino, Motoki; Nagai, Hiroshi; Bach, Uwe; Blackburn, Lindy; Bouman, Katherine L.; Chael, Andrew; Crew, Geoffrey B.; Doeleman, Sheperd S.; Fish, Vincent L.; Goddi, Ciriaco; Issaoun, Sara; Johnson, Michael D.; Jorstad, Svetlana; Koyama, Shoko; Lonsdale, Colin J.; Lu, Ru-Sen; Marti-Vidal, Ivan; Matthews, Lynn D.; Mizuno, Yosuke; Moriyama, Kotaro; Nakamura, Masanori; Pu, Hung-Yi; Ros, Eduardo; Savolainen, Tuomas; Tazaki, Fumie; Wagner, Jan; Wielgus, Maciek; Zensus, Anton.--This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited., 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. © 2022. The Author(s). Published by the American Astronomical Society., This work was financially supported by grants from the National Science Foundation (NSF; AST-1440254; AST-1614868; AST-2034306) and the Japan Society for Promotion of Science (JSPS; JP21H01137; JP20H01951; JP19KK0081; JP18H03721.). This work has been partially supported by the Generalitat Valenciana GenT Project CIDEGENT/2018/021 and by the MICINN Research Project PID2019-108995GB-C22. This work was partially supported by FAPESP (Fundação de Amparo à Pesquisa do Estado de São Paulo) under grant 2021/01183-8. K. Akiyama is also financially supported by the following NSF grants (AST-1935980, AST-2107681, AST-2132700, OMA-2029670). A.C. was supported by Hubble Fellowship grant HST-HF2-51431.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. A.C. also gratefully acknowledges support from the Princeton Gravity Initiative. S.I. is supported by NASA Hubble Fellowship grant HST-HF2-51482.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. R.-S. L. is supported by the Max Planck Partner Group of the MPG and the CAS and acknowledges support by the Key Program of the National Natural Science Foundation of China (grant No. 11933007), the Key Research Program of Frontier Sciences, CAS (grant No. ZDBS-LY-SLH011), and the Shanghai Pilot Program for Basic Research—Chinese Academy of Science, Shanghai Branch (JCYJ-SHFY-2022-013). Y.M. acknowledges the support of the National Natural Science Foundation of China (grant No. 12273022). The black hole Initiative at Harvard University is financially supported by a grant from the John Templeton Foundation., With funding from the Spanish government through the Severo Ochoa Centre of Excellence accreditation SEV-2017-0709.

We study the effects of including a nonzero positron-to-electron fraction in emitting plasma on the polarized SEDs and sub-millimeter images of jet and accretion flow models for near-horizon emission from M87* and Sgr A*. For M87*, we consider a semi-analytic fit to the force-free plasma regions of a general relativistic magnetohydrodynamic jet simulation which we populate with power-law leptons with a constant electron-to-magnetic pressure ratio. For Sgr A*, we consider a standard self-similar radiatively inefficient accretion flow where the emission is predominantly from thermal leptons with a small fraction in a power-law tail. In both models, we fix the positron-to-electron ratio throughout the emission region. We generate polarized images and spectra from our models using the general-relativistic ray tracing and radiative transfer from GRTRANS. We find that a substantial positron fraction reduces the circular polarization fraction at infrared and higher frequencies. However, in sub-millimeter images higher positron fractions increase polarization fractions due to strong effects of Faraday conversion. We find a M87* jet model that best matches the available broadband total intensity and 230 GHz polarization data is a sub-equipartition, with positron fraction of $\simeq$ 10%. We show that jet models with significant positron fractions do not satisfy the polarimetric constraints at 230 GHz from the Event Horizon Telescope (EHT). Sgr A* models show similar trends in their polarization fractions with increasing pair fraction. Both models suggest that resolved, polarized EHT images are useful to constrain the presence of pairs at 230 GHz emitting regions of M87* and Sgr A*., Comment: 32 pages, 12 figures, accepted to ApJ

We present the results of nine simulations of radiatively-inefficient magnetically arrested disks (MADs) across different values of the black hole spin parameter $a_*$: $-0.9$, $-0.7$, $-0.5$, $-0.3$, 0, 0.3, 0.5, 0.7, and 0.9. Each simulation was run up to $t \gtrsim 100,000\,GM/c^3$ to ensure disk inflow equilibrium out to large radii. We find that the saturated magnetic flux level, and consequently also jet power, of MAD disks depends strongly on the black hole spin, confirming previous results. Prograde disks saturate at a much higher relative magnetic flux and have more powerful jets than their retrograde counterparts. MADs with spinning black holes naturally launch jets with generalized parabolic profiles whose widths vary as a power of distance from the black hole. For distances up to $100\;GM/c^2$, the power-law index is $k \approx 0.27-0.42$. There is a strong correlation between the disk-jet geometry and the dimensionless magnetic flux, resulting in prograde systems displaying thinner equatorial accretion flows near the black hole and wider jets, compared to retrograde systems. Prograde and retrograde MADs also exhibit different trends in disk variability: accretion rate variability increases with increasing spin for $a_*>0$ and remains almost constant for $a_*\lesssim 0$, while magnetic flux variability shows the opposite trend. Jets in the MAD state remove more angular momentum from black holes than is accreted, effectively spinning down the black hole. If powerful jets from MAD systems in Nature are persistent, this loss of angular momentum will notably reduce the black hole spin over cosmic time., Accepted by MNRAS 2022 January 26. Received 2022 January 26; in original form 2021 August 27

We present measurements of the properties of the central radio source in M87 using Event Horizon Telescope data obtained during the 2017 campaign. We develop and fit geometric crescent models (asymmetric rings with interior brightness depressions) using two independent sampling algorithms that consider distinct representations of the visibility data. We show that the crescent family of models is statistically preferred over other comparably complex geometric models that we explore. We calibrate the geometric model parameters using general relativistic magnetohydrodynamic (GRMHD) models of the emission region and estimate physical properties of the source. We further fit images generated from GRMHD models directly to the data. We compare the derived emission region and black hole parameters from these analyses with those recovered from reconstructed images. There is a remarkable consistency among all methods and data sets. We find that >50% of the total flux at arcsecond scales comes from near the horizon, and that the emission is dramatically suppressed interior to this region by a factor >10, providing direct evidence of the predicted shadow of a black hole. Across all methods, we measure a crescent diameter of 42 ± 3 μas and constrain its fractional width to be, The authors of this Letter thank the following organizations and programs: the Academy of Finland (projects 274477, 284495, 312496); the Advanced European Network of E-infrastructures for Astronomy with the SKA (AENEAS) project, supported by the European Commission Framework Programme Horizon 2020 Research and Innovation action under grant agreement 731016; the Alexander von Humboldt Stiftung; the Black Hole Initiative at Harvard University, through a grant (60477) from the John Templeton Foundation; the China Scholarship Council; Comisión Nacional de Investigación Científica y Tecnológica (CONICYT, Chile, via PIA ACT172033, Fondecyt 1171506, BASAL AFB170002, ALMA-conicyt 31140007); Consejo Nacional de Ciencia y Tecnología (CONACYT, Mexico, projects 104497, 275201, 279006, 281692); the Delaney Family via the Delaney Family John A. Wheeler Chair at Perimeter Institute; Dirección General de Asuntos del Personal Académico-Universidad Nacional 9 The Astrophysical Journal Letters, 875:L1 (17pp), 2019 April 10 The EHT Collaboration et al. Autónoma de México (DGAPA-UNAM, project IN112417); the European Research Council (ERC) Synergy Grant “BlackHoleCam: Imaging the Event Horizon of Black Holes” (grant 610058); the Generalitat Valenciana postdoctoral grant APOSTD/2018/177; the Gordon and Betty Moore Foundation (grants GBMF-3561, GBMF-5278); the Istituto Nazionale di Fisica Nucleare (INFN) sezione di Napoli, iniziative specifiche TEONGRAV; the International Max Planck Research School for Astronomy and Astrophysics at the Universities of Bonn and Cologne; the Jansky Fellowship program of the National Radio Astronomy Observatory (NRAO); the Japanese Government (Monbukagakusho: MEXT) Scholarship; the Japan Society for the Promotion of Science (JSPS) Grant-in-Aid for JSPS Research Fellowship (JP17J08829); JSPS Overseas Research Fellowships; the Key Research Program of Frontier Sciences, Chinese, Academy of Sciences (CAS, grants QYZDJ-SSW-SLH057, QYZDJ-SSW-SYS008); the Leverhulme Trust Early Career Research Fellowship; the Max-Planck-Gesellschaft (MPG); the Max Planck Partner Group of the MPG and the CAS; the MEXT/JSPS KAKENHI (grants 18KK0090, JP18K13594, JP18K03656, JP18H03721, 18K03709, 18H01245, 25120007); the MIT International Science and Technology Initiatives (MISTI) Funds; the Ministry of Science and Technology (MOST) of Taiwan (105-2112-M-001-025-MY3, 106-2112-M001-011, 106-2119-M-001-027, 107-2119-M-001-017, 107- 2119-M-001-020, and 107-2119-M-110-005); the National Aeronautics and Space Administration (NASA, Fermi Guest Investigator grant 80NSSC17K0649); the National Institute of Natural Sciences (NINS) of Japan; the National Key Research and Development Program of China (grant 2016YFA0400704, 2016YFA0400702); the National Science Foundation (NSF, grants AST-0096454, AST-0352953, AST-0521233, AST0705062, AST-0905844, AST-0922984, AST-1126433, AST1140030, DGE-1144085, AST-1207704, AST-1207730, AST1207752, MRI-1228509, OPP-1248097, AST-1310896, AST1312651, AST-1337663, AST-1440254, AST-1555365, AST1715061, AST-1615796, AST-1614868, AST-1716327, OISE1743747, AST-1816420); the Natural Science Foundation of China (grants 11573051, 11633006, 11650110427, 10625314, 11721303, 11725312, 11873028, 11873073, U1531245, 11473010); the Natural Sciences and Engineering Research Council of Canada (NSERC, including a Discovery Grant and the NSERC Alexander Graham Bell Canada Graduate Scholarships-Doctoral Program); the National Youth Thousand Talents Program of China; the National Research Foundation of Korea (grant 2015-R1D1A1A01056807, the Global PhD Fellowship Grant: NRF-2015H1A2A1033752, and the Korea Research Fellowship Program: NRF-2015H1D3A1066561); the Netherlands Organization for Scientific Research (NWO) VICI award (grant 639.043.513) and Spinoza Prize (SPI 78-409); the New Scientific Frontiers with Precision, Radio Interferometry Fellowship awarded by the South African Radio Astronomy Observatory (SARAO), which is a facility of the National Research Foundation (NRF), an agency of the Department of Science and Technology (DST) of South Africa; the Onsala Space Observatory (OSO) national infrastructure, for the provisioning of its facilities/observational support (OSO receives funding through the Swedish Research Council under grant 2017-00648); the Perimeter Institute for Theoretical Physics (research at Perimeter Institute is supported by the Government of Canada through the Department of Innovation, Science and Economic Development Canada and by the Province of Ontario through the Ministry of Economic Development, Job Creation and Trade); the Russian Science Foundation (grant 17-12-01029); the Spanish Ministerio de Economía y Competitividad (grants AYA2015-63939-C2-1-P, AYA2016-80889-P); the State Agency for Research of the Spanish MCIU through the “Center of Excellence Severo Ochoa” award for the Instituto de Astrofísica de Andalucía (SEV-2017-0709); the Toray Science Foundation; the US Department of Energy (USDOE) through the Los Alamos National Laboratory (operated by Triad National Security, LLC, for the National Nuclear Security Administration of the USDOE (Contract 89233218CNA000001)); the Italian Ministero dell’Istruzione Università e Ricerca through the grant Progetti Premiali 2012-iALMA (CUP C52I13000140001); the European Unionʼs Horizon 2020 research and innovation programme under grant agreement No 730562 RadioNet; ALMA North America Development Fund; Chandra TM6-17006X

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

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., Comment: 31 pages, 17 figures

Current simulations of hot accretion flows around black holes assume either a single-temperature gas or, at best, a two-temperature gas with thermal ions and electrons. However, processes like magnetic reconnection and shocks can accelerate electrons into a nonthermal distribution, which will not quickly thermalise at the very low densities found in many systems. Such nonthermal electrons have been invoked to explain the infrared and X-ray spectra and strong variability of Sagittarius A* (Sgr A*), the black hole at the Galactic Center. We present a method for self-consistent evolution of a nonthermal electron population in the GRMHD code KORAL. The electron distribution is tracked across Lorentz factor space and is evolved in space and time, in parallel with thermal electrons, thermal ions, and radiation. In the present study, for simplicity, energy injection into the nonthermal distribution is taken as a fixed fraction of the local electron viscous heating rate. Numerical results are presented for a model with a low mass accretion rate similar to that of Sgr A*. We find that the presence of a nonthermal population of electrons has negligible effect on the overall dynamics of the system. Due to our simple uniform particle injection prescription, the radiative power in the nonthermal simulation is enhanced at large radii. The energy distribution of the nonthermal electrons shows a synchrotron cooling break, with the break Lorentz factor varying with location and time, reflecting the complex interplay between the local viscous heating rate, magnetic field strength, and fluid velocity., Comment: 21 pages, 10 figures. Accepted to MNRAS

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~\mu\rm{as}) \times (120\pm12~\mu\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., Comment: 10 pages, 8 figures, accepted for publication in ApJ

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., Comment: 19 pages, 9 figures. Accepted for publication in ApJ. Imaging code available at https://github.com/achael/eht-imaging/

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., Comment: 8 Pages, 4 Figures, accepted for publication in ApJ