An active stage of relativistic astrophysics started in 1963 since in this year, quasars were discovered, Kerr solution has been found and the first Texas Symposium on Relativistic Astrophysics was organized in Dallas. Five years later, in 1967--1968 pulsars were discovered and their model as rotating neutron stars has been proposed, meanwhile J. A. Wheeler claimed that Kerr and Schwarzschild vacuum solutions of Einstein equations provide an efficient approach for astronomical objects with different masses. Wheeler suggested to call these objects black holes. Neutron stars were observed in different spectral band of electromagnetic radiation. In addition, a neutrino signal has been found for SN1987A. Therefore, multi-messenger astronomy demonstrated its efficiency for decades even before observations of the first gravitational radiation sources. However, usually, one has only manifestations of black holes in a weak gravitational field limit and sometimes a model with a black hole could be substituted with an alternative approach which very often looks much less natural, however, it is necessary to find observational evidences to reject such an alternative model. After two observational runs the LIGO-- Virgo collaboration provided a confirmation for an presence of mergers for ten binary black holes and one binary neutron star system where gravitational wave signals were found. In addition, in last years a remarkable progress has been reached in a development of observational facilities to investigate a gravitational potential, for instance, a number of telescopes operating in the Event Horizon Telescope network is increasing and accuracy of a shadow reconstruction near the Galactic Center is improving, meanwhile largest VLT, Keck telescopes with adaptive optics and especially, GRAVITY facilities observe bright IR stars at the Galactic Center with a perfecting accuracy.

As it was pointed out recently in Hees et al. (2017), observations of stars near the Galactic Center with current and future facilities provide an unique tool to test general relativity (GR) and alternative theories of gravity in a strong gravitational field regime. In particular, the authors showed that the Yukawa gravity could be constrained with Keck and TMT observations. Some time ago, Dadhich et al. (2001) showed that the Reissner -- Nordstr\"om metric with a tidal charge is naturally appeared in the framework of Randall -- Sundrum model with an extra dimension ($Q^2$ is called tidal charge and it could be negative in such an approach). Astrophysical consequences of of presence of black holes with a tidal charge are considerered, in particular, geodesics and shadows in Kerr -- Newman braneworld metric are analyzed in (Schee and Stuchlik, 2009a), while profiles of emission lines generated by rings orbiting braneworld Kerr black hole are considered in (Schee and Stuchlik, 2009b). Possible observational signatures of gravitational lensing in a presence of the Reissner -- Nordstr\"om black hole with a tidal charge at the Galactic Center are discussed in papers by Bin-Nun (2010a, 2010b, 2011). Here we are following such an approach and we obtain analytical expressions for orbital precession for Reissner -- Nordstrom -- de-Sitter solution in post-Newtonian approximation and discuss opportunities to constrain parameters of the metric from observations of bright stars with current and future astrometric observational facilities such as VLT, Keck, GRAVITY, E-ELT and TMT.

To evaluate a potential usually one analyzes trajectories of test particles. For the Galactic Center case astronomers use bright stars or photons, so there are two basic observational techniques to investigate a gravitational potential, namely, (a) monitoring the orbits of bright stars near the Galactic Center as it is going on with 10m Keck twin and four 8m VLT telescopes equipped with adaptive optics facilities (in addition, recently the IR interferometer GRAVITY started to operate with VLT); (b) measuring the size and shape of shadows around black hole with VLBI-technique using telescopes operating in mm-band. At the moment, one can use a small relativistic correction approach for stellar orbit analysis, however, in the future the approximation will not be precise enough due to enormous progress of observational facilities and recently the GRAVITY team found that the first post-Newtonian correction has to be taken into account for the gravitational redshift in the S2 star orbit case. Meanwhile for smallest structure analysis in VLBI observations one really needs a strong gravitational field approximation. We discuss results of observations and their interpretations.

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.

One of the most interesting astronomical objects is the Galactic Center. We concentrate our discussion on a theoretical analysis of observational data of bright stars in the IR-band obtained with large telescopes. We also discuss the importance of VLBI observations of bright structures which could characterize the shadow at the Galactic Center. There are attempts to describe the Galactic Center with alternative theories of gravity and in this case one can constrain parameters of such theories with observational data for the Galactic Center. In particular, theories of massive gravity are intensively developing and theorists have overcome pathologies presented in initial versions of these theories. In theories of massive gravity, a graviton is massive in contrast with GR where a graviton is massless. Now these theories are considered as an alternative to GR. For example, the LIGO-Virgo collaboration obtained the graviton mass constraint of about $1.2 \times 10^{-22}$ eV in their first publication about the discovery of the first gravitational wave detection event that resulted of the merger of two massive black holes. Surprisingly, one could obtain a consistent and comparable constraint of graviton mass at a level around $m_{g} < 2.9 \times 10^{-21}$ eV from analysis of observational data on the trajectory of the star S2 near the Galactic Center. Therefore, observations of bright stars with existing and forthcoming telescopes such as the European Extremely Large Telescope (E-ELT) and the Thirty Meter Telescope (TMT) are extremely useful for investigating the structure of the Galactic Center in the framework of GR, but these observations also give a tool to confirm, rule out or constrain alternative theories of gravity.

The most likely cause of levitation of dust above the surface of atmosphereless planets is the electrostatic mechanism. However, the crucial problem in the explanation of this effect is a determination of the reason why a large electric charge (units or even dozens of elementary charges) required for take-off can be accumulated on the smallest dust particles. Due to the photoeffect the charge of such value could be easily accumulated on a solitary dust particle, but if a dust particle has not yet taken off, the average value of its charge is several orders of magnitude lower because of the extremely low charge density on the planet's surface. The paper shows that surface charge density is really small only for averaging over regions of macroscopic size, and on a submicron scale the surface appear to be a collection of chaotic "spots" with charges of different signs. The positively charged "spots" are only slightly larger than the negatively charged spots, which explains the small value of the charge density averaged over macroscopic regions. "Spots" arise due to fluctuations in the fluxes of the photoelectrons taking off and falling back the surface, and the charge density inside the "spots" is sufficient to allow a takeoff of particles with the dimensions less than or on the order of 0.1 micron in the field of a double layer.

Neglecting the effects associated with the solar wind plasma, the photoelectrons are the only elementary particles which create an electrical current through sunlit surface of the moon. They are knocked off of the surface soil, rise above the surface, and then fall back. Therefore, on average, on any unit of surface area there is a positive charge, equal in magnitude to the charge of photoelectrons flying over this area. However, the charge of any small dust particle can strongly fluctuate discretely: a photoelectron can be either knocked off of the or be reacquired by the particle. The result is a "random walk" in sign and magnitude of the charge of grains. In a few minutes after sunrise, almost every dust particle on the surface has at least one extra or missing electron, and the average modulus of the charge accumulated on a particle is proportional to the square root of the number of "steps" (knocking off /returns of photoelectrons). Therefore, the average value of the modulus of the charge of a fine dust particle exceeds by several orders of magnitude the proportion of the average surface charge attributable to the particle. So dust particles that have ejected a sufficient number of photoelectrons can take off from the surface because of the electric field of the near surface charge double layer. It is shown that: (i) almost half of all dust particles on the illuminated lunar surface are missing at least one electron; (ii) a significant portion of particles up to 100 nm in size emits several photoelectrons and acquire a positive charge, sufficient to take off from the surface; (iii) above the surface there is a "boiling layer" of dust with a maximum thickness of several hundred meters where the average size of the particles and their density non-monotonically depend on the altitude.

The principal goal of this paper is to specify conditions of detection of unresolved binaries by multicolor photometry. We have developed a method for estimating the critical distance at which an unresolved binary of given mass and age can be detected. The method is applied to the photometric system of the planned Lyra-B spaceborne experiment. We have shown that some types of unresolved binary stars can be discovered and distinguished from single stars solely by means of photometric observations.

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.

We propose an improved method for the atmospheric extinction reduction within optical photometry. Our method is based on the simultaneous multicolor observations of photometric standards. Such data are now available within the modern wide-field sky surveys and contain a large amount of information about instant atmospheric conditions. So, it became possible to estimate the extinction parameters on the basis of a quite short observational dataset and, hence, to trace the rapid stars twinkling accurately. Having been developed for a new MiniMegaTORTORA observational system, the proposed method can be adopted for a wide range of modern observational programs.

Using an algebraic condition of vanishing discriminant for multiple roots of fourth degree polynomials we derive an analytical expression of a shadow size as a function of a charge in the Reissner -- Nordstr\"om (RN) metric \cite{Reissner_16,Nordstrom_18}. We consider shadows for negative tidal charges and charges corresponding to naked singularities $q=\mathcal{Q}^2/M^2 > 1$, where $\mathcal{Q}$ and $M$ are black hole charge and mass, respectively, with the derived expression. An introduction of a negative tidal charge $q$ can describe black hole solutions in theories with extra dimensions, so following the approach we consider an opportunity to extend RN metric to negative $\mathcal{Q}^2$, while for the standard RN metric $\mathcal{Q}^2$ is always non-negative. We found that for $q > 9/8$ black hole shadows disappear. Significant tidal charges $q=-6.4$ (suggested by Bin-Nun (2010)) are not consistent with observations of a minimal spot size at the Galactic Center observed in mm-band, moreover, these observations demonstrate that a Reissner -- Nordstr\"om black hole with a significant charge $q \approx 1$ provides a better fit of recent observational data for the black hole at the Galactic Center in comparison with the Schwarzschild black hole.

Using Schwarzschild -- de-Sitter (Kottler) metric we derive a simple analytical relation between a shadow size and $\Lambda$-term. Current observations of the smallest spot to evaluate shadow size at the Galactic Center do not reach an accuracy comparable with cosmological $\Lambda$-term $\sim 10^{-52}{\rm m}^{-2}$, however, if in reality we have dark energy instead of a constant $\Lambda$-term then dark energy may be a function depending on time and space and it could be approximated with a local constant $\Lambda$-term near the Galactic Center and it is important to introduce a procedure to evaluate the $\Lambda$-term. We suggest such a procedure based on a black hole shadow evaluation. Surprisingly, current observational estimates of shadows are in agreement with anti-de-Sitter spacetimes corresponding to a negative $\Lambda$-term which is about $-0.4\times 10^{-20}{\rm m}^{-2}$. A negative $\Lambda$-term has been predicted in the framework of a some class of multidimensional string models.

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.

It is shown that the appearance of lunar horizon glow and streamers observed above the lunar terminator described by the dynamic fountain mode byl Stubbs et al. (2006) requires a value of dust particles charge several orders exceed what they obtain on the Lunar surface. To obtain a sufficient charge due to a departure of photoelectrons separated submicron particles have to flow over the surface during tens of seconds or even several minutes. Therefore for emergence of dust streamers at the lunar sunrise it is necessary that dust particles do not lose the positive charge during a night.

There are different methods for finding exoplanets such as radial spectral shifts, astrometrical measurements, transits, timing etc. Gravitational microlensing (including pixel-lensing) is among the most promising techniques with the potentiality of detecting Earth-like planets at distances about a few astronomical units from their host star or near the so-called snow line with a temperature in the range $0-100^0$ C on a solid surface of an exoplanet. We emphasize the importance of polarization measurements which can help to resolve degeneracies in theoretical models. In particular, the polarization angle could give additional information about the relative position of the lens with respect to the source.

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.

The light received by source stars in microlensing events may be significantly polarized if both an efficient photon scattering mechanism is active in the source stellar atmosphere and a differential magnification is therein induced by the lensing system. The best candidate events for observing polarization are highly magnified events with source stars belonging to the class of cool, giant stars {in which the stellar light is polarized by photon scattering on dust grains contained in their envelopes. The presence in the stellar atmosphere of an internal cavity devoid of dust produces polarization profiles with a two peaks structure. Hence, the time interval between them gives an important observable quantity directly related to the size of the internal cavity and to the model parameters of the lens system.} We show that {during a microlensing event} the expected polarization variability can solve an ambiguity, that arises in some cases, related to the binary or planetary lensing interpretation of the perturbations observed near the maximum of the event light-curve. We consider a specific event case for which the parameter values corresponding to the two solutions are given. Then, assuming a polarization model for the source star, we compute the two expected polarization profiles. The position of the two peaks appearing in the polarization curves and the characteristic time interval between them allow us to distinguish between the binary and planetary lens solutions.

Gravitational microlensing, in case of relevant finite source size effects, provides an unique tool for the study of stellar atmospheres through the enhancement of a characteristic polarization signal. Here, we consider a set of highly magnified events and show that for different types of source stars (as hot, late type main sequence and cool giants) showing that the polarization strength may be of $\simeq 0.04$ percent for late type stars and up to a few percent for cool giants.

Gravitational microlensing, when finite size source effects are relevant, provides an unique tool for the study of source star stellar atmospheres through an enhancement of a characteristic polarization signal. This is due to the differential magnification induced during the crossing of the source star. In this paper we consider a specific set of reported highly magnified, both single and binary exoplanetary systems, microlensing events towards the Galactic bulge and evaluate the expected polarization signal. To this purpose, we consider several polarization models which apply to different types of source stars: hot, late type main sequence and cool giants. As a result we compute the polarization signal P,which goes up to P=0.04% for late type stars and up to a few percent for cool giants, depending on the underlying physical polarization processes and atmosphere model parameters. Given a I band magnitude at maximum magnification of about 12, and a typical duration of the polarization signal up to 1 day, we conclude that the currently available technology, in particular the polarimeter in FORS2 on the VLT, potentially may allow the detection of such signals. This observational programme may take advantage of the currently available surveys plus follow up strategy already routinely used for microlensing monitoring towards the Galactic bulge (aimed at the detection of exoplanets). In particular, this allows one to predict in advance for which events and at which exact time the observing resources may be focused to make intensive polarization measurements.

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

This paper on ASTROD I is based on our 2010 proposal submitted for the ESA call for class-M mission proposals, and is a sequel and an update to our previous paper [Experimental Astronomy 23 (2009) 491-527; designated as Paper I] which was based on our last proposal submitted for the 2007 ESA call. In this paper, we present our orbit selection with one Venus swing-by together with orbit simulation. In Paper I, our orbit choice is with two Venus swing-bys. The present choice takes shorter time (about 250 days) to reach the opposite side of the Sun. We also present a preliminary design of the optical bench, and elaborate on the solar physics goals with the radiation monitor payload. We discuss telescope size, trade-offs of drag-free sensitivities, thermal issues and present an outlook. ASTROD I is a planned interplanetary space mission with multiple goals. The primary aims are: to test General Relativity with an improvement in sensitivity of over 3 orders of magnitude, improving our understanding of gravity and aiding the development of a new quantum gravity theory; to measure key solar system parameters with increased accuracy, advancing solar physics and our knowledge of the solar system; and to measure the time rate of change of the gravitational constant with an order of magnitude improvement and the anomalous Pioneer acceleration, thereby probing dark matter and dark energy gravitationally. It is envisaged as the first in a series of ASTROD missions. ASTROD I will consist of one spacecraft carrying a telescope, four lasers, two event timers and a clock. Two-way, two-wavelength laser pulse ranging will be used between the spacecraft in a solar orbit and deep space laser stations on Earth, to achieve the ASTROD I goals.

The measurement of relativistic effects around the galactic center may allow in the near future to strongly constrain the parameters of the supermassive black hole likely present at the galactic center (Sgr A*). As a by-product of these measurements it would be possible to severely constrain, in addition, also the parameters of the mass-density distributions of both the innermost star cluster and the dark matter clump around the galactic center.

Assuming that supernovae type Ia (SNe Ia) are standard candles one could use them to test cosmological theories. The Hubble Space Telescope team analyzed 186 SNe Ia\cite{Riess_04} to test the Standard Cosmological model (SC) associated with expanded lengths in the Universe and evaluate its parameters. We use the same sample to determine parameters of Conformal Cosmological model (CC) with relative reference units of intervals, so that conformal quantities of General Relativity are interpreted as observables. We concluded, that really the test is extremely useful and allows to evaluate parameters of the model. From a formal statistical point of view the best fit of the CC model is almost the same quality approximation as the best fit of SC model with $\Omega_\Lambda=0.72, \Omega_m=0.28$. As it was noted earlier, for CC models, a rigid matter component could substitute the $\Lambda$-term (or quintessence) existing in the SC model. We note that a free massless scalar field can generate such a rigid matter. We describe results of our analysis for more recent "gold" data (for 192 SNe Ia).

The space photometric project "Lyra" is developed now in Russia.The project purpose is determination of the photometric information and coordinates of the natural and artificial space objects, from the brightest ones to 16^m in visual lights. It is supposed to obtain the data for about 40-400 million objects from board of International Space Station, using an astronomical telescope with diameter of the main mirror of 0.5 m. The observations will be carried out in a scanning mode. Photometry will be obtained in 10 spectral bands. The expected uncertainty of magnitudes for objects of 16^m in the V-band is 0.001^m. The main results of experiment should become: 1) creation of spatial model of the Galaxy on distances to 3 kpc from the Sun; 2) specification of physical parameters of stars and models of star evolution; 3) discovering of a huge number (to several millions) of variable stars and determination of their variability parameters. The central wavelengths of the 10 bands of "Lyra" photometric system will be at 195, 218, 270, 350, 440, 550, 700, 825, 930 and 1000 nm. It is shown that combinations of various colour indices will allow to determine confidently both effective temperature and metallicity of stars. The presence of a 218 nm band allows to determine confidently interstellar extinction on stars of O-F spectral classes. The photometric system will make it possible to separate halo stars from disk stars and to define physical parameters of their atmospheres.

Several exoplanets have been detected towards the Galactic bulge with the microlensing technique. We show that exoplanets in M31 may also be detected with the pixel-lensing method, if telescopes making high cadence observations of an ongoing microlensing event are used. Using a Monte Carlo approach we find that the mean mass for detectable planetary systems is about $2 M_{\rm {J}}$. However, even small mass exoplanets ($M_{\rm P} < 20 M_{\oplus}$) can cause significant deviations, which are observable with large telescopes. We reanalysed the POINT-AGAPE microlensing event PA-99-N2. First, we test the robustness of the binary lens conclusion for this light curve. Second, we show that for such long duration and bright microlensing events, the efficiency for finding planetary-like deviations is strongly enhanced with respect to that evaluated for all planetary detectable events.

We show that exoplanets in the M31 galaxy may be detected with the pixel-lensing method by using telescopes making high cadence observations of an ongoing microlensing event. Although the mean mass for detectable exoplanets is about $2 M_{\rm {J}}$, even small mass exoplanets ($M_{\rm P} < 20 M_{\oplus}$) can cause significant deviations, which are observable with large telescopes. We reanalysed the POINT-AGAPE microlensing event PA-99-N2. First, we test the robustness of the binary lens conclusion for this light curve. Second, we show that for such long duration and bright microlensing events, the efficiency for finding planetary-like deviations is strongly enhanced.

We study the possibility to detect extrasolar planets in M31 through pixel-lensing observations. Using a Monte Carlo approach, we select the physical parameters of the binary lens system, a star hosting a planet, and we calculate the pixel-lensing light curve taking into account the finite source effects. Indeed, their inclusion is crucial since the sources in M31 microlensing events are mainly giant stars. Light curves with detectable planetary features are selected by looking for significant deviations from the corresponding Paczy\'{n}ski shapes. We find that the time scale of planetary deviations in light curves increase (up to 3-4 days) as the source size increases. This means that only few exposures per day, depending also on the required accuracy, may be sufficient to reveal in the light curve a planetary companion. Although the mean planet mass for the selected events is about $2 M_{\rm {Jupiter}}$, even small mass planets ($M_{\rm P} < 20 M_{\oplus}$) can cause significant deviations, at least in the observations with large telescopes. However, even in the former case, the probability to find detectable planetary features in pixel-lensing light curves is at most a few percent of the detectable events, and therefore many events have to be collected in order to detect an extrasolar planet in M31. Our analysis also supports the claim that the anomaly found in the candidate event PA-99-N2 towards M31 can be explained by a companion object orbiting the lens star.

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.

The existence of dark matter (DM) at scales of few pc down to $\simeq 10^{-5}$ pc around the centers of galaxies and in particular in the Galactic Center region has been considered in the literature. Under the assumption that such a DM clump, principally constituted by non-baryonic matter (like WIMPs) does exist at the center of our galaxy, the study of the $\gamma$-ray emission from the Galactic Center region allows us to constrain both the mass and the size of this DM sphere. Further constraints on the DM distribution parameters may be derived by observations of bright infrared stars around the Galactic Center. Hall and Gondolo \cite{hallgondolo} used estimates of the enclosed mass obtained in various ways and tabulated by Ghez et al. \cite{Ghez_2003,Ghez_2005}. Moreover, if a DM cusp does exist around the Galactic Center it could modify the trajectories of stars moving around it in a sensible way depending on the DM mass distribution. Here, we discuss the constraints that can be obtained with the orbit analysis of stars (as S2 and S16) moving inside the DM concentration with present and next generations of large telescopes. In particular, consideration of the S2 star apoastron shift may allow improving limits on the DM mass and size.

It is well established that certain detached eclipsing binary stars exhibit apsidal motions whose value is in disagreement with with calculated deviations from Keplerian motion based on tidal effects and the general theory of relativity. Although many theoretical senarios have been demonstrated to bring calculations into line with observations, all have seemed unlikely for various reasons. In particular, it has been established that the hypothesis of a third star in an orbit almost perpendicular to the orbital plane of the close binary system can explain the anomalous motion in at least some cases. The stability of triple star systems with highly inclined orbits has been in doubt, however. We have found conditions which allow the long term stability of such systems so that the third body hypothesis now seems a likely resolution of the apsidal motion problem. We apply our stability criteria to the cases of AS Cam and DI Her and recommend observations at the new Keck interferometer which should be able to directly observe the third bodies in these systems.

Several stars orbit around a black hole candidate of mass $3.7\times 10^6$ M$_{\odot}$, in the region of the Galactic Center (GC). Looking for General Relativistic (GR) periastron shifts is limited by the existence of a stellar cluster around the black hole that would modify the orbits due to classical effects that might mask the GR effect. Only if one knows the cluster parameters (its mass and core radius) it is possible to unequivocally deduce the GR effects expected and then test them. In this paper it is shown that the observation of the proper motion of Sgr A$^*$, $v_{Sgr A^*} = (0.4\pm 0.9)$ km s$^{-1}$ (\citealt{reid2004}), could help us to constrain the cluster parameters significantly and that future measurements of the periastron shifts for at least three stars may adequately determine the cluster parameters and the mass of the black hole.

Basics of standard theory of microlensing are introduced. Results of microlensing observations toward Magellanic Clouds and relations with dark matter (DM) problem in our Galaxy are described. Pixel microlensing observations and recent discoveries of planets with microlensing observations are listed.

Assuming that supernovae type Ia (SNe Ia) are standard candles one could use them to test cosmological theories. The Hubble Space Telescope team analyzed 186 SNe Ia (Riess et al. (2004)) to test the standard cosmological model (SC) and evaluate its parameters. We use the same sample to determine parameters of Conformal Cosmological models (CC). We concluded, that really the test is extremely useful and allows to evaluate parameters of the model. From a formal statistical point of view the best fit of the CC model is almost the same quality approximation as the best fit of SC model with $\Omega_\Lambda=0.72, \Omega_m=0.28$. As it was noted earlier, for CC models, a rigid matter component could substitute the $\Lambda$-term (or quintessence) existing in the SC model.

Recently, gravitational microlensing has been investigated in the framework of the weak field limit of fourth order gravity theory. However, solar system data (i.e. planetary periods and light bending) can be used to put strong constraints on the parameters of this class of gravity theories. We find that these parameters must be very close to those corresponding to the Newtonian limit of the theory.

The Hamiltonian formulation of general relativity (GR) is considered in finite space-time and a specific frame of reference given by the diffeo-invariant components of the Fock simplex in terms of the Dirac -- ADM variables. The evolution parameter and energy invariant with respect to the time-coordinate transformations are constructed by the separation of the cosmological scale factor and its identification with the spatial averaging of the metric determinant, so that the dimension of the kinemetric group of diffeomorphisms coincides with the dimension of a set of variables whose velocities are removed by the Gauss-type constraints in accordance with the second N\"other theorem. This coincidence allows us to solve the energy constraint, fulfil Dirac's Hamiltonian reduction, and to describe the potential perturbations in terms of the Lichnerowicz scale-invariant variables distinguished by the absence of the time derivatives of the spatial metric determinant. It was shown that the Hamiltonian version of the cosmological perturbation theory acquires attributes of the theory of superfluid liquid, and it leads to the generalization of the Schwarzschild solution. The astrophysical application of this approach to GR is considered under supposition that the Dirac -- ADM Hamiltonian frame is identified with that of the Cosmic Microwave Background radiation distinguished by its dipole component in the frame of an Earth observer.

For microlensing case angular distances between images or typical astrometric shifts due to microlensing are about $10^{-5}-10^{-6} \mu as$. Such an angular resolution will be reached with the space space--ground interferometer Radioastron. The basic targets for microlensing searches should be bright point-like radio sources at cosmological distances. In this case, an analysis of their variability and a solid determination of microlensing could lead to an estimation of their cosmological mass density, moreover, in this case one could not exclude a possibility that non-baryonic dark matter also form microlenses if the corresponding optical depth will be high enough. To search for microlensing the most perspective objects are gravitational lensed systems as usually, like CLASS gravitational lens B1600+434, for instance. However, for direct resolving these images and direct detection of apparent motion of the knots, a Radioastron sensitivity have to be improved, since an estimated flux density is too low and to observe the phenomena one should improve sensitivity in 10 times at 6 cm wavelength, for instance, otherwise, it is necessary to increase an integration time (assuming that a radio source has a typical core -- jet structure and microlensing phenomenon is caused superluminal apparent motion of knots). Therefore, a launch of space interferometer Radioastron will give new excellent facilities to investigate microlensing in radio band, since in this case there is a possibility not only to resolve microimages but also observe astrometric microlensing.

A version of the cosmological perturbation theory in general relativity (GR) is developed, where the cosmological scale factor is identified with spatial averaging of the metric determinant logarithm and the cosmic evolution acquires the pattern of a superfluid motion: the absence of "friction-type" interaction, the London-type wave function, and the Bogoliubov condensation of quantum universes. This identification keeps the number of variables of GR and leads to a new type of potential perturbations. A set of arguments is given in favor of that this "superfluid" version of GR is in agreement with the observational data.

Recently, Zakharov et al. (2005) considered the possibility of evaluating the spin parameter and the inclination angle for Kerr black holes in nearby galactic centers by using future advanced astrometrical instruments. A similar approach which uses the characteristic properties of gravitational retro-lensing images can be followed to measure the charge of Reissner-Nordstrom black hole. Indeed, in spite of the fact that their formation might be problematic, charged black holes are objects of intensive investigations. From the theoretical point of view it is well-known that a black hole is described by only three parameters, namely, its mass M, angular momentum J and charge Q. Therefore, it would be important to have a method for measuring all these parameters, preferably by model independent way. In this paper, we propose a procedure to measure the black hole charge by using the size of the retro-lensing images that can be revealed by future astrometrical missions. A discussion of the Kerr-Newmann black hole case is also offered.

Recent X-ray observations of microquasars and Seyfert galaxies reveal the broad emission lines in their spectra, which can arise in the innermost parts of accretion disks. A theoretical analysis of observations and their interpretations were discussed in a number of papers. We consider a radiating annulus model to simulate spectral line shapes. That is a natural approximation for narrow emitting circular rings without extra astrophysical assumptions about emissivity laws. Recently Muller & Camenzind (2004) presented results of their calculations and classified different types of spectral line shapes and described their origin. We clarified their hypothesis about an origin of doubled peaked and double horned line shapes. Based on results of numerical simulations we showed that double peaked spectral lines arise almost for any location of narrow emission rings (annuli) although Muller & Camenzind (2004) suggested that such profiles arise for relatively flat space-times and typical radii for emission region about 25 r_g. We showed that triangular spectral lines could arise for nearest annuli and high inclination angles. We discuss a possibility of appearance of narrow spectral line shapes as a result of spiralling evolution of matter along quasi-circular orbits which could be approximated by narrow annuli.

The observed enhancement of the Fe K$\alpha$ line in three gravitationally lensed QSOs (MG J0414+0534, QSO 2237+0305, H1413+117) is interpreted in terms of microlensing, even when equivalent X-ray continuum amplification is not observed. In order to interpret these observations, first we studied the effects of microlensing on quasars spectra, produced by straight fold caustic crossing over standard relativistic accretion disk. The disk emission was analyzed using the ray tracing method, considering Schwarzschild and Kerr metrics. When the emission is separated in two regions (an inner disk corresponding to the Fe K$\alpha$ line and an outer annulus corresponding to the continuum, or vice-versa) we find microlensing events which enhance the Fe K$\alpha$ line without noticeable amplification of the X-ray continuum, but only during a limited time interval. Continuum amplification is expected if a complete microlensing event is monitored. Second, we studied a more realistic case of amplification by caustic magnification pattern. In this case we could satisfactorily explain the observations if the Fe K$\alpha$ line is emitted from the innermost part of the accretion disk, while the continuum is emitted from a larger region. Also, we studied the chromatic effects of microlensing, finding that the radial distribution of temperature in the accretion disk, combined with microlensing itself, can induce wavelength dependent variability of $\sim$ 30% for microlenses with very small masses. All these results show that X-ray monitoring of gravitational lenses is a well suited method for studying of the innermost structure of AGN accretion disks.

The version of the cosmological perturbation theory based on exact resolution of energy constraint is developed in accordance with the diffeomorphisms of general relativity in the Dirac Hamiltonian approach. Such exact resolution gives one a possibility to fulfil the Hamiltonian reduction and to explain the ``CMBR primordial power spectrum'' and other topical problems of modern cosmology by quantization of the energy constraint and quantum cosmological origin of matter.

Pixel lensing is the gravitational microlensing of light from unresolved stars contributing to the luminosity flux collected by a single pixel. A star must be sufficiently magnified, that is, the lens impact parameter must be less than a threshold value $u_T$ if the excess photon flux in a pixel is to be detected over the background. Assuming the parameters of the Isaac Newton Telescope and typical observing conditions, we present maps in the sky plane towards M31 of threshold impact parameter, optical depth, event number and event time scale, analyzing in particular how these quantities depend on $u_T$ in pixel lensing searches. We use an analytical approach consisting of averaging on $u_T$ and the star column density the optical depth, microlensing rate and event duration time scale. An overall decrease in the expected optical depth and event number with respect to the classical microlensing results is found, particularly towards the high luminosity M31 inner regions. As expected, pixel lensing events towards the inner region of M31 are mostly due to self-lensing, while in the outer region dark events dominate even for a 20% MACHO halo fraction. We also find a far-disk/near-disk asymmetry in the expected event number, smaller than that found by \cite{Kerins04}. Both for self and dark lensing events, the pixel lensing time scale we obtain is $\simeq 1 - 7$ days, dark events lasting roughly twice as long as self-lensing events. The shortest events are found to occur towards the M31 South Semisphere. We also note that the pixel lensing results depend on $<u_T>$ and $<u_T^{~2}>$ values and ultimately on the observing conditions and telescope capabilities.

Observations of AGNs and microquasars by ASCA, RXTE, Chandra and XMM-Newton indicate the existence of broad X-ray emission lines of ionized heavy elements in their spectra. Such spectral lines were discovered also in X-ray spectra of neutron stars and X-ray afterglows of GRBs. Recently, Zakharov et al. (MNRAS, 2003, 342, 1325) described a procedure to estimate an upper limit of the magnetic fields in regions from which X-ray photons are emitted. The authors simulated typical profiles of the iron $K_\alpha$ line in the presence of magnetic field and compared them with observational data in the framework of the widely accepted accretion disk model. Here we further consider typical Zeeman splitting in the framework of a model of non-flat accretion flows, which is a generalization of previous consideration into non-equatorial plane motion of particles emitting X-ray photons. Using perspective facilities of space borne instruments (e.g. Constellation-X mission) a better resolution of the blue peak structure of iron $K_\alpha$ line will allow to evaluate the magnetic fields with higher accuracy.

Recent X-ray observations of microquasars and Seyfert galaxies reveal broad emission lines in their spectra, which can arise in the innermost parts of accretion disks. Simulations indicate that at low inclination angle the line is measured by a distant observer as characteristic two-peak profile. However, at high inclination angles ($> 85^0$) two additional peaks arise. This phenomenon was discovered by Matt et al. (1993) using the Schwarzschild black hole metric to analyze such effect. They assumed that the effect is applicable to a Kerr metric far beyond the range of parameters that they exploited. We check and confirm their hypothesis about such a structure of the spectral line shape for the Kerr metric case. We use no astrophysical assumptions about the physical structure of the emission region except the assumption that the region should be narrow enough. Positions and heights of these extra peaks drastically depend on both the radial coordinate of the emitting region (annuli) and the inclination angle. It was found that these extra peaks arise due to gravitational lens effect in the strong gravitational field, namely they are formed by photons with some number of revolutions around black hole. This conclusion is based only on relativistic calculations without any assumption about physical parameters of the accretion disc like X-ray surface emissivity etc. We discuss how analysis of the iron spectral line shapes could give an information about an upper limit of magnetic field near black hole horizon.

Recently Holz & Wheeler (2002) considered a very attracting possibility to detect retro-MACHOs, i.e. retro-images of the Sun by a Schwarzschild black hole. In this paper we discuss glories (mirages) formed near rapidly rotating Kerr black hole horizons and propose a procedure to measure masses and rotation parameters analyzing these forms of mirages. In some sense that is a manifestation of gravitational lens effect in the strong gravitational field near black hole horizon and a generalization of the retro-gravitational lens phenomenon. We analyze the case of a Kerr black hole rotating at arbitrary speed for some selected positions of a distant observer with respect to the equatorial plane of a Kerr black hole. We discuss glories (mirages) formed near rapidly rotating Kerr black hole horizons and propose a procedure to measure masses and rotation parameters analyzing these forms of mirages. Some time ago Falcke et al. (2000) suggested to search shadows at the Galactic Center. In this paper we present the boundaries for shadows calculated numerically. We also propose to use future radio interferometer RADIOASTRON facilities to measure shapes of mirages (glories) and to evaluate the black hole spin as a function of the position angle of a distant observer.

We consider a contribution of microlensing to the X-ray variability of high-redshifted QSOs. Cosmologically distributed gravitational microlenses could be localized in galaxies (or even in bulge or halo of gravitational macrolenses) or could be distributed in a uniform way. We have analyzed both cases of such distributions. We found that the optical depth for gravitational microlensing caused by cosmologically distributed deflectors could be significant and could reach $10^{-2} - 0.1$ at $z\sim 2$. This means that cosmologically distributed deflectors may contribute significantlly to the X-ray variability of high-redshifted QSOs ($z>2$). Considering that the upper limit of the optical depth ($\tau\sim 0.1$) corresponds to the case where dark matter forms cosmologically distributed deflectors, observations of the X-ray variations of unlensed QSOs can be used for the estimation of the dark matter fraction of microlenses.

We consider a contribution of microlensing to the X-ray variability of high-redshifted QSOs. Such an effect could be caused by stellar mass objects (SMO) located in a bulge or/and in a halo of this quasar as well as at cosmological distances between an observer and a quasar. Cosmologically distributed gravitational microlenses could be localized in galaxies (or even in bulge or halo of gravitational macrolenses) or could be distributed in a uniform way. We have analyzed both cases of such distributions. As a result of our analysis, we obtained that the optical depth for microlensing caused by stellar mass objects is usually small for quasar bulge and quasar halo gravitational microlens distributions. On the other hand, the optical depth for gravitational microlensing caused by cosmologically distributed deflectors could be significant and could reach 0.01 - 0.1 at z ~ 2. This means that cosmologically distributed deflectors may contribute significantlly to the X-ray variability of high-redshifted QSOs (z > 2). Considering that the upper limit of the optical depth (~ 0.1) corresponds to the case where dark matter forms cosmologically distributed deflectors, observations of the X-ray variations of unlensed QSOs can be used for the estimation of the dark matter fraction of microlenses.

Foundations of standard theory of microlensing are described, namely we consider microlensing stars in Galactic bulge, the Magellanic Clouds or other nearby galaxies. We suppose that gravitational microlenses lie between an Earth observer and these stars. Criteria of an identification of microlensing events are discussed. We also consider such microlensing events which do not satisfy these criteria (non-symmetrical light curves, chromatic effects, polarization effects). We describe results of MACHO collaboration observations towards the Large Magellanic Cloud (LMC) and the Galactic bulge. Results of EROS observations towards the LMC and OGLE observations towards the Galactic bulge are also presented. Future microlensing searches are discussed.

We consider a contribution of microlensing in X-ray variability of high-redshifted QSOs. Such an effect could be caused by stellar mass objects (SMO) located in a bulge or/and in a halo of this quasar as well as at cosmological distances between an observer and a quasar. Here, we not considerability of microlensing caused by deflectors in our Galaxy since it is well-known from recent MACHO, EROS and OGLE observations that the corresponding optical depth for the Galactic halo and the Galactic bulge is lower than $10^{-6}$. Cosmologically distributed gravitational microlenses could be localized in galaxies (or even in bulge or halo of gravitational macrolenses) or could be distributed in a uniform way. We have analyzed both cases of such distributions. As a result of our analysis, we obtained that an optical depth for microlensing caused by stellar mass objects is usually small for quasar bulge and quasar halo gravitational microlens distributions ($\tau\sim 10^{-4}$). On the other hand, the optical depth for gravitational microlensing caused by cosmologically distributed deflectors could be significant and could reach $10^{-2} - 0.1$ at $z\sim 2$. It means that cosmologically distributed deflectors may significantlly contribute to the X-ray variability of high-redshifted QSOs ($z>2$). Considering that upper limit of the optical depth ($\tau\sim 0.1$) corresponds to the case when dark matter forms cosmologically distributed deflectors, therefore observations of X-ray variations of unlensend QSOs can be used for the estimation of the dark matter fraction of microlenses.

The possibility that some Markarian objects (e.g. Mkn 501, Mkn 421 and Mkn 766) host massive binary black hole systems with eccentric orbits at their centers has been considered. These systems could be sources of gravitational radiation for space-based gravitational wave interferometers like LISA and ASTROD. In the framework of the Lincoln -- Will approximation we simulate coalescence of such systems, calculate gravitational wave templates and discuss parameters of these binary black hole systems corresponding to the facilities of LISA and ASTROD. We discuss also the possibility to extract information about parameters of the binary black hole systems (masses, of components, distances between them, eccentricity and orbit inclination angle with respect to line of sight) from future gravitational wave measurements.

Recent X-ray observations of Galactic Black Hole Candidates and AGNs and Seyfert galaxies reveal the broad emission lines in their spectra, which can arise in the innermost parts of accretion disks. Simulation indicates that at low inclination angle the line is registered by distant observer with characteristic two-peak profile. However, at high inclination angles ($> 85 \degr$) two additional peaks arise. This phenomenon was discovered by Matt, Perola, Stella (1993) using the Schwarzschild black hole metric to analyze such effect. They assumed that the effect is applicable to a Kerr metric far beyond of a range of parameters that they exploited. We check and confirm their hypothesis about such structure of the spectral line shape for a Kerr metric case. We use no astrophysical assumptions about physical structure of the emission region except assumption that the region should be narrow enough. Positions and heights of these extra peaks drastically depends on both the radial coordinate of the emitting region (circular hot spot) and the inclination angle. We find that best conditions to observe this effect are realized at $\theta > 85\degr$ and $r > 5r_g$ and may exist in microquasars or low-mass black holes in X-ray binary systems, because there is some precession (and nutation of accretion disks) with not very long time periods (see, for example, SS433 binary system). The line profiles for different inclination angles and radial coordinates are presented. To analyze an influence of disk models on the spectral line shapes we simulate the line profiles for Shakura -- Sunyaev disk model for accretion disks with the high inclination.

The original idea to show the spacetime geometry using few geodesics was developed by Johnson and Ruffini (1974). We used this idea to interpret the observational data for rotating BH's. We developed the imitation approach to simulate a propagation of radiation near BH's. An important problem for this approach is the diagnostics of a black hole metric using X-ray observational data of the iron $K_\alpha$-line. Observations of Seyfert galaxies in X-ray region reveal the broad emissiion lines in their spectra, which can arise in inner parts of accretion disks, where the effects of General Relativity (GR) must be counted. A spectrum of a solitary emission line (the $K_\alpha$-line of iron, for example) of a hot spot in Kerr accretion disk is simulated, depending on the radial coordinate $r$ and the angular momentum $a=J/M$ of a black hole, under the assumption of an equatorial circular motion of a hot spot. Using results of numerical simulations it is shown that the characteristic two-peak line profile with the sharp edges arises at a large distance, (about $r \approx (3-10)r_g$). The inner regions emit the line, which is observed with one maximum and extremely broad red wing. High accuracy future spectral observations, being carried out, could detect the angular momentum $a$ of the black hole. We analyzed the different parameters of problems on the observable shape of this line and discussed some possible kinds of these shapes. The total number of geodesics is about $10^9$ (to simulate possible shapes of the $K_\alpha$-line), so the number is great enough, especially in comparison with few geodesics in the original paper by Johnson and Ruffini (1974).

Observations of AGNs and microquasars by ASCA, RXTE, Chandra and XMM-Newton indicate the existence of wide X-ray emission lines of heavy ionized elements in their spectra. The emission can arise in the inner parts of accretion discs where the effects of General Relativity (GR) must be counted, moreover such effects can dominate. We describe a procedure to estimate an upper limit of the magnetic fields in the regions where X-ray photons are emitted. We simulate typical profiles of the iron $K_\alpha$ line in the presence of a magnetic field and compare them with observational data. As an illustration we find $H < 10^{10} - 10^{11}$ Gs for Seyfert galaxy MCG--6--30--15. Using the perspective facilities of measurement devices (e.g. Constellation-X mission) a better resolution of the blue peak structure of iron $K_\alpha$ line will allow to find the value of the magnetic fields if the latter are high enough.

Foundations of standard theory of microlensing are described, namely we consider microlensing stars in Galactic bulge, the Magellanic Clouds or other nearby galaxies. We suppose that gravitational microlenses lie between an Earth observer and these stars. Criteria of an identification of microlensing events are discussed. We also consider such microlensing events which do not satisfy these criteria (non-symmetrical light curves, chromatic effects, polarization effects). We describe results of MACHO collaboration observations towards the Large Magellanic Cloud (LMC) and the Galactic bulge in detail. Results of EROS observations towards the LMC and OGLE observations towards the Galactic bulge are also presented. A comparison of the microlensing theory and observations is discussed in full.