Expanding space, quasars and St. Augustine's fireworks

An attempt is made to explain time non-dilation allegedly observed in quasar light curves. The explanation is based on the assumption that quasar black holes are, in some sense, foreign for our Friedmann-Robertson-Walker universe and do not participate in the Hubble flow. Although at first sight such a weird explanation requires unreasonably fine-tuned Big Bang initial conditions, we find a natural justification for it using the Milne cosmological model as an inspiration.


INTRODUCTION
"Quasar light curves do not show the effects of time dilation"-this result of the paper [2] seems incredible. The cosmological time dilation is a very basic phenomenon predicted by the Friedmann-Robertson-Walker (FRW) metric. How can quasars ignore this fundamental effect of relativistic cosmology?
Of course, extraordinary claims require extraordinary evidence. Until the puzzling results of [2] are firmly confirmed by independent observations, we cannot be certain that what is observed in [2] is a real observation of time non-dilation in quasar light curves and not a some unaccounted systematic error, for example, in data analysis.
Nevertheless, we will assume in this article that the results of [2] are correct and try to explain why quasars behave like strangers in our FRW world. Strangely enough, our main idea came from Fyodor Tyutchev's poetry (see the epigraph). The image of capricious Hebe tipping a quasar-foaming goblet on our universe was a poetic metaphor guiding this investigation and what follows can be considered as an attempt of more or less scientific incarnation of this irresistible metaphor.
The paper is organized in five parts. In the first section we consider caveats of the standard picture of expanding space as a rather subjective notion. On the other hand, in contrast to expanding space, which can be considered as a kind of coordinate artifact, cosmological time dilation is argued to be an unavoidable, coordinate independent prediction for any internal time variability of any object participating in the Hubble flow. Therefore, if quasars really manage to escape bonds of time dilation, a natural inference will be that they somehow do not participate in the Hubble flow. This quasar mystery is accompanied with some others that are discussed in the second section.
The following section gives a detailed description of the Milne model. In the ideal case, the Milne universe is a part of the Minkowski space-time with the future light-cone of the Big Bang event as its impenetrable boundary. However, the density of the Milne fundamental observers, tending to infinity on this boundary, makes unrealistic the realization of the ideal Milne universe. This conclusion is further strengthened if we consider a quantum scalar field in the Milne universe: if the boundary is assumed to be impenetrable, the corresponding choice of the initial vacuum state is not an adiabatic but a conformal vacuum that leads to the pair production phenomenon indicating the presence of physically unrealistic infinite power sources on the Milne universe's singular boundary.
In the forth section it is argued that in any realistic incarnation of the Milne universe it is expected that objects from "outside" can penetrate inside the Milne universe. Then we postulate that the existence of analogous objects in our universe is possible and we comment how these objects which do not participate in the Hubble flow due to their "otherworldly" origin, can offer an explanation of some quasar mysteries.
In the final section we provide concluding remarks and discuss how the result of this work can be related to other ideas described in the literature.
We have tried to make the paper as self-contained as possible by providing enough details of all calculations, and we hope that "the reader will not need to have his fingers at eleven places to follow an argument" [3].

SPACE EXPANSION AND ENIGMA OF TIME NON-DILATION IN QUASAR LIGHT CURVES
Expansion of space is, probably, the most familiar concept in modern cosmology. However, this beguilingly simple idea harbors many dangers of misunderstanding and misuse even for professional physicists, "most scientists think they understand it, but few agree on what it really means" [4]. It is not surprising that the subtle notion of expanding space became a subject of continuing debates [5][6][7][8][9][10][11][12][13][14][15], largely triggered by a beautiful Scientific American article by Lineweaver and Davis [4].
We can identify at least two main reasons why it is dangerous and misleading to speak about expanding space without explicitly clarifying what is really meant by this combination of words.
First of all, as special relativity teaches us, space by itself does not constitute an objective reality-only space-time does. As eloquently expressed by Minkowski "Henceforth, space by itself and time by itself are doomed to fade away into mere shadows, and only a kind of union of the two will preserve an independent reality" [16]. It does not make sense, therefore, to speak about expanding space without clarifying what foliation of space-time we have in mind and why we have chosen this specific foliation.
In general relativity understanding is further complicated by the fact that coordinates loose their usual meaning as temporal and spatial labels of events. We can use any coordinate system in general relativity to describe a given physical phenomenon and we should be very careful to separate the real physical events from mere coordinate artifacts. Moreover, when we speak about temporal development in space of some physical process, we should specify what coordinate system was chosen and for what reason. Without doing so, "to some degree we mislead both our students and ourselves when we calculate, for instance, the mercury perihelion motion without explaining how our coordinate system is fixed in space, what defines it in such a way that it cannot be rotated, by a few seconds a year, to follow the perihelion's apparent motion . . . Expressing our results in terms of the values of coordinates became a habit to such an extend, that we adhere to this habit also in general relativity, where values of coordinates are not meaningful per se" [17].
It is the symmetry of space-time, namely isotropy and homogeneity of space, which determines its preferred foliation. More precisely, the cosmological or Copernican principle states that we do not occupy a special place in the universe, which, at large, is very much the same everywhere at any given instance. Mathematically this means (see, for example, [18]) that the space-time can be foliated into spacelike slices and has a form R × Σ, where every spacelike slice Σ is a maximally symmetric space (homogeneous and isotropic) and R represents the cosmic time (a succession of instances of "now"). Geometrically the cosmic time is essentially the parameter that labels spacelike slices and each slice hypersurface is an orbit of the symmetry group of spatial homogeneity and isotropy [19].
The metric which corresponds to this preferred foliation is the celebrated Friedmann-Robertson-Walker metric (in units where c = 1): dr 2 1 − kr 2 + r 2 (dθ 2 + sin 2 θdφ 2 ) , with k > 0, k = 0 or k < 0 corresponding to the close, flat or open universes respectively. The expansion of space has a clear meaning in FRW coordinates: the scale factor a(t) depends on cosmic time. The particular form of this dependence is determined from Einstein equations. However, this fact does not mean much without connecting it to real physical observables: as was mentioned above, we should distinguish real physical effects from coordinate artifacts.
How can we be sure that the time dependence of the scale factor cannot be similarly removed by a suitable change of coordinates? For example, let us consider the following cosmological model. The metric Eq. (1) has two characteristic length scales both of which in general evolve in time: the spatial curvature radius R c = a/ 6|k| and the Hubble length R H = H −1 , where H =ȧ/a is the Hubble parameter that sets the time scale of the cosmic expansion. An attractive possibility is the perfect cosmological principle that the universe at large is not only homogeneous in space but also in time. Bondi's derisive formulation [20] "Geography does not matter, and history does not matter either" expresses the perfect cosmological principle more eloquently. It is immediately obvious that for the non-dependence on time of the curvature radius and of the Hubble length, it is necessary to have k = 0 and a(t) = a 0 e Ht with some constant a 0 which can be absorbed into the exponential part and then eliminated by a time translation (by suitable redefinition of the origin of time scale). Thus, the metric takes the form: This beautiful theory, based on the perfect cosmological principle, used to be a base for the steady state theory (see, for example, [23]) before the latter being killed by some ugly observational facts [24]. However, nowadays this de Sitter cosmology surprisingly resurrected as a possible description of reality as it may correspond to universe's initial inflationary phase or its terminal state. Let us introduce the new coordinates T, R instead of t, r through relations (see, for example, [10]): In this so called static coordinates the metric Eq. (5) becomes The static coordinates determine another slicing (foliation) of the de Sitter space, which doesn't show any sign of space expansion. What about other foliations?
Let us introduce the following fancy coordinatization of the de Sitter space: x 1 = re Ht sin θ cos φ, x 2 = re Ht sin θ sin φ, x 3 = re Ht cos θ.
This set of five coordinates is obviously not independent as the de Sitter space-time is four dimensional. From Eq. (8) we easily find that Therefore, these coordinates define a four-dimensional hypersurface in the five-dimensional space. The really remarkable fact about these coordinates is however that the de Sitter expanding metric Eq. (5) on this hypersurface is induced by the five-dimensional Minkowski metric in the ambient space: This can be easily checked by explicit calculations starting from relations Eq. (8). But these coordinates are in fact based on the symmetry considerations [25].
Thus, the de Sitter universe can be considered as a highly symmetric four-dimensional hypersurface (a Minkowskian sphere) in circumambient five-dimensional Minkowski space-time. Each foliation, which endows this Minkowskian sphere with a sense of time, are on equal footing according to the spirit of relativity. However, the cosmological circumstances can dictate the preferred foliation. For example, the exponentially expanding universe Eq. (5) corresponds to the sense of time of comoving observers for which the cosmic microwave background is isotropic and the universe itself is spatially flat. However, other foliations exist and can give closed or open universes [25,26]. Namely, if we define comoving FRW coordinates as: x 0 = R H sinh (Ht), x 1 = r cosh (Ht) sin θ cos φ, x 2 = r cosh (Ht) sin θ sin φ, x 3 = r cosh (Ht) cos θ, the ambient Minkowski metric induces on the four-dimensional Minkowskian sphere the metric: which corresponds to the close, cosmological constant dominated, FRW universe. While using the coordinatization x 2 = r sinh (Ht) sin θ sin φ, x 3 = r sinh (Ht) cos θ, we get an open, cosmological constant dominated, FRW universe with the metric As we see, the notion of expanding space is rather subjective. Only comoving observers (galaxies) would find it natural, and our insistence on using this notion reflects our belief that there is no other natural set of observers which would define a different preferred foliation for cosmology. Nevertheless, it seems that quasars dispute this common wisdom and do not find expanding space natural.
Cosmological redshift is usually associated with the expanding space. However, in contrast to the expanding space, redshift is a directly observable objective property of light signals emitted by distant sources.
Suppose two light signals are emitted at cosmological times t e and t e + ∆t e in a distant galaxy at the comoving coordinate ψ. We at ψ = 0 will receive signals at cosmological times t r and t r + ∆t r . How are ∆t e and ∆t r related to each other? We can find this relation as follows [20].
Light propagates along null-geodesics or with ds = 0 along its world-line. Then, for radially propagating light-signals we have from Eq. (4) dt = −a(t)dψ and, therefore, Hence, we obtain the cosmological time dilation relation As its derivation indicates, it is an elementary and general consequence of general relativity and space-time geometry. Usually this time dilation relation is applied to redshifts: if ∆t represents the period (inverse frequency) of light-wave with the wavelength λ ∼ ∆t, then the redshift The derivation of the cosmological time dilation relation Eq. (16) is the most simple in the comoving coordinates where it can be interpreted as a result of the space expansion. But we can still use any other coordinates. Let us consider, for example, the de Sitter space-time in static coordinates where nothing expands or shrinks and rederive Eq. (16).
For radial light propagation, we have from Eq. (7) that (R decreases when the light signal moves towards us from a distant galaxy. Hence the minus sign). We obtain then the following relations: Assuming infinitesimal ∆T r and ∆T e , we can easily get from Eq. (19) The coordinate r of a comoving galaxy does not change with time; taking it into account, the second relation in Eq. (6) indicates that Using this expression and R r = 0, we get from Eq. (20) But the coordinate time ∆T is not the time that clocks of the comoving observers actually measure. The latter according to Eq. (5) is just the cosmic time ∆t, and the first equation in Eq. (6) together with Eq. (21) indicate the following relation between them Therefore, The fact that Eq. (24) is the same relation as Eq. (16) is not obvious, but can be shown as follows [10]. Along a null radial geodesic from (t e , r e ) to (t r , 0), we have according to Eq. (5) dt = −e Ht dr and respectively Then, the second equation in Eq. (6) indicates that and, therefore, Which proves that Eq.s (16) and (24) relations are equivalent. However, the physical interpretation of the cosmological time dilation in static coordinates is different. The relation following from Eq. (7): between the cosmic time t (the proper time of comoving observers) and the coordinate time T , where can be interpreted [10] as resulting from both the gravitational redshift due to "potential" Φ(R), and the kinematical effect (the Doppler shift) due to the Hubble recession velocity dR dT of the light-source. As the discussion above shows, in contrast to the expanding space, which is a subjective notion depending upon the coordinate choice, the cosmological time dilation is an objective and very basic phenomenon with the important observable consequence that the rate of any time variation of the radiation emitted by a distant source with the cosmological redshift z should be proportional to (1 + z) −1 (see Eq.s (16) and (17)).
There is no way to escape this conclusion under the standard belief that all distant astrophysical objects participate in the Hubble flow. Therefore, it is not surprising that observations of distant type Ia supernova confirms this time dilation prediction [27] (it was suggested long ago to use supernova to detect the time dilation [28]).
However, a recent study of over 800 quasar light curves found no evidence of the expected effects of time dilation [2]. This is a very surprising result with no explanation within the conventional cosmological framework (note that some studies report the lack of time dilation signatures also in gamma-ray burst light curves [29,30]. See however [31,32]).
Some non-conventional explanations were proposed in [2], for example, the idea that the expected effects of the quasar time dilation are exactly compensated by an increase in the intrinsic timescale of quasar variability due to the growth of the quasar black holes, or that the observed variations are not intrinsic for quasars but are caused by microlensing events at much lower redshifts. Nonetheless, these explanations do not seem plausible.
We will try to argue in this note that this surprising result, if confirmed without any doubt, will imply a radical change of our cosmological perspective, namely that (at least some) quasars are related to the pre-Big-Bang objects, which we call St. Augustine's objects for reasons explained later. These objects would define a preferred foliation different from the preferred foliation of the Hubble flow comoving observables and this new preferred foliation implies a pre-existing space-time predating the Big Bang.
Perhaps some comments are appropriate here. Application of the cosmological principle assumes that the cosmic time and the corresponding three dimensional space, defined as a set of simultaneous events according to this cosmic time, are well defined concepts in cosmology. It may seem that this is obvious and should be taken for granted. However, Gödel's discovery [33] of cosmological space-time with closed timelike curves passing through every event made it clear that in general, the existence of global cosmic time is not at all guaranteed (the first example of spacetime with closed timelike curves was in fact given by Lanczos [34] in 1923 and then rediscovered by van Stockum [35] in 1937). The existence of cosmic time requires a condition which is stronger than causality. Namely, a space-time admits a cosmic time if and only if it is stably causal [36]. Space-time is stably causal if closed timelike curves are not only absent but they do not appear even after a deformation of the space-time metric which slightly opens up the light cones all over the space-time [37][38][39].
Friedmann-Robertson-Walker metric is stably causal and therefore admits the cosmic time. However, it is the matter content of the universe (cosmic substratum) and motion pattern of this substratum which actually endow the cosmic time and Friedmann-Robertson-Walker metric with the real physical meaning via celebrated Weyl's principle [40,41]. Weyl's principle is an assumption about the nature of the Hubble flow and it states that world lines of the "fundamental particles" of the cosmic substratum (usually assumed to be galaxies or clusters of galaxies) form, after averaging over the peculiar motion, a space-time-filling family (congruence) of non-crossing geodesics converging towards the common past [41]. Friedmann-Robertson-Walker metric holds a privileged position as it is co-moving with the fundamental observers. The above mentioned congruence of world lines is unique in the FRW space-time and thus defines its preferred foliation into spatial hypersurfaces that are ordered in cosmic time.
Unlike the FRW space-time, in the de Sitter space-time there is no unique choice of congruence and that's why the de Sitter space-time can be considered either as a static or as an expanding universe [41].
When we speak about St. Augustine's objects it is in fact assumed that the FRW space-time is only a part of a bigger entity (imagine, for example, an embedding of a spatially flat FRW space-time into the five-dimensional Minkowski space considered at the end of the paper) and that motion pattern of the cosmic substratum in other parts of this encompassing space are quite different from the Weyl congruence defining the Hubble flow.

OTHER QUASAR MYSTERIES
Quasars (Quasi-Stellar Objects) are fascinating objects discovered half a century ago [42,43]. Shortly after their discovery, Zel'dovich [44] and Salpeter [45] independently suggested accretion on a black hole as an energy source of quasars. This hypothesis became universally accepted after the publication of Lynden-Bell [46] showing that the main properties of quasars could indeed be understood in terms of release of gravitational binding energy of matter accreting onto a supermassive black hole. Nevertheless, quasars largely remain the most mysterious and enigmatic objects in the Universe and even the most basic and major questions about them are far from being settled [47][48][49]. Let us briefly outline some problems in quasar research with no definite answers [49].

The Origin of Supermassive Black Holes
If we adopt the standard view that the observed quasar redshifts are of cosmological origin, when the unavoidable consequence is that the quasars are immensely bright: brighter than thousands of the brightest galaxies at low redshifts. At that strong variability of quasars in short time scales indicates that the majority of this immense energy is emitted from a tiny region not exceeding the Solar System in size. Supermassive black holes are indeed the only powerhouses that can be imagined to meet such extreme energy and size requirements.
A simple argument why the central black hole in quasars should be supermassive goes as follows [50][51][52]. Let us assume a steady accretion so that the inward gravitational force −ρdV ∇Φ acting on a small element dV of accreting gas of density ρ in the gravitational potential Φ of the black hole is balanced everywhere by the radiation pressure force d F on this gas element. If the gas is fully ionized hydrogen then radiation pressure force acts mainly on free electrons through Thomson scattering (for protons the scattering cross section is smaller by the factor (m e /m p ) 2 ≈ 2.5 · 10 −8 ). However, due to Coulomb attraction, the electrons drag the protons with them so that the charge separation does not occur). The Thomson cross section σ T gives the effective area for an electron to absorb radiation when it is illuminated. Therefore, if the radiation energy flux is d Σ, each electron absorbs momentum at a rate (d Σ/c)σ T . The radiation pressure force d F equals the rate at which the volume element dV absorbs momentum. In this volume element we have N = ρdV /m p protons and hence, the same amount of electrons, as the gas is assumed to be fully ionized. Therefore, Then the condition of hydrostatic equilibrium, −ρdV ∇Φ + d F = 0, gives the equation Using this relation, the Gauss's theorem and the Poisson equation ∇ 2 Φ = 4πGρ, G being the Newton constant, we can calculate the total luminosity as a surface integral over some closed surface S around the quasar: Therefore, finally, where M is the quasar mass. This so called Eddington luminosity is an estimate of the maximal luminosity a quasar of mass M can have if it is powered by a steady accretion. Substituting σ T ≈ 6.7 · 10 −25 cm 2 and other physical constants, Eq. (30) can be rewritten as follows [50] L E ≈ 1.3 · 10 38 M M ⊙ erg s −1 (31) where M ⊙ is the solar mass. Typical quasar luminosity is L Q ≈ 10 46 erg s −1 [52]. Then Eq. (31)indicates that a central mass of the order of 10 8 M ⊙ is required. Various methods were developed in past decades to estimate quasar masses from observations [53] and the results confirm the above conclusion that quasars are supermassive. Moreover, the results indicate [53] that quasars with masses 10 9 M ⊙ were probably already present at z ∼ 7 when the universe was less than 1 Gyr old. This raises a potential problem: how these supermassive black holes have had grown up given the limited time they had?
If the accretion rate is µ, then the associated luminosity can be written as L = ǫµc 2 , with ǫ representing the radiative efficiency that equals the fraction of the accreted rest energy converted into radiation. The growth rate of the black hole under such accretion isṀ = (1 − ǫ)µ. Therefore, from these two relations, using the Eddington luminosity Eq. (30), we getṀ where is the so called Salpeter time. Equation (32) indicates that under the steady Eddington accretion the black hole mass grows exponentially with time: where M 0 is the mass of the black hole seed. What can we say about the radiative efficiency ǫ? To estimate ǫ, we can use a Newtonian argument. In fact not quite Newtonian, because the pure Newtonian estimate, assuming that the accretion disk extends up to the Schwarzschild radius, is a factor of two higher than the more appropriate general relativistic result. However, it was shown be Paczyńsky and Wiita that the following simple pseudo-Newtonian potential reproduces the correct general relativistic results surprisingly accurately as far as the black hole accretion flaw is concerned [54,55]. Here, r S is the Schwarzschild radius of the black hole. For a test body of mass m moving on a circular orbit of radius r in the potential Eq. (35), we find for its velocity and for its total energy As we see, only for r > 2r S the bound circular orbits are possible. But not all of them are stable. The stability condition is the following: where U is the effective potential corresponding to the orbital momentum l. The orbital radius r is determined from the condition This enables us to exclude l and rewrite the stability condition in the form Therefore, only circular orbits with r > 3r S are stable and hence, the steady accretion disk ends at r = 3r S . The binding energy at this innermost stable orbit is, according to Eq. (36), When the particle of mass m migrates from the outskirts of the accretion disk at r ≫ r S , where its binding energy −E ≈ 0, to the inner edge of the accretion disk, this mc 2 /16 energy is released as heat and eventually radiation. Therefore, the radiative efficiency of such a Schwarzschild black hole engine is ǫ = 1/16 ≈ 6.2%. More correct treatment of general relativistic effects gives ǫ ≈ 5.7% [56]. Rotating black holes can lead to even higher radiative efficiencies: up to about 30% in realistic case [56]. For comparison, nuclear burning of hydrogen to helium in stars has the efficiency of only 0.7%.
Usually, in quasar research one assumes a canonical radiative efficiency of ǫ = 0.1 [53]. For M = 10 9 M ⊙ at z ∼ 7 and t − t 0 ∼ 0.7 Gyr that corresponds to z 0 ∼ 30 when the first stars were formed, we get from Eq. (34) that the black hole seeds of mass M 0 ∼ 10 3 M ⊙ is needed. This seems somewhat high as it is expected that the black hole remnants from the first generation stars had typical values less than tens of solar masses [53].
Although this is not considered as an immediate problem, because we could imagine ways how to produce such massive black hole seeds (for example, direct collapse of primordial gas clouds [53]), but it will be fair to say that we do not completely understand the origin of supermassive black holes at present. Even if we would have a sufficiently massive black hole seed, it is not clear how a steady accretion at Eddington rates can be maintained for such a long period of time. The accreting gas has to lose about 99.9% of its angular momentum to migrate from the host galaxy to the vicinity of the black hole [57]. During this time, various processes can intervene to make the accretion not continuous. For example, star formation which inevitably occur in the accretion disks on large scales, can reduce the amount of gas which can be directly accreted onto the black hole [50,57]. In fact, our present day arguments about how quasars ignite and produce their stupendous energies are of the kind of hand-waving lacking necessary details [58] and so we still do not have good predictive physical models [48].

Evolution of Quasars and Age Problems
Another mysterious thing about quasars which is poorly understood is the evolution of quasar luminosity over redshifts: quasars at high redshifts are far more luminous than quasars at low redshifts [49]. But, strangely, quasar spectra show little evidence, if any, of this evolution: quasar spectral features at low and high redshifts are very similar [49]. It seems that quasars were already mature at z ∼ 6 − 7 and the physical properties of their accretion disks have not changed since then despite strong evolution of their luminosities.
Not only quasars at high redshifts seem to be mature but their surroundings too. It is assumed that it requires many generations of supernovae to produce sizable amount of metals. So, the natural expectation is that the accreting material for high redshift quasars should be metal-free. Observations, on the contrary, indicate that metallicities are even higher than the solar metallicity [49]. Recent discovery of z = 7.51 galaxy with significant metallicity and very high star formation rate [59] also indicates more rapid chemical evolution in the early universe than it is usually assumed.
Some quasars were found [60] at z ∼ 6 without any detectable infrared emission from the hot dust which is expected in the standard accretion disk model and is indeed universally present at low redshifts. This may be a sign of the long evolution of quasars absent in the ultraviolet, optical and X-ray bands of their spectra. However, the rarity of such dust-free quasars indicate a very rapid mass accretion and dust formation at their alleged early evolutionary stages [60].
Another question is what ignites a quasar. It is usually assumed that the major merger events of gas-rich galaxies trigger the accretion on the central black hole and ignite a quasar. However about half of host galaxies have undistorted disks and show no signs of distortions induced by merger [49]. Besides the problem of explaining the appearance at a surprisingly early times after the Big Bang of chemically mature galaxies hosting a supermassive black holes, there is still another age problem related to quasars. Namely, a quasar was discovered at z = 3.91 whose age is 2-3 Gyr, while according to the standard cosmological model the Universe was only 1.63 Gyr old at this redshift [61]. Therefore, this quasar appears to be older than the Universe-an obvious paradox requiring explanation.
While the so called α-process elements (O, N e, M g, Si, S) are mainly produced in type II supernovae, the so called iron-peak elements (V , Cr, M n, F e, Co, N i) are produced in type Ia supernovae [62,63]. The corresponding characteristic time scales are vastly different which enables to estimate the age of astrophysical objects by measuring their F e/α abundance ratios.
Type II supernovae are explosions of massive stars induced by the collapse of their inert iron-nickel cores when the core exceeds the Chandrasekhar limit of about 1.4 solar masses and the electron degeneracy pressure in the core is no longer sufficient to maintain the equilibrium against gravity. The lifetime of such massive stars (M > 8M ⊙ ) is quite small, about 10 6 − 10 7 years.
Type Ia supernovae are thermonuclear explosions of accreting carbon-oxygen white dwarfs in close binaries. The lifetime of such binaries, estimated from the chemical evolution of our galaxy, is about 1 Gyr [63], although a lifetime distribution can span within a wide range of 0.1-20 Gyr [64].
As it was already mentioned, this difference in timescales gives us a specific "clock" to measure the age of astrophysical objects. In the case of the quasar APM 08279+5255 at z = 3.91, using the measured F e/O abundance ratio and a detailed chemodynamical model, it was estimated [65] that the age of this quasar is t Q = 2.1 ± 0.3 Gyr.
The age of the Universe at this redshift can be obtained in the following way. From Eq. (17), with the identifications t e = t, t r = t 0 , where t 0 is the present age, we get On the other hand, the Hubble parameter H =ȧ/a is given by the Friedmann equation which for the flat Universe (k = 0) takes the form (see, for example, [20]): ρ m , ρ r and ρ v being the energy densities of matter, radiation and vacuum, respectively. The vacuum energy density ρ v (the so called dark energy) does not depend on the scale factor. The matter energy density ρ m scales as a −3 and the radiation energy density ρ r as a −4 (because in this case not only the photon number density decreases as a −3 upon space expansion, but also each photon experiences a redshift and its energy decreases as a −1 ). Therefore, we get from Eq. (38), where are the present fractional energy densities of matter and radiation given in terms of the critical density and we have taken into account that for the flat Universe ρ m0 + ρ r0 + ρ v0 = ρ c0 . Equation (39) gives the Hubble parameter as a function of redshift. Therefore, we can separate variables in Eq. (37) and get the following integral (note that at the Big Bang, t = 0, the scale factor is zero and hence the redshift is infinity) [65]: This is an elliptic integral. However, if we ignore the radiation energy density, which is very small at present epoch (Ω r0 < 10 −4 ) and becomes important only at very high redshifts hence contributing little in the integral Eq. (40), then the integral can be calculated in terms of elementary functions by a sequence of substitutions [66] x and we finally get [66,67] For z = 3.91, H 0 = 70.4 kms −1 Mpc −1 and Ω m0 = 0.27, Eq. (41) gives t ≈ 1.63 Gyr and we have the above mentioned age problem t < t Q .

Correlation of Quasars with Galaxies of Lower Redshifts
A long standing controversy is an alleged association of high-redshift quasars with low-redshift galaxies [49,[68][69][70]. Some quasars are surprisingly close, only a few arcseconds away, to the galaxies with widely different redshift. In some amazing cases filaments apparently connecting the pairs of objects with different redshifts have been observed. Maybe one of the most impressive cases is the Seyfert galaxy NGC 7603 [71,72]. Seyfert galaxies, less energetic cousins of quasars, are spiral or irregular galaxies containing an extremely bright quasar-like nucleus. The redshift of NGC 7603 is z = 0.029. There is another smaller galaxy, NGC 7603B, with z = 0.057 nearby and the two are apparently connected by a narrow luminous filament, NGC 7603B lying at the end of the filament. The redshift of the filament itself is z = 0.030 (the same as for NGC 7603 within the measurement errors which are of the order of 0.001 [71]). There are two more galaxies lying exactly on the filament: one with the redshift z = 0.243 positioned near NGC 7603B, and the another one with z = 0.391 positioned near NGC 7603 [71,72] (see Figure 1a in [71]).
Mainstream cosmologists simply ignore this kind of evidence supporting association of objects with different redshifts considering all these cases as a just random projections of background objects like constellations. They might be right in some cases but the statistics of anomalies seems to be too high and indicates that we do not completely understand these phenomena [49].

Apparent Superluminal Motion
The last conundrum related to quasars which we will discuss in this section is the problem of apparent superluminal motion [49,69,74,75]. In fact this phenomenon was predicted [73] before it was discovered following the development of Very Large Base Interferometry (VLBI) that allowed a very accurate mapping of the radio sources morphologies. It was found that in many radio galaxies and quasars there are several compact sources of radio emission that appear to move apart superluminaly in successive VLBI images [76].
Canonical explanation of this apparent superluminal motion (see, for example, [52]) is ingenious but simple enough to be used as an undergraduate problem for the first year physics students [77]. The standard derivation, however, is valid only for nearby sources as it ignores cosmological space expansion, and this fact is not sufficiently stressed in existing literature [78]. Therefore, to avoid a confusion, we provide, following [78], a derivation valid for all redshifts.
Suppose two radio emitting blobs coincide at a point A (see Figure 1) at cosmic time t e . The observer at a point C notices this coincidence at a later time t r because of the finite speed of light. One blob remains stationary while the other moves at a relativistic speed V towards the observer under a small angle θ = ∠BAC to the line-of-sight of the observer. After infinitesimal time ∆t e , the moving blob will reach the point B. The observer notices this event at time t r + ∆t r with some ∆t r . As was already explained, in FRW universe the emission and reception times along radial null-geodesics are related as follows (in units where c = 1, so below we will use V and β = V /c interchangeably) where ψ is the co-moving radial coordinate of the point A, ψ + ∆ψ-of the point B and it is assumed that the observer is at the origin ψ = 0. From Eq. (42) we have On the other hand, the point B has, at the first order in ∆t e , the same radial coordinate as the point D. Therefore, Combining Eq.s (43) and (44) we get, after using Eq. (17), Therefore, the observed apparent transverse velocity of the moving blob with respect to the stationary one is [78] The quantity which is actually measured on Earth is the observed angular proper motion of the moving blob on the sky per time unit µ = ∆α/∆t r and to transform it into the apparent transverse velocity V a = µD A we need to know the so called angular diameter distance D A . It is defined in such a way [20] that for a small solid angle ∆Ω = 2π(1 − cos (∆α)) ≈ π(∆α) 2 the corresponding cross-sectional area ∆S at a distance D A is On the other hand, for the FRW metric Eq. (4), the total area of a sphere ψ = const at cosmic time t e is 4πa 2 (t e )S 2 (ψ; k) and this area corresponds to the solid angle 4π. Therefore, From Eq.s (46) and (47), we get D A = a(t e )S(ψ; k). In particular, for the spatially flat universe, To express Eq. (48) in terms of observed quantities, we note again that for radial null-geodesics dψ = −dt/a(t) and therefore, by using Eq.s (17) and (39), we get [78] As we see, the apparent transverse velocity can be calculated from the observed quantities as follows To concentrate on the intrinsic parameters of the source, it is a common practice to eliminate the redshift dependence of the apparent transverse velocity by considering (we have recovered the light velocity c for a moment) Rather misleadingly,V a is widely called "apparent velocity" despite the fact that it is the apparent velocity only for a co-moving observer near the source [78]. The true apparent velocity for a terrestrial observer at a cosmological distance from the source is V a , as our derivation above indicates, and it is suppressed by a factor 1 + z [78].
As a function of the angle θ,V a has a maximumV am = γV when θ = θ m , where cos θ m = β. As we see,V am can be arbitrarily large and at first sight relativistic beaming model gives a natural solution of the apparent superluminal motion problem. However, matters are not as simple as they seem.
We need γ ≫ 1, if we wantV am ≫ c. But then β ∼ 1 and θ m is very small. If quasar jets are isotropically distributed in space, then the probability of observing an approaching jet with θ < 10 • is (assuming that the quasar jets are double-sided and symmetrical) p = 2π(1 − cos θ)/(2π) ≈ θ 2 /2 ≈ 0.015. This quite small probability indicates that the observation by chance of an apparent superluminal motion in quasars should be relatively rare phenomenon. It was clear from the beginning that real observations do not confirm this conclusion [79]. In fact apparent superluminal motion is common and ubiquitous in quasars [69,80]. For example, the MOJAVE (MOnitoring of Jets in AGN with VLBA Experiments) survey contains 101 quasars with 354 observed radio components, 95% of which move with apparent superluminal velocities with respect to the core. The maximum observed velocity equals to 60c and about 45% of radio emitting blobs have velocities larger than 10c [80].
A close inspection of the formula (50) indicates that it requires roughly θ < 10 • to obtainV a > 10c. Then randomly, using the above mentioned probability p ≈ 0.015, one would expect only 354 * p ≈ 5 blobs, out of 354, to haveV a > 10c. In reality the observed number is 158-a totally improbable event by chance [80].
At this point, the relativistic beaming model is clearly in trouble. Fortunately for the standard paradigm, the Doppler boosting which is intrinsically related to the relativistic beaming comes at rescue. In particular, the observed flux density, F ν , from an emitting blob which moves with a relativistic velocity towards observer making an angle θ with the line of sight is enhanced compared to the intrinsic flux from the blob at rest, F 0,ν , according to the formula [79,81] where α is the spectral index of the radiating source (in other words, it is assumed that F ν ∼ ν −α ). This relation follows from the relativistic invariance of I ν /ν 3 and the Doppler blue-shift formula for the frequency ν 0 = γ(1 − β cos θ)ν [82]. Here I ν is the specific intensity related to the flux density as follows F ν = I ν cos θ ′ dΩ ′ [83]. Indeed, since I ν /ν 3 = I 0,ν0 /ν 3 0 , we have But (k 0 is some constant) and, therefore, which is equivalent to Eq. (51).
To prove that I ν /ν 3 is indeed relativistic invariant, we first show, following [83] (see also [84]), that a phase space volume element d 3 r d 3 p is Lorentz invariant. Let a group of particles occupy an infinitesimal volume d 3 r around some radius-vector and have a slight spread d p of momenta around a mean momentum in some inertial reference frame S. In the frame S ′ , comoving with the particles, the mean momentum equals to zero and all particles have the same energy mc 2 to the first order because the infinitesimally small momenta of individual particles contribute quadratically to their energy and therefore this contribution vanishes to the first order. Let us imagine some particle from the group with a momentum p ′ in the frame S ′ . Then in the frame S we have and therefore the momentum spreads in the frames S and S ′ are related as follows dp x = γ dp ′ x , dp y = dp ′ y , dp z = dp ′ z On the other hand, due to Lorentz contraction, There is one subtlety here which sometimes can lead to confusion (see, for example, [85,86]). We need particles to occupy a given volume d 3 r element simultaneously. But simultaneity is relative: if a group of particles are all inside a volume element d 3 r at time t in reference frame S, in the comoving frame S ′ they are inside the corresponding volume element d 3 r ′ = d 3 r/γ at different times and according to the Lorentz transformation the spread of these times around the mean time t ′ is −γβdx/c. In the frame S ′ the particles are not in general at rest but move (in the x-direction) with infinitesimal velocities of the order of dp ′ x /m = dp x /(mγ). Therefore, to bring all particles at the time instant t ′ will require a shift of their x-coordinates of the order of ∆x = γβ dx c dp x mγ = βdx dp x mc Analogously, the required shifts in the y and z directions are ∆y = βγdx dp y mc , ∆z = βγdx dp z mc Fortunately, for finite m and γ, these quantities are of the second order and hence, can safely be ignored. Let us stress, however, that for this argument to be valid it is essential for S ′ to be the instantaneous comoving frame because only in this frame particle velocities are infinitesimal. If S ′ is not the comoving frame, we can no longer in general ignore the effects of desynchronization, and if we do, we get an erroneous conclusion, as in [85], that the phase space volume element is not Lorentz invariant.
Relations Eq.s (52) and (53) do indeed show that dx dy dz dp x dp y dp z = dx ′ dy ′ dz ′ dp ′ x dp ′ y dp ′ z Of course this reasoning is strictly valid only for particles of finite mass. However, the final result Eq. (56) has no reference altogether to the particle mass and it is reasonable to propose that it should be true also in the m → 0 limiting case of photons. Note, however, that Eq.s (54) and (55) indicate that the m → 0 limit is rather subtle. More formal and rigorous derivations of Eq. (56) for photons can be found in [87,88]. Now let us consider an energy content dE of a phase space volume element d 3 r d 3 p occupied by dN monochromatic photons of frequency ν. It can be calculated in two different ways. First, from the definition of specific intensity, we have dE = I ν dν dΩ dt dA, where dA is the cross-sectional area of d 3 r in the direction perpendicular to the photons mean momentum. On the other hand, where we have used the relation d 3 r = cdt dA and f = dN d 3 r d 3 p is the photon distribution function which is Lorentzinvariant because both dN and d 3 r d 3 p are Lorentz-invariant. Comparing these two expressions for dE, we see that [88] This completes the proof that I ν /ν 3 is Lorentz-invariant. According to Eq. (51), relativistic beaming has a very strong effect on the observed luminosity. For γ ∼ 10, α ∼ 0 and θ ∼ 10 • , the enhancement is about two orders of magnitude (this is true for a simple ballistic model of one radiating approaching blob for which Eq. (51) applies. In the case of more realistic jet geometries the enhancement factor may be smaller [89]). This Doppler boosting effect allows to explain the above mentioned huge discrepancy in the statistical properties of observed superluminal sources as a kind of the Malmquist bias [80]. Simply, due to the Doppler boosting, it is natural to expect that flux-limited surveys of radio emitting sources will indeed preferentially pick up brighter sources that are ejecting material with small θ. However, this argument assumes that the main part of the observed flux, even the one that comes from the core, originates from relativistic jets, and not everybody is ready to except this assumption [90].
To conclude this section, at present we have no compelling reason to reject the standard cosmological concepts about quasars. However, it seems that there is a growing number of alleged observational facts which are difficult to reconcile with the standard paradigm. It is not a good strategy to simply ignore these facts, in our opinion. It would be much better to remember the nineteenth century humorist Josh Billings' aphorism "It ain't ignorance causes so much trouble; it's folks knowing so much that ain't so" [91] and keep eyes wide open to new insights.

THE MILNE MODEL
Let us consider a FRW universe with vanishing total energy density ρ and vanishing cosmological constant Λ. The Friedmann equation [20] then gives (for not static universe) k < 0 andȧ = √ −k. By rescaling the radial coordinate and time, we may assume that the units are such that the curvature constant is k = −1 together with c = 1. Therefore, in this case the scale factor a(t) = t linearly grows with time. The possible integration constant can always be eliminated by a suitable change of the time origin. The metric Eq. (4) takes the form Such a space-time is called the Milne universe. It was Milne who realized that such a model must have an alternative, special relativistic description because it has the gravity "turned off" (vanishing energy density) and no cosmological constant [92][93][94][95]. Indeed, if we introduce new coordinates the metric Eq. (57) takes the form which is just the Minkowski line element in spherical coordinates. Let us note that the Milne coordinates t and ψ cover only a part of the Minkowski space-time. Indeed, equations Eq. (58) indicate that T > 0 (because the cosmic time t > 0, t = 0 corresponding to the Big Bang) and R < T . Therefore, the Milne universe corresponds to the future light cone of the Big Bang event. To cover the past light cone, we can redefine the Milne coordinates ψ ≥ 0, t ≥ 0 in such a way that This leads to the same metric Eq. (57) with a singularity at t = 0. For the expanding Milne universe Eq. (58), the singularity corresponds to the Big Bang, and for the contracting Milne universe (with −t ≤ 0 as a cosmic time), the singularity corresponds to the Big Crunch.
Equation (58) give R/T = tanh ψ, which for constant ψ represents a world-line of an observer uniformly moving radially with rapidity ψ. This clarifies the physical meaning of the Rindler coordinate ψ and suggests the following interesting interpretation of the Milne universe [96], originally due to Milne [92][93][94][95].
In the Minkowski space there is a point-like explosive event, the Big Bang, resulting in infinitely many debris of point particles (fundamental observers) shot out at all speeds up to the speed of light in all directions. It is possible to arrange a special kind of velocity distribution of fundamental observers for which the cosmological principle will remain true: all fundamental observers will see the identical velocity distributions. Indeed, let us show that this is possible [93]. Let f (u) be such a distribution. Because of the assumed isotropy of space, this distribution can depend only on the magnitude u of the velocity u, not on its direction. The distribution function f (u) does not depend on time because after the Big Bang no new fundamental particles are created or destroyed and each of them moves with a constant velocity. The cosmological principle demands where (according to our conventions, the light velocity c = 1) is given by the relativistic velocity addition formula (see, for example, [97]). The Jacobian of this transformation is (to avoid involved calculations by hand, REDUCE Computer Algebra System [98] can be used to get this result) where at the last step we have used the so called gamma identity [99,100] Therefore, Eq. (60) takes the form and we see that (we have restored the light velocity c in the last formula for a moment) with some positive constant B. This velocity distribution function has several remarkable properties [93]. The total number of particles N = f (u) du x du y du z is infinite and the velocity-centroid (the mean velocity) < u >= 1 N u f (u) du x du y du z cannot be defined. The density of particles (in velocity-space) increases and tends to infinity as the velocity u approaches the light velocity. Each fundamental observer can be regarded as the stationary center of an expanding ball which looks one and the same for every choice of the fundamental observer. Not every explosion readily produces such kind of distribution. In fact, "it must have been an incredibly delicately tuned Big Bang to achieve this!" [96].
The Milne distribution Eq. (62) is not of statistical but of hydrodynamic nature: the velocity u of every fundamental particle at time T is uniquely determined by its position through the Hubble-like law This fact allows to rewrite Eq. (62) as a spatiotemporal distribution with the particle density function The density tends to infinity at the unattainable boundary R = T of the expanding ball which moves away at the light velocity.
The density Eq. (63) is inhomogeneous and therefore the Copernican principle is not evident in the Minkowski coordinates R and T . Of course, indirectly it is implemented through the equivalence of fundamental observers (cosmological principle). However, for the Milne universe there exists a natural foliation of space-time under which the spatial distribution of matter becomes homogeneous and thus the Copernican principle becomes obvious.
The symmetry group of special relativity (Minkowski space-time) is the ten parameter Poincaré group of inhomogeneous (including space-time translations) Lorentz transformations. The symmetry group of the Milne universe is a six parameter subgroup of the Poincaré group of homogeneous Lorentz transformations which leave invariant a special event-the Big Bang. As a result, there is a preferred class of inertial observers, namely those whose world-lines pass through the Big Bang event. The Milne universe is a description of the Minkowski space-time from the point of view of members of this preferred class [101]. The invariant proper time after the Big Bang measured by the clocks of this fundamental observers, can be used for a natural definition of simultaneity t = const in the Milne universe, alternative to the usual Einsteinian simultaneity T = const (interestingly, Einstein's definition of simultaneity by means of the light signal exchanges coincides, in principle, with the St. Augustine's criterion of simultaneity given by him 1500 years earlier [102]). This cosmic time t leads to the slicing (of a part) of the Minkowski space-time into space-like hypersurfaces Eq. (64) (two-sheeted hyperboloids T 2 − X 2 − Y 2 − Z 2 = t 2 ) which Milne calls the public space to distinguish it from the private space of fundamental observers defined through the T = const slicing, T being the usual Minkowski time. Note that immediately after the Big Bang (for every t > 0) the public space is infinite, while the private space for every T > 0 is a finite sphere of radius R = T , increasing in size with light velocity as T goes on. Thus, the Milne model provides an excellent demonstration of a seemingly paradoxical situation that a point-like Big Bang can give birth to a universe with infinite spatial extension [96,103]. Interestingly, the Milne foliation of Minkowski space-time is exactly the foliation which leads to the Dirac's point-form of relativistic dynamics [104,105].
Letñ(ψ, t) be the density of particles in public space so that the number of particles in the volume element dV = |g| dψ dθ dφ at cosmic time t is dN =ñ(ψ, t)dV . For the metric Eq. (57), the spatial components of the metric tensor areg ψψ = t 2 ,g θθ = t 2 sinh 2 ψ,g φφ = t 2 sinh 2 ψ sin 2 θ Therefore, the square root of the absolute value of the determinant of the metric tensor equals to |g| = t 3 sinh 2 ψ sin θ and dN =ñ(ψ, t) t 3 sinh 2 ψ sin θ dψ dθ dφ On the other hand, let us calculate the same number of particles in the private space. For fixed cosmic time t, the radial coordinates lay in-between R = t sinh ψ and R + dR for the selected bunch of particles in the private space, where dR = t cosh ψ dψ. But, according to Eq. (58), this particles have different Minkowski time coordinates T = t cosh ψ. Namely, the particles with the radial coordinate R + dR have the Milne coordinate ψ + dψ and hence the T -time coordinate T + dT with dT = t sinh ψ dψ. However, in private space the fundamental particles move radially outward with velocity V = tanh ψ. Therefore, at time T the bunch of particles have a reduced radial spread dR = dR − V dT = t cosh ψ dψ − t tanh ψ sinh ψ dψ = t cosh ψ dψ and hence occupy dV = R 2 sin θdR dθ dφ volume element in the private space. The corresponding number of particles is From Eq.s (65) and (66), we getñ where the last step follows from equations Eq.s (58), (63) and (64). As we see, the densityñ depends only on cosmic time t and not on the position in the public space-the distribution of fundamental particles in the public space is not only isotropic but also homogeneous and the Copernican principle becomes explicit. In more formal way, the same result can be obtained as follows [101]. The density-flow four-vector where the density n is given by Eq. (63), transforms under the change of coordinates Differentiating Other portion of the partial derivatives can be obtained by differentiating R/T = tanh ψ. The results are and correspondingly we get from Eq. (69) j ψ = 0.
Analogously, differentiating tan 2 θ = (X 2 + Y 2 )/Z 2 , we get Substituting these partial derivatives in Eq. (69), we find that j θ = 0 and j φ = 0. Therefore, in the Milne universe the density-flow four-vector has the form We see that matter (fundamental particles) is at rest and uniformly distributed in the Milne Universe. As was already mentioned, the Milne distribution Eq. (62) does not allow to determine the mean velocity of the whole system of fundamental particles and therefore the global standard of rest. However, using the mean motion of matter in the limited portion of space-time, we can define a local standard of rest in the neighborhood of a given point by defining rest and velocity relative to this mean motion. Introduction of the cosmic time and public space (co-moving coordinates) is just a realization of this strategy [101].
In the following, for simplicity, we will use a two-dimensional space-time to illustrate our arguments. The twodimensional Milne metric has the form (t ≥ 0, −∞ < ψ < ∞) How can imaginary intelligent observers in the Milne universe discover that their space-time is just a part of the more vast Minkowski space-time? The coordinate transformation Eq. (58) is not simple enough to be trivially guessed. Instead we can envisage the following line of thoughts (adapted from [105]).
Null geodesics (light rays) can be used to construct a natural coordinate grid in the Milne universe. Coordinates with respect to this grid (null coordinates) are expected to convey and reveal the intrinsic geometry of the Milne space-time. According to Eq. (72), null geodesics are determined by the condition Integrating, we get that along null geodesics where g > 0 is an integration constant such that gt is dimensionless. Accordingly, we introduce null coordinates The inverse transformation is and after substituting it into Eq. (72), we get The form of Eq. (74) suggests to make further transformation which brings the metric into the very simple form ds 2 = − 1 g 2 dU dV (76) Equation (75) indicate that the Milne space-time corresponds to the coordinate ranges U < 0 and V > 0. However the metric Eq. (76) is no longer singular and we can extend our space-time by allowing the whole ranges −∞ < U < ∞ and −∞ < V < ∞. The identification puts the metric in the canonical Minkowskian form ds 2 = dT 2 − dX 2 and thus shows that our extended space-time is just the Minkowski space-time. The coordinate transformation, inverse to Eq. (77), is which implies the following relations between the Minkowski coordinates (T, X) and the Milne coordinates (t, ψ) (compare with Eq. (58)) If we choose instead U = e −u , V = −e v , we get the contacting Milne universe with which correspond to the U > 0, V < 0 patch of the Minkowski space-time. What about U > 0, V > 0 and U < 0, V < 0 patches of the extended space-time? Let us take Substituting into Eq. (76), we get Therefore, α is timelike and β is spacelike. Thus, we identify α = gt R = ψ, β = gx R which brings the metric into the form Such a metric is known as Rindler metric [20]. According to Eq. (78), the Minkowski coordinates (T, X) and Rindler coordinates (t R , x R ) are related as follows It is clear from Eq. (81) that we must have x R > 0. Therefore the Rindler space-time corresponds to the X > 0, |T | < X portion of the Minkowski space-time (right Rindler wedge). The choice U = −e −ũ , V = −eṽ leads to the left Rindler wedge with the same metric Eq. (82) but with Figure 2 summarizes our findings. It is interesting to note that the above discussion of the Milne (or Rindler) spacetime extension closely resembles the celebrated Kruskal extension of the Schwarzschild space-time [20,106,107]. More precisely, the expanding Milne universe corresponds to the black hole region of the Kruskal diagram, the contracting Milne universe corresponds to the white hole region, and the Rindler space-time corresponds to outer Schwarzschild regions [20]. What is the physical meaning of the Rindler coordinates? Inspired by the Milne interpretation outlined above, let us consider a set of Rindler fundamental observers which are at rest in the Rindler space-time having constant x R spatial coordinates. The line of simultaneity for an arbitrary event P = (T, X) at the worldline of a Rindler observer is the set of events that are simultaneous with P in the inertial instantaneous rest frame of the Rindler observer at P . As x R is fixed, this instantaneous rest frame moves with the the Minkowski velocity (in the private space of the Big Bang event) Therefore, according to the Lorentz transformations, in the instantaneous rest frame the event P has the time coordinate T ′ = γ v (T − vX) = 0, what means that it is simultaneous with the Big Bang event. As both the Rindler observer and the point on its worldline were arbitrary, we come to an amusing conclusion that the Big Bang event is simultaneous with every event on the worldline of every Rindler observer. In other words, the space axes of the fundamental Rindler observers all pass through the Big Bang event. In comparison, the defining characteristic property of the Milne observers is that all their time axes pass through the Big Bang event. This interchange of space and time makes a big difference: the Rindler observers cannot be inertial because the space axes of the inertial observer are all parallel to each other and thus cannot cross at the origin.
From Eq. (84) the acceleration of the Rindler observer is On the other hand, differentiating the velocity addition formula For the instantaneous rest frame V = v, so Eq.s (85) and (86) will give As we see from Eq. (87), the Rindler observers in the private space of the Big Bang event move with the constant proper acceleration which is inversely proportional to their Rindler coordinate x R and approaches the infinity for observers near the Big Bang, x R → 0.
Equation (84) shows that at Rindler time t R all Rindler observers have the same rapidity ψ = gt R . In fact the rapidity is the most natural (dimensionless) time coordinate in the Rindler space-time. The Rindler metric expressed in terms of the rapidity ψ does not contain an arbitrary constant g. Similarly, in the Milne space-time the rapidity is the most natural (dimensionless) spatial coordinate.
What is the origin of the constant g? The fact that Rindler observers share their lines of simultaneity enables us to define global synchronization in the Rindler space-time: we can simply pick out an arbitrary Rindler observer with the spatial coordinate x R0 and declare its proper time to be the coordinate time t R . But the proper time of the observer with the spatial coordinate x R is Comparing with Eq. (83), we see that We have chosen to define synchronization τ = t R for the selected observer, and hence, gx R0 = 1, so that g is the proper acceleration of this observer. Therefore, the freedom to use an arbitrary constant g = 1/x R0 in description of the Rindler space-time is the freedom to select any Rindler observer to define global synchronization and hence, the coordinate time t R . For Rindler observers other than the selected one, the proper time will not coincide with the coordinate time but will be proportional to it, as Eq. (88) shows.
As we see, the Rindler space-time is a part of the Minkowski space-time as viewed by a set of non-inertial observers whose constant proper accelerations are adjusted in such a way that their lines of simultaneity all pass through one event (Big Bang). Further details and a lucid discussion of the physical meaning of Rindler coordinates can be found in [108].
Note that in the four-dimensional space-time radially moving Rindler observers define spherical Rindler coordinates (t R , r R , θ, φ) as follows [109] T = r R sinh (gt R ), X = r R cosh (gt R ) sin θ cos φ Y = r R cosh (gt R ) sin θ sin φ, Z = r R cosh (gt R ) cos θ They cover the whole external part of the Big Bang's light-cone (subdivision into Rindler wedges arises only in the (1 + 1)-dimensional space-time) where the induced metric is

ST. AUGUSTINE'S OBJECTS
Both Milne and Rindler fundamental observers may not suspect that their space-times are just some parts of a bigger entity -the Minkowski space-time. However, we know that this is the case. A legitimate question then is whether Milne and Rindler parts of the Minkowski space-time can exchange material objects. As far as the space-time itself is concerned, the only restriction is causality, because the Minkowski space-time is not singular. Thus, naively, material objects from other parts of the Minkowski space-time can penetrate expanding Milne universe. One timelike world line of this kind is shown on Figure 3. We will call such objects, maybe rather frivolously, the St. Augustine's objects.
St. Augustine of Hippo was the most influential philosopher who pondered on the nature of time. "What then is time? If no one asks me, I know what it is. If I wish to explain it to him who asks, I do not know" -ruminates he in his Confessions [110]. St. Augustine's version of the modern question what was going on before the big bang sounds like this: "What was God doing before He made heaven and earth?". St. Augustine overcomes his own temptation to laugh at the question and give a jesting answer: "He was preparing hell for those who pry into mysteries." Instead he comes to an answer which strikingly resonates with the modern day ideas about space-time and that the universe was created by quantum tunneling from a state with no classical space time [111,112]. Namely, in City of God [113] he comes to the conclusion that "time does not exist without some movement and transition" and this is what distinguishes eternity and time. His main inference then is that "the world was made, not in time, but simultaneously with time". In modern words, time itself was created by the Big Bang and the question what was before the Big Bang is meaningless, it is like asking what is to the north to the North Pole [114].
However St. Augustine's brilliant insight that the time was created simultaneously with the material world should be taken with a grain of salt, because the use of the word "create" may imply that there was some concept of time even at the epoch when time was still not existing [115,116]. The Milne model provides a ready resolution of this dilemma.
As we know, there are two notions of time in the Milne model. The public time is determined by an expanding substratum of fundamental particles. Space-like hypersurfaces orthogonal to the world lines of these fundamental particles define a public space. The public time is the proper time along world lines of the fundamental particles and this cosmic time in a sense labels the series of the space-like hypersurfaces. Therefore, the public time is created by the Big Bang and it is really meaningless to ask what happened before the public time zero.
However, to any inertial fundamental particle in the Milne model we can attach the usual Minkowski frame and thus define the private time of this observer. If there are fundamental observers whose world lines begin not at the Big Bang but in the preceding contracting phase, we can use the private time of such observers supposing that they survived during the Big Bang, to make the question of what happened before the Big Bang obviously meaningful.
The central question, therefore, is whether material objects can cross the boundaries (the future light-cone of the Big Bang event) of the expanding Milne universe. This question is not a trivial one, because the density of the fundamental observers, as we have seen above, becomes infinite at these boundaries, as a result of the cosmological principle.
Milne himself was well aware of this problem. However, he considers the question what is outside the sphere R = T in the private space as meaningless based on the following arguments [93].
Milne observers inside the sphere R = T cannot observe any object outside, because infinite density of fundamental particles on this sphere serves as a curtain obscuring the vision and precludes any window into outer space in any direction. Therefore, only objects which were overtaken by the expanding frontier and appeared inside the sphere R = T can be observed. For inhabitants of the Milne universe the St. Augustine objects would suddenly appear to observation in an act of creation. More precisely, from the point of view of public time, such objects are just a part of the Big Bang initial conditions. However, Milne thinks that the infinite density of fundamental particles on the expanding frontier can lead to the infinite number of collisions between the St. Augustine object and the fundamental particles, thus, not allowing the St. Augustine object to penetrate the boundary. In fact, not the infinite density by itself is crucial in this argument, but a particular form of the density function Eq. (63), which implies that the total number of particles in any thin layer including the boundary is infinite. Milne concludes: "such objects never come, and never can come, into interaction with observable members of the given system. They may, mathematically described, tend to be swept up by the expanding frontier, but they never penetrate it. Penetration would be a logical self-inconsistency of description of the whole system" [93].
But can the ideal Milne universe actually be realized? The infinite number of needed fundamental particles and the initial (in cosmic time) singularity of infinite density, implied by Eq. (67), speak against such a possibility. Another argument against the realizability of the ideal Milne universe can be obtained by considering quantum fields in the Milne universe.
In analogy with Eq. (78), let us define the dimensionfull temporal and spatial (t M , x M ) Milne coordinates: They are related to the Minkowski coordinates (T, X) through: The Milne metric Eq. (72) in terms of these coordinates takes conformally flat form In contrast to the cosmic time t whose range is 0 ≤ t < ∞, the range of the conformal Milne time t M (as well as of the spatial coordinate x M ) is the whole real axis from −∞ to ∞.
Free scalar field Φ(t M , x M ) with mass m = 0 in the Milne universe satisfies the covariant Klein-Gordon equation where the covariant d'Alembertian (Laplace-Beltrami operator) has the form For scalar Φ and vector V µ = g µν ∂ ν Φ fields the covariant derivative ∇ µ is defined as follows [20] where the connection coefficients (Christoffel symbols) are and For the metric Eq. (92), the only non-zero components of the metric tensor are g 00 = e 2gtM , g 11 = −e 2gtM , g 00 = e −2gtM , g 11 = −e −2gtM Then Γ µ µ0 = 2g, Γ µ µ1 = 0 and it can be easily checked that the second term in Eq. (95) equals to zero. Finally, and the Klein-Gordon Eq. (93) takes the form We can separate variables, taking Φ(t M , x M ) = φ(x M )ψ(t M ), and geẗ This equation is formally identical to the s-wave radial Schrödinger equation for an exponential potential with well known solution [117]. Namely, in terms of the cosmic time t = 1 g e gtM , Eq. (97) takes the form which is Bessel's equation (of pure imaginary order) with the general solution [117] where ν = |k|/g and c 1 , c 2 are arbitrary constants. For small arguments [118], Therefore, to get the positive-frequency mode (which behaves as e −i|k|tM ) in the asymptotic past (in Milne conformal time t M ), we must take c 1 = 0. The second constant, c 2 , can be fixed by normalizing the mode using the covariant norm [119,120] where g = det(g µν ) is the determinant of the metric tensor, Σ is a spacelike hypersurface with n µ being the unit timelike vector normal to it, and √ −g dΣ is the invariant (proper) volume element in this hypersurface. Using the Gauss' theorem for curved manifolds [121], it can be proved [119,120] that the scalar product is independent of the choice of spacelike surface Σ provided Ψ 1 and Ψ 2 are solutions of the Klein-Gordon Eq. (93) that vanish sufficiently quickly at spatial infinity. Indeed, let V be the four-volume whose boundary ∂V consists of non-intersecting spacelike hypersurfaces Σ 1 , Σ 2 and, possibly, by timelike boundaries at spatial infinity on which Ψ 1 = Ψ 2 = 0. Then we have But as Ψ 1 and Ψ 2 are solutions of the Klein-Gordon Eq. (93), then: Therefore, Let us specify the scalar product Eq. (100) for our case of two-dimensional Milne space-time with the metric Eq. (92). The convenient choice of Σ is the spacelike hyperbola defined by the condition t M = const. Then the proper line element in Σ is √ −g dΣ = t dψ = e gtM dx M , and the unit vector (normalized according to the metric Eq. (92)) is n µ = (e −gtM , 0). Therefore, the scalar product takes the form Now we can determine the constant c 2 . In fact, using J * iν (mt) = J −iν (mt), the Wronskian relation [118] W and the normalization (for the positive-frequency modes) we find that the normalized mode, which behaves as of pure positive frequency in the asymptotic past, has the form In the asymptotic future, t M → ∞, the cosmic time t also tends to infinity and we should use the following asymptotic behavior [118]: when |z| → ∞, | arg z| < π, then Using Eq. (103), it can be verified that the combinations Eq. (99) that correspond to the positive and negative frequency behavior (∼ e ∓imt ) in asymptotic future are, respectively, the Hankel functions H (2) iν (mt) and H (1) iν (mt), defined by Using the Wronskian relation [118] W we find that the normalized modes, corresponding to the positive and negative frequencies in the asymptotic future, have the form Note that the negative frequency mode has the negative norm: Equations (102), (104) and (105) imply the relation which shows that the pure positive-frequency mode Ψ k,out of the asymptotic future [122]. The standard interpretation [119] of this fact is that the time-dependent background metric leads to a pair production with the averaged number density of produced pairs in the k-mode The problem, however, is that, if we consider the Milne universe as just a re-parametrization of the part of Minkowski space-time, the time dependence of the background metric is just a coordinate artifact, not related to the real presence of a variable gravitational field, and, hence, no particle production is expected [122]. The crux of the problem can be traced back to the definition of the initial vacuum state [123].
Suppose we have a nonzero rate R of particle production due to time-variable background metric. To measure precisely the number of particles in some volume, the measurement duration ∆t must be sufficiently small, namely such that |R|∆t ≪ 1. However, the smaller is the measurement time ∆t, the greater is the uncertainty in energy ∆E ∼ 1/∆t, and due to the possibility of temporary creation of virtual particle pairs, the number of particles will become uncertain by the amount ∆E/m ∼ 1/(m∆t). Therefore, over a time interval ∆t, the total uncertainty in the number of particles is [119,124] The optimum is reached at ∆t = 1 m|R| which gives the following minimum inherent uncertainty of the number of particles: This simple computation indicates that in an adiabatic region, where the Hubble parameter H =ȧ/a is small (the scale factor a(t) varies slowly and hence |R| is expected to be small), the concept of particles, and hence of the vacuum state, is a well defined concept. For the Milne universe H = 1/t. Therefore, the vacuum state is well defined in remote future. However, the region of the Milne space-time near the initial singularity is not an adiabatic region and we do not have a natural well defined vacuum state there. We can use (as we have, in fact, done above) the conformal time t M and quantum states that contain only positive frequencies with respect to this conformal time, when t M → −∞, to define the initial vacuum state in this region. Such conformal vacuum state is, of course, not adiabatic and its use in the role of the initial vacuum state may seem as a somewhat arbitrary choice [123]. Nevertheless, this choice is justified in the path-integral formulation, where it corresponds to the natural restriction of the paths summed to those that lie inside the future light-cone of the initial singularity (the Big Bang event) [125].
However, this restriction is natural only when the initial singularity is real and not a coordinate artifact, as in the case of the naive interpretation of the Milne model. When x > 0, we have the following integral representations of the Hankel functions [126]: Using these integral representations in Eq. (105), we get After changing the integration variable u in the above integrals respectively toψ = ±(gx m − u) = ±(ψ − u), Eq.s (109) take the form where E = m coshψ = √ m 2 + P 2 and P = m sinhψ. As we see, Ψ k,out is a superposition of positive frequency Minkowski plane waves with different rapiditiesψ, while Ψ (−) k,out is a superposition of negative frequency Minkowski plane waves [119,127]. Of course, instead of the rapidityψ, we can use the momentum P = m sinhψ as an integration variable [123]: The integral representations Eq. k,out on the whole Minkowski space-time [123,127] and, of course, these analytic continuations pick up a privileged vacuum state corresponding to no particle production. If we consider the Milne universe as a part of Minkowski space-time, this choice of the vacuum state is the most natural. However, from the point of view of the inside observers of the Milne universe, it may appear rather contrived and will require subtle arguments to justify it from the inside perspective [123].
An other possibility is given by considering the Milne universe as a limiting case of some nonsingular space-time whose metric differs from the Milne-Rindler metric only in the narrow transition region around the Big-Bang light-cone [128]. In this case the particle production can be related to the nonzero time-variable curvature of the space-time (a real time dependent gravitational field) in the transition region. One can imagine that this effect of particle production can take place even in the limit of zero width of the transition region, however such a situation will require a physically unrealistic sources of infinite power on the Milne universe singular boundary [128]. Of course, this picture is in accord with the infinite density of the Milne observers at the boundary, but these infinities make the construction of the ideal Milne universe with impenetrable boundaries (for example, in hydrodynamical laboratory experiments to mimic an arbitrary FRW space-time by relativistic acoustic geometry [129]) impossible. Therefore, any realistic approximation to the ideal Milne universe is expected to have boundaries penetrable for St. Augustine's objects. Now we assume that the same is true to our real universe: the initial singularity is only an unrealistic idealization and in realistic settings it is traversable for St. Augustine's objects. As we expect that the initial density at the Big Bang was nevertheless very high, we assume that only black holes as St. Augustine's objects can survive such a dramatic event. Therefore, we conclude that it is possible to have some amount of black hole population, moving relativistically with respect to the nearby Hubble flow at high redshifts, as a part of the Big Bang initial conditions. If the amount of such St. Augustine's objects is small, they cannot significantly change the dynamics of the universe and spoil the successful predictions of the standard cosmological model.
Let us sketch, from the Milne model perspective, how St. Augustine's objects in the role of quasars can offer an explanation of some above mentioned mysteries.
Let us begin with the time non-dilation mystery. It is expected that the most radiation from quasars comes from the ambient matter heated by shock waves when St. Augustine's black hole moves relativistically through this matter. Therefore, the associated red-shifts will be truly cosmological as originated from the matter which participates in the Hubble flow. On the other hand, any internal variability related to the black hole itself (for example, due to instabilities in the accretion disc) will not be time-dilated if the St. Augustine's object is nearly motionless with respect to us in our private space.
The mystery related to the origin of supermassive black holes also becomes less acute. Indeed, St. Augustine's black holes might already have been sufficiently massive upon entering our universe. Moreover, initially St. Augustine's black hole moves ultrarelativistically through ambient matter and extremely high value of its gamma-factor drives it into the fast growing mode [130].
Such high gamma-factors are considered as unrealistic for primordial black holes in the early universe [130], but for St. Augustine's black holes this is just what is expected most naturally from anthropic principle. Indeed, if we assume that the relativistic velocity of the St. Augustine's black hole with respect to ambient Hubble flow matter is essential for triggering the quasar activity, an immediate anthropic principle type consequence would be that we are not a typical observer but rather a special one, namely the one which is nearly at rest with respect to potential St. Augustine's objects population, because otherwise the existence of nearby quasars would make development of intelligent life impossible. At that, of course, we have implicitly assumed that St. Augustine's objects population indeed define a reasonable standard of rest. To make sure that this is a plausible supposition, imagine an analog FRW space-time created in a hydrodynamical laboratory experiment, as advocated in [129]. Then the vessel in which an expanding relativistic fluid mimics the FRW universe provides a natural standard of rest and various residual gas molecules, as candidate St. Augustine's objects for this simulated FRW universe, approximately can be considered as nearly at rest with respect to the vessel from the point of view of participants of the relativistic flow.
As relative velocities of St. Augustine's black holes with respect to ambient Hubble flow matter increase with redshift, it is expected that high redshift quasars will be more luminous than low redshift quasars. Therefore, the observed evolution of quasar luminosity over redshifts is no longer mysterious but finds its explanation if quasars are related to St. Augustine's objects.
The shock wave, which surely accompanies the motion of the St. Augustine's black holes through ambient matter, is expected to trigger intense star formation in this ambient matter and hence increase the expected supernovae rates. Thus, it becomes possible to explain the appearance of unusually young chemically mature galaxies around high-redshift quasars.
Probably, shock waves caused by St. Augustine's objects can give a clue also to the most mysterious of paradoxes related to quasars: alleged correlations of objects that have vastly different redshifts. Schematically this is explained on Figure 4. At a space-time point A a St. Augustine's object enters into the Milne universe and at the present time we, situated at the point E, observe the light emitted by the corresponding quasar at point C. The accompanying shock wave propagates along the world line AB and at a point B triggers the formation of the seed of future galaxy which sends to us the light from the point D when it is already formed. As a result, we observe a distant quasar and much closer galaxy aligned, and this alignment is really not by a chance, but due to above described fancy causal relation between them.
It can be expected that the shock wave will give birth not to one but to several galaxies along its world line. Therefore, we expect chains of galaxies with different redshifts all aligned with a quasar. Remarkably, according to Arp [68], there are actually many examples of such chains of galaxies throughout extragalactic space. In fact this effect is so pronounced that it led Victor Ambartsumian in 1958 to the exotic idea that new galaxies are formed through the ejection from older active galaxies [70]. He came to this idea by simply looking at pictures of galaxies. Eight years later Arp came to the same conclusion by inspecting much more detailed and significantly better quality images of peculiar galaxies [131]. Ambartsumian presented his conclusions at the prestigious Solvay conference. However, "this select group of the best known scientists in the world had either been completely baffled or laughed privately at these crazy ideas" [131]. Interestingly enough, the paradigm of St. Augustine's objects reconciles Ambartsumian's and Arp's exotic idea that galaxies beget galaxies with the standard cosmological interpretation of redshifts in a manner which does not undermine the standard theory of galaxy formation for the bulk of galaxy population, if the number of the St. Augustine's objects is small enough.
As the last bonus, the paradigm of St. Augustine's objects also explains why the observation of an apparent superluminal motion in quasars is not a relatively rare phenomenon. In public space, St. Augustine's objects move relativistically nearly towards us so that the associated blue-shift nearly compensates the cosmological redshift for any internal variation of intensity. As a by product, relative motion of any quasar-related two radio emitting blobs will be with a high probability radial and hence can lead to the apparent superluminal motion.

CONCLUDING REMARKS
We have heavily used the Milne model in this article to justify our suppositions. It is time to ask what are the relations of this model with reality. We introduced the Milne metric as a metric of an empty FRW universe. The present density of matter is rather low in terms of the critical density. Therefore, probably, it is not surprising that the Milne metric describes some aspects of the present day universe quite satisfactorily. Namely, curiously the Milne model explains supernovae observations rather well without requiring any dark energy [132,133]. The standard concordance ΛCDM model with positive cosmological constant, of course, provides a better fit and, as expected, an empty universe is not its viable alternative. Nevertheless, "the Milne model has great pedagogical value, elucidating the kinematic aspect of the universe's expansion" [133].
However, in fact, Milne's aspiration was significantly more ambitious than just providing late time approximate description of FRW cosmology in terms of special relativity. Milne believed that some very general rational and at the same time aesthetic principles were at the foundations of the observed order and regularities in the universe. His approach, named kinematic relativity [94], was conceived as a rival theory to the standard relativistic cosmology. Milne's aim was "to deduce as much as possible merely from the cosmological principle and the basic properties of space, time and the propagation of light" [134] revealed by special relativity.
Milne's unorthodox methodology triggered a heated debate involving the most eminent cosmologists of the time [135][136][137]. When the kinematic relativity was first developed "it met great hostility and was criticized very severely, often unjustly, and sometimes frivolously" [134].
In retrospect, the influence of Milne on the development of modern cosmology was very substantial and he inspired others, in particular Robertson and Walker, in shaping fundamental concepts of modern cosmology [135][136][137]. However, as for the kinematic relativity, very little work has been done on it since the death of Milne and the theory has been left "in a curiously unfinished state" [134]. Walker showed [138,139] that Milne's cosmological construction, if described in terms of geometry, was fundamentally different from that of general relativity and, in general, such a description requires Finsler geometry, a generalization of Riemannian geometry (for applications of Finsler geometry see, for example, [140,141]). Although we feel that the Finslerian perspective of the Milne model was not sufficiently explored, at present the Milne model cannot be considered as a viable alternative to the standard relativistic cosmology, of course. Therefore, in the spirit of this paper we can use it as only a source of inspiration and not a solid ground to justify the introduction of the St. Augustine's objects. Can we find a support for them in more mainstream theories?
In non-empty FRW cosmological models we cannot define global inertial frames and thus, we cannot introduce private time like Milne did. However, a plausible generalization of the concept of private time can be envisaged in the general case as well thanks to the fact that Weyl Tensor for all FRW Cosmological metrics vanishes and thus all FRW space-times are conformally flat [142,143]. Suitably defined conformal time can serve as a substitute of private time. For example, in the case of the flat FRW metric we can introduce conformal time-coordinate T and conformal radial coordinate R in the following way [143] where η = t t0 dt a(t) with arbitrary t 0 . Substituting the inverse transformations into Eq. (4) and taking into account that according to Eq. (113) dt = a(t)dη, we find that the metric Eq. (4) (with k = 0) indeed takes the conformally flat form: ds 2 = a(t(T, R)) 2 (T 2 − R 2 ) 2 dT 2 − dR 2 − R 2 (dθ 2 + sin 2 θdφ 2 ) According to Eq. (112), 0 < ψ < η region is mapped onto the region 0 < R < T , while the region 0 < η < ψ-onto the region R < T < 0. As in the case of private space and time in the Milne model, the FRW universe coordinates occupy only a quarter of the (T, R) plane and this circumstance hints at a possibility that the FRW universe might be just a part of some bigger entity. Let us note, however, that the conformal coordinates considered are not exact analogues of the Milne's private time and private space. In particular, transformations Eq. (112) are singular at the light-cone η = ψ.
Another analogue of the Milne foliation is the Rindler's foliation of a non-empty FRW space-time [96] by spacelike slices of finite volume. Rindler showed that, just like the Milne universe, "every open big-bang FRW universe can be regarded as an expanding finite-volume ball of matter, springing from zero volume" [96]. This feature of FRW space-time becomes particularly clear from the multi-dimensional perspective [96]. It is well known that any fourdimensional space-time can be locally embedded (immersed) in a ten-dimensional flat pseudo-Euclidean space [144]. Because of the high degree of symmetry of the FRW space-time, its embedding is possible even in a five-dimensional pseudo-Euclidean space [145]. Namely, the embedding of a spatially flat FRW spacetime ds 2 = dt 2 − a 2 (t)(dx 2 + dy 2 + dz 2 ) into the five-dimensional Minkowski space is defined by relations [145]: T = a(t) 2α (α 2 + r 2 ) + 1 2α t t0 dτ a(τ ) , X = a(t)x, Y = a(t)y Z = a(t)z, Z 5 = a(t) 2α (α 2 − r 2 ) − 1 2α t t0 dτ a(τ ) = αa(t) − T where r 2 = x 2 + y 2 + x 2 and α is a constant with the dimension of length. It follows from Eq. (117) that [145] In the right part of Eq. (118), the cosmic time t can be considered as a function of T + Z 5 because, according to Eq. (117), Therefore, Eq.s (118) and (119) define a hyper-surface in the five-dimensional Minkowski space. For example, in the radiation dominated universe a(t) = t/τ with some τ , and the Eq. (118), defining the hyper-surface, takes the form (in the case t 0 = 0) In fact, the FRW space-time corresponds to only a portion of this hyper-surface, because according to Eq. (119) we must have T + Z 5 ≥ 0.
As we see, a multi-dimensional perspective allows to consider a non-empty FRW universe as a part of a bigger entity and the overall picture is similar to the one in the Milne model that led us to the St. Augustine's objects hypothesis.
It is worth mentioning that extra-dimensional models were very popular in the last decade. In such models, an ambient space is not the Minkowski space but has more complex geometry. The universe as a thin shell expanding in a five-dimensional space with so called non-factorizable warped geometry was first considered by Gogberashvili in [146], and such models with warped extra dimensions became very popular after the Randall and Sundrum paper [147]. Other popular extra-dimensional models are Kaluza-Klein type models with compactified extra-dimensions [148]. It is assumed that all Standard Model particles are confined to the four-dimensional hyper-surface, a brane, in a higher dimensional space. Gravity, on the contrary, can freely propagate in the bulk thus explaining the weakness of the gravitational interactions between the brane confined particles by accompanying spreading of gravitational flux into the large volume of the extra dimensions. In our terminology, gravitons are St. Augustine's objects in these models, and who knows what other type of St. Augustine's objects are lurking in the bulk.
Gravitons also play the role of St. Augustine's objects in Penrose's conformal cyclic cosmology [149]. Gravitational waves emitted in close encounters between super-massive black holes in the last contracting phase of the previous aeon can propagate into the next aeon giving rise to fancy concentric circular rings with slightly different temperature on the Cosmic microwave background sky [149].
Finally, let us mention the so called Ekpyrotic cosmology in which the Big Bang is associated not with a singularity but with a collision of two branes in the extra-dimensional bulk [150]. As a result, cyclic model of the universe can be constructed in which there is an endless sequence of cosmic epochs of expansion and contraction [151]. Black holes as possible St. Augustine's objects in such cyclic scenarios, surviving a bounce of the Big Crunch-Big Bang transition, are discussed in [152]. Interestingly, modern attempts to discuss sailing through the big crunch-big bang transition involves the (compactified) Milne universe in an essential way [122,153]. Milne's influence on the modern cosmology is ongoing!