skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Stochastic Averaging of Radiative Transfer Coefficients for Relativistic Electrons
Abstract Synchrotron emissivities, absorptivities, and Faraday rotation and conversion coefficients are needed in modeling a variety of astrophysical sources, including Event Horizon Telescope (EHT) sources. We develop a method for estimating transfer coefficients that exploits their linear dependence on the electron distribution function, decomposing the distribution function into a sum of parts each of whose emissivity can be calculated easily. We refer to this procedure as stochastic averaging and apply it in two contexts. First, we use it to estimate the emissivity of an isotropicκdistribution function with a high-energy cutoff. The resulting coefficients can be evaluated efficiently enough to be used directly in ray-tracing calculations, and we provide an example calculation. Second, we use stochastic averaging to assess the effect of subgrid turbulence on the volume-averaged emissivity and along the way provide a prescription for a turbulent emissivity. We find that for parameters appropriate to EHT sources turbulence reduces the emissivity slightly. In the infrared, turbulence can dramatically increase the emissivity.  more » « less
Award ID(s):
2034306
PAR ID:
10563626
Author(s) / Creator(s):
;
Publisher / Repository:
Astrophysical Journal
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
968
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
6
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The gradient technique is a promising tool with theoretical foundations based on the fundamental properties of MHD turbulence and turbulent reconnection. Its various incarnations use spectroscopic, synchrotron, and intensity data to trace the magnetic field and measure the media magnetization in terms of Alfvén Mach number. We provide an analytical theory of gradient measurements and quantify the effects of averaging gradients along the line of sight and over the plane of the sky. We derive analytical expressions that relate the properties of gradient distribution with the Alfvén Mach numberMA. We show that these measurements can be combined with measures of sonic Mach number or line broadening to obtain the magnetic field strength. The corresponding technique has advantages to the Davis–Chandrasekhar–Fermi way of obtaining the magnetic field strength. 
    more » « less
  2. Abstract Turbulence quantities in geophysical fluids roughly follow a lognormal distribution. Consequently, extreme values dominate arithmetic means, and it is challenging for regression algorithms to accurately predict point-to-point and mean values. We train neural networks (NNs) to predict logarithmic values of turbulence diffusivityKTfrom multiyear observational records of turbulence in the deep-cycle layer (DCL) of the upper ocean at Earth’s equator, with the objective of accurately predicting long-term statistics ofKT. Depending on prescribed input variables, temporal averaging, and the number of internal parametersNnn, NNs can predict instantaneous values of log10KTwith correlation coefficientsRas high as 0.6–0.65, root-mean-square error less than an order of magnitude, and long-term averages ofKTwithin a factor of 2 between predictions and observations. Prescribing smallNnncompared to the training dataset size results in poor representation of the distribution’s tails. Conversely, prescribing largeNnncauses overfitting degrading the instantaneous predictability. NNs reproduce the observed spread inKTof multiple orders of magnitude at given gradient Richardson number Ri, unlike commonly used physics-based parameterizations which are single-valued functions of Ri. Predictions of the log10of vertical turbulent heat fluxJqare qualitatively similar to those of log10KTbut with poorer correlation because of differences between the observed distributions. Tests for spatial generalizability show that when training on two of three equatorial locations, each having DCLs, with multiyear records (140°, 23°, and 10°W), predictions at the third location are less accurate than when training from the same site. Significance StatementBy enhancing thermodynamic mixing, ocean turbulence transports heat from the surface to the ocean interior. Directly quantifying this transport requires careful, small-scale, long-term observations. Since such observations are rare in both space and time, it is necessary to infer turbulence parameters from conventional measurements like temperature or current velocity. Machine learning predictive algorithms, trained using existing long-term observations, show promise as a means to achieve this goal. A challenge to overcome is how to characterize the lognormal-like distribution of turbulence levels, in which values vary over orders of magnitude and a few extreme values dominate the arithmetic mean. Reproducing this distribution with neural networks is not trivial, and identifying how to do so is a focus of this paper. 
    more » « less
  3. Abstract The Event Horizon Telescope (EHT) has produced horizon-resolving images of Sagittarius A* (Sgr A*). Scattering in the turbulent plasma of the interstellar medium distorts the appearance of Sgr A* on scales only marginally smaller than the fiducial resolution of EHT. The scattering process both diffractively blurs and adds stochastic refractive substructures that limits the practical angular resolution of EHT images of Sgr A*. We explore the ability of a novel recurrent neural network machine learning framework to mitigate these scattering effects, after training on sample data that are agnostic to general relativistic magnetohydrodynamics (GRMHD). We demonstrate that if instrumental limitations are negligible, it is possible to nearly completely mitigate interstellar scattering at a wavelength of 1.3 mm. We validate and quantify the fidelity of this scattering mitigation scheme with physically relevant GRMHD simulations. We find that we can accurately reconstruct resolved structures at the scale of 3μas, well below the nominal instrumental resolution of EHT, 24μas. 
    more » « less
  4. Context. In a series of publications, we describe a comprehensive comparison of Event Horizon Telescope (EHT) data with theoretical models of the observed Sagittarius A* (Sgr A*) and Messier 87* (M87*) horizon-scale sources. Aims. In this article, we report on improvements made to our observational data reduction pipeline and present the generation of observables derived from the EHT models. We make use of ray-traced general relativistic magnetohydrodynamic simulations that are based on different black hole spacetime metrics and accretion physics parameters. These broad classes of models provide a good representation of the primary targets observed by the EHT. Methods. We describe how we combined multiple frequency bands and polarization channels of the observational data to improve our fringe-finding sensitivity and stabilization of atmospheric phase fluctuations. To generate realistic synthetic data from our models, we took the signal path as well as the calibration process, and thereby the aforementioned improvements, into account. We could thus produce synthetic visibilities akin to calibrated EHT data and identify salient features for the discrimination of model parameters. Results. We have produced a library consisting of an unparalleled 962 000 synthetic Sgr A*and M87*datasets. In terms of baseline coverage and noise properties, the library encompasses 2017 EHT measurements as well as future observations with an extended telescope array. Conclusions. We differentiate between robust visibility data products related to model features and data products that are strongly affected by data corruption effects. Parameter inference is mostly limited by intrinsic model variability, which highlights the importance of long-term monitoring observations with the EHT. In later papers in this series, we will show how a Bayesian neural network trained on our synthetic data is capable of dealing with the model variability and extracting physical parameters from EHT observations. With our calibration improvements, our newly reduced EHT datasets have a considerably better quality compared to previously analyzed data. 
    more » « less
  5. Abstract Models of the resolved Event Horizon Telescope (EHT) sources Sgr A* and M87* are constrained by observations at multiple wavelengths, resolutions, polarizations, and time cadences. In this paper, we compare unresolved circular polarization (CP) measurements to a library of models, where each model is characterized by a distribution of CP over time. In the library, we vary the spin of the black hole, the magnetic field strength at the horizon (i.e., both SANE and magnetically arrested disk or MAD models), the observer inclination, a parameter for the maximum ion–electron temperature ratio assuming a thermal plasma, and the direction of the magnetic field dipole moment. We find that Atacama Large Millimeter/submillimeter Array (ALMA) observations of Sgr A* are inconsistent with all edge-on (i= 90°) models. Restricting attention to the MAD models favored by earlier EHT studies of Sgr A*, we find that only models with magnetic dipole moment pointing away from the observer are consistent with ALMA data. We also note that in 26 of the 27 passing MAD models, the accretion flow rotates clockwise on the sky. We provide a table of the means and standard deviations of the CP distributions for all model parameters, along with their trends. 
    more » « less