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

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 𝐹₁₂ 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. [15–18] 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 [22–31]. 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 [34–38]. 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. [40–44] 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 [61–64], 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 [66–68] 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. [69–73], 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 [84–87] 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.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]

𝐸^𝑗𝑖=𝜇ℎ_𝑙ℎ^𝑙𝑢^𝑖𝑢^𝑗+12𝑔^𝑗𝑖−ℎ^𝑖ℎ^𝑗,(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 𝑢¹ = 0 = 𝑢² = 𝑢³ and 𝑢⁴=1/𝐴. If we consider that the current flows along the 𝑧-axis, then 𝐹₁₂ is the only nonvanishing component of 𝐹_𝑖𝑗. The Maxwell's equations

𝐹1[𝑖𝑗;𝑘]=0,𝜇𝐹^𝑖𝑗_;𝑗=0,(2.7) require that 𝐹₁₂ 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

𝑅_𝑖𝑗−12𝑔^𝑖𝑗𝑅+32𝜙^𝑖𝜙_𝑗−34𝑔^𝑖𝑗𝜙_𝑘𝜙^𝑘=−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𝐴²−𝐵𝐵̈−𝐶𝐶̈+𝐴𝐴̇𝐵𝐵̇+𝐶𝐶̇−𝐵𝐵̇𝐶𝐶̇+𝐵^𝐶^𝐵𝐶−34𝛽²=8𝜋𝑝+𝐹²¹²12𝜇𝐴²𝐵²,(2.10)𝐴²𝐴𝐴̇²²−𝐴̈𝐴−𝐶𝐶̈+𝐶𝐶^−34𝛽²=8𝜋𝑝+𝐹²¹²12𝜇𝐴²𝐵²,(2.11)𝐴²𝐴𝐴̇²²−𝐴̈𝐴−𝐵𝐵̈+𝐵𝐵^−34𝛽²=8𝜋𝑝−𝐹²¹²12𝜇𝐴²𝐵²,(2.12)𝐴²−𝐵𝐵^−𝐶𝐶^+𝐴𝐴̇𝐵𝐵̇+𝐶𝐶̇−𝐵^𝐶^𝐵𝐶+𝐵𝐵̇𝐶𝐶̇+34𝛽²=8𝜋𝜌+𝐹²¹²𝐵̇2𝜇𝐴²𝐵²,(2.13)_𝐵+𝐶𝐶̇^−𝐴𝐴̇𝐵𝐵^+𝐶𝐶^=0.(2.14) 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.15)_𝑗𝑖−12𝑔^𝑗𝑖𝑅_;𝑗+32𝜙^𝑖𝜙^𝑗_;𝑗−34𝑔^𝑗𝑖𝜙_𝑘𝜙^𝑘_;𝑗=0.(2.16) Equation (2.16) leads to 32𝜙^𝑖𝜕𝜕𝜙𝑥^𝑗𝑗+𝜙^𝑙Γ^𝑗𝑙𝑗+23𝜙^𝑗𝜕𝜕𝑥𝜙^𝑖𝑗−𝜙_𝑙Γ^𝑙𝑖𝑗−43𝑔^𝑗𝑖𝜙_𝑘𝜕𝜕𝜙𝑥^𝑘𝑗+𝜙^𝑙Γ^𝑘𝑙𝑗−43𝑔^𝑗𝑖𝜙^𝑘𝜕𝜕𝜙𝑥^𝑘𝑗+𝜙_𝑙Γ^𝑙𝑘𝑗=0.(2.17) Equation (2.17) is identically satisfied for 𝑖=1,2,3. For 𝑖=4, (2.17) reduces to

32𝛽𝛽̇+23𝛽²2𝐴𝐴̇+𝐵𝐵̇+𝐶𝐶̇=0.(2.18)

3. Solution of Field Equations in Presence of Magnetic Field

Equations (2.10)–(2.14) are five independent equations in seven unknowns 𝐴, 𝐵, 𝐶, 𝜌, 𝑝, 𝛽, and 𝐹₁₂. 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 ≈30 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

𝐴₄₄𝐴−𝐴𝐴²⁴²+𝐴⁴𝐵₄𝐴𝐵+𝐴⁴𝐶₄𝐴𝐶−𝐵𝐵⁴⁴−𝐵⁴𝐶₄𝐵𝐶=𝐶𝐶¹¹−𝐵¹𝐶₁8𝜋𝐹₂₁₂𝐵𝐶=K(constant),(3.3)𝜇𝐵²=−𝐶⁴⁴𝐶+𝐶𝐶¹¹+𝐵𝐵⁴⁴−𝐵𝐵¹¹.(3.4) We also assume that

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

𝑘₄(𝑛−1)𝑔𝑘₄₄=((22𝑛𝑛−+11))𝑔𝑔⁴,(3.6)𝑓𝑔𝑓₁₁−𝑛𝑘𝑘⁴⁴−𝑔𝑔⁴𝑘𝑘⁴=𝐾,(3.7)−𝑓²¹=𝐾𝑓².(3.8) Equation (3.6) leads to

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

𝑔₄₄𝑔+ℓ𝑔𝑔²⁴²=𝑁,(3.10) where

ℓ=𝑛𝛼(𝛼−1)+𝛼𝑛(𝛼−1)+1,𝑁=𝑛(1−𝐾𝛼)−1.(3.11) Equation (3.8) leads to

𝑓=exp12𝐾𝑥+𝑥⁰²,(3.12) where 𝑥₀ is an integrating constant. Equation (3.10) leads to

𝑔=𝑐₁𝑒^𝑏𝑡+𝑐₂𝑒^−𝑏𝑡^1/(ℓ+1),(3.13) where 𝑏=√(ℓ+1)𝑁 and 𝑐₁, 𝑐₂ are integrating constants. Hence from (3.9) and (3.13), we have

𝑘=𝑐₀𝑐₁𝑒^𝑏𝑡+𝑐₂𝑒^−𝑏𝑡^𝛼/(ℓ+1).(3.14) Therefore, we obtain

𝐵=exp12𝐾𝑥+𝑥⁰²𝑐¹𝐶=exp1𝑒^𝑏𝑡+𝑐₂𝑒^−𝑏𝑡^1/(ℓ+1),2𝐾𝑥+𝑥⁰²𝑐⁰𝑐₁𝐴=𝑎𝑐₁𝑒^𝑏𝑡+𝑐₂𝑒^−𝑏𝑡^𝛼/(ℓ+1),𝑒^𝑏𝑡+𝑐₂𝑒^−𝑏𝑡^{𝑛(1−𝛼)/(ℓ+1)},(3.15) where 𝑎=𝑐₃/𝑐₀, 𝑐₃ being a constant of integration.

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

𝑑𝑠²=+𝑎²𝑒𝑐₁𝑒^𝑏𝑇_𝐾𝑋2+𝑐₂𝑒^−𝑏𝑇^{2𝑛(1−𝛼)/(ℓ+1)}𝑑𝑋²−𝑑𝑇²𝑐₁𝑒^𝑏𝑇+𝑐₂𝑒^−𝑏𝑇^2/(ℓ+1)𝑑𝑌²+𝑒^𝐾𝑋2𝑐₁𝑒^𝑏𝑇+𝑐₂𝑒^−𝑏𝑇^{2𝛼/(ℓ+1)}𝑑𝑍²,(3.16) where 𝑥+𝑥₀=𝑋, 𝑡=𝑇, 𝑦=𝑌, and 𝑐₀𝑧=𝑍.

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

Equation (2.18) gives

𝛽̇𝛽=−2𝐴𝐴̇+𝐵𝐵̇+𝐶𝐶̇,as𝛽≠0,(4.1) which leads to

𝛽̇𝛽=−𝑏{2𝑛(1(−ℓ+𝛼)1+)1+𝛼}𝑐¹𝑒^𝑏𝑇−𝑐₂𝑒^−𝑏𝑇𝑐₁𝑒^𝑏𝑇+𝑐₂𝑒^−𝑏𝑇.(4.2) Equation (4.2) on integration gives

𝛽=𝐾₀𝑐₁𝑒^𝑏𝑇+𝑐₂𝑒^−𝑏𝑇^𝜅,(4.3) where 𝐾₀ 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

8𝜋𝑝=1𝑎²𝜓^{22𝑛(1−𝛼)/(ℓ+1)}𝐾²𝑋²−2(3+𝛼)𝑏²𝑐₁𝑐₂(ℓ+1)𝜓²²−2𝑛𝛼²+𝛼²+2𝛼−2𝑛+3𝑏²2(ℓ+1)²𝜓²¹𝜓²²−34𝐾²⁰8𝜋𝜌=1𝜓^22𝜅,𝑎²𝜓^{22𝑛(1−𝛼)/(ℓ+1)}−3𝐾²𝑋²−2𝐾+2𝑏²(𝛼−1)𝑐₁𝑐₂(ℓ+1)𝜓²²−2𝑛𝛼²−𝛼²−2𝛼−2𝑛+1𝑏²2(ℓ+1)²𝜓²¹+3𝜓²²4𝐾²⁰𝜓^22𝜅,(4.5) where

𝜓𝜓₁₂=𝑐₁𝑒^𝑏𝑇−𝑐₂𝑒^−𝑏𝑇,=𝑐₁𝑒^𝑏𝑇+𝑐₂𝑒^−𝑏𝑇.(4.6) From (3.4) the nonvanishing component 𝐹₁₂ of the electromagnetic field tensor is obtained as

𝐹²¹²=𝜇8𝜋𝑏²(1−𝛼)(ℓ+1)²𝑒^𝐾𝑋2𝜓^22/(ℓ+1)4(ℓ+1)𝑐¹𝑐₂+(1+𝛼)𝜓²¹𝜓²².(4.7) From the above, equation it is observed that the electromagnetic field tensor increases with time.

The reality conditions (Ellis [94])

((iii))𝜌𝜌++𝑝3𝑝>>00,,(4.8) lead to

𝑏²𝑛−𝑛𝛼²−1(ℓ+1)²𝜓²¹𝜓²²−4𝑏²𝑐₁𝑐₂(ℓ+1)𝜓²²𝑏>𝐾𝐾𝑋²+1,(4.9)₂4𝑛−4𝑛𝛼²−𝛼²−2𝛼−5(ℓ+1)²𝜓²¹𝜓²²−4𝑏²(𝛼+5)𝑐₁𝑐₂(ℓ+1)𝜓²²>2𝐾+32𝛽²𝑎²𝜓^{22𝑛(1−𝛼)/(ℓ+1)},(4.10) respectively.

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

((iii))𝜌𝜌+−𝑝𝑝≥≥00,,(4.11) lead to

𝑏²(𝛼+1)²(ℓ+1)²𝜓²¹𝜓²²+4𝑏²(𝛼+1)𝑐₁𝑐₂(ℓ+1)𝜓²²+32𝛽²𝑏𝑎²𝜓^{22𝑛(1−𝛼)/(ℓ+1)}≥2𝐾2𝐾𝑋²+1,(4.12)₂𝑛−𝑛𝛼²−1(ℓ+1)²𝜓²¹𝜓²²−4𝑏²𝑐₁𝑐₂(ℓ+1)𝜓²²≥𝐾𝐾𝑋²+1,(4.13) respectively. The conditions (4.10) and (4.12) impose a restriction on displacement vector 𝛽(𝑡).

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

𝐻=3𝜃=3𝑏{𝑛(1−𝛼)+(1+𝛼)}(ℓ+1)𝑎𝜓^{2𝑛(1−𝛼)/(ℓ+1)}𝜓₁𝜓₂𝜎,(4.14)₂=𝑏²{𝑛(1−𝛼)+(1+𝛼)}²−3𝑛(1−𝛼)(1+𝛼)−3𝛼3(ℓ+1)²𝑎²𝜓^{22𝑛(1−𝛼)/(ℓ+1)}𝜓²¹𝜓²²𝑞=−1−6𝑐₁,(4.15)𝑐₂(ℓ+1)𝑛1−𝛼²𝜓²¹𝑉,(4.16)₃=√−𝑔=𝑎²𝜓^{22𝑛(1+𝛼)(1−𝛼)/(ℓ+1)}𝑒^𝐾𝑋2.(4.17) From (4.14) and (4.15), we obtain

𝜎²𝜃²={𝑛(1−𝛼)+(1+𝛼)}²−3𝑛1−𝛼²−3𝛼3{𝑛(1−𝛼)+(1+𝛼)}²=constant.(4.18) The rotation 𝜔 is identically zero.

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

𝐻_𝑥=𝐴⁴𝐴=𝑛𝑏(ℓ(1+−1𝛼))𝜓𝜓¹₂𝐻_𝑦,=𝐵⁴𝐵=(ℓ+𝑏1)𝜓𝜓¹₂𝐻_𝑧,=𝐶⁴𝐶=(ℓ𝑏+𝛼1)𝜓𝜓¹₂.(4.19) 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 𝐹₁₂ 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 𝑐₁>0 and 𝑐₂>0 which implies an accelerating model of the universe. Recent observations of type Ia supernovae [96–100] 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 𝑐₁=0 or 𝑐₂=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𝐴²−𝐵𝐵⁴⁴−𝐶𝐶⁴⁴+𝐴𝐴⁴𝐵𝐵⁴+𝐶𝐶⁴−𝐵⁴𝐶₄𝐵𝐶+𝐵¹𝐶₁𝐵𝐶=8𝜋𝑝+34𝛽²1,(5.1)𝐴²𝐴𝐴²⁴²−𝐴⁴⁴𝐴−𝐶𝐶⁴⁴+𝐶𝐶¹¹=8𝜋𝑝+43𝛽²1,(5.2)𝐴²𝐴𝐴²⁴²−𝐴⁴⁴𝐴−𝐵𝐵⁴⁴+𝐵𝐵¹¹=8𝜋𝑝+34𝛽²1,(5.3)𝐴²−𝐵𝐵¹¹−𝐶𝐶¹¹+𝐴𝐴⁴𝐵𝐵⁴+𝐶𝐶⁴−𝐵¹𝐶₁𝐵𝐶+𝐵⁴𝐶₄𝐵𝐶=8𝜋𝜌−43𝛽²𝐵,(5.4)₁₄𝐵+𝐶𝐶¹⁴−𝐴𝐴⁴𝐵𝐵¹+𝐶𝐶¹=0.(5.5) Equations (5.2) and (5.3) lead to

𝐵₄₄𝐵−𝐵𝐵¹¹−𝐶𝐶⁴⁴+𝐶𝐶¹¹=0.(5.6) Equations (3.5) and (5.6) lead to

𝑔₄₄𝑔−𝑘𝑘⁴⁴=0.(5.7) Equations (3.9) and (5.7) lead to

𝑔₄₄𝑔+𝛼𝑔𝑔²⁴²=0,(5.8) which on integration gives

𝑔=𝑐₄𝑡+𝑐₅^1/(𝛼+1),(5.9) where 𝑐₄ and 𝑐₅ are constants of integration. Hence from (3.9) and (5.9), we have

𝑘=𝑐₀𝑐₄𝑡+𝑐₅^{𝛼/(𝛼+1)}.(5.10) In this case (3.8) also leads to the same as (3.12).

Therefore, in absence of magnetic field, we have

𝐵=exp12𝐾𝑥+𝑥⁰²𝑐⁴𝐶=exp1𝑡+𝑐₅^1/(𝛼+1),2𝐾𝑥+𝑥⁰²𝑐𝑐⁴𝐴=𝑎𝑐₄𝑡+𝑐₅^{𝛼/(𝛼+1)},𝑡+𝑐₅^{𝑛(1−𝛼)/(1+𝛼)},(5.11) where 𝑎 is already defined in previous section.

After using suitable transformation of the coordinates, the metric (2.1) reduces to the form 𝑑𝑠²=𝑎²𝑐₄𝑇^{2𝑛(1−𝛼)/(1+𝛼)}𝑑𝑋²−𝑑𝑇²+𝑒^𝐾𝑋2𝑐₄𝑇^2/(𝛼+1)𝑑𝑌²+𝑒^𝐾𝑋2𝑐₄𝑇^{2𝛼/(𝛼+1)}𝑑𝑍²,(5.12) where 𝑥+𝑥₀=𝑋, 𝑦=𝑌, 𝑐₀𝑧=𝑍, and 𝑡+𝑐₅/𝑐₄=𝑇.

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),as𝛽≠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

8𝜋𝑝=1𝑎²𝑐₄𝑇^{2𝑛(1−𝛼)/(1+𝛼)}𝑛1−𝛼²+𝛼(𝛼+1)²𝑇1²+𝐾²𝑋²−34ℜ²8𝜋𝜌=1𝑇^{2((2𝑛(𝛼−1)−(𝛼+1))/(𝛼+1))},𝑎²𝑐₄𝑇^{2𝑛(1−𝛼)/(1+𝛼)}𝑛1−𝛼²+𝛼(𝛼+1)²𝑇1²−𝐾2+3𝐾𝑋²+34ℜ²𝑇^{2((2𝑛(𝛼−1)−(𝛼+1))/(𝛼+1))}.(6.3) The dominant energy conditions (Hawking and Ellis [95])

((iii))𝜌𝜌+−𝑝𝑝≥≥00,,(6.4) lead to

34𝛽²𝑛𝑎²1𝑐−₄𝑇𝛼^{2𝑛(1−𝛼)/(1+𝛼)}₂≥𝐾1+2𝐾𝑋²,(6.5)+𝛼(1+𝛼)²𝑇1²≥𝐾1+𝐾𝑋²,(6.6) respectively.

The reality conditions (Ellis [94])

((iii))𝜌𝜌++𝑝3𝑝>>00,,(6.7) lead to

𝑛1−𝛼²+𝛼(1+𝛼)²𝑇1²2𝑛1−𝛼₂>𝐾1+𝐾𝑋²,(6.8)+𝛼(1+𝛼)²1𝑇²>𝐾+34𝛽²𝑐₄𝑇^{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 𝜎², deceleration parameter 𝑞, and proper volume 𝑉³ 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.10)₂={𝑛(1−𝛼)+(1+𝛼)}²−3𝑛1−𝛼²−3𝛼3𝑎²(1+𝛼)²𝑐^{4𝑛(1−𝛼)/(1+𝛼)}1𝑇^{(2𝑛(1−𝛼)+2(1+𝛼))/(1+𝛼)}𝑞=−1+3(𝛼+1),(6.11)𝑉2𝑛(1−𝛼)+2(1+𝛼),(6.12)₃=√−𝑔=𝑎²𝑒^𝐾𝑋2𝑐₄𝑇^{(2𝑛(1−𝛼)+(1+𝛼))/(1+𝛼)}.(6.13) From (6.10) and (6.11), we obtain

𝜎²𝜃²={𝑛(1−𝛼)+(1+𝛼)}²−3𝑛1−𝛼²−3𝛼3{𝑛(1−𝛼)+(1+𝛼)}²=constant.(6.14) The rotation 𝜔 is identically zero.

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

𝐻_𝑥=𝐴⁴𝐻_𝑦𝐴=𝑛((11+−𝛼𝛼))𝑇1,=𝐵⁴𝐻_𝑧𝐵=(1+1𝛼)𝑇1,=𝐶⁴𝐶=(1+𝛼𝛼)𝑇1.(6.15) The model (5.12) starts expanding with a big bang at 𝑇=0 and it stops expanding at 𝑇=∞. It should be noted that the universe exhibits initial singularity of the Point-type at 𝑇=0. The space-time is well behaved in the range 0<𝑇<𝑇₀. In absence of magnetic field, the model represents a shearing and nonrotating universe in which the flow vector is geodetic. At the initial moment 𝑇=0, the parameters 𝜌, 𝑝, 𝛽, 𝜃, 𝜎² and 𝐻 tend to infinity. So the universe starts from initial singularity with infinite energy density, infinite internal pressure, infinitely large gauge function, infinite rate of shear and expansion. Moreover, 𝜌, 𝑝, 𝛽, 𝜃, 𝜎² and 𝐻 are monotonically decreasing toward a nonzero finite quantity for 𝑇 in the range 0<𝑇<𝑇₀ in absence of magnetic field. Since 𝜎/𝜃 = constant, the model does not approach isotropy. As 𝑇 increases the proper volume also increases. It is observed that for the derived model, the displacement vector 𝛽(𝑡) is a decreasing function of time and therefore it behaves like cosmological term Λ. It is observed from (6.12) that 𝑞<0 when 𝛼<(2𝑛−1)/(2𝑛+1) which implies an accelerating model of the universe. When 𝛼=−1, the deceleration parameter 𝑞 approaches the value (−1) as in the case of de-Sitter universe. Thus, also in absence of magnetic field, our models of the universe are consistent with recent observations.

7. Discussion and Concluding Remarks

In this paper, we have obtained a new class of exact solutions of Einstein's modified field equations for cylindrically symmetric space-time with perfect fluid distribution within the framework of Lyra's geometry both in presence and absence of magnetic field. The solutions are obtained using the functional separability of the metric coefficients. The source of the magnetic field is due to an electric current produced along the 𝑧-axis. 𝐹₁₂ is the only nonvanishing component of electromagnetic field tensor. The electromagnetic field tensor is given by (4.7), 𝜇 remains undetermined as function of both 𝑥 and 𝑡. The electromagnetic field tensor does not vanish if 𝑏≠0 and 𝛼≠1. It is observed that in presence of magnetic field, the rate of expansion of the universe is faster than that in absence of magnetic field. The idea of primordial magnetism is appealing because it can potentially explain all the large-scale fields seen in the universe today, specially those found in remote proto-galaxies. As a result, the literature contains many studies examining the role and the implications of magnetic fields for cosmology. In presence of magnetic field, the model (3.16) represents an expanding, shearing and nonrotating universe in which the flow vector is geodetic. But in the absence of magnetic field, the model (5.12) found that in the universe all the matter and radiation are concentrated at the big bang epoch and the cosmic expansion is driven by the big bang impulse. The universe has singular origin and it exhibits power-law expansion after the big bang impulse. The rate of expansion slows down and finally stops at 𝑇→∞. In absence of magnetic field, the pressure, energy density and displacement field become zero whereas the spatial volume becomes infinitely large as 𝑇→∞.

It is possible to discuss entropy in our universe. In thermodynamics the expression for entropy is given by

𝑇𝑑𝑆=𝑑𝜌𝑉³+𝑝𝑑𝑉³,(7.1) where 𝑉³=𝐴²𝐵𝐶 is the proper volume in our case. To solve the entropy problem of the standard model, it is necessary to treat 𝑑𝑆>0 for at least a part of evolution of the universe. Hence (7.1) reduces to

𝑇𝑑𝑆=𝜌₄+(𝜌+𝑝)2𝐴⁴𝐴+𝐵𝐵⁴+𝐶𝐶⁴>0.(7.2) The conservation equation 𝑇^{𝑗𝑖∶𝑗}=0 for (2.1) leads to

𝜌₄+(𝜌+𝑝)𝐴⁴𝐴+𝐵𝐵⁴+𝐶𝐶⁴+23𝛽𝛽⁴+32𝛽²2𝐴𝐴⁴+𝐵𝐵⁴+𝐶𝐶⁴=0.(7.3) Therefore, (7.1) and (7.2) lead to

32𝛽𝛽⁴+32𝛽²2𝐴𝐴⁴+𝐵𝐵⁴+𝐶𝐶⁴<0,(7.4) which gives to 𝛽<0. Thus, the displacement vector 𝛽(𝑡) affects entropy because for entropy 𝑑𝑆>0 leads to 𝛽(𝑡)<0.

In spite of homogeneity at large scale, our universe is inhomogeneous at small scale, so physical quantities being position-dependent are more natural in our observable universe if we do not go to super high scale. This result shows this kind of physical importance. It is observed that the displacement vector 𝛽(𝑡) coincides with the nature of the cosmological constant Λ which has been supported by the work of several authors as discussed in the physical behaviour of the model in Sections 4 and 6. In the recent time Λ-term has attracted theoreticians and observers for many a reason. The nontrivial role of the vacuum in the early universe generates a Λ-term that leads to inflationary phase. Observationally, this term provides an additional parameter to accommodate conflicting data on the values of the Hubble constant, the deceleration parameter, the density parameter and the age of the universe (see, e.g., [101, 102]). In recent past there is an upsurge of interest in scalar fields in general relativity and alternative theories of gravitation in the context of inflationary cosmology [103–105]. Therefore the study of cosmological models in Lyra's geometry may be relevant for inflationary models. There seems a good possibility of Lyra's geometry to provide a theoretical foundation for relativistic gravitation, astrophysics, and cosmology. However, the importance of Lyra's geometry for astrophysical bodies is still an open question. In fact, it needs a fair trial for experiment.

Acknowledgments

The authors would like to thank the Harish-Chandra Research Institute, Allahabad, India for local hospitality where this work is done. The authors also thank the referee for his fruitful comments.