#### Abstract

We determine the energy-momentum tensor of nonperfect fluids in thermodynamic equilibrium and, respectively, near to it. To this end, we derive the constitutive equations for energy density and isotropic and anisotropic pressure as well as for heat-flux from the corresponding propagation equations and by drawing on Einstein’s equations. Following Obukhov on this, we assume the corresponding space-times to be conform-stationary and homogeneous. This procedure provides these quantities in closed form, that is, in terms of the structure constants of the three-dimensional isometry group of homogeneity and, respectively, in terms of the kinematical quantities expansion, rotation, and acceleration. In particular, we find a generalized form of the Friedmann equations. As special cases we recover Friedmann and Gödel models as well as nontilted Bianchi solutions with anisotropic pressure. All of our results are derived without assuming any equations of state or other specific thermodynamic conditions a priori. For the considered models, results in literature are generalized to rotating fluids with dissipative fluxes.

#### 1. Introduction

In this paper, we consider systems described by Einstein’s equations:with an energy-momentum tensor and equations of state, neither of which are specified by any* ad hoc* assumptions. Instead, we discuss the whole question from a thermodynamic perspective (we emphasize that we approach this without any specific thermodynamic conditions as done, e.g., in [1]; i.e., we refrain from applying linear or extended thermodynamics). This consideration is discussed for a class of cosmologically interesting metrics (introduced and shown to be observationally admissible in [2–5], see (5) below). In terms of the temperature and the kinematic invariants characterizing the matter, our consideration provides a class of general equations of state (“matter equations”) which are compatible with Einstein’s equations and correspond to generalized Friedmann equations. This framework can find (and has found) applications in relativistic cosmology and astrophysics. Basically, it allows to go beyond the standard phase cosmology (governed by phases with certain equations of state like inflation, radiation, and dust that have to be fitted by fine-tuning) and to describe the cosmological state transitions from phase to phase by intermediate stages. However, in [6–9] a fine-tuned sudden passage from the decelerated to the accelerated regime, as observed today, produces inadequacies. These are then avoided by* ad hoc* introduced equations of state where viscosity originates from geometry (e.g., , ). Our calculations can provide a theoretical foundation of such equations. Furthermore, this framework also contributes to a physical discussion of no-go theorems like the shear-free fluid conjecture [10]. For instance, this thermodynamic approach enables one to sharpen the theorem (proved in [11], without explicitly referring to thermodynamics) which states for nonvanishing acceleration that rotation and expansion cannot simultaneously be equal to zero: in [11] it has been shown that models with vanishing acceleration do not allow for nonvanishing rotation.

Ehlers et al. have proven [12] that the high isotropy of the cosmic microwave background (CMB) and the vanishing of shear of a congruence of curves (for the definition of the quantities, shear, rotation, acceleration, and expansion, see [1] or see also Section 2) are closely related to the requirement that there exists a conformal Killing vector field (in [13], it is shown that while for fluids the CKV property is essential, for fluids this property has to be generalized to conformal collineation) being parallel to the tangents of the curves (i.e., to the velocity of a streaming fluid). In detail, it was shown (Lemma 3 in [12]) that a space-time admits a time-like conformal Killing vector field (CKV)with , if and only if there is a velocity field with satisfyingwhere is the acceleration and the expansion. Oliver Jr. and Davis [14] showed that (3) is a necessary and sufficient condition for the existence of a CKV in the case of rotating space-times, too. In the following, we consider such conform-stationary space-times (according to [15, 16] this is equivalent to parallax-free cosmological models).

In particular, the second condition allows implying a parameter which can be identified as the inverse temperature, , so that can be interpreted as temperature vector. This parameter occurs if the second equation in (3) is rewritten as (Theorem 2.1 in [14])

Additionally, we assume that the considered space-times are spatially homogeneous. This reduction to Bianchi-type models still allows for the matter distribution to be anisotropic, while the CMB is isotropic.

Altogether we are led to the subclass of tilted Bianchi models constructed by Obukhov [2–5] that admit a CKV,Thereby, rotating and expanding models with acceleration and isotropic CMB are considered in [2–5]. In contrast to that, in [17–19], Ehlers-Geren-Sachs theorems were, partly in a generalized version, used to study and determine a class of space-times containing also inhomogeneous cosmological models, with nontrivial acceleration but zero rotation.

To complete notation used in (5), we definewith and as arbitrary constants. Here the triad components form a basis which is invariant under the spatial isometries of the Bianchi models. Accordingly, their Lie derivative with respect to the generating Killing vector fields (KV) vanishes (for details see [20]). The components are supposed to be functions of the space-like canonic coordinates only and determine the metric (5) as to the Bianchi-type. The coordinate denotes the proper time with respect to a fluid particle and is the scale factor.

Furthermore, Latin indices are used for the coordinate while Greek indices are used for the tetrad description of the tensor components. All indices with hat denote the three spatial dimensions, for example, and , whereas those without hat run through all four space-time dimensions, and .

Regarding the thermodynamic proposition, we follow the Eckart approach and assume that the model under consideration is in thermodynamic (near)equilibrium [21]. The Eckart approach to the relativistic Theory of Irreversible Processes [22] (see also [23–25]) is based on the balance equations for the particle number (where represents mass density and the specific volume)the energy-momentum (in particular, regarding (7), the null-component can be interpreted as the first law of thermodynamics [26, 27], i.e., as the conservation of internal energy), and the entropywhere denotes the entropy vector and the density of the nonnegative entropy production. In the case of thermodynamic equilibrium (vanishing of entropy production), appropriated supplementary conditions have to be added by hand. In the general-relativistic version of this theory the framework is completed by Einstein’s gravitational equations (1). However, in the general-relativistic Theory of Irreversible Processes no further assumptions have to be introduced* ad hoc* in order to yield thermodynamic equilibrium [21].

Now, if the entropy vector is defined according to [23, 24, 28],the entropy production can be reformulated asHere denotes the heat-flux, the energy density, the pressure, and the temperature vector and . Finally, by decomposing the energy-momentum tensor (27) this yields [28]As shown in [21], regarding the conformal Killing equation (2), the second term in brackets turns out to be traceless which results in a vanishing entropy production This shows that nonperfect fluids are not necessarily incompatible with reversible thermodynamics [28–30]. However, for space-times without a CKV or KV (13) can only be solved by assuming a perfect fluid. That this CKV property is not purely mathematical but has also a physical meaning is supported by the following arguments.

Firstly, the derivations of (12) and (13), given in [28], show that the quantities , , and are thermodynamically well determined. Indeed, it is assumed that the specific entropy is given as a function of the specific internal energy and the specific volume ; that is,so that (Gibbs equation) Moreover, for a comoving observer the relation holds. For the thermodynamic quantities defined in this way (13) is valid. Secondly, (13) has the following solutions: Either the fluid is perfect or the temperature vector has to be a CKV (containing the special case of a Killing vector field). Therefore, the CKV property is justified by defining equilibrium or near-equilibrium states in the framework of reversible thermodynamics [21].

This is confirmed by the fact that , being a CKV, leads to some well-known models like Friedmann’s and Gödel’s space-times with the corresponding equations of state (see Section 3.2). Furthermore, it should be emphasized that the justification of the thermodynamic meaning of the CKV condition given via (13), that is, in the context of phenomenological continuum theory, is supported by considerations in the framework of kinetic theory, where the CKV property of in combination with related equations of state for some special cases is derived from Boltzmann’s equation [31–34].

Based on the existence of such a CKV one can derive a set of four propagation equations for nonperfect fluids (see [11]), which link the propagation of the matter content to the kinematic description of the space-time (see Section 3).

The paper is organized as follows: In Section 2, we introduce a suitable tetrad frame that allows us to establish manageable equations. In addition, the decomposition of the energy-momentum tensor with respect to kinematic invariants is shortly reviewed and their form in tetrads for the space-times (5) is derived. Subsequently, by solving the propagation equations, we deduce in Section 3 general expressions for the whole matter content depending on the structure constants and the kinematic invariants, respectively. After checking the consistency with Einstein’s equations the general case of a nonperfect fluid and particular cases like nontilted [35] and stationary models are discussed. Among the special cases that can be recovered are the Friedmann and the Gödel models. In Section 4 we discuss the results and provide alternative formulations, relevant for further observational and thermodynamic considerations.

#### 2. Tetrad Formulation and Kinematic Invariants

In the following, we introduce tetrads (see, e.g., [36]) that allow for a convenient separation of the variable objects, and , and the constants, and in (5). Definingwiththe tetrads can be chosen asTo fulfill the relationsthe constant and symmetric matrix has to take the formThe structure constants of the isometry groups acting on the space-like hypersurfaces and specific to the Bianchi models can be expressed by a 4-dimensional representation:such that , , and . Expressions for the curvature tensors and scalars in terms of these newly introduced tetrads are derived in Appendix B.

On the basis of these preliminaries, we now introduce the kinematic invariants and the respective decomposition of the energy-momentum tensor.

Assuming a one-component fluid with the four-velocity , such that , the gradient of can be decomposed kinematically [1]:Thus, rotation, shear, acceleration, and expansion as well as the scalar quantities of rotation and acceleration readChoosing , these kinematic quantities are rewritten in the tetrad representation and for the space-times (5) as follows (the subscript denotes the covariant tetrad derivative)where .

According to [1], the energy-momentum tensor can be decomposed with respect to the timelike velocity field :Here the quantities can be identified with the appropriate projections, for the energy density, for the isotropic pressure, for the heat-flux, and for the anisotropic pressure.

#### 3. Matter Equations

The conditions for the temperature vector , being a CKV, areBoth can be found in [14] (the latter reproduces the -component of (4)). The two results (28) and (29) are obtained by inserting into (2) and multiplying this equation by and , respectively.

Furthermore, integration of (29) leads to an expression for the temperature scalar,and the conformal factor,where is the constant of integration.

The existence of the CKV has far-reaching consequences for the geometry of the space-time and, factoring in Einstein’s field equations, for the matter. By drawing on the Ricci identity for the CKV and the Bianchi identity subsequently, we deduce a set of four propagation equations [11, 21]. The first two describe the evolution of the energy density and the isotropic pressure :The other two equations describe the change of the heat-flux ,and the anisotropic pressure ,

##### 3.1. Solutions of the Propagation Equations

The reformulation of the dynamic equations (32)–(34) in terms of the space-times (5) and some tedious algebra brings us to a set of ordinary differential equations which can be solved analytically.

To this end, we decouple (32) toBy means of (18), (35) can be rewritten asIntegrating (36) yieldswhere the objects and represent the summed constants of integration of the additive antiderivatives; see Appendix C. For concrete cosmological or astrophysical models, for example, stars, the boundaries of the respective integrals are specified.

With the identityand (24), the tetrad formulations for the propagation equations of the heat-flux (33) and the anisotropic pressure (34) together with the kinematic quantities (24)–(26) yield the following integrable partial differential equations: for the heat-flux and the anisotropic pressure, respectively.

Integration and reorganization of terms bring the wanted solutions where, similarly to the case of the energy density (37) and the isotropic pressure (38) above, the objects and represent the constants of integration (see Appendix C).

With the help of the kinematic quantities (24)–(26) and (A.1)–(A.4) of Appendix A the solutions (37), (38), (41), and (42) can be rewritten as follows: or, in terms of purely kinematic quantities, The expressions (43) and (44) can be understood as generalized Friedmann equations.

According to (27), we can now reconstruct the energy-momentum tensor by inserting the four solutions above: or whereIn order to verify the consistency of solutions (43)–(46) with Einstein’s field equations and in order to recover special cases, the constants of integration are determined by the calculations of Appendix C. The matter equations (43)–(46) then take the exclusively kinematic forms:in which for the latter, , relation (A.3) was used in addition. By multiplying these expressions with the tetrads, one obtains the coordinate representation without any additional terms.

Notice that (51) and (52) satisfy the Raychaudhuri equation:

##### 3.2. Special Cases

###### 3.2.1. Nontilted Models

The nontilted limit () leads to purely expanding models, that is, those with vanishing rotation and acceleration. In this case also the coefficients of (5) become zero, so that the space-times are conformally static [37, 38].

Following up on this premise, (51) and (52) turn into the equationsThese correspond to the Friedmann equationsif the curvature parameter and the Ricci scalar of the 3-dimensional Bianchi spaces are related by This result is in accordance with [39, p. 474]. The two constants of integration are then related by

Furthermore, the heat-flux (53) is identically zero, while for the anisotropic pressure (54) one getsAccording to Section 4 of [37], the nonperfect fluid models investigated here can be subdivided into three further classes. In detail, this amounts to determining the number of distinct eigenvalues of the anisotropic pressure (60): If has three different eigenvalues, the space-time is of Petrov-type I. In the case of two different eigenvalues, one obtains Petrov-type D. Finally, if there is only one eigenvalue, it can only be zero and results in vanishing identically. Therefore, only this latter case is in general a sufficient condition for obtaining perfect fluid Friedmann models.

An example for a nontilted space-time subclass of (5) which does not contain Friedmann models is provided by the Bianchi-type IV metric with canonic coordinates (, , , ) and undetermined constants introduced in (5). Since here one finds three distinct eigenvalues for (60), models (61) are of Petrov-type I.

For the particular case , metric (61) takes the formwhich coincides exactly with the example discussed in [37] (see (4.9) therein).

###### 3.2.2. Stationary Models

We check the consistency of our results with the stationary limit as to [3, 4]. Accordingly, for vanishing expansion and generally nontrivial rotation, (51) and (52) reduce toso thatIn the perfect fluid limit with vanishing heat-flux and anisotropic pressure, one obtains from (53) the condition and from (54) As a more concrete* ansatz* we choose a Bianchi-type III subclass of space-times (5): where and and with , , and being constant. Admitting in general nonvanishing rotation and expansion, this metric is also denoted as the Gödel-type model (see [3, 4]). By this choice, the heat-flux (65) vanishes identically, while the anisotropic pressure condition (66) holds only for at least either of the two relations:Furthermore, one has or, by (68), respectively, This yields and thus, for expression (64), According to, for example, [28], this is just the equation of state of the classical Gödel space-time. Indeed, in [3] it is stated that yields closed timelike curves.

#### 4. Discussion

In this paper, we considered homogeneous and conform-stationary space-times (5) with Bianchi group invariance and an arbitrary matter source, which allows for generally tilted models.

By solving the propagation equations (32)–(34), we deduce explicit expressions for the energy density (37), the isotropic pressure (38), the heat-flux (41), and the anisotropic pressure (42) in terms of the scale factor, the tetrad components (18), and the structure constants. These results are rewritten in terms of the kinematic quantities, as to be found in (43), (44), (45), and (47), and are combined to the energy-momentum tensor, (48) or (49). Similar equations are* ad hoc* assumed in [6–9] in order to solve problems arising during cosmological evolution for different reasons.

In addition to the Raychaudhuri equation and the other propagation and constraint equations (see, e.g., [1, 40]), we obtain equations in which the expressions for the matter content are decoupled and independent of higher derivatives of the kinematic quantities (except for the expansion and acceleration) or depending on the electric part of the Weyl tensor. Particularly, no equations of state or further thermodynamic relations have to be assumed to arrive at these results. Here it should be emphasized that the vanishing shear does not necessarily imply a zero anisotropic pressure as required by linear thermodynamics. It should also be pointed out that more-component fluids or a cosmological constant can easily be included.

Equations (37), (38), (41), and (42) represent a class of models which does not only contain physically relevant space-times. To take into account well-motivated (energy) conditions or global aspects (as considered in [41]) which should provide broader restrictions is therefore a subject of future research.

Moreover, inspection of (43), (44), (45), and (47) underlines that further thermodynamic assumptions like an equation of state, Fourier’s law, Cauchy’s law, or expressions from extended thermodynamics will further restrict possible solutions. This becomes manifest, if one rewrites (37), (38), (41), and (42) with the help of (29): These equations describe the temperature dependence of the matter content which has to be fulfilled for the considered class of models.

Expressions (72) and (73) can be used to construct equations of state. For instance, one can combine the two in such a way that the outcome does not contain the structure constants:which is a possible equation of state for the considered space-time class. This relation clearly shows that the pressure has a difficult dependence on the temperature and its first and second derivatives. Of course, has an explicit temperature dependence as given in (72), but assuming the validity of simple equations of state, like , an effective fine-tuning has to be done in order to prevent an additional temperature dependence of .

It becomes obvious from (74) and (75) that the assumption of Fourier’s or Cauchy’s laws consequently generates additional strong restrictions on the space-time and its matter content. The same is true for other* ad hoc* introduced constitutive equations. This includes nonlinear ones like the heat-flux law of Israel-Stewart-type [23] which is physically motivated by the fact that it overcomes stability and causality problems arising in the linear case. Our point, however, is to ask for those constitutive equations and equations of state, respectively, which follow from the conservation laws in a prescribed geometry and a given temperature field. Thus, it is not in the sense of the present consideration, to additionally impose* ad hoc* constitutive equations on (74) and (75). In nonrelativistic continuum thermodynamics the situation is different. There, one has to complete the system of basic equations resting on the conservation or balance equations by adding such* ad hoc* relations manually. If (74) and (75) differ from those* ansatzes* made by hand, this can have a variety of reasons and implications. To call only one, if there were severe thermodynamic arguments for one of the linear or nonlinear* ad hoc ansatzes*, for example, for the heat-flux, one was obliged to ask under which condition it is compatible with (74).

This view is reinforced by the results obtained in [42]. There, it was shown that in conform-stationary models the heat-flux must vanish for zero anisotropic pressure and under the assumption of a heat-flux law of the Israel-Stewart type. An example for physical processes in which this does not hold (Landau damping) is also provided in [42].

Moreover, the form of expressions (74) and (75), which denotes the modified laws of Fourier and Cauchy, is pointing in a direction that is to be found in various formulations of extended thermodynamics [23, 43]. This becomes evident if one rewrites (75) with the help of (72) and (74):As a constitutive quantity is a function which is linear and quadratic in the heat-flow and linear in the energy density, while the temperature is also included with its first and second time derivative.

The consideration of simple models like nontilted or stationary ones leads back to, for example, the well-known Friedmann or the Gödel space-times (in both cases the constants of integration are determined; see Section 3) and similar anisotropic models as discussed in [37]. In this context expressions (43) and (44) or (51) and (52), respectively, can be understood as generalized Friedmann equations.

By rewriting (43), (44), (45), and (47) in terms of the observational quantities for the Hubble function and for the deceleration parameter one receives limits on acceleration, rotation, heat-flux, and anisotropic pressure.

The corresponding equations take the formso that the matter content can be described by the observable quantities and and the model-dependent constants , , and as well as the constants of integration, eventually given by initial or boundary conditions.

In analogy to the calculations which lead to (76), one obtains from (78) and (79)which can again be regarded as an equation of state given by observational quantities. The class of models we consider here may have an anisotropic behavior of the Hubble flow and the galaxy distribution function [3, 4]. In this context, the observation of a large-scale flow of galaxies, called “dark flow,” with respect to the CMB is remarkable (see [44] for a review). A detailed discussion of this and other possibly observable effects in nonrotating models can be found in [17–19].

Refraining from possible further restrictions on relations (78)–(81), one finds the following hypothetical scenario.

For a large scale factor , the structure constants and the constants of integration are negligible, such that, for the behavior of the matter content, the expansion rate and the deceleration rate are most important. Moreover, one sees that has critical values at which the behavior of the matter variables changes. For instance, in the case of large accelerations (), which for cosmological models means a strongly increasing expansion and for local models (like stars) a strongly increasing collapse, most matter variables change the sign. All matter variables display generally the same dependence on the expansion rate and are therefore of likewise importance. A more detailed discussion can only be achieved if the dependence on and is fixed for specified Bianchi models.

For small values of the scale factor , that is, in the early cosmological phase or for objects which become very dense, the structure and the integration constants become much more important in comparison to and . Besides, all matter variables show the same behavior and are therefore of equal importance. When the scale factor increases, the heat-flux and the anisotropic pressure essentially behave like the energy and the pressure; they dilute.

As a final remark, we would like to stress out that our results are obtained using the phenomenological approach to relativistic thermodynamics. As indicated in the introduction, in this context the vanishing entropy production (13) is a necessary condition for equilibrium, whereas it is a necessary and sufficient condition for the description of reversible processes. To consider the compatibility of these results with relativistic kinetic theory would be an interesting task for further investigations.

#### Appendix

#### A. Kinematic Relations

The following relations between the kinematic invariants are used to obtain and simplify results of the Sections 3 and 4 in kinematic terms.

#### B. Tetrad Formulation of Curvature

The connection coefficients in the tetrad formulation (see Section 2), the so-called Ricci rotation coefficients, can be expressed in terms of the Christoffel symbols ,or, due to (18) by the structure constants (21), respectively, Determining the Riemannian curvature tensor by the Ricci-identity and the tetrads, the Ricci tensor can be brought to the formwithAccordingly, the Ricci scalar becomeswithThe expressions of (B.5) and equivalently of (B.7) can be identified with the Ricci tensor and the Ricci scalar of 3-dimensional Bianchi spaces [45].

This results in the following shape of the Einstein tensor:

#### C. Constants of Integration

From the field equations, , in terms of the tetrad formulation from Section 2 and together with (5), one finds Then, because of (50) the constants of integration, that is, , , , and , for the energy density, the isotropic pressure, the heat-flux, and the anisotropic pressure (43)–(46), become, in this order,where in (C.4) and (C.5) it was made use of relation (A.4). If one reinserts the constants of integration (C.2)–(C.5) into the matter equations (43)–(46), they take the purely kinematic forms (51)–(54).

The* summarized* constants of integration, , , , and , in this paper are pieced together as follows:The occurring objects , , , , , , and are the* actual* constants of integration yielded by the following integrals, which are to be calculated in Section 3.1:

#### Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.