We investigated realistic neutron stars in axion $R^{2}$ gravity. The coupling between curvature and axion field $\phi$ is assumed in the simple form $\sim R^2\phi$. For the axion mass in the range $m_{a}\sim 10^{-11}-10^{-10}$ eV the solitonic core within neutron star and corresponding halo with size $\sim 100$ km can exist. Therefore the effective contribution of $R^2$ term grows inside the star and it leads to change of star parameters (namely, mass and radius). We obtained the increase of star mass independent from central density for wide range of masses. Therefore, maximal possible mass for given equation of state grows. At the same time, the star radius increases not so considerably in comparison with GR. {Hence, our model may predict possible existence of supermassive compact stars with masses $M\sim 2.2-2.3M_\odot$ and radii $R_{s}\sim 11$ km for realistic equation of state (we considered APR equation of state). In General Relativity one can obtain neutron stars with such characteristics only for unrealistic, extremely stiff equations of state.} Note that this increase of mass occurs due to change of solution for scalar curvature outside the star. In GR curvature drops to zero on star surface where $\rho=p=0$. In the model under consideration the scalar curvature dumps more slowly in comparison with vacuum $R^2$ gravity due to axion "galo" around the star.

The article is dedicated to one of the most undeservedly overlooked properties of the cosmological models: the behaviour at, near and due to a jump discontinuity. It is most interesting that while the usual considerations of the cosmological dynamics deals heavily in the singularities produced by the discontinuities of the second kind (a.k.a. the essential discontinuities) of one (or more) of the physical parameters, almost no research exists to date that would turn to their natural extension/counterpart: the singularities induced by the discontinuities of the first kind (a.k.a. the jump discontinuities). It is this oversight that this article aims to amend. In fact, it demonstrates that the inclusion of such singularities allows one to produce a number of very interesting scenarios of cosmological evolution. For example, it produces the cosmological models with a finite value of the equation of state parameter $w=p/\rho$ even when both the energy density and the pressure diverge, while at the same time keeping the scale factor finite. Such a dynamics is shown to be possible only when the scale factor experiences a finite jump at some moment of time. Furthermore, if it is the first derivative of the scale factor that experiences a jump, then a whole new and different type of a sudden future singularity appears. Finally, jump discontinuities suffered by either a second or third derivatives of a scale factor lead to cosmological models experiencing a sudden dephantomization -- or avoiding the phantomization altogether. This implies that theoretically there should not be any obstacles for extending the cosmological evolution beyond the corresponding singularities; therefore, such singularities can be considered a sort of a cosmological phase transition.

Neutron stars are extremely strong cosmic magnets which fields are expected to decay with time. Here we report on the simple test of this process. Adopting a novel approach, we have estimated surface magnetic fields $B$ for 76 radiopulsars (the most numerous subclass of the known isolated neutron stars) which ages $t$ are known independently. Focusing on the accurate evaluation of the precision of both quantities, we determined a significant power-law trend $B(t) \propto t^{-\beta}$ with index $\beta = 0.19^{+0.05}_{-0.06}$ at 95\% C.L. The effects of the observational selection turn this value into the upper limit for the intrinsic field decay rate. If so, then neutron star crusts are close to the "impurity-free crystals", which results in a relatively slow magnetic fields decay.

In the context of f(R)=R + alpha R^2 gravity, we study the existence of neutron and quark stars with no intermediate approximations in the generalised system of Tolman-Oppenheimer-Volkov equations. Analysis shows that for positive alpha's the scalar curvature does not drop to zero at the star surface (as in General Relativity) but exponentially decreases with distance. Also the stellar mass bounded by star surface decreases when the value alpha increases. Nonetheless distant observers would observe a gravitational mass due to appearance of a so-called gravitational sphere around the star. The non-zero curvature contribution to the gravitational mass eventually is shown to compensate the stellar mass decrease for growing alpha's. We perform our analysis for several equations of state including purely hadronic configurations as well as hyperons and quark stars. In all cases, we assess that the relation between the parameter $\alpha$ and the gravitational mass weakly depend upon the chosen equation of state. Another interesting feature is the increase of the star radius in comparison to General Relativity for stars with masses close to maximal, whereas for intermediate masses around 1.4-1.6 solar masses, the radius of star depends upon alpha very weakly. Also the decrease in the mass bounded by star surface may cause the surface redshift to decrease in R^2-gravity when compared to Einsteinian predictions. This effect is shown to hardly depend upon the observed gravitational mass. Finally, for negative values of alpha our analysis shows that outside the star the scalar curvature has damped oscillations but the contribution of the gravitational sphere into the gravitational mass increases indefinitely with radial distance putting into question the very existence of such relativistic stars.

Numerical simulations of realistic non-vacuum magnetospheres of isolated neutron stars have shown that pulsar spin-down luminosities depend weakly on the magnetic obliquity. This result provides the opportunity to estimate the surface magnetic field for a given radiopulsar with quite a high accuracy. In the current work, we present a refinement of the classical `magneto-dipolar' formula for pulsar magnetic fields $B_{\rm md} = (3.2\times 10^{19}\mbox{ G})\sqrt{P\dot P}$, where $P$ is the neutron star spin period. The new, robust timing-based estimator is introduced as $\log B = \log B_{\rm md} + \Delta_{\rm B}(M, \alpha)$, where the correction $\Delta_{\rm B}$ depends on the equation of state (EOS) of dense matter, the individual pulsar obliquity $\alpha$ and the mass $M$. Adopting state-of-the-art statistics for $M$ and $\alpha$ we calculate the distributions of $\Delta_{\rm B}$ for a representative subset of 22 EOSs that do not contradict observations. It has been found that $\Delta_{\rm B}$ is distributed nearly normally, with the average in the range -0.5 to -0.25 dex and standard deviation $\sigma[\Delta_{\rm B}] \approx$ 0.06 to 0.09 dex, depending on the adopted EOS. The latter quantity represents a formal uncertainty of the corrected estimation of $\log B$ because $\Delta_{\rm B}$ is weakly correlated with $\log B_{\rm md}$. At the same time, if it is assumed that every considered EOS has the same chance of occurring in nature, then another, more generalized, estimator $B^* \approx 3B_{\rm md}/7$ can be introduced providing an unbiased value of the pulsar surface magnetic field with $\sim$30 per cent uncertainty with 68 per cent confidence. Finally, we discuss the possible impact of pulsar timing irregularities on the timing-based estimation of $B$ and review the astrophysical applications of the obtained results.

In this paper we propose and extensively study mimetic $f({\cal G})$ modified gravity models, with various scenarios of cosmological evolution, with or without extra matter fluids. The easiest formulation is based on the use of Lagrange multiplier constraint. In certain versions of this theory, it is possible to realize accelerated expansion of the Universe or even unified evolution which includes inflation with dark energy, and at the same time in the same theoretical framework, dark matter is described by the theory. This is achieved by the re-parametrization of the metric tensor, which introduces a new degree of freedom in the cosmological equations and leads to appearance of the mimetic "dark matter" component. In the context of mimetic $f({\cal G})$ theory, we also provide some quite general reconstruction schemes, which enable us to find which $f({\cal G})$ gravity generates a specific cosmological evolution. In addition, we also provide the general reconstruction technique for Lagrange multiplier $f({\cal G})$ gravity. All our results are accompanied by illustrative examples, with special emphasis on bouncing cosmologies.

Quark star models with realistic equation of state in nonperturbative $f(R)$ gravity are considered. The mass-radius relation for $f(R)=R+\alpha R^2$ model is obtained. Considering scalar curvature $R$ as an independent function, one can find out, for each value of central density, the unique value of central curvature for which one has solutions with the required asymptotic $R\rightarrow 0$ for $r\rightarrow\infty$. In another words, one needs a fine-tuning for $R$ to achieve quark stars in $f(R)$ gravity. We consider also the analogue description in corresponding scalar-tensor gravity. The fine-tuning on $R$ is equivalent to the fine-tuning on the scalar field $\phi$ in this description. For distant observers, the gravitational mass of the star increases with increasing $\alpha$ ($\alpha>0$) but the interpretation of this fact depends on frame where we work. Considering directly $f(R)$ gravity, one can say that increasing of mass occurs by the "gravitational sphere" outside the star with some "effective mass". On the other hand, in scalar-tensor theory, we also have a dilaton sphere (or "disphere") outside the star but its contribution to gravitational mass for distant observer is negligible. We show that it is possible to discriminate modified theories of gravity from General Relativity due to the gravitational redshift of the thermal spectrum emerging from the surface of the star.

We discuss neutron stars with strong magnetic mean fields in the framework of Extended Theories of Gravity. In particular, we take into account models derived from $f(R)$ and $f(\cal G)$ extensions of General Relativity where functions of the Ricci curvature invariant $R$ and the Gauss-Bonnet invariant ${\cal G}$ are respectively considered. Dense matter in magnetic mean field, generated by magnetic properties of particles, is described by assuming a model with three meson fields and baryons octet. As result, the considerable increasing of maximal mass of neutron stars can be achieved by cubic corrections in $f(R)$ gravity. In principle, massive stars with $M> 4 M_{\odot}$ can be obtained. On the other hand, stable stars with high strangeness fraction (with central densities $\rho_{c}\sim 1.5-2.0$ GeV/fm$^{3}$) are possible considering quadratic corrections of $f(\cal {G})$ gravity. The magnetic field strength in the star center is of order $6-8\times 10^{18}$ G. In general, we can say that other branches of massive neutron stars are possible considering the extra pressure contributions coming from gravity extensions. Such a feature can constitute both a probe for alternative theories and a way out to address anomalous self-gravitating compact systems.

We consider the well-known Boltzmann brains problem in frames of simple phantom energy models with little rip and big rip singularity. It is showed that these models (i) satisfy to observational data and (ii) may be free from Boltzmann brains problem. The human observers in phantom models can exist only in during for a certain period $t<t_{f}$ ($t_{f}$ is lifetime of universe) via Bekenstein bound. If fraction of unordered observers in this part of universe history is negligible in comparison with ordered observers than Boltzmann brains problem doesn't appear. The bounds on model parameters derived from such requirement don't contradict to allowable range from observational data.

Neutron stars with strong magnetic fields are considered in the framework of f(R) gravity. In order to describe dense matter in magnetic field, the model with baryon octet interacting through $\sigma$$\rho$$\omega$-fields is used. The hyperonization process results in softening the equation of state (EoS) and in decreasing the maximal mass. We investigate the effect of strong magnetic field in models involving quadratic and cubic corrections in the Ricci scalar $R$ to the Hilbert-Einstein action. For large fields, the Mass-Radius relation differs considerably from that of General Relativity only for stars with masses close to the maximal one. Another interesting feature is the possible existence of more compact stable stars with extremely large magnetic fields ($\sim 6\times 10^{18}$ G instead of $\sim 4\times 10^{18}$ G as in General Relativity) in the central regions of the stars. Due to cubic terms, a significant increasing of the maximal mass is possible.

Neutron star models in perturbative $f(R)$ gravity are considered with realistic equations of state. In particular, we consider the FPS, SLy and other equations of state and a case of piecewise equation of state for stars with quark cores. The mass-radius relations for $f(R)=R+R(e^{-R/R_{0}}-1)$ model and for $R^2$ models with logarithmic and cubic corrections are obtained. In the case of $R^2$ gravity with cubic corrections, we obtain that at high central densities ($\rho>10\rho_{ns}$, where $\rho_{ns}=2.7\times 10^{14}$ g/cm$^{3}$ is the nuclear saturation density), stable star configurations exist. The minimal radius of such stars is close to $9$ km with maximal mass $\sim 1.9 M_{\odot}$ (SLy equation). A similar situation takes place for AP4 and BSK20 EoS. Such an effect can give rise to more compact stars than in General Relativity. If observationally identified, such objects could constitute a formidable signature for modified gravity at astrophysical level. Another interesting result can be achieved in modified gravity with only a cubic correction. For some EoS, the upper limit of neutron star mass increases and therefore these EoS can describe realistic star configurations (although, in General Relativity, these EoS are excluded by observational constraints).

We develop a method for constructing exact cosmological solutions in brane world cosmology. New classes of cosmological solutions on Randall-Sandrum brane are obtained. The superpotential and Hubble parameter are represented in quadratures. These solutions have inflationary phases under general assumptions and also describe an exit from the inflationary phase without a fine tuning of the parameters. Another class solutions can describe the current phase of accelerated expansion with or without possible exit from it.

Several dark energy models on the brane are investigated. They are compared with corresponding theories in the frame of 4d Friedmann-Robertson-Walker cosmology. To constrain the parameters of the models considered, recent observational data, including SNIa apparent magnitude measurements, baryon acoustic oscillation results, Hubble parameter evolution data and matter density perturbations are used. Explicit formulas of the so-called {\it state-finder} parameters in teleparallel theories are obtained that could be useful to test these models and to establish a link between Loop Quantum Cosmology and Brane Cosmology. It is concluded that a joint analysis as the one developed here allows to estimate, in a very convenient way, possible deviation of the real universe cosmology from the standard Friedmann-Robertson-Walker one.

We confront dark energy models which are currently similar to $\Lambda$CDM theory with observational data which include the SNe data, matter density perturbations and baryon acoustic oscillations data. DE cosmology under consideration may evolve to Big Rip, type II or type III future singularity, or to Little Rip or Pseudo-Rip universe. It is shown that matter perturbations data define more precisely the possible deviation from $\Lambda$CDM model than consideration of SNe data only. The combined data analysis proves that DE models under consideration are as consistent as $\Lambda$CDM model. We demonstrate that growth of matter density perturbations may occur at sufficiently small background density but still before the possible disintegration of bound objects (like clusters of galaxies, galaxies, etc) in Big Rip, type III singularity, Little Rip or Pseudo-Rip universe. This new effect may bring the future universe to chaotic state well before disintegration or Rip.

It is observed that one of Einstein-Friedmann's equations has formally the aspect of a Sturm-Liouville problem, and that the cosmological constant, $\Lambda$, plays thereby the role of spectral parameter (what hints to its connection with the Casimir effect). The subsequent formulation of appropriate boundary conditions leads to a set of admissible values for $\Lambda$, considered as eigenvalues of the corresponding linear operator. Simplest boundary conditions are assumed, namely that the eigenfunctions belong to $L^2$ space, with the result that, when all energy conditions are satisfied, they yield a discrete spectrum for $\Lambda>0$ and a continuous one for $\Lambda<0$. A very interesting situation is seen to occur when the discrete spectrum contains only one point: then, there is the possibility to obtain appropriate cosmological conditions without invoking the anthropic principle. This possibility is shown to be realized in cyclic cosmological models, provided the potential of the matter field is similar to the potential of the scalar field. The dynamics of the universe in this case contains a sudden future singularity.

Brane dark energy cosmologies, leading to various possible evolutions of our universe, are investigated. The discussion shows that while all these models can be made arbitrarily close to the standard $\Lambda$CDM cosmology at present, their future evolutions can be very different, even diverge with time in a number of ways. This includes asymptotic de-Sitter evolution, Little Rip with dissolution of bound structures, and various possible singularities, as the Big Rip, a sudden future singularity (Type II), and Type III and Type IV cases. Specifically, some interesting effects coming from the brane tension are investigated. It is shown, in particular, that the Little Rip occurs faster on the brane model than in usual FRW cosmology. And in the asymptotic de-Sitter regime the influence of the brane tension leads to a deviation of the value of the effective cosmological constant from that corresponding to ordinary dark energy. As a consequence, the value of the inertial force from the accelerating expansion can greatly exceed the corresponding inertial force in ordinary cosmological models.

Dark energy models with various scenarios of evolution are considered from the viewpoint of the formalism for the equation of state. It is shown that these models are compatible with current astronomical data. Some of the models presented here evolve arbitrarily close to $\Lambda$CDM up to the present, but diverge in the future into a number of different possible asymptotic states, including asymptotic de-Sitter (pseudo-rip) evolution, little rips with disintegration of bound structures, and various forms of finite-time future singularities. Therefore it is impossible from observational data to determine whether the universe will end in a future singularity or not. We demonstrate that the models under consideration are stable for a long period of time (billions of years) before entering a Little Rip/Pseudo-Rip induced dissolution of bound structures or before entering a soft finite-time future singularity. Finally, the physical consequences of Little Rip, Type II, III and Big Crush singularities are briefly compared.

In this paper, we give a possible solution to the cosmological constant problem. It is shown hat the traditional approach, based on volume weighting of probabilities, leads to an incoherent onclusion: the probability that a randomly chosen observer measures $\Lambda=0$ is exactly equal to 1. Using an alternative, volume averaging measure, instead of volume weighting can explain why the osmological constant is non-zero.

An anthropic principle has made it possible to answer the difficult question of why the observable value of cosmological constant ($\Lambda\sim 10^{-47}$ GeV${}^4$) is so disconcertingly tiny compared to predicted value of vacuum energy density $\rho_{SUSY}\sim 10^{12}$ GeV${}^4$. Unfortunately, there is a darker side to this argument, as it consequently leads to another absurd prediction: that the probability to observe the value $\Lambda=0$ for randomly selected observer exactly equals to 1. We'll call this controversy an infrared divergence problem. It is shown that the IRD prediction can be avoided with the help of a Linde-Vanchurin {\em singular runaway measure} coupled with the calculation of relative Bayesian probabilities by the means of the {\em doomsday argument}. Moreover, it is shown that while the IRD problem occurs for the {\em prediction stage} of value of $\Lambda$, it disappears at the {\em explanatory stage} when $\Lambda$ has already been measured by the observer.

We construct phantom energy models with the equation-of-state parameter $w$ such that $w<-1$, but finite-time future singularity does not occur. Such models can be divided into two classes: (i) energy density increases with time ("phantom energy" without "Big Rip" singularity) and (ii) energy density tends to constant value with time ("cosmological constant" with asymptotically de Sitter evolution). The disintegration of bound structure is confirmed in Little Rip cosmology. Surprisingly, we find that such disintegration (on example of Sun-Earth system) may occur even in asymptotically de Sitter phantom universe consistent with observational data. We also demonstrate that non-singular phantom models admit wormhole solutions as well as possibility of big trip via wormholes.

We develop a method for constructing exact cosmological solutions of the Einstein equations based on representing them as a second-order linear differential equation. In particular, the method allows using an arbitrary known solution to construct a more general solution parameterized by a set of 3\textit{N} constants, where \textit{N} is an arbitrary natural number. The large number of free parameters may prove useful for constructing a theoretical model that agrees satisfactorily with the results of astronomical observations. Cosmological solutions on the Randall-Sundrum brane have similar properties. We show that three-parameter solutions in the general case already exhibit inflationary regimes. In contrast to previously studied two-parameter solutions, these three-parameter solutions can describe an exit from inflation without a fine tuning of the parameters and also several consecutive inflationary regimes.

Recently it was found that dark energy in the form of phantom generalized Chaplygin gas may lead to a new form of the cosmic doomsday, the big freeze singularity. Like the big rip singularity, the big freeze singularity would also take place at a finite future cosmic time, but unlike the big rip singularity it happens for a finite scale factor.Our goal is to test if a universe filled with phantom generalized Chaplygin gas can conform to the data of astronomical observations. We shall see that if the universe is only filled with generalized phantom Chaplygin gas with equation of state $p=-c^2s^2/\rho^{\alpha}$ with $\alpha<-1$, then such a model cannot be matched to the data of astronomical observations. To construct matched models one actually need to add dark matter. This procedure results in cosmological scenarios which do not contradict the data of astronomical observations and allows one to estimate how long we are now from the future big freeze doomsday.