In this work, we report some results on the numerical exploration of the model of Guendelman-Nissimov-Pacheva. This model has been previously applied to cosmology, but there were open questions regarding its parameters. Here we demonstrate the existence of families of solutions on the slope of the effective potential which preserve the duration of the inflation and its power. For this solutions, one can see the previously reported phenomenon of the inflaton scalar field climbing up the slope, with the effect more pronounced when starting lower on the potential slope. Finally we compare the dynamical and the potential slow-roll parameters for the model and we find that the latter describe the numerically observed inflationary period better.

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.

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.

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.

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.

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.