Changing-look phenomenon observed now in a growing number of active galaxies challenges our understanding of the accretion process close to a black hole. We propose a simple explanation for the sources where multiple semi-periodic outbursts are observed, and the sources are operating at a few per cent of the Eddington limit. The outburst are caused by the radiation pressure instability operating in the narrow ring between the standard gas-dominated outer disk and the hot optically thin inner Advection-Dominated Accretion Flow. The corresponding limit cycle is responsible for periodic outbursts, and the timescales are much shorter than the standard viscous timescale due to the narrowness of the unstable radial zone. Our toy model gives quantitative predictions and works well for multiple outbursts like those observed in NGC 1566, NGC 4151, NGC 5548 and GSN 069, although the shapes of the outbursts are not yet well modeled, and further development of the model is necessary.

For about 500 intermediate-redshift sources ($0.04<z<0.4$), whose radio flux densities at $1.4\,{\rm GHz}$ are larger than 10 mJy, we performed additional observations at 4.85 GHz and 10.45 GHz using 100-m Effelsberg telescope. Our radio-optical galaxies are located preferentially in the composite and AGN spectral classes in the narrow line optical diagnostic diagrams (ODD). In the analysis, we focused on the distribution of radio spectral indices of radio synchrotron power-law profiles, $S_{\nu}\propto \nu^{+\alpha}$, in the ODDs. Using different analysis techniques, both observationally motivated and machine-learning based, we found three distinct groups--clusters in the radio loudness, [OIII]/H$\beta$ ratio, and spectral index volume: (1) sources with steep radio spectral index, large radio loudness and large [OIII]/H$\beta$ ratio; (2) sources with flat radio spectral index, intermediate radio loudness and lower [OIII]/H$\beta$ ratio; (3) sources with inverted radio spectral index, low radio loudness and low [OIII]/H$\beta$. The groups (1), (2), and (3) are located along the Seyfert-LINER spectral classes towards lower ionization ratios in the ODDs and hence can represent different activity cycles/accretion modes of AGNs, which could be in some cases associated with different merger stages.

Recently, near-infrared GRAVITY@ESO observations at $2.2\,\mu{\rm m}$ have announced the detection of three bright "flares" in the vicinity of the Galactic center supermassive black hole (SMBH) that exhibited orbital motion at a distance of about $6 - 11$ gravitational radii from an $\sim 4\times 10^6\, M_{\odot}$ black hole. There are indications of the presence of a large-scale, organized component of the magnetic field at the Galactic center. Electromagnetic effects on the flare dynamics were previously not taken into account despite the relativistic motion of a plasma in magnetic field leading to the charge separation and nonnegligible net charge density in the plasma. Applying various approaches, we find the net charge number density of the flare components of the order of $10^{-3} - 10^{-4}$ cm$^{-3}$, while the particles' total number density is of the order of $10^{6} - 10^{8}$ cm$^{-3}$. However, even such a tiny excess of charged particles in the quasi-neutral plasma can significantly affect the dynamics of flare components, which can then lead to the degeneracy in the measurements of spin of the SMBH. Analyzing the dynamics of recent flares in the case of the rapidly rotating black hole, we also constrain the inclination angle between the magnetic field and spin axis to $\alpha < 50^{\circ}$, as for larger angles, the motion of the hot spot is strongly chaotic.

In Failed Radiatively Accelerated Dusty Outflow (FRADO) model which provides the source of material above the accretion disk (AD) as an option to explain the formation mechanism of Broad Line Region (BLR) in AGNs, the BLR inner radius ($\rm{BLR}_{in}$ hereafter) is set by the condition that the dust evaporates immediately upon departure from the AD surface. On the other hand, the location of BLR clouds obtained observationaly via reverberation mapping shows some scaling with the source luminosity, so-called RL relation. We assume $\rm{BLR}_{in}$ to be the location of BLR clouds, then using a realistic expression for the radiation pressure of an AD, and having included the proper values of dust opacity, and shielding effect as well, we report our numerical results on calculation of $\rm{BLR}_{in}$ based on FRADO model. We investigate how it scales with monochromatic luminosity at 5100 angstrom for a grid of blackhole masses and Eddington ratios to compare along with the FRADO analytically predicted RL directly to observational data.

We study the multi-phase medium in the innermost regions of the Centaurus A (Cen A) galaxy. Combined high-resolution observations with the ALMA and Chandra observatories indicate that the hot X-ray emitting plasma coexists with the warm and cold medium in Cen A. This complex environment is a source of CO lines with great impact for its diagnostics. We show the images from the two mentioned instruments covering the nuclear region (diameter of 10'' or ~184 pc), and we study the conditions for plasma thermal equilibrium and possible coexistence of cool clouds embedded within the hot X-ray emitting gas. Further, we demonstrate that the multi-phase medium originates naturally by the thermal instability arising due to the interaction of the high-energy radiation field from the nucleus with the ambient gas and dust. We show that cold gas clouds can coexist in the mutual contact with hot plasma, but even colder dusty molecular clouds have to be distanced by several hundred pc from the hot region. Finally, we propose a 3-D model of the appearance of the hot plasma and the CO line-emitting regions consistent with the Chandra image and we derive the integrated emissivity in specific molecular lines observed by ALMA from this model. To reproduce the observed images and the CO line luminosity the dusty shell has to be 419 pc thick and located at 1020 pc from the centre.

We present the monitoring of the AGN continuum and MgII broad line emission for the quasar HE 0413-4031 ($z=1.38$) based on the six-year monitoring by the South African Large Telescope (SALT). We managed to estimate a time-delay of $302.6^{+28.7}_{-33.1}$ days in the rest frame of the source using seven different methods: interpolated cross-correlation function (ICCF), discrete correlation function (DCF), $z$-transformed DCF, JAVELIN, two estimators of data regularity (Von Neumann, Bartels), and $\chi^2$ method. This time-delay is below the value expected from the standard radius-luminosity relation. However, based on the monochromatic luminosity of the source and the SED modelling, we interpret this departure as the shortening of the time-delay due to the higher accretion rate of the source, with the inferred Eddington ratio of $\sim 0.4$. The MgII line luminosity of HE 0413-4031 responds to the continuum variability as $L_{\rm line}\propto L_{\rm cont}^{0.43\pm 0.10}$, which is consistent with the light-travel distance of the location of MgII emission at $R_{\rm out} \sim 10^{18}\,{\rm cm}$. Using the data of 10 other quasars, we confirm the radius-luminosity relation for broad MgII line, which was previously determined for broad H$\beta$ line for lower-redshift sources. In addition, we detect a general departure of higher-accreting quasars from this relation in analogy to H$\beta$ sample. After the accretion-rate correction of the light-travel distance, the MgII-based radius-luminosity relation has a small scatter of only $0.10$ dex.

Optical Fe II emission is a strong feature of quasar spectra originating in the broad-line regions (BLRs). The difficulty in understanding the complex Fe II pseudo-continuum has led us in search of other reliable, simpler ionic species such as Ca II. In this first part of the series, we confirm the strong correlation between the strengths of two emission features, the optical Fe II and the NIR Ca II, both from observations and photoionization modelling. With the inclusion of an up-to-date compilation of observations with both optical Fe II and NIR Ca II measurements, we span a wider and more extended parameter space and confirm the common origin of these two species with our photoionization models using CLOUDY. Taking into account the effect of dust into our modelling, we retrieve the log U - log n parameter space as a function of the strengths of Fe II and Ca II emission. Incorporating a broad range of metallicity and extended cloud columns we reinstate this correlation and explain the emission from the strong Fe II emitting Narrow-Line Seyfert 1 galaxies (NLS1s).

The effective size of Broad Line Region (BLR), so-called the BLR radius, in galaxies with active galactic nuclei (AGN) scales with the source luminosity. Therefore by determining this location either observationally through reverberation mapping or theoretically, one can use AGNs as an interesting laboratory to test cosmological models. In this article we focus on the theoretical side of BLR based on the Failed Radiatively Accelerated Dusty Outflow (FRADO) model. By simulating the dynamics of matter in BLR through a realistic model of radiation of accretion disk (AD) including the shielding effect, as well as incorporating the proper values of dust opacities, we investigate how the radial extension and geometrical height of the BLR depends on the Eddington ratio [and blackhole mass], and modeling of shielding effect. We show that assuming a range of Eddington ratios and shielding we are able to explain the measured time-delays in a sample of reverberation-measured AGNs.

The main sequence offers a method for the systematization of quasar spectral properties. Extreme FeII emitters (or extreme Population A, xA) are believed to be sources accreting matter at very high rates. They are easily identifiable along the quasar main sequence, in large spectroscopic surveys over a broad redshift range. The very high accretion rate makes it possible that massive black holes hosted in xA quasars radiate at a stable, extreme luminosity-to-mass ratio. After reviewing the basic interpretation of the main sequence, we report on the possibility of identifying virial broadening estimators from low-ionization line widths, and provide evidence of the conceptual validity of redshift-independent luminosities based on virial broadening for a known luminosity-to-mass ratio.

A soft X-ray excess above the 2-10 keV power law extrapolation is generally observed in AGN X-ray spectra. Presently there are two competitive models to explain it: blurred ionized reflection and warm Comptonisation. In the latter case, observations suggest a corona temperature $\sim$ 1 keV and a corona optical depth $\sim$ 10. Moreover, radiative constraints from spectral fits with Comptonisation models suggest that most of the accretion power should be released in the warm corona. The disk below is basically non-dissipative, radiating only the reprocessed emission from the corona. The true radiative properties of such a warm and optically thick plasma are not well-known, however. For instance, the importance of the Comptonisation process, the potential presence of strong absorption/emission features or the spectral shape of the output spectrum have been studied only very recently. We present in this paper simulations of warm and optically thick coronae using the TITAN radiative transfer code coupled with the NOAR Monte-Carlo code, the latter fully accounting for Compton scattering of continuum and lines. Illumination from above by a hard X-ray emission and from below by an optically thick accretion disk is taken into account as well as (uniform) internal heating. Our simulations show that for a large part of the parameter space, the warm corona with sufficient internal mechanical heating is dominated by Compton cooling and neither strong absorption nor emission lines are present in the outgoing spectra. In a smaller part of the parameter space, the calculated emission agrees with the spectral shape of the observed soft X-ray excess. Remarkably, this also corresponds to the conditions of radiative equilibrium of an extended warm corona covering almost entirely a non-dissipative accretion disk. These results confirm the warm Comptonisation as a valuable model for the soft X-ray excess.

We address the effect of orientation of the accretion disk plane and the geometry of the broad-line region (BLR) in the context of understanding the distribution of quasars along their Main Sequence. We utilize the photoionization code CLOUDY to model the BLR, incorporating the `un-constant' virial factor. We show the preliminary results of the analysis to highlight the co-dependence of the Eigenvector 1 parameter, $\mathrm{R_{FeII}}$ on the broad H$\beta$ FWHM (i.e. the line dispersion) and the inclination angle ($\theta$), assuming fixed values for the Eddington ratio ($\mathrm{L_{bol}}/\mathrm{L_{Edd}}$), black hole mass ($\mathrm{M_{BH}}$), spectral energy distribution (SED) shape, cloud density (n$\rm{_{H}}$) and composition.

The interpretation of the main sequence of quasars has become a frontier subject in the last years. The consider the effect of a highly flattened, axially symmetric geometry for the broad-line region (BLR) on the parameters related to the distribution of quasars along their Main Sequence. We utilize the photoionization code CLOUDY to model the BLR, assuming `un-constant' virial factor with a strong dependence on the viewing angle. We show the preliminary results of the analysis to highlight the co-dependence of the Eigenvector 1 parameter, R$_{\mathrm{FeII}}$ on the broad H$\beta$ FWHM (i.e. the line dispersion) and the inclination angle ($\theta$), assuming fixed values for the Eddington ratio ($\mathrm{L_{bol}}/\mathrm{L_{Edd}}$), black hole mass ($\mathrm{M_{BH}}$) and spectral energy distribution (SED) shape. We consider four cases with changing cloud density (n$\rm{_{H}}$) and composition. Understanding the FeII emitting region is crucial as this knowledge can be extended to the use of quasars as distance indicators for Cosmology.

In the near future, the arrival of new surveys promises a significant increment in the number of quasars (QSO) at large redshift ranges. This will help to constrain the dark energy models using quasars. The Large Synoptic Survey Telescope (LSST) will cover over 10 million QSO in six photometric bands during its 10-year run. QSO will be monitored and subsequently analysed using the photometric reverberation mapping technique. Oncoming LSST data will decrease the uncertainties in the dark energy determination using reverberating--mapped sources, particularly at high redshifts. We present the first steps in modelling of light curves for H$\beta$ and MgII and discuss the quasar selection in the context of photometric reverberation mapping with LSST. With the onset of the LSST era, we expect a huge rise in the overall quasar counts and redshift range covered ($z\leq7.0$), which will provide a better constraint of AGN properties with cosmological purposes.

A quasar main sequence is widely believed to reveal itself through objects represented in a plane spanned by two parameters: the full-width at half-maximum (FWHM) of H$\beta{}$, and the ratio of Fe II to H$\beta{}$ equivalent width. This sequence is related to the application to quasar properties of principal component analysis (PCA), which reveals that the main axis of variance (Eigenvector 1) is co-directional with a strong anti-correlation between these two measurements. We aim to determine whether the dominance of two Eigenvectors, originally discovered over two decades ago, is replicated in newer high-quality quasar samples. If so, we aim to test if a non-linear approach is an improvement on the linear PCA method by finding two new parameters which represent a more accurate projection of the variances than the Eigenvectors recovered from PCA. We selected quasars from the XSHOOTER archive and a major quasar catalog to build high quality samples. These samples were tested with PCA. We find that the new high-quality samples do indeed have two dominant Eigenvectors as originally discovered. Subsequently we find that the fitting of a non-linear decay curve to the main sequence allows a new plane spanned by linearly independent axes to be defined, based on the distance along the decay curve as the main axis and the distance of each quasar datapoint from the curve as the secondary axis respectively. The results show that it is possible to define a new plane based on the quasar main sequence which accounts for the majority of the variance. The most likely candidate for the new main axis is an anti-correlation with black hole mass. In this case the secondary axis likely represents luminosity. However, given the results of previous studies, inclination angle likely plays a role in H$\beta{}$ width.

In the near future, new surveys promise a significant increase in the number of quasars (QSO) at large redshifts. This will help to constrain the dark energy models using quasars. The Large Synoptic Survey Telescope (LSST) will cover over 10 million QSO in six photometric bands during its 10-year run. QSO will be monitored and subsequently analyzed using the photometric reverberation mapping (RM) technique. In low-redshift quasars, the combination of reverberation-mapped and spectroscopic results have provided important progress. However, there are still some facts which have to be taken into account for future results. It has been found that super-Eddington sources show time delays shorter than the expected from the well-known Radius-Luminosity ($R\mathrm{_{H\beta}}-L_{5100}$) relation. Using a sample of 117 H$\beta$ reverberation-mapped AGN with $0.02<z<0.9$, we propose a correction by the accretion rate effect recovering the classical $R\mathrm{_{H\beta}}-L_{5100}$ relation. We determined the cosmological constants, which are in agreement with $\Lambda$-Cold Dark Matter model within 2$\sigma$ confidence level, which is still not suitable for testing possible departures from the standard model. Upcoming LSST data will decrease the uncertainties in the dark energy determination using reverberation-mapped sources, particularly at high redshifts. We show the first steps in the modeling of the expected light curves for H$\beta$ and MgII.

Broad Emission Lines are the most characteristic features of Active Galaxies, but the mechanism of creating a medium able to emit these intense lines is not quite clear. Observations clearly indicate that the motion of the material is predominantly Keplerian, with traces of inflow and clear signatures of outflow, but this still does not point out whether the lines partially come from the disk surface, or exclusively from the circumnuclear material, and whether this material originates from the disk as a wind, or comes, at least partially, from outside. I review the basic scenarios for the formation of the Broad Line Region (BLR), and the recent progress in modelling the physical conditions in the emitting medium. The current state is the the outer radius of the BLR is fixed by the dust sublimation temperature in the medium exposed to the irradiation from the central source, the inner radius is likely fixed by the dust sublimation temperature in the atmosphere of the non-illuminated accretion disk, and the local cloud density is a universal number fixed by the radiation pressure confinement. The time-dependent aspects of the BLR formation, however, still wait for serious modelling effort matching the quality of the observational data.

Quasars have been proposed as a new class of standard candles analogous to Supernovae, as their large redshift range and high luminosities make them excellent candidates. Reverberation mapping (RM) method enables to estimate the distance to the source from the time delay measurement of the emission lines with respect to the continuum, since the time delay depends on the absolute luminosity of the source. This radius-luminosity (RL) relation exhibits a low scatter and it was thus proposed for cosmological purposes. However, in recent years the increase in the studied sample, and in particular the inclusion of highly accreting QSO has increased the dispersion in the RL relation, with many objects showing time delays shorter than the expected. Using H$\beta$ RM measurements for 117 sources with $0.002<z<0.9$ and $41.5<\mathrm{log}~L_{5100}<45.9$, we find a correction for the time delay based on the dimensionless accretion rate considering the virial factor anti-correlated with the FWHM of H$\beta$. This correction decreases the scattering in the accretion parameters compared with the typical value used, which is directly reflected by suppressing the RL relation dispersion. We also confirm the anti-correlation between the excess of variability and the accretion parameters. With this correction we are able to build the Hubble diagram and estimate the cosmological constants $\Omega_m$ and $\Omega_\Lambda$, which are consistent with the standard cosmological model at 2$\sigma$ confidence level. Therefore, RM results can be used to constrain cosmological models in the future.

We address the effect of the viewing angle of the accretion disk plane and the geometry of the broad-line region (BLR) with the goal of interpreting the distribution of quasars along the main sequence (MS). We utilize photoionization code CLOUDY to model the BLR FeII emission, incorporating the grossly underestimated role of the form factor (f). We recover the dependence of the strength of the FeII emission in the optical (R$_{\rm{FeII}}$) on L$_{\rm{bol}}$/L$_{\rm{Edd}}$ ratio and related observational trends - as a function of the spectral energy distribution (SED) shape, cloud density, composition and intra-cloud dynamics, assumed following prior observational constraints. With this approach, we are now able to explain the diversity of quasars and the change of the quasar properties along the Main Sequence (MS). Our approach also explains the rarity of the highest FeII emitters known as the extreme xA sources and can be used as a predictive tool in future reverberation mapping studies of Type-1 AGNs. This approach further justifies the use of quasars as `cosmological probes'.

We address the effect of orientation of the accretion disk plane and the geometry of the broad-line region (BLR) as part of an effort to understand the distribution of quasars in the optical plane of the quasar main sequence. We utilize the photoionization code CLOUDY to model the BLR incorporating the grossly underestimated form factor ($f$). Treating the aspect of viewing angle appropriately, we confirm the dependence of the $R_{\mathrm{FeII}}$ sequence on Eddington ratio and on the related observational trends - as a function of the SED shape, cloud density and composition, verified from prior observations. Sources with $R_{\mathrm{FeII}}$ in the range 1 -- 2 (about 10\% of all quasars, the so-called extreme Population A [xA] quasars) are explained as sources of high, and possibly extreme Eddington ratio along the $R_{\mathrm{FeII}}$ sequence. This result has important implication for the exploitation of xA sources as distance indicators for Cosmology. $\mathrm{FeII}$ emitters with $R_{\mathrm{FeII}} > 2$ are very rare (<1\% of all type 1 quasars). Our approach also explains the rarity of these highest $\mathrm{FeII}$ emitters as extreme xA sources and constrains the viewing angle ranges with increasing H$\beta$ FWHM.

By applying different statistically robust methods, we analyze the time-lag between the continuum and ionized line-emission (Mg II line) light curves for the distant bright quasar CTS C$30.10$ (redshift $z \sim 0.9$). The data were obtained by the 10-meter SALT telescope in South Africa. In detail, we demonstrate the application of several methods using interpolated cross-correlation function (ICCF), discrete correlation function (DCF), z-transformed discrete correlation function (zDCF), von Neumann estimator, and the JAVELIN code package. In particular, we discuss the uncertainties of these methods. In conclusion, we find that the quasar lies on the broad-line region (BLR) size -- monochromatic luminosity power-law scaling, $R_{\rm BLR}\propto L_{5100}^{1/2}$, which was already confirmed for low-redshift sources. In case the BLR size-luminosity relation holds for other distant sources, quasars could be used for probing cosmological models as ``standard candles'' complementary to supernovae Ia.

The Absorption Measure Distribution (AMD) in the X-ray outflows of Seyfert active galactic nuclei (AGN) describes the distribution of absorbing column density as a function of ionization parameter. Up to now, AMD has been measured only for seven objects with high-resolution X-ray data that contain absorption lines from ionized heavy elements. Even though the number of measured AMDs is not large, they display a universal broad shape containing a prominent dip, for which the absorbing column drops by around two orders of magnitude. In this paper, we tested a range of photoionization models against the overall shape of the AMD as observed in Seyferts. In particular, we demonstrate that the shape of the AMD depends both on the spectral energy distribution (SED) of radiation which enters the outflow, and the density of the warm absorber (WA). The model that best reproduces the observed shape of the AMD is one wherein the gas density of the WA is of the order of $10^{12}$ cm$^{-3}$, irradiated by an SED whose optical/UV luminosity is 100 times higher than the X-ray luminosity. When the cloud density is higher than $\sim 10^{11}$ cm$^{-3}$, free-free heating dominates the entire absorber, and only one instability zone occurs, which is in agreement with observations.

Slim accretion disks idea emerged over 30 years ago as an answer to several unsolved problems. Since that time there was a tremendous increase in the amount of observational data where this model applies. However, many critical issues on the theoretical side remain unsolved, as they are inherently difficult. This is the issue of the disk stability under the radiation pressure, the role of the magnetic field in the energy transfer inside the disk and the formation (or not) of a warm corona, and outflows. Thus the progress has to be done both through further developments of the model and through the careful comparison to the observational data.

Bright active galaxies show a range of properties but many of these properties are correlated which has led to the concept of the Quasar Main Sequence. We test whether our current understanding of the quasar structure allows to reproduce the pattern observed in the optical plane formed by the kinematic line width of H$\beta$ and the relative importance of the Fe II optical emission. We performed simulations of the H$\beta$ and Fe II production using the code CLOUDY and well justified assumptions about the broad band spectra, distance of the emission line region, and the cloud properties. We show that the presence of the warm corona is an important element of the broad band spectrum which decreases the dependence of the relative Fe II emissivity on the Eddington ratio, and allows to reproduce the rare cases of the particularly strong Fe II emitters. Results are sensitive to the adopted cloud distance, and strong Fe II emission can be obtain either by adopting strongly super-solar metallicity, or much shorter distance than traditionally obtained from reverberation mapping. We modeled in a similar way the UV plane defined by the Mg II line and Fe II UV pseudo-continuum, but here our approach is less successful, in general overproducing the Fe II strength. We found that the Fe II optical and UV emissivity depend in a different way on the turbulent velocity and metallicity, and the best extension of the model in order to cover both planes is to allow very large turbulent velocities in the Broad Line Region clouds.

Changing-look phenomenon observed now in a growing number of active galaxies challenges our understanding of the accretion process close to a black hole. We propose a simple explanation for periodic outbursts in sources operating at a few per cent of the Eddington limit. The mechanism is based on two relatively well understood phenomena: radiation pressure instability and formation of the inner optically thin Advection-Dominated Accretion Flow. The limit cycle behaviour takes place in a relatively narrow transition zone between the standard disk and optically thin flow. Large changes in the cold disk are due to the irradiation by the hot flow with accretion rate strongly varying during the cycle. The model gives quantitative predictions and works well for multiple outbursts of NGC 1566.

In this contribution, we summarize our results concerning the observational constraints on the electric charge associated with the Galactic centre black hole - Sgr A*. According to the no-hair theorem, every astrophysical black hole, including supermassive black holes, is characterized by at most three classical, externally observable parameters - mass, spin, and the electric charge. While the mass and the spin have routinely been measured by several methods, the electric charge has usually been neglected, based on the arguments of efficient discharge in astrophysical plasmas. From a theoretical point of view, the black hole can attain charge due to the mass imbalance between protons and electrons in fully ionized plasmas, which yields about $\sim 10^8\,{\rm C}$ for Sgr A*. The second, induction mechanism concerns rotating Kerr black holes embedded in an external magnetic field, which leads to electric field generation due to the twisting of magnetic field lines. This electric field can be associated with the induced Wald charge, for which we calculate the upper limit of $\sim 10^{15}\,{\rm C}$ for Sgr A*. Although the maximum theoretical limit of $\sim 10^{15}\,{\rm C}$ is still 12 orders of magnitude smaller than the extremal charge of Sgr A*, we analyse a few astrophysical consequences of having a black hole with a small charge in the Galactic centre. Two most prominent ones are the effect on the X-ray bremsstrahlung profile and the effect on the position of the innermost stable circular orbit.

Context: Giant radio galaxies (with projected linear size of radio structure larger than 0.7 Mpc) are very rare and unusual objects. Only $\sim$5% of extended radio sources reach such sizes. Understanding of the processes responsible for their large sizes is crucial to further our knowledge about the radio source's evolution. Aims: We check the hypothesis that giants become extremely large due to the specific history of their host galaxy formation, as well as in the context of the cluster or group of galaxies where they evolve. Therefore we study the star formation histories in their host galaxies and in galaxies located in their neighbourhood. Methods: We studied 41 giant-size radio galaxies as well as galaxies located within a radius of 5 Mpc around giants to verify whether the external conditions of the intergalactic medium somehow influence the internal evolution of galaxies in the group/cluster. We compared the results with a control sample of smaller-sized Fanaroff--Riley type II radio galaxies and their neighbouring galaxies. We fit stellar continua in all galaxy spectra using the spectral synthesis code STARLIGHT and provide statistical analysis of the results. Results: We find that hosts of giant radio galaxies have a larger amount of intermediate age stellar populations compared with smaller-sized FRII radio sources. The same result is also visible when we compare neighbouring galaxies located up to 1.5 Mpc around giants and FRIIs. This may be evidence that star formation in groups with giants was triggered due to global processes occurring in the ambient intergalactic medium. These processes may also contribute to mechanisms responsible for the extremely large sizes of giants.

We write in response to the call from the 2020 Decadal Survey to submit white papers illustrating the most pressing scientific questions in astrophysics for the coming decade. We propose exploration as the central question for the Decadal Committee's discussions.The history of astronomy shows that paradigm changing discoveries are not driven by well formulated scientific questions, based on the knowledge of the time. They were instead the result of the increase in discovery space fostered by new telescopes and instruments. An additional tool for increasing the discovery space is provided by the analysis and mining of the increasingly larger amount of archival data available to astronomers. Revolutionary observing facilities, and the state of the art astronomy archives needed to support these facilities, will open up the universe to new discovery. Here we focus on exploration for compact objects and multi messenger science. This white paper includes science examples of the power of the discovery approach, encompassing all the areas of astrophysics covered by the 2020 Decadal Survey.

We review the current status of the understanding how the accretion onto the central black hole proceeds in Active Galaxies. Standard accretion disk is a key element in all relatively bright active galaxies like Seyferts and quasars, although it is not present in very low luminosity sources, like Sgr A*. However, the standard disk does not explain the broad band spectrum, so the disk has to be supplemented with a number of additional components, and our deeper understanding of these components is still far from being complete. These additional elements are: compact hard X-ray corona, inner hot flow, warm corona, disk wind and the Broad Line Region, and finally dusty/molecular torus. All these elements seem to be needed in various proportions, depending predominantly on the Eddington ratio of a given source. These elements also interact with each other which is not yet fully taken into account.

We report 6 yr monitoring of a distant bright quasar CTS C30.10 (z = 0.90052) with the Southern African Large Telescope (SALT). We measured the rest-frame time-lag of $562\pm 2$ days between the continuum variations and the response of the Mg II emission line, using the Javelin approach. More conservative approach, based on five different methods, imply the time delay of $564^{+109}_{-71}$ days. This time delay, combined with other available measurements of Mg II line delay, mostly for lower redshift sources, shows that the Mg II line reverberation implies a radius-luminosity relation very similar to the one based on a more frequently studied H$\beta$ line.

The work initially started as a test to retrace the Shen & Ho (2014) Quasar Main Sequence diagram where they claimed that the parameter RFeII, which defines the Eigenvector 1 (EV1) is driven by the Eddington ratio alone. We subsequently construct a refined (error and redshift limited) sample from the original Shen et al. 2011 QSO catalogue. Based on our hypothesis - the main driver of the Quasar Main Sequence is the maximum of the accretion disk temperature defined by the Big Blue Bump on the Spectral Energy Distribution. We select the four extreme sources that have RFeII $\geq$ 4.0 and use the SED modelling code CIGALE to fit the multi-band photometric data for these sources. We also perform detailed spectral fitting including the FeII pseudo-continuum to estimate and compare the value of RFeII for them. We show the metallicity dependent FeII strength in the context of this study.

Most results of the reverberation monitoring of active galaxies showed a universal scaling of the time delay of the Hbeta emission region with the monochromatic flux at 5100 A, with very small dipersion. Such a scaling favored the dust-based formation mechanism of the Broad Line Region (BLR). Recent reverberation measurements showed that actually a significant fraction of objects exhibits horter lags than the previously found scaling. Here we demonstrate that these shorter lags can be explained by the old concept of scaling of the BLR size with the ionization parameter. Assuming a universal value of this parameter and universal value of the cloud density reproduces the distribution of observational points in the time delay vs. monochromatic flux plane, provided that a range of black hole spins is allowed. However, a confirmation of the new measurements for low/moderate Eddington ratio sources is strongly needed before the dust-based origin of the BLR can be excluded.

Active galaxies form a clear pattern in the optical plane showing the correlation between the Full Width Half Maximum (FWHM) of the H$\beta$ line and the ratio of the Equivalent Width (EW) of the optical FeII emission and the EW(H$\beta$). This pattern is frequently refereed to as the quasar main sequence. In this paper we study the UV plane showing FWHM of the MgII line against the ratio of EW of UV FeII emission to EW(MgII). We show that the UV plane trends are different, with the underlying strong correlation between FWHM(MgII) and EW(MgII). This correlation is entirely driven by the choice of the continuum used to measure EW(MgII). If we correct for the strong curvature in the UV continuum related to the character of emission expected from accretion disk, the behavior of the FWHM vs. EW for MgII becomes similar to the behaviour of H$\beta$. Such a similarity is expected since both lines belong to the Low Ionization group of emission lines and come from a similar region. Such a similarity and weak coupling between FWHM and EW of the corresponding line is then also consistent with the Broad Line Region model based on the presence of dust in the accretion disk atmosphere.

The concept of the quasar main sequence is very attractive since it stresses correlations between various parameters and implies the underlying simplicity. In the optical plane defined by the width of the H\{beta} line and the ratio of the equivalent width of the Fe II to H\{beta} observed objects form a characteristic pattern. In this paper, we use a physically motivated model to explain the distribution of quasars in the optical plane. Continuum is modelled as an accretion disk with a hard X-ray power law uniquely tight to the disk at the basis of observational scaling, and the Broad Line Region distance is determined also from observational scaling. We perform the computations of the FeII and H\{beta} line production with the code CLOUDY. We have only six free parameters for an individual source: maximum temperature of the accretion disk, Eddington ratio, cloud density, cloud column density, microturbulence, and iron abundance, and only the last four remain as global parameters in our modelling of the whole sequence. Our theoretically computed points cover well the optical plane part populated with the observed quasars, particularly if we allow for super-Solar abundance of heavy elements. Explanation of the exceptionally strong Fe II emitter requires a stronger contribution from the dark sides of the clouds. Analyzing the way how our model covers the optical plane we conclude that there is no single simple driver behind the sequence, as neither the Eddington ratio nor broad band spectrum shape plays the dominant role. Also, the role of the viewing angle in providing the dispersion of the quasar main sequence is apparently not as strong as expected.

Narrow-line Seyfert 1 galaxies are defined on the basis of their line widths, and they are generally considered to be high Eddington ratio sources. But in the context of the Quasar Main Sequence, high Eddington rate sources are those which have weak [O III] lines and strong Fe II lines. There is an overlap between the two populations, but they are not identical. Thus these two selection criteria give a different view on which objects are actually high Eddington ratio sources. We discuss this issue in the context of the broad band spectral energy density, emission line shape modeling, Fe II pseudo-continuum strength, and the level of X-ray variability. We also discuss the issue of the viewing angle and the insight one can gain from spectropolarimetric observations. We conclude that there is no single driver behind the Quasar Main Sequence, and the Eddington ratio of the source cannot be determined from the source location on the optical plane alone. On the other hand, an expected range of AGN parameters combined with a simple model of the Fe II production represent well the observed coverage pattern of the plane, with not much effect needed from the dispersion due to the viewing angle.

The work initially started as a test to retrace the Shen & Ho (2014) Quasar Main Sequence diagram where they (and references therein) claimed that the parameter RFeII, which defines the Eigenvector 1 (EV1) is driven by the Eddington ratio alone. We subsequently construct a refined (error and redshift limited) sample from the original Shen et. al (2011) QSO catalogue. Based on our hypothesis - the main driver of the Quasar Main Sequence is the maximum of the accretion disk temperature defined by the Big Blue Bump on the Spectral Energy Distribution. We select five extreme sources that have RFeII >= 4.0 and use the SED modelling code CIGALE to fit the multi-band photometric data for these sources. Incorporating the prescription for broad emission-line species dependent virial factors, we derive the black hole masses for the entire Shen+11 catalogue and compare their distribution with the black hole masses in the original catalogue. Here, we show the results of our analysis for one of the 5 extreme sources, SDSSJ082358.30+213545.2. We also show the detailed modelling, including the Fe II pseudo-continuum to estimate and compare the value of RFeII for this object.

We investigate the relationships between Mg ii {\lambda}2798 emission line properties, as well as between these properties and inclination angle and Fe ii strength, in a lobe-dominated quasar sample. We find no correlation between Mg ii line width and inclination, unlike previous studies of the general quasar population. This suggests that the Mg ii emission region in these objects is not of a thin disk geometry, however the average equivalent width of the line negates a spherical alternative. A positive correlation between Mg ii equivalent width and inclination cannot be ruled out, meaning there is no strong evidence that Mg ii emission is anisotropic. Since thin disk emission would be highly directional, the geometric implications of these two findings are compatible. The lack of line width correlation with inclination may also indicate that Mg ii is useful for estimating black hole masses in lobe-dominated quasar samples, since it is unbiased by orientation. Some quasars in our sample have almost edge-on viewing angles and therefore cannot have a smooth toroidal obscurer co-planar with the accretion disk. Alternatives may be a distorted dusty disk or a clumpy obscurer. This could result from the sample selection bias towards high inclination objects, rather than intrinsic differences between lobe-dominated and typical quasars. Five objects have visible [O iii] allowing equivalent width calculation, revealing it to be higher than in typical quasars. Since these objects are of high inclination, this finding supports the positive correlation between [O iii] equivalent width and inclination found in a previous study.

The requirement of intermediate line component in the recently observed spectra of several AGNs points to possibility of the existence of a physically separate region between broad line region (BLR) and narrow line region (NLR). In this paper we explore the emission from intermediate line region (ILR) by using the photoionization simulations of the gas clouds distributed radially from the AGN center. The gas clouds span distances typical for BLR, ILR and NLR, and the appearance of dust at the sublimation radius is fully taken into account in our model. Single cloud structure is calculated under the assumption of the constant pressure. We show that the slope of the power law cloud density radial profile does not affect the existence of ILR in major types of AGN. We found that the low ionization iron line, Fe~II, appears to be highly sensitive for the presence of dust and therefore becomes potential tracer of dust content in line emitting regions. We show that the use of disk-like cloud density profile computed at the upper part of the accretion disc atmosphere reproduces the observed properties of the line emissivities. In particular, the distance of H${\beta}$ line inferred from our model agrees with that obtained from the reverberation mapping studies in Sy1 galaxy NGC 5548.

Eigenvector 1 (EV1) is the formal parameter which allows the introduction of some order in the properties of the unobscured type 1 active galaxies. We aim to understand the nature of this parameter by analyzing the most extreme examples of quasars with the highest possible values of the corresponding eigenvalues $R_{Fe}$. We selected the appropriate sources from the Sloan Digital Sky Survey (SDSS) and performed detailed modeling, including various templates for the Fe II pseudo-continuum and the starlight contribution to the spectrum. Out of 27 sources with $R_{Fe}$ larger than 1.3 and with the measurement errors smaller than 20\% selected from the SDSS quasar catalog, only six sources were confirmed to have a high value of $R_{Fe}$, defined as being above 1.3. All other sources have $an R_{Fe}$ of approximately 1. Three of the high $R_{Fe}$ objects have a very narrow H$\beta$ line, below 2100 km s$^{-1}$ but three sources have broad lines, above 4500 km s$^{-1}$, that do not seem to form a uniform group, differing considerably in black hole mass and Eddington ratio; they simply have a very similar EW([OIII]5007) line. Therefore, the interpretation of the EV1 remains an open issue.

The formal division of the distance indicators into primary and secondary leads to difficulties in description of methods which can actually be used in two ways: with, and without the support of the other methods for scaling. Thus instead of concentrating on the scaling requirement we concentrate on all methods of distance determination to extragalactic sources which are designated, at least formally, to use for individual sources. Among those, the Supernovae Ia is clearly the leader due to its enormous success in determination of the expansion rate of the Universe. However, new methods are rapidly developing, and there is also a progress in more traditional methods. We give a general overview of the methods but we mostly concentrate on the most recent developments in each field, and future expectations.

Eigenvector 1 (EV1) was found to be the dominant component behind the significant correlations for the measured parameters in quasar spectra (Boroson & Green, 1992). The parameter R$_{\mathrm{FeII}}$, which strongly correlates to the EV1, is the ${\mathrm{FeII}}$ strength, defined to be the ratio of the equivalent width of ${\mathrm{FeII}}$ to the equivalent width of ${\mathrm{H\beta}}$. This allows to construct a quasar main sequence analogous to the stellar properties driven HR diagram (Sulentic et al. 2001). We try to find the main driver behind the EV1 among the basic (theoretically motivated) parameters of an active nucleus (Eddington ratio, black hole mass, accretion rate, spin, and viewing angle). Based on theoretical modeling using the photoionization code CLOUDY (Ferland et al. 2013), we test the hypothesis that the physical driver of EV1 is the maximum of the accretion disk temperature ($\mathrm{T_{BBB}}$), reflected in the shape of the spectral energy distribution (SED). We have assumed that both H$\mathrm{\beta}$ and Fe${\mathrm{II}}$ emission come from the Broad Line Region represented as a constant density cloud in a plane-parallel geometry. We test the effect of changing Eddington ratio on the $\mathrm{R_{FeII} - T_{BBB}}$ trends with varying mean hydrogen densities. We also test the effect of adding microturbulence that affect the line intensities on the overall $\mathrm{R_{FeII} - T_{BBB}}$ picture.

The most common observed feature of the X-ray outflows in active galactic nuclei (AGN) is their broad ionization distribution spanning up to $\sim$ 5 orders of magnitude in ionization levels. This feature is quantified in terms of absorption measure distribution (AMD), defined as the distribution of column density with the ionization parameter. Recently, the photoionization models with constant pressure assumption are shown to well reproduce the observed shape of AMD. However, there exist inconsistencies in the normalization and the position of discontinuities presented in the AMD shape between the observation and model. In this work, we show that AMD normalization differs by one order of magnitude depending on the shape of ionizing spectral energy distribution (SED).

Quasars are among the most luminous sources characterized by their broad band spectra ranging from radio through optical to X-ray band, with numerous emission and absorption features. Using the Principal Component Analysis (PCA), Boroson & Green (1992) were able to show significant correlations between the measured parameters. Among the significant correlations projected, the leading component, related to Eigenvector 1 (EV1) was dominated by the anti-correlation between the Fe${\mathrm{II}}$ optical emission and [OIII] line where the EV1 alone contained 30% of the total variance. This introduced a way to define a quasar main sequence, in close analogy to the stellar main sequence in the Hertzsprung-Russel (HR) diagram (Sulentic et. al 2001). Which of the basic theoretically motivated parameters of an active nucleus (Eddington ratio, black hole mass, accretion rate, spin, and viewing angle) is the main driver behind the EV1 yet remains to be answered. We currently limit ourselves to the optical waveband, and concentrate on theoretical modelling the Fe${\mathrm{II}}$ to H$\mathrm{\beta}$ ratio, and test the hypothesis that the physical driver of EV1 is the maximum of the accretion disk temperature, reflected in the shape of the spectral energy distribution (SED). We performed computations of the H$\mathrm{\beta}$ and optical Fe${\mathrm{II}}$ for a broad range of SED peak position using CLOUDY photoionisation code. We assumed that both H$\mathrm{\beta}$ and Fe${\mathrm{II}}$ emission come from the Broad Line Region represented as a constant density cloud in a plane-parallel geometry. We compare the results for two different approaches: (1) considering a fixed bolometric luminosity for the SED; (2) considering $\mathrm{L_{bol}/L_{Edd}}$ = 1.

Time-series photometry from space in the ultraviolet can be presently done with only a few platforms, none of which is able to provide wide-field long-term high-cadence photometry. We present a concept of UVSat, a twin space telescope which will be capable to perform this kind of photometry, filling an observational niche. The satellite will host two telescopes, one for observations in the ultraviolet, the other for observations in the optical band. We also briefly show what science can be done with UVSat.

Several observed spectral properties of quasars are believed to be influenced by quasar orientation. In this investigation we examine the effect of orientation on the Mg II line located at 2798 {\AA} in a sample of 36 radio-loud quasars, with orientation angles having been obtained in a previous study using radio observations. We find no significant relationship between orientation angle and either Mg II line full-width at half-maximum or equivalent width. The lack of correlation with inclination angle contradicts previous studies which also use radio data as a proxy for inclination angle and suggests the Mg II emission region does not occupy a disk-like geometry. The lack of correlation with Mg II equivalent width, however, is reported in at least one previous study. Although the significance is not very strong (86 percent), there is a possible negative relationship between inclination angle and Fe II strength which, if true, could explain the Fe II anti-correlation with [O III ] strength associated with Eigenvector 1. Interestingly, there are objects having almost edge-on inclinations while still exhibiting broad lines. This could be explained by a torus which is either clumpy (allowing sight lines to the central engine) or mis-aligned with the accretion disk.

Quasars are complex sources, characterized by broad band spectra from radio through optical to X-ray band, with numerous emission and absorption features. However, Boroson & Green (1992) used Principal Component Analysis (PCA), and with this analysis they were able to show significant correlations between the measured parameters. The leading component, related to Eigenvector 1 (EV1) was dominated by the anticorrelation between the Fe${\mathrm{II}}$ optical emission and [OIII] line and EV1 alone contained 30% of the total variance. It opened a way in defining a quasar main sequence, in close analogy to the stellar main sequence on the Hertzsprung-Russel (HR) diagram (Sulentic et al. 2001). The question still remains which of the basic theoretically motivated parameters of an active nucleus (Eddington ratio, black hole mass, accretion rate, spin, and viewing angle) is the main driver behind the EV1. Here we limit ourselves to the optical waveband, and concentrate on theoretical modelling the Fe${\mathrm{II}}$ to H$\mathrm{\beta}$ ratio, and we test the hypothesis that the physical driver of EV1 is the maximum of the accretion disk temperature, reflected in the shape of the spectral energy distribution (SED). We performed computations of the H$\mathrm{\beta}$ and optical Fe${\mathrm{II}}$ for a broad range of SED peak position using CLOUDY photoionisation code. We assumed that both H$\mathrm{\beta}$ and Fe${\mathrm{II}}$ emission come from the Broad Line Region represented as a constant density cloud in a plane-parallel geometry. We expected that a hotter disk continuum will lead to more efficient production of Fe${\mathrm{II}}$ but our computations show that the Fe${\mathrm{II}}$ to H$\mathrm{\beta}$ ratio actually drops with the rise of the disk temperature. Thus either hypothesis is incorrect, or approximations used in our paper for the description of the line emissivity is inadequate.

The mid-infrared (IR) range contains many spectral features associated with large molecules and dust grains such as polycyclic aromatic hydrocarbons (PAHs) and silicates. These are usually very strong compared to fine-structure gas lines, and thus valuable in studying the spectral properties of faint distant galaxies. In this paper, we evaluate the capability of low-resolution mid-IR spectroscopic surveys of galaxies that could be performed by SPICA. The surveys are designed to address the question how star formation and black hole accretion activities evolved over cosmic time through spectral diagnostics of the physical conditions of the interstellar/circumnuclear media in galaxies. On the basis of results obtained with Herschel far-IR photometric surveys of distant galaxies and Spitzer and AKARI near- to mid-IR spectroscopic observations of nearby galaxies, we estimate the numbers of the galaxies at redshift z > 0.5, which are expected to be detected in the PAH features or dust continuum by a wide (10 deg^2) or deep (1 deg^2) blind survey, both for a given observation time of 600 hours. As by-products of the wide blind survey, we also expect to detect debris disks, through the mid-IR excess above the photospheric emission of nearby main-sequence stars, and we estimate their number. We demonstrate that the SPICA mid-IR surveys will efficiently provide us with unprecedentedly large spectral samples, which can be studied further in the far-IR with SPICA.

In this paper we explore the intermediate line region (ILR) by using the photoionisation simulations of the gas clouds present at different radial distances from the center, corresponding to the locations from BLR out to NLR in four types of AGNs. We let for the presence of dust whenever conditions allow for dust existence. All spectral shapes are taken from the recent multi-wavelength campaigns. The cloud density decreases with distance as a power law. We found that the slope of the power law density profile does not affect the line emissivity radial profiles of major emission lines: H${\beta}$, He~II, Mg~II, C~III] ~and [O~III]. When the density of the cloud at the sublimation radius is as high as 10$^{11.5}$ cm$^{-3}$, the ILR should clearly be seen in the observations independently of the shape of the illuminating radiation. Moreover, our result is valid for low ionization nuclear emission regions of active galaxies.

The physical origin of the Broad Line Region in Active Galactic Nuclei is still unclear despite many years of observational studies. The reason is that the region is unresolved and the reverberation mapping results imply complex velocity field. We adopt a theory-motivated approach to identify the principal mechanism responsible for this complex phenomenon. We consider the possibility that the role of dust is essential. We assume that the local radiation pressure acting on the dust in the accretion disk atmosphere launches the outflow of material, but higher above the disk the irradiation from the central parts cause the dust evaporation and a subsequent fall back. This failed radiatively accelerated dusty outflow (FRADO) is expected to represent the material forming low ionization lines. In this paper we formulate simple analytical equations describing the cloud motion, including the evaporation phase. The model is fully described just by the basic parameters: black hole mass, accretion rate, black hole spin and the viewing angle. We study how the spectral line generic profiles correspond to this dynamics. We show that the virial factor calculated from our model strongly depends on the black hole mass in case of enhanced dust opacity, and thus it then correlates with the line width. This could explain why the virial factor measured in galaxies with pseudo-bulges differs from that obtained from objects with classical bulges although the trend predicted by the current version of the model is opposite to the observed trend.

The thermal stability of accretion disk and the possibility to see a limit-cycle behaviour strongly depends on the ability of the disk plasma to cool down. Various processes connected with radiation-matter interaction appearing in hot accretion disk plasma contribute to opacity. For the case of geometrically thin and optically thick accretion disk, we can estimate the influence of several different components of function \kappa, given by the Roseland mean. In the case of high temperatures, the electron Thomson scattering is dominant. At lower temperatures atomic processes become important. The slope d log \kappa/d log T can have locally stabilizing or destabilizing effect on the disk. Although the local MHD simulation postulate the stabilizing influence of the atomic processes, only the global time-dependent model can reveal the global disk stability range estimation. This is due to global diffusive nature of that processes. In this paper, using previously tested GLADIS code with modified prescription of the viscous dissipation, we examine the stabilizing effect of the Iron Opacity Bump.

We develope a self-consistent description of the Broad Line Region based on the concept of the failed wind powered by the radiation pressure acting on dusty accretion disk atmosphere in Keplerian motion. The material raised high above the disk is illuminated, dust evaportes, and the matter falls back towards the disk. This material is the source of emission lines. The model predicts the inner and outer radius of the region, the cloud dynamics under the dust radiation pressure and, subsequently, just the gravitational field of the central black hole, which results in assymetry between the rise and fall. Knowledge of the dynamics allows to predict the shapes of the emission lines as functions of the basic parameters of an active nucleus: black hole mass, accretion rate, black hole spin (or accretion efficiency) and the viewing angle with respect to the symmetry axis. Here we show preliminary results based on analytical approximations to the cloud motion.

Black holes (BHs) surrounded by accretion disks are present in the Universe at different scales of masses, from microquasars up to the active galactic nuclei (AGNs). The current picture remains ad hoc due to the complexity of the magnetic field action. In addition, accretion disks at high Eddington rates can be radiation-pressure dominated and, according to some of the heating prescriptions, thermally unstable. The observational verification of their resulting variability patterns may shed light on both the role of radiation pressure and magnetic fields. We compute the structure and time evolution of an accretion disk, using the code GLADIS (which models the global accretion disk instability). We supplement this model with a modified viscosity prescription, which can to some extent describe the magnetisation of the disk. We study the results for a large grid of models, to cover the whole parameter space, and we derive conclusions separately for different scales of black hole masses, which are characteristic for various types of cosmic sources. We show the dependencies between the flare or outburst duration, its amplitude, and period, on the accretion rate and viscosity scaling. We show that if the heating rate in the accretion disk grows more rapidly with the total pressure and temperature, the instability results in longer and sharper flares. In general, we confirm that the disks around the supermassive black holes are more radiation-pressure dominated and present relatively brighter bursts. Our method can also be used as an independent tool for the black hole mass determination, which we confront now for the intermediate black hole in the source HLX-1.

The Mg II emission line is visible in the optical band for intermediate redshift quasars (0.4 < z < 1.6) and it is thus an extremely important tool to measure the black hole mass and to understand the structure of the Broad Line Region. We aim to determine the substructure and the variability of the Mg II line with the aim to identify which part of the line comes from a medium in Keplerian motion. Using the Southern African Large Telescope (SALT) with the Robert Stobie Spectrograph (RSS) we performed ten spectroscopic observations of quasar HE 0435-4312 (z = 1.2231) over a period of three years (Dec 23/24, 2012 to Dec 7/8, 2015). Both the Mg II line and the Fe II pseudo-continuum increase with time. We clearly detect the systematic shift of the Mg II line with respect to the Fe II over the years, corresponding to the acceleration of 104 pm 14 km/s/year in the quasar rest frame. The Mg II line shape is clearly non-Gaussian but single-component, and the increase in line equivalent width and line shift is not accompanied with significant evolution of the line shape. We analyse the conditions in the Mg II and Fe II formation region and we note that the very large difference in the covering factor and the turbulent velocity also support the conclusion that the two regions are spatially separated. The measured acceleration of the line systematic shift is too large to connect it with the orbital motion at a distance of the Broad Line Region (BLR) in this source. It may imply a precessing inner disk illuminating the BLR. Further monitoring is still needed to better constrain the variability mechanism.

The hyperluminous X-ray source (HLX-1, the peak X-ray luminosity $\sim 10^{42}\rm erg\ s^{-1}$) near the spiral galaxy ESO 243-49 is possibly the best candidate for intermediate mass black hole (IMBH), which underwent recurrent outbursts with a period of $\sim 400$ days. The physical reason for this quasi-periodic variability is still unclear. We explore the possibility of radiation-pressure instability in accretion disk by modeling the light curve of HLX-1, and find that it can roughly reproduce the duration, period and amplitude of the recurrent outbursts HLX-1 with an IMBH of ~10^5Msun. Our result provides a possible mechanism to explain the recurrent outbursts in HLX-1. We further find a universal correlation between the outburst duration and the bolometric luminosity for the BH sources with a very broad mass range (e.g., X-ray binaries, XRBs, HLX-1 and active galactic nuclei, AGNs), which is roughly consistent with the prediction of radiation-pressure instability of the accretion disk. These results imply that "heartbeat" oscillations triggered by radiation-pressure instability may appears in different-scale BH systems.

The origin of the Broad Line Region (BLR) in active galaxies remains unknown. It seems to be related to the underlying accretion disk but an efficient mechanism is required to rise the material from the disk surface without giving too strong signatures of the outflow in the case of the low ionization lines. We discuss in detail two proposed mechanisms: (i) radiation pressure acting on dust in the disk atmosphere creating a failed wind (ii) the gravitational instability of the underlying disk. We compare the predicted location of the inner radius of the BLR in those two scenarios with the observed position obtained from the reverberation studies of several active galaxies. The failed dusty outflow model well represents the observational data while the predictions of the self-gravitational instability are not consistent with observations. The issue remains why actually we do not see any imprints of the underlying disk instability in the BLR properties.

We study the conditions for the onset of Thermal Instability in the innermost regions of compact galactic nuclei, where the properties of the interstellar environment are governed by the interplay of quasi-spherical accretion onto a supermassive black hole (SMBH) and the heating/cooling processes of gas in a dense nuclear star cluster. Stellar winds are the source of material for radiatively inefficient (quasi-spherical, non-magnetised) inflow/outflow onto the central SMBH, where a stagnation point develops within the Bondi type accretion. We study the local thermal equilibrium to determine the parameter space which allows cold and hot phases in mutual contact to co-exist. We include the effects of mechanical heating by stellar winds and radiative cooling/heating by the ambient field of the dense star cluster. We consider two examples: the Nuclear Star Cluster (NSC) in the Milky Way central region (including the gaseous Mini-spiral of Sgr~A*), and the Ultra-Compact Dwarf galaxy M60-UCD1. We find that the two systems behave in different ways because they are placed in different areas of parameter space in the instability diagram: gas temperature vs. dynamical ionization parameter. In the case of Sgr~A*, stellar heating prevents the spontaneous formation of cold clouds. The plasma from stellar winds joins the hot X-ray emitting phase and forms an outflow. In M60-UCD1 our model predicts spontaneous formation of cold clouds in the inner part of the galaxy. These cold clouds may survive since the cooling timescale is shorter than the inflow/outflow timescale.

We show that the recently observed suppression of the gap between the broad line region (BLR) and the narrow line region (NLR) in some AGN can be fully explained by an increase of the gas density in the emitting region. Our model predicts the formation of the intermediate line region (ILR) that is observed in some Seyfert galaxies by the detection of emission lines with intermediate velocity full width half maximum (FWHM) $\sim$ 700 - 1200 km s$^{-1}$. These lines are believed to be originating from an ILR located somewhere between the BLR and NLR. As it was previously proved, the apparent gap is assumed to be caused by the presence of dust beyond the sublimation radius. Our computations with the use of {\sc cloudy} photoionization code, show that the differences in the shape of spectral energy distribution (SED) from the central region of AGN, do not diminish the apparent gap in the line emission in those objects. A strong discontinuity in the line emission vs radius exists for all lines at the dust sublimation radius. However, increasing the gas density to $\sim$ 10$^{11.5}$ cm$^{-3}$ at the sublimation radius provides the continuous line emission vs radius and fully explains the recently observed lack of apparent gap in some AGN. We show that such a high density is consistent with the density of upper layers of an accretion disk atmosphere. Therefore, the upper layers of the disk atmosphere can give rise to the formation of observed emission line clouds.

The black hole mass measurement in active galaxies is a challenge, particularly in sources where the reverberation method cannot be applied. We aim to determine the black hole mass in a very special object, RE J1034+396, one of the two AGN with QPO oscillations detected in X-rays, and a single bright AGN with optical band totally dominated by starlight. We fit the stellar content using the code starlight, and the broad band disk contribution to optical/UV/X-ray emission is modeled with optxagnf. We also determine the black hole mass using several other independent methods. Various methods give contradictory results. Most measurements of the blacc hole mass are in the range 1.e6-1.e7 Msun, and the measurements based on dynamics give higher values than measurements based on Hbeta and Mg II emission lines.

The nature of dark energy, driving the accelerated expansion of the Universe, is one of the most important issues in modern astrophysics. In order to understand this phenomenon, we need precise astrophysical probes of the universal expansion spanning wide redshift ranges. Quasars have recently emerged as such a probe, thanks to their high intrinsic luminosities and, most importantly, our ability to measure their luminosity distances independently of redshifts. Here we report our ongoing work on observational reverberation mapping using the time delay of the Mg II line, performed with the South African Large Telescope (SALT).

The continuum-fitting method is one of the two most advanced methods of determining the black hole spin in accreting X-ray binary systems. There are, however, still some unresolved issues with the underlying disk models. One of them manifests as an apparent decrease in spin for increasing source luminosity. Here, we perform a few simple tests to establish whether outflows from the disk close to the inner radius can address this problem. We employ four different parametric models to describe the wind and compare these to the apparent decrease in spin with luminosity measured in the sources LMC~X-3 and GRS~1915+105. Wind models in which parameters do not explicitly depend on the accretion rate cannot reproduce the spin measurements. Models with mass accretion rate dependent outflows, however, have spectra that emulate the observed ones. The assumption of a wind thus effectively removes the artifact of spin decrease. This solution is not unique; the same conclusion can be obtained with a truncated inner disk model. To distinguish among valid models, high resolution X-ray data and a realistic description of the Comptonization in the wind will be needed.

Warm absorber (WA) is an ionised gas present in the line of sight to the AGN central engine. The effect of the absorber is imprinted in the absorption lines observed in X-ray spectra of AGN. In this work, we model the WA in Seyfert 1 galaxy Mrk 509 using its recently published shape of broad band spectral energy distribution (SED) as a continuum illuminating the absorber. Using the photoionization code {\sc Titan}, recently we have shown that the absorption measure distribution (AMD) found for this object can be successfully modelled as a single slab of gas in total pressure (radiation+gas) equilibrium, contrary to the usual models of constant density multiple slabs. We discuss the transmitted spectrum that would be recorded by an observer after the radiation from the nucleus passes through the WA.

In this paper we model the observed absorption measure distribution (AMD) in Mrk 509, which spans three orders of magnitude in ionization level with a single-zone absorber in pressure equilibrium. AMD is usually constructed from observations of narrow absorption lines in radio-quiet active galaxies with warm absorbers. We study the properties of the warm absorber in Mrk~509 using recently published broad-band spectral energy distribution observed with different instruments. This spectrum is an input in radiative transfer computations with full photoionisation treatment using {\sc titan} code. We show that the simplest way to fully reproduce the shape of AMD is to assume that the warm absorber is a single zone under constant total pressure. With this assumption we found theoretical AMD which matches the observed one determined on the basis of 600 ks RGS XMM-Newton spectrum of Mrk~509. The softness of the source spectrum and the important role of the free-free emission breaks the usual degeneracy in the ionization state calculations, and the explicit dependence of the depths of AMD dips on density open a new path to density diagnostic for the warm absorber. In Mrk 509 the implied density is of the order of $10^8$ cm$^{-3}$.

We review the current state of understanding how accretion onto a black hole proceeds and what are they key elements needed to form relativistic jets. Since the theoretical progress is severely halted by the lack of deep understanding of the microphysics involved, particularly in the presence of strong magnetic fields, all proposed solutions have a status of models, and their selection has to be based directly on observational constraints.

In the past years, several observations of AGN and X-ray binaries have suggested the existence of a warm T around 0.5-1 keV and optically thick, \tau ~ 10-20, corona covering the inner parts of the accretion disk. These properties are directly derived from spectral fitting in UV to soft-X-rays using Comptonization models. However, whether such a medium can be both in radiative and hydrostatic equilibrium with an accretion disk is still uncertain. We investigate the properties of such warm, optically thick coronae and put constraints on their existence. We solve the radiative transfer equation for grey atmosphere analytically in a pure scattering medium, including local dissipation as an additional heating term in the warm corona. The temperature profile of the warm corona is calculated assuming it is cooled by Compton scattering, with the underlying dissipative disk providing photons to the corona. Our analytic calculations show that a dissipative thick, (\tau_{cor} ~ 10-12) corona on the top of a standard accretion disk can reach temperatures of the order of 0.5-1 keV in its upper layers provided that the disk is passive. But, in absence of strong magnetic fields, the requirement of a Compton cooled corona in hydrostatic equilibrium in the vertical direction sets an upper limit on the Thomson optical depth \tau_{cor} < 5 . We show this value cannot be exceeded independently of the accretion disk parameters. However, magnetic pressure can extend this result to larger optical depths. Namely, a dissipative corona might have an optical depth up to ~ 20 when the magnetic pressure is 100 times higher that the gas pressure. The observation of warm coronae with Thomson depth larger than ~ 5 puts tights constraints on the physics of the accretion disk/corona systems and requires either strong magnetic fields or vertical outflows to stabilize the system.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of active galactic nuclei. For a summary, we refer to the paper.

We propose a very simple physical mechanism responsible for the formation of the Low Ionization Line part of the Broad Line Region in Active Galactic Nuclei. It explains the scaling of the Broad Line Region size with the monochromatic luminosity, including the exact slope and the proportionality constant, seen in the reverberation studies of nearby sources. The scaling is independent from the mass and accretion rate of an active nucleus. The mechanism predicts the formation of a dust-driven wind in the disk region where the local effective temperature of a non-illuminated accretion disk drops below 1000 K and allows for dust formation. We explore now the predictive power of the model with the aim to differentiate between this model and the previously proposed mechanisms of the formation of the Broad Line Region. We discuss the expected departures from the universal scaling at long wavelength, and the role of the inclination angle of the accretion disk in the source. We compare the expected line profiles with Mg II line profiles in the quasars observed by us with the SALT telescope. We also discuss the tests based on the presence or absence of the broad emission lines in low luminosity active galaxies. Finally, we discuss the future tests of the model to be done with expected ground-based observations and satellite missions.

We introduce a novel formalism to investigate the role of the spin angular momentum of astrophysical black holes in influencing the behaviour of low angular momentum general relativistic accretion. We propose a metric independent analysis of axisymmetric general relativistic flow, and consequently formulate the space and time dependent equations describing the general relativistic hydrodynamic accretion flow in the Kerr metric. The associated stationary critical solutions for such flow equations are provided and the stability of the stationary transonic configuration is examined using an elegant linear perturbation technique. We examine the properties of infalling material for both prograde and retrograde accretion as a function of the Kerr parameter at extremely close proximity to the event horizon. Our formalism can be used to identify a new spectral signature of black hole spin, and has the potential of performing the black hole shadow imaging corresponding to the low angular momentum accretion flow.

Quasars, as the most luminous persistent sources in the Universe, have broad applications for cosmological studies. In particular, they can be employed to directly measure the expansion history of the Universe, similarly to SNe Ia. The advantage of quasars is that they are numerous, cover a broad range of redshifts, up to $z = 7$, and do not show significant evolution of metallicity with redshift. The idea is based on the relation between the time delay of an emission line and the continuum, and the absolute monochromatic luminosity of a quasar. For intermediate redshift quasars, the suitable line is Mg II. Between December 2012 and March 2014, we performed five spectroscopic observations of the QSO CTS C30.10 ($z = 0.900$) using the South African Large Telesope (SALT), supplemented with photometric monitoring, with the aim of determining the variability of the line shape, changes in the total line intensity and in the continuum. We show that the method is very promising.

We explore the conditions for the thermal instability to operate in the mini-spiral region of the Galactic centre (Sgr A*), where both the hot and cold media are known to coexist. The photoionisation Cloudy calculations are performed for different physical states of plasma. We neglect the dynamics of the material and concentrate on the study of the parameter ranges where the thermal instability may operate, taking into account the past history of Sgr A* bolometric luminosity. We show that the thermal instability does not operate at the present very low level of the Sgr A* activity. However, Sgr A* was much more luminous in the past. For the highest luminosity states the two-phase medium can be created up to 1.4 pc from the centre. The presence of dust grains tends to suppress the instability, but the dust is destroyed in the presence of strong radiation field and hot plasma. The clumpiness is thus induced in the high activity period, and the cooling/heating timescales are long enough to preserve later the past multi-phase structure. The instability enhances the clumpiness of the mini-spiral medium and creates a possibility of episodes of enhanced accretion of cold clumps towards Sgr A*. The mechanism determines the range of masses and sizes of clouds; under the conditions of Sgr A*, the likely values come out $1$ - $10^2M_{\oplus}$ for the cloud typical mass.

Quasars can be used as a complementary tool to SN Ia to probe the distribution of dark energy in the Universe by measuring the time delay of the emission line with respect to the continuum. The understanding of the Mg II emission line structure is important for cosmological application and for the black hole mass measurements of intermediate redshift quasars. The knowledge of the shape of Mg II line and its variability allows to determine which part of the line should be used to measure the time delay and the black hole mass. We thus aim at determination of the structure and the variability of the Mg II line as well as of the underlying Fe II pseudo-continuum. We performed five spectroscopic observations of a quasar CTS C30.10 (z = 0.9000) with the SALT telescope between December 2012 and March 2014, and we studied the variations of the spectral shape in the 2700 - 2900 A rest frame. We showed that the Mg II line in this source consists of two kinematic components, which makes the source representative of type B quasars. Both components were well modeled with a Lorentzian shape. Both components vary in a similar way. The Fe II contribution seems to be related only to the first (blue) Mg II component}. Broad band spectral fitting rather favor the use of the whole line profile. The contribution of the Narrow Line Region to Mg II is very low, below 2%. The Mg II variability is lower than the variability of the continuum, which is consistent with the simple reprocessing scenario. The variability level of CTS C30.10 and the measurement accuracy of the line and continuum is high enough to expect that further monitoring will allow to measure the time delay between the Mg II line and continuum.

The Mg II line is of extreme importance in intermediate redshift quasars since it allows us to measure the black hole mass in these sources and to use these sources as probes of the distribution of dark energy in the Universe, as a complementary tool to SN Ia. Reliable use of Mg II requires a good understanding of all the systematic effects involved in the measurement of the line properties, including the contamination by Fe II UV emission. We performed three spectroscopic observations of a quasar LBQS 2113-4538 (z = 0.956) with the SALT telescope separated in time by several months and we analyze in detail the mean spectrum and the variability in the spectral shape. We show that even in our good-quality spectra the Mg II doublet is well fit by a single Lorentzian shape. We tested several models of the Fe II pseudo-continuum and showed that one of them well represents all the data. The amplitudes of both components vary in time, but the shapes do not change significantly. The measured line width of LBQS 2113-4538 identifies this object as a class A quasar. The upper limit of $3\%$ for the contribution of the Narrow Line Region (NLR) to Mg II may suggest that the separation of the Broad Line Region (BLR) and NLR disappears in this class of objects.

We study some general properties of accretion disc variability in the context of stationary random processes. In particular, we are interested in mathematical constraints that can be imposed on the functional form of the Fourier power-spectrum density (PSD) that exhibits a multiply broken shape and several local maxima. We develop a methodology for determining the regions of the model parameter space that can in principle reproduce a PSD shape with a given number and position of local peaks and breaks of the PSD slope. Given the vast space of possible parameters, it is an important requirement that the method is fast in estimating the PSD shape for a given parameter set of the model. We generated and discuss the theoretical PSD profiles of a shot-noise-type random process with exponentially decaying flares. Then we determined conditions under which one, two, or more breaks or local maxima occur in the PSD. We calculated positions of these features and determined the changing slope of the model PSD. Furthermore, we considered the influence of the modulation by the orbital motion for a variability pattern assumed to result from an orbiting-spot model. We suggest that our general methodology can be useful in for describing non-monotonic PSD profiles (such as the trend seen, on different scales, in exemplary cases of the high-mass X-ray binary Cygnus X-1 and the narrow-line Seyfert galaxy Ark 564). We adopt a model where these power spectra are reproduced as a superposition of several Lorentzians with varying amplitudes in the X-ray-band light curve. Our general approach can help in constraining the model parameters and in determining which parts of the parameter space are accessible under various circumstances.

Gas clouds are present in the Galactic centre, where they orbit around the supermassive black hole. Collisions between these clumps reduce their angular momentum, and as a result some of the clumps are set on a plunging trajectory. Constraints can be imposed on the nature of past accretion events based on the currently observed X-ray reflection from the molecular clouds. We discuss accretion of clouds in the context of enhanced activity of Sgr A* during the past few hundred years. We put forward a scenario according to which gas clouds bring material close to the horizon of the black hole on <~0.1 parsec scale. We have modelled the source intrinsic luminosity assuming that multiple events occur at various moments in time. These events are characterized by the amount of accreted material and the distribution of angular momentum. We parameterized the activity in the form of a sequence of discrete events, followed the viscous evolution, and calculated the luminosity of the system from the time-dependent accretion rate across the inner boundary. Accreting clumps settle near a circularization radius, spread there during the viscous time, and subsequently feed the black hole over a certain period. A significant enhancement (by factor of ten) of the luminosity is only expected if the viscous timescale of the inflow is very short. On the other hand, the increase in source activity is expected to be much less prominent if the latter timescale is longer and a considerable fraction of the material does not reach the centre. A solution is obtained under two additional assumptions: (i) the radiative efficiency is a decreasing function of the Eddington ratio; (ii) the viscous decay of the luminosity proceeds somewhat faster than the canonical L(t)~t^{-5/3} profile. We applied our scheme to the case of G2 cloud in the Galactic centre to obtain constraints on the core-less gaseous cloud model.

High redshift quasars can be used to deduce the distribution of dark energy in the Universe, as a complementary tool to SN Ia. The method is based on determination of the size of the Broad Line Region from the emission line delay, determination of the absolute monochromatic luminosity either from the observed statistical relation or from a model of the formation of the Broad Line Region, and determination of the observed monochromatic flux from photometry. This allows to obtain the luminosity distance to a quasar independently from its redshift. The accuracy of the measurements is however, a key issue. We model the expected accuracy of the measurements by creating artificial quasar monochromatic lightcurves and responses from the Broad Line Region under various assumptions about the variability of a quasar, Broad Line Region extension, distribution of the measurements in time, accuracy of the measurements and the intrinsic line variability. We show that the five year monitoring based on Mg II line should give the accuracy of 0.06 - 0.32 magnitude in the distance modulus which allows to put interesting constraints on the cosmological models. Monitoring of higher redshift quasars based on CIV lines is problematic due to much higher level of the intrinsic variability of CIV in comparison with Mg II.

We present photoionisation modelling of the intrinsic absorber in the bright quasar HS 1603+3820. We construct spectral energy distribution using optical/UV/X-ray observations as an input to photoionisation calculations. Spectra from Keck telescope show extremely high ratio of CIV to HI, for the first absorber in system A, named A1. This value, together with high column density of CIV ion, puts strong constraints on photoionisation model. We use photoionisation modeling to derive hydrogen number density at the cloud illuminated surface. Estimating bolometric luminosity of HS 1603+3820, from typical formula for quasars, we calculate the distance to A1. The derived location is as close as 0.1 pc, and situates an absorber even closer to the nucleus than the possible location of the Broad Line Region in this object

We discuss the issues of stability of accretion disks that may undergo the limit-cycle oscillations due to the two main types of thermal-viscous instabilities. These are induced either by the domination of radiation pressure in the innermost regions close to the central black hole, or by the partial ionization of hydrogen in the zone of appropriate temperatures. These physical processes may lead to the intermittent activity in AGN on timescales between hundreds and millions of years. We list a number of observational facts that support the idea of the cyclic activity in high accretion rate sources. We conclude however that the observed features of quasars may provide only indirect signatures of the underlying instabilities. Also, the support from the sources with stellar mass black holes, whose variability timescales are observationally feasible, is limited to a few cases of the microquasars. Therefore we consider a number of plausible mechanisms of stabilization of the limit cycle oscillations in high accretion rate accretion disks. The newly found is the stabilizing effect of the stochastic viscosity fluctuations.

The most characteristic property of active galaxies, including quasars, are prominent broad emission lines. I will discuss an interesting possibility that dust is responsible for this phenomenon. The dust is known to be present in quasars in the form of a dusty/molecular torus which results in complexity of the appearance of active galaxies. However, this dust is located further from the black hole than the Broad Line Region. We propose that the dust is present also closer in and it is actually responsible for formation of the broad emission lines. The argument is based on determination of the temperature of the disk atmosphere underlying the Broad Line Region: it is close to 1000 K, independently from the black hole mass and accretion rate of the object. The mechanism is simple and universal but leads to a considerable complexity of the active nucleus surrounding. The understanding the formation of BLR opens a way to use it reliably - in combination with reverberation measurement of its size - as standard candles in cosmology.

The question of the origin of the gas supplying the accretion process is pertinent especially in the context of enhanced activity of Galactic Center during the past few hundred years, seen now as echo from the surrounding molecular clouds, and the currently observed new cloud approaching Sgr A*. We discuss the so-called Galactic Center mini-spiral as a possible source of material feeding the supermassive black hole on a 0.1 parsec scale. The collisions between individual clumps reduce their angular momentum. and set some of the clumps on a plunging trajectory. We conclude that the amount of material contained in the mini-spiral is sufficient to sustain the luminosity of Sgr A* at the required level. The accretion episodes of relatively dense gas from the mini-spiral passing through a transient ring mode at $\sim 10^4$ Rg provide a viable scenario for the bright phase of Galactic Center.

Gamma-ray catalogs contain a considerable amount of unidentified sources. Many of these are located out of the Galactic plane and therefore may have extragalactic origin. Here we assume that the formation of massive black holes in galactic nuclei proceeds through a quasi-star stage and consider the possibility of jet production by such objects. Those jets would be the sources of collimated synchrotron and Compton emission, extending from radio to gamma rays. The expected lifetimes of quasi-stars are of the order of million of years while the jet luminosities, somewhat smaller than that of quasar jets, are sufficient to account for the unidentified gamma-ray sources. The jet emission dominates over the thermal emission of a quasi-star in all energy bands, except when the jet is not directed towards an observer. The predicted synchrotron emission peaks in the IR band, with the flux close to the limits of the available IR all sky surveys. The ratio of the $\gamma$-ray flux to the IR flux is found to be very large ($\sim 60$), much larger than in BL Lac objects but reached by some radio-loud quasars. On the other hand, radio-loud quasars show broad emission lines while no such lines are expected from quasi-stars. Therefore the differentiation between various scenarios accounting for the unidentified gamma-ray sources will be possible at the basis of the photometry and spectroscopy of the IR/optical counterparts.

Strong broad emission lines are the most important signatures of active galactic nuclei. These lines allowed to discover the cosmological nature of quasars, and at present these lines allow for convenient method of weighting the black holes residing in their nuclei. However, a question remains why such strong lines form there in the first place. Specifically, in the case of Low Ionization Lines, there must be a mechanism which leads to an efficient rise of the material from the surface of the accretion disk surrounding a black hole but at the same time should not give a strong signature of the systematic outflow, as the Balmer lines are not significantly shifted with respect to the Narrow Line Region. We determine the effective temperature of the accretion disk underlying the H$\beta$ line at the basis of the time delay measured from reverberation and the simple Shakura-Sunyaev theory of accretion disks. We obtain that this temperature is universal, and equal $995 \pm 74$ K, independently from the black hole mass and accretion rate of the source. This result suggests to us that the dust formation in the disk atmosphere is responsible for the strong rise of the material. However, as the material gains height above the disk it becomes irradiated, the dust evaporates, the radiation pressure force suddently drops and the material fall back again at the disk. Therefore, a failed wind forms. In the simple version of the model the disk irradiation is neglected, but in the present paper we also discuss this irradiation and we use the observed variation of the Broad Line Region in NGC 5548 to constrain the character of this non-local non-stationary phenomenon. The current instruments cannot resolve the Broad Line Region but future instrumentation may allow to test the model directly.

We analyze the time variability of the X-ray emission of RE J1034+396, an active galactic nucleus with the first firm detection of a quasi-periodic oscillations (QPO). Based on the results of a wavelet analysis, we find a drift in the QPO central frequency. The change inthe QPO frequency correlates with the change in the X-ray flux with a short time delay. Linear structures such as shocks, spiral waves, orvery distant flares seem to be a favored explanation for this particular QPO event.

Blandford-Znajek (BZ) mechanism has usually been studied in the literature for accretion with considerably high angular momentum leading either to the formation of a cold Keplerian disc, or a hot and geometrically thick sub-Keplerian flow as described within the framework of ADAF/RIAF. However, in nearby elliptical galaxies, as well as for our own Galactic centre, accretion with very low angular momentum is prevalent. Such quasi-spherical strongly sub-Keplerian accretion has complex dynamical features and can accommodate stationary shocks. In this letter, we present our calculation for the maximum efficiency obtainable through the BZ mechanism for complete general relativistic weakly rotating axisymmetric flow in the Kerr metric. Both shocked and shock free flow has been studied in detail for rotating and counter rotating accretion. Such study has never been done in the literature before. We find that the energy extraction efficiency is low, about 0.1%, and increases by a factor 15 if the ram pressure is included. Such an efficiency is still much higher than the radiative efficiency of such optically thin flows. For BZ mechanism, shocked flow produces higher efficiency than the shock free solutions and retrograde flow provides a slightly larger value of the efficiency than that for the prograde flow.

We present the state of current research of Very Short Gamma Ray Bursts (VSGRBs) from seven GRB detectors. We found that VSGRBs form distinct class of GRBs, which in our opinion, in most cases can originate from the evaporating Primordial Black Holes (PBHs). Arguments supporting our opinion: 1. GRBs with time duration (T90) < 100 ms form distinct class: VSGRBs. 2. We observe significant anisotropy in the galactic angular distribution of BATSE VSGRB events. 3. V/Vmax distribution for BATSE VSGRB events indicates the local distance production. 4. VSGBBs have more energetic {\gamma}-ray burst than other GRBs with longer duration (KONUS). 5. We observe small number of afterglows in SWIFT VSGRB sample (25%), in contrast with the noticeable afterglow frequency in SGRB sample (78%). 6. Time profile of rising part BATSE VSGRBs is in agreement with the evaporation PBH model.

We study the variability mechanism of active galactic nuclei (AGN) within the framework of the flare model. To this end we examine the case of Seyfert/LINER galaxy NGC 4258, which is observed at high inclination angle and exhibits rapid fluctuations of the X-ray light curve. We construct a model light curve based on the assumption of magnetic flares localized in the equatorial plane and orbiting with Keplerian speed at each given radius. We calculate the level of variability as a function of the inclination of an observer, taking into account all effects of general relativity near a rotating supermassive black hole. The variability level is a monotonic function of the source inclination. It rises more rapidly for larger values of the black hole spin (Kerr parameter) and for steeper emissivity (index beta of the radial profile). We compare the expected level of variability for the viewing angle 81.6 deg, as inferred for NGC 4258, with the case of moderate viewing angles about 30 deg, typical for Seyfert type-1 galaxies. Highly inclined sources such as this one are particularly suitable to test the flare model because the effects of orbital motion, Doppler boosting and light bending are all expected to have maximum when the accretion disk is seen almost edge-on. The model is consistent with the NGC 4258 variability, where the obscuring material is thought to be localized mainly towards the equatorial plane rather than forming a geometrically thick torus. Once the intrinsic time-scales of the flare duration are determined to better precision, this kind of highly inclined objects with a precisely known mass of the black hole can be used to set independent constraints on the spin parameter.

The spin of the black hole is an important parameter which may be responsible for the properties of the inflow and outflow of the material surrounding a black hole. Broad band IR/optical/UV spectrum of the quasar SDSS J094533.99+100950.1 is clearly disk-dominated, with the spectrum peaking up in the observed frequency range. Therefore, disk fitting method usually used for Galactic black holes can be used in this object to determine the black hole spin. We develop the numerical code for computing disk properties, including radius-dependent hardening factor, and we apply the ray-tracing method to incorporate all general relativity effects in light propagation. We show that the simple multicolor disk model gives a good fit, without any other component required, and the disk extends down to the marginally stable orbit. The best fit accretion rate is 0.13, well below the Eddington limit, and the black hole spin is moderate, 0.3. The contour error for the fit combined with the constraints for the black hole mass and the disk inclination gives a constraint that the spin is lower than 0.8. We discuss the sources of possible systematic errors in the parameter determinations.

We discuss two important instability mechanisms that may lead to the limit-cycle oscillations of the luminosity of the accretion disks around compact objects: ionization instability and radiation-pressure instability. Ionization instability is well established as a mechanism of X-ray novae eruptions in black hole binary systems but its applicability to AGN is still problematic. Radiation pressure theory has still very weak observational background in any of these sources. In the present paper we attempt to confront the parameter space of these instabilities with the observational data. At the basis of this simple survey of sources properties we argue that the radiation pressure instability is likely to be present in several Galactic sources with the Eddington ratios above 0.15, and in AGN with the Eddington ratio above 0.025. Our results favor the parameterization of the viscosity through the geometrical mean of the radiation and gas pressure both in Galactic sources and AGN. More examples of the quasi-regular outbursts in the timescales of 100 seconds in Galactic sources, and hundreds of years in AGN are needed to formulate firm conclusions. We also show that the disk sizes in the X-ray novae are consistent with the ionization instability. This instability may also considerably influence the lifetime cycle and overall complexity in the supermassive black hole environment.

Observation of the bright Seyfert 1 galaxy RE J1034+396 is believed to demonstrate a drift of the central period of the Quasi Periodic Oscillation (QPO) linearly correlated with the temporary X-ray luminosity. We show, using a specific scenario of the oscillation mechanism in black hole accretion disc, that modeling such correlated trends puts very strong constraints on the nature of this oscillation and the characteristic features of the hot flow in Active Galactic Nuclei (AGN). In our model, QPO oscillations are due to the oscillations of the shock formed in the low angular momentum hot accretion flow, and the variation of the shock location corresponds to the observed changes in the QPO period and the X-ray flux. In this scenario, change in the shock location caused by perturbation of the flow angular momentum is compatible with the trends observed in RE J1034+396, whereas the perturbation of the specific flow energy results in too strong flux response to the change of the oscillation period. Using a complete general relativistic framework to study the accretion flow in the Kerr metric, we discuss the role of the black hole spin in the period drift. Future missions are expected to bring more active galaxies with time-resolved quasi-periodic oscillations so similar quantitative study for other QPO scenarios will be necessary.

Although broad emission lines are the most reliable signature of the nuclear activity of a galaxy and the location of the emitting material is well measured by the reverberation method, the physical cause of the formation of the broad line region remains unclear. We attempt to place some constraints on its origin. We study the properties of the accretion disk underlying the broad line region. We find that the effective temperature at the disk radius corresponding to the location of the broad line region, as inferred from the H$\beta$ line, is universal in all monitored sources and equal to 1000 K. This value is close to the limiting value that permits for the existence of the dust. The likely origin of the low ionization part of the broad line region is the strong local dusty wind from the disk. This wind becomes exposed to the irradiation by the central regions when moving higher above the disk surface and subsequently behaves like a failed wind, thus leading to a local mixture of inflow and outflow. This may provide the physical explanation of the turbulence needed both to smooth the line profiles as well as provide additional mechanical heating.

We present various instability mechanisms in the accreting black hole systems which might indicate at the connection between the accretion disk and jet. The jets observed in microquasars can have a persistent or blobby morphology. Correlated with the accretion luminosity, this might provide a link to the cyclic outbursts of the disk. Such duty-cycle type of behavior on short timescales results from the thermal instability caused by the radiation pressure domination. The same type of instability may explain the cyclic radioactivity of the supermassive black hole systems. The somewhat longer timescales are characteristic for the instability caused by the partial hydrogen ionization. The distortions of the jet direction and complex morphology of the sources can be caused by precession of the disk-jet axis.

We analyze the time variability of the X-ray emission of RE J1034+396 -- an active galactic nucleus with the first firm detection of a quasi-periodic oscillations (QPO). Based on the results of a wavelet analysis, we find a drift in the QPO central frequency. The change in the QPO frequency correlates with the change in the X-ray flux with a short time delay. The data specifically suggest a linear dependence between the QPO period and the flux, and this gives important constraints on the QPO models. In particular, it excludes explanation in terms of the orbiting hot spot model close to a black hole. Linear structures such as shocks, spiral waves, or very distant flares are favored.

We discuss a very peculiar subgroup of gamma-ray bursts among the BATSE sources. These bursts are very short ($T_{90} \le $0.1 s), hard, and came predominantly from a restricted direction of the sky (close to the Galactic anti-center). We analyze their arrival times and possible correlations, as well as the profiles of individual bursts. We find no peculiarities in the arrival times of Very Short Bursts (VSBs) despite their highly non-uniform spatial distribution. There is no dependence in the burst shapes on location. Bursts coming both from the burst-enhancement Galactic Anticenter region and from all other directions show considerable dispersion in their rise and fall times. Significant fraction of VSBs have multiple peaks despite their extremely short duration. Burst time properties are most likely to be consistent with two origin mechanisms: either with binary NS-NS mergers with low total masses passing through a phase of hypermassive neutron star, or with evaporation of the primordial black holes in the scenario of no photosphere formation.

Large scale X-ray jets that extend to >100 kpc distances from the host galaxy indicate the importance of jets interactions with the environment on many different physical scales. Morphology of X-ray clusters indicate that the radio-jet activity of a cD galaxy is intermittent. This intermittency might be a result of a feedback and/or interactions between galaxies within the cluster. Here we consider the radiation pressure instability operating on short timescales (<10^5 years) as the origin of the intermittent behaviour. We test whether this instability can be responsible for short ages (< 10^4 years) of Compact Symmetric Objects measured by hot spots propagation velocities in VLBI observations. We model the accretion disk evolution and constrain model parameters that may explain the observed compact radio structures and over-abundance of GPS sources. We also describe effects of consequent outbursts.

Weak emission line quasars are a rare and puzzling group of objects. In this paper we present one more object of this class found in the Sloan Digital Sky Survey (SDSS). The quasar SDSS J094533.99+100950.1, lying at z = 1.66, has practically no C IV emission line, a red continuum very similar to the second steepest of the quasar composite spectra of Richards et al., is not strongly affected by absorption and the Mg II line, although relatively weak, is strong enough to measure the black hole mass. The Eddington ratio in this object is about 0.45, and the line properties are not consistent with the trends expected at high accretion rates. We propose that the most probable explanation of the line properties in this object, and perhaps in all weak emission line quasars, is that the quasar activity has just started. A disk wind is freshly launched so the low ionization lines which form close to the disk surface are already observed but the wind has not yet reached the regions where high ionization lines or narrow line components are formed. The relatively high occurrence of such a phenomenon may additionally indicate that the quasar active phase consists of several sub-phases, each starting with a fresh build-up of the Broad Line Region.

In this small review we present the actual state the knowledge about weighting black holes. Black holes can be found in stellar binary systems in our Galaxy and in other nearby galaxies, in globular clusters, which we can see in our and nearby galaxies, and in centres of all well-developed galaxies. Range of values of their masses is wide and cover about ten orders of magnitude (not taking into account the hypothetic primordial black holes). Establishing the presence of black holes, and in particular the measurement of their mass is one on the key issues for many branches of astronomy, from stellar evolution to cosmology.

The secular evolution of the purely general relativistic low angular momentum accretion flow around a spinning black hole is shown to exhibit hysteresis effects. This confirms that a stationary shock is an integral part of such an accretion disc in the Kerr metric. The equations describing the space gradient of the dynamical flow velocity of the accreting matter have been shown to be equivalent to a first order autonomous dynamical systems. Fixed point analysis ensures that such flow must be multi-transonic for certain astrophysically relevant initial boundary conditions. Contrary to the existing consensus in the literature, the critical points and the sonic points are proved not to be isomorphic in general. Homoclinic orbits for the flow flow possessing multiple critical points select the critical point with the higher entropy accretion rate, confirming that the entropy accretion rate is the degeneracy removing agent in the system. However, heteroclinic orbits are also observed for some special situation, where both the saddle type critical points of the flow configuration possesses identical entropy accretion rate. Topologies with heteroclinic orbits are thus the only allowed non removable degenerate solutions for accretion flow with multiple critical points, and are shown to be structurally unstable. Depending on suitable initial boundary conditions, a homoclinic trajectory can be combined with a standard non homoclinic orbit through an energy preserving Rankine-Hugoniot type of stationary shock. An effective Lyapunov index has been proposed to analytically confirm why certain class of transonic flow can not accommodate shock solutions even if it produces multiple critical points. (Abridged)

Black hole mass determination in active galaxies is a key issue in understanding various luminosity states. In the present paper we try to generalise the mass determination method based on the X-ray excess variance, successfully used for typical broad line Seyfert 1 galaxies (BLS1) to Narrow Line Seyfert 1 (NLS1) galaxies. NLS1 galaxies differ from BLS1 with respect to several properties. They are generally more variable in 2-10 keV energy band so the natural expectation is the need to use a different scaling coefficient between the mass and the variance in these two types of sources. However, we find that such a simple approach is not enough. Although for majority of the 21 NLS1 galaxies in our sample a single scaling coefficient (larger by a factor 20) provided us with a satisfactory method of mass determination, in a small subset of NLS1 galaxies this approach failed. Variability of those objects appeared to be at the intermediate level between NLS1 and BLS1 galaxies. These exceptional NLS1 galaxies have much harder soft X-ray spectra than majority of NLS1 galaxies. We thus postulate that the division of Seyfert 1 galaxies into BLS1 and NLS1 according to the widths of the Hbeta line is less generic than according to the soft X-ray slope.

We associate the existence of short-lived compact radio sources with the intermittent activity of the central engine caused by a radiation pressure instability within an accretion disk. Such objects may constitute a numerous sub-class of Giga-Hertz Peaked Spectrum sources, in accordance with the population studies of radio-loud active galaxies, as well as detailed investigations of their radio morphologies. We perform the model computations assuming the viscosity parametrization as proportional to a geometrical mean of the total and gas pressure. The implied timescales are consistent with the observed ages of the sources. The duration of an active phase for a moderate accretion rate is short enough (< 10^3-10^4 years) that the ejecta are confined within the host galaxy and thus these sources cannot evolve into large size radio galaxies unless they are close to the Eddington limit.

A photoionized gas under constant pressure can display a thermal instability, with three or more solutions for possible thermal equilibrium. A unique solution of the structure of the irradiated medium is obtained only if electron conduction is considered. The subject of our study is to estimate how the effect of thermal conduction affects the structure and transmitted spectrum of the warm absorber computed by solving radiative transfer with the code TITAN. We developed a new computational mode for the code TITAN to obtain several solutions for a given external conditions and we test a posteriori which solution is the closest one to the required integral condition based on conduction. We demonstrate that the automatic mode of the code TITAN provides the solution to the radiative transfer which is generally consistent with the estimated exact solution within a few per cent accuracy, with larger errors for some line intensities (up to 20 per cent) for iron lines at intermediate ionization state.

We analyze the spectral energy distributions (SEDs) and emission line properties of the red (J-K$_S$ > 2) 2MASS AGN observed by Chandra using principle component analysis. The sample includes 44 low redshift AGN with low or moderate obscuration (N_H < 10^{23} cm^{-2}) as indicated by X-rays and SED modeling. The obscuration of the AGN allows us to see weaker emission components (host galaxy emission, AGN scattered light) which are usually outshone by the AGN. The first four eigenvectors explain 70% of the variance in the sample. The dominant cause of variance in the sample (eigenvector 1) is the L/Ledd ratio strengthened by intrinsic absorption. Eigenvector 2 is related to host galaxy (relative to the observed AGN) emission and eigenvectors 3 and 4 distinguish between two sources of obscuration: host galaxy absorption and circumnuclear absorption. Although our analysis is consistent with unification schemes where inclination dependent obscuration is important in determining the AGN SEDs, the L/Ledd ratio is the most important factor, followed by host galaxy emission.

Sgr A* is the closest massive black hole and can be observed with the highest angular resolution. Nevertheless, our current understanding of the accretion process in this source is very poor. The inflow is almost certainly of low radiative efficiency and it is accompanied by a strong outflow and the flow is strongly variable but the details of the dynamics are unknown. Even the amount of angular momentum in the flow is an open question. Here we argue that low angular momentum scenario is better suited to explain the flow variability. We present a new hybrid model which describes such a flow and consists of an outer spherically symmetric Bondi flow and an inner axially symmetric flow described through MHD simulations. The assumed angular momentum of the matter is low, i.e. the corresponding circularization radius in the equatorial plane of the flow is just above the innermost stable circular orbit in pseudo-Newtonian potential. We compare the radiation spectrum from such a flow to the broad band observational data for Sgr A*.

Theory of random processes provides an attractive mathematical tool to describe the fluctuating signal from accreting sources, such as active galactic nuclei and Galactic black holes observed in X-rays. These objects exhibit featureless variability on different timescales, probably originating from an accretion disc. We study the basic features of the power spectra in terms of a general framework, which permits semi-analytical determination of the power spectral density (PSD) of the resulting light curve. We consider the expected signal generated by an ensemble of spots randomly created on the accretion disc surface. Spot generation is governed by Poisson or by Hawkes processes. We include general relativity effects shaping the signal on its propagation to a distant observer. We analyse the PSD of a spotted disc light curve and show the accuracy of our semi-analytical approach by comparing the obtained PSD with the results of Monte Carlo simulations. The asymptotic slopes of PSD are 0 at low frequencies and they drop to -2 at high frequencies, usually with a single frequency break. More complex two-peak solutions also occur. The amplitude of the peaks and their frequency difference depend on the inherent timescales of the model. At intermediate frequencies, the intrinsic PSD is influenced by the individual light curve profile as well as by the type of the underlying process. However, even in cases when two Lorentzians seem to dominate the PSD, it does not necessarily imply that two single oscillation mechanisms operate simultaneously. Instead, it may well be the manifestation of the avalanche mechanism. The main advantage of our approach is an insight in the model functioning and the fast evaluation of the PSD.

We select a sample of 10 radio-quiet quasars with confirmed intranight optical variability and with available X-ray data. We compare the variability properties and the broad band spectral constraints to the predictions of intranight variability by three models: (i) irradiation of an accretion disk by a variable X-ray flux (ii) an accretion disk instability (iii) the presence of a weak blazar component. We concluded that the third model, e.g. the blazar component model, is the most promising if we adopt a cannonball model for the jet variable emission. In this case, the probability of detecting the intranight variability is within 20-80%, depending on the ratio of the disk to the jet optical luminosity. Variable X-ray irradiation mechanism is also possible but only under additional requirement: either the source should have a very narrow Hbeta line or occasional extremely strong flares should appear at very large disk radii.

Sgr A* is a source of strongly variable emission in several energy bands. It is generally agreed that this emission comes from the material surrounding the black hole which is either falling in or flowing out. The activity must be driven by accretion but the character of accretion flow in this object is an open question. We suggest that the inflow is dominated by the relatively low angular momentum material originating in one of the nearby group of stars. Such material flows in directly towards the black hole up to the distance of order of ten Schwarzschild radii or less, where it hits the angular momentum barrier which leads naturally to a flow variability. We study both the analytical and the numerical solutions for the flow dynamics, and we analyze the radiation spectra in both cases using the Monte Carlo code to simulate the synchrotron, bremsstrahlung and the Compton scattering. Our model roughly reproduces the broad band spectrum of Sgr A* and its variability if we allow for a small fraction of energy to be converted to non-thermal population of electrons. It is also consistent (for a range of viewing angles) with the strong constraints on the amount of circumnuclear material imposed by the measurements of the Faraday rotation.

There is increasing evidence that supermassive black holes in active galactic nuclei (AGN) are scaled-up versions of Galactic black holes. We show that the amplitude of high-frequency X-ray variability in the hard spectral state is inversely proportional to the black hole mass over eight orders of magnitude. We have analyzed all available hard-state data from RXTE of seven Galactic black holes. Their power density spectra change dramatically from observation to observation, except for the high-frequency (>10 Hz) tail, which seems to have a universal shape, roughly represented by a power law of index -2. The amplitude of the tail, C_M (extrapolated to 1 Hz), remains approximately constant for a given source, regardless of the luminosity, unlike the break or QPO frequencies, which are usually strongly correlated with luminosity. Comparison with a moderate-luminosity sample of AGN shows that the amplitude of the tail is a simple function of black hole mass, C_M = C/M, where C = 1.25 M_Sol / Hz. This makes C_M a robust estimator of the black hole mass which is easy to apply to low- to moderate-luminosity supermassive black holes. The high-frequency tail with its universal shape is an invariant feature of a black hole and, possibly, an imprint of the last stable orbit.

We present the analysis of multi-wavelength observations of bright quasar HS1603+3820: the optical data taken with the MMT and Keck telescopes, and X-ray data obtained with the Chandra X-ray Observatory. The optical spectra contain a very large number of absorption lines from numerous heavy elements. We derived X-ray properties of HS1603. The quasar has the optical-to-X-ray slope index alpha_ox of 1.70, which is on the high end of the typical range for radio quiet QSOs. We found 49 individual heavy element absorption clouds, which can be grouped into eleven distinct systems. We determined column densities and redshifts of the individual components. Absorbers from the associated system which is likely spatially closest to the QSO show large CIV to HI column density ratio, reaching ca.20.

[abridged] We model the X-ray reprocessing from a strong co-rotating flare above an accretion disk in active galactic nuclei. We explore the horizontal structure and evolution of the underlying hot spot. To obtain the spectral evolution seen by a distant observer, we apply a general relativity ray-tracing technique. We concentrate on the energy band around the iron K-line, where the relativistic effects are most pronounced. Persistent flares lasting for a significant fraction of the orbital time scale and short, transient flares are considered. In our time-resolved analysis, the spectra recorded by a distant observer depend on the position of the flare/spot with respect to the central black hole. If the flare duration significantly exceeds the light travel time across the spot, then the spot horizontal stratification is unimportant. On the other hand, if the flare duration is comparable to the light travel time across the spot radius, the lightcurves exhibit a typical asymmetry in their time profiles. The sequence of dynamical spectra proceeds from more strongly to less strongly ionized re-emission. At all locations within the spot the spectral intensity increases towards edge-on emission angles, revealing the limb brightening effect. Future X-ray observatories with significantly larger effective collecting areas will enable to spectroscopically map out the azimuthal irradiation structure of the accretion disk and to localize persistent flares. If the hot spot is not located too close to the marginally stable orbit of the black hole, it will be possible to probe the reflecting medium via the sub-structure of the iron K-line. Indications for transient flares will only be obtained from analyzing the observed lightcurves on the gravitational time scale of the accreting supermassive black hole.

Numerical simulations of MHD accretion flows in the vicinity of a supermasssive black hole provide important insights to the problem of why and how systems -- such as the Galactic Center -- are underluminous and variable. To access applicability of such flows to real objects, we examine the dynamical MHD studies with computations of the time dependent radiation spectra predicted by the simulations. We apply Monte Carlo methods to calculate spectra predicted by the time-dependent model of an axisymmetic MHD flow accreting onto a black hole presented by Proga and Begelman. Our calculations show that variability in an accretion flow is not always reflected in the corresponding spectra, at least not in all wavelengths. We find no one-to-one correspondence between the accretion state and the predicted spectrum

We show the locations of the SWIFT short hard bursts (SHB) with afterglows on the galactic map and compare with the VSB BATSE events. As we have pointed out before, there is an excess of events in the galactic map of BATSE VSB events. We not that none of VSB SWIFT era events fall into this cluster. More SWIFT events are needed to check this claim. We also report a new study with KONUS data of the VSB sample with an average energy above 90 keV showing a clear excess of events below 100 ms duration (T90) that have large mean energy protons. We suggest that VSB themselves consist of two subclasses: a fraction of events have peculiar distribution properties and have no detectable counterparts, as might be expected for exotic sources such as primordial black holes. We show how GLAST could add key new information to the study of VSB bursts and could help test the black hole concept.

We present modeling results for the reprocessed radiation expected from magnetic flares above AGN accretion disks. Relativistic corrections for the orbital motion of the flare and for the curved space-time in the vicinity of the black hole are taken into account. We investigate the local emission spectra, as seen in a frame co-orbiting with the disk, and the observed spectra at infinity. We investigate long-term flares at different orbital phases and short-term flares for various global parameters of the accreting black hole. Particular emphasis is put on the relation between the iron Kalpha line and the Compton hump as these two features can be simultaneously observed by the Suzaku satellite and later by Simbol-X.

We propose a reflection model of the time delays detected during an exceptionally bright, single flare in MCG-6-30-15. We consider a scenario in which the delays of the hard X-rays with respect to the soft X-rays are caused by the presence of the delayed reflection component. We employ a model of the flare, which is accompanied by reprocessed emission. We consider two geometries/thermal states of the reprocessing medium: a partially ionized accretion disk surface and a distribution of magnetically confined, cold blobs. The reprocessing by cold blobs predicts positive time delays and a saturation in the time delay -- energy relation, which is likely present in the data. The model requires a strong reflection component and relies on the apparent pivoting of the combined primary and reflected spectrum. The reflection by the ionized disk surface does not reproduce the observed delays. We discuss the relation between the two reflection scenarios and argue that they are both present in MCG-6-30-15.

We develop a two-flow model of accretion onto a black hole which incorporates the effect of the local magneto-rotational instability. The flow consists of an accretion disk and an accreting corona, and the local dynamo affects the disk/corona mass exchange. The model is aimed to explain the power spectrum density of the sources in their soft, disk-dominated states. The local perturbations of the magnetic field in the disk are described as in King et al. (2004) and Mayer and Pringle (2005), but the time-dependent local magnetic field is assumed to affect the local supply of the material to the corona. Since the viscous timescale in the corona is much shorter than in the disk, the local perturbations are not smeared in the corona. Simple analytical estimates and full time-dependent computations of the disk-corona system are performed. The accreting corona model can reproduce the broad power spectra of Soft State X-ray binaries and AGN. The model, however, predicts that (i) sources undergoing radiation pressure instability (GRS 1915+105) should have systematically steeper power spectra than other sources, (ii) AGN should have systematically steeper power spectra than GBH, even if their disks are described using viscosity proportional to the gas pressure. More precise measurements of power spectra of Soft State sources are clearly needed.

We discuss implications of a strong flare event observed in the Seyfert galaxy MCG-6-30-15 assuming that the emission is due to localized magnetic reconnection. We conduct detailed radiative transfer modeling of the reprocessed radiation for a primary source that is elevated above the disk. The model includes relativistic effects and Keplerian motion around the black hole. We show that for such a model setup the observed time-modulation must be intrinsic to the primary source. Using a simple analytical model we then investigate time delays between hard and soft X-rays during the flare. The model considers an intrinsic delay between primary and reprocessed radiation, which measures the geometrical distance of the flare source to the reprocessing sites. The observed time delays are well reproduced if one assumes that the reprocessing happens in magnetically confined, cold clouds.

Using Monte-Carlo simulations of X-ray flare distributions across the accretion disk of active galactic nuclei (AGN), we obtain modeling results for the energy-dependent fractional variability amplitude. Referring to previous results of this model, we illustrate the relation between the shape of the point-to-point fractional variability spectrum, F_pp, and the time-integrated spectral energy distribution, F_E. The results confirm that the spectral shape and variability of the iron Kalpha line are dominated by the flares closest to the disk center.

Recent modeling of the warm absorber in active galactic nuclei has proved the usefulness of constant total (gas plus radiation) pressure models, which are highly stratified in temperature and density. We explore the consistency of those models when the typical variation of the flux from the central source is taken into account. We perform a variability study of the warm absorber response, based on timescales and our photoionization code TITAN. We show that the ionization and recombination timescales are much shorter than the dynamical timescale. Clouds very close to the central black hole will maintain their equilibrium since the characteristic variability timescales of the nuclear source are longer than cloud timescales. For more distant clouds, the density structure has no time to vary, in response to the variations of the temperature or ionization structure, and such clouds will show the departure from the constant pressure equilibrium. We explore the impact of this departure on the observed properties of the transmitted spectrum and soft X-ray variability: (i) non uniform velocities, of the order of sound speed, appear due to pressure gradients, up to typical values of 100 km/s. These velocities lead to the broadening of lines. This broadening is usually observed and very difficult to explain otherwise. (ii) Energy-dependent fractional variability amplitude in soft X-ray range has a broader hump around ~ 1-2 keV, and (iv) the plot of the equivalent hydrogen column density vs. ionization parameter is steeper than for equilibrium clouds. The results have the character of a preliminary study and should be supplemented in the future with full time-dependent radiation transfer and dynamical computations.

The presence of the dust in the circumnuclear region strongly affects our view of the nucleus itself. The effect is strong in type 2 objects but weaker effect is likely to be present in type 1 objects as well. In these objects a correction to the observed optical/UV spectrum must be done in order to recover the intrinsic spectrum of a nucleus. The approach based on the extinction curve is convenient for that purpose so significant effort has been recently done in order to determine the extinction curve for the circumnuclear material. It seems clear that the circumnuclear dust is different from the average properties of the dust in the Interstellar Medium in our galaxy: the well known 2175 A feature is weak or absent in AGN nuclear dust, and the extinction curve at shorter wavelength does not seem to be rising as steeply. The circumnuclear dust is therefore more similar to SMC dust, or more likely, to the dust in very dense molecular clouds in our Galaxy. However, the exact shape of the extinction curve in the far UV is still a matter of debate, and various effects are difficult to disentangle.

We estimate black hole masses in Narrow Line Seyfert 1 (NLS1) galaxies at the basis of their X-ray excess variance. We apply the standard approach appropriate for Broad Line Seyfert 1 (BLS1) galaxies. In general, we find that the obtained masses are by a factor \sim 20 too small to agree with values obtained from other methods (reverberation, stellar dispersion). However, a small subset of our NLS1 objects does not require that multiplication, or the correction factor is less than 4. We find that this subset have a soft X-ray photon index, \Gamma_{0.1-2.4 keV}, smaller than 2. We thus postulate that this subclass of NLS1 actually belongs to BLS1.

We develop a model of time delays between the continuum bands in the Narrow Line Seyfert 1 galaxy Mrk 335 to explain the observed delays measured in this source. We consider two geometries: an accretion disk with fully ionized warm absorber of considerable optical depth, located close to the symmetry axis, and an accretion disk with a hot corona. Both media lead to significant disk irradiation but the disk/corona geometry gives lower values of the time delays. Only the disk/corona models give results consistent with measurements of Sergeev et al., and a low value of the disk inclination is favored. The presence of an optically thick, fully ionized outflow is ruled out at the 2-sigma level.

The X-ray spectra of Active Galactic Nuclei (AGN) are complex and vary rapidly in time as seen in recent observations. Magnetic flares above the accretion disk can account for the extreme variability of AGN. They also explain the observed iron Kalpha fluorescence lines. We present radiative transfer modeling of the X-ray reflection due to emission from magnetic flares close to the marginally stable orbit. The hard X-ray primary radiation coming from the flare source illuminates the accretion disk. A Compton reflection/reprocessed component coming from the disk surface is computed for different emission directions. We assume that the density structure remains adjusted to the hydrostatic equilibrium without external illumination because the flare duration is only a quarter-orbit. The model takes into account the variations of the incident radiation across the hot spot underneath the flare source. The integrated spectrum seen by a distant observer is computed for flares at different orbital phases close to the marginally stable orbit of a Schwarzschild black hole and of a maximally rotating Kerr black hole. The calculations include relativistic and Doppler corrections of the spectra using a ray tracing technique. We explore the practical possibilities to map out the azimuthal irradiation pattern of the inner accretion disks and conclude that the next generation of X-ray satellites should reveal this structure from iron Kalpha line profiles and X-ray lightcurves.

We show the locations of the SWIFT short hard bursts (SHB) with afterglows on the Galactic map and compare with the very short bursts (VSB) BATSE events. As we have pointed out before, there is an excess of events in the Galactic map of BATSE VSB events. We note, that none of VSB SWIFT era events fall into this cluster. More SWIFT events are needed to check this claim. We also report a new study with KONUS data of the VSB sample with an average energy above 90 keV showing a clear excess of events below 100 ms duration (T_90) that have large mean energy photons. We suggest that VSB themselves consists of two subclasses: a fraction of events have peculiar distribution properties and have no detectable counter parts, as might be expected for exotic sources such as Primordial Black Holes.

Following our previous work of Lachowicz & Czerny (2005), we explore further the application of the continuous wavelet transform to X-ray astronomical signals. Using the public archive of the XMM-Newton satellite, we analyze all available EPIC-pn observations for nearby Seyfert 1 galaxies MCG-6-30-15 and NGC 4051. We confine our analysis to 0.002-0.007 Hz frequency band in which, on the way of theoretically motivated premises, some quasi-periodic oscillations (QPOs) are expected to be found. We find that indeed wavelet power histogram analysis reveals such QPOs centered at two frequencies of ~2.5E-3 Hz and 4-6E-3 Hz, respectively. We show that these quasi-periodic features can be disentangled from the Poisson noise contamination level what is hardly to achieve with the standard Fourier analysis. Interestingly, we find some of them to be in 2:1 or 3:2 ratio. If real, our finding may be considered as a link between QPOs observed in AGN and kHz QPOs seen in X-ray binary systems.

Values of black hole masses are frequently determined with the help of the reverberation method. This method requires a specific geometrical factor related to the distribution of the orbits of the Broad Line Region clouds. Onken et al. determined the value f^2= 1.37+/-0.45 from the black hole mass - dispersion relation. In this paper we determine this factor using an independent mass determination from the X-ray variance method for a number of Seyfert 1 galaxies and comparing them with the reverberation results by Peterson et al. We obtain mean value f^2 = 1.12 +/- 0.54, consistent with Onken et al. Both values are larger than the value 0.75 corresponding to a spherical geometry. It indicates that most probably all values of the black hole masses obtained with the use of the Kaspi et al. formulae should be multiplied by a factor of \sim 1.7. This also shows that the Broad Line Region is rather flat, and hints for a dependence of the factor f^2 on a source inclination seem to be present in the data.

[abridged] Extending the work of Czerny et al. (2004), we model the fractional variability amplitude due to distributions of hot spots co-orbiting on the accretion disk around a supermassive black hole. From defined radial distributions, our code samples random positions for the hot spots across the disk. The local spot emission is computed as reprocessed radiation coming from a compact primary source above the disk. The structure of the hot spot and the anisotropy of the re-emission are taken into account. We compute the fractional variability spectra expected from such spot ensembles and investigate dependencies on the parameters describing the radial spot distribution. We consider the fractional variability F_var with respect to the spectral mean and also the so-called point-to-point F_pp. Our method includes relativistic corrections for the curved space-time; the black hole angular momentum is a free parameter and subject to the fitting procedure. We confirm that the rms-variability spectra involve intrinsic randomness at a significant level when the number of flares appearing during the total observation time is too small. Furthermore, F_var is not always compatible with F_pp. For MCG -6-30-15, we can reproduce the short-timescale variability and model the suppressed variability in the energy range of the Kalpha line without any need to postulate reprocessing farther away from the center. An increasing rate of energy production by the flares toward the center of the disk, a fast rotation of the central black hole, and considerable suppression of the primary flare emission are required. The modeled line remains consistent with the measured equivalent width of the iron Kalpha line complex.

The material accreting onto Sgr A* most probably comes from the nearby stars. We analyze the pattern of this flow at distances of a fraction of a parsec and we argue that the net angular momentum of this material is low but non-negligible, and the initially supersonic disk accretion changes into subsonic flow with constant angular momentum. Next we estimate the flow parameters at a distance $R_{BHL}$ from the black hole and we argue that for the plausible parameter range the accretion flow is non-stationary. The inflow becomes supersonic at distance of $\sim 10^4 R_g$ but the solution does not continue below the horizon and the material piles up forming a torus, or a ring, at a distance of a few up to tens of Schwarzchild radii. Such a torus is known to be unstable and may explain strong variability of the flow in Sgr A*. Our considerations show that the temporary formation of such a torus seems to be unavoidable. Our best fitting model predicts a rather large accretion rate of around $4 \cdot 10^{-6} M_{\odot}/yr$ directly on Sgr A*. We argue that magnetic fields in the flow are tangled and this allows our model to overcome the disagreement with the Faraday rotation limits.

The 2-10 keV continuum of AGN is generally well represented by a single power law. However, at smaller energies the continuum displays an excess with respect to the extrapolation of this power law, called the ''soft X-ray excess''. Until now this soft X-ray excess was attributed, either to reflection of the hard X-ray source by the accretion disk, or to the presence of an additional comptonizing medium, giving a steep spectrum. An alternative solution proposed by Gierlinski and Done (2004) is that a single power law well represents both the soft and the hard X-ray emission and the impression of the soft X-ray excess is due to absorption of a primary power law by a relativistic wind. We examine the advantages and drawbacks of reflection versus absorption models, and we conclude that the observed spectra can be well modeled, either by absorption (for a strong excess), or by reflection (for a weak excess). However the physical conditions required by the absorption models do not seem very realistic: we would prefer an ''hybrid model''.

The modified version of a bullet model for gamma ray bursts is studied. The central engine of the source produces multiple sub-jets that are contained within a cone. The emission of photons in the source frame of a sub-jet either takes part in an infinitesimally thin shell, or during its expansion for a finite time. The analysis of the observed profiles of GRBs taken by BATSE leads us to the conclusion that the latter possibility is much more favored. We also study the statistical distribution of GRBs, in the context of their bimodality of durations, taking into account the detector's capability of observing the signal above a certain flux limit. The model with shells emitting for a finite time is able to reproduce only one class of bursts, short or long, depending on the adopted physical parameters. Therefore we suggest that the GRB bimodality is intrinsically connected with two separate classes of sources.

We present a wide grid of models for the structure and transmission properties of warm absorbers in active galactic nuclei (AGN). Contrary to commonly used constant density models, our absorbing cloud is assumed to be under constant total (gas plus radiation) pressure. This assumption implies the coexistence of material at different temperatures and ionization states, which is a natural consequence of pressure and thermal equilibrium. Our photoionization code allows us to compute the profiles of the density, the temperature, the gas pressure, the radiation pressure and the ionization state across the cloud, and to calculate the radiative transfer of continuum and lines including Compton scattering. Therefore, equivalent widths of both saturated and unsaturated lines are properly modeled. For each pair of the incident spectrum slope and the ionization parameter at the cloud surface there is a natural upper limit to the total column densities of the cloud due to thermal instabilities. These maximum values are comparable to the observational constraints on the column density of warm absorbers which may give support to constant total pressure models. In all models we note considerable absorption around 6.4 keV which modifies the intrinsic relativistically broadened iron line profile originating in an accretion disk illuminated atmosphere. Our models can be applied to fitting the spectroscopic data from the {\it XMM-Newton} and {\it Chandra} satellites.

The 2-10 keV continuum of AGN is well represented by a single power law, generally attributed to a hot comptonizing medium, such as a corona above the accretion disk. At smaller energies the continuum displays an excess with respect to the extrapolation of this power law, called the ``soft X-ray excess". Until now it was attributed, either to reflection of the hard X-ray source by the accretion disk, or to the presence of an additional comptonizing medium. An alternative solution proposed by Gierli\'nski & Done (2004) is that a single power law represents correctly both the soft and the hard X-ray emission, and the soft X-ray excess is an artefact due to the absorption of the primary power law by a relativistic wind. We examine the advantages and drawbacks of the reflection versus absorption models. We argue that in the absorption hypothesis, the absorbing medium should be in total pressure equilibrium, to constrain the spectral distribution which otherwise would be too strongly variable in time and from one object to the other, as compared to observations. We conclude that some X-ray spectra, in particular those with strong soft X-ray excesses, can be modelled by absorption in the 0.3-10 keV range. However, due to the lack of a complete grid of models and good data extending above 10 keV, we are not able to conclude presently that all objects can be accommodated with such models. These absorption models imply either strong relativistic outflowing winds with mass rates of the order of the Eddington value (or even larger), or quasi-spherical inhomogeneous accretion flows. Only weak excesses can be modelled by reflection, unless the primary continuum is not directly seen. Finally, a reflection model absorbed by a modest relativistic wind could be the best solution to the problem.

Our recent studies of Very Short Gamma Ray Bursts (VSB) of duration T_{90} \le 100 ms have indicated a significant angular asymmetry and a uniform V/V_{max} distribution from the BATSE data. Here we update these studies, and we extend our research to events observed with KONUS satellite which gives a new insight into the spectra not possible with the BATSE data. KONUS has observed 18 events with T_{90} \le 100 ms duration. These events display considerable numbers of photons with energies above 1 MeV and have photons above 5 MeV in all cases. These appear to be some of the most energetic photons observed from any classes of GRB to date.

Active Galactic Nuclei are powered by accretion onto massive black holes. Although radio-quiet objects are not as spectacular sources of very high energy photons as radio-loud ones this class of objects also represents a challenge for modeling high energy processes close to a black hole. Both a hot optically thin plasma and a cooler optically thick accretion disk are usually thought to be present in the vicinity of a black hole although the details of the accretion flow are still under discussion. The role of the disk seems to decrease with a drop in the Eddington ratio: in sources like quasars and Narrow Line Seyfert 1 galaxies disk flow dominates while in Seyfert galaxies the disk retreats, and in sources like LINERS or Sgr A* a disk is most likely absent. Shocks and reconnections are possibly taking place in an inner hot flow and in the magnetic corona above the cold disk. Uncollimated outflow is also present and it may carry significant fraction of available mass and energy.

Short-term variability of X-ray continuum spectra has been reported for several Active Galactic Nuclei. Significant X-ray flux variations are observed within time scales down to 10^3-10^5 seconds. We discuss short variability time scales in the frame of the X-ray flare model, which assumes the release of a large hard X-ray flux above a small portion of the accretion disk. The resulting observed X-ray spectrum is composed of the primary radiation and of a reprocessed Compton reflection component that we model with numerical radiative transfer simulations. The incident hard X-rays of the flare will heat up the atmosphere of the accretion disk and hence induce thermal expansion. Eventually, the flare source will be surrounded by an optically thick medium, which should modify the observed spectra.

Application of multi-resolution analysis adopting wavelets in order to investigate millisecond aperiodic X-ray variability of Cyg X-1 is presented. This relatively new approach in a time-series analysis allows us to study both the significance of any flux oscillations localized in time-frequency space, as well as the stability of its period. Using the observations of the RXTE/PCA we analyze archival data from Cyg X-1 during its failed state transition. The power spectrum of such a state is strongly dominated by a single Lorentzian peak corresponding to a damped oscillator. Our wavelet analysis presented in this paper shows the existence of short lasting oscillations without an obvious trend in the frequency domain. Based on this result, we suggest the interpretation of the dominant components of Cyg X-1 power spectrum in all spectral states: (i) a power law component dominating a soft state is due to accretion rate fluctuations related to MRI instability in the cold disk, (ii) the low frequency Lorentzian is related to the same fluctuations but in the inner ion torus, and (iii) the high frequency Lorentzian (dominating the power spectrum during the failed state transition) is due to the dynamical pulsations of the inner ion torus. The usage of the wavelet analysis as a potential and attractive tool in order to spot the semi-direct evidences of accretion process onto black-holes is proposed.

Short-term variability of X-ray spectra has been reported for several Active Galactic Nuclei. Significant X-ray flux variations are observed within time scales down to 10^4 -- 10^5 seconds. We discuss short variability time scales in the framework of the X-ray flare model, which assumes the release of a large hard X-ray flux above a small portion of the accretion disk. The resulting observed X-ray spectrum is composed of the primary radiation as well as a reprocessed/Compton reflection component that we model performing numerical radiative transfer simulations. We conduct Monte-Carlo simulations of large flare distributions across the disk including relativistic corrections of the observed radiation and discuss general implications on AGN spectra and their variability.

We analyze the consequences of the hypothesis that the formation of the Broad Line Region is intrinsically connected with the existence of the cold accretion disk. We assume that the Broad Line Region radius is well estimated by the formula of Kaspi et al. (2000). We consider three models of the disappearance of the inner disk which limit the existence of the Broad Line Region: (i) classical ADAF approach, i.e. the inner hot flow develops whenever it can exist (ii) disk evaporation model of Meyer & Meyer-Hofmeister (2002) (iii) generalized disk evaporation model of Rozanska & Czerny (2000b). For each of the models, we determine the minimum value of the Eddington ratio and the maximum value of the broad line widths as functions of the viscosity parameter alpha and the magnetic field parameter beta. We compare the predicted parameter space with observations of several AGN. Weak dependence of the maximum value of the FWHM and minimum value of the Eddington ratio on the black hole mass in our sample is noticeable. It seems to favor the description of the cold disk/hot inner flow transition as in the classical ADAF approach than with the model of disk evaporation due to conduction between the disk and accreting corona.

The hard X-ray lightcurves exhibit delays of $\sim 1$ s with respect to the soft X-ray lightcurves when the microquasar GRS 1915+105 is in the state of frequent, regular outbursts (states $\rho$ and $\kappa$ of Belloni et al. 2000). Such outbursts are supposed to be driven by the radiation pressure instability of the inner disc parts. The hard X-ray delays are then caused by the time needed for the adjustment of the corona to changing conditions in the underlying disc. We support this claim by the computation of the time evolution of the disc, including a non-stationary evaporation of the disc and mass exchange with the corona.

Optically thick accretion disks are considered to be important ingredients of luminous AGN. The claim of their existence is well supported by observations and recent years brought some progress in understanding of their dynamics. However, the role of accretion disks in optical/UV/X-ray variability of AGN is not quite clear. Most probably, in short timescales the disk reprocesses the variable X-ray flux but at longer timescales the variations of the disk structure lead directly to optical/UV variations as well as affect, or even create, the X-ray variability pattern. We urgently need a considerable progress in time-dependent disk models to close the gap between the theory and the stream of data coming from the AGN monitoring.

The presence of the warm absorber of considerable optical depth is seen in many AGN. We show that this medium may affect significantly the optical/UV spectrum of an AGN by backscattering a fraction of the total radiation flux towards the disc surface. We consider in detail the case when the disk extends down to a marginally stable orbit, all the emission comes from the disk surface and the scattering medium forms a cone around the symmetry axis. Disc irradiation results in much flatter optical/UV continuum than predicted by standard disc models. The effect depends both on the total optical depth of the warm absorber and on the specific density distribution of this medium so the analysis of the optical/UV continuum allows to obtain constraints for the warm absorber complementary to those obtained from the soft X-ray data analysis. We give results for two exemplary sources - RE J1034+396 and PG1211+143, and for the bluest composite quasar spectrum of Richards et al. obtained from SDSS.

Rapid, hyper-Eddington accretion is likely to power the central engines of gamma-ray bursts (GRBs). In the extreme conditions of densities and temperatures the accreting torus is cooled by neutrino emission rather than by radiation. Another important cooling mechanism is the advection of energy into the central black hole. We compute the time evolution of a neutrino-dominated disc that proceeds during the burst and investigate the changes in its density and temperature. The discrimination between short and long bursts is made on the basis of the different rates of material inflow to the outer parts of the disc, thus favoring the binary merger scenario for the short GRBs and the collapsar scenario for the long ones. Within the context of the collapsar model, we also study the evolution of the photon luminosity of the remnant disc up to times of the order of 1 day, and discuss its implications for the production of emission lines in GRB spectra.

We discuss the local physical changes at the surface of an AGN accretion disk after the onset of a magnetic flare. The X-ray irradiation by a flare creates a hot spot at the disk surface where the plasma both heats up and expands in the vertical direction in order to regain the hydrostatic equilibrium. Assuming that the magnetic loop causing the flare is anchored deeply within the disk interior, we derive analytical estimates for the vertical dimension H_hot and the optical depth tau_es of the heated atmosphere as a function of the position within the spot. We perform computations for various values of the accretion rate dm/dt, the fraction f_cor of radiation dissipated within the disk corona, and the covering factor f_cover of the disk surface with flare-illuminated patches. It turns out that generally we can distinguish three characteristic radial zones within the disk showing a qualitatively different behavior of the heated material. In the innermost regions of the disk (inner zone) the expansion of the disk material is restricted by strong gravitational forces. Further out, the flare source, initially above the disk, soon becomes embedded by the expanding disk atmosphere. At these intermediate disk radii (middle zone) the material is optically thick thus greatly modifying the observed radiation by multiple Compton scattering. We show exemplary spectra models obtained from Monte Carlo simulations illustrating the trends. In the outermost regions of the disk (outer zone) the expanding material is optically thin and its influence on the observed spectra is smaller but pressure gradients in radial directions should cause the development of a fountain-like dynamical structure around the flare source. We discuss the observational consequences of our results.

We employ a low angular momentum accretion-outflow scenario to model the flares emanating out from the central region of Sgr A$^*$. The primary donor for matter accreting onto the central SMBH of Sgr A$^*$ is assumed to be the WR star ISR 13 E3. We analytically calculate the specific energy and angular momentum density of stellar wind originating from ISR 13 E3 and study the dynamics of that wind down to the very close vicinity of the central SMBH of Sgr A$^*$. We show that on the way to the Galactic centre, such wind-fed accretion may encounter standing shocks and such shock drives outflow from the close vicinity of the SMBH. Matter content of such outflow is computed and it is argued that such outflow is responsible for production of the Galactic centre flares. We then self-consistently compute the luminosity ${\cal L}_j$ and the duration time scale $\tau$ of such flares, as a function of fundamental accretion parameters. Our theoretical calculation of ${\cal L}_j$ and $\tau$ are in good agreement with observational results.

We propose a new method of estimation of the black hole masses in AGN based on the normalized excess variance, sigma_{nxs}^2. We derive a relation between sigma_{nxs}^2, the length of the observation, T, the light curve bin size, Delta t, and the black hole mass, assuming that (i) the power spectrum above the high frequency break, f_{bf}, has a slope of -2, (ii) the high frequency break scales with black hole mass, (iii) the power spectrum amplitude (in 'frequency x power' space) is universal and (iv) sigma_{nxs}^2 is calculated from observations of length T < 1/f_{bf}. Values of black hole masses in AGN obtained with this method are consistent with estimates based on other techniques such as reverberation mapping or the Mbh-stellar velocity dispersion relation. The method is formally equivalent to methods based on power spectrum scaling with mass but the use of the normalized excess variance has the big advantage of being applicable to relatively low quality data.

We discuss a model of X-ray variability of active galactic nuclei (AGN). We consider multiple spots which originate on the surface of an accretion disk following intense irradiation by coronal flares. The spots move with the disk around the central black hole and eventually decay while new spots continuously emerge. We construct time sequences of the spectra of the spotted disk and compute the corresponding energy-dependent fractional variability amplitude. We explore the dependence on the disk inclination and other model parameters. AGN seen at higher inclination with respect to the observer, such as Seyfert 2 galaxies, are expected to have fractional variability amplitude of the direct emission by a factor of a few higher than objects seen face on, such as the Seyfert 1s.

We construct a quasar extinction curve based on the blue and red composite quasar spectra of Richards et al. (2003) prepared from the SDSS survey. This extinction curve does not show any traces of the 2200 A feature characteristic of the Interstellar Medium, and this indicates that graphite grains are likely absent close to quasar nuclei. The extinction is best modeled by AC amorphous carbon grains, assuming a standard distribution of grain sizes (p=3.5) but slightly larger minimum grain size (a_{min} = 0.016 microns) and lower maximum grain size (a_{min} = 0.12 microns) than the respective canonical values for the interstellar medium. The dust composition is thus similar to that of the dust in AGB stars. Since graphite grains form from amorphous carbon exposed to strong UV irradiation the results indicate that either the dust forms surprisingly far from the active nucleus or in a wind that leaves the nucleus quickly enough to avoid crystallization into graphite.

We consider effects of the MHD turbulence on the viscosity during the evolution of the thermal-viscous ionization instability in the standard $\alpha$-accretion disks. We consider the possibility that the accretion onto a supermassive black hole proceeds through an outer standard accretion disk and inner, radiatively inefficient and advection dominated flow. In this scenario we follow the time evolution of the accretion disk in which the viscosity parameter $\alpha$ is constant throughout the whole instability cycle, as implied by the strength of MHD turbulence. We conclude that the hydrogen ionization instability is a promising mechanism to explain the intermittent activity in AGN.

We present new interpretation of Ton S180 spectrum obtained by {\it Chandra} Spectrometer (Low Energy Transmission Grating). Several narrow absorption lines and a few emission disk lines have been successfully fitted to the data. We have not found any significant edges accompanying line emission. We propose the interpretation of narrow lines consistent with the paper recently written by Krolik (2002), where warm absorber is strongly inhomogeneous. Such situation is possible in so called multi-phase medium, where regions with different ionization states, densities and temperatures may coexist in thermal equilibrium under constant pressure. We illustrate this scenario with theoretical spectra of radiation transfered through a stratified cloud with constant pressure (instead of constant density) computed by code {\sc titan} in plane parallel approximation. Detected spectral features are faint and their presence do not alter the broad band continuum. We model the broad band continuum of Ton S180 assuming an irradiated accretion disk with a dissipative warm skin. The set of parameters appropriate for the data cannot be determined uniquely but models with low values of the black hole mass have too hot and radially extended warm skin to explain the formation of soft X-ray disk lines seen in the data.

The spectra of radio quiet and NLS1 galaxies show suprising similarity in their shape. They seem to scale only with the accretion rate but not with central black hole mass. We consider two mechanisms modifying the disk spectrum. First, the outer parts of the disk are irradiated by the flux emerging from the inner parts. This is due to the scattering of the flux by the extended hot medium (warm absorber). Second process is connected with the development of the disk warm Comptonizing skin above the disk and/or coronae. Our scenario applies only to object with relatively high luminosity to the Eddington luminosity ratio for which disk evaporation is inefficient.

The spectra of quasars and NLS1 galaxies show surprising similarity in their spectral shape. They seem to scale only with the accretion rate. This is in contradiction with the simple expectations from the standard disk model which predicts lower disk temperature for higher black hole mass. Here we consider two mechanisms modifying the disk spectrum: the irradiation of the outer disk due to the scattering of the flux by the extended ionized medium (warm absorber and the development of the warm Comptonizing disk skin under the effect of the radiation pressure instability. Those two mechanisms seem to lead to a spectrum which indeed roughly scales, as observed, only with the accretion rate. The scenario applies only to objects with relatively high luminosity to the Eddington luminosity ratio for which disk evaporation is inefficient.

We analyze observations of the Seyfert galaxy NGC 4151 covering 90 years in the optical band and 27 years in the 2-10 keV X-ray band. We compute the Normalized Power Spectrum Density (NPSD), the Structure Function (SF) and the Autocorrelation Function (ACF) for these data. The results show that the optical and X-ray variability properties are significantly different. X-ray variations are predominantly in the timescale range of 5 - 1000 days. The optical variations have also a short timescale component which may be related to X-ray variability but the dominant effect is the long timescale variability, with timescales longer than $\sim$ 10 years. We compare our results with observations of NGC 5548 and Cyg X-1. We conclude that the long timescale variability may be caused by radiation pressure instability in the accretion disk, although the observed timescale in NGC 4151 is by a factor of few longer than expected. X-ray variability of this source is very similar to what is observed in Cyg X-1 but scaled with the mass of the black hole, which suggests that the radiation pressure instability does not affect considerably the X-ray production.

The time dependent evolution of the accretion disk around black hole is computed. The classical description of the $\alpha$-viscosity is adopted so the evolution is driven by the instability operating in the innermost radiation-pressure dominated part of the accretion disk. We assume that the optically thick disk always extends down to the marginally stable orbit so it is never evacuated completely. We include the effect of the advection, coronal dissipation and vertical outflow. We show that the presence of the corona and/or the outflow reduce the amplitude of the outburst. If only about half of the energy is dissipated in the disk (with the other half dissipated in the corona and carried away by the outflow) the outburst amplitude and duration are consistent with observations of the microquasar GRS 1915+105. Viscous evolution explains in a natural way the lack of direct transitions from the state C to the state B in color-color diagram of this source. Further reduction of the fraction of energy dissipated in the optically thick disk switches off the outbursts which may explain why they are not seen in all high accretion rate sources being in the Very High State.

We observed 17 optically-selected, radio-quiet high-redshift quasars with the Chandra Observatory ACIS, and detected 16 of them. The quasars have redshift between 3.70 and 6.28 and include the highest redshift quasars known. When compared to low-redshift quasars observed with ROSAT, these high redshift quasars are significantly more X-ray quiet. We also find that the X-ray spectral index of the high redshift objects is flatter than the average at lower redshift. These trends confirm the predictions of models where the accretion flow is described by a cold, optically-thick accretion disk surrounded by a hot, optically thin corona, provided the viscosity parameter alpha >= 0.02. The high redshift quasars have supermassive black holes with masses ~10^{10} M_{sun}, and are accreting material at ~0.1 the Eddington limit. We detect 10 X-ray photons from the z=6.28 quasar SDS 1030+0524, which may have a Gunn-Peterson trough and be near the redshift of reionization of the intergalactic medium. The X-ray data place an upper limit on the optical depth of the intergalactic medium tau(IGM) < 10^6, compared to the lower limit from the spectrum of Lyalpha and Lybeta, which implies tau(IGM) > 20.

The subject of the paper is the role of the massive stellar cluster in the activity phenomenon and in the structure of active galactic nuclei. We introduce a simple model of stellar dynamics in the internal part of the cluster, which allows us to include both the star-disk and the star-star interactions. It is shown that the properties of the distribution of stars in the vicinity of the black hole are determined both by the interaction of the stars with the accretion disk and by the pair gravitational and contact interaction between the stars. We calculate the distribution of stars in the central parts of the cluster and we discuss possible effects of stellar mass-loss due to the star-disk interaction. Finally, we study the implications of the central cluster for active galactic nuclei activity. We model the broad line region assuming that the gaseous wakes, following stars after each disk crossing, play the role of the broad line region clouds, and we calculate the corresponding line profiles. We also analyze the contribution of star-star and star-disk collisions to active galactic nuclei variability.

We present a model of the vertical stratification and the spectra of an irradiated medium under the assumption of constant pressure. Such a solution has properties intermediate between constant density models and hydrostatic equilibrium models, and it may represent a flattened configuration of gas clumps accreting onto the central black hole. Such a medium develops a hot skin, thicker than hydrostatic models, but thinner than constant density models, under comparable irradiation. The range of theoretical values of the alpha_ox index is comparable to those from hydrostatic models and both are close to the observed values for Seyfert galaxies but lower than in quasars. The amount of X-ray Compton reflection is consistent with the observed range. The characteristic property of the model is a frequently multicomponent iron K alpha line.

We present detailed computations of the vertical structure of an accretion disc illuminated by hard X-ray radiation with the code {\sc titan-noar} suitable for Compton thick media. The energy generated via accretion is dissipated partially in the cold disc as well as in the X-ray source. We study the differences between the case where the X-ray source is in the form of a lamp post above the accretion disc and the case of a heavy corona. We consider radiative heating via Comptonization together with heating via photo-absorption on numerous heavy elements as carbon, oxygen, silicon, iron. The transfer in lines is precisely calculated. A better description of the heating/cooling through the inclusion of line transfer, a correct description of the temperature in the deeper layers, a correct description of the entire disc vertical structure, as well as the study of the possible coronal pressure effect, constitute an improvement in comparison to previous works. We show that exact calculations of hydrostatic equilibrium and determination of the disc thickness has a crucial impact on the optical depth of the hot illuminated zone. We assume a moderate illumination where the viscous flux equals the X-ray radiation flux. A highly ionized skin is created in the lamp post model, with the outgoing spectrum containing many emission lines and ionization edges in emission or absorption in the soft X-ray domain, as well as an iron line at $\sim 7 $ keV consisting of a blend of low ionization line from the deepest layers and hydrogen and helium like resonance line from the upper layers, and almost no absorption edge, contrary to the case of a slab of constant density.A full heavy corona completely suppresses the highly ionized zone on the top of the accretion disc and in such case the spectrum is featureless.

We propose a theoretical model to explain the spectrum of the quasar PG1211+143 emitted in the Optical/X-ray bands. In particular, we suggest that the inner accretion disk may develop a warm, optically thick skin, which produces a profound emission feature observed in the soft X-ray band. This is well modelled with the Comptonized black body emission. The same warm, mildly ionized medium may also be responsible for the hard X-ray reflection and the presence of the iron $K_{\alpha}$ line. However, in our model it still remains an open question, whether the seed photons for Comptonization come from the cold accretion disc or from the hotter plasma. High resolution spectroscopy available through the Chandra and XMM data may provide now an independent test of the physical conditions in the Comptonizing and reflecting warm skin.

The central engine of the the Narrow Line Seyfert 1 galaxies is being probed. We use the ASCA and RXTE data to model the X-ray primary continuum as well as the reflected component and iron $K\alpha$ line. Since these are strongly coupled, we obtain independent measurements of the disc ionization level and the orientation dependent reflection amplitude. Using the available Optical/UV data we also estimate the black hole masses and the $L/L_{Edd}$ ratios, which are probably related to Boroson & Green eigenvector~1.

We present the results of computations performed with a new photoionisation-transfer code designed for hot Compton thick media. Presently the code solves the transfer of the continuum with the Accelerated Lambda Iteration method (ALI) and that of the lines in a two stream Eddington approximation, without using the local escape probability formalism to approximate the line transfer. We show the influence that the approximations made in the treatment of the transfer and of the atomic data can have on the broad band UV-X spectrum and on the detailed spectral features (X-ray lines and ionization edges) emitted/absorbed by an X-ray irradiated medium. This transfer code is coupled with a Monte Carlo code which allows to take into account direct and inverse Compton diffusions, and to compute the spectrum emitted up to hundreds of keV energies. The influence of a few physical parameters is shown, and the importance of the density and pressure distribution (constant density, pressure equilibrium) is discussed.

We present the study of the emission properties of the quasar PG1211+143, which belongs to the class of Narrow Line Seyfert 1 galaxies. On the basis of observational data analyzed by us and collected from the literature, we study the temporal and spectral variability of the source in the optical/UV/X-ray bands and we propose a model that explains the spectrum emitted in this broad energy range. In this model, the intrinsic emission originating in the warm skin of the accretion disk is responsible for the spectral component that is dominant in the softest X-ray range. The shape of reflected spectrum as well as Fe K line detected in hard X-rays require the reflecting medium to be mildly ionized (xi~500). We identify this reflector with the warm skin of the disk and we show that the heating of the skin is consistent with the classical alpha P_{tot} prescription, while alpha P_{gas} option is at least two orders of magnitude too low to provide the required heating. We find that the mass of the central black hole is relatively small (M_BH~10^7- 10^8 Msun, which is consistent with the Broad Line Region mapping results and characteristic for NLS1 class.

Mathur (2000) has proposed that the narrow line Seyfert 1 galaxies (NLS1s) are likely to be the active galaxies in early stage of evolution. To test this hypothesis, we have calculated the black hole (BH) masses and the host galaxy bulge masses for a sample of NLS1s. We find that the mean BH to bulge mass ratio of NLS1s is significantly smaller than that for normal Seyfert galaxies. We also find that the ratio of BH mass to bulge velocity dispersion is also significantly smaller for NLS1s. A scenario of BH growth is our preferred interpretation, though alternative explanations are discussed. Assuming that the BHs grow with accretion with a radiative efficiency of 0.1, it will take them about $3.3\times 10^8$ years to become as massive as in normal Seyfert 1s. These timescales are consistent with the theoretical estimates of quasar timescale $t_Q$ less than about $4.5\times 10^8$ years calculated by Haehnelt, Natarajan & Rees (1998). Studies of low redshift NLS1s thus provide a powerful, and due to their proximity relatively easy, way to understand the high redshift quasars and their evolution.

We analyze the scaling of the X-ray power density spectra with the mass of the black hole on the example of Cyg X-1 and Seyfert 1 galaxy NGC 5548. We show that the high frequency tail of the power density spectrum can be successfully used for determination of the black hole mass. We determine the masses of the black holes in 6 Broad Line Seyfert 1 galaxies, 5 Narrow Line Seyfert 1 galaxies and two QSOs using available power density spectra. The proposed scaling is clearly appropriate for other Seyfert galaxies and QSOs. In all but 1 normal Seyferts the resulting luminosity to the Eddington luminosity ratio is smaller than 0.15, with a source MCG -6-15-30 being an exception. The applicability of the same scaling to Narrow Line Seyfert 1 is less clear and there may be a systematic shift between the power spectra of NLS1 and S1 galaxies of the same mass, leading to underestimation of the black hole mass. However, both the method based on variability and the method based on spectral fitting show that those galaxies have relatively low masses and high luminosity to the Eddington luminosity ratio, supporting the view of those objects as analogs of galactic sources in their high/soft or very high state based on the overall spectral shape. Bulge masses of their host galaxies are similar to normal Seyfert galaxies so they do not follow the black hole mass-bulge mass relation for Seyfert galaxies, being evolutionary less advanced, as suggested by Mathur (2000). The bulge mass-black hole mass relation in our sample is consistent with being linear, with black hole to bulge ratio $\sim$ 0.03 %, similar to Wandel (1999) and Laor (1998, 2001) for low mass objects but significantly shifted from the relation of Magorrian et al. (1998) and McLure & Dunlop (2000).

We test the disk/corona model of Janiuk et al. (2000) against the composite mean optical/UV/soft X-ray spectrum of radio loud and radio quiet quasars from Laor et al. (1997) which applies to faint quasars. The model well represents the optical/UV continuum if the hardening of the locally emitted spectrum is included, with the color to effective temperature ratio $f \sim 2$ in the inner 20-30 Schwarzschild radii. Comptonization seen in the soft X-ray band is well explained by the adopted corona model in the case of radio loud objects. However, this Comptonization is much stronger in radio quiet objects and additional ad hoc assumed Comptonizing medium must be present. We speculate that perhaps in both types of quasars there is a strong outflow of the hot plasma; this plasma is collimated in radio loud objects but not collimated in the radio quiet objects so it is present along the line of sight, serving as the required Comptonizing medium. Hard X-ray power law is not explained by the model and it may come from non-thermal component of electron plasma, as it is the case for galactic black holes in their soft states.

Physical mechanism responsible for high viscosity in accretion disks is still under debate. Parameterization of the viscous stress as $\alpha P$ proved to be a successful representation of this mechanism in the outer parts of the disk, explaining the dwarf novae and X-ray novae outbursts as due to ionization instability. We show that this parameterization can be also adopted in the innermost part of the disk where the adoption of the $\alpha$-viscosity law implies the presence of the instability in the radiation pressure dominated region. We study the time evolution of such disks. We show that the time-dependent behavior of GRS 1915+105 can be well reproduced if $\alpha$-viscosity disk model is calculated accurately (with proper numerical coefficients in vertically averaged equations and with advection included), and if the model is supplemented with (i) moderate corona dissipating 50% of energy (ii) jet carrying luminosity-dependent fraction of energy. These necessary modifications in the form of the presence of a corona and a jet are well justified observationally. The model predicts outbursts at luminosity larger than 0.16$\dot M_{Edd}$, as required, correct outburst timescales and amplitudes, including the effect of increasing outburst timescale with mean luminosity. This result strongly suggests that the $\alpha$-viscosity law is a good description of the actual mechanism responsible for angular momentum transfer also in the innermost, radiation pressure dominated part of the disk around a black hole.

The standard accretion discs are known to be thermally and viscously unstable over a certain range of temperatures. In the inner disc regions there may develop radiation pressure driven instability, which is possibly related to the rapid variability detected in AGNs in the UV range. In the outer disc develops the ionization instability, similar to that in the cataclysmic variables, but operating on much longer timescales. Due to this process the spectrum of the accretion disc differs from that of a stationary one.

We present the results of both analytical and numerical calculations of the amplitude of the reflection component in X-ray spectra of galactic black hole systems. We take into account the anisotropy of Compton scattering and the systematic relativistic bulk motion of the hot plasma. In the case of single scattering approximation the reflection from the disc surface is significantly enhanced due to the anisotropy of Compton scattering. On the other hand the calculations of multiple scattering obtained using the Monte Carlo method show that the anisotropy effect is much weaker in that case. Therefore, the enhanced backscattered flux may affect the observed spectra only if the disc surface is highly ionized, which reduces the absorption in the energy band corresponding to the first Compton scattering.

We study the UV spectra of narrow-line Seyfert~1 (NLSy1) galaxies and compare them with ``normal'' AGN. Similar to their optical lines, the NLSy1s show narrower UV lines. They are also characterized by weaker CIV 1549, CIII] 1909, and stronger AlIII 1857 emission. These UV line properties add to the optical and X-ray properties known to be part of the Boroson and Green eigenvector 1. We show that the steep soft-X-rays, which characterize the NLSy1s SEDs, change the equilibrium of the two phase cloud-intercloud medium resulting in somewhat higher BLR cloud densities, lower ionization parameter and larger BLR radii. These modified conditions can explain the unusual emission line properties we find in NLSy1. Using a specific model of an accretion disk with corona presented by Witt, Czerny and Zycki (1997), we also show that the steep soft and hard-X-ray continua can be explained if the $L/L_{Edd}$ ratios are larger than in ``normal'' Seyfert1/QSO strengthening earlier suggestions that the $L/L_{Edd}$ is the physical parameter driving this eigenvector.

We investigate the model of the disc/corona accretion flow around the black hole. The model, parameterized by the total (i.e. disk plus corona) accretion rate, $\dot m$, mass of the black hole, $M$ and the viscosity parameter, $\alpha$, uniquely determines the fraction of energy released in the corona as a function of radius and, in particular, the transition radius to the single-phase ADAF flow. Self-consistent solutions with the mass exchange between phases display radial dependence of the parameters. This model reproduces all characteristic luminosity states of accretion black hole without any additional ad hoc assumptions. In particular. the mechanism of the disk evaporation leads to a new, almost horizontal branch on the accretion flow's stability curve (i.e. the dependence of accretion rate on surface density) at the critical accretion rate.

The origin and the profile of the iron K-alpha line from active galactic nuclei is being studied. The model with a quasi-spherical distribution of clouds around a compact galactic nucleus is explored in which the intrinsic spectral feature around 6-7 keV is substantially affected by Comptonization of primary X-rays. The predicted spectral profiles are influenced mainly by motion of the clouds in strong gravitational field of the centre, and by orientation of the source with respect to the observer. The case of the Seyfert galaxy MCG-6-30-15 is examined in more detail and free parameters of the model are determined. A robust result is obtained: acceptable fits to the observed profile correspond to the spherical distribution of the clouds with strong self-obscuration effects and with rather high value of ionization parameter. This is a different situation from that considered in the model of a cold illuminated disc. Nevertheless, geometrically flat (disc-type) configurations are not rejected.

We study the UV spectra of NLS1 galaxies and compare them with typical Seyfert 1 galaxies and quasars. The NLS1 spectra show narrower UV lines as well as weaker CIV lambda 1549 and CIII] lambda 1909 emission. We show that these line properties are due to a lower ionization parameter and somewhat higher BLR cloud densities. These modified conditions can be explained by the hotter big blue bumps observed in NLS1s, which are in turn due to higher L/L_Edd ratios, as shown by our accretion disk and corona modeling of the NLS1 continua. We also present evidence that the Boroson & Green eigenvector 1, which is correlated with the optical and UV emission-line properties, is not driven by orientation and hence NLS1s, which have extreme eigenvector 1 values, are not viewed from an extreme viewing angle.

We present the results of spectral analysis of the ASCA data for the Narrow-Line Seyfert 1 galaxy (NLS1) Ton S180 and simultaneous ASCA and RXTE data modelling for the NLS1 Ark 564. We model both the primary and reflected continuum as well as the iron K alpha line, the energy of which depends on the ionization state of the reprocessor. We show that the reprocessing matter is mildly ionized, and we find the soft to hard luminosity ratio to be about 2.5. The accretion rate approximately corresponds to the Eddington limit value.

We propose a model of the accretion flow onto a black hole consisting of the accretion disk with an accreting two-temperature corona. The model is based on assumptions about the radiative and conductive energy exchange between the two phases and the pressure equilibrium. The complete model is determined by the mass, the accretion rate, and the viscosity parameter. We present the radial dependencies of parameters of such a two-phase flow, with advection in the corona and the disk/corona mass exchange due to evaporation/condensation included, and we determine the transition radius from a two-phase disk/corona accretion to a single-phase optically thin flow (ADAF) in the innermost part of the disk as a function of accretion rate. We identify the NLS1 galaxies with objects accreting at a rate close to the Eddington accretion rate. The strong variability of these objects may be related to the limit cycle behaviour expected in this luminosity range, as the disk, unstable due to the dominance by the radiation pressure, oscillates between the two stable branches: the advection-dominated optically thick branch and the evaporation branch.

We analyse the radiation pressure instability of the standard, optically thisc disc and propose it as an explanation of the obseved variability of the superluminal microquasar GRS 1915+105. The quantitave agreement with the observed outburst amplitudes is obtained under the assumption that substantial part of the energy is carried away by the jet.

The main characteristics of the average spectrum of radio quiet AGN in the UV and X-ray range are reviewed, and the emission mechanisms are discussed in the framework of accretion disk models, in particular the ``irradiated cold relativistic disk". It is shown that some problems arise in confronting the predictions of the model to the observations. We propose an alternative model in terms of a hot disk surrounded by a cold Compton thick medium. Finally we mention problems with remote regions of the accretion disk.

We propose an explanation of the variability of AGN based on a cloud model of accretion onto a black hole. These clouds, present at a distance of 10-100 Schwarzschild radii, possibly come from the disrupted innermost disk and they can partially obscure the central X-ray source due to their extreme Compton optical depth. We consider the implications of this scenario using a toy model of AGN spectra and support the results with full radiative transfer computations from the codes TITAN and NOAR. We show that small random rearrangements of the cloud distribution can happen on timescales of the order of 10^2 - 10^6 s and may lead to relatively high variability amplitude in the X-ray band if the mean covering factor is large. The normalized variability amplitude in the UV band is either the same as in the X-ray band, if the contribution from the dark sides of the clouds to the UV band is negligible, or it is smaller. The X-ray spectrum should basically preserve its spectral shape even during the high amplitude variability, within the frame of our model.

We study in detail the properties of the accreting corona model for active galactic nuclei and galactic black holes. In this model the fraction of the energy liberated in the corona at a given radius is calculated from the global parameters of the model (mass of the central object, accretion rate and viscosity parameter) and it appears to be a strong function of the radius. The model predicts the relative decrease of the coronal hard X-ray emission with an increase of the accretion rate. The presented description of disc/corona interaction forms a basis for further studies of disc disruption mechanism.

We consider the role of advection in the two-temperature accreting corona with an underlying optically thick disc. The properties of coronal solutions depend significantly on the description of advection. Local parameterization of advection by a constant coeficient $\delta$ replacing the radial derivatives lead to complex topology of solutions, similar to some extent to other advection-dominated accretion flow solutions. One, radiatively cooled branch exists for low accertion rates. For higher accretion rates two solutions exist in a broad range of radii: one is radiatively cooled and the other one is advection-dominated. With further increase of accretion rate the radial extensions of the two solutions shrink and no solutions are found above certain critical value. However, these trends change if the local parameterization of advection is replaced by proper radial derivatives computed iteratively from the model. Only one, radiatively cooled solution remains, and it exists even for high accretion rates. The advection-dominated branch disappears during the iteration process which means that a self-consistently described advection-dominated flow cannot exist in the presence of an underlying cold disc.

There are strong suggestions that the disk-like accretion flow onto massive black hole in AGN is disrupted in its innermost part (10-100 Rg), possibly due to the radiation pressure instability. It may form a hot optically thin quasi spherical (ADAF) flow surrounded by or containing denser clouds due to the disruption of the disk. Such clouds might be optically thick, with a Thompson depth of order of 10 or more. Within the frame of this cloud scenario (Collin-Souffrin et al. 1996, Czerny & Dumont 1998), obscuration events are expected and the effect would be seen as a variability. We consider expected random variability due to statistical dispersion in location of clouds along the line of sight for a constant covering factor. We discuss a simple analytical toy model which provides us with the estimates of the mean spectral properties and variability amplitude of AGN, and we support them with radiative transfer computations done with the use of TITAN code of Dumont, Abrassart & Collin (1999) and NOAR code of Abrassart (1999).

We study vertical structure of radiation pressure dominated disc with a hot corona. We include all the relevant processes like bound--free opacity and convection. We show that the presence of the corona modifies considerably the density and the opacity of the disc surface layers which are important from the point of view of spectrum formation. The surface of the disc with a corona is much denser and less ionized than the surface of a bare disc. Such a disc is likely to produce a neutral reflection and a local spectrum close to a black body. This effect will help to reconcile the predictions of accretion disc models with the observational data since a neutral reflection and a lack of Lyman edge are generally seen in AGN.

Recently Ghosh (1998) reported a new regime of instability in Keplerian accretion disks which is caused by relativistic effects. This instability appears in the gas pressure dominated region when all relativistic corrections to the disk structure equations are taken into account. We show that he uses the stability criterion in completely wrong way leading to inappropriate conclusions. We perform a standard stability analysis to show that no unstable region can be found when the relativistic disk is gas pressure dominated.

We derive power density spectra in the optical and X-ray band in the timescale range from several years down to $\sim $ a day. We suggest that the optical power density spectrum consists of two separate components: long timescale variations and short timescale variations, with the dividing timescale around 100 days. The shape of the short timescale component is similar to X-ray power density spectrum which is consistent with the interpretation of short timescale optical variations being caused by X-ray reprocessing. We show that the observed optical long timescale variability is consistent with thermal pulsations of the accretion disc.

We assume a distribution of clouds optically thick for electron scattering (OTCM) which are moderately optically thin for absorption and we consider them as a model of the mean quasar spectrum of Laor et al. (1997). We show that the model is particularly sensitive to the value of the ionization parameter $\xi$ and that for $\xi \sim 500$ the model well reproduces the optical/UV/X-ray mean quasar spectrum, in agreement with the estimates of the ionization parameter based on the energy of the iron $K_{\alpha}$ line. We cannot definitively reject synchrotron emission as a source of primary radiation but we favor the model in which the hard X-ray emission is produced by Compton scattering of soft photons in a central hot medium surrounded by cool clouds. In such a model clouds are located typically at the distance of $\sim 12 R_{Schw}$, with the covering factor about 0.88 and the radius of hot plasma is $\sim 9 R_{Schw}$. The model explains optical/UV emission as predominantly due to the dark sides of the clouds and the soft X-ray emission as due to the reflection by the irradiated sides of the clouds. Therefore, atomic features are expected in these bands although they are hardly present in the observational data. The kinematical effects connected with the cloud motion affect those features but do not remove them. The level of the primary emission required to model the mean quasar spectrum is too low to reproduce the equivalent width of the iron K$_{\alpha}$ line correctly but more detailed computations may resolve this problem.

We study the extensively monitored Seyfert galaxy NGC 5548 and compare its nuclear emission with models of accretion disk with accreting corona. To obtain the intrinsic nuclear spectra from the observed spectra we had to estimate and subtract the contribution from circumnuclear components such as stars, the Balmer continuum and blended FeII lines, and the FC2 extended, featureless continuum. The true nuclear spectra were compared with a two parameter model of the accreting disk with an accreting corona, described by the mass of the central black hole and viscosity. The model that best fitted the data was for M_{BH}=1.4e8 solar masses and the viscosity parameter alpha=0.033. Such a low viscosity parameter was necessary to produce the sufficient amount of X-rays. The vertical outflow of mass from corona in the form of wind had to be neglected in our model in order to fit into high and low states that NGC 5548 underwent. The model also predicts the behavior of the overall opt/UV/X continuum of NGC 5548 during the whole five year campaign.

The accurate position determination of GRB 970228 by the Beppo/SAX Satellite led to the discovery of a fading X-ray and optical counterparts of the burst. About a month after the GRB event the proper motion of the ejecta was detected and the extended optical source has faded below the Keck detection level. We analyze these observations in the framework of the most popular gamma ray burst models. We find a lower limit on the distance to GRB 970228 d>6.6 pc. We estimate the amount of energy required to explain the motion of a point source and the variability of the extended source. We find that the cosmological models are ruled out. The Galactic corona models suffer difficulties if the distance to the burst is > 100 kpc or the variability of the extended source is directly connected with the burst No constraints are found for GRBs originating between 6.6 pc and $\sim 10$ kpc.

We study a model of an accretion disc surrounded by an extended corona of the temperature 10^7 - 10^8 K. This corona modifies the disk spectrum since it redirects a significant fraction of the emission from the central parts towards the more distant parts of the disk. The same corona is indirectly the source of the broad emission lines because we expect the formation of cool clouds at the basis of the corona due to thermal instabilities. We model the number of the clouds and their motion through the corona adopting a few different physically sound assumptions. Comparing the predicted optical/UV continua and emission line ratios and profiles with the observed values we favor a particular model of a typical quasar. It radiates at ~0.5 of the Eddington luminosity, the corona surrounding the disk is additionally heated in excess of the Inverse Compton heating, and the broad line clouds are most probably destroyed when accelerated vertically above the sound speed within the corona.

The analysis of eight EXOSAT X-ray lightcurves of six active galactic nuclei by nonlinear prediction methods indicates that the observed short time-scale variability is truly stochastic and is not caused by deterministic chaos. This result favors X-ray emission models with mutliple centers of production of hot, possibly relativistic electrons.

We discuss the optical/UV/Xray spectra of AGN within the frame of the corona model of Witt, Czerny and Zycki (1996). In this model both the disk and the corona accrete and release energy through the viscosity. The relative strength the the disk and the corona is therefore determined by the model. Translated into the spectra, this model allows to predict the spectral index alpha_{ox}, measuring the relative strenth of the big blue bump with respect to the hard X-ray power law. Comparison of the predicted and observed distributions of alpha_{ox} allow us to conclude that the luminosity to the Eddington luminosity ratio in quasars cover typically the range 0.01 - 0.1 and it is broader in Seyferts, 0.001 - 0.3, or even higher. We identify Narrow Line Seyfert galaxies as objects having high luminosity to the Eddington luminosity ratio although, in principle, the second branch of high big blue bump objects for this ratio being below 0.001 is also predicted.

We analyse the data from the optical/IUE observational campaign of the Seyfert galaxy NGC 5548 in the context of 10 phenomenological models. On the basis of the optical/UV data as well as constraints from the X-ray observations we can favour one model of the nucleus: an accretion disc with an inner radius cutoff surrounded by a hot corona. The second acceptable model for optical/UV data is a distribution of optically thin clouds. However, X-ray constraints which were crucial in the analysis of disc type models could not have been applied; further development of this model is necessary.

We study an accretion disc with a hot continuous corona. We assume that the corona itself accretes and therefore it is powered directly by the release of the gravitational energy and cooled by radiative interaction with the disc. We consider the vertical structure of such a corona and show that the radial infall is accompanied by strong vertical outflow. Such a model has two important consequences: (i) at a given radius the corona forms only for accretion rate larger than a limiting value and the fraction of energy dissipated in the corona decreases with increasing accretion rate, and (ii) the disc spectra are significantly softer in the optical/UV band in comparison with the predictions of standard accretion discs due to the mass loss and the decrease of internal dissipation in the disc. Both trends correspond well to the mean spectra of radio quiet AGN and observed luminosity states in galactic black hole candidates.

We discuss a model for the central region of radio-quiet AGN in which coronal accretion is the source of energy for X-ray emission. We start from the Haardt & Maraschi model but supplement it with pressure balance condition, two possibilities of the angular momentum transfer and three models of the disc-corona transition. We employ the alpha-description of viscosity within the two layers. Both descriptions of the angular momentum transport lead to similar conclusions about the fraction of energy liberated in the disc although they predict strongly different properties of the disc interiors. The coupled disc/corona system has one family of solutions at low accretion rates, two families of solutions at intermediate accretion rates and no solution for high accretion rate. Of the two families of solutions, the first one has weakly dissipating coronae and the second one has coronae which dominate energetically the system. The two solutions merge at sub-Eddington accretion rate because a corona in hydrostatic equilibrium can dissipate only a limited amount of energy, which depends on the viscosity parameter.