The object G2 was recently discovered descending into the gravitational potential of the supermassive black hole (BH) Sgr A*. We test the photoionized cloud scenario, determine the cloud properties, and estimate the emission during the pericenter passage. The incident radiation is computed starting from the individual stars at the locations of G2. The radiative transfer calculations are conducted with CLOUDY code and $2011$ broadband and line luminosities are fitted. The spherically symmetric, tidally distorted, and magnetically arrested cloud shapes are tested with both the interstellar medium dust and $10$~nm graphite dust. The best-fitting magnetically arrested model has the initial density $n_{\rm init}=1.8\times10^5{\rm cm}^{-3}$, initial radius $R_{\rm init}=2.2\times10^{15}{\rm cm}=17 {\rm mas}$, mass $m_{\rm cloud}=4M_{\rm Earth}$, and dust relative abundance $A=0.072$. It provides a good fit to $2011$ data, is consistent with the luminosities in $2004$ and $2008$, and reaches an agreement with the observed size. We revise down the predicted radio and X-ray bow shock luminosities to be below the quiescent level of Sgr A*, which readily leads to non-detection in agreement to observations. The magnetic energy dissipation in the cloud at the pericenter coupled with more powerful irradiation may lead to an infrared source with an apparent magnitude $m_{L'}\approx13.0$. No shock into the cloud and no X-rays are expected from cloud squeezing by the ambient gas pressure. Larger than previously estimated cloud mass $m_{\rm cloud}=(4-20)M_{\rm Earth}$ may produce a higher accretion rate and a brighter state of Sgr A* as the debris descend onto the BH.

Based on the dynamical black hole (BH) mass estimates, NGC3115 hosts the closest billion solar mass BH. Deep studies of the center revealed a very underluminous active galactic nucleus (AGN) immersed in an old massive nuclear star cluster. Recent $1$~Ms \textit{Chandra} X-ray visionary project observations of the NGC3115 nucleus resolved hot tenuous gas, which fuels the AGN. In this paper we connect the processes in the nuclear star cluster with the feeding of the supermassive BH. We model the hot gas flow sustained by the injection of matter and energy from the stars and supernova explosions. We incorporate electron heat conduction as the small-scale feedback mechanism, the gravitational pull of the stellar mass, cooling, and Coulomb collisions. Fitting simulated X-ray emission to the spatially and spectrally resolved observed data, we find the best-fitting solutions with $\chi^2/dof=1.00$ for $dof=236$ both with and without conduction. The radial modeling favors a low BH mass $<1.3\times10^{9}M_\odot$. The best-fitting supernova rate and the best-fitting mass injection rate are consistent with their expected values. The stagnation point is at $r_{\rm st}\lesssim1$arcsec, so that most of gas, including the gas at a Bondi radius $r_B=2-4$arcsec, outflows from the region. We put an upper limit on the accretion rate at $2\times10^{-3}M_\odot{\rm yr}^{-1}$. We find a shallow density profile $n\propto r^{-\beta}$ with $\beta\approx1$ over a large dynamic range. This density profile is determined in the feeding region $0.5-10$arcsec as an interplay of four processes and effects: (1) the radius-dependent mass injection, (2) the effect of the galactic gravitational potential, (3) the accretion flow onset at $r\lesssim1$arcsec, and (4) the outflow at $r\gtrsim1$arcsec. The gas temperature is close to the virial temperature $T_v$ at any radius.

An object called G2 was recently discovered moving towards the supermassive black hole in the Galactic Center. G2 emits infrared (IR) lines and continuum, which allows constraining its properties. The question is still unresolved whether G2 has a central windy star or it is a coreless cloud. Assuming the object is a cloud originating near the apocenter I perform line/continuum IR diagnostics, revisit estimates of non-thermal emission from pericenter passage, and speculate about future observational prospects. This work is partially reported in arXiv:1309.2282 and partially consists of new ideas discussed at the conference.

Observational confirmation of hot accretion model predictions has been hindered by the challenge to resolve spatially the Bondi radii of black holes with X-ray telescopes. Here, we use the Megasecond Chandra X-ray Visionary Project (XVP) observation of the NGC~3115 supermassive black hole to place the first direct observational constraints on the spatially and spectroscopically resolved structures of the X-ray emitting gas inside the Bondi radius of a black hole. We measured temperature and density profiles of the hot gas from a fraction out to tens of the Bondi radius (R_B = 2.4-4.8 arcsec = 112-224 pc). The projected temperature jumps significantly from ~0.3 keV beyond 5 arcsec to ~0.7 keV within ~4-5 arcsec, but then abruptly drops back to ~0.3 keV within ~3 arcsec. This is contrary to the expectation that the temperature should rise toward the center for a radiatively inefficient accretion flow. A hotter thermal component of ~1 keV inside 3 arcsec (~150 pc) is revealed using a two component thermal model, with the cooler ~0.3 keV thermal component dominating the spectra. We argue that the softer emission comes from diffuse gas physically located within $\sim 150$~pc from the black hole. The density profile is broadly consistent with rho ~ r^{-1} within the Bondi radius for either the single temperature or the two-temperature model. The X-ray data alone with physical reasoning argue against the absence of a black hole, supporting that we are witnessing the onset of the gravitational influence of the supermassive black hole.

Most supermassive black holes (SMBHs) are accreting at very low levels and are difficult to distinguish from the galaxy centers where they reside. Our own Galaxy's SMBH provides a uniquely instructive exception, and we present a close-up view of its quiescent X-ray emission based on 3 mega-second of Chandra observations. Although the X-ray emission is elongated and aligns well with a surrounding disk of massive stars, we can rule out a concentration of low-mass coronally active stars as the origin of the emission based on the lack of predicted Fe Kalpha emission. The extremely weak H-like Fe Kalpha line further suggests the presence of an outflow from the accretion flow onto the SMBH. These results provide important constraints for models of the prevalent radiatively inefficient accretion state.

High-frequency quasi-periodic oscillations (QPOs) appear in general-relativistic magnetohydrodynamic simulations of magnetically choked accretion flows around rapidly rotating black holes (BHs). We perform polarized radiative transfer calculations with ASTRORAY code to explore the manifestations of these QPOs for SgrA*. We construct a simulation-based model of a radiatively inefficient accretion flow and find model parameters by fitting the mean polarized source spectrum. The simulated QPOs have a total sub-mm flux amplitude up to 5% and a linearly polarized flux amplitude up to 2%. The oscillations reach high levels of significance 10-30\sigma and high quality factors Q~5. The oscillation period T~100M~35min corresponds to the rotation period of the BH magnetosphere that produces a trailing spiral in resolved disk images. The total flux signal is significant over noise for all tested frequencies 87GHz, 230GHz, and 857GHz and inclination angles 10deg, 37deg, and 80deg. The non-detection in the 230GHz Sub-Millimeter Array light curve is consistent with a low signal level and a low sampling rate. The presence of sub-mm QPOs in SgrA* will be better tested with the Atacama Large Millimeter Array.

A highly unusual pair of a gamma-ray burst (GRB) GRB060218 and an associated supernova SN2006aj has puzzled theorists for years. A supernova shock breakout and a jet from a newborn stellar mass compact object were put forward to explain its multiwavelength signature. We propose that the source is naturally explained by another channel, a tidal disruption of a white dwarf (WD) by an intermediate mass black hole (IMBH). The tidal disruption is accompanied by a tidal pinching, which leads to the ignition of a WD and a supernova. Some debris falls back onto the IMBH, forms a disk, which quickly amplifies the magnetic field, and launches a jet. We successfully fit soft X-ray spectrum with the Comptonized blackbody emission from a jet photosphere. The optical/UV emission is consistent with self-absorbed synchrotron from the expanding jet front. The accretion rate temporal dependence Mdot(t) in a tidal disruption provides a good fit to soft X-ray lightcurve. The IMBH mass is found to be about 10^4Msun in three independent estimates: (1) fitting tidal disruption Mdot(t) to soft X-ray lightcurve; (2) computing the jet base radius in a jet photospheric emission model; (3) inferring the central BH mass based on a host dwarf galaxy stellar mass. The supernova position is consistent with the center of the host galaxy, while low supernova ejecta mass is consistent with a WD mass. High expected rate of tidal disruptions in dwarf galaxies is consistent with one source observed by Swift satellite over several years at GRB060218 distance of 150Mpc. The encounters with the WDs provide a lot of fuel for IMBH growth.

We present a qualitative picture of prompt emission from tidal disruptions of white dwarfs (WD) by intermediate mass black holes (IMBH). The smaller size of an IMBH compared to a supermassive black hole and a smaller tidal radius of a WD disruption lead to a very fast event with high peak luminosity. Magnetic field is generated in situ following the tidal disruption, which leads to effective accretion. Since large-scale magnetic field is also produced, geometrically thick super-Eddington inflow leads to a relativistic jet. The dense jet possesses a photosphere, which emits quasi-thermal radiation in soft X-rays. The source can be classified as a long low-luminosity gamma-ray burst (ll-GRB). Tidal compression of a WD causes nuclear ignition, which is observable as an accompanying supernova. We suggest that GRB060218 and SN2006aj is such a pair of ll-GRB and supernova. We argue that in a flux-limited sample the disruptions of WDs by IMBHs are more frequent then the disruptions of other stars by IMBHs.

The unmatched X-ray resolution of Chandra allows probing the gas flow near quiescent supermassive black holes (BHs). The radius of BH gravitational influence on gas, called the Bondi radius, is resolved in Sgr A* and NGC 3115. Shallow accretion flow density profiles n \sim r^{-beta} with beta=0.7-1.0 were found for Sgr A* and NGC 3115 with the help of Chandra. We construct self-consistent models with gas feeding and dynamics from near the Bondi radius to the event horizon to explain the observations. Gas is mainly supplied to the region by hot colliding stellar winds. Small-scale feedback such as conduction effectively flattens the density profile from steep beta=1.5 in a Bondi flow. We further constrain density and temperature profiles using the observed radio/sub-mm radiation emitted near the event horizon. We discuss the present state of our numerical model and its qualitative features, such as the role of the galactic gravitational potential and the random motion of wind-emitting stars.

We briefly summarize the method of simulating Sgr A* polarized sub-mm spectra from the accretion flow and fitting the observed spectrum. The dynamical flow model is based on three-dimensional general relativistic magneto hydrodynamic simulations. Fully self-consistent radiative transfer of polarized cyclo-synchrotron emission is performed. We compile a mean sub-mm spectrum of Sgr A* and fit it with the mean simulated spectra. We estimate the ranges of inclination angle theta=42-75deg, mass accretion rate Mdot=(1.4-7.0)*10^{-8}Msun/yr, and electron temperature Te=(3-4)*10^{10}K at 6M. We discuss multiple caveats in dynamical modeling, which must be resolved to make further progress.

We obtain estimates of Sgr A* accretion flow and black hole parameters by fitting polarized sub-mm observations with spectra computed using three-dimensional (3D) general relativistic (GR) magnetohydrodynamical (MHD) (GRMHD) simulations. Observations are compiled from averages over many epochs from reports in 29 papers for estimating the mean fluxes Fnu, linear polarization (LP) fractions, circular polarization (CP) fractions, and electric vector position angles (EVPAs). GRMHD simulations are computed with dimensionless spins a_*=0,0.5,0.7,0.9,0.98 over a 20,000M time interval. We perform fully self-consistent GR polarized radiative transfer using our new code to explore the effects of spin a_*, inclination angle \theta, position angle (PA), accretion rate Mdot, and electron temperature Te (Te is reported for radius 6M). By fitting the mean sub-mm fluxes and LP/CP fractions, we obtain estimates for these model parameters and determine the physical effects that could produce polarization signatures. Our best bet model has a_*=0.5, \theta=75deg, PA=115deg, Mdot=4.6*10^{-8}M_Sun/year, and Te=3.1*10^10K at 6M. The sub-mm CP is mainly produced by Faraday conversion as modified by Faraday rotation, and the emission region size at 230GHz is consistent with the VLBI size of 37microas. Across all spins, model parameters are in the ranges \theta=42deg-75deg, Mdot=(1.4-7.0)*10^{-8}M_Sun/year, and Te=(3-4)*10^10K. Polarization is found both to help differentiate models and to introduce new observational constraints on the effects of the magnetic field that might not be fit by accretion models so-far considered.

We present numerical relativity results of tidal disruptions of white dwarfs from ultra-close encounters with a spinning, intermediate mass black hole. These encounters require a full general relativistic treatment of gravity. We show that the disruption process and prompt accretion of the debris strongly depend on the magnitude and orientation of the black hole spin. However, the late-time accretion onto the black hole follows the same decay, $\dot{M}$ ~ t^{-5/3}, estimated from Newtonian gravity disruption studies. We compute the spectrum of the disk formed from the fallback material using a slim disk model. The disk spectrum peaks in the soft X-rays and sustains Eddington luminosity for 1-3 yrs after the disruption. For arbitrary black hole spin orientations, the disrupted material is scattered away from the orbital plane by relativistic frame dragging, which often leads to obscuration of the inner fallback disk by the outflowing debris. The disruption events also yield bursts of gravitational radiation with characteristic frequencies of ~3.2 Hz and strain amplitudes of ~10^{-18} for galactic intermediate mass black holes. The optimistic rate of considered ultra-close disruptions is consistent with no sources found in ROSAT all-sky survey. The future missions like Wide-Field X-ray Telescope (WFXT) could observe dozens of events.

We provide precise fitting formulae for Faraday conversion and rotation coefficients in uniformly magnetized relativistic plasma. The formulae are immediately applicable to Rotation Measure and Circular Polarization (CP) production in jets and hot accretion flows. We show the recipe and results for arbitrary isotropic particle distributions, in particular thermal and power-law. The exact Faraday conversion coefficient is found to approach zero with the increasing particle energy. The non-linear corrections of Faraday conversion and rotation coefficients are found essential for reliable CP inter- pretation of Sgr A*.

The X-ray spectra of accretion discs of eight stellar-mass black holes have been analyzed to date using the thermal continuum fitting method, and the spectral fits have been used to estimate the spin parameters of the black holes. However, the underlying model used in this method of estimating spin is the general relativistic thin-disc model of Novikov & Thorne, which is only valid for razor-thin discs. We therefore expect errors in the measured values of spin due to inadequacies in the theoretical model. We investigate this issue by computing spectra of numerically calculated models of thin accretion discs around black holes, obtained via three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations. We apply the continuum fitting method to these computed spectra to estimate the black hole spins and check how closely the values match the actual spin used in the GRMHD simulations. We find that the error in the dimensionless spin parameter is up to about 0.2 for a non-spinning black hole, depending on the inclination. For black holes with spins of 0.7, 0.9 and 0.98, the errors are up to about 0.1, 0.03 and 0.01 respectively. These errors are comparable to or smaller than those arising from current levels of observational uncertainty. Furthermore, we estimate that the GRMHD simulated discs from which these error estimates are obtained correspond to effective disc luminosities of about 0.4-0.7 Eddington, and that the errors will be smaller for discs with luminosities of 0.3 Eddington or less, which are used in the continuum-fitting method. We thus conclude that use of the Novikov-Thorne thin-disc model does not presently limit the accuracy of the continuum-fitting method of measuring black hole spin.

The rising amount of polarized observations of relativistic sources requires the correct theory for proper model fitting. The equations for general relativistic (GR) polarized radiative transfer are derived starting from the Boltzmann equation and basic ideas of general relativity. The derivation is aimed at providing a practical guide to reproducing the synchrotron part of radio & sub-mm emission from low luminosity active galactic nuclei (LLAGNs), in particular Sgr A*, and jets. The recipe for fast exact calculation of cyclo-synchrotron emissivities, absorptivities, Faraday rotation and conversion coefficients is given for isotropic particle distributions. The multitude of physical effects influencing simulated spectrum is discussed. The application of the prescribed technique is necessary to determine the black hole (BH) spin in LLAGNs, constraining it with all observations of total flux, linear and circular polarization fractions, and electric vector position angle as functions of the observed frequency.

We propose a two-temperature radial inflow-outflow model near Sgr A* with self-consistent feeding and conduction. Stellar winds from individual stars are considered to find the rates of mass injection and energy injection. These source terms help to partially eliminate the boundary conditions on the inflow. Electron thermal conduction is crucial for inhibiting the accretion. Energy diffuses out from several gravitational radii, unbinding more gas at several arcseconds and limiting the accretion rate to <1% of Bondi rate. We successfully fit the X-Ray surface brightness profile found from the extensive Chandra observations and reveal the X-Ray point source in the center. The super-resolution technique allows us to infer the presence and estimate the unabsorbed luminosity $L\approx4\cdot10^{32}{\rm erg s^{-1}}$ of the point source. The employed relativistic heat capacity and direct heating of electrons naturally lead to low electron temperature $T_e\approx 4\cdot10^{10}$ K near the black hole. Within the same model we fit 86 GHz optically thick emission and obtain the order of magnitude agreement of Faraday rotation measure, thus achieving a single accretion model suitable at all radii.

We propose a 2-temperature radial dynamical model of plasma flow near Sgr A* and fit the bremsstrahlung emission to extensive quiescent X-Ray Chandra data. The model extends from several arcseconds to black hole (BH) gravitational radius, describing the outer accretion flow together with the infalling region. The model incorporates electron heat conduction, relativistic heat capacity of particles and feeding by stellar winds. Stellar winds from each star are considered separately as sources of mass, momentum and energy. Self-consistent search for the stagnation and sonic points is performed. Most of gas is found to outflow from the region. The accretion rate is limited to below 1% of Bondi rate due to the effect of thermal conduction enhanced by entropy production in a turbulent flow. The X-Ray brightness profile proves too steep near the BH, thus a synchrotron self-Compton point source is inferred with luminosity L=3x10^32erg/s. We fit the sub-mm emission from the inner flow, thus aiming at a single model of Sgr A* accretion suitable at any radius.

We present the combination of dynamical accretion model based on 3D GRMHD simulations and general relativistic (GR) polarized radiative transfer. We write down the formalism of and perform the GR ray-tracing of cyclo-synchrotron radiation through the model of accretion flow in Sagittarius A*. GR polarimetric imaging is presented as well as the results for spectrum for a probable set of spins and orientations. Precise fitting formulae for Faraday rotation and Faraday conversion coefficients are employed for thermal plasma. The axisymmetic flow pattern and the magnetic field geometry correspond to averaged 3D GRMHD simulations near the black hole, whereas the analytic model was used far from the black hole. The density scaling is found by fitting the sub-mm flux. Spin $a=0.7$ and inclination angle $\theta=0.6$ produce the best fit to sub-mm flux and linear polarization fraction.

The transfer of polarized radiation in magnetized and non-magnetized relativistic plasmas is an area of research with numerous flaws and gaps. The present paper is aimed at filling some gaps and eliminating the flaws. Starting from a Trubnikov's linear response tensor for a vacuum wave with ${\bf k}=\omega/c$ in thermal plasma, the analytic expression for the dielectric tensor is found in the limit of high frequencies. The Faraday rotation and Faraday conversion measures are computed in their first orders in the ratio of the cyclotron frequency $\Omega_0$ to the observed frequency $\omega$. The computed temperature dependencies of propagation effects bridge the known non-relativistic and ultra-relativistic limiting formulas. The fitting expressions are found for high temperatures, where the higher orders in $\Omega_0/\omega$ cannot be neglected. The plasma eigenmodes are found to become linearly polarized at much larger temperatures than thought before. The results are applied to the diagnostics of the hot ISM, hot accretion flows, and jets.

The first spherical accretion model was developed 55 years ago, but the theory is yet far from being complete. The real accretion flow was found to be time-dependent and turbulent. This paper presents the minimal MHD spherical accretion model that separately deals with turbulence. Treatment of turbulence is based on simulations of several regimes of collisional MHD. The effects of freezing-in amplification, dissipation, dynamo action, isotropization, and constant magnetic helicity are self-consistently included. The assumptions of equipartition and magnetic field isotropy are released. Correct dynamics of magnetized flow is calculated. Diffusion, convection, and radiation are not accounted for. Two different types of Radiatively Inefficient accretion flows are found: a transonic non-rotating flow (I), a flow with effective transport of angular momentum outward (II). Non-rotating flow has an accretion rate several times smaller than Bondi rate, because turbulence inhibits accretion. Flow with angular momentum transport has accretion rate about 10-100 times smaller than Bondi rate. The effects of highly helical turbulence, states of outer magnetization, and different equations of state are discussed. The flows were found to be convectively stable on average, despite gas entropy increases inward. The proposed model has a small number of free parameters and the following attractive property. Inner density in the non-rotating magnetized flow was found to be several times lower than density in a non-magnetized accretion. Still several times lower density is required to explain the observed low IR luminosity and low Faraday rotation measure of accretion onto Sgr A*.

Transonic accretion flow with self-consistent treatment of random magnetic field is presented.

We investigate the properties of an axisymmetric non-magnetized gas flow without angular momentum on a small compact object, in particular, on a Schwarzschild black hole in the supersonic region near the object; the velocity of the object itself is assumed to be low compared to the speed of sound at infinity. First of all, we see that the streamlines intersect (i.e., a caustic forms) on the symmetry axis at a certain distance $r_x$ from the center on the front side if the pressure gradient is neglected. The characteristic radial size of the region, in which the streamlines emerging from the sonic surface at an angle no larger than $\theta_0$ to the axis intersect, is $\Delta r= r_x\theta^2_0/3.$ To refine the flow structure in this region, we numerically compute the system in the adiabatic approximation without ignoring the pressure. We estimate the parameters of the inferred region with anomalously high matter temperature and density accompanied by anomalously high energy release.

We derive an analytical expression of the inter-arrival time distribution for a non-homogeneous Poisson process (NHPP). This expression is exact and is applicable to any time interval, finite or infinite. As an illustration, we present simulation results for three different intensity functions.