Long-term observations have shown that black hole X-ray binaries exhibit strong, aperiodic variability on time-scales of a few milliseconds to seconds. The observed light curves display various characteristic features like a log-normal distribution of flux and a linear rms-flux relation, which indicate that the underlying variability process is stochastic in nature. It is also thought to be intrinsic to accretion. This variability has been modelled as inward propagating fluctuations of mass accretion rate, although the physical process driving the fluctuations remains puzzling. In this work, we analyse five exceptionally long duration general relativistic magnetohydrodynamic (GRMHD) simulations of optically thin, geometrically thick, black hole accretion flows to look for hints of propagating fluctuations in the simulation data. We find that the accretion profiles from these simulations do show evidence for inward propagating fluctuations below the viscous frequency by featuring strong radial coherence and positive time lags when comparing smaller to larger radii, although these time lags are generally shorter than the viscous time-scale and frequency independent. Our simulations also support the notion that the fluctuations in $\dot{M}$ build up in a multiplicative manner, as the simulations exhibit linear rms-mass flux relations, as well as log-normal distributions of their mass fluxes. When combining the mass fluxes from the simulations with an assumed emissivity profile, we additionally find broad agreement with observed power spectra and time lags, including a recovery of the frequency dependency of the time lags.

We perform a set of numerical experiments studying the interaction of Type I X-ray bursts with thin, Shakura-Sunyaev type accretion discs. Careful observations of X-ray spectra during such bursts have hinted at changes occurring in the inner regions of the disc. We now clearly demonstrate a number of key effects that take place simultaneously, including: evidence for weak, radiation-driven outflows along the surface of the disc; significant levels of Poynting-Robertson (PR) drag, leading to enhanced accretion; and prominent heating in the disc, which increases the height, while lowering the density and optical depth. The PR drag causes the inner edge of the disc to retreat from the neutron star surface toward larger radii and then recover on the timescale of the burst. We conclude that the rich interaction of an X-ray burst with the surrounding disc provides a novel way to study the physics of accretion onto compact objects.

We analyze the two brightest Chandra X-ray flares detected from Sagittarius A*, with peak luminosities more than 600 x and 245 x greater than the quiescent X-ray emission. The brightest flare has a distinctive double-peaked morphology --- it lasts 5.7 ksec ($\sim 2$ hours), with a rapid rise time of 1500 sec and a decay time of 2500 sec. The second flare lasts 3.4 ksec, with rise and decay times of 1700 sec and 1400 sec. These luminous flares are significantly harder than quiescence: the first has a power law spectral index $\Gamma = 2.06\pm 0.14$ and the second has $\Gamma = 2.03\pm 0.27$, compared to $\Gamma = 3.0\pm0.2$ for the quiescent accretion flow. These spectral indices (as well as the flare hardness ratios) are consistent with previously-detected Sgr A* flares, suggesting that bright and faint flares arise from similar physical processes. Leveraging the brightest flare's long duration and high signal-to-noise, we search for intraflare variability and detect excess X-ray power at a frequency of $\nu \approx 3$ mHz, but show that it is an instrumental artifact and not of astrophysical origin. We find no other evidence (at the 95% confidence level) for periodic or quasi-periodic variability in either flares' time series. We also search for non-periodic excess power but do not find compelling evidence in the power spectrum. Bright flares like these remain our most promising avenue for identifying Sgr A*'s short timescale variability in the X-ray, which may probe the characteristic size scale for the X-ray emission region.

Across black hole (BH) and neutron star (NS) low-mass X-ray binaries (LMXBs), there appears to be some correlation between certain high- and low-frequency quasi-periodic oscillations (QPOs). In a previous paper, we showed that for BH LMXBs, this could be explained by the simultaneous oscillation and precession of a hot, thick, torus-like corona. In the current work, we extend this idea to NS LMXBs by associating the horizontal branch oscillations (HBO) with precession and the upper-kiloHertz (ukHz) QPO with vertical epicyclic motion. For the Atoll source 4U 1608-52, the model can match many distinct, simultaneous observations of the HBO and ukHz QPO by varying the inner and outer radius of the torus, while maintaining fixed values for the mass (M_{NS}) and spin (a_*) of the neutron star. The best fit values are M_{NS} = 1.38 \pm 0.03 M_\odot and a_* = 0.325 \pm 0.005. By combining these constraints with the measured spin frequency, we are able to obtain an estimate for the moment of inertia of I_{NS} = 1.40 \pm 0.02 \times 10^{45} g cm^2, which places constraints on the equation of state. The model is unable to fit the lower-kHz QPO, but evidence suggests that QPO may be associated with the boundary layer between the accretion flow and the neutron star surface, which is not treated in this work.

We present results from general relativistic calculations of nuclear ignition in white dwarf stars triggered by near encounters with rotating intermediate mass black holes with different spin and alignment parameters. These encounters create thermonuclear environments characteristic of Type Ia supernovae capable of producing both calcium and iron group elements in arbitrary ratios, depending primarily on the proximity of the interaction which acts as a strong moderator of nucleosynthesis. We explore the effects of black hole spin and spin-orbital alignment on burn product synthesis to determine whether they might also be capable of moderating reactive flows. When normalized to equivalent impact penetration, accounting for frame dragging corrections, the influence of spin is weak, no more than 25% as measured by nuclear energy release and mass of burn products, even for near maximally rotating black holes. Stars on prograde trajectories approach closer to the black hole and produce significantly more unbound debris and iron group elements than is possible by encounters with nonrotating black holes or by retrograde orbits, at more than 50% mass conversion efficiency. The debris contains several radioisotopes, most notably Ni56, made in amounts that produce sub-luminous (but still observable) light curves compared to branch-normal SNe Ia.

We perform detailed variability analysis of two-dimensional viscous, radiation hydrodynamic numerical simulations of Shakura-Sunyaev thin disks around a stellar mass black hole. Disk models are initialized on both the gas-, as well as radiation-, pressure-dominated branches of the thermal equilibrium curve, with mass accretion rates spanning the range from $\dot{M} = 0.01 L_\mathrm{Edd}/c^2$ to $10 L_\mathrm{Edd}/c^2$. An analysis of temporal variations of the numerically simulated disk reveals multiple robust, coherent oscillations. Considering the local mass flux variability, we find an oscillation occurring at the maximum radial epicyclic frequency, $3.5\times 10^{-3}\,t_\mathrm{g}^{-1}$, a possible signature of a trapped fundamental ${\it g}$-mode. Although present in each of our simulated models, the trapped ${\it g}$-mode feature is most prominent in the gas-pressure-dominated case. The total pressure fluctuations in the disk suggest strong evidence for standing-wave ${\it p}$-modes, some trapped in the inner disk close to the ISCO, others present in the middle/outer parts of the disk. Knowing that the trapped ${\it g}$-mode frequency and maximum radial epicyclic frequency differ by only $0.01\%$ in the case of a non-rotating black hole, we simulated an additional initially gas-pressure-dominated disk with a dimensionless black hole spin parameter $a_* = 0.5$. The oscillation frequency in the spinning black hole case confirms that this oscillation is indeed a trapped ${\it g}$-mode. All the numerical models we report here also show a set of high frequency oscillations at the vertical epicyclic and breathing mode frequencies. The vertical oscillations show a 3:2 frequency ratio with oscillations occurring approximately at the radial epicyclic frequency, which could be of astrophysical importance in observed twin peak, high-frequency quasi-periodic oscillations.

Identifying the compact object in ultraluminous X-ray sources (ULXs) has to-date required detection of pulsations or a cyclotron resonance scattering feature (CRSF), indicating a magnetised neutron star. However, pulsations are observed to be transient and it is plausible that accretion onto the neutron star may have suppressed the surface magnetic field such that pulsations and CRSFs will be entirely absent. We may therefore lack direct means to identify neutron star systems whilst we presently lack an effective means by which to identify black hole ULXs. Here we present a possible method for separating the ULX population by assuming the X-ray, mHz quasi-periodic oscillations (QPOs) and day timescale periods/QPOs are associated with Lense-Thirring precession of the inflow and outflowing wind respectively. The precession timescales combined with the temperature of the soft X-ray component produce planes where the accretor mass enters as a free parameter. Depending on the properties of the wind, use of these planes may be robust to a range in the angular momentum (spin) and, for high accretion rates, essentially independent of the neutron star's surface dipole field strength. Our model also predicts the mHz QPO frequency and magnitude of the phase-lag imprinted due to propagation through the optically thick wind; in the case of NGC 5408 X-1 we subsequently infer a black hole mass and high spin. Finally, we note that observing secular QPO evolution over sufficient baselines may indicate a neutron star, as the precession responds to spin-up which is not readily observable for black hole primaries.

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.

The immense tidal forces of massive black holes can rip apart stars that come too close to them. As the resulting stellar debris spirals inwards, it heats up and emits x-rays when near the black hole. Here, we report the discovery of an exceptionally stable 131-second x-ray quasi-periodicity from a black hole after it disrupted a star. Using a black hole mass indicated from host galaxy scaling relations implies that, (1) this periodicity originates from very close to the black hole's event horizon, and (2) the black hole is rapidly spinning. Our findings suggest that other disruption events with similar highly sensitive observations likely also exhibit quasi-periodicities that encode information about the fundamental properties of their black holes.

We use two-dimensional, general relativistic, viscous, radiation hydrodynamic simulations to study the impact of a Type I X-ray burst on a hot and geometrically thick accretion disk surrounding an unmagnetized, non-rotating neutron star. The disk is initially consistent with a system in its low/hard spectral state, and is subject to a burst which rises to a peak luminosity of $10^{38}$ erg s$^{-1}$ in $2.05$ s. At the peak of the burst, the temperature of the disk has dropped by more than three orders of magnitude and its scale height has gone down by more than one order of magnitude. The simulations show that these effects predominantly happen due to Compton cooling of the hot plasma, and clearly illustrate the potential cooling effects of bursts on accretion disk coronae. In addition, we demonstrate the presence of Poynting-Robertson drag, though it only enhances the mass accretion rate onto the neutron star by a factor of $\sim 3$-$4$ compared to a simulation with no burst. Simulations such as these are important for building a general understanding of the response of an accretion disk to an intense X-ray impulse, which, in turn, will be crucial for deciphering burst spectra. Detailed analysis of such spectra offers the potential to measure neutron star radii, and hence constrain the neutron star equation of state, but only if the contributions coming from the impacted disk and its associated corona can be understood.

We study the oscillations of an axisymmetric, viscous, radiative, general relativistic hydrodynamical simulation of a geometrically thin disk around a non-rotating, $6.62\,M_\odot$ black hole. The numerical setup is initialized with a Novikov-Thorne, gas-pressure-dominated accretion disk, with an initial mass accretion rate of $\dot{m} = 0.01\,L_\mathrm{Edd}/c^2$ (where $L_\mathrm{Edd}$ is the Eddington luminosity and $c$ is the speed of light). Viscosity is treated with the $\alpha$-prescription. The simulation was evolved for about $1000$ Keplerian orbital periods at three Schwarzschild radii (ISCO radius). Power density spectra of the radial and vertical fluid velocity components, the total (gas $+$ radiation) midplane pressure, and the vertical component of radiative flux from the photosphere, all reveal strong power at the local breathing oscillation frequency. The first, second and third harmonics of the breathing oscillation are also clearly seen in the data. We quantify the properties of these oscillations by extracting eigenfunctions of the radial and vertical velocity components and total pressure. This confirms that these oscillations are associated with breathing motion.

We present results from general relativistic calculations of the tidal disruption of white dwarf stars from near encounters with intermediate mass black holes. We follow the evolution of 0.2 and $0.6 M_\odot$ stars on parabolic trajectories that approach $10^3$ - $10^4 M_\odot$ black holes as close as a few Schwarzschild radii at periapsis, paying particular attention to the effect tidal disruption has on thermonuclear reactions and the synthesis of intermediate to heavy ion elements. These encounters create diverse thermonuclear environments characteristic of Type I supernovae and capable of producing both intermediate and heavy mass elements in arbitrary ratios, depending on the strength (or proximity) of the interaction. Nuclear ignition is triggered in all of our calculations, even at weak tidal strengths $\beta \sim 2.6$ and large periapsis radius $R_P \sim 28$ Schwarzschild radii. A strong inverse correlation exists between the mass ratio of calcium to iron group elements and tidal strength, with $\beta \lesssim 5$ producing predominately calcium-rich debris. At these moderate to weak interactions, nucleosynthesis is not especially efficient, limiting the total mass and outflows of calcium group elements to $< 15$\% of available nuclear fuel. Iron group elements however continue to be produced in greater quantity and ratio with increasing tidal strength, peaking at $\sim 60$\% mass conversion efficiency in our closest encounter cases. These events generate short bursts of gravitational waves with characteristic frequencies 0.1-0.7 Hz and strain amplitudes $0.5\times10^{-22}$ - $3.5\times10^{-22}$ at 10 Mpc source distance.

We present results from two-dimensional, general relativistic, viscous, radiation hydrodynamic numerical simulations of Shakura-Sunyaev thin disks accreting onto stellar mass Schwarzschild black holes. We consider cases on both the gas- and radiation-pressure-dominated branches of the thermal equilibrium curve, with mass accretion rates spanning the range from $\dot{M} = 0.01 L_\mathrm{Edd}/c^2$ to $10 L_\mathrm{Edd}/c^2$. The simulations directly test the stability of this standard disk model on the different branches. We find clear evidence of thermal instability for all radiation-pressure-dominated disks, resulting universally in the vertical collapse of the disks, which in some cases then settle onto the stable, gas-pressure-dominated branch. Although these results are consistent with decades-old theoretical predictions, they appear to be in conflict with available observational data from black hole X-ray binaries. We also find evidence for a radiation-pressure-driven instability that breaks the unstable disks up into alternating rings of high and low surface density on a timescale comparable to the thermal collapse. Since radiation is included self-consistently in the simulations, we are able to calculate lightcurves and power density spectra (PDS). For the most part, we measure radiative efficiencies (ratio of luminosity to mass accretion rate) close to 6%, as expected for a non-rotating black hole. The PDS appear as broken power laws, with a break typically around 100 Hz. There is no evidence of significant excess power at any frequencies, i.e. no quasi-periodic oscillations are observed.

The presence of neutron stars in at least three ultraluminous X-ray sources is now firmly established and offers an unambiguous view of super-critical accretion. All three systems show long-timescale periods (60-80 days) in the X-rays and/or optical, two of which are known to be super-orbital in nature. Should the flow be classically super critical, i.e. the Eddington limit is reached locally in the disc (implying surface dipole fields that are sub-magnetar in strength), then the large scale-height flow can precess through the Lense-Thirring effect which could provide an explanation for the observed super-orbital periods. By connecting the details of the Lense-Thirring effect with the observed pulsar spin period, we are able to infer the moment-of-inertia and therefore equation-of-state of the neutron star without relying on the inclination of, or distance to the system. We apply our technique to the case of NGC 7793 P13 and demonstrate that stronger magnetic fields imply stiffer equations of state. We discuss the caveats and uncertainties, many of which can be addressed through forthcoming radiative magnetohydrodynamic (RMHD) simulations and their connection to observation.

We present results of three-dimensional, multi-physics simulations of an AGN jet colliding with an intergalactic cloud. The purpose of these simulations is to assess the degree of "positive feedback," i.e. jet-induced star formation, that results. We have specifically tailored our simulation parameters to facilitate comparison with recent observations of Minkowski's Object (M.O.), a stellar nursery located at the termination point of a radio jet coming from galaxy NGC 541. As shown in our simulations, such a collision triggers shocks which propagate around and through the cloud. These shocks condense the gas and under the right circumstances may trigger cooling instabilities, creating runaway increases in density, to the point that individual clumps can become Jeans unstable. Our simulations provide information about the expected star formation rate, total mass converted to \ion{H}{1}, H$_2$, and stars, and the relative velocity of the stars and gas. Our results confirm the possibility of jet-induced star formation, and agree well with the observations of M.O.

Monitoring of Sagittarius A* from X-ray to radio wavelengths has revealed structured variability --- including X-ray flares --- but it is challenging to establish correlations between them. Most studies have focused on variability in the X-ray and infrared, where variations are often simultaneous, and because long time series at sub-millimeter and radio wavelengths are limited. Previous work on sub-mm and radio variability hints at a lag between X-ray flares and their candidate sub-millimeter or radio counterparts, with the long wavelength data lagging the X-ray. However, there is only one published time lag between an X-ray flare and a possible radio counterpart. Here we report 9 contemporaneous X-ray and radio observations of Sgr A*. We detect significant radio variability peaking $\gtrsim$176 minutes after the brightest X-ray flare ever detected from Sgr A*. We also report other potentially associated X-ray and radio variability, with the radio peaks appearing $\lesssim$80 minutes after these weaker X-ray flares. Taken at face value, these results suggest that stronger X-ray flares lead to longer time lags in the radio. However, we also test the possibility that the variability at X-ray and radio wavelengths is not temporally correlated. We cross-correlate data from mismatched X-ray and radio epochs and obtain comparable correlations to the matched data. Hence, we find no overall statistical evidence that X-ray flares and radio variability are correlated, underscoring a need for more simultaneous, long duration X-ray--radio monitoring of Sgr A*.

We have extended Cosmos++, a multi-dimensional unstructured adaptive mesh code for solving the covariant Newtonian and general relativistic radiation magnetohydrodynamic (MHD) equations, to accommodate both discrete finite volume and arbitrarily high order finite element structures. The new finite element implementation, called CosmosDG, is based on a discontinuous Galerkin (DG) formulation, using both entropy-based artificial viscosity and slope limiting procedures for regularization of shocks. High order multi-stage forward Euler and strong stability preserving Runge-Kutta time integration options complement high order spatial discretization. We have also added flexibility in the code infrastructure allowing for both adaptive mesh and adaptive basis order refinement to be performed separately or simultaneously in a local (cell-by-cell) manner. We discuss in this report the DG formulation and present tests demonstrating the robustness, accuracy, and convergence of our numerical methods applied to special and general relativistic MHD, though we note an equivalent capability currently also exists in CosmosDG for Newtonian systems.

Strong magnetization in accretion discs could resolve a number of outstanding issues related to stability and state transitions in low-mass X-ray binaries. However, it is unclear how real discs become strongly magnetized and, even if they do, whether they can remain in such a state. In this paper, we address the latter issue through a pair of global disc simulations. Here, we only consider cases of initially purely toroidal magnetic fields contained entirely within a compact torus. We find that, over only a few tens of orbital periods, the magnetization of an initially strongly magnetized disc, $P_\mathrm{mag}/P_\mathrm{gas} \ge 10$, drops to $\lesssim 0.1$, similar to the steady-state value reached in initially weakly magnetized discs. This is consistent with recent shearing box simulations with initially strong toroidal fields, the robust conclusion being that strongly magnetized toroidal fields can not be locally self-sustaining. These results appear to leave net poloidal flux or extended radial fields as the only avenues for establishing strongly magnetized discs, ruling out the thermal collapse scenario.

We present results of a set of three-dimensional, general relativistic radiation magnetohydrodynamics simulations of thin accretion discs around a non-rotating black hole to test their thermal stability. We consider two cases, one that is initially radiation pressure dominated and expected to be thermally unstable and another that is initially gas-pressure dominated and expected to remain stable. Indeed, we find that cooling dominates over heating in the radiation pressure dominated model, causing the disc to collapse vertically on roughly the local cooling timescale. We also find that heating and cooling within the disc have a different dependence on the mid-plane pressure, a prerequisite of thermal instability. Comparison of our data with the relevant thin-disc thermal equilibrium curve suggests that our disc may be headed for the thermally stable, gas-pressure-dominated branch. However, because the disc collapses to the point that we are no longer able to resolve it, we had to terminate the simulation. On the other hand, the gas pressure dominated model, which was run for twice as long as the radiation pressure dominated one, remains stable, with heating and cooling roughly in balance. Finally, the radiation pressure dominated simulation shows some evidence of viscous instability. The strongest evidence is in plots of surface density, which show the disc breaking up into rings.

Motivated by recent studies showing an apparent correlation between the high-frequency quasi-periodic oscillations (QPOs) and the low-frequency, type-C QPO in low-mass, black hole X-ray binaries (LMXBs), we explore a model that explains all three QPOs in terms of an oscillating, precessing hot flow in the truncated-disk geometry. Our model favors attributing the two high-frequency QPOs, often occurring in a near 3:2 frequency ratio, to the breathing and vertical epicyclic frequency modes of the hot, thick flow, although we can not rule out the Keplerian and m=-1 radial epicyclic modes. In either case, the type-C QPO is attributed to precession. The correlation of the QPOs comes from the fact that all three frequencies are associated with the same geometrical structure. While the exact QPO frequencies are sensitive to the black hole mass and spin, their evolution over the course of an outburst is mainly tied to the truncation radius between the geometrically thin, optically thick disk and the inner, hot flow. We show that, in the case of the LMXB GRO J1655-40, this model can explain the one simultaneous observation of all three QPOs and that an extrapolation of the model appears to match lower frequency observations where only two of the three components are seen. Thus, this model may be able to unify multiple QPO observations using the properties of a single, simple, geometrical model.

We simulate an oscillating purely hydrodynamical torus with constant specific angular mo- mentum around a Schwarzschild black hole. The goal is to search for quasi-periodic oscil- lations (QPOs) in the light curve of the torus. The initial torus setup is subjected to radial, vertical and diagonal (combination of radial and vertical) velocity perturbations. The hydro- dynamical simulations are performed using the general relativistic magnetohydrodynamics code Cosmos++ and ray-traced using the GYOTO code. We found that a horizontal velocity perturbation triggers the radial and plus modes, while a vertical velocity perturbation trig- gers the vertical and X modes. The diagonal perturbation gives a combination of the modes triggered in the radial and vertical perturbations.

We report new observations with the Very Large Array, Atacama Large Millimeter Array, and Submillimeter Array at frequencies from 1.0 to 355 GHz of the Galactic Center black hole, Sagittarius A*. These observations were conducted between October 2012 and November 2014. While we see variability over the whole spectrum with an amplitude as large as a factor of 2 at millimeter wavelengths, we find no evidence for a change in the mean flux density or spectrum of Sgr A* that can be attributed to interaction with the G2 source. The absence of a bow shock at low frequencies is consistent with a cross-sectional area for G2 that is less than $2 \times 10^{29}$ cm$^2$. This result fits with several model predictions including a magnetically arrested cloud, a pressure-confined stellar wind, and a stellar photosphere of a binary merger. There is no evidence for enhanced accretion onto the black hole driving greater jet and/or accretion flow emission. Finally, we measure the millimeter wavelength spectral index of Sgr A* to be flat; combined with previous measurements, this suggests that there is no spectral break between 230 and 690 GHz. The emission region is thus likely in a transition between optically thick and thin at these frequencies and requires a mix of lepton distributions with varying temperatures consistent with stratification.

We investigate the local stability of strongly magnetized relativistic tori orbiting Kerr black holes, for the case of a purely toroidal magnetic field topology. Our approach encompasses both solving the full dispersion relation numerically and solving a simplified analytic treatment. Both approaches indicate that such tori are subject to an unstable non-axisymmetric magnetorotational mode, regardless of such background properties as the gas-to-magnetic pressure ratio and specific angular momentum distribution. We demonstrate that our modal analysis matches well with what is seen in global, three-dimensional, general relativistic magnetohydrodynamic simulations.

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.

In this paper we introduce the first results that use data extracted directly from numerical simulations as inputs to the analytic twisted disk model of Zhuravlev & Ivanov. In both numerical and analytic approaches, fully relativistic models of tilted and twisted disks having a moderate effective viscosity around a slowly rotating Kerr black hole are considered. Qualitatively, the analytic model demonstrates the same dynamics as the simulations, although with some quantitative offset. Namely, the GRMHD simulations generally give smaller variations of tilt and twist across the disk. When the black hole and the disk rotate in the same sense, the simulated tilted disk and analytic model show no sign of Bardeen-Petterson alignment, even in the innermost parts of the disk, where the characteristic time for relaxation to a quasi-stationary configuration is of the same order as the computation time. In the opposite case, when the direction of the disk's rotation is opposite to that of black hole, a partial alignment is observed, in agreement with previous theoretical estimates. Thus, both fully numerical and analytic schemes demonstrate that the Bardeen-Petterson effect may not be possible for the case of prograde rotation provided that disk's effective viscosity is sufficiently small. This may have implications in modeling of different astrophysical phenomena, such as disk spectra and jet orientation.

This paper presents our latest numerical simulations of accretion disks that are misaligned with respect to the rotation axis of a Kerr black hole. In this work we use a new, fully conservative version of the Cosmos++ general relativistic magnetohydrodynamics (GRMHD) code, coupled with an ad hoc cooling function designed to control the thickness of the disk. Together these allow us to simulate the thinnest tilted accretion disks ever using a GRMHD code. In this way, we are able to probe the regime where the dimensionless stress and scale height of the disk become comparable. We present results for both prograde and retrograde cases. The simulated prograde tilted disk shows no sign of Bardeen-Petterson alignment even in the innermost parts of the disk. The simulated retrograde tilted disk, however, does show modest alignment. The implication of these results is that the parameter space associated with Bardeen-Petterson alignment for prograde disks may be rather small, only including very thin disks. Unlike our previous work, we find no evidence for standing shocks in our simulated tilted disks. We ascribe this to the combination of small black hole spin, small tilt angle, and small disk scale height in these simulations. We also add to the growing body of literature pointing out that the turbulence driven by the magnetorotational instability in global simulations of accretion disks is not isotropic. Finally, we provide a comparison between our moderately thin, untilted reference simulation and other numerical simulations of thin disks in the literature.

In this paper we report on recent upgrades to our general relativistic radiation magnetohydrodynamics code, Cosmos++, including the development of a new primitive inversion scheme and a hybrid implicit-explicit solver with a more general closure relation for the radiation equations. The new hybrid solver helps stabilize the treatment of the radiation source terms, while the new closure allows for a much broader range of optical depths to be considered. These changes allow us to expand by orders of magnitude the range of temperatures, opacities, and mass accretion rates, and move a step closer toward our goal of performing global simulations of radiation-pressure-dominated black hole accretion disks. In this work we test and validate the new method against an array of problems. We also demonstrate its ability to handle super-Eddington, quasi-spherical accretion. Even with just a single proof-of-principle simulation, we already see tantalizing hints of the interesting phenomenology associated with the coupling of radiation and gas in super-Eddington accretion flows.

We show that the low frequency QPO seen in the power density spectra of black hole binaries (and neutron stars) can be explained by Lense-Thirring precession. This has been proposed many times in the past, and simple, single radius models can qualitatively match the observed increase in QPO frequency by decreasing a characteristic radius, as predicted by the truncated disc models. However, this also predicts that the frequency is strongly dependent on spin, and gives a maximum frequency at the last stable orbit which is generally much higher than the remarkably constant maximum frequency at ~10Hz observed in all black hole binaries. The key aspect of our model which makes it match these observations is that the precession is of a radially extended region of the hot inner flow. The outer radius is set by the truncation radius of the disc as above, but the inner radius lies well outside of the last stable orbit at the point where numerical simulations show that the density drops off sharply for a misaligned flow. Physically motivated analytic estimates for this inner radius show that it increases with a_*, decreasing the expected frequency in a way which almost completely cancels the expected increase with spin. This ties the maximum predicted frequency to around 10Hz irrespective of a_*, as observed. This is the first QPO model which explains both frequencies and spectrum in the context of a well established geometry for the accretion flow.

Using three-dimensional, moving-mesh simulations, we investigate the future evolution of the recently discovered gas cloud G2 traveling through the galactic center. From our simulations we expect an average feeding rate onto Sgr A* in the range of $(5-19) \times 10^{-8} M_\odot\mathrm{~yr}^{-1}$ beginning in 2014. This accretion varies by less than a factor of three on timescales of about 1 month, and shows no more than a factor of 10 difference between the maximum and minimum observed rates within any given model. These rates are comparable to the current estimated accretion rate in the immediate vicinity of Sgr A*, although they represent only a small (<10%) increase over the current expected feeding rate at the effective inner boundary of our simulations $(r_\mathrm{acc} = 750 R_S \sim 10^{15} \mathrm{cm})$. We also explore multiple possible equations of state to describe the gas. In examining the Br-$\gamma$ light curves produced from our simulations, we find that all of our isothermal models predict significant (factor of 10) enhancements in the luminosity of G2 as it approaches pericenter, in conflict with observations. Models that instead allow the cloud to heat as it is compressed do better at matching observations.

Using three-dimensional, moving-mesh simulations, we investigate the future evolution of the recently discovered gas cloud G2 traveling through the galactic center. We consider the case of a spherical cloud initially in pressure equilibrium with the background. Our suite of simulations explores the following parameters: the equation of state, radial profiles of the background gas, and start times for the evolution. Our primary focus is on how the fate of this cloud will affect the future activity of Sgr A*. From our simulations we expect an average feeding rate in the range of $5-19 \times 10^{-8}$ solar masses per year beginning in 2013 and lasting for at least 7 years (our simulations stop in year 2020). The accretion varies by less than a factor of three on timescales <1 month, and shows no more than a factor of 10 difference between the maximum and minimum. These rates are comparable to the current estimated accretion rate in the immediate vicinity of Sgr A*, although they represent only a small (<5%) increase over the current expected feeding rate at the effective inner boundary of our simulations. Therefore, the break up of cloud G2 may have only a minimal effect on the brightness and variability of Sgr A* over the next decade. This is because current models of the galactic center predict that most of the gas will be caught up in outflows. However, if the accreted G2 material can remain cold, it may not mix well with the hot, diffuse background gas, and instead accrete efficiently onto Sgr A*. Further observations of G2 will give us an unprecedented opportunity to test this idea. The break up of the cloud itself may also be observable. By tracking the amount of cloud energy that is dissipated during our simulations, we are able to get a rough estimate of the luminosity associated with its tidal disruption; we find values of a few $10^{36}$ erg/s.

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.

Simulations of hot, pressure supported, tilted black hole accretion flows, in which the angular momentum of the flow is misaligned with the black hole spin axis, can exhibit two non-axisymmetric shock structures in the inner regions of the flow. We analyze the strength and significance of these shock structures in simulations with tilt angles of 10 and 15 degrees. By integrating fluid trajectories in the simulations through the shocks, and tracking the variations of fluid quantities along these trajectories, we show that these shocks are strong, with substantial compression ratios, in contrast to earlier claims. However, they are only moderately relativistic. We also show that the two density enhancements resembling flow streams in their shape are in fact merely post-shock compressions, as fluid trajectories cut across, rather than flow along, them. The dissipation associated with the shocks is a substantial fraction ($\simeq3-12$ percent) of the rest mass energy advected into the hole, and therefore comparable to the dissipation expected from turbulence. The shocks should therefore make order unity changes in the observed properties of black hole accretion flows that are tilted.

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.

As the title suggests, the purpose of this chapter is to review the current status of numerical simulations of black hole accretion disks. This chapter focuses exclusively on global simulations of the accretion process within a few tens of gravitational radii of the black hole. Most of the simulations discussed are performed using general relativistic magnetohydrodynamic (MHD) schemes, although some mention is made of Newtonian radiation MHD simulations and smoothed particle hydrodynamics. The goal is to convey some of the exciting work that has been going on in the past few years and provide some speculation on future directions.

This review covers the main aspects of black hole accretion disk theory. We begin with the view that one of the main goals of the theory is to better understand the nature of black holes themselves. In this light we discuss how accretion disks might reveal some of the unique signatures of strong gravity: the event horizon, the innermost stable circular orbit, and the ergosphere. We then review, from a first-principles perspective, the physical processes at play in accretion disks. This leads us to the four primary accretion disk models that we review: Polish doughnuts (thick disks), Shakura-Sunyaev (thin) disks, slim disks, and advection-dominated accretion flows (ADAFs). After presenting the models we discuss issues of stability, oscillations, and jets. Following our review of the analytic work, we take a parallel approach in reviewing numerical studies of black hole accretion disks. We finish with a few select applications that highlight particular astrophysical applications: measurements of black hole mass and spin, black hole vs. neutron star accretion disks, black hole accretion disk spectral states, and quasi-periodic oscillations (QPOs).

We study the spatial and temporal behavior of fluid in fully three-dimensional, general relativistic, magnetohydrodynamical simulations of both tilted and untilted black hole accretion flows. We uncover characteristically greater variability in tilted simulations at frequencies similar to those predicted by the formalism of trapped modes, but ultimately conclude that its spatial structure is inconsistent with a modal interpretation. We find instead that previously identified, transient, over-dense clumps orbiting on roughly Keplerian trajectories appear generically in our global simulations, independent of tilt. Associated with these fluctuations are acoustic spiral waves interior to the orbits of the clumps. We show that the two nonaxisymmetric standing shock structures that exist in the inner regions of these tilted flows effectively amplify the variability caused by these spiral waves to markedly higher levels than in untilted flows, which lack standing shocks. Our identification of clumps, spirals, and spiral-shock interactions in these fully general relativistic, magnetohydrodynamical simulations suggests that these features may be important dynamical elements in models which incorporate tilt as a way to explain the observed variability in black hole accretion flows.

We present the first spectral energy distributions produced self-consistently by 2.5D general relativistic magneto-hydrodynamical (GRMHD) numerical simulations, where radiative cooling is included in the dynamical calculation. As a case study, we focus on the accretion flow around the supermassive black hole in the Galactic Centre, Sagittarius A* (Sgr A*), which has the best constrained physical parameters. We compare the simulated spectra to the observational data of Sgr A* and explore the parameter space of our model to determine the effect of changing the initial magnetic field configuration, ion to electron temperature ratio T_i/T_e and the target accretion rate. We find the best description of the data for a mass accretion rate of ~ 1e-9 Msun/yr, and rapid spin (0.7 < a_* < 0.9). The submillimeter peak flux seems largely independent of initial conditions, while the higher energies can be very sensitive to the initial magnetic field configuration. Finally, we also discuss flaring features observed in some simulations, that may be due to artifacts of the 2D configuration.

We present general relativistic magnetohydrodynamic (GRMHD) numerical simulations of the accretion flow around the supermassive black hole in the Galactic centre, Sagittarius A* (Sgr A*). The simulations include for the first time radiative cooling processes (synchrotron, bremsstrahlung, and inverse Compton) self-consistently in the dynamics, allowing us to test the common simplification of ignoring all cooling losses in the modeling of Sgr A*. We confirm that for Sgr A*, neglecting the cooling losses is a reasonable approximation if the Galactic centre is accreting below ~10^{-8} Msun/yr i.e. Mdot < 10^{-7} Mdot_Edd. But above this limit, we show that radiative losses should be taken into account as significant differences appear in the dynamics and the resulting spectra when comparing simulations with and without cooling. This limit implies that most nearby low-luminosity active galactic nuclei are in the regime where cooling should be taken into account. We further make a parameter study of axisymmetric gas accretion around the supermassive black hole at the Galactic centre. This approach allows us to investigate the physics of gas accretion in general, while confronting our results with the well studied and observed source, Sgr A*, as a test case. We confirm that the nature of the accretion flow and outflow is strongly dependent on the initial geometry of the magnetic field. For example, we find it difficult, even with very high spins, to generate powerful outflows from discs threaded with multiple, separate poloidal field loops.

It is now well established that changes in the X-ray spectral state of black hole low-mass X-ray binaries are correlated with changes in the radio properties of those systems. Assuming radio power is a proxy for jet power, we can say that the jet is continuously present in the hard state and undetectable (and therefore weaker) in the soft state. Since the different accretion states are also generally assumed to be associated with different disc geometries -- the hard state with a hot, thick flow, and the soft state with a cold, thin disc -- we investigate the possibility that these two phenomena are linked; i.e., that the difference in disc geometry is the cause of the difference in observed jet power. We do this by comparing various measures of jet power in numerical simulations of accretion discs of differing temperatures and thicknesses. We perform these simulations using the general relativistic magnetohydrodynamic code Cosmos++ and a newly added cooling function, which allows us to regulate the disc scale height H/r at different radii. We find no apparent correlation between the disc scale height and jet power whenever we normalize the latter by the mass accretion history of each simulation. We attribute this result to the role that the "corona" plays in confining and accelerating the jet (our corona may also be considered a failed MHD "wind"). The properties of the corona do not vary significantly from one simulation to another, even though the scale heights of the discs vary by up to a factor of four. If this holds true in nature, then it suggests that the correlation between spectral state and jet power must be attributable to some other property, possibly the topology of the magnetic field. Alternatively, it could be that the corona disappears altogether in the soft state, which would be consistent with observations, but has so far not been seen in simulations.

Modeling the radiation generated by accreting matter is an important step towards realistic simulations of black hole accretion disks, especially at high accretion rates. To this end, we have recently added radiation transport to the existing general relativistic magnetohydrodynamic code, Cosmos++. However, before attempting to model radiative accretion disks, we have tested the new code using a series of shock tube and Bondi (spherical inflow) problems. The four radiative shock tube tests, first presented by Farris et al. (2008), have known analytic solutions, allowing us to calculate errors and convergence rates for our code. The Bondi problem only has an analytic solution when radiative processes are ignored, but is pertinent because it is closer to the physics we ultimately want to study. In our simulations, we include Thomson scattering and thermal bremsstrahlung in the opacity, focusing exclusively on the super-Eddington regime. Unlike accretion onto bodies with solid surfaces, super-Eddington accretion onto black holes does not produce super-Eddington luminosity. In our examples, despite accreting at up to 300 times the Eddington rate, our measured luminosity is always several orders of magnitude below Eddington.

High-resolution, multi-wavelength, and time-domain observations of the Galactic centre black hole candidate, Sgr A*, allow for a direct test of contemporary accretion theory. To date, all models have assumed alignment between the accretion disc and black hole angular momentum axes, but this is unjustified for geometrically thick accretion flows like that onto Sgr A*. Instead, we calculate images and spectra from a set of simulations of accretion flows misaligned ('tilted') by 15 degrees from the black hole spin axis and compare them with millimetre (mm) to near-infrared (NIR) observations. Non-axisymmetric standing shocks from eccentric fluid orbits dominate the emission, leading to a wide range of possible image morphologies. These effects invalidate previous parameter estimates from model fitting, including estimates of the dimensionless black hole spin, except possibly at low values of spin or tilt. At 1.3mm, the images have crescent morphologies, and the black hole shadow may still be accessible to future mm-VLBI observations. Shock heating leads to high energy electrons (T > 10^12 K), which can naturally produce the observed NIR flux, spectral index, and rapid variability ('flaring'). This NIR emission is uncorrelated with that in the mm, which also agrees with observations. These are the first models to self-consistently explain the time-variable mm to NIR emission of Sgr A*. Predictions of the model include significant structural changes observable with mm-VLBI on both the dynamical (hour) and Lense-Thirring precession (day-year) timescales; and ~30-50 microarcsecond changes in centroid position from extreme gravitational lensing events during NIR flares, detectable with the future VLT instrument GRAVITY. If the observed NIR emission is caused by shock heating in a tilted accretion disc, then the Galactic centre black hole has a positive, non-zero spin parameter (a > 0).

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.

Geometrically thick accretion flows may be present in black hole X-ray binaries observed in the low/hard state and in low-luminosity active galactic nuclei. Unlike in geometrically thin disks, the angular momentum axis in these sources is not expected to align with the black hole spin axis. We compute images from three-dimensional general relativistic magnetohydrodynamic simulations of misaligned (tilted) accretion flows using relativistic radiative transfer, and compare the estimated locations of the radiation edge with expectations from their aligned (untilted) counterparts. The radiation edge in the tilted simulations is independent of black hole spin for a tilt of 15 degrees, in stark contrast to the results for untilted simulations, which agree with the monotonic dependence on spin expected from thin accretion disk theory. Synthetic emission line profiles from the tilted simulations depend strongly on the observer's azimuth, and exhibit unique features such as broad "blue wings." Coupled with precession, the azimuthal variation could generate time fluctuations in observed emission lines, which would be a clear "signature" of a tilted accretion flow. Finally, we evaluate the possibility that the observed low- and high-frequency quasi-periodic oscillations (QPOs) from black hole binaries could be produced by misaligned accretion flows. Although low-frequency QPOs from precessing, tilted disks remains a viable option, we find little evidence for significant power in our light curves in the frequency range of high-frequency QPOs.

Recent high resolution observations of the Galactic center black hole allow for direct comparison with accretion disk simulations. We compare two-temperature synchrotron emission models from three dimensional, general relativistic magnetohydrodynamic simulations to millimeter observations of Sgr A*. Fits to very long baseline interferometry and spectral index measurements disfavor the monochromatic face-on black hole shadow models from our previous work. Inclination angles \le 20 degrees are ruled out to 3 \sigma. We estimate the inclination and position angles of the black hole, as well as the electron temperature of the accretion flow and the accretion rate, to be i=50+35-15 degrees, \xi=-23+97-22 degrees, T_e=(5.4 +/- 3.0)x10^10 K and Mdot=(5+15-2)x10^-9 M_sun / yr respectively, with 90% confidence. The black hole shadow is unobscured in all best fit models, and may be detected by observations on baselines between Chile and California, Arizona or Mexico at 1.3mm or .87mm either through direct sampling of the visibility amplitude or using closure phase information. Millimeter flaring behavior consistent with the observations is present in all viable models, and is caused by magnetic turbulence in the inner radii of the accretion flow. The variability at optically thin frequencies is strongly correlated with that in the accretion rate. The simulations provide a universal picture of the 1.3mm emission region as a small region near the midplane in the inner radii of the accretion flow, which is roughly isothermal and has \nu/\nu_c ~ 1-20, where \nu_c is the critical frequency for thermal synchrotron emission.

This paper presents a numerical study over a wide parameter space of the likelihood of the dynamical bar-mode instability in differentially rotating magnetized neutron stars. The innovative aspect of this study is the incorporation of magnetic fields in such a context, which have thus far been neglected in the purely hydrodynamical simulations available in the literature. The investigation uses the Cosmos++ code which allows us to perform three dimensional simulations on a cylindrical grid at high resolution. A sample of Newtonian magneto-hydrodynamical simulations starting from a set of models previously analyzed by other authors without magnetic fields has been performed, providing estimates of the effects of magnetic fields on the dynamical bar-mode deformation of rotating neutron stars. Overall, our results suggest that the effect of magnetic fields are not likely to be very significant in realistic configurations. Only in the most extreme cases are the magnetic fields able to suppress growth of the bar mode.

One of the primary means of determining the spin of an astrophysical black hole is by actually measuring the inner radius of a surrounding accretion disk and using that to infer the spin. By comparing a number of different estimates of the inner radius from simulations of tilted accretion disks with differing black-hole spins, we show that such a procedure can give quite wrong answers. Over the range 0 <= a/M <= 0.9, we find that, for moderately thick disks (H/r ~ 0.2) with modest tilt (15 degrees), the inner radius is nearly independent of spin. This result is likely dependent on tilt, such that for larger tilts, it may even be that the inner radius would increase with increasing spin. In the opposite limit, we confirm through numerical simulations of untilted disks that, in the limit of zero tilt, the inner radius recovers approximately the expected dependence on spin.

We analyze the time dependence of fluid variables in general relativistic, magnetohydrodynamic simulations of accretion flows onto a black hole with dimensionless spin parameter a/M=0.9. We consider both the case where the angular momentum of the accretion material is aligned with the black hole spin axis (an untilted flow) and where it is misaligned by 15 degrees (a tilted flow). In comparison to the untilted simulation, the tilted simulation exhibits a clear excess of inertial variability, that is, variability at frequencies below the local radial epicyclic frequency. We further study the radial structure of this inertial-like power by focusing on a radially extended band at 118 (M/10Msol)^-1 Hz found in each of the three analyzed fluid variables. The three dimensional density structure at this frequency suggests that the power is a composite oscillation whose dominant components are an over dense clump corotating with the background flow, a low order inertial wave, and a low order inertial-acoustic wave. Our results provide preliminary confirmation of earlier suggestions that disk tilt can be an important excitation mechanism for inertial waves.

The recent VLBI observation of the Galactic center black hole candidate Sgr A* at 1.3mm shows source structure on event-horizon scales. This detection enables a direct comparison of the emission region with models of the accretion flow onto the black hole. We present the first results from time-dependent radiative transfer of general relativistic MHD simulation data, and compare simulated synchrotron images at black hole spin a=0.9 with the VLBI measurements. After tuning the accretion rate to match the millimeter flux, we find excellent agreement between predicted and observed visibilities, even when viewed face-on (i < 30 degrees). VLBI measurements on 2000-3000km baselines should constrain the inclination. The data constrain the accretion rate to be (1.0-2.3)x10^-9 M_sun / yr with 99% confidence, consistent with but independent of prior estimates derived from spectroscopic and polarimetric measurements. Finally, we compute light curves, which show that magnetic turbulence can directly produce flaring events with .5 hour rise times, 2-3.5 hour durations and 40-50% flux modulation, in agreement with observations of Sgr A* at millimeter wavelengths.

We construct a new family of analytic models of black hole accretion disks in dynamical equilibria. Our construction is based on assuming distributions of angular momentum and entropy. For a particular choice of the distribution of angular momentum, we calculate the shapes of equipressure surfaces. The equipressure surfaces we find are similar to those in thick, slim and thin disks, and to those in ADAFs.

Correlations among pulse properties in the prompt emission of long GRBs can potentially be used as cosmological distance indicators to estimate redshifts of GRBs to which these pulses belong. We demonstrate application of this technique to a sample of GRBs for which redshifts are not known. We also study the scatter of predicted redshifts of pulses found within individual bursts. We explore the characteristics of this scatter in hopes of identifying systematic corrections and/or pulse subsets that can be used to increase the technique's reliability.

(Abridged) We present one of the first physically-motivated two-dimensional general relativistic magnetohydrodynamic (GRMHD) numerical simulations of a radiatively-cooled black-hole accretion disk. The fiducial simulation combines a total-energy-conserving formulation with a radiative cooling function, which includes bremsstrahlung, synchrotron, and Compton effects. By comparison with other simulations we show that in optically thin advection-dominated accretion flows, radiative cooling can significantly affect the structure, without necessarily leading to an optically thick, geometrically thin accretion disk. We further compare the results of our radiatively-cooled simulation to the predictions of a previously developed analytic model for such flows. For the very low stress parameter and accretion rate found in our simulated disk, we closely match a state called the "transition" solution between an outer advection-dominated accretion flow and what would be a magnetically-dominated accretion flow (MDAF) in the interior. The qualitative and quantitative agreement between the numerical and analytic models is quite good, with only a few well-understood exceptions. According to the analytic model then, at significantly higher stress or accretion, we would expect a full MDAF to form. The collection of simulations in this work also provide important data for interpreting other numerical results in the literature, as they span the most common treatments of thermodynamics, including simulations evolving: 1) the internal energy only; 2) the internal energy plus an explicit cooling function; 3) the total energy without cooling; and 4) total energy including cooling. We find that the total energy formulation is a necessary prerequisite for proper treatment of radiative cooling in MRI accretion flows.

In this talk I review the current status of jet formation in direct numerical simulations of black-hole accretion disks and magnetospheres. I address the following critical questions: What constitutes the jet? What is the launching mechanism? Where is the launching point of the jet? What is the Lorentz factor? What is the opening angle? How is the jet collimated? Just as importantly, I also discuss how dependent the answers to the above questions are on factors such as the initial conditions of the simulation. I end by discussing possible future directions for this research.

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.

This work presents a detailed analysis of the overall flow structure and unique features of the inner region of the tilted disk simulations described in Fragile et al. (2007). The primary new feature identified in the main disk body is a latitude-dependent radial epicyclic motion driven by pressure gradients attributable to the gravitomagnetic warping of the disk. The induced motion of the gas is coherent over the scale of the entire disk and is fast enough that it could be observable in features such as relativistic iron lines. The eccentricity of the associated fluid element trajectories increases with decreasing radius, leading to a crowding of orbit trajectories near their apocenters. This results in a local density enhancement akin to a compression. These compressions are sufficiently strong to produce a pair of weak shocks in the vicinity of the black hole. These shocks are roughly aligned with the line-of-nodes between the black-hole symmetry plane and disk midplane, with one shock above the line-of-nodes on one side of the black hole and the other below on the opposite side. These shocks enhance angular momentum transport and energy dissipation near the hole, forcing some material to plunge toward the black hole from well outside the innermost stable circular orbit. A new, extended simulation, which was evolved for more than a full disk precession period, allows us to confirm that these shocks and the previously identified ``plunging streams'' precess with the disk in such a way as to remain aligned relative to the line-of-nodes, as expected based on our physical understanding of these phenomena. Such a precessing structure would likely present a strong quasi-periodic signal.

We derive a new peak lag vs. peak luminosity relation in gamma-ray burst (GRB) pulses. We demonstrate conclusively that GRB spectral lags are pulse rather than burst properties and show how the lag vs. luminosity relation determined from CCF measurements of burst properties is essentially just a rough measure of this newly derived relation for individual pulses. We further show that most GRB pulses have correlated properties: short-lag pulses have shorter durations, are more luminous, and are harder within a burst than long-lag pulses. We also uncover a new pulse duration vs. pulse peak luminosity relation, and indicate that long-lag pulses often precede short-lag pulses. Although most pulse behaviors are supportive of internal shocks (including long-lag pulses), we identify some pulse shapes that could result from external shocks.

This paper presents a continuation of our efforts to numerically study accretion disks that are misaligned (tilted) with respect to the rotation axis of a Kerr black hole. Here we present results of a global numerical simulation which fully incorporates the effects of the black hole spacetime as well as magnetorotational turbulence that is the primary source of angular momentum transport in the flow. This simulation shows dramatic differences from comparable simulations of untilted disks. Accretion onto the hole occurs predominantly through two opposing plunging streams that start from high latitudes with respect to both the black-hole and disk midplanes. This is due to the aspherical nature of the gravitational spacetime around the rotating black hole. These plunging streams start from a larger radius than would be expected for an untilted disk. In this regard the tilted black hole effectively acts like an untilted black hole of lesser spin. Throughout the duration of the simulation, the main body of the disk remains tilted with respect to the symmetry plane of the black hole; thus there is no indication of a Bardeen-Petterson effect in the disk at large. The torque of the black hole instead principally causes a global precession of the main disk body. In this simulation the precession has a frequency of $3 (M_\odot/M)$ Hz, a value consistent with many observed low-frequency quasi-periodic oscillations. However, this value is strongly dependent on the size of the disk, so this frequency may be expected to vary over a large range.

We posit that accreting compact objects, including stellar mass black holes and neutron stars as well as supermassive black holes, may undergo extended periods of accretion during which the angular momentum of the disk at large scales is misaligned with that of the compact object. In such a scenario, Lense-Thirring precession caused by the rotating compact object can dramatically affect the disk. In this presentation we describe results from a three-dimensional relativistic magnetohydrodynamic simulation of an MRI turbulent disk accreting onto a tilted rapidly rotating black hole. For this case, the disk does not achieve the commonly described Bardeen-Petterson configuration; rather, it remains nearly planar, undergoing a slow global precession. Accretion from the disk onto the hole occurs predominantly through two opposing plunging streams that start from high latitudes with respect to both the black-hole and disk midplanes. This is a consequence of the non-sphericity of the gravitational spacetime of the black hole.

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.

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.

A new code and methodology are introduced for solving the general relativistic magnetohydrodynamic (GRMHD) equations in fixed background spacetimes using time-explicit, finite-volume discretization. The code has options for solving the GRMHD equations using traditional artificial-viscosity (AV) or non-oscillatory central difference (NOCD) methods, or a new extended AV (eAV) scheme using artificial-viscosity together with a dual energy-flux-conserving formulation. The dual energy approach allows for accurate modeling of highly relativistic flows at boost factors well beyond what has been achieved to date by standard artificial viscosity methods. It provides the benefit of Godunov methods in capturing high Lorentz boosted flows but without complicated Riemann solvers, and the advantages of traditional artificial viscosity methods in their speed and flexibility. Additionally, the GRMHD equations are solved on an unstructured grid that supports local adaptive mesh refinement using a fully threaded oct-tree (in three dimensions) network to traverse the grid hierarchy across levels and immediate neighbors. A number of tests are presented to demonstrate robustness of the numerical algorithms and adaptive mesh framework over a wide spectrum of problems, boosts, and astrophysical applications, including relativistic shock tubes, shock collisions, magnetosonic shocks, Alfven wave propagation, blast waves, magnetized Bondi flow, and the magneto-rotational instability in Kerr black hole spacetimes.

In this work we show that Bardeen-Petterson accretion disks can exhibit unique, detectable features in relativistically broadened emission line profiles. Some of the unique characteristics include inverted line profiles with sharper red horns and softer blue horns and even profiles with more than 2 horns from a single rest-frame line. We demonstrate these points by constructing a series of synthetic line profiles using simple two-component disk models. We find that the resultant profiles are very sensitive to the two key parameters one would like to constrain, namely the Bardeen-Petterson transition radius r_{BP} and the relative tilt \beta between the two disk components over a range of likely values [10 < r_{BP}/(GM/c^2) < 40 ; 15deg < \beta < 45deg]. We use our findings to show that some of the ``extra'' line features observed in the spectrum of the Seyfert-I galaxy MCG--6-30-15 may be attributable to a Bardeen-Petterson disk structure. Similarly, we apply our findings to two likely Bardeen-Petterson candidate Galactic black holes - GRO J1655-40 and XTE J1550-564. We provide synthetic line profiles of these systems using observationally constrained sets of parameters. Although we do not formally fit the data for any of these systems, we confirm that our synthetic spectra are consistent with current observations.

Motivated by recent interesting work on p-mode oscillations in axisymmetric hydrodynamic black-hole tori by Rezzolla, Zanotti, and collaborators, I explore the robustness of these oscillations by means of two and three-dimensional relativistic hydrodynamic and MHD simulations. The primary purpose of this investigation is to determine how the amplitudes of these oscillations are affected by the presence of known instabilities of black-hole tori, including the Papaloizou-Pringle instability (PPI) and the magneto-rotational instability (MRI). Both instabilities drive accretion at rates above those considered in Rezzolla's work. The increased accretion can allow wave energy to leak out of the torus into the hole. Furthermore, with the MRI, the presence of turbulence, which is absent in the hydrodynamic simulations, can lead to turbulent damping (or excitation) of modes. The current numerical results are preliminary, but suggest that the PPI and MRI both significantly damp acoustic oscillations in tori.

We present results from fully general relativistic three-dimensional numerical studies of thick-disk accretion onto a rapidly-rotating (Kerr) black hole with a spin axis that is tilted (not aligned) with the angular momentum vector of the disk. We initialize the problem with the solution for an aligned, constant angular momentum, accreting thick disk, which is then allowed to respond to the Lense-Thirring precession of the tilted black hole. The precession causes the disk to warp, beginning at the inner edge and moving out on roughly the Lense-Thirring precession timescale. The propagation of the warp stops at a radius in the disk at which other dynamical timescales, primarily the azimuthal sound-crossing time, become shorter than the precession time. At this point, the warp effectively freezes into the disk and the evolution becomes quasi-static, except in cases where the sound-crossing time in the bulk of the disk is shorter than the local precession timescale. We see evidence that such disks undergo near solid-body precession after the initial warping has frozen in. Simultaneous to the warping of the disk, there is also a tendency for the midplane to align with the symmetry plane of the black hole due to the preferential accretion of the most tilted disk gas. This alignment is not as pronounced, however, as it would be if more efficient angular momentum transport (e.g. from viscosity or magneto-rotational instability) were considered.

We present results from two-dimensional numerical simulations of the interactions between magnetized shocks and radiative clouds. Our primary goal is to characterize the dynamical evolution of the shocked clouds. We perform runs in both the strong and weak magnetic field limits and consider three different field orientations. For the geometries considered, we generally find that magnetic fields external to, but concentrated near, the surface of the cloud suppress the growth of destructive hydrodynamic instabilities. External fields also increase the compression of the cloud by effectively acting as a confinement mechanism driven by the interstellar flow and local field stretching. This can have a dramatic effect on both the efficiency of radiative cooling, which tends to increase with increasing magnetic field strength, and on the size and distribution of condensed cooled fragments. In contrast, fields acting predominately internally to the cloud tend to resist compression, thereby inhibiting cooling. We observe that, even at modest strengths, internal fields can completely suppress low-temperature cooling.

We examine the efficiency with which supernova-enriched gas may be ejected from dwarf disk galaxies, using a methodology previously employed to study the self-enrichment efficiency of dwarf spheroidal systems. Unlike previous studies that focused on highly concentrated starbursts, in the current work we consider discrete supernova events spread throughout various fractions of the disk. We model disk systems having gas masses of 10^8 and 10^9 solar masses with supernova rates of 30, 300, and 3000 per Myr. The supernova events are confined to the midplane of the disk, but distributed over radii of 0, 30, and 80% of the disk radius, consistent with expectations for Type II supernovae. In agreement with earlier studies, we find that the enriched material from supernovae is largely lost when the supernovae are concentrated near the nucleus, as expected for a starburst event. In contrast, however, we find the loss of enriched material to be much less efficient when the supernovae occur over even a relatively small fraction of the disk. The difference is due to the ability of the system to relax following supernova events that occur over more extended regions. Larger physical separations also reduce the likelihood of supernovae going off within low-density "chimneys" swept out by previous supernovae. We also find that, for the most distributed systems, significant metal loss is more likely to be accompanied by significant mass loss. A comparison with theoretical predications indicates that, when undergoing self-regulated star formation, galaxies in the mass range considered shall efficiently retain the products of Type II supernovae.

Accumulating observational evidence for a number of radio galaxies suggests an association between their jets and regions of active star formation. The standard picture is that shocks generated by the jet propagate through an inhomogeneous medium and trigger the collapse of overdense clouds, which then become active star-forming regions. In this contribution, we report on recent hydrodynamic simulations of radiative shock-cloud interactions using two different cooling models: an equilibrium cooling-curve model assuming solar metallicities and a non-equilibrium chemistry model appropriate for primordial gas clouds. We consider a range of initial cloud densities and shock speeds in order to quantify the role of cooling in the evolution. Our results indicate that for moderate cloud densities (>1 cm^{-3}) and shock Mach numbers (<20), cooling processes can be highly efficient and result in more than 50% of the initial cloud mass cooling to below 100 K. We also use our results to estimate the final H_2 mass fraction for the simulations that use the non-equilibrium chemistry package. This is an important measurement, since H_2 is the dominant coolant for a primordial gas cloud. We find peak H_2 mass fractions of >0.01 and total H_2 mass fractions of >10^{-5} for the cloud gas. Finally, we compare our results with the observations of jet-induced star formation in ``Minkowski's Object.'' We conclude that its morphology, star formation rate (~ 0.3M_solar/yr) and stellar mass (~ 1.2 x 10^7 M_solar) can be explained by the interaction of a 90,000 km/s jet with an ensemble of moderately dense (~ 10 cm^{-3}), warm (10^4 K) intergalactic clouds in the vicinity of its associated radio galaxy at the center of the galaxy cluster.

We explore several models which might be proposed to explain recent possible detections of high-energy (TeV) gamma rays in association with low-energy gamma-ray bursts (GRBs). Likely values (and/or upper limits) for the source energies in low- and high-energy gamma rays and hadrons are deduced for the burst sources associated with possible TeV gamma-ray detections by the Project GRAND array. Possible spectra for energetic gammas are deduced for three models: 1) inverse-Compton scattering of ambient photons from relativistic electrons; 2) proton-synchrotron emission; and 3) inelastic scattering of relativistic protons from ambient photons creating high-energy neutral pions, which decay into high-energy photons. These models rely on some basic assumptions about the GRB properties, e.g. that: the low- and high-energy gamma rays are produced at the same location; the time variability of the high-energy component can be estimated from the FWHM of the highest peak in the low-energy gamma ray light curve; and the variability-luminosity relation of Fenimore & Ramirez-Ruiz (2000) gives a reliable estimate of the redshifts of these bursts. We also explore the impact of each of these assumptions upon our models. We conclude that the energetic requirements are difficult to satisfy for any of these models unless, perhaps, either the photon beaming angle is much narrower for the high-energy component than for the low-energy GRB or the bursts occur at very low redshifts (z<0.01). Nevertheless, we find that the energetic requirements are most easily satisfied if TeV gamma rays are produced predominantly by inverse-Compton scattering with a magnetic field strength well below equipartition or by proton-synchrotron emission with a magnetic field strength near equipartition.

(Abridged) Many dwarf galaxies exhibit sub-Solar metallicities, with some star-to-star variation, despite often containing multiple generations of stars. The total metal content in these systems is much less than expected from the heavy element production of massive stars in each episode of star formation. Such a deficiency implies that a substantial fraction of the enriched material has been lost from these small galaxies. Mass ejection from dwarf galaxies may have important consequences for the evolution of the intergalactic medium and for the evolution of massive galaxies, which themselves may have formed via the merger of smaller systems. We report here the results of three-dimensional simulations of the evolution of supernova-enriched gas within dwarf spheroidal galaxies (dSph's), with the aim of determining the retention efficiency of supernova ejecta. We consider two galaxy models, selected to represent opposite ends of the dSph sequence. For each model galaxy we investigate a number of scenarios, ranging from a single supernova in smooth gas distributions to more complex multiple supernovae in highly disturbed gas distributions. The results of these investigations suggest that, for low star-formation efficiencies, it is difficult to completely expel the enriched material from the galaxy. Most of the enriched gas is, however, lost from the core of the galaxy following multiple supernovae, especially if the interstellar medium is already highly disturbed by processes such as photo-ionization and stellar winds. If subsequent star formation occurs predominantly within the core where most of the residual gas is concentrated, then these results could explain the poor self-enrichment efficiency observed in dwarf galaxies.

We have developed a new massively-parallel radiation-hydrodynamics code (Cosmos) for Newtonian and relativistic astrophysical problems that also includes radiative cooling, self-gravity, and non-equilibrium, multi-species chemistry. Several numerical methods are implemented for the hydrodynamics, including options for both internal and total energy conserving schemes. Radiation is treated using flux-limited diffusion. The chemistry incorporates 27 reactions, including both collisional and radiative processes for atomic hydrogen and helium gases, and molecular hydrogen chains. In this paper we discuss the equations and present results from test problems carried out to verify the robustness and accuracy of our code in the Newtonian regime. An earlier paper presented tests of the relativistic capabilities of Cosmos.

Recent results from linear perturbation theory suggest that first-order cosmological quark-hadron phase transitions occurring as deflagrations may be ``borderline'' unstable, and those occurring as detonations may give rise to growing modes behind the interface boundary. However, since nonlinear effects can play important roles in the development of perturbations, unstable behavior cannot be asserted entirely by linear analysis, and the uncertainty of these recent studies is compounded further by nonlinearities in the hydrodynamics and self-interaction fields. In this paper we investigate the growth of perturbations and the stability of quark-hadron phase transitions in the early Universe by solving numerically the fully nonlinear relativistic hydrodynamics equations coupled to a scalar field with a quartic self-interaction potential regulating the transitions. We consider single, perturbed, phase transitions propagating either by detonation or deflagration, as well as multiple phase and shock front interactions in 1+2 dimensional spacetimes.

We report on a study of sub-TeV ($E_\gamma > 10$ GeV) gamma-ray-induced muon secondaries in coincidence with BATSE gamma ray bursts (GRBs). Each TeV gamma ray striking the atmosphere produces $\approx$0.2 muons whose identity and angle can be measured by the Project GRAND array. Eight GRB candidates were studied; seven were selected based upon the optimum product of (detected BATSE fluence) $\times$ (GRAND's acceptance). One candidate was added because it was reported as a detection by the Milagrito collaboration. Seven candidates show a positive, though not statistically significant, muon excess. The most significant possible coincidence shows an excess of $466 \pm 171$ muons during the BATSE T90 time interval. It is also of interest that the event with the second most significant signal to noise (1.2$\sigma$) is the Milagrito event. If our most significant event is real, the implied fluence of energetic ($> 10$ GeV) gammas necessary to account for the observed muon excess would require that most of the GRB fluence arrived in the form of energetic gammas.

High resolution, non-oscillatory, central difference (NOCD) numerical schemes are introduced as alternatives to more traditional artificial viscosity (AV) and Godunov methods for solving the fully general relativistic hydrodynamics equations. These new approaches provide the advantages of Godunov methods in capturing ultra-relativistic flows without the cost and complication of Riemann solvers, and the advantages of AV methods in their speed, ease of implementation, and general applicability without explicitly using artificial viscosity for shock capturing. Shock tube, wall shock, and dust accretion tests, all with adiabatic equations of state, are presented and compared against equivalent solutions from both AV and Godunov based codes. In the process we address the accuracy of time-explicit NOCD and AV methods over a wide range of Lorentz factors.

We present an analysis of the longest timescale microlensing events discovered by the MACHO Collaboration during a 7 year survey of the Galactic bulge. We find 6 events that exhibit very strong microlensing parallax signals due, in part, to accurate photometric data from the GMAN and MPS collaborations. The microlensing parallax fit parameters are used in a likelihood analysis, which is able to estimate the distance and masses of the lens objects based upon a standard model of the Galactic velocity distribution. This analysis indicates that the most likely masses of 5 of the 6 lenses are > 1 Msun, which suggests that a substantial fraction of the Galactic lenses may be massive stellar remnants. This could explain the observed excess of long timescale microlensing events. The lenses for events MACHO-96-BLG-5 and MACHO-98-BLG-6 are the most massive, with mass estimates of M/Msun = 6 +10/-3 and M/Msun = 6 +7/-3, respectively. The observed upper limits on the absolute brightness of main sequence stars for these lenses are < 1 Lsun, so both lenses are black hole candidates. The black hole interpretation is also favored by a likelihood analysis with a Bayesian prior using a conventional model for the lens mass function. We consider the possibility that the source stars for some of these 6 events may lie in the foreground or background of the Galactic bulge, but we find that this is unlikely. Future HST observations of these events can either confirm the black hole lens hypothesis or detect the lens stars and provide a direct measurement of their masses. Future observations of similar events by SIM or the Keck or VLTI interferometers will allow direct measurements of the lens masses for stellar remnant lenses as well.

The Bardeen-Petterson effect around a rapidly-rotating compact object causes a tilted accretion disk to warp into the equatorial plane of the rotating body. Viscous forces cause the accretion flow to divide into two distinct regions - an inner aligned accretion disk and an outer tilted accretion disk. The transition between these two occurs at a characteristic radius that depends on the mass and angular momentum of the central object and possibly on the accretion rate through the disk. We propose that accreting material passing through the transition region may generate quasi-periodic brightness oscillations (QPOs) such as have been observed in a number of X-ray binaries. We show that this effect may be present in the black-hole X-ray binary GRO J1655-40. We also argue that the QPO frequency range predicted by this model is consistent with observed QPO frequencies in both black-hole and neutron-star low-mass X-ray binaries.

We have reconstructed possible orbits for a collection of stars located within 0.5 arcsec of Sgr A*. These orbits are constrained by observed stellar positions and angular proper motions. The construction of such orbits serves as a baseline from which to search for possible deviations due to the unseen mass distribution in the central 1000 AU of the Galaxy. We also discuss the likelihood that some of these stars may eventually exhibit detectable relativistic effects, allowing for interesting tests of general relativity around the 2.6 x 10^6 solar mass central object.

We present observations of microlensing event MACHO-98-BLG-35 which reached a peak magnification factor of almost 80. These observations by the Microlensing Planet Search (MPS) and the MOA Collaborations place strong constraints on the possible planetary system of the lens star and show intriguing evidence for a low mass planet with a mass fraction $4\times 10^{-5} \leq \epsilon \leq 2\times 10^{-4}$. A giant planet with $\epsilon = 10^{-3}$ is excluded from 95% of the region between 0.4 and 2.5 $R_E$ from the lens star, where $R_E$ is the Einstein ring radius of the lens. This exclusion region is more extensive than the generic "lensing zone" which is $0.6 - 1.6 R_E$. For smaller mass planets, we can exclude 57% of the "lensing zone" for $\epsilon = 10^{-4}$ and 14% of the lensing zone for $\epsilon = 10^{-5}$. The mass fraction $\epsilon = 10^{-5}$ corresponds to an Earth mass planet for a lensing star of mass $\sim 0.3 \msun$. A number of similar events will provide statistically significant constraints on the prevalence of Earth mass planets. In order to put our limits in more familiar terms, we have compared our results to those expected for a Solar System clone averaging over possible lens system distances and orientations. We find that such a system is ruled out at the 90% confidence level. A copy of the Solar System with Jupiter replaced by a second Saturn mass planet can be ruled out at 70% confidence. Our low mass planetary signal (few Earth masses to Neptune mass) is significant at the $4.5\sigma$ confidence level. If this planetary interpretation is correct, the MACHO-98-BLG-35 lens system constitutes the first detection of a low mass planet orbiting an ordinary star without gas giant planets.

We present the observations of the binary lensing event MACHO-98-SMC-1 conducted at the Mt.~Stromlo 74" telescope by the Microlensing Planet Search (MPS) collaboration. The MPS data constrain the first caustic crossing to have occurred after 1998 June 5.55 UT and thus directly rule out one of the two fits presented by the PLANET collaboration (model II). This substantially reduces the uncertainty in the the relative proper motion estimations of the lens object. We perform joint binary microlensing fits of the MPS data together with the publicly available data from the EROS, MACHO/GMAN and OGLE collaborations. We also study the binary lens fit parameters previously published by the PLANET and MACHO/GMAN collaborations by using them as initial values for $\chi^2$ minimization. Fits based on the PLANET model I appear to be in conflict with the GMAN-CTIO data. From our best fit, we find that the lens system has a proper motion of $\mu = 1.3\pm 0.2 \kmsk$ with respect to the source, which implies that the lens system is most likely to be located in the Small Magellanic Cloud strengthening the conclusion of previous reports.