Source author record

Charles F. Gammie

Charles F. Gammie appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

34works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

34 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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

preprint2016arXiv

A hot big bang theory: magnetic fields and the early evolution of the protolunar disk

The leading theory for the formation of the Earth's moon invokes a collision between a Mars-sized body and the proto-Earth to produce a disk of orbiting material that later condenses to form the Moon. Here we study the early evolution of the protolunar disk. First, we show that the disk opacity is large and cooling is therefore inefficient (t_{cool} Ω>> 1). In this regime angular momentum transport in the disk leads to steady heating unless α< (t_{cool} Ω)^{-1} << 1. Following earlier work by Charnoz and Michaut, and Carballido et al., we show that once the disk is completely vaporized it is well coupled to the magnetic field. We consider a scenario in which turbulence driven by magnetic fields leads to a brief, hot phase where the disk is geometrically thick, wit h strong turbulent mixing. The disk cools by spreading until it decouples from the field. We point out that approximately half the accretion energy is dissipated in the boundary layer where the disk meets the Earth's surface. This creates high entropy material close to the Earth, driving convection and mixing. Finally, a hot, magnetized disk could drive bipolar outflows that remove mass and angular momentum from the Earth-Moon system.

preprint2016arXiv

An Extension of the Athena++ Code Framework for GRMHD Based on Advanced Riemann Solvers and Staggered-Mesh Constrained Transport

We present a new general relativistic magnetohydrodynamics (GRMHD) code integrated into the Athena++ framework. Improving upon the techniques used in most GRMHD codes, ours allows the use of advanced, less diffusive Riemann solvers, in particular HLLC and HLLD. We also employ a staggered-mesh constrained transport algorithm suited for curvilinear coordinate systems in order to maintain the divergence-free constraint of the magnetic field. Our code is designed to work with arbitrary stationary spacetimes in one, two, or three dimensions, and we demonstrate its reliability in a number of tests. We also report on its promising performance and scalability.

preprint2016arXiv

Polarized Synchrotron Emissivities and Absorptivities for Relativistic Thermal, Power-Law, and Kappa Distribution Functions

Synchrotron emission and absorption determine the observational appearance of many astronomical systems. In this paper, we describe a numerical scheme for calculating synchrotron emissivities and absorptivities in all four Stokes parameters for arbitrary gyrotropic electron distribution functions, building on earlier work by Leung, Gammie, and Noble. We use this technique to evaluate the emissivities and the absorptivities for a thermal (Maxwell-Jüttner), isotropic power-law, and isotropic kappa distribution function. The latter contains a power-law tail at high particle energies that smoothly merges with a thermal core at low energies, as is characteristic of observed particle spectra in collisionless plasmas. We provide fitting formulae and error bounds on the fitting formulae for use in codes that solve the radiative transfer equation. The numerical method and the fitting formulae are implemented in a compact C library called ${\tt symphony}$. We find that: the kappa distribution has a source function that is indistinguishable from a thermal spectrum at low frequencies and transitions to the characteristic self-absorbed synchrotron spectrum, $\propto ν^{5/2}$, at high frequency; the linear polarization fraction for a thermal spectrum is near unity at high frequency; and all distributions produce $O(10\%)$ circular polarization at low frequency for lines of sight sufficiently close to the magnetic field vector.

preprint2015arXiv

An Extended Magnetohydrodynamics Model for Relativistic Weakly Collisional Plasmas

Black holes that accrete far below the Eddington limit are believed to accrete through a geometrically thick, optically thin, rotationally supported plasma that we will refer to as a radiatively inefficient accretion flow (RIAF). RIAFs are typically collisionless in the sense that the Coulomb mean free path is large compared to $GM/c^2$, and relativistically hot near the event horizon. In this paper we develop a phenomenological model for the plasma in RIAFs, motivated by the application to sources such as Sgr A* and M87. The model is derived using Israel-Stewart theory, which considers deviations up to second order from thermal equilibrium, but modified for a magnetized plasma. This leads to thermal conduction along magnetic field lines and a difference in pressure, parallel and perpendicular to the field lines (which is equivalent to anisotrotropic viscosity). In the non-relativistic limit, our model reduces to the widely used Braginskii theory of magnetized, weakly collisional plasmas. We compare our model to the existing literature on dissipative relativistic fluids, describe the linear theory of the plasma, and elucidate the physical meaning of the free parameters in the model. We also describe limits of the model when the conduction is saturated and when the viscosity implies a large pressure anisotropy. In future work, the formalism developed in this paper will be used in numerical models of RIAFs to assess the importance of non-ideal processes for the dynamics and radiative properties of slowly accreting black holes.

preprint2015arXiv

bhlight: General Relativistic Radiation Magnetohydrodynamics with Monte Carlo Transport

We present bhlight, a numerical scheme for solving the equations of general relativistic radiation magnetohydrodynamics (GRRMHD) using a direct Monte Carlo solution of the frequency-dependent radiative transport equation. bhlight is designed to evolve black hole accretion flows at intermediate accretion rate, in the regime between the classical radiatively efficient disk and the radiatively inefficient accretion flow (RIAF), in which global radiative effects play a sub-dominant but non-negligible role in disk dynamics. We describe the governing equations, numerical method, idiosyncrasies of our implementation, and a suite of test and convergence results. We also describe example applications to radiative Bondi accretion and to a slowly accreting Kerr black hole in axisymmetry.

preprint2015arXiv

Electron Thermodynamics in GRMHD Simulations of Low-Luminosity Black Hole Accretion

Simple assumptions made regarding electron thermodynamics often limit the extent to which general relativistic magnetohydrodynamic (GRMHD) simulations can be applied to observations of low-luminosity accreting black holes. We present, implement, and test a model that self-consistently evolves an electron entropy equation and takes into account the effects of spatially varying electron heating and relativistic anisotropic thermal conduction along magnetic field lines. We neglect the back-reaction of electron pressure on the dynamics of the accretion flow. Our model is appropriate for systems accreting at $\ll 10^{-5}$ of the Eddington rate, so radiative cooling by electrons can be neglected. It can be extended to higher accretion rates in the future by including electron cooling and proton-electron Coulomb collisions. We present a suite of tests showing that our method recovers the correct solution for electron heating under a range of circumstances, including strong shocks and driven turbulence. Our initial applications to axisymmetric simulations of accreting black holes show that (1)~physically-motivated electron heating rates yield electron temperature distributions significantly different from the constant electron to proton temperature ratios assumed in previous work, with higher electron temperatures concentrated in the coronal region between the disc and the jet; (2)~electron thermal conduction significantly modifies the electron temperature in the inner regions of black hole accretion flows if the effective electron mean free path is larger than the local scale-height of the disc (at least for the initial conditions and magnetic field configurations we study). The methods developed in this work are important for producing more realistic predictions for the emission from accreting black holes such as Sagittarius A* and M87; these applications will be explored in future work.

preprint2015arXiv

Evolution of Accretion Discs around a Kerr Black Hole using Extended Magnetohydrodynamics

Black holes accreting well below the Eddington rate are believed to have geometrically thick, optically thin, rotationally supported accretion discs in which the Coulomb mean free path is large compared to $GM/c^2$. In such an environment, the disc evolution may differ significantly from ideal magnetohydrodynamic predictions. We present non-ideal global axisymmetric simulations of geometrically thick discs around a rotating black hole. The simulations are carried out using a new code ${\rm\it grim}$, which evolves a covariant extended magnetohydrodynamics model derived by treating non-ideal effects as a perturbation of ideal magnetohydrodynamics. Non-ideal effects are modeled through heat conduction along magnetic field lines, and a difference between the pressure parallel and perpendicular to the field lines. The model relies on an effective collisionality in the disc from wave-particle scattering and velocity-space (mirror and firehose) instabilities. We find that the pressure anisotropy grows to match the magnetic pressure, at which point it saturates due to the mirror instability. The pressure anisotropy produces outward angular momentum transport with a magnitude comparable to that of MHD turbulence in the disc, and a significant increase in the temperature in the wall of the jet. We also find that, at least in our axisymmetric simulations, conduction has a small effect on the disc evolution because (1) the heat flux is constrained to be parallel to the field and the field is close to perpendicular to temperature gradients, and (2) the heat flux is choked by an increase in effective collisionality associated with the mirror instability.

preprint2015arXiv

Imaging an Event Horizon: Mitigation of Source Variability of Sagittarius A*

The black hole in the center of the Galaxy, associated with the compact source Sagittarius A* (Sgr A*), is predicted to cast a shadow upon the emission of the surrounding plasma flow, which encodes the influence of general relativity in the strong-field regime. The Event Horizon Telescope (EHT) is a Very Long Baseline Interferometry (VLBI) network with a goal of imaging nearby supermassive black holes (in particular Sgr A* and M87) with angular resolution sufficient to observe strong gravity effects near the event horizon. General relativistic magnetohydrodynamic (GRMHD) simulations show that radio emission from Sgr A* exhibits vari- ability on timescales of minutes, much shorter than the duration of a typical VLBI imaging experiment, which usually takes several hours. A changing source structure during the observations, however, violates one of the basic assumptions needed for aperture synthesis in radio interferometry imaging to work. By simulating realistic EHT observations of a model movie of Sgr A*, we demonstrate that an image of the average quiescent emission, featuring the characteristic black hole shadow and photon ring predicted by general relativity, can nonetheless be obtained by observing over multiple days and subsequent processing of the visibilities (scaling, averaging, and smoothing) before imaging. Moreover, it is shown that this procedure can be combined with an existing method to mitigate the effects of interstellar scattering. Taken together, these techniques allow the black hole shadow in the Galactic center to be recovered on the reconstructed image.

preprint2014arXiv

CARMA Large Area Star Formation Survey: Structure and Kinematics of Dense Gas in Serpens Main

We present observations of N2H+(1-0), HCO+(1-0), and HCN(1-0) toward the Serpens Main molecular cloud from the CARMA Large Area Star Formation Survey (CLASSy). We mapped 150 square arcminutes of Serpens Main with an angular resolution of 7 arcsecs. The gas emission is concentrated in two subclusters (the NW and SE subclusters). The SE subcluster has more prominent filamentary structures and more complicated kinematics compared to the NW subcluster. The majority of gas in the two subclusters has subsonic to sonic velocity dispersions. We applied a dendrogram technique with N2H+(1-0) to study the gas structures; the SE subcluster has a higher degree of hierarchy than the NW subcluster. Combining the dendrogram and line fitting analyses reveals two distinct relations: a flat relation between nonthermal velocity dispersion and size, and a positive correlation between variation in velocity centroids and size. The two relations imply a characteristic depth of 0.15 pc for the cloud. Furthermore, we have identified six filaments in the SE subcluster. These filaments have lengths of 0.2 pc and widths of 0.03 pc, which is smaller than a characteristic width of 0.1 pc suggested by Herschel observations. The filaments can be classified into two types based on their properties. The first type, located in the northeast of the SE subcluster, has larger velocity gradients, smaller masses, and nearly critical mass-per-unit-length ratios. The other type, located in the southwest of the SE subcluster, has the opposite properties. Several YSOs are formed along two filaments which have supercritical mass per unit length ratios, while filaments with nearly critical mass-per-unit-length ratios are not associated with YSOs, suggesting that stars are formed on gravitationally unstable filaments.

preprint2014arXiv

Observational appearance of inefficient accretion flows and jets in 3D GRMHD simulations: Application to Sgr~A*

Radiatively inefficient accretion flows (RIAFs) are believed to power supermassive black holes (SMBH) in the underluminous cores of galaxies. Such black holes are typically accompanied by flat-spectrum radio cores indicating the presence of moderately relativistic jets. One of the best constrained RIAFs is associated with the SMBH in the Galactic center, Sgr A*. Since the plasma in RIAFs is only weakly collisional, the dynamics and the radiative properties of these systems are very uncertain. Here we want to study the impact of varying electron temperature on the appearance of accretion flows and jets. Using 3-D GRMHD accretion flow simulations, we use ray tracing methods to predict spectra and radio images of RIAFs allowing for different electron heating mechanisms in the in- and outflowing parts of the simulations. We find that small changes in the electron temperature can result in dramatic differences in the relative dominance of jets and accretion flows. Application to Sgr A* shows that radio spectrum and size of this source can be well reproduced with a model where electrons are more efficiently heated in the jet. The X-ray emission is sensitive to the electron heating mechanism in the jets and disk and therefore X-ray observations put strong constraints on electron temperatures and geometry of the accretion flow and jet. For Sgr A*, the jet model also predicts a significant frequency-dependent core shift which could place independent constraints on the model once measured accurately. We conclude that more sophisticated models for electron distribution functions are crucial for constraining GRMHD simulations with actual observations. For Sgr A*, the radio appearance may well be dominated by the outflowing plasma. Nonetheless, at the highest radio frequencies, the shadow of the event horizon should still be detectable with future Very Long Baseline Interferometric observations.

preprint2012arXiv

Near-Infrared and X-ray Quasi-Periodic Oscillations in Numerical Models of Sgr A*

We report transient quasi-periodic oscillations (QPOs) on minute timescales in relativistic, radiative models of the galactic center source Sgr A*. The QPOs result from nonaxisymmetric $m=1$ structure in the accretion flow excited by MHD turbulence. Near-infrared (NIR) and X-ray power spectra show significant peaks at frequencies comparable to the orbital frequency at the innermost stable circular orbit (ISCO) $f_o$. The excess power is associated with inward propagating magnetic filaments inside the ISCO. The amplitudes of the QPOs are sensitive to the electron distribution function. We argue that transient QPOs appear at a range of frequencies in the neighborhood of $f_o$ and that the power spectra, averaged over long times, likely show a broad bump near $f_o$ rather than distinct, narrow QPO features.

preprint2012arXiv

The Galactic Center Weather Forecast

In accretion-based models for Sgr A* the X-ray, infrared, and millimeter emission arise in a hot, geometrically thick accretion flow close to the black hole. The spectrum and size of the source depend on the black hole mass accretion rate $\dot{M}$. Since Gillessen et al. have recently discovered a cloud moving toward Sgr A* that will arrive in summer 2013, $\dot{M}$ may increase from its present value $\dot{M}_0$. We therefore reconsider the "best-bet" accretion model of Moscibrodzka et al., which is based on a general relativistic MHD flow model and fully relativistic radiative transfer, for a range of $\dot{M}$. We find that for modest increases in $\dot{M}$ the characteristic ring of emission due to the photon orbit becomes brighter, more extended, and easier to detect by the planned Event Horizon Telescope submm VLBI experiment. If $\dot{M} \gtrsim 8 \dot{M}_0$ this "silhouette of the black hole will be hidden beneath the synchrotron photosphere at 230 GHz, and for $\dot{M} \gtrsim 16 \dot{M}_0$ the silhouette is hidden at 345 GHz. We also find that for $\dot{M} > 2 \dot{M}_0$ the near-horizon accretion flow becomes a persistent X-ray and mid-infrared source, and in the near-infrared Sgr A* will acquire a persistent component that is brighter than currently observed flares.

preprint2011arXiv

Challenges in Forming Planets by Gravitational Instability: Disk Irradiation and Clump Migration, Accretion & Tidal Destructio

We present two-dimensional hydrodynamic simulations of self-gravitating protostellar disks subject to axisymmetric infall from envelopes and irradiation from the central star, to explore disk fragmentation due to gravitational instability (GI), and the fragmented clump evolution. We assume that the disk is built gradually and smoothly by the infall, resulting in good numerical convergence. We confirm that for disks around solar-mass stars, infall at high rates at radii beyond ~50 AU leads to disk fragmentation. At lower infall rates <1e-5 Msun/yr, however, irradiation suppresses fragmentation. We find that, once formed, the fragments or clumps migrate inward on typical type-I time scales of ~2e3 yr initially, but later migration deviates from the type-I time scale when the clump becomes more massive than the local disk mass, and/or when they starts to open gaps. As they migrate, the clumps accrete from the disk at a rate 1e-3 to 1e-1 MJupiter/yr, consistent with analytic estimates that assume a 1-2 Hill radii cross section. Most clumps can grow to their isolation masses >0.1 Msun quickly. The eventual fates of these clumps, however, diverges depending on the migration speed: 3 out of 13 clumps become massive enough (brown dwarf mass) to open gaps in the disk and essentially stop migrating; 4 are tidally destroyed during inward migration; 6 migrate across the inner simulated disk boundary. A simple analytic model for clump evolution is derived to explain these different fates. Overall, our results indicate that fast migration, accretion, and tidal destruction of the clumps pose challenges to the scenario of giant planet formation by GI in situ, but may provide a formation mechanism for close binary systems.

preprint2011arXiv

Global GRMHD Simulations of Black Hole Accretion Flows: a Convergence Study

Global, general relativistic magnetohydrodynamic (GRMHD) simulations of nonradiative, magnetized disks are widely used to model accreting black holes. We have performed a convergence study of GRMHD models computed with HARM3D. The models span a factor of 4 in linear resolution, from 96x96x64 to 384x384x256. We consider three diagnostics of convergence: (1) dimensionless shell-averaged quantities such as plasma β; (2) the azimuthal correlation length of fluid variables; and (3) synthetic spectra of the source including synchrotron emission, absorption, and Compton scattering. Shell-averaged temperature is, except for the lowest resolution run, nearly independent of resolution; shell-averaged plasma β decreases steadily with resolution but shows signs of convergence. The azimuthal correlation lengths of density, internal energy, and temperature decrease steadily with resolution but show signs of convergence. In contrast, the azimuthal correlation length of magnetic field decreases nearly linearly with grid size. We argue by analogy with local models, however, that convergence should be achieved with another factor of 2 in resolution. Synthetic spectra are, except for the lowest resolution run, nearly independent of resolution. The convergence behavior is consistent with that of higher physical resolution local model (shearing box) calculations and with the recent nonrelativistic global convergence studies of Hawley et al. (2011).

preprint2011arXiv

Pair Production in Low Luminosity Galactic Nuclei

Electron-positron pairs may be produced near accreting black holes by a variety of physical processes, and the resulting pair plasma may be accelerated and collimated into a relativistic jet. Here we use a self-consistent dynamical and radiative model to investigate pair production by γγcollisions in weakly radiative accretion flows around a black hole of mass M and accretion rate \dot{M}. Our flow model is drawn from general relativistic magnetohydrodynamic simulations, and our radiation field is computed by a Monte Carlo transport scheme assuming the electron distribution function is thermal. We argue that the pair production rate scales as r^{-6} M^{-1} \dot{M}^{6}. We confirm this numerically and calibrate the scaling relation. This relation is self-consistent in a wedge in M, \dot{M} parameter space. If \dot{M} is too low the implied pair density over the poles of the black hole is below the Goldreich-Julian density and γγpair production is relatively unimportant; if \dot{M} is too high the models are radiatively efficient. We also argue that for a power-law spectrum the pair production rate should scale with the observables L_X \equiv X-ray luminosity and M as L_X^2 M^{-4}. We confirm this numerically and argue that this relation likely holds even for radiatively efficient flows. The pair production rates are sensitive to black hole spin and to the ion-electron temperature ratio which are fixed in this exploratory calculation. We finish with a brief discussion of the implications for Sgr A* and M87.

preprint2010arXiv

Long-term Evolution of Protostellar and Protoplanetary Disks. I. Outbursts

As an initial investigation into the long-term evolution of protostellar disks, we explore the conditions required to explain the large outbursts of disk accretion seen in some young stellar objects. We use one-dimensional time-dependent disk models with a phenomenological treatment of the magnetorotational instability (MRI) and gravitational torques to follow disk evolution over long timescales. Comparison with our previous two-dimensional disk model calculations (Zhu et al. 2009b, Z2009b) indicates that the neglect of radial effects and two-dimensional disk structure in the one-dimensional case makes only modest differences in the results; this allows us to use the simpler models to explore parameter space efficiently. We find that the mass infall rates typically estimated for low-mass protostars generally result in AU-scale disk accretion outbursts, as predicted by our previous analysis (Zhu et al. 2009a,Z2009a). We also confirm quasi-steady accretion behavior for high mass infall rates if the values of $α$-parameter for the magnetorotational instability is small, while at this high accretion rate convection from the thermal instability may lead to some variations. We further constrain the combinations of the $α$-parameter and the MRI critical temperature, which can reproduce observed outburst behavior. Our results suggest that dust sublimation may be connected with full activation of the MRI. This is consistent with the idea that small dust captures ions and electrons to suppress the MRI. In a later paper we will explore both long-term outburst and disk evolution with this model, allowing for infall from protostellar envelopes with differing angular momenta.

preprint2010arXiv

Radially Extended, Stratified, Local Models of Isothermal Disks

We consider local, stratified, numerical models of isothermal accretion disks. The novel feature of our treatment is that radial extent L_x and azimuthal extent L_y satisfy H << L_x, L_y << R, where H is the scale height and R is the local radius. This enables us to probe mesoscale structure in stratified thin disks. We evolve the model at several resolutions, sizes, and initial magnetic field strengths. Consistent with earlier work, we find that the saturated, turbulent state consists of a weakly magnetized disk midplane coupled to a strongly magnetized corona, with a transition at |z| ~ 2H. The saturated α~ 0.01 - 0.02. A two-point correlation function analysis reveals that the central 4H of the disk is dominated by small scale turbulence that is statistically similar to unstratified disk models, while the coronal magnetic fields are correlated on scales ~ 10 H. Nevertheless angular momentum transport through the corona is small. A study of magnetic field loops in the corona reveals few open field lines and predominantly toroidal loops with a characteristic distance between footpoints that is ~ H. Finally we find quasi-periodic oscillations with characteristic timescale ~ 30 Ω^{-1} in the magnetic field energy density. These oscillations are correlated with oscillations in the mean azimuthal field; we present a phenomenological, alpha-dynamo model that captures most aspects of the oscillations.

preprint2010arXiv

The Jet in the Galactic Center: An Ideal Laboratory for Magnetohydrodynamics and General Relativity

In this paper we review and discuss some of the intriguing properties of the Galactic Center supermassive black hole candidate Sgr A*. Of all possible black hole sources, the event horizon of Sgr A*, subtends the largest angular scale on the sky. It is therefore a prime candidate to study and image plasma processes in strong gravity and it even allows imaging of the shadow cast by the event horizon. Recent mm-wave VLBI and radio timing observations as well as numerical GRMHD simulations now have provided several breakthroughs that put Sgr A* back into the focus. Firstly, VLBI observations have now measured the intrinsic size of Sgr A* at multiple frequencies, where the highest frequency measurements have approached the scale of the black hole shadow. Moreover, measurements of the radio variability show a clear time lag between 22 GHz and 43 GHz. The combination of size and timing measurements, allows one to actually measure the flow speed and direction of magnetized plasma at some tens of Schwarzschild radii. This data strongly support a moderately relativistic outflow, consistent with an accelerating jet model. This is compared to recent GRMHD simulation that show the presence of a moderately relativistic outflow coupled to an accretion flow Sgr A*. Further VLBI and timing observations coupled to simulations have the potential to map out the velocity profile from 5-40 Schwarzschild radii and to provide a first glimpse at the appearance of a jet-disk system near the event horizon. Future submm-VLBI experiments would even be able to directly image those processes in strong gravity and directly confirm the presence of an event horizon.

preprint2009arXiv

grmonty: a Monte Carlo Code for Relativistic Radiative Transport

We describe a Monte Carlo radiative transport code intended for calculating spectra of hot, optically thin plasmas in full general relativity. The version we describe here is designed to model hot accretion flows in the Kerr metric and therefore incorporates synchrotron emission and absorption, and Compton scattering. The code can be readily generalized, however, to account for other radiative processes and an arbitrary spacetime. We describe a suite of test problems, and demonstrate the expected $N^{-1/2}$ convergence rate, where $N$ is the number of Monte Carlo samples. Finally we illustrate the capabilities of the code with a model calculation, a spectrum of the slowly accreting black hole Sgr A* based on data provided by a numerical general relativistic MHD model of the accreting plasma.

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

preprint2009arXiv

Radiative Models of Sgr A* from GRMHD Simulations

Using flow models based on axisymmetric general relativistic magnetohydrodynamics (GRMHD) simulations, we construct radiative models for sgr A*. Spectral energy distributions that include the effects of thermal synchrotron emission and absorption, and Compton scattering, are calculated using a Monte Carlo technique. Images are calculated using a ray-tracing scheme. All models are scaled so that the 230 GHz flux density is 3.4 Jy. The key model parameters are the dimensionless black hole spin a*, the inclination i, and the ion-to-electron temperature ratio Ti/Te. We find that: (1) models with Ti/Te=1 are inconsistent with the observed submillimeter spectral slope; (2) the X-ray flux is a strongly increasing function of a*; (3) the X-ray flux is a strongly increasing function of i; (4) 230 GHz image size is a complicated function of i, a*, and Ti/Te, but the Ti/Te = 10 models are generally large and at most marginally consistent with the 230 GHz VLBI data; (5) for models with Ti/Te=10 and i=85 deg the event horizon is cloaked behind a synchrotron photosphere at 230 GHz and will not be seen by VLBI, but these models overproduce NIR and X-ray flux; (6) in all models whose SEDs are consistent with observations the event horizon is uncloaked at 230 GHz; (7) the models that are most consistent with the observations have a* \sim 0.9. We finish with a discussion of the limitations of our model and prospects for future improvements.

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.

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.

preprint1996arXiv

Can Nonlinear Hydromagnetic Waves Support a Self-Gravitating Cloud?

Using self-consistent magnetohydrodynamic (MHD) simulations, we explore the hypothesis that nonlinear MHD waves dominate the internal dynamics of galactic molecular clouds. We employ an isothermal equation of state and allow for self-gravity. We adopt ``slab-symmetry,'' which permits motions $\bf v_\perp$ and fields $\bf B_\perp$ perpendicular to the mean field, but permits gradients only parallel to the mean field. The Alfvén speed $v_A$ exceeds the sound speed $c_s$ by a factor $3-30$. We simulate the free decay of a spectrum of Alfvén waves, with and without self-gravity. We also perform simulations with and without self-gravity that include small-scale stochastic forcing. Our major results are as follows: (1) We confirm that fluctuating transverse fields inhibit the mean-field collapse of clouds when the energy in Alfvén- like disturbances remains comparable to the cloud's gravitational binding energy. (2) We characterize the turbulent energy spectrum and density structure in magnetically-dominated clouds. The spectra evolve to approximately $v_{\perp,\,k}^2\approx B_{\perp,\,k}^2/4πρ\propto k^{-s}$ with $s\sim 2$, i.e. approximately consistent with a ``linewidth-size'' relation $σ_v(R) \propto R^{1/2}$. The simulations show large density contrasts, with high density regions confined in part by the fluctuating magnetic fields. (3) We evaluate the input power required to offset dissipation through shocks, as a function of $c_s/v_A$, the velocity dispersion $σ_v$, and the scale $λ$ of the forcing. In equilibrium, the volume dissipation rate is $5.5(c_s/v_a)^{1/2} (λ/L)^{-1/2}\times ρσ_v^3/L$, for a cloud of linear size $L$ and density $ρ$. (4) Somewhat speculatively, we apply our results to a ``typical'' molecular cloud. The mechanical power input required