The model of minimal dilatonic gravity (MDG), called also the massive Branse-Dicke model with $\omega =0$, is an alternative model of gravitation, which uses one Branse-Dicke gravitation-dilaton field $\Phi$ and offers a simultaneous explanation of the effects of dark energy (DE) and dark matter (DM). Here we present an extensive research of non-rotating neutron star models in MDG with four different realistic equations of state (EOS), which are in agreement with the latest observational data. The equations describing static spherically symmetric stars in MDG are solved numerically. The effects corresponding to DE and DM are clearly seen and discussed.

Using an effective one body approach we describe in detail gravitational waves from classical three body problem on a non-rotating straight line and derive their basic physical characteristics. Special attention is paid to the irregular motions of such systems and to the significance of double and triple collisions. The conclusive role of the collinear solutions is also discussed in short. It is shown that the residuals may contain information about irregular motion of the source of gravitational waves.

The proper understanding of the electromagnetic counterpart of gravity-waves emitters is of serious interest to the multimessenger astronomy. In this article, we study the numerical stability of the quasinormal modes (QNM) and quasibound modes (QBM) obtained as solutions of the Teukolsky Angular Equation and the Teukolsky Radial Equation with appropriate boundary conditions. We use the epsilon-method for the system featuring the confluent Heun functions to study the stability of the spectra with respect to changes in the radial variable. We find that the QNM and QBM are stable in certain regions of the complex plane, just as expected, while the third "spurious" spectrum was found to be numerically unstable and thus unphysical. This analysis shows the importance of understanding the numerical results in the framework of the theory of the confluent Heun functions, in order to be able to distinguish the physical spectra from the numerical artifacts.

We present a derivation of the basic equations and boundary conditions for relativistic static spherically symmetric stars (SSSS) in the model of minimal dilatonic gravity (MDG) which offers an alternative and simultaneous description of the effects of dark matter (DM) and dark energy (DE) using one dilaton field $\Phi$. The numerical results for a realistic equation of state (EOS) MPA1 of neutron matter are presented for the first time. The three very different scales, the Compton length of the scalar field $\lambda_\Phi$, the star's radius $r^*$, and the finite radius of the MDG Universe $r_{U}$ are a source of numerical difficulties. Owing to the introduction of a new dark scalar field $\varphi=\ln(1+\ln\Phi)$, we have been able to study numerically an unprecedentedly large interval of $\lambda_\Phi$ and have discovered the existence of $\lambda_\Phi^{crit}\approx 2.1$\ km for a neutron star with MPA1 EOS. This is related to the bifurcation of the physical domain in the phase space of the system. Some novel physical consequences are discussed.

The perturbations of the Kerr metric and the miracle of their exact solutions play a critical role in the comparison of predictions of general relativity with astrophysical observations of compact massive objects. The differential equations governing the late-time ring-down of the perturbations of the Kerr metric, the Teukolsky Angular Equation and the Teukolsky Radial Equation, can be solved analytically in terms of confluent Heun functions. In this article, we solve numerically the spectral system formed by those exact solutions and we obtain the electromagnetic (EM) spectra of the Kerr black hole. Because of the novel direct way of imposing the boundary conditions, one is able to discern three different types of spectra: the well-known quasinormal modes (QNM), the symmetric with respect to the real axis quasibound modes (QBM) and a spurious spectrum who is radially unstable. This approach allows clearer justification of the term "spurious" spectrum, which may be important considering the recent interest in the spectra of the electromagnetic counterparts of events producing gravitational waves.

We present solution of the equations for relativistic static spherically symmetric stars (SSSS) in the model of minimal dilatonic gravity (MDG) using the polytropic equation of state. A polytropic equation of state, which has a good fitting with a more realistic one, is used. Results are obtained for all variables of a single neutron star in the model of MDG. The maximum mass about two solar masses is in accordance with the latest observations of pulsars. Several new effects are observed for the variables related with the dilaton $\Phi$ and the cosmological constant $\Lambda$. The mass-radius relation is also obtained. Special attention is paid to the behavior of the quantities which describe the effects analogous to those of dark energy and dark matter in MDG. The results of the present paper confirm the conclusion that the dilaton $\Phi$ is able to play simultaneously the role of dark energy and dark matter.

We present the basic equations and relations for the relativistic static spherically symmetric stars (SSSS) in the model of minimal dilatonic gravity (MDG) which is {\em locally} equivalent to the f(R) theories of gravity and gives an alternative description of the effects of dark matter and dark energy. The results for the simplest form of the relativistic equation of state (EOS) of neutron matter are represented. Our approach overcomes the well-known difficulties of the physics of SSSS in the f(R) theories of gravity introducing two novel EOS for cosmological energy-pressure densities and dilaton energy-pressure densities, as well as proper boundary conditions.

In the present article we introduce and study a novel type of solutions of the general Heun's equation. Our approach is based on the symmetric form of the Heun's differential equation yielded by development of the Felix Klein symmetric form of the Fuchsian equations with an arbitrary number $N\geq 4$ of regular singular points. We derive the symmetry group of these equations which turns to be a proper extension of the Mobius group. We also introduce and study new series solution of symmetric form of the general Heun's differential equation (N=4) which treats simultaneously and on an equal footing all singular points. Hopefully, this new form will simplify the resolution of the existing open problems in the theory of general Heun's functions and can be used for development of new effective computational methods.

We study the relation between Minimal Dilatonic Gravity (MDG) and f(R) theories of gravity and establish strict conditions for their {\em global} equivalence. Such equivalence takes place only for a certain class of cosmological potentials, dubbed here {\em withholding potentials}, since they prevent change of the sign of dilaton $\Phi$. The withholding property ensures the attractive character of gravity, as well as absence of ghosts and a tachyon in the gravi-dilaton sector and yields certain asymptotic of the admissible functions $f(R)$. Large classes of withholding cosmological potentials and functions $f(R)$ are found and described in detail. It is shown that the popular choices of $f(R)$ functions are not withholding ones. The particle content of the gravi-dilaton sector is found using perturbation theory around de Sitter vacuum of MDG. The graviton remains massless, since it obeys conformal invariant field equation in the de Sitter space-time. The $R/6$ term in the conformal invariant wave operator introduces a very small mass scale $m_{R} \approx 1.5 \times 10^{-38} m_e$, $m_e$ being the mass of the electron. The mass of the dilaton is much larger: $m\gtrapprox 10^{29} m_{R}$. Two new phenomena: scalaron waves and induction of gravitational waves by the scalaron field are discussed using the derived wave equations for scalaron and graviton. The MDG and f(R) theories are shown to predict physical deviations from GR. Seemingly, the MDG and f(R) theories, when globally equivalent, offer a unified description of dark energy and dark matter.

The differential equations governing the late-time ring-down of the perturbations of the Kerr metric, the Teukolsky Angular Equation and the Teukolsky Radial Equation, can be solved analytically in terms of confluent Heun functions. In this article, for the first time, we use those exact solutions to obtain the electromagnetic (EM) quasinormal spectra of the Kerr black hole . This is done by imposing the appropriate boundary conditions on the solutions and solving numerically the so-obtained two-dimensional transcendental system. The EM QNM spectra are compared with already published results, evaluated through the continued fraction method. The comparison shows that the modes with lower $n$ coincide for both methods, while those with higher $n$ may demonstrate significant deviations. To study those deviations, we employ the $\epsilon$-method which introduces small variations in the argument of the complex radial variable. Using the $\epsilon$-method one can move in the complex $r-$plane the branch cuts of the solutions of the radial equation and examine the dependence of the spectrum on their position. For different values of $\epsilon$, one can obtain both the frequencies evaluated through the well-established continued fraction method or somewhat different spectra calculated here for the first time. This result raises the question what spectrum should be compared with the observational data and why. This choice should come from better understanding of the physics of the problem and it may become particularly important considering the recent interest in the spectra of the electromagnetic counterparts of events producing gravitational waves.

The Heun functions have wide application in modern physics and are expected to succeed the hypergeometrical functions in the physical problems of the 21st century. The numerical work with those functions, however, is complicated and requires filling the gaps in the theory of the Heun functions and also, creating new algorithms able to work with them efficiently. We propose a new algorithm for solving a system of two nonlinear transcendental equations with two complex variables based on the M\"uller algorithm. The new algorithm is particularly useful in systems featuring the Heun functions and for them, the new algorithm gives distinctly better results than Newton's and Broyden's methods. As an example for its application in physics, the new algorithm was used to find the quasi-normal modes (QNM) of Schwarzschild black hole described by the Regge-Wheeler equation. The numerical results obtained by our method are compared with the already published QNM frequencies and are found to coincide to a great extent with them. Also discussed are the QNM of the Kerr black hole, described by the Teukolsky Master equation.

The paper presents a generalization and further development of our recent publications where solutions of the Klein-Fock-Gordon equation defined on a few particular $D=(2+1)$-dim static space-time manifolds were considered. The latter involve toy models of 2-dim spaces with axial symmetry, including dimension reduction to the 1-dim space as a singular limiting case. Here the non-static models of space geometry with axial symmetry are under consideration. To make these models closer to physical reality, we define a set of "admissible" shape functions $\rho(t,z)$ as the $(2+1)$-dim Einstein equations solutions in the vacuum space-time, in the presence of the $\Lambda$-term, and for the space-time filled with the standard "dust". It is curious that in the last case the Einstein equations reduce to the well-known Monge-Amp\`{e}re equation, thus enabling one to obtain the general solution of the Cauchy problem, as well as a set of other specific solutions involving one arbitrary function. A few explicit solutions of the Klein-Fock-Gordon equation in this set are given. An interesting qualitative feature of these solutions relates to the dimension reduction points, their classification, and time behavior. In particular, these new entities could provide us with novel insight into the nature of P- and T-violation, and of Big Bang. A short comparison with other attempts to utilize dimensional reduction of the space-time is given.

Although finding numerically the quasinormal modes of a nonrotating black hole is a well-studied question, the physics of the problem is often hidden behind complicated numerical procedures aimed at avoiding the direct solution of the spectral system in this case. In this article, we use the exact analytical solutions of the Regge-Wheeler equation and the Teukolsky radial equation, written in terms of confluent Heun functions. In both cases, we obtain the quasinormal modes numerically from spectral condition written in terms of the Heun functions. The frequencies are compared with ones already published by Andersson and other authors. A new method of studying the branch cuts in the solutions is presented -- the epsilon-method. In particular, we prove that the mode $n=8$ is not algebraically special and find its value with more than 6 firm figures of precision for the first time. The stability of that mode is explored using the $\epsilon$ method, and the results show that this new method provides a natural way of studying the behavior of the modes around the branch cut points.

We propose a new algorithm for solving a system of two nonlinear transcendental equations with two complex variables based on the Muller algorithm. The two-dimensional Muller algorithm is tested on systems of different type and is found to work comparably to Newton's method and Broyden's method in many cases. The new algorithm is particularly useful in systems featuring the Heun functions whose complexity may make the already known algorithms not efficient enough or not working at all. In those specific cases, the new algorithm gives distinctly better results than the other two methods. As an example for its application in physics, the new algorithm was used to find the quasi-normal modes (QNM) of Schwarzschild black hole described by the Regge-Wheeler equation. The numerical results obtained by our method are compared with the already published QNM frequencies and are found to coincide to a great extent with them. Also discussed are the QNM of the Kerr black hole, described by the Teukolsky Master equation.

We give a new theoretical basis for examination of the presence of the Kerr black hole (KBH) or the Kerr naked singularity (KNS) in the central engine of different astrophysical objects around which astrophysical jets are typically formed: X-ray binary systems, gamma ray bursts (GRBs), active galactic nuclei (AGN), etc. Our method is based on the study of the exact solutions of the Teukolsky master equation for electromagnetic perturbations of the Kerr metric. By imposing original boundary conditions on the solutions so that they describe a collimated electromagnetic outflow, we obtain the spectra of possible {\em primary jets} of radiation, introduced here for the first time. The theoretical spectra of primary electromagnetic jets are calculated numerically. Our main result is a detailed description of the qualitative change of the behavior of primary electromagnetic jet frequencies under the transition from the KBH to the KNS, considered here as a bifurcation of the Kerr metric. We show that quite surprisingly the novel spectra describe linearly stable primary electromagnetic jets from both the KBH and the KNS. Numerical investigation of the dependence of these primary jet spectra on the rotation of the Kerr metric is presented and discussed.

Weak gravitational, electromagnetic, neutrino and scalar fields, considered as perturbations on Kerr background satisfy Teukolsky Master Equation. The two non-trivial equations obtained after separating the variables are the polar angle equation and the radial equation. We solve them by transforming each one into the form of a confluent Heun equation. The transformation depends on a set of parameters, which can be chosen in a such a way, so the resulting equations have simple polynomial solutions for neutrino, electromagnetic, and gravitational perturbations, provided some additional conditions are satisfied. Remarkably there exists a class of solutions for which these additional conditions are the same for the two different equations for $|s|=1/2$ and $|s|=1$. As a result the additional conditions fix the dependence of the separation constant on the angular frequency but the frequency itself remains unconstrained and belongs to a continuous spectrum.

The Teukolsky Master Equation is the basic tool for study of perturbations of the Kerr metric in linear approximation. It admits separation of variables, thus yielding the Teukolsky Radial Equation and the Teukolsky Angular Equation. We present here a unified description of all classes of exact solutions to these equations in terms of the confluent Heun functions. Large classes of new exact solutions are found and classified with respect to their characteristic properties. Special attention is paid to the polynomial solutions which are singular ones and introduce collimated one-way-running waves. It is shown that a proper linear combination of such solutions can present bounded one-way-running waves. This type of waves may be suitable as models of the observed astrophysical jets.

We present new developments of the simple model of the central engine of GRB, proposed recently. The model is based on minimal assumptions: some rotating compact relativistic object at the center and stable perturbations of its rotating gravitational field, described by Teukolsky Master Equation. We show that using nonstandard polynomial solutions to the angular Teukolsky equation we can describe the formation of collimated jets of various forms. Appearance of imaginary part of the superradiance-like frequency is established for the first time for pure vacuum black hole jet solutions of Teukolsky equation.

Despite all the already existing observational data, current models still cannot explain completely the excessive energy output and the time variability of GRB. One of the reasons for this is the lack of a good model of the central engine of GRB. A major problem in the proposed models with a black hole (BH) in the center is that they don't explain the observed evidences of late time activity of the central engine. In this paper we are starting the search for a possible model of that central engine as a rotating compact body of still unknown nature. The formation of jets in the new model lies entirely on the fundamental Teukolsky Master Equation. We demonstrate that this general model can describe the formation of collimated GRB-jets of various forms. Some preliminary results are presented.

We describe a novel class of geometrical models of relativistic stars. Our approach to the static spherically symmetric solutions of Einstein equations is based on a careful physical analysis of radial gauge conditions. It brings us to a two parameter family of relativistic stars without stiff functional dependence between the stelar radius and stelar mass. As a result, a point particle idealization -- a limiting case of bodies with finite dimension, becomes possible in GR, much like in Newtonian gravity. We devote this article to detailed mathematical study of this limit.

In a series of articles we describe a novel class of geometrical models of relativistic stars. Our approach to the static spherically symmetric solutions of Einstein equations is based on a careful physical analysis of radial gauge conditions. It turns out that there exist heavy black dwarfs: relativistic stars with arbitrary large mass, which are to have arbitrary small radius and arbitrary small luminosity. In the present article we mathematically prove this new phenomena, using a detailed consideration of incompressible GR stars. We study the whole two parameter family of solutions of extended TOV equations for incompressible stars. This example is used to illustrate most of the basic features of the new geometrical models of relativistic stars. Comparison with newest observational data is discussed

In a series of articles we describe a novel class of geometrical models of relativistic stars. Our approach to the static spherically symmetric solutions of Einstein equations is based on a careful physical analysis of radial gauge conditions. It brings us to a two parameter family of relativistic stars without stiff functional dependence between the stelar radius and stelar mass. It turns out that within this family there do exist relativistic stars with arbitrary large mass, which are to have arbitrary small radius and arbitrary small luminosity. In addition, point particle idealization, as a limiting case of bodies with finite dimension, becomes possible in GR, much like in Newton gravity.

We present a class of simple scalar-tensor models of gravity with one scalar field (dilaton $\Phi$) and only one unknown function (cosmological potential $U(\Phi)$). These models might be considered as a stringy inspired ones with broken SUSY. They have the following basic properties: 1) Positive dilaton mass, $m_\Phi$, and positive cosmological constant $\Lambda$, define two extremely different scales. The models under consideration are consistent with the known experimental facts if $m_\Phi > 10^{-3} eV$ and $\Lambda=\Lambda^{obs}\sim 10^{-56} cm^{-2}$. 2) Einstein week equivalence principle is strictly satisfied and extended to scalar-tensor theories of gravity using a novel form of principle of "constancy of fundamental constants". 3) The dilaton plays simultaneously role of inflation field and quintessence field and yields a sequential hyper-inflation with graceful exit to asymptotic de Sitter space-time which is an attractor, and is approached as $\exp(-\sqrt{3\Lambda^{obs}} ct/2)$. The time duration of inflation is $\Delta t_{infl} \sim m_\Phi^{-1}$. 4) Ultra-high frequency ($\omega_\Phi \sim m_\Phi$) dilatonic oscillations take place in asymptotic regime. 5) No fine tuning. (The Robertson-Walker solutions of general type have the above properties.) 6) A novel adjustment mechanism for cosmological constant problem seems to be possible: the huge value of cosmological constant in the stringy frame is re-scaled to its observed value by dilaton after transition to phenomenological frame.

We numerically construct static and spherically symmetric electrically charged black hole solutions in Einstein-Born-Infeld gravity with massive dilaton. The numerical solutions show that the dilaton potential allows many more black hole causal structures than the massless dilaton. We find that depending on the black hole mass and charge and the dilaton mass the black holes can have either one, two, or three horizons. The extremal solutions are also found out. As an interesting peculiarity we note that there are extremal black holes with an inner horizon and with triply degenerated horizon.

We investigate numerically models of the static spherically symmetric boson-fermion stars in the scalar-tensor theory of gravity with massive dilaton field. The proper mathematical model of such stars is interpreted as a nonlinear two-parametric eigenvalue problem with unknown internal boundary. To solve this problem the Continuous Analogue of Newton Method is used.

We consider: minimal scalar-tensor model of gravity with Brans-Dicke factor $\omega(\Phi)\equiv 0$ and cosmological factor $\Pi(\Phi)$; restrictions on it from gravitational experiments; qualitative analysis of new approach to cosmological constant problem based on the huge amount of action in Universe; determination of $\Pi(\Phi)$ using time evolution of scale factor of Universe.

We study a new minimal scalar-tensor model of gravity with Brans-Dicke factor $\omega(\Phi)\equiv 0$ and cosmological factor $\Pi(\Phi)$. The constraints on $\Pi(\Phi)$ from known gravitational experiments are derived. We show that almost any time evolution of the scale factor in a homogeneous isotropic Universe can be obtained via properly chosen $\Pi(\Phi)$ and discuss the general properties of models of this type.

In the framework of a model of minimal of dilatonic gravity (MDG) with cosmological potential we consider: the relations of MDG with nonlinear gravity and string theory; natural cosmological units, defined by cosmological constant; the properties of cosmological factor, derived from solar system and Earth-surface gravitational experiments; universal anty-gravitational interactions, induced by positive cosmological constant and by Nordtved effect; a new formulation of cosmological constant problem using the ratio of introduced cosmological action and Planck constant $\sim 10 ^{122}$;qualitative analysis of this huge number based on classical action of effective Bohr hydrogen atoms; inverse cosmological problem: to find cosmological potential which yields given evolution of the RW Universe; and comment other general properties of MDG.

We develop the general theory of stars in Saa's model of gravity with propagating torsion and study the basic stationary state of neutron star. Our numerical results show that the torsion force decreases the role of the gravity in the star configuration leading to significant changes in the neutron star masses depending on the equation of state of star matter. The inconsistency of the Saa's model with Roll-Krotkov-Dicke and Braginsky-Panov experiments is discussed.

We propose a novel self consistent minimal coupling principle in presence of torsion dilaton field. This principle yields a new local dilatation symmetry and predicts the interactions of torsion dilaton with the real matter and with metric. The soft violation of this symmetry yields a physical dilaton and a simple relation between Cartan scalar curvature and cosmological constant in this new model of gravity with propagating torsion. Its relation with scalar-tensor theories of gravity and a possible use of torsion dilaton in the inflation scenario is discussed. \noindent{PACS number(s): 04.50.+h, 04.40.Nr, 04.62.+v}

We present a review of some recent models of gravitation theory with propagating torsion based on the use of a torsion-dilaton field and propose one more model of this type which promises to be more realistic. A proper universal self-consistent minimal action principle yields the properties of this model and predicts the interactions of torsion-dilaton field with the real matter. The new model may be compatible with the string models with dilaton field and gives a novel interpretation of the dilaton as a part of the space-time torsion. A relation with some recent models of dilatonic gravity is also possible.

It is shown that the recently proposed interpretation of the transposed equi-affine theory of gravity as a theory with variable Plank "constant" is inconsistent with basic solar system gravitational experiments.