Anisotropic Generalization of the ΛCDM Universe Model with Application to the Hubble Tension

: I deduce an exact and analytic Bianchi type I solution of Einstein’s field equations, which generalizes the isotropic ΛCDM universe model to a corresponding model with anisotropic expansion. The main point of the article is to present the anisotropic generalization of the ΛCDM universe model in a way suitable for investigating how anisotropic expansion modifies observable properties of the ΛCDM universe model. Although such generalizations of the isotropic ΛCDM universe model have been considered earlier, they have never been presented in this form before. Several physical properties of the model are pointed out and compared with properties of special cases, such as the isotropic ΛCDM universe model. The solution is then used to investigate the Hubble tension. It has recently been suggested that the cosmic large-scale anisotropy may solve the Hubble tension. I consider those earlier suggestions and find that the formulae of these papers lead to the result that the anisotropy of the cosmic expansion is too small to solve the Hubble tension. Then, I investigate the problem in a new way, using the exact solution of the field equations. This gives the result that the cosmic expansion anisotropy is still too small to solve the Hubble tension in the general Bianchi type I universe with dust and LIVE (Lorentz Invariant Vacuum Energy with a constant energy density, which is represented by the cosmological constant) and anisotropic expansion in all three directions—even if one neglects the constraints coming from the requirement that the anisotropy should be sufficiently small so that it does not have any significant effect upon the results coming from the calculations of the comic nucleosynthesis during the first ten minutes of the universe. If this constraint is taken into account, the cosmic expansion anisotropy is much too small to solve the Hubble tension.


Introduction
In the present article, I solve Einstein's field equations and calculate exact analytic expressions for the co-moving volume, the Hubble parameter, and the shear scalar as functions of time for the general Bianchi type I universe model with anisotropy in all three directions.The solution is presented in a form suitable for investigating how anisotropic expansion modifies observable properties of corresponding isotropic universe models.
I subsequently apply these solutions to the question of whether the upper limit of the expansion anisotropy of the universe, as restricted by different types of observations, is large enough to permit a solution of the Hubble tension.
The Hubble tension is the fact that early universe measurements and calculations to determine the value of Hubble's constant and late time measurements have given results with a difference that is larger than the uncertainties of the determined values.A recent review of late time measurements has been given by Adam Riess and co-workers [1] with the result that the Hubble constant is H0 = 73.04 ± 1.04 km s −1 Mpc −1 .The most recent result using supernovae and quasars was announced on 30 August 2023 by T. Liu et al. [2].They found H0 = 73.51± 0.67 km s −1 Mpc −1 .A review of the early universe measurements has been given by the Planck team [3] with the result H0 = (67.4± 0.5) km s −1 Mpc −1 .Further references are found in these articles.
A large number of articles have been published with proposed solutions to the Hubble tension, but so far, no solution has been generally accepted.A review [4] with 709 references has recently been published by Maria Dainotti and co-workers.
I here discuss the proposal by V. Yadav [5] that cosmic anisotropy can solve the problem.Ö. Akarsu et al. [6,7] have given constraints on the anisotropy of a Bianchi type I spacetime extension of the standard ΛCDM universe model.They found that the present value of the anisotropy parameter is restricted to 18 0 10    (95% C.L.) from CMB+Lense data and that the introduction of spatial curvature or anisotropic expansion, or both, on a generalized ΛCDM universe model, does not offer a possible relaxation to the H0 tension.The results of the present analysis, based upon an exact solution of Einstein's field equations, confirm this result.
In Figure 1 and Table 3 of ref. [5] it is indicated that the Hubble tension is solved in the context of an anisotropic Bianchi type I universe model.However it is not explained where this solution comes from.It is tempting to think that it is the anisotropy of the expansion which solves the Hubble tension.In the present paper, I investigate whether this is the case.
The method which is used here to determine the effect of the anisotropy of the expansion upon the value of the Hubble's constant, i.e., the present value of the Hubble parameter, is the following.First a formula giving the evolution of the Hubble parameter for the considered anisotropic universe model is deduced.Then, the anisotropy is assumed to vanish, and the corresponding formula with vanishing anisotropy is written down.These formulae connect the Hubble parameter at an arbitrary point of time with its value at the present time, i.e., the Hubble constant, for an anisotropic and the corresponding isotropic universe model.The initial value is chosen as the value determined from the Planck observations of the temperature fluctuations of the cosmic microwave background radiation.These measurements determined the value of the Hubble parameter at the point of time 380,000 years after the Big Bang.This value is an observed quantity and is independent of whether the cosmic expansion is assumed to be anisotropic or isotropic.Hence, this is the initial value in both the evolution equation of the Hubble parameter in the anisotropic and the isotropic case.Then, these equations are used to calculate the present values of the Hubble parameter, i.e., the values of the Hubble constant, in an anisotropic-and an isotropic universe.Finally, the difference of the values of the Hubble constant with anisotropic and isotropic expansion is calculated and compared with the difference of the values from early and late time measurements, which makes up the Hubble tension.

Anisotropic Generalization of the CDM  Universe Model
The Bianchi universes are universe models which are homogeneous but not necessarily isotropic.There exist nine Bianchi-type universe models.The Bianchi type I universe models are the only universe models of the Bianchi type that reduce to flat, isotropic FLRW universe models in the case of vanishing anisotropy.Hence they are candidates for generalizing the CDM  universe models to universe models permitting anisotropic expansion.
Concerning the possibility that cosmic expansion anisotropy should be able to solve the Hubble tension, the investigations in the present paper show that the anisotropy of our universe, as described in terms of the most general Bianchitype I universe model, is much too small to be able to solve the Hubble tension.Hence describing the expansion anisotropy in terms of any other type of the Bianchi universe models is not expected to lead to a solution of the Hubble tension.
In particular, the Bianchi type V universe models do not contain the CDM  universe models because they are spatially curved.They are not expected to be more realistic universe models than those of the Bianchi type I and will, therefore, not be considered in the present work.

Field Equations for the Bianchi Type I Universe Models
A 55-year-old article [8] by P. T. Saunders deserves to be mentioned as a pioneering work.He considered a general Bianchi type I universe model, allowing for different expansion factors in three orthogonal directions with scale factors , , a t a t a t .Us- ing units with the velocity of light ds dt a t dx a t dy a t dz The average scale factor,   at, and co-moving volume,   The scale factors are normalized so that the co-moving volume at the present time is   0 1 Vt  .The redshift of observed radiation from a source emitting the radiation at a point of time The directional Hubble parameters and the average Hubble parameter are The average deceleration parameter is As usual 0 q  means deceleration of the cosmic expansion and 0 q  means ac- celeration.
It has also been usual to introduce the shear scalar, which is a kinematic quantity representing the anisotropy of the cosmic expansion.Yadav [5] has shown that I will here deduce a relationship first shown by Saunders, following [9], but giving more detailed explanations.With the line element (1) Einstein's field equations can be written in the form [10]   and furthermore, as shown by P. Sarmah and U. D. Goswami [11], where  and p are the sum of the mass-energy densities and pressure, respectively, of the matter, radiation, and dark energy contents of the universe.
Due to Equation (11), it is natural in the Bianchi type I universe models to define an effective gravitational mass density, G  , by There is attractive gravity for Adding Equation (10) gives or Equations ( 10) and (13) give The last two equations lead to Integration gives where i K are integration constants representing the deviation of the directional Hubble parameters from isotropy, and K is defined in the last quation.Integration of this Equation and use of Equation (3) gives Physicists investigating the Bianchi type I universe models have often defined an average expansion anisotropy by Using Equation ( 16), one then obtains Hence, the constant K is a measure of the expansion anisotropy of the universe.The It has lately become more usual to introduce an anisotropy parameter: As seen from Equation (24), the quantities   and A have the same physical in- terpretation, which is not obvious from the definitions ( 6) and (21).In this paper, we shall use the parameter   and not A. However both of these parameters have been much in use and still are, so I have demonstrated their physical identity for the benefit of those who want to read more on the Bianchi type I universe models.
For the present value of   , I shall use , which is the upper value permitted by the Planck observations of the temperature fluctuations of the CMB radiation, as suggested by Akarsu et al. [6,7].
Using Equation ( 6), the anisotropic generalization of Friedmann's first Equation takes the form where , , , are the densities of LIVE, cold matter, radiation, and Zeldovich fluid, respectively?A consequence of Einstein's field equations is the laws of energy and momentum conservation in the form Assuming that there is no transition between matter and dark energy, this gives the equations of continuity for LIVE, cold matter, radiation, and stiff Zeldovich fluid with the Inserting Equation (27) into Equation (24) and using that 3 Note that the anisotropy appears with the same dependency upon the scale factor as a Zeldovich fluid.Hence, the effect upon the expansion of the universe of the Zeldovich fluid is equivalent to that of the anisotropy of the expansion.In the present article, the Zeldovich fluid will, therefore, be neglected.
If the expansion of the universe is anisotropic, the anisotropy would have a dominating influence on the expansion of the universe very early in the history of the universe.
Saunders was interested in those epochs of the universe that it was possible to observe in the sixties-Rather late in the history of the universe.Then, the radiation term was much less than the matter term in the expression (31).Therefore, he chose to integrate this Equation for an era of the universe where the contribution of radiation could be neglected.He found three solutions: one for 0  , one for 0  , and one for 0  .
Since  represents the density of LIVE we are here only interested in the case 0  .In this case, Equation (31) reduces to For later comparison, we shall consider four special cases before this Equation is integrated into the general case.We start by presenting the isotropic ΛCDM universe model, following [9], which is presently used as the standard model of the universe.

The ΛCDM Universe Model
where It follows from Equations ( 4) and (33) that the average Hubble parameter is Using Equations ( 5) and (39), the deceleration parameter is There is a transition from decelerated expansion to accelerated expansion at the point of time, or Hence, the corresponding redshift is   . According to Equations ( 25) and (30), the density of the cold matter decreases with time as There is a transition from a matter-dominated era to a LIVE-dominated era at a point of time . Due to the strong repulsive gravitational effect of the negative pressure of LIVE, as seen from the relativistic expression (13) of the effective gravitational mass density, the transition to accelerated cosmic expansion happens before the universe becomes LIVE-dominated.

The Kasner Universe
The empty anisotropic Bianchi type I universe has 0 0 0 M      and is called the Kasner universe.For this universe 0 1   , and Equation (32) reduces to Integrating with   where is the present age of this universe.Using that Equation ( 20) gives the directional scale factors i a in the form The Hubble parameter of the Kasner universe is The Hubble horizon is a surface around an observer separating an internal region where the expansion velocity is less than the velocity of light relative to the observer from an external region where the expansion velocity is larger than that of light.The radius of the Hubble horizon is giving for the Kasner universe.It follows from Equations ( 5) and ( 49) that the deceleration parameter of the Kasner universe is 2 K q  , which means cosmic expansion deceleration.
This may be somewhat surprising since this universe is empty and the corresponding isotropic universe, the Milne universe, has deceleration equal to zero and hence vanishing cosmic deceleration as expected for an empty universe.Hence, the anisotropy induces a deceleration of the cosmic expansion.
In this connection, it may be noted that the isotropic Milne universe is not a special case of the Bianchi type I universe models because these models have vanishing spatial curvature, while the Milne universe has negative spatial curvature.

The LIVE-Dominated Bianchi Type I Universe
Let us first consider the LIVE-dominated isotropic case with 0 0 0 M      .This is the De Sitter universe.Then, Equation (32) reduces to This universe model does not permit the initial condition   and constant Hubble parameter 0 DS HH  .Next we shall consider the anisotropic LIVE-dominated Bianchi type I universe.This universe model is particularly relevant as a model of the early part of the inflationary era since the relationship (7) implies that even a very small anisotropy at the present time means that the anisotropy may have been great, and even dominating, during the first part of the inflationary era.It was thoroughly described in ref. [9], and here I shall only recapitulate the main properties of this model.

In this model
It follows from Equations ( 4) and (56) that in the LIVE-dominated universe, the Hubble parameter is which is similar to the corresponding formula (39) for the Hubble parameter in the In the calculation of the directional scale factors, i a , from Equation (20), we need the integral The resulting expression for i a can be simplified by introducing the constants Using that and Equation ( 13) with 0 1 V  we find that the constants i p fulfil the relationships This leads to the following expression for the directional scale factors: In agreement with Equation (22) in ref. [9], this LIVE-dominated Bianchi type I universe model has been generalized to include viscosity by Mostafapoor and Grøn [12].
The relationships (61) mean that There are two cases, The second model has scale factors The models initially behave rather differently.Two scale factors of the first model initially vanish, and the third is infinite.Hence, this model starts as a needle singularity and approaches isotropy at late times.One scale factor of the second model initially vanishes, and two are infinite.Thus, this model begins as a plane singularity with vanishing thickness and with infinitely great extension.But like the first model it approaches isotropy at late times.This indicates that very different initial behavior at the beginning of the inflationary era would not lead to corresponding late time differences in the models.The anisotropic models evolve in an exponential way towards isotropy.
It follows from Equations ( 24), (57), and (58) that in this universe model, the anisotropy parameter varies with time as   was determined using the isotropic ΛCDM universe model.It will be shown below that the upper bound on the anisotropy of the cosmic expansion is so small that the effect of the anisotropy upon the calculated value of the recombination time is negligible.Hence, we shall use the value above in the main part of this article.At the time of the recombination the anisotropy parameter was still For the present time, which is much larger than allowed by the observational restriction found by Akarsü et al. [7].We shall later see how the inclusion of mass in the universe model modifies this result.
In the early era with t t   , we can make the approximations 57) and ( 58).This gives in agreement with Equations ( 49) and (52) with 0 1   .At late times, tt   , the anisotropy decreases exponentially, and the universe model approaches the de Sitter universe.

The Anisotropic Generalization of the ΛCDM Universe Model
We now go back to the general case with cold matter, LIVE, and anisotropy.Introducing a new variable, Equation (32) takes the form Integrating this Equation with the initial condition   where t is given in equation (34), and the cosmological constant is The value of the Hubble constant, as given by the Planck project, is H0 = (67.4± 0.5) km s −1 Mpc −1 .This gives Using that the age of the universe is 0 t  13.8•10 9 years, the present value of t is 0 ˆ1.25 t  . The universe became transparent to the CMB background radiation at a point of time Equation (70) may also be written as This is the most useful expression for the time-dependence of the co-moving volume in the Bianchi type I anisotropic generalization of the ΛCDM -universe model since it separates nicely the effects of anisotropy and matter upon the time evolution of the co-moving volume.Formula (71) shows that the co-moving volume of the anisotropic generalization of the ΛCDM-model is the sum of the co-moving volume of the isotropic ΛCDM-model as given in Equation ( 33) and the anisotropic LIVE-dominated Bianchi type I universe as given in Equation (57).Expression (71) shows that the universe comes from an initial singularity with Differentiating Equation (71) and using Equation (4) gives the average Hubble parameter of the anisotropic generalization of the CDM  universe as This expression can be written as These expressions, together with the expression for t in Equation (70), show that the universe started from a Kasner-like era with Hubble parameter as given in Equation (52).The average deceleration parameter, q, of this universe model is In the early era with ˆ1 t  , this expression approaches that of the Kasner era, From Equations ( 7) and ( 71), the time-evolution of the shear scalar is This shows that for ˆ1 t  , the shear scalar decreases exponentially with time.Us- ing Equations ( 71) and (72), expression (24) for the anisotropy parameter can be written as Again, we see that the initial value of the anisotropy parameter is We proceed to calculate the directional scale factors in the same way as above for the LIVE-dominated Bianchi type I universe.We then need the integral Again, the resulting expression for i a can be simplified by introducing the con- stants in Equation (60), fulfilling the relationships (61).From Equations (20), (71), and (79), we obtain Let us here, too, consider the two cases permitted by Equation (63) with two equal scale factors.The first one is . In this case, the directional scale factors are The evolution of this model can be separated into three periods.There is a transition between the first two periods at a point of time The transition time between the two first periods is the same as in the first case, but the behavior is the opposite.The universe starts from a needle singularity; there is a transition at  through an isotropic moment to a plane-like era, and for t t   , the universe approaches isotropy exponentially.Also, in the general case with different scale factors in all directions, the evolution of this universe model can be separated into three periods with different behaviors.
1. Anisotropy dominated era.Expression (76) for the anisotropy parameter shows that the universe started from a state with maximal anisotropy,  , which means a large deceleration.Hence, the uni- verse must have started by a process not described by this solution, which has given the universe an initial expansion velocity.At the end of this singular process, which marks the beginning of the evolution described by the anisotropic ΛCDM universe model, the Hubble parameter had an infinitely large value.The universe then entered an anisotropy-dominated era with Kasner-like behavior.In this era ˆ1 t  and we can use the approximations , cosh 1 t t t . To first order in t Equations ( 71) and (73) then give in agreement with Equations (49) and (52); 2. Matter dominated era.The transition to a matter-dominated era happened at a time giving   .This happened quite early in the matter-dominated era; 3. LIVE dominated era.The transition from a matter-dominated to a LIVE-dominated era is here defined by the condition that cosmic deceleration due to the attractive gravity of the matter changes to accelerated expansion due to the repulsive gravity of the LIVE.Hence, the point of time, 3 t , of this transition is given by the condition that the deceleration parameter vanishes,   3 0 qt  .
In the case of the isotropic universe, this transition happens at a point of time, .This gives where is given in Equation (40).It is shown in Appendix A that this condition leads to the result that the change of the transition time due to the expansion anisotropy is not greater than Inserting the values of the constants gives   3 0.67 zt  like the transition redshift in the isotropic universe.
In the anisotropic universe the age of the universe is given by applying the normal- To 1. order in 0   this gives In the case of the isotropic ΛCDM-universe, this reduces to In accordance with the standard expression, (35), for the age of the isotropic ΛCDM universe model.
It follows from Equations (90) and (91) that Inserting the values for the constants as given above we obtain

Review of Some Papers on the Bianchi Type I Universe Models
In this Section I shall review some selected articles discussing Bianchi type I universe models in light of the results in the previous sections.

Anisotropic Brane Universe Models
C. M. Chen, T. Harko, and M. K. Mak [13] have discussed anisotropic brane universe models.They obtained solutions of Einstein's field equations describing several types of Universe models.I shall here focus on the model most closely related to the present work.They argued that the equation of state most appropriate to describe the high-density regime of the early Universe is the stiff Zeldovich one, with p   .In the case of ordinary general theory of relativity with no brane their constant 5 0 k  .Then their solutions of the field equations reduce to those in Section 2.4 above.

Anisotropic Universe Models with Scalar Field and Phantom Field
In 2008 B. C. Paul and D. Paul published a paper [14] where they investigated some Bianchi type I universe models with non-vanishing cosmological constant and either a scalar field or a phantom field.In the case where these fields can be neglected compared to the LIVE represented by a cosmological constant, their model reduces to the one considered in [9].
On 8 January 2024, Mark P. Herzberg and Abraham Loeb published a preprint [15] with the title Constraints on an Anisotropic Universe.They first considered a matter dominated Bianchi type I universe with vanishing cosmological constant.Then, Equation (32) reduces to Integration with   This Equation shows that there is a transition from an early period where the anisotropy dominates to a late period where matter dominates at the point of time  28) and (94) show that in the early anisotropy-dominated era, the matter density decreases as 1/t , and in the later matter-dominated era, the density decreases as 2 1/t , as it does in the isotropic universe.

Inserting Equation (94) into Equation (20) and introducing the constants
fulfilling Equation (49) leads to Herzberg and Abraham Loeb considered a model with two equal scale factors.Again, Equation (63) shows that there are two possibilities with two equal scale factors: The first one is where The other case with two equal scale factors is In this model, The age of this universe model is found either from   To 1. order in Hence, the anisotropy of the expansion reduces the age of the universe by Akarsu et al. [6,7] have found that The directional Hubble parameters of the anisotropic, mass-dominated universe model are It may be noted that 3 H can be written as Hence, the difference between 1 H and 3 H decreases as 1/t in agreement with a result found by Herzberg and Loeb [15].
In the fully asymmetric case, they write: "Now let us consider the more general case in which all the scale factors are different.In this case, we do not have an analytical solution to the above equations since the equations are all coupled.Nevertheless, we can solve the equations numerically."The reason for this is an unfortunate way of writing the field equations.As shown above, there are analytic solutions in this general case, too, for a matter-dominated universe, both with and without a cosmological constant.These were originally found by Saunders [8], although in another form than in the present paper.
Herzberg and Loeb [15] also considered a radiation-dominated universe with The solution of this Equation with   The age of this universe is given by where is the age of the corresponding radiation-dominated isotropic universe.Hence, the expansion anisotropy increases the age a little.With   which is 23 orders of magnitude greater than that per- mitted by the effect of the expansion anisotropy upon the cosmic nucleosynthesis.

A Simple Bianchi Type I Universe Applied to the Hubble Tension
In this Section, we shall first consider the most simple universe model of this type, where only one direction expands differently from the others, following M. Le Delliou et al. [16].
A model with a general perfect fluid, including pressure, was considered.However, it should be noted that the pressure of the radiation appeared only in the Einstein Equations (3)-( 5) in ref. [7], not in the integrated equations.Hence, putting p = 0, i.e., neglecting pressure, makes no changes in their results.
The line element has the form (using units so that the velocity of light where the scale factor   at is normalized so that its present value is   The authors have given the following formula for the evolution of the average Hubble parameter of this anisotropic universe model in terms of the scale factor and the anisotropy parameter and its rate of change, Hence, the Hubble constant of this anisotropic universe model is The last term is much less than 1, so we can, with sufficient accuracy, use a series expansion to first order in this term, giving Thus, the difference between the Hubble constant in this universe model with anisotropic expansion in one direction and in a universe with isotropic expansion is with good accuracy The present value of 0  coming from observations can be estimated from the for- mulae in the appendix of [16].According to their Equation (A8), , we can approximate this expression with Hence, since we are only interested in an order of magnitude estimate, we can calculate the value of the present rate of change of the anisotropy parameter from Inserting this into Equation (102) gives Using the values from ref. [16], . This is a factor 10 3.2 10   smaller than that needed to solve the Hubble tension.Hence, when the formulae of ref. [16] is applied to the Hubble tension, we find that the anisotropy of the expansion velocity, as deduced from the Planck data, is less than a hundred million times too small to solve the Hubble tension for this universe model.
I have a question from a referee: In Equation (112), only a small deviation from isotropy along one axis is considered.Will adding more general perturbations along all axes help to at least alleviate the Hubble tension?In order to answer this question I shall now present the above result in a different way.
According to Equation (112) where , It then follows from Equation (6) that for this universe model Equation ( 24) then gives the anisotropy parameter for this universe model, With the normalization   0 1 Vt  , the present value of the anisotropy parameter is It follows from Equations ( 115) and (125) that in a universe model with only a small deviation (perturbation) from isotropy in one direction, the present magnitude of the anisotropy parameter is where 0 H  is the difference between the Hubble constant in this universe model with anisotropic expansion in one direction and in a universe with isotropic expansion.
In order for the expansion anisotropy to have a possibility of solving the Hubble tension, .This result is in conflict with that of ref. [5].This may be a result of the special character of the universe model used to deduce Equation (120)-A model with anisotropy in only one direction, i.e., a model with minimal deviation from the isotropic models.Hence, in order to shed new light upon the question of whether the expansion anisotropy can solve the Hubble tension, I shall, in the next Section, consider this question from a new point of view and use the exact solution of Einstein's field equations given in Section 2 with arbitrarily large anisotropy in all three directions to calculate how much the expansion anisotropy changes the value of the Hubble constant.

The Hubble Tension Analysed by Means of the Anisotropic ΛCDM Universe Model
The value of the Hubble constant, as determined from measurements of the temperature fluctuations in the microwave background radiation by the Planck team, is H0 = (67.4± 0.5) km s −1 Mpc −1 = 1/1.45•10 10 years.With the above values of the density parameters and the anisotropy parameter, we obtain We shall now use the exact solution (71) of Einstein's field equations and the corresponding expression (73) for the Hubble parameter to investigate whether the expansion isotropy is large enough to solve the Hubble tension.We define the difference between the Hubble constant in an isotropic-and anisotropic universe as . In order to solve the Hubble tension, the present value of 0 H  must be at least as large as the difference between the late time and early time measurements of the Hubble constant, The method takes as a point of departure that the value of the Hubble parameter determined from the temperature fluctuations 380,000 years after the Big Bang, at the time of recombination, 5 ˆ3. 3 This is the relationship between the Hubble constant in the anisotropic-and the isotropic ΛCDM -universe as determined from the Planck measurements determining the value of the Hubble parameter at the point of time RC t .Hence, the difference between the Hubble constant in an isotropic-and anisotropic ΛCDM universe can be expressed as Compared with Equation (76), this may be written as It was shown above that   .This is too small to be able to solve the Hubble tension.In appendix B it is shown that if we make a series expansion of the expression (31) for the Hubble parameter to 1. order in 0   and then calculates 0 H  , the result is that 0 H  seems to be proportional to 0   .However, eq.( 129) and the expression for K  in eq.( 70), shows that Next, we consider whether the presence of radiation is of greater importance in this connection.From Equation (31), it follows that the dimensionless point of time of the cosmic recombination for an isotropic, flat universe with dust, radiation, and LIVE is The last scattering redshift is given by Akarsu et al. [6]  I will finally answer two questions that I have received from one of the referees: 1.What level of anisotropy would be necessary to solve the Hubble tension? 2. What other observations would then be contradicted?
Concerning the first question, it follows from Equation (123) that in order to have a possibility of solving the Hubble tension, the value of the anisotropy parameter at the recombination time must obey  (76)   this requires that the present value of the anisotropy parameter is Hence, so that the expansion anisotropy should be able to solve the Hubble tension, the present value of anisotropy parameter must obey .This value of the anisotropy parameter is in conflict with several measurements of the properties of the early universe.I will here refer to the analysis of Akarsu and co-workers [6,7].From observations of baryon acoustic oscillations (BAO) and the Planck measurements of CMB temperature fluctuations, they found the following restriction on the present value of the anisotropy parameter: Furthermore, they considered the restriction upon the expansion anisotropy coming from the Big Bang nucleosynthesis (BBN).It follows from Equations ( 24) and (27) that From this and Equation (3) for the redshift, we have The calculations of the amounts of the lightest element from the BBN in a universe with isotropic cosmic expansion are in agreement with observations.A large anisotropy would change the calculated quantities and would hence give calculated amounts of helium and lithium conflicting with observations.In order to avoid such a conflict, Akarsu et al. required that the anisotropy parameter at the point of time of the BBN, around three minutes after the Big Bang, should not be larger than the radiation density at the same time, . One can, with good accuracy, use the approximation . Hence, they concluded that in order not to destroy the agreement of the predictions of the BBN calculations with observations, the present anisotropy parameter is restricted to values obeying In Section 4, a universe model with only a small deviation from isotropy along one axis was considered.We can now answer the question: "Will adding more general perturbations along all axes help to at least alleviate the Hubble tension?"It was found that for a universe model with only a small deviation from isotropy along one axis, the present value of the expansion anisotropy parameter must be rather large, .Therefore, "adding more general perturbations along all axes" does indeed alleviate the Hubble tension', in fact, by 10 orders of magnitude.But even that is not sufficient to solve the Hubble tension.

Conclusions
The main topic of the present paper has been the presentation of the anisotropic Bianchi type I generalization of the isotropic ΛCDM-universe model.Considering the Equations ( 59)-(80), the evolution of this universe model can be separated into three periods.(1) Initially, the universe was in a Kasner-like era, where anisotropy dominated over matter and LIVE.The universe came from a singular state with vanishing co-moving volume inside the Hubble horizon and infinitely great Hubble parameter.The anisotropy parameter originally had its maximal value 1  .In this period, the decel- eration parameter vanished.(2) There was a transition from an early anisotropy-dominated period to a matter-dominated era at the point of time 1 t given in Equa- tion (73).The recombination time when the universe became transparent to the CMB-radiation took place early in this era.(3) Transition from a matter-dominated era with cosmic deceleration due to the attractive gravity of the matter to an era with repul-sive gravity of the LIVE took place at a point of time given in Equation (34).The predicted redshift at the transition is zTR = 0.67.
It has recently been argued that anisotropy of the universe can solve the Hubble tension.In this connection, the Bianchi type I models are of unique importance because they are the only Bianchi models that contain the ΛCDM universe model as a special case.Hence, they are the proper generalization of our standard model of the universe.
Einstein's field equations have therefore here been solved for a general Bianchi type I universe with cold matter and LIVE with a constant density, which can be represented by a cosmological constant, and the most general model of this type has been presented in a way suitable for investigations of the effects of expansion anisotropy upon observable properties of the models.In particular two models have been applied to an investigation of the Hubble tension, one with deviation from isotropic expansion in only one direction, and a general Bianchi type I model.
In this way, I have shown that the expansion anisotropy of the universe model with deviation from isotropy in only one direction is much too small to solve the Hubble tension.In the case of the general Bianchi type I, a calculation based upon the exact solution of Einstein's field equations for a general Bianchi type I universe with dust and LIVE leads to a larger effect of the expansion anisotropy upon the value of the Hubble constant than in the previous cases.However, even if radiation is included in the universe model when the point of time of the recombination is calculated, and the constraints come from the requirement that the cosmic expansion anisotropy shall not have a significant effect upon the cosmic nucleosynthesis during the first ten minutes of the history of the universe is neglected, the permitted anisotropy is still not large enough to solve the Hubble tension.
Akarsu and co-workers have shown that if one takes into account the constraint that the expansion anisotropy shall have no significant effect on the standard Big Bang nucleosynthesis, the largest permitted value of the anisotropy parameter is Ωσ0 = 10 −23 [6], which is a hundred thousand times smaller than if these constraints are neglected.Hence, if this last requirement is not neglected, the permitted anisotropy is much too small to solve the Hubble tension.

Appendix B. Significance of a Universe Model, Which Is an Exact Solution of Einstein's Field Equations
In order to illustrate the point of the heading of this Appendix, it is useful to consider the anisotropic extension of the ΛCDM universe, i.e., a flat universe with LIVE and dust, neglecting radiation, since the field equations can be solved exactly in terms of elementary functions for such a model.In this case, Equation (31) takes the form Hence, Putting     in Equations (A9) and (A10) we obtain,

where 0 
is the present value of the shear scalar.This relationship is valid inde- pendently of the energy and matter contents of the universe.It is a kinematical relationship in the Bianchi type I universe models.Saunders first considered a Bianchi type I universe filled with cold matter in the form of dust with density M  and vanishing pressure, radiation with density R  , and dark energy in the form of LIVE having a constant density,   , which can be repre- sented by the cosmological constant  .
cause attractive gravity.On the other hand, LIVE has negative pressure p strong repulsive gravity corresponding to twice the gravity of the mass density.


are the present values of the average density of LIVE, cold matter, radiation, and Zeldovich fluid.
It has become usual to express Equation (25) in terms of the present values of the density parameters of LIVE, cold matter, radiation, the stiff fluid with

H
is the present value of the av- erage Hubble parameter, i.e., the Hubble constant of the anisotropic universe model.From Equation (23) and the expression (24) for giving follows from Equation (42) that

BL
two equal scale factors in this very early LIVE and anisotropy-dominated era.The first one has scale factors that of the empty Kasner universe.The time of recombination,

3 M
70) shows that for ˆ1 t  , this universe goes into an era with eternal expo- nential expansion.The present values of the mass parameters of matter and dark energy are  .From the baryonic acoustic oscillations and cosmic micro- wave background (CMB) data, Akarsu and co-workers[6,7] obtained the constraint Ωσ0 < 10 −18 .Furthermore, they wrote: "Demanding that the expansion anisotropy has no significant effect on the standard big bang nucleosynthesis (BBN), we find the constraint Ωσ0  10 −23 ."I shall here use Ωσ0 ≈ 10 −18 .
to the constant anisotropy of an empty Kasner universe with vanishing cosmological constant, and ends in an isotropic state with 0  .The initial value of the decelera- tion parameter was   02 q

2 t
when the two terms of the expression (71) for the co-moving volume had the same size, i.e., values of the co-moving volume, Hubble parameter, and the anisotropy parameter at this point of time was universe became transparent for the CMB background radiation, as given by Equation (76), was the transition time is close to that of the corresponding isotropic universe given in Equation (46).The corresponding redshift is

and 2 aaaa
M H is the Hubble constant with two equal scale factors.Although the initial co-moving volume vanishes, the scale factors 1 a are non-vanishing at the initial moment and then form a plane while 3 (0) 0 a  .Hence, this universe model starts from an initial plane-like singularity.There is no transition to a needle-like period with  for all values of t .

1 M
Hence, there is an initial singularity with a needle-like character.Similarly, as in the first case, there is no transition to a plane-like period since this mass-dominated Bianchi type I universe model reduces to the Einstein de Sitter-universe with mass parameter of the cold matter 0   , and the formulae for the co-moving volume and the scale factor reduce to finds the age of the corresponding anisotropic universe model by in- anisotropy decreases the age by less than 28 years.The average Hubble parameter of the matter-dominated universe model with two equal scale factors is found by differentiation of Equation (94): few minutes.It should be noted, however, that this model cannot give a realistic description of the universe.It requires that


departure from isotropy is measured by the anisotropic perturbation parameter  .In accordance with observational constraints, Delliou et al. applied the initial condition of time of the cosmic recombination, .It is shown below that  is a decreasing func- tion of time, so its present value is

0 H
must be at least as large as the difference between the early type and late type determinations of the Hubble constant,
difference between the late time and early time determinations of the Hubble constant, are in conflict with these measurements.
needed to solve the Hubble tension would lead to a severe conflict between BBN calculations in such an anisotropic universe and the observed cosmic amounts of the lightest element.
for the anisotropy shall be able to solve the Hubble tension.In the context of a universe model with anisotropic expansion in all three directions it is sufficient that

)
Comment [M1]: Please add citation in the main text for Appendix B.
From then on the anisotropy decreases at an increasing rate, and the value of the anisotropy parameter at the point of time of the recombination, 10 to the square root of /s per Mpc.It is still too small to have any potential for solving the Hubble tension. km Since the effect of the expansion anisotropy is small, we can, with good accuracy, use a series expansion of the deceleration parameter about the point of time