Magnetogravitational instability in strongly coupled rotating clumpy molecular clouds including heating and cooling functions

Molecular cloud (MC) formation is caused by the gravitational collapse mechanism and is significantly affected by radiative heating and cooling processes. This paper analyzes the gravitational instability in strongly coupled clumpy molecular clouds (MCs), under the effects of uniform rotation, magnetic field, and heat-loss functions. The generalized hydrodynamic equations coupled with the modified energy equation (which incorporates the heating and cooling effects due to cloud-cloud collisions) are used to describe the mathematical model. Following Jeans stability analysis, it is found that the value of the critical Jeans wavenumber decreases due to the strong coupling between the plasma particles (coupling parameter) and clump stirring processes (heating rate), so both have a stabilizing influence on the onset of gravitational collapse in clumpy MCs. The influence of various parameters on the growth rate of the instability is discussed numerically, and it is found that the cooling rate parameter that describes cloud-cloud collisions shows a destabilizing effect. The region of instability is observed to be smaller in the strongly coupled clumps (kinetic limit) than in the weakly coupled (hydrodynamic limit) clumps. The results are helpful in understanding the role of heating and cooling mechanisms in the MC formation.


Introduction
The radiative transfer mechanism in astrophysical systems plays a decisive role in the pre-stage star formation process (Juvela 2011).In the interstellar medium (ISM), the thermal processes allow the gas to become denser and hotter, providing the necessary conditions for gravitational collapse.The ISM observations report that the radiative heat-loss mechanism via thermal instability plays an important role in the star formation process in molecular clouds (MCs) (Fukue and Kamaya 2007).Wurster et al. (2018) studied the radiation and non-ideal magnetohydrodynamic (MHD) effects, e.g., ambipolar diffusion, Ohmic resistivity, and the Hall effect, on the collapse of rotating and magnetized molecular cloud cores.Hegmann and Kegel (1996) presented a radia-tive transfer model in the form of a Fokker-Planck equation, considering the density fluctuations in the system.Patrick and Shu-ichiro (2019) explored the influence of magnetic fields on the formation of pre-stellar dense cores and stellar clusters.Molaro et al. (2016) suggested a new method for probing global properties of the clump and core population in giant molecular clouds (GMCs), based on the study of their overall effect on the reflected X-ray signal.Bethell et al. (2007) calculated the photoionization rates in clumpy MCs.Micic et al. (2013) studied the influence of cooling function on the formation of MCs, using high-resolution three-dimensional simulations of converging flows.
MCs are the main sources of star formation in which the equilibrium temperature is maintained due to the balance of heating and cooling mechanisms (Hollenbach 1988).In the ongoing radiative processes in the star formation mechanism, Jeans instability plays a fundamental role in the gravitational collapse of a gas cloud.Jeans (1902) described the conditions under which a region of gas or plasma can collapse gravitationally and form a new structure, such as a star or a planet.The minimum perturbation wavenumber for which Jeans instability is excited in a system is expressed by the Jeans instability criterion as k < k j = (4πGρ 0 ) 1/2 /c s , where G is the gravitational constant, ρ 0 is the gas density, and c s is the sound speed.If the wavenumber of perturbations in the system is smaller than the critical Jeans wavenumber, then the system will be unstable and collapse under its own gravity, leading to the formation of structures such as stars and galaxies.In other words, if the size of the perturbation is larger than some critical wavelength of the order λ j = 2π/k j , then the enhanced self-gravity can overpower the excess pressure so that the perturbation grows.Motivated by the works of Jeans (1902), many authors have shown keen interest in exploring the mechanism of star formation and studied Jeans instability in various astrophysical systems, such as gaseous clouds (Chandrasekhar 1961), clumpy MCs (Elmegreen 1989), strongly coupled medium (Janaki et al. 2011;Prajapati and Chhajlani 2013;Dhiman and Sharma 2014;Dasgupta and Karmakar 2019;Dolai and Prajapati 2020), radiative MCs (Prajapati 2022), and self-gravitating optically thick gas-and-dust medium (Kolesnichenko 2023).
A strongly coupled medium refers to a state of matter where the interactions between charged particles, such as ions and electrons, are so strong that they begin to interact with each other through long-range electric and magnetic forces, which can result in the formation of complex structures and patterns.The present work aims to examine the role of heating and cooling effects on the growth rate of the Jeans instability in strongly coupled fluid (for which j 1, i.e., when average potential energy dominates the kinetic energy of the particles).Previous works include the effects of viscoelastic coefficients (Janaki et al. 2011), electrical resistivity (Prajapati and Chhajlani 2013), rotation (Dhiman and Sharma 2014) and quantum corrections (Sharma and Chhajlani 2014) on the Jeans instability in strongly coupled viscoelastic fluids.Sharma et al. (2015) investigated the Jeans instability of rotating viscoelastic fluids including magnetic fields.Dhiman and Mahajan (2023) studied the gravitational instability in strongly coupled viscoelastic clumpy MCs considering dissipative effects.However, the effects of the magnetic field and rotation were ignored in the considered configuration.Recently, Yang et al. (2023) studied Jeans instability in viscoelastic astrofluids using Eddington-Inspired-Born-Infeld (EiBI) gravity and found that the positive EiBI gravity parameter and effective generalized fluid viscosity act as stabilizing agents whereas the negative EiBI parameter and viscoelastic relaxation time act as destabilizing factors.From the above studies, it is evident that the strong coupling effects in terms of viscoelastic coefficients play a dominant role in investigating the gravitational instability in viscoelastic fluids.
The presence of magnetic fields in MCs is essential for understanding their role in the evolution of dense clouds and star formation process (Braine et al. 2020;Crutcher 2012).Körtgen and Banerjee (2015) showed that dense cores of MCs could build up under all the conditions; however, the star formation process in these cores is either delayed or completely suppressed if the initial field strength is B > 3μG.Ibáñez-Mejía et al. (2022) discussed the competition between gravity versus the magnetic fields in the starforming region of the MCs.Recently, Patidar et al. (2023) investigated the influence of the magnetic field on the Jeans criteria for quantum plasma and found that highly magnetized white dwarfs have a different mass-radius relation than their non-magnetic counterparts, which results in a modified super-Chandrashekhar mass limit.In addition, rotation is a unique characteristic of astrophysical systems such as accretion disks, circumstellar disks, and MCs.The spin angular momenta of the clouds may be due to the orbital rotation of the gas.The observational data show that GMCs have a spin rotation with angular momenta of the order of 100 km s −1 pc per unit mass (Chernin and Efremov 1995).In the gravitational collapse and subsequent star formation processes, the role of the magnetic field and rotation has been studied previously by many authors.However to the best of our knowledge, none of them has studied the combined effects of heating and cooling functions, magnetic field and rotation on the gravitational instability of strongly coupled clumpy MCs.
In this paper, we have analyzed the gravitational instability of strongly coupled MCs and have investigated the effects of rotation, magnetic field, and heat-loss functions on the onset of gravitational collapse.Here, we have considered the modified energy equation for an ideal fluid with heating and cooling functions (Elmegreen 1989) along with standard fluid equations.The paper is organized as follows: Sec. 1 is the introductory in nature and also includes a review of the literature.In Sec. 2, the mathematical model is constructed for the rotating and magnetized clumpy MCs, using the generalized hydrodynamic equations, which include the effect of the Coulomb coupling parameter j .In Sec. 3, the linearized perturbed equations are obtained, the general dispersion relation is derived using the plane wave solutions, and the results are discussed in parallel wave propagation mode for both the strongly (kinetic) and weakly (hydrodynamic) coupled limits.In Sec. 4, the outcomes of the present work have been summarized.

Mathematical model
Consider an infinite homogeneous, self-gravitating clumpy MCs in which the uniform rotation ( (0, 0, z )) and uniform magnetic field (H(0, 0, H 0 )) act simultaneously.The considered medium is strongly coupled in nature and hence exhibits both viscous and elastic behaviour in terms of the viscoelastic coefficients.Following Kaw and Sen (1998), the basic equations for the considered system using the generalized hydrodynamic (GH) model are given by (Elmegreen 1989;Janaki et al. 2011;Dhiman and Sharma 2014): The equation of continuity; The equation of motion; (2) The magnetic induction equation; The Poisson equation for the gravitational field; Here, g = −∇φ, φ is the gravitational potential, u, P and ρ represent the fluid velocity, fluid pressure, and fluid density, respectively.τ m is the relaxation time or memory parameter, which depends on the viscoelastic coefficients There has been a great deal of concern among researchers about analyzing the role of the radiative heat-loss mechanism in MCs due to its significance in structure formation.Field (1965) studied the effects of temperature and density-dependent radiative heat-loss functions on the thermal stability of a dilute, non-gravitating gas.Later, this problem was extensively studied by many authors, considering the effects of finite Larmor frequency and electrical resistivity (Bora and Talwar 1993), uniform rotation and magnetic field (Aggarwal and Talwar 1969), ion-neutral collisions (Fukue and Kamaya 2007), Hall effect and electron inertia (Prajapati et al. 2010), dust temperature and viscoelastic effects (Sharma et al. 2020) and porosity of the medium (Kaothekar et al. 2022).In these works, the temperature and density-dependent radiative heat-loss functions are considered and it is found that the fundamental thermal instability criterion remains unchanged due to the considered effects.However, when we consider cloud-cloud collisions, the heating and cooling functions are considered to be independent of time (Elmegreen 1989) and these functions are different than those considered by Field (1965) in his work.In the present work, we shall focus on the physical conditions of molecular cloud clumps to analyze the Jeans instability of strongly coupled fluids in the presence of rotation and magnetic fields, considering the heating and cooling effects.
Thus, the energy equation is therefore considered as (Elmegreen 1989); where, γ is the ratio of specific heat.The term represents the cooling function, which is attributed to cloud-cloud collisions with an isotropic, Maxwellian distribution of cloud velocities and is defined as; Here, M c , R c and σ c respectively denote the average masses, radii, and column densities of the colliding clouds.Further, = e ρ( 0 ρ 0 ) describes the heating rate where e , 0 are the values of heating and cooling rate in the equilibrium states.

Dispersion relation
The dispersion relation for the considered configuration is obtained using the first-order perturbation method.The basic equations are linearized, considering infinitesimal small perturbations; where, ρ 1 , g 1 , P 1 , u 1 (v x , v y , v z ), and H 1 (h x , h y , h z ) denote the perturbations in fluid density, gravitational field, fluid pressure, fluid velocity and magnetic field, respectively.All the above quantities with subscript "0" represent the value of that quantity at equilibrium.In the stability analysis of gravitating fluid, the streaming effects are neglected; hence we put u 0 = 0. To check the stability of the perturbed system, the arbitrary perturbations can be decomposed into a complete set of normal modes and are analyzed individually.Therefore, assuming the plane wave solutions of each of the perturbed quantity as; where, ω = −iω is the perturbation frequency of the harmonic disturbance and k z is the wavenumber.This allows us to substitute ∇ = ik z ẑ and ∂/∂t = ω in the perturbed equations obtained after substituting the perturbed quantities (7) into equations ( 1)-( 5), we obtained the following linearized equations; (1 where ω c = (γ − 1) 0 /(2ρ 0 c 2 ) represents the cooling rate (Elmegreen 1989) and is approximated as the inverse of the cloud-cloud collision time.Equation ( 13) can be simplified as Using the values of ρ 1 , h x , h y , g 1 and P 1 from equations ( 10)-( 12) and ( 14) in equation ( 9), we obtain the following dispersion relation after some simplifications and setting k z = k; (1 Equation ( 15) is the modified form of the dispersion relation for the clumpy cloudy medium in the presence of uniform rotation, magnetic fields, and viscoelastic effects.
If the effects of rotation and magnetic field are neglected, we obtain the same dispersion relation as obtained by Dhiman and Mahajan (2023) neglecting dissipative effects in that case.Further, if we neglect the viscoelastic coefficients (ζ = μ = 0), roational frequency ( z = 0), magnetic field intensity (H 0 = 0), the equation ( 15) reduces to the same dispersion relation obtained by Elmegreen (1989).Also, if the clump stirring processes and cloud-cloud collisions are neglected (i.e.= e = ω c = 0), then we obtain the same dispersion relation ( 17) obtained by Prajapati and Chhajlani (2013) for an infinitely conducting medium without rotation.
The dispersion relation ( 15) is useful to discuss the Jeans instability in a molecular cloud clump undergoing heating and cooling mechanisms.To obtain the instability criterion and effects of the various considered parameters on the growth rate, let us analyze these two factors, separately.

Kinetic limit
The first factor of the dispersion relation ( 15) under the kinetic limit (known as the strongly coupled limit) i.e. τ m ω 1 reduces to where 3 )/ρ 0 τ m is defined as the square of the velocity of the compressional viscoelastic mode.Equation ( 16) clearly represents non-gravitating, rotating, damped viscoelastic, Alfvén wave mode that is dynamically stable.In the absence of rotation (i.e.z = 0), equation ( 16) reduces to The above equation represents the magnetosonic mode in the compressional viscoelastic fluid.It is clear that in the absence of rotation, the usual magnetosonic mode is modified due to the presence of viscoelastic coefficients.Thus, due to the presence of uniform rotation, the phase speed of the MHD mode will be significantly affected.
The second factor of the dispersion relation (16) under the kinetic limit reduces to Equation ( 18) is the modified form of the dispersion relation obtained by Elmegreen (1989) in the presence of viscoelastic effects.If the parameters; e , γ, , and ω c representing the clumpy cloudy medium are neglected, we obtain the same dispersion relation obtained by Janaki et al. (2011).It is very difficult to measure the real values of the viscoelastic coefficients for astrophysical systems.Thus, to study the effects of various parameters on the growth rate of the Jeans instability, we replace v c in terms of the Coulomb coupling parameter j , the ion mass and the ion temperature.Following Dolai and Prajapati (2017), the simplified expression of compressional wave velocity (v c ) in terms of coupling parameter ( j ) is given by Here, we have used the following relation between the memory parameter, coupling parameter and compressibility of the medium Substituting the value of compressional velocity in terms of the coupling parameter from ( 19) in ( 18) and after simplifying the resulting equation, we obtain The dispersion relation ( 21) shows the influence of the Coulomb coupling parameter and heating and cooling functions on the Jeans instability of clumpy MCs.Using the Guillemin (1949) criterion for the sign of roots, we obtain the following condition for the onset of gravitational instability from the constant term of equation ( 21) and is given by Condition ( 22) is termed as Jeans instability criterion.The expression for the critical Jeans wavenumber is given by The above condition of Jeans instability and expression of the Jeans wavenumber show the modified form of the Chandrasekhar (1961) instability criterion c 2 s k 2 − 4πGρ 0 < 0 (or k J < ω j /c s ) and is modified due to the presence of the Coulomb coupling parameter ( j ) and heating rate parameter ( e ).It should be noted that the cooling rate parameter does not affect the Jeans instability criterion or the expression of the Jeans wavenumber.Therefore, it is clear from expression (23) that the Jeans wavenumber decreases with increasing values of j and e in the clumpy cloudy medium, which yields that the coupling effects and heating rate postpone the gravitational collapse.
The above theoretical expressions are used to calculate the length scale and Jeans frequency for the considered systems.We now consider the real physical conditions of the MCs for estimating the fundamental plasma parameters.In the region of the molecular cloud clumps, we consider the following parameters; number density of H 2 molecules n = 10 3 cm −3 (Bergin and Tafalla 2007), mass of H 2 molecules m i = 3.32 × 10 −24 g, ion temperature T i = 10K and Coulomb coupling parameter j = 2.0 (Prajapati 2022).With these values of the parameters, the calculated value of the Jeans frequency is measured to be ω j 0.52 × 10 −13 s −1 , and the corresponding Jeans wavenumber (using expression ( 23)) k j 1 0.46 × 10 −17 cm −1 .Hence, the Jeans length of the molecular clouds λ j 5 0.44 pc which is comparable to the length scale of cores within the cloud clumps (Bergin and Tafalla 2007).It is also observed that the calculated value of the Jeans wavelength is larger than the value calculated by Dhiman and Mahajan (2023) for the clumpy MCs with dissipative effects.
The growth rate of the Jeans instability describes how quickly small perturbations or fluctuations within the cloud grow and amplify over time, leading to the eventual collapse of the cloud.The exact behavior of the growth rate depends on the balance between gravitational forces and pressure forces within the cloud.The growth rate of the Jeans instability depends on several factors, including the density of the gas, heating, and cooling functions, its temperature, and the wavelength of the perturbations.To study the impact of the Coulomb coupling parameter and the heating-cooling rate parameter on the growth rate of instability, let us use the following non-dimensional parameters defined in terms of the Jeans frequency; the dispersion relation ( 21) assumes the following nondimensional form; The normalized growth rate of Jeans instability is calculated from ( 25) for different values of wavenumber and is depicted in the following figures.Figure 1.represents the normalized growth rate of the Jeans instability (Re( ω)) versus the normalized wavenumber ( k), for different values of the Coulomb coupling parameter j = 0.5 for weakly coupled plasma (WCP), 1.0 (for viscoelastic fluid), and 2.0 for strongly coupled plasma (SCP).It is clear from the curve corresponding to j = 0.5 that due to the weak coupling between the plasma particles, the growth rate for WCP is observed to be larger than that of SCP ( j = 2.0) and viscoelastic fluid ( j = 1.0).Further, the growth rate of the instability region is decreased due to an increase in the coupling parameter, and the region is further decreased in the case of SCP with j = 2.0.This is because it represents high-frequency oscillations with strong coupling between the plasma particles.The Coulomb coupling parameter thus stabilizes the growth rate of the Jeans instability.
Figure 2. depicts the change of the normalized growth rate with normalized wavenumber for the different values of the cooling rate parameter; ωc = 0.3, 0.6, 1.0, in the presence ( e = 1) of cloud stirring of the molecular cloud.The variation reveals that the growth rate increases with increasing values of the cooling rate ωc , which means that the onset of gravitational instability sets early in the diffused molecular cloud, as compared to dense clouds, which can be validated from the analysis of Elmegreen (1989).Hence, the cooling rate parameter ( ωc ) shows destabilizing character against the self-gravitating collapse dynamics by increasing the rate at which structure fomation takes place.It is also clear that the cutoff wavenumber at which the growth rate becomes zero remains the same for each value of the cooling parameter ( ωc ).
In Fig. 3., the growth rate of Jeans instability is depicted for the different values of the heating rate parameter, e = 0.0, 1.0, 2.0.From the figure, it is observed that the values of the normalized growth rate of gravitational insta-Fig.3 The normalized growth rate of Jeans instability (Re( ω)) versus normalized wavenumber k for different values of heating rate e = 0, 1.0, 2.0 for fixed values of j = 2.0, ωc = 0.2 and γ = 5/3 bility are higher in the absence of heating ( e = 0.0) rate parameter as compared to the values in the presence of heating ( e = 1.0, 2.0) rate parameter.Thus, it is concluded that the heating rate, which arises due to dense cloud stirring, has a stabilizing influence on gravitational instability.In this case, the cutoff wavenumbers at which the growth rate becomes zero are different for the different values of e .In the absence of the heating rate parameter ( e = 0.0) the cutoff wavenumber is much larger than that in the remaining cases.The suppression of the growth rate with e = 0 lies in the range of finite perturbation wavenumbers.Thus, we conclude that both coupling and heating rate parameters delay the onset of collapse of MCs and the structure formation.

Hydrodynamic limit
Let us now discuss the dispersion properties in the hydrodynamic limit or weakly coupled limit corresponding to τ m ω 1.The weakly coupled limit refers to a regime where the interactions between individual particles in the plasma are relatively weak compared to their kinetic energies.In this limit, the first factor of the dispersion relation (15) reduces to Equation ( 26) represents a non-gravitating stable mode that includes the effect of rotational frequency, Alfvén wave velocity, and viscoelastic coefficients.It is independent of the heating and cooling parameters of the clumpy MCs.The stability of the system represented by the above equation can be discussed using the Routh-Hurwitz criterion.Accordingly, all the coefficients of the polynomial must be positive, which satisfies the necessary condition of stability.To satisfy the sufficient condition, all the principal diagonal minors must be positive.We calculate the principal diagonal minors, which are given by; We found that the values all the are positive, hence the system represented by equation ( 26) is a stable system.
The second factor of dispersion relation (15) under the weakly coupled limit reduces to Clearly, equation ( 27) represents the dispersion relation for weakly coupled plasma in a rotating and magnetized clumpy molecular cloud, which is modified due to the presence of viscoelastic coefficients, and heating and cooling functions.The onset of gravitational instability occurs when, Furthermore, the value of the critical Jeans wavenumber is given by We find that the above Jeans instability criterion is independent of the cooling and Coulomb coupling parameters and is different from the instability criterion obtained in the kinetic limit.In order to observe the Jeans instability and get real values of the critical wavenumber, 2 e − 1 > 0 or e > 0.5.Elmegreen (1989) suggested the significance of the heating rate parameter in the cloud clump formation.He mentioned that regions without much star formation activity (i.e. e = 0) will spontaneously clump into cloud complexes on a variety of scales.This scale ranges from the collisional mean free path up to the scale of the conventional Jeans instability in the interstellar gas disk.The smallscale clumping occurs at approximately the same rate as the large-scale clumping if the component clouds are mildly self-gravitating with ω c ω j ∼ 1.Once star formation begins and energy is available to stir the cloud population, thermal equilibrium becomes possible ( e = 1).Under these circumstances, the Jeans instability is observed, which grows continuously for sufficiently large molecular cloud clumps.
Let us now illustrate the growth rate of the Jeans instability in the hydrodynamic limit.The dimensionless form of dispersion relation ( 27) is derived using the non-dimensional parameters; ( ω = ω/ω j , k = kc/ω j , τm = τ m ω j and ωc = ω c /ω j ), as; ω3 + 3 ωc + 0.16 j τm k2 ω2 + k2 (γ + 0.48 j τm ωc ) Figure 4. shows the variation in normalized growth rate (Re( ω)) against normalized wave number k for various values of cooling rate parameter; ωc = 0.3, 0.6, 1.0, in the presence of heating function e = 1.0 and relaxation time τm = 1.0.The nature of the curves indicates that the growth rate increases with increasing values of the cooling rate ωc , which is the same as that obtained in Fig. 2. Notably, in the kinetic limit (Fig. 2.), the cutoff wavenumber at which the growth rates become zero is kc1 = 1.25, while in the hydrodynamic limit, it has been increased to kc2 = 1.75.Thus, the cooling function ( ωc ) shows destabilizing behaviour on the growth rate of self-gravitational instability and encourages the gravitational collapse of clouds to start the structure formation process.
Figure 5. represents the variation in the normalized growth rate (Re( ω)) versus the normalized wavenumber k for different values of heating rate e = 0.0, 1.0, 2.0 in the presence of the memory parameter ( τm = 1.0) and cooling Fig. 6 Comparison of the normalized growth rate of Jeans instability (Re( ω)) versus normalized wave number k between kinetic limit (with j = 2.0) and hydrodynamic limit (with j = 0.5) for fixed values of ωc = 0.2, e = 1.0, τm = 1.0 and γ = 5/3 rate parameter ( ωc = 0.2).The curves show that the growth rate of gravitational instability decreases with the increasing values of heating rate and wavenumber.Hence, it stabilizes the collapse dynamics of the clumpy MCs.Additionally, in the presence and absence of cloud stirring, the cutoff wavenumber is very large in the hydrodynamic regime as compared to the kinetic regime.
The normalized growth rate of Jeans instability in the hydrodynamic regime is calculated from (30) for different values of wavenumber and depicted in figures.In Fig. 6., the normalized growth rate of the Jeans instability (Re( ω)) is plotted against the normalized wavenumber ( k) using ( 25) and (30) for kinetic and hydrodynamic limits, respectively.From the curves, it is clear that the growth rate is larger for WCP with the Coulomb coupling parameter j = 0.5 as compared to SCP with the Coulomb coupling parameter j = 2.0 in the presence of cooling and heating func-tions.This might be because weakly coupled limits exhibit a slower decrease in the growth rate of unstable Jeans modes than strongly coupled limits.The kinetic limit corresponds to the high-frequency limit (ωτ m 1) thus, the growth rate decreases faster than the low-frequency hydrodynamic limit.

Conclusions
The radiative heating and cooling mechanisms play a crucial role in the energy transfer and subsequent gravitational collapse in dense clumps of MCs.This paper studies the influence of heating and cooling parameters on the gravitational instability of rotating and magnetized strongly coupled plasma in the clumps of MCs.A dispersion relation for the considered configuration has been derived analytically using the normal mode analysis and discussed in hydrodynamic and kinetic limits.In both cases, the Jeans criteria are modified due to the presence of the Coloumb coupling parameter and heating rate, but uniform rotation and magnetic field do not affect the instability criterion.However, the presence of magnetic field and uniform rotation significantly modify the compressional Alfvén viscoelastic mode.Furthermore, we calculated the growth rate of instability against the wavenumber for both cases and found that the normalized growth rate decreased with increasing wavenumber.
The graphical illustrations show that the Coulomb coupling parameter ( j ) and heating rate parameter ( e ) stabilizes the growth rate of Jeans instability.The cooling rate parameter (ω c ) shows a destabilizing character against the self-gravitating collapse dynamics.The dynamical stability of the system has been discussed using the Routh-Hurwitz criterion.The variation of growth rate is also plotted in the hydrodynamic limit to study the effects of heating, and cooling parameter and to compare the growth rates with kinetic limits.We found that the growth rate decreases faster in the kinetic limits than in the hydrodynamic limit.Thus, the collapse rate is slower in the hydrodynamic or weakly coupled limit.
The role of strong coupling effects is ubiquitous in many astrophysical systems, such as white dwarfs, neutron stars and in some regions of MCs, as they exist in strongly coupled states.The numerical parameters are calculated for clumpy MCs and measure the Jeans frequency to be ω j 0.52 × 10 −13 s −1 , Jeans wavenumber k j 1 0.46 × 10 −17 cm −1 .The Jeans length for the MCs is calculated to be λ j 5 0.44 pc which is comparable to the length scale of cores within cloud clumps.The present results are helpful to discuss the gravitational collapse in dense clumps of MCs consisting of strongly coupled magnetized plasma with rotation.The inclusion of finite electrical resistivity and the Hall effect is a possible future scope of the present work.