Hotaka Shiokawa

RADIATIVE MODELS OF SGR A* FROM GRMHD SIMULATIONS

Using flow models based on axisymmetric general relativistic magnetohydrodynamics simulations, we construct radiative models for Sgr A*. Spectral energy distributions (SEDs) 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 T i/T e. We find that (1) models with T i/T e = 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 T i/T e, but the T i/T e = 10 models are generally large and at most marginally consistent with the 230 GHz very long baseline interferometry (VLBI) data; (5) for models with T i/T e = 10 and i = 85° the event horizon is cloaked behind a synchrotron photosphere at 230 GHz and will not be seen by VLBI, but these models overproduce near-infrared 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 * ~ 0.9. We finish with a discussion of the limitations of our model and prospects for future improvements.

General Relativistic Hydrodynamic Simulation of Accretion Flow from a Stellar Tidal Disruption

We study how the matter dispersed when a supermassive black hole tidally disrupts a star joins an accretion flow. Combining a relativistic hydrodynamic simulation of the stellar disruption with a relativistic hydrodynamics simulation of the tidal debris motion, we track such a system until ~80% of the stellar mass bound to the black hole has settled into an accretion flow. Shocks near the stellar pericenter and also near the apocenter of the most tightly-bound debris dissipate orbital energy, but only enough to make the characteristic radius comparable to the semi-major axis of the most-bound material, not the tidal radius as previously thought. The outer shocks are caused by post-Newtonian effects, both on the stellar orbit during its disruption and on the tidal forces. Accumulation of mass into the accretion flow is non-monotonic and slow, requiring ~3--10x the orbital period of the most tightly-bound tidal streams, while the inflow time for most of the mass may be comparable to or longer than the mass accumulation time. Deflection by shocks does, however, remove enough angular momentum and energy from some mass for it to move inward even before most of the mass is accumulated into the accretion flow. Although the accretion rate rises sharply and then decays roughly as a power-law, its maximum is ~0.1x the previous expectation, and the duration of the peak is ~5x longer than previously predicted. The geometric mean of the black hole mass and stellar mass inferred from a measured event timescale is therefore ~0.2x the value given by classical theory.

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

Context. Radiatively ineffi cient accretion flows (RIAFs) are believed to power supermas sive black holes 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 supermassive black hole in the Galactic center, Sgr A*. Aims. 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 an d jets. Methods. Using three-dimensional general relativistic magnetohydrodynamics 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 p arts of the simulations. Results. We find that small changes in the electron temperature can res ult in dramatic differences in the relative dominance of jets and accretion flows. Applic ation to Sgr A* shows that radio spectrum and size of this source can be well reproduced with a model where electrons are more effi ciently 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 m odel also predicts a significant

′Disk Formation Versus Disk Accretion—What Powers Tidal Disruption Events?

A tidal disruption event (TDE) takes place when a star passes near enough to a massive black hole to be disrupted. About half the star’s matter is given elliptical trajectories with large apocenter distances, and the other half is unbound. To form an accretion flow, the bound matter must lose a significant amount of energy, with the actual amount depending on the characteristic scale of the flow measured in units of the black hole’s gravitational radius ( erg). Recent numerical simulations have revealed that the accretion flow scale is close to the scale of the most bound initial orbits, cm from the black hole, and the corresponding energy dissipation rate is erg s−1. We suggest that the energy liberated during the formation of the accretion disk, rather than the energy liberated by subsequent accretion onto the black hole, powers the observed optical TDE candidates. The observed rise times, luminosities, temperatures, emission radii, and line widths seen in these TDEs are all more readily explained in terms of heating associated with disk formation rather than in terms of accretion.

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{sub 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{sub o} and that the power spectra, averaged over long times, likely show a broad bump near f{sub o} rather than distinct, narrow QPO features.

GLOBAL GENERAL RELATIVISTIC MAGNETOHYDRODYNAMIC SIMULATIONS OF BLACK HOLE ACCRETION FLOWS: A CONVERGENCE STUDY

Global, general relativistic magnetohydrodynamic (GRMHD) simulations of non-radiative, 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 four in linear resolution, from 96 × 96 × 64 to 384 × 384 × 256. 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 two 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 non-relativistic global convergence studies of Hawley et al.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 96× 96× 64 to 384×384×256. 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). Department of Astronomy, University of Illinois at Urbana-Champaign, 1002West Green Street, Urbana, IL 61801 Department of Astronomy, University of Illinois at Urbana-Champaign, 1002West Green Street, Urbana, IL 61801; Current address: Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544 Department of Physics, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, IL 61801 Center for Computational Relativity and Gravitation, School of Mathematical Sciences, Rochester Institute of Technology, Rochester, NY 14623

General relativistic magnetohydrodynamical simulations of the jet in M 87

(abridged) The connection between black hole, accretion disk, and radio jet can be best constrained by fitting models to observations of nearby low luminosity galactic nuclei, in particular the well studied sources Sgr~A* and M87. There has been considerable progress in modeling the central engine of active galactic nuclei by an accreting supermassive black hole coupled to a relativistic plasma jet. However, can a single model be applied to a range of black hole masses and accretion rates? Here we want to compare the latest three-dimensional numerical model, originally developed for Sgr A* in the center of the Milky Way, to radio observations of the much more powerful and more massive black hole in M87. We postprocess three-dimensional GRMHD models of a jet-producing radiatively inefficient accretion flow around a spinning black hole using relativistic radiative transfer and ray-tracing to produce model spectra and images. As a key new ingredient to these models, we allow the proton-electron coupling in these simulations depend on the magnetic properties of the plasma. We find that the radio emission in M87 is well described by a combination of a two-temperature accretion flow and a hot single-temperature jet. The model fits the basic observed characteristics of the M87 radio core. The best fit model has a mass-accretion rate of Mdot approx 9x10^{-3} MSUN/YR and a total jet power of P_j \sim 10^{43} erg/s. Emission at 1.3mm is produced by the counter jet close to the event horizon. Its characteristic crescent shape surrounding the black hole shadow could be resolved by future millimeter-wave VLBI experiments. The model was successfully derived from one for the supermassive black hole in center of the Milky Way by appropriately scaling mass and accretion rate. This suggests the possibility that this model could also apply to a larger range of low-luminosity black holes.

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 . 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 . We confirm this numerically and calibrate the scaling relation. This relation is self-consistent in a wedge in parameter space. If 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 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 LX ? X-ray luminosity and M as L 2 X 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.

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.

The Intrinsic Shape of Sagittarius A* at 3.5-mm Wavelength

The radio emission from Sgr A