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.

Recently, the LIGO-Virgo collaboration discovered gravitational waves and in their first publication on the subject the authors also presented a graviton mass constraint as $m_g < 1.2 \times 10^{-22}$ eV (Abbott et al., 2016). In the paper we analyze a potential to reduce upper bounds for graviton mass with future observational data on trajectories of bright stars near the Galactic Center. Since gravitational potentials are different for these two cases, expressions for relativistic advance for general relativity and Yukawa potential are different functions on eccentricity and semimajor axis, it gives an opportunity to improve current estimates of graviton mass with future observational facilities. In our considerations of an improvement potential for a graviton mass estimate we adopt a conservative strategy and assume that trajectories of bright stars and their apocenter advance will be described with general relativity expressions and it gives opportunities to improve graviton mass constraints. In contrast with our previous studies, where we present current constraints on parameters of Yukawa gravity (Borka et al., 2013) and graviton mass (Zakharov et al., 2016) from observations of S2 star, in the paper we express expectations to improve current constraints for graviton mass, assuming the GR predictions about apocenter shifts will be confirmed with future observations. We concluded that if future observations of bright star orbits during around fifty years will confirm GR predictions about apocenter shifts of bright star orbits it give an opportunity to constrain a graviton mass at a level around $5 \times 10^{-23}$ eV or slightly better than current estimates obtained with LIGO observations.

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.

Recently LIGO collaboration discovered gravitational waves \cite{Abbott_16} predicted 100 years ago by A. Einstein. Moreover, in the key paper reporting about the discovery, the joint LIGO \& VIRGO team presented an upper limit on graviton mass such as $m_g < 1.2 \times 10^{-22} eV$ (Abbott et al. (LIGO collaboration) PRL 116 (2016) 061102). Since the graviton mass limit is so small the authors concluded that their observational data do not show violations of classical general relativity. We consider another opportunity to evaluate a graviton mass from phenomenological consequences of massive gravity and show that an analysis of bright star trajectories could bound graviton mass with a comparable accuracy with accuracies reached with gravitational wave interferometers and expected with forthcoming pulsar timing observations for gravitational wave detection. It gives an opportunity to treat observations of bright stars near the Galactic Center as a wonderful tool not only for an evaluation specific parameters of the black hole but also to obtain constraints on the fundamental gravity law such as a modifications of Newton gravity law in a weak field approximation. In particular, we obtain bounds on a graviton mass based on a potential reconstruction at the Galactic Center.

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 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).

We use observations of the continuum radio emission at 1420, 820, 408, 34.5 and 22 MHz to estimate the mean brightness temperatures of the HB 21 supernova remnant (SNR) at five frequencies. We also presented mean spectral index of HB 21. The spectra of HB 21 are estimated for mean temperatures versus frequency for 1420, 820, 408, 34.5 and 22 MHz. We also presented $T-T$ plots of three frequency pairs: between 1420--34.5, 1420--22, 34.5--22 MHz. We noticed flatter spectral indices at frequencies below 408 MHz. Probably this is due to the absorption by thermal plasma at low frequencies.