We report on a new class of solutions of black hole accretion disks that we have found through three-dimensional, global, radiative magnetohydrodynamic simulations in general relativity. It combines features of the canonical thin, slim and thick disk models but differs in crucial respects from each of them. We expect these new solutions to provide a more realistic description of black hole disks than the slim disk model. We are presenting a disk solution for a non-spinning black hole at a sub-Eddington mass accretion rate, $\dot M=0.6\,\dot M_{\rm Edd}$. By the density scale-height measure the disk appears to be thin, having a high density core near the equatorial plane of height $h_{\rho} \sim 0.1 \,r$, but most of the inflow occurs through a highly advective, turbulent, optically thick, Keplerian region that sandwiches the core and has a substantial geometrical thickness comparable to the radius, $H \sim r$. The accreting fluid is supported above the midplane in large part by the magnetic field, with the gas and radiation to magnetic pressure ratio $\beta \sim 1$, this makes the disk thermally stable, even though the radiation pressure strongly dominates over gas pressure. A significant part of the radiation emerging from the disk is captured by the black hole, so the disk is less luminous than a thin disk would be at the same accretion rate.

Neutron stars with near-Eddington observable luminosities were shown to harbor levitating atmospheres, suspended above their surface. We report a new method to simultaneously measure the mass and radius of a neutron star based on oscillations of such atmospheres. In this paper, we present an analytic derivation of a family of relativistic, oscillatory, spherically symmetric eigenmodes of the optically and geometrically thin levitating atmospheres, including the damping effects induced by the radiation drag. We discover characteristic maxima in the frequencies of the damped oscillations and show that using the frequency maxima, one can estimate mass and radius of the neutron star, given the observed frequency and the corresponding luminosity of the star during the X-ray burst. Thus, our model provides a new way to probe the stellar parameters. We also show that the ratio of any two undamped eigenfrequencies depends only on the adiabatic index of the atmosphere, while for the damped eigenfrequencies, this ratio varies with the luminosity. The damping coefficient is independent of the mode number of the oscillations. Signatures of these atmospheres' dynamics will be reflected in the source's X-ray light curves.

We run two GRRMHD simulations of super-Eddington accretion disks around a black hole and a non-magnetized, non-rotating neutron star. The neutron star was modeled using a reflective inner boundary condition. We observe the formation of a transition layer in the inner region of the disk in the neutron star simulation which leads to a larger mass outflow rate and a lower radiative luminosity over the black hole case. Sphereization of the flow leads to an observable luminosity at infinity around the Eddington value when viewed from all directions for the neutron star case, contrasting to the black hole case where collimation of the emission leads to observable luminosities about an order of magnitude higher when observed along the disk axis. We find the outflow to be optically thick to scattering, which would lead to the obscuring of any neutron star pulsations observed in corresponding ULXs.

We describe some recent results on the evolution of accretion disks around black holes and neutron stars obtained in magnetohydrodynamic (MHD) simulations in general relativity (GR) with the inclusion of radiation (GRRMHD).

We present a numerical method which evolves a two-temperature, magnetized, radiative, accretion flow around a black hole, within the framework of general relativistic radiation magnetohydrodynamics. As implemented in the code KORAL, the gas consists of two sub-components -- ions and electrons -- which share the same dynamics but experience independent, relativistically consistent, thermodynamical evolution. The electrons and ions are heated independently according to a standard prescription from the literature for magnetohydrodynamical turbulent dissipation. Energy exchange between the particle species via Coulomb collisions is included. In addition, electrons gain and lose energy and momentum by absorbing and emitting synchrotron and bremsstrahlung radiation, and through Compton scattering. All evolution equations are handled within a fully covariant framework in the relativistic fixed-metric spacetime of the black hole. Numerical results are presented for five models of low luminosity black hole accretion. In the case of a model with a mass accretion rate $\dot{M}\sim10^{-8} \dot M_{\rm Edd}$, we find that radiation has a negligible effect on either the dynamics or the thermodynamics of the accreting gas. In contrast, a model with a larger $\dot{M}\sim 4\times 10^{-4} \dot M_{\rm Edd}$ behaves very differently. The accreting gas is much cooler and the flow is geometrically less thick, though it is not quite a thin accretion disk.

In general relativity, it has been shown that radiation-supported atmospheres exist outside the surface of a radiating spherical body close to a radius where the gravitational and radiative forces balance each other. We calculate the frequency of oscillation of the incompressible radial mode of such a thin atmospheric shell and show that in the optically thin case, this particular mode is overdamped by radiation drag.

We perform hydro- and magnetohydrodynamical general relativistic simulations of a tidal disruption of a $0.1\,M_\odot$ red dwarf approaching a $10^5\,M_\odot$ non-rotating massive black hole on a close (impact parameter $\beta=10$) elliptical (eccentricity $e=0.97$) orbit. We track the debris self-interaction, circularization, and the accompanying accretion through the black hole horizon. We find that the relativistic precession leads to the formation of a self-crossing shock. The dissipated kinetic energy heats up the incoming debris and efficiently generates a quasi-spherical outflow. The self-interaction is modulated because of the feedback exerted by the flow on itself. The debris quickly forms a thick, almost marginally bound disc that remains turbulent for many orbital periods. Initially, the accretion through the black hole horizon results from the self-interaction, while in the later stages it is dominated by the debris originally ejected in the shocked region, as it gradually falls back towards the hole. The effective viscosity in the debris disc stems from the original hydrodynamical turbulence, which dominates over the magnetic component. The radiative efficiency is very low because of low energetics of the gas crossing the horizon and large optical depth that results in photon trapping.

We present a sub-grid model that emulates the magnetic dynamo operating in magnetized accretion disks. We have implemented this model in the general relativisic radiation magnetohydrodynamic (GRRMHD) code \koral, using results from local shearing sheet simulations of the magnetorotational instability to fix the parameters of the dynamo. With the inclusion of this dynamo, we are able to run 2D axisymmetric GRRMHD simulations of accretion disks for arbitrarily long times. The simulated disks exhibit sustained turbulence, with the poloidal and toroidal magnetic field components driven towards a state similar to that seen in 3D studies. Using this dynamo code, we present a set of long-duration global simulations of super-Eddington, optically-thick disks around non-spinning and spinning black holes. Super-Eddington disks around non-rotating black holes exhibit a surprisingly large efficiency, $\eta\approx0.04$, independent of the accretion rate, where we measure efficiency in terms of the total energy output, both radiation and mechanical, flowing out to infinity. Super-Eddington disks around spinning black holes are even more efficient, and appear to extract black hole rotational energy through a process similar to the Blandford-Znajek mechanism. All the simulated models are characterized by highly super-Eddington radiative fluxes collimated along the rotation axis. We also present a set of simulations that were designed to have Eddington or slightly sub-Eddington accretion rates ($\dot{M} \lesssim 2\dot M_{\rm Edd}$). None of these models reached a steady state. Instead, the disks collapsed as a result of runaway cooling, presumably because of a thermal instability.

We study the dynamical interaction of the G2 cloud with the accretion flow around Sgr A* by means of three-dimensional, hydrodynamic simulations. We show the effects of the rotating accretion flow on the evolution of G2 by projecting the cloud density onto the plane of the sky, and extracting position-velocity diagrams. We study a number of possible orientations of the cloud orbit with respect to the disk. We find that once the center of mass of the cloud has crossed the pericenter, the differences between models becomes significant. Models with the cloud counter-rotating with respect to the disk are expected to reach higher blue-shifted line of sight velocities. The spatial extent of the emission depends strongly on the cloud-to-disk inclination angle. Future imaging and spectroscopy of G2 emission will shed light both on the structure of the Sgr A* disk and on the properties of the cloud.

We calculate radio light curves produced by the bow shock that is likely to form in front of the G2 cloud when it penetrates the accretion disk of Sgr A*. The shock acceleration of the radio-emitting electrons is captured self-consistently by means of first-principles particle-in-cell simulations. We show that the radio luminosity is expected to reach maximum in early 2013, roughly a month after the bow shock crosses the orbit pericenter. We estimate the peak radio flux at 1.4 GHz to be 1.4 - 22 Jy depending on the assumed orbit orientation and parameters. We show that the most promising frequencies for radio observations are in the 0.1<nu<1 GHz range, for which the bow shock emission will be much stronger than the intrinsic radio flux for all the models considered.