## Computation of Partially Invariant Solutions for the Einstein Walker Manifolds' Identifying Equations

Mehdi Nadjafikhah a, ∗ , Mehdi Jafari b

a b

- School of Mathematics, Iran University of Science and Technology, Narmak, Tehran 1684613114, Iran
- Department of Complementary Education, Payame Noor University, PO BOX 19395-3697, Tehran, Iran

## Abstract

In this paper, partially invariant solutions (PISs) method is applied in order to obtain new four-dimensional Einstein Walker manifolds. This method is based on subgroup classification for the symmetry group of partial differential equations (PDEs) and can be regarded as the generalization of the similarity reduction method. For this purpose, those cases of PISs which have the defect structure δ = 1 and are resulted from two-dimensional subalgebras are considered in the present paper. Also it is shown that the obtained PISs are distinct from the invariant solutions that obtained by similarity reduction method.

Keywords: Einstein Walker manifolds, Lie symmetry group, Optimal system of Lie subalgebras, Partially invariant solutions (PISs)

2000 MSC:

70G65, 34C14, 53C50

## 1. Introduction

The idea of analyzing the differential equations by applying the transformation groups implied a new theory: the symmetry group theory, also called Lie group analysis. This method was originated at the end of nineteenth century from the pioneering work of Sophus Lie [1]. Since that time, several books have been dedicated to this concept and its generalizations [2-5]. Classification of the group invariant solutions and reduction of the original system can be regarded as significant applications of the Lie group method in the theory of differential equations. The fact that symmetry reductions for many PDEs can not be obtained via the classical symmetry method, motivated the creation of several generalizations of the classical Lie group method for symmetry reductions. Consequently, several alternative reduction methods have been proposed, going beyond Lie's classical procedure and providing further solutions. Partially invariant solutions (PISs) method is one of these procedures. This algorithmic method is a powerful tool for the reduction of PDEs and is based on classifying the subgroups of the symmetry group. The notion of PISs was first developed by Ovsiannikov [5] and can be considered as the extension of invariant solutions. The algorithm of constructing PISs is similar to that of invariant solutions. Indeed, obtaining the invariant solutions by applying the PISs method is easier than similarity reduction method whenever we deal with low-dimensional groups. One significant concept which appears while constructing PISs is the defect structure. This quantity is determined by the dimension of orbits and is denoted by δ .

In this paper, the PISs method will be applied in order to construct some new classes of four-dimensional Einstein Walker manifolds. Those Manifolds which admit null parallel distributions are called Walker manifolds . A Walker manifold is called Einstein Walker manifold if its Ricci tensor is a scaler multiple of the metric at each point. Four-dimensional Einstein Walker manifolds form the underling structure of many geometric and physical models such as; hh-space in general relativity, pp-wave model and other areas [6-12].

∗ Corresponding author. Tel.:+98 21 73225426; Fax: +98 21 73228426 Email addresses: m nadjafikhah@iust.ac.ir (Mehdi Nadjafikhah), m.jafari@phd.pnu.ac.ir (Mehdi Jafari)

Preprint submitted to Commun Nonlinear Sci Numer Simulat

The general form of the metric tensor of four-dimensional walker manifolds is expressed as follows [13]:

$$g _ { a, b, c } & \coloneqq 2 ( d x \circ d y$$

where a , b and c are smooth functions with respect to ( x, t, y, z ). Let M a,b,c := ( O , g a,b,c ), where O be an open subset of R 4 . We can see that M a,b,c is Einstein if and only if the functions a , b and c verify the following system of PDEs ([13], page 81).

$$\begin{array} { c } a _ { 1 1 } - b _ { 2 2 } =$$

where the index 1, 2, 3 and 4 for functions a , b and c represent the derivatives of these functions with respect to x , t , y and z , respectively. The system (2) is hard to handle, so we consider a spacial case in this paper; where a , b and c only depend on x and t . Therefore the following system must be solved.

$$\begin{array} {$$

In system (3), x and t are independent and a , b and c are dependent variables. It is worthwhile to say that some other special cases have been considered in many references which yield some results about the structures admitted by these manifolds [14-16].

In [17], we have comprehensively analyzed the problem of symmetries of the system (3). By applying the basic Lie symmetry method, we have obtained the classical Lie point symmetry operators of the system (3) and proved the following result (refer to [17] for more details):

Corollary 1.1. The Lie group of point symmetries of the PDE system (3) has a seven-dimensional Lie subalgebra generated by the following vector fields:

$$\begin{array} { c c c } X _ { 1 } = \partial _ { x }$$

$$( \partial _ { x } \equiv \frac { \partial } { \partial x },$$

Mainly, we have constructed an optimal system of one-dimensional subalgebras in [17] which provides the preliminary classification of group invariant solutions for the system (3). Also, we have obtained the corresponding invariant solutions of this system via the method of similarity reduction.

It is worth considering that some of the partially invariant solutions are not invariant with respect to the subalgebras of lower dimensions. They are known as non-reducible PISs. By determining non-reducible PISs for the system (3), we can obtain some new four-dimensional Einstein Walker manifolds.

This paper is organized as follows: In section 2, we recall the general procedure of determining PISs. In section 3, we classify the two-dimensional subalgebras of the symmetry Lie algebra and construct an optimal system. Section 4 is devoted to the computation of some of the PISs associated to the system (3). In section 5, the relation between PISs and invariant solutions is investigated and the condition for obtaining the non-reducible PISs is stated. Meanwhile, some non-reducible PISs for the system (3) are presented. Some concluding remarks are declared at the end of the paper.

## 2. Partially invariant solutions method

In this section, we recall the general procedure for determining PISs for an arbitrary system of PDEs. First, we present a brief review of the concept of PISs [18, 5]. Consider a system of PDEs of n th order with

p independent variables ( x = ( x i ) ∈ X,i = 1 , ..., p ) and q dependent variables ( u = ( u j ) ∈ U, j = 1 , ..., q ) such as:

$$\Delta = \Delta _ { \mu } ( x, u ^ { ( n ) } ) = 0, \quad \ \mu = 1, \dots, r.$$

Let G be a local symmetry group of the above system which acts on the total space X × U with r -dimensional orbits. If u = f ( x ) is a solution of the system (5) with graph Γ f , then the orbit space of Γ f can be defined as follows:

$$G \Gamma _ { f } = \{ g. ( x, u ) \ | \ ( x, u ) \in \Gamma _ { f }, g \in G \}$$

which is the union of the orbits of the Γ f -elements.

The defect structure of the solution u = f ( x ) with respect to the group G is computed by the matrix of generators' characteristics and is defined as:

$$\delta = \dim ( G \Gamma _ { f } ) - \dim ( \Gamma _ { f } ) = \dim ( G \Gamma _ { f } ) - p.$$

Also, we have 0 ≤ δ ≤ min { r, q } ([5], p 276-277).

If δ = 0, then u = f ( x ) is an invariant solution and if 0 &lt; δ &lt; min { r, q } , then u = f ( x ) is a partially invariant solution.

In order to calculate the PISs, first of all, it is necessary to classify the symmetry group into conjugacy classes. For obtaining the PISs with the defect structure δ , those subgroups H ⊂ G which have the property that if the dimension of the orbits of H on the space X × U is r , then the dimension of the orbits is r -δ on the space X , must be selected [19]. Let H be a subgroup with this property that mentioned above and h be its Lie algebra with infinitesimal generators { v , ..., v 1 s } . Hence, we can obtain a complete set of functionally independent invariants of the form

$$\{ \xi _ { i } ( x ), I _ { j } ( x, u ) \},$$

where i = 1 , ..., p + δ -s and j = 1 , ..., q -δ . Then we have

$$\text{rank} \left ( \frac { \partial I _ { j } ( x, u ) } { \partial u } \right ) = q - \delta = q ^ { \prime }. \quad \text{--} \cdots \text{--} \text{rrc} \quad \text{--} \cdots \cdots \cdots \quad \text{--} \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \ddots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdets \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \collections \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots$$

If u = f ( x ) is a function, then the manifold H Γ f can be expressed in terms of the invariants (8). So, we have:

$$I _ { j } ( x, u ) = f _ { j } ( \xi _ { i } ( x ) )$$

where the functions f j are arbitrary. Now, by applying the implicit function theorem, we conclude that

$$u ^ { i _ { \alpha } } = U ^ { i _ { \alpha } } ( x, u ^ { j \beta }, f _ { j } ( \xi _ { i } ( x ) ) )$$

where α = 1 , ..., q ′ and β = 1 , ..., δ . The remaining dependent variables only depend on the original independent variables:

$$u ^ { j _ { \beta } } = U ^ { j _ { \beta } } ( x _ { 1 }, \dots, x _ { p } ), \quad \ \beta = 1, \dots, \delta.$$

Now, the derivatives of the functions u , ..., u 1 q with respect to the new variables which are obtained from equations (11) and (12) must be calculated. Hence, by substituting these quantities into the original system, a new system is obtained, involving the q ′ functions f j and the invariants ξ i . The resulted equations are not generally consistent, so that the compatibility conditions must be computed. Consequently, a system of PDEs is deduced from these constraints which is denoted by ∆ /H . On the other hand, a system of PDEs is resulted from (12) denoted by ∆ . ′ Now, the system ∆ /H must be solved first. Then corresponding to each of the solution of this system, the system ∆ ′ will be solved. Finally, the partially invariant solutions are obtained by substituting the resulted solutions into equations (11) and (12). For more details about the methods of determining the PISs refer to [20, 5].

## 3. Classification of Subalgebras for the system (3)

In this section, we want to classify the subgroups of symmetry group of the system (3), into conjugacy classes. Searching for the invariant solutions can be regarded as the main motivation of computing the symmetries of a differential equation. As it is well known, the problem of classifying invariant solutions is equivalent to the problem of classifying the subgroups of the full symmetry group under conjugation. Let H and ˜ be two connected, H s -dimensional Lie subgroups of the Lie group G with corresponding Lie subalgebras h and ˜ of the Lie algebra h g of G . Let g ∈ G , then ˜ H = gHg -1 are conjugate subgroups if and only if ˜ = Ad( h g ) · h are conjugate subalgebras, where Ad( g ) is adjoint represen Hence, the problem of determining an optimal system of subgroups is equivalent to that of obtaining an optimal system of subalgebras, and so we focus on it [4, 5]. The latter problem, tends to obtain a list of conjugacy inequivalent subalgebras with the property that any other subalgebra is equivalent to only a unique member of the list under some element of the adjoint representation for some element of the investigated Lie group.

## 3.1. Optimal system of one-dimensional subalgebras for the system (3)

Indeed, for one-dimensional subalgebras, the classification problem is necessarily the same as the problem of classifying the orbits of the adjoint representation. Thus, an optimal set of subalgebras is constructed if we select just one representative from each family of equivalent subalgebras. Consequently, the associated set of invariant solutions is then the minimal list from which we can obtain all other invariant solutions of one-dimensional subalgebras simply via transformations.

In [17], we have presented a comprehensive analysis of this problem and have constructed an optimal system of one-dimensional subalgebras for the system (3) as follows :

Theorem 3.1. An optimal system of one-dimensional Lie subalgebras of the system (3) is provided by the following generators:

$$\text{now many genericours} \colon \\ 1 ) \, X ^ { 1 } & = X _ { 7 }, \quad & \, 8 ) \, X ^ { 8 } \ = X _ { 4 } + a X _ { 5 } + b X _ { 6 } + c X _ { 7 }, \\ 2 ) \, X ^ { 2 } & = X _ { 1 } + a X _ { 7 }, \quad & \, 9 ) \, X ^ { 9 } \ = \varepsilon X _ { 1 } + X _ { 4 } + a X _ { 5 } + b X _ { 6 } + c X _ { 7 }, \\ 3 ) \, X ^ { 3 } & = X _ { 2 } + a X _ { 7 }, \quad & \, 1 0 ) \, X ^ { 1 0 } & = X _ { 3 } + a X _ { 5 } + b X _ { 6 } + c X _ { 7 }, \\ 4 ) \, X ^ { 4 } & = X _ { 6 } + a X _ { 7 } \quad & \, 1 1 ) \, X ^ { 1 1 } & = \varepsilon X _ { 2 } + X _ { 3 } + a X _ { 5 } + b X _ { 6 } + c X _ { 7 }, \quad & \, ( 1 3 ) \\ 5 ) \, X ^ { 5 } & = \varepsilon X _ { 1 } + X _ { 6 } + a X _ { 7 }, \quad & \, 1 2 ) \, X ^ { 1 2 } & = X _ { 3 } + \varepsilon X _ { 4 } + a X _ { 5 } + b X _ { 6 } + c X _ { 7 }, \\ 6 ) \, X ^ { 6 } & = X _ { 5 } + a X _ { 6 } + b X _ { 7 }, \quad & \, 1 3 ) \, X ^ { 1 3 } & = \varepsilon X _ { 2 } + X _ { 3 } + \varepsilon ^ { \prime } X _ { 4 } + a X _ { 5 } + b X _ { 6 } + c X _ { 7 }. \\ 7 ) \, X ^ { 7 } & = \varepsilon X _ { 2 } + X _ { 5 } + a X _ { 6 } + b X _ { 7 },$$

- 7) X = εX 2 + X 5 + aX 6 + bX , 7

where ε and ε ′ are ± 1 and a, b, c ∈ R are arbitrary numbers [17].

## 3.2. Optimal system of two-dimensional subalgebras for the system (3)

In this paper, we need to classify the two-dimensional subalgebras. Because, we want to calculate those PISs which have the defect structure δ = 1 and the reduced system ∆ /H is a system of ordinary differential equations. Since p = 2, δ = 1 and p + δ -s = 1, then we have s = 2. So, we should consider the twodimensional subgroups.

Consequently, constructing the two-dimensional optimal system, i.e., classification of the two-dimensional subalgebras of g is our next step. This process is performed by choosing one of the vector fields as stated in theorem (3.1). Let us consider X 1 (or X i , i = 2 , · · · , 13). Corresponding to it, an optional vector field Y = b 1 X 1 + · · · + b 7 X 7 is selected, so we must have

$$[ \mathbf X ^ { 1 }, Y ] = \lambda \mathbf X ^ { 1 } + \mu Y.$$

Equation (14) leads us to the system

$$C _ { j k } ^ { i } \alpha _ { j } a _ { k } = \lambda a _ { i } + \mu \alpha _ { i } & & ( i = 1, \cdots, 7 ). & & ( 1 5 )$$

where the constant coefficients C i jk are the structure constants. The solutions of the system (15), give one of the two-dimensional generators and the second generator is X 1 (or X i , i = 2 , · · · , 13) if selected.

After the construction of all two-dimensional subalgebras, for every vector fields of theorem (3.1), they need to be simplified by the action of adjoint transformations in the manner analogous to the way of onedimensional optimal system. Hence, we can state the following theorem:

## Theorem 3.2. Any two-dimensional subalgebras of (4) is conjugate to precisely one of the following subalgebras:

| A 1 1 : 〈 X 1 , X 3 + αX 6 + βX 7 〉 , A 4 1 : 〈 X 1 , X 3 + /epsilon1X 5 + X 6 + αX 7 〉 , A 7 1 : 〈 X 1 , X 7 〉 ,   | A 2 1 : 〈 X 1 , X 2 + αX 5 + βX 7 〉 , A 5 1 : 〈 X 1 , X 2 + αX 3 + βX 7 〉 ,        | A 3 1 : 〈 X 1 , X 5 + αX 7 〉 , A 6 1 : 〈 X 1 , X 6 + αX 7 〉 ,   |
|---------------------------------------------------------------------------------------------------------------------------|----------------------------------------------------------------------------------------|---------------------------------------------------------------------|
| A 1 2 : 〈 X 2 , X 3 + αX 6 + βX 7 〉 , A 4 2 : 〈 X 2 , X 3 + /epsilon1X 4 + X 6 + αX 7 〉 , A 7 2 : 〈 X 2 , X 7 〉 ,   | A 2 2 : 〈 X 2 , X 1 + /epsilon1X 4 + βX 7 〉 , A 5 2 : 〈 X 2 , X 1 + X 6 + αX 7 〉 , | A 3 2 : 〈 X 2 , X 4 + αX 7 〉 , A 6 2 : 〈 X 2 , X 6 + αX 7 〉 ,   |
| A 1 3 : 〈 X 6 , X 3 + αX 7 〉 , A 4 3 : 〈 X 6 , X 1 + αX 7 〉 ,                                                         | A 2 3 : 〈 X 6 , X 4 〉 , A 5 3 : 〈 X 6 , X 2 〉 ,                                    | A 3 3 : 〈 X 6 , X 5 〉 , A 6 3 : 〈 X 6 , X 7 〉 ,                 |
| A 1 4 : 〈 εX 1 + X 6 , X 2 〉 ,                                                                                          | A 2 4 : 〈 εX 1 + X 6 , X 5 〉 ,                                                       | A 3 4 : 〈 εX 1 + X 6 , X 7 〉 ,                                    |
| A 1 5 : 〈 X 5 , X 3 + αX 6 + βX 7 〉 , A 4 5 : 〈 X 5 , X 7 〉 ,                                                         | A 2 5 : 〈 X 5 , X 1 + αX 6 + βX 7 〉 ,                                                | A 3 5 : 〈 X 5 , X 6 + αX 7 〉 ,                                    |
| A 1 6 : 〈 εX 2 + X 5 , X 3 + 1 2 X 6 + αX 7 〉 ,                                                                         | A 2 6 : 〈 εX 2 + X 5 , X 1 + αX 7 〉 ,                                                | A 3 6 : 〈 εX 2 + X 5 , X 7 〉 ,                                    |
| A 1 7 : 〈 X 4 , X 3 + αX 6 + βX 7 〉 , A 4 7 : 〈 X 4 , X 7 〉 ,                                                         | A 2 7 : 〈 X 4 , X 2 + αX 3 + βX 7 〉 ,                                                | A 3 7 : 〈 X 4 , X 6 + αX 7 〉 ,                                    |
| A 1 8 : 〈 εX 1 + X 4 , X 3 +2 X 6 + αX 7 〉 ,                                                                            | A 2 8 : 〈 εX 1 + X 4 , X 2 + αX 7 〉 ,                                                | A 3 8 : 〈 εX 1 + X 4 , X 7 〉 ,                                    |
| A 1 9 : 〈 X 3 , X 2 + αX 7 〉 , A 4 9 : 〈 X 3 , X 6 + αX 7 〉 ,                                                         | A 2 9 : 〈 X 3 , X 5 〉 , A 5 9 : 〈 X 3 , X 7 〉 ,                                    | A 3 9 : 〈 X 3 , X 1 〉 , A 6 9 : 〈 X 3 , X 4 〉 ,                 |
| A 1 10 : 〈 εX 2 + X 3 , X 1 〉 , A 4 10 : 〈 X 2 , X 3 + αX 7 〉 ,                                                       | A 2 10 : 〈 εX 2 + X 3 , X 4 〉 ,                                                      | A 3 10 : 〈 εX 2 + X 3 , X 7 〉 ,                                   |
| A 1 11 : 〈 X 3 + εX 4 , εX 5 + X 6 - 2 X 7 〉 , A 4 11 : 〈 X 3 + εX 4 , X 3 + X 6 + αX 7 〉 ,                           | A 2 11 : 〈 X 3 + εX 4 , X 2 + αX 7 〉 , A 5 11 : 〈 X 3 + εX 4 , X 1 + εX 2 〉 ,      | A 3 11 : 〈 X 3 + εX 4 , X 7 〉 ,                                   |
| 1 12 : εX 2 + X 3 + ε ′ X 4 , X 1 + ε ′ X 2 ,                                                                             | 2 12 : εX 2 + X 3 + ε ′ X 4 , X 7 ,                                                    |                                                                     |

A 1 12 : 〈 εX 2 + X 3 + ε X , X ′ 4 1 + ε X ′ 2 〉 , A 2 12 : 〈 εX 2 + X 3 + ε X , X ′ 4 7 〉 , where α and β are arbitrary constants, ε and ε ′ are ± 1 and /epsilon1 is ± 1 or 0 .

Proof. Each of the two-dimensional subalgebra has two generators. For classifying two-dimensional subalgebras, we must select one of the generators from the list of one-dimensional optimal system (13) and another generator must be taken optionally.

Suppose that h = span { X,Y } is a two-dimensional subalgebra of g where X is a one-dimensional subalgebra which is chosen from the list (13) and Y is an optional vector defined by: Y = b 1 X 1 + · · · + b 7 X 7 . Now, we must simplify h as much as possible by imposing various adjoint transformations on it [4].

Each adjoint transformation is a linear map F s i : g → g defined by X → Ad(exp( sX i ) .X ), for i = 1 , · · · , 7.

↦

Since this procedure is so lengthy, we only explain one of the cases in the following.

Case1 : If X = X 1 + aX 7 then

$$\mathfrak { h } = \langle X, Y \rangle & = \langle X _ { 1 } + a X _ { 7 }, \sum _ { i = 1 } ^ { 7 } b _ { i } X _ { i } \rangle \\ & \quad = \langle X _ { 1 }, b _ { 2 } X _ { 2 } + b _ { 3 } X _ { 3 } + b _ { 4 } X _ { 4 } + b _ { 5 } X _ { 5 } + b _ { 6 } X _ { 6 } + b _ { 7 } X _ { 7 } \rangle \\.$$

So, we have:

- a) If b 2 = b 3 = b 4 = b 5 = b 6 = 0, then we have h = 〈 X ,b X 1 7 7 〉 = 〈 X ,X 1 7 〉 . Since [ X ,X 1 7 ] = 0 = 0 X 1 +0 X 7 , then h is closed under the Lie bracket. So, h is reduced to the case A 7 1 .

/negationslash

- b) If b 2 = b 3 = b 4 = b 5 = 0 and b 6 = 0, then h = 〈 X ,b X 1 6 6 + b 7 X 7 〉 = 〈 X ,X 1 6 + αX 7 〉 . Since [ X ,X 1 6 + αX 7 ] = 0, then h is closed under the Lie bracket. So, h is reduced to the case A 6 1 .

/negationslash

- c) If b 2 = b 3 = b 4 = b 6 = 0 and b 5 = 0, then h = 〈 X ,b X 1 5 5 + b 7 X 7 〉 = 〈 X ,X 1 5 + αX 7 〉 . Since [ X ,X 1 5 + αX 7 ] = 0, then h is closed under the Lie bracket. So, h is reduced to the case A 3 1 .

/negationslash

- d) If b 2 = b 3 = b 4 = 0 and b 6 = 0, then we can make the coefficient of X 5 vanish by F s 5 5 ; By setting s 5 = 1 b 6 . So, we have h = 〈 X ,b X 1 6 6 + b 7 X 7 〉 which is reduced to the case A 6 1 like (b).

/negationslash

- e) If b 2 = b 3 = 0 and b 4 = 0, then h = 〈 X ,b X 1 4 4 + · · · + b 7 X 7 〉 and [ X ,b X 1 4 4 + · · · + b 7 X 7 ] = b 4 X 2 = rX 1 + ( s b 4 X 4 + · · · + b 7 X 7 ) for any r, s ∈ R . So, h is not closed under the Lie bracket and we have not any two-dimensional subalgebras in this case.

/negationslash

/negationslash

- f) If b 2 = b 4 = 0, b 3 = 0 and b 6 = 1, then we can make the coefficient of X 5 vanish by F s 5 5 ; By setting s 5 = b 5 -1+ b 6 . Also, by scaling if necessary, we can assume that b 3 = 1. Then we have h = 〈 X ,X 1 3 + b 6 X 6 + b 7 X 7 〉 and [ X ,X 1 3 + b 6 X 6 + b 7 X 7 ] = X 1 . So, h is closed under the Lie bracket and the case A 1 1 is concluded.

/negationslash

/negationslash

- g) If b 2 = b 4 = 0, b 3 = 0 and b 6 = 1, then by scaling if necessary, we can assume that b 3 = 1. Also, the coefficient of X 5 can be vanished or be ± 1 by F s 6 6 ; By setting s 6 = -ln | b 5 | . Then we have h = 〈 X ,X 1 3 + εX 5 + X 6 + b 7 X 7 〉 and [ X ,X 1 3 + εX 5 + X 6 + b 7 X 7 ] = X 1 . So, h is closed under the Lie bracket and the case A 4 1 is deduced.

/negationslash

- h) If b 3 = b 4 = b 6 = 0 and b 2 = 0, then by scaling if necessary, we can assume that b 2 = 1. Also, we have h = 〈 X ,X 1 2 + b 5 X 5 + b 7 X 7 〉 and [ X ,X 1 2 + b 5 X 5 + b 7 X 7 ] = 0. So, h is closed under the Lie bracket and the case A 2 1 is resulted.

/negationslash

- i) If b 4 = b 6 = 0, b 2 = 0 and b 3 = 0, then by scaling if necessary, we can assume that b 2 = 1. Also we can make the coefficient of X 5 vanish by F s 5 5 ; By setting s 5 = -b 5 b 3 . Then we have h = 〈 X ,X 1 2 + b 3 X 3 + b 7 X 7 〉 and [ X ,X 1 2 + b 3 X 3 + b 7 X 7 ] = b 3 X 1 . So, h is closed under the Lie bracket and the case A 5 1 is resulted .

/negationslash

/negationslash

- j) If b 4 = 0, b 6 = 0 and b 3 = 1, then by scaling if necessary, we can assume that b 6 = 1. We can make the coefficient of X 2 vanish by F s 2 2 ; By setting s 2 = -1 b 6 . Also, the coefficient of X 5 can be vanished or be ± 1 by F s 6 6 ; By setting s 6 = -ln b | 5 | . So, we have h = 〈 X ,X 1 3 + εX 5 + X 6 + b 7 X 7 〉 which is reduced to the case A 4 1 similar to (g).

/negationslash

- k) If b 4 = 0, b 6 = 0 and b 3 = 1, then we can make the coefficient of X 2 and X 5 vanish by F s 2 2 and F s 5 5 ; By setting s 2 = -1 b 6 and s 5 = -b 5 b 3 -1 , respectively. By scaling if necessary, we have h = 〈 X ,X 1 3 + b 6 X 6 + b 7 X 7 〉 that is reduced to the case A 1 1 similar to (f).

/negationslash

In each case, h can not be simplified any more by F s i i , i = 1 , ..., 7. By a similar method, we can find the two-dimensional subalgebras for the other 11 cases.

## 4. Computation of the partially invariant solutions for the system (3)

In this section, we will calculate the PISs of the system (3). For example, consider the Lie subalgebra A 7 1 : 〈 X ,X 1 7 〉 . The set of functionally independent invariants for A 7 1 is a set of functions I with the following property: X 1 ( I ) = X 7 ( I ) = 0. By calculating these functions, a set of functionally independent invariants is obtained as follows:

$$\{ t, \frac { b } { a }, \frac { c } { a } \}$$

So, we have:

$$\text{ive} \colon & & \text{rank} \left ( \frac { \partial ( \frac { b } { a }, \frac { c } { a } ) } { \partial ( a, b, c ) } \right ) = \text{rank} \left ( \begin{array} { c c } - b / a ^ { 2 } & - c / a ^ { 2 } \\ 1 / a & 0 \\ 0 & 1 / a \end{array} \right ) = 2 = q - \delta. & & ( 1 8 ) \\ \colon & & 3, \text{then} \, \delta = 1. \text{ The equations corresponding to the relation } ( 1 0 ) \, \text{ are}$$

Since q = 3, then δ = 1. The equations corresponding to the relation (10) are

$$\frac { b } { a } = f ( t ), \quad \frac { c } { a } = g ( t ),$$

and the solutions corresponding to the equations (11) and (12) are expressed as follows:

$$b = a f ( t ), \quad c = a g ( t ), \quad a = a ( x, t ).$$

Note that a x, t ( ) is an arbitrary function.

Now, we can compute the derivatives of the functions a, b and c from relations (20):

$$\begin{array} { c c c } b _ { 1 } = a _ { 1 } f, & b _ { 2 } = a _ { 2 } f + a f ^ { \prime }$$

Hence, by substituting the above relations into the system (3), we obtain:

$$\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cd$$

The consistency conditions, obtained from the system (22), conclude the following equations

$$\begin{array} { c c c c c c c } \cdot & \cdot & \cdot & \cdot & \cdot \\ & a ^ { 2 } f ^ { \$$

and these inequations

$$f \neq 0, \quad g \neq 0, \quad a \neq 0, \quad f - g ^ { 2 } \neq 0. \quad \quad \quad$$

/negationslash

/negationslash

/negationslash

/negationslash

Equations (23) form a system of ODEs. By solving this system, we obtain four types of solutions as follows:

$$\text{ons} \, ( 2 3 ) \, \text{form} \, a \$$

Now, by applying the relations (20), the partially invariant solutions associated to the subalgebra A 7 1 can be obtained for the system (3).

$$\begin{array} {$$

  These partially invariant solutions seem to be trivial. In a similar way, we can calculate nontrivial PISs by applying other subalgebras listed in theorem (3.2). Some other of these PISs are presented in Table 1.

## 5. Non-reducible partially invariant solutions

In this section, we will deal with those partially invariant solutions which are not invariant with respect to some of the subgroups of G . They are called non-reducible PISs and they are usually rare. For a reducible partially invariant solution u = f ( x ), we can determine a subgroup H ′ ⊂ H which u is a H ′ -invariant solution and dim( H ′ Γ ) f /greaterorequalslant dim( H ) -δ = s -δ ([5], page 290). So, reducible PISs can be obtained via the method of similarity reduction from the reduced system involving p -dim( H ′ Γ ) f /lessorequalslant p + δ -s independent variables, which is indeed easier than obtaining them by PISs method. For example, the PISs that we computed in the previous section are non-reducible PISs. Because, reducible PISs must be invariant with respect to the one-parameter subgroups of its Lie group. These subgroups have Lie subalgebras of A 7 1 generated by an infinitesimal generator of the form αX 1 + βX 7 where α, β ∈ R . In [17], we obtained the invariant solutions of the system (3) with respect to this generator by applying the method of similarity reduction:

$$\text{use system} \, ( \mathcal { J } \, \text{with respect$$

 Since, none of these solutions is similar to the PISs in (26), we conclude that PISs in (26) are non-reducible PISs. Several other non-reducible partially invariant solutions of the system (3) are listed in the following table.

Table 1: Non-reducible partially invariant solutions of (3).

| Subalgebras           | Invariants                                      | Dependent variables                                                               | PISs                                                                                                                     |
|-----------------------|-------------------------------------------------|-----------------------------------------------------------------------------------|--------------------------------------------------------------------------------------------------------------------------|
| 〈 X 2 , X 7 〉       | { x, b a , c a }                                | a = a ( x, t ) b = af ( x ) c = ag ( x )                                          | a = c 1 x + c 2 b = c 1 c 2 3 x +( ∗ ) c = c 3 ( c 1 x + c 2 ) + c 1 c 4                                                 |
| 〈 X 5 , X 7 〉       | { t, c 2 - ba b 2 , - ct + bx b }               | a = b ( x - g ( t )) 2 t 2 - bf ( t ) b = b ( x, t ) c = b t ( x - g ( t ))       | a = ( c 1 + c 2 t 3 )( x - c 3 ) 2 t 3 - ( ∗∗ ) b = c 1 + c 2 t 3 t c = ( c 1 + c 2 t 3 )( x - c 3 ) t 2                 |
| 〈 X 2 , X 6 + X 7 〉 | { x, b a 3 , c a 2 }                            | a = a ( x, t ) b = a 3 f ( x ) c = a 2 g ( x )                                    | a = 4( t + c 1 ) ( c 2 x + c 3 ) 2 b = 4 c 2 2 ( t + c 1 ) 3 ( c 2 x + c 3 ) 4 c = 4 c 2 ( t + c 1 ) 2 ( c 2 x + c 3 ) 3 |
| 〈 X 4 , X 7 〉       | { x, - at + cx xa , at 2 - 2 xtc + bx 2 x 2 a } | a = a ( x, t ) b = a ( g ( x ) +2 t x f ( x ) + t 2 x 2 ) c = a ( f ( x ) + t x ) | a = c 1 + c 2 x 3 x b = ( c 1 + c 2 x 3 )( t + c 3 ) 2 x 3 +( ∗ ∗ ∗ ) c = ( c 1 + c 2 x 3 )( t + c 3 ) x 2               |

where c i 's are arbitrary constants and

- ∗ = ( c 5 c 1 -c 2 3 c 2 )ln( c 1 x + c 2 ) + c 6 , ∗∗ = ( t c 4 + c 5 )(ln( t +( c 1 c 2 ) ( 1 3 ) ) 2 -ln( t 2 -t ( c 1 c 2 ) ( 1 3 ) +( c 1 c 2 ) ( 2 3 ) ) + 2 √ 3tan -1 ( 2 c 2 t √ 3 c 1 ( c 1 c 2 ) ( 2 3 ) -1 √ 3 )) , √

$$* * * = x ( c _ { 4 } + c _ { 5 } ) ( \ln ($$

## Conclusion

Partially invariant solutions can be regarded as the natural extension of the invariant ones and the method of obtaining them is algorithmic. Some of the PISs are not invariant with respect to the lower dimensional subalgebras. They are called non-reducible PISs and can be applied in order to obtain some new solutions for an arbitrary system of PDEs. For example, in the present paper, we have obtained a non-reducible PISs for the system (3) as follows:

$$a = \frac { 4 ( t + c _ { 1 } ) } { ( c _ { 2$$

Consequently, according to relation (1), we could determine a general form for the metric associated to a set of four-dimensional Einstein Walker manifolds as follows:

$$g _ { a, b, c } & \coloneqq 2 ( d x \circ d y + d t \circ d z ) + \frac { 4 ( t + c _ { 1 } ) } { ( c _ { 2 } x + c _ { 3 } ) ^ { 2 } } d y \circ d y \\ & \quad + \frac { 4 c _ { 2 } ^ { 2 } ( t + c _ { 1 } ) ^ { 3 } } { ( c _ { 2 } x + c _ { 3 } ) ^ { 4 } } d z \circ d z + \frac { 8 c _ { 2 } ( t + c _ { 1 } ) ^ { 2 } } { ( c _ { 2 } x + c _ { 3 } ) ^ { 3 } } d y \circ d z. \\ \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot$$

This set of solutions can not be obtained from the similarity reduction method. In a similar manner, we can deduce some other four-dimensional Einstein Walker manifolds from Table 1.

## Acknowledgements

It is a pleasure to thank the anonymous referees for their constructive suggestions and helpful comments which have improved the presentation of the paper. The authors wish to express their sincere gratitude to Fatemeh Ahangari for her useful advise and suggestions.

## References

- [1] S. Lie, On integration of a class of linear partial differential equations by means of definite integrals, translation by N.H. Ibragimov, Arch Math. 6 (1881) 328-368.
- [2] G. Bluman, A.F. Cheviakov, S. Anco, Application of symmetry methods to partial differential equations, Springer, New York, 2010.
- [3] A. Kushner, V. Lychagin, V. Rubstov, Contact geometry and nonlinear differential equations, Cambridge University Press, Cambridge, 2007.
- [4] P.J. Olver, Applications of Lie Groups to Differential Equations, Springer, New York, 1986.
- [5] L.V. Ovsiannikov, Group Analysis of Differential Equations, Academic Press. New York, 1982.
- [6] C.P. Boyer, J. D. Finley III, and J. F. Pleba´ski, Complex general relativity, n h and hh spaces - a survey of one approach, Gen. Relativity Gravitation, 2 (1980) 241-281, Plenum, New York-London.
- [7] J.D. Finley III, and J. F. Pleba´ski, The intrinsic spinorial structure of hyperheavens, J. Math. Phys. 17 (1976), 2207-2214. n
- [8] R. Abounasr, A. Belhaj, J. Rasmussen, and E.H. Saidi, Superstring theory on pp-waves with ADE geometries, J. Phys. A 39 (2006), 2797-2841.
- [9] J. Kerimo, AdS pp-waves, J. High Energy Phys. (2005), 025, 18 pp.
- [10] J. Michelson, and X. Wu, Dynamics of antimembranes in the maximally supersymmetric eleven-dimensional pp-wave, J. High Energy Phys. (2006), 028, 37 pp.
- [11] V. Pravda, A. Pravdova, A. Coley, and R. Milson, All spacetimes with vanishing curvature invariants, Classical Quantum Gravity 19 (2002), 6213-6236.
- [12] Th. Leistner, Screen bundles of Lorentzian manifolds and some generalizations of pp-waves, J. Geom. Phys. 56 (2006), 2117-2134.

- [13] M. Brozos-V´ azquez, E. Garcia-Rio, P. Gilkey, S. Nikˇ cevi´, c R. V´ azquez-Lorenzo, The Geometry of Walker Manifolds. Synthesis Lectures on Mathematics and Statistics, Vol. 5, Morgan and Claypool Publ, 2009.
- [14] Y. Matsushita, Walker 4-manifolds with proper almost complex structures, J. Geom. Phys. 55 (2005), 385-398.
- [15] M. Brozos-V´ azquez, E. Garcia-Rio, and R. V´zquez-Lorenzo, Conformally Osserman four-dimensional manifolds whose a conformal Jacobi operators have complex eigenvalues, Proc. R. Soc. A. 462 (2006), 1425-1441.
- [16] P. Gilkey, The geometry of curvature homogeneous pseudo-Riemannian manifolds. ICP Advanced Texts in Mathematics, 2. Imperial College Press, London, 2007.

[17] M. Nadjafikhah, M. Jafari, Some general new Einstein Walker manifolds, arXiv:1206.3730v1 [math.DG] 17 Jun 2012.

- [18] S.V. Meleshko, Methods for Constructing Exact Solutions of Partial Differential Equations, Springer Science+Business Media, Inc, NewYork, 2005.
- [19] A.M. Grundland, L. Lalague, Invariant and partially-invariant solutions of the equations describing a non-stationary and isentropic flow for an ideal and compressible fluid in (3 + 1) dimensions, J. Phys. A: Math. Gen. 29 (1996) 1723-1739.
- [20] A.M. Grundland, P. Tempesta, P. Winternitz, Weak transversality and partially invariant solutions, J. Math. Phys. 44(6) (2003) 2704-2722.
Mehdi Nadjafikhah Mehdi Jafari
Computation of Partially Invariant Solutions for the Einstein Walker Manifolds' Identifying Equations
In this paper, partially invariant solutions (PISs) method is applied in
order to obtain new four-dimensional Einstein Walker manifolds. This method is
based on subgroup classification for the symmetry group of partial differential
equations (PDEs) and can be regarded as the generalization of the similarity
reduction method. For this purpose, those cases of PISs which have the defect
structure delta=1 and are resulted from two-dimensional subalgebras are
considered in the present paper. Also it is shown that the obtained PISs are
distinct from the invariant solutions that obtained by similarity reduction
method.