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.

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.

The observable quantities that carry the most information regarding the structures of the images of black holes in the interferometric observations with the Event Horizon Telescope are the closure phases along different baseline triangles. We use long time span, high cadence, GRMHD+radiative transfer models of Sgr A$^*$ to investigate the expected variability of closure phases in such observations. We find that, in general, closure phases along small baseline triangles show little variability, except in the cases when one of the triangle vertices crosses one of a small regions of low visibility amplitude. The closure phase variability increases with the size of the baseline triangle, as larger baselines probe the small-scale structures of the images, which are highly variable. On average, the jet-dominated MAD models show less closure phase variability than the disk-dominated SANE models, even in the large baseline triangles, because the images from the latter are more sensitive to the turbulence in the accretion flow. Our results suggest that image reconstruction techniques need to explicitly take into account the closure phase variability, especially if the quality and quantity of data allow for a detailed characterization of the nature of variability. This also implies that, if image reconstruction techniques that rely on the assumption of a static image are utilized, regions of the $u-v$ space that show a high level of variability will need to be identified and excised.

In general relativity static gaseous atmospheres may be in hydrostatic balance in the absence of a supporting stellar surface, provided that the luminosity is close to the Eddington value. We construct analytic models of optically thin, spherically symmetric shells supported by the radiation pressure of a luminous central body in the Schwarzschild metric. Opacity is assumed to be dominated by Thomson scattering. The inner parts of the atmospheres, where the luminosity locally has supercritical values, are characterized by a density and pressure inversion. The atmospheres are convectively and Rayleigh-Taylor stable, and there is no outflow of gas.

We study the dynamics of radiation pressure supported tori around Schwarzschild black holes, focusing on their oscillatory response to an external perturbation. Using KORAL, a general relativistic radiation hydrodynamics code capable of modeling all radiative regimes from the optically thick to the optically thin, we monitor a sample of models at different initial temperatures and opacities, evolving them in two spatial dimensions for $\sim 165$ orbital periods. The dynamics of models with high opacity is very similar to that of purely hydrodynamics models, and it is characterized by regular oscillations which are visible also in the light curves. As the opacity is decreased, the tori quickly and violently migrate towards the gas-pressure dominated regime, collapsing towards the equatorial plane. When the spectra of the $L_2$ norm of the mass density are considered, high frequency inertial-acoustic modes of oscillations are detected (with the fundamental mode at a frequency $68 M_{\rm BH}^{-1}\,\rm Hz$), in close analogy to the phenomenology of purely hydrodynamic models. An additional mode of oscillation, at a frequency $129 M_{\rm BH}^{-1}\,\rm Hz$, is also found, which can be unambiguously attributed to the radiation. The spectra extracted from the light curves are typically more noisy, indicating that in a real observation such modes would not be easily detected.

For the first time ever, we derive equations governing the time-evolution of fully relativistic slim accretion disks in the Kerr metric, and numerically construct their detailed non-stationary models. We discuss applications of these general results to a possible limit-cycle behavior of thermally unstable disks. Our equations and numerical method are applicable in a wide class of possible viscosity prescriptions, but in this paper we use a diffusive form of the "standard alpha prescription" that assumes the viscous torque is proportional to the total pressure. In this particular case, we find that the parameters which dominate the limit-cycle properties are the mass-supply rate and the value of the alpha-viscosity parameter. Although the duration of the cycle (or the outburst) does not exhibit any clear dependence on the black hole spin, the maximal outburst luminosity (in the Eddington units) is positively correlated with the spin value. We suggest a simple method for a rough estimate of the black hole spin based on the maximal luminosity and the ratio of outburst to cycle durations. We also discuss a temperature-luminosity relation for the Kerr black hole accretion discs limit-cycle. Based on these results we discuss the limit-cycle behavior observed in microquasar GRS 1915+105. We also extend this study to several non-standard viscosity prescriptions, including a "delayed heating" prescription recently stimulated by the recent MHD simulations of accretion disks.