The recent observations of luminosity of type Ia Supernovae (SNe Ia) [1] -[8] indicate that the universe is currently undergoing an accelerated expansion. The dark energy (DE) with negative pressure is responsible for this scenario. Researchers have developed number of DE models of the universe in general relativity [9] [10] and also in different theories of gravitation [11] -[15] . These models have a lot of significance in theoretical study of the structure of the universe. It is considered that DE is to be the best candidate to explain cosmic acceleration. It is known from the history that our present universe consists of about DE and about dark matter (DM) [16] [17] . DE is usually characterized by the EoS parameter defined by the equation of state, where is the fluid pressure and is energy density. Usually EoS parameter is assumed to be a constant with the values −1, 0, 1/3 and +1 for vacuum, dust, radiation and stiff matter dominated universe respectively. However, latest observations from SNe Ia data [18] -[20] indicate that is not constant. Many authors [21] -[28] have obtained DE models in general relativity with variable EoS parameter.

Magnetic field plays a vital role in the description of the energy distribution in the universe as it contains highly ionized matter. Strong magnetic fields can be created due to adiabatic compression in cluster of galaxies. Large scale magnetic fields give rise to anisotropies in the universe. The magnetic field has the significant role in the dynamics of the universe depending on the direction of the field lines [29] [30] . Many research works have studied the importance of magnetic field for various astrophysical phenomena. Misra and Radhakrishna [31] have obtained the analogous relation between some components of metric and electromagnetic potentials for source free Einstein-Maxwell field described by Einstein-Rosen metric. Thorne [32] , Jacobs [33] [34] , Collins [35] , Roy and Prakash [36] have investigated magnetized cosmological models for perfect fluid distributions in general relativity. Milaneschhi and Fabbri [37] studied the anisotropy and polarization properties of CMB radiations in homogeneous Bianchi type-I cosmological model. Roy et al. [38] investigated Bianchi type-I cosmological models containing perfect fluid and magnetic field directed along x-axis. Adhav et al. [39] have investigated Bianchi type-III model with DE from a wet dark fluid [WDF] in presence and absence of magnetic field in general theory of relativity. Katore et al. [40] have explored Einstein-Rosen cosmological model with magnetized anisotropic DE.

Bianchi type cosmological models are important in the sense that, these are homogeneous and anisotropic, from which the process of isotropization of the universe is studied through the passage of time. Moreover, from the theoretical point of view anisotropic universe has a greater generality than isotropic models. The simplicity of the field equations made Bianchi space-times useful in constructing models of spatially homogeneous and anisotropic cosmologies. Reddy et al. [41] have studied Bianchi-II, VIII and IX models in scale-covariant theory of gravitation. Chakraborty [42] , Raj Bali and Dave [43] , Raj Bali and Yadav [44] have studied Bianchi type-IX string as well as viscous fluid models in general relativity. Pradhan [45] has studied some homogeneous Bianchi type-IX viscous fluid cosmological models with varying. Tyagi et al. [46] have obtained Bianchi type-IX string cosmological models for perfect fluid distribution in general relativity. Rao et al. [47] have obtained Bianchi type-II, VIII and IX DE cosmological models in Saez-Ballester theory of gravitation. Recently, Ghate and Sontakke [48] -[50] have studied Bianchi type-IX cosmological model with anisotropic DE in Lyra geometry, model with binary mixture of perfect fluid and dark energy and DE model in a Brans-Dicke theory of gravitation repectively.

In the last few decades, alternative theories to Einstein’s theory of gravitation have developed mainly scalar tensor theories proposed by Brans and Dicke [51] , Nordtvedt [52] , Wagoner [53] , Ross [54] , Dun [55] , Saez and Ballester [56] , Barber [57] and La & Steinhardt [58] . Among them Brans-Dicke and Saez-Ballester theories are considered to be viable alternatives to general relativity. Brans-Dicke theory includes a long range scalar field interacting equally with all forms of matter (with the exception of electromagnetism) while in Saez-Ballester theory, metric is coupled with a dimensionless scalar field in a simple manner. This coupling gives satisfactory description of weak fields. This theory suggests a possible way to solve the missing-matter problem in non-flat FRW cosmologies. In earlier literature, cosmological models in Saez-Ballester theory of gravitation have been studied by Singh and Agrawal [59] , Shri Ram and Tiwari [60] , Singh and Shri Ram [61] . In recent years, Reddy and Naidu [62] , Adhav et al. [63] , Katore et al. [64] Pradhan and Singh [65] are some of the authors who have obtained the solutions in Saez-Ballester theory in different contexts.

In this paper, we have studied the solutions of Bianchi type-IX universe with variable in Saez-Ballester theory of gravitation in the presence and absence of magnetic field of energy density together with constant deceleration parameter. A special law of variation of Hubble’s parameter proposed by Berman [66] is used to solve the field equations. The physical and kinematical properties of the model are also discussed. The out-line of the paper is as follows: In section 2, the model and field equations are described. The solution of field equations are presented in section 3 and section 4 concludes the findings.

Bianchi type-IX metric is considered in the form

where are scale factors and are functions of cosmic time.

The model has one transverse direction and two equivalent longitudinal directions and.

The field equations in Saez-Ballester theory of gravitation are

and the scalar field satisfies the equation

where is the Einstein tensor, R is the scalar curvature, and n are constants, is the stress tensor of matter, is a dimensionless coupling constant; comma and semicolon denote partial and covariant differentiation respectively.

Also, we have energy-conservation equation

is a consequence of the field equations.

King and Coles [30] and Jacobs [34] used the magnetized perfect fluid energy momentum tensor to discuss the effects of magnetic field on the evolution of the universe filled with perfect fluid.

The energy-momentum tensor for the magnetized anisotropic DE fluid is in the form

where is the energy density of magnetic fluid, is the energy density of the fluid and, , are pressures on x, y, z axes respectively.

The equation of state for an anisotropic fluid is taken of the form, where is not necessarily constant (Carroll et al. [67] );

Here is the deviation free parameter., , are the directional EoS parameters on, , axes respectively.

For Bianchi type-IX metric, using Equation (6), the field Equations (2) and (3) takes the form

where over dot denotes the differentiation with respect to.

We have the following equation from the Bianchi identity:

The field Equations (7)-(10) are a system of four equations with six unknown parameters, , , , and. The system is thus initially undetermined and we need additional constraints to obtain the solution of field equations.

We assumed that the magnetized DE is minimally interacting, hence the Bianchi identity has been split into two separately additive conserved components namely, the conservation of the energy-momentum tensor for the anisotropic fluid and for the magnetic field (King and Coles [30] ).

Finally, we constrain the system of equations with a law of variation for the average Hubble’s parameter that yields a constant value of deceleration parameter. Such types of relation have already been considered by Berman [66] , Berman and Gomide [68] for solving FRW models. Later on many authors (Singh et al. [69] -[71] , Singh and Baghel [72] ) have studied flat FRW and Bianchi type models by using the special law of Hubble parameter that yields constant value of deceleration parameter.

The average scale factor R of Bianchi type-IX metric is given by

The proper volume V is defined by

We defined the generalized mean Hubble’s parameter H as

where, and are the directional Hubble parameters in the direction of x, y, z axes respectively.

From Equations (14) and (16), we obtain

Since, the line element (1) is completely characterized by Hubble’s parameter H. Therefore, let us consider that, the mean Hubble parameter H is related to the average scale factor by the relation

where and are constants.

An important observational quantity is the deceleration parameter, which is defined as

From Equations (17) and (18), we obtain

From Equations (19), (20), (21) we get constant value for the deceleration parameter for the mean scale factor as

, for

, for

The sign of indicates whether the model accelerates or not. The positive sign of (i.e.) corresponds to decelerating models whereas the negative sign of for indicates acceleration and for corresponds to expansion with constant velocity.

Using, Equation (20), we obtain the law of average scale factor as

, for

. for

where, and are constants of integration.

Case (i): Model for: We consider that, the scalar expansion for is proportional to the shear scalar for i.e. which leads to

where m is a positive constant.

From Equations (14) and (24), we get

Thus, the Bianchi type-IX magnetized anisotropic DE model in Saez-Ballester scalar-tensor theory can be written as

The average Hubble’s parameter, expansion scalar, anisotropic parameter of the expansion, shear scalar for the model (29) are given by

From Equations (31) and (33), we get

Using Equations (13) and (27), we obtain

Using Equations (27) and (28), Equation (10) takes the form

Using Equations (7), (27), (28), (35) and (36), we obtain energy density for fluid as

It is observed that the Hubble parameter, expansion scalar, mean anisotropic parameter of the expansion, shear scalar, magnetized DE density and energy density are decreasing functions of time and approaches to 0 as.

Since, constant, the model is not isotropic for large values of t.

Using Equations (8), (27), (28), (35), (36) and (37), the equation of state parameter is obtained as

It is observed that the equation of state is time dependent, it can be function of red shift or scale factor R as well.

Using Equations (9), (27), (28), (35), (36) and (37), the skew-ness parameter is given by

In the absence of magnetic field i.e., the values of the Hubble parameter, expansion scalar, mean anisotropic parameter of the expansion, shear scalar remains as it is whereas magnetized DE density and energy density, the EoS parameter, and the skewness parameter are given by

Case (ii): Model for:

Using Equations (14) (25) and (26), we get

Thus, the Bianch Type-IX magnetized anisotropic DE model in Saez-Ballester scalar-tensor theory for can be written as

The expression for kinematical parameters i.e. the average Hubble’s parameter, expansion scalar, anisotropic parameter of the expansion, shear scalar for the model (46) are given by

From Equations (48) and (49),

Using Equations (13) and (45), we obtain

Using Equations (44) and (45), Equation (10) takes the form

Using Equations (7), (44), (45), (52) and (53), we obtain energy density for fluid as

It is observed that the Hubble parameter, expansion scalar, mean anisotropic parameter of the expansion, shear scalar, magnetized DE density and energy density is decreasing  function of time and approaches 0 as.

Since, the model is not isotropic for large values of t.

Using Equations (8), (44), (45), (52), (53) and (54), the equation of state parameter is obtained as

It is observed that the equation of state is time dependent, it can be function of red shift z or scale factor R as well.

Using Equations (9), (44), (45), (52), (53) and (54), the skew-ness parameter is given by

In the absence of magnetic field i.e. , the values of the Hubble parameter, expansion scalar, mean anisotropic parameter of the expansion, shear scalar remains as it is whereas magnetized DE density and energy density, the EoS parameter, and the skew-ness parameter are given by

In this paper, we have obtained Bianchi type-IX cosmological model with variable EoS parameter in Saez-Ballester theory of gravitation in the presence and absence of magnetic field of energy density. In this model, the magnetic field used is as in King and Coles [30] . The solution of the field equations is obtained by using special law of variation for Hubble’s parameter proposed by Bermann [66] . We observed that, the model is free from big-bang singularity and the values of Hubble parameter, expansion scalar, shear scalar are constant at the initial epoch and decreases with time approaching to zero as. implies the model does not approach to isotropy i.e. the model is anisotropic throughout the evolution. The values of Hubble parameter, expansion scalar, mean anisotropic parameter of the expansion and the shear scalar remain same in the presence and absence of magnetic field while the component of magnetic field reduces energy density of the anisotropic fluid. It is noted that our investigations resemble with the investigations of Katore [40] .