International Journal of Mathematics and Mathematical Sciences
Volume 2009 (2009), Article ID 471938, 20 pages
doi:10.1155/2009/471938
Research Article

Accelerated Lyra's Cosmology Driven by Electromagnetic Field in Inhomogeneous Universe

1Department of Mathematics, Hindu P. G. College, Zamania, Ghazipur, Uttar Pardesh 232 331, India
2Department of Mathematics, P. G. College, Ghazipur 233 001, India

Received 13 July 2009; Accepted 11 December 2009

Academic Editor: Christian Corda

Copyright © 2009 Anirudh Pradhan and Padmini Yadav. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

A new class of cylindrically symmetric inhomogeneous cosmological models for perfect fluid distribution with electromagnetic field is obtained in the context of Lyra's geometry. We have obtained solutions by considering the time dependent displacement field. The source of the magnetic field is due to an electric current produced along the z-axis. Only 𝐹 1 2 is a nonvanishing component of electromagnetic field tensor. To get the deterministic solution, it has been assumed that the expansion 𝜃 in the model is proportional to the shear 𝜎 . It has been found that the solutions are consistent with the recent observations of type Ia supernovae, and the displacement vector 𝛽 ( 𝑡 ) affects entropy. Physical and geometric aspects of the models are also discussed in presence and absence of magnetic field.

1. Introduction and Motivations

The inhomogeneous cosmological models play a significant role in understanding some essential features of the universe such as the formation of galaxies during the early stages of evolution and process of homogenization. The early attempts at the construction of such models have been done by Tolman [1] and Bondi [2] who considered spherically symmetric models. Inhomogeneous plane-symmetric models were considered by Taub [3, 4] and later by Tomimura [5], Szekeres [6], Collins and Szafron [7, 8], and Szafron and Collins [9]. Senovilla [10] obtained a new class of exact solutions of Einstein's equations without big bang singularity, representing a cylindrically symmetric, inhomogeneous cosmological model filled with perfect fluid which is smooth and regular everywhere satisfying energy and causality conditions. Later, Ruiz and Senovilla [11] have examined a fairly large class of singularity-free models through a comprehensive study of general cylindrically symmetric metric with separable function of 𝑟 and 𝑡 as metric coefficients. Dadhich et al. [12] have established a link between the FRW model and the singularity-free family by deducing the latter through a natural and simple in-homogenization and anisotropization of the former. Recently, Patel et al. [13] have presented a general class of inhomogeneous cosmological models filled with nonthermalized perfect fluid assuming that the background space-time admits two space-like commuting Killing vectors and has separable metric coefficients. Singh et al. [14] obtained inhomogeneous cosmological models of perfect fluid distribution with electromagnetic field. Recently, Pradhan et al. [1518] have investigated cylindrically-symmetric inhomogeneous cosmological models in various contexts.

The occurrence of magnetic field on galactic scale is a well-established fact today, and its importance for a variety of astrophysical phenomena is generally acknowledged as pointed out by Zeldovich et al. [19]. Also Harrison [20] suggests that magnetic field could have a cosmological origin. As natural consequences, we should include magnetic fields in the energy-momentum tensor of the early universe. The choice of anisotropic cosmological models in Einstein system of field equations leads to the cosmological models more general than Robertson-Walker model [21]. The presence of primordial magnetic field in the early stages of the evolution of the universe is discussed by many [2231]. Strong magnetic field can be created due to adiabatic compression in clusters of galaxies. Large-scale magnetic field gives rise to anisotropies in the universe. The anisotropic pressure created by the magnetic fields dominates the evolution of the shear anisotropy and decays slowly as compared to the case when the pressure is held isotropic [32, 33]. Such fields can be generated at the end of an inflationary epoch [3438]. Anisotropic magnetic field models have significant contribution in the evolution of galaxies and stellar objects. Bali and Ali [39] obtained a magnetized cylindrically symmetric universe with an electrically neutral perfect fluid as the source of matter. Pradhan et al. [4044] have investigated magnetized cosmological models in various contexts.

In 1917 Einstein introduced the cosmological constant into his field equations of general relativity in order to obtain a static cosmological model since, as is well known, without the cosmological term his field equations admit only nonstatic solutions. After the discovery of the red-shift of galaxies and explanation thereof Einstein regretted for the introduction of the cosmological constant. Recently, there has been much interest in the cosmological term in context of quantum field theories, quantum gravity, super-gravity theories, Kaluza-Klein theories and the inflationary-universe scenario. Shortly after Einstein's general theory of relativity Weyl [45] suggested the first so-called unified field theory based on a generalization of Riemannian geometry. With its backdrop, it would seem more appropriate to call Weyl's theory a geometrized theory of gravitation and electromagnetism (just as the general theory was a geometrized theory of gravitation only), instead of a unified field theory. It is not clear as to what extent the two fields have been unified, even though they acquire (different) geometrical significance in the same geometry. The theory was never taken seriously inasmuchas it was based on the concept of nonintegrability of length transfer; and, as pointed out by Einstein, this implies that spectral frequencies of atoms depend on their past histories and therefore have no absolute significance. Nevertheless, Weyl's geometry provides an interesting example of nonRiemannian connections, and recently Folland [46] has given a global formulation of Weyl manifolds clarifying considerably many of Weyl's basic ideas thereby.

In 1951 Lyra [47] proposed a modification of Riemannian geometry by introducing a gauge function into the structureless manifold, as a result of which the cosmological constant arises naturally from the geometry. This bears a remarkable resemblance to Weyl's geometry. But in Lyra's geometry, unlike that of Weyl, the connection is metric preserving as in Remannian; in other words, length transfers are integrable. Lyra also introduced the notion of a gauge and in the “normal” gauge the curvature scalar in identical to that of Weyl. In consecutive investigations Sen [48], Sen and Dunn [49] proposed a new scalar-tensor theory of gravitation and constructed an analog of the Einstein field equations based on Lyra's geometry. It is, thus, possible [48] to construct a geometrized theory of gravitation and electromagnetism much along the lines of Weyl's “unified” field theory, however, without the inconvenience of nonintegrability length transfer.

Halford [50] has pointed out that the constant vector displacement field 𝜙 𝑖 in Lyra's geometry plays the role of cosmological constant Λ in the normal general relativistic treatment. It is shown by Halford [51] that the scalar-tensor treatment based on Lyra's geometry predicts the same effects within observational limits as the Einstein's theory. Several authors Sen and Vanstone [52], Bhamra [53], Karade and Borikar [54], Kalyanshetti and Waghmode [55], Reddy and Innaiah [56], Beesham [57], Reddy and Venkateswarlu [58], Soleng [59], have studied cosmological models based on Lyra's manifold with a constant displacement field vector. However, this restriction of the displacement field to be constant is merely one for convenience, and there is no a priori reason for it. Beesham [60] considered FRW models with time-dependent displacement field. He has shown that by assuming the energy density of the universe to be equal to its critical value, the models have the 𝑘 = 1 geometry. T. Singh and G. P. Singh [6164], Singh and Desikan [65] have studied Bianchi-type I, III, Kantowaski-Sachs and a new class of cosmological models with time-dependent displacement field and have made a comparative study of Robertson-Walker models with constant deceleration parameter in Einstein's theory with cosmological term and in the cosmological theory based on Lyra's geometry. Soleng [59] has pointed out that the cosmologies based on Lyra's manifold with constant gauge vector 𝜙 will either include a creation field and be equal to Hoyle's creation field cosmology [6668] or contain a special vacuum field, which together with the gauge vector term, may be considered as a cosmological term. In the latter case the solutions are equal to the general relativistic cosmologies with a cosmological term.

Recently, Pradhan et al. [6973], Casana et al. [74], Rahaman et al. [75, 76], Bali and Chandnani [77], Kumar and Singh [78], Singh [79], Rao, Vinutha et al. [80] and Pradhan [81] have studied cosmological models based on Lyra's geometry in various contexts. Rahaman et al. [82, 83] have evaluated solutions for plane-symmetric thick domain wall in Lyra geometry by using the separable form for the metric coefficients. Rahaman [8487] has also studied some topological defects within the framework of Lyra geometry. With these motivations, in this paper, we have obtained exact solutions of Einstein's modified field equations in cylindrically symmetric inhomogeneous space-time within the frame work of Lyra's geometry in the presence and absence of magnetic field for time varying displacement vector. This paper is organized as follows. In Section 1 the motivation for the present work is discussed. The metric and the field equations are presented in Section 2. In Section 3 the solutions of field equations are derived for time varying displacement field 𝛽 ( 𝑡 ) in presence of magnetic field. Section 4 contains the physical and geometric properties of the model in presence of magnetic field. The solutions in absence of magnetic field are given in Section 5. The physical and geometric properties of the model in absence of magnetic field are discussed in Section 6. Finally, in Section 7 discussion and concluding remarks are given.

2. The Metric and Field Equations

We consider the cylindrically symmetric metric in the form

𝑑 𝑠 2 = 𝐴 2 𝑑 𝑥 2 𝑑 𝑡 2 + 𝐵 2 𝑑 𝑦 2 + 𝐶 2 𝑑 𝑧 2 , ( 2 . 1 ) where 𝐴 is the function of 𝑡 , alone and 𝐵 and 𝐶 are functions of 𝑥 and 𝑡 . The energy momentum tensor is taken as having the form

𝑇 𝑗 𝑖 = ( 𝜌 + 𝑝 ) 𝑢 𝑖 𝑢 𝑗 + 𝑝 𝑔 𝑗 𝑖 + 𝐸 𝑗 𝑖 , ( 2 . 2 ) where 𝜌 and 𝑝 are, respectively, the energy density and pressure of the cosmic fluid, and 𝑢 𝑖 is the fluid four-velocity vector satisfying the condition

𝑢 𝑖 𝑢 𝑖 = 1 , 𝑢 𝑖 𝑥 𝑖 = 0 . ( 2 . 3 ) In (2.2), 𝐸 𝑗 𝑖 is the electromagnetic field given by Lichnerowicz [88]

𝐸 𝑗 𝑖 = 𝜇 𝑙 𝑙 𝑢 𝑖 𝑢 𝑗 + 1 2 𝑔 𝑗 𝑖 𝑖 𝑗 , ( 2 . 4 ) where 𝜇 is the magnetic permeability and 𝑖 the magnetic flux vector defined by

𝑖 = 1 𝜇 𝐹 𝑗 𝑖 𝑢 𝑗 , ( 2 . 5 ) where the dual electromagnetic field tensor 𝐹 𝑖 𝑗 is defined by Synge [89]

𝐹 𝑖 𝑗 = 𝑔 2 𝜖 𝑖 𝑗 𝑘 𝑙 𝐹 𝑘 𝑙 . ( 2 . 6 ) Here 𝐹 𝑖 𝑗 is the electromagnetic field tensor and 𝜖 𝑖 𝑗 𝑘 𝑙 is the Levi-Civita tensor density.

The coordinates are considered to be comoving so that 𝑢 1 = 0 = 𝑢 2 = 𝑢 3 and 𝑢 4 = 1 / 𝐴 . If we consider that the current flows along the 𝑧 -axis, then 𝐹 1 2 is the only nonvanishing component of 𝐹 𝑖 𝑗 . The Maxwell's equations

𝐹 [ ] 1 𝑖 𝑗 ; 𝑘 = 0 , 𝜇 𝐹 𝑖 𝑗 ; 𝑗 = 0 , ( 2 . 7 ) require that 𝐹 1 2 is the function of 𝑥 -alone. We assume that the magnetic permeability is the functions of 𝑥 and 𝑡 both. Here the semicolon represents a covariant differentiation.

The field equations (in gravitational units 𝑐 = 1 , 𝐺 = 1 ), in normal gauge for Lyra's manifold, were obtained by Sen [48] as

𝑅 𝑖 𝑗 1 2 𝑔 𝑖 𝑗 3 𝑅 + 2 𝜙 𝑖 𝜙 𝑗 3 4 𝑔 𝑖 𝑗 𝜙 𝑘 𝜙 𝑘 = 8 𝜋 𝑇 𝑖 𝑗 , ( 2 . 8 ) where 𝜙 𝑖 is the displacement field vector defined as

𝜙 𝑖 = ( 0 , 0 , 0 , 𝛽 ( 𝑡 ) ) , ( 2 . 9 ) where other symbols have their usual meaning as in Riemannian geometry.

For the line-element (2.1), the field (2.8) with (2.2) and (2.9) leads to the following system of equations:

1 𝐴 2 ̈ 𝐵 𝐵 ̈ 𝐶 𝐶 + ̇ 𝐴 𝐴 ̇ 𝐵 𝐵 + ̇ 𝐶 𝐶 ̇ 𝐵 ̇ 𝐶 + 𝐵 𝐵 𝐶 𝐶 3 𝐵 𝐶 4 𝛽 2 𝐹 = 8 𝜋 𝑝 + 2 1 2 2 𝜇 𝐴 2 𝐵 2 1 , ( 2 . 1 0 ) 𝐴 2 ̇ 𝐴 2 𝐴 2 ̈ 𝐴 𝐴 ̈ 𝐶 𝐶 + 𝐶 𝐶 3 4 𝛽 2 𝐹 = 8 𝜋 𝑝 + 2 1 2 2 𝜇 𝐴 2 𝐵 2 1 , ( 2 . 1 1 ) 𝐴 2 ̇ 𝐴 2 𝐴 2 ̈ 𝐴 𝐴 ̈ 𝐵 𝐵 + 𝐵 𝐵 3 4 𝛽 2 𝐹 = 8 𝜋 𝑝 2 1 2 2 𝜇 𝐴 2 𝐵 2 1 , ( 2 . 1 2 ) 𝐴 2 𝐵 𝐵 𝐶 𝐶 + ̇ 𝐴 𝐴 ̇ 𝐵 𝐵 + ̇ 𝐶 𝐶 𝐵 𝐶 + ̇ 𝐵 ̇ 𝐶 𝐵 𝐶 + 3 𝐵 𝐶 4 𝛽 2 𝐹 = 8 𝜋 𝜌 + 2 1 2 2 𝜇 𝐴 2 𝐵 2 ̇ 𝐵 , ( 2 . 1 3 ) 𝐵 + ̇ 𝐶 𝐶 ̇ 𝐴 𝐴 𝐵 𝐵 + 𝐶 𝐶 = 0 . ( 2 . 1 4 ) Here, and also in the following expressions, a dot and a dash indicate ordinary differentiation with, respect to 𝑡 and 𝑥 respectively.

The energy conservation equation 𝑇 𝑖 𝑖 ; 𝑗 = 0 leads to

2 ̇ 𝐴 ̇ 𝜌 + ( 𝜌 + 𝑝 ) 𝐴 + ̇ 𝐵 𝐵 + ̇ 𝐶 𝐶 𝑅 = 0 , ( 2 . 1 5 ) 𝑗 𝑖 1 2 𝑔 𝑗 𝑖 𝑅 ; 𝑗 + 3 2 𝜙 𝑖 𝜙 𝑗 ; 𝑗 3 4 𝑔 𝑗 𝑖 𝜙 𝑘 𝜙 𝑘 ; 𝑗 = 0 . ( 2 . 1 6 ) Equation (2.16) leads to 3 2 𝜙 𝑖 𝜕 𝜙 𝑗 𝜕 𝑥 𝑗 + 𝜙 𝑙 Γ 𝑗 𝑙 𝑗 + 3 2 𝜙 𝑗 𝜕 𝜙 𝑖 𝜕 𝑥 𝑗 𝜙 𝑙 Γ 𝑙 𝑖 𝑗 3 4 𝑔 𝑗 𝑖 𝜙 𝑘 𝜕 𝜙 𝑘 𝜕 𝑥 𝑗 + 𝜙 𝑙 Γ 𝑘 𝑙 𝑗 3 4 𝑔 𝑗 𝑖 𝜙 𝑘 𝜕 𝜙 𝑘 𝜕 𝑥 𝑗 + 𝜙 𝑙 Γ 𝑙 𝑘 𝑗 = 0 . ( 2 . 1 7 ) Equation (2.17) is identically satisfied for 𝑖 = 1 , 2 , 3 . For 𝑖 = 4 , (2.17) reduces to

3 2 𝛽 ̇ 3 𝛽 + 2 𝛽 2 2 ̇ 𝐴 𝐴 + ̇ 𝐵 𝐵 + ̇ 𝐶 𝐶 = 0 . ( 2 . 1 8 )

3. Solution of Field Equations in Presence of Magnetic Field

Equations (2.10)–(2.14) are five independent equations in seven unknowns 𝐴 , 𝐵 , 𝐶 , 𝜌 , 𝑝 , 𝛽 , and 𝐹 1 2 . For the complete determinacy of the system, we need two extra conditions which are narrated hereinafter. The research on exact solutions is based on some physically reasonable restrictions used to simplify the field equations.

To get determinate solution, we assume that the expansion 𝜃 in the model is proportional to the shear 𝜎 . This condition leads to

𝐵 𝐴 = 𝐶 𝑛 , ( 3 . 1 ) where 𝑛 is a constant. The motive behind assuming this condition is explained with reference to Thorne [90]; the observations of the velocity-red-shift relation for extragalactic sources suggest that Hubble expansion of the universe is isotropic today within 3 0 percent [91, 92]. To put more precisely, red-shift studies place the limit

𝜎 𝐻 0 . 3 ( 3 . 2 ) on the ratio of shear, 𝜎 , to Hubble constant, 𝐻 , in the neighbourhood of our Galaxy today. Collins et al. [93] have pointed out that for spatially homogeneous metric, the normal congruence to the homogeneous expansion satisfies that the condition 𝜎 / 𝜃 is constant.

From (2.15)–(2.17), we have

𝐴 4 4 𝐴 𝐴 2 4 𝐴 2 + 𝐴 4 𝐵 4 + 𝐴 𝐴 𝐵 4 𝐶 4 𝐵 𝐴 𝐶 4 4 𝐵 𝐵 4 𝐶 4 = 𝐶 𝐵 𝐶 1 1 𝐶 𝐵 1 𝐶 1 𝐵 𝐶 = K ( c o n s t a n t ) , ( 3 . 3 ) 8 𝜋 𝐹 2 1 2 𝜇 𝐵 2 𝐶 = 4 4 𝐶 + 𝐶 1 1 𝐶 + 𝐵 4 4 𝐵 𝐵 1 1 𝐵 . ( 3 . 4 ) We also assume that

𝐵 = 𝑓 ( 𝑥 ) 𝑔 ( 𝑡 ) , 𝐶 = 𝑓 ( 𝑥 ) 𝑘 ( 𝑡 ) . ( 3 . 5 ) Using (3.1) and (3.5) in (2.14) and (3.3) leads to

𝑘 4 𝑘 = ( 2 𝑛 1 ) 𝑔 ( 2 𝑛 + 1 ) 4 𝑔 , ( 𝑔 ( 3 . 6 ) 𝑛 1 ) 4 4 𝑔 𝑘 𝑛 4 4 𝑘 𝑔 4 𝑔 𝑘 4 𝑘 = 𝐾 , ( 3 . 7 ) 𝑓 𝑓 1 1 𝑓 2 1 = 𝐾 𝑓 2 . ( 3 . 8 ) Equation (3.6) leads to

𝑘 = 𝑐 0 𝑔 𝛼 , ( 3 . 9 ) where 𝛼 = ( 2 𝑛 1 ) / ( 2 𝑛 + 1 ) and 𝑐 0 is the constant of integration. From (3.7) and (3.9), we have

𝑔 4 4 𝑔 𝑔 + 2 4 𝑔 2 = 𝑁 , ( 3 . 1 0 ) where

= 𝑛 𝛼 ( 𝛼 1 ) + 𝛼 𝐾 𝑛 ( 𝛼 1 ) + 1 , 𝑁 = . 𝑛 ( 1 𝛼 ) 1 ( 3 . 1 1 ) Equation (3.8) leads to

1 𝑓 = e x p 2 𝐾 𝑥 + 𝑥 0 2 , ( 3 . 1 2 ) where 𝑥 0 is an integrating constant. Equation (3.10) leads to

𝑐 𝑔 = 1 𝑒 𝑏 𝑡 + 𝑐 2 𝑒 𝑏 𝑡 1 / ( + 1 ) , ( 3 . 1 3 ) where 𝑏 = ( + 1 ) 𝑁 and 𝑐 1 , 𝑐 2 are integrating constants. Hence from (3.9) and (3.13), we have

𝑘 = 𝑐 0 𝑐 1 𝑒 𝑏 𝑡 + 𝑐 2 𝑒 𝑏 𝑡 𝛼 / ( + 1 ) . ( 3 . 1 4 ) Therefore, we obtain

1 𝐵 = e x p 2 𝐾 𝑥 + 𝑥 0 2 𝑐 1 𝑒 𝑏 𝑡 + 𝑐 2 𝑒 𝑏 𝑡 1 / ( + 1 ) , 1 𝐶 = e x p 2 𝐾 𝑥 + 𝑥 0 2 𝑐 0 𝑐 1 𝑒 𝑏 𝑡 + 𝑐 2 𝑒 𝑏 𝑡 𝛼 / ( + 1 ) , 𝑐 𝐴 = 𝑎 1 𝑒 𝑏 𝑡 + 𝑐 2 𝑒 𝑏 𝑡 𝑛 ( 1 𝛼 ) / ( + 1 ) , ( 3 . 1 5 ) where 𝑎 = 𝑐 3 / 𝑐 0 , 𝑐 3 being a constant of integration.

After using suitable transformation of the coordinates, the model (2.1) reduces to the form

𝑑 𝑠 2 = 𝑎 2 𝑐 1 𝑒 𝑏 𝑇 + 𝑐 2 𝑒 𝑏 𝑇 2 𝑛 ( 1 𝛼 ) / ( + 1 ) 𝑑 𝑋 2 𝑑 𝑇 2 + 𝑒 𝐾 𝑋 2 𝑐 1 𝑒 𝑏 𝑇 + 𝑐 2 𝑒 𝑏 𝑇 2 / ( + 1 ) 𝑑 𝑌 2 + 𝑒 𝐾 𝑋 2 𝑐 1 𝑒 𝑏 𝑇 + 𝑐 2 𝑒 𝑏 𝑇 2 𝛼 / ( + 1 ) 𝑑 𝑍 2 , ( 3 . 1 6 ) where 𝑥 + 𝑥 0 = 𝑋 , 𝑡 = 𝑇 , 𝑦 = 𝑌 , and 𝑐 0 𝑧 = 𝑍 .

4. Some Physical and Geometric Properties of the Model in Presence of Magnetic Field

Equation (2.18) gives

̇ 𝛽 𝛽 2 ̇ 𝐴 = 𝐴 + ̇ 𝐵 𝐵 + ̇ 𝐶 𝐶 , a s 𝛽 0 , ( 4 . 1 ) which leads to

̇ 𝛽 𝛽 = 𝑏 { 2 𝑛 ( 1 𝛼 ) + 1 + 𝛼 } 𝑐 ( + 1 ) 1 𝑒 𝑏 𝑇 𝑐 2 𝑒 𝑏 𝑇 𝑐 1 𝑒 𝑏 𝑇 + 𝑐 2 𝑒 𝑏 𝑇 . ( 4 . 2 ) Equation (4.2) on integration gives

𝛽 = 𝐾 0 𝑐 1 𝑒 𝑏 𝑇 + 𝑐 2 𝑒 𝑏 𝑇 𝜅 , ( 4 . 3 ) where 𝐾 0 is a constant of integration and

𝜅 = 𝑏 { 2 𝑛 ( 𝛼 1 ) ( 𝛼 + 1 ) } . ( + 1 ) ( 4 . 4 ) Using (3.15) and (4.3) in (2.10) and (2.13), the expressions for pressure 𝑝 and density 𝜌 for the model (3.16) are given by

1 8 𝜋 𝑝 = 𝑎 2 𝜓 2 2 𝑛 ( 1 𝛼 ) / ( + 1 ) 𝐾 2 𝑋 2 2 ( 3 + 𝛼 ) 𝑏 2 𝑐 1 𝑐 2 ( + 1 ) 𝜓 2 2 2 𝑛 𝛼 2 + 𝛼 2 𝑏 + 2 𝛼 2 𝑛 + 3 2 2 ( + 1 ) 2 𝜓 2 1 𝜓 2 2 3 4 𝐾 2 0 𝜓 2 2 𝜅 , 1 8 𝜋 𝜌 = 𝑎 2 𝜓 2 2 𝑛 ( 1 𝛼 ) / ( + 1 ) 3 𝐾 2 𝑋 2 2 𝐾 + 2 𝑏 2 ( 𝛼 1 ) 𝑐 1 𝑐 2 ( + 1 ) 𝜓 2 2 2 𝑛 𝛼 2 𝛼 2 𝑏 2 𝛼 2 𝑛 + 1 2 2 ( + 1 ) 2 𝜓 2 1 𝜓 2 2 + 3 4 𝐾 2 0 𝜓 2 2 𝜅 , ( 4 . 5 ) where

𝜓 1 = 𝑐 1 𝑒 𝑏 𝑇 𝑐 2 𝑒 𝑏 𝑇 , 𝜓 2 = 𝑐 1 𝑒 𝑏 𝑇 + 𝑐 2 𝑒 𝑏 𝑇 . ( 4 . 6 ) From (3.4) the nonvanishing component 𝐹 1 2 of the electromagnetic field tensor is obtained as

𝐹 2 1 2 = 𝜇 𝑏 8 𝜋 2 ( 1 𝛼 ) ( + 1 ) 2 𝑒 𝐾 𝑋 2 𝜓 2 2 / ( + 1 ) 4 ( + 1 ) 𝑐 1 𝑐 2 + ( 1 + 𝛼 ) 𝜓 2 1 𝜓 2 2 . ( 4 . 7 ) From the above, equation it is observed that the electromagnetic field tensor increases with time.

The reality conditions (Ellis [94])

( i ) 𝜌 + 𝑝 > 0 , ( i i ) 𝜌 + 3 𝑝 > 0 , ( 4 . 8 ) lead to

𝑏 2 𝑛 𝑛 𝛼 2 1 ( + 1 ) 2 𝜓 2 1 𝜓 2 2 4 𝑏 2 𝑐 1 𝑐 2 ( + 1 ) 𝜓 2 2 > 𝐾 𝐾 𝑋 2 , 𝑏 + 1 ( 4 . 9 ) 2 4 𝑛 4 𝑛 𝛼 2 𝛼 2 2 𝛼 5 ( + 1 ) 2 𝜓 2 1 𝜓 2 2 4 𝑏 2 ( 𝛼 + 5 ) 𝑐 1 𝑐 2 ( + 1 ) 𝜓 2 2 3 > 2 𝐾 + 2 𝛽 2 𝑎 2 𝜓 2 2 𝑛 ( 1 𝛼 ) / ( + 1 ) , ( 4 . 1 0 ) respectively.

The dominant energy conditions (Hawking and Ellis [95])

( i ) 𝜌 𝑝 0 , ( i i ) 𝜌 + 𝑝 0 , ( 4 . 1 1 ) lead to

𝑏 2 ( 𝛼 + 1 ) 2 ( + 1 ) 2 𝜓 2 1 𝜓 2 2 + 4 𝑏 2 ( 𝛼 + 1 ) 𝑐 1 𝑐 2 ( + 1 ) 𝜓 2 2 + 3 2 𝛽 2 𝑎 2 𝜓 2 2 𝑛 ( 1 𝛼 ) / ( + 1 ) 2 𝐾 2 𝐾 𝑋 2 𝑏 + 1 , ( 4 . 1 2 ) 2 𝑛 𝑛 𝛼 2 1 ( + 1 ) 2 𝜓 2 1 𝜓 2 2 4 𝑏 2 𝑐 1 𝑐 2 ( + 1 ) 𝜓 2 2 𝐾 𝐾 𝑋 2 , + 1 ( 4 . 1 3 ) respectively. The conditions (4.10) and (4.12) impose a restriction on displacement vector 𝛽 ( 𝑡 ) .

The expressions for the expansion 𝜃 , Hubble parameter 𝐻 , shear scalar 𝜎 2 , deceleration parameter 𝑞 , and proper volume 𝑉 3 for the model (3.16) are given by

𝐻 = 3 𝜃 = 3 𝑏 { 𝑛 ( 1 𝛼 ) + ( 1 + 𝛼 ) } ( + 1 ) 𝑎 𝜓 2 𝑛 ( 1 𝛼 ) / ( + 1 ) 𝜓 1 𝜓 2 , 𝜎 ( 4 . 1 4 ) 2 = 𝑏 2 { 𝑛 ( 1 𝛼 ) + ( 1 + 𝛼 ) } 2 3 𝑛 ( 1 𝛼 ) ( 1 + 𝛼 ) 3 𝛼 3 ( + 1 ) 2 𝑎 2 𝜓 2 2 𝑛 ( 1 𝛼 ) / ( + 1 ) 𝜓 2 1 𝜓 2 2 , ( 4 . 1 5 ) 𝑞 = 1 6 𝑐 1 𝑐 2 ( + 1 ) 𝑛 1 𝛼 2 𝜓 2 1 𝑉 , ( 4 . 1 6 ) 3 = 𝑔 = 𝑎 2 𝜓 2 2 𝑛 ( 1 + 𝛼 ) ( 1 𝛼 ) / ( + 1 ) 𝑒 𝐾 𝑋 2 . ( 4 . 1 7 ) From (4.14) and (4.15), we obtain

𝜎 2 𝜃 2 = { 𝑛 ( 1 𝛼 ) + ( 1 + 𝛼 ) } 2 3 𝑛 1 𝛼 2 3 𝛼 3 { 𝑛 ( 1 𝛼 ) + ( 1 + 𝛼 ) } 2 = c o n s t a n t . ( 4 . 1 8 ) The rotation 𝜔 is identically zero.

The rate of expansion 𝐻 𝑖 in the direction of 𝑥 , 𝑦 , and 𝑧 is given by

𝐻 𝑥 = 𝐴 4 𝐴 = 𝑛 𝑏 ( 1 𝛼 ) 𝜓 ( + 1 ) 1 𝜓 2 , 𝐻 𝑦 = 𝐵 4 𝐵 = 𝑏 ( 𝜓 + 1 ) 1 𝜓 2 , 𝐻 𝑧 = 𝐶 4 𝐶 = 𝑏 𝛼 𝜓 ( + 1 ) 1 𝜓 2 . ( 4 . 1 9 ) Generally the model (3.16) represents an expanding, shearing, and nonrotating universe in which the flow vector is geodetic. The model (3.16) starts expanding at 𝑇 > 0 and goes on expanding indefinitely when 𝑛 ( 1 𝛼 ) / ( + 1 ) < 0 . Since 𝜎 / 𝜃 = constant, the model does not approach isotropy. As 𝑇 increases the proper volume also increases. The physical quantities 𝑝 and 𝜌 decrease as 𝐹 1 2 increases. However, if 𝑛 ( 1 𝛼 ) / ( 𝛽 + 1 ) > 0 , the process of contraction starts at 𝑇 > 0 , and at 𝑇 = the expansion stops. The electromagnetic field tensor does not vanish when 𝑏 0 , and 𝛼 1 . It is observed from (4.16) that 𝑞 < 0 when 𝑐 1 > 0 and 𝑐 2 > 0 which implies an accelerating model of the universe. Recent observations of type Ia supernovae [96100] reveal that the present universe is in accelerating phase and deceleration parameter lies somewhere in the range 1 < 𝑞 0 . It follows that our models of the universe are consistent with recent observations. Either when 𝑐 1 = 0 or 𝑐 2 = 0 , the deceleration parameter 𝑞 approaches the value ( 1 ) as in the case of de-Sitter universe.

5. Field Equations and Their Solution in Absence of Magnetic Field

In absence of magnetic field, the field (2.8) with (2.2) and (2.9) for metric (2.1) reads as

1 𝐴 2 𝐵 4 4 𝐵 𝐶 4 4 𝐶 + 𝐴 4 𝐴 𝐵 4 𝐵 + 𝐶 4 𝐶 𝐵 4 𝐶 4 + 𝐵 𝐵 𝐶 1 𝐶 1 3 𝐵 𝐶 = 8 𝜋 𝑝 + 4 𝛽 2 1 , ( 5 . 1 ) 𝐴 2 𝐴 2 4 𝐴 2 𝐴 4 4 𝐴 𝐶 4 4 𝐶 + 𝐶 1 1 𝐶 3 = 8 𝜋 𝑝 + 4 𝛽 2 1 , ( 5 . 2 ) 𝐴 2 𝐴 2 4 𝐴 2 𝐴 4 4 𝐴 𝐵 4 4 𝐵 + 𝐵 1 1 𝐵 3 = 8 𝜋 𝑝 + 4 𝛽 2 1 , ( 5 . 3 ) 𝐴 2 𝐵 1 1 𝐵 𝐶 1 1 𝐶 + 𝐴 4 𝐴 𝐵 4 𝐵 + 𝐶 4 𝐶 𝐵 1 𝐶 1 + 𝐵 𝐵 𝐶 4 𝐶 4 3 𝐵 𝐶 = 8 𝜋 𝜌 4 𝛽 2 𝐵 , ( 5 . 4 ) 1 4 𝐵 + 𝐶 1 4 𝐶 𝐴 4 𝐴 𝐵 1 𝐵 + 𝐶 1 𝐶 = 0 . ( 5 . 5 ) Equations (5.2) and (5.3) lead to

𝐵 4 4 𝐵 𝐵 1 1 𝐵 𝐶 4 4 𝐶 + 𝐶 1 1 𝐶 = 0 . ( 5 . 6 ) Equations (3.5) and (5.6) lead to

𝑔 4 4 𝑔 𝑘 4 4 𝑘 = 0 . ( 5 . 7 ) Equations (3.9) and (5.7) lead to

𝑔 4 4 𝑔 𝑔 + 𝛼 2 4 𝑔 2 = 0 , ( 5 . 8 ) which on integration gives

𝑐 𝑔 = 4 𝑡 + 𝑐 5 1 / ( 𝛼 + 1 ) , ( 5 . 9 ) where 𝑐 4 and 𝑐 5 are constants of integration. Hence from (3.9) and (5.9), we have

𝑘 = 𝑐 0 𝑐 4 𝑡 + 𝑐 5 𝛼 / ( 𝛼 + 1 ) . ( 5 . 1 0 ) In this case (3.8) also leads to the same as (3.12).

Therefore, in absence of magnetic field, we have

1 𝐵 = e x p 2 𝐾 𝑥 + 𝑥 0 2 𝑐 4 𝑡 + 𝑐 5 1 / ( 𝛼 + 1 ) , 1 𝐶 = e x p 2 𝐾 𝑥 + 𝑥 0 2 𝑐 𝑐 4 𝑡 + 𝑐 5 𝛼 / ( 𝛼 + 1 ) , 𝑐 𝐴 = 𝑎 4 𝑡 + 𝑐 5 𝑛 ( 1 𝛼 ) / ( 1 + 𝛼 ) , ( 5 . 1 1 ) where 𝑎 is already defined in previous section.

After using suitable transformation of the coordinates, the metric (2.1) reduces to the form 𝑑 𝑠 2 = 𝑎 2 𝑐 4 𝑇 2 𝑛 ( 1 𝛼 ) / ( 1 + 𝛼 ) 𝑑 𝑋 2 𝑑 𝑇 2 + 𝑒 𝐾 𝑋 2 𝑐 4 𝑇 2 / ( 𝛼 + 1 ) 𝑑 𝑌 2 + 𝑒 𝐾 𝑋 2 𝑐 4 𝑇 2 𝛼 / ( 𝛼 + 1 ) 𝑑 𝑍 2 , ( 5 . 1 2 ) where 𝑥 + 𝑥 0 = 𝑋 , 𝑦 = 𝑌 , 𝑐 0 𝑧 = 𝑍 , and 𝑡 + 𝑐 5 / 𝑐 4 = 𝑇 .

6. Some Physical and Geometric Properties of the Model in Absence of Magnetic Field

With the use of (5.11), equation (2.18) leads to

̇ 𝛽 𝛽 = 1 𝑇 2 𝑛 ( 𝛼 1 ) ( 𝛼 + 1 ) ( 𝛼 + 1 ) , a s 𝛽 0 , ( 6 . 1 ) which upon integration leads to

𝛽 = 𝑇 ( ( 2 𝑛 ( 𝛼 1 ) ( 𝛼 + 1 ) ) / ( 𝛼 + 1 ) ) , ( 6 . 2 ) where is an integrating constant.

Using (5.11) and (6.2) in (5.1) and (5.4), the expressions for pressure 𝑝 and density 𝜌 for the model (5.12) are given by

1 8 𝜋 𝑝 = 𝑎 2 𝑐 4 𝑇 2 𝑛 ( 1 𝛼 ) / ( 1 + 𝛼 ) 𝑛 1 𝛼 2 + 𝛼 ( 𝛼 + 1 ) 2 1 𝑇 2 + 𝐾 2 𝑋 2 3 4 2 𝑇 2 ( ( 2 𝑛 ( 𝛼 1 ) ( 𝛼 + 1 ) ) / ( 𝛼 + 1 ) ) , 1 8 𝜋 𝜌 = 𝑎 2 𝑐 4 𝑇 2 𝑛 ( 1 𝛼 ) / ( 1 + 𝛼 ) 𝑛 1 𝛼 2 + 𝛼 ( 𝛼 + 1 ) 2 1 𝑇 2 𝐾 2 + 3 𝐾 𝑋 2 + 3 4 2 𝑇 2 ( ( 2 𝑛 ( 𝛼 1 ) ( 𝛼 + 1 ) ) / ( 𝛼 + 1 ) ) . ( 6 . 3 ) The dominant energy conditions (Hawking and Ellis [95])

( i ) 𝜌 𝑝 0 , ( i i ) 𝜌 + 𝑝 0 , ( 6 . 4 ) lead to

3 4 𝛽 2 𝑎 2 𝑐 4 𝑇 2 𝑛 ( 1 𝛼 ) / ( 1 + 𝛼 ) 𝐾 1 + 2 𝐾 𝑋 2 , 𝑛 ( 6 . 5 ) 1 𝛼 2 + 𝛼 ( 1 + 𝛼 ) 2 1 𝑇 2 𝐾 1 + 𝐾 𝑋 2 , ( 6 . 6 ) respectively.

The reality conditions (Ellis [94])

( i ) 𝜌 + 𝑝 > 0 , ( i i ) 𝜌 + 3 𝑝 > 0 , ( 6 . 7 ) lead to

𝑛 1 𝛼 2 + 𝛼 ( 1 + 𝛼 ) 2 1 𝑇 2 > 𝐾 1 + 𝐾 𝑋 2 2 𝑛 , ( 6 . 8 ) 1 𝛼 2 + 𝛼 ( 1 + 𝛼 ) 2 1 𝑇 2 3 > 𝐾 + 4 𝛽 2 𝑐 4 𝑇 2 𝑛 ( 1 𝛼 ) / ( 1 + 𝛼 ) . ( 6 . 9 ) The conditions (6.5) and (6.9) impose a restriction on 𝛽 ( 𝑡 ) .

The expressions for the expansion 𝜃 , Hubble parameter 𝐻 , shear scalar 𝜎 2 , deceleration parameter 𝑞 , and proper volume 𝑉 3 for the model (5.12) in absence of magnetic field are given by

𝐻 = 3 𝜃 = 𝑛 ( 1 𝛼 ) + ( 1 + 𝛼 ) 𝑎 ( 1 + 𝛼 ) 𝑐 4 𝑛 ( 1 𝛼 ) / ( 1 + 𝛼 ) 1 𝑇 ( 𝑛 ( 1 𝛼 ) + ( 1 + 𝛼 ) ) / ( 1 + 𝛼 ) , 𝜎 ( 6 . 1 0 ) 2 = { 𝑛 ( 1 𝛼 ) + ( 1 + 𝛼 ) } 2 3 𝑛 1 𝛼 2 3 𝛼 3 𝑎 2 ( 1 + 𝛼 ) 2 𝑐 4 𝑛 ( 1 𝛼 ) / ( 1 + 𝛼 ) 1 𝑇 ( 2 𝑛 ( 1 𝛼 ) + 2 ( 1 + 𝛼