Dispersive Friedmann universes and synchronization

We introduce consideration of dispersive aspects of standard perfect fluid Friedmann cosmology and study the new qualitative behaviours of cosmological solutions that emerge as the fluid parameter changes and zero eigenvalues appear in the linear part of the Friedmann equations. We find that due to their insufficient degeneracy, the Milne, flat, Einstein-static, and de Sitter solutions cannot properly bifurcate. However, the dispersive versions of Milne and flat universes contained in the versal unfolding of the standard Friedmann equations possess novel long-term properties not met in their standard counterparts. We apply these results to the horizon problem and show that unlike their hyperbolic versions, the dispersive Milne and flat solutions completely synchronize in the future, hence offering a solution to the homogeneity, isotropy, and causal connectedness puzzles.


Dispersive dynamics and bifurcations
It is a mainstream assumption in cosmology that the history of the evolving universe is made up of a number of different epochs such as the vacuum, inflationary, radiation, matter, and other periods, coherently joined to produce the standard cosmological scenario from some early phase to the future. This scenario usually assumes a cosmic fluid with an equation of state given by p = (γ −1)ρ, relating the fluid density ρ to its pressure through the fluid parameter γ, specific values of which are associated with special epochs of cosmological evolution, cf. e.g., [1], [2].
An efficient way to understand and unify many of the dynamical aspects of a variety of cosmological models is through the application of dynamical systems techniques, cf.
e.g., [3]. In this way, using only the simplest themes and ideas from a rich underline dynamical theory, one arrives at a reliable picture of many of the possibilities available in cosmological dynamics, which are difficult to unravel by other methods.
In introductory dynamical systems treatments ( This is, however, a very strong restriction on a vector field. Systems that describe physical phenomena, in particular cosmological dynamical systems, always depend on parameters that are never known exactly, and, in addition, have fixed points that are not hyperbolic, i.e., the linearized jacobian has some eigenvalues on the imaginary axis.
In this case, however, contrary to what happens in structurally stable systems, radically new dynamical behaviour may arise.
Bifurcation theory is the broad framework that deals with this problem, where one is interested not in the properties of an individual hyperbolic system but in a family of dynamical systems that depend on one or more parameters, such that the fixed points are non-hyperbolic for some values of the parameters. It comprises the study of all sorts of topological changes and qualitative reorganizations of the dynamics that result from a change of the parameters on which the system depends. In distinction to structurally stable systems, no complete analysis exists of bifurcating or dispersive systems. Instead, progress has been made in two fronts: local bifurcation theory which studies the orbit structure near a non-hyperbolic fixed point of a parametrized vector field, and global (non-local) bifurcation theory that examines dispersive phenomena which cannot be explained by local bifurcations.
In this case, nearby vector fields can and do have very different orbit structure than that of the individual vector field that corresponds to the value of the parameter associated with the non-hyperbolic fixed point. In an important sense, bifurcation theory is realized by a program of study initiated and emphasized by H. Poincaré to describe changes in the behaviour of the solutions of an individual system (with non-hyperbolic behaviour) that result from its suitable embedding in a parametric family of dynamical systems, these changes being unavoidable and unremovable and occur as the parameters vary.
Dispersive methods include bifurcation theory, but also singularity theory and its applications, the latter are sometimes called 'catastrophe theory' [5]. Whereas the theory of bifurcations describes qualitative changes in the phase portraits of dynamical systems, singularity theory is concerned with such changes for more general objects (such as mappings, manifolds, etc) as the parameters change, and because of this it is very closely related to bifurcation theory. These methods are also ideal for the description of the emerging behaviours in complex systems showing some type of collective behaviour, resulting from interactions between their different parts. Such behaviours are again due to the joint effect of two factors, namely, various degeneracies because of the presence, in the linear part of the vector field, of eigenvalues on the imaginary axis, and secondly because of the presence of parameters.

Dispersive problems in cosmology and nonlinear dynamics
There are various important situations where such dispersive behaviour may arise in a cosmological context. In fact, a fluid parameter, a cosmological constant, the mass of a scalar field, etc, can all be considered as separate parameters in the sense of bifurcation theory. As such, their presence in the cosmological equations is enough to easily produce dynamical systems with an overwhelming dispersive behaviour, leading to novel effects which are totally absent in the standard approach in which these parameters are taken to be constants.
Like in general nonlinear dynamics problems, in cosmology one expects such phenomena to arise when a change in a parameter of the system occurs, and will be always accompanied by dynamical behaviour in which the linear part of the vector field has eigenvalues on the imaginary axis corresponding to a 'bifurcating value' of the parameter. In such a case, the structure of the linear part of the system determines the degree of 'degeneracy', or complexity, of the new equilibria, and leads to possibilities of renewed significance: the vector field is then expected to 'unfold', and the problem becomes one to capture its extended dynamics.
This approach has to our knowledge not been seriously applied before in a cosmological context, and in this paper we make a start by looking for dispersive behaviour in the simplest, nontrivial cosmological model, a Robertson-Walker spacetime with a perfect fluid with an equation of state p = (γ − 1)ρ, where γ is the fluid parameter. This leads to a different situation than that which one is accustomed with in the context of standard Friedmann cosmology, and below we explore the new possibilities to their full extent. To compare our findings with the standard 'hyperbolic' approach to cosmological dynamics, we have also chosen to apply our results to the problem of synchronization, a problem that is closely linked the cosmological horizon problem.
Such 'collective' behaviour of dynamical synchronization ('sync' hereafter) is in fact very common among different areas of nonlinear dynamics, and eventually wins over possible inherent difficulties in communication [6]- [15]. From cyberspace security and social networks, to biological and laser oscillators, to chaotically synchronized behaviour, one witnesses an arbitrary number of identical, coupled subsystems gradually approaching a synced state to remain there forever. The almost magical common property of sync in these systems is that it emerges spontaneously through self-organization, in the sense that although the 'receiver' has received only part of the information of the 'transmitter', it somehow manages to reconstruct the remaining piece.
In this paper, we consider the evolution of a single Friedmannian domain defined by a Robertson-Walker spacetime satisfying the standard cosmological evolution equations with a single perfect fluid with the equation of state p = (γ − 1)ρ. This represents a single completely homogeneous and isotropic spatial domain evolving in (proper) time.
However, the known causality properties of spacetime lead to it containing causal subregions, i.e., parts of any spatial slice of it, that were never in causal contact with each other, yet they now appear to an observer as being completely synchronized -the horizon problem.
The horizon problem of having synchronized but causally disjoint subregions in a homogeneous and isotropic Friedmannian domain, has been known for quite some time, and our description of it completely adheres to the standard one as the following quotations show: In [16], p. 815, we read: . . . one concludes that the foundations for the homogeneity and isotropy of the universe were laid long before the universe became approximately Friedmann, for if statistical homogeneity and isotropy of the universe had not already been achieved at the longest wavelengths earlier, these horizon limitations would have prevented any further synchronization of conditions over large scales while the universe was in a nearly Friedmann state, and small amplitude (10%) deviations from isotropy should be observed now. . . , in [17], p. 506, . . . By comparing radiation background intensities across the sky it is also found that the temperature and expansion rate are precisely synchronized across the visible universe.
Even though the separate parts of the visible universe are not visible to each other they are evolving in precise unison. . .
while in inflation one postulates that such regions were initially 'unsynchronized' but were brought to become causally connected very early, possibly around Planck time, due to a rapid thermalization process, [18], p. 54: . . . expansion began practically simultaneously in different regions of the observable part of the universe with a size l ≤ 10 −33 cm, since they all came into being as a result of inflation of a region of the universe no bigger than 10 −43 sec which started simultaneously to within ∆t ∼ t P ∼ 10 −43 sec. The exponential expansion of the universe makes it causally connected at scales many orders of magnitude greater than the horizon size in a hot universe, R p ∼ ct. . . .
In addition, a lack of proper causal synchronization of spatial subdomains was generally believed to lead to inhomogeneities and anisotropies, as the following passage from [19], pp. 525-6 shows: . . . If the homogeneity of the universe is achieved by the physical transport of energy and momentum . . . we should expect the universe . . . [at last scattering] . . . to be inhomogeneous over distances larger than twice the 'particle horizon' . . . the microwave background ought to exhibit large anisotropies . . . However, there is no sign of any appreciable anisotropy . . . on the contrary, the microwave radiation appears to be highly isotropic [therefore] . . . it is difficult to understand how such a high degree of isotropy could be produced by any physical process occurring at any time since the initial singularity. . . .
In this paper, the possibility of explaining this problem by a non-instantaneous but gradual, or dynamical, synchronization mechanism is proposed as a kind of temporal analogue of phase transitions 1 . From the above quotations, we are naturally led to ask the following two questions presently: Can a single Friedmann domain that has initially become completely synchronized and homogenized by some process (e.g., inflation) maintain its synchronization property during its future evolution? Secondly: can causally disconnected Friedmann subdomains of the single Friedmann domain synchronize during their evolution even if they were not so initially (so that the single Friedmann domain was also not synced initially)?
Perhaps somewhat surprisingly, in this paper we show that the answer to both questions is 'no', if one restricts to the simplest Friedmann-fluid models with hyperbolic equilibria. However, the examination of the dispersive analogues of the simplest Robertson-Walker cosmologies reveals that the models synchronize completely. This we believe provides a new approach to the horizon problem.

Outline of this paper
Our work develops bifurcation phenomena in cosmology, and also points to a possible connection between bifurcations and sync in a cosmological context.
The structure of this paper is as follows. In the next Section, we provide a more detailed guide to the main results of this paper. Section 3 contains basic background material in bifurcation theory. This is used in an essential way throughout the remaining of this work. Sections 4-7 contain results and developments about the dispersive nature of Friedmann cosmologies. In Section 8, we formulate the sync problem in Friedmann cosmology, and then apply our previous results on bifurcating Friedmann universes to this problem. We discuss our findings in the last Section. 1 A totally different kind of synchronization problem based on Mixmaster domains all chaotically oscillating and interacting with weaker couplings and cooperating to sync on approach to the initial singularity, was recently considered in [20,21].
In this Section, we describe the most important results and developments of this work, to guide the reader to the essential points of our arguments. We may tabulate some essential points and results as follows: 1. We compute the center manifolds of all dispersive equilibria of the Friedmann-fluid system.
2. We prove that the unfoldings of the dispersive Milne and flat states are versal.
3. We discuss the implications of the existence of a versal unfolding for the long-term behaviour of the cosmological solutions. 4. We show that the only systems that synchronize in the future are the dispersive Milne and flat unfoldings (horizon problem).
A more detailed description of the results of this paper is given below. These theorems are valid for the non-parametric and parametric center manifolds respectively, and reduce the dimensionality of the original system (3.1) to another system of lower dimensionality holding on the center manifold. These theorems also guide us to compute the center manifold efficiently using the so-called 'tangency condition' (parametric or not).
In subsection 3.5, we introduce the notions of bifurcation, codimension, and unfolding.
We also discuss the important case where one is able to find one unfolding that includes them all, the versal unfolding.
In subsection 3.6, we develop these ideas in the simplest context of one-zero-eigenvalue, codimension-1 bifurcations, which correspond to those found in in later Sections.
Sections 4-7 are the main body of this work. After a preliminary discussion and motivation in subsection 4.1 of how a bifurcation theory approach arises in cosmology for the Friedmann system (4.1)-(4.2), we present the extra dispersive (i.e., non-hyperbolic) equilibria for single-fluid Friedmann cosmology in Section 4.2. We call these fixed points, type-I, IIa, IIb-1, and IIb-2, and show that they correspond to dispersive versions of the Milne, flat, Einstein-static, and de Sitter universes respectively.
Then for each one of these basic solutions, we show that the dynamics is suitably reduced to the corresponding parametric center manifold in each case, namely, Eqns. the Milne state, and (6.14) for the flat state, and we discuss their properties.
It is interesting that no such result can hold for the Einstein static and de Sitter spaces in the present context, and we show that these represent non-generic solutions in this sense (cf. Section 7).
In Section 8, we revisit the two sync questions stated in the Introduction. In subsection 8.1, we show that with respect to sync, the problem of a distribution of causally disjoint domains is reduced to that of the dynamics of a single Friedmann domain.

Center manifolds and bifurcation theory
In this Section, we describe basic results of bifurcation theory necessary for the devel-  [5], [25]. We assume that the reader has some background in flows, stability, invariant manifolds, and asymptotic sets, as the subject is given in standard sources, e.g., [4,26,27]. We exclusively deal below with vector fields, although everything we discuss is directly applicable to maps.

Generalities
We consider the family of dynamical systems, a 'parametrized vector field', where G is a C r function on some open subset of R n × R p , and suppose that The simplest way in which (Y 0 , 0 ) can be non-hyperbolic is when D Y G(Y 0 , 0 ) has a single zero eigenvalue and the remaining eigenvalues have nonzero real parts, this belongs to the so-called 'codimension-1 bifurcation theory', and is the only case we consider here.
(Our treatment does not develop codimension-1 bifurcation theory with a pure imaginary pair of eigenvalues, nor higher-codimension bifurcations, but in fact also covers the case where D Y G(Y 0 , 0 ) has more than one but not all of its eigenvalues equal to zero and the remaining ones having nonzero real parts.) Another name for vector fields with some zero or purely imaginary eigenvalues is dispersive, borrowing from similar usage in partial differential equations (cf. [28]). Hence, we introduce 'dispersive' as another word for 'structurally unstable' vector fields. In this case, a 0-eigenvector defines a dispersive direction, while for a nonzero eigenvalue λ, we call the λ-eigendirection with λ negative (resp. positive) a diffusive, or stable (resp. anti-diffusive, or unstable) direction. The set of all bifurcating dynamical systems is a subset of the dispersive ones; in fact, it is a proper subset: for example, the equatioṅ

Normal form of the linearization
Given the system (3.1), there is some preliminary work that needs to be done in order to bring it to a form suitable for subsequent treatment.
We may assume without loss of generality that the fixed point in (3.1) is located at  having eigenvalues with zero real parts, B is a (s × s)-matrix, having eigenvalues with negative real parts, and C is a (u × u)-matrix, having eigenvalues with positive real parts (a basic example is obtained when A = 0, B = −1, C = 0, this is in fact a case that often appears in the present work). Then the system (3.2) becomes,  In addition, we shall regard terms of the form µ i x j , µ i y j , µ i z j (with the indices i, j in their appropriate ranges each time) as nonlinear terms, so that they will not be present in the linearized Jacobian.
The original system (3.1) is now in suitable form in order to examine how bifurcation considerations enter in the study of parametrized systems.

The center manifold
Let us first consider the case where there are no unstable directions, that is C = 0, h = 0 in Eq. (3.5), and also no parameters. In this case, the system (3.2) instead of having the form (3.3)-(3.5), becomes,ẋ where f, g are C r functions, and A, B as before. We assume that this system has a fixed point, which is dispersive (that is non-hyperbolic), We give the following definition.
Definition 1 (Center manifold) A center manifold is an invariant manifold for the system (3.7)-(3.8) that can be locally represented as the graph,
We note that here x denotes the dispersive dimension, while y is the diffusive (stable) dimension of the system. (We also use the letter h to denote the center manifold function, without a danger of confusing it with the same letter we used before for the nonlinear function associated with the unstable dimensions.) The fundamental properties of (non-parametric) center manifolds are contained in the following theorem.
Theorem 1 For the system (3.7)-(3.8) there exists a C r center manifold having the following properties: 1. Reduction: The dynamics of the system on the center manifold is given by the following c-dimensional 'reduced' system (c is the number of zero eigenvalues), 2. Stability: The stability of the center manifold solution x(t) is the same as the stability of the full solution x = x(t), y = h(x(t)) of the original system, and, for sufficiently small x, the two solutions (that captured by the center manifold and the original one) are identical up to exponentially small terms.

Approximation:
The center manifold can be computed to any finite order of x for small solutions by using the tangency condition,ẏ − Dh(x)ẋ = 0.
We refer to [29] for a proof.

Center manifolds with parameters
The theory developed in the previous subsection has a simple but important extension when we include the parameter µ ∈ R p . In this case, instead of the system (3.7)-(3.8) we have the parametric system,ẋ We assume that this system has a dispersive fixed point at the origin, and that similar assumptions hold as before for A, B, f, g. Then Theorem 1 still holds on the parametric center manifold, for the reduced vector field,ẋ that is, we have, Theorem 2 For the system (3.13)-(3.15) there exists a C r parametric center manifold given by (3.17), such that the reduced dynamics on W c loc is given by (3.18). The stability properties of the solutions of (3.13)-(3.15) are described by those of the reduced vector field (3.18)- (3.19) up to exponentially small terms. The parametric center manifold can be computed to any finite order of x for small solutions by using the tangency condition, We refer to [A83], [GH83], and [W3] for further comments and for the original references.
We note that for the two center manifold theorems above, the inclusion of the remaining unstable (anti-diffusive) dimensions for both the parametric and the non-parametric cases, is direct, but of course the fixed point would then be unstable. However, all this behaviour will again be captured (or, 'realized') on the center manifold (cf.

Degeneracy, codimensions, and unfoldings
In this subsection we provide a general, intuitive discussion of basic notions of bifurcation theory.
We shall say that (the dynamics of) two (possibly parametric) vector fields f, g (or simply their flows) are C k -equivalent (k = 0, 1, 2, . . . ) if there is a C k -map that takes the orbits of the flow of f to the orbits of the flow of g preserving their orientation. If in addition the map preserves parametrization by time, then we say that f, g are C kconjugate. When the flows of f, g are C k -equivalent, the fixed points of f are mapped to the fixed points of g and the same is true for their periodic orbits, but their periods may not be equal. However, C k -conjugate flows have the same mapped eigenvalues and equal periods of their mapped periodic orbits. For instance, the famous Hartman-Grobman theorem says that every hyperbolic flow is C 0 -conjugate to its linearization (cf. e.g., [W3], p. 350).
We now give the following basic definition. We consider a system of the form (3.18)- that is a general system after center manifold reduction. We shall also assume that (3.20) is one-dimensional. As we remarked above, the condition that the fixed point is non-hyperbolic is a necessary but not a sufficient one for it to be a bifurcation point. Intuitively speaking, the nature of a bifurcation of the system is such that when the parameter µ is varied through the bifurcation point, the system (i.e., the vector field) itself is 'deformed' and new dynamical phenomena occur, for example, fixed points can be created or annihilated.
The most distinct characteristic of a bifurcation is the existence of one or more bifurcation sets, that is smooth submanifolds of fixed points in the phase-parameter Such questions are studied in detail in the sister fields of bifurcation theory, singularity theory, and catastrophe theory, and their consideration becomes very essential for bifurcation problems [5,25,30,31,32].
For a given parametric family of systems such as (3.20), as the parameter of the system varies through its bifurcation value and the vector field changes, we say that the system 'bifurcates', or 'deforms', or 'unfolds', and the problem becomes one to classify all these different 'unfoldings', that is all those families which contain the given bifurcation in a 'persistent' way (cf. [GH83], p. 123). As the parameter changes, the unfoldings contain those behaviours of the system which are unremovable and so should not be missed. This is so because the removal of some degeneracy at one individual value of the parameter in one unfolding, appears anew at a nearby value of another (one cannot remove all degeneracies simultaneously, e.g., [25], p. 16).
To study the unfoldings of a given family, we must perturb it in suitable ways. Not In some cases, the study of all possible unfoldings of a system can be obtained from a single one that yields all other possible bifurcations. This is called the versal deformation, or the universal unfolding of the system, and was first studied by Poincaré in his Thesis. Bifurcation theory provides general methods to compute the codimension of a bifurcations and, in particular, to construct the versal unfolding (we use this term collectively to describe both versal deformations and universal unfoldings (although there are generally subtle differences between the two concepts, these are discussed in [33,34]).

The simplest bifurcations with a single-zero eigenvalue
The remaining problem is to describe the dynamics of (3.1) in the special case of the reduced family (3.18)-(3.19) on the center manifold. We shall assume that A = 0 (but (B = 0), this way we have a 1-dimensional system with a single zero eigenvalue on the center manifold and the remaining eigenvalues have nonzero real parts 2 .
The problem is therefore reduced to describing the dynamics of the 1-dimensional system (similar to the general system (3.20)-(3.21)), near the origin, when, so that it has a non-hyperbolic fixed point at the origin (in both phase and parameter space) with a single-zero eigenvalue, and no terms linear in x (we note that a term of the form µx is nonlinear of order 2). We shall also assume that p = 1 (when more than one parameters are present we regard all of them except one as fixed). 2 Here we do not discuss the Hopf bifurcation, the more complicated case where there is a pair of imaginary eigenvalues instead of a single-zero eigenvalue. However, we note the interesting fact that dynamical considerations restrict the codimension of the Hopf bifurcation (and also others) to a smaller value than that found by using singularity theory considerations alone (1 instead of 2). This explains the word 'simplest' in the title of this subsection.
The simplest bifurcations of the system (3.23) are given below. In each case, we give the name of the bifurcation and the 'normal form' of the vector field (the function f (x, µ)): 1. The normal form of the saddle-node bifurcation: 2. The normal form of the transcritical bifurcation: 3. The normal form of the pitchfork bifurcation: There are bifurcation conditions on f for the behaviour of the system (3.23) to be one of these prototypical three types on the center manifold, but we shall not be repeated The term 'normal form of the saddle-node bifurcation' introduced above means that the dynamics of any system of the form (3.23) with f (x, µ) = a 0 µ + a 1 x 2 + a 2 µx + a 3 µ 2 and having a saddle-node bifurcation is qualitatively the same as that of the systeṁ All O(3)-terms or higher could be neglected without qualitatively affecting the dynamics.
In fact, one may show that the saddle-node and the transcritical are codimension-1, but the pitchfork is a codimension-2 bifurcation. In terms of versal (unfolding) families, one may better understand the significance, or typicality, of these bifurcations, in particular that of the saddle-node's: The systemẋ = ax 2 , a : const, has the saddle-node bifurcationẋ = µ+ax 2 as its versal unfolding, whereas the fixed point ofẋ = ax 3 , a : const is of codimension-2, with versal unfolding the pitchfork bifurcationẋ = µ 1 + µ 2 x + ax 3

Dispersive Friedmann domains
In this Section, after some remarks mainly for motivating the idea of cosmological bifurcations, we discuss the dispersive fixed points of the Friedmann fluid system.

A bifurcation approach to Friedmannian evolution
As discussed in the Introduction, we consider a single Friedmannian causal domain A, that is a spatial region with all 'points' (or subregions inside it) 'homogenized'. (For the moment we shall restrict our attention only to this case, leaving to subsection 8.1 the more general situation where we also allow for subdomains of A to be causally disconnected with respect to others.) The evolution of A is governed by the usual dynamical equations (cf. [3], Sec. 2.3), for the Hubble parameter H =ȧ/a, defined using the scale factor a as a function of the proper time t (a dot is derivative with respect to t), and the density parameter , with the fluid density ρ related to the pressure by p = (γ − 1)ρ, γ being the fluid parameter. Then using a dimensionless time variable τ defined by dt/dτ = 1/H, the evolution equations for A are given by, Here q = µΩ/2 denotes the deceleration parameter, and we have set µ = 3γ − 2.
Below we shall be interested in the question of whether or not the system (4.1)-(4.2) is dispersive, and, in particular, whether or not it admits a bifurcation set of the form (3.22). However, keeping µ constant, as is the case in the usual approach to cosmological dynamics described by (4.1)-(4.2), cannot provide a suitable framework to address this question, because in this way one would restrict the evolution of Friedmann domain A to a particular era in its evolution.
We shall assume instead that the (4.1)-(4.2) defines a family of differential equations parametrized by µ which is now promoted to be a true continuous parameter of the problem 3 . This is a crucial point that plays an important role below, and we shall be interested in the properties of the solutions of the system (4. In order to understand the various possibilities, it is useful to distinguish between two broad categories of dynamical systems: generic and degenerate. For the former, it was shown by Poincaré that the only behaviour of the phase curves that occur in a neighborhood of an equilibrium is either a focus, or a node, or a saddle. This behaviour associated with generic systems is also called hyperbolic, to distinguish it from more complicated ones, like say a (non-hyperbolic) center, or periodic orbit, as in the case of the simple harmonic oscillator, (y, −ω 2 0 x) for each frequency, at the origin. This degenerate behaviour as is well-known is, however, unstable with respect to small generic perturbations of the vector field, such as (y, −ω 2 0 x − y), for constant: the center is altered into a (hyperbolic) sink or a source and is consequently destroyed into stable (for > 0) or unstable (for < 0) foci. That is, more complicated cases turn into generic ones under small generic perturbations of the system.
However, this situation changes dramatically if we are interested not in an individual system but in a family of systems that depend on one or more varying parameters. Let us give one last example to illustrate the distinction between constants vs. continuous varying parameters. This is a Lorenz dynamical system, a simplified model for fluid convection in a 2-dimensional layer heated from below. This system is described by the vector field (σ(y − x), ρx + x − y − xz, −βz + xy), with (x, y, z) ∈ R 3 , σ, β positive constants, and ρ a parameter (usually it is typical to put ρ = r−1 in the standard version of the Lorenz equations). Below we shall consider only the possibility where two of the parameters are fixed and only one remains a true varying parameter (the ρ). If for the moment we consider ρ also as a positive constant, then it is easy to see that the system has a hyperbolic fixed point at the origin, a saddle (when ρ > 0), and further hyperbolic equilibria which are kinds of spiral saddles (one positive and two complex eigenvalues with nonzero real parts). For ρ < 0, the origin is a hyperbolic sink, and becomes the unique attractor of the system. This is the overall hyperbolic behaviour of the Lorenz system. However, as is well known, the Lorenz system contains much more structure than the simplified hyperbolic analysis suggests, even in the simple case we consider, with only one of the parameters as varying (ρ = 0). Indeed, for ρ = 0, a pitchfork bifur-

Dispersive Friedmann equilibria
In sharp contrast to the usual approach to cosmological dynamics but in line with what has been discussed above, we shall take the system (4.1)-(4.2) to play the same role as that which the primary system (3.1) played for subsequent developments in Section 3, and look for the possible existence of dispersive fixed points.
Then the Jacobian J Y G for (4.1)-(4.2) is, The behaviour of the corresponding solutions near the equilibria EQ-I, EQ-IIa, EQ-IIb is analyzed qualitatively below, as orbits in the corresponding phase spaces.
We shall refer to them as Type-I, IIa, and IIb solutions. We emphasize that these solutions correspond to non-hyperbolic (i.e., dispersive) behaviour, and consequently their properties cannot be deduced from the linearized parts of the systems, the same way as in hyperbolic problems.
The main tool for a complete analysis of these issues is bifurcation theory, and in this work we shall use this tool to obtain a detailed picture of the dynamics. Once this is done, we shall see that the system (4.1)-(4.2), albeit being non-hyperbolic, cannot properly bifurcate.
In our ensuing investigations to explain this remarkable conclusion for the Friedmann system (4.1)-(4.2), we shall discover two novel situations: Firstly, that the failure to bifurcate lies in a deeper effect, that of violation of the fundamental bifurcating conditions.
Secondly, we shall be led to the existence and construction of a unique perturbation of the Friedmann equations (4.1)-(4.2) that has the sought for bifurcation property. This in turn will have interesting implications for the structure of the universe. As we saw earlier, the Jacobian of the system (5.2)-(5.4) at (Ω, H; µ) = (0, 0; 0) is, with eigenvalues 0, −1. The center eigenspace is then the set H = 0, that is the Ωdimension, while the diffusive eigenspace (i.e., that corresponding to the eigenvalue −1) is the set Ω = 0, that is the H-axis. Both sets Ω = 0, H = 0 are invariant.

Center manifold reduction
It follows from the parametric center manifold theorem (Theorem 2 of Section 3.4) that for the system (5.2)-(5.4) a center manifold exists for µ sufficiently close to µ = 0, given by, The result is that the center manifold W c loc from Eq. (5.9) is the Ω-axis, and the reduced dynamics is governed by, We note the important fact that in Eq. (5.12), quadratic terms include the Ω 2 term. The significance of this will be seen more clearly later. Since the other eigenvalue is negative, it follows from the center manifold theorem that the unstable manifold W u loc is empty, and hence, all solutions of the full system are stably attracted by the center manifold (i.e., rapidly decay to it).

Do dispersive Friedmann type-I models bifurcate?
In Fig. 1, we see three examples of phase portraits for the family of systems in (5.2)-(5.4), for µ = 0 and two values of µ near zero. Since the diagrams for µ = 0 have no zero eigenvalue, they are qualitative different than that at µ = 0, and one may conclude that the system apparently satisfies the general conditions of a bifurcation in the Definition 2 of Section 3.5.
In practice, however, we need the existence of a curve of fixed points for any system to undergo a bifurcation at the bifurcating value µ = 0 5 . We can use the center manifold reduction theorem of Section 3.4 and focus on the reduced vector field (5.12) to examine more carefully the question of whether the field bifurcates and, in particular, how this relates to the codimension-1 bifurcations (3.25)-(3.27).
We can see the behaviour of f (Ω, µ) more clearly, if we draw the phase-parameter diagram for the vector field (5.12), as in Fig. 2.
The global evolution of the system is characterized by the different signs of the µparameter: when µ < 0, the equilibrium (0, 0) is unstable and (1, 0) is stable, when µ > 0, (0, 0) is stable and (1, 0) is unstable, while at µ = 0 we have the fixed point EQ-I. 5 As we shall see below, this will illustrate the somewhat subtle property that the non-hyperbolicity of a fixed point is a necessary but not a sufficient condition to warrant a proper bifurcation.  However, since f (Ω, µ) violates the transversality condition, namely, it satisfies, ∂f ∂µ (0, 0) = 0, (5.15) it follows that f (Ω, µ) does not satisfy the conditions for a saddle-node bifurcation. Also because it violates the nondegeneracy condition, that is it satisfies, it follows that it cannot bifurcate transcritically. Finally, because it satisfies the condition, it follows that f (Ω, µ) cannot satisfy the conditions of the pitchfork bifurcation either.
The conclusion from these results is that because of the lack of bifurcations curves, and the consequent violation of the transversality and the nondegeneracy conditions, the Friedmann system (5.2)-(5.4) cannot properly bifurcate.

The versal unfolding
We now show that there is a uniquely defined perturbation of the system (5.12)-(5.13) (and so of the original Friedmann equations (5.2)-(5.4)) that satisfies the conditions of the typical saddle-node bifurcation. We reiterate the fact that in this paper we do not consider perturbations of the Friedmann metric, but instead study the issue of structural instability of the Friedmann equations.
Instead of the one-dimensional field f (Ω, µ) = µΩ(Ω − 1) considered above, we now introduce the following perturbation of it, where σ is a new parameter in a Taylor expansion of the field. We call σ the unfolding parameter of (5.2)-(5.4). Although it looks as if (5.18) has two parameters, we can eliminate µ by defining a new time, 19) to combine the two parameters µ = 0 and σ into one, which we call the unfolding-fluid parameter, ν = σ µ , (5.20) to get dΩ dT = Ω(Ω − 1) + ν. Therefore if we could also show that (5.21) is also a generic family, it would follow that (5.18) is a versal unfolding. The proof proceeds in two steps. The first step is to show that (5.21) can be written as the versal deformation of the generic quadratic family.
The second step, which is completed at the end of Section 6, is to suitably shift (5.12) to show that it matches with the generic quadratic family.
We introduce the new variable, There is a very great difference in the implications of the two laws, (5.12) and (5.25), for the structure and evolution of the universe. These become apparent in form of the two bifurcation diagrams in Figs. 2, 6, and this has direct implications for the long term behaviour of the cosmological solutions as we now discuss.
In the case of the standard system (5.12)-(5.13), Fig. 2 shows that for all nonzero values of µ there are always the two equilibria at Ω = 0, 1. In particular, the equilibrium corresponding to Ω = 1 implies the presence of an all-encompassing big bang singularity in the past, that exists for all models.
On the other hand, the bifurcation diagram for the system (5.25) is the standard one for the normal form of the saddle-node bifurcation, cf. Fig. 6 6 . Whenμ > 0, or ν < 1/4, there are two new fixed points for the versal unfolding (5.25) at, that is when, (so that in this definition, Ω ± ≷ 1/2), but no fixed points forμ < 0 (or, ν > 1/4).
This means that when the system (5.12)-(5.13) unfolds to become ( Consequently they experience no big bang or big crunch singularities as the asymptotic limits of the endpoints of their neighborhoods are approached. In Fig. 6, each vertical line is the phase space of one particular 'universe', for instance the standard Friedmann evolution given by (5.12)-(5.13) in general relativity is on thē µ = 1/4-phase line. When the parameterμ passes through its bifurcation value zero, the parabola is created or disappears (if zero is approached by going the other way). In this sense, all vertical phase lines suddenly appear together with their 'prescribed' distances between their fixed points (one stable, one unstable) depending smoothly onμ, or the other way, they disappear because the equilibria combine with one another.

Dispersive type-IIa states
We now turn to the study of the flat equilibrium solution EQ-IIa of the system (4.1)- At the origin of the (v, H, µ) phase-parameter space, the jacobian of the linearized system is again (5.8). This leads to the (3.2)-form of the equations (6.2)-(6.3), Proceeding as before, we find that the center manifold is the v-line, that is the horizontal line Ω = 1, with the evolution equation on it given by the reduced vector field on W c loc , v = µv + µv 2 (6.8) From the dynamics on the center manifold it then follows that when µ < 0 we have v < 0, whereas when µ > 0 we have v > 0, and we arrive at the phase diagram given in Fig. 7. This shows that the flat state F is a past attractor of all models in the sense that, The remaining bifurcation analysis for the vector field (6.8) on the center manifold is analogous to that in Sections 5.2 and will not be repeated here. We find that Eq. (6.8) does not bifurcate.
However, passing on to the unfolding, with the unfolding-fluid parameter ν given again by (5.20), and defining we find the saddle-node bifurcation, Specifically, the unfolding (6.11) is versal. Same comments apply here as those at the end of Section 5.3, in particular, for the motion of the system near the stable orbit of the parabola, that is near W = − √ −μ. We also note that the previous discussion about the long term behaviour and singularities in the type-I unfolding of the previous Section apply here as well.
We end this Section by completing the proof of the second step that (5.18) and also (6.11) are generic families, by reducing (5.12) and (6.2) to purely quadratic vector fields.
Both (5.12) and (6.2) are of the form, We note that here, in distinction to the previous analyses in Sections 5, 6 , the dispersive dimension is H. Proceeding as before, we find that the reduced dynamics on the center manifold is given by, When µ 1 > 0 (or, µ < −2), and H > 0, we have H < 0, and we have stability of the fixed point, whereas when H < 0 we have H > 0, and the solution is unstable. These results describe the stability properties of the dispersive Einstein static universe.
It is also easily checked that due to the violations of the transversality and nondegeneracy conditions, the system (7.5) cannot bifurcate.
We shall comment on the question of unfolding (7.5) at the end of this Section.

The equilibrium (1, 1, −2)
We now come to the stability of the dispersive de Sitter space (1, 1, −2). The linear change of variables, brings the fixed point (1, 1, −2) of the system (4.1)-(4.2) to the origin, where the system assumes the form, For this system, we are interested in the behaviour of solutions near the origin (where also µ 1 = 0). We note that the dispersive direction (i.e., the 0-eigenspace) is the w-axis, while the (−2)-eigenspace is the line v = −2w in the (w, v) phase space.
In this new formulation, we can bring the system (7.10) to having its linear part in canonical form as follows. The eigenvector matrix for the linearized Jacobian of (7.8)-(7.9) is, 14) and so introducing the new variables, we find that the system (7.10) takes the form, where, Eq. (7.16) is a block diagonal system for its linear part, with the zero eigenvalues coming first and the non zero eigenvalue after the zero ones. Therefore we can now proceed with the calculation of the parametric center manifold and set, After some algebra, we find that the dynamics on the center manifold is described by, 20) or, using the definition (7.15), and the fact that on the center manifold, v = 0, we find 21) or using (7.7), This equation is to be compared with Eq. (7.5), the corresponding one for the other type II-b state, they are the same, and so similar stability conclusions apply presently.
We end this Section with a comment about the two families (7.5), (7.22). The vector field f (µ 1 , H) = µ 1 H is an unfolding of the zero vector field f (0, H) = 0, for any H.
However, it is not a versal unfolding of f (0, H), for instance the unfolding g(µ 1 , H) = µ 1 does not have a fixed point for µ 1 = 0, and therefore is not conjugate to f (µ 1 , H) which has an equilibrium at zero. This argument shows that unlike the type I and IIa cases, the type II-b solutions are rather special because they cannot be versally deformed to give generic unfoldings.

Applications to synchronization
We now return to the problem stated in the Introduction, namely,

Common evolution of a pair of uncorrelated subdomains
We imagine that the Friedmannian domain A is the union of an arbitrary number of Let us first suppose that we have two causally disjoint spatial domains A 1 , A 2 that evolve according to different dynamical laws given by: and we are back to the single domain problem.
To examine the two domains for sync (see also below for the general case of an arbitrary number of domains), we assume that the spatial domains

General synchronization
The problem of Question 2 for a distribution of an arbitrary number of causally dis- To examine the distribution of domains A i , i ∈ I for sync, we proceed as in the next four steps: Step 1: Split the transmitter in two subsystems a, b, with the indexing set written as Here we are assuming a common time for all subsystems.
Step 2: The receiver at G receives only the 'v-part' of the signals from the transmitter, and at G we set, Step 3: Using the remaining equations, the 'u-part' of the system (8.11), (8.12), namely, u a = f a (u, v), the receiver then checks the difference in the flows of the remaining avariables, namely, the sync function, Step 4: The single domain A is said to sync in the future (resp. past) provided the orbits of the system (8.11) satisfy, ω → 0, as τ → ∞ (resp. − ∞). (8.15) In this case, A evolves keeping its causally disjoint subdomains in perfect unison.
In the presence of parameters, the above procedure extends analogously to parameterdependent flows.
The almost magical property of dynamical systems that have the ability to sync with each other is that although the receiver has received only part of the information of the transmitters (the 'v-part'), it somehow manages to reconstruct the remaining 'u-part'.

Friedmann synchronization
As we already discussed earlier, the above procedure complies perfectly with our center manifold reduction, through which we arrived at the reduced equation (5.12) holding on the center manifold. To introduce a suitable sync function for our problem, we therefore proceed as follows.
Let us start with the general system (3.20)-(3.21) and set x = Ω, and φ τ (Ω, µ) for the flow of the reduced system (this covers subsequent cases like the system (5.12)-(5.13), or the various unfoldings on the center manifold). We define the sync function ω(Ω, µ) to be, ω(Ω, µ) = |φ τ (Ω + , µ + ) − φ τ (Ω − , µ − )|, (8.16) where Ω ± , µ ± denote any pair of states of Ω and any two values of µ respectively. Usually, in bifurcation problems, µ ± denote values of the bifurcation parameter on either side of its bifurcation value, and similarly for the states Ω ± . The sync function checks for traces of any overall difference in the density function solutions as functions of µ during the evolution of the domain, usually as µ passes through its bifurcation value µ = 0.
We now have the following definition. is stated in future terms, but sync in the past evolution simply corresponds to instead taking the τ → −∞ limit in the above definitions.

Sync and proper-time sync
This assumption i.e., common proper time for all uncorrelated domains -subdomains of A), is a rather strong one because in this way any observer field U on A (that is any timelike, future-pointing unit vector field, parametrized by proper time t), becomes geodesic and irrotational, and such that U = −grad t, across all subdomains of A. This assumption thus implies a fine-tuning of all observer fields across A, and so it is rather special.
However, this problem can be easily rectified by defining t to be not the proper time but a kind of average time over all observers in A, so that there exists a smooth function p > 0 on A with U = −p grad t. In this way, the elapsed proper time between two spacelike hypersurfaces normal to U will be different in different A i 's (that is different U observers). In this case, U is just called 'synchronizable', rather than proper time synchronizable (cf. [37], p. 54 for this terminology where our 'average time' is called a 'compromise time').
Dynamical sync as defined earlier requires the validity of the further condition (8.17) that appears in the sync definition above, and so is more subtle than the standard proper time sync. The former refers to a property of the flow of the Einstein equations (4.1), (4.2), whereas the latter to the flow of the ∂ t vector field being geodesic and irrotational.
In the remaining of this Section, we shall examine how dynamical sync applies to the dispersive states of Type I, II.

Do dispersive type-I states sync?
Looking back at the diagram in Fig. 2  From (5.12), it follows that for any µ = µ − < 0, the dynamics will be qualitatively the same as that of the system, and for any µ = µ + > 0 the dynamics of (5.12) will be equivalent to that of the vector field, (We note however, that the two vector fields in each pair will not be C k -conjugate, k ≥ 1, to each other except when µ ± = ±1.) Here, we used the same symbol x to denote Ω in these ranges of µ (cf. the diagram in Fig. 2).
This result implies that there are states for which the sync function ω(µ − , µ + ) cannot satisfy the sync definition (8.17), and therefore the Friedmann domain A desynchronizes during its future evolution in this case.

Sync for unfoldings of the dispersive type-I states
Let us now examine the versal unfolding (5.25) for sync. Since the system bifurcates at µ = 0, we need to examine each of the two branches of the bifurcation curve in Fig. 6 separately.
In the unstable branch, the flows on each phase line for each µ have the property of moving all nearby phase points away from the branch, and so the difference of any two such flows may not tend to a well-defined limit during evolution. Therefore the sync function will not tend to a definite limit, and so the system will desync in this case.
We therefore find that ω(μ 1 ,μ 2 ) satisfies the definition (8.17), when the sign(µτ ) > 0, and it can happen either in the past or in the future in τ -time.
Since the system bifurcates atμ = 0, it will either follow the stable or the unstable branch depending on whether the system finds itself near the one or the other of its two fixed points.

Sync and dispersive type-IIa states
Using a similar approach like in the previous subsection, it is evident how to proceed in order to examine the reduced vector field (6.8) as well as for its unfolding (6.13) for sync, and so we shall be brief. The results are similar to those achieved for the type-I states sync problem.
The original vector field (6.8) does not sync and the proof is completely analogous to that given in Section 8.3, but now we are interested in the behaviour of states near to the Ω = 1 line, instead of those on the Ω = 0 line. For µ > 0, the vector field is qualitatively equivalent toẋ = x + x 2 , with flow given by the form, while for µ < 0, we find These results imply that the function φ τ − ψ τ is increasing for states with x > 0. We conclude that the limit condition (8.17) is violated and therefore Friedmann domains cannot sync during evolution in this case.
It is not difficult to see using the method of the previous subsection that for the unfolding (6.13), sync is achieved during evolution near the stable branch of the orbit which is now located at the third quadrant of the (W,μ)-plane.

Sync and dispersive type-IIb states
For reasons completely analogous to those discussed in previous subsections of this Section, we conclude that in the case of the type-IIb solutions the model cannot sync. Also for reasons described in Section 7, they cannot unfold to describe synced states. Therefore we may regard type II-b solutions as non-generic manifestations that cannot unfold the degenerate solution H = 0. The main reason is the absence of quadratic terms of the form H 2 in the center manifold reduced vector fields.

Discussion
In this paper we introduced and studied two new factors that affect the dynamics of cosmological models. These are the consideration of parameters -as opposed to constants -in the basic equations that govern the evolution of the universe, and the appearance of zero eigenvalues in the linear part of the vector field that defines these equations.
These factors challenge basic aspects of standard cosmology, and prove that restricting attention to only hyperbolic fixed points is inadequate for the complete treatment of cosmology.
Consequently, we are lead by necessity to the wider consideration of dispersive methods in cosmological dynamics, where the newer fields of bifurcation theory, singularity theory, and their applications come into play and need to be taken seriously. At a most basic level, this view requires the consideration of not the usual 'individual' system of The question of synchronization required for a resolution of the horizon problem was discussed in Section 8. We showed that there is a synchronization function which asymptotically vanishes in the future for the versal unfoldings of the Milne and the flat equilibria. Therefore in these two cases, Friedmann universes completely synchronize in the future. This constitutes a new approach to the horizon problem independently of any consideration of an initial inflationary phase.
There are various directions towards which this work may be usefully extended.
Throughout this work we have neglected the possible effects of a cosmological constant as it leads to a drastic change in the complexity of the problem because of the presence of a second parameter (besides the fluid parameter) in the dynamical equations (the Λ).
Keeping the fluid parameter fixed, a cosmological constant would dramatically change the structure of the Jordan canonical form of the linear part of the vector field, thus probably extending the present results to the next level of degeneracy. Of course other cosmological problems would similarly lead to more degenerate dynamical systems by increasing the codimension of the bifurcation.
An further interesting question related to the present results is about the nature of the underlined gravity theory that contains a perfect fluid as matter source and leads directly to the unfolding (5.18). An answer to this question is required in order to understand the physical meaning of the unfolding parameter σ, or its dimensionless version ν.