We investigate inflationary dynamics in the framework of the Einstein-Gauss-Bonnet gravity. In the model under consideration, the inflaton field is non-minimally coupled to the Gauss-Bonnet curvature invariant, so that the latter appears to be dynamically important. We consider a quartic potential for the inflaton field, in particular the one asymptotically connected to the Higgs inflation, and a wider class of coupling functions not considered in the earlier work. Keeping in mind the observational bounds on the parameters - the amplitude of scalar perturbations $A_s$, spectral index $n_s$ and tensor-to-scalar ratio $r$, we demonstrate that the model a quartic potential and the proposed coupling function is in agreement with observation.

We construct the $F(R)$ gravity models with exact particular solutions using the conformal transformation and the superpotential method for the corresponding models in the Einstein frame. The functions $F(R)$ are obtained explicitly. We consider exact solutions for the obtained $R^2$ gravity model with the cosmological constant in detail.

We investigate the scalar field dynamics of models with nonminimally coupled scalar fields in the presence of the Gauss-Bonnet term and derive the structure of the effective potential and conditions for stable de Sitter solutions in general. Specializing to specific couplings, we explore the possibility of realizing the stable de Sitter configurations which may have implications for both the early Universe and late time evolution.

We investigate the cosmological dynamics of nonlocally corrected gravity involving a function of the inverse d'Alembertian acting on the Gauss-Bonnet term. Casting the dynamical equations in local form, we derive the reconstruction procedure. We find conditions on the model parameters that are sufficient for the existence of de Sitter solutions and obtain these solutions explicitly.

We consider the Chiral Cosmological Models (CCMs) and modified gravity theories associated with them. Generalization of the superpotential method for a general CCM with several scalar fields is performed, and the method of construction CCMs admitting exact solutions is developed. New classes of exact solutions in the two-component CCM connected with an $f(R)$ gravity model with an additional scalar field have been constructed. We construct new cosmological solutions for a diagonal metric of the target space, including modified power-law solutions. In particular, we propose the reconstruction procedure based on the superpotential method and present examples of kinetic part reconstruction for periodic and hyperbolic Hubble parameters. We also focus on a cyclic type of Universe dubbed the Quasi-Steady State (QSS) model, with the aim of constructing single- and double-field potentials for one and the same behaviour of the Hubble parameter using the developed superpotential method for the CCM. The realization of this task includes a new set of solutions for a CCM with a scale factor characterized by the QSS theory. We also propose a method for reducing the two-field CCM to the single scalar field model.

The cosmological dynamics of a non-locally corrected gravity theory, involving a power of the inverse d'Alembertian, is investigated. Casting the dynamical equations into local form, the fixed points of the models are derived, as well as corresponding de Sitter and power-law solutions. Necessary and sufficient conditions on the model parameters for the existence of de Sitter solutions are obtained. The possible existence of power-law solutions is investigated, and it is proven that models with de Sitter solutions have no power-law solutions. A model is found, which allows to describe the matter-dominated phase of the Universe evolution.

Inflationary scenarios motivated by the Minimal Supersymmetric Standard Model (MSSM) where five scalar fields are non-minimally coupled to gravity are considered. The potential of the model and the function of non-minimal coupling are polynomials of two Higgs doublet convolutions. We show that the use of the strong coupling approximation allows to obtain inflationary parameters in the case when a combination of the four scalar fields plays a role of inflaton. Numerical calculations show that the cosmological evolution leads to inflationary scenarios fully compatible with observational data for different values of the MSSM mixing angle $\beta$.

We study the cosmological evolution and singularity crossing in the Bianchi-I universe filled with a conformally coupled scalar field and compare them with those of the Bianchi-I universe filled with a minimally coupled scalar field. We also write down the solution for the Bianchi-I Universe in the induced gravity cosmology.

Despite the fact that experimentally with a high degree of statistical significance only a single Standard Model--like Higgs boson is discovered at the LHC, extended Higgs sectors with multiple scalar fields not excluded by combined fits of the data are more preferable theoretically for internally consistent realistic models of particle physics. We analyze the inflationary scenarios which could be induced by the two-Higgs doublet potential of the Minimal Supersymmetric Standard Model (MSSM) where five scalar fields have nonminimal couplings to gravity. Observables following from such MSSM-inspired multifield inflation are calculated and a number of consistent inflationary scenarios are constructed. Cosmological evolution with different initial conditions for the multifield system leads to consequences fully compatible with observational data on the spectral index and the tensor-to-scalar ratio. It is demonstrated that the strong coupling approximation is precise enough to describe such inflationary scenarios.

We study dynamics of induced gravity cosmological models with the sixth degree polynomial potentials, that have been constructed using the superpotential method. We find conditions on the potential under which exact bounce solutions exist and study the stability of these solutions.

We consider a rather simple method for the description of the Big Bang - Big Crunch cosmological singularity crossing. For the flat Friedmann universe this method gives the same results as more complicated methods, using Weyl symmetry or the transitions between the Jordan and Einstein frames. It is then easily generalized for the case of a Bianchi-I anisotropic universe. We also present early-time and late-time asymptotic solutions for a Bianchi-I universe, filled with a conformally coupled massless scalar field.

We explore dynamics of cosmological models with bounce solutions evolving on a spatially flat Friedmann-Lemaitre-Robertson-Walker background. We consider cosmological models that contain the Hilbert-Einstein curvature term, the induced gravity term with a negative coupled constant, and even polynomial potentials of the scalar field. Bounce solutions with non-monotonic Hubble parameters have been obtained and analyzed. The case when the scalar field has the conformal coupling and the Higgs-like potential with an opposite sign is studied in detail. In this model the evolution of the Hubble parameter of the bounce solution essentially depends on the sign of the cosmological constant.

We remind the way to obtain integrable models with non-minimally coupled scalar fields. We are interesting to models with bounce solutions and compare bounce solutions in two known integrable models. We show that only one model has a bounce solution that tends to a stable de Sitter solution.

We study the relation between the Jordan-Einstein frame transition and the possible description of the crossing of singularities in flat Friedmann universes, using the fact that the regular evolution in one frame can correspond to crossing singularities in the other frame. We show that some interesting effects arise in simple models such as one with a massless scalar field or another wherein the potential is constant in the Einstein frame. The dynamics in these models and in their conformally coupled counterparts are described in detail, and a method for the continuation of such cosmological evolutions beyond the singularity is developed. We compare our approach with some other, recently developed, approaches to the problem of the crossing of singularities.

We study integrable models in the Bianchi I metric case with scalar fields minimally and non-minimally coupled with gravity and the correspondence between their general solutions. Using the model with a minimally coupled scalar field and a constant potential as an example, we demonstrate how to obtain the general solutions of the corresponding models in the Jordan frame.

The possibility to construct an inflationary scenario for renormalization-group improved potentials corresponding to the Higgs sector of quantum field models is investigated. Taking into account quantum corrections to the renormalization-group potential which sums all leading logs of perturbation theory is essential for a successful realization of the inflationary scenario, with very reasonable parameters values. The scalar electrodynamics inflationary scenario thus obtained are seen to be in good agreement with the most recent observational data.

We study the integrable model with minimally and non-minimally coupled scalar fields and the correspondence of their general solutions. Using the model with a minimally coupled scalar field and a the constant potential as an example we demonstrate the difference between the general solutions of the corresponding models in the Jordan and the Einstein frames.

The possibility to construct an inflationary scenario for renormalization-group improved potentials corresponding to the Higgs sector of finite gauge models is investigated. Taking into account quantum corrections to the renormalization-group potential which sums all leading logs of perturbation theory is essential for a successful realization of the inflationary scenario, with very reasonable parameter values. The inflationary models thus obtained are seen to be in good agreement with the most recent and accurate observational data. More specifically, the values of the relevant inflationary parameters, $n_s$ and $r$, are close to the corresponding ones in the $R^2$ and Higgs-driven inflation scenarios. It is shown that the model here constructed and Higgs-driven inflation belong to the same class of cosmological attractors.

We consider the relation between exact solutions of cosmological models having minimally and non-minimally coupled scalar fields. This is done for a particular class of solvable models which, in the Einstein frame, have potentials depending on hyperbolic functions and in the Jordan frame, where the non-minimal coupling is conformal, possess a relatively simple dynamics. We show that a particular model in this class can be generalized to the cases of closed and open Friedmann universes and still exhibits a simple dynamics. Further we illustrate the conditions for the existences of bounces in some sub-classes of the set of integrable models we have considered.

We explore dynamics of cosmological models with a nonminimally coupled scalar field evolving on a spatially flat Friedmann-Lemaitre-Robertson-Walker background. We consider cosmological models including the Hilbert-Einstein curvature term and the $N$ degree monomial of the scalar field nonminimally coupled to gravity. The potential of the scalar field is the $n$ degree monomial or polynomial. We describe several qualitatively different types of dynamics depending on values of power indices $N$ and $n$. We identify that three main possible pictures correspond to $n<N$, $N<n<2N$ and $n>2N$ cases. Some special features connected with the important cases of $N=n$ (including the quadratic potential with quadratic coupling) and $n=2N$ (which shares its asymptotics with the potential of the Higgs-driven inflation) are described separately. A global qualitative analysis allows us to cover the most interesting cases of small $N$ and $n$ by a limiting number of phase-space diagrams. The influence of the cosmological constant to the global features of dynamics is also studied.

The possibility to construct inflationary models for the renormalization-group improved potentials corresponding to scalar electrodynamics and to $SU(2)$ and $SU(5)$ models is investigated. In all cases, the tree-level potential, which corresponds to the cosmological constant in the Einstein frame, is seen to be non-suitable for inflation. Rather than adding the Hilbert-Einstein term to the action, quantum corrections to the potential, coming from to the RG-equation, are included. The inflationary scenario is analyzed with unstable de Sitter solutions which correspond to positive values of the coupling function, only. We show that, for the finite $SU(2)$ model and $SU(2)$ gauge model, there are no de Sitter solutions suitable for inflation, unless exit from it occurs according to some weird, non-standard scenarios. Inflation is realized both for scalar electrodynamics and for $SU(5)$ RG-improved potentials, and the corresponding values of the coupling function are seen to be positive. It is shown that, for quite reasonable values of the parameters, the inflationary models obtained both from scalar electrodynamics and from the $SU(5)$ RG-improved potentials, are in good agreement with the most recent observational data coming from the Planck2013 and BICEP2 collaborations.

A non-local modified gravity model with an analytic function of the d'Alembert operator that has been proposed as a possible way of resolving the singularities problems in cosmology is considered. We show that the anzats that is usually used to obtain exact solutions in this model provides a connection with $f(R)$ gravity models.

We obtain general solutions for some flat Friedmann universes filled with a scalar field in induced gravity models and models including the Hilbert-Einstein curvature term plus a scalar field conformally coupled to gravity. As is well known, these models are connected to minimally coupled models through the combination of a conformal transformation and a transformation of the scalar field. The explicit forms of the self-interaction potentials for six exactly solvable models are presented here. We obtain the general solution for one of the integrable models, namely, the induced gravity model with a power-law potential for the self-interaction of the scalar field. We argue that although being mathematically in a one-to-one correspondence with the solutions in the minimally coupled models, the solutions in the corresponding non-minimally coupled models are physically different. This is because the cosmological evolutions seen by an internal observer connected with the cosmic time can be quite different. The study of a few induced gravity models with particular potentials gives us an explicit example of such a difference.

We study dynamics of induced gravity cosmological models with sixth degree potential, that have found using the superpotential method. The important property of these models are existence of exact cosmological solutions that tend to fixed points. The stability of these cosmological solutions have been obtained. In particular, we find conditions under which solutions with a non-monotonic Hubble parameter that tend to a fixed point are attractors.

We study dynamics of non-minimally coupled scalar field cosmological models with Higgs-like potentials and a negative cosmological constant. In these models the inflationary stage of the Universe evolution changes into a quasi-cyclic stage of the Universe evolution with oscillation behaviour of the Hubble parameter from positive to negative values. Depending on the initial conditions the Hubble parameter can perform either one or several cycles before to become negative forever.

A nonlocal gravity model which does not assume the existence of a new dimensional parameter in the action and includes a function $f(\Box^{-1} R)$, with $\Box$ the d'Alembertian operator, is studied. By specifying an exponential form for the function f and including a matter sector with a constant equation of state parameter, all available power-law solutions in the Jordan frame are obtained. New power-law solutions in the Einstein frame are also probed. Furthermore, the relationship between power-law solutions in both frames, established through conformal transformation, is substantially clarified. The correspondence between power-law solutions in these two frames is proven to be a very useful tool in order to obtain new solutions in the Einstein frame.

A nonlocal gravity model is considered which does not assume the existence of a new dimensional parameter in the action and includes a function $f(\Box^{-1} R)$, with $\Box$ the d'Alembertian operator. Using a reconstruction procedure for the local scalar-tensor formulation of this model, a function f is obtained for which the model exhibits power-law solutions with the Hubble parameter $H=n/t$, for any value of the constant n. For generic n - namely except for a few special values which are characterized and also specifically studied - the corresponding function f is a sum of exponential functions. Corresponding power-law solutions are found explicitly. Also the case is solved in all detail of a function f such that the model contains both de Sitter and power-law solutions.

We employ the superpotential technique for the reconstruction of cosmological models with a non-minimally coupled scalar field evolving on a spatially flat Friedmann-Robertson-Walker background. The key point in this method is that the Hubble parameter is considered as a function of the scalar field and this allows one to reconstruct the scalar field potential and determine the dynamics of the field itself, without a priori fixing the Hubble parameter as a function of time or of the scale factor. The scalar field potentials that lead to de Sitter or asymptotic de Sitter solutions, and those that reproduce the cosmological evolution given by Einstein-Hilbert action plus a barotropic perfect fluid, have been obtained.

A non-local modified gravity model with an analytic function of the d'Alembert operator is considered. This model has been recently proposed as a possible way of resolving the singularities problem in cosmology. We present an exact bouncing solution, which is simpler compared to the already known one in this model in the sense it does not require an additional matter to satisfy all the gravitational equations.

One of the greatest problems of primordial inflation is that the inflationary space-time is past-incomplete. This is mainly because Einstein's GR suffers from a space-like Big Bang singularity. It has recently been shown that ghost-free, non-local higher-derivative ultra-violet modifications of Einstein's gravity may be able to resolve the cosmological Big Bang singularity via a non-singular bounce. Within the framework of such non-local cosmological models, we are going to study both sub- and super-Hubble perturbations around an inflationary trajectory which is preceded by the Big Bounce in the past, and demonstrate that the inflationary trajectory has an ultra-violet completion and that perturbations do not suffer from any pathologies.

We study cosmological perturbations in models with a single non-local scalar field originating from the string field theory description of the rolling tachyon dynamics. We construct the equation for the energy density perturbations of the non-local scalar field and explicitly prove that for the free field it is identical to a system of local cosmological perturbation equations in a particular model with multiple (maybe infinitely many) local free scalar fields.

A nonlocal gravity model with a function $f(\Box^{-1} R)$, where $\Box$ is the d'Alembert operator, is considered. The algorithm, allowing to reconstruct $f(\Box^{-1} R)$, corresponding to the given Hubble parameter and the state parameter of the matter, is proposed. Using this algorithm, we find the functions $f(\Box^{-1} R)$, corresponding to de Sitter solutions.

Cosmological models with an SU(2) Yang-Mills field are studied. For a specific model with a minimally coupled Yang-Mills Lagrangian, which includes an arbitrary function of the second-order term and a fourth-order term, a corresponding reconstruction program is proposed. It is shown that the model with minimal coupling has no de Sitter solutions, for any nontrivial function of the second-order term. To get de Sitter solutions, a gravitational model with nonminimally coupled Yang-Mills fields is then investigated. It is shown that the model with non-minimal coupling has in fact a de Sitter solution, even in absence of the cosmological constant term.

A nonlocal gravity model, which does not assume the existence of a new dimensional parameter in the action and includes a function $f(\Box^{-1} R)$, with $\Box$ the d'Alembertian operator, is considered. The model is proven to have de Sitter solutions only if the function f satisfies a certain second-order linear differential equation. The de Sitter solutions corresponding to the simplest case, an exponential function f, are explored, without any restrictions on the parameters. If the value of the Hubble parameter is positive, the de Sitter solution is stable at late times, both for negative and for positive values of the cosmological constant. Also, the stability of the solutions with zero cosmological constant is discussed and sufficient conditions for it are established in this case. New de Sitter solutions are obtained which correspond to the model with dark matter and stability is proven in this situation for nonzero values of the cosmological constant.

A non-local gravity model, which includes a function $f(\Box^{-1} R)$, where $\Box$ is the d'Alembert operator, is considered. For the model with an exponential $f(\Box^{-1} R)$ de Sitter solutions are explored, without any restrictions on the parameters. Using Hubble-normalized variables, the stability of the de Sitter solutions is investigated, with respect to perturbations in the Bianchi I metric, in the case of zero cosmological constant, and sufficient conditions for stability are obtained.

A general class of f(R) gravity models with minimally coupling a nonlocal scalar field is considered. The Ostrogradski representation for nonlocal gravitational models with a quadratic potential and the way of its localization are proposed. We study the action with an arbitrary analytic function $F(\Box)$, which has both simple and double roots. The way of localization allows to find particular solutions of nonlocal equations of gravity.

We develop the cosmological perturbations formalism in models with a single non-local scalar field originating from the string field theory description of the rolling tachyon dynamics. We construct the equation for the energy density perturbations of the non-local scalar field in the presence of the arbitrary potential and formulate the local system of equations for perturbations in the linearized model when both simple and double roots of the characteristic equation are present. We carry out the general analysis related to the curvature and entropy perturbations and consider the most specific example of perturbations when important quantities in the model become complex.

A general class of gravitational models driven by a nonlocal scalar field with a linear or quadratic potential is considered. We study the action with an arbitrary analytic function $F(\Box)$, which has both simple and double roots. The way of localization of nonlocal Einstein equations is generalized on models with linear potentials. Exact solutions in the Friedmann-Robertson-Walker and Bianchi I metrics are presented.

A method for the search of exact solutions for equation of a nonlocal scalar field in a non-flat metric is considered. In the Friedmann-Robertson-Walker metric the proposed method can be used in the case of an arbitrary potential, with the exception of linear and quadratic potentials, and allows to get in quadratures solutions, which depend on two arbitrary parameters. Exact solutions have been found for an arbitrary cubic potential, which consideration is motivated by the string field theory, as well as for exponential, logarithmic and power potentials. It has been shown that one can add the k-essence field to the model to get exact solutions for all Einstein equations.

The stability of isotropic cosmological solutions for two-field models in the Bianchi I metric is considered. We prove that the sufficient conditions for the Lyapunov stability in the Friedmann-Robertson-Walker metric provide the stability with respect to anisotropic perturbations in the Bianchi I metric and with respect to the cold dark matter energy density fluctuations. Sufficient conditions for the Lyapunov stability of the isotropic fixed points of the system of the Einstein equations have been found. We use the superpotential method to construct stable kink-type solutions and obtain sufficient conditions on the superpotential for the Lyapunov stability of the corresponding exact solutions. We analyze the stability of isotropic kink-type solutions for string field theory inspired cosmological models.

Nonlocal cosmological models with quadratic potentials are considered. We study the action with an arbitrary analytic function F(\Box_g), which has both double and simple roots. The formulae for nonlocal energy-momentum tensor, which correspond to double roots, have been obtained. The way to find particular solutions for nonlocal Einstein equations in the case when F(\Box_g) has both simple and double roots has been proposed. One and the same functions solve the initial nonlocal Einstein equations and the obtained local Einstein equations.

The stability of isotropic cosmological solutions in the Bianchi I model is considered. We prove that the stability of isotropic solutions in the Bianchi I metric for a positive Hubble parameter follows from their stability in the Friedmann-Robertson-Walker metric. This result is applied to models inspired by string field theory, which violate the null energy condition. Examples of stable isotropic solutions are presented. We also consider the k-essence model and analyse the stability of solutions of the form $\Phi(t)=t$.

A general class of cosmological models driven by a nonlocal scalar field inspired by the string field theory is studied. Using the fact that the considering linear nonlocal model is equivalent to an infinite number of local models we have found an exact special solution of the nonlocal Friedmann equations. This solution describes a monotonically increasing Universe with the phantom dark energy.

A general class of cosmological models driven by a nonlocal scalar field inspired by the string field theory is studied. Using the fact that the considering linear nonlocal model is equivalent to an infinite number of local models we have found an exact special solution of the nonlocal Friedmann equations. This solution describes a monotonically increasing Universe with the phantom dark energy.

A general class of cosmological models driven by a non-local scalar field inspired by string field theories is studied. In particular cases the scalar field is a string dilaton or a string tachyon. A distinguished feature of these models is a crossing of the phantom divide. We reveal the nature of this phenomena showing that it is caused by an equivalence of the initial non-local model to a model with an infinite number of local fields some of which are ghosts. Deformations of the model that admit exact solutions are constructed. These deformations contain locking potentials that stabilize solutions. Bouncing and accelerating solutions are presented.

A dark energy model with a phantom scalar field, an usual scalar field and the string field theory inspired polynomial potential has been constructed. A two-parameter set of exact solutions to the Friedmann equations has been found. We have constructed such string inspired potential that some exact solutions correspond to the state parameter $w_{DE}>-1$ at large time, whereas other ones correspond to $w_{DE}<-1$ at large time. We demonstrate that the superpotential method is very effective to seek new exact solutions. We also present a two-fields model with a polynomial potential and the state parameter, which crosses the cosmological barrier infinitely often.

The Painleve test is very useful to construct not only the Laurent series solutions of systems of nonlinear ordinary differential equations but also the elliptic and trigonometric ones. The standard methods for constructing the elliptic solutions consist of two independent steps: transformation of a nonlinear polynomial differential equation into a nonlinear algebraic system and a search for solutions of the obtained system. It has been demonstrated by the example of the generalized Henon-Heiles system that the use of the Laurent series solutions of the initial differential equation assists to solve the obtained algebraic system. This procedure has been automatized and generalized on some type of multivalued solutions. To find solutions of the initial differential equation in the form of the Laurent or Puiseux series we use the Painleve test. This test can also assist to solve the inverse problem: to find the form of a polynomial potential, which corresponds to the required type of solutions. We consider the five-dimensional gravitational model with a scalar field to demonstrate this.

An exact solution to the Friedmann equations with a string inspired phantom scalar matter field is constructed and the absence of the "Big Rip" singularity is shown explicitly. The notable features of the concerned model are a ghost sign of the kinetic term and a special polynomial form of the effective tachyon potential. The constructed solution is stable with respect to small fluctuations of the initial conditions and special deviations of the form of the potential.

Cosmological consequences of adding the Cold Dark Matter (CDM) to the exactly solvable stringy Dark Energy (DE) model are investigated. The model is motivated by the consideration of our Universe as a slowly decaying D3-brane. The decay of this D-brane is described in the String Field Theory framework. Stability conditions of the exact solution with respect to small fluctuations of the initial value of the CDM energy density are found. Solutions with large initial value of the CDM energy density attracted by the exact solution without CDM are constructed numerically. In contrast to the LambdaCDM model the Hubble parameter in the model is not a monotonic function of time. For specific initial data the DE state parameter w_{DE} is also not monotonic function of time. For these cases there are two separate regions of time where w_{DE} being less than -1 is close to -1.

We explore a possibility for the Universe to cross the w=-1 cosmological constant barrier for the dark energy state parameter. We consider the Universe as a slowly decaying D3-brane. The D3-brane dynamics is approximately described by a nonlocal string tachyon interaction and a back reaction of gravity is incorporated in the closed string tachyon dynamics. In a local effective approximation this model contains one phantom component and one usual field with a simple polynomial interaction. To understand cosmological properties of this system we study toy models with the same scalar fields but with modified interactions. These modifications admit polynomial superpotentials. We find restrictions on these interactions under which it is possible to reach w=-1 from below at large time. Explicit solutions with the dark energy state parameter crossing/non-crossing the barrier w=-1 at large time are presented.

The Painleve test is very useful to construct not only the Laurent-series solutions but also the elliptic and trigonometric ones. Such single-valued functions are solutions of some polynomial first order differential equations. To find the elliptic solutions we transform an initial nonlinear differential equation in a nonlinear algebraic system in parameters of the Laurent-series solutions of the initial equation. The number of unknowns in the obtained nonlinear system does not depend on number of arbitrary coefficients of the used first order equation. In this paper we describe the corresponding algorithm, which has been realized in REDUCE and Maple.

The generalized Henon-Heiles system with an additional nonpolynomial term is considered. In two nonintegrable cases new two-parameter solutions have been obtained in terms of elliptic functions. These solutions generalize the known one-parameter solutions. The singularity analysis shows that it is possible that three-parameter single-valued solutions exist in these two nonintegrable cases. The knowledge of the Laurent series solutions simplifies search of the elliptic solutions and allows to automatize it.

The generalized Henon-Heiles system has been considered. In two nonintegrable cases with the help of the Painleve test new special solutions have been found as Laurent series, depending on three parameters. The obtained series converge in some ring. One of parameters determines the singularity point location, other parameters determine coefficients of series. For some values of these parameters the obtained Laurent series coincide with the Laurent series of the known exact solutions. The Painleve test can be used not only to construct local solutions as the Laurent series but also to find elliptic solutions.

The Henon-Heiles system in the general form has been considered. In a nonintegrable case with the help of the Painleve test new solutions have been found as formal Laurent or Puiseux series, depending on three parameters. One of parameters determines a location of the singularity point, other parameters determine coefficients of series. It has been proved, that if absolute values of these two parameters are less or equal to unit, then obtained series converge in some ring. For some values of these parameters the obtained Laurent series coincide with the Laurent series of the known exact solutions.

The H\'enon--Heiles system in the general form is studied. In a nonintegrable case new solutions have been found as formal Laurent series, depending on three parameters. One of parameters determines a location of the singularity point, other parameters determine coefficients of the Laurent series. For some values of these two parameters the obtained Laurent series coincide with the Laurent series of the known exact solutions.