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 consider properties of the trajectory with the zero momentum inside a spherically symmetric black hole. We work mostly in the Painlev\'{e}% -Gullstrand frame and use the concept of the "river model of black hole". This consept allows us to decompose (in a "cosmological manner") the geodesic motion of a test particle into a "flow" of the frame and a peculiar motion with respect to this frame. After this decomposition the application of standard formulae of special relativity for kinematic processes becomes possible. The present paper expands the notion of peculiar velocities to the region under the horzion and exploits it for the description of two physical processess - high energy collisions and redshift. Using this approach we (i) present a novel description of particle collisions occuring near black hole horizons inside the event horizon. In particular, we show that the trajectory under discussion is relevant for ultra-high energy collisions. (ii) In the framework of the river model, we derive a simple formula (both outside and inside the horizon) for the redshift in the case of radial motion. It represents the product of two factors. One of them is responsible for pure gravitational part whereas the other one gives the Doppler shift due to peculiar motion with respect to the "flow".

We revisit the problem of generality of Starobinsky and Higgs inflation. The known results obtained in the Einstein frame are generalized for the case of an arbitrary initial energy of the scalar field. These results are compared with the results obtained directly in the Jordan frame. We demonstrate that the qualitative picture of initial conditions zone in the $(\phi, \dot \phi)$ plane, which leads to sufficient amount of inflation, is quite similar for both the frames in the case of Higgs inflation. For Starobinsky inflation, the conformal transformation between the frames relates the geometrical variables in the Jordan frame with the properties of an effective scalar field in the Einstein frame. We show that the transformation $(H, R) \to (\phi, \dot \phi)$ is not regular everywhere, leading to some peculiarities in the zone of good initial conditions in the $(H, R)$ plane.

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 examine the class of initial conditions which give rise to inflation. Our analysis is carried out for several popular models including: Higgs inflation, Starobinsky inflation, chaotic inflation, axion monodromy inflation and non-canonical inflation. In each case we determine the set of initial conditions which give rise to sufficient inflation, with at least $60$ e-foldings. A phase-space analysis has been performed for each of these models and the effect of the initial inflationary energy scale on inflation has been studied numerically. This paper discusses two scenarios of Higgs inflation: (i) the Higgs is coupled to the scalar curvature, (ii) the Higgs Lagrangian contains a non-canonical kinetic term. In both cases we find Higgs inflation to be very robust since it can arise for a large class of initial conditions. One of the central results of our analysis is that, for plateau-like potentials associated with the Higgs and Starobinsky models, inflation can be realised even for initial scalar field values which lie close to the minimum of the potential. This dispels a misconception relating to plateau potentials prevailing in the literature. We also find that inflation in all models is more robust for larger values of the initial energy scale.

In this paper we investigate the past evolution of an anisotropic Bianchi I universe in $R+R^2$ gravity. Using the dynamical system approach we show that there exists a new two-parameters set of solutions that includes both an isotropic "false radiation" solution and an anisotropic generalized Kasner solution, which is stable. We derive the analytic behaviour of the shear from a specific property of $f(R)$ gravity and the analytic asymptotic form of the Ricci scalar when approaching the initial singularity. Finally we numerically checked our results.

In this paper we address two important issues which could affect reaching the exponential and Kasner asymptotes in Einstein-Gauss-Bonnet cosmologies -- spatial curvature and anisotropy in both three- and extra-dimensional subspaces. In the first part of the paper we consider cosmological evolution of spaces being the product of two isotropic and spatially curved subspaces. It is demonstrated that the dynamics in $D=2$ (the number of extra dimensions) and $D \geqslant 3$ is different. It was already known that for the $\Lambda$-term case there is a regime with "stabilization" of extra dimensions, where the expansion rate of the three-dimensional subspace as well as the scale factor (the "size") associated with extra dimensions reach constant value. This regime is achieved if the curvature of the extra dimensions is negative. We demonstrate that it take place only if the number of extra dimensions is $D \geqslant 3$. In the second part of the paper we study the influence of initial anisotropy. Our study reveals that the transition from Gauss-Bonnet Kasner regime to anisotropic exponential expansion (with expanding three and contracting extra dimensions) is stable with respect to breaking the symmetry within both three- and extra-dimensional subspaces. However, the details of the dynamics in $D=2$ and $D \geqslant 3$ are different. Combining the two described affects allows us to construct a scenario in $D \geqslant 3$, where isotropisation of outer and inner subspaces is reached dynamically from rather general anisotropic initial conditions.

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 consider dynamics of a flat anisotropic Universe filled by a perfect fluid near a cosmological singularity in quadratic gravity. Two possible regimes are described -- the Kasner anisotropic solution and an isotropic "vacuum radiation" solution which has three sub cases depending on whether the equation of state parameter $w$ is bigger, smaller or equals to $1/3$. Initial conditions for numerical integrations have been chosen near General Relativity anisotropic solution with matter (Jacobs solution). We have found that for such initial conditions there is a range of values of coupling constants so that the resulting cosmological singularity is isotropic.

In this paper cosmological dynamics in Einstein-Gauss-Bonnet gravity with a perfect fluid source in arbitrary dimension is studied. A systematic analysis is performed for the case that the theory does not admit maximally symmetric solutions. Considering two independent scale factors, namely one for the three dimensional space and one for the extra dimensional space, is found that a regime exists where the two scale factors tend to a constant value via damped oscillations for not too negative pressure of the fluid, so that asymptotically the evolution of the $(3+1)$-dimensional Friedmann model with perfect fluid is recovered. At last, it is worth emphasizing that the present numerical results strongly support a 't Hooft-like interpretation of the parameter $1/D$ (where $D$ is the number of extra dimensions) as a small expansion parameter in very much the same way as it happens in the large $N$ expansion of gauge theories with $1/N$. Indeed, the dependence on $D$ of many of the relevant physical quantities computed here manifests a clear WKB-like pattern, as expected on the basis of large $N$ arguments.

We study flat Friedmann-Lemaitre-Robertson-Walker cosmological models for a scalar field coupled nonminimally to teleparallel gravity with generic coupling and potential functions. The goal of this paper is to determine the conditions under which the cosmological evolution tends to the limit where the variation of the gravitational "constant" ceases and the system evolves close to general relativity. These conditions can be read off from the approximate analytical solutions describing the process in matter and potential domination eras. Only those models where the GR limit exists and is an attractor can be considered viable. We expect the results to hold in the original "pure tetrad" formulation as well as in the recently suggested covariant formulation of the teleparallel theory. In the former case the GR attractor simultaneously provides a mechanism how cosmological evolution suppresses the problematic degrees of freedom stemming from the lack of local Lorentz invariance.

It is generally believed that a cosmological arrow of time must be associated with entropy production. Indeed, in his seminal work on cyclic cosmology, Tolman introduced a viscous fluid in order to make successive expansion/contraction cycles larger than previous ones, thereby generating an arrow of time. However, as we demonstrate in this letter, the production of entropy is not the only means by which a cosmological arrow of time may emerge. Remarkably, systems which are dissipationless may nevertheless demonstrate a preferred direction of time provided they possess attractors. An example of a system with well defined attractors is scalar-field driven cosmology. In this case, for a wide class of potentials (especially those responsible for inflation), the attractor equation of state during expansion can have the form $p \simeq -\rho$, and during contraction $p \simeq \rho$. If the resulting cosmology is cyclic, then the presence of cosmological hysteresis, $\oint p~dV \neq 0$ during successive cycles, causes an arrow of time to emerge in a system which is formally dissipationless. An important analogy is drawn between the arrow of time in cyclic cosmology and an arrow of time in an $N$-body system of gravitationally interacting particles. We find that, like the $N$-body system, a cyclic universe can evolve from a single past into two futures with oppositely directed arrows of time.

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.

In this paper we investigate the constant volume exponential solutions (i.e. the solutions with the scale factors change exponentially over time so that the comoving volume remains the same) in the Einstein-Gauss-Bonnet gravity. We find conditions for these solutions to exist and show that they are compatible with any perfect fluid with the equation of state parameter $\omega<1/3$ if the matter density of the Universe exceeds some critical value. We write down some exact solutions which generalize ones found in our previous paper for models with a cosmological constant.

We study cosmological solutions in $R + \beta R^{N}$-gravity for an isotropic Universe filled with ordinary matter with the equation of state parameter $\gamma$. Using the Bogolyubov-Krylov-Mitropol'skii averaging method we find asymptotic oscillatory solutions in terms of new functions, which have been specially introduced by us for this problem and appeared as a natural generalization of the usual sine and cosine. It is shown that the late-time behaviour of the Universe in the model under investigation is determined by the sign of the difference $\gamma-\gamma_{crit}$ where $\gamma_{crit}=2N/(3N-2)$. If $\gamma < \gamma_{crit}$, the Universe reaches the regime of small oscillations near values of Hubble parameter and matter density, corresponding to General Relativity solution. Otherwise higher-curvature corrections become important at late times. We also study numerically basins of attraction for the oscillatory and phantom solutions, which are present in the theory for $N>2$. Some important differences between $N=2$ and $N>2$ cases are discussed.

In this methodological note we discuss some details and peculiarities of the cosmic expansion as viewed by a realistic observer. We show that the velocity $v_\Theta$ related to a change (measured by observer's clock) of the angular distance, plays an important role in formation of a meaningful observed picture of the expansion of the universe. Usage of this velocity and the angular distance (in addition to the standard approach --- proper distance and corresponding velocity) allows to present the cosmic expansion in a more illustrative manner. These parameters play a key role in visualization of the expansion.

We investigate the cosmological dynamics of non-minimally coupled scalar field system described by $F(\phi)R$ coupling with $F(\phi)=(1-\xi\phi^N)R$($N\ge2$) and the field potential, $V(\phi)=V_0\phi^n$. We use a generic set of dynamical variables to bring out new asymptotic regimes of the underlying dynamics. However, our dynamical variables miss the most important fixed point$-$ the de Sitter solution. We make use of the original form of system of equations to investigate the issues related to this important solution. In particular, we show that the de-Sitter solution which is a dynamical attractor of the system lies in the region of negative effective gravitational constant $G_N$ thereby leading to a ghost dominated universe in future and a transient quintessence(phantom) phase with $G_N>0 $ around the present epoch (however, as demonstrated by Starobinsky in 1981, the ghost dominated universe, if exists, can not be accessed from the Universe we live in, we shall say more about this important result in the last section). We also carry out comparison of the model with other competing models of dark energy such as galileon modified gravity and others.

This paper revisits the Inflationary scenario within the framework of scalar field models possessing a non-canonical kinetic term. We obtain closed form solutions for all essential quantities associated with chaotic inflation including slow roll parameters, scalar and tensor power spectra, spectral indices, the tensor-to-scalar ratio, etc. We also examine the Hamilton-Jacobi equation and demonstrate the existence of an inflationary attractor. Our results highlight the fact that non-canonical scalars can significantly improve the viability of inflationary models. They accomplish this by decreasing the tensor-to-scalar ratio while simultaneously increasing the value of the scalar spectral index, thereby redeeming models which are incompatible with the cosmic microwave background (CMB) in their canonical version. For instance, the non-canonical version of the chaotic inflationary potential, $V(\phi) \sim \lambda\phi^4$, is found to agree with observations for values of $\lambda$ as large as unity ! The exponential potential can also provide a reasonable fit to CMB observations. A central result of this paper is that {\em steep potentials} (such as $V \propto \phi^{-n}$) usually associated with dark energy, can drive inflation in the non-canonical setting. Interestingly, non-canonical scalars violate the consistency relation $r = -8n_T$, which emerges as a {\em smoking gun} test for this class of models.

A Universe filled with a homogeneous scalar field exhibits `Cosmological hysteresis'. Cosmological hysteresis is caused by the asymmetry in the equation of state during expansion and contraction. This asymmetry results in the formation of a hysteresis loop: $\oint pdV$, whose value can be non-vanishing during each oscillatory cycle. For flat potentials, a negative value of the hysteresis loop leads to the increase in amplitude of consecutive cycles and to a universe with older and larger successive cycles. Such a universe appears to possess an arrow of time even though entropy production is absent and all of the equations respect time-reversal symmetry ! Cosmological hysteresis appears to be widespread and exists for a large class of scalar field potentials and mechanisms for making the universe bounce. For steep potentials, the value of the hysteresis loop can be positive as well as negative. The expansion factor in this case displays quasi-periodic behaviour in which successive cycles can be both larger as well as smaller than previous ones. This quasi-regular pattern resembles the phenomenon of BEATS displayed by acoustic systems. Remarkably, the expression relating the increase/decrease in oscillatory cycles to the quantum of hysteresis appears to be model independent. The cyclic scenario is extended to spatially anisotropic models and it is shown that the anisotropy density decreases during successive cycles if the hysteresis loop is negative.

We have investigated the decay statistics of triple systems with different masses in Newtonial dynamics. We demonstrate that in a broad interval of mass ratios this statistics has good approximation by power-law tails. The power indices do not show any significant dependence on mass ratios.

We discuss the cosmological reconstruction in modified Gauss-Bonnet and F(R) gravities. Two alternative representations of the action (with and without auxiliary scalar) are considered. The approximate description of deceleration-acceleration transition cosmologies is reconstructed. It is shown that cosmological solution containing Big Bang and Big Rip singularities may be reconstructed only using the representation with the auxiliary field. The analytical description of the deceleration-acceleration transition cosmology in modified Gauss-Bonnet gravity is demonstrated to be impossible at sufficiently general conditions.

We investigate the cosmological properties of an "induced gravity" brane scenario in the absence of mirror symmetry with respect to the brane. We find that brane evolution can proceed along one of four distinct branches. By contrast, when mirror symmetry is imposed, only two branches exist, one of which represents the self-accelerating brane, while the other is the so-called normal branch. This model incorporates many of the well-known possibilities of brane cosmology including phantom acceleration (w < -1), self-acceleration, transient acceleration, quiescent singularities, and cosmic mimicry. Significantly, the absence of mirror symmetry also provides an interesting way of inducing a sufficiently small cosmological constant on the brane. A small (positive) Lambda-term in this case is induced by a small asymmetry in the values of bulk fundamental constants on the two sides of the brane.

In this paper, we analyze the viability of a vacuum Gauss-Bonnet cosmology by examining the dynamics of the homogeneous and anisotropic background in 4+1 dimensions. The trajectories of the system either originate from the standard singularity or from non-standard type, the later is characterized by the divergence of time derivative of the Hubble parameters for its finite value. At the onset, the system should relax to Einstein phase at late times as the effect of Gauss-Bonnet term becomes negligible in the low energy regime. However, we find that most of the trajectories emerging from the standard big-bang singularity lead to future re-collapse whereas the system beginning its evolution from the non-standard singularity enters the Kasner regime at late times. This leads to the conclusion that the measure of trajectories giving rise to a smooth evolution from a standard singularity to the Einstein phase is negligibly small for generic initial conditions.

We examine a class of braneworld models in which the expanding universe encounters a "quiescent" future singularity. At a quiescent singularity, the energy density and pressure of the cosmic fluid as well as the Hubble parameter remain finite while all derivatives of the Hubble parameter diverge (i.e., ${\dot H}$, ${\ddot H}$, etc. $\to \infty$). Since the Kretschmann invariant diverges ($R_{iklm}R^{iklm} \to \infty$) at the singularity, one expects quantum effects to play an important role as the quiescent singularity is approached. We explore the effects of vacuum polarization due to massless conformally coupled fields near the singularity and show that these can either cause the universe to recollapse or, else, lead to a softer singularity at which $H$, ${\dot H}$, and ${\ddot H}$ remain finite while ${\dddot H}$ and higher derivatives of the Hubble parameter diverge. An important aspect of the quiescent singularity is that it is encountered in regions of low density, which has obvious implications for a universe consisting of a cosmic web of high and low density regions -- superclusters and voids. In addition to vacuum polarization, the effects of quantum particle production of non-conformal fields are also likely to be important. A preliminary examination shows that intense particle production can lead to an accelerating universe whose Hubble parameter shows oscillations about a constant value.

We study the evolution of (phantom) dark energy universe by taking into account the higher-order string corrections to Einstein-Hilbert action with a fixed dilaton. While the presence of a cosmological constant gives stable de-Sitter fixed points in the cases of heterotic and bosonic strings, no stable de-Sitter solutions exist when a phantom fluid is present. We find that the universe can exhibit a Big Crunch singularity with a finite time for type II string, whereas it reaches a Big Rip singularity for heterotic and bosonic strings. Thus the fate of dark energy universe crucially depends upon the type of string theory under consideration.

The Einstein-Gauss-Bonnet equations projected from the bulk to brane lead to a complicated Friedmann equation which simplifies to $H^2 \sim \rho^q$ in the asymptotic regimes. The Randall-Sundrum (RS) scenario corresponds to $q=2$ whereas $q=2/3$ $&$ $q=1$ give rise to high energy Gauss-Bonnet (GB) regime and the standard GR respectively. Amazingly, while evolving from RS regime to high energy GB limit, one passes through a GR like region which has important implications for brane world inflation. For tachyon GB inflation with potentials $V(\phi) \sim \phi^p$ investigated in this paper, the scalar to tensor ratio of perturbations $R$ is maximum around the RS region and is generally suppressed in the high energy regime for the positive values of $p$. The ratio is very low for $p>0$ at all energy scales relative to GB inflation with ordinary scalar field. The models based upon tachyon inflation with polynomial type of potentials with generic positive values of $p$ turn out to be in the $1 \sigma$ observational contour bound at all energy scales varying from GR to high energy GB limit. The spectral index $n_S$ improves for the lower values of $p$ and approaches its scale invariant limit for $p=-2$ in the high energy GB regime. The ratio $R$ also remains small for large negative values of $p$, however, difference arises for models close to scale invariance limit. In this case, the tensor to scale ratio is large in the GB regime whereas it is suppressed in the intermediate region between RS and GB. Within the frame work of patch cosmologies governed by $H^2 \sim \rho^q$, the behavior of ordinary scalar field near cosmological singularity and the nature of scaling solutions are distinguished for the values of $q < 1$ and $q > 1$.

It is well known that a closed universe with a minimally coupled massive scalar field always collapses to a singularity unless the initial conditions are extremely fine tuned. We show that the corrections to the equations of motion for the massive scalar field, given by loop quantum gravity in high curvature regime, always lead to a bounce independently of the initial conditions. In contrast to the previous works in loop quantum cosmology, we note that the singularity can be avoided even at the semi-classical level of effective dynamical equations with non-perturbative quantum gravity modifications, without using a discrete quantum evolution.

In this paper we analyze the cosmological dynamics of phantom field in a variety of potentials unbounded from above. We demonstrate that the nature of future evolution generically depends upon the steepness of the phantom potential and discuss the fate of Universe accordingly.

We study nature of singularities in anisotropic string-inspired cosmological models in the presence of a Gauss-Bonnet term. We analyze two string gravity models-- dilaton-driven and modulus-driven cases-- in the Bianchi type-I background without an axion field. In both scenarios singularities can be classified in two ways- the determinant singularity where the main determinant of the system vanishes and the ordinary singularity where at least one of the anisotropic expansion rates of the Universe diverges. In the dilaton case, either of these singularities inevitably appears during the evolution of the system. In the modulus case, nonsingular cosmological solutions exist both in asymptotic past and future with determinant $D=+\infty$ and D=2, respectively. In both scenarios nonsingular trajectories in either future or past typically meet the determinant singularity in past/future when the solutions are singular, apart from the exceptional case where the sign of the time-derivative of dilaton is negative. This implies that the determinant singularity may play a crucial role to lead to singular solutions in an anisotropic background.

The cosmic microwave background (CMB) polarization, whose experimental detection must become the next natural step in the study of the CMB radiation, depends on the dynamics of hydrogen recombination in the Universe, on the relation between the power of scalar perturbations and gravitational waves, and on the presence or absence of secondary ionization. The spectrum of the CMB anisotropy is currently known with an accuracy that is high enough to draw some conclusions about the cosmological parameters and about the presence of secondary ionization of matter. The reduction of observational data strongly suggests the presence of secondary ionization. We consider the effect of this ionization at redshifts $z<100$ on the CMB polarization generated by the scalar mode on various angular scales.