ISRN Astronomy and AstrophysicsVolume 2012 (2012), Article ID 420938, 14 pagesdoi:10.5402/2012/420938

Research Article

Effect of Radiative Heat-Loss Function and Finite Larmor Radius Corrections on Jeans Instability of Viscous Thermally Conducting Self-Gravitating Astrophysical Plasma

Sachin Kaothekar1,2 and R. K. Chhajlani¹

¹School of Studies in Physics, Vikram University, Madhya Pradesh, Ujjain 456010, India²Department of Physics, Mahakal Institute of Technology, Madhya Pradesh, Ujjain 456664, India

Received 3 April 2012; Accepted 15 June 2012

Academic Editors: C. W. Engelbracht, A. Ferrari, F. Fraschetti, I. Goldman, and C. Meegan

Copyright © 2012 Sachin Kaothekar and R. K. Chhajlani. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The effect of radiative heat-loss function and finite ion Larmor radius (FLR) corrections on the self-gravitational instability of infinite homogeneous viscous plasma has been investigated incorporating the effects of thermal conductivity and finite electrical resistivity for the formation of a star in astrophysical plasma. The general dispersion relation is derived using the normal mode analysis method with the help of relevant linearized perturbation equations of the problem. Furthermore the wave propagation along and perpendicular to the direction of external magnetic field has been discussed. Stability of the medium is discussed by applying Routh Hurwitz’s criterion. We find that the presence of radiative heat-loss function and thermal conductivity modify the fundamental Jeans criterion of gravitational instability into radiative instability criterion. From the curves we see that temperature dependent heat-loss function, FLR corrections and viscosity have stabilizing effect, while density dependent heat-loss function has destabilizing effect on the growth rate of self-gravitational instability. Our result shows that the FLR corrections and radiative heat-loss functions affect the star formation.

1. Introduction

The problem of self-gravitational instability is widely investigated due to its relevance to the fragmentation of interstellar medium and its role in star formation. Also, the self-gravitational instability of molecular clouds is connected to the cloud collapse and star formation. Hayashi [1] has discussed the problem of evolution of protostars and discusses the different phase of formation of protostar in connection with variation of temperature and density. Shu et al. [2] have investigated the problem of star formation in molecular clouds and concluded that the star formation occurs mainly in four phases. In this problem, the self-gravitational instability of dust and gas plays important role. Draine and Mckee [3] have studied the theory of interstellar shocks. Mckee and Ostriker [4] have discussed the theory of star formation. They concluded that the key dynamical processes involved in star formation are turbulence, magnetic fields, and self-gravity. In this connection, the gravitational instability of infinite homogeneous self-gravitating magnetized and rotating plasma is also discussed by Chandrasekhar [5]. Several authors (Pacholczyk and Stodolkiewicz [6], Nayyar [7], and Shaikh et al. [8]) have investigated the problem of gravitational instability of plasma with different physical parameters such as viscosity, finite electrical conductivity, Hall current, thermal conductivity, magnetic field, and rotation. Yang et al. [9] have investigated the problem of large-scale gravitational instability and star formation. Borah and Sen [10] have investigated the gravitational instability of partially ionized molecular clouds considering the effects of electrons, ions, and charged dust grains. Avinash et al. [11] have studied the dynamics of self-gravitating dust clouds and the formation of planetesimals. Thus we find that a large number of problems are discussed for self-gravitating dusty and nondusty plasma with different parameters under various assumptions due to its importance in star formation and in many astrophysical situations.

Along with this, in above discussed problems, the effect of finite ion Larmor radius is not considered. In many astrophysical situations such as in interstellar and interplanetary plasmas the approximation of zero Larmor radius is not valid. Several authors (Roberts and Taylor [12], Jeffery and Taniuti [13], Jukes [14], and Vandakurov [15]) have pointed out the importance of finite ion Larmor radius (FLR) effects in the form of magnetic resistivity, on the plasma instability. Sharma [16] has shown the stabilizing effect of FLR on gravitational instability of rotating plasma. Bhatia and Chhonkar [17] have investigated the stabilizing effect of FLR on the instability of a rotating and self-gravitating plasma. Herrnegger [18] has studied the effects of collision and gyroviscosity on gravitational instability in a two-component plasma and concluded that the critical wave number becomes smaller with increasing gyroviscosity for finite Alfven numbers and showed that Jeans criterion is changed by FLR for wave propagating perpendicular to magnetic field. Vaghela and Chhajlani [19] have investigated the stabilizing effect of FLR on magneto-gravitational stability of resistive plasma through porous medium with thermal conduction, but they neglect the effect of radiative heat-loss function on gravitational instability. Chhajlani and Parihar [20] have carried out the stabilizing effect of FLR on magnetogravitational instability of anisotropic plasma with generalized polytrope laws. Ferraro [21] has shown the stabilizing effect of FLR on magneto-rotational instability. Recently Sandberg et al. [22] have investigated the stabilizing effect of FLR on the coupled trapped electron and ion temperature-gradient modes. More recently Devlen and Pekünlü [23] have studied the effect of FLR on weakly magnetized, dilute plasma. Thus FLR effect is an important factor in discussion of self-gravitational instability and other hydrodynamic instability.

In addition to this it is well known that thermal and radiative effects do play an important role in the stability investigations. The thermal instability arising due to various heat-loss mechanisms may be the cause of astrophysical condensation and the formation of large and small objects. Several authors (Field [24], Hunter [25], Ibanez [26], Kim and Narayan [27], Radwan [28], Menou et al. [29], and Inutsuka et al. [30]) have investigated the phenomenon of thermal instability arising due to heat-loss mechanism in plasma. Shadmehri and Dib [31] have discussed the thermal instability in magnetized partially ionized plasma with charged dust particles and radiative cooling function. Shaikh et al. [32] have investigated the Jeans gravitational instability of thermally conducting plasma in a variable magnetic field with Hall current, finite conductivity, and viscosity, but they neglect the effect of FLR corrections and radiative heat-loss function on gravitational instability. Aggarwal and Talwar [33] have discussed magnetothermal instability in a rotating gravitating fluid, but they neglect the effect of FLR correction on radiative instability. Bora and Talwar [34] have investigated the magnetothermal instability with finite electrical resistivity and Hall current, but they neglect the effect of FLR corrections and viscosity on radiative instability. Recently El-Sayed and Mohamed [35] have discussed the gravitational instability of rotating viscoelastic partially ionized plasma in the presence of oblique magnetic field and Hall current. More recently Kaothekar and Chhajlani [36] have carried out the problem of gravitational instability of radiative plasma with FLR corrections.

In the light of above work, we find that all these authors have studied the problem of gravitational instability with different combinations of these parameters, but none studied the joint effect of all the parameters together with radiative effects. We also find that none of the above authors has tried to explore Jeans condition for viscous and nonviscous plasma with FLR corrections, radiative heat-loss function, and thermal conductivity. The problem given by Kaothekar and Chhajlani [36] is briefly discussed with permeability parameter, but in present paper we have made a detailed analysis of the problem of self-gravitational instability with all the considered physical parameters, without including the effect of permeability to give a better insight in real problem. Therefore in the present work self-gravitational instability of magnetized plasma with FLR corrections, radiative heat-loss function, viscosity, thermal conductivity, and finite electrical resistivity for self-gravitating configuration is studied. We also wish to explore the importance of viscous and nonviscous system and its impact on the self-gravitational instability of plasma in connection with FLR correction, radiative heat-loss function, thermal conductivity, and finite electrical resistivity. The stability of the system is discussed by applying Routh-Hurwitz criterion. The above work is applicable to formation of stars in astrophysical plasma.

2. Basic Equations of the Problem and Perturbation

Let us consider an infinite homogeneous, self-gravitating, radiating, thermally conducting, and viscous plasma of finite electrical resistivity, including the effect of finite ion Larmor radius (FLR) corrections in the presence of magnetic field  . The equations of the problem with these effects are written as where  , and denote the fluid pressure, density, kinematic viscosity, temperature, electrical resistivity, thermal conductivity, gravitational potential, gravitational constant, gas constant, and ratio of two specific heats respectively. , the radiative heat-loss function, depends on local values of density and temperature of the fluid. The operator is the substantial derivative given as  .

The components of pressure tensor  , considering the finite ion gyration radius for the magnetic field along  -axis as given by Roberts and Taylor [12], are The parameter has the dimensions of the kinematic viscosity and defined as , where is the ion-Larmor radius, and is the ion gyration frequency. The perturbation in fluid pressure, density, temperature, velocity, magnetic field, heat-loss function, and gravitational potential is given as , , , , , , and , respectively. The perturbation state is given as Suffix “0” represents the initial equilibrium state, which is independent of space and time.

3. Dispersion Relation

We linearize the basic set of MHD equations by substituting (3) in (1)-(2), and we assume the variation of perturbed quantities as where is the frequency, and and are the wave numbers of the perturbations along and axes. Using (4) in linearized perturbed form of (1) with (2), we get four equations governing the perturbations of velocity and the condensation of medium. These four equations can be written in the following determinant form: We have made following substitutions: where   is the adiabatic velocity of sound in the medium,   is the isothermal velocity of sound in the medium, and is the condensation of the medium.

The general dispersion relation can be obtained from the determinant of matrix of (5) is

The dispersion relation (7) represents the simultaneous inclusion of FLR corrections, radiative heat-loss function, thermal conductivity, viscosity, magnetic field, and finite electrical conductivity on self-gravitational instability of plasma. In absence of radiative heat-loss function, the general dispersion relation (7) is identical to that of Vaghela and Chhajlani [19] for and in their case. In absence of radiative heat-loss function, thermal conductivity, finite electrical resistivity, and viscosity the general dispersion relation (7) is identical to Sharma [16] for nonrotational case. In absence of FLR corrections and radiative heat-loss function the general dispersion relation (7) is identical to Shaikh et al. [32] neglecting Hall current in their case. In absence of FLR corrections and viscosity dispersion relation (7) is identical to Bora and Talwar [34] neglecting Hall current and electron inertia in their case. Also in absence of FLR corrections, viscosity and finite conductivity dispersion relation (7) reduces to that obtained by Field [24] for nongravitating medium. The dispersion relation (7) is different from (9) of Kaothekar and Chhajlani [36].

Thus we have obtained the modified dispersion relation of self-gravitational instability including the combined effect of FLR corrections, radiative heat-loss function, viscosity, thermal conductivity, finite electrical resistivity, and magnetic field. Now we discuss the general dispersion relation (7) for longitudinal and transverse wave propagation.

4. Discussion

4.1. Longitudinal Propagation

For this case, we assume all the perturbations are longitudinal to the direction of magnetic field (i.e.,  ,  ). The dispersion relation (7) reduces to

The above dispersion relation shows the combined effect of FLR corrections, radiative heat-loss function, thermal conductivity, viscosity, magnetic field strength, and finite electrical resistivity on self-gravitation instability of plasma. On multiplying all the components of (8), we get the dispersion relation, which is an equation of degree eight in , and it is cumbersome to write such a lengthy equation. If we remove the effect of FLR corrections and viscosity in the above relation, then we recover the relation given by Bora and Talwar [34] excluding Hall current and electron inertia in their case. Hence the above dispersion relation is the modified form of equation of Bora and Talwar [34] due to the inclusion of FLR corrections and viscosity in our case and by neglecting Hall current and electron inertia in their case for longitudinal propagation in dimensional form. But the condition of instability is unaffected by the presence of FLR correction and viscosity. Thus we conclude that FLR corrections and viscosity have no effect on the condition of instability, but presence of these parameters modifies the growth rate of instability in the present case. Hence these are the new findings in our case than that of Bora and Talwar [34]. Also on comparing our dispersion relation (8) with dispersion relation of Vaghela and Chhajlani [19] we find that two factors are the same, but the third factor is different and gets modified because of radiative terms.

The dispersion relation (8) has three different components, and our aim is to take out the physics involved in each component. So we discuss each component separately. The first component of the dispersion relation (8) gives

This represents a damped mode modified by the presence of viscosity of the medium. Thus viscous force is capable of stabilizing the growth rate of the considered system. The above mode is unaffected by the presence of FLR correction, magnetic field strength, finite electrical resistivity, thermal conductivity, radiative heat-loss function, and self-gravitation.

The second factor of (8) on simplification gives

The above dispersion relation represents the Alfven mode modified by the presence of FLR corrections, viscosity, and finite electrical resistivity. But it is independent of thermal conductivity, radiative heat-loss function, and self-gravitation. Equation (10) is a four-degree equation in power of having its all coefficients positive which is a necessary condition for the stability of the system. To achieve the sufficient condition the principal diagonal minors of Hurwitz matrix must be positive. On calculating we get all the principal diagonal minors positive. Hence (10) always represents stability.

The third component of the dispersion relation (8) on simplifying gives

The above equation represents the combined influence of radiative heat-loss function, thermal conductivity, and viscosity on the self-gravitational instability of plasma. But there is no effect of FLR corrections, finite electrical conductivity, and magnetic field on the self-gravitational instability of the considered system. If the constant term of cubic equation (11) is less than zero, this allows at least one positive real root which corresponds to the instability of the system. The condition of instability obtained from constant term of (11) is given as This condition of instability is modified form of Jeans condition by inclusion of thermal conductivity and radiative heat-loss functions. The above condition of instability is independent of FLR corrections, finite electrical conductivity, magnetic field strength, and viscosity. The above inequality (12) is same as obtained by Bora and Talwar [34], and by Kaothekar and Chhajlani [36]. But Kaothekar and Chhajlani [36] have also considered the effect of permeability in their analysis. Here in the present problem we neglect the effect of permeability to give better insight in real problem. Thus in this mode the dispersion relation and the growth rate is modified due to FLR corrections and viscosity, but the condition of radiative instability is unaffected by the presence of FLR corrections and viscosity. The condition of instability (12) can be solved to yield the following modified critical Jeans wavelength:

The medium is unstable for wave length . Here it may be noted that the modified critical wave length involves the derivatives of temperature-dependent, density-dependent heat-loss functions and thermal conductivity of the medium.

In dispersion relation (8) first and second factors show wave propagation, but third factor shows instability. Thus to discuss the effect of each parameter on the growth rate of instability, we solve (11) numerically by introducing the following dimensionless quantities: The nondimensional form of (11) in terms of self-gravitation is given as

Here for calculations we take the numerical values corresponding to the conditions in interstellar molecular clouds for star formation as

And the parameters are taken as nondimensional. Numerical calculations were performed to determine the roots of from dispersion relation (15), as a function of wave number for several values of the different parameters involved taking . Out of the three modes only one mode is unstable for which the calculations are presented in Figures 1–3, where the growth rate (positive real value of , after multiplying by ) has been plotted against the wave number (after multiplying by ) to show the dependence of the growth rate on the different physical parameters such as temperature-dependent heat-loss function (after multiplying by), density-dependent heat-loss function (after multiplying by), and viscosity (after multiplying by ).

Figure 1: Growth rate (positive values of ) against wave number for three values of parameter , keeping the other parameters fixed ,  ,  .

Figure 2: Growth rate (positive values of ) against wave number for three values of parameter , keeping the other parameters fixed ,  ,  .

Figure 3: Growth rate (positive values of ) against wave number for three values of parameter , keeping the other parameters fixed ,  ,  .

It is clear from Figure 1 that growth rate decreases with increasing temperature-dependent heat-loss function. Thus the effect of temperature-dependent heat-loss function is stabilizing. From Figure 2 we conclude that growth rate increases with increasing density-dependent heat-loss function. Thus the effect of density-dependent heat-loss function is destabilizing. From Figure 3 we conclude that growth rate decreases with increasing viscosity of the medium. Thus the effect of viscosity is stabilizing.

To discuss the stability of the system, if constant term of cubic equation (11) is greater than zero, then all the coefficients of (11) must be positive. Equation (11) is a third-degree equation in the power of having its coefficients positive, which is a necessary condition for the stability of the system. To achieve the sufficient condition the principal diagonal minors of Hurwitz matrix must be positive. The principal diagonal minors are Since ,   and , it is clear that all the    are positive hence system represented by (11) is stable system.

Now we wish to examine the effect of radiative heat-loss function in the considered system with some simplifications, and at the same time we wish to investigate the physics involved in such simplifications in the present problem.

4.1.1. For Viscous Medium

(1)Effect of Radiative Heat-Loss Function

For nonradiative medium (11) becomes The condition of instability from constant term of (18) is , or if

The medium is unstable for wave length , where is the modified Jeans wave length for thermally conducting medium. On comparing (18) and (11) we see that there is no new mode coming due to inclusion of radiative heat-loss functions, but the condition of instability gets modified in form of (11). Also it is interesting to note that the thermal conductivity term is separated from condition of instability when we neglect the effect of heat-loss function in (11).

For nonthermal and nonradiating medium the condition of instability obtained from (11) is

The medium is unstable for wave length , where is the Jeans wave length. The condition of Jeans instability given by (21) is identical to Chandrasekhar [5]. On comparing (12) and (21) it is clear that due to the presence of radiative heat-loss functions and thermal conductivity the fundamental criterion of Jeans gravitational instability changes into radiative instability criterion. Also, we conclude that one mode is increased because of inclusion of radiative heat-loss functions.

4.1.2. For Nonviscous Medium

(1) Effect of Radiative Heat-Loss Function

For inviscid and nonradiating medium (11) becomes The condition of instability is , it is discussed in (19) and (20). On comparing (23) and (11) we see that there is no new mode which comes due to inclusion of radiative heat-loss functions, but the condition of instability gets modified in form of (12). In both the cases whether the medium is viscous or nonviscous, the condition of instability is the same; hence, we conclude that viscosity of the medium has no effect on the condition of instability, and it only modifies the growth rate of the instability for viscous case.

For inviscid, thermally nonconducting and nonradiating medium the condition of instability obtained from (11) is , and it is discussed in (21) and (22). We see that inclusion of viscosity cannot change the condition of instability, but it modifies the growth rate of instability.

Thus we conclude that for longitudinal wave propagation as given by (8) the system is unstable only for modified Jeans condition, and else it is stable. Also the modified Jeans criterion remains unaffected by FLR corrections, viscosity, magnetic field, and finite electrical resistivity, but radiative heat-loss function and thermal conductivity modifies the Jeans expression and the fundamental Jeans instability criterion becomes radiative instability criterion. The growth rate of the system is modified by the presence of radiative heat-loss functions, thermal conductivity, and viscosity. The growth rate of Alfven mode is modified by the presence of viscosity, finite electrical resistivity, and FLR corrections. From the curves we find that temperature-dependent heat-loss function and viscosity have stabilizing influence, whereas density-dependent heat-loss function has a destabilizing influence on the self-gravitational instability of plasma.

4.2. Transverse Propagation

For this case, we assume that all the perturbations are transverse to the direction of magnetic field (i.e.,  ). The dispersion relation (7) reduces to

We find that for transverse mode of propagation the dispersion relation is modified due to the presence of FLR corrections, radiative heat-loss function, thermal conductivity, viscosity, finite electrical resistivity, and magnetic field. The dispersion relation (24) has two different components. The first component of the dispersion relation (24) represents a damped mode modified by the presence of viscosity of the medium and is discussed in (9).

The second component of the dispersion relation (24) on simplifying gives

The above fifth degree equation in represents the combined influence of FLR corrections, radiative heat-loss function, thermal conductivity, finite electrical resistivity, viscosity, and magnetic field on self-gravitational instability of plasma. When constant term of (25) is less than zero, this allows at least one positive real root which corresponds to the instability of the system. The condition of instability obtained from constant term of (25) is given as The above condition of instability is independent of FLR correction, finite electrical resistivity, viscosity, and magnetic field strength and is discussed in (12). Thus from (25), we find that the modified Jeans condition of instability remains the same whether we consider the effect of FLR corrections or not, but, owing to the presence of FLR corrections, the growth rate of radiative instability is modified. This result can be easily explained by graphical representation as shown in Figure 6; by increasing values of FLR, we see that the growth rate of radiative instability decreases. Hence FLR shows stabilizing influence on the growth rate of radiative instability.

In several astrophysical plasma situations, such as in interstellar molecular clouds the formation of objects is mainly due to the unstable modes. Thus to discuss the effect of each parameter on the growth rate of unstable modes, we solve (25) numerically by using the relation given in (14) with the following additional dimensionless quantities:

For calculations we take the numerical values corresponding to the conditions in interstellar molecular clouds which are shown in (16), and the rest of values are taken as nondimensional. Numerical calculations were performed to determine the roots of as a function of wave number for several values of the different parameters involved taking . Out of the five modes, only one mode is unstable for which the calculations are presented in Figures 4–6, where the growth rate (positive real value of , after multiplying by ) has been plotted against the wave number (after multiplying by ) to show the dependence of the growth rate on the different physical parameters such as temperature-dependent heat-loss function (after multiplying by), density-dependent heat-loss function (after multiplying by ), viscosity (after multiplying by ), and FLR corrections (after multiplying by ).

Figure 4: Growth rate (positive values of ) against wave number for three values of parameter , keeping the other parameters fixed ,  ,  ,  ,  ,  .

Figure 5: Growth rate (positive values of ) against wave number for three values of parameter , keeping the other parameters fixed ,  ,  ,  ,  ,  .

Figure 6: Growth rate (positive values of ) against wave number   for several values of parameter   keeping the other parameters fixed ,  ,  ,  ,  .  :   , (2) (3) .  : (4) (5) .

It is clear from Figure 4 that growth rate decreases with increasing temperature-dependent heat-loss function. Thus the effect of temperature-dependent heat-loss function is stabilizing. From Figure 5 we conclude that growth rate increases with increasing density-dependent heat-loss function. Thus the effect of density-dependent heat-loss function is destabilizing. One can observe from Figure 6 that the growth rate decreases with increasing FLR corrections. Thus the effect of FLR corrections is stabilizing. If the viscous effect of plasma is removed, then the stabilizing rate of FLR corrections increases and the peak value of the curve also increases. Hence we conclude that both FLR corrections and viscosity stabilize the system.

To discuss the stability of the system given by dispersion relation (25), if constant term of equation (25) is greater than zero, then all the coefficients of the equation (25) must be positive. Equation (25) is a fifth degree equation in the power of ; hence, it will involve a lengthy algebra to calculate the principal diagonal minors of Hurwitz matrix. Since (25) has positive real coefficients, having ,   and , so it is sufficient to show that alternate    are positive; hence, system represented by (25) is a stable system.

Now we wish to examine the effect of FLR corrections and radiative heat-loss function on the considered system with some simplifications such as viscous and nonviscous medium, and at the same time we wish to investigate the physics involved in such simplifications in the present problem.

4.2.1. For Viscous Medium

(1)Effect of Radiative Heat-Loss Function with FLR Corrections

For nonradiative medium, but in the presence of FLR corrections the condition of instability obtained from the constant term of (25) is given as

The above (28) is identical to that of Vaghela and Chhajlani [19]. We see that there is no new mode due to inclusion of radiative heat-loss function, but the growth rate of instability gets modified due to inclusion of radiative heat-loss functions and FLR corrections. Also it is interesting to note that, on neglecting the effect of radiative heat-loss function in (25), the thermal conductivity parameter (i.e., ) is separated from instability condition as shown in (28). So we conclude that on neglecting the effect of radiative heat-loss function the condition of instability is independent of both the radiative heat-loss function and thermal conductivity of the medium, but the growth rate of instability is affected by thermal conductivity and FLR corrections.

For nonradiative and nonthermal medium, but in presence of FLR corrections (,  ) the condition of instability obtained from (25) is given as

We find that one mode is increased due to inclusion of radiative heat-loss functions. On comparing (26) and (29), it is clear that owing to the presence of thermal conductivity and radiative heat-loss functions the fundamental criterion of Jeans gravitational instability changes into radiative instability criterion. In case of transverse mode of propagation, for viscous medium, FLR corrections have no effect on the condition of instability for radiating or nonradiating medium, but the growth rate of instability is modified by FLR corrections; that is, FLR stabilizes the system in all the cases. The joint stabilizing effect of FLR corrections and viscosity decreases the peak value of growth rate of the system shown in Figure 6.

4.2.2. For Nonviscous Medium

(1) Effect of FLR Corrections

For inviscid medium, but in presence of FLR corrections (25) becomes The above equation (30) is a modified form of Chhajlani and Parihar [20] by inclusion of radiative heat-loss function, thermal conductivity, and finite electrical resistivity. The condition of instability obtained from constant term of (30) is given as From the above condition of instability given by (31) we conclude that FLR correction tries to stabilize the radiative instability. Also on comparing (25) and (30) we see that inclusion of viscosity removes the effects of FLR correction from condition of instability. So in both the cases either the system is viscous or nonviscous, FLR correction stabilizes the growth rate of radiative instability, as shown in Figure 6. Thus in presence of FLR correction the condition of instability given by Bora and Talwar [34] is modified as given by (31). In absence of FLR corrections (30) is identical to Bora and Talwar [34]; thus, the dispersion relation is modified due to the inclusion of effects of FLR in the present system. Thus we find that due to the presence of FLR corrections Jeans criterion of instability is modified, and also the growth rate of instability in this case is modified with the inclusion of FLR correction in the present MHD set of equations. The above condition of instability (31) can be written in nondimensional form as

In absence of FLR corrections the above condition of instability is the same as given by Bora and Talwar [34] and by Hunter [25] for an electrically nonconducting gas. We conclude that condition of instability given by Bora and Talwar [34] is modified by inclusion of FLR corrections in our case. Also the condition of instability given by (31) is the modified form of Field [24] by inclusion of FLR correction and self-gravitation. Thus the present result is the improvement of Field [24], Hunter [25], and Bora and Talwar [34]. Condition (32) involves the derivatives and of the heat-loss function with respect of local temperature and density, and is FLR correction factor in the system. It may be noted that for a density-independent heat-loss function , which increases with temperature , condition (32) gives instability if . However if instead, the heat-loss function decreases with temperature , and the instability comes to lying between and .

For a purely density-dependent heat-loss function , condition (32) gives instability if The above equation (33) applies for or if the heat-loss function, respectively, increases or decreases with increase in density.(2) Effect of Radiative Heat-Loss Function with FLR Corrections

For nonradiating and inviscid medium, but in presence of FLR the condition of instability obtained from (25) is given as The critical Jeans wave number and critical Jeans wave length are given as Or

Condition (34) is identical of Vaghela and Chhajlani [19] for in their case. From (34) we conclude that FLR have stabilizing effect on the instability of the system. We find that one mode is increased due to inclusion of viscosity, and we conclude that inclusion of viscosity removes the effect of FLR corrections from the condition of instability. Also on comparing (31) and (34), we conclude that condition of instability given by Vaghela and Chhajlani [19] is modified by inclusion of radiative effects and the growth rate of instability given by (31) is modified by radiative heat-loss function and FLR corrections.

For inviscid, nonradiating and thermally nonconducting medium, but in presence of FLR , the condition of instability obtained from (25) is given as The critical Jeans wave number or critical Jeans wave length is given as Or

From (37) we find that Jeans condition of gravitational instability gets modified by FLR corrections. On comparing (37) and (31) we conclude that inclusion of thermal conductivity and radiative heat-loss function modifies both the Jeans condition of instability and growth rate of instability in presence of FLR corrections. Hence the fundamental Jeans criterion of gravitational instability gets modified into radiative instability criterion by inclusion of thermal conductivity and radiative heat-loss functions in the presence of FLR corrections.(3) Effect of Electrical Conductivity with FLR Corrections

For inviscid and infinitely conducting medium, but in presence of FLR (), the condition of instability obtained from constant term of (25) is given as

From the above condition of instability given by (40) we conclude that FLR correction and magnetic field strength try to stabilize the radiative instability. Also on comparing (26) and (40) we see that inclusion of viscosity and finite conductivity removes the effect of FLR correction and magnetic field from condition of instability. So in both the cases whether the system is finitely conducting or infinitely conducting magnetic field and FLR correction stabilizes the growth rate of radiative instability. The above condition of instability is identical to Kaothekar and Chhajlani [36]. Kaothekar and Chhajlani [36] have considered the effect of permeability, but, in the present problem, we have not considered the effect of permeability, so we conclude that the conditions of instability for permeable or nonpermeable medium are same, but the growth rate of the system is modified by the presence of permeability. The above condition of instability (40) can be written in nondimensional form as In absence of FLR corrections the above condition of instability is identical to Bora and Talwar [34] excluding electron inertia in that case. We conclude that condition of instability given by Bora and Talwar [34] is modified by inclusion of FLR corrections in our case. Thus the present results are the improvement of Bora and Talwar [34] due to FLR corrections.

From condition of instability (41) it may be noted that for a density-independent heat-loss function , which increases with temperature , it gives instability if . However if instead, the heat-loss function decreases with temperature , the instability comes to lying between and .

For a purely density-dependent heat-loss function , condition (41) gives instability if The above equation (42) applies for or , if the heat-loss function, respectively, increases or decreases with increase in density.(4) Effect of Nongravitation and FLR Corrections

For inviscid, nongravitating, infinitely conducting, radiating, thermally conducting, magnetized medium with FLR corrections (25) becomes The condition of instability obtained from constant term of (43) is given as

From above equation (44) we see that FLR correction and magnetic field strength try to stabilize the system. If we neglect the effect of FLR and magnetic field the above dispersion relation (43) is identical to Field [24]. In the present case we have considered the effects of FLR correction and magnetic field, but Field [24] has not considered these effects. Thus in the present analysis both the dispersion relation and the condition of instability get modified due to the presence of FLR corrections and magnetic field strength. On comparing (40) and (44) we see that consideration of self-gravitation modifies the thermal instability criterion into radiative instability criterion. Also from (25) it is clear that the growth rate of the dispersion relation given by Field [24] is modified due to the presence of FLR correction, viscosity, magnetic field, finite electrical conductivity, and self-gravitation in our present case. Hence these are the new results in our present problem than that of Field [24].

Thus we conclude that for transverse wave propagation the Jeans criterion is affected by FLR corrections, radiative heat-loss functions, thermal conductivity, viscosity, magnetic field strength, and finite electrical resistivity. From curves we find that FLR corrections, temperature-dependent heat-loss function, and viscosity have stabilizing influence, whereas density-dependent heat-loss function has destabilizing influence on the self-gravitational instability of plasma.

5. Conclusion

We have dealt with the self-gravitation instability of an infinite homogeneous viscous electrically and thermally conducting fluid including the effects of FLR corrections and radiative heat-loss function for star formation. The general dispersion relation is obtained, which is modified due to the presence of considered physical parameters. This dispersion relation is reduced for longitudinal and transverse mode of propagation to the direction of the magnetic field. We find that Jeans criterion remains valid and gets modified because of FLR corrections, radiative heat-loss function, thermal conductivity, and magnetic field. We also find that the presence of thermal conductivity and radiative heat-loss function modifies the fundamental Jeans criterion of gravitational instability into radiative instability criterion. The effect of viscosity parameter is found to stabilize the system in both the longitudinal and transverse mode of propagation. For longitudinal wave propagation magnetic field, viscosity, finite electrical resistivity, and FLR correction have no effect on Jeans criterion, but FLR corrections, viscosity, and finite electrical resistivity modify the growth rate of Alfven mode. For transverse wave propagation FLR corrections stabilize the growth rate of the system in all the cases, but it modifies the condition of instability only for the case of nonviscous medium. Also, magnetic field stabilizes the system but finite conductivity removes the effect of magnetic field there by destabilizing the system. Numerical calculation shows stabilizing effect of temperature-dependent heat-loss function, FLR corrections, and viscosity, where as destabilizing effect of density-dependent heat-loss function is on the self-gravitational instability.

Acknowledgments

The authors are highly thankful to Prof. and Head Dr. S. K. Ghosh, S. S. in Physics and Dr. T. R. Thapak, Hon’ble V. C. Vikram University Ujjain (MP) for their constant encouragements. S. Kaothekar is grateful to Er. Praveen Vashishtha, Chairman of Mahakal Institute of Technology for continuous support and guidance.

References

C. Hayashi, “Evolution of protostars,” Annual Review of Astronomy and Astrophysics, vol. 4, pp. 171–192, 1966.

F. H. Shu, F. C. Adams, and S. Lizano, “Star formation in molecular clouds,” Annual Review of Astronomy and Astrophysics, vol. 25, pp. 23–81, 1987.

B. T. Draine and C. F. McKee, “Theory of interstellar shocks,” Annual Review of Astronomy and Astrophysics, vol. 31, pp. 373–432, 1993.

C. F. McKee and E. C. Ostriker, “Theory of star formation,” Annual Review of Astronomy and Astrophysics, vol. 45, pp. 565–687, 2007.

S. Chandrasekhar, Hydrodynamics and Hydromagnetic Stability, Clarendon Press, Oxford, UK, 1961.

A. G. Pacholczyk and J. S. Stodolkiewicz, “On the gravitational instability of some magnetohydrodynamical systems of astrophysical interest,” Acta Astronomica, vol. 10, pp. 1–29, 1960.

N. K. Nayyar, “Magnetogravitational instability of an infinite medium with finite electrical and thermal conductivity,” Zeitschrift fur Astrophysik, vol. 52, pp. 266–271, 1961.

S. Shaikh, A. Khan, and P. K. Bhatia, “Stability of thermally conducting plasma in a variable magnetic field,” Astrophysics and Space Science, vol. 312, pp. 35–40, 2007.

C.-C. Yang, R. A. Gruendl, Y.-H. Chu, M.-M. Mac Low, and Y. Fukui, “Large-scale gravitational instability and star formation in the large magellanic cloud,” The Astrophysical Journal, vol. 671, pp. 374–379, 2007.

A. C. Borah and A. K. Sen, “Gravitational instability of partially ionized molecular clouds,” Journal of Plasma Physics, vol. 73, no. 6, pp. 831–838, 2007.

K. Avinash, B. Eliasson, and P. K. Shukla, “Dynamics of self-gravitating dust clouds and the formation of planetesimals,” Physics Letters, Section A, vol. 353, pp. 105–108, 2006.

K. V. Roberts and J. B. Taylor, “Magnetohydrodynamic equations for finite Larmor radius,” Physical Review Letters, vol. 8, pp. 197–198, 1962.

A. Jeffery and T. Taniuti, MHD Stability and Thermonuclear Containment, Academic Press, New York, NY, USA, 1966.

J. D. Jukes, “Gravitational resistive instabilities in plasma with finite Larmor radius,” Physics of Fluids, vol. 7, no. 1, pp. 52–58, 1964.

L. V. Vandakurov, “On the theory of stability of a plasma in a strong longitudinal magnetic field, variable along the axis of symmetry,” Journal of Applied Mathematics and Mechanics, vol. 28, no. 1, pp. 77–89, 1964.

R. C. Sharma, “Gravitational instability of rotating plasma,” Astrophysics and Space Science, vol. 29, pp. L1–L4, 1974.

P. K. Bhatia and R. P. S. Chhonkar, “Larmor radius effects on the instability of a rotating layer of a self-gravitating plasma,” Astrophysics and Space Science, vol. 115, no. 2, pp. 327–344, 1985.

F. Herrnegger, “Effect of collisions and gyroviscosity on gravitational instability in a two-component plasma,” Journal of Plasma Physics, vol. 8, pp. 393–400, 1972.

D. S. Vaghela and R. K. Chhajlani, “Magnetogravitational stability of resistive plasma through porous medium with thermal conduction and FLR corrections,” Contributions to Plasma Physics, vol. 29, pp. 77–89, 1989.

R. K. Chhajlani and A. K. Parihar, “Magnetogravitational instability of anisotropic plasma with finite ion Larmor and generalized polytrope law,” Contributions to Plasma Physics, vol. 34, pp. 669–681, 1994.

N. M. Ferraro, “Finite Larmor radius effects on the magnetorotational instability,” The Astrophysical Journal, vol. 662, pp. 512–516, 2007.

I. Sandberg, H. Isliker, and V. P. Pavlenko, “Finite Larmor radius effects on the coupled trapped electron and ion temperature gradient modes,” Physics of Plasmas, vol. 14, Article ID 092504, 2007.

E. Devlen and E. R. Pekünlü, “Finite Larmor radius effects on weakly magnetized, dilute plasmas,” Monthly Notices of the Royal Astronomical Society, vol. 404, pp. 830–836, 2010.

G. B. Field, “Thermal instability,” The Astrophysical Journal, vol. 142, pp. 531–567, 1965.

J. H. Hunter, “The role of thermal instabilities in star formation,” Monthly Notices of the Royal Astronomical Society, vol. 133, pp. 239–245, 1966.

M. H. Ibanez S., “Sound and thermal waves in a fluid with an arbitrary heat-loss function,” The Astrophysical Journal, vol. 290, pp. 33–46, 1985.

W.-T. Kim and R. Narayan, “Thermal instability in clusters of galaxies with conduction,” The Astrophysical Journal, vol. 596, pp. 889–902, 2003.

A. E. Radwan, “Variable streams self-gravitating instability of radiating rotating gas cloud,” Applied Mathematics and Computation, vol. 148, no. 2, pp. 331–339, 2004.

K. Menou, S. A. Balbus, and H. C. Spruit, “Local axisymmetric diffusive stability of weakly magnetized, differentially rotating, stratified fluids,” The Astrophysical Journal, vol. 607, no. 1, pp. 564–574, 2004.

S. I. Inutsuka, H. Koyama, and T. Inoue, “The role of thermal instability in interstellar medium,” in Proceedings of the AIP Conference on Magnetic Fields in the Universe: from Laboratory and Stars to Primordial Structures, vol. 784, pp. 318–328, Angra dos Reis, Brazil, December 2004.

M. Shadmehri and S. Dib, “Magnetothermal condensation modes including the effects of charged dust particles,” Monthly Notices of the Royal Astronomical Society, vol. 395, pp. 985–990, 2009.

S. Shaikh, A. Khan, and P. K. Bhatia, “Jeans' gravitational instability of a thermally conducting plasma,” Physics Letters, Section A, vol. 372, no. 9, pp. 1451–1457, 2008.

M. Aggarwal and S. P. Talwar, “Magnetothermal instability in a rotating gravitating fluid,” Monthly Notices of the Royal Astronomical Society, vol. 146, pp. 235–242, 1969.

M. P. Bora and S. P. Talwar, “Magnetothermal instability with generalized Ohm's law,” Physics of Fluids B, vol. 5, no. 3, pp. 950–955, 1993.

M. F. El-Sayed and R. A. Mohamed, “Gravitational instability of rotating viscoelastic partially ionized plasma in the presence of an oblique magnetic field and Hall current,” ISRN Mechanical Engineering, vol. 2011, Article ID 597172, 8 pages, 2011.

S. Kaothekar and R. K. Chhajlani, “Gravitational instability of radiative plasma with finite Larmor radius corrections,” Journal of Physics: Conference Series, vol. 365, no. 1, Article ID 012053, 2012.