Next Article in Journal
Genuine Dilatons in Gauge Theories
Previous Article in Journal
Excess of Soft Dielectrons and Photons
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spontaneous Color Polarization as A Modus Originis of the Dynamic Aether

by
Alexander Balakin
*,† and
Gleb Kiselev
Department of General Relativity and Gravitation, Institute of Physics, Kazan Federal University, Kremlevskaya str. 16a, Kazan 420008, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Universe 2020, 6(7), 95; https://doi.org/10.3390/universe6070095
Submission received: 1 June 2020 / Revised: 30 June 2020 / Accepted: 6 July 2020 / Published: 10 July 2020
(This article belongs to the Section Cosmology)

Abstract

:
We suggest the phenomenological model of emergence of the dynamic aether as a result of decay of the SU(N) symmetric field configuration containing the multiplet of vector fields. The scenario of the transition to the dynamic aether, which is characterized by one unit timelike vector field that is associated with the aether velocity, is based on the idea of spontaneous color polarization analogous to the spontaneous electric polarization in ferroelectric materials. The mechanism of spontaneous color polarization is described in the framework of anisotropic cosmological model of the Bianchi-I type; it involves consideration of the idea of critical behavior of the eigenvalues of the tensor of color polarization in the course of the Universe accelerated expansion. The interim stage of transition from the color aether to the canonic dynamic aether takes the finite period of time, the duration of which is predetermined by the phenomenologically introduced critical value of the expansion scalar.

1. Introduction

The term dynamic aether was introduced in the framework of the Einstein-aether theory, which was formulated in the beginning of 21 century [1,2,3]. This theory operates with a timelike unit dynamic vector field U i , which is interpreted as the velocity four-vector of some cosmic substratum indicated as the aether. The color dynamic aether, which is considered in the presented work, appears in a SU(N) symmetric generalization of the dynamic aether model. The motives of such an extension include physical and mathematical aspects; we briefly consider them in the context of our study.

1.1. Basic Elements of the Einstein-Aether Theory and the Problem of Identification of the Coupling Constants

The Einstein–aether theory belongs to the vector-tensor branch of the science, known as Modified Gravity (see, e.g., [4,5] for details of its classification); it is based on the action functional
S ( EA ) = d 4 x g 2 κ [ R + 2 Λ + λ ( g i k U i U k 1 ) + K ] ,
where R is the Ricci scalar; Λ is the cosmological constant. The function λ is the Lagrange multiplier and the term λ ( g i k U i U k 1 ) has to provide the normalization condition U k U k = 1 . One of the consequences of this normalization condition is the differential constraint U k j U k = 0 . The scalar K = K m n a b a U m b U n is the kinetic term quadratic in the covariant derivative a U m . According to the idea of the authors of the works [1,2,3], the constitutive tensor K m n a b has the form
K m n a b = C 1 g a b g m n + C 2 δ m a δ n b + C 3 δ n a δ m b + C 4 U a U b g m n ,
i.e., it contains not only the basic geometric objects: the metric tensor g m n and the Kronecker tensor δ k i , but the four-vector U i also. The decomposition (2) is equipped by four coupling constants C 1 , C 2 , C 3 , C 4 introduced phenomenologically. When the term K is considered to be free of the four-vector U i , one has to put C 4 = 0 . When the four-vector U i is the necessary element of the constitutive tensor K m n a b , one can consider, formally speaking, the terms C 5 U m U n g a b , C 6 U a U n δ m b , C 7 U b U n δ m a , C 8 U a U m δ n b , C 9 U b U m δ n a , C 10 U a U b U m U n , as the additional elements of the Lagrangian; however, due to the consequence U k j U k = 0 all of these terms disappear from the master equations, and the coupling constants C 5 , , C 10 become hidden. In other words, the canonic variant of the Einstein–aether theory contains four parameters. The constraints on the quartet of coupling constants C 1 , C 2 , C 3 , and C 4 were formulated based on Parameterized Post-Newtonian (PPN) formalism [6]; the details can be found, e.g., in the status report [3]).
In the context of the coupling constant estimation, we would like to mention an interesting example, which is indicated as the Maxwell-like sub-model. It appears when C 1 = C 3 , C 2 = 0 , C 4 = 0 , so that the scalar K contains the square of the skew-symmetric term F i k = i U k k U i only, i.e., K = 1 2 C 1 F i k F i k . Of course, the covariant derivative i U k has both symmetric and skew-symmetric parts, but the only skew-symmetric element appears in the scalar K due to the specific choice of the coupling parameters. In 2017, the events, encoded as GW170817 and GRB 170817A [7], changed the status of the constraint C 1 = C 3 from a hypothesis to the experimental fact. Based on the observation of the binary neutron star merger, it was established that the ratio of the velocities of the gravitational and electromagnetic waves differs from one by the quantity about 10 15 (to be more precise, 1 3 × 10 15 < v gw c < 1 + 7 × 10 16 ). Because the square of the velocity of the tensorial spin-2 aether mode was calculated in [2] to be S ( 2 ) 2 = 1 1 ( C 1 + C 3 ) , one can estimate the sum of the parameters C 1 + C 3 , as follows: 6 × 10 15 < C 1 + C 3 < 1 . 4 × 10 15 . This result admits three variants of interpretation. First, when C 1 + C 3 0 and thus S ( 2 ) 1 , we can additionally require only that the squares of the velocities of the vectorial spin-1 and of the scalar spin-0 aether modes are non-negative, thus providing the inequalities S ( 1 ) 2 = C 1 C 1 + C 4 0 and S ( 0 ) 2 = [ 2 C 1 + C 4 1 ] / [ 3 + 2 C 2 ] 0 , respectively. Second, if we assume that 1 3 × 10 15 < v gw c < 1 , i.e., the spin-2 aether wave mode is subluminal, we have to take into account the possible gravitational Cherenkov effect [8,9]). As it was shown in [10,11], the corresponding anomalous acceleration of cosmic particles is absent, so the spin-1 and spin-0 aether modes should be superluminal, and we have to consider additional inequalities S ( 1 ) > 1 and S ( 0 ) > 1 , providing the constraint 0 < C 2 < 0 . 095 [10]. Third, if we assume that 1 < v gw c < 1 + 7 × 10 16 , i.e., the spin-2 mode is superluminal, the argument concerning the gravitational Cherenkov effect can not be used, and the requirements that the phase velocities of the spin-1 and spin-0 modes are bigger than 1, are not necessary. In this case, the constraint for C 2 is 2 27 < C 2 < 2 21 (see [10,12]).
We see that the recent success in the gravitational wave experiments has provided the progress in the constraining of the Einstein–aether theory. We can hope that the predictions based on the solutions describing black holes [13,14,15], static spherically symmetric objects, neutron stars and binary systems [16,17,18], universal horizons and thermodynamics of the Einstein–aether black holes (see, e.g., [19,20,21]), will also be tested in the nearest future and will introduce new constraints on the coupling constants.

1.2. Einstein–Aether Theory: Historical Motives and Related Models

The Einstein-aether theory [1,2,3] realizes two remarkable ideas. The first idea is connected with the assumption that there exists a preferred frame of reference advocated by Nordtvedt and Will [6,22,23]. Indeed, if one has a timelike global unit continuous vector field U i , one can reconstruct the preferred coordinate system, for which d x i d s = U i ; in this context, the vector field U i can be interpreted as the velocity four-vector associated with the flow of some medium (e.g., the aether). The second remarkable idea, which was suggested by Dirac in [24,25,26], is that the unit timelike vector field can play the dual role: of the vector potential, which satisfies a special gauge condition, and of the velocity four-vector of a charged particle flow. This Dirac’s idea can be extended as follows: let a timelike vector field potential U i as a source of the gravity field exist, but let its norm be not equal to one. In this case, one can introduce the scalar field as the modulus of this vector field, U = g m n U m U n , and the normalized unit four-vector V i U i U . Thus, we deal with a version of the scalar-vector-tensor theory of gravity, based on the interacting trio { U , V j , g i k } .
There exist a class of models, which can be indicated as the related to the Einstein-aether theory. As an example of such model, we would like to consider the so-called bumblebee model (see, e.g., [27,28,29,30]). Such a model differs from the Einstein–aether theory by one term only: it is the term λ ( g i k B i B k ± b 2 ) 2 introduced into the Lagrangian instead of the term λ ( g i k U i U k 1 ) . Now, λ is not a Lagrange multiplier, it is a phenomenological constant; and the new term gives the fourth order potential of the Higgs type. The vector field B k in this theory is not normalized, and the constant | b | is not a norm, but it is the vacuum expectation value corresponding to the minimum of the potential.

1.3. Generalizations and Extensions of the Einstein-Aether Theory

With any approach, the interpretation of the vector field U i as the aether velocity four-vector inevitably revives the discussions concerning the taboo that is associated with the Lorentz invariance violation (see, e.g., [31,32,33]). Breaking this taboo, many authors consider the modifications of the Einstein–aether theory; we distinguish between generalizations and extensions of this theory. The generalization of the first type is based on the introduction of the nonlinear term F ( K ) instead of the linear term K into the Lagrangian, K F ( K ) (see, e.g., [34,35]). Clearly, the idea of such nonlinear generalization is inspired by the f ( R ) -extension of the theory of gravity, in the framework of which the Ricci scalar R is replaced by the nonlinear function R f ( R ) (see, e.g., [4,5]). The generalization of the second type is predetermined by modifications in the gravitational part of the Lagrangian; there are modifications of the Horava–Lifshitz type (see, e.g., [36,37]), the modifications based on the metric-affine gravity [38], on the f ( R ) , f ( G ) , f ( T ) , etc., models of gravity (see, e.g., [39]). The generalizations of the third type involve consideration of extra-dimensions [40] and the ideas of supersymmetry [41,42,43].
We speak about extensions of the Einstein–aether theory, when one already has the unit vector field U i , and one adds new terms into the action functional, which describe the scalar, pseudoscalar, electromagnetic, gauge, etc., fields. The scalar extension of the Einstein–aether theory was motivated by cosmological applications, (see, e.g., [44,45,46,47,48,49]), in particular, to precise some details of inflation, dynamics of perturbations, and the large-scale structure formation. In order to describe the impact of the scalar field φ on the unit vector field U k , the authors of [44] suggested to modify the constitutive tensor K m n a b (2) by introduction of the functions β 1 ( φ ) , β 2 ( φ ) , β 3 ( φ ) , β 4 ( φ ) instead of constants C 1 , C 2 , C 3 , C 4 , respectively. The Einstein–Maxwell–aether theory as an electromagnetic extension of the Einstein–aether theory, was established in [50] and applied to the problem of birefringence in [51]. The Einstein–Yang–Mills–aether theory as a SU(N) generalization of the U(1) symmetric Einstein–Maxwell–aether theory is considered in the work [52]. Pseudoscalar (axionic) extension of the Einstein–aether theory was presented in the work [53]. The next natural step in the development of the Einstein–aether theory was the combination of the electromagnetic and pseudoscalar modifications indicated as axionic extension of the Einstein–Maxwell–aether theory [54,55,56].

1.4. What Is the Aim of This Work?

In the paper [52] we proposed an idea of the color aether, which is based on the introduction of the SU(N) symmetric multiplet of vector fields U ( a ) i . The subscript ( a ) denotes the group index in the adjoint representation of the SU(N) group and this Latin index runs over the values { ( 1 ) , ( 2 ) , ( N 2 1 ) } . If all of the vectors U ( a ) i from this multiplet become parallel, i.e., when U ( a ) i = q ( a ) U i , and, if U i is the unit timelike vector field, we happen to be faced with the version of the dynamic aether, supplemented by the SU(N) symmetric multiplet of scalar fields { q ( a ) } . Clearly, q ( a ) are the scalars from the point of view of spacetime transformations, but the corresponding object with upper index, q ( a ) = G ( a ) ( b ) q ( b ) , can be considered as vectors in the N 2 1 dimensional group space with the metric G ( a ) ( b ) (color vectors, for short). When we deal with the SU(N) symmetric multiplet of vector fields, we have to use the self-consistent formalism, which includes the Yang–Mills gauge fields with the potentials A k ( a ) , the gauge covariant derivative D ^ m ( a ) , etc. This full-format SU(N) symmetric formalism is elaborated and presented in the work [52] as the extension of the U(1) symmetric formalism that is discussed in [50] in the framework of Einstein–Maxwell–aether theory.
In order to establish the phenomenological model of the transition from the color aether to the dynamic aether, we exploit the known analogy with the spontaneous polarization in the ferroelectric materials [57,58]; the term spontaneous means that the electric polarization is induced by variations of temperature or pressure, not by the external electric field. We keep in mind three principal details that appeared in both models. The first detail is connected with the character of evolution. In the theory of ferroelectricity, one can associate the evolution with the temperature decrease: high temperature relates to the symmetric phase without polarization, low temperature corresponds to the polarized dissymmetric phase. When we consider the evolution of the color system in the expanding Universe, the effective temperature is also dropping, providing the transition from the symmetric phase to the dissymmetric one. The second detail is connected with the Curie temperature T C : below this critical value, T < T C , the spontaneous electric polarization appears, and, thus, T C indicates the temperature of the ferroelectric phase transition. Similarly, we can introduce the critical value of the expansion scalar Θ * , which relates to the cosmological phase transition taking place at the cosmological time t * . The third formal accordance is based on the analogy between the electric domain structure in ferroelectric (chaotic in the symmetric phase and aligned in the dissymmetric one) and the anisotropy of the group space that is associated with the color vector fields (multi-axial structure in the symmetric phase and uni-axial in the dissymmetric one). In other words, we assume that, in the early Universe, the phase transition took place, which transformed the SU(N) symmetric multiplet of vector fields to the bundle of vectors, which happen to be parallel in the group space. Below, we describe the mathematical details of the phenomenological model of the spontaneous color polarization.
This paper is organized, as follows. In Section 2, we recall the basic elements of the SU(N) symmetric field formalism and the corresponding master equations. In Section 3 we discuss the possible mechanism of transformation of the color aether, equipped by the SU(N) symmetric multiplet of the vector fields { U ( a ) i } , to the canonic dynamic aether, which possesses one unit vector field U i . Section 4 contains the discussion and conclusion.

2. The Formalism

2.1. Basic Elements of the Theory

We follow the book [59] and use the Hermitian traceless generators of the SU(N) group, t ( a ) , which satisfy the commutation relations
[ t ( a ) , t ( b ) ] = i f ( a ) ( b ) ( c ) t ( c ) ,
where f ( a ) ( b ) ( c ) are the structure constants of the SU(N) gauge group. The scalar product of the generators t ( a ) and t ( b ) introduces the metric in the group space, G ( a ) ( b ) :
G ( a ) ( b ) = ( t ( a ) , t ( b ) ) 2 Tr t ( a ) t ( b ) .
The quantities
f ( c ) ( a ) ( b ) G ( c ) ( d ) f ( a ) ( b ) ( d ) = 2 i Tr [ t ( a ) , t ( b ) ] t ( c )
are antisymmetric with respect to transposition of any two indices. The Yang-Mills field potential A m and the Yang-Mills field strength F i k are defined as follows:
A m = i G t ( a ) A m ( a ) ,     F m n = i G t ( a ) F m n ( a ) ,
where G is the coupling constant, and the multiplets of real fields A i ( a ) and F i k ( a ) are connected as follows:
F m n ( a ) = m A n ( a ) n A m ( a ) + G f ( b ) ( c ) ( a ) A m ( b ) A n ( c ) .
Here and below m is a covariant derivative. The tensors F ( a ) i k has its dual elements
* F ( a ) i k = 1 2 ϵ i k l s F l s ( a ) ,
with universal Levi-Civita tensor ϵ i k l s . The multiplet of vector fields U m ( a ) appears as the decomposition of the Hermitian quantity
U m = t ( a ) U m ( a ) .
The extended (gauge covariant) derivative D ^ m is defined as follows:
D ^ m Q ( d ) ( a ) m Q ( d ) ( a ) + G f · ( b ) ( c ) ( a ) A m ( b ) Q ( d ) ( c ) G f · ( b ) ( d ) ( c ) A m ( b ) Q ( c ) ( a ) + ,
where Q ( d ) ( a ) is arbitrary tensor in the group space (see, e.g., [60]). For the vector fields this formula gives
D ^ m U n ( a ) m U n ( a ) + G f ( b ) ( c ) ( a ) A m ( b ) U n ( c ) .
The dual tensor * F ( a ) i k , due to the definitions (7) and (8), satisfies the relations
D ^ k * F ( a ) i k = 0 .
Finally, we introduce the multiplet of scalar fields Ω ( a ) , which forms the color vector
Ω ( a ) D ^ m U ( a ) m = m U ( a ) m + G f ( b ) ( c ) ( a ) A m ( b ) U ( c ) m ,
and the scalar Ω based on the relationship
Ω 2 = Ω ( a ) Ω ( a ) .
In the Abelian model of the aether with U(1) symmetry the scalar Ω coincides with the expansion scalar Θ m U m .

2.2. The Ansatz

Our additional ansatz is that the vector fields satisfy the condition
G ( a ) ( b ) U m ( a ) U n ( b ) g m n = 1 .
This is the direct generalization of the normalization condition g m n U m U n = 1 , appeared in the canonic Einstein-aether theory, however, we have to keep in mind two new essential details. First, the set of vector fields U ( a ) i can include s 0 timelike vectors with positive norms, s 1 spacelike vectors with negative norms, and N 2 1 s 0 s 1 null vectors. Clearly, the number s 0 can not be equal to zero, s 0 1 , but this is the only general condition for the numbers s 0 and s 1 . The second detail is the following: in contrast to the condition U k j U k = 0 , we see that generally U k ( a ) j U ( b ) k 0 . This means, in particular, that the terms of the type C 5 U m ( c ) U n ( d ) g a b in the generalized constitutive tensor will not be hidden.

2.3. Action Functional and Master Equations

We use the action functional, which describes interaction of gauge, vector, and gravitational fields
S = d 4 x g { 1 2 κ [ R + 2 Λ + λ ( U ( a ) m U m ( a ) 1 ) + K ( a ) ( b ) i j m n D ^ i U m ( a ) D ^ j U n ( b ) ] + 1 4 F m n ( a ) F ( a ) m n } .
The full-format version of the constitutive tensor K ( a ) ( b ) i j m n is presented in the paper [52]; that version contains 75 coupling constants, since now U k ( a ) j U ( b ) k 0 , and the phenomenological approach requires to list all of the possible terms. In this work, based on the direct analogy with the canonic Einstein–aether theory and following the idea to study some illustrative toy model, we postulate the constitutive tensor to contain only four coupling constants:
K ( a ) ( b ) i j m n = C 4 U ( a ) i U ( b ) j g m n + G ( a ) ( b ) [ C 1 g i j g m n + C 2 g i m g j n + C 3 g i n g j m ] .
Clearly, for the U(1) - symmetric theory the group indices disappear, and this constitutive tensor coincides with the Jacobson’s one (2). On the other hand, this example explicitly shows that the constitutive tensor K ( a ) ( b ) i j m n is not generally proportional to the metric G ( a ) ( b ) , and in this sense it does not possess the isotropy with respect to the group indices.

2.3.1. Master Equations for the Gauge Fields

Master equations appear as a result of variation of the action functional (16) with respect to four quantities: the Lagrange multiplier λ , the vector fields U ( a ) j , the gauge potential four-vector A i ( b ) , and space-time metric g p q . The variation with respect to λ gives the normalization condition (15). Variation with respect to A i ( a ) gives the extended Yang–Mills equations
D ^ k F ( a ) i k = Γ ( a ) i ,
where the color current Γ ( a ) i is given by the formula
Γ ( a ) i = G κ f ( c ) ( a ) ( d ) U k ( c ) K ( d ) ( b ) i m k n D ^ m U n ( b ) .
This quantity appears, when we fulfill the variation of the second term in the right-hand side of (11) with respect to A i ( d ) . From the physical point of view, the color current Γ ( a ) i is induced by the interaction between the color vector field and non-Abelian gauge field; clearly, it disappears in the U(1) case, since the structure constants are equal to zero in that model.

2.3.2. Master Equations for the Vector Fields U ( a ) k

The variation with respect to U ( a ) j yields the standard balance equation
D ^ i J ( a ) i j = λ U ( a ) j + I ( a ) j ,
where the color objects J ( a ) i j and I ( a ) j are introduced in analogy with the canonic aether theory [1]:
J ( a ) i j = K ( a ) ( b ) i m j n D ^ m U n ( b ) = C 1 D ^ i U ( a ) j + C 2 g i j D ^ m U ( a ) m + C 3 D ^ j U ( a ) i + C 4 U ( a ) i U ( b ) m D ^ m U j ( b ) ,
I ( a ) j = C 4 D ^ j U m ( a ) U ( b ) n D ^ n U m ( b ) .
Finally, the Lagrange multiplier can be found as
λ = U j ( a ) [ D ^ i J ( a ) i j I ( a ) j ] ,
i.e., the set (20) contains 4 N 2 5 independent equations for the vector fields U ( a ) k (instead of three ones in the U(1) model).

2.3.3. Master Equations for the Gravitational Field

The gravitational field is described by the set of equations
R p q 1 2 R g p q = Λ g p q + λ U p ( a ) U ( a ) q + T p q + T p q ( YM ) .
The stress-energy tensor of the color vector fields T p q is of the following form:
T p q = 1 2 g p q K ( a ) ( b ) i j m n D ^ i U m ( a ) D ^ j U n ( b ) + C 1 [ D ^ m U p ( a ) D ^ m U q ( a ) D ^ p U m ( a ) D ^ q U m ( a ) ] +
+ C 4 U ( a ) m D ^ m U p ( a ) U ( b ) n D ^ n U q ( b ) + G ( a ) ( b ) D ^ m [ U ( p ( b ) J q ) m ( a ) J m ( p ( a ) U q ) ( b ) J ( p q ) ( a ) U m ( b ) ] .
Here, the parentheses ( p q ) denote the symmetrization with respect to the coordinate indices p and q. The symbol T p q ( YM ) relates to the stress-energy tensor of the Yang–Mills field
T p q ( YM ) = 1 4 g p q F m n ( a ) F ( a ) m n F p m ( a ) F ( a ) q n g m n .
Thus, the total set of master equations is composed by (12), (18), (20) and (24).

2.4. The Tensor of Color Polarization

Based on the multiplet of the color fields U k ( a ) , we introduce an auxiliary quantity
H ( a ) ( b ) g p q U p ( a ) U q ( b ) .
It presents a set of scalar functions from the point of view of spacetime transformations, and is a symmetric real tensor in the group space. The quantity H ( a ) ( b ) can be indicated as the tensor of color polarization, the N 2 1 dimensional analog of the two-dimensional polarization tensor in optics (see, e.g., [61]). The components of this tensor depend on time and spatial coordinates. If we work with the matrix H ( b ) ( a ) G ( c ) ( b ) H ( a ) ( c ) as with the algebraic quantity calculated in the fixed point of the spacetime, we can find N 2 1 eigenvalues, σ { α } , and the corresponding N 2 1 eigenvectors q { α } ( a ) (the subscript { α } indicates the serial number of the eigenvalue):
H ( b ) ( a ) q { α } ( b ) = σ { α } q { α } ( a ) .
Because the trace H ( a ) ( a ) is equal to one due to the relationship
H ( a ) ( a ) = G ( a ) ( b ) g m n U m ( a ) U n ( b ) = 1 ,
the sum of eigenvalues is equal to one,
α = 1 N 2 1 σ { α } = 1 .
One can decompose the tensor H ( a ) ( b ) into the series of products of eigenvectors:
H ( a ) ( b ) = α = 1 N 2 1 σ { α } q { α } ( a ) q { α } ( b ) .
There are two physically motivated limiting cases that are associated with the tensor of color polarization.
  • When all of the eigenvalues are equal to one another, and thus σ { α } = σ = 1 / ( N 2 1 ) , we deal with the color analog of the so-called natural light [61].
  • When all of the eigenvalues, except one, are equal to zero, we deal with the color analog of the linearly polarized light; in this case only one eigenvalue is non-vanishing, say, σ { 1 } , and it is equal to one. We obtain now H ( a ) ( b ) = q { 1 } ( a ) q { 1 } ( b ) , where q { 1 } ( a ) is the corresponding eigenvector; the determinant of the matrix of this color tensor is equal to zero, the matrix is degenerated and, thus, the equality det H ( a ) ( b ) = 0 can be considered as a criterion for recognition of the color polarization.
The crucial point of the presented model is the spontaneous process, which transforms the multiplet of color vector fields U ( a ) i into the bundle of vectors q ( a ) U i , which are parallel in the group space. We can monitor this process while using the behavior of the eigenvalues σ { α } of the matrix of color polarization H ( a ) ( b ) . Our ansatz is that N 2 2 eigenvalues have the stepwise form
σ { α } = η [ Ω { α } * Ω ] σ { α } 0 [ Ω { α } * Ω ] μ .
Here, η [ z ] is the Heaviside function that is equal to zero, when the argument z is negative. The quantities Ω { α } * , where { α } = 2 , 3 N 2 1 , are the critical values of the scalar Ω introduced in (13) and (14). These critical values can coincide or be different. The quantity σ { α } 0 is a constant, and the parameter μ plays the role of critical index; we assume that μ > 1 providing the function σ { α } and its derivative to be equal to zero at the critical moment. This stepwise representation means that there are critical moments of time t { α } * , or equivalently, critical values of the scalar Ω , for which N 2 2 eigenvalues vanish and only one of them happens to be equal to one, since the total sum has to be unit. It may be a cascade of the eigenvalue disappearances (with many different critical times), or synchronized process (with only one critical time t * ), however the final of this process can be characterized by the degenerated matrix H ( a ) ( b ) with the rank equal to one. Thus, when Ω > max { Ω { α } * } , the color vector q ( a ) appears, which forms the polarization tensor H ( a ) ( b ) = q ( a ) q ( b ) according to the decomposition (31). The evolution of the color vector q ( a ) is described in the next Section.

3. Phenomenology of the Spontaneous Color Polarization

3.1. Basic Ansatz and Scheme of Analysis

3.1.1. On the Universe Evolution Scenario

We consider the toy model of the Universe expansion with three epochs, which we indicate as symmetric, interim transition, and dissymmetric, respectively, while keeping in mind the analogy with the theory of ferroelectricity.
(a)
The evolution of the Universe during the symmetric epoch is considered to be predetermined by the coupling of the multiplet of vector fields U ( a ) i , non-Abelian Yang–Mills field with the potentials A k ( a ) , and the gravitational field. Objectively, a self-consistent description of this stage is very complicated because of a few reasons. First of all, in order to properly describe the non-Abelian system on this stage we have to use the formalism of the Quantum Field Theory. Second, because some of the vector fields from the set { U ( a ) i } can be the spacelike or null ones, the theory cannot ignore the appearance of ghosts. Third, generally, the model with spacelike and null vector fields is not stable. The presence of instability complicates the mathematical description of the model, but facilitates the understanding of the idea, that a phase transition is possible in the Universe history, designed to stabilize the situation. In principle, in this stage the idea of Dirac could be developed (many thanks to the first Referee for the hint!). Indeed, because the vector fields U ( a ) i are not normalized, one can consider the coincidence (or at least the proportionality) of the gauge field potentials A i ( a ) with U i ( a ) .
(b)
The scenario of the transition stage includes three processes: first, the Yang–Mills field becomes quasi-Abelian, A i ( a ) Q ( a ) A i , decouples from the interaction with the color vector fields and degrades A i 0 ; second, the procedure of color polarization reduces the color multiplet U ( a ) i to one timelike unit vector U ( a ) i q ( a ) U i ; and third, the color vector q ( a ) precesses and lines up along the constant color vector Q ( a ) . Of course, such a scenario raises a lot of questions. One of them is connected with the origin of the parallel alignment in the group space of the gauge and vector fields. Another question requires to precise the term spontaneous symmetry breaking. We will try to discuss them below.
(c)
The dissymmetric epoch is assumed to be characterized by unit vector field U i , associated with the canonic aether velocity four-vector; some aspects of the corresponding model are already studied in the framework of isotropic Friedmann type cosmology.
In this paper we focus on the description of the transition stage only. The anisotropic spatially homogeneous Bianchi-I model is considered to be the spacetime platform for three indicated cosmological epochs. We monitor the state of the Universe using the color polarization tensor H ( a ) ( b ) ( t ) as a function of cosmological time. In the symmetric epoch, the rank of the matrix H ( a ) ( b ) is assumed to be equal to N 2 1 , and the determinant H = det H ( a ) ( b ) to be non-vanishing. In the dissymmetric epoch, the corresponding rank is equal to one, the determinant is equal to zero. We assume that during the transition epoch N 2 2 eigenvalues of the color polarization tensor (from the total N 2 1 ones) take zero values.

3.1.2. On the Parallel Fields in the Group Space

We assume that the mechanism of the color polarization can be adequately described in terms of fields “parallel” in the group space. This term appeared in the theory of gauge fields for description of the quasi-Abelian gauge field configuration, for which A i ( a ) = Q ( a ) A i with constants Q ( a ) satisfying the condition Q ( a ) Q ( a ) = 1 (see the work [62]). Initially, this ansatz was chosen to facilitate the analysis of the Einstein–Yang–Mills equations, and for analysis of colored particle motion (see, e.g., [63,64]. Indeed, for the parallel potentials, all of the nonlinear terms in (7) and in the Yang–Mills equations disappear, thus simplifying the model equations. However, we would like to look at the problem from the other side. Let the gauge field be strong, so that the term G f ( b ) ( c ) ( a ) A m ( b ) A n ( c ) makes a significant contribution into the field strength (7). When A i ( a ) = Q ( a ) A i , this term converts into G f ( b ) ( c ) ( a ) Q ( b ) Q ( c ) A m A n 0 , i.e., it vanishes, even if the potential is huge. Because the energy density of the gauge field contains the square of the F m n ( a ) , one can see that the energy of the non-Abelian configuration has to be bigger than the energy of the one with the parallel fields. From the standard physical point of view, in the process of evolution the system tends to go into the state, for which the energy is less. In other words, the state with parallel potentials may be more preferable for the system and, thus, the vector field alignment is a good mechanism for the minimization of its energy.
Now we apply this ansatz to the vector field assuming that it converts into the product
U ( a ) i = q ( a ) U i ,
but generally the color vector q ( a ) is not constant and it does not coincide with the quantity Q ( a ) . The gauge covariant derivative can be now transformed into
D ^ m U n ( a ) = q ( a ) m U n + U n m q ( a ) + G f ( b ) ( c ) ( a ) Q ( b ) q ( c ) A m U n .
Clearly, when q ( a ) = Q ( a ) , the potential four-vector A i drops out from the derivative D ^ m U n ( a ) , since the structure constant are skew-symmetric. Because the derivative (34) enters the expression of the energy density of the vector field, we can use the same argument: the last term in (34) vanishes, when q ( a ) = Q ( a ) , even if the potential A i is huge. In other words, in the process of evolution the color vector q ( a ) has to tend to the color vector Q ( a ) in order to minimize the total energy of the system.
Keeping in mind the relationship (15) and assuming that the four-vector U i is unit and time-like
g i k U i U k = 1 ,
we obtain that the vector q ( a ) in the group space is also unit
G ( a ) ( b ) q ( a ) q ( b ) = 1 .
Clearly, the gauge-covariant derivative D ^ m q ( a ) is orthogonal to the color vector q ( a ) , i.e.,
0 = 1 2 D ^ m [ q ( a ) q ( a ) ] = q ( a ) D ^ m q ( a ) = q ( a ) [ m q ( a ) + G f ( b ) ( c ) ( a ) A m Q ( b ) q ( c ) ] = q ( a ) m q ( a ) = 0 .
The tensor of color polarization H ( a ) ( b ) now takes the form
H ( a ) ( b ) = q ( a ) q ( b ) .
Clearly, the rank of the matrix H ( b ) ( a ) is equal to one, and the determinant is equal to zero. The color vector q ( a ) is an eigenvector of H ( b ) ( a ) with the unit eigenvalue:
H ( b ) ( a ) q ( b ) = q ( a ) .
In other words, the mechanism of parallelization in the group space applied to the color vectors U ( a ) i gives the result that is analogous to the polarization of light in optics. Indeed, one can say that a special direction in the group space appears, which is associated with the color vector q ( a ) , along which all of the color vectors U ( a ) i are lined up. Generally, q ( a ) ( t ) precesses in the group space, but when the color vector q ( a ) becomes constant, one can transform the matrix H ( a ) ( b ) into diag { 1 , 0 , 0 0 } for arbitrary time moment, using the admissible rotation in the group space. The description of the evolution of the color vector q ( a ) ( t ) is the crucial point of our analysis.

3.1.3. On the Spacetime Platform

For the illustration of the suggested idea, we consider the class of homogeneous spacetimes with the metric
d s 2 = d t 2 a 2 ( t ) d x 2 b 2 ( t ) d y 2 c 2 ( t ) d z 2 .
Generally, this metric describes the Bianchi-I model. When a ( t ) = b ( t ) = c ( t ) , we obtain the Friedmann type model. We assume that all of the unknown model state functions inherit the spacetime symmetry and depend on the cosmological time only.
The global unit timelike vector U i is assumed to be of the form U i = δ t i ; below, we show that the model admits this construction. The metric (40) provides the covariant derivative to be symmetric, the acceleration four-vector a i to vanish, and the scalar of expansion Θ to have very simple form
m U n = n U m = 1 2 g ˙ m n , a i U k k U i = 0 , Θ ( t ) k U k = a ˙ a + b ˙ b + c ˙ c .
Here and below the dot denotes the ordinary derivative with respect to time. In addition, we use two differential consequences of the normalization conditions: U k i U k = 0 and q ( c ) U k k q ( c ) = 0 . The listed assumptions simplify the model essentially, and we start to discuss its details.

3.2. Reduced Master Equations

3.2.1. The Yang–Mills Field in the Transition Epoch

As we already noticed, we assume that the first manifestation of the phase transition in the model under discussion is connected with the parallelization of the Yang–Mills potentials, i.e., A i ( a ) Q ( a ) A i , where Q ( a ) is a constant vector in the group space; the second manifestation is the parallelization of the color vector fields U i ( a ) q ( a ) U i . For such fields, parallel in the group space, the strength of the Yang–Mills field (7) has the quasi-Maxwellian form
F m n ( a ) = Q ( a ) [ m A n n A m ] .
Based on the Bianchi-I spacetime platform and taking into account the antisymmetric properties of the structure constants, we find the color current Γ ( a ) i that is produced by the color vector fields:
Γ ( a ) i = G κ ( C 1 + C 2 + C 3 ) U i f ( d ) ( c ) ( a ) q ( c ) [ q ˙ ( d ) G f ( d ) ( b ) ( h ) Q ( b ) q ( h ) U m A m ]
G 2 κ C 1 A k Δ k i f ( c ) ( a ) ( d ) q ( c ) f ( d ) ( b ) ( h ) Q ( b ) q ( h ) .
Here, Δ k i δ k i U i U k is the projector. Clearly, as long as the color current is not vanishing, Γ ( a ) i 0 , or is not linear in the potential A i , the Yang–Mills Equation (18) do not admit the trivial solutions A i = 0 . In principle, this color current disappears, when q ( a ) = Q ( a ) = c o n s t ; we consider this state as the final one. When q ( a ) depends on time and, thus, q ˙ ( a ) 0 , one can obtain the trivial solution A i = 0 , if to use the special condition for the Jacobson’s coupling constants
C 1 + C 2 + C 3 = 0 ,
which was motivated, e.g., in [54]. For this case, the color current is proportional to the potential A k , and the trivial exact solution A i = 0 exists. This means that the system can choose this trivial solution for the Yang–Mills field in the bifurcation associated with the end of the first phase transition. Of course, the condition (44) narrows down the frameworks of the canonic Einstein–aether theory; however, this restricted model is very illustrative for our purpose. By the way, if we accept that C 1 + C 3 0 (see the Section 1.1 for motivation), the condition (44) gives C 2 0 , i.e., S ( 0 ) 2 = 0 and the spin-0 aether modes can not propagate.

3.2.2. Solutions to the Equations for the Color Vector Fields

Now we can consider the reduced master equations for the multiplet of vector fields with the conditions A ( a ) i = 0 and C 1 + C 2 + C 3 = 0 . Clearly, the vectorial quantity I ( a ) j (22) vanishes, since the following vector is equal to zero:
U ( b ) n D ^ n U m ( b ) = q ( b ) q ( b ) U n n U m + U m q ( b ) q ˙ ( b ) .
The tensor J ( a ) i j (21) happens to be proportional to the coupling constant C 2 :
J ( a ) i j = C 2 { q ( a ) [ g i j Θ ( i U j ) ] + q ˙ ( a ) Δ i j } .
Subsequently, we see that 3 ( N 2 1 ) equations from the 4 ( N 2 1 ) ones, written in (20) convert into the trivial identities 0 = 0 ; only N 2 1 equations for j = 0 have the nontrivial form
q ( a ) C 2 ( a ¨ a + b ¨ b + c ¨ c ) = λ q ( a ) .
All of these equations are satisfied simultaneously, if the Lagrange parameter λ ( t ) is chosen to be of the form:
λ = C 2 ( a ¨ a + b ¨ b + c ¨ c ) .
We have to emphasize that, in the presented model, the master equations for the multiplet of vector fields do not impose any restrictions on the color vector q ( a ) ( t ) . This means that we can suggest our version of the model of the color vector q ( a ) ( t ) evolution.

3.2.3. Gravity Field Equations

To explain the details of the structure of the reduced aether stress-energy tensor (25) we have to make four preliminary remarks.
  • Using (46) and the relationships
    i ( q ( a ) U j ) = 1 2 q ( a ) g ˙ i j + q ˙ a U i U j
    we can calculate the scalar
    K ( a ) ( b ) i j m n D ^ i U m ( a ) D ^ j U n ( b ) = C 2 [ Θ 2 i U j i U j ] = 2 C 2 [ a ˙ a b ˙ b + a ˙ a c ˙ c + b ˙ b c ˙ c ] .
  • Because the tensor i U j is symmetric, the term proportional to C 1 in (25) vanishes.
  • Because U k k U j = 0 and q ( a ) q ˙ ( a ) = 0 , the term with C 4 in (25) also vanishes.
  • Because J i j ( a ) = J j i ( a ) , the last term in (25) converts into
    C 2 [ g p q ( Θ ˙ + Θ 2 ) + 1 2 Θ g ˙ p q + U m m ( p U q ) ] .
Thus, for the tensor T q p we obtain the formula, which surprisingly contains neither the color vector q ( a ) nor its derivative q ˙ ( a ) :
T q p = C 2 { δ q p ( a ¨ a + b ¨ b + c ¨ c + a ˙ a b ˙ b + a ˙ a c ˙ c + b ˙ b c ˙ c ) + δ 1 p δ q 1 ( a ¨ a + a ˙ a b ˙ b + a ˙ a c ˙ c ) + δ 2 p δ q 2 ( b ¨ b + b ˙ b c ˙ c + a ˙ a b ˙ b ) + δ 3 p δ q 3 ( c ¨ c + a ˙ a c ˙ c + b ˙ b c ˙ c ) } .
Taking into account the term for λ (48), we can now represent four non-trivial gravity field equations in the following form:
[ a ˙ a b ˙ b + a ˙ a c ˙ c + b ˙ b c ˙ c ] ( 1 + C 2 ) = Λ , [ b ¨ b + c ¨ c + b ˙ b c ˙ c ] ( 1 + C 2 ) = Λ , [ a ¨ a + c ¨ c + a ˙ a c ˙ c ] ( 1 + C 2 ) = Λ , [ b ¨ b + a ¨ a + b ˙ b a ˙ a ] ( 1 + C 2 ) = Λ .
It is the well known system of equations, the novelty is only in the structure of the constant: we deal now with the modified cosmological constant Λ 1 + C 2 instead of Λ . The symmetry of the Equation (53) allows us to explicitly find the expansion scalar Θ ( t ) as the function of time. Indeed, if we calculate the auxiliary function Θ ˙ + Θ 2 , we obtain
Θ ˙ + Θ 2 = a ¨ a + b ¨ b + c ¨ c + 2 [ a ˙ a b ˙ b + a ˙ a c ˙ c + b ˙ b c ˙ c ] .
The appropriate linear combination of the Equation (53) shows that the right-hand side of (54) is equal to 3 Λ 1 + C 2 . In other words, we are faced with the equation
Θ ˙ + Θ 2 = 3 Λ 1 + C 2 ,
the solution to which is the following:
Θ ( t ) = 3 H 0 { Θ ( t 0 ) + 3 H 0 th [ 3 H 0 ( t t 0 ) ] Θ ( t 0 ) th [ 3 H 0 ( t t 0 ) ] + 3 H 0 } ,
where the parameter H 0 is introduced as
H 0 Λ 3 ( 1 + C 2 ) .
The function Θ ( t ) starts from the value Θ ( t 0 ) at t = t 0 and finishes with the value 3 H 0 at t . When | Θ ( t 0 ) | < 3 H 0 , the function Θ ( t ) is regular in the interval t 0 < t < , and it is monotonic with Θ ˙ > 0 . When Θ ( t 0 ) = 3 H 0 , the solution is constant, Θ ( t ) = 3 H 0 .
Additionally, one can find the explicit formula for the unit volume in the three-dimensional space
V ( t ) = a ( t ) · b ( t ) · c ( t ) a ( t 0 ) b ( t 0 ) c ( t 0 ) ,
for which V ( t 0 ) = 1 , and
V ˙ V = a ˙ a + b ˙ b + c ˙ c = Θ ( t ) = k U k .
Using (56) and (59), we now obtain
V ( t ) = Θ ( t 0 ) 3 H 0 sh [ 3 H 0 ( t t 0 ) ] + ch [ 3 H 0 ( t t 0 ) ] .
Finally, it is clear that the differences of logarithmic derivatives have the form
b ˙ b a ˙ a = K 1 V ( t ) , c ˙ c b ˙ b = K 2 V ( t ) , a ˙ a c ˙ c = K 3 V ( t ) ,
where K 1 , K 2 , K 3 are the integration constants, which satisfy the condition K 1 + K 2 + K 3 = 0 . To make sure that (61) take place, we can check, e.g., due to the second and third equations from the set (53)
d d t ( b ˙ b a ˙ a ) = ( b ˙ b a ˙ a ) Θ ,
and then integrate this equation. Asymptotically, at V , the rates of expansion a ˙ a , b ˙ b and c ˙ c coincide, i.e., the Universe tends to be spatially isotropic. If we are interested to find the functions a ( t ) , b ( t ) , c ( t ) , we can use the relationships
a ˙ a = 1 3 Θ + K 3 K 1 3 V , b ˙ b = 1 3 Θ + K 1 K 2 3 V , c ˙ c = 1 3 Θ + K 2 K 3 3 V ,
supplemented by (56) and (60), and then integrate them. Finally, we can state the following.
(i)
When the cosmological constant is equal to zero, Λ = 0 , and C 2 1 , the Equation (53) give the known anisotropic Kasner solutions, and the parameter C 2 , associated with the aether, becomes hidden.
(ii)
When Λ 0 , we deal with the spacetime isotropization, and we asymptotically obtain the de Sitter spacetime with a ( t ) exp { Λ 3 ( 1 + C 2 ) t } ; in this case the coupling constant C 2 modifies the rate of expansion.
(iii)
When Λ = 0 and C 2 = 1 simultaneously, the gravity field Equation (53) are satisfied identically, so that the metric functions a ( t ) , b ( t ) and c ( t ) are arbitrary (in the presented toy model).
(iv)
When C 3 = C 1 and C 2 = 0 , the parameters C 1 and C 4 happen to be hidden, and there are no imprints of the dynamic aether in the equations for the gravity field describing the Universe evolution (dynamic aether becomes invisible?).

3.3. Dynamics of the Color Vector q ( a ) ( t )

In the scenario proposed above, the minimization of the energy of the gauge field can be reached, if A i ( a ) = Q ( a ) A i ; the minimization of the interaction energy between the aether and the gauge field requires the color vector q ( a ) to coincide with Q ( a ) . However, during the first stage of the phase transition, we deal with the relationship q ( a ) Q ( a ) , and one needs to describe the mechanism of re-orientation of the color vector q ( a ) . In this subsection we propose our version of the color vector evolution; we indicate it as the precession mechanism.

3.3.1. Exact Solution to the Precession Equation

According to our scenario, after the phase transition, when the specific direction Q ( a ) in the group space appeared, and the color vector fields have become parallel, U ( a ) i = q ( a ) U i , the process of precession of the color vector q ( a ) around the color vector Q ( a ) was started up. This process can be described by the equation, which is the generalization of the Wong equation for the color charges [65]
U m D ^ m q ( a ) = F ( a ) , F ( a ) = ν [ δ ( b ) ( a ) q ( a ) q ( b ) ] Q ( b ) .
In the Bianchi-I spacetime platform, these equations can be presented in more detailed form:
q ˙ ( a ) + G f ( b ) ( c ) ( a ) Q ( b ) q ( c ) ( U m A m ) = ν ( t ) [ Q ( a ) q ( a ) q ( b ) Q ( b ) ] ,
where ν ( t ) is some function of time. The force-like term F ( a ) in (64) is orthogonal to the color vector q ( a ) : i.e., q ( a ) F ( a ) = 0 . This property of the force-like term F ( a ) guaranties that the normalization condition q ( a ) q ( a ) = 1 is supported. The second term in the left-hand side of the Equation (65) is also orthogonal to the color vector q ( a ) ; it can be eliminated using the so-called Landau’s gauge condition for the Abelian potential four-vector, U i A i = 0 . The equations, modified correspondingly,
q ˙ ( a ) = ν ( t ) [ Q ( a ) q ( a ) q ( b ) Q ( b ) ] ,
look like the equations of motion of the macroscopic particle with unit mass under the influence of the Stokes-type force in the fluid flow, which is characterized by the velocity Q ( a ) . In order to solve the Equation (66), we consider the scalar product X ( t ) q ( a ) ( t ) Q ( a ) ; since q ( a ) ( t ) q ( a ) ( t ) = 1 and Q ( a ) Q ( a ) = 1 , the scalar product X is equal to the cosine of the angle between the vectors q ( a ) and Q ( b ) in the group space. Then we find its derivative X ˙ = Q ( a ) F ( a ) , and obtain the equation, which is very famous in the theory of differential equations:
X ˙ = ν ( t ) ( 1 X 2 ) .
This equation admits two special constant solutions X = ± 1 , which relate to the cases q ( a ) ( t ) = ± Q ( a ) , i.e., q ( a ) and Q ( a ) are parallel or anti-parallel, respectively, for arbitrary moment of the cosmological time. For arbitrary initial value X ( t 0 ) , the solution to the Equation (67) is known to be of the following form:
X ( t ) = e 2 T ( t ) ( 1 + X ( t 0 ) ) ( 1 X ( t 0 ) ) e 2 T ( t ) ( 1 + X ( t 0 ) ) + ( 1 X ( t 0 ) ) , T ( t ) t 0 t d τ ν ( τ ) .
Returning to the Equation (66) rewritten in the form
d q ( a ) d T = Q ( a ) q ( a ) X ( T ) ,
we readily find the solution
q ( a ) ( t ) = q ( a ) ( t 0 ) [ ch T + X ( t 0 ) sh T ] + Q ( a ) [ sh T + X ( t 0 ) ( ch T 1 ) ] [ ch T + X ( t 0 ) sh T ] .
For illustration, it is convenient to consider the case with X ( t 0 ) = 0 ; it corresponds to the configuration with initially orthogonal color vectors q ( a ) ( t 0 ) and Q ( a ) ; now we obtain
X ( t ) = th T , q ( a ) ( t ) = q ( a ) ( t 0 ) ch T + Q ( a ) th T ,
H ( a ) ( b ) ( t ) = H ( a ) ( b ) ( t 0 ) ch 2 T + Q ( a ) Q ( b ) th 2 T + [ Q ( a ) q ( b ) ( t 0 ) + Q ( b ) q ( a ) ( t 0 ) ] sh T ch 2 T .
Asymptotically, at T we obtain q ( a ) Q ( a ) . In other words, during the evolution process, the color vector q ( a ) forgets its initial value q ( a ) ( t 0 ) and it is forced to be lined up along the direction Q ( a ) . The question arises: what is the rate of asymptotic parallelization of the vectors q ( a ) and Q ( d ) . Clearly, this process is predetermined by the properties of the function ν ( t ) ; let us consider an example of its modeling.

3.3.2. Model Function ν ( t )

For the modeling of the function ν ( t ) , we use the following construction:
ν ( t ) = 2 γ T * Θ ( t ) V 2 [ V * 2 V 2 ( t ) ] 1 + γ ,
and calculate the function T ( t ) (68):
T ( t ) = T * { 1 [ V * 2 V 2 ( t ) ] γ 1 [ V * 2 1 ] γ } .
Here T * and γ are positive parameters. Let us make three remarks concerning this formula.
  • When the function V ( t ) grows, it reaches the value V * , so, at the finite moment of the cosmological time t = t * the function T ( t ) reaches infinity, T ( t * ) = .
  • According to (70), at T ( t * ) = , the color vector q ( a ) ( t * ) coincides with Q ( a ) , i.e., the process of their parallelization finished during the finite time interval.
  • At the moment t = t * , the derivative of the color vector q ( a ) takes zero value, q ˙ ( a ) ( t * ) = 0 .
In order to illustrate the last statement, we calculate the derivative of the function q ( a ) ( t ) presented in (71) for the case X ( t 0 ) = 0 :
q ˙ ( a ) ( t ) = ν ( t ) ch T [ Q ( a ) ch T q ( a ) ( t 0 ) th T ] .
For big values of the function T, one can use the following replacement (compare (73) and (74)):
ν ( t ) 2 γ T * 1 γ Θ V 2 ( t * ) T 1 + 1 γ ,
where Θ (56) is the regular finite function. Because, for arbitrary α , the limit lim T [ T α ch T ] is equal to zero, we can confirm the condition q ˙ ( a ) ( t * ) = 0 . The critical time moment t * , for which V ( t * ) = V * can be found as
t * = t 0 + 1 3 H 0 ln { 3 H 0 V * 3 H 0 + Θ ( t 0 ) [ 1 + 1 + Θ 2 ( t 0 ) 9 H 0 2 9 H 0 2 V * 2 ] } .
One can check that t * > t 0 , when V * > 1 . For the special case, when 3 H 0 = Θ ( t 0 ) = Θ ( t ) , we obtain t * = t 0 + 1 3 H 0 ln V * .

3.3.3. Model Threshold Function Ω

Let us now consider the function Ω introduced by (13) and (14), the value of which indicates the nearness of the event of the spontaneous color polarization. We have now that
Ω ( a ) = q ˙ ( a ) + q ( a ) Θ , Ω 2 = q ˙ ( a ) q ˙ ( a ) + Θ 2 ,
thus, for the illustrative function (71) we obtain
Ω = ν 2 ch 2 T + Θ 2 .
Taking into account (73) and (74), we finally reconstruct the function Ω as
Ω = Θ 1 + 4 γ 2 T * 2 V 4 ch 2 T [ T T * + 1 ( V * 2 1 ) γ ] 2 ( 1 + γ ) γ .
Near the critical point t = t * , when T , the threshold function Ω behaves as the expansion scalar Θ , i.e., the summarized rate of the Universe expansion Θ = a ˙ a + b ˙ b + c ˙ c predetermines the rate of the onset of the spontaneous color polarization.

4. Discussion and Conclusions

We share the point of view of many cosmologists, that a plethora of fields and particles in the early hot Universe were lost, destroyed, or transformed as a result of a series of phase transitions; now, we observe a limited number of fundamental fields and stable elementary particles. Our goal is to reconstruct the history of emergence of the cosmic substratum indicated as dynamic aether; to be more precise, we are interested to establish a model of formation of the unit timelike vector field as a splinter of the family of vector fields from the SU(N) symmetric multiplet. Our phenomenological model is based on three elements.
(1)
The first element of the model is the procedure of the spontaneous polarization of the SU(N) multiplet of vector fields. This procedure is described in terms of critical behavior of the eigenvalues of the color polarization tensor (see (27), (31) and (32)). When the Universe expands, and the scalar Ω reaches the critical value Ω { α } * , the corresponding eigenvalue with the serial number α vanishes. When Ω > max { Ω { α } * } , we obtain the situation with only one non-vanishing eigenvalue; this last eigenvalue is equal to one, because the trace of the polarization tensor is unit. Based on the analogy with optics, we have to state now, that we deal with linearly polarized system, i.e., the color polarization took place. Mathematically, this procedure means that all of the vectors from the SU(N) symmetric multiplet have become parallel in the group space to one vector U i , namely, U ( a ) i = q ( a ) U i .
(2)
The second element of the model is the mechanism of precession of the color vector q ( a ) around the constant color vector Q ( a ) formed in course of parallelization of the Yang–Mills potential A i ( a ) Q ( a ) A i . Evolution of the color vector q ( a ) was described by the equation of the Wong type (64); the exact solution to this equation is presented in (70). Asymptotically, q ( a ) ( t ) tends to the constant Q ( a ) , and the system as a whole becomes quasi-Abelian, since all of the structure constants f ( a ) ( b ) ( c ) fall out from the basic formulas.
(3)
The third element of the model is the calculation of the rate of the precession of the color vector q ( a ) . We have chosen the guiding function ν ( t ) in (65) in the critical form (73) providing the process of parallelization of q ( a ) and Q ( a ) to be finished during the finite time interval (see (70) and (74)).
In the context of our study, we can not elude the question concerning the breaking of a symmetry; what is the symmetry, which we mean? In our scenario we assume that a phase transition of the second kind takes place in the early Universe. According to the Landau’s theory of the phase transitions in crystals, one can indicate explicitly the spatial symmetry groups, which characterize the system before and after the transition, and one can point out a specific order parameter (the scalar, spatial vector, or tensor). In our case, the spacetime symmetry in the Universe is not broken, and the restructuring takes place in the field conglomerate and concerns the internal symmetry of the self-interacting system. Before such an effective phase transition the set of interacting fields could be characterized by the SU(N) symmetry group, and after the transition these fields display the symmetry, which relates to the U(1) group of internal symmetry. In other words, the group space is the scene of the transition; the transformations A i ( a ) Q ( a ) A i and U i ( a ) q ( a ) U i relate to first act of the spectacle, as a result the group space becomes bi-axial, because two color vectors do not coincide, q ( a ) Q ( a ) . During the second act, the color vector q ( a ) precesses and tends to be equal to Q ( a ) ; at the final of the second act the group space becomes uni-axial; the corresponding color vector plays the role of the order parameter of this effective phase transition. What is the driven force of this phase transition? We think that it is the natural tendency of the physical system to find the state with minimal energy level; since the energy of the quasi-Abelian configuration is less than the energy of the non-Abelian one, the transition is inevitable.
The conclusion of our study is the following: the aether velocity four-vector can be considered as a rudiment, which remained after the spontaneous color polarization accompanying the phase transition in the early Universe.

Author Contributions

Conceptualization, A.B.; formal analysis, A.B. and G.K.; investigation, A.B. and G.K.; writing—original draft preparation, A.B.; writing—review and editing, A.B.; All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Russian Foundation for Basic Research (Project No. 20-02-00280).

Acknowledgments

The work was supported by Russian Foundation for Basic Research (Project No. 20-02-00280), and, partially, by the Program of Competitive Growth of Kazan Federal University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jacobson, T.; Mattingly, D. Gravity with a dynamical preferred frame. Phys. Rev. D 2001, 64, 024028. [Google Scholar] [CrossRef] [Green Version]
  2. Jacobson, T.; Mattingly, D. Einstein-aether waves. Phys. Rev. D 2004, 70, 024003. [Google Scholar] [CrossRef] [Green Version]
  3. Jacobson, T. Einstein-aether gravity: A status report. arXiv 2007, arXiv:0801.1547. [Google Scholar]
  4. Nojiri, S.; Odintsov, S.D. Unified cosmic history in modified gravity: From F(R) theory to Lorentz non-invariant models. Phys. Rep. 2011, 505, 59–144. [Google Scholar] [CrossRef] [Green Version]
  5. Nojiri, S.; Odintsov, S.D.; Oikonomou, V.K. Modified gravity theories on a nutshell: Inflation, bounce and late-time evolution. Phys. Rep. 2017, 692, 1–104. [Google Scholar] [CrossRef] [Green Version]
  6. Will, C.M. Theory and Experiment in Gravitational Physics; Cambridge University Press: Cambridge, UK, 1993. [Google Scholar]
  7. Abbott, B.P.; Abbott, R.; Abbott, T.D.; Acernese, F.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R.X.; Adya, V.B.; et al. [LIGO Scientific Collaboration and Virgo Collaboration], [Fermi Gamma-Ray Burst Monitor], [INTEGRAL] Gravitational Waves and Gamma-rays from a Binary Neutron Star Merger: GW170817 and GRB 170817A. Astrophys. J. Lett. 2017, 848, L13. [Google Scholar] [CrossRef]
  8. Elliott, J.W.; Moore, G.D.; Stoica, H. Constraining the New Aether: Gravitational Cherenkov Radiation. J. High Energy Phys. 2005, 0508, 066. [Google Scholar] [CrossRef]
  9. Kostelecky, A.; Tasson, J.D. Constraints on Lorentz violation from gravitational Cherenkov radiation. Phys. Lett. B 2015, 749, 551–559. [Google Scholar] [CrossRef] [Green Version]
  10. Oost, J.; Mukohyama, S.; Wang, A. Constraints on Einstein-aether theory after GW170817. Phys. Rev. D 2018, 97, 124023. [Google Scholar] [CrossRef] [Green Version]
  11. Trinh, D.; Pace, F.; Battye, R.A.; Bolliet, B. Cosmologically viable generalized Einstein-Aether theories. Phys. Rev. D 2019, 99, 043515. [Google Scholar]
  12. Balakin, A.B.; Shakirzyanov, A.F. Axionic extension of the Einstein-aether theory: How does dynamic aether regulate the state of axionic dark matter? Phys. Dark Univ. 2019, 24, 100283. [Google Scholar] [CrossRef] [Green Version]
  13. Eling, C.; Jacobson, T. Black holes in Einstein-aether theory. Class. Quantum Gravity 2006, 23, 5643. [Google Scholar] [CrossRef] [Green Version]
  14. Foster, B.Z. Noether charges and black hole mechanics in Einstein-aether theory. Phys. Rev. D 2006, 73, 024005. [Google Scholar] [CrossRef] [Green Version]
  15. Tamaki, T.; Miyamoto, U. Generic features of Einstein-Aether black holes. Phys. Rev. D 2008, 77, 024026. [Google Scholar] [CrossRef] [Green Version]
  16. Eling, C.; Jacobson, T. Spherical solutions in Einstein-aether theory: Static aether and stars. Class. Quantum Gravity 2006, 23, 5625–5642. [Google Scholar] [CrossRef] [Green Version]
  17. Eling, C.; Jacobson, T.; Miller, M.C. Neutron stars in Einstein-aether theory. Phys. Rev. D 2007, 76, 042003. [Google Scholar] [CrossRef] [Green Version]
  18. Foster, B.Z. Strong field effects on binary systems in Einstein-aether theory. Phys. Rev. D 2007, 76, 084033. [Google Scholar] [CrossRef] [Green Version]
  19. Ding, C.; Wang, A.; Wang, X.; Zhu, T. Hawking radiation of charged Einstein-aether black holes at both Killing and universal horizons. Nucl. Phys. B 2016, 913, 694–715. [Google Scholar] [CrossRef] [Green Version]
  20. Pacilio, C.; Liberati, S. An improved derivation of the Smarr Formula for Lorentz-breaking gravity. Phys. Rev. D 2017, 95, 124010. [Google Scholar] [CrossRef] [Green Version]
  21. Bhattacharjee, M.; Mukohyama, S.; Wan, M.-B.; Wang, A. Gravitational collapse and formation of universal horizons in Einstein-aether theory. Phys. Rev. D 2018, 98, 064010. [Google Scholar] [CrossRef] [Green Version]
  22. Will, C.M.; Nordtvedt, K. Conservation laws and preferred frames in relativistic gravity. I. Preferred-frame theories and an extended PPN formalism. Astrophys. J. 1972, 177, 757–774. [Google Scholar] [CrossRef]
  23. Nordtvedt, K.; Will, C.M. Conservation laws and preferred frames in relativistic gravity. II. Experimental evidence to rule out preferred-frame theories of gravity. Astrophys. J. 1972, 177, 775–792. [Google Scholar] [CrossRef]
  24. Dirac, P.A.M. A New classical theory of electrons (Part I). Proc. Roy. Soc. A 1951, 209, 291–296. [Google Scholar]
  25. Dirac, P.A.M. A New classical theory of electrons (Part II). Proc. Roy. Soc. A 1952, 212, 330–339. [Google Scholar]
  26. Dirac, P.A.M. A New classical theory of electrons (Part III). Proc. Roy. Soc. A 1954, 223, 438–445. [Google Scholar]
  27. Bertolami, O.; Paramos, J. The flight of the bumblebee: Vacuum solutions of a gravity model with vector-induced spontaneous Lorentz symmetry breaking. Phys. Rev. D 2005, 72, 044001. [Google Scholar] [CrossRef] [Green Version]
  28. Seifert, M.D. Generalized bumblebee models and Lorentz-violating electrodynamics. Phys. Rev. D 2010, 81, 065010. [Google Scholar] [CrossRef] [Green Version]
  29. Escobar, C.A.; Martín-Ruiz, A. On the equivalence between bumblebee models and electrodynamics in a non-linear gauge. Phys. Rev. D 2017, 95, 095006. [Google Scholar] [CrossRef] [Green Version]
  30. Assunção, J.F.; Mariz, T.; Nascimento, J.R.; Petrov, A.Y. Dynamical Lorentz symmetry breaking in a tensor bumblebee model. Phys. Rev. D 2019, 100, 085009. [Google Scholar] [CrossRef] [Green Version]
  31. Kostelecky, A.; Mewes, M. Electrodynamics with Lorentz-violating operators of arbitrary dimension. Phys. Rev. D 2009, 80, 015020. [Google Scholar] [CrossRef] [Green Version]
  32. Liberati, S.; Maccione, L. Lorentz violation: Motivation and new constraints. Ann. Rev. Nucl. Part. Sci. 2009, 59, 245. [Google Scholar] [CrossRef]
  33. Kostelecky, A.; Russell, N. Data tables for Lorentz and CPT violation. Rev. Mod. Phys. 2011, 83, 11. [Google Scholar] [CrossRef] [Green Version]
  34. Zlosnik, T.G.; Ferreira, P.G.; Starkman, G.D. Modifying gravity with the Aether: An alternative to Dark Matter. Phys. Rev. D 2007, 75, 044017. [Google Scholar] [CrossRef] [Green Version]
  35. Battye, R.A.; Pace, F.; Trinh, D. Cosmological perturbation theory in Generalized Einstein-Aether models. Phys. Rev. D 2017, 96, 064041. [Google Scholar] [CrossRef] [Green Version]
  36. Jacobson, T. Extended Horava gravity and Einstein-aether theory. Phys. Rev. D 2010, 81, 101502. [Google Scholar] [CrossRef]
  37. Barausse, E.; Jacobson, T.; Sotiriou, T.P. Black holes in Einstein-aether and Horava-Lifshitz gravity. Phys. Rev. D 2011, 83, 124043. [Google Scholar] [CrossRef] [Green Version]
  38. Heinicke, C.; Baekler, P.; Hehl, F.W. Einstein-aether theory, violation of Lorentz invariance, and metric-affine gravity. Phys. Rev. D 2005, 72, 025012. [Google Scholar] [CrossRef] [Green Version]
  39. Ranjit, C.; Debnath, U. Reconstruction of Einstein-Aether Gravity from other Modified Gravity Models. Eur. Phys. J. Plus 2014, 129, 235. [Google Scholar] [CrossRef]
  40. Carroll, S.M.; Tam, H. Aether Compactification. Phys. Rev. D 2008, 78, 044047. [Google Scholar] [CrossRef] [Green Version]
  41. Pujolas, O.; Sibiryakov, S. Supersymmetric Aether. J. High Energy Phys. 2012, 01, 062. [Google Scholar] [CrossRef] [Green Version]
  42. Farias, C.F.; Lehum, A.C.; Nascimento, J.R.; Petrov, A.Y. On the superfield supersymmetric aether-like Lorentz-breaking models. Phys. Rev. D 2012, 86, 065035. [Google Scholar] [CrossRef] [Green Version]
  43. Lehum, A.C.; Nascimento, J.R.; Petrov, A.Y.; da Silva, A.J. Supergauge theories in aether superspace. Phys. Rev. D 2013, 88, 045022. [Google Scholar] [CrossRef] [Green Version]
  44. Kanno, S.; Soda, J. Lorentz violating inflation. Phys. Rev. D 2006, 74, 063505. [Google Scholar] [CrossRef] [Green Version]
  45. Donnelly, W.; Jacobson, T. Coupling the inflaton to an expanding aether. Phys. Rev. D 2010, 82, 064032. [Google Scholar] [CrossRef] [Green Version]
  46. Barrow, J.D. Some Inflationary Einstein-Aether Cosmologies. Phys. Rev. D 2012, 85, 047503. [Google Scholar] [CrossRef] [Green Version]
  47. Solomon, A.R.; Barrow, J.D. Inflationary Instabilities of Einstein-Aether Cosmology. Phys. Rev. D 2014, 89, 024001. [Google Scholar] [CrossRef] [Green Version]
  48. Wei, H.; Yan, X.-P.; Zhou, Y.-N. Cosmological Evolution of Einstein-Aether Models with Power-law-like Potential. Gen. Rel. Grav. 2014, 46, 1719. [Google Scholar] [CrossRef] [Green Version]
  49. Alhulaimi, B.; van den Hoogen, R.J.; Coley, A.A. Spatially Homogeneous Einstein-Aether Cosmological Models: Scalar Fields with a Generalized Harmonic Potential. Phys. Rev. D 2017, 96, 104021. [Google Scholar] [CrossRef] [Green Version]
  50. Balakin, A.B.; Lemos, J.P.S. Einstein-aether theory with a Maxwell field: General formalism. Ann. Phys. 2014, 350, 454–484. [Google Scholar] [CrossRef] [Green Version]
  51. Alpin, T.Y.; Balakin, A.B. Birefringence induced by pp-wave modes in an electromagnetically active dynamic aether. Eur. Phys. J. C 2017, 77, 699. [Google Scholar] [CrossRef] [Green Version]
  52. Balakin, A.B.; Andreyanov, A.V. SU(N)—symmetric dynamic aether: General formalism and a hypothesis on spontaneous color polarization. Space Time Fund. Int. 2017, 4, 36–58. [Google Scholar]
  53. Balakin, A.B. Axionic extension of the Einstein-aether theory. Phys. Rev. D 2016, 94, 024021. [Google Scholar] [CrossRef] [Green Version]
  54. Alpin, T.Y.; Balakin, A.B. The Einstein-Maxwell-aether-axion theory: Dynamo-optical anomaly in the electromagnetic response. Int. J. Mod. Phys. D 2016, 25, 1650048. [Google Scholar] [CrossRef]
  55. Balakin, A.B. Electrodynamics of a cosmic dark fluid. Symmetry 2016, 8, 56. [Google Scholar] [CrossRef]
  56. Balakin, A.B. The extended Einstein-Maxwell-aether-axion model: Exact solutions for axionically controlled pp-wave aether modes. Mod. Phys. Lett. A 2018, 33, 1850050. [Google Scholar] [CrossRef] [Green Version]
  57. Lines, M.E.; Glass, A.M. Principles and Application of Ferroelectrics and Related Materials; Clarendon Press: Oxford, UK, 1977. [Google Scholar]
  58. Landau, L.D.; Lifshitz, E.M.; Pitaevskii, L.P. Electrodynamics of Continuous Media; Butterworth Heinemann: Oxford, UK, 1960. [Google Scholar]
  59. Rubakov, V. Classical Theory of Gauge Fields; Princeton University Press: Princeton, NJ, USA; Oxford, UK, 2002. [Google Scholar]
  60. Akhiezer, A.I.; Peletminsky, S.V. Fields and Fundamental Interactions; Taylor and Francis: London, UK, 2002. [Google Scholar]
  61. Landau, L.D.; Lifshitz, E.M. The Classical Theory of Fields; Pergamon Press: Oxford, UK, 1971; Volume 2. [Google Scholar]
  62. Yasskin, P.B. Solutions for gravity coupled to massless gauge fields. Phys. Rev. D 1975, 12, 2212–2217. [Google Scholar] [CrossRef]
  63. Galtsov, D.V. Particles and Fields in the Environment of Black Holes; Moscow State University Press: Moscow, Russia, 1986. [Google Scholar]
  64. Balakin, A.B.; Kurbanova, V.R.; Zimdahl, W. Parametric phenomena of the particle dynamics in a periodic gravitational wave field. J. Math. Phys. 2003, 44, 5120–5140. [Google Scholar] [CrossRef] [Green Version]
  65. Wong, S.K. Field and particle equations for the classical Yang-Mills field and particles with isotopic spin. Il Nuovo Cimento 1970, 65, 689–694. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Balakin, A.; Kiselev, G. Spontaneous Color Polarization as A Modus Originis of the Dynamic Aether. Universe 2020, 6, 95. https://doi.org/10.3390/universe6070095

AMA Style

Balakin A, Kiselev G. Spontaneous Color Polarization as A Modus Originis of the Dynamic Aether. Universe. 2020; 6(7):95. https://doi.org/10.3390/universe6070095

Chicago/Turabian Style

Balakin, Alexander, and Gleb Kiselev. 2020. "Spontaneous Color Polarization as A Modus Originis of the Dynamic Aether" Universe 6, no. 7: 95. https://doi.org/10.3390/universe6070095

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop