Next Article in Journal
Synthesis of 2,6-Diaminotriptycene Conjugates with Chiral Auxiliaries: Towards the Scalable Resolution of Crucial Triptycene Intermediates
Next Article in Special Issue
Wave Surface Symmetry and Petrov Types in General Relativity
Previous Article in Journal
Modeling of Particle Size Distribution in the Presence of Flocculant
Previous Article in Special Issue
Dynamics of Particles with Electric Charge and Magnetic Dipole Moment near Schwarzschild-MOG Black Hole
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Noether Symmetry Analysis of the Klein–Gordon and Wave Equations in Bianchi I Spacetime

Department of Chemistry and Physics, Roger Williams University, One Old Ferry Road, Bristol, RI 02809, USA
Symmetry 2024, 16(1), 115; https://doi.org/10.3390/sym16010115
Submission received: 1 December 2023 / Revised: 3 January 2024 / Accepted: 16 January 2024 / Published: 18 January 2024
(This article belongs to the Special Issue Noether and Space-Time Symmetries in Physics—Volume Ⅱ)

Abstract

:
We investigate the Noether symmetries of the Klein–Gordon Lagrangian for Bianchi I spacetime. This is accomplished using a set of new Noether symmetry relations for the Klein–Gordon Lagrangian of Bianchi I spacetime, which reduces to the wave equation in a special case. A detailed Noether symmetry analysis of the Klein–Gordon and the wave equations for Bianchi I spacetime is presented, and the corresponding conservation laws are derived.

1. Introduction

A conformal Killing vector (CKV) K has to satisfy
£ K g i j = 2 σ ( x k ) g i j ,
where g i j is the metric tensor, £ K is the Lie derivative operator along K , and σ ( x k ) is a conformal factor. When σ ; i j 0 , the CKV field is said to be proper [1]. The vector field K is called a special conformal Killing vector (SCKV) field if σ ; i j = 0 ; a homothetic Killing vector (HKV) field if σ , i = 0 , e.g., σ is a constant on the manifold; and a Killing vector (KV) field if σ = 0 , which is also called the isometry of spacetime. The set of all CKVs (respectively SCKV, HKV, and KV) forms a finite-dimensional Lie algebra. The maximum dimensions of the CKV algebra on the manifold M is fifteen if M is conformally flat, and it is seven if the spacetime is not conformally flat. The physical features of differential equations, in terms of the conservation laws admitted by them, are directly associated with the Noether symmetries, which is facilitated using a Lagrangian of the corresponding dynamical system.As we show in the following section, there is a direct relation between the conformal symmetries and Noether symmetries.
If a Lagrangian L for a given dynamical system exhibits symmetry, this property is strongly related to Noether symmetries, which describe the physical characteristics of differential equations associated with a Lagrangian L , in terms of the first integrals they possess [2,3]. This relationship can be viewed from two different perspectives. First, one can take a strict Noether symmetry approach [4,5,6], which results in £ X L = 0 , where £ X is the Lie derivative operator along X . On the other hand, one can employ the Noether symmetry approach with a gauge term [7,8,9,10,11,12,13], a generalization of the strict Noether symmetry approach where the Noether symmetry equation includes the gauge term. Noether symmetries with a gauge term are equally valuable in addressing a variety of problems in physics and applied mathematics. In the following section, we will discuss the relationship between Noether symmetries with a gauge term of the Klein–Gordon Lagrangian and the geometric symmetries of spacetimes. The study of differential equations involving geometry is an active area of research. Recent literature [14,15] has delved into the connection between geometrical structures and conserved quantities. Noether symmetries are directly linked to conserved quantities or conservation laws [3], which naturally emerge in a wide range of applications.
The cosmological principle assumes that the universe is homogeneous and isotropic at large scales, and the geometrical model that satisfies these properties is Friedmann–Lemaître–Robertson–Walker (FLRW) spacetime. However, there has been a suggestion from studies of the cosmic microwave background (CMB) temperature anisotropies that the assumption of statistical isotropy is violated at the largest angular scales, leading to some intriguing anomalies [16]. In order to make predictions for the CMB anisotropies, one can explore cosmological models that are homogeneous but anisotropic, such as the Bianchi-type spacetimes, which encompass both isotropic and homogeneous FLRW models.
Bianchi classified all three-dimensional real Lie algebras and demonstrated that there are nine possible simply transitive groups of motions, denoted as G 3 . These groups are generated using Killing vectors K α , where α = 1 , 2 , 3 , with structure constants C α β γ , defined by K α , K β = C α β γ K γ . The structure constants C α β γ can be decomposed into irreducible parts, as follows:
C α β γ = ϵ α β σ N σ γ + A α δ β γ A β δ α γ
here N α β is symmetric. It follows from the Jacobi identities that N α β A β = 0 . One can diagonalize N α β without loss of generality, choosing A α in the 1-direction when it is non-zero. The real Lie algebras can then be classified into nine types of spatially homogeneous Bianchi spacetimes, distinguished by the particular form of the structure constants C α β γ , based on whether A α = 0 (Class A) or A α 0 (Class B). This classification results in the Bianchi types, such as I, I I , V I 0 , V I I 0 , V I I I , and I X for Class A models, and I I I , I V , V, V I h , and V I I h for Class B models. It is important to note that the Bianchi models include FLRW models as special cases, with Bianchi I (flat), Bianchi V (open), and Bianchi IX (closed) representing the FLRW models. In this study, we consider the Bianchi I metric as the background spacetime in the Klein–Gordon equation. The line element for Bianchi I spacetime can be written as follows [17]:
d s 2 = d t 2 + A ( t ) 2 d x 2 + B ( t ) 2 d y 2 + C ( t ) 2 d z 2 .
The above spacetime yields a flat FLRW metric if A = B = C . For any form of the metric coefficients A , B , and C, it can easily be computed from Equation (1) that the KVs for Bianchi I spacetime are K 1 = x , K 2 = y , and K 3 = z . In this study, all vector fields are written with bold letters.
The rest of the paper is organized as follows: In the next section, Section 2, we present an analysis of Noether symmetries with a gauge term for the Klein–Gordon Lagrangian in the context of the Bianchi I spacetime model. In Section 3, we apply the Noether symmetry approach to the Bianchi I spacetime. In Section 4, we study the field equations of the Bianchi I spacetime with an imperfect fluid. Conclusions and discussions are presented in the final section, Section 5.

2. Klein–Gordon Lagrangian and Noether Symmetry Equations

The Klein–Gordon equation in the Riemannian space with metric g i k is a second-order partial differential equation of the form
ψ = G ( x k , ψ ) ,
where ψ = g i k ψ ; i k and □ refers to the de’Alembertian or Laplace operator defined by = 1 g x i g g i k x k in terms of the Riemannian space. The Klein–Gordon equation given in (3) follows from the first-order Lagrangian, which is called the Klein–Gordon Lagrangian,
L = g 1 2 g i j ψ i ψ j F ( x k , ψ ) ,
where ψ i ( x k ) ψ ( x k ) x i and F ( x k , ψ ) = G ( x k , ψ ) . Here, the prime represents the derivative with respect to the scalar field ψ . The above Lagrangian reduces to the Lagrangian of the wave equation when F = 0 .
The wave equation of spacetimes is one of the most important equations in physics, and it is common to study this equation in terms of the Lie and Noether symmetry generators they admit [18,19,20,21,22,23]. When symmetry generators exist, they play a crucial role in finding exact solutions. Jamal et al. [21] investigated the wave equation of Bianchi III spacetime, calculating and classifying Noether symmetries and constructing corresponding conservation laws. They also obtained reductions of the wave equation and identified some invariant solutions. In the papers [22,23], the authors utilized the invariance and multiplier method [24,25] to conduct conservation law classifications of the wave equation for the Bianchi I spacetime, using power law metric coefficients.
Noether symmetries of the Klein–Gordon equations for well-known spacetimes have been calculated and classified according to their symmetry generators [15,26,27,28,29,30]. In [26], the symmetry analysis of the Klein–Gordon equation in de Sitter spacetime was classified, and the obtained symmetries were utilized to find exact solutions using quadratures. Paliathanasis et al. [27] conducted a classification of Lie and Noether symmetries for the Klein–Gordon and wave equations in pp-wave spacetimes. The symmetry properties and conservation laws of wave and Gordon-type equations in Milne and Bianchi III spacetimes were investigated in the papers [28,29]. In [15], a geometric procedure was employed for the symmetry classification of the Klein–Gordon equation in Bianchi I spacetime, connecting the Noether symmetries with a gauge function of the Klein–Gordon Lagrangian to the conformal symmetries of the metric tensor. Their study extended the results of Bozhkov and Freire [31] for the Klein–Gordon equation with a constant potential, i.e., V ( x k ) = V 0 . More recently, Paliathanasis [30] considered the Klein–Gordon equations for the conformal forms of Bianchi I, Bianchi III, and Bianchi V spacetimes, deriving closed-form expressions for potential functions that admit Lie and Noether symmetries of the Klein–Gordon equations. In this work, we aim to derive Noether symmetries with a gauge term for the dynamical Lagrangian L of the Klein–Gordon equation within the background of Bianchi I spacetimes. To obtain the Klein–Gordon equation in this spacetime, we construct a Lagrangian model. Using the obtained Lagrangian for the Klein–Gordon equation in Bianchi I spacetime, we calculate and classify Noether symmetry generators with gauge terms. Furthermore, we identify conservation laws provided by the Lagrangian for representing the Klein–Gordon equation.
In a general program of research into the problem of integrating the classical and quantum equations of motion of a test particle in external fields of different nature in spaces with symmetry following the sets of Killing fields, Obukhov found all admissible electromagnetic fields for the case, when the groups of motions G 3 act simply transitively on the hypersurfaces of spacetime V 4 [32,33,34,35]. In [32], he found all external electromagnetic fields in which the Klein–Gordon–Fock equation admits the first-order symmetry operators and completed the classification of admissible electromagnetic fields in which the Hamilton–Jacobi and Klein–Gordon–Fock equations admit algebras of motion integrals that are isomorphic to the algebras of operators of the r-parametric groups of motions, G r , of spacetime manifolds if r 4 [33]. In the paper [34], the case when the groups G 4 act on V 3 was considered. The remaining case in the latter article, when the groups G 4 act simply transitively on the space V 4 , was studied in the paper [35].
The Noether symmetry generator for the Klein–Gordon Lagrangian (4) is
X = ξ i ( x k , ψ ) x i + Φ ( x k , ψ ) ψ
if there a gauge function f i ( x k , ψ ) exists and the Noether symmetry condition
X [ 1 ] L + L ( D i ξ i ) = D i f i
is satisfied, where
D i = x i + ψ i ψ + ψ i k ψ k +
is the total derivative operator and X [ 1 ] is the first prolongation of Noether symmetry generator X , i.e.,
X [ 1 ] = X + ζ i ( x k , ψ , ψ k ) ψ i
where ζ i ( x k , ψ , ψ k ) = D i Φ ψ j D i ξ j . The corresponding Noether flow T i is defined by the expression
T i = ξ i L + ξ k ψ k Φ L ψ i + f i ,
which is called the conserved vector  T = ( T 1 , , T n ) , where i = 1 , , n , or if n = 1 , then it is called the conserved quantity. The Noether flow (7) satisfies the local conservation law
D i T i = 0 .
It is crucial to discover conservation laws when studying physical systems. For studying differential equations, conservation laws are useful for integrability, linearization, analyzing solutions, and understanding constants of motion. Noether’s theorem enables the derivation of all local conservation laws for a (system of) differential equation(s) derived from a Lagrangian. This method helps resolve the problem of calculating conservation laws for given differential equation(s), as Noether’s theorem provides a formula utilizing symmetries of the action to derive these local conservation laws. One can find a detailed discussion about the physical significance of Noether symmetries in [36].
For the Klein–Gordon Lagrangian (4), we obtain the first prolongation of the Noether symmetry generator X as
X [ 1 ] L = g { Φ F , ψ F ξ i Γ i k k F , i ξ i + 1 2 £ ξ g m n + ξ i Γ i k k + 2 Φ , ψ g m n ψ m ψ n + g i m Φ , i ψ m ξ , ψ j ψ m ψ i ψ j } ,
where £ ξ is the Lie derivative operator along ξ = ξ k / x k . Putting (9) into (5) together with D i ξ i = ξ , i i + ξ , ψ i ψ i and D i f i = f , i i + f , ψ i ψ i , the Noether symmetry condition (5) gives rise to
ξ , ψ i = 0 , g g i j Φ , j f , ψ i = 0 ,
£ ξ g i j = ξ ; k k + 2 Φ , ψ g i j ,
g £ ξ F + F ξ ; i i + Φ F , ψ + f , i i = 0 .
Thus, we find the geometrical form of Noether symmetry Equation (5) in terms of Lie derivatives of the metric tensor. Here, Equation (10) yields ξ i = ξ i ( x k ) and f i = g g i j Φ , j d ψ + A i ( x k ) , where A i ( x k ) is an integration function. If ξ i is a CKV with conformal factor σ ( x k ) , then Equations (1) and (11) imply ξ ; i i = 2 σ 2 Φ , ψ . Finally, the Noether symmetry condition (12) becomes
F , i ξ i + 2 F ( σ Φ , ψ ) + Φ F , ψ + ( Φ ) d ψ + A , i i g = 0 .
For the Bianchi I spacetime (2), the Klein–Gordon Lagrangian (4) has the form
L = A B C 2 ψ t 2 + B C 2 A ψ x 2 + A C 2 B ψ y 2 + A B 2 C ψ z 2 A B C F ( t , x , y , z , ψ ) ,
which is the the Lagrangian of the wave equation when F = 0 . The Klein–Gordon Equation (3) is obtained through variation of this Lagrangian with respect to the scalar field ψ , as the following:
1 A 2 ψ x x + 1 B 2 ψ y y + 1 C 2 ψ z z ψ t t ( A B C ) A B C ψ t = G ( t , x , y , z , ψ ) ,
where G ( t , x , y , z , ψ ) = F ψ and the dot represents the derivative with respect to time t.
Let us consider the Noether symmetry generator for the Klein–Gordon Lagrangian (14) as follows:
X = ξ 0 t + ξ 1 x + ξ 2 y + ξ 3 z + Φ ψ ,
where the components of ξ = ( ξ 0 , ξ 1 , ξ 2 , ξ 3 ) and Φ are dependent on t , x , y , z and ψ . Now, we seek the dependent variables ξ 0 , ξ 1 , ξ 2 , ξ 3 , Φ that will be solved from the geometrical Noether symmetry conditions (10)–(12), in order that the Lagrangian (14) would admit any Noether symmetry. For the Bianchi I spacetime (2), the geometrical Noether symmetry conditions (10)–(12) yield 19 PDEs:
ξ , ψ 0 = 0 , ξ , ψ 1 = 0 , ξ , ψ 2 = 0 , ξ , ψ 3 = 0 , A B C Φ , t + f , ψ 0 = 0 , B C A Φ , x f , ψ 1 = 0 , A C B Φ , y f , ψ 2 = 0 , A B C Φ , z f , ψ 3 = 0 , A 2 ξ , t 1 ξ , x 0 = 0 , B 2 ξ , t 2 ξ , y 0 = 0 , C 2 ξ , t 3 ξ , z 0 = 0 , A 2 ξ , y 1 + B 2 ξ , x 2 = 0 , A 2 ξ , z 1 + C 2 ξ , x 3 = 0 , B 2 ξ , z 2 + C 2 ξ , y 3 = 0 , ξ , t 0 + ξ , x 1 + ξ , y 2 + ξ , z 3 + A ˙ A + B ˙ B + C ˙ C ξ 0 + 2 Φ , ψ = 0 , ξ , t 0 ξ , x 1 + ξ , y 2 + ξ , z 3 + A ˙ A + B ˙ B + C ˙ C ξ 0 + 2 Φ , ψ = 0 , ξ , t 0 + ξ , x 1 ξ , y 2 + ξ , z 3 + A ˙ A B ˙ B + C ˙ C ξ 0 + 2 Φ , ψ = 0 , ξ , t 0 + ξ , x 1 + ξ , y 2 ξ , z 3 + A ˙ A + B ˙ B C ˙ C ξ 0 + 2 Φ , ψ = 0 , F , t ξ 0 + F , x ξ 1 + F , y ξ 2 + F , z ξ 3 + Φ F , ψ + F ξ , t 0 + ξ , x 1 + ξ , y 2 + ξ , z 3 + A ˙ A + B ˙ B + C ˙ C ξ 0 + 1 A B C f , t 0 + f , x 1 + f , y 2 + f , z 3 = 0 .
It is noted here that the set of all Noether symmetries with the gauge functions f i form a finite dimensional Lie algebra.

3. Noether Symmetries and Conservation Laws

It is easily seen from Equation (1) that for arbitrary forms of metric functions A ( t ) , B ( t ) and C ( t ) , the background spacetime of Bianchi I metric admit the three KVs, the generators of translations in x, y and z directions which implies momentum conservation,
K 1 = x , K 2 = y , K 3 = z .
It is obvious that these KVs are also Noether, and so Lie, symmetries of the Klein–Gordon Equation (15) in the background of Bianchi I spacetime. Hence, applying the expression (7) of conservation law, the resulting conserved flow vector components related to the Klein–Gordon Equation (15) are as follows:
T t = A B C ψ t ψ x , T x = A B C 2 ψ t 2 + 1 A 2 ψ x 2 1 B 2 ψ y 2 1 C 2 ψ z 2 + 2 V 1 ( t , y , z , ψ ) , T y = A C B ψ x ψ y , T z = A B C ψ y ψ z ,
for K 1 = x ,
T t = A B C ψ t ψ y , T x = B C A ψ x ψ y , T y = A B C 2 ψ t 2 1 A 2 ψ x 2 + 1 B 2 ψ y 2 1 C 2 ψ z 2 + 2 V 2 ( t , x , z , ψ ) , T z = A B C ψ y ψ z .
for K 2 = y ,
T t = A B C ψ t ψ z , T x = B C A ψ x ψ z , T y = A C B ψ y ψ z T z = A B C 2 ψ t 2 1 A 2 ψ x 2 1 B 2 ψ y 2 + 1 C 2 ψ z 2 + 2 V 3 ( t , x , y , ψ ) .
for K 3 = z . Here, the conserved flow vector is T = ( T t , T x , T y , T z ) , and V 1 ( t , y , z , ψ ) , V 2 ( t , x , z , ψ ) , and V 3 ( t , x , y , ψ ) are integration functions. We give a complete solution of Noether symmetry conditions (17) for the Bianchi I spacetime in the following.
  • Case (i): First, let us consider the Klein–Gordon equation that requires that F 0 . Taking F = U 0 + U 1 ψ + 1 2 U 2 2 ψ 2 for any metric coefficients of Bianchi I spacetime, where U 0 , U 1 and U 2 are constants, the components of the Noether symmetry generator (16) are found from Equations (17) as
ξ 0 = 0 , ξ 1 = c 1 , ξ 2 = c 2 , ξ 3 = c 3 ,
Φ = c 4 e a 1 x + c 5 e a 1 x c 6 e a 2 y + c 7 e a 2 y c 8 e a 3 z + c 9 e a 3 z Y ( t ) ,
along with the gauge vector components
f 0 = A B C c 4 e a 1 x + c 5 e a 1 x c 6 e a 2 y + c 7 e a 2 y c 8 e a 3 z + c 9 e a 3 z ψ Y ˙ ( t ) + F 0 ( t , x , y , z ) , f 1 = a 1 B C A c 4 e a 1 x c 5 e a 1 x c 6 e a 2 y + c 7 e a 2 y c 8 e a 3 z + c 9 e a 3 z ψ Y ( t ) + F 1 ( t , x , y , z ) , f 2 = a 2 A C B c 4 e a 1 x + c 5 e a 1 x c 6 e a 2 y c 7 e a 2 y c 8 e a 3 z + c 9 e a 3 z ψ Y ( t ) + F 2 ( t , x , y , z ) , f 3 = A B a 3 ψ C U 1 C a 3 c 4 e a 1 x + c 5 e a 1 x c 6 e a 2 y + c 7 e a 2 y c 8 e a 3 z c 9 e a 3 z Y ( t ) ( F 0 , t + F 1 , x + F 2 , y ) d z + F 3 ( t , x , y ) ,
where c 1 , , c 9 and a 1 , a 2 , a 3 are constant parameters, F 0 , , F 3 are integration functions, and Y ( t ) solves the following second-order ordinary differential equation:
Y ¨ + A ˙ A + B ˙ B + C ˙ C Y ˙ a 1 2 A 2 + a 2 2 B 2 + a 3 2 C 2 + U 2 2 Y ( t ) = 0 .
Therefore, there are eleven Noether symmetries, such that
X 1 = K 1 , X 2 = K 2 , X 3 = K 3 ,
X 4 = e a 1 x + a 2 y + a 3 z Y ( t ) ψ , X 5 = e a 1 x + a 2 y + a 3 z Y ( t ) ψ , X 6 = e a 1 x a 2 y + a 3 z Y ( t ) ψ ,
X 7 = e a 1 x + a 2 y a 3 z Y ( t ) ψ , X 8 = e a 1 x a 2 y a 3 z Y ( t ) ψ , X 9 = e a 1 x + a 2 y a 3 z Y ( t ) ψ ,
X 10 = e a 1 x a 2 y + a 3 z Y ( t ) ψ , X 11 = e a 1 x a 2 y a 3 z Y ( t ) ψ .
with the corresponding non-zero gauge vectors:
f 4 , 5 = e ± a 1 x + a 2 y + a 3 z A B C ψ Y ˙ , ± a 1 Y A 2 , a 2 Y B 2 , ( a 3 C 2 U 1 a 3 ψ ) Y , f 6 , 7 = e a 1 x a 2 y ± a 3 z A B C ψ Y ˙ , a 1 Y A 2 , a 2 Y B 2 , ± ( a 3 C 2 U 1 a 3 ψ ) Y , f 8 , 9 = e ± a 1 x a 2 y a 3 z A B C ψ Y ˙ , ± a 1 Y A 2 , a 2 Y B 2 , ( a 3 C 2 U 1 a 3 ψ ) Y , f 10 , 11 = e a 1 x a 2 y ± a 3 z A B C ψ Y ˙ , a 1 Y A 2 , a 2 Y B 2 , ± ( a 3 C 2 U 1 a 3 ψ ) Y .
The conserved vector fields associated with X 1 , , X 11 given in (26)–(29) are obtained as
T 1 = 1 2 A B C W K 1 ψ x T 0 , T 2 = 1 2 A B C W K 2 ψ y T 0 , T 3 = 1 2 A B C W K 3 ψ z T 0 , T 4 = e a 1 x + a 2 y + a 3 z Y T 0 + f 4 , T 5 = e a 1 x + a 2 y + a 3 z Y T 0 + f 5 , T 6 = e a 1 x a 2 y + a 3 z Y T 0 + f 6 , T 7 = e a 1 x + a 2 y a 3 z Y T 0 + f 7 , T 8 = e a 1 x a 2 y a 3 z Y T 0 + f 8 , T 9 = e a 1 x + a 2 y a 3 z Y T 0 + f 9 , T 10 = e a 1 x a 2 y + a 3 z Y T 0 + f 10 , T 11 = e a 1 x a 2 y a 3 z Y T 0 + f 11 ,
where the gauge vectors f 5 , , f 11 are the same as in (30), and W and T 0 are defined by
W = ψ t 2 + 1 A 2 ψ x 2 + 1 B 2 ψ y 2 + 1 C 2 ψ z 2 2 F , T 0 = A B C ψ t , 1 A 2 ψ x , 1 B 2 ψ y , 1 C 2 ψ z ,
where F = U 0 + U 1 ψ + 1 2 U 2 ψ 2 . The conserved vector components (7) for the integration functions F 0 , , F 3 of the gauge functions which have the property of T i = f i are
T = F 0 ( t , x , y , z ) , F 1 ( t , x , y , z ) , F 2 ( t , x , y , z ) , ( F 0 , t + F 1 , x + F 2 , y ) d z + F 3 ( t , x , y ) .
We note, here, that the above conserved quantities will appear in each of the possible cases. Therefore, we will not mention these quantities again.
Through to the end of this section, following Ref. [21], we will take into account the power-law form of metric functions, such that
A ( t ) = t L , B ( t ) = t p , C ( t ) = t q ,
where L , p , and q are constant parameters. For the latter forms of the metric functions, we find that there exists an additional HKV, which is a scaling transformation or a dilation,
K 4 = t t + ( L 1 ) x x + ( p 1 ) y y + ( q 1 ) z z ,
where σ = const . = 1 , in addition to the KVs K 1 , K 2 and K 3 in (18). Furthermore, we will consider some subcases in which we obtain the symmetry generators for the wave ( F = U 0 = c o n s t . , where one can take F = 0 without loss of generality) and Klein–Gordon equations ( F c o n s t . ) of Bianchi I spacetime.
Obviously, other choices for the metric functions will lead to a different solution for the function of Y ( t ) from Equation (25). We give some examples of these choices: (i.1)  A ( t ) = sin ( L t ) , B ( t ) = cos ( p t ) , C ( t ) = cos ( q t ) ; (i.2)  A ( t ) = sin ( L t ) , B ( t ) = sin ( p t ) , C ( t ) = sin ( q t ) ; (i.3)  A ( t ) = cos ( L t ) , B ( t ) = cos ( p t ) , C ( t ) = cos ( q t ) ; and (i.4)  A ( t ) = sech ( L t ) , B ( t ) = sech ( p t ) , C ( t ) = sech ( q t ) . In subcase (i.1), if L = p = 1 and q = 0 , then the line element reduces to the conformally flat Bianchi I spacetime, and point symmetries and potentials for the Klein–Gordon equation with F = V ( t , x , y , z ) ψ 2 / 2 in this metric were studied by Ref. [15]. The solution Y ( t ) of (25) for L = p and q = 0 is found as follows:
Y ( t ) = g 2 ( t ) p + a 1 2 p sin ( 2 p t ) [ b 1 g 1 ( t ) d 1 2 2 F 1 n 1 d 1 , n 2 d 1 ; 2 d 1 ; g 1 ( t ) + b 2 g 1 ( t ) d 2 2 2 F 1 n 1 + d 2 , n 2 + d 2 ; 2 d 2 ; g 1 ( t ) ] ,
where p 0 , d 1 = 1 a 2 / p , d 2 = 1 + a 2 / p , n 1 = ( a 1 + p 2 a 3 2 U 2 2 ) / ( 2 p ) , n 2 = ( a 1 p 2 a 3 2 U 2 2 ) / ( 2 p ) , and the functions g 1 ( t ) , g 2 ( t ) are defined as
g 1 ( t ) = 1 2 cos ( 2 p t ) + 1 , g 2 ( t ) = 1 2 cos ( 2 p t ) 1 .
The solution (36) is a new one that was not mentioned in Ref. [15] and generalizes the case 4.4 of this reference, where they were taken as L = p = 1 . For L = p = 1 and q = 3 , we have the following solution for Equation (25):
Y ( t ) = sin ( 2 t ) g 2 ( t ) a 1 2 g 3 ( t ) a 3 36 d 2 ( 4 a 3 + 3 25 4 U 2 2 ) [ 2 g 1 ( t ) ] 3 4 2 cos 2 ( 2 t ) 3 cos ( 2 t ) + 1 [ b 1 ( cos t ) d 1 H G 3 4 , α 1 , β 1 , γ 1 , δ 1 , 1 + a 1 , g 1 ( t ) + b 2 ( cos t ) d 1 H G 3 4 , α 2 , β 2 , γ 2 , δ 2 , 1 + a 1 , g 1 ( t ) ] ,
where H G is the Heun general function, the functions g 1 ( t ) , g 2 ( t ) are the same ones as in (37) with p = 1 , and we define that g 3 ( t ) = [ cos ( 2 t ) 1 / 2 ] / 2 , β 1 = 1 + ( a 1 + d 1 ) / 2 + d 2 , β 2 = 1 + ( a 1 d 1 ) / 2 + d 2 , δ 1 = 1 + d 1 , δ 2 = 1 d 1 , and
α 1 = d 1 8 ( 3 a 1 + 7 ) + a 1 2 + a 2 2 + 11 27 a 3 2 + 2 a 1 + U 2 2 7 3 + a 3 ( 1 + d 1 ) 9 d 2 ( 4 a 3 + 3 25 4 U 2 2 ) ,
α 2 = d 1 8 ( 3 a 1 + 7 ) + a 1 2 + a 2 2 + 11 27 a 3 2 + 2 a 1 + U 2 2 7 3 + a 3 ( 1 d 1 ) 9 d 2 ( 4 a 3 + 3 25 4 U 2 2 ) ,
in which d 1 and d 2 are defined as
d 1 = 1 6 9 + 4 ( 9 a 2 2 + a 3 2 ) , d 2 = 1 12 225 + 16 a 3 2 36 U 2 2 + 24 a 3 25 4 U 2 2 .
One can find other solutions of (25) for different choices of the parameters L , p and q. Some solutions of the Equation (25) for the remaining subcases (i.2)(i.4) are included in Table 1.
  • Case (ii): If F = U 0 + U 1 ψ + 1 2 U 2 2 ψ 2 , it is found for the power-law form of Bianchi I spacetime that twelve Noether symmetries appear:
X 1 1 = K 1 , X 2 1 = K 2 , X 3 1 = K 3 ,
X 4 1 = K 4 + ψ ψ , X 5 1 = X 4 , X 6 1 = X 5 , X 7 1 = X 6 ,
X 8 1 = X 7 , X 9 1 = X 8 , X 10 1 = X 9 , X 11 1 = X 10 , X 12 1 = X 11 ,
with the corresponding non-zero gauge vectors:
f 5 , 6 1 = e ± a 1 x + a 2 y + a 3 z f 0 Y ˙ , ± a 1 t 2 L Y , a 2 t 2 p Y , a 3 t 2 q U 1 a 3 ψ Y , f 7 , 8 1 = e a 1 x a 2 y ± a 3 z f 0 Y ˙ , a 1 t 2 L Y , a 2 t 2 p Y , ± a 3 t 2 q U 1 a 3 ψ Y , f 9 , 10 1 = e ± a 1 x a 2 y a 3 z f 0 Y ˙ , ± a 1 t 2 L Y , a 2 t 2 p Y , a 3 t 2 q U 1 a 3 ψ Y , f 11 , 12 1 = e a 1 x a 2 y ± a 3 z f 0 Y ˙ , a 1 t 2 L Y , a 2 t 2 p Y , ± a 3 t 2 q U 1 a 3 ψ Y ,
where f 0 = t L + p + q ψ and X 4 , , X 11 are the same as in (27)–(29). Meanwhile, in this subcase, the second-order ordinary differential Equation (25) for Y ( t ) becomes
Y ¨ + ( L + p + q ) t Y ˙ a 1 2 t 2 L + a 2 2 t 2 p + a 3 2 t 2 q + U 2 2 Y ( t ) = 0 .
The conserved flow vectors associated with the KVs X 1 1 , X 2 1 , X 3 1 and non-Killing Noether symmetries X 4 1 , , X 12 1 given in (40) and (41) are obtained as
T 1 1 = 1 2 t L + p + q W K 1 ψ x T 0 , T 2 1 = 1 2 t L + p + q W K 2 ψ y T 0 ,
T 3 1 = 1 2 t L + p + q W K 3 ψ z T 0 , T 4 1 = 1 2 t L + p + q W K 4 + ( ψ Q 0 ) T 0 ,
T 5 1 = e a 1 x + a 2 y + a 3 z Y T 0 + f 5 1 , T 6 1 = e a 1 x + a 2 y + a 3 z Y T 0 + f 6 1 , T 7 1 = e a 1 x a 2 y + a 3 z Y T 0 + f 7 1 ,
T 8 1 = e a 1 x + a 2 y a 3 z Y T 0 + f 8 1 , T 9 1 = e a 1 x a 2 y a 3 z Y T 0 + f 9 1 , T 10 1 = e a 1 x + a 2 y a 3 z Y T 0 + f 10 1 ,
T 11 1 = e a 1 x a 2 y + a 3 z Y T 0 + f 11 1 , T 12 1 = e ( a 1 x + a 2 y + a 3 z ) Y T 0 + f 12 1 ,
where W , T 0 are given in Equation (32) by taking A = t L , B = t p , C = t q , Q 0 is defined as
Q 0 = t ψ t + ( L 1 ) x ψ x + ( p 1 ) y ψ y + ( q 1 ) z ψ z ,
and the gauge vectors f 5 1 , , f 12 1 have the same forms as in (42). Furthermore, this subcase yields a thirteenth Noether symmetry, in addition to the twelve obtained above, which is X = y x x y if L = p , and q is an arbitrary constant; X = z x x z if L = q , and p is an arbitrary constant, and finally X = z y y z if p = q , and L is an arbitrary constant.
In this case, let us examine some solutions of Equation (43) for specific values of the constant powers L , p , and q of the metric functions. If we assume L = 2 , p = q = 0 , the solution of the Equation (43) is expressed in terms of double confluent Heun functions:
Y ( t ) = 1 t b 1 H D + b 2 H D d t t H D 2 ,
where b 1 , b 2 are constant parameters, and H D = H e u n D α , β , γ , δ , t 2 + 1 t 2 1 with α = 0 , β = 1 4 ( a 1 2 + a 2 2 + a 3 2 ) U 2 2 , γ = 2 ( a 1 2 a 2 2 a 3 2 U 2 2 ) , and δ = 1 4 ( a 1 2 + a 2 2 + a 3 2 ) U 2 2 . For L = 3 , p = q = 0 , we have the solution
Y ( t ) = t 4 e a 1 t 4 4 b 1 H B + b 2 H B e a 1 t 4 2 d t t 5 H B 2 ,
where H B is the Heun biconfluent function such that H B = H e u n B α , β , γ , δ , a 1 2 t 2 with α = 2 , β = 0 , γ = 0 , and δ = 1 2 a 1 ( a 3 2 + U 2 2 ) . Further analysis of Equation (43) with respect to values of L , p , and q reveals non-trivial solutions in terms of some familiar special functions. Table 2 contains some of these solutions corresponding to specific values for L , p , and q.
  • Subcase (ii.1). If we take F = U 0 + U 1 ψ in the Klein–Gordon Equation (15) for the power law form of Bianchi I spacetime (34) with L p q , where U 0 , U 1 are constants, it is found that there are fifteen Noether symmetries, as follows:
X 1 5 = K 1 , X 2 5 = K 2 , X 3 5 = K 3 ,
X 4 , 5 5 = m e a 1 x + a 2 y ± a 3 z Y ( t ) ψ , X 6 , 7 5 = m e a 1 x a 2 y ± a 3 z Y ( t ) ψ ,
X 8 , 9 5 = ( e a 1 x + m e a 1 x ) e a 2 y ± a 3 z Y ( t ) ψ , X 10 , 11 5 = ( e a 1 x + m e a 1 x ) e a 2 y ± a 3 z Y ( t ) ψ ,
X 12 , 13 5 = K 4 + X 4 , 5 5 + ψ 3 V 1 t 2 2 ( 1 + L + p + q ) 3 V 1 t 1 L p q 2 m ( 1 + L + p + q ) ψ ,
X 14 , 15 5 = K 4 + X 6 , 7 5 + ψ 3 V 1 t 2 2 ( 1 + L + p + q ) ψ ,
and the corresponding non-zero gauge vectors are
f 4 , 5 5 = m e a 1 x + a 2 y ± a 3 z f 0 Y ˙ , a 1 t 2 L Y , a 2 t 2 p Y , ± ( a 3 t 2 q U 1 a 3 ψ ) Y , f 6 , 7 5 = m e a 1 x a 2 y ± a 3 z f 0 Y ˙ , a 1 t 2 L Y , a 2 t 2 p Y , ± ( a 3 t 2 q U 1 a 3 ψ ) Y , f 8 , 9 5 = ( e a 1 x + m e a 1 x ) e a 2 y ± a 3 z f 0 Y ˙ , a 1 ( e a 1 x m e a 1 x ) ( e a 1 x + m e a 1 x ) Y , a 2 t 2 p Y , ± ( a 3 t 2 q U 1 a 3 ψ ) Y , f 10 , 11 5 = ( e a 1 x + m e a 1 x ) e a 2 y ± a 3 z f 0 Y ˙ , a 1 ( e a 1 x m e a 1 x ) ( e a 1 x + m e a 1 x ) Y , a 2 t 2 p Y , ± ( a 3 t 2 q U 1 a 3 ψ ) Y , f 12 , 13 5 = f 4 , 5 5 + 3 U 1 ψ 1 + L + p + q t 1 + L + p + q + ( 1 L p q ) 2 m t + 3 U 1 t 2 m ( 1 + L + p + q ) + t L + p + q 4 U 0 + 3 U 1 2 t 2 2 ( 1 + L + p + q ) z z , f 14 , 15 5 = f 6 , 7 5 + 3 U 1 ψ t 1 + L + p + q 1 + L + p + q t + t L + p + q 4 U 0 + 3 U 1 2 t 2 2 ( 1 + L + p + q ) z z ,
where f 0 = t L + p + q ψ , m is a non-zero constant parameter, and K 4 is the HKV given in (35). Thus, one can write the conserved flow vectors for the Noether symmetries X 4 5 , , X 15 5 as
T 4 , 5 5 = m Y e a 1 x + a 2 y ± a 3 z T 0 + f 4 , 5 5 , T 6 , 7 5 = m Y e a 1 x a 2 y ± a 3 z T 0 + f 6 , 7 5 T 8 , 9 5 = Y ( e a 1 x + m e a 1 x ) e a 2 y ± a 3 z T 0 + f 8 , 9 5 , T 10 , 11 5 = Y ( e a 1 x + m e a 1 x ) e a 2 y ± a 3 z T 0 + f 10 , 11 5 T 12 , 13 5 = m Y e a 1 x a 2 y ± a 3 z 3 U 1 t 2 2 ( 1 + L + p + q ) 3 U 1 t 1 L p q 2 m ( 1 + L + p + q ) + ψ T 0 1 2 t L + p + q W K 4 + f 12 , 13 5 T 14 , 15 5 = m Y e a 1 x + a 2 y ± a 3 z 3 U 1 t 2 2 ( 1 + L + p + q ) + ψ T 0 1 2 t L + p + q W K 4 + f 14 , 15 5 ,
where W = ψ t 2 + t 2 L ψ x 2 + t 2 p ψ y 2 + t 2 q ψ z 2 2 ( U 0 + U 1 ψ ) .
  • Subcase (ii.2). For F = 0 , which requires U 0 = U 1 = U 2 = 0 , Equation (15) reduces to the wave equation for Bianchi I spacetime. If we take L = p = q in (34), the Bianchi I spacetime yields the flat FLRW spacetime. For the latter assumption of metric coefficients, and assuming F = U 0 , we find twenty six Noether symmetries, which are the KVs X 1 1 , X 2 1 , X 3 1 given in (39), and
X 4 2 = K 4 + ψ ψ , X 5 2 = K 5 , X 6 2 = K 6 , X 7 2 = K 7 ,
X 8 2 = K 8 1 + 2 ( q 1 ) x ψ ψ , X 9 2 = K 9 1 + 2 ( q 1 ) y ψ ψ , X 10 2 = K 10 1 + 2 ( q 1 ) z ψ ψ ,
X 11 2 = M e b 1 x + b 2 y sin ( z ) ψ , X 12 2 = M e b 1 x + b 2 y sin ( z ) ψ ,
X 13 2 = M e b 1 x b 2 y sin ( z ) ψ , X 14 2 = M e ( b 1 x + b 2 y ) sin ( z ) ψ ,
X 15 2 = M e b 1 x + b 2 y cos ( z ) ψ , X 16 2 = M e b 1 x + b 2 y cos ( z ) ψ ,
X 17 2 = M e b 1 x b 2 y cos ( z ) ψ , X 18 2 = M e ( b 1 x + b 2 y ) cos ( z ) ψ ,
X 19 2 = N e b 1 x + b 2 y sin ( z ) ψ , X 20 2 = N e b 1 x + b 2 y sin ( z ) ψ ,
X 21 2 = N e b 1 x b 2 y sin ( z ) ψ , X 22 2 = N e ( b 1 x + b 2 y ) sin ( z ) ψ ,
X 23 2 = N e b 1 x + b 2 y cos ( z ) ψ , X 24 2 = N e b 1 x + b 2 y cos ( z ) ψ ,
X 25 2 = N e b 1 x b 2 y cos ( z ) ψ , X 26 2 = N e ( b 1 x + b 2 y ) cos ( z ) ψ ,
in which the vector fields K 5 , K 6 , K 7 , K 8 1 , K 9 1 and K 10 1 have the form
K 5 = y x x y , K 6 = z x x z , K 7 = z y y z ,
K 8 1 = 2 ( q 1 ) x t t + ( q 1 ) ( y y + z z ) + ( q 1 ) 2 ( x 2 y 2 z 2 ) + t 2 ( 1 q ) x ,
K 9 1 = 2 ( q 1 ) y t t + ( q 1 ) ( x x + z z ) + ( q 1 ) 2 ( x 2 + y 2 z 2 ) + t 2 ( 1 q ) y ,
K 10 1 = 2 ( q 1 ) z t t + ( q 1 ) ( x x + y y ) + ( q 1 ) 2 ( x 2 y 2 + z 2 ) + t 2 ( 1 q ) z ,
and M = M ( t ) and N = N ( t ) are defined by
M ( t ) = ( q + 1 ) b 3 t ( q + 1 ) 2 J ν 1 ( τ ) + t ( 1 3 q ) 2 J ν 2 ( τ ) , N ( t ) = ( q + 1 ) b 3 t ( q + 1 ) 2 Y ν 1 ( τ ) + t ( 1 3 q ) 2 Y ν 2 ( τ ) ,
where b 1 , b 2 and b 3 ( 0 ) are constant parameters, q 0 , 1 , 1 / 2 , ν 1 = ( q + 1 ) / ( 2 ( q 1 ) ) , ν 2 = ( q 3 ) / ( 2 ( q 1 ) ) , = b 1 2 + b 2 2 + b 3 2 , τ = ( b 3 / ( q 1 ) ) t 1 q , { J ν 1 ( τ ) , J ν 2 ( τ ) } and { Y ν 1 ( τ ) , Y ν 2 ( τ ) } are first and second kind Bessel functions, respectively. We note here that the vector fields K 5 , K 6 , K 7 in (69) are the KVs, and the other ones K 8 1 , K 9 1 , and K 10 1 in (70)–(72) are the SCKVs with the conformal factors σ = 2 ( q 1 ) x , σ = 2 ( q 1 ) y , and σ = 2 ( q 1 ) z , respectively. Furthermore, the gauge vectors for X 8 2 , , X 26 2 are
f 8 2 = ( q 1 ) t q ψ 2 x , f 9 2 = ( q 1 ) t q ψ 2 y , f 10 2 = ( q 1 ) t q ψ 2 z , f 11 , 12 2 = ψ e ± b 1 x + b 2 y sin ( z ) b 3 t q + 1 2 J ν 1 ( τ ) , t q M ( ± b 1 , b 2 , cot ( z ) ) , f 13 , 14 2 = ψ e ± b 1 x b 2 y sin ( z ) b 3 t q + 1 2 J ν 1 ( τ ) , t q M ( ± b 1 , b 2 , cot ( z ) ) , f 15 , 16 2 = ψ e ± b 1 x + b 2 y cos ( z ) b 3 t q + 1 2 J ν 1 ( τ ) , t q M ( ± b 1 , b 2 , tan ( z ) ) , f 17 , 18 2 = ψ e ± b 1 x b 2 y cos ( z ) b 3 t q + 1 2 J ν 1 ( τ ) , t q M ( ± b 1 , b 2 , tan ( z ) ) , f 19 , 20 2 = ψ e ± b 1 x + b 2 y sin ( z ) b 3 t q + 1 2 Y ν 1 ( τ ) , t q N ( ± b 1 , b 2 , cot ( z ) ) , f 21 , 22 2 = ψ e ± b 1 x b 2 y sin ( z ) b 3 t q + 1 2 Y ν 1 ( τ ) , t q N ( ± b 1 , b 2 , cot ( z ) ) , f 23 , 24 2 = ψ e ± b 1 x + b 2 y cos ( z ) b 3 t q + 1 2 Y ν 1 ( τ ) , t q N ( ± b 1 , b 2 , tan ( z ) ) , f 25 , 26 2 = ψ e ± b 1 x b 2 y cos ( z ) b 3 t q + 1 2 Y ν 1 ( τ ) , t q N ( ± b 1 , b 2 , tan ( z ) ) .
Then, in addition to the conserved flow vectors for X 1 , X 2 , X 3 , and X 4 2 obtained in (44) and (45), the conserved vector fields for the remaining Noether symmetries of this subcase become
T 5 2 = 1 2 t 3 q W K 5 + ( y ψ x x ψ y ) T 0 , T 6 2 = 1 2 t 3 q W K 6 + ( z ψ x x ψ z ) T 0 , T 7 2 = 1 2 t 3 q W K 7 + ( z ψ y y ψ z ) T 0 , T 8 2 = 1 2 t 3 q W K 8 1 + [ 2 ( q 1 ) x ψ Q 1 ] T 0 + f 8 2 , T 9 2 = 1 2 t 3 q W K 9 1 + [ 2 ( q 1 ) y ψ Q 2 ] T 0 + f 9 2 , T 10 2 = 1 2 t 3 q W K 10 1 + [ 2 ( q 1 ) z ψ Q 3 ) ] T 0 + f 10 2 , T 11 , 12 2 = M e ± b 1 x + b 2 y sin ( z ) T 0 + f 11 , 12 2 , T 13 , 14 2 = M e ± b 1 x b 2 y sin ( z ) T 0 + f 13 , 14 2 T 15 , 16 2 = M e ± b 1 x + b 2 y cos ( z ) T 0 + f 15 , 16 2 , T 17 , 18 2 = M e ± b 1 x b 2 y cos ( z ) T 0 + f 17 , 18 2 T 19 , 20 2 = N e ± b 1 x + b 2 y sin ( z ) T 0 + f 19 , 20 2 , T 21 , 22 2 = N e ± b 1 x b 2 y sin ( z ) T 0 + f 21 , 22 2 T 23 , 24 2 = N e ± b 1 x + b 2 y cos ( z ) T 0 + f 23 , 24 2 , T 25 , 26 2 = N e ± b 1 x b 2 y cos ( z ) T 0 + f 25 , 26 2
where W , T 0 are the same as given in (32), in which F = 0 and A = B = C = t q , and Q 1 , Q 2 and Q 3 are defined as the following:
Q 1 = 2 ( q 1 ) x t ψ t + ( q 1 ) ( y ψ y + z ψ z ) + ( q 1 ) 2 ( x 2 y 2 z 2 ) + t 2 ( 1 q ) ψ x ,
Q 2 = 2 ( q 1 ) y t ψ t + ( q 1 ) ( x ψ x + z ψ z ) + ( q 1 ) 2 ( x 2 + y 2 z 2 ) + t 2 ( 1 q ) ψ y ,
Q 3 = 2 ( q 1 ) z t ψ t + ( q 1 ) ( x ψ x + y ψ y ) + ( q 1 ) 2 ( x 2 y 2 + z 2 ) + t 2 ( 1 q ) ψ z .
For some values of q in this subcase, one can obtain different Noether symmetry generators from those above. Now, we will give some examples of this situation, as seen below.
If L = p = q = 1 and F = 0 , then there are again twenty six Noether symmetries, such that the KVs K 1 , K 2 , K 3 , and
X 4 3 = t t + ψ ψ , X 5 3 = K 5 , X 6 3 = K 6 , X 7 2 = K 7 ,
X 8 3 = K 8 2 x ψ ψ , X 9 3 = K 9 2 y ψ ψ , X 10 3 = K 10 2 z ψ ψ ,
X 11 3 = t 1 h 1 ( t ) e b 1 x + b 2 y sin ( z ) ψ , X 12 3 = t 1 h 1 ( t ) e b 1 x + b 2 y sin ( z ) ψ ,
X 13 3 = t 1 h 1 ( t ) e b 1 x b 2 y sin ( z ) ψ , X 14 3 = t 1 h 1 ( t ) e ( b 1 x + b 2 y ) sin ( z ) ψ ,
X 15 3 = t 1 h 1 ( t ) e b 1 x + b 2 y cos ( z ) ψ , X 16 3 = t 1 h 1 ( t ) e b 1 x + b 2 y cos ( z ) ψ ,
X 17 3 = t 1 h 1 ( t ) e b 1 x b 2 y cos ( z ) ψ , X 18 3 = t 1 h 1 ( t ) e ( b 1 x + b 2 y ) cos ( z ) ψ ,
X 19 3 = t 1 h 2 ( t ) e b 1 x + b 2 y sin ( z ) ψ , X 20 3 = t 1 h 2 ( t ) e b 1 x + b 2 y sin ( z ) ψ ,
X 21 3 = t 1 h 2 ( t ) e b 1 x b 2 y sin ( z ) ψ , X 22 3 = t 1 h 2 ( t ) e ( b 1 x + b 2 y ) sin ( z ) ψ ,
X 23 3 = t 1 h 2 ( t ) e b 1 x + b 2 y cos ( z ) ψ , X 24 3 = t 1 h 2 ( t ) e b 1 x + b 2 y cos ( z ) ψ ,
X 25 3 = t 1 h 2 ( t ) e b 1 x b 2 y cos ( z ) ψ , X 26 3 = t 1 h 2 ( t ) e ( b 1 x + b 2 y ) cos ( z ) ψ ,
where K 8 2 , K 9 2 and K 10 2 are given by
K 8 2 = x t t + ln t x , K 9 2 = y t t + ln t y , K 10 2 = z t t + ln t z ,
which are SCKVs with the conformal factors σ = x for K 8 2 , σ = y for K 9 2 , and σ = z for K 10 2 . Here, the functions h 1 ( t ) and h 2 ( t ) are of the forms h 1 ( t ) = t 1 b 3 2 and h 2 ( t ) = t 1 b 3 2 . The corresponding gauge vector fields of the above Noether symmetries are
f 8 3 = t 2 ψ 2 x , f 9 3 = t 2 ψ 2 y , f 10 3 = t 2 ψ 2 z , f 11 , 12 3 = h 1 ( t ) e ± b 1 x + b 2 y sin ( z ) ψ ( 1 1 b 3 2 ) t , ± b 1 , b 2 , cot ( z ) , f 13 , 14 3 = h 1 ( t ) e ± b 1 x b 2 y sin ( z ) ψ ( 1 1 b 3 2 ) t , ± b 1 , b 2 , cot ( z ) f 15 , 16 3 = h 1 ( t ) e ± b 1 x + b 2 y cos ( z ) ψ ( 1 1 b 3 2 ) t , ± b 1 , b 2 , tan ( z ) , f 17 , 18 3 = h 1 ( t ) e ± b 1 x b 2 y cos ( z ) ψ ( 1 1 b 3 2 ) t , ± b 1 , b 2 , tan ( z ) , f 19 , 20 3 = h 2 ( t ) e ± b 1 x + b 2 y sin ( z ) ψ ( 1 1 b 3 2 ) t , ± b 1 , b 2 , cot ( z ) , f 21 , 22 3 = h 2 ( t ) e ± b 1 x b 2 y sin ( z ) ψ ( 1 1 b 3 2 ) t , ± b 1 , b 2 , cot ( z ) , f 23 , 24 3 = h 2 ( t ) e ± b 1 x + b 2 y cos ( z ) ψ ( 1 1 b 3 2 ) t , ± b 1 , b 2 , tan ( z ) , f 25 , 26 3 = h 2 ( t ) e ± b 1 x b 2 y cos ( z ) ψ ( 1 1 b 3 2 ) t , ± b 1 , b 2 , tan ( z ) .
Using the conservation law relation (7), the conserved flow vectors for X 4 3 , , X 26 3 become
T 4 3 = 1 2 t 3 W t + ( t ψ t + ψ ) T 0 , T 5 3 = 1 2 t 3 W K 5 + ( x ψ y y ψ x ) T 0 , T 6 3 = 1 2 t 3 W K 6 + ( x ψ z ψ x ) T 0 , T 7 3 = 1 2 t 3 W K 7 + ( y ψ z z ψ y ) T 0 , T 8 3 = 1 2 t 3 W K 8 2 x ( t ψ t + ψ ) + ln t ψ x T 0 + f 8 3 , T 9 3 = 1 2 t 3 W K 9 2 y ( t ψ t + ψ ) + ln t ψ y T 0 + f 9 3 , T 10 3 = 1 2 t 3 W K 10 2 z ( t ψ t + ψ ) + ln t ψ z T 0 + f 10 3 , T 11 , 12 3 = t 1 h 1 ( t ) e ± b 1 x + b 2 y sin ( z ) T 0 + f 11 , 12 3 , T 13 , 14 3 = t 1 h 1 ( t ) e ± b 1 x b 2 y sin ( z ) T 0 + f 13 , 14 3 T 15 , 16 3 = t 1 h 1 ( t ) e ± b 1 x + b 2 y cos ( z ) T 0 + f 15 , 16 3 , T 17 , 18 3 = t 1 h 1 ( t ) e ± b 1 x b 2 y cos ( z ) T 0 + f 17 , 18 3 T 19 , 20 3 = t 1 h 2 ( t ) e ± b 1 x + b 2 y sin ( z ) T 0 + f 19 , 20 3 , T 21 , 22 3 = t 1 h 2 ( t ) e ± b 1 x b 2 y sin ( z ) T 0 + f 21 , 22 3 T 23 , 24 3 = t 1 h 2 ( t ) e ± b 1 x + b 2 y cos ( z ) T 0 + f 23 , 24 3 , T 25 , 26 3 = t 1 h 2 ( t ) e ± b 1 x b 2 y cos ( z ) T 0 + f 25 , 26 3
where T 0 and W are the same as given in (32) by taking A = B = C = t and F = 0 .
When L = p = q = 1 / 2 and F = 0 , we obtain thirty one Noether symmetries, which is the maximum number of symmetries, and these include six KVs K 1 , K 2 , K 3 , K 5 , K 6 , K 7 and the following:
X 7 4 = X 4 2 , X 8 4 = X 8 2 , X 9 4 = X 9 2 , X 10 4 = X 10 2 ,
X 11 4 = K 11 ψ 2 t ψ , X 12 4 = K 12 x ψ 2 t ψ , X 13 4 = K 13 y ψ 2 t ψ ,
X 14 4 = K 14 z ψ 2 t ψ , X 15 4 = K 15 1 2 t [ 3 t + 1 4 ( x 2 + y 2 + z 2 ) ] ψ ψ ,
X 16 4 = h 3 ( t ) e b 1 x + b 2 y sin ( β z ) ψ , X 17 4 = h 3 ( t ) e b 1 x + b 2 y sin ( β z ) ψ ,
X 18 4 = h 3 ( t ) e b 1 x b 2 y sin ( β z ) ψ , X 19 4 = h 3 ( t ) e ( b 1 x + b 2 y ) sin ( β z ) ψ ,
X 20 4 = h 7 ( t ) e b 1 x + b 2 y cos ( β z ) ψ , X 21 4 = h 3 ( t ) e b 1 x + b 2 y cos ( β z ) ψ ,
X 22 4 = h 3 ( t ) e b 1 x b 2 y cos ( β z ) ψ , X 23 4 = h 3 ( t ) e ( b 1 x + b 2 y ) cos ( β z ) ψ ,
X 24 4 = h 4 ( t ) e b 1 x + b 2 y sin ( β z ) ψ , X 25 4 = h 4 ( t ) e b 1 x + b 2 y sin ( β z ) ψ ,
X 26 4 = h 4 ( t ) e b 1 x b 2 y sin ( β z ) ψ , X 27 4 = h 4 ( t ) e ( b 1 x + b 2 y ) sin ( β z ) ψ ,
X 28 4 = h 4 ( t ) e b 1 x + b 2 y cos ( β z ) ψ , X 29 4 = h 4 ( t ) e b 1 x + b 2 y cos ( β z ) ψ ,
X 30 4 = h 4 ( t ) e b 1 x b 2 y cos ( β z ) ψ , X 31 4 = h 4 ( t ) e ( b 1 x + b 2 y ) cos ( β z ) ψ ,
where β = b 1 2 + b 2 2 b 3 2 , h 3 ( t ) = sinh ( 2 b 3 t ) / t , and h 4 ( t ) = cosh ( 2 b 3 t ) / t . Furthermore, K 11 , K 12 , K 13 , K 14 , and K 15 are defined as follows:
K 11 = t t , K 12 = t ( x t + 2 x ) , K 13 = t ( y t + 2 y ) , K 14 = t ( z t + 2 z ) ,
K 15 = t ( t + 1 4 ( x 2 + y 2 + z 2 ) ) t + x x + y y + z z ,
which are CKVs of Bianchi I spacetime and the conformal factors are σ = 1 / ( 2 t ) for K 11 , σ = x / ( 2 t ) for K 12 , σ = y / ( 2 t ) for K 13 , σ = z / ( 2 t ) for K 14 , and σ = [ 3 t + ( x 2 + y 2 + z 2 ) / 4 ] / t for K 15 . Here, the gauge vector fields for the Noether symmetries obtained in (92)–(102) are
f 8 4 = 1 2 t ψ 2 x , f 9 4 = 1 2 t ψ 2 y , f 10 4 = 1 2 t ψ 2 z , f 11 4 = 1 8 ψ 2 t , f 12 4 = 1 8 ψ 2 ( x t + 2 x ) , f 13 4 = 1 8 ψ 2 ( y t + 2 y ) , f 14 4 = 1 8 ψ 2 ( z t + 2 z ) , f 15 4 = 1 8 ψ 2 [ 3 t 1 4 ( x 2 + y 2 + z 2 ) ] t + x x + y y + z z , f 16 , 17 4 = t h 3 ( t ) e ± b 1 + b 2 y sin ( β z ) ψ 1 2 b 3 t coth ( 2 b 3 t ) , ± b 1 , b 2 , β cot ( β z ) , f 18 , 19 4 = t h 3 ( t ) e ± b 1 b 2 y sin ( β z ) ψ 1 2 b 3 t coth ( 2 b 3 t ) , ± b 1 , b 2 , β cot ( β z ) , f 20 , 21 4 = t h 3 ( t ) e ± b 1 + b 2 y cos ( β z ) ψ 1 2 b 3 t coth ( 2 b 3 t ) , ± b 1 , b 2 , β tan ( β z ) , f 22 , 23 4 = t h 3 ( t ) e ± b 1 b 2 y cos ( β z ) ψ 1 2 b 3 t coth ( 2 b 3 t ) , ± b 1 , b 2 , β tan ( β z ) , f 24 , 25 4 = t h 4 ( t ) e ± b 1 + b 2 y sin ( β z ) ψ 1 2 b 3 t coth ( 2 b 3 t ) , ± b 1 , b 2 , β cot ( β z ) , f 26 , 27 4 = t h 4 ( t ) e ± b 1 b 2 y sin ( β z ) ψ 1 2 b 3 t coth ( 2 b 3 t ) , ± b 1 , b 2 , β cot ( β z ) , f 28 , 29 4 = t h 4 ( t ) e ± b 1 + b 2 y cos ( β z ) ψ 1 2 b 3 t coth ( 2 b 3 t ) , ± b 1 , b 2 , β tan ( β z ) , f 30 , 31 4 = t h 4 ( t ) e ± b 1 b 2 y cos ( β z ) ψ 1 2 b 3 t coth ( 2 b 3 t ) , ± b 1 , b 2 , β tan ( β z ) .
Then, the corresponding conserved flow vectors for the Noether symmetries given in (92)–(102) yield
T 8 4 = T 8 2 with q = 1 / 2 , T 9 4 = T 9 2 with q = 1 / 2 , T 10 4 = T 10 2 with q = 1 / 2 , T 11 4 = 1 2 t 3 / 2 W K 11 t ( ψ t + ψ 2 t ) T 0 + f 11 4 , T 12 4 = 1 2 t 3 / 2 W K 12 t ( x ( ψ t + ψ 2 t ) + 2 ψ x ) T 0 + f 12 4 , T 13 4 = 1 2 t 3 / 2 W K 13 t ( y ( ψ t + ψ 2 t ) + 2 ψ y ) T 0 + f 13 4 , T 14 4 = 1 2 t 3 / 2 W K 14 t ( z ( ψ t + ψ 2 t ) + 2 ψ z ) T 0 + f 14 4 , T 15 4 = 1 2 t 3 / 2 W K 15 t [ Q 4 + 1 2 t ( 3 t + 1 4 ( x 2 + y 2 + z 2 ) ) ψ ] T 0 + f 15 4 , T 16 , 17 4 = h 3 ( t ) e ± b 1 x + b 2 y sin ( β z ) T 0 + f 16 , 17 4 , T 18 , 19 4 = h 3 ( t ) e ± b 1 x b 2 y sin ( β z ) T 0 + f 18 , 19 4 T 20 , 21 4 = h 3 ( t ) e ± b 1 x + b 2 y cos ( β z ) T 0 + f 20 , 21 4 , T 22 , 23 4 = h 3 ( t ) e ± b 1 x b 2 y cos ( β z ) T 0 + f 22 , 23 4 T 24 , 25 4 = h 4 ( t ) e ± b 1 x + b 2 y sin ( β z ) T 0 + f 24 , 25 4 , T 26 , 27 4 = h 4 ( t ) e ± b 1 x b 2 y sin ( β z ) T 0 + f 26 , 27 4 T 28 , 29 4 = h 4 ( t ) e ± b 1 x + b 2 y cos ( β z ) T 0 + f 28 , 29 4 , T 30 , 31 4 = h 4 ( t ) e ± b 1 x b 2 y cos ( β z ) T 0 + f 30 , 31 4
where T 0 and W are given in (32) through F = 0 , A = B = C = t 1 / 2 , and Q 4 is defined by
Q 4 = t + 1 4 ( x 2 + y 2 + z 2 ) ψ t + x ψ x + y ψ y + z ψ z .
  • Subcase (ii.3). If A ( t ) = B ( t ) = C ( t ) = t q through F = ( b ψ + c ) 4 (see Ref. [11] for selection of this function to study the Noether symmetries in the flat FLRW spacetime), where b and c are constants, then we find that there are ten Noether symmetries:
X 1 6 = K 1 , X 2 6 = K 2 , X 3 6 = K 3 , X 4 6 = K 4 + 1 b ( b ψ + c ) ψ ,
X 5 6 = K 5 , X 6 6 = K 6 , X 7 6 = K 7 ,
X 8 6 = K 8 1 + 2 b ( q 1 ) ( b ψ + c ) x ψ , with f 8 6 = 1 b ( q 1 ) t q ψ ( b ψ + c ) x ,
X 9 6 = K 9 1 + 2 b ( q 1 ) ( b ψ + c ) y ψ , with f 9 6 = 1 b ( q 1 ) t q ψ ( b ψ + c ) y ,
X 10 6 = K 10 1 + 2 b ( q 1 ) ( b ψ + c ) z ψ , with f 10 6 = 1 b ( q 1 ) t q ψ ( b ψ + c ) z ,
where K 8 1 , K 9 1 , and K 10 1 are the same as given in (70), (71), and (72), respectively. Then, the conserved vectors of X 4 6 , , X 10 6 are
T 4 6 = 1 2 t 3 q W K 4 + ( ψ + c b Q ) T 0 , T 5 6 = T 5 2 , T 6 6 = T 6 2 , T 7 6 = T 7 2 , T 8 6 = 1 2 t 3 q W K 8 1 + 2 b ( q 1 ) x ( b ψ + c ) Q 1 T 0 + f 8 6 , T 9 6 = 1 2 t 3 q W K 9 1 + 2 b ( q 1 ) y ( b ψ + c ) Q 2 T 0 + f 9 6 , T 10 6 = 1 2 t 3 q W K 10 1 + 2 b ( q 1 ) z ( b ψ + c ) Q 3 T 0 + f 10 6 ,
where W = ψ t 2 + t 2 q ( ψ x 2 + ψ y 2 + ψ z 2 ) 2 ( b ψ + c ) 4 , Q = t ψ t + ( q 1 ) ( x ψ x + y ψ y + z ψ z ) , Q i ’s ( i = 1 , 2 , 3 ) are the same as in (76)–(78), and q 0 , 1 , 1 / 2 . For certain values of q, one can obtain either ten or fifteen Noether symmetry generators. Now, we will provide some examples illustrating this situation below.
  • (ii.3.1): For q = 1 , there are again ten Noether symmetries, such that X 1 7 = K 1 , X 2 7 = K 2 , X 3 7 = K 3 and
X 4 7 = t t + 1 b ( b ψ + c ) ψ , X 5 7 = K 5 , X 6 7 = K 6 , X 7 7 = K 7 ,
X 8 7 = K 8 2 x b ( b ψ + c ) ψ with f 8 7 = t 2 b ψ ( b ψ + c ) x ,
X 9 7 = K 9 2 y b ( b ψ + c ) ψ with f 9 7 = t 2 b ψ ( b ψ + c ) y ,
X 10 7 = K 10 2 z b ( b ψ + c ) ψ with f 10 7 = t 2 b ψ ( b ψ + c ) z ,
where K 8 2 , K 9 2 , K 10 2 are the same vector fields defined in (89). The conserved flow vectors for the above Noether symmetries are
T 5 7 = T 5 3 , T 6 7 = T 6 3 , T 7 7 = T 7 3 , T 8 7 = 1 2 t 3 W K 8 2 [ x ( t ψ t + 1 b ( b ψ + c ) ) + ln t ψ x ] T 0 + f 8 7 , T 9 7 = 1 2 t 3 W K 9 2 [ y ( t ψ t + 1 b ( b ψ + c ) ) + ln t ψ y ] T 0 + f 9 7 , T 10 7 = 1 2 t 3 W K 10 2 [ z ( t ψ t + 1 b ( b ψ + c ) ) + ln t ψ z ] T 0 + f 10 7 ,
where T 5 3 , T 6 3 , T 7 3 are of the same form as obtained in (91), and W = ψ t 2 + t 2 ( ψ x 2 + ψ y 2 + ψ z 2 ) ( b ψ + c ) 4 .
  • (ii.3.2): For q = 1 / 2 , we find that the number of Noether symmetries is fifteen, which are the KVs K 1 , K 2 , K 3 , K 5 , K 6 , K 7 and
X 4 8 = K 4 + 1 b ( b ψ + c ) ψ , X 8 8 = K 8 1 x b ( b ψ + c ) ψ ,
X 9 8 = K 9 1 y b ( b ψ + c ) ψ , X 10 8 = K 10 1 z b ( b ψ + c ) ψ ,
X 11 8 = K 11 ( b ψ + c ) 2 b t ψ , X 12 8 = K 12 x ( b ψ + c ) 2 b t ψ , X 13 8 = K 13 y ( b ψ + c ) 2 t ψ ,
X 14 8 = K 14 z ( b ψ + c ) 2 b t ψ , X 15 8 = K 15 1 2 b t [ 3 t + 1 4 ( x 2 + y 2 + z 2 ) ] ( b ψ + c ) ψ ,
where K 11 , , K 15 are of the form given in Equations (103) and (104). Then, one can obtain the corresponding non-zero gauge vectors, such that
f 8 8 = t 2 b ψ ( b ψ + 2 c ) x , f 9 8 = t 2 b ψ ( b ψ + 2 c ) y , f 10 8 = t 2 b ψ ( b ψ + 2 c ) z , f 11 8 = 1 8 b ψ ( b ψ + 2 c ) x , f 12 8 = 1 8 b ψ ( b ψ + 2 c ) ( x t + 2 x ) , f 13 8 = 1 8 b ψ ( b ψ + 2 c ) ( y t + 2 y ) , f 14 8 = 1 8 b ψ ( b ψ + 2 c ) ( z t + 2 z ) , f 15 8 = 1 4 b ψ ( b ψ + 2 c ) 3 t + 1 4 ( x 2 + y 2 + z 2 ) + x x + y y + z z .
The conserved vectors for the Noether symmetries obtained in (119)–(122) yield
T 4 8 = 1 2 t 3 / 2 W K 4 + ( ψ + c b Q ) T 0 , T 8 8 = 1 2 t 3 / 2 W K 8 1 ( Q 1 + x b ( b ψ + c ) ) T 0 + f 8 8 , T 9 8 = 1 2 t 3 / 2 W K 9 1 ( Q 2 + y b ( b ψ + c ) ) T 0 + f 19 8 , T 10 8 = 1 2 t 3 / 2 W K 10 1 ( Q 3 + z b ( b ψ + c ) ) T 0 + f 10 8 , T 11 4 = 1 2 t 3 / 2 W K 11 t ( ψ t + ψ 2 t ) T 0 + f 11 4 , T 12 8 = 1 2 t 3 / 2 W K 12 t [ x ( ψ t + 1 2 b t ( b ψ + c ) ) + 2 ψ x ] T 0 + f 12 8 , T 13 8 = 1 2 t 3 / 2 W K 13 t [ y ( ψ t + 1 2 b t ( b ψ + c ) ) + 2 ψ y ] T 0 + f 13 8 , T 14 8 = 1 2 t 3 / 2 W K 14 t [ z ( ψ t + 1 2 b t ( b ψ + c ) ) + 2 ψ z ] T 0 + f 14 8 , T 15 8 = 1 2 t 3 / 2 W K 15 t [ Q 4 + 1 2 b t ( 3 t + 1 4 ( x 2 + y 2 + z 2 ) ) ( b ψ + c ) ] T 0 + f 15 8 ,
where Q 0 , Q 1 , Q 2 , Q 3 , and Q 4 are given in Equations (49), (76), (77), (78), and (107), respectively.

4. The Field Equations

In this section, we briefly examine the gravitational field equations of the Bianchi I spacetime when influenced by an imperfect fluid source. The general expression for the energy-momentum tensor of an imperfect fluid is given by
T i j = ( ρ + P ) u i u j + P g i j + ( P x P ) χ i χ j + ( P y P ) y i y j + ( P z P ) z i z j ,
where ρ is the energy density; and P x , P y , and P z are the anisotropic pressures in the x , y , and z directions, respectively. P represents the isotropic pressure and is defined as P = ( P x + P y + P z ) / 3 . The energy density ρ is measured using a comoving observer with the timelike unit four-velocity u i , satisfying the normalization condition u i u i = 1 . Additionally, χ i , y i , z i are spacelike vectors with orthonormality conditions: χ i χ i = 1 , y i y i = 1 , z i z i = 1 , χ i y i = 0 , χ i z i = 0 , y i z i = 0 , χ i u i = 0 , y i u i = 0 and z i u i = 0 . When P x = P y = P z = P , the energy-momentum tensor T i j for the imperfect fluid, as given in (125), reduces to that of a perfect fluid. The pressure is treated as a vectorial quantity, consistent with the nature of anisotropic fluids. The equation of state (EoS) parameter for the imperfect fluid in the Bianchi I spacetime may be determined separately along each spatial axis, denoted as w i = P i / ρ , where P i = { P x , P y , P z } . The average EoS parameter for the imperfect fluid in the Bianchi I spacetime is defined as w = P / ρ = ( w x + w y + w z ) / 3 . A well-known example of dark energy occurs within the range w < 1 3 , indicating an accelerating phase for the universe. This form of dark energy is termed quintessence when 1 3 > w > 1 , and labeled a phantom fluid if w < 1 .
In the geometry described by the Bianchi I metric given in (2), the timelike four-velocity in comoving coordinates is u i = δ 0 i . From the properties of the spacelike vectors, it follows that χ i = A 1 δ 1 i , y i = B 1 δ 2 i , and z i = C 1 δ 3 i . Subsequently, the Einstein field equations G i j = 8 π T i j in natural units ( G = 1 and c = 1 ) can be expressed as follows:
A ˙ B ˙ A B + A ˙ C ˙ A C + B ˙ C ˙ B C = 8 π ρ ,
A ¨ A + B ¨ B + A ˙ B ˙ A B = 8 π P x ,
A ¨ A + C ¨ C + A ˙ C ˙ A C = 8 π P y ,
B ¨ B + C ¨ C + B ˙ C ˙ B C = 8 π P z ,
where G i j = R i j 1 2 R g i j is the Einstein tensor, R i j is the Ricci tensor, and R is the Ricci scalar, which has the form
R = 2 A ¨ A + B ¨ B + C ¨ C + A ˙ B ˙ A B + A ˙ C ˙ A C + B ˙ C ˙ B C .
A curvature scalar characterizes the spacetime curvature. One of the important curvature scalars is the Kretschmann scalar, defined by K = R i j k l R i j k l , which is a quadratic scalar invariant of the Riemann tensor R i j k l . It is crucial for measuring curvature in a vacuum. This curvature scalar also characterizes the spacetime curvature of a realistic rotating black hole, allowing us to mathematically perceive the black hole. Additionally, we consider an important curvature scalar, the so-called Gauss–Bonnet (GB) invariant G , defined as G = R 2 4 R i j R i j + R i j k l R i j k l . Apart from emerging in the context of defining quantum fields in curved spacetimes, the GB invariant G can encapsulate all the curvature information stemming from the Riemann tensor in dynamical equations. For the Bianchi I spacetime, both the K and G take the following forms
K = 4 A ¨ 2 A 2 + B ¨ 2 B 2 + C ¨ 2 C 2 + ( A ˙ B ˙ ) 2 ( A B ) 2 + ( A ˙ C ˙ ) 2 ( A C ) 2 + ( B ˙ C ˙ ) 2 ( B C ) 2 ,
G = 8 A ¨ A + B ¨ B + C ¨ C .
In the previous section, we investigated the power law form of metric functions, such as A ( t ) = t L , B ( t ) = t p , and C ( t ) = t q in specific cases. With this form of the Bianchi I metric, analysis of the field Equations (126) to (128) reveals the expressions for the physical variables ρ , P x , P y and P z to be
ρ = ( L p + L q + p q ) 8 π t 2 , P x = ( p 2 + p q + q 2 p q ) 8 π t 2 ,
P y = ( L 2 + L q + q 2 L q ) 8 π t 2 , P z = ( L 2 + L p + p 2 L p ) 8 π t 2 ,
As a result, the directional EoS parameters w x , w y , and w z along the x , y , and z axes, respectively, are determined as follows:
w x = p + q p 2 p q q 2 L p + L q + p q , w y = L + q L 2 L q q 2 L p + L q + p q , w z = L + p L 2 L p p 2 L p + L q + p q .
Subsequently, the average EoS parameter w is expressed as
w = 1 + 2 3 ( L + p + q ) ( 1 L p q ) L p + L q + p q ,
where L p + L q + p q 0 . Finally, the Ricci, Kretschmann, and the GB scalars for the power law form of metric coefficients are obtained as follows:
R = 4 t 2 L 2 + p 2 + q 2 + L p + L q + p q L p q , G = 8 L p q ( L + p + q 3 ) t 4 ,
K = 8 t 4 L 4 + p 4 + q 4 2 ( L 3 + p 3 + q 3 ) + L 2 p 2 + L 2 q 2 + p 2 q 2 + L 2 + p 2 + q 2 .
If L p + L q + p q = 0 , both p and q should be zero, and there is no need for vanishing of the power L. This means that the energy density ρ is zero, while the anisotropic pressures in x , y and z directions are P x = 0 , P y = L ( 1 L ) / ( 8 π t 2 ) and P z = P y , respectively. In theoretical physics and cosmology, scenarios where the energy density vanishes (or is extremely low) while the pressure is different from zero are less common. However, there are still certain contexts where this can occur, often involving exotic fields or conditions. In certain models of the early universe, there is a concept known as the false vacuum. The vacuum state of a field can be thought of as the state of lowest energy. In a false vacuum, the energy density is very close to zero, but the field is in a metastable state rather than the true vacuum state.
When considering a perfect fluid source ( P x = P y = P z = P ), Equations (133) to (134) imply the following constraint relations:
L ( L 1 ) p ( p 1 ) + q ( L p ) = 0 , L ( L 1 ) q ( q 1 ) + p ( L q ) = 0 .
For cases where L = p = q , as seen in cases (ii.2) and (ii.3), these constraint equations are precisely satisfied. Consequently, the physical quantities ρ and P take the forms
ρ = 3 q 2 8 π t 2 , P = q ( 2 3 q ) 8 π t 2 ,
yielding the EoS parameter w = ( 2 3 q ) / ( 3 q ) where q 0 . So, dark energy occurs when q > 1 , which mentions that a well-known example of dark energy ( w = 1 ), the cosmological constant, is not possible for the power law Bianchi I spacetime with a perfect fluid. In addition, we have dust fluid if q = 1 / 3 , and a stiff fluid if q = 1 / 2 . Additionally, the curvature scalars R , K and G for cases (ii.1) and (ii.3) are given by
R = 6 q ( 2 q 1 ) t 2 , K = 12 q 2 ( 2 q 2 2 q + 1 ) t 4 , G = 24 q 3 ( q 1 ) t 4 .
Moreover, when L = p = q = 2 / 3 , the spacetime corresponds to the Einstein–de Sitter model, a solution of field Equations (126) to (128) with P = 0 and ρ = 1 / ( 6 π t 2 ) , i.e., w = 0 (representing a dust fluid).
Kasner spacetimes, non-trivial Bianchi I solutions satisfying Einstein’s vacuum equations ( R i j = 0 , or ρ = 0 and P = 0 ), are characterized using the Kasner metric. In this metric, the three constant parameters L , p , and q are termed the Kasner components. These components adhere to the following relations:
L + p + q = 1 , L 2 + p 2 + q 2 = 1 ,
resulting in L p + L q + p q = 0 . Kasner spacetimes exhibit spatial anisotropy, potentially expanding along one direction while contracting along another. In the Kasner scenario, the three parameters cannot all be equal. Thus, at least one of L , p , q must satisfy 0 L , p , q 1 . If any of the parameters L , p , or q equals one, the other two must be zero, resulting in a flat Kasner spacetime. Conversely, if L , p , q 1 are all not equal to one, indicating all are non-zero, the Kasner spacetime is non-flat [17].

5. Discussions and Concluding Remarks

Though the standard procedure for determining Lie and Noether symmetry generators can be cumbersome, it is feasible to consider reduction or obtaining conservation laws through Noether’s theorem. This paper delved into the geometric nature of the Klein–Gordon/wave equation within the framework of the Lagrangian for the Bianchi I spacetime. We demonstrated that computing Noether symmetries of the first-order Lagrangian (14) for the Klein–Gordon Equation (15) is simplified to solving a set of differential conditions outlined by (17). For the Bianchi I spacetime, employing this method to find Noether symmetries yields solutions encompassing scenarios where F ( t , x , y , z , ψ ) = U 0 (let us say, U 0 = 0 for convenience), reducing Equation (15) to the wave equation. Conversely, when F ( t , x , y , z , ψ ) const . , Equation (15) represents the Klein–Gordon equation. Exploring various functional forms of metric coefficients, we derived exact solutions of the Noether symmetry Equations (17) for the Bianchi I spacetime.
In case (i), for F = U 0 + U 1 ψ 2 + 1 2 U 2 2 ψ 2 , we obtained the Noether symmetry generators and gauge vector components for arbitrary metric coefficients. Additionally, we constructed conserved vector fields corresponding to these Noether symmetries. Selected solutions of Equation (25) are presented in Table 1 for specific trigonometric or hyperbolic metric functions other than power-law forms, involving the unknown function Y ( t ) in the components of Noether symmetry generators.
In case (ii), utilizing power-law forms of metric coefficients, we derived Noether symmetries in terms of constant powers L , p , and q. In Table 2, a set of nontrivial solutions of Equation (43), which solve Y ( t ) akin to Equation (25), is provided for a specific class of Bianchi I metrics. In subcases of (ii), particularly (ii.1), (ii.2), and (ii.3), we discovered special values of powers L , p , and q, resulting in various dimensions of Noether symmetry groups. Additionally, Einstein field equations incorporating an imperfect fluid were derived for arbitrary metric functions within the context of Bianchi I spacetime and particularly studied for the power-law form of that geometry.
Let us briefly discuss a comparison and contrast between this study and some related references. In [15], the study focused on the Lie symmetries applied to the Klein–Gordon and wave equations within Bianchi I spacetime in terms of the potential function V ( x i ) . They considered the function F as F = V ( x i ) ψ 2 / 2 , which gives G ( x i , ψ ) = V ( x i ) ψ , specifically for the Klein–Gordon equation. In this context, the investigation in [14] revealed that for F = V ( x i ) ψ 2 / 2 , the Noether symmetry component of the scalar field ψ is Φ = 2 σ ψ + b ( x k ) , with b ( x k ) being a solution of (3). Then, it follows from the Noether symmetry condition (13) that the equation A , i i = 0 and the relationship
V , i ξ i + 2 σ V 2 σ = 0 ,
were established, incorporating the potential function V ( x i ) and the conformal factor σ ( x i ) . In contrast, our study directly delved into the Noether symmetry Equations (17) employing an arbitrary function F ( x i , ψ ) , bypassing the need to solve condition (143) to find potentials V ( x i ) concerning point symmetries generated by the spacetime symmetries of Bianchi I metric (2). We obtained solutions for these equations, including components of Noether symmetry generators, employing specific solvable forms of the function F within Bianchi I spacetime. Consequently, we not only acquired conserved Noether currents but also derived non-trivial gauge vectors for each Noether symmetry. Moreover, in the paper [30], the author presented solutions for the classification problem concerning the Klein–Gordon equation with a non-constant potential function V ( x i ) within the framework of Bianchi I, Bianchi III, and Bianchi V spacetimes. This investigation showcased that the Klein–Gordon equation admits non-trivial Lie symmetries derived from conformal algebra. However, the metric tensor provided for Bianchi I spacetime in that reference, given as
g i j = diag C 2 ( t ) , e 2 c t C 2 ( t ) , e 2 α 1 c t C 2 ( t ) , C 2 ( t ) ,
is less general compared to the one employed in our study, described by (2). This study investigated the Klein–Gordon equation for Bianchi I spacetime, employing the Noether symmetry approach as a means to derive conservation laws. The novelty of this research lies in utilizing the Noether symmetry method to discern the symmetries inherent in the Lagrangian, thereby establishing the associated conservation laws for the Klein–Gordon equation of Bianchi I spacetime.
It is widely recognized that Lie point symmetries serve as a tool for reducing the Klein–Gordon/wave equation, thereby identifying corresponding invariant solutions. The reduction based on the Lie invariants of the symmetry vector leads to a solution of the Klein–Gordon/wave equation, coupled with a constraint equation related to an arbitrary integration function within the solution. In the paper [22], invariant solutions of the wave equation on Bianchi I spacetime using the Lie symmetry method were investigated. In our study, we used the power-law form of metric coefficients outlined in this reference. However, we extended the investigation by exploring both the Noether symmetries of the wave and Klein–Gordon equations, obtaining conserved quantities for each Noether symmetry. Considering that every Noether symmetry should be a Lie symmetry, we aimed to utilize Noether symmetries for reducing both the wave equation and Klein–Gordon equation. Meanwhile, the common Lie symmetry vector for Bianchi I spacetime is ψ ψ [15], which is not a Noether symmetry vector. Based on this observation, we inferred that the scalar field ψ , representing the solution of the Klein–Gordon equation, should be equal to the component Φ of the Noether symmetry generator X ; that is, ψ = Φ . Further exploration of this relationship will be the focus of another paper.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The author would like to thank unknown referees for their valuable comments and suggestions.

Conflicts of Interest

The author declares no conflicts of interest.

References

  1. Katzin, G.H.; Levine, J.; Davis, W.R. Curvature Collineations: A Fundamental Symmetry Property of the Space-Times of General Relativity Defined by the Vanishing Lie Derivative of the Riemann Curvature Tensor. J. Math. Phys. 1969, 10, 617. [Google Scholar]
  2. Stephani, H. Differential Equations: Their Solution Using Symmetries; Cambridge University Press: Cambridge, UK, 1989. [Google Scholar]
  3. Ibragimov, N.H. CRC Handbook of Lie Group Analysis of Differential Equations: Symmetries, Exact Solutions and Conservation Laws; CRC Press: Boca Raton, FL, USA, 1994. [Google Scholar]
  4. Capozziello, S.; Lambiase, G. Selection rules in minisuperspace quantum cosmology. Gen. Relativ. Gravit. 2000, 32, 673. [Google Scholar]
  5. Camci, U.; Kucukakca, Y. Noether symmetries of Bianchi I, Bianchi III, and Kantowski-Sachs spacetimes in scalar-coupled gravity theories. Phys. Rev. D 2007, 76, 084023. [Google Scholar] [CrossRef]
  6. Kucukakca, Y.; Camci, U.; Semiz, I. LRS Bianchi type I universes exhibiting Noether symmetry in the scalar-tensor Brans-Dicke theory. Gen. Relativ. Gravit. 2012, 44, 1893. [Google Scholar]
  7. Feroze, T.; Mahomed, F.M.; Qadir, A. The connection between isometries and symmetries of geodesic equations of underlying spaces. Nonlinear Dyn. 2006, 45, 65. [Google Scholar] [CrossRef]
  8. Tsamparlis, T.; Paliathanasis, A. Lie and Noether symmetries of geodesic equations and collineations. Gen. Relativ. Gravit. 2010, 42, 2957. [Google Scholar] [CrossRef]
  9. Tsamparlis, T.; Paliathanasis, A. The geometric nature of Lie and Noether symmetries. Gen. Relativ. Gravit. 2011, 43, 1861. [Google Scholar]
  10. Kucukakca, K.; Camci, U. Noether gauge symmetry for f(R) gravity in Palatini formalism. Astrophys. Space Sci. 2012, 338, 211. [Google Scholar] [CrossRef]
  11. Camci, U.; Jamal, S.; Kara, A.H. Invariances and Conservation Laws Based on Some FRW Universes. Int. J. Theor. Phys. 2014, 53, 1483. [Google Scholar] [CrossRef]
  12. Camci, U.; Yildirim, A. Lie and Noether symmetries in some classes of pp-wave spacetimes. Phys. Scr. 2014, 89, 084003. [Google Scholar]
  13. Camci, U.; Yildirim, A. Noether gauge symmetry classes for pp-wave spacetimes. Int. J. Geom. Meth. Mod. Phys. 2015, 12, 1550120. [Google Scholar] [CrossRef]
  14. Paliathanasis, A.; Tsamparlis, M. The geometric orijin of Lie point symmetries of the Schrödinger and the Klein-Gordon equations. Int. J. Geom. Methods Mod. Phys. 2014, 11, 1450037. [Google Scholar] [CrossRef]
  15. Paliathanasis, A.; Tsamparlis, M.; Mustafa, M.T. Symmetry analysis of the Klein-Gordon equation in Bianchi I spacetimes. Int. J. Geom. Meth. Mod. Phys. 2015, 12, 1550033. [Google Scholar]
  16. Ellis, G.F.R.; Maartens, R.; MacCallum, M.A.H. Relativistic Cosmology; Cambridge University Press: Cambridge, UK, 2012. [Google Scholar]
  17. Stephani, H.; Kramer, D.; MacCallum, M.A.H.; Hoenselaers, C.; Herlt, E. Exact Solutions of Einstein Field Equations; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar]
  18. Sharif, S.; Jhangeer, A. Conserved quantities for (1+ 2)-dimensional non-linear wave equation on curved surfaces. Commun. Nonlinear Sci. Numer. Simulat. 2013, 18, 1684. [Google Scholar] [CrossRef]
  19. Jhangeer, A.; Sharif, S. Conserved quantities and group classification of wave equation on hyperbolic space. Commun. Nonlinear Sci. Numer. Simulat. 2013, 18, 236. [Google Scholar]
  20. Jamal, S.; Shabbir, G. Noether symmetries of vacuum classes of pp-waves and the wave equation. Int. J. Geom. Meth. Mod. Phys. 2016, 13, 1650109. [Google Scholar] [CrossRef]
  21. Jamal, S.; Kara, A.H.; Narain, R. Wave Equations in Bianchi Space-Times. J. App. Math. 2012, 2012, 765361. [Google Scholar]
  22. Abdulwahhab, M.A.; Jhangeer, A. Symmetries and generalized higher order conserved vectors of the wave equation on Bianchi I spacetime. Int. J. Geom. Meth. Mod. Phys. 2017, 14, 1750028. [Google Scholar] [CrossRef]
  23. Abdulwahhab, M.A. Conservation laws and invariant solutions of the wave equation on Bianchi I space-time: A comprehensive analysis. Optik 2021, 3231, 166364. [Google Scholar] [CrossRef]
  24. Olver, P.J. Applications of Lie Groups to Differential Equations, 2nd ed.; Springer: New York, NY, USA, 1993. [Google Scholar]
  25. Anco, S.C.; Bluman, G.W. Direct construction method for conservation laws of partial differential equations. Part I. Examples of conservation law classifications. Eur. J. Appl. Math. 2002, 13, 545. [Google Scholar]
  26. Jamal, S.; Kara, A.H.; Bokhari, A.H. Symmetries, conservation laws, reduction, and exact solutions for the Klein-Gordon equation in de Sitter spacetimes. Can. J. Phys. 2012, 90, 667. [Google Scholar]
  27. Paliathanasis, A.; Tsamparlis, M.; Mustafa, M.T. Classification of the Lie and Noether point symmetries for the Wave and the Klein-Gordon equations in pp-wave spacetimes. Commun. Nonlinear Sci. Numer. Simul. 2018, 55, 68. [Google Scholar] [CrossRef]
  28. Jamal, S.; Kara, A.H.; Bokhari, A.H. Symmetries, conservation laws and reduction of wave and Gordon-type equations on Riemannian manifolds. In Proceedings of the International Conference on Applied Mathematics and Engineering (WASET’11), Paris, France, 14–16 November 2011; World Academy of Science Engineering and Technology: Puket, Thailand, 2011. [Google Scholar]
  29. Jamal, S.; Kara, A.H.; Bokhari, A.H.; Zaman, F.D. The symmetries and conservation laws of some Gordon-type equations in Milne spacetime. Paramana J. Phys. 2013, 80, 739. [Google Scholar] [CrossRef]
  30. Paliathanasis, A. Classification of the Lie and Noether Symmetries for the Klein-Gordon Equation in Anisotropic Cosmology. Symmetry 2023, 15, 306. [Google Scholar] [CrossRef]
  31. Bozhkov, Y.; Freire, I.L. Special conformal groups of a Riemannian manifolds and Lie point symmetries of the nonlinear Poisson equation. J.Differ. Equ. 2010, 249, 872. [Google Scholar] [CrossRef]
  32. Obukhov, V.V. Algebra of Symmetry Operators for Klein-Gordon-Fock Equation. Symmetry 2021, 13, 727. [Google Scholar] [CrossRef]
  33. Obukhov, V.V. Algebra of the Symmetry Operators of the Klein-Gordon-Fock Equation for the Case When Groups of Motions G3 Act Transitively on Null Subsurfaces of Spacetime. Symmetry 2022, 14, 346. [Google Scholar] [CrossRef]
  34. Obukhov, V.V. Algebras of integrals of motion for the Hamilton-Jacobi and Klein-Gordon-Fock equations in spacetime with four-parameter groups of motions in the presence of an external electromagnetic field. J. Math. Phys. 2022, 63, 023505. [Google Scholar]
  35. Obukhov, V.V. Hamilton-Jacobi and Klein-Gordon-Fock equations for a charged test particle in space-time with simply transitive four-parameter groups of motions. J. Math. Phys. 2023, 64, 093507. [Google Scholar] [CrossRef]
  36. Qadir, A.; Camci, U. Physical Significance of Noether Symmetries. Symmetry 2022, 14, 476. [Google Scholar] [CrossRef]
Table 1. Some solutions of the Equation (25) for the metric functions given in subcases (i.2), (i.3), and (i.4). Here, it is defined that k = a 1 2 + a 2 2 + a 3 2 , k 1 = 1 + a 2 2 + a 3 2 , k 2 = k 2 + q 2 / ( 2 q ) , k 3 = k 2 + 4 / 4 , k 4 = a 2 2 + a 3 2 / q , k 5 = a 1 2 + q 2 / q , k 6 = k 2 U 2 2 / q , u 1 = 9 q 2 4 U 2 2 / ( 4 q ) , u 2 = 4 U 2 2 / 2 , u 3 = 25 4 U 2 2 and u 4 = 26 4 U 2 2 + 2 u 3 . We also defined that 1 ( t ) = cos ( 2 q t ) + 1 , 2 ( t ) = cos ( 2 q t ) 1 , and 3 ± ( t ) = ( cos 2 t ) 1 2 ( k 1 ± 1 ) . The special functions H G = H e u n G 1 , α , β , γ , δ , τ ( t ) and H C = H e u n C α , β , γ , δ , η , τ ( t ) are the Heun general function and the Heun confluent function, respectively. Furthermore, 2 F 1 n 1 , n 2 ; d ; ( t ) represents the hypergeometric function.
Table 1. Some solutions of the Equation (25) for the metric functions given in subcases (i.2), (i.3), and (i.4). Here, it is defined that k = a 1 2 + a 2 2 + a 3 2 , k 1 = 1 + a 2 2 + a 3 2 , k 2 = k 2 + q 2 / ( 2 q ) , k 3 = k 2 + 4 / 4 , k 4 = a 2 2 + a 3 2 / q , k 5 = a 1 2 + q 2 / q , k 6 = k 2 U 2 2 / q , u 1 = 9 q 2 4 U 2 2 / ( 4 q ) , u 2 = 4 U 2 2 / 2 , u 3 = 25 4 U 2 2 and u 4 = 26 4 U 2 2 + 2 u 3 . We also defined that 1 ( t ) = cos ( 2 q t ) + 1 , 2 ( t ) = cos ( 2 q t ) 1 , and 3 ± ( t ) = ( cos 2 t ) 1 2 ( k 1 ± 1 ) . The special functions H G = H e u n G 1 , α , β , γ , δ , τ ( t ) and H C = H e u n C α , β , γ , δ , η , τ ( t ) are the Heun general function and the Heun confluent function, respectively. Furthermore, 2 F 1 n 1 , n 2 ; d ; ( t ) represents the hypergeometric function.
Case L , p , q Y ( t )
( i . 2 ) L = p = q 2 ( t ) / 2 k 2 + 1 2 2 ( t ) [ b 1 2 F 1 n , n ; 1 2 ; 1 2 1 ( t ) + b 2 1 ( t ) 2 F 1 n , n ; 3 2 ; 1 2 1 ( t ) ]
where n = k 2 + 1 4 + u 1
L = 2 , p = q = 1 sin ( 2 t ) ( cos 2 t 1 ) k 3 2 ( t ) 1 ( t ) [ b 1 ( cos t ) a 1 2 2 F 1 n 1 , n 2 ; 2 d ; 1 2 1 ( t )
+ b 2 ( cos t ) a 1 2 2 F 1 n 3 , n 4 ; 2 d + ; 1 2 1 ( t ) ]
where d ± = 1 2 ( 1 ± a 1 2 ) , n 1 = k 3 + d u 2 , n 2 = k 3 + d + u 2
n 3 = k 3 + d + + u 2 , n 4 = k 3 + d + u 2
( i . 3 ) L = p = q 2 ( t ) 1 ( t ) [ b 1 [ 2 2 ( t ) ] 1 2 + k 2 2 F 1 n 1 , n 2 ; 1 + 2 k 2 ; 1 2 1 ( t )
+ b 2 [ 2 ( t ) ] 1 2 k 2 2 F 1 n 3 , n 4 ; 1 2 k 2 ; 1 2 1 ( t ) ]
where n 1 = k 2 + 3 4 u 1 , n 2 = k 2 + 3 4 + u 1 , n 3 = k 2 + 3 4 + u 1 ,
n 4 = k 2 + 3 4 u 1
L = 1 , p = q = 2 ( cos t ) a 1 2 ( t ) 1 + u 3 + u 4 4 u 4 [ b 1 3 H G 1 , α 1 , β 1 , γ 1 , δ 1 , η , cos 2 t
+ b 2 3 + H G 1 , α 2 , β 2 , γ 2 , δ 2 , η , cos 2 t ]
where λ 1 ± = 1 + 1 2 ( a 1 ± k 1 ) , δ 1 = 1 + k 1 , δ 2 = 1 k 1 , η = 1 + a 1
α 1 = 1 2 ( 1 a 1 δ 1 ) + k 1 2 u 4 k 1 ( 1 k 1 ) + 1 2 ( 1 + u 3 ) + 1 + u 3 4 u 4 ,
α 2 = 1 2 ( 1 a 1 δ 2 ) k 1 2 u 4 k 1 ( 1 k 1 ) 1 2 ( 1 + u 3 ) + 1 + u 3 4 u 4 ,
β 1 = λ 1 + + 1 4 u 4 , β 2 = λ 1 + 1 4 u 4 , γ 1 = λ 1 + + U 2 2 6 u 4 , γ 2 = λ 1 + U 2 2 6 u 4
( i . 4 ) L = p = q cosh 4 ( q t ) [ b 1 H C 0 , 1 2 , 2 , k 2 4 q 2 , 3 q 2 k 2 U 2 2 4 q 2 , sin 2 t
+ b 2 sinh ( q t ) H C 0 , 1 2 , 2 , k 2 4 q 2 , 3 q 2 k 2 U 2 2 4 q 2 , sin 2 t ]
L = 0 , p = q b 1 H C 0 , 3 2 , 1 2 , k 4 2 4 , 5 q 2 2 a 1 2 2 U 2 2 8 q 2 , cosh 2 ( q t )
+ b 2 cosh 3 ( q t ) H C 0 , 3 2 , 1 2 , k 4 2 4 , 5 q 2 2 a 1 2 2 U 2 2 8 q 2 , cosh 2 ( q t )
Table 2. Some solutions of Equation (43) for specific values of L , p , and q. Here, it is defined that ϵ 1 = a 2 2 + a 3 2 + U 2 2 , ϵ 2 = a 3 2 + U 2 2 , η 1 = 1 + 4 a 1 2 and η 2 = 1 + 16 a 1 2 . Furthermore, the special functions K M = K u m m e r M ( μ , ν , z ) and K U = K u m m e r U ( μ , ν , z ) are the Kummer functions, W M = W h i t t a k e r M ( μ , ν , z ) and W W = W h i t t a k e r W ( μ , ν , z ) are the Whittaker functions, respectively.
Table 2. Some solutions of Equation (43) for specific values of L , p , and q. Here, it is defined that ϵ 1 = a 2 2 + a 3 2 + U 2 2 , ϵ 2 = a 3 2 + U 2 2 , η 1 = 1 + 4 a 1 2 and η 2 = 1 + 16 a 1 2 . Furthermore, the special functions K M = K u m m e r M ( μ , ν , z ) and K U = K u m m e r U ( μ , ν , z ) are the Kummer functions, W M = W h i t t a k e r M ( μ , ν , z ) and W W = W h i t t a k e r W ( μ , ν , z ) are the Whittaker functions, respectively.
Case L , p , q Y ( t )
1 L = 1 2 , p = q = 0 t e ϵ 1 t b 1 K M ( μ , ν , 2 ϵ 1 t ) + b 2 K U ( μ , ν , 2 ϵ 1 t )
where μ = a 1 2 2 ϵ 1 + 3 4 , ν = 3 2 .
2 L = 1 , p = q = 0 t 2 e a 1 t 2 2 b 1 K M ( μ , ν , a 1 t 2 ) + b 2 K U ( μ , ν , a 1 t 2 )
where μ = 1 + ϵ 1 2 4 a 1 , ν = 2 .
3 L = 1 , p = 1 , q = 0 1 t b 1 W M ( μ , ν , a 2 t 2 ) + b 2 W W ( μ , ν , a 2 t 2 )
where μ = ϵ 2 2 4 a 2 , ν = 1 4 η 1 .
4 L = 1 , p = 1 2 , q = 0 t 3 4 b 1 W M ( μ , ν , 2 ϵ 2 t ) + b 2 W W ( μ , ν , 2 ϵ 2 t )
where μ = a 2 2 2 ϵ 2 , ν = 1 4 η 2 .
5 L = 1 , p = 1 , q = 3 t 2 e a 3 t 4 4 [ b 1 t α H B ( α , β , γ , δ , a 3 2 t 2 ) + b 2 t α H B ( α , β , γ , δ , a 3 2 t 2 ) ]
where α = 4 + a 1 2 , β = 0 , γ = a 2 2 2 a 3 , δ = U 2 2 2 a 3
6 L = 3 , p = 2 , q = 1 t 3 2 [ b 1 e ( a 3 t 4 a 1 ) 2 t 2 H D ( α , β , γ , δ , α t 2 4 a 1 α t 2 + 4 a 1 )
+ b 2 e ( a 3 t 4 a 1 ) 2 t 2 H D ( α , β , γ , δ , α t 2 4 a 1 α t 2 + 4 a 1 ) ]
where α = 4 a 1 a 3 , β = α a 2 2 4 a 1 + α 2 8 9 4 + U 2 2 a 1 a 3 ,
γ = α a 2 2 2 a 1 + 2 U 2 2 a 1 a 3 , δ = α a 2 2 4 a 1 α 2 8 + 9 4 + U 2 2 a 1 a 3 .
7 L = 1 , p = 3 , q = 1 1 t [ b 1 e ( a 1 t 4 a 2 ) 2 t 2 H D ( α , β , γ , δ , α t 2 4 a 2 α t 2 + 4 a 2 )
+ b 2 e ( a 1 t 4 a 2 ) 2 t 2 H D ( α , β , γ , δ , α t 2 4 a 2 α t 2 + 4 a 2 ) ]
α = 4 a 1 a 2 , β = α 2 8 a 3 2 1 + U 2 2 a 2 a 1
γ = 2 U 2 2 a 2 a 1 , δ = α 2 8 + a 3 2 + 1 + U 2 2 a 2 a 1 .
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Camci, U. Noether Symmetry Analysis of the Klein–Gordon and Wave Equations in Bianchi I Spacetime. Symmetry 2024, 16, 115. https://doi.org/10.3390/sym16010115

AMA Style

Camci U. Noether Symmetry Analysis of the Klein–Gordon and Wave Equations in Bianchi I Spacetime. Symmetry. 2024; 16(1):115. https://doi.org/10.3390/sym16010115

Chicago/Turabian Style

Camci, Ugur. 2024. "Noether Symmetry Analysis of the Klein–Gordon and Wave Equations in Bianchi I Spacetime" Symmetry 16, no. 1: 115. https://doi.org/10.3390/sym16010115

APA Style

Camci, U. (2024). Noether Symmetry Analysis of the Klein–Gordon and Wave Equations in Bianchi I Spacetime. Symmetry, 16(1), 115. https://doi.org/10.3390/sym16010115

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

Article Metrics

Back to TopTop