Showing 1–50 of 73 results for author: Wong, G

The Black Hole Explorer: Motivation and Vision

Authors: Michael D. Johnson, Kazunori Akiyama, Rebecca Baturin, Bryan Bilyeu, Lindy Blackburn, Don Boroson, Alejandro Cardenas-Avendano, Andrew Chael, Chi-kwan Chan, Dominic Chang, Peter Cheimets, Cathy Chou, Sheperd S. Doeleman, Joseph Farah, Peter Galison, Ronald Gamble, Charles F. Gammie, Zachary Gelles, Jose L. Gomez, Samuel E. Gralla, Paul Grimes, Leonid I. Gurvits, Shahar Hadar, Kari Haworth, Kazuhiro Hada , et al. (43 additional authors not shown)

Abstract: We present the Black Hole Explorer (BHEX), a mission that will produce the sharpest images in the history of astronomy by extending submillimeter Very-Long-Baseline Interferometry (VLBI) to space. BHEX will discover and measure the bright and narrow "photon ring" that is predicted to exist in images of black holes, produced from light that has orbited the black hole before escaping. This discovery… ▽ More We present the Black Hole Explorer (BHEX), a mission that will produce the sharpest images in the history of astronomy by extending submillimeter Very-Long-Baseline Interferometry (VLBI) to space. BHEX will discover and measure the bright and narrow "photon ring" that is predicted to exist in images of black holes, produced from light that has orbited the black hole before escaping. This discovery will expose universal features of a black hole's spacetime that are distinct from the complex astrophysics of the emitting plasma, allowing the first direct measurements of a supermassive black hole's spin. In addition to studying the properties of the nearby supermassive black holes M87* and Sgr A*, BHEX will measure the properties of dozens of additional supermassive black holes, providing crucial insights into the processes that drive their creation and growth. BHEX will also connect these supermassive black holes to their relativistic jets, elucidating the power source for the brightest and most efficient engines in the universe. BHEX will address fundamental open questions in the physics and astrophysics of black holes that cannot be answered without submillimeter space VLBI. The mission is enabled by recent technological breakthroughs, including the development of ultra-high-speed downlink using laser communications, and it leverages billions of dollars of existing ground infrastructure. We present the motivation for BHEX, its science goals and associated requirements, and the pathway to launch within the next decade. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: Proceedings for SPIE Astronomical Telescopes and Instrumentation

Photon Ring Interferometric Signatures Beyond The Universal Regime

Authors: He Jia, Eliot Quataert, Alexandru Lupsasca, George N. Wong

Abstract: We calculate the interferometric signatures of black hole photon rings beyond the universal regime by perturbatively including the effects of finite ring width. Our approach first slices a thick ring into a series of thin rings, each of which falls within the universal regime. We thus calculate the visibility of the thick ring by aggregating the contributions from each thin ring, and then perturba… ▽ More We calculate the interferometric signatures of black hole photon rings beyond the universal regime by perturbatively including the effects of finite ring width. Our approach first slices a thick ring into a series of thin rings, each of which falls within the universal regime. We thus calculate the visibility of the thick ring by aggregating the contributions from each thin ring, and then perturbatively expand the result into polynomials of the baseline length $u$. We show that the visibility amplitude of a thick ring depends on its "center-of-light" diameter; it also includes additional higher-order corrections due to the width of the ring, with the leading correction terms proportional to $u^2$ for the envelope and $u^3$ for the phase. We apply our method to images ray traced from general-relativistic magnetohydrodynamic (GRMHD) simulations and demonstrate that incorporating the higher-order corrections is crucial for accurately modeling the visibility of the first photon ring around M87*. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 10+6 pages, 7+3 figures, to be submitted

Ordered magnetic fields around the 3C 84 central black hole

Authors: G. F. Paraschos, J. -Y. Kim, M. Wielgus, J. Röder, T. P. Krichbaum, E. Ros, I. Agudo, I. Myserlis, M. Moscibrodzka, E. Traianou, J. A. Zensus, L. Blackburn, C. -K. Chan, S. Issaoun, M. Janssen, M. D. Johnson, V. L. Fish, K. Akiyama, A. Alberdi, W. Alef, J. C. Algaba, R. Anantua, K. Asada, R. Azulay, U. Bach , et al. (258 additional authors not shown)

Abstract: 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures a… ▽ More 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures and understand the physical conditions in the compact region of 3C84. We used EHT 228GHz observations and, given the limited (u,v)-coverage, applied geometric model fitting to the data. We also employed quasi-simultaneously observed, multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure. We report the detection of a highly ordered, strong magnetic field around the central, SMBH of 3C84. The brightness temperature analysis suggests that the system is in equipartition. We determined a turnover frequency of $ν_m=(113\pm4)$GHz, a corresponding synchrotron self-absorbed magnetic field of $B_{SSA}=(2.9\pm1.6)$G, and an equipartition magnetic field of $B_{eq}=(5.2\pm0.6)$G. Three components are resolved with the highest fractional polarisation detected for this object ($m_\textrm{net}=(17.0\pm3.9)$%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017-2018. We report a steeply negative slope of the spectrum at 228GHz. We used these findings to test models of jet formation, propagation, and Faraday rotation in 3C84. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C84. However, systematic uncertainties due to the limited (u,v)-coverage, however, cannot be ignored. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 15 pages, 6 figures, published in A&A

Journal ref: Issue: A&A Volume 682, February 2024; Article number: L3; Number of pages: 15

Balanced Turbulence and the Helicity Barrier in Black Hole Accretion

Authors: George N. Wong, Lev Arzamasskiy

Abstract: Horizon-scale observations from the Event Horizon Telescope (EHT) have enabled precision study of supermassive black hole accretion. Contemporary accretion modeling often treats the inflowing plasma as a single, thermal fluid, but microphysical kinetic effects can lead to significant deviations from this idealized picture. We investigate how the helicity barrier influences EHT-accessible electroma… ▽ More Horizon-scale observations from the Event Horizon Telescope (EHT) have enabled precision study of supermassive black hole accretion. Contemporary accretion modeling often treats the inflowing plasma as a single, thermal fluid, but microphysical kinetic effects can lead to significant deviations from this idealized picture. We investigate how the helicity barrier influences EHT-accessible electromagnetic observables by employing a simple model for electron heating based on kinetic physics and the cascade of energy and helicity in unbalanced turbulence. Although the helicity barrier plays only a minor role in regions with high plasma-beta, like in SANE disks, it may substantially impact in regions with more ordered magnetic fields, such as the jet and its surrounding wind in SANE flows as well as throughout the entire domain in MAD flows. In SANE flows, emission shifts from the funnel wall towards the lower-magnetization disk region; in MAD flows the emission morphology remains largely unchanged. Including the helicity barrier leads to characteristically lower electron temperatures, and neglecting it can lead to underestimated accretion rates and inferred jet powers. The corresponding higher plasma densities result in increased depolarization and Faraday depths thereby decreasing the amplitude of the beta_2 coefficient while leaving its angle unchanged. Both the increased jet power and lower |beta_2| may help alleviate outstanding tensions between modeling and EHT observations. We also find that the estimated ring diameter may be underestimated when the helicity barrier is neglected. Our results underscore the significance of the helicity barrier in shaping black hole observables and inferred accretion system parameters. △ Less

Submitted 26 December, 2023; originally announced December 2023.

Comments: 19 pages, 10 figures, accepted for publication in ApJ

The $230$ GHz Variability of Numerical Models of Sagittarius~A* I. Parameter Surveys on Varying Ion-electron Temperature Ratios Under Strongly Magnetized Conditions

Authors: Ho-Sang Chan, Chi-kwan Chan, Ben S. Prather, George N. Wong, Charles Gammie

Abstract: The $230$ GHz lightcurves of Sagittarius~A* (Sgr~A*) predicted by general relativistic magnetohydrodynamics (GRMHD) and ray-tracing (GRRT) models in Event Horizon Telescope Collaboration et al. (2022) have higher variability $M_{ΔT}$ compared to observations. In this series of papers, we explore the origin of such large brightness variability. In this first paper, we performed large GRRT parameter… ▽ More The $230$ GHz lightcurves of Sagittarius~A* (Sgr~A*) predicted by general relativistic magnetohydrodynamics (GRMHD) and ray-tracing (GRRT) models in Event Horizon Telescope Collaboration et al. (2022) have higher variability $M_{ΔT}$ compared to observations. In this series of papers, we explore the origin of such large brightness variability. In this first paper, we performed large GRRT parameter surveys that span from the optically thin to the optically thick regimes, covering the ion-electron temperature ratio under strongly magnetized conditions, $R_{\rm Low}$, from $1$ to $60$. We find that increasing $R_{\rm Low}$ can lead to either an increase or a reduction in $M_{ΔT}$ depending on other model parameters, making it consistent with the observed variability of Sgr~A* in some cases. Our analysis of GRRT image snapshots finds that the major contribution to the large $M_{ΔT}$ for the $R_{\rm Low} = 1$ models comes from the photon rings. However, secondary contributions from the accretion flow are also visible depending on the spin parameter. Our work demonstrates the importance of the electron temperature used for modelling radiatively inefficient accretion flows and places new constraints on the ion-electron temperature ratio. A more in-depth analysis for understanding the dependencies of $M_{ΔT}$ on $R_{\rm Low}$ will be performed in subsequent papers. △ Less

Submitted 4 February, 2024; v1 submitted 7 November, 2023; originally announced November 2023.

Comments: 15 Pages, 9 Figures

On the Comparison of AGN with GRMHD Simulations: II. M87

Authors: Richard Anantua, Angelo Ricarte, George Wong, Razieh Emami, Roger Blandford, Lani Oramas, Hayley West, Joaquin Duran, Brandon Curd

Abstract: Horizon-scale observations of the jetted active galactic nucleus M87 are compared with simulations spanning a broad range of dissipation mechanisms and plasma content in three-dimensional general relativistic flows around spinning black holes. Observations of synchrotron radiation from radio to X-ray frequencies can be compared with simulations by adding prescriptions specifying the relativistic e… ▽ More Horizon-scale observations of the jetted active galactic nucleus M87 are compared with simulations spanning a broad range of dissipation mechanisms and plasma content in three-dimensional general relativistic flows around spinning black holes. Observations of synchrotron radiation from radio to X-ray frequencies can be compared with simulations by adding prescriptions specifying the relativistic electron-plus-positron distribution function and associated radiative transfer coefficients. A suite of time-varying simulations with various spins, plasma magnetizations and turbulent heating and equipartition-based emission prescriptions (and piecewise combinations thereof) is chosen to represent distinct possibilities for the M87 jet/accretion flow/black hole (JAB) system. Simulation jet morphology, polarization and variation are then "observed" and compared with real observations to infer the rules that govern the polarized emissivity. Our models support several possible spin/emission model/plasma composition combinations supplying the jet in M87, whose black hole shadow has been observed down to the photon ring at 230 GHz by the Event Horizon Telescope (EHT). Net linear polarization and circular polarization constraints favor magnetically arrested disk (MAD) models whereas resolved linear polarization favors standard and normal evolution (SANE) in our parameter space. We also show that some MAD cases dominated by intrinsic circular polarization have near-linear V/I dependence on unpaired electron or positron content while SANE polarization exhibits markedly greater positron-dependent Faraday effects - future probes of the SANE/MAD dichotomy and plasma content with the EHT. This is the second work in a series also applying the "observing" simulations methodology to near-horizon regions of supermassive black holes in Sgr A* and 3C 279. △ Less

Submitted 22 December, 2023; v1 submitted 11 September, 2023; originally announced September 2023.

Comments: 27 pages, 27 figures

A search for pulsars around Sgr A* in the first Event Horizon Telescope dataset

Authors: Pablo Torne, Kuo Liu, Ralph P. Eatough, Jompoj Wongphechauxsorn, James M. Cordes, Gregory Desvignes, Mariafelicia De Laurentis, Michael Kramer, Scott M. Ransom, Shami Chatterjee, Robert Wharton, Ramesh Karuppusamy, Lindy Blackburn, Michael Janssen, Chi-kwan Chan, Geoffrey B. Crew, Lynn D. Matthews, Ciriaco Goddi, Helge Rottmann, Jan Wagner, Salvador Sanchez, Ignacio Ruiz, Federico Abbate, Geoffrey C. Bower, Juan J. Salamanca , et al. (261 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission… ▽ More The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission spectra - are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic Center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most-sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the Fast-Folding-Algorithm and single pulse search targeting both pulsars and burst-like transient emission; using the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction ($\lesssim$2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

Black Hole Polarimetry I: A Signature of Electromagnetic Energy Extraction

Authors: Andrew Chael, Alexandru Lupsasca, George N. Wong, Eliot Quataert

Abstract: 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… ▽ More 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 -- the phase $\angleβ_2$ of the second azimuthal Fourier mode of the linear polarization 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^φB^r$. The phase $\angleβ_2$ likewise depends on the ratio $B^φ/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β_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*. △ Less

Submitted 14 November, 2023; v1 submitted 12 July, 2023; originally announced July 2023.

Comments: 35 pages, 5 figures. Published in ApJ

Demonstrating Photon Ring Existence with Single-Baseline Polarimetry

Authors: Daniel C. M. Palumbo, George N. Wong, Andrew A. Chael, Michael D. Johnson

Abstract: Images of supermassive black hole accretion flows contain features of both curved spacetime and plasma structure. Inferring properties of the spacetime from images requires modeling the plasma properties, and vice versa. The Event Horizon Telescope Collaboration has imaged near-horizon millimeter emission from both Messier 87* (M87*) and Sagittarius A* (Sgr A*) with very-long-baseline interferomet… ▽ More Images of supermassive black hole accretion flows contain features of both curved spacetime and plasma structure. Inferring properties of the spacetime from images requires modeling the plasma properties, and vice versa. The Event Horizon Telescope Collaboration has imaged near-horizon millimeter emission from both Messier 87* (M87*) and Sagittarius A* (Sgr A*) with very-long-baseline interferometry (VLBI) and has found a preference for magnetically arrested disk (MAD) accretion in each case. MAD accretion enables spacetime measurements through future observations of the photon ring, the image feature composed of near-orbiting photons. The ordered fields and relatively weak Faraday rotation of MADs yield rotationally symmetric polarization when viewed at modest inclination. In this letter, we utilize this symmetry along with parallel transport symmetries to construct a gain-robust interferometric quantity that detects the transition between the weakly lensed accretion flow image and the strongly lensed photon ring. We predict a shift in polarimetric phases on long baselines and demonstrate that the photon rings in M87* and Sgr A* can be unambiguously detected {with sensitive, long-baseline measurements. For M87* we find that photon ring detection in snapshot observations requires $\sim1$ mJy sensitivity on $>15$ G$λ$ baselines at 230 GHz and above, which could be achieved with space-VLBI or higher-frequency ground-based VLBI. For Sgr A*, we find that interstellar scattering inhibits photon ring detectability at 230 GHz, but $\sim10$ mJy sensitivity on $>12$ G$λ$ baselines at 345 GHz is sufficient, which is accessible from the ground. For both sources, these sensitivity requirements may be relaxed by repeated observations and averaging. △ Less

Submitted 11 July, 2023; originally announced July 2023.

Comments: 14 pages, 7 figures, Accepted to ApJL

Mahakala: a Python-based Modular Ray-tracing and Radiative Transfer Algorithm for Curved Space-times

Authors: Aniket Sharma, Lia Medeiros, Chi-kwan Chan, Goni Halevi, Patrick D. Mullen, James M. Stone, George N. Wong

Abstract: We introduce Mahakala, a Python-based, modular, radiative ray-tracing code for curved space-times. We employ Google's JAX framework for accelerated automatic differentiation, which can efficiently compute Christoffel symbols directly from the metric, allowing the user to easily and quickly simulate photon trajectories through non-Kerr metrics. JAX also enables Mahakala to run in parallel on both C… ▽ More We introduce Mahakala, a Python-based, modular, radiative ray-tracing code for curved space-times. We employ Google's JAX framework for accelerated automatic differentiation, which can efficiently compute Christoffel symbols directly from the metric, allowing the user to easily and quickly simulate photon trajectories through non-Kerr metrics. JAX also enables Mahakala to run in parallel on both CPUs and GPUs and achieve speeds comparable to C-based codes. Mahakala natively uses the Cartesian Kerr-Schild coordinate system, which avoids numerical issues caused by the "pole" of spherical coordinates. We demonstrate Mahakala's capabilities by simulating the 1.3 mm wavelength images (the wavelength of Event Horizon Telescope observations) of general relativistic magnetohydrodynamic simulations of low-accretion rate supermassive black holes. The modular nature of Mahakala allows us to easily quantify the relative contribution of different regions of the flow to image features. We show that most of the emission seen in 1.3 mm images originates close to the black hole. We also quantify the relative contribution of the disk, forward jet, and counter jet to 1.3 mm images. △ Less

Submitted 7 April, 2023; originally announced April 2023.

Comments: 15 pages, 11 figures

$κ$monty: a Monte Carlo Compton Scattering code including non-thermal electrons

Authors: Jordy Davelaar, Benjamin R. Ryan, George N. Wong, Thomas Bronzwaer, Hector Olivares, Monika Mościbrodzka, Charles F. Gammie, Heino Falcke

Abstract: Low-luminosity active galactic nuclei are strong sources of X-ray emission produced by Compton scattering originating from the accretion flows surrounding their supermassive black holes. The shape and energy of the resulting spectrum depend on the shape of the underlying electron distribution function (DF). In this work, we present an extended version of the grmonty code, called $κ$monty. The grmo… ▽ More Low-luminosity active galactic nuclei are strong sources of X-ray emission produced by Compton scattering originating from the accretion flows surrounding their supermassive black holes. The shape and energy of the resulting spectrum depend on the shape of the underlying electron distribution function (DF). In this work, we present an extended version of the grmonty code, called $κ$monty. The grmonty code previously only included a thermal Maxwell Jütner electron distribution function. We extend the gromty code with non-thermal electron DFs, namely the $κ$ and power-law DFs, implement Cartesian Kerr-Schild coordinates, accelerate the code with MPI, and couple the code to the non-uniform AMR grid data from the GRMHD code BHAC. For the Compton scattering process, we derive two sampling kernels for both distribution functions. Finally, we present a series of code tests to verify the accuracy of our schemes. The implementation of non-thermal DFs opens the possibility of studying the effect of non-thermal emission on previously developed black hole accretion models. △ Less

Submitted 2 October, 2023; v1 submitted 27 March, 2023; originally announced March 2023.

Comments: 12 pages, 9 figures, accepted for publication in MNRAS

Comparison of Polarized Radiative Transfer Codes used by the EHT Collaboration

Authors: Ben S. Prather, Jason Dexter, Monika Moscibrodzka, Hung-Yi Pu, Thomas Bronzwaer, Jordy Davelaar, Ziri Younsi, Charles F. Gammie, Roman Gold, George N. Wong, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, Anne-Kathrin Baczko, David Ball, Mislav Baloković, John Barrett, Michi Bauböck, Bradford A. Benson, Dan Bintley , et al. (248 additional authors not shown)

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

Submitted 21 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ

Using Machine Learning to Link Black Hole Accretion Flows with Spatially Resolved Polarimetric Observables

Authors: Richard Qiu, Angelo Ricarte, Ramesh Narayan, George N. Wong, Andrew Chael, Daniel Palumbo

Abstract: 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, an… ▽ More 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. △ Less

Submitted 9 February, 2023; v1 submitted 9 December, 2022; originally announced December 2022.

Comments: 24 pages, 27 figures

Probing plasma composition with the next generation Event Horizon Telescope (ngEHT)

Authors: Razieh Emami, Richard Anantua, Angelo Ricarte, Sheperd S. Doeleman, Avery Broderick, George Wong, Lindy Blackburn, Maciek Wielgus, Ramesh Narayan, Grant Tremblay, Charles Alcock, Lars Hernquist, Randall Smith, Matthew Liska, Priyamvada Natarajan, Mark Vogelsberger, Brandon Curd, Joana A. Kramer

Abstract: We explore the plasma matter content in the innermost accretion disk/jet in M87* as relevant for an enthusiastic search for the signatures of anti-matter in the next generation of the Event Horizon Telescope (ngEHT). We model the impact of non-zero positron-to-electron ratio using different emission models including a constant electron to magnetic pressure (constant $β_e$ model) with a population… ▽ More We explore the plasma matter content in the innermost accretion disk/jet in M87* as relevant for an enthusiastic search for the signatures of anti-matter in the next generation of the Event Horizon Telescope (ngEHT). We model the impact of non-zero positron-to-electron ratio using different emission models including a constant electron to magnetic pressure (constant $β_e$ model) with a population of non-thermal electrons as well as a R-beta model populated with thermal electrons. In the former case, we pick a semi-analytic fit to the force-free region of a general relativistic magnetohydrodynamic (GRMHD) simulation, while in the latter case, we analyze the GRMHD simulations directly. In both cases, positrons are being added at the post-processing level. We generate polarized images and spectra for some of these models and find out that at the radio frequencies, both of the linear and the circular polarizations get enhanced per adding pairs. On the contrary, we show that at higher frequencies a substantial positron fraction washes out the circular polarization. We report strong degeneracies between different emission models and the positron fraction, though our non-thermal models show more sensitivities to the pair fraction than the thermal models. We conclude that a large theoretical image library is indeed required to fully understand the trends probed in this study, and to place them in the context of large set of parameters which also affect polarimetric images, such as magnetic field strength, black hole spin, and detailed aspects of the electron temperature and the distribution function. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: 14 pages, 9 figures

Emission Modeling in the EHT-ngEHT Age

Authors: Richard Anantua, Joaquín Dúran, Nathan Ngata, Lani Oramas, Razieh Emami, Angelo Ricarte, Brandon Curd, Jan Röder, Avery Broderick, Jeremy Wayland, George N. Wong, Sean Ressler

Abstract: This work proposes a methodology to test phenomenologically-motivated emission processes that account for the flux and polarization distribution and global structure of the 230 GHz sources imaged by the Event Horizon Telescope (EHT): Messier (M)87* and Sagittarius (Sgr) A*. We introduce to general relativistic magnetohydrodynamic (GRMHD) simulations some novel models to bridge the largely uncertai… ▽ More This work proposes a methodology to test phenomenologically-motivated emission processes that account for the flux and polarization distribution and global structure of the 230 GHz sources imaged by the Event Horizon Telescope (EHT): Messier (M)87* and Sagittarius (Sgr) A*. We introduce to general relativistic magnetohydrodynamic (GRMHD) simulations some novel models to bridge the largely uncertain mechanisms by which high-energy particles in jet/accretion flow/black hole (JAB) system plasmas attain billion degree temperatures and emit synchrotron radiation. The "Observing" JAB Systems methodology then partitions the simulation to apply different parametric models to regions governed by different plasma physics -- an advance over methods where one parametrization is used over simulation regions spanning thousands of gravitational radii from the central supermassive black hole. We present several classes of viewing-angle dependent morphologies, and highlight signatures of piecewise modeling and positron effects -- including a MAD/SANE dichotomy in which polarized maps appear dominated by intrinsic polarization in the MAD case and by Faraday effects in the SANE case. The library of images thus produced spans a wide range of morphologies awaiting discovery by the groundbreaking EHT instrument and its yet more sensitive, higher resolution next-generation counterpart ngEHT. △ Less

Submitted 11 November, 2022; originally announced November 2022.

Comments: 21 pages, 9 figures

Unraveling Twisty Linear Polarization Morphologies in Black Hole Images

Authors: Razieh Emami, Angelo Ricarte, George N. Wong, Daniel Palumbo, Dominic Chang, Sheperd S. Doeleman, Avery Broaderick, Ramesh Narayan, Maciek Wielgus, Lindy Blackburn, Ben S. Prather, Andrew A. Chael, Richard Anantua, Koushik Chatterjee, Ivan Marti-Vidal, Jose L. Gomez, Kazunori Akiyama, Matthew Liska, Lars Hernquist, Grant Tremblay, Mark Vogelsberger, Charles Alcock, Randall Smith, James Steiner, Paul Tiede , et al. (1 additional authors not shown)

Abstract: 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 magnet… ▽ More 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. △ Less

Submitted 28 March, 2023; v1 submitted 3 October, 2022; originally announced October 2022.

Comments: 26 pages, 20 figures, accepted for publication in ApJ

Effects of Hydrogen vs. Helium on Electromagnetic Black Hole Observables

Authors: George N. Wong, Charles F. Gammie

Abstract: The centers of our galaxy and the nearby Messier 87 are known to contain supermassive black holes, which support accretion flows that radiate across the electromagnetic spectrum. Although the composition of the accreting gas is unknown, it is likely a mix of ionized hydrogen and helium. We use a simple analytic model and a suite of numerical general relativistic magnetohydrodynamic accretion simul… ▽ More The centers of our galaxy and the nearby Messier 87 are known to contain supermassive black holes, which support accretion flows that radiate across the electromagnetic spectrum. Although the composition of the accreting gas is unknown, it is likely a mix of ionized hydrogen and helium. We use a simple analytic model and a suite of numerical general relativistic magnetohydrodynamic accretion simulations to study how polarimetric images and spectral energy distributions of the source are influenced by the hydrogen/helium content of the accreting matter. We aim to identify general trends rather than make quantitatively precise predictions, since it is not possible to fully explore the parameter space of accretion models. If the ion-to-electron temperature ratio is fixed, then increasing the helium fraction increases the gas temperature; to match the observational flux density constraints, the number density of electrons and magnetic field strengths must therefore decrease. In our numerical simulations, emission shifts from regions of low to high plasma beta -- both altering the morphology of the image and decreasing the variability of the light curve -- especially in strongly magnetized models with emission close to the midplane. In polarized images, we find that the model gas composition influences the degree to which linear polarization is (de)scrambled and therefore affects estimates for the resolved linear polarization fraction. We also find that the spectra of helium-composition flows peak at higher frequencies and exhibit higher luminosities. We conclude that gas composition may play an important role in predictive models for black hole accretion. △ Less

Submitted 27 July, 2022; originally announced July 2022.

Comments: 17 pages, 8 figures, accepted for publication in ApJ

A deep spectromorphological study of the $γ$-ray emission surrounding the young massive stellar cluster Westerlund 1

Authors: F. Aharonian, H. Ashkar, M. Backes, V. Barbosa Martins, Y. Becherini, D. Berge, B. Bi, M. Böttcher, M. de Bony de Lavergne, F. Bradascio, R. Brose, F. Brun, T. Bulik, C. Burger-Scheidlin, F. Cangemi, S. Caroff, S. Casanova, M. Cerruti, T. Chand, S. Chandra, A. Chen, O. Chibueze, P. Cristofari, J. Damascene Mbarubucyeye, A. Djannati-Ataï , et al. (134 additional authors not shown)

Abstract: Young massive stellar clusters are extreme environments and potentially provide the means for efficient particle acceleration. Indeed, they are increasingly considered as being responsible for a significant fraction of cosmic rays (CRs) accelerated within the Milky Way. Westerlund 1, the most massive known young stellar cluster in our Galaxy is a prime candidate for studying this hypothesis. While… ▽ More Young massive stellar clusters are extreme environments and potentially provide the means for efficient particle acceleration. Indeed, they are increasingly considered as being responsible for a significant fraction of cosmic rays (CRs) accelerated within the Milky Way. Westerlund 1, the most massive known young stellar cluster in our Galaxy is a prime candidate for studying this hypothesis. While the very-high-energy $γ$-ray source HESS J1646-458 has been detected in the vicinity of Westerlund 1 in the past, its association could not be firmly identified. We aim to identify the physical processes responsible for the $γ$-ray emission around Westerlund 1 and thus to better understand the role of massive stellar clusters in the acceleration of Galactic CRs. Using 164 hours of data recorded with the High Energy Stereoscopic System (H.E.S.S.), we carried out a deep spectromorphological study of the $γ$-ray emission of HESS J1646-458. We furthermore employed H I and CO observations of the region to infer the presence of gas that could serve as target material for interactions of accelerated CRs. We detected large-scale ($\sim 2^\circ$ diameter) $γ$-ray emission with a complex morphology, exhibiting a shell-like structure and showing no significant variation with $γ$-ray energy. The combined energy spectrum of the emission extends to several tens of TeV, and is uniform across the entire source region. We did not find a clear correlation of the $γ$-ray emission with gas clouds as identified through H I and CO observations. We conclude that, of the known objects within the region, only Westerlund 1 can explain the bulk of the $γ$-ray emission. Several CR acceleration sites and mechanisms are conceivable, and discussed in detail. (abridged) △ Less

Submitted 10 November, 2022; v1 submitted 22 July, 2022; originally announced July 2022.

Comments: 15 pages, 9 figures. Corresponding authors: L. Mohrmann, S. Ohm, R. Rauth, A. Specovius. v2: corrected affiliation of M. Vecchi (still incorrect in journal version)

Journal ref: A&A 666, A124 (2022)

Photon Ring Symmetries in Simulated Linear Polarization Images of Messier 87*

Authors: Daniel C. M. Palumbo, George N. Wong

Abstract: The Event Horizon Telescope (EHT) recently released the first linearly polarized images of the accretion flow around the supermassive black hole Messier 87*, hereafter \m{}. The spiraling polarization pattern found in EHT images favored magnetically arrested disks (MADs) as the explanation for the EHT image. With next-generation improvements to very long baseline interferometry (VLBI) on the horiz… ▽ More The Event Horizon Telescope (EHT) recently released the first linearly polarized images of the accretion flow around the supermassive black hole Messier 87*, hereafter \m{}. The spiraling polarization pattern found in EHT images favored magnetically arrested disks (MADs) as the explanation for the EHT image. With next-generation improvements to very long baseline interferometry (VLBI) on the horizon, understanding similar polarized features in the highly lensed structure known as the "photon ring," where photons make multiple half-orbits about the black hole before reaching the observer, will be critical to analysis of future images. Recent work has indicated that this image region may be depolarized relative to more direct emission. We expand this observation by decomposing photon half-orbits in the EHT library of simulated images of the \m{} accretion system and find that images of MAD simulations show a relative depolarization of the photon ring attributable to destructive interference of oppositely spiraling electric field vectors; this antisymmetry, which arises purely from strong gravitational lensing, can produce up to ${\sim}50\%$ depolarization in the photon ring region with respect to the direct image. In systems that are not magnetically arrested and with the exception of systems with high spin and ions and electrons of equal temperature, we find that highly lensed indirect sub-images are almost completely depolarized, causing a modest depolarization of the photon ring region in the complete image. We predict that next-generation EHT observations of \m{} polarization should jointly constrain the black hole spin and the underlying emission and magnetic field geometry. △ Less

Submitted 1 March, 2022; originally announced March 2022.

Comments: Accepted to ApJ. 11 pages, 4 figures

PATOKA: Simulating Electromagnetic Observables of Black Hole Accretion

Authors: George N. Wong, Ben S. Prather, Vedant Dhruv, Benjamin R. Ryan, Monika Moscibrodzka, Chi-kwan Chan, Abhishek V. Joshi, Ricardo Yarza, Angelo Ricarte, Hotaka Shiokawa, Joshua C. Dolence, Scott C. Noble, Jonathan C. McKinney, Charles F. Gammie

Abstract: The Event Horizon Telescope (EHT) has released analyses of reconstructed images of horizon-scale millimeter emission near the supermassive black hole at the center of the M87 galaxy. Parts of the analyses made use of a large library of synthetic black hole images and spectra, which were produced using numerical general relativistic magnetohydrodynamics fluid simulations and polarized ray tracing.… ▽ More The Event Horizon Telescope (EHT) has released analyses of reconstructed images of horizon-scale millimeter emission near the supermassive black hole at the center of the M87 galaxy. Parts of the analyses made use of a large library of synthetic black hole images and spectra, which were produced using numerical general relativistic magnetohydrodynamics fluid simulations and polarized ray tracing. In this article, we describe the PATOKA pipeline, which was used to generate the Illinois contribution to the EHT simulation library. We begin by describing the relevant accretion systems and radiative processes. We then describe the details of the three numerical codes we use, iharm, ipole, and igrmonty, paying particular attention to differences between the current generation of the codes and the originally published versions. Finally, we provide a brief overview of simulated data as produced by PATOKA and conclude with a discussion of limitations and future directions. △ Less

Submitted 11 March, 2022; v1 submitted 23 February, 2022; originally announced February 2022.

Comments: 42 pages, 17 figures, accepted for publication in ApJS

Modelling the Gamma-Ray Morphology of HESS J1804-216 from Two Supernova Remnants in a Hadronic Scenario

Authors: K. Feijen, S. Einecke, G. Rowell, C. Braiding, M. G. Burton, G. F. Wong

Abstract: HESS J1804-216 is one of the brightest yet most mysterious TeV gamma-ray sources discovered to date. Previous arc-minute scale studies of the interstellar medium (ISM) surrounding this TeV gamma-ray source revealed HESS J1804-216 is likely powered by a mature supernova remnant (SNR) or pulsar, hence its origin remains uncertain. In this paper, we focus on the diffusive escape of cosmic-ray protons… ▽ More HESS J1804-216 is one of the brightest yet most mysterious TeV gamma-ray sources discovered to date. Previous arc-minute scale studies of the interstellar medium (ISM) surrounding this TeV gamma-ray source revealed HESS J1804-216 is likely powered by a mature supernova remnant (SNR) or pulsar, hence its origin remains uncertain. In this paper, we focus on the diffusive escape of cosmic-ray protons from potential SNR accelerators. These cosmic rays interact with the ISM to produce TeV gamma-rays. We utilise the isotropic diffusion equation solution for particles escaping from a shell, to model the energy-dependent escape and propagation of protons into the ISM. This work is the first attempt at modelling the spatial morphology of gamma-rays towards HESS J1804-216, using arc-minute ISM observations from both Mopra and the Southern Galactic Plane Survey. The spectral and spatial distributions of gamma-rays for the two nearby potential SNR counterparts, SNR G8.7-0.1 and the progenitor SNR of PSR J1803-2137, are presented here. We vary the diffusion parameters and particle spectrum and use a grid search approach to find the best combination of model parameters. We conclude that moderately slow diffusion is required for both candidates. The most promising candidate to be powering the TeV gamma-rays from HESS J1804-216 in a hadronic scenario is the progenitor SNR of PSR J1803-213. △ Less

Submitted 27 January, 2022; originally announced January 2022.

Comments: 13 pages, 13 figures, 4 tables

Event Horizon Telescope observations of the jet launching and collimation in Centaurus A

Authors: Michael Janssen, Heino Falcke, Matthias Kadler, Eduardo Ros, Maciek Wielgus, Kazunori Akiyama, Mislav Baloković, Lindy Blackburn, Katherine L. Bouman, Andrew Chael, Chi-kwan Chan, Koushik Chatterjee, Jordy Davelaar, Philip G. Edwards, Christian M. Fromm, José L. Gómez, Ciriaco Goddi, Sara Issaoun, Michael D. Johnson, Junhan Kim, Jun Yi Koay, Thomas P. Krichbaum, Jun Liu, Elisabetta Liuzzo, Sera Markoff , et al. (215 additional authors not shown)

Abstract: Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supe… ▽ More Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supermassive black holes (SMBHs) in Messier 87 and our galactic center. A large southern declination of $-43^{\circ}$ has however prevented VLBI imaging of Centaurus A below $λ1$cm thus far. Here, we show the millimeter VLBI image of the source, which we obtained with the Event Horizon Telescope at $228$GHz. Compared to previous observations, we image Centaurus A's jet at a tenfold higher frequency and sixteen times sharper resolution and thereby probe sub-lightday structures. We reveal a highly-collimated, asymmetrically edge-brightened jet as well as the fainter counterjet. We find that Centaurus A's source structure resembles the jet in Messier 87 on ${\sim}500r_g$ scales remarkably well. Furthermore, we identify the location of Centaurus A's SMBH with respect to its resolved jet core at $λ1.3$mm and conclude that the source's event horizon shadow should be visible at THz frequencies. This location further supports the universal scale invariance of black holes over a wide range of masses. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 27 pages, 9 figures. This is a post-peer-review, pre-copyedit version of an article published in Nature Astronomy. The final authenticated version is available online at: https://doi.org/10.1038/s41550-021-01417-w

Journal ref: Nature Astronomy, July 2021, Volume 5, p. 1017-1028

The Variability of the Black-Hole Image in M87 at the Dynamical Time Scale

Authors: Kaushik Satapathy, Dimitrios Psaltis, Feryal Ozel, Lia Medeiros, Sean T. Dougall, Chi-kwan Chan, Maciek Wielgus, Ben S. Prather, George N. Wong, Charles F. Gammie, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David R. Ball, Mislav Baloković, John Barrett, Bradford A. Benson, Dan Bintley, Lindy Blackburn, Raymond Blundell , et al. (213 additional authors not shown)

Abstract: The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expect… ▽ More The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure phase measurements on all six linearly independent non-trivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of $\sim3-5^\circ$. The only triangles that exhibit substantially higher variability ($\sim90-180^\circ$) are the ones with baselines that cross visibility amplitude minima on the $u-v$ plane, as expected from theoretical modeling. We used two sets of General Relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black-hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black-hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: Accepted for Publication in ApJ

Iharm3D: Vectorized General Relativistic Magnetohydrodynamics

Authors: Ben S. Prather, George N. Wong, Vedant Dhruv, Benjamin R. Ryan, Joshua C. Dolence, Sean M. Ressler, Charles F. Gammie

Abstract: Iharm3D is an open-source C code for simulating black hole accretion systems in arbitrary stationary spacetimes using ideal general-relativistic magnetohydrodynamics (GRMHD). It is an implementation of the HARM ("High Accuracy Relativistic Magnetohydrodynamics") algorithm outlined in Gammie et al. (2003) with updates as outlined in McKinney & Gammie (2004) and Noble et al. (2006). The code is most… ▽ More Iharm3D is an open-source C code for simulating black hole accretion systems in arbitrary stationary spacetimes using ideal general-relativistic magnetohydrodynamics (GRMHD). It is an implementation of the HARM ("High Accuracy Relativistic Magnetohydrodynamics") algorithm outlined in Gammie et al. (2003) with updates as outlined in McKinney & Gammie (2004) and Noble et al. (2006). The code is most directly derived from Ryan et al. (2015) but with radiative transfer portions removed. HARM is a conservative finite-volume scheme for solving the equations of ideal GRMHD, a hyperbolic system of partial differential equations, on a logically Cartesian mesh in arbitrary coordinates. △ Less

Submitted 19 October, 2021; originally announced October 2021.

Comments: 4 pages, 2 figures. Associated code at https://github.com/AFD-Illinois/iharm3d

Journal ref: Journal of Open Source Software, 6(66), 3336 (2021)

VLBInet: Radio Interferometry Data Classification for EHT with Neural Networks

Authors: Joshua Yao-Yu Lin, Dominic W. Pesce, George N. Wong, Ajay Uppili Arasanipalai, Ben S. Prather, Charles F. Gammie

Abstract: The Event Horizon Telescope (EHT) recently released the first horizon-scale images of the black hole in M87. Combined with other astronomical data, these images constrain the mass and spin of the hole as well as the accretion rate and magnetic flux trapped on the hole. An important question for the EHT is how well key parameters, such as trapped magnetic flux and the associated disk models, can be… ▽ More The Event Horizon Telescope (EHT) recently released the first horizon-scale images of the black hole in M87. Combined with other astronomical data, these images constrain the mass and spin of the hole as well as the accretion rate and magnetic flux trapped on the hole. An important question for the EHT is how well key parameters, such as trapped magnetic flux and the associated disk models, can be extracted from present and future EHT VLBI data products. The process of modeling visibilities and analyzing them is complicated by the fact that the data are sparsely sampled in the Fourier domain while most of the theory/simulation is constructed in the image domain. Here we propose a data-driven approach to analyze complex visibilities and closure quantities for radio interferometric data with neural networks. Using mock interferometric data, we show that our neural networks are able to infer the accretion state as either high magnetic flux (MAD) or low magnetic flux (SANE), suggesting that it is possible to perform parameter extraction directly in the visibility domain without image reconstruction. We have applied VLBInet to real M87 EHT data taken on four different days in 2017 (April 5, 6, 10, 11), and our neural networks give a score prediction 0.52, 0.4, 0.43, 0.76 for each day, with an average score 0.53, which shows no significant indication for the data to lean toward either the MAD or SANE state. △ Less

Submitted 14 October, 2021; originally announced October 2021.

Comments: 10 pages, 7 figures

Radiation GRMHD simulations of M87: funnel properties and prospects for gap acceleration

Authors: Philippe Z. Yao, Jason Dexter, Alexander Y. Chen, Benjamin R. Ryan, George N. Wong

Abstract: We use the public code ebhlight to carry out 3D radiative general relativistic magnetohydrodynamics (GRMHD) simulations of accretion onto the supermassive black hole in M87. The simulations self-consistently evolve a frequency-dependent Monte Carlo description of the radiation field produced by the accretion flow. We explore two limits of accumulated magnetic flux at the black hole (SANE and MAD),… ▽ More We use the public code ebhlight to carry out 3D radiative general relativistic magnetohydrodynamics (GRMHD) simulations of accretion onto the supermassive black hole in M87. The simulations self-consistently evolve a frequency-dependent Monte Carlo description of the radiation field produced by the accretion flow. We explore two limits of accumulated magnetic flux at the black hole (SANE and MAD), each coupled to several sub-grid prescriptions for electron heating that are motivated by models of turbulence and magnetic reconnection. We present convergence studies for the radiation field and study its properties. We find that the near-horizon photon energy density is an order of magnitude higher than is predicted by simple isotropic estimates from the observed luminosity. The radially dependent photon momentum distribution is anisotropic and can be modeled by a set of point-sources near the equatorial plane. We draw properties of the radiation and magnetic field from the simulation and feed them into an analytic model of gap acceleration to estimate the very high energy (VHE) gamma-ray luminosity from the magnetized jet funnel, assuming that a gap is able to form. We find luminosities of $\rm \sim 10^{41} \, erg \, s^{-1}$ for MAD models and $\rm \sim 2\times 10^{40} \, erg \, s^{-1}$ for SANE models, which are comparable to measurements of M87's VHE flares. The time-dependence seen in our calculations is insufficient to explain the flaring behavior. Our results provide a step towards bridging theoretical models of near-horizon properties seen in black hole images with the VHE activity of M87. △ Less

Submitted 5 September, 2021; v1 submitted 26 August, 2021; originally announced August 2021.

Comments: 16 pages, 17 figures, Accepted for publication by MNRAS

The Polarized Image of a Synchrotron Emitting Ring of Gas Orbiting a Black Hole

Authors: Ramesh Narayan, Daniel C. M. Palumbo, Michael D. Johnson, Zachary Gelles, Elizabeth Himwich, Dominic O. Chang, Angelo Ricarte, Jason Dexter, Charles F. Gammie, Andrew A. Chael, The Event Horizon Telescope Collaboration, :, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David Ball, Mislav Balokovic, John Barrett, Bradford A. Benson, Dan Bintley , et al. (215 additional authors not shown)

Abstract: Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equ… ▽ More Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov (2002) and conservation of the Walker-Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images. △ Less

Submitted 13 May, 2021; v1 submitted 4 May, 2021; originally announced May 2021.

Comments: 29 pages, 11 figures, published in ApJ on May 3

Journal ref: ApJ 912 35 (2021)

The Jet-Disk Boundary Layer in Black Hole Accretion

Authors: George N. Wong, Yufeng Du, Ben S. Prather, Charles F. Gammie

Abstract: Magnetic fields lines are trapped in black hole event horizons by accreting plasma. If the trapped field lines are lightly loaded with plasma, then their motion is controlled by their footpoints on the horizon and thus by the spin of the black hole. In this paper, we investigate the boundary layer between lightly loaded polar field lines and a dense, equatorial accretion flow. We present an analyt… ▽ More Magnetic fields lines are trapped in black hole event horizons by accreting plasma. If the trapped field lines are lightly loaded with plasma, then their motion is controlled by their footpoints on the horizon and thus by the spin of the black hole. In this paper, we investigate the boundary layer between lightly loaded polar field lines and a dense, equatorial accretion flow. We present an analytic model for aligned prograde and retrograde accretion systems and argue that there is significant shear across this "jet-disk boundary" at most radii for all black hole spins. Specializing to retrograde aligned accretion, where the model predicts the strongest shear, we show numerically that the jet-disk boundary is unstable. The resulting mixing layer episodically loads plasma onto trapped field lines where it is heated, forced to rotate with the hole, and permitted to escape outward into the jet. In one case we follow the mass loading in detail using Lagrangian tracer particles and find a time-averaged mass-loading rate ~ 0.01 Mdot. △ Less

Submitted 14 April, 2021; originally announced April 2021.

Comments: 17 pages, 13 figures, accepted for publication in ApJ

The Role of Adaptive Ray Tracing in Analyzing Black Hole Structure

Authors: Z. Gelles, B. S. Prather, D. C. M. Palumbo, M. D. Johnson, G. N. Wong, B. Georgiev

Abstract: The recent advent of the Event Horizon Telescope (EHT) has made direct imaging of supermassive black holes a reality. Simulated images of black holes produced via general relativistic ray tracing and radiative transfer provide a key counterpart to these observational efforts. Black hole images have a wide range of physically interesting image structures, ranging from extremely fine scales in their… ▽ More The recent advent of the Event Horizon Telescope (EHT) has made direct imaging of supermassive black holes a reality. Simulated images of black holes produced via general relativistic ray tracing and radiative transfer provide a key counterpart to these observational efforts. Black hole images have a wide range of physically interesting image structures, ranging from extremely fine scales in their lensed "photon rings" to the very large scales in their relativistic jets. The multi-scale nature of the black hole system is therefore suitable for a multi-scale approach to generating simulated images that capture all key elements of the system. Here, we present a prescription for adaptive ray tracing, which enables efficient computation of extremely high resolution images of black holes. Using the polarized ray-tracing code ipole, we image a combination of semi-analytic and GRMHD models, and we show that images can be reproduced with mean squared error of less than 0.1% even after tracing 12x fewer rays. We then use adaptive ray tracing to explore properties of the photon ring. We illustrate the behavior of individual subrings in GRMHD simulations, and we explore their signatures in interferometric visibilities. △ Less

Submitted 12 May, 2021; v1 submitted 12 March, 2021; originally announced March 2021.

Comments: 17 pages, 12 figures. Published in ApJ on May 3, 2021

Journal ref: ApJ 912 39 (2021)

Self-interacting Dark Matter with Scalar Dilepton Mediator

Authors: Chung Kao, Yue-Lin Sming Tsai, Gwo-Guang Wong

Abstract: The cold dark matter (CDM) candidate with weakly interacting massive particles can successfully explain the observed dark matter relic density in cosmic scale and the large-scale structure of the Universe. However, a number of observations at the satellite galaxy scale seem to be inconsistent with CDM simulation. This is known as the small-scale problem of CDM. In recent years, it has been demonst… ▽ More The cold dark matter (CDM) candidate with weakly interacting massive particles can successfully explain the observed dark matter relic density in cosmic scale and the large-scale structure of the Universe. However, a number of observations at the satellite galaxy scale seem to be inconsistent with CDM simulation. This is known as the small-scale problem of CDM. In recent years, it has been demonstrated that self-interacting dark matter (SIDM) with a light mediator offers a reasonable explanation for the small-scale problem. We adopt a simple model with SIDM and focus on the effects of Sommerfeld enhancement. In this model, the dark matter candidate is a leptonic scalar particle with a light mediator. We have found several regions of the parameter space with proper masses and coupling strength generating a relic density that is consistent with the observed CDM relic density. Furthermore, this model satisfies the constraints of recent direct searches and indirect detection for dark matter as well as the effective number of neutrinos and the observed small-scale structure of the Universe. In addition, this model with the favored parameters can resolve the discrepancies between astrophysical observations and $N$-body simulations. △ Less

Submitted 3 March, 2021; v1 submitted 30 December, 2020; originally announced December 2020.

Comments: REVTEX, 32 pages, 11 figures. To be published in Phys. Rev. D

Report number: OU-HEP-201228

Journal ref: Phys. Rev. D 103, 055021 (2021)

Pair Drizzle around Sub-Eddington Supermassive Black Holes

Authors: George N. Wong, Benjamin R. Ryan, Charles F. Gammie

Abstract: Electron-positron pair creation near sub-Eddington accretion rate black holes is believed to be dominated by the Breit-Wheeler process (photon-photon collisions). The interacting high energy photons are produced when unscreened electric fields accelerate leptons either in coherent, macroscopic gaps or in incoherent structures embedded in the turbulent plasma flow. The latter type of acceleration r… ▽ More Electron-positron pair creation near sub-Eddington accretion rate black holes is believed to be dominated by the Breit-Wheeler process (photon-photon collisions). The interacting high energy photons are produced when unscreened electric fields accelerate leptons either in coherent, macroscopic gaps or in incoherent structures embedded in the turbulent plasma flow. The latter type of acceleration results in a drizzle of pair production sourced by photons from the background radiation field whose energies are near the pair-production threshold. In this work, we use radiation GRMHD simulations to extend an earlier study of pair drizzle by Moscibrodzka et al. We focus on low-magnetization (SANE) accretion onto supermassive Kerr black holes and consider radiation due to synchrotron, bremsstrahlung, and Compton upscattering processes. We confirm that pair drizzle in M87 is sufficient to keep the magnetospheric charge density orders of magnitude above the Goldreich-Julian density. We also find that pair production peaks along the jet-disk boundary. △ Less

Submitted 8 December, 2020; originally announced December 2020.

Comments: 14 pages, 8 figures, accepted for publication in ApJ

Arc-minute-scale studies of the interstellar gas towards HESS$\,$J1804$-$216: Still an unidentified TeV $γ$-ray source

Authors: K. Feijen, G. Rowell, S. Einecke, C. Braiding, M. G. Burton, N. Maxted, F. Voisin, G. F. Wong

Abstract: The Galactic TeV $γ$-ray source HESS$\,$J1804$-$216 is currently an unidentified source. In an attempt to unveil its origin, we present here the most detailed study of interstellar gas using data from the Mopra Southern Galactic Plane CO Survey, 7 and 12$\,$mm wavelength Mopra surveys and Southern Galactic Plane Survey of HI. Several components of atomic and molecular gas are found to overlap HESS… ▽ More The Galactic TeV $γ$-ray source HESS$\,$J1804$-$216 is currently an unidentified source. In an attempt to unveil its origin, we present here the most detailed study of interstellar gas using data from the Mopra Southern Galactic Plane CO Survey, 7 and 12$\,$mm wavelength Mopra surveys and Southern Galactic Plane Survey of HI. Several components of atomic and molecular gas are found to overlap HESS$\,$J1804$-$216 at various velocities along the line of sight. The CS(1-0) emission clumps confirm the presence of dense gas. Both correlation and anti-correlation between the gas and TeV $γ$-ray emission have been identified in various gas tracers, enabling several origin scenarios for the TeV $γ$-ray emission from HESS$\,$J1804$-$216. For a hadronic scenario, SNR$\,$G8.7$-$0.1 and the progenitor SNR of PSR$\,$J1803$-$2137 require cosmic ray (CR) enhancement factors of $\mathord{\sim} 50$ times the solar neighbour CR flux value to produce the TeV $γ$-ray emission. Assuming an isotropic diffusion model, CRs from both these SNRs require a slow diffusion coefficient, as found for other TeV SNRs associated with adjacent ISM gas. The morphology of gas located at 3.8$\,$kpc (the dispersion measure distance to PSR$\,$J1803$-$2137) tends to anti-correlate with features of the TeV emission from HESS$\,$J1804$-$216, making the leptonic scenario possible. Both pure hadronic and pure leptonic scenarios thus remain plausible. △ Less

Submitted 17 November, 2020; originally announced November 2020.

Comments: 29 pages, 23 figures, 5 tables, accepted for publication in PASA

Photon Ring Autocorrelations

Authors: Shahar Hadar, Michael D. Johnson, Alexandru Lupsasca, George N. Wong

Abstract: In the presence of a black hole, light sources connect to observers along multiple paths. As a result, observed brightness fluctuations must be correlated across different times and positions in black hole images. Photons that execute multiple orbits around the black hole appear near a critical curve in the observer sky, giving rise to the photon ring. In this paper, a novel observable is proposed… ▽ More In the presence of a black hole, light sources connect to observers along multiple paths. As a result, observed brightness fluctuations must be correlated across different times and positions in black hole images. Photons that execute multiple orbits around the black hole appear near a critical curve in the observer sky, giving rise to the photon ring. In this paper, a novel observable is proposed: the two-point correlation function of intensity fluctuations on the photon ring. This correlation function is analytically computed for a Kerr black hole surrounded by stochastic equatorial emission, with source statistics motivated by simulations of a turbulent accretion flow. It is shown that this two-point function exhibits a universal, self-similar structure consisting of multiple peaks of identical shape: while the profile of each peak encodes statistical properties of fluctuations in the source, the locations and heights of the peaks are determined purely by the black hole parameters. Measuring these peaks would demonstrate the existence of the photon ring without resolving its thickness, and would provide estimates of black hole mass and spin. With regular monitoring over sufficiently long timescales, this measurement could be possible via interferometric imaging with modest improvements to the Event Horizon Telescope. △ Less

Submitted 18 June, 2023; v1 submitted 7 October, 2020; originally announced October 2020.

Comments: 31 pages, 3 figures. v2: published version. v3: minor corrections, summarized in an appendix; conclusions essentially unchanged

Journal ref: Phys. Rev. D 103, 104038 (2021)

Gravitational Test Beyond the First Post-Newtonian Order with the Shadow of the M87 Black Hole

Authors: Dimitrios Psaltis, Lia Medeiros, Pierre Christian, Feryal Ozel, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Keiichi Asada, Rebecca Azulay, David Ball, Mislav Balokovic, John Barrett, Dan Bintley, Lindy Blackburn, Wilfred Boland, Geoffrey C. Bower, Michael Bremer, Christiaan D. Brinkerink, Roger Brissenden, Silke Britzen, Dominique Broguiere, Thomas Bronzwaer, Do-Young Byun, John E. Carlstrom, Andrew Chael , et al. (163 additional authors not shown)

Abstract: The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the p… ▽ More The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the predicted black-hole shadows that are inconsistent with even the current EHT measurements. We use numerical calculations of regular, parametric, non-Kerr metrics to identify the common characteristic among these different parametrizations that control the predicted shadow size. We show that the shadow-size measurements place significant constraints on deviation parameters that control the second post-Newtonian and higher orders of each metric and are, therefore, inaccessible to weak-field tests. The new constraints are complementary to those imposed by observations of gravitational waves from stellar-mass sources. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: Physical Review Letters

Monitoring the Morphology of M87* in 2009-2017 with the Event Horizon Telescope

Authors: Maciek Wielgus, Kazunori Akiyama, Lindy Blackburn, Chi-kwan Chan, Jason Dexter, Sheperd S. Doeleman, Vincent L. Fish, Sara Issaoun, Michael D. Johnson, Thomas P. Krichbaum, Ru-Sen Lu, Dominic W. Pesce, George N. Wong, Geoffrey C. Bower, Avery E. Broderick, Andrew Chael, Koushik Chatterjee, Charles F. Gammie, Boris Georgiev, Kazuhiro Hada, Laurent Loinard, Sera Markoff, Daniel P. Marrone, Richard Plambeck, Jonathan Weintroub , et al. (4 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) has recently delivered the first resolved images of M87*, the supermassive black hole in the center of the M87 galaxy. These images were produced using 230 GHz observations performed in 2017 April. Additional observations are required to investigate the persistence of the primary image feature - a ring with azimuthal brightness asymmetry - and to quantify the imag… ▽ More The Event Horizon Telescope (EHT) has recently delivered the first resolved images of M87*, the supermassive black hole in the center of the M87 galaxy. These images were produced using 230 GHz observations performed in 2017 April. Additional observations are required to investigate the persistence of the primary image feature - a ring with azimuthal brightness asymmetry - and to quantify the image variability on event horizon scales. To address this need, we analyze M87* data collected with prototype EHT arrays in 2009, 2011, 2012, and 2013. While these observations do not contain enough information to produce images, they are sufficient to constrain simple geometric models. We develop a modeling approach based on the framework utilized for the 2017 EHT data analysis and validate our procedures using synthetic data. Applying the same approach to the observational data sets, we find the M87* morphology in 2009-2017 to be consistent with a persistent asymmetric ring of ~40 uas diameter. The position angle of the peak intensity varies in time. In particular, we find a significant difference between the position angle measured in 2013 and 2017. These variations are in broad agreement with predictions of a subset of general relativistic magnetohydrodynamic simulations. We show that quantifying the variability across multiple observational epochs has the potential to constrain the physical properties of the source, such as the accretion state or the black hole spin. △ Less

Submitted 24 September, 2020; originally announced September 2020.

Journal ref: The Astrophysical Journal, 901:67 (2020)

Black Hole Glimmer Signatures of Mass, Spin, and Inclination

Authors: George N. Wong

Abstract: Gravitational lensing near a black hole is strong enough that light rays can circle the event horizon multiple times. Photons emitted in multiple directions at a single event, perhaps because of localized, impulsive heating of accreting plasma, take multiple paths to a distant observer. In the Kerr geometry, each path is associated with a distinct light travel time and a distinct arrival location… ▽ More Gravitational lensing near a black hole is strong enough that light rays can circle the event horizon multiple times. Photons emitted in multiple directions at a single event, perhaps because of localized, impulsive heating of accreting plasma, take multiple paths to a distant observer. In the Kerr geometry, each path is associated with a distinct light travel time and a distinct arrival location in the image plane, producing black hole glimmer. This sequence of arrival times and locations uniquely encodes the mass and spin of the black hole and can be understood in terms of properties of bound photon orbits. We provide a geometrically motivated treatment of Kerr glimmer and evaluate it numerically for simple hotspot models to show that glimmer can be measured in a finite-resolution observation. We discuss potential measurement methods and implications for tests of the Kerr hypothesis. △ Less

Submitted 28 January, 2021; v1 submitted 14 September, 2020; originally announced September 2020.

Comments: 12 pages, 7 figures; added more detailed discussion of inclination; accepted to ApJ

Decomposing the Internal Faraday Rotation of Black Hole Accretion Flows

Authors: Angelo Ricarte, Ben S. Prather, George N. Wong, Ramesh Narayan, Charles Gammie, Michael Johnson

Abstract: Faraday rotation has been seen at millimeter wavelengths in several low luminosity active galactic nuclei, including Event Horizon Telescope (EHT) targets M87* and Sgr A*. The observed rotation measure (RM) probes the density, magnetic field, and temperature of material integrated along the line of sight. To better understand how accretion disc conditions are reflected in the RM, we perform polari… ▽ More Faraday rotation has been seen at millimeter wavelengths in several low luminosity active galactic nuclei, including Event Horizon Telescope (EHT) targets M87* and Sgr A*. The observed rotation measure (RM) probes the density, magnetic field, and temperature of material integrated along the line of sight. To better understand how accretion disc conditions are reflected in the RM, we perform polarized radiative transfer calculations using a set of general relativistic magneto-hydrodynamic (GRMHD) simulations appropriate for M87*. We find that in spatially resolved millimetre wavelength images on event horizon scales, the RM can vary by orders of magnitude and even flip sign. The observational consequences of this spatial structure include significant time-variability, sign-flips, and non-$λ^2$ evolution of the polarization plane. For some models, we find that internal rotation measure can cause significant bandwidth depolarization even across the relatively narrow fractional bandwidths observed by the EHT. We decompose the linearly polarized emission in these models based on their RM and find that emission in front of the mid-plane can exhibit orders of magnitude less Faraday rotation than emission originating from behind the mid-plane or within the photon ring. We confirm that the spatially unresolved (i.e., image integrated) RM is a poor predictor of the accretion rate, with substantial scatter stemming from time variability and inclination effects. Models can be constrained with repeated observations to characterise time variability and the degree of non-$λ^2$ evolution of the polarization plane. △ Less

Submitted 4 September, 2020; originally announced September 2020.

Comments: Accepted for publication by MNRAS

Feature Extraction on Synthetic Black Hole Images

Authors: Joshua Yao-Yu Lin, George N. Wong, Ben S. Prather, Charles F. Gammie

Abstract: The Event Horizon Telescope (EHT) recently released the first horizon-scale images of the black hole in M87. Combined with other astronomical data, these images constrain the mass and spin of the hole as well as the accretion rate and magnetic flux trapped on the hole. An important question for EHT is how well key parameters such as spin and trapped magnetic flux can be extracted from present and… ▽ More The Event Horizon Telescope (EHT) recently released the first horizon-scale images of the black hole in M87. Combined with other astronomical data, these images constrain the mass and spin of the hole as well as the accretion rate and magnetic flux trapped on the hole. An important question for EHT is how well key parameters such as spin and trapped magnetic flux can be extracted from present and future EHT data alone. Here we explore parameter extraction using a neural network trained on high resolution synthetic images drawn from state-of-the-art simulations. We find that the neural network is able to recover spin and flux with high accuracy. We are particularly interested in interpreting the neural network output and understanding which features are used to identify, e.g., black hole spin. Using feature maps, we find that the network keys on low surface brightness features in particular. △ Less

Submitted 1 July, 2020; originally announced July 2020.

Comments: Accepted by ML Interpretability for Scientific Discovery Workshop at ICML 2020

Bremsstrahlung in GRMHD models of accreting black holes

Authors: R. Yarza, G. N. Wong, B. R. Ryan, C. F. Gammie

Abstract: The role of bremsstrahlung in the emission from hot accretion flows around slowly accreting supermassive black holes is not thoroughly understood. In order to appraise the importance of bremsstrahlung relative to other radiative processes, we compute spectral energy distributions (SEDs) of accretion disks around slowly accreting supermassive black holes including synchrotron radiation, inverse Com… ▽ More The role of bremsstrahlung in the emission from hot accretion flows around slowly accreting supermassive black holes is not thoroughly understood. In order to appraise the importance of bremsstrahlung relative to other radiative processes, we compute spectral energy distributions (SEDs) of accretion disks around slowly accreting supermassive black holes including synchrotron radiation, inverse Compton scattering, and bremsstrahlung. We compute SEDs for (i) four axisymmetric radiative general relativistic magnetohydrodynamics (RadGRMHD) simulations of $10^{8}M_{\odot}$ black holes with accretion rates between $10^{-8}\dot{M}_{\text{Edd}}$ and $10^{-5}\dot{M}_{\text{Edd}}$, (ii) four axisymmetric RadGRMHD simulations of M87$^\ast$ with varying dimensionless spin $a_\ast$ and black hole mass, and (iii) a 3D GRMHD simulation scaled for Sgr A$^\ast$. At $10^{-8}\dot{M}_{\text{Edd}}$, most of the luminosity is synchrotron radiation, while at $10^{-5}\dot{M}_{\text{Edd}}$ the three radiative processes have similar luminosities. In most models, bremsstrahlung dominates the SED near $512\text{ keV}$. In the M87$^\ast$ models, bremsstrahlung dominates this part of the SED if $a_{\ast} = 0.5$, but inverse Compton scattering dominates if $a_{\ast}= 0.9375$. Since scattering is more variable than bremsstrahlung, this result suggests that $512\text{ keV}$ variability could be a diagnostic of black hole spin. In the appendix, we compare some bremsstrahlung formulae found in the literature. △ Less

Submitted 1 June, 2020; originally announced June 2020.

Comments: 10 pages, 5 figures, accepted for publication in ApJ

Triggered high-mass star formation in the HII region W28A2: A cloud-cloud collision scenario

Authors: Katsuhiro Hayashi, Satoshi Yoshiike, Rei Enokiya, Shinji Fujita, Rin Yamada, Hidetoshi Sano, Kazufumi Torii, Mikito Kohno, Atsushi Nishimura, Akio Ohama, Hiroaki Yamamoto, Kengo Tachihara, Graeme Wong, Nigel Maxted, Catherine Braiding, Gavin Rowell, Michael Burton, Yasuo Fukui

Abstract: We report on a study of the high-mass star formation in the the HII region W28A2 by investigating the molecular clouds extended over ~5-10 pc from the exciting stars using the 12CO and 13CO (J=1-0) and 12CO (J=2-1) data taken by the NANTEN2 and Mopra observations. These molecular clouds consist of three velocity components with the CO intensity peaks at V_LSR ~ -4 km s$^{-1}$, 9 km s$^{-1}$ and 16… ▽ More We report on a study of the high-mass star formation in the the HII region W28A2 by investigating the molecular clouds extended over ~5-10 pc from the exciting stars using the 12CO and 13CO (J=1-0) and 12CO (J=2-1) data taken by the NANTEN2 and Mopra observations. These molecular clouds consist of three velocity components with the CO intensity peaks at V_LSR ~ -4 km s$^{-1}$, 9 km s$^{-1}$ and 16 km s$^{-1}$. The highest CO intensity is detected at V_LSR ~ 9 km s$^{-1}$, where the high-mass stars with the spectral types of O6.5-B0.5 are embedded. We found bridging features connecting these clouds toward the directions of the exciting sources. Comparisons of the gas distributions with the radio continuum emission and 8 um infrared emission show spatial coincidence/anti-coincidence, suggesting physical associations between the gas and the exciting sources. The 12CO J=2-1 to 1-0 intensity ratio shows a high value (> 0.8) toward the exciting sources for the -4 km s$^{-1}$ and +9 km s$^{-1}$ clouds, possibly due to heating by the high-mass stars, whereas the intensity ratio at the CO intensity peak (V_LSR ~ 9 km s$^{-1}$) lowers down to ~0.6, suggesting self absorption by the dense gas in the near side of the +9 km s$^{-1}$ cloud. We found partly complementary gas distributions between the -4 km s$^{-1}$ and +9 km s$^{-1}$ clouds, and the -4 km s$^{-1}$ and +16 km s$^{-1}$ clouds. The exciting sources are located toward the overlapping region in the -4 km s$^{-1}$ and +9 km s$^{-1}$ clouds. Similar gas properties are found in the Galactic massive star clusters, RCW 38 and NGC 6334, where an early stage of cloud collision to trigger the star formation is suggested. Based on these results, we discuss a possibility of the formation of high-mass stars in the W28A2 region triggered by the cloud-cloud collision. △ Less

Submitted 16 May, 2020; originally announced May 2020.

Comments: 26 pages, 16 figures, Accepted for publication in PASJ

Discriminating Accretion States via Rotational Symmetry in Simulated Polarimetric Images of M87

Authors: Daniel C. M. Palumbo, George N. Wong, Ben S. Prather

Abstract: In April 2017, the Event Horizon Telescope observed the shadow of the supermassive black hole at the core of the elliptical galaxy Messier 87. While the original image was constructed from measurements of the total intensity, full polarimetric data were also collected, and linear polarimetric images are expected in the near future. We propose a modal image decomposition of the linear polarization… ▽ More In April 2017, the Event Horizon Telescope observed the shadow of the supermassive black hole at the core of the elliptical galaxy Messier 87. While the original image was constructed from measurements of the total intensity, full polarimetric data were also collected, and linear polarimetric images are expected in the near future. We propose a modal image decomposition of the linear polarization field into basis functions with varying azimuthal dependence of the electric vector position angle. We apply this decomposition to images of ray traced general relativistic magnetohydrodynamics simulations of the Messier 87 accretion disk. For simulated images that are physically consistent with previous observations, the magnitude of the coefficient associated with rotational symmetry, $β_2$, is a useful discriminator between accretion states. We find that at 20 $μ$as resolution, $|β_2|$ is greater than 0.2 only for models of disks with horizon-scale magnetic pressures large enough to disrupt steady accretion. We also find that images with a more radially directed electric vector position angle correspond to models with higher black hole spin. Our analysis demonstrates the utility of the proposed decomposition as a diagnostic framework to improve constraints on theoretical models. △ Less

Submitted 3 April, 2020; originally announced April 2020.

Comments: 11 pages, 7 figures. Accepted into ApJ

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

Authors: F. Roelofs, M. Janssen, I. Natarajan, R. Deane, J. Davelaar, H. Olivares, O. Porth, S. N. Paine, K. L. Bouman, R. P. J. Tilanus, I. M. van Bemmel, H. Falcke, K. Akiyama, A. Alberdi, W. Alef, K. Asada, R. Azulay, A. Baczko, D. Ball, M. Baloković, J. Barrett, D. Bintley, L. Blackburn, W. Boland, G. C. Bower , et al. (183 additional authors not shown)

Abstract: Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabili… ▽ More Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images. △ Less

Submitted 2 April, 2020; originally announced April 2020.

Comments: 20 pages, 15 figures, accepted for publication in A&A

Radio Observations of Supernova Remnant G1.9+0.3

Authors: Kieran J. Luken, Miroslav D. Filipović, Nigel I. Maxted, Roland Kothes, Ray P. Norris, James R. Allison, Rebecca Blackwell, Catherine Braiding, Robert Brose, Michael Burton, Ain Y. De Horta, Tim J. Galvin, Lisa Harvey-Smith, Natasha Hurley-Walker, Denis Leahy, Nicholas O. Ralph, Quentin Roper, Gavin Rowell, Iurii Sushch, Dejan Urošević, Graeme F. Wong

Abstract: We present 1 to 10GHz radio continuum flux density, spectral index, polarisation and Rotation Measure (RM) images of the youngest known Galactic Supernova Remnant (SNR) G1.9+0.3, using observations from the Australia Telescope Compact Array (ATCA). We have conducted an expansion study spanning 8 epochs between 1984 and 2017, yielding results consistent with previous expansion studies of G1.9+0.3.… ▽ More We present 1 to 10GHz radio continuum flux density, spectral index, polarisation and Rotation Measure (RM) images of the youngest known Galactic Supernova Remnant (SNR) G1.9+0.3, using observations from the Australia Telescope Compact Array (ATCA). We have conducted an expansion study spanning 8 epochs between 1984 and 2017, yielding results consistent with previous expansion studies of G1.9+0.3. We find a mean radio continuum expansion rate of ($0.78 \pm 0.09$) per cent year$^{-1}$ (or $\sim8900$ km s$^{-1}$ at an assumed distance of 8.5 kpc), although the expansion rate varies across the SNR perimeter. In the case of the most recent epoch between 2016 and 2017, we observe faster-than-expected expansion of the northern region. We find a global spectral index for G1.9+0.3 of $-0.81\pm0.02$ (76 MHz$-$10 GHz). Towards the northern region, however, the radio spectrum is observed to steepen significantly ($\sim -$1). Towards the two so called (east & west) "ears" of G1.9+0.3, we find very different RM values of 400-600 rad m$^{2}$ and 100-200 rad m$^{2}$ respectively. The fractional polarisation of the radio continuum emission reaches (19 $\pm$ 2)~per~cent, consistent with other, slightly older, SNRs such as Cas~A. △ Less

Submitted 3 December, 2019; originally announced December 2019.

Comments: 17 Pages, 13 Figures

Radio Emission from Interstellar Shocks: Young Type Ia Supernova Remnants and the Case of N 103B in the Large Magellanic Cloud

Authors: R. Z. E. Alsaberi, L. A. Barnes, M. D. Filipovic, N. I. Maxted, H. Sano, G. Rowell, L. M. Bozzetto, S. Gurovich, D. Urovsevic, D. Onic, B. Q. For, P. Manojlovic, G. Wong, T. Galvin, P. Kavanagh, N. Ralph, E. J. Crawford, M. Sasaki, F. Haberl, P. Maggi, N. F. H. Tothil, Y. Fukui

Abstract: Here we present a radio continuum study based on new and archival data from the Australia Telescope Compact Array towards N 103B, a young (<=1000 yrs) spectroscopically confirmed type Ia SNR in the Large Magellanic Cloud (LMC) and proposed to have originated from a single degenerate progenitor. The radio morphology of this SNR is asymmetrical with two bright regions towards the north-west and sout… ▽ More Here we present a radio continuum study based on new and archival data from the Australia Telescope Compact Array towards N 103B, a young (<=1000 yrs) spectroscopically confirmed type Ia SNR in the Large Magellanic Cloud (LMC) and proposed to have originated from a single degenerate progenitor. The radio morphology of this SNR is asymmetrical with two bright regions towards the north-west and south-west of the central location as defined by radio emission. N 103B identified features include: a radio spectral index of -0.75+-0.01 (consistent with other young type Ia SNRs in the Galaxy); a bulk SNR expansion rate as in X-rays; morphology and polarised electrical field vector measurements where we note radial polarisation peak towards the north-west of the remnant at both 5500 and 9000 MHz. The spectrum is concave-up and the most likely reason is the non-linear diffusive shock acceleration effects or presence of two different populations of ultra-relativistic electrons. We also note unpolarized clumps near the south-west region which is in agreement with this above scenario. We derive a typical magnetic field strength for N 103B, of 16.4 microG for an average rotation measurement of 200 rad m^-2. However, we estimate the equipartition field to be of the order of ~235 microG with an estimated minimum energy of Emin=6.3*10^48 erg. The close (~0.5 degree) proximity of N 103B to the LMC mid-plane indicates that an early encounter with dense interstellar medium may have set an important constrain on SNR evolution. Finally, we compare features of N 103B, to six other young type Ia SNRs in the LMC and Galaxy, with a range of proposed degeneracy scenarios to highlight potential differences due to a different models. We suggest that the single degenerate scenario might point to morphologically asymmetric type Ia supernova explosions. △ Less

Submitted 10 November, 2019; originally announced November 2019.

The ASKAP-EMU Early Science Project:Radio Continuum Survey of the Small Magellanic Cloud

Authors: T. D. Joseph, M. D. Filipović, E. J. Crawford, I. Bojičić, E. L. Alexander, G. F. Wong, H. Andernach, H. Leverenz, R. P. Norris, R. Z. E. Alsaberi, C. Anderson, L. A. Barnes, L. M. Bozzetto, F. Bufano, J. D. Bunton, F. Cavallaro, J. D. Collier, H. Dénes, Y. Fukui, T. Galvin, F. Haberl, A. Ingallinera, A. D. Kapinska, B. S. Koribalski, R. Kothes , et al. (31 additional authors not shown)

Abstract: We present two new radio continuum images from the Australian Square Kilometre Array Pathfinder (ASKAP) survey in the direction of the Small Magellanic Cloud (SMC). These images are part of the Evolutionary Map of the Universe (EMU) Early Science Project (ESP) survey of the Small and Large Magellanic Clouds. The two new source lists produced from these images contain radio continuum sources observ… ▽ More We present two new radio continuum images from the Australian Square Kilometre Array Pathfinder (ASKAP) survey in the direction of the Small Magellanic Cloud (SMC). These images are part of the Evolutionary Map of the Universe (EMU) Early Science Project (ESP) survey of the Small and Large Magellanic Clouds. The two new source lists produced from these images contain radio continuum sources observed at 960 MHz (4489 sources) and 1320 MHz (5954 sources) with a bandwidth of 192 MHz and beam sizes of 30.0"x30.0" and 16.3"x15.1", respectively. The median Root Mean Squared (RMS) noise values are 186$μ$Jy beam$^{-1}$ (960 MHz) and 165$μ$Jy beam$^{-1}$ (1320 MHz). To create point source catalogues, we use these two source lists, together with the previously published Molonglo Observatory Synthesis Telescope (MOST) and the Australia Telescope Compact Array (ATCA) point source catalogues to estimate spectral indices for the whole population of radio point sources found in the survey region. Combining our ASKAP catalogues with these radio continuum surveys, we found 7736 point-like sources in common over an area of 30 deg$^2$. In addition, we report the detection of two new, low surface brightness supernova remnant candidates in the SMC. The high sensitivity of the new ASKAP ESP survey also enabled us to detect the bright end of the SMC planetary nebula sample, with 22 out of 102 optically known planetary nebulae showing point-like radio continuum emission. Lastly, we present several morphologically interesting background radio galaxies. △ Less

Submitted 23 September, 2019; originally announced September 2019.

Comments: 18 pages, 23 figues, 9 tables, Accepted for publication in the Monthly Notices of the Royal Astronomical Society

A Supernova Remnant Counterpart for HESS J1832-085

Authors: Nigel I. Maxted, M. D. Filipovic, N. Hurley-Walker, I. Bojicic, G. P. Rowell, F. Haberl, A. J. Ruiter, I. R. Seitenzahl, F. Panther, G. F. Wong, C. Braiding, M. Burton, G. Puhlhofer, H. Sano, Y. Fukui, M. Sasaki, W. Tian, H. Su, X. Cui, D. Leahy, P. J. Hancock

Abstract: We examine the new Galactic supernova remnant (SNR) candidate, G23.11+0.18, as seen by the Murchison Widefield Array (MWA) radio telescope. We describe the morphology of the candidate and find a spectral index of -0.63+/-0.05 in the 70-170MHz domain. A coincident TeV gamma-ray detection in High-Energy Stereoscopic System (HESS) data supports the SNR nature of G23.11+0.18 and suggests that G23.11+0… ▽ More We examine the new Galactic supernova remnant (SNR) candidate, G23.11+0.18, as seen by the Murchison Widefield Array (MWA) radio telescope. We describe the morphology of the candidate and find a spectral index of -0.63+/-0.05 in the 70-170MHz domain. A coincident TeV gamma-ray detection in High-Energy Stereoscopic System (HESS) data supports the SNR nature of G23.11+0.18 and suggests that G23.11+0.18 is accelerating particles beyond TeV energies, thus making this object a promising new cosmic ray hadron source candidate. The remnant cannot be seen in current optical, infrared and X-ray data-sets. We do find, however, a dip in CO-traced molecular gas at a line-of-sight velocity of ~85 km/s, suggesting the existence of a G23.11+0.18 progenitor wind-blown bubble. Furthermore, the discovery of molecular gas clumps at a neighbouring velocity towards HESS J1832-085 adheres to the notion that a hadronic gamma-ray production mechanism is plausible towards the north of the remnant. Based on these morphological arguments, we propose an interstellar medium association for G23.11+0.18 at a kinematic distance of 4.6+/-0.8 kpc. △ Less

Submitted 25 August, 2019; originally announced August 2019.

Comments: 22 pages, 12 figures, Accepted for publication in The Astrophysical Journal

Universal Interferometric Signatures of a Black Hole's Photon Ring

Authors: Michael D. Johnson, Alexandru Lupsasca, Andrew Strominger, George N. Wong, Shahar Hadar, Daniel Kapec, Ramesh Narayan, Andrew Chael, Charles F. Gammie, Peter Galison, Daniel C. M. Palumbo, Sheperd S. Doeleman, Lindy Blackburn, Maciek Wielgus, Dominic W. Pesce, Joseph R. Farah, James M. Moran

Abstract: 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… ▽ More 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. △ Less

Submitted 27 March, 2020; v1 submitted 9 July, 2019; originally announced July 2019.

Comments: 11 pages, 8 figures. Published in Science Advances

Journal ref: Sci. Adv. 6, eaaz1310 (2020)

Murchison Widefield Array and XMM-Newton observations of the Galactic supernova remnant G5.9+3.1

Authors: D. Onić, M. D. Filipović, I. Bojičić, N. Hurley-Walker, B. Arbutina, T. G. Pannuti, C. Maitra, D. Urošević, F. Haberl, N. Maxted, G. F. Wong, G. Rowell, M. E. Bell, J. R. Callingham, K. S. Dwarakanath, B. -Q. For, P. J. Hancock, L. Hindson, M. Johnston-Hollitt, A. D. Kapińska, E. Lenc, B. McKinley, J. Morgan, A. R. Offringa, L. E. Porter , et al. (5 additional authors not shown)

Abstract: In this paper we discuss the radio continuum and X-ray properties of the so-far poorly studied Galactic supernova remnant (SNR) G5.9+3.1. We present the radio spectral energy distribution (SED) of the Galactic SNR G5.9+3.1 obtained with the Murchison Widefield Array (MWA). Combining these new observations with the surveys at other radio continuum frequencies, we discuss the integrated radio contin… ▽ More In this paper we discuss the radio continuum and X-ray properties of the so-far poorly studied Galactic supernova remnant (SNR) G5.9+3.1. We present the radio spectral energy distribution (SED) of the Galactic SNR G5.9+3.1 obtained with the Murchison Widefield Array (MWA). Combining these new observations with the surveys at other radio continuum frequencies, we discuss the integrated radio continuum spectrum of this particular remnant. We have also analyzed an archival XMM-Newton observation, which represents the first detection of X-ray emission from this remnant. The SNR SED is very well explained by a simple power-law relation. The synchrotron radio spectral index of G5.9+3.1, is estimated to be 0.42$\pm$0.03 and the integrated flux density at 1GHz to be around 2.7Jy. Furthermore, we propose that the identified point radio source, located centrally inside the SNR shell, is most probably a compact remnant of the supernova explosion. The shell-like X-ray morphology of G5.9+3.1 as revealed by XMM-Newton broadly matches the spatial distribution of the radio emission, where the radio-bright eastern and western rims are also readily detected in the X-ray while the radio-weak northern and southern rims are weak or absent in the X-ray. Extracted MOS1+MOS2+PN spectra from the whole SNR as well as the north, east, and west rims of the SNR are fit successfully with an optically thin thermal plasma model in collisional ionization equilibrium with a column density N_H~0.80x$10^{22}$ cm$^{-2}$ and fitted temperatures spanning the range kT~0.14-0.23keV for all of the regions. The derived electron number densities n_e for the whole SNR and the rims are also roughly comparable (ranging from ~$0.20f^{-1/2}$ cm$^{-3}$ to ~$0.40f^{-1/2}$ cm$^{-3}$, where f is the volume filling factor). We also estimate the swept-up mass of the X-ray emitting plasma associated with G5.9+3.1 to be ~$46f^{-1/2}M_{\odot}$. △ Less

Submitted 13 April, 2019; originally announced April 2019.

Comments: Accepted for publication in A&A

Journal ref: A&A 625, A93 (2019)

Probing ISM Structure in Trumpler 14 & Carina I Using The Stratospheric Terahertz Observatory 2

Authors: Young Min Seo, Paul F. Goldsmith, Chris Walker, David J. Hollenbach, Mark G. Wolfire, Craig Kulesa, Volker Tolls, Pietro N. Bernasconi, Umit Kavak, Floris F. S. van der Tak, Russ Shipman, Jian Rong Gao, Alexander Tielens, Michael G. Burton, Harold Yorke, Erick Young, William L. Peters, Abram Young, Christopher Groppi, Kristina Davis, Jorge L. Pineda, William D. Langer, Jonathan H. Kawamura, Antony Stark, Gary Melnick , et al. (4 additional authors not shown)

Abstract: We present observations of the Trumpler 14/Carina I region carried out using the Stratospheric Terahertz Observatory 2 (STO2). The Trumpler 14/Carina I region is in the west part of the Carina Nebula Complex, which is one of the most extreme star-forming regions in the Milky Way. We observed Trumpler 14/Carina I in the 158 $μ$m transition of [C\,{\sc ii}] with a spatial resolution of 48$''$ and a… ▽ More We present observations of the Trumpler 14/Carina I region carried out using the Stratospheric Terahertz Observatory 2 (STO2). The Trumpler 14/Carina I region is in the west part of the Carina Nebula Complex, which is one of the most extreme star-forming regions in the Milky Way. We observed Trumpler 14/Carina I in the 158 $μ$m transition of [C\,{\sc ii}] with a spatial resolution of 48$''$ and a velocity resolution of 0.17 km s$^{-1}$. The observations cover a 0.25$^\circ$ by 0.28$^\circ$ area with central position {\it l} = 297.34$^\circ$, {\it b} = -0.60$^\circ$. The kinematics show that bright [C\,{\sc ii}] structures are spatially and spectrally correlated with the surfaces of CO clouds, tracing the photodissociation region and ionization front of each molecular cloud. Along 7 lines of sight that traverse Tr 14 into the dark ridge to the southwest, we find that the [C\,{\sc ii}] luminosity from the HII region is 3.7 times that from the PDR. In same los we find in the PDRs an average ratio of 1:4.1:5.6 for the mass in atomic gas:dark-CO gas: molecular gas traced by CO. Comparing multiple gas tracers including HI 21cm, [C\,{\sc ii}], CO, and radio recombination lines, we find that the HII regions of the Carina Nebula Complex are well-described as HII regions with one-side freely expanding towards us, consistent with the champagne model of ionized gas evolution. The dispersal of the GMC in this region is dominated by EUV photoevaporation; the dispersal timescale is 20-30 Myr. △ Less

Submitted 22 May, 2019; v1 submitted 22 March, 2019; originally announced March 2019.

Comments: ApJ accepted

Discovery of a Pulsar-powered Bow Shock Nebula in the Small Magellanic Cloud Supernova Remnant DEMS5

Authors: Rami Z. E. Alsaberi, C. Maitra, M. D. Filipovi'c, L. M. Bozzetto, F. Haberl, P. Maggi, M. Sasaki, P. Manjolovi'c, V. Velovi'c, P. Kavanagh, N. I. Maxted, D. Urovsevi'c, G. P. Rowell, G. F. Wong, B. -Q. For, A. N. O'Brien, T. J. Galvin, L. Staveley-Smith, R. P. Norris, T. Jarrett, R. Kothes, K. J. Luken, N. Hurley-Walker, H. Sano, D. Oni'c , et al. (10 additional authors not shown)

Abstract: We report the discovery of a new Small Magellanic Cloud Pulsar Wind Nebula (PWN) at the edge of the Supernova Remnant (SNR)-DEM S5. The pulsar powered object has a cometary morphology similar to the Galactic PWN analogs PSR B1951+32 and 'the mouse'. It is travelling supersonically through the interstellar medium. We estimate the Pulsar kick velocity to be in the range of 700-2000 km/s for an age b… ▽ More We report the discovery of a new Small Magellanic Cloud Pulsar Wind Nebula (PWN) at the edge of the Supernova Remnant (SNR)-DEM S5. The pulsar powered object has a cometary morphology similar to the Galactic PWN analogs PSR B1951+32 and 'the mouse'. It is travelling supersonically through the interstellar medium. We estimate the Pulsar kick velocity to be in the range of 700-2000 km/s for an age between 28-10 kyr. The radio spectral index for this SNR PWN pulsar system is flat (-0.29 $\pm$ 0.01) consistent with other similar objects. We infer that the putative pulsar has a radio spectral index of -1.8, which is typical for Galactic pulsars. We searched for dispersion measures (DMs) up to 1000 cm/pc^3 but found no convincing candidates with a S/N greater than 8. We produce a polarisation map for this PWN at 5500 MHz and find a mean fractional polarisation of P $\sim 23$ percent. The X-ray power-law spectrum (Gamma $\sim 2$) is indicative of non-thermal synchrotron emission as is expected from PWN-pulsar system. Finally, we detect DEM S5 in Infrared (IR) bands. Our IR photometric measurements strongly indicate the presence of shocked gas which is expected for SNRs. However, it is unusual to detect such IR emission in a SNR with a supersonic bow-shock PWN. We also find a low-velocity HI cloud of $\sim 107$ km/s which is possibly interacting with DEM S5. SNR DEM S5 is the first confirmed detection of a pulsar-powered bow shock nebula found outside the Galaxy. △ Less

Submitted 11 April, 2019; v1 submitted 7 March, 2019; originally announced March 2019.

Comments: 20 pages, 15 figures, accepted for publication in MNRAS