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

Authors: Sachin Kaothekar; R. K. Chhajlani
Journal: ISRN Astronomy and Astrophysics
 (2012)

---

## 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.


---

## Body

## 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 fieldH(0,0,H). The equations of the problem with these effects are written as
(1)dudt=-1ρ∇p+∇U-∇·Pρ+14πρ(∇×H)×H+υ∇2u,1γ-1dpdt-γγ-1pρdρdt+ρL-∇·(λ∇T)=0,p=ρRT,dρdt+ρ∇·u=0,∇2U=-4πGρ,∂H∂t=∇×(u×H)+η∇2H,∇·H=0,
where  p,ρ,υ,T,η,λ,U,G,R, 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. L(ρ,T), the radiative heat-loss function, depends on local values of density and temperature of the fluid. The operator (d/dt) is the substantial derivative given as  (d/dt)=(∂/∂t+u·∇).The components of pressure tensorP, considering the finite ion gyration radius for the magnetic field along  z-axis as given by Roberts and Taylor [12], are
(2)Pxx=-ρυ0(∂uy∂x+∂ux∂y),Pyy=ρυ0(∂uy∂x+∂ux∂y),Pzz=0,Pxy=Pyx=ρυ0(∂ux∂x-∂uy∂y),Pxz=Pzx=-2ρυ0(∂uy∂z+∂uz∂y),Pyz=Pzy=2ρυ0(∂uz∂x+∂ux∂z).
The parameter υ0 has the dimensions of the kinematic viscosity and defined as υ0=ΩLRL2/4, where RL is the ion-Larmor radius, and ΩL is the ion gyration frequency. The perturbation in fluid pressure, density, temperature, velocity, magnetic field, heat-loss function, and gravitational potential is given as δp, δρ, δT, u(ux,uy,uz), h(hx,hy,hz), L, and δU, respectively. The perturbation state is given as
(3)p=p0+δp,ρ=ρ0+δρ,T=T0+δT,u=u0+u(withu0=0),H=H0+h,L=L0+L(withL0=0),U=U0+δU.
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
(4)exp(ikxx+ikzz+iσt),
where σ is the frequency, and kx and kz are the wave numbers of the perturbations along x and z 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:
(5)[PF0ikxk2ΩT2-FQ-2υ0kxkz002υ0kxkzMikzk2ΩT2ikxV2k2dikxυ0(kx2+4kz2)0-R]×[uxuyuzs]=0.
We have made following substitutions:
(6)iσ=ω,R=ω2+ωυk2+ΩT2,F=υ0(kx2+2kz2),d=(ω+ηk2),M=ω+υk2,P=M+V2k2d,Q=M+V2kz2d,Lρ=(∂L∂ρ)T,V2=H24πρ,ΩT2=ΩI2+ωΩj2ω+B,Ωj2=c2k2-4πGρ,ΩI2=k2A-4πGρB,LT=(∂L∂T)ρ,A=(γ-1)(TLT-ρLρ+λk2Tρ),B=(γ-1)(TρLTp+λk2Tp),
where c=(γp/ρ)1/2  is the adiabatic velocity of sound in the medium, c'=(p/ρ)1/2  is the isothermal velocity of sound in the medium, and s=δρ/ρ is the condensation of the medium.The general dispersion relation can be obtained from the determinant of matrix of (5) is (7)[ω+υk2+V2k2ω+ηk2](ω2+ωυk2+ΩI2+ωΩj2ω+B)×{4υ02kx2kz2+(ω+υk2)×[ω+υk2V2kz2ω+ηk2]}-(kx2+4kz2)[ω+υk2+V2k2ω+ηk2]×(2υ02kx2kz2k2)(ΩI2+ωΩj2ω+B)+υ02(kx2+2kz2)2[ω2+ωυk2+ΩI2+ωΩj2ω+B](ω+υk2)+(ΩI2+ωΩj2ω+B)(V2ω+ηk2)×2υ02kx2kz2(kx2+2kz2)-(ω+υk2)υ02kx2k2(ΩI2+ωΩj2ω+B)(kx2+2kz2)(kx2+4kz2)-(ω+υk2)×[ω+υk2+V2kz2ω+ηk2](V2kx2ω+ηk2)(ΩI2+ωΩj2ω+B)-4υ02kx4kz2(V2ω+ηk2)(ΩI2+ωΩj2ω+B)=0.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 (ε=1 and K1=∞) 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(kx=0,kz=k)

For this case, we assume all the perturbations are longitudinal to the direction of magnetic field (i.e.,kz=k,  kx=0). The dispersion relation (7) reduces to(8)(ω+υk2)×{4υ02k4+[ω+υk2+V2k2ω+ηk2]2}×[(γ-1){k2(TLT-ρLρ+λk2T/p)-4πGρ(TρLT/p+λk2T/p)}ω+(γ-1)(TρLT/p+λk2T/p)ω2+ωυk2+ω(c2k2-4πGρ)ω+(γ-1)(TρLT/p+λk2T/p)+(γ-1){k2(TLT-ρLρ+λk2T/p)-4πGρ(TρLT/p+λk2T/p)}ω+(γ-1)(TρLT/p+λk2T/p)]=0.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 (21) 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 (20) 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
(9)ω+υk2=0.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
(10)ω4+2(υk2+ηk2)ω3+[(υk2+ηk2)2+2(V2k2+ηk2υk2)+4υ02k4]ω2+[2(υk2+ηk2)×(V2k2+ηk2υk2)+8ηk2υ02k4]ω+(V2k2+ηk2υk2)2+4η2k4υ02k4=0.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
(11)ω3+[υk2+(γ-1)(TρLTp+λk2Tp)]ω2+[(γ-1)(TρLTp+λk2Tp)υk2+(c2k2-4πGρ)]ω+{k2(γ-1)(TLT-ρLρ+λk2Tρ)-4πGρ(γ-1)(TρLTp+λk2Tp)}=0.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
(12)k2(TLT-ρLρ+λk2Tρ)<4πGρ(TρLTp+λk2Tp).
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: (13)λJ1=[8π{(4Gρ/c′2+ρ2Lρ/πλT-ρLT/πλ)±[(4Gρ/c′2+ρ2Lρ/πλT-ρLT/πλ)2+(16Gρ2/πλc′2)LT]1/2}]1/2.The medium is unstable for wave lengthλ>λJ1. 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:
(14)k*=kc(4πGρ)1/2,ω*=ω(4πGρ)1/2,υ*=υ(4πGρ)1/2c2,LT*=(γ-1)ρTLTp(4πGρ)1/2,Lρ*=(γ-1)ρLρc2(4πGρ)1/2,λ*=(γ-1)Tλ(4πGρ)1/2pc2.
The nondimensional form of (11) in terms of self-gravitation is given as
(15)ω*3+(υ*k*2+LT*+λ*k*2)ω*2+[υ*k*2(LT*+λ*k*2)+k*2-1]ω*+{k*2[1γ(LT*+λ*k*2)-Lρ*]-(LT*+λ*k*2)}=0.Here for calculations we take the numerical values corresponding to the conditions in interstellar molecular clouds for star formation as(16)ρ=1.7×10-21kgm-3G=6.658×10-11(kg)-1m3s-2c2=2.25×108m2s-2V2=5.0×108m2s-2.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 k* for several values of the different parameters involved taking γ=5/3. 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 1016) has been plotted against the wave number k* (after multiplying by 1020) to show the dependence of the growth rate on the different physical parameters such as temperature-dependent heat-loss function LT* (after multiplying by108), density-dependent heat-loss function Lρ* (after multiplying by108), and viscosity υ* (after multiplying by 1024).Figure 1
Growth rate (positive values ofω*) against wave number k* for three values of parameter LT*=0.0,2.0,5.0, keeping the other parameters fixed Lρ*=1.0,  λ*=1.0,  υ*=1.0.Figure 2
Growth rate (positive values ofω*) against wave number k* for three values of parameter Lρ*=0.0,2.0,5.0, keeping the other parameters fixed λ*=1.0,  LT*=1.0,  υ*=1.0.Figure 3
Growth rate (positive values ofω*) against wave number k* for three values of parameter υ*=0.0,2.0,5.0, keeping the other parameters fixed Lρ*=1.0,  λ*=1.0,  LT*=1.0.It is clear from Figure1 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
(17)Δ1={υk2+(γ-1)(TρLTp+λk2Tp)}>0,Δ2={υk2[Δ1(γ-1)(TρLTp+λk2Tp)+c2k2-4πGρ]+(γ-1)k2ρLρ(TρLTp+λk2Tp)}>0,Δ3=Δ2{k2(γ-1)(TLT-ρLρ+λk2Tρ)-4πGρ(γ-1)(TρLTp+λk2Tp)}>0.
Since Ωj2>0,  ΩI2>0 and γ>1, it is clear that all the  Δs  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 FunctionFor nonradiative medium(LT,ρ=0,υ=λ≠0) (11) becomes
(18)ω3+(υk2+γλk2ρcp)ω2+(υk2γλk2ρcp+c2k2-4πGρ)ω+γλk2ρcp(c'2k2-4πGρ)=0.
The condition of instability from constant term of (18) is (c'2k2-4πGρ)<0, or
(19)k2<4πGρc'2
if
(20)kJ22=4πGρc'2,orλJ2=c′[πGρ]1/2.The medium is unstable for wave lengthλ>λJ2, where λJ2 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(λ=LT,ρ=0,υ≠0) the condition of instability obtained from (11) is
(21)(c2k2-4πGρ)<0,(22)λJ=c[πGρ]1/2.The medium is unstable for wave lengthλ>λJ, where λJ 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 FunctionFor inviscid and nonradiating medium(υ=LT,ρ=0,λ≠0) (11) becomes
(23)ω3+γλk2ρcpω2+(c2k2-4πGρ)ω+γλk2ρcp(c′2k2-4πGρ)=0.
The condition of instability is (c'2k2-4πGρ)<0, 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(υ=λ=LT,ρ=0) the condition of instability obtained from (11) is (c2k2-4πGρ)<0, 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(kx=k,kz=0)

For this case, we assume that all the perturbations are transverse to the direction of magnetic field (i.e.,kx=k,kz=0). The dispersion relation (7) reduces to(24)(ω+υk2)2{(γ-1){k2(TLT-ρLρ+λk2T/p)-4πGρ(TρLT/p+λk2T/p)}ω+(γ-1)(TρLT/p+λk2T/p)]{ωυ02k4+(ω+υk2)[ω2+ωυk2+ωV2k2ω+ηk2+ω(c2k2-4πGρ)ω+(γ-1)(TρLT/p+λk2T/p)}ωυ02k4+(ω+υk2)[(γ-1){k2(TLT-ρLρ+λk2T/p)-4πGρ(TρLT/p+λk2T/p)}ω+(γ-1)(TρLT/p+λk2T/p)ω2+ωυk2+ωV2k2ω+ηk2+ω(c2k2-4πGρ)ω+(γ-1)(TρLT/p+λk2T/p)+(γ-1){k2(TLT-ρLρ+λk2T/p)-4πGρ(TρLT/p+λk2T/p)}ω+(γ-1)(TρLT/p+λk2T/p)]}=0.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
(25)ω5+[2υk2+ηk2+(γ-1)(TρLTp+λk2Tp)]ω4+{2υk2[ηk2+(γ-1)(TρLTp+λk2Tp)]+υ2k4+ηk2(γ-1)(TρLTp+λk2Tp)+V2k2+υ02k4+c2k2-4πGρ[ηk2+(γ-1)(TρLTp+λk2Tp)]}ω3+{(TρLTp+λk2Tp)2ηk2υk2(γ-1)×(TρLTp+λk2Tp)+[ηk2+(γ-1)(TρLTp+λk2Tp)](υ2k4+υ02k4)+V2k2[υk2+(γ-1)×(TρLTp+λk2Tp)]+(c2k2-4πGρ)(υk2+ηk2)+(γ-1)[k2(TLT-ρLρ+λk2Tρ)-4πGρ×(TρLTp+λk2Tp)]}ω2+{υk2[ηk2υk2(γ-1)(TρLTp+λk2Tp)+V2k2(γ-1)×(TρLTp+λk2Tp)+ηk2(c2k2-4πGρ)]+ηk2υ02k4(γ-1)(TρLTp+λk2Tp)+(γ-1)[k2(TLT-ρLρ+λk2Tρ)-4πGρ(TρLTp+λk2Tp)](υk2+ηk2)}ω+ηk2υk2×[k2(γ-1)(TLT-ρLρ+λk2Tρ)-4πGρ(γ-1)(TρLTp+λk2Tp)]=0.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
(26)(γ-1)[k2(TLT-ρLρ+λk2Tρ)-4πGρ(TρLTp+λk2Tp)]<0.
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:
(27)V*2=V2c2,η*=η(4πGρ)1/2c2,υ0*=υ0(4πGρ)1/2c2.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 k* for several values of the different parameters involved taking γ=5/3. 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 1016) has been plotted against the wave number k* (after multiplying by 1020) to show the dependence of the growth rate on the different physical parameters such as temperature-dependent heat-loss function LT* (after multiplying by108), density-dependent heat-loss function Lρ* (after multiplying by 108), viscosity  υ* (after multiplying by 1024), and FLR corrections υ0* (after multiplying by 1024).Figure 4
Growth rate (positive values ofω*) against wave number k* for three values of parameter LT*=0.0,1.5,3.0, keeping the other parameters fixed υ*=0.5,  Lρ*=0.0,  λ*=1.0,  υ0*=1.0,  V*=1.0,  η*=0.5.Figure 5
Growth rate (positive values ofω*) against wave number k* for three values of parameter Lρ*=0.0,1.5,3.0, keeping the other parameters fixed υ0*=1.0,  λ*=1.0,  υ*=0.5,  V*=1.0,  η*=0.5,  LT*=0.5.Figure 6
Growth rate (positive values ofω*) against wave number k*  for several values of parameter υ0*  keeping the other parameters fixed Lρ*=0.5,  LT*=0.5,  λ*=1.0,  V*=1.0,  η*=0.5.  υ*=0.0: (1)  υ0*=0.0, (2) υ0*=1.5 (3) υ0*=3.0.  υ*=0.9: (4) υ0*=1.5, (5) υ0*=3.0.It is clear from Figure4 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 Ωj2>0,  ΩI2>0 and γ>1, so it is sufficient to show that alternate  Δs  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 CorrectionsFor nonradiative medium, but in the presence of FLR corrections(LT,ρ=0,υ=λ=V=η=υ0≠0) the condition of instability obtained from the constant term of (25) is given as
(28)(c'2k2-4πGρ)<0.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 (LT,ρ=λ=0,  υ=V=η=υ0≠0) the condition of instability obtained from (25) is given as
(29)(c2k2-4πGρ)<0.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 CorrectionsFor inviscid medium, but in presence of FLR corrections(υ=0,η=LT,ρ=V=λ=υ0≠0) (25) becomes
(30)ω4+[ηk2+(γ-1)(TρLTp+λk2Tp)]ω3+[ηk2(γ-1)(TρLTp+λk2Tp)+V2k2+υ02k4+c2k2-4πGρ(TρLTp+λk2Tp)]ω2+{υ02k4[ηk2+(γ-1)(TρLTp+λk2Tp)]+V2k2(γ-1)×(TρLTp+λk2Tp)+ηk2(c2k2-4πGρ)+[k2(γ-1)(TLT-ρLρ+λk2Tρ)-4πGρ(γ-1)(TρLTp+λk2Tp)]}ω+ηk2{υ02k4(γ-1)(TρLTp+λk2Tp)+k2(γ-1)×(TLT-ρLρ+λk2Tρ)-4πGρ(γ-1)(TρLTp+λk2Tp)}=0.
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
(31){k2(TLT-ρLρ+λk2Tρ)-(4πGρ-υ02k4)(TρLTp+λk2Tp)}<0.
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
(32){(λ*k*2+LT*)[υ0*2k*4+k*2γ-1]-k*2Lρ*}<0.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 LT* and Lρ* of the heat-loss function with respect of local temperature and density, and υ0* is FLR correction factor in the system. It may be noted that for a density-independent heat-loss function (Lρ*=0), which increases with temperature (LT*>0), condition (32) gives instability if k*2<(γ-1/γυ0*2). However if instead, the heat-loss function decreases with temperature (LT*<0), and the instability comes to k*2 lying between (γ-1/γυ0*2) and {|LT*|/λ*}.For a purely density-dependent heat-loss function(LT*=0), condition (32) gives instability if
(33)k*2<12{±[(1γυ0*2)2+4υ0*2(1+Lρ*λ*)]1/2-(1γυ0*2)}.
The above equation (33) applies for Lρ*>0 or |Lρ*|<λ* if the heat-loss function, respectively, increases or decreases with increase in density.(2)
Effect of Radiative Heat-Loss Function with FLR CorrectionsFor nonradiating and inviscid medium, but in presence of FLR(LT,ρ=υ=0,η=λ=υ0=V≠0) the condition of instability obtained from (25) is given as
(34)k2(c'2+υ02k2)<4πGρ.
The critical Jeans wave number and critical Jeans wave length are given as
(35)kJ32={-c'2±(c'4+16πGρυ02)1/2}2υ02.
Or
(36)λJ3=[8π2υ02{-c'2±(c'4+16πGρυ02)1/2}]1/2.Condition (34) is identical of Vaghela and Chhajlani [19] for (ε=1) 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(υ=LT,ρ=λ=0,η=υ0≠0), the condition of instability obtained from (25) is given as
(37){(c2+υ02k2)k2-4πGρ}<0.
The critical Jeans wave number or critical Jeans wave length is given as
(38)kJ42={-c2±(c4+16πGρυ02)1/2}2υ02.
Or
(39)λJ4=[8π2υ02{-c2±(c4+16πGρυ02)1/2}]1/2.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 CorrectionsFor inviscid and infinitely conducting medium, but in presence of FLR (υ=η=0,LT,ρ=V=λ=υ0≠0), the condition of instability obtained from constant term of  (25) is given as
(40){k2(TLT-ρLρ+λk2Tρ)-(4πGρ-υ02k4-V2k2)(TρLTp+λk2Tp)}<0.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
(41){(λ*k*2+LT*)[υ0*2k*4+k*2(1γ+V*2)-1]-k*2Lρ*(1γ+V*2)}<0.
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 (Lρ*=0), which increases with temperature (LT*>0), it gives instability if k*2<(1/υ0*2){1-[(1/γ)+V*2]}. However if instead, the heat-loss function decreases with temperature (LT*<0), the instability comes to k*2 lying between (1/υ0*2){1-[(1/γ)+V*2]} and {|LT*|/λ*}.For a purely density-dependent heat-loss function(LT*=0), condition (41) gives instability if
(42)k*2<12{±[1υ0*4(1γ+V*2)2+4υ0*2(1+Lρ*λ*)]1/2-1υ0*2(1γ+V*2)[1υ0*4(1γ+V*2)2+4υ0*2(1+Lρ*λ*)]1/2}.
The above equation (42) applies for Lρ*>0 or |Lρ*|<λ*, if the heat-loss function, respectively, increases or decreases with increase in density.(4)
Effect of Nongravitation and FLR CorrectionsFor inviscid, nongravitating, infinitely conducting, radiating, thermally conducting, magnetized medium with FLR corrections(υ=η=G=0,υ0=V=LT,ρ=λ≠0) (25) becomes
(43)ω3+[(γ-1)(TρLTp+λk2Tp)]ω2+[V2k2+υ02k4+c2k2]ω+(γ-1)k2[(TρLTp+λk2Tp)(V2+υ02k2)×(TρLTp+λk2Tp)+(TLT-ρLρ+λk2Tρ)]=0.
The condition of instability obtained from constant term of (43) is given as
(44)[(υ02k2+V2)(TρLTp+λk2Tp)+(TLT+λk2Tρ-ρLρ)]<0.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.


## 4.1. Longitudinal Propagation(kx=0,kz=k)

For this case, we assume all the perturbations are longitudinal to the direction of magnetic field (i.e.,kz=k,  kx=0). The dispersion relation (7) reduces to(8)(ω+υk2)×{4υ02k4+[ω+υk2+V2k2ω+ηk2]2}×[(γ-1){k2(TLT-ρLρ+λk2T/p)-4πGρ(TρLT/p+λk2T/p)}ω+(γ-1)(TρLT/p+λk2T/p)ω2+ωυk2+ω(c2k2-4πGρ)ω+(γ-1)(TρLT/p+λk2T/p)+(γ-1){k2(TLT-ρLρ+λk2T/p)-4πGρ(TρLT/p+λk2T/p)}ω+(γ-1)(TρLT/p+λk2T/p)]=0.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 (21) 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 (20) 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
(9)ω+υk2=0.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
(10)ω4+2(υk2+ηk2)ω3+[(υk2+ηk2)2+2(V2k2+ηk2υk2)+4υ02k4]ω2+[2(υk2+ηk2)×(V2k2+ηk2υk2)+8ηk2υ02k4]ω+(V2k2+ηk2υk2)2+4η2k4υ02k4=0.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
(11)ω3+[υk2+(γ-1)(TρLTp+λk2Tp)]ω2+[(γ-1)(TρLTp+λk2Tp)υk2+(c2k2-4πGρ)]ω+{k2(γ-1)(TLT-ρLρ+λk2Tρ)-4πGρ(γ-1)(TρLTp+λk2Tp)}=0.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
(12)k2(TLT-ρLρ+λk2Tρ)<4πGρ(TρLTp+λk2Tp).
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: (13)λJ1=[8π{(4Gρ/c′2+ρ2Lρ/πλT-ρLT/πλ)±[(4Gρ/c′2+ρ2Lρ/πλT-ρLT/πλ)2+(16Gρ2/πλc′2)LT]1/2}]1/2.The medium is unstable for wave lengthλ>λJ1. 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:
(14)k*=kc(4πGρ)1/2,ω*=ω(4πGρ)1/2,υ*=υ(4πGρ)1/2c2,LT*=(γ-1)ρTLTp(4πGρ)1/2,Lρ*=(γ-1)ρLρc2(4πGρ)1/2,λ*=(γ-1)Tλ(4πGρ)1/2pc2.
The nondimensional form of (11) in terms of self-gravitation is given as
(15)ω*3+(υ*k*2+LT*+λ*k*2)ω*2+[υ*k*2(LT*+λ*k*2)+k*2-1]ω*+{k*2[1γ(LT*+λ*k*2)-Lρ*]-(LT*+λ*k*2)}=0.Here for calculations we take the numerical values corresponding to the conditions in interstellar molecular clouds for star formation as(16)ρ=1.7×10-21kgm-3G=6.658×10-11(kg)-1m3s-2c2=2.25×108m2s-2V2=5.0×108m2s-2.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 k* for several values of the different parameters involved taking γ=5/3. 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 1016) has been plotted against the wave number k* (after multiplying by 1020) to show the dependence of the growth rate on the different physical parameters such as temperature-dependent heat-loss function LT* (after multiplying by108), density-dependent heat-loss function Lρ* (after multiplying by108), and viscosity υ* (after multiplying by 1024).Figure 1
Growth rate (positive values ofω*) against wave number k* for three values of parameter LT*=0.0,2.0,5.0, keeping the other parameters fixed Lρ*=1.0,  λ*=1.0,  υ*=1.0.Figure 2
Growth rate (positive values ofω*) against wave number k* for three values of parameter Lρ*=0.0,2.0,5.0, keeping the other parameters fixed λ*=1.0,  LT*=1.0,  υ*=1.0.Figure 3
Growth rate (positive values ofω*) against wave number k* for three values of parameter υ*=0.0,2.0,5.0, keeping the other parameters fixed Lρ*=1.0,  λ*=1.0,  LT*=1.0.It is clear from Figure1 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
(17)Δ1={υk2+(γ-1)(TρLTp+λk2Tp)}>0,Δ2={υk2[Δ1(γ-1)(TρLTp+λk2Tp)+c2k2-4πGρ]+(γ-1)k2ρLρ(TρLTp+λk2Tp)}>0,Δ3=Δ2{k2(γ-1)(TLT-ρLρ+λk2Tρ)-4πGρ(γ-1)(TρLTp+λk2Tp)}>0.
Since Ωj2>0,  ΩI2>0 and γ>1, it is clear that all the  Δs  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 FunctionFor nonradiative medium(LT,ρ=0,υ=λ≠0) (11) becomes
(18)ω3+(υk2+γλk2ρcp)ω2+(υk2γλk2ρcp+c2k2-4πGρ)ω+γλk2ρcp(c'2k2-4πGρ)=0.
The condition of instability from constant term of (18) is (c'2k2-4πGρ)<0, or
(19)k2<4πGρc'2
if
(20)kJ22=4πGρc'2,orλJ2=c′[πGρ]1/2.The medium is unstable for wave lengthλ>λJ2, where λJ2 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(λ=LT,ρ=0,υ≠0) the condition of instability obtained from (11) is
(21)(c2k2-4πGρ)<0,(22)λJ=c[πGρ]1/2.The medium is unstable for wave lengthλ>λJ, where λJ 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 FunctionFor inviscid and nonradiating medium(υ=LT,ρ=0,λ≠0) (11) becomes
(23)ω3+γλk2ρcpω2+(c2k2-4πGρ)ω+γλk2ρcp(c′2k2-4πGρ)=0.
The condition of instability is (c'2k2-4πGρ)<0, 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(υ=λ=LT,ρ=0) the condition of instability obtained from (11) is (c2k2-4πGρ)<0, 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.1.1. For Viscous Medium

(1)
Effect of Radiative Heat-Loss FunctionFor nonradiative medium(LT,ρ=0,υ=λ≠0) (11) becomes
(18)ω3+(υk2+γλk2ρcp)ω2+(υk2γλk2ρcp+c2k2-4πGρ)ω+γλk2ρcp(c'2k2-4πGρ)=0.
The condition of instability from constant term of (18) is (c'2k2-4πGρ)<0, or
(19)k2<4πGρc'2
if
(20)kJ22=4πGρc'2,orλJ2=c′[πGρ]1/2.The medium is unstable for wave lengthλ>λJ2, where λJ2 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(λ=LT,ρ=0,υ≠0) the condition of instability obtained from (11) is
(21)(c2k2-4πGρ)<0,(22)λJ=c[πGρ]1/2.The medium is unstable for wave lengthλ>λJ, where λJ 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 FunctionFor inviscid and nonradiating medium(υ=LT,ρ=0,λ≠0) (11) becomes
(23)ω3+γλk2ρcpω2+(c2k2-4πGρ)ω+γλk2ρcp(c′2k2-4πGρ)=0.
The condition of instability is (c'2k2-4πGρ)<0, 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(υ=λ=LT,ρ=0) the condition of instability obtained from (11) is (c2k2-4πGρ)<0, 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(kx=k,kz=0)

For this case, we assume that all the perturbations are transverse to the direction of magnetic field (i.e.,kx=k,kz=0). The dispersion relation (7) reduces to(24)(ω+υk2)2{(γ-1){k2(TLT-ρLρ+λk2T/p)-4πGρ(TρLT/p+λk2T/p)}ω+(γ-1)(TρLT/p+λk2T/p)]{ωυ02k4+(ω+υk2)[ω2+ωυk2+ωV2k2ω+ηk2+ω(c2k2-4πGρ)ω+(γ-1)(TρLT/p+λk2T/p)}ωυ02k4+(ω+υk2)[(γ-1){k2(TLT-ρLρ+λk2T/p)-4πGρ(TρLT/p+λk2T/p)}ω+(γ-1)(TρLT/p+λk2T/p)ω2+ωυk2+ωV2k2ω+ηk2+ω(c2k2-4πGρ)ω+(γ-1)(TρLT/p+λk2T/p)+(γ-1){k2(TLT-ρLρ+λk2T/p)-4πGρ(TρLT/p+λk2T/p)}ω+(γ-1)(TρLT/p+λk2T/p)]}=0.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
(25)ω5+[2υk2+ηk2+(γ-1)(TρLTp+λk2Tp)]ω4+{2υk2[ηk2+(γ-1)(TρLTp+λk2Tp)]+υ2k4+ηk2(γ-1)(TρLTp+λk2Tp)+V2k2+υ02k4+c2k2-4πGρ[ηk2+(γ-1)(TρLTp+λk2Tp)]}ω3+{(TρLTp+λk2Tp)2ηk2υk2(γ-1)×(TρLTp+λk2Tp)+[ηk2+(γ-1)(TρLTp+λk2Tp)](υ2k4+υ02k4)+V2k2[υk2+(γ-1)×(TρLTp+λk2Tp)]+(c2k2-4πGρ)(υk2+ηk2)+(γ-1)[k2(TLT-ρLρ+λk2Tρ)-4πGρ×(TρLTp+λk2Tp)]}ω2+{υk2[ηk2υk2(γ-1)(TρLTp+λk2Tp)+V2k2(γ-1)×(TρLTp+λk2Tp)+ηk2(c2k2-4πGρ)]+ηk2υ02k4(γ-1)(TρLTp+λk2Tp)+(γ-1)[k2(TLT-ρLρ+λk2Tρ)-4πGρ(TρLTp+λk2Tp)](υk2+ηk2)}ω+ηk2υk2×[k2(γ-1)(TLT-ρLρ+λk2Tρ)-4πGρ(γ-1)(TρLTp+λk2Tp)]=0.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
(26)(γ-1)[k2(TLT-ρLρ+λk2Tρ)-4πGρ(TρLTp+λk2Tp)]<0.
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:
(27)V*2=V2c2,η*=η(4πGρ)1/2c2,υ0*=υ0(4πGρ)1/2c2.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 k* for several values of the different parameters involved taking γ=5/3. 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 1016) has been plotted against the wave number k* (after multiplying by 1020) to show the dependence of the growth rate on the different physical parameters such as temperature-dependent heat-loss function LT* (after multiplying by108), density-dependent heat-loss function Lρ* (after multiplying by 108), viscosity  υ* (after multiplying by 1024), and FLR corrections υ0* (after multiplying by 1024).Figure 4
Growth rate (positive values ofω*) against wave number k* for three values of parameter LT*=0.0,1.5,3.0, keeping the other parameters fixed υ*=0.5,  Lρ*=0.0,  λ*=1.0,  υ0*=1.0,  V*=1.0,  η*=0.5.Figure 5
Growth rate (positive values ofω*) against wave number k* for three values of parameter Lρ*=0.0,1.5,3.0, keeping the other parameters fixed υ0*=1.0,  λ*=1.0,  υ*=0.5,  V*=1.0,  η*=0.5,  LT*=0.5.Figure 6
Growth rate (positive values ofω*) against wave number k*  for several values of parameter υ0*  keeping the other parameters fixed Lρ*=0.5,  LT*=0.5,  λ*=1.0,  V*=1.0,  η*=0.5.  υ*=0.0: (1)  υ0*=0.0, (2) υ0*=1.5 (3) υ0*=3.0.  υ*=0.9: (4) υ0*=1.5, (5) υ0*=3.0.It is clear from Figure4 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 Ωj2>0,  ΩI2>0 and γ>1, so it is sufficient to show that alternate  Δs  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 CorrectionsFor nonradiative medium, but in the presence of FLR corrections(LT,ρ=0,υ=λ=V=η=υ0≠0) the condition of instability obtained from the constant term of (25) is given as
(28)(c'2k2-4πGρ)<0.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 (LT,ρ=λ=0,  υ=V=η=υ0≠0) the condition of instability obtained from (25) is given as
(29)(c2k2-4πGρ)<0.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 CorrectionsFor inviscid medium, but in presence of FLR corrections(υ=0,η=LT,ρ=V=λ=υ0≠0) (25) becomes
(30)ω4+[ηk2+(γ-1)(TρLTp+λk2Tp)]ω3+[ηk2(γ-1)(TρLTp+λk2Tp)+V2k2+υ02k4+c2k2-4πGρ(TρLTp+λk2Tp)]ω2+{υ02k4[ηk2+(γ-1)(TρLTp+λk2Tp)]+V2k2(γ-1)×(TρLTp+λk2Tp)+ηk2(c2k2-4πGρ)+[k2(γ-1)(TLT-ρLρ+λk2Tρ)-4πGρ(γ-1)(TρLTp+λk2Tp)]}ω+ηk2{υ02k4(γ-1)(TρLTp+λk2Tp)+k2(γ-1)×(TLT-ρLρ+λk2Tρ)-4πGρ(γ-1)(TρLTp+λk2Tp)}=0.
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
(31){k2(TLT-ρLρ+λk2Tρ)-(4πGρ-υ02k4)(TρLTp+λk2Tp)}<0.
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
(32){(λ*k*2+LT*)[υ0*2k*4+k*2γ-1]-k*2Lρ*}<0.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 LT* and Lρ* of the heat-loss function with respect of local temperature and density, and υ0* is FLR correction factor in the system. It may be noted that for a density-independent heat-loss function (Lρ*=0), which increases with temperature (LT*>0), condition (32) gives instability if k*2<(γ-1/γυ0*2). However if instead, the heat-loss function decreases with temperature (LT*<0), and the instability comes to k*2 lying between (γ-1/γυ0*2) and {|LT*|/λ*}.For a purely density-dependent heat-loss function(LT*=0), condition (32) gives instability if
(33)k*2<12{±[(1γυ0*2)2+4υ0*2(1+Lρ*λ*)]1/2-(1γυ0*2)}.
The above equation (33) applies for Lρ*>0 or |Lρ*|<λ* if the heat-loss function, respectively, increases or decreases with increase in density.(2)
Effect of Radiative Heat-Loss Function with FLR CorrectionsFor nonradiating and inviscid medium, but in presence of FLR(LT,ρ=υ=0,η=λ=υ0=V≠0) the condition of instability obtained from (25) is given as
(34)k2(c'2+υ02k2)<4πGρ.
The critical Jeans wave number and critical Jeans wave length are given as
(35)kJ32={-c'2±(c'4+16πGρυ02)1/2}2υ02.
Or
(36)λJ3=[8π2υ02{-c'2±(c'4+16πGρυ02)1/2}]1/2.Condition (34) is identical of Vaghela and Chhajlani [19] for (ε=1) 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(υ=LT,ρ=λ=0,η=υ0≠0), the condition of instability obtained from (25) is given as
(37){(c2+υ02k2)k2-4πGρ}<0.
The critical Jeans wave number or critical Jeans wave length is given as
(38)kJ42={-c2±(c4+16πGρυ02)1/2}2υ02.
Or
(39)λJ4=[8π2υ02{-c2±(c4+16πGρυ02)1/2}]1/2.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 CorrectionsFor inviscid and infinitely conducting medium, but in presence of FLR (υ=η=0,LT,ρ=V=λ=υ0≠0), the condition of instability obtained from constant term of  (25) is given as
(40){k2(TLT-ρLρ+λk2Tρ)-(4πGρ-υ02k4-V2k2)(TρLTp+λk2Tp)}<0.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
(41){(λ*k*2+LT*)[υ0*2k*4+k*2(1γ+V*2)-1]-k*2Lρ*(1γ+V*2)}<0.
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 (Lρ*=0), which increases with temperature (LT*>0), it gives instability if k*2<(1/υ0*2){1-[(1/γ)+V*2]}. However if instead, the heat-loss function decreases with temperature (LT*<0), the instability comes to k*2 lying between (1/υ0*2){1-[(1/γ)+V*2]} and {|LT*|/λ*}.For a purely density-dependent heat-loss function(LT*=0), condition (41) gives instability if
(42)k*2<12{±[1υ0*4(1γ+V*2)2+4υ0*2(1+Lρ*λ*)]1/2-1υ0*2(1γ+V*2)[1υ0*4(1γ+V*2)2+4υ0*2(1+Lρ*λ*)]1/2}.
The above equation (42) applies for Lρ*>0 or |Lρ*|<λ*, if the heat-loss function, respectively, increases or decreases with increase in density.(4)
Effect of Nongravitation and FLR CorrectionsFor inviscid, nongravitating, infinitely conducting, radiating, thermally conducting, magnetized medium with FLR corrections(υ=η=G=0,υ0=V=LT,ρ=λ≠0) (25) becomes
(43)ω3+[(γ-1)(TρLTp+λk2Tp)]ω2+[V2k2+υ02k4+c2k2]ω+(γ-1)k2[(TρLTp+λk2Tp)(V2+υ02k2)×(TρLTp+λk2Tp)+(TLT-ρLρ+λk2Tρ)]=0.
The condition of instability obtained from constant term of (43) is given as
(44)[(υ02k2+V2)(TρLTp+λk2Tp)+(TLT+λk2Tρ-ρLρ)]<0.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.


## 4.2.1. For Viscous Medium

(1)
Effect of Radiative Heat-Loss Function with FLR CorrectionsFor nonradiative medium, but in the presence of FLR corrections(LT,ρ=0,υ=λ=V=η=υ0≠0) the condition of instability obtained from the constant term of (25) is given as
(28)(c'2k2-4πGρ)<0.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 (LT,ρ=λ=0,  υ=V=η=υ0≠0) the condition of instability obtained from (25) is given as
(29)(c2k2-4πGρ)<0.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 CorrectionsFor inviscid medium, but in presence of FLR corrections(υ=0,η=LT,ρ=V=λ=υ0≠0) (25) becomes
(30)ω4+[ηk2+(γ-1)(TρLTp+λk2Tp)]ω3+[ηk2(γ-1)(TρLTp+λk2Tp)+V2k2+υ02k4+c2k2-4πGρ(TρLTp+λk2Tp)]ω2+{υ02k4[ηk2+(γ-1)(TρLTp+λk2Tp)]+V2k2(γ-1)×(TρLTp+λk2Tp)+ηk2(c2k2-4πGρ)+[k2(γ-1)(TLT-ρLρ+λk2Tρ)-4πGρ(γ-1)(TρLTp+λk2Tp)]}ω+ηk2{υ02k4(γ-1)(TρLTp+λk2Tp)+k2(γ-1)×(TLT-ρLρ+λk2Tρ)-4πGρ(γ-1)(TρLTp+λk2Tp)}=0.
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
(31){k2(TLT-ρLρ+λk2Tρ)-(4πGρ-υ02k4)(TρLTp+λk2Tp)}<0.
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
(32){(λ*k*2+LT*)[υ0*2k*4+k*2γ-1]-k*2Lρ*}<0.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 LT* and Lρ* of the heat-loss function with respect of local temperature and density, and υ0* is FLR correction factor in the system. It may be noted that for a density-independent heat-loss function (Lρ*=0), which increases with temperature (LT*>0), condition (32) gives instability if k*2<(γ-1/γυ0*2). However if instead, the heat-loss function decreases with temperature (LT*<0), and the instability comes to k*2 lying between (γ-1/γυ0*2) and {|LT*|/λ*}.For a purely density-dependent heat-loss function(LT*=0), condition (32) gives instability if
(33)k*2<12{±[(1γυ0*2)2+4υ0*2(1+Lρ*λ*)]1/2-(1γυ0*2)}.
The above equation (33) applies for Lρ*>0 or |Lρ*|<λ* if the heat-loss function, respectively, increases or decreases with increase in density.(2)
Effect of Radiative Heat-Loss Function with FLR CorrectionsFor nonradiating and inviscid medium, but in presence of FLR(LT,ρ=υ=0,η=λ=υ0=V≠0) the condition of instability obtained from (25) is given as
(34)k2(c'2+υ02k2)<4πGρ.
The critical Jeans wave number and critical Jeans wave length are given as
(35)kJ32={-c'2±(c'4+16πGρυ02)1/2}2υ02.
Or
(36)λJ3=[8π2υ02{-c'2±(c'4+16πGρυ02)1/2}]1/2.Condition (34) is identical of Vaghela and Chhajlani [19] for (ε=1) 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(υ=LT,ρ=λ=0,η=υ0≠0), the condition of instability obtained from (25) is given as
(37){(c2+υ02k2)k2-4πGρ}<0.
The critical Jeans wave number or critical Jeans wave length is given as
(38)kJ42={-c2±(c4+16πGρυ02)1/2}2υ02.
Or
(39)λJ4=[8π2υ02{-c2±(c4+16πGρυ02)1/2}]1/2.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 CorrectionsFor inviscid and infinitely conducting medium, but in presence of FLR (υ=η=0,LT,ρ=V=λ=υ0≠0), the condition of instability obtained from constant term of  (25) is given as
(40){k2(TLT-ρLρ+λk2Tρ)-(4πGρ-υ02k4-V2k2)(TρLTp+λk2Tp)}<0.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
(41){(λ*k*2+LT*)[υ0*2k*4+k*2(1γ+V*2)-1]-k*2Lρ*(1γ+V*2)}<0.
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 (Lρ*=0), which increases with temperature (LT*>0), it gives instability if k*2<(1/υ0*2){1-[(1/γ)+V*2]}. However if instead, the heat-loss function decreases with temperature (LT*<0), the instability comes to k*2 lying between (1/υ0*2){1-[(1/γ)+V*2]} and {|LT*|/λ*}.For a purely density-dependent heat-loss function(LT*=0), condition (41) gives instability if
(42)k*2<12{±[1υ0*4(1γ+V*2)2+4υ0*2(1+Lρ*λ*)]1/2-1υ0*2(1γ+V*2)[1υ0*4(1γ+V*2)2+4υ0*2(1+Lρ*λ*)]1/2}.
The above equation (42) applies for Lρ*>0 or |Lρ*|<λ*, if the heat-loss function, respectively, increases or decreases with increase in density.(4)
Effect of Nongravitation and FLR CorrectionsFor inviscid, nongravitating, infinitely conducting, radiating, thermally conducting, magnetized medium with FLR corrections(υ=η=G=0,υ0=V=LT,ρ=λ≠0) (25) becomes
(43)ω3+[(γ-1)(TρLTp+λk2Tp)]ω2+[V2k2+υ02k4+c2k2]ω+(γ-1)k2[(TρLTp+λk2Tp)(V2+υ02k2)×(TρLTp+λk2Tp)+(TLT-ρLρ+λk2Tρ)]=0.
The condition of instability obtained from constant term of (43) is given as
(44)[(υ02k2+V2)(TρLTp+λk2Tp)+(TLT+λk2Tρ-ρLρ)]<0.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.



---

*Source: 420938-2012-09-29.xml*