BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Charles F. Gammie

Charles F. Gammie contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEastro-ph.GAgr-qcastro-ph.IMastro-phhep-th

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

13 published item(s)

preprint2022arXiv

Blandford-Znajek Process in Quadratic Gravity

The Blandford-Znajek process, which uses a magnetized plasma to extract energy from a rotating black hole, is one of the leading candidates for powering relativistic jets. In this work, we investigate the Blandford-Znajek process in two well-motivated quadratic gravity theories: scalar Gauss-Bonnet and dynamical Chern-Simons gravity. We solve analytically for a split-monopole magnetosphere to first order in the small-coupling approximation and second relative order in the slow-rotation approximation. The extracted power at fixed spin and magnetic flux is enhanced in scalar Gauss-Bonnet and reduced in dynamical Chern-Simons gravity, compared to general relativity. We find that there is a degeneracy between spin and the coupling constants of the theories at leading order in the slow rotation approximation that is broken at higher orders.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

PATOKA: Simulating Electromagnetic Observables of Black Hole Accretion

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Spherical accretion in alternative theories of gravity

The groundbreaking image of the black hole at the center of the M87 galaxy has raised questions at the intersection of observational astronomy and black hole physics. How well can the radius of a black hole shadow can be measured, and can this measurement be used to distinguish general relativity from other theories of gravity? We explore these questions using a simple spherical flow model in general relativity, scalar Gauss--Bonnet gravity, and the Rezzolla and Zhidenko parameterized metric. We assume an optically thin plasma with power-law emissivity in radius. Along the way we present a generalized Bondi flow as well as a piecewise-analytic model for the brightness profile of a cold inflow. We use the second moment of a synthetic image as a proxy for EHT observables and compute the ratio of the second moment to the radius of the black hole shadow. We show that corrections to this ratio from modifications to general relativity are subdominant compared to corrections to the critical impact parameter, and argue that this is generally true. We find that astrophysical model parameters are the dominant source of uncertainty in this calculation, emphasizing the importance of understanding the astrophysical model. Given a sufficiently accurate astrophysical model, however, it is possible using measurements of the black hole shadow to distinguish between general relativity and other theories of gravity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Disks as Inhomogeneous, Anisotropic Gaussian Random Fields

We model astrophysical disk surface brightness fluctuations as an inhomogeneous, anisotropic, time-dependent Gaussian random field. The field locally obeys the stochastic partial differential equation of a Matérn field, which has a power spectrum that is flat at large scales and falls off as a power law at small scales. We provide a series of pedagogical examples and along the way provide a convenient parameterization for the local covariance. We then consider two applications to disks. In the first we generate a movie of a disk. In the second, by integrating over a movie of a disk, we generate synthetic light curves and show that the high frequency slope of the resulting power spectrum depends on the local covariance model. We finish with a summary and a brief discussion of other possible astrophysical applications.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Feature Extraction on Synthetic Black Hole Images

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Pair Drizzle around Sub-Eddington Supermassive Black Holes

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Structure of Radiatively Inefficient Black Hole Accretion Flows

We run three long-timescale general-relativistic magnetohydrodynamic simulations of radiatively inefficient accretion flows onto non-rotating black holes. Our aim is to achieve steady-state behavior out to large radii and understand the resulting flow structure. A simulation with adiabatic index Gamma = 4/3 and small initial alternating poloidal magnetic field loops is run to a time of 440,000 GM/c^3, reaching inflow equilibrium inside a radius of 370 GM/c^2. Variations with larger alternating field loops and with Gamma = 5/3 are run to 220,000 GM/c^3, attaining equilibrium out to 170 GM/c^2 and 440 GM/c^2. There is no universal self-similar behavior obtained at radii in inflow equilibrium: the Gamma = 5/3 simulation shows a radial density profile with power law index ranging from -1 in the inner regions to -1/2 in the outer regions, while the others have a power-law slope ranging from -1/2 to close to -2. Both simulations with small field loops reach a state with polar inflow of matter, while the more ordered initial field has polar outflows. However, unbound outflows remove only a factor of order unity of the inflowing material over a factor of ~300 in radius. Our results suggest that the dynamics of radiatively inefficient accretion flows are sensitive to how the flow is fed from larger radii, and may differ appreciably in different astrophysical systems. Millimeter images appropriate for Sgr A* are qualitatively (but not quantitatively) similar in all simulations, with a prominent asymmetric image due to Doppler boosting.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Universal Interferometric Signatures of a Black Hole's Photon Ring

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

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Covariant Radiative Transfer for Black Hole Spacetimes

It has now become possible to study directly, via numerical simulation, the evolution of relativistic, radiation-dominated flows around compact objects. With this in mind we set out explicitly covariant forms of the radiative transfer equation that are suitable for numerical integration in curved spacetime or flat spacetime in curvilinear coordinates. Our work builds on and summarizes in consistent form earlier work by Lindquist, Thorne, Morita and Kaneko, and others. We give explicitly the basic equations in spherical-polar coordinates for Minkowski space and the Kerr spacetime in Kerr-Schild coordinates.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Locality of MHD Turbulence in Isothermal Disks

We numerically evolve turbulence driven by the magnetorotational instability (MRI) in a 3D, unstratified shearing box and study its structure using two-point correlation functions. We confirm Fromang and Papaloizou's result that shearing box models with zero net magnetic flux are not converged; the dimensionless shear stress $α$ is proportional to the grid scale. We find that the two-point correlation of the magnetic field shows that it is composed of narrow filaments that are swept back by differential rotation into a trailing spiral. The correlation lengths along each of the correlation function principal axes decrease monotonically with the grid scale. For mean azimuthal field models, which we argue are more relevant to astrophysical disks than the zero net field models, we find that: $α$ increases weakly with increasing resolution at fixed box size; $α$ increases slightly as the box size is increased; $α$ increases linearly with net field strength, confirming earlier results; the two-point correlation function of the magnetic field is resolved and converged, and is composed of narrow filaments swept back by the shear; the major axis of the two-point increases slightly as the box size is increased; these results are code independent, based on a comparison of ATHENA and ZEUS runs. The velocity, density, and magnetic fields decorrelate over scales larger than $\sim H$, as do the dynamical terms in the magnetic energy evolution equations. We conclude that MHD turbulence in disks is localized, subject to the limitations imposed by the absence of vertical stratification, the use of an isothermal equation of state, finite box size, finite run time, and finite resolution

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Turbulent Magnetic Prandtl Number of MHD Turbulence in Disks

The magnetic Prandtl number Pr_M is the ratio of viscosity to resistivity. In astrophysical disks the diffusion of angular momentum (viscosity) and magnetic fields (resistivity) are controlled by turbulence. Phenomenological models of the evolution of large scale poloidal magnetic fields in disks suggest that the turbulent magnetic Prandtl number Pr_{M,T} controls the rate of escape of vertical field from the disk; for Pr_{M,T} leq R/H vertical field diffuses outward before it can be advected inward by accretion. Here we measure field diffusion and angular momentum transport due to MHD turbulence in a shearing box, and thus Pr_{M,T}, by studying the evolution of a sinusoidal perturbation in the magnetic field that is injected into a turbulent background. We show that the perturbation is always stable, decays approximately exponentially, has decay rate proportional to k^2, and that the implied Pr_{M,T} ~ 1.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Black Hole Spin Evolution

We consider a subset of the physical processes that determine the spin j = a/M of astrophysical black holes. These include: (1) Initial conditions. Recent models suggest that the collapse of supermassive stars are likely to produce black holes with j ~ 0.7. (2) Major mergers. The outcome of a nearly equal mass black hole-black hole merger is not yet known, but we review the current best guesses and analytic bounds. (3) Minor mergers. We recover the result of Blandford & Hughes that accretion of small companions with isotropically distributed orbital angular momenta results in spindown, with j ~ M^{-7/3}. (4) Accretion. We present new results from fully relativistic magnetohydrodynamic accretion simulations. These show that, at least for one sequence of flow models, spin equilibrium (dj/dt = 0) is reached for j ~ 0.9, far less than the canonical value 0.998 of Thorne that was derived in the absence of MHD effects. This equilibrium value may not apply to all accretion flows, particularly thin disks. Nevertheless, it opens the possibility that black holes that have grown primarily through accretion are not maximally rotating.

0 citations0 reviewsNo code linkedOpen access