Application of the Cubed-Sphere Grid to Tilted Black-Hole Accretion Disks

In recent work we presented the first results of global general relativistic magnetohydrodynamic (GRMHD) simulations of tilted (or misaligned) accretion disks around rotating black holes. The simulated tilted disks showed dramatic differences from comparable untilted disks, such as asymmetrical accretion onto the hole through opposing "plunging streams" and global precession of the disk powered by a torque provided by the black hole. However, those simulations used a traditional spherical-polar grid that was purposefully underresolved along the pole, which prevented us from assessing the behavior of any jets that may have been associated with the tilted disks. To address this shortcoming we have added a block-structured "cubed-sphere" grid option to the Cosmos++ GRMHD code, which will allow us to simultaneously resolve the disk and polar regions. Here we present our implementation of this grid and the results of a small suite of validation tests intended to demonstrate that the new grid performs as expected. The most important test in this work is a comparison of identical tilted disks, one evolved using our spherical-polar grid and the other with the cubed-sphere grid. We also demonstrate an interesting dependence of the early-time evolution of our disks on their orientation with respect to the grid alignment. This dependence arises from the differing treatment of current sheets within the disks, especially whether they are aligned with symmetry planes of the grid or not.Comment: 15 pages, 11 figures, submitted to Ap

Oscillation modes of relativistic slender tori

Accretion flows with pressure gradients permit the existence of standing waves which may be responsible for observed quasi-periodic oscillations (QPO's) in X-ray binaries. We present a comprehensive treatment of the linear modes of a hydrodynamic, non-self-gravitating, polytropic slender torus, with arbitrary specific angular momentum distribution, orbiting in an arbitrary axisymmetric spacetime with reflection symmetry. We discuss the physical nature of the modes, present general analytic expressions and illustrations for those which are low order, and show that they can be excited in numerical simulations of relativistic tori. The mode oscillation spectrum simplifies dramatically for near Keplerian angular momentum distributions, which appear to be generic in global simulations of the magnetorotational instability. We discuss our results in light of observations of high frequency QPO's, and point out the existence of a new pair of modes which can be in an approximate 3:2 ratio for arbitrary black hole spins and angular momentum distributions, provided the torus is radiation pressure dominated. This mode pair consists of the axisymmetric vertical epicyclic mode and the lowest order axisymmetric breathing mode.Comment: submitted to MNRA

Radiative Models of Sagittarius A* and M87 from Relativistic MHD Simulations

Ongoing millimeter VLBI observations with the Event Horizon Telescope allow unprecedented study of the innermost portion of black hole accretion flows. Interpreting the observations requires relativistic, time-dependent physical modeling. We discuss the comparison of radiative transfer calculations from general relativistic MHD simulations of Sagittarius A* and M87 with current and future mm-VLBI observations. This comparison allows estimates of the viewing geometry and physical conditions of the Sgr A* accretion flow. The viewing geometry for M87 is already constrained from observations of its large-scale jet, but, unlike Sgr A*, there is no consensus for its millimeter emission geometry or electron population. Despite this uncertainty, as long as the emission region is compact, robust predictions for the size of its jet launching region can be made. For both sources, the black hole shadow may be detected with future observations including ALMA and/or the LMT, which would constitute the first direct evidence for a black hole event horizon.Comment: 8 pages, 2 figures, submitted to the proceedings of AHAR 2011: The Central Kiloparse

Observing Lense-Thirring Precession in Tidal Disruption Flares

When a star is tidally disrupted by a supermassive black hole (SMBH), the streams of liberated gas form an accretion disk after their return to pericenter. We demonstrate that Lense-Thirring precession in the spacetime around a rotating SMBH can produce significant time evolution of the disk angular momentum vector, due to both the periodic precession of the disk and the nonperiodic, differential precession of the bound debris streams. Jet precession and periodic modulation of disk luminosity are possible consequences. The persistence of the jetted X-ray emission in the Swift J164449.3+573451 flare suggests that the jet axis was aligned with the spin axis of the SMBH during this event.Comment: 4 pages, 4 figures. Accepted for publication in Physical Review Letters. Minor changes made to match proof

Numerical Modeling of the Radio Nebula from the 2004 December 27 Giant Flare of SGR 1806-20

We use the relativistic hydrodynamics code Cosmos++ to model the evolution of the radio nebula triggered by the Dec. 27, 2004 giant flare event of soft gamma repeater 1806-20. We primarily focus on the rebrightening and centroid motion occurring subsequent to day 20 following the flare event. We model this period as a mildly relativistic (gamma ~ 1.07 - 1.67) jetted outflow expanding into the interstellar medium (ISM). We demonstrate that a jet with total energy ~ 10^46 ergs confined to a half opening angle ~ 20 degrees fits the key observables of this event, e.g. the flux lightcurve, emission map centroid position, and aspect ratio. In particular, we find excellent agreement with observations if the rebrightening is due to the jet, moving at 0.5c and inclined ~ 0 - 40 degrees toward the observer, colliding with a density discontinuity in the ISM at a radius of several 10^16 cm. We also find that a jet with a higher velocity, >~ 0.7c, and larger inclination, >~ 70 degrees, moving into a uniform ISM can fit the observations in general, but tends to miss the details of rebrightening. The latter, uniform ISM model predicts an ISM density more than 100 times lower than that of the former model, and thus suggests an independent test which might discriminate between the two. One of the strongest constraints of both models is that the data seems to require a non-uniform jet in order to be well fit

Impact of neutron star spin on Poynting-Robertson drag during a Type I X-ray burst

External irradiation of a neutron star (NS) accretion disc induces Poynting-Robertson (PR) drag, removing angular momentum and increasing the mass accretion rate. Recent simulations show PR drag significantly enhancing the mass accretion rate during Type I X-ray bursts, which could explain X-ray spectral features such as an increase in the persistent emission and a soft excess. However, prograde spin of the NS is expected to weaken PR drag, challenging its importance during bursts. Here, we study the effect of spin on PR drag during X-ray bursts. We run four simulations, with two assuming a non-spinning NS and two using a spin parameter of $a_*=0.2$, corresponding to a rotation frequency of 500 Hz. For each scenario, we simulate the disc evolution subject to an X-ray burst and compare it to the evolution found with no burst. PR drag drains the inner disc region during a burst, moving the inner disc radius outward by $\approx1.6$ km in the $a_*=0$ and by $\approx2.2$ km in the $a_*=0.2$ simulation. The burst enhances the mass accretion rate across the innermost stable circular orbit $\approx7.9$ times when the NS is not spinning and $\approx11.2$ times when it is spinning. The explanation for this seemingly contradictory result is that the disc is closer to the NS when $a_*=0.2$, and the resulting stronger irradiating flux offsets the weakening effect of spin on the PR drag. Hence, PR drag remains a viable explanation for the increased persistent emission and soft excess observed during X-ray bursts in spinning NS systems.Comment: 9 pages, 8 figures, accepted for publication in MNRA