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.

The role of bremsstrahlung in the emission from hot accretion flows around slowly accreting supermassive black holes is not thoroughly understood. In order to appraise the importance of bremsstrahlung relative to other radiative processes, we compute spectral energy distributions (SEDs) of accretion disks around slowly accreting supermassive black holes including synchrotron radiation, inverse Compton scattering, and bremsstrahlung. We compute SEDs for (i) four axisymmetric radiative general relativistic magnetohydrodynamics (RadGRMHD) simulations of $10^{8}M_{\odot}$ black holes with accretion rates between $10^{-8}\dot{M}_{\text{Edd}}$ and $10^{-5}\dot{M}_{\text{Edd}}$, (ii) four axisymmetric RadGRMHD simulations of M87$^\ast$ with varying dimensionless spin $a_\ast$ and black hole mass, and (iii) a 3D GRMHD simulation scaled for Sgr A$^\ast$. At $10^{-8}\dot{M}_{\text{Edd}}$, most of the luminosity is synchrotron radiation, while at $10^{-5}\dot{M}_{\text{Edd}}$ the three radiative processes have similar luminosities. In most models, bremsstrahlung dominates the SED near $512\text{ keV}$. In the M87$^\ast$ models, bremsstrahlung dominates this part of the SED if $a_{\ast} = 0.5$, but inverse Compton scattering dominates if $a_{\ast}= 0.9375$. Since scattering is more variable than bremsstrahlung, this result suggests that $512\text{ keV}$ variability could be a diagnostic of black hole spin. In the appendix, we compare some bremsstrahlung formulae found in the literature.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

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.

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.

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.

High-resolution imaging of supermassive black holes is now possible, with new applications to testing general relativity and horizon-scale accretion and relativistic jet formation processes. Over the coming decade, the EHT will propose to add new strategically placed VLBI elements operating at 1.3mm and 0.87mm wavelength. In parallel, development of next-generation backend instrumentation, coupled with high throughput correlation architectures, will boost sensitivity, allowing the new stations to be of modest collecting area while still improving imaging fidelity and angular resolution. The goal of these efforts is to move from imaging static horizon scale structure to dynamic reconstructions that capture the processes of accretion and jet launching in near real time.

We explore a simple spherical model of optically thin accretion on a Schwarzschild black hole, and study the properties of the image as seen by a distant observer. We show that a dark circular region in the center --- a shadow --- is always present. The outer edge of the shadow is located at the photon ring radius $b_{\rm ph} \equiv \sqrt{27}r_g$, where $r_g=GM/c^2$ is the gravitational radius of the accreting mass $M$. The location of the shadow edge is independent of the inner radius at which the accreting gas stops radiating. The size of the observed shadow is thus a signature of the spacetime geometry and it is hardly influenced by accretion details. We briefly discuss the relevance of these results for the Event Horizon Telescope image of the supermassive black hole in M87.

Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL. Agreement between the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems.

This white paper highlights compact object and fundamental physics science opportunities afforded by high-throughput broadband (0.1-60 keV) X-ray polarization observations. X-ray polarimetry gives new observables with geometric information about stellar remnants which are many orders of magnitude too small for direct imaging. The X-ray polarimetric data also reveal details about the emission mechanisms and the structure of the magnetic fields in and around the most extreme objects in the Universe. Whereas the Imaging X-ray Polarimetry Explorer (IXPE) to be launched in 2021 will obtain first results for bright objects, a follow-up mission could be one order of magnitude more sensitive and would be able to use a broader bandpass to perform physics type experiments for representative samples of sources.

As the closest example of a galactic nucleus, the Galactic center (GC) presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environment. We describe several exciting new research directions that, over the next 10 years, hold the potential to answer some of the biggest scientific questions raised in recent decades: Is General Relativity (GR) the correct description for supermassive black holes? What is the nature of star formation in extreme environments? How do stars and compact objects dynamically interact with the supermassive black hole? What physical processes drive gas accretion in low-luminosity black holes? We describe how the high sensitivity, angular resolution, and astrometric precision offered by the next generation of large ground-based telescopes with adaptive optics will help us answer these questions. First, it will be possible to obtain precision measurements of stellar orbits in the Galaxy's central potential, providing both tests of GR in the unexplored regime near a SMBH and measurements of the extended dark matter distribution that is predicted to exist at the GC. Second, we will probe stellar populations at the GC to significantly lower masses than are possible today, down to brown dwarfs. Their structure and dynamics will provide an unprecedented view of the stellar cusp around the SMBH and will distinguish between models of star formation in this extreme environment. This increase in depth will also allow us to measure the currently unknown population of compact remnants at the GC by observing their effects on luminous sources. Third, uncertainties on the mass of and distance to the SMBH can be improved by a factor of $\sim$10. Finally, we can also study the near-infrared accretion onto the black hole at unprecedented sensitivity and time resolution, which can reveal the underlying physics of black hole accretion.

This white paper, submitted for the Plasma 2020 Decadal Survey, concerns the physics of weakly collisional, high-beta plasmas -- plasmas in which the thermal pressure dominates over the magnetic pressure and in which the inter-particle collision time is comparable to the characteristic timescales of bulk motions. This state of matter, although widespread in the Universe, remains poorly understood: we lack a predictive theory for how it responds to perturbations, how it transports momentum and energy, and how it generates and amplifies magnetic fields. Such topics are foundational to the scientific study of plasmas, and are of intrinsic interest to those who regard plasma physics as a fundamental physics discipline. But these topics are also of extrinsic interest: addressing them directly informs upon our understanding of a wide variety of space and astrophysical systems, including accretion flows around supermassive black holes, the intracluster medium (ICM) between galaxies in clusters, and regions of the near-Earth solar wind. Specific recommendations to advance this field of study are discussed.

We present axisymmetric two-temperature general relativistic radiation magnetohydrodynamic (GRRMHD) simulations of the inner region of the accretion flow onto the supermassive black hole M87. We address uncertainties from previous modeling efforts through inclusion of models for (1) self-consistent dissipative and Coulomb electron heating (2) radiation transport (3) frequency-dependent synchrotron emission, self-absorption, and Compton scattering. We adopt a distance $D=16.7$ Mpc, an observer angle $\theta = 20^{\circ}$, and consider black hole masses $M/M_{\odot} = (3.3\times10^{9}, 6.2\times10^{9})$ and spins $a_{\star} = (0.5, 0.9375)$ in a four-simulation suite. For each $(M, a_{\star})$, we identify the accretion rate that recovers the 230 GHz flux from VLBI measurements. We report on disk thermodynamics at these accretion rates ($\dot{M}/\dot{M}_{\mathrm{Edd}} \sim 10^{-5}$). The disk remains geometrically thick; cooling does not lead to a thin disk component. While electron heating is dominated by Coulomb rather than dissipation for $r \gtrsim 10 GM/c^2$, the accretion disk remains two-temperature. Radiative cooling of electrons is not negligible, especially for $r \lesssim 10 GM/c^2$. The Compton $y$ parameter is of order unity. We then compare derived and observed or inferred spectra, mm images, and jet powers. Simulations with $M/M_{\odot} = 3.3\times10^{9}$ are in conflict with observations. These simulations produce mm images that are too small, while the low-spin simulation also overproduces X-rays. For $M/M_{\odot} = 6.2\times10^{9}$, both simulations agree with constraints on radio/IR/X-ray fluxes and mm image sizes. Simulation jet power is a factor $10^2-10^3$ below inferred values, a possible consequence of the modest net magnetic flux in our models.

Polarized models of relativistically hot astrophysical plasmas require transport coefficients as input: synchrotron absorption and emission coefficients in each of the four Stokes parameters, as well as three Faraday rotation coefficients. Approximations are known for all coefficients for a small set of electron distribution functions, such as the Maxwell-Juttner relativistic thermal distribution, and a general procedure has been obtained by Huang & Shcherbakov for an isotropic distribution function. Here we provide an alternative general procedure, with a full derivation, for calculating absorption and rotation coefficients for an arbitrary isotropic distribution function. Our method involves the computation of the full plasma susceptibility tensor, which in addition to absorption and rotation coefficients may be used to determine plasma modes and the dispersion relation. We implement the scheme in a publicly available library with a simple interface, thus allowing for easy incorporation into radiation transport codes. We also provide a comprehensive survey of the literature and comparison with earlier results.

Sagittarius A* (Sgr A*) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We present an analysis of the most comprehensive NIR variability dataset of Sgr A* to date: eight 24-hour epochs of continuous monitoring of Sgr A* at 4.5 $\mu$m with the IRAC instrument on the Spitzer Space Telescope, 93 epochs of 2.18 $\mu$m data from Naos Conica at the Very Large Telescope, and 30 epochs of 2.12 $\mu$m data from the NIRC2 camera at the Keck Observatory, in total 94,929 measurements. A new approximate Bayesian computation method for fitting the first-order structure function extracts information beyond current Fast Fourier Transformation (FFT) methods of power spectral density (PSD) estimation. With a combined fit of the data of all three observatories, the characteristic coherence timescale of Sgr A* is $\tau_{b} = 243^{+82}_{-57}$ minutes ($90\%$ credible interval). The PSD has no detectable features on timescales down to 8.5 minutes ($95\%$ credible level), which is the ISCO orbital frequency for a dimensionless spin parameter $a = 0.92$. One light curve measured simultaneously at 2.12 and 4.5 $\mu$m during a low flux-density phase gave a spectral index $\alpha_s = 1.6 \pm 0.1$ ($F_\nu \propto \nu^{-\alpha_s}$). This value implies that the Sgr A* NIR color becomes bluer during higher flux-density phases. The probability densities of flux densities of the combined datasets are best fit by log-normal distributions. Based on these distributions, the Sgr A* spectral energy distribution is consistent with synchrotron radiation from a non-thermal electron population from below 20 GHz through the NIR.

Sgr A*, the supermassive black hole (SMBH) at the center of our Milky Way Galaxy, is known to be a variable source of X-ray, near-infrared (NIR), and submillimeter (submm) radiation and therefore a prime candidate to study the electromagnetic radiation generated by mass accretion flow onto a black hole and/or a related jet. Disentangling the power source and emission mechanisms of this variability is a central challenge to our understanding of accretion flows around SMBHs. Simultaneous multiwavelength observations of the flux variations and their time correlations can play an important role in obtaining a better understanding of possible emission mechanisms and their origin. This paper presents observations of two flares that both apparently violate the previously established patterns in the relative timing of submm/NIR/X-ray flares from Sgr A*. One of these events provides the first evidence of coeval structure between NIR and submm flux increases, while the second event is the first example of the sequence of submm/X-ray/NIR flux increases all occurring within ~1 hr. Each of these two events appears to upend assumptions that have been the basis of some analytic models of flaring in Sgr A*. However, it cannot be ruled out that these events, even though unusual, were just coincidental. These observations demonstrate that we do not fully understand the origin of the multiwavelength variability of Sgr A*, and show that there is a continued and important need for long-term, coordinated, and precise multiwavelength observations of Sgr A* to characterize the full range of variability behavior.

We describe ${\tt ipole}$, a new public ray-tracing code for covariant, polarized radiative transport. The code extends the ${\tt ibothros}$ scheme for covariant, unpolarized transport using two representations of the polarized radiation field: in the coordinate frame, it parallel transports the coherency tensor; in the frame of the plasma it evolves the Stokes parameters under emission, absorption, and Faraday conversion. The transport step is implemented to be as spacetime- and coordinate- independent as possible. The emission, absorption, and Faraday conversion step is implemented using an analytic solution to the polarized transport equation with constant coefficients. As a result, ${\tt ipole}$ is stable, efficient, and produces a physically reasonable solution even for a step with high optical depth and Faraday depth. We show that the code matches analytic results in flat space, and that it produces results that converge to those produced by Dexter's ${\tt grtrans}$ polarized transport code on a complicated model problem. We expect ${\tt ipole}$ will mainly find applications in modeling Event Horizon Telescope sources, but it may also be useful in other relativistic transport problems such as modeling for the IXPE mission.

The goal of the Event Horizon Telescope (EHT) is to provide spatially resolved images of Sgr A*, the source associated with the Galactic Center black hole. Because Sgr A* varies on timescales short compared to an EHT observing campaign, it is interesting to ask whether variability contains information about the structure and dynamics of the accretion flow. In this paper, we introduce "time-domain filtering", a technique to filter time fluctuating images with specific temporal frequency ranges, and demonstrate the power and usage of the technique by applying it to mock millimeter wavelength images of Sgr A*. The mock image data is generated from General Relativistic Magnetohydrodynamic (GRMHD) simulation and general relativistic ray-tracing method. We show that the variability on each line of sight is tightly correlated with a typical radius of emission. This is because disk emissivity fluctuates on a timescale of order the local orbital period. Time-domain filtered images therefore reflect the model dependent emission radius distribution, which is not accessible in time-averaged images. We show that, in principle, filtered data have the power to distinguish between models with different black hole spins, different disk viewing angles, and different disk orientations in the sky.

We present axisymmetric numerical simulations of radiatively inefficient accretion flows onto black holes combining general relativity, magnetohydrodynamics, self-consistent electron thermodynamics, and frequency-dependent radiation transport. We investigate a range of accretion rates up to $10^{-5} \dot{M}_{\mathrm{Edd}}$ onto a $10^8 M_{\odot}$ black hole with spin $a_{\star} = 0.5$. We report on averaged flow thermodynamics as a function of accretion rate. We present the spectra of outgoing radiation and find that it varies strongly with accretion rate, from synchrotron-dominated in the radio at low $\dot{M}$ to inverse Compton-dominated at our highest $\dot{M}$. In contrast to canonical analytic models, we find that by $\dot{M} \approx 10^{-5} \dot{M}_{\mathrm{Edd}}$, the flow approaches $\sim 1\%$ radiative efficiency, with much of the radiation due to inverse Compton scattering off Coulomb-heated electrons far from the black hole. These results have broad implications for modeling of accreting black holes across a large fraction of the accretion rates realized in observed systems.

Black holes with accretion rates well below the Eddington rate are expected to be surrounded by low-density, hot, geometrically thick accretion disks. This includes the two black holes being imaged at sub-horizon resolution by the Event Horizon Telescope. In these disks, the mean free path for Coulomb interactions between charged particles is large, and the accreting matter is a nearly collisionless plasma. Despite this, numerical simulations have so far modeled these accretion flows using ideal magnetohydrodynamics. Here, we present the first global, general relativistic, 3D simulations of accretion flows onto a Kerr black hole including the non-ideal effects most likely to affect the dynamics of the disk: the anisotropy between the pressure parallel and perpendicular to the magnetic field, and the heat flux along magnetic field lines. We show that for both standard and magnetically arrested disks, the pressure anisotropy is comparable to the magnetic pressure, while the heat flux remains dynamically unimportant. Despite this large pressure anisotropy, however, the time-averaged structure of the accretion flow is strikingly similar to that found in simulations treating the plasma as an ideal fluid. We argue that these similarities are largely due to the interchangeability of the viscous and magnetic shear stresses as long as the magnetic pressure is small compared to the gas pressure, and to the sub-dominant role of pressure/viscous effects in magnetically arrested disks. We conclude by highlighting outstanding questions in modeling the dynamics of low collisionality accretion flows.

We calculate the radiative properties of Sagittarius A* -- spectral energy distribution, variability, and radio-infrared images -- using the first 3D, physically motivated black hole accretion models that directly evolve the electron thermodynamics in general relativistic MHD simulations. These models reproduce the coupled disc-jet structure for the emission favored by previous phenomenological analytic and numerical works. More specifically, we find that the low frequency radio emission is dominated by emission from a polar outflow while the emission above 100 GHz is dominated by the inner region of the accretion disc. The latter produces time variable near infrared (NIR) and X-ray emission, with frequent flaring events (including IR flares without corresponding X-ray flares and IR flares with weak X-ray flares). The photon ring is clearly visible at 230 GHz and 2 microns, which is encouraging for future horizon-scale observations. We also show that anisotropic electron thermal conduction along magnetic field lines has a negligible effect on the radiative properties of our model. We conclude by noting limitations of our current generation of first-principles models, particularly that the outflow is closer to adiabatic than isothermal and thus underpredicts the low frequency radio emission.

Hot, diffuse, relativistic plasmas such as sub-Eddington black hole accretion flows are expected to be collisionless, yet are commonly modeled as a fluid using ideal general relativistic magnetohydrodynamics (GRMHD). Dissipative effects such as heat conduction and viscosity can be important in a collisionless plasma and will potentially alter the dynamics and radiative properties of the flow from that in ideal fluid models; we refer to models that include these processes as Extended GRMHD. Here we describe a new conservative code, grim, that enables all the above and additional physics to be efficiently incorporated. grim combines time evolution and primitive variable inversion needed for conservative schemes into a single step using an algorithm that only requires the residuals of the governing equations as inputs. This algorithm enables the code to be physics agnostic as well as flexibility regarding time-stepping schemes. grim runs on CPUs, as well as on GPUs, using the same code. We formulate a performance model, and use it to show that our implementation runs optimally on both architectures. grim correctly captures classical GRMHD test problems as well as a new suite of linear and nonlinear test problems with anisotropic conduction and viscosity in special and general relativity. As tests and example applications, we resolve the shock substructure due to the presence of dissipation, and report on relativistic versions of the magneto-thermal instability and heat flux driven buoyancy instability, which arise due to anisotropic heat conduction, and of the firehose instability, which occurs due to anisotropic pressure (i.e. viscosity). Finally, we show an example integration of an accretion flow around a Kerr black hole, using Extended GRMHD.

Magnetohydrodynamic (MHD) turbulence driven by the magnetorotational instability can provide diffusive transport of angular momentum in astrophysical disks, and a widely studied computational model for this process is the ideal, stratified, isothermal shearing box. Here we report results of a convergence study of such boxes up to a resolution of $N = 256$ zones per scale height, performed on blue waters at NCSA with ramses-gpu. We find that the time and vertically integrated dimensionless shear stress $\overline{\alpha} \sim N^{-1/3}$, i.e. the shear stress is resolution dependent. We also find that the magnetic field correlation length decreases with resolution, $\lambda \sim N^{-1/2}$. This variation is strongest at the disk midplane. We show that our measurements of $\alpha$ are consistent with earlier studies. We discuss possible reasons for the lack of convergence.

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} \Omega >> 1). In this regime angular momentum transport in the disk leads to steady heating unless \alpha < (t_{cool} \Omega)^{-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.

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.

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\"uttner), 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 \nu^{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.

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.

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.

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.

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.

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.

Radio and mm-wavelength observations of Sagittarius A* (Sgr A*), the radio source associated with the supermassive black hole at the center of our Galaxy, show that it behaves as a partially self-absorbed synchrotron-emitting source. The measured size of Sgr A* shows that the mm-wavelength emission comes from a small region and consists of the inner accretion flow and a possible collimated outflow. Existing observations of Sgr A* have revealed a time lag between light curves at 43 GHz and 22 GHz, which is consistent with a rapidly expanding plasma flow and supports the presence of a collimated outflow from the environment of an accreting black hole. Here we wish to measure simultaneous frequency-dependent time lags in the light curves of Sgr A* across a broad frequency range to constrain direction and speed of the radio-emitting plasma in the vicinity of the black hole. Light curves of Sgr A* were taken in May 2012 using ALMA at 100 GHz using the VLA at 48, 39, 37, 27, 25.5, and 19 GHz. As a result of elevation limits and the longitude difference between the stations, the usable overlap in the light curves is approximately four hours. Although Sgr A* was in a relatively quiet phase, the high sensitivity of ALMA and the VLA allowed us to detect and fit maxima of an observed minor flare where flux density varied by ~10%. The fitted times of flux density maxima at frequencies from 100 GHz to 19 GHz, as well as a cross-correlation analysis, reveal a simple frequency-dependent time lag relation where maxima at higher frequencies lead those at lower frequencies. Taking the observed size-frequency relation of Sgr A* into account, these time lags suggest a moderately relativistic (lower estimates: 0.5c for two-sided, 0.77c for one-sided) collimated outflow.

We present a statistical analysis of the X-ray flux distribution of Sgr A* from the Chandra X-ray Observatory's 3 Ms Sgr A* X-ray Visionary Project (XVP) in 2012. Our analysis indicates that the observed X-ray flux distribution can be decomposed into a steady quiescent component, represented by a Poisson process with rate $Q=(5.24\pm0.08)\times10^{-3}$ cts s$^{-1},$ and a variable component, represented by a power law process ($dN/dF\propto F^{-\xi},$ $\xi=1.92_{-0.02}^{+0.03}$). This slope matches our recently-reported distribution of flare luminosities. The variability may also be described by a log-normal process with a median unabsorbed 2-8 keV flux of $1.8^{+0.9}_{-0.6}\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$ and a shape parameter $\sigma=2.4\pm0.2,$ but the power law provides a superior description of the data. In this decomposition of the flux distribution, all of the intrinsic X-ray variability of Sgr A* (spanning at least three orders of magnitude in flux) can be attributed to flaring activity, likely in the inner accretion flow. We confirm that at the faint end, the variable component contributes ~10% of the apparent quiescent flux, as previously indicated by our statistical analysis of X-ray flares in these Chandra observations. Our flux distribution provides a new and important observational constraint on theoretical models of Sgr A*, and we use simple radiation models to explore the extent to which a statistical comparison of the X-ray and infrared can provide insights into the physics of the X-ray emission mechanism.

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.

We present details of the CARMA Large Area Star Formation Survey (CLASSy), while focusing on observations of Barnard 1. CLASSy is a CARMA Key Project that spectrally imaged N2H+, HCO+, and HCN (J=1-0 transitions) across over 800 square arcminutes of the Perseus and Serpens Molecular Clouds. The observations have angular resolution near 7" and spectral resolution near 0.16 km/s. We imaged ~150 square arcminutes of Barnard 1, focusing on the main core, and the B1 Ridge and clumps to its southwest. N2H+ shows the strongest emission, with morphology similar to cool dust in the region, while HCO+ and HCN trace several molecular outflows from a collection of protostars in the main core. We identify a range of kinematic complexity, with N2H+ velocity dispersions ranging from ~0.05-0.50 km/s across the field. Simultaneous continuum mapping at 3 mm reveals six compact object detections, three of which are new detections. A new non-binary dendrogram algorithm is used to analyze dense gas structures in the N2H+ position-position-velocity (PPV) cube. The projected sizes of dendrogram-identified structures range from about 0.01-0.34 pc. Size-linewidth relations using those structures show that non-thermal line-of-sight velocity dispersion varies weakly with projected size, while rms variation in the centroid velocity rises steeply with projected size. Comparing these relations, we propose that all dense gas structures in Barnard 1 have comparable depths into the sky, around 0.1-0.2 pc; this suggests that over-dense, parsec-scale regions within molecular clouds are better described as flattened structures rather than spherical collections of gas. Science-ready PPV cubes for Barnard 1 molecular emission are available for download.

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.

Planets appear to form in environments shaped by the gas flowing through protostellar disks to the central young stars. The flows in turn are governed by orbital angular momentum transfer. In this chapter we summarize current understanding of the transfer processes best able to account for the flows, including magneto-rotational turbulence, magnetically-launched winds, self-gravitational instability and vortices driven by hydrodynamical instabilities. For each in turn we outline the major achievements of the past few years and the outstanding questions. We underscore the requirements for operation, especially ionization for the magnetic processes and heating and cooling for the others. We describe the distribution and strength of the resulting flows and compare with the long-used phenomenological $\alpha$-picture, highlighting issues where the fuller physical picture yields substantially different answers. We also discuss the links between magnetized turbulence and magnetically-launched outflows, and between magnetized turbulence and hydrodynamical vortices. We end with a summary of the status of efforts to detect specific signatures of the flows.

Over the last decade, X-ray observations of Sgr A* have revealed a black hole in a deep sleep, punctuated roughly once per day by brief flares. The extreme X-ray faintness of this supermassive black hole has been a long-standing puzzle in black hole accretion. To study the accretion processes in the Galactic Center, Chandra (in concert with numerous ground- and space-based observatories) undertook a 3 Ms campaign on Sgr A* in 2012. With its excellent observing cadence, sensitivity, and spectral resolution, this Chandra X-ray Visionary Project (XVP) provides an unprecedented opportunity to study the behavior of the closest supermassive black hole. We present a progress report from our ongoing study of X-ray flares, including the brightest flare ever seen from Sgr A*. Focusing on the statistics of the flares and the quiescent emission, we discuss the physical implications of X-ray variability in the Galactic Center.

An international consortium is presently constructing a beamformer for the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile that will be available as a facility instrument. The beamformer will aggregate the entire collecting area of the array into a single, very large aperture. The extraordinary sensitivity of phased ALMA, combined with the extremely fine angular resolution available on baselines to the Northern Hemisphere, will enable transformational new very long baseline interferometry (VLBI) observations in Bands 6 and 7 (1.3 and 0.8 mm) and provide substantial improvements to existing VLBI arrays in Bands 1 and 3 (7 and 3 mm). The ALMA beamformer will have impact on a variety of scientific topics, including accretion and outflow processes around black holes in active galactic nuclei (AGN), tests of general relativity near black holes, jet launch and collimation from AGN and microquasars, pulsar and magnetar emission processes, the chemical history of the universe and the evolution of fundamental constants across cosmic time, maser science, and astrometry.

We present the first systematic analysis of the X-ray variability of Sgr A* during the Chandra X-ray Observatory's 2012 Sgr A* X-ray Visionary Project (XVP). With 38 High Energy Transmission Grating Spectrometer (HETGS) observations spaced an average of 7 days apart, this unprecedented campaign enables detailed study of the X-ray emission from this supermassive black hole at high spatial, spectral and timing resolution. In 3 Ms of observations, we detect 39 X-ray flares from Sgr A*, lasting from a few hundred seconds to approximately 8 ks, and ranging in 2-10 keV luminosity from ~1e34 erg/s to 2e35 erg/s. Despite tentative evidence for a gap in the distribution of flare peak count rates, there is no evidence for X-ray color differences between faint and bright flares. Our preliminary X-ray flare luminosity distribution dN/dL is consistent with a power law with index -1.9 (+0.3 -0.4); this is similar to some estimates of Sgr A*'s NIR flux distribution. The observed flares contribute one-third of the total X-ray output of Sgr A* during the campaign, and as much as 10% of the quiescent X-ray emission could be comprised of weak, undetected flares, which may also contribute high-frequency variability. We argue that flares may be the only source of X-ray emission from the inner accretion flow.

We perform calculations of our one-dimensional, two-zone disk model to study the long-term evolution of the circumstellar disk. In particular, we adopt published photoevaporation prescriptions and examine whether the photoevaporative loss alone, coupled with a range of initial angular momenta of the protostellar cloud, can explain the observed decline of the frequency of optically-thick dusty disks with increasing age. In the parameter space we explore, disks have accreting and/or non-accreting transitional phases lasting of $\lesssim20 %$ of their lifetime, which is in reasonable agreement with observed statistics. Assuming that photoevaporation controls disk clearing, we find that initial angular momentum distribution of clouds needs to be weighted in favor of slowly rotating protostellar cloud cores. Again, assuming inner disk dispersal by photoevaporation, we conjecture that this skewed angular momentum distribution is a result of fragmentation into binary or multiple stellar systems in rapidly-rotating cores. Accreting and non-accreting transitional disks show different evolutionary paths on the $\dot{M}-R_{\rm wall}$ plane, which possibly explains the different observed properties between the two populations. However, we further find that scaling the photoevaporation rates downward by a factor of 10 makes it difficult to clear the disks on the observed timescales, showing that the precise value of the photoevaporative loss is crucial to setting the clearing times. While our results apply only to pure photoevaporative loss (plus disk accretion), there may be implications for models in which planets clear disks preferentially at radii of order 10 AU.

We extend the one-dimensional, two-zone models of long-term protostellar disk evolution with infall of Zhu et al. to consider the potential effects of a finite viscosity in regions where the ionization is too low for the magnetorotational instability (MRI) to operate (the "dead zone"). We find that the presence of a small but finite dead zone viscosity, as suggested by simulations of stratified disks with MRI-active outer layers, can trigger inside-out bursts of accretion, starting at or near the inner edge of the disk, instead of the previously-found outside-in bursts with zero dead zone viscosity, which originate at a few AU in radius. These inside-out bursts of accretion bear a qualitative resemblance to the outburst behavior of one FU Ori object, V1515 Cyg, in contrast to the outside-in burst models which more closely resemble the accretion events in FU Ori and V1057 Cyg. Our results suggest that the type and frequency of outbursts are potentially a probe of transport efficiency in the dead zone. Simulations must treat the inner disk regions, $R \lesssim 0.5$ AU, to show the detailed time evolution of accretion outbursts in general and to observe the inside-out bursts in particular.

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.

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.

Measurements of stellar orbits provide compelling evidence that the compact radio source Sagittarius A* at the Galactic Centre is a black hole four million times the mass of the Sun. With the exception of modest X-ray and infrared flares, Sgr A* is surprisingly faint, suggesting that the accretion rate and radiation efficiency near the event horizon are currently very low. Here we report the presence of a dense gas cloud approximately three times the mass of Earth that is falling into the accretion zone of Sgr A*. Our observations tightly constrain the cloud's orbit to be highly eccentric, with an innermost radius of approach of only ~3,100 times the event horizon that will be reached in 2013. Over the past three years the cloud has begun to disrupt, probably mainly through tidal shearing arising from the black hole's gravitational force. The cloud's dynamic evolution and radiation in the next few years will probe the properties of the accretion flow and the feeding processes of the super-massive black hole. The kilo-electronvolt X-ray emission of Sgr A* may brighten significantly when the cloud reaches pericentre. There may also be a giant radiation flare several years from now if the cloud breaks up and its fragments feed gas into the central accretion zone.

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.

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 \beta; (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 \beta\ 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).

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 \gamma\gamma 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 \gamma\gamma 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.

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 \alpha ~ 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 \Omega^{-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.

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.

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 $\alpha$-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 $\alpha$-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.

We use one-dimensional two-zone time-dependent accretion disk models to study the long-term evolution of protostellar disks subject to mass addition from the collapse of a rotating cloud core. Our model consists of a constant surface density magnetically coupled active layer, with transport and dissipation in inactive regions only via gravitational instability. We start our simulations after a central protostar has formed, containing ~ 10% of the mass of the protostellar cloud. Subsequent evolution depends on the angular momentum of the accreting envelope. We find that disk accretion matches the infall rate early in the disk evolution because much of the inner disk is hot enough to couple to the magnetic field. Later infall reaches the disk beyond ~10 AU, and the disk undergoes outbursts of accretion in FU Ori-like events as described in Zhu et al. 2009c. If the initial cloud core is moderately rotating most of the central star's mass is built up by these outburst events. Our results suggest that the protostellar "luminosity problem" is eased by accretion during these FU Ori-like outbursts. After infall stops the disk enters the T Tauri phase. An outer, viscously evolving disk has structure that is in reasonable agreement with recent submillimeter studies and its surface density evolves from $\Sigma \propto R^{-1}$ to $R^{-1.5}$. An inner, massive belt of material-- the "dead zone" -- would not have been observed yet but should be seen in future high angular resolution observations by EVLA and ALMA. This high surface density belt is a generic consequence of low angular momentum transport efficiency at radii where the disk is magnetically decoupled, and would strongly affect planet formation and migration.

We review results from general relativistic axisymmetric magnetohydrodynamic simulations of accretion in Sgr A*. We use general relativistic radiative transfer methods and to produce a broad band (from millimeter to gamma-rays) spectrum. Using a ray tracing scheme we also model images of Sgr A* and compare the size of image to the VLBI observations at 230 GHz. We perform a parameter survey and study radiative properties of the flow models for various black hole spins, ion to electron temperature ratios, and inclinations. We scale our models to reconstruct the flux and the spectral slope around 230 GHz. The combination of Monte Carlo spectral energy distribution calculations and 230 GHz image modeling constrains the parameter space of the numerical models. Our models suggest rather high black hole spin ($a_*\approx 0.9$), electron temperatures close to the ion temperature ($T_i/T_e \sim 3$) and high inclination angles ($i \approx 90 \deg$).

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.

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.

A long standing goal in astrophysics is to directly observe the immediate environment of a black hole with angular resolution comparable to the event horizon. Realizing this goal would open a new window on the study of General Relativity in the strong field regime, accretion and outflow processes at the edge of a black hole, the existence of an event horizon, and fundamental black hole physics (e.g., spin). Steady long-term progress on improving the capability of Very Long Baseline Interferometry (VLBI) at short wavelengths has now made it extremely likely that this goal will be achieved within the next decade. The most compelling evidence for this is the recent observation by 1.3mm VLBI of Schwarzschild radius scale structure in SgrA*, the compact source of radio, submm, NIR and xrays at the center of the Milky Way. SgrA* is thought to mark the position of a ~4 million solar mass black hole, and because of its proximity and estimated mass presents the largest apparent event horizon size of any black hole candidate in the Universe. Over the next decade, existing and planned mm/submm facilities will be combined into a high sensitivity, high angular resolution "Event Horizon Telescope" that will bring us as close to the edge of black hole as we will come for decades. This white paper describes the science case for mm/submm VLBI observations of both SgrA* and M87 (a radio loud AGN of a much more luminous class that SgrA*). We emphasize that while there is development and procurement involved, the technical path forward is clear, and the recent successful observations have removed much of the risk that would normally be associated with such an ambitious project.

We have developed time-dependent models of FU Ori accretion outbursts to explore the physical properties of protostellar disks. Our two-dimensional, axisymmetric models incorporate full vertical structure with a new treatment of the radiative boundary condition for the disk photosphere. We find that FU Ori-type outbursts can be explained by a slow accumulation of matter due to gravitational instability. Eventually this triggers the magnetorotational instability, which leads to rapid accretion. The thermal instability is triggered in the inner disk but this instability is not necessary for the outburst. An accurate disk vertical structure, including convection, is important for understanding the outburst behavior. Large convective eddies develop during the high state in the inner disk. The models are in agreement with Spitzer IRS spectra and also with peak accretion rates and decay timescales of observed outbursts, though some objects show faster rise timescale. We also propose that convection may account for the observed mild-supersonic turbulence and the short-timescale variations of FU Orionis objects.

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.

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 $\alpha$ 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: $\alpha$ increases weakly with increasing resolution at fixed box size; $\alpha$ increases slightly as the box size is increased; $\alpha$ 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

Observations indicate that mass accretion rates onto low-mass protostars are generally lower than the rates of infall to their disks; this suggests that much of the protostellar mass must be accreted during rare, short outbursts of rapid accretion. We explore when protostellar disk accretion is likely to be highly variable. While constant $\alpha$ disks can in principle adjust their accretion rates to match infall rates, protostellar disks are unlikely to have constant $\alpha$. In particular we show that neither models with angular momentum ransport due solely to the magnetorotational instability (MRI) nor ravitational instability (GI) are likely to transport disk mass at rotostellar infall rates over the large range of radii needed to move infalling envelope material down to the central protostar. We show that the MRI and GI are likely to combine to produce outbursts of rapid accretion starting at a few AU. Our analysis is consistent with the time-dependent models of Armitage, Livio, & Pringle (2001) and agrees with our observational study of the outbursting object FU Ori.

We describe three-dimensional general relativistic magnetohydrodynamic simulations of a geometrically thin accretion disk around a non-spinning black hole. The disk has a thickness $h/r\sim0.05-0.1$ over the radial range $(2-20)GM/c^2$. In steady state, the specific angular momentum profile of the inflowing magnetized gas deviates by less than 2% from that of the standard thin disk model of Novikov & Thorne (1973). Also, the magnetic torque at the radius of the innermost stable circular orbit (ISCO) is only $\sim2%$ of the inward flux of angular momentum at this radius. Both results indicate that magnetic coupling across the ISCO is relatively unimportant for geometrically thin disks.

In numerical models of thin astrophysical disks that use an Eulerian scheme, gas orbits supersonically through a fixed grid. As a result the time step is sharply limited by the Courant condition. Also, because the mean flow speed with respect to the grid varies with position, the truncation error varies systematically with position. For hydrodynamic (unmagnetized) disks an algorithm called FARGO has been developed that advects the gas along its mean orbit using a separate interpolation substep. This relaxes the constraint imposed by the Courant condition, which now depends only on the peculiar velocity of the gas, and results in a truncation error that is more nearly independent of position. This paper describes a FARGO-like algorithm suitable for evolving magnetized disks. Our method is second order accurate on a smooth flow and preserves the divergence-free constraint to machine precision. The main restriction is that the magnetic field must be discretized on a staggered mesh. We give a detailed description of an implementation of the code and demonstrate that it produces the expected results on linear and nonlinear problems. We also point out how the scheme might be generalized to make the integration of other supersonic/super-fast flows more efficient. Although our scheme reduces the variation of truncation error with position, it does not eliminate it. We show that the residual position dependence leads to characteristic radial variations in the density over long integrations.

The radio source Sagittarius A* (Sgr A*) is believed to be a hot, inhomogeneous, magnetized plasma flowing near the event horizon of the 3 million solar mass black hole at the galactic center. At a distance of 8000 parsecs the black hole would be among the largest black holes as judged by angular size. Recent observations are consistent with the idea that the millimeter and sub-millimeter photons are dominated by optically thin, thermal synchrotron emission. Anticipating future Very Long Baseline Interferometry (VLBI) observations of Sgr A* at these wavelengths, we present here the first dynamically self-consistent models of millimeter and sub-millimeter emission from Sgr A* based on general relativistic numerical simulations of the accretion flow. Angle-dependent spectra are calculated assuming a thermal distribution of electrons at the baryonic temperature dictated by the simulation and the accretion rate, which acts as a free parameter in our model. The effects of varying model parameters (black hole spin and inclination of the spin to the line of sight) and source variability on the spectrum are shown. We find that the accretion rate value needed to match our calculated millimeter flux to the observed flux is consistent with constraints on the accretion rate inferred from detections of the rotation measure. We also describe the relativistic jet that is launched along the black hole spin axis by the accretion disk and evolves to scales of 1000 gravitational radii.

We consider the formation and evolution of vortices in a hydrodynamic shearing-sheet model. The evolution is done numerically using a version of the ZEUS code. Consistent with earlier results, an injected vorticity field evolves into a set of long-lived vortices, each of which has a radial extent comparable to the local scale height. But we also find that the resulting velocity field has a positive shear stress, <Sigma delta v_r delta v_phi>. This effect appears only at high resolution. The transport, which decays with time as t^-1/2, arises primarily because the vortices drive compressive motions. This result suggests a possible mechanism for angular momentum transport in low-ionization disks, with two important caveats: a mechanism must be found to inject vorticity into the disk, and the vortices must not decay rapidly due to three-dimensional instabilities.

Dynamically significant magnetic fields are routinely observed in molecular clouds, with mass-to-flux ratio lambda = (2 pi sqrt{G}) (Sigma/B) ~ 1 (here Sigma is the total column density and B is the field strength). It is widely believed that ``subcritical'' clouds with lambda < 1 cannot collapse, based on virial arguments by Mestel and Spitzer and a linear stability analysis by Nakano and Nakamura. Here we confirm, using high resolution numerical models that begin with a strongly supersonic velocity dispersion, that this criterion is a fully nonlinear stability condition. All the high-resolution models with lambda <= 0.95 form ``Spitzer sheets'' but collapse no further. All models with lambda >= 1.02 collapse to the maximum numerically resolvable density. We also investigate other factors determining the collapse time for supercritical models. We show that there is a strong stochastic element in the collapse time: models that differ only in details of their initial conditions can have collapse times that vary by as much as a factor of 3. The collapse time cannot be determined from just the velocity dispersion; it depends also on its distribution. Finally, we discuss the astrophysical implications of our results.

Conservative numerical schemes for general relativistic magnetohydrodynamics (GRMHD) require a method for transforming between ``conserved'' variables such as momentum and energy density and ``primitive'' variables such as rest-mass density, internal energy, and components of the four-velocity. The forward transformation (primitive to conserved) has a closed-form solution, but the inverse transformation (conserved to primitive) requires the solution of a set of five nonlinear equations. Here we discuss the mathematical properties of the inverse transformation and present six numerical methods for performing the inversion. The first method solves the full set of five nonlinear equations directly using a Newton-Raphson scheme and a guess from the previous timestep. The other methods reduce the five nonlinear equations to either one or two nonlinear equations that are solved numerically. Comparisons between the methods are made using a survey over phase space, a two-dimensional explosion problem, and a general relativistic MHD accretion disk simulation. The run-time of the methods is also examined. Code implementing the schemes is available for download on the web.

We perform local numerical experiments to investigate the nonlinear stability of thin, radially-stratified disks. We demonstrate the presence of radial convective instability when the disk is nearly in uniform rotation, and show that the net angular momentum transport is slightly inwards, consistent with previous investigations of vertical convection. We then show that a convectively-unstable equilibrium is stabilized by differential rotation. Convective instability is determined by the Richardson number Ri = N_r^2/(q\Omega)^2, where N_r is the radial Brunt-Vaisala frequency and q\Omega is the shear rate. Classical convective instability in a nonshearing medium (Ri -> -infinity) is suppressed when Ri > -1, i.e. when the shear rate becomes greater than the growth rate. Disks with a nearly-Keplerian rotation profile and radial gradients on the order of the disk radius have Ri > -0.01 and are therefore stable to local nonaxisymmetric disturbances. One implication of our results is that the ``baroclinic'' instability recently claimed by Klahr & Bodenheimer is either global or nonexistent. We estimate that our simulations would detect any genuine growth rate > 0.0025\Omega.

We consider the nonaxisymmetric linear theory of radially-stratified disks. We work in a shearing-sheet-like approximation, where the vertical structure of the disk is neglected, and develop equations for the evolution of a plane-wave perturbation comoving with the shear flow (a shearing wave, or ``shwave''). We calculate a complete solution set for compressive and incompressive short-wavelength perturbations in both the stratified and unstratified shearing-sheet models. We develop expressions for the late-time asymptotic evolution of an individual shwave as well as for the expectation value of the energy for an ensemble of shwaves that are initially distributed isotropically in k-space. We find that: (i) incompressive, short-wavelength perturbations in the unstratified shearing sheet exhibit transient growth and asymptotic decay, but the energy of an ensemble of such shwaves is constant with time (consistent with Afshordi, Mukhopadhyay & Narayan 2004); (ii) short-wavelength compressive shwaves grow asymptotically in the unstratified shearing sheet, as does the energy of an ensemble of such shwaves; (iii) incompressive shwaves in the stratified shearing sheet have density and azimuthal velocity perturbations \delta \Sigma,\delta v_y ~ t^{-Ri} (for |Ri| << 1), where Ri = N_x^2/(q \Omega)^2 is the Richardson number, N_x^2 is the square of the radial Brunt-Vaisala frequency and q\Omega is the effective shear rate; (iv) the energy of an ensemble of incompressive shwaves in the stratified shearing sheet behaves asymptotically as Ri t^{1-4Ri} for |Ri| << 1. For Keplerian disks with modest radial gradients, |Ri| is expected to be << 1, and there will therefore be weak growth in a single shwave for Ri < 0 and near-linear growth in the energy of an ensemble of shwaves, independent of the sign of Ri.

The magnetorotational instability (MRI) is the leading candidate for driving turbulence, angular momentum transport, and accretion in astrophysical disks. I consider the linear theory of the MRI in a thin, equatorial disk in the Kerr metric. I begin by analyzing a mechanical model for the MRI that consists of two point masses on nearly circular orbits connected by a spring. I then develop a local Cartesian coordinate system for thin, equatorial Kerr disks. In this local model general relativistic effects manifest themselves solely through changes in the Coriolis parameter and in the tidal expansion of the effective potential. The MRI can be analyzed in the context of the local model using nonrelativistic magnetohydrodynamics, and the growth rates agree with those found in the mechanical model. The maximum growth rate measured by a circular orbit observer differs from a naive estimate using Newtonian gravity by a factor that varies between 1 and 4/3 for all radii and for all a/M.

Some active galactic nuclei, microquasars, and gamma ray bursts may be powered by the electromagnetic braking of a rapidly rotating black hole. We investigate this possibility via axisymmetric numerical simulations of a black hole surrounded by a magnetized plasma. The plasma is described by the equations of general relativistic magnetohydrodynamics, and the effects of radiation are neglected. The evolution is followed for $2000 G M/c^3$, and the computational domain extends from inside the event horizon to typically $40 G M/c^2$. We compare our results to two analytic steady state models, including the force-free magnetosphere of Blandford & Znajek. Along the way we present a self-contained rederivation of the Blandford-Znajek model in Kerr-Schild (horizon penetrating) coordinates. We find that (1) low density polar regions of the numerical models agree well with the Blandford-Znajek model; (2) many of our models have an outward Poynting flux on the horizon in the Kerr-Schild frame; (3) none of our models have a net outward energy flux on the horizon; and (4) one of our models, in which the initial disk has net magnetic flux, shows a net outward angular momentum flux on the horizon. We conclude with a discussion of the limitations of our model, astrophysical implications, and problems to be addressed by future numerical experiments.

We consider the nonlinear outcome of gravitational instability in optically thick disks with a realistic cooling function. We use a numerical model that is local, razor-thin, and unmagnetized. External illumination is ignored. Cooling is calculated from a one-zone model using analytic fits to low temperature Rosseland mean opacities. The model has two parameters: the initial surface density Sigma_0 and the rotation frequency Omega. We survey the parameter space and find: (1) The disk fragments when t_c,eff Omega = 1, where t_c,eff is an effective cooling time defined as the average internal energy of the model divided by the average cooling rate. This is consistent with earlier results that used a simplified cooling function. (2) The initial cooling time t_c0 or a uniform disk with Q = 1 can differ by orders of magnitude from t_c,eff in the nonlinear outcome. The difference is caused by sharp variations in the opacity with temperature. The condition t_c0 Omega = 1 therefore does not necessarily indicate where fragmentation will occur. (3) The largest difference between t_c,eff and t_c0 is near the opacity gap, where dust is absent and hydrogen is largely molecular. (4) In the limit of strong illumination the disk is isothermal; we find that an isothermal version of our model fragments for Q < 1.4. Finally, we discuss some physical processes not included in our model, and find that most are likely to make disks more susceptible to fragmentation. We conclude that disks with t_c,eff Omega < 1 do not exist.

We consider astrophysical maser radiation that is created in the presence of mildly supersonic, magnetohydrodynamic (MHD) turbulence. The focus is on the OH masers for which the magnetic field is strong enough that the separations of the Zeeman components are greater than the spectral linebreadths. A longstanding puzzle has been the absence of the Zeeman $\pi $ components and the high circular polarization in the observed spectra of these masers. We first argue that the elongation of eddies along the field that has recently been recognized in MHD turbulence will enhance the optical depth parallel to the magnetic field in comparison with that perpendicular to the magnetic field. We then simulate maser emission with a numerical model of MHD turbulence to demonstrate quantitatively how the intensities of the linearly polarized $\pi $ components are suppressed and the intensities of the nearly circularly polarized $\sigma $ components are enhanced. This effect is also generic in the sense that most spectral lines in MHD turbulence with Mach number $M\sim 1$ should have larger optical depths parallel to the magnetic field than perpendicular. The effect is reduced considerably when $M<1$. The simulations also demonstrate that the velocity and magnetic field variations due to the turbulence can (but do not necessarily) cause one of the $\sigma $ components to be much more intense than the other, as is often observed for mainline OH masers.

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.

Observations reveal concentrations of molecular line emission on the sky, called ``clumps,'' in dense, star-forming molecular clouds. These clumps are believed to be the eventual sites of star formation. We study the three-dimensional analogs of clumps using a set of self-consistent, time-dependent numerical models of molecular clouds. The models follow the decay of initially supersonic turbulence in an isothermal, self-gravitating, magnetized fluid. We find the following. (1) Clumps are intrinsically triaxial. This explains the observed deficit of clumps with a projected axis ratio near unity, and the apparent prolateness of clumps. (2) Simulated clump axes are not strongly aligned with the mean magnetic field within clumps, nor with the large-scale mean fields. This is in agreement with observations. (3) The clump mass spectrum has a high-mass slope that is consistent with the Salpeter value. There is a low-mass break in the slope at $\sim 0.5 \msun$, although this may depend on model parameters including numerical resolution. (4) The typical specific spin angular momentum of clumps is $4 \times 10^{22} {\rm cm^2 s^{-1}}$. This is larger than the median specific angular momentum of binary stars. Scaling arguments suggest that higher resolution simulations may soon be able to resolve the scales at which the angular momentum of binary stars is determined.

We describe a conservative, shock-capturing scheme for evolving the equations of general relativistic magnetohydrodynamics. The fluxes are calculated using the Harten, Lax, and van Leer scheme. A variant of constrained transport, proposed earlier by T\'oth, is used to maintain a divergence free magnetic field. Only the covariant form of the metric in a coordinate basis is required to specify the geometry. We describe code performance on a full suite of test problems in both special and general relativity. On smooth flows we show that it converges at second order. We conclude by showing some results from the evolution of a magnetized torus near a rotating black hole.

We report on a numerical study of viscous fluid accretion onto a black hole. The flow is axisymmetric and uses a pseudo-Newtonian potential to model relativistic effects near the event horizon. The numerical method is a variant of the ZEUS code. As a test of our numerical scheme, we are able to reproduce results from earlier, similar work by Igumenshchev and Abramowicz and Stone et al. We consider models in which mass is injected onto the grid as well as models in which an initial equilibrium torus is accreted. In each model we measure three ``eigenvalues'' of the flow: the accretion rate of mass, angular momentum, and energy. We find that the eigenvalues are sensitive to r_{in}, the location of the inner radial boundary. Only when the flow is always supersonic on the inner boundary are the eigenvalues insensitive to small changes in r_{in}. We also report on the sensitivity of the results to other numerical parameters.

Thin, Keplerian accretion disks generically become gravitationally unstable at large radius. I investigate the nonlinear outcome of such instability in cool disks using razor-thin, local, numerical models. Cooling, characterized by a constant cooling time t_c, drives the instability. I show analytically that, if the disk can reach a steady state in which heating by dissipation of turbulence balances cooling, then the dimensionless angular momentum flux density \alpha = ((9/4) \gamma (\gamma-1) \Omega t_c)^{-1}. Numerical experiments show that: (1) if t_c \gtrsim 3\Omega^{-1} then the disk reaches a steady, gravito-turbulent state in which Q \sim 1 and cooling is balanced by heating due to dissipation of turbulence; (2) if t_c \lesssim 3\Omega^{-1}, then the disk fragments, possibly forming planets or stars; (3) in a steady, gravito-turbulent state, surface density structures have a characteristic physical scale \sim 64 G \Sigma/\Omega^2 that is independent of the size of the computational domain.

We use 3D numerical MHD simulations to follow the evolution of cold, turbulent, gaseous systems with parameters representing GMC conditions. We study three cloud simulations with varying mean magnetic fields, but identical initial velocity fields. We show that turbulent energy is reduced by a factor two after 0.4-0.8 flow crossing times (2-4 Myr), and that the magnetically supercritical cloud models collapse after ~6 Myr, while the subcritical cloud does not collapse. We compare density, velocity, and magnetic field structure in three sets of snapshots with matched Mach numbers. The volume and column densities are both log-normally distributed, with mean volume density a factor 3-6 times the unperturbed value, but mean column density only a factor 1.1-1.4 times the unperturbed value. We use a binning algorithm to investigate the dependence of kinetic quantities on spatial scale for regions of column density contrast (ROCs). The average velocity dispersion for the ROCs is only weakly correlated with scale, similar to the mean size-linewidth relation for clumps within GMCs. ROCs are often superpositions of spatially unconnected regions that cannot easily be separated using velocity information; the same difficulty may affect observed GMC clumps. We analyze magnetic field structure, and show that in the high density regime, total magnetic field strengths increase with density with logarithmic slope 1/3 -2/3. Mean line-of-sight magnetic field strengths vary widely across a projected cloud, and do not correlate with column density. We compute simulated interstellar polarization maps at varying orientations, and determine that the Chandrasekhar-Fermi formula multiplied by a factor ~0.5 yields a good estimate of the plane-of sky magnetic field strength provided the dispersion in polarization angles is < 25 degrees.

We consider the stability of warping modes in Keplerian discs. We find them to be parametrically unstable using two lines of attack, one based on three-mode couplings and the other on Floquet theory. We confirm the existence of the instability, and investigate its nonlinear development in three dimensions, via numerical experiment. The most rapidly growing non-axisymmetric disturbances are the most nearly axisymmetric (low m) ones. Finally, we offer a simple, somewhat speculative model for the interaction of the parametric instability with the warp. We apply this model to the masing disc in NGC 4258 and show that, provided the warp is not forced too strongly, parametric instability can fix the amplitude of the warp.

The efficiency of thin disk accretion onto black holes depends on the inner boundary condition, specifically the torque applied to the disk at the last stable orbit. This is usually assumed to vanish. I estimate the torque on a magnetized disk using a steady magnetohydrodynamic inflow model originally developed by Takahashi et al., 1990. I find that the efficiency epsilon can depart significantly from the classical thin disk value. In some cases epsilon > 1, i.e. energy is extracted from the black hole.

We review recent results from numerical simulations and related models of MHD turbulence in the interstellar medium (ISM) and in molecular clouds. We discuss the implications of turbulence for the processes of cloud formation and evolution, and the determination of clouds' physical properties. Numerical simulations of the turbulent ISM to date have included magnetic fields, self-gravity, parameterized heating and cooling, modeled star formation and other turbulent inputs. The structures which form reproduce well observed velocity-size scaling properties, while predicting the non-existence of a general density-size scaling law. Criteria for the formation of gravitationally-bound structures by turbulent compression are summarized. For flows with equations of state $P\propto \rho^\gamma$, the statistics of the density field depend on the exponent $\gamma$. Numerical simulations of both forced and decaying MHD compressible turbulence have shown that the decay rate is comparable to the non-magnetic case. For virialized clouds, the turbulent decay time is shorter than the gravitational free-fall time, so wholesale cloud collapse is only prevented by ongoing turbulent inputs and/or a strong mean magnetic field. Finally, perspectives for future work in this field are briefly discussed.

The molecular component of the Galaxy is comprised of turbulent, magnetized clouds, many of which are self-gravitating and form stars. To understand how these clouds' evolution may depend on their level of turbulence, mean magnetization, and degree of self-gravity, we perform a survey of direct numerical MHD simulations in 2.5 dimensions. Two of our independent survey parameters allow us to model clouds which either meet or fail conditions for magneto-Jeans stability and magnetic criticality. Our third survey parameter allows us to initiate turbulence of either sub- or super-Alfv\'enic amplitude. We evaluate the times for each cloud model to become gravitationally bound, and measure each model's kinetic energy loss over the fluid flow crossing time. We compare the evolution of density and magnetic field structural morphology, and quantify the differences in the density contrast generated by internal stresses, for models of differing mean magnetization. We find that the values of (c_s/v_A)^2 and L/L_Jeans, but not the initial Mach number v/c_s, determine the time for cloud gravitational binding and collapse. We find, contrary to some previous expectations, less than a factor of two difference between turbulent decay times for models with varying magnetic field strength. In all models, we find turbulent amplification in the magnetic field strength, with the turbulent magnetic energy between 25-60% of the turbulent kinetic energy after one flow crossing time. We find that for non-self-gravitating stages of evolution, the mass-averaged density contrast magnitudes <log(\rho/\bar\rho)> are in the range 0.2-0.5; we note that only the more strongly-magnetized models appear consistent with estimates of clump/interclump density contrasts inferred in Galactic GMCs.

We report results of a three dimensional, high resolution (up to 512^3) numerical investigation of supersonic compressible magnetohydrodynamic turbulence. We consider both forced and decaying turbulence. The model parameters are appropriate to conditions found in Galactic molecular clouds. We find that the dissipation time of turbulence is of order the flow crossing time or smaller, even in the presence of strong magnetic fields. About half the dissipation occurs in shocks. Weak magnetic fields are amplified and tangled by the turbulence, while strong fields remain well ordered.

We consider the implications of recent infrared and radio observations of the nucleus of NGC 4258. There is no direct evidence that the nucleus has been steadily accreting on the viscous timescale of the outer masing disk, which is >~ 10^9 \yr. Thus the mass accretion rate in the outer disk need not be the same as in the inner accretion flow where most of the gravitational binding energy is released. We show that an advection-dominated flow model with a transition radius of \sim (10-100) M (where M is the mass of the hole) and \mdot \approx 10^{-2}\msun \yr^{-1} is consistent with the observed spectrum from radio to X-rays. We also show that a thin (flat or warped) disk can fit the observed fluxes outside the X-ray band. The X-rays can be explained by means of a corona in such a model, but the absence of radio emission from the location of the putative central black hole provides a serious constraint on the properties of the corona. A wide range of accretion rates, 10^{-4} <~ \mdot <~ 10^{-2} \msun\yr^{-1}, can be made to fit the data, but the most ``natural'' models have 10^{-3} <~ \mdot <~ 10^{-2} \msun\yr^{-1}. Physical conditions in the observed VLBI jet features can also be related to conditions in the inner accretion flow. We conclude with a list of future observations that might help to constrain the accretion rate.

In a previous paper we have written down equations describing steady-state, optically thin, advection-dominated accretion onto a Kerr black hole (Gammie and Popham 1997, hereafter Paper I). In this paper we survey the numerical solutions to these equations. We find that the temperature and density of the gas in the inner part of the accretion flow depend strongly on the black hole spin parameter $a$. The rate of angular momentum accretion is also shown to depend on $a$; for $a$ greater than an equilibrium spin parameter $a_{eq}$ the black hole is de-spun by the accretion flow. We also investigate the dependence of the flow on the angular momentum transport efficiency $\alpha$, the advected fraction of the dissipated energy $f$, and the adiabatic index $\gamma$. We find solutions for $-1 < a < 1$, $10^{-4} \le \alpha \le 0.44$, $0.01 \le f \le 1$, and $4/3 < \gamma < 5/3$. For low values of $\alpha$ and $f$ the inner part of the flow exhibits a pressure maximum and appears similar to equilibrium thick disk solutions.

We show that radiation dominated accretion discs are likely to suffer from a ``photon bubble'' instability similar to that described by Arons in the context of accretion onto neutron star polar caps. The instability requires a magnetic field for its existence. In an asymptotic regime appropriate to accretion discs, we find that the overstable modes obey the remarkably simple dispersion relation \omega^2 = -i g k F(B,k). Here g is the vertical gravitational acceleration, B the magnetic field, and F is a geometric factor of order unity that depends on the relative orientation of the magnetic field and the wavevector. In the nonlinear outcome it seems likely that the instability will enhance vertical energy transport and thereby change the structure of the innermost parts of relativistic accretion discs.

I review recent numerical studies of accretion disks, focusing on measurements of the turbulent shear stress, or $\alpha$, in the shearing box model. I conclude with a list of astronomically relevant open questions that can be settled via future numerical experiments.

We show that dwarf nova disks in quiescence have rather low magnetic Reynolds number, of order 10^3. Numerical simulations of magnetized accretion disks suggest that under these conditions magnetohydrodynamic turbulence and the associated angular momentum transport is sharply reduced. This could be the physical origin of episodic accretion in dwarf nova disks. If so, the standard disk instability model needs to be revised.

Sgr A* at the Galactic Center is a puzzling source. It has a mass M=(2.5+/-0.4) x 10^6 solar masses which makes it an excellent black hole candidate. Observations of stellar winds and other gas flows in its vicinity suggest a mass accretion rate approximately few x 10^{-6} solar masses per year. However, such an accretion rate would imply a luminosity > 10^{40} erg/s if the radiative efficiency is the usual 10 percent, whereas observations indicate a bolometric luminosity <10^{37} erg/s. The spectrum of Sgr A* is unusual, with emission extending over many decades of wavelength. We present a model of Sgr A* which is based on a two-temperature optically-thin advection-dominated accretion flow. The model is consistent with the estimated mass and accretion rate, and fits the observed fluxes in the cm/mm and X-ray bands as well as upper limits in the sub-mm and infrared bands; the fit is less good in the radio below 86 GHz and in gamma-rays above 100 MeV. The very low luminosity of Sgr A* is explained naturally in the model by means of advection. Most of the viscously dissipated energy is advected into the central mass by the accreting gas, and therefore the radiative efficiency is extremely low, approximately 5 x 10^{-6}. A critical element of the model is the presence of an event horizon at the center which swallows the advected energy. The success of the model could thus be viewed as confirmation that Sgr A* is a black hole.

We write down and solve equations describing steady state, optically thin, advection-dominated accretion onto a Kerr black hole. The mean flow, described by the relativistic fluid equations, is axisymmetric and vertically averaged. The effect of turbulence in the flow is represented by a viscous shear stress. Our treatment differs in several important ways from earlier work: we use a causal prescription for the shear stress, we do not assume the relativistic enthalpy is unity (this is important for rapidly rotating holes), and we use a relativistic equation of state. We present several representative solutions and use them to evaluate the importance of relativistic effects, to check our approximations, and to evaluate the robustness of the input physics. Detailed properties of the solutions are described in an accompanying paper.

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\'en speed $v_A$ exceeds the sound speed $c_s$ by a factor $3-30$. We simulate the free decay of a spectrum of Alfv\'en 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\'en- 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\pi\rho\propto k^{-s}$ with $s\sim 2$, i.e. approximately consistent with a ``linewidth-size'' relation $\sigma_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 $\sigma_v$, and the scale $\lambda$ of the forcing. In equilibrium, the volume dissipation rate is $5.5(c_s/v_a)^{1/2} (\lambda/L)^{-1/2}\times \rho \sigma_v^3/L$, for a cloud of linear size $L$ and density $\rho$. (4) Somewhat speculatively, we apply our results to a ``typical'' molecular cloud. The mechanical power input required