Observations of the redshift z=7.085 quasar J1120+0641 have been used to search for variations of the fine structure constant, alpha, over the redshift range 5.5 to 7.1. Observations at z=7.1 probe the physics of the universe when it was only 0.8 billion years old. These are the most distant direct measurements of alpha to date and the first measurements made with a near-IR spectrograph. A new AI analysis method has been employed. Four measurements from the X-SHOOTER spectrograph on the European Southern Observatory's Very Large Telescope (VLT) directly constrain any changes in alpha relative to the value measured on Earth (alpha_0). The weighted mean strength of the electromagnetic force over this redshift range in this location in the universe is da/a = (alpha_z - alpha_0)/alpha_0 = (-2.18 +/- 7.27) X 10^{-5}, i.e. we find no evidence for a temporal change from the 4 new very high redshift measurements. When the 4 new measurements are combined with a large existing sample of lower redshift measurements, a new limit on possible spatial variation of da/a is marginally preferred over a no-variation model at the 3.7 sigma level.

We study the conformal structure of exotic (non-big-bang) singularity universes using the hybrid big-bang/exotic singularity/big-bang and big-rip/exotic singularity/big-rip models by investigating their appropriate Penrose diagrams. We show that the diagrams have the standard structure for the big-bang and big-rip and that exotic singularities appear just as the constant time hypersurfaces for the time of a singularity and because of their geodesic completeness are potentially transversable. We also comment on some applications and extensions of the Penrose diagram method in studying exotic singularities.

This paper evaluates some important aspects of the multiverse concept. Firstly, the most realistic opportunity for it which is the spacetime variability of the physical constants and may deliver worlds with different physics, hopefully fulfilling the conditions of the anthropic principles. Then, more esoteric versions of the multiverse being the realisation of some abstract mathematics or even logic. Finally, it evaluates the big challenge of getting any signal from "other universes" using recent achievements of the quantum theory.

We study the spontaneous baryogenesis scenario in the early universe for three different frameworks of varying constants theories. We replace the constants by dynamical scalar fields playing the role of thermions. We first obtain the results for baryogenesis driven by the varying gravitational constant, $G$, as in the previous literature, then challenge the problem for varying fine structure constant $\alpha$ models as well as for varying speed of light $c$ models. We show that in each of these frameworks the current observational value of the baryon to entropy ratio, $\eta_B \sim 8.6 \cdot 10^{-11}$, can be obtained for large set of parameters of dynamical constants models as well as the decoupling temperature, and the characteristic cut-off length scale.

We modify the standard relativistic dispersion relation in a way which breaks Lorentz symmetry - the effect is predicted in a high-energy regime of some modern theories of quantum gravity. We show that it is possible to realise this scenario within the framework of Rainbow Gravity which introduces two new energy-dependent functions $f_1(E)$ and $f_2(E)$ into the dispersion relation. Additionally, we assume that the gravitational constant $G$ and the cosmological constant $\Lambda$ also depend on energy $E$ and introduce the scaling function $h(E)$ in order to express this dependence. For cosmological applications we specify the functions $f_1$ and $f_2$ in order to fit massless particles which allows us to derive modified cosmological equations. Finally, by using Hubble+SNIa+BAO(BOSS+Lyman $\alpha$)+CMB data, we constrain the energy scale $E_{LV}$ to be at least of the order of $10^{16}$ GeV at $1\sigma$ which is the GUT scale or even higher $10^{17}$ GeV at $3\sigma$. Our claim is that this energy can be interpreted as the decoupling scale of massless particles from spacetime Lorentz violating effects.

The models of cyclic universes and cyclic multiverses based on the alternative gravity theories of varying constants are considered.

Variation of the speed of light is quite a debated issue in cosmology with some benefits, but also with some controversial concerns. Many approaches to develop a consistent varying speed of light (VSL) theory have been developed recently. Although a lot of theoretical debate has sprout out about their feasibility and reliability, the most obvious and straightforward way to discriminate and check if such theories are really workable has been missed out or not fully employed. What is meant here is the comparison of these theories with observational data in a fully comprehensive way. In this paper we try to address this point i.e., by using the most updated cosmological probes, we test three different candidates for a VSL theory (Barrow \& Magueijo, Avelino \& Martins, and Moffat) signal. We consider many different ans\"{a}tze for both the functional form of $c(z)$ (which cannot be fixed by theoretical motivations) and for the dark energy dynamics, in order to have a clear global picture from which we extract the results. We compare these results using a reliable statistical tool such as the Bayesian Evidence. We find that the present cosmological data is perfectly compatible with any of these VSL scenarios, but in one case (Moffat model) we have a higher Bayesian Evidence ratio in favour of VSL than in the standard $c=$ constant $\Lambda$CDM scenario. Moreover, in such a scenario, the VSL signal can help to strengthen constraints on the spatial curvature (with indication toward an open universe), to clarify some properties of dark energy (exclusion of a cosmological constant at $2\sigma$ level) and is also falsifiable in the nearest future due to some peculiar issues which differentiate this model from the standard model. Finally, we have applied some priors which come from cosmology and, in particular, from information theory and gravitational thermodynamics.

Modified gravity theories with a screening mechanism have acquired much interest recently in the quest for a viable alternative to General Relativity on cosmological scales, given their intrinsic property of being able to pass Solar System scale tests and, at the same time, to possibly drive universe acceleration on much larger scales. Here, we explore the possibility that the same screening mechanism, or its breaking at a certain astrophysical scale, might be responsible of those gravitational effects which, in the context of general relativity, are generally attributed to Dark Matter. We consider a recently proposed extension of covariant Galileon models in the so-called "beyond Horndeski" scenario, where a breaking of the Vainshtein mechanism is possible and, thus, some peculiar observational signatures should be detectable and make it distinguishable from general relativity. We apply this model to a sample of clusters of galaxies observed under the \textit{CLASH} survey, using both new data from gravitational lensing events and archival data from X-ray intra-cluster hot gas observations. In particular, we use the latter to model the gas density, and then use it as the only ingredient in the matter clusters' budget to calculate the expected lensing convergence map. Results show that, in the context of this extended Galileon, the assumption of having only gas and no Dark Matter at all in the clusters is able to match observations. We also obtain narrow and very interesting bounds on the parameters which characterize this model. In particular, we find that, at least for one of them, the general relativity limit is excluded at $2\sigma$ confidence level, thus making this model clearly statistically different and competitive with respect to general relativity.

Using the idea of regularisation of singularities due to the variability of the fundamental constants in cosmology we study the cyclic universe models. We find two models of oscillating and non-singular mass density and pressure ("non-singular" bounce) regularised by varying gravitational constant $G$ despite the scale factor evolution is oscillating and having sharp turning points ("singular" bounce). Both violating (big-bang) and non-violating (phantom) null energy condition models appear. Then, we extend this idea onto the multiverse containing cyclic individual universes with either growing or decreasing entropy though leaving the net entropy constant. In order to get an insight into the key idea, we consider the doubleverse with the same geometrical evolution of the two "parallel" universes with their physical evolution (physical coupling constants $c(t)$ and $G(t)$) being different. An interesting point is that there is a possibility to exchange the universes at the point of maximum expansion -- the fact which was already noticed in quantum cosmology. Similar scenario is also possible within the framework of Brans-Dicke theory where varying $G(t)$ is replaced by the dynamical Brans-Dicke field $\phi(t)$ though these theories are slightly different.

In this paper we extend a new method to measure possible variation of the speed of light by using Baryon Acoustic Oscillations and the Hubble function presented in our earlier paper [V. Salzano, M. P. D\c{a}browski, and R. Lazkoz, Phys. Rev. D93, 063521 (2016)] onto an inhomogeneous model of the universe. The method relies on the fact that there is a simple relation between the angular diameter distance $(D_{A})$ maximum and the Hubble function $(H)$ evaluated at the same maximum-condition redshift, which includes speed of light $c$. One limit of such method was the assumption of null spatial curvature (even if we showed that even a non-zero curvature would have negligible effects). Here, we move one step further: we explicitly assume a model with intrinsic non-null curvature, and calculate the exact relation between $D_{A}$ and $H$ in this case. Then, we evaluate if current or future missions such as SKA can be sensitive enough to detect any such kind of spatial variation of $c$ which can perhaps be related to the recently observed spatial variation of the fine structure constant (an effect known as $\alpha$-dipole).

In [1] a new method to measure the speed of light through Baryon Acoustic Oscillations (BAO) was introduced. Here, we describe in much more detail the theoretical basis of that method, its implementation, and give some newly updated results about its application to the forecast data. In particular, we will show that SKA will be able to detect a 1% variation (if any) in the speed of light at 3$\sigma$ level. Smaller signals will be hardly detected by already-planned future galaxy surveys, but we give indications about what sensitivity requirements should a survey fulfill in order to be successful.

We formulate the basic framework of thermodynamical entropic force cosmology which allows variation of the gravitational constant $G$ and the speed of light $c$. Three different approaches to the formulation of the field equations are presented. Some cosmological solutions for each framework are given and one of them is tested against combined observational data (supernovae, BAO, and CMB). From the fit of the data it is found that the Hawking temperature numerical coefficient $\gamma$ is two to four orders of magnitude less than usually assumed on the geometrical ground theoretical value of $O(1)$ and that it is also compatible with zero. Besides, in the entropic scenario we observationally test that the fit of the data is allowed for the speed of light $c$ growing and the gravitational constant $G$ diminishing during the evolution of the universe. We also obtain a bound on the variation of $c$ to be $\Delta c/c \propto 10^{-5} >0$ which is at least one order of magnitude weaker than the quasar spectra observational bound.

We find the series of example theories for which the relativistic limit of maximum tension $F_{max} = c^4/4G$ represented by the entropic force can be abolished. Among them the varying constants theories, some generalized entropy models applied both for cosmological and black hole horizons as well as some generalized uncertainty principle models.

We consider spherically symmetric inhomogeneous pressure Stephani universes, the center of symmetry being our location. The main feature of these models is that comoving observers do not follow geodesics. In particular, comoving perfect fluids have necessarily a radially dependent pressure. We consider a subclass of these models characterized by some inhomogeneity parameter $\beta$. We show that also the velocity of sound, like the (effective) equation of state parameter, of comoving perfect fluids acquire away from the origin a time and radial dependent change proportional to $\beta$. In order to produce a realistic universe accelerating at late times without dark energy component one must take $\beta < 0$. The redshift gets a modified dependence on the scale factor $a(t)$ with a relative modification of $-9\%$ peaking at $z\sim 4$ and vanishing at the big-bang and today on our past lightcone. The equation of state parameter and the speed of sound of dustlike matter (corresponding to a vanishing pressure at the center of symmetry $r=0$) behave in a similar way and away from the center of symmetry they become negative -- a property usually encountered for the dark energy component only. In order to mimic the observed late-time accelerated expansion, the matter component must significantly depart from standard dust, presumably ruling this subclass of Stephani models out as a realistic cosmology. The only way to accept these models is to keep all standard matter components of the universe including dark energy and take an inhomogeneity parameter $\beta$ small enough.

We discuss minisuperspace models within the framework of varying physical constants theories including $\Lambda$-term. In particular, we consider the varying speed of light (VSL) theory and varying gravitational constant theory (VG) using the specific ans\"atze for the variability of constants: $c(a) = c_0 a^n$ and $G(a)=G_0 a^q$. We find that most of the varying $c$ and $G$ minisuperspace potentials are of the tunneling type which allows to use WKB approximation of quantum mechanics. Using this method we show that the probability of tunneling of the universe "from nothing" ($a=0)$ to a Friedmann geometry with the scale factor $a_t$ is large for growing $c$ models and is strongly suppressed for diminishing $c$ models. As for $G$ varying, the probability of tunneling is large for $G$ diminishing, while it is small for $G$ increasing. In general, both varying $c$ and $G$ change the probability of tunneling in comparison to the standard matter content (cosmological term, dust, radiation) universe models.

In this letter we describe a new method to use Baryon Acoustic Oscillations (BAO) to derive a constraint on the possible variation of the speed of light. The method relies on the fact that there is a simple relation between the angular diameter distance $(D_{A})$ maximum and the Hubble function $(H)$ evaluated at the same maximum-condition redshift, which includes speed of light $c$. We note the close analogy of the BAO probe with a laboratory experiment: here we have $D_{A}$ which plays the role of a standard (cosmological) ruler, and $H^{-1}$, with the dimension of time, as a (cosmological) clock. We evaluate if current or future missions such as Euclid can be sensitive enough to detect any variation of $c$.

Exact luminosity distance and apparent magnitude formulas are applied to Union2 557 supernovae sample in order to constrain possible position of an observer outside of the center of symmetry in spherically symmetric inhomogeneous pressure Stephani universes which are complementary to inhomogeneous density Lema\^itre-Tolman-Bondi (LTB) void models. Two specific models are investigated. The first which allows a barotropic equation of state at the center of symmetry with no scale factor function being specified (model IIA), and the second which has no barotropic equation of state at the center, but has an explicit dust-like scale factor evolution (model IIB). It is shown that even at $3\sigma$ CL, an off-center observer cannot be further than about 4.4 Gpc away from the center of symmetry which is comparable to the reported size of a void in LTB models with the most likely value of the distance from the center about 341 Mpc for model IIA and 68 Mpc for model IIB. The off-center observer cannot be farther away from the center than about 577 Mpc for model IIB at $3\sigma$ CL. It is evaluated that the best-fit parameters which characterize inhomogeneity are: $\Omega_{inh} = 0.77$ (dimensionless - model IIA) and $\alpha = 7.31 \cdot 10^{-9}$ $(s/km)^{2/3} Mpc^{-4/3}$ (model IIB).

Various classes of exotic singularity models have been studied as possible mimic models for the observed recent acceleration of the universe. Here we further study one of these classes and, under the assumption that they are phenomenological toy models for the behavior of an underlying scalar field which also couples to the electromagnetic sector of the theory, obtain the corresponding behavior of the fine-structure constant $\alpha$ for particular choices of model parameters that have been previously shown to be in reasonable agreement with cosmological observations. We then compare this predicted behavior with available measurements of $\alpha$, thus constraining this putative coupling to electromagnetism. We find that values of the coupling which would provide a good fit to spectroscopic measurements of $\alpha$ are in more than three-sigma tension with local atomic clock bounds. Future measurements by ESPRESSO and ELT-HIRES will provide a definitive test of these models.

We derive a redshift drift formula within the framework of varying speed of light (VSL) theory using the specific ansatz for the variability of $c(t) = c_0 a^n(t)$. We show that negative values of the parameter $n$, which correspond to diminishing value of the speed of light during the evolution of the universe, effectively rescales dust matter to become little negative pressure matter, and the cosmological constant to became phantom. Positive values of $n$ (growing $c(t)$) make VSL model to become more like Cold Dark Matter (CDM) model. Observationally, there is a distinction between the VSL model and the $\Lambda$CDM model for the admissible values of the parameter $n \sim - 10^{-5}$, though it will be rather difficult to detect by planned extremely large telescopes (E-ELT, TMT, GMT) within their accuracy.

We derive a luminosity distance formula for the varying speed of light (VSL) theory which involves higher order characteristics of expansion such as jerk, snap and lerk which can test the impact of varying $c$ onto the evolution of the universe. We show that the effect of varying $c$ is possible to be isolated due to the relations connecting observational parameters already by measuring the second-order term in redshift $z$ unless there is a redundancy between the curvature and an exotic fluid of cosmic strings scaling the same way as the curvature.

We discuss how dynamical dark energy universes with exotic singularities may be distinguished from the standard $\Lambda$CDM model on the basis of their redshift drift signal, for which measurements both in the acceleration phase and in the deep matter era will be provided by forthcoming astrophysical facilities. Two specific classes of exotic singularity models are studied: sudden future singularity models and finite scale factor singularity models. In each class we identify the models which can mimic $\Lambda$CDM and play the role of dark energy as well as models for which redshift drift signals are significantly different from $\Lambda$CDM and the test can differentiate between them.

We review the variety of new singularities in homogeneous and isotropic FRW cosmology which differ from standard Big-Bang and Big-Crunch singularities and suggest how the nature of these singularities can be influenced by the varying fundamental constants.

We derive a redshift drift formula for the spherically symmetric inhomogeneous pressure Stephani universes which are complementary to the spherically symmetric inhomogeneous density Lema\^itre-Tolman-Bondi models. We show that there is a clear difference between redshift drift predictions for these two models as well as between the Stephani models and the standard $\Lambda$CDM Friedmann models. The Stephani models have positive drift values at small redshift and behave qualitatively (but not quantitatively) as the $\Lambda$CDM models at large redshift, while the drift for LTB models is always negative. This prediction may perhaps be tested in future telescopes such as European Extremely Large Telescope (EELT), Thirty Meter Telescope (TMT), Giant Magellan Telescope (GMT), and especially, in gravitational wave interferometers DECi-Hertz Interferometer Gravitational Wave Observatory and Big Bang Observer (DECIGO/BBO), which aim at low redshift.

Varying physical constant cosmologies were claimed to solve standard cosmological problems such as the horizon, the flatness and the $\Lambda$-problem. In this paper, we suggest yet another possible application of these theories: solving the singularity problem. By specifying some examples we show that various cosmological singularities may be regularized provided the physical constants evolve in time in an appropriate way.

We discuss combined constraints, coming from the cosmic microwave background shift parameter $\mathcal{R}$, baryon acoustic oscillations (BAO) distance parameter $\mathcal{A}$, and from the latest type Ia supernovae data, imposed on cosmological models which allow sudden future singularities of pressure. We show that due to their weakness such sudden singularities may happen in the very near future and that at present they can mimic standard dark energy models.

By setting some special boundary conditions in the variational principle we obtain junction conditions for the five-dimensional $f(R)$ gravity which in the Einstein limit $f(R)\rightarrow R$ transform into the standard Randall-Sundrum junction conditions. We apply these junction conditions to a particular model of a Friedmann universe on the brane and show explicitly that the limit gives the standard Randall-Sundrum model Friedmann equation.

Taking a spacetime average as a measure of the strength of singularities we show that big-rips (type I) are stronger than big-bangs. The former have infinite spacetime averages while the latter have them equal to zero. The sudden future singularities (type II) and $w-$singularities (type V) have finite spacetime averages. The finite scale factor (type III) singularities for some values of the parameters may have an infinite average and in that sense they may be considered stronger than big-bangs.

Current observational evidence does not yet exclude the possibility that dark energy could be in the form of phantom energy. A universe consisting of a phantom constituent will be driven toward a drastic end known as the `Big Rip' singularity where all the matter in the universe will be destroyed. Motivated by this possibility, other evolutionary scenarios have been explored by Barrow, including the phenomena which he called Sudden Future Singularities (SFS). In such a model it is possible to have a blow up of the pressure occurring at sometime in the future evolution of the universe while the energy density would remain unaffected. The particular evolution of the scale factor of the universe in this model that results in a singular behaviour of the pressure also admits acceleration in the current era. In this paper we will present the results of our confrontation of one example class of SFS models with the available cosmological data from high redshift supernovae, baryon acoustic oscillations (BAO) and the cosmic microwave background (CMB). We then discuss the viability of the model in question as an alternative to dark energy.

By the application of the generalized Israel junction conditions we derive cosmological equations for the fourth-order $f(R)$ brane gravity and study their cosmological solutions. We show that there exists a non-static solution which describes a four-dimensional de-Sitter $(dS_4)$ brane embedded in a five-dimensional anti-de-Sitter $(AdS_5)$ bulk for a vanishing Weyl tensor contribution. On the other hand, for the case of a non-vanishing Weyl tensor contribution, there exists a static brane only. We claim that in order to get some more general non-static $f(R)$ brane configurations, one needs to admit a dynamical matter energy-momentum tensor in the bulk rather than just a bulk cosmological constant.

We discuss various types of exotic (non-standard) singularities in the Universe: a Big-Rip (BR or type I), a Sudden Future Singularity (SFS or type II), a Generalized Sudden Future Singularity, a Finite Scale Factor singularity (FSF or type III), a Big-Separation (BS or type IV) and a $w$-singularity. They are characterized by violation of all or some of the energy conditions which results in a blow-up of all or some of the physical quantities: the scale factor, the energy density, the pressure, and the barotropic index. We relate the emergence of these singularities with physical theories (superstring, brane, higher-order gravity, loop quantum cosmology). We show how the models involving exotic singularities may serve as dark energy by applying the observational data. In particular, we show that some of these exotic singularities (though being of a weak type according to relativistic definitions) may occur in the near future of the universe.

We discuss a very general theory of gravity, of which Lagrangian is an arbitrary function of the curvature invariants, on the brane. In general, the formulation of the junction conditions (except for Euler characteristics such as Gauss-Bonnet term) leads to the powers of the delta function and requires regularization. We suggest the way to avoid such a problem by imposing the metric and its first derivative to be regular at the brane, the second derivative to have a kink, the third derivative of the metric to have a step function discontinuity, and no sooner as the fourth derivative of the metric to give the delta function contribution to the field equations. Alternatively, we discuss the reduction of the fourth-order gravity to the second order theory by introducing extra scalar and tensor fields: the scalaron and the tensoron. In order to obtain junction conditions we apply two methods: the application of the Gauss-Codazzi formalism and the application of the generalized Gibbons-Hawking boundary terms which are appended to the appropriate actions. In the most general case we derive junction conditions without assuming the continuity of the scalaron and the tensoron on the brane. The derived junction conditions can serve studying the cosmological implications of the higher-order brane gravity models.

The superenergy of the universe is a tensorial quantity and it is a general relativistic analogue of the Appell's energy of acceleration in classical mechanics. We propose the measurement of this quantity by the observational parameters such as the Hubble parameter, the deceleration parameter, the jerk and the snap (kerk) known as statefinders. We show that the superenergy of gravity requires only the Hubble and deceleration parameter to be measured, while the superenergy of matter requires also the measurement of the higher-order characteristics of expansion: the jerk and the snap. In such a way, the superenergy becomes another parameter characterizing the evolution of the universe.

We find an explicit cosmological model which allows a special type of cosmological singularity which we call a $w$-singularity. This singularity has the scale factor finite, the energy density and pressure vanishing, and the only singular behaviour appears in a time-dependent barotropic index $w(t)$. It is different from the type IV cosmological singularity in that it does not exhibit the divergence of the higher derivatives of the Hubble parameter and from the big-brake since it does not fulfill the anti-Chaplygin gas equation of state. We also find an interesting duality between the $w$-singularities and the big-bang singularities. Physical examples of $w$-singularities appear in $f(R)$, scalar field and brane cosmologies.

A redshift-magnitude relation for the two exact non-uniform pressure spherically symmetric Stephani universes is presented. The Kristian-Sachs method expanding the relativistic quantities in series is used, but only first order terms in redshift $z$ are considered. The numerical results are given both for centrally placed and non-centrally placed observers. In the former case the redshift-magnitude relation does not depend on the direction in the sky and the Friedman limit can be easily performed. It appears that the effect of spatial dependence of pressure is similar to the effect of the deceleration parameter in Friedman models. In the latter case the angular dependence of the relation is important. This may serve as another possible explanation of the non-compatibility of the theoretical curve of the redshift-magnitude relation with observations for large redshift objects in the Friedman universe. On the other hand, comparing the magnitudes of equal redshifts objects in different directions in the sky, one can test the reliability of these models.

We derive the most general junction conditions for the fourth-order brane gravity constructed of arbitrary functions of curvature invariants. We reduce these fourth-order theories to second order theories at the expense of introducing new scalar and tensor fields - the scalaron and the tensoron. In order to obtain junction conditions we apply the method of generalized Gibbons-Hawking boundary terms which are appended to the appropriate actions. After assuming the continuity of the scalaron and the tensoron on the brane, we recover junction conditions for such general brane universe models previously obtained by different methods. The derived junction conditions can serve studying the cosmological implications of the higher-order brane gravity models.

Conformal transformations are frequently used tools in order to study relations between various theories of gravity and the Einstein relativity. In this paper we discuss the rules of these transformations for geometric quantities as well as for the matter energy-momentum tensor. We show the subtlety of the matter energy-momentum conservation law which refers to the fact that the conformal transformation "creates" an extra matter term composed of the conformal factor which enters the conservation law. In an extreme case of the flat original spacetime the matter is "created" due to work done by the conformal transformation to bend the spacetime which was originally flat. We discuss how to construct the conformally invariant gravity theories and also find the conformal transformation rules for the curvature invariants $R^2$, $R_{ab}R^{ab}$, $R_{abcd}R^{abcd}$ and the Gauss-Bonnet invariant in a spacetime of an arbitrary dimension. Finally, we present the conformal transformation rules in the fashion of the duality transformations of the superstring theory. In such a case the transitions between conformal frames reduce to a simple change of the sign of a redefined conformal factor.

We consider an alternative mechanism for the production of the cosmic microwave background (CMB) radiation. It is basically due to vacuum pair creation (VPC) of vector bosons (W and Z) as a consequence of a rapid W and Z mass generation during the electroweak phase transition in the early Universe. The mechanism is as follows: after their pair crreation, the vector bosons may either annihilate directly into photons or decay into leptons and quarks which subsequently annihilate as lepton-antilepton and quark-antiquark pairs into photons. Preliminary estimates show that the number of CMB photons obtained this way can be sufficient to explain the presently observed CMB photon density. In this contribution we present an exactly soluble model for vacuum pair creation kinetics.

We find the rules of the conformal transformation for the energetic quantities such as the Einstein energy-momentum complex, the Bergmann-Thomson angular momentum complex, the superenergy tensor, and the angular supermomentum tensor of gravitation and matter. We show that the conformal transformation rules for the matter parts of both the Einstein complex and the Bergmann-Thomson complex are fairly simple, while the transformation rules for their gravitational parts are more complicated. We also find that the transformational rules of the superenergy tensor of matter and the superenergy tensor of gravity are quite complicated except for the case of a pure gravity. In such a special case the superenergy density as well as the sum of the superenergy density and the matter energy density are invariants of the conformal transformation. Besides, in that case, a conformal invariant is also the Bel-Robinson tensor which is a part of the superenergy tensor. As for the angular supermomentum tensor of gravity - it emerges that its transformational rule even for a pure gravity is quite complicated but this is not the case for the angular supermomentum tensor of matter. Having investigated some technicalities of the conformal transformations, we also find the conformal transformation rule for the curvature invariants and, in particular, for the Gauss-Bonnet invariant in a spacetime of arbitrary dimension.

We discuss a general fourth-order theory of gravity on the brane. In general, the formulation of the junction conditions (except for Euler characteristics such as Gauss-Bonnet term) leads to the higher powers of the delta function and requires regularization. We suggest the way to avoid such a problem by imposing the metric and its first derivative to be regular at the brane, while the second derivative to have a kink, the third derivative of the metric to have a step function discontinuity, and no sooner as the fourth derivative of the metric to give the delta function contribution to the field equations. Alternatively, we discuss the reduction of the fourth-order gravity to the second-order theory by introducing an extra tensor field. We formulate the appropriate junction conditions on the brane. We prove the equivalence of both theories. In particular, we prove the equivalence of the junction conditions with different assumptions related to the continuity of the metric along the brane.

We discuss the constraints coming from current observations of type Ia supernovae on cosmological models which allow sudden future singularities of pressure (with the scale factor and the energy density regular). We show that such a sudden singularity may happen in the very near future (e.g. within ten million years) and its prediction at the present moment of cosmic evolution cannot be distinguished, with current observational data, from the prediction given by the standard quintessence scenario of future evolution. Fortunately, sudden future singularities are characterized by a momentary peak of infinite tidal forces only; there is no geodesic incompletness which means that the evolution of the universe may eventually be continued throughout until another ``more serious'' singularity such as Big-Crunch or Big-Rip.

I discuss the dark energy characterized by the violation of the null energy condition ($\varrho + p \geq 0$), dubbed phantom. Amazingly, it is admitted by the current astronomical data from supernovae. We discuss both classical and quantum cosmological models with phantom as a source of matter and present the phenomenon called phantom duality.

We discuss exotic singularities in the evolution of the universe motivated by the progress of observations in cosmology. Among them there are: Big-Rip (BR), Sudden Future Singularities (SFS), Generalized Sudden Future Singularities (GSFS), Finite Density Singularities (FD), type III, and type IV singularities. We relate some of these singularities with higher-order characteristics of expansion such as jerk and snap. We also discuss the behaviour of pointlike objects and classical strings on the approach to these singularities.

We discuss observational constrains coming from supernovae imposed on the behaviour of the Randall-Sundrum models. The data indicates that, under the assumption that we admit zero pressure dust matter on the brane, the cosmological constant is still needed to explain current observations. The observations from supernovae give a large value of the density parameter for brane matter $\Omega_{\lambda,0} \simeq 0.01$ as the best fit. For high redshifts $z > 1.2$, the difference between the brane model and the $\Lambda$CDM (Perlmutter) model becomes detectable observationally. From the maximum likelihood method we obtained the favored value of $\Omega_{\lambda,0} = 0.004 \pm 0.016$ for $\Omega_{k,0}=0$ and $\Omega_{\mathrm{m},0}=0.3$. This gives the limit $\Omega_{\lambda,0} < 0.02$ at $1 \sigma$ level. While the model with brane effects is preferred by the supernovae type Ia data, the model without brane fluid is still statistically admissible. For comparison the limit on brane effects which comes from CMB anisotropies and BBN is also obtained. The uncertainty in the location of the first peak gives a stronger limit $\Omega_{\lambda,0} < 1.0 \cdot 10^{-12}$, whereas from BBN we obtain that $\Omega_{\lambda,0} < 1.0 \cdot 10^{-27}$. However, both very strict limits are obtained with the assumption that brane effects do not change the physics in the pre-recombination era, while the SN Ia limit is model independent.We demonstrate that the fit to supernovae data can also be obtained if we admit the phantom matter $p = - (4/3) \varrho$ on the brane, where this matter mimics the influence of the cosmological constant. Finally, we propose to check for dark radiation and brane tension by the application of the angular diameter of galaxies minimum value test.

Following the observational evidence for cosmic acceleration which may exclude a possibility for the universe to recollapse to a second singularity, we review alternative scenarios of its future evolution. Although the de Sitter asymptotic state is still an option, some other asymptotic states which allow new types of singularities such as Big-Rip (due to a phantom matter) and sudden future singularities are also admissible and are reviewed in detail. The reality of these singularities which comes from the relation to observational characteristics of the universe expansion are also revealed and widely discussed.

We apply the formalism of quantum cosmology to models containing a phantom field. Three models are discussed explicitly: a toy model, a model with an exponential phantom potential, and a model with phantom field accompanied by a negative cosmological constant. In all these cases we calculate the classical trajectories in configuration space and give solutions to the Wheeler-DeWitt equation in quantum cosmology. In the cases of the toy model and the model with exponential potential we are able to solve the Wheeler-DeWitt equation exactly. For comparison, we also give the corresponding solutions for an ordinary scalar field. We discuss in particular the behaviour of wave packets in minisuperspace. For the phantom field these packets disperse in the region that corresponds to the Big Rip singularity. This thus constitutes a genuine quantum region at large scales, described by a regular solution of the Wheeler-DeWitt equation. For the ordinary scalar field, the Big-Bang singularity is avoided. Some remarks on the arrow of time in phantom models as well as on the relation of phantom models to loop quantum cosmology are given.

We discuss the behaviour of strings propagating in spacetimes which allow future singularities of either a sudden future or a Big-Rip type. We show that in general the invariant string size remains finite at sudden future singularities while it grows to infinity at a Big-Rip. This claim is based on the discussion of both the tensile and null strings. In conclusion, strings may survive a sudden future singularity, but not a Big-Rip where they are infinitely stretched.

Following recent fit of supernovae data to Brans-Dicke theory which favours the model with $\omega = - 3/2$ \cite{fabris} we discuss the status of this special case of Brans-Dicke cosmology in both isotropic and anisotropic framework. It emerges that the limit $\omega = -3/2$ is consistent only with the vacuum field equations and it makes such a Brans-Dicke theory conformally invariant. Then it is an example of the conformal relativity theory which allows the invariance with respect to conformal transformations of the metric. Besides, Brans-Dicke theory with $\omega = -3/2$ gives a border between a standard scalar field model and a ghost/phantom model. In this paper we show that in $\omega = -3/2$ Brans-Dicke theory, i.e., in the conformal relativity there are no isotropic Friedmann solutions of non-zero spatial curvature except for $k=-1$ case. Further we show that this $k=-1$ case, after the conformal transformation into the Einstein frame, is just the Milne universe and, as such, it is equivalent to Minkowski spacetime. It generally means that only flat models are fully consistent with the field equations. On the other hand, it is shown explicitly that the anisotropic non-zero spatial curvature models of Kantowski-Sachs type are admissible in $\omega = -3/2$ Brans-Dicke theory. It then seems that an additional scale factor which appears in anisotropic models gives an extra deegre of freedom and makes it less restrictive than in an isotropic Friedmann case.

We discuss a more general class of phantom ($p < -\varrho$) cosmologies with various forms of both phantom ($w < -1$), and standard ($w > -1$) matter. We show that many types of evolution which include both Big-Bang and Big-Rip singularities are admitted and give explicit examples. Among some interesting models, there exist non-singular oscillating (or "bounce") cosmologies, which appear due to a competition between positive and negative pressure of variety of matter content. From the point of view of the current observations the most interesting cosmologies are the ones which start with a Big-Bang and terminate at a Big-Rip. A related consequence of having a possibility of two types of singularities is that there exists an unstable static universe approached by the two asymptotic models - one of them reaches Big-Bang, and another reaches Big-Rip. We also give explicit relations between density parameters $\Omega$ and the dynamical characteristics for these generalized phantom models, including higher-order observational characteristics such as jerk and "kerk". Finally, we discuss the observational quantities such as luminosity distance, angular diameter, and source counts, both in series expansion and explicitly, for phantom models. Our series expansion formulas for the luminosity distance and the apparent magnitude go as far as to the fourth-order in redshift $z$ term, which includes explicitly not only the jerk, but also the "kerk" (or "snap") which may serve as an indicator of the curvature of the universe.

We link observational parameters such as the deceleration parameter, the jerk, the kerk (snap) and higher-order derivatives of the scale factor, called statefinders, to the conditions which allow to develop sudden future singularities of pressure with finite energy density. In this context, and within the framework of Friedmann cosmology, we also propose higher-order energy conditions which relate time derivatives of the energy density and pressure which may be useful in general relativity.

The Friedman equation for the universe with arbitrary curvature $(k = 0, \pm 1)$, filled with mutually noninteracting pressureless dust, radiation, cosmological constant, and strings is considered. We assume the string domination scenario for the evolution of the latter component. Moreover, we discuss the simplest possibility for the scaling of the string energy density: $\varrho \propto R^{-2}$. For such models we write down the explicit solution of the Friedman equation. We realize that corresponding cosmological models do not essentially differ from those without strings. We find an analytic formula for the radial coordinate $\chi$ of a galaxy with a redshift z and express it in terms of astronomical parameters. This relation is then used for the derivation of the astrophysical formulas such as luminosity distance, angular diameter, and source counts, which may serve for testing the string-dominated universe. It seems that the most sensitive test, at least from the formal point of view, is the formula for the number of galaxies $N(z)$ corresponding to a given value of the redshift. We show that the maximum of $N(z)$ strongly depends on the density of strings, especially if the density is large enough to explain the $\Omega$ problem. Other tests are also proposed.

Conformal relativity theory which is also known as Hoyle-Narlikar theory has recently been given some new interest. It is an extended relativity theory which is invariant with respect to conformal transformations of the metric. In this paper we show how conformal relativity is related to the Brans-Dicke theory and to the low-energy-effective superstring theory. We show that conformal relativity action is equaivalent to a transformed Brans-Dicke action for Brans-Dicke parameter $\omega = -3/2$ in contrast to a reduced (graviton-dilaton) low-energy-effective superstring action which corresponds to a Brans-Dicke action with Brans-Dicke parameter $\omega = -1$. In fact, Brans-Dicke parameter $\omega =-3/2$ gives a border between a standard scalar field evolution and a ghost. We also present basic cosmological solutions of conformal relativity in both Einstein and string frames. The Eintein limit for flat conformal cosmology solutions is unique and it is flat Minkowski space. This requires the scalar field/mass evolution instead of the scale factor evolution in order to explain cosmological redshift. It is interesting that like in ekpyrotic/cyclic models, a possible transition through a singularity in conformal cosmology in the string frame takes place in the weak coupling regime.

We discuss the Einstein energy-momentum complex and the Bergmann-Thomson angular momentum complex in general relativity and calculate them for space-time homogeneous Goedel universes. The calculations are performed for a dust acausal model and for a scalar-field causal model. It is shown that the Einstein pseudotensor is traceless, not symmetric, the gravitational energy is "density" is negative and the gravitational Poynting vector vanishes. Significantly, the total (gravitational and matter) energy "density" fro the acausal model is zero while for the casual model it is negative.The Bergmann-Thomson angular momentum complex does not vanish for both G\"odel models.

We discuss a class of phantom ($p < - \varrho$) cosmological models. Except for phantom we admit various forms of standard types of matter and discuss the problem of singularities for these cosmologies. The singularities are different from those of standard matter cosmology since they appear for infinite values of the scale factor. We also find an interesting relation between the phantom models and standard matter models which is like the duality symmetry of string cosmology.

We discuss observational constrains coming from supernovae Ia \cite{Perlmutter99} imposed on the behaviour of the Randall-Sundrum models. In the case of dust matter on the brane, the difference between the best-fit general relativistic model with a $\Lambda$-term \cite{Perlmutter99} and the best-fit brane models becomes detectable for redshifts $z > 0.6$. It is interesting that brane models predict brighter galaxies for such redshifts which is in agreement with the measurement of the $z = 1.7$ supernova \cite{Riess01} and with the New Data from the High Z Supernovae Search Team \cite{schmit02}. We also demonstrate that the fit to supernovae data can also be obtained, if we admit the "super-negative" dark energy $p = - (4/3) \varrho$ on the brane, where the dark energy in a way mimics the influence of the cosmological constant. It also appears that the dark energy enlarges the age of the universe which is demanded in cosmology. Finally, we propose to check for dark radiation and brane tension by the application of the angular diameter of galaxies minimum value test.

We apply methods of dynamical systems to study the behaviour of the Randall-Sundrum models. We determine evolutionary paths for all possible initial conditions in a 2-dimensional phase space and we investigate the set of accelerated models. The simplicity of our formulation in comparison to some earlier studies is expressed in the following: our dynamical system is a 2-dimensional Hamiltonian system, and what is more advantageous, it is free from the degeneracy of critical points so that the system is structurally stable. The phase plane analysis of Randall-Sundrum models with isotropic Friedmann geometry clearly shows that qualitatively we deal with the same types of evolution as in general relativity, although quantitatively there are important differences.

We discuss the problem of the motion of classical strings in some black hole and cosmological spacetimes. In particular, the null string limit (zero tension) of tensile strings is considered. We present some new exact string solutions in Reissner-Nordstr\"om black hole background as well as in the Einstein Static Universe and in the Einstein-Schwarzschild (a black hole in the Einstein Static Universe) spacetime. These solutions can give some insight into a general nature of propagation of strings (cosmic and fundamental) in curved backgrounds.

We review the canonical superenergy tensor and the canonical angular supermomentum tensors in general relativity and calculate them for space-time homogeneous G\"odel universes to show that both of these tensors do not, in general, vanish. We consider both an original dust-filled pressureless acausal G\"odel model of 1949 and a scalar-field-filled causal G\"odel model of Rebou\c cas and Tiomno. For the acausal model, the non-vanishing components of superenergy of matter are different from those of gravitation. The angular supermomentum tensors of matter and gravitation do not vanish either which simply reflects the fact that G\"odel universe rotates. However, the axial (totally antisymmetric) and vectorial parts of supermomentum tensors vanish. It is interesting that superenergetic quantities are {\it sensitive} to causality in a way that superenergy density $_g S_{00}$ of gravitation in the acausal model is {\it positive}, while superenergy density $_g S_{00}$ in the causal model is {\it negative}. That means superenergetic quantities might serve as criterion of causality in cosmology and prove useful.

We discuss various superstring effective actions and, in particular, their common sector which leads to the so-called pre-big-bang cosmology (cosmology in a weak coupling limit of heterotic superstring). Then, we review the main ideas of the Horava-Witten theory which is a strong coupling limit of heterotic superstring theory. Using the conformal relationship between these two theories we present Kasner asymptotic solutions of Bianchi type IX geometries within these theories and make predictions about possible emergence of chaos. Finally, we present a possible method of generating Horava-Witten cosmological solutions out of the well-known general relativistic pre-big-bang solutions.

We discuss boson stars as possible gravitational lenses and study the lensing effect by these objects made of scalar particles. The mass and the size of a boson star may vary from an individual Newtonian object similar to the Sun to the general relativistic size and mass of a galaxy close to its Schwarzschild radius. We assume boson stars to be transparent which allows the light to pass through them though the light is gravitationally deflected. We assume boson stars of the mass $M = 10^{10}M_\odot$ to be on non-cosmological distance from the observer. We discuss the lens equation for these stars as well as the details of magnification. We find that there are typically three images of a star but the deflection angles may vary from arcseconds to even degrees. There is one tangential critical curve (Einstein ring) and one radial critical curve for tangential and radial magnification, respectively. Moreover, the deflection angles for the light passing in the gravitational field of boson stars can be very large (even of the order of degrees) which reflects the fact they are very strong relativistic objects. We also propose a suitable formula for the lens equation for such large deflection angles, and with the reservation that large deflection angle images are highly demagnified but in the area of the tangential critical curve, their existence may help in observational detection of suitable lenses possessing characteristic features of boson stars which could also serve as a direct evidence for scalar fields in the universe.

Bianchi type I and type IX ('Mixmaster') geometries are investigated within the framework of Ho\v{r}ava-Witten cosmology. We consider the models for which the fifth coordinate is a $S^1/Z_2$ orbifold while the four coordinates are such that the 3-space is homogeneous and has geometry of Bianchi type I or IX while the rest six dimensions have already been compactified on a Calabi-Yau space. In particular, we study Kasner-type solutions of the Bianchi I field equations and discuss Kasner asymptotics of Bianchi IX field equations. We are able to recover the isotropic 3-space solutions found by Lukas {\it et al}. Finally, we discuss if such Bianchi IX configuration can result in chaotic behaviour of these Ho\v{r}ava-Witten cosmologies.

Most of the observational claims in cosmology are based on the assumption that the universe is isotropic and homogeneous so they essentially test different types of Friedmann models. This also refers to recent observations of supernovae Ia, which, within the framework of Friedmann cosmologies give strong support to negative pressure matter and also weaken the age conflict. In this essay we drop the assumption of homogeneity, though temporarily leaving the assumption of isotropy with respect to one point, and show that supernovae data can be consistent with a model of the universe with inhomogeneous pressure known as the Stephani model. Being consistent with supernovae data we are able to get the age of the universe in this model to be about 3.8 Gyr more than in its Friedmann counterpart.

String dynamics near the photon sphere in Schwarzschild spacetime is considered on the basis of a perturbative approach with respect to a rescaled string tension as a small parameter. The perturbative string solution in the zeroth and first approximation is presented. The perturbative solution corresponds to a small deformations of the photon sphere in Schwarzschild spacetime.

We investigate Bianchi type IX ''Mixmaster'' universes within the framework of the low-energy tree-level effective action for string theory, which (when the ''stringy'' 2-form axion potential vanishes) is formally the same as the Brans-Dicke action with $\omega =-1$. We show that, unlike the case of general relativity in vacuum, there is no Mixmaster chaos in these string cosmologies. In the Einstein frame an infinite sequence of chaotic oscillations of the scale factors on approach to the initial singularity is impossible, as it was in general relativistic Mixmaster universes in the presence of stiff -fluid matter (or a massless scalar field). A finite sequence of oscillations of the scale factors approximated by Kasner metrics is possible, but it always ceases when all expansion rates become positive. In the string frame the evolution through Kasner epochs changes to a new form which reflects the duality symmetry of the theory. Again, we show that chaotic oscillations must end after a finite time. The need for duality symmetry appears to be incompatible with the presence of chaotic behaviour as $t\to 0$. We obtain our results using the Hamiltonian qualitative cosmological picture for Mixmaster models. We mention possible relations of this picture to diffeomorphism-independent methods of measuring chaos in general relativity. Finally, we discuss the effect of inhomogeneities and higher dimensions on the possible emergence of chaos within the string cosmology.

We show that homogeneous G\"odel spacetimes need not contain closed timelike curves in low-energy-effective string theories. We find exact solutions for the G\"odel metric in string theory for the full $O(\alpha ^{\prime})$ action including both dilaton and axion fields. The results are valid for bosonic, heterotic and super-strings. To first order in the inverse string tension $\alpha ^{\prime}$, these solutions display a simple relation between the angular velocity of the G\"odel universe, $\Omega ,$ and the inverse string tension of the form $\alpha ^{\prime}=1/\Omega ^2$ in the absence of the axion field. The generalization of this relationship is also found when the axion field is present.

We use the redshift-magnitude relation, as derived by D\c{a}browski (1995), for the two exact non-uniform pressure spherically symmetric Stephani universes with the observer positioned at the center of symmetry, to test the agreement of these models with recent observations of high redshift type Ia supernovae (SNIa), as reported in Perlmutter et al. (1997). By a particular choice of model parameters, we show that these models give an excellent fit to the observed redshifts and (corrected) B band apparent magnitudes of the SNIa data, but for an age of the Universe which is typically about two Gyr greater than in the corresponding Friedmann model. Based on a value of $H_0 \sim 65$ and assuming $\Lambda \geq 0$, the P97 data implies a Friedmann age of at most 13 Gyr and in fact a best-fit (for $q_0 = 0.5$) age of only 10 Gyr. Our Stephani models, on the other hand, can give a good fit to the P97 data with an age of up to 15 Gyr and could, therefore, significantly alleviate the conflict between recent cosmological and astrophysical age predictions. The choice of model parameters is quite robust: one requires only that the non-uniform pressure parameter, $a$, in one of the models is negative and satisfies $|a| \lte$ 3 km$^2$ s$^{-2}$ Mpc$^{-1}$. By allowing slightly larger, negative, values of $a$ one may `fine tune' the model to give an even better fit to the P97 data.

It is shown that growing-entropy stiff-fluid Kantowski-Sachs universes become time-symmetric (if they start with time-asymmetric phase) and isotropize. Isotropization happens without any inflationary era during the evolution since there is no cosmological term here. It seems that this approach is an alternative to inflation since the universe gets bigger and bigger approaching 'flatness'.

It is shown that the tunneling effect in quantum cosmology is possible not only at the very beginning or the very end of the evolution, but also at the moment of maximum expansion of the universe. A positive curvature expanding Friedmann universe changes its state of evolution spontaneously and completely, {\it without} any changes in the matter content, avoiding recollapse, and falling into oscillations between the nonzero values of the scale factor. On the other hand, an oscillating nonsingular universe can tunnel spontaneously to a recollapsing regime. The probability of such kind of tunneling is given explicitly. It is inversely related to the amount of nonrelativistic matter (dust), and grows from a certain fixed value to unity if the negative cosmological constant approaches zero.