We show that fundamental plane of elliptical galaxies can be used to obtain observational constraints on metric theories of gravity. Being it connected to global properties of ellipticals, it can fix parameters of modified gravity. Specifically, we use fundamental plane to constrain modified theories of gravity with Yukawa-like corrections which commonly emerge in the post-Newtonian limit. After giving examples on how these corrections are derived, we first analyze the velocity distribution of elliptical galaxies comparing theoretical results of modified gravity with Yukawa-like corrections with astronomical data. According to these results, it is possible to constrain the parameters of the corrections discriminating among classes of models compatible with astronomical observations. We conclude that fundamental plane can be used as a standard tool to probe different theories of gravity in the weak field limit.

Here we study the potential observational signatures of supermassive black hole binaries (SMBHBs) in the Fe K$\alpha$ line profiles emitted from the accretion disks around their components. We simulated the Fe K$\alpha$ line emission from the relativistic accretion disks using ray tracing method in Kerr metric. The obtained profiles from the SMBHBs are then compared with those in the case of the single supermassive black holes (SMBHs). We considered two models of the SMBHBs: a model when the secondary SMBH is embedded in the accretion disk around the primary, causing an empty gap in the disk, and a model with clearly separated components, where the accretion disks around both primary and secondary give a significant contribution to the composite Fe K$\alpha$ line emission of a such SMBHB. The obtained results showed that both models of SMBHBs can leave imprints in the form of ripples in the cores of the emitted Fe K$\alpha$ line profiles, which may look like an absorption component in the line profile. However, in the case of the composite line profiles emitted from two accretion disks, these ripples could have much higher amplitudes and strongly depend on orbital phase of the system, while for those emitted from a disk with an empty gap, the corresponding ripples mostly have lower amplitudes and do not vary significantly with orbital phase. The present day X-ray telescopes are not able to detect such signatures in the observed X-ray spectra of SMBHBs. However this will be possible with the next generation of X-ray observatories, which will also enable application of such effects as a tool for studying the properties of these objects.

The aim of our investigation is to derive a particular theory among the class of scalar-tensor(ST) theories of gravity, and then to test it by studying kinematics and dynamics of S-stars around supermassive black hole (BH) at Galactic Center (GC). We also discuss the Newtonian limit of this class of ST theories of gravity, as well as its parameters. We compare the observed orbit of S2 star with our simulated orbit which we obtained theoretically with the derived ST potential and constrained the parameters. Using the obtained best fit parameters we calculated orbital precession of S2 star in ST gravity,and found that it has the same direction as in General Relativity (GR), but causes much larger pericenter shift.

Non-local theories of gravity have recently gained a lot of interest because they can suitably represent the behavior of gravitational interaction in the ultraviolet regime. Furthermore, at infrared scales, they give rise to notable cosmological effects which could be important to describe the dark energy behavior. In particular, exponential forms of the distortion function seem particularly useful for this purpose. Using Noether Symmetries, it can be shown that the only non-trivial form of the distortion function is the exponential one, which is working not only for cosmological mini-superspaces, but also in a spherically symmetric spacetime. Taking this result into account, we study the weak field approximation of this type of non-local gravity, and comparing with the orbits of S2 stars around the Galactic center (NTT/VLT data), we set constraints on the parameters of the theory. Non-local effects do not play a significant role on the orbits of S2 stars around Sgr A*, but give richer phenomenology at cosmological scales than the $\Lambda$CDM model. Also, we show that non-local gravity model gives better agreement between theory and astronomical observations than Keplerian orbits.

The global properties of elliptical galaxies are connected through the so-called fundamental plane of ellipticals, which is an empirical relation between their parameters: effective radius, central velocity dispersion and mean surface brightness within the effective radius. We investigated the relation between the parameters of the fundamental plane equation and the parameters of modified gravity potential $f(R)$. With that aim, we compared theoretical predictions for circular velocity in $f(R)$ gravity with the corresponding values from a large sample of observed elliptical galaxies. Besides, we consistently reproduced the values of coefficients of the fundamental plane equation as deduced from observations, showing that the photometric quantities like mean surface brightness are related to gravitational parameters. We show that this type of modified gravity, especially its power-law version - $R^n$, is able to reproduce the stellar dynamics in elliptical galaxies. Also, it is shown that $R^n$ gravity fits the observations very well, without need for a dark matter.

The relativistically broadened Fe K$\alpha$ line, originating from the accretion disc in a vicinity of a super massive black hole, is observed in only less than 50\% of type 1 Active Galactic Nuclei (AGN). In this study we investigate could this lack of detections be explained by the effects of certain parameters of the accretion disc and black hole, such as the inclination, the inner and outer radius of disc and emissivity index. In order to determine how these parameters affect the Fe K $\alpha$ line shape, we simulated about 60,000 Fe K $\alpha$ line profiles emitted from the relativistic disc. Based on simulated line profiles, we conclude that the lack of the Fe K$\alpha$ line detection in type 1 AGN could, be caused by the specific emitting disc parameters, but also by the limits in the spectral resolution and sensitivity of the X-ray detectors.

We demonstrate that the existence of a Noether symmetry in $f(R)$ theories of gravity gives rise to a further gravitational radius, besides the standard Schwarzschild one, determining the dynamics at galactic scales. By this feature, it is possible to explain the baryonic Tully-Fisher relation and the rotation curve of gas-rich galaxies without the dark matter hypothesis.

The fundamental plane (FP) of galaxies can be recovered in the framework of $f(R)$ gravity avoiding the issues related to dark matter to fit the observations. In particular, the power-law version $f(R)\propto R^n$, resulting from the existence of Noether symmetries for $f(R)$, is sufficient to implement the approach. In fact, relations between the FP parameters and the corrected Newtonian potential, coming from $R^n$, can be found and justified from a physical point of view. Specifically, we analyze the velocity distribution of elliptical galaxies and obtain that $r_c$, the scale-length depending on the gravitational system properties, is proportional to $r_e$, the galaxy effective radius. This fact points out that the gravitational corrections induced by $f(R)$ can lead photometry and dynamics of the system. Furthermore, the main byproduct of such an approach is that gravity could work in different ways depending on the scales of self-gravitating systems.

We use radio-continuum all-sky surveys at 1420 and 408 MHz with the aim to investigate properties of the Galactic radio source Lupus Loop. The survey data at 1435 MHz, with the linear polarization of the southern sky, is also used. We calculate properties of this supernova remnant: the brightness temperature, surface brightness and radio spectral index. For determining borders and calculation of its properties, we use the method we have developed. The non-thermal nature of its radiation is confirmed. The distribution of spectral index over its area is also given. A significant correlation between the radio spectral index distribution and the corresponding polarized intensity distribution inside the loop borders is found, indicating that the polarization maps could provide us information about the distribution of interstellar medium, and thus could represent one additional way to search for new Galactic loops.

We present an analysis of 43 years (1972 to 2015) of spectroscopic observations of the Seyfert 1 galaxy NGC 5548. This includes 12 years of new unpublished observations (2003 to 2015). We compiled about 1600 H$\beta$ spectra and analyzed the long-term spectral variations of the 5100 \AA\ continuum and the H$\beta$ line. Our analysis is based on standard procedures including the Lomb-Scargle method, which is known to be rather limited to such heterogeneous data sets, and new method developed specifically for this project that is more robust and reveals a $\sim$5700 day periodicity in the continuum light curve, the H$\beta$ light curve and the radial velocity curve of the red wing of the H$\beta$ line. The data are consistent with orbital motion inside the broad emission line region of the source. We discuss several possible mechanisms that can explain this periodicity, including orbiting dusty and dust-free clouds, a binary black hole system, tidal disruption events, and the effect of an orbiting star periodically passing through an accretion disc.

The detection of gravitationally redshifted optical emission lines has been reported just for a few active galaxies. In this paper we give a short overview of studies that analyzed or exploited the detection of the gravitational redshift in optical AGN spectra. In addition, we tested the consistency of gravitational redshift as the physical origin of the redward shifts observed in their spectra using a sample of $\approx$ 50 Hamburg-ESO intermediate to high redshift quasars that are among the most luminous quasars known ($10^{47} \lesssim L \lesssim 10^{48}$ erg/s), and are expected to host very massive black holes. To this aim we modeled the line profile with accretion disk models around a black hole.

Here we present a short overview and main results of our investigations of several effects which can induce shifts in the broad Fe K$\alpha$ line emitted from relativistic accretion disks around single and binary supermassive black holes. We used numerical simulations based on ray-tracing method in the Kerr metric to study the role of classical Doppler shift, special relativistic transverse Doppler shift and Doppler beaming, general relativistic gravitational redshift, and perturbations of the disk emissivity in the formation of the observed Fe K$\alpha$ line profiles. Besides, we also investigated whether the observed line profiles from the binary systems of supermassive black holes could be affected by the Doppler shifts due to dynamics of such systems. The presented results demonstrate that all these effects could have a significant influence on the observed profiles of the broad Fe K$\alpha$ line emitted from relativistic accretion disks around single and binary supermassive black holes.

We investigate the possibility to explain theoretically the observed deviations of S2 star orbit around the Galactic Centre using gravitational potentials derived from modified gravity models in absence of dark matter. To this aim, an analytic fourth-order theory of gravity, non-minimally coupled with a massive scalar field is considered. Specifically, the interaction term is given by analytic functions $f(R)$ and $f(R,\phi)$ where $R$ is the Ricci scalar and $\phi$ is a scalar field whose meaning can be related to further gravitational degrees of freedom. We simulate the orbit of S2 star around the Galactic Centre in $f(R)$ (Yukawa-like) and $f(R,\phi)$ (Sanders-like) gravity potentials and compare it with NTT/VLT observations. Our simulations result in strong constraints on the range of gravity interaction. In the case of analytic functions $f(R)$, we are not able to obtain reliable constraints on the derivative constants $f_1$ and $f_2$, because the current observations of S2 star indicated that they may be highly mutually correlated. In the case of analytic functions $f(R,\phi)$, we are able to obtain reliable constraints on the derivative constants $f_0$, $f_R$, $f_{RR}$, $f_{\phi}$, $f_{\phi\phi}$ and $f_{\phi R}$. The approach we are proposing seems to be sufficiently reliable to constrain the modified gravity models from stellar orbits around Galactic Centre.

Here we investigate possible applications of observed stellar orbits around Galactic Center for constraining the R$^n$ gravity at Galactic scales. For that purpose, we simulated orbits of S2-like stars around the massive black hole at Galactic Center, and study the constraints on the R$^n$ gravity which could be obtained by the present and next generations of large telescopes. Our results show that R$^n$ gravity affects the simulated orbits in the qualitatively similar way as a bulk distribution of matter (including a stellar cluster and dark matter distributions) in Newton's gravity. In the cases where the density of extended mass is higher, the maximum allowed value of parameter $\beta$ in R$^n$ gravity is noticeably smaller, due to the fact that the both extended mass and $R^n$ gravity cause the retrograde orbital precession.

We consider possible signatures for Yukawa gravity within the Galactic Central Parsec, based on our analysis of the S2 star orbital precession around the massive compact dark object at the Galactic Centre, and on the comparisons between the simulated orbits in Yukawa gravity and two independent sets of observations. Our simulations resulted in strong constraints on the range of Yukawa interaction $\Lambda$ and showed that its most probable value in the case of S2 star is around 5000 - 7000 AU. At the same time, we were not able to obtain reliable constrains on the universal constant $\delta$ of Yukawa gravity, because the current observations of S2 star indicated that it may be highly correlated with parameter $\Lambda$ in the range $(0 <\delta < 1)$. For $\delta > 2$ they are not correlated. However, the same universal constant which was successfully applied to clusters of galaxies and rotation curves of spiral galaxies ($\delta=1/3$) also gives a satisfactory agreement with the observed orbital precession of the S2 star, and in that case the most probable value for the scale parameter is $\Lambda \approx 3000 \pm 1500$ AU. Also, the Yukawa gravity potential induces precession of S2 star orbit in the same direction as General Relativity for $\delta > 0$ and for $\delta < -1$, and in the opposite direction for $-1 <\delta < 0$. The future observations with advanced facilities, such as GRAVITY or/and European Extremely Large Telescope, are needed in order to verify these claims.

We used a model of a relativistic accretion disk around a supermassive black hole (SMBH), based on ray-tracing method in the Kerr metric, to study the variations of the composite Fe K$\alpha$ line emitted from two accretion disks around SMBHs in a binary system. We assumed that the orbit of such a binary is approximately Keplerian, and simulated the composite line shapes for different orbital elements, accretion disk parameters and mass ratios of the components. The obtained results show that, if observed in the spectra of some SMBH binaries during their different orbital phases, such composite Fe K$\alpha$ line profiles could be used to constrain the orbits and several properties of such SMBH binaries.

One of the most intriguing scenarios proposed to explain how active galactic nuclei are triggered involves the existence of a supermassive binary black hole system in their cores. Here we present an observational evidence for the first spectroscopically resolved sub-parsec orbit of a such system in the core of Seyfert galaxy NGC 4151. Using a method similar to those typically applied for spectroscopic binary stars we obtained radial velocity curves of the supermassive binary system, from which we calculated orbital elements and made estimates about the masses of components. Our analysis shows that periodic variations in the light and radial velocity curves can be accounted for an eccentric, sub-parsec Keplerian orbit of a 15.9-year period. The flux maximum in the lightcurve correspond to the approaching phase of a secondary component towards the observer. According to the obtained results we speculate that the periodic variations in the observed H{\alpha} line shape and flux are due to shock waves generated by the supersonic motion of the components through the surrounding medium. Given the large observational effort needed to reveal this spectroscopically resolved binary orbital motion we suggest that many such systems may exist in similar objects even if they are hard to find. Detecting more of them will provide us with insight into black hole mass growth process.

We investigated the gravitational microlensing of active galactic nucleus dusty tori in the case of lensed quasars in the infrared domain. The dusty torus is modeled as a clumpy two-phase medium. To obtain spectral energy distributions and images of tori at different wavelengths, we used the 3D Monte Carlo radiative transfer code SKIRT. A ray-shooting technique has been used to calculate microlensing magnification maps. We simulated microlensing by the stars in the lens galaxy for different configurations of the lensed system and different values of the torus parameters, in order to estimate (a) amplitudes and timescales of high magnification events, and (b) the influence of geometrical and physical properties of dusty tori on light curves in the infrared domain. We found that, despite their large size, dusty tori could be significantly affected by microlensing in some cases, especially in the near-infrared domain (rest-frame). The very long time-scales of such events, in the range from several decades to hundreds of years, are limiting the practical use of this method to study the properties of dusty tori. However, our results indicate that, when studying flux ratios between the images in different wavebands of lensed quasars, one should not disregard the possibility that the near and mid-infrared flux ratios could be under the influence of microlensing.

We study some possible observational signatures of $R^n$ gravity at Galactic scales and how these signatures could be used for constraining this type of $f(R)$ gravity. For that purpose, we performed two-body simulations in $R^n$ gravity potential and analyzed the obtained trajectories of S2-like stars around Galactic center, as well as resulting parameter space of $R^n$ gravity potential. Here, we discuss the constraints on the $R^n$ gravity which can be obtained from the observations of orbits of S2-like stars with the present and next generations of large telescopes. We make comparison between the theoretical results and observations. Our results show that the most probable value for the parameter $r_c$ in $R^n$ gravity potential in the case of S2-like stars is $\sim$100 AU, while the universal parameter $\beta$ is close to 0.01. Also, the $R^n$ gravity potential induces the precession of S2-like stars orbit in opposite direction with respect to General Relativity, therefore, such a behavior of orbits qualitatively is similar to a behavior of Newtonian orbits with a bulk distribution of matter (including a stellar cluster and dark matter distributions).

Here we present an overview of some of the most significant observational and theoretical studies of the broad Fe K{\alpha} spectral line, which is believed to originate from the innermost regions of relativistic accretion disks around central supermassive black holes of galaxies. The most important results of our investigations in this field are also listed. All these investigations indicate that the broad Fe K{\alpha} line is a powerful tool for studying the properties of the supermassive black holes (such as their masses and spins), space-time geometry (metric) in their vicinity, their accretion physics, probing the effects of their strong gravitational fields, and for testing the certain predictions of General Relativity.

We report the new spectroscopic observations of the gravitational lens RXJ 021+4529 with the multi-mode focal reducer SCORPIO of the SAO RAS 6-m telescope. The new spectral observations were compared with the previously observed spectra of components A and B of RXJ 0921+4529, i.e. the same components observed in different epochs. We found a significant difference in the spectrum between the components that cannot be explained with microlensing and/or spectral variation. We conclude that RXJ 0921+4529 is a binary quasar system, where redshifts of quasars A and B are 1.6535 +/- 0.0005 and 1.6625 +/- 0.0015, respectively.

Here we show that in the case when double peaked emission lines originate from outer parts of accretion disk, their variability could be caused by perturbations in the disk emissivity. In order to test this hypothesis, we introduced a model of disk perturbing region in the form of a single bright spot (or flare) by a modification of the power law disk emissivity in appropriate way. The disk emission was then analyzed using numerical simulations based on ray-tracing method in Kerr metric and the corresponding simulated line profiles were obtained. We applied this model to the observed H-beta line profiles of 3C 390.3 (observed in the period 1995-1999), and estimated the parameters of both, accretion disk and perturbing region. Our results show that two large amplitude outbursts of the H-beta line observed in 3C 390.3 could be explained by successive occurrences of two bright spots on approaching side of the disk. These bright spots are either moving, originating in the inner regions of the disk and spiralling outwards by crossing small distances during the period of several years, or stationary. In both cases, their widths increase with time, indicating that they most likely decay.

Here we give a brief overview of some investigations of the gravitational microlensing influence on broad absorption spectral lines of lensed QSOs. Especially, we consider the microlensing influence on the Fe K$\alpha$ broad absorption lines using a model of an accretion disk covered by an absorption region. Gravitational microlensing is modeled by ray shooting method which enables us to obtain realistic microlensing patterns. We obtain that microlensing can affect both emission and absorption component of line that depends on dimensions on emission and absorption line regions. Here we give detailed analysis of emission and absorption line shape variations due to gravitational microlensing.

In this chapter we discuss the X-ray radiation from relativistic accretion disks around supermassive black holes, supposed to exist in the centers of Active Galactic Nuclei (AGN). Our focus is on the X-ray radiation, especially in the Fe K$\alpha$ line which originates in the innermost parts of an accretion disk. Moreover, here we discuss some effects which can disturb the Fe K$\alpha$ profile and cause its rapid and irregular variability, observed in the X-ray spectra of some AGN. We will pay attention to three such effects: perturbations in the disk emissivity, absorbtion by warm absorbers and gravitational microlensing. The X-ray emission from accretion disks around non-rotating (Schwarzschild metric), as well as rotating (Kerr metric) supermassive black holes, is discussed. The X-ray radiation of AGN is probably produced in a compact region near their central supermassive black holes, and can provide us some essential information about the plasma conditions and the space-time geometry in these regions. The goal of this chapter is mainly to present a short overview of some important and recent investigations in this field.

Here we discuss the effects of strong gravity that can be observed in electromagnetic spectra of active galactic nuclei (AGN). According to the unification model of an AGN, there is a supermassive black hole ($10^7 - 10^9 M_\odot$) in its center, surrounded by an accretion disk that radiates in the X-ray band. Accretion disks could have different forms, dimensions, and emission, depending on the type of central black hole (BH), whether it is rotating (Kerr metric) or nonrotating (Schwarzschild metric). We modeled the emission of an accretion disk around supermassive BH using numerical simulations based on a ray-tracing method in the Kerr metric. A broad emission line Fe K$\alpha$ at 6.4 keV with asymmetric profile (narrow bright blue peak and a wide faint red wing) has been observed in a number of type 1 AGN. The effects of strong gravitational field are investigated by comparison between the modeled and observed iron K$\alpha$ line profiles. The results of our modeling show that the parameters of the Fe K$\alpha$ line emitting region have significant influence on the line profile and thus, allow us to determine the space-time geometry (metric) in vicinity of the central BH of AGN, and also can give us information about the plasma conditions in these regions.

We consider the influence of microlensing on different spectral bands of lensed QSOs. We assumed that the emitting X-ray, UV and optical regions are different in size, but that the continuum emission in these spectral bands is originating from an accretion disc. Estimations of the time scales for microlensing and flux amplification in different bands are given. We found that the microlensing duration should be shorter in the X-ray (several months) than in UV/optical emitting region (several years). This result indicates that monitoring of the X-ray variations in lensed QSOs that show a 'flux anomaly' can clarify the source of this anomaly.