Search for: All records

Abstract In the collisionless plasmas of radiatively inefficient accretion flows, heating and acceleration of ions and electrons are not well understood. Recent studies in the gyrokinetic limit revealed the importance of incorporating both the compressive and Alfvénic cascades when calculating the partition of dissipated energy between the plasma species. In this paper, we use a covariant analytic model of the accretion flow to explore the impact of compressive and Alfvénic heating, Coulomb collisions, compressional heating, and radiative cooling on the radial temperature profiles of ions and electrons. We show that, independent of the partition of heat between the plasma species, even a small fraction of turbulent energy dissipated to the electrons makes their temperature scale with a virial profile and the ion-to-electron temperature ratio smaller than in the case of pure Coulomb heating. In contrast, the presence of compressive cascades makes this ratio larger because compressive turbulent energy is channeled primarily into the ions. We calculate the ion-to-electron temperature in the inner accretion flow for a broad range of plasma properties, mass accretion rates, and black hole spins and show that it ranges between 5 ≲T_i/T_e≲ 40. We provide a physically motivated expression for this ratio that can be used to calculate observables from simulations of black hole accretion flows for a wide range of conditions. 

Abstract Observations of the S stars, the cluster of young stars in the inner 0.1 pc of the Galactic center, have been crucial in providing conclusive evidence for a supermassive black hole at the center of our galaxy. Since some of the stars have orbits less than that of a typical human lifetime, it is possible to observe multiple orbits and test the weak-field regime of general relativity. Current calculations of orbits require relatively slow and expensive computations in order to perform numerical integrations for the position and momentum of each star at each observing time. In this paper, we present a computationally efficient, first-order post-Newtonian model for the astrometric and spectroscopic data gathered for the S stars. We find that future, 30 m class telescopes—and potentially even current large telescopes with very high spectroscopic resolution—may be able to detect the Shapiro effect for an S star in the next decade or so. 

Abstract The sparse interferometric coverage of the Event Horizon Telescope (EHT) poses a significant challenge for both reconstruction and model fitting of black hole images.PRIMOis a new principal components analysis-based algorithm for image reconstruction that uses the results of high-fidelity general relativistic, magnetohydrodynamic simulations of low-luminosity accretion flows as a training set. This allows the reconstruction of images that are consistent with the interferometric data and that live in the space of images that is spanned by the simulations.PRIMOfollows Monte Carlo Markov Chains to fit a linear combination of principal components derived from an ensemble of simulated images to interferometric data. We show thatPRIMOcan efficiently and accurately reconstruct synthetic EHT data sets for several simulated images, even when the simulation parameters are significantly different from those of the image ensemble that was used to generate the principal components. The resulting reconstructions achieve resolution that is consistent with the performance of the array and do not introduce significant biases in image features such as the diameter of the ring of emission. 

Abstract The Event Horizon Telescope recently captured images of the supermassive black hole in the center of the M87 galaxy, which shows a ring-like emission structure with the south side only slightly brighter than the north side. This relatively weak asymmetry in the brightness profile along the ring has been interpreted as a consequence of the low inclination of the observer (around 17° for M87), which suppresses the Doppler beaming and boosting effects that might otherwise be expected due to the nearly relativistic velocities of the orbiting plasma. In this work, we use a large suite of general relativistic magnetohydrodynamic simulations to reassess the validity of this argument. By constructing explicit counterexamples, we show that low inclination is a sufficient but not necessary condition for images to have low brightness asymmetry. Accretion flow models with high accumulated magnetic flux close to the black hole horizon (the so-called magnetically arrested disks) and low black hole spins have angular velocities that are substantially smaller than the orbital velocities of test particles at the same location. As a result, such models can produce images with low brightness asymmetry even when viewed edge on. 

Abstract We introduce a new Markov Chain Monte Carlo (MCMC) algorithm with parallel tempering for fitting theoretical models of horizon-scale images of black holes to the interferometric data from the Event Horizon Telescope (EHT). The algorithm implements forms of the noise distribution in the data that are accurate for all signal-to-noise ratios. In addition to being trivially parallelizable, the algorithm is optimized for high performance, achieving 1 million MCMC chain steps in under 20 s on a single processor. We use synthetic data for the 2017 EHT coverage of M87 that are generated based on analytic as well as General Relativistic Magnetohydrodynamic (GRMHD) model images to explore several potential sources of biases in fitting models to sparse interferometric data. We demonstrate that a very small number of data points that lie near salient features of the interferometric data exert disproportionate influence on the inferred model parameters. We also show that the preferred orientations of the EHT baselines introduce significant biases in the inference of the orientation of the model images. Finally, we discuss strategies that help identify the presence and severity of such biases in realistic applications. 

The Event Horizon Telescope recently captured images of the supermassive black hole in the center of the M87 galaxy, which show a ring-like emission structure with the South side only slightly brighter than the North side. This relatively weak asymmetry in the brightness profile along the ring has been interpreted as a consequence of the low inclination of the observer (around 17 deg for M87), which suppresses the Doppler beaming and boosting effects that might otherwise be expected due to the nearly relativistic velocities of the orbiting plasma. In this work, we use a large suite of general relativistic magnetohydrodynamic simulations to reassess the validity of this argument. By constructing explicit counter examples, we show that low-inclination is a sufficient but not necessary condition for images to have low brightness asymmetry. Accretion flow models with high accumulated magnetic flux close to the black hole horizon (the so-called magnetically arrested disks) and low black-hole spins have angular velocities that are substantially smaller than the orbital velocities of test particles at the same location. As a result, such models can produce images with low brightness asymmetry even when viewed edge on. 

Thermal X-ray emission from rotation-powered pulsars is believed to originate from localized “hotspots” on the stellar surface occurring where large-scale currents from the magnetosphere return to heat the atmosphere. Lightcurve modeling has primarily been limited to simple models, such as circular antipodal emitting regions with constant temperature. We calculate more realistic temperature distributions within the polar caps, taking advantage of recent advances in magnetospheric theory, and we consider their effect on the predicted lightcurves. The emitting regions are non-circular even for a pure dipole magnetic field, and the inclusion of an aligned magnetic quadrupole moment introduces a north-south asymmetry. As the quadrupole moment is increased, one hotspot grows in size before becoming a thin ring surrounding the star. For the pure dipole case, moving to the more realistic model changes the lightcurves by $$5-10\%$$ for millisecond pulsars, helping to quantify the systematic uncertainty present in current dipolar models. Including the quadrupole gives considerable freedom in generating more complex lightcurves. We explore whether these simple dipole+quadrupole models can account for the qualitative features of the lightcurve of PSR J0437−4715.