Abstract
We theoretically study the nonlinear pulsational mode dynamics in strongly correlated (viscoelastic) self-gravitating complex charge-fluctuating dust molecular clouds (DMCs) on the astrophysical spatiotemporal scales. A nonlinear normal mode (local) analysis results in a unique pair of extended Kortweg de-Vries-Burgers (KdV-B) equations on the conjugational gravito-electrostatic potential fluctuations in a mixed form. The KdV-B system is numerically analysed in a wide-range parametric window relevant to realistic astronomical DMC circumstances. It is found that the fluctuation dynamics evolves as solitary-chain patterns. The electrostatic fluctuation amplitude increases with the referral frame velocity; whereas, the gravitational fluctuations are insensitive to this velocity. The fluctuation dynamics are found to be independent of variation in the equilibrium dust density, equilibrium dust charge, dust mass, and so forth. The validation and reliability checkup of our results is highlighted in fair corroboration with the literature. Our results could be useful in understanding the mechanism behind bounded structure formation via the non-local self gravitational collapse dynamics.
1. Introduction
The fact that intense locations of the interstellar dense molecular clouds, which are basically self-gravitating partially ionized plasmas, are active sites triggering the formation of stars, planets and other bounded structures in galaxies, is now widely acknowledged [Citation1–4]. These relatively complex formation processes are prompted by a wide variety of collective waves, instabilities and oscillations [Citation5–10]. Various studies [Citation11–14] have brought to light that partially charged clouds display a comparatively new type of collective mode, christened as the pulsational mode.
The dusty plasma plays a pivotal role in astrophysical bodies and space environments like planetary ring systems, cometary tails, interstellar and circumstellar clouds and asteroid zones [Citation15] apart from laboratory plasmas [Citation16]. The ubiquitous nature of the dusty plasmas and its vital role in the space environment results in new interesting salient features effected by dust-charged grains [Citation15,Citation17–21]. Astronomical and satellite observations [Citation17,Citation22,Citation23] have brought to the fore the presence of the massive dust grains in the interstellar dust molecular clouds, in the vicinity of comets and other astrophysical objects. The presence of massive grains in astrophysical plasmas results in the considerable importance of the gravitational effects. The self-gravitational and electrostatic repulsive effects balance each other [Citation12,Citation24,Citation25] in plasmas with medium-sized charged grains.
Numerous studies on nonlinear waves in multispecies self-gravitating dusty plasmas have been carried out over the last few decades in the light of the significant role played by them in the formation processes of varied galactic structures [Citation15,Citation26,Citation27]. As a consequence of the nonlinear and dispersive properties of the considered self-gravitating medium [Citation27], Korteweg–de Vries (K-dV) type of equations and their derivative forms usually govern the evolution of such waves in a nonlinear regime. It is a well-known fact that there exists a wide spectrum of different nonlinear waves and collective eigen modes in star-forming dusty molecular clouds. The excitation of perturbations and waves in such interstellar fluid media saturates into various eigen mode structures like shock, solitons and vertices [Citation9,Citation12–14,Citation28,Citation29]. Karmakar et al. [Citation30] developed a theoretical model to investigate the formation of nonlinear conjugational eigen structure in complex magnetized collisional astrofluid and showed that the fluid perturbation dynamics is administered by a unique bimodal pair forced Korteweg-de Vries (f-KdV) equations. Gohain and Karmakar [Citation31] have come up with an analytical model to investigate the weakly nonlinear gravito-electrostatic waves in a field-free polytropic dust molecular cloud. Later, a classical formalism for the weakly nonlinear instability analysis of a gravitating rotating viscoelastic gaseous cloud in the presence of gyratory dark matter [Citation32].
Of late, Dutta and Karmakar [Citation33] have come up with an evolutionary model analysis with a view to studying the linear pulsational mode dynamics in viscoelastic complex self-gravitating clouds. Various stabilizing and destabilizing agencies have been identified and characterized systematically. However, it has completely ignored any kind of quassi-linear or nonlinear gravito-electrostatic mode-mode coupling and evolution thereof in such visco-elastic media. Subsequently, the nonlinear fluctuation dynamics in the case of self-gravitating nonthermal viscoelastic complex fluid was also investigated in a semi-analytical manner [Citation34]. Later El-Taibany et al. [Citation35] have carried out dust-acoustic solitary and periodic wave investigations in a magnetized self-gravitating dusty plasma. In the recent past, Yang et al.[Citation36] have studied the behaviours of Jeans modes for gravitational systems composed of dark and baryonic matters.
Motivated by this property of the star-forming nonlinear molecular clouds, we herein propose a theoretical model formalism to investigate the nonlinear dynamics of collective pulsational mode excitation and propagation contributing to the triggered initiation of structure formalism in such complex viscoelastic fluids with full incorporation of dust charge fluctuations. Application of multiple scaling techniques in the molecular cloud results in a conjugational pair extended KdV-B equation of complex non-integrable mathematical construct. The KdV-B system is numerically integrated for the exact solutions evolving as diversified eigen – patterns of broadband K-dV family. The key factors, stabilizing and destabilizing the perturbed cloud dynamics, are investigated in detail followed by a realistic comparison with the real astronomical circumstances as already reported in the literature.
2. Basic governing equations
The model under consideration encompasses a four-component isothermal strongly coupled dusty plasma consisting of thermal electrons, ions and inertial dust grains with partial ionization. The comparatively lighter electrons and ions, with the non-consideration of inertial effects, are treated to be Boltzmannian in case of the low-frequency fluctuations on the gravitational scales of space and time; on the contrary, the heavier dust species under the impact of high inertia are taken as fluids. The dust grains are deemed to be viscoelastic fluids because of their low thermal velocities and as a consequence, the mean-free path of the gravitationally interacting grains is practically negligible. It is evident from observation that the massive grains in the cloud which interact gravitationally under the impact of long-range Newtonian coupling processes, add to the net mass of the cloud. The random interaction of the electron and ion thermal currents on the surface of the dust grains results in their charge fluctuation [Citation12], which is taken into consideration. Moreover, the effects of collisional momentum transfer among all the constituent species are also taken care of.
The classical dynamics of the dust cloud constituents, i.e. thermal electrons and ions are controlled by the equations of continuity and momentum respectively [Citation9–12] in flat space-time under non-relativistic conditions as (1) (1) (2) (2) where , , , and represent respectively the equilibrium number density, non-equilibrium number density, mass, velocity and temperature of the jth particles (j = e, i where “e” and “i” stands for electrons and ions respectively). stand for the collision frequencies of the j-th particles with the charged dust grains and ϕ denotes the electrostatic plasma potential.
The dynamics of the viscoelastic fluid comprising of uncharged dust grains is taken care of by the equation of continuity and viscoelastic momentum equation which respectively are (3) (3) (4) (4) On the other hand, the dynamics of the viscoelastic fluid consisting of the charged dust grains is described by the equation of continuity and viscoelastic momentum equations as follows (5) (5) (6) (6) Here , , and stand for population density, mass, flow velocity and temperature of neutral (charged) dust grains respectively [Citation33]. is the charge of the dust grain. Moreover, , , represent respectively viscoelastic relaxation time, shear viscosity co-efficient and bulk viscosity co-efficient of neutral (charged) dust particles. and designate the neutral-charged and charged-neutral collision frequencies respectively. ϕ and ψ respectively denote the electrostatic potential and self-gravitating potential generated by fluid density fields.
Poisson equations for these potential distributions are respectively given by (7) (7) (8) (8) where is the universal gravitational constant.
Based on the assumption that , Equation (Equation8(8) (8) ) takes the form (9) (9) The charge fluctuation equation which dictates the dynamics of the fluctuating dust charge is (10) (10) It may be noted that Equation (Equation1(1) (1) ) represents the equation of continuity for the flux conservation of the thermal species (electrons and ions). It is not a conservation law for mass but for population density. As this model considers charge fluctuating dust molecular clouds (DMCs), therefore the net flux conservation is in correlation with the fluctuating charge dynamical Equation (Equation10(10) (10) ). As a consequence, the net results revealed by Equation (Equation1(1) (1) ) in association with Equation (Equation10(10) (10) ) can represent the conservation laws of thermal particle population density and dust charge conservation laws, respectively. If the dust charge is conserved to be static, then the RHS of Equation (Equation1(1) (1) ) would be exactly zero, thereby representing incompressible fluids associated with the thermal species.
It is clearly evident from Equations (Equation1(1) (1) )–(Equation10(10) (10) ) that the plasma thermal species (hot, mutually thermalized) constitute here inertialess fluids (Equations (Equation1(1) (1) )–(Equation2(2) (2) )) and cold dusty species form the inertial fluids (Equations (Equation3(3) (3) )–(Equation6(6) (6) )). In limiting case of , and , all the species can collectively be treated in a common footing of Boltzmannian equilibrium. Now, we adopt a standard astrophysical normalization procedure [Citation33] as follows
, , (), , (, , ,
where is the Jeans length, is the Jeans time, is the Jeans frequency, is the equilibrium density, is the acoustic phase speed, is the plasma thermal potential and is the equilibrium dust charge. The typical values of the various normalizing parameters adopted here are the same as that already reported in the literature [Citation33].
The normalized form of Equations (Equation1(1) (1) )–(Equation10(10) (10) ) are (11) (11) (12) (12) (13) (13) (14) (14) (15) (15) (16) (16) (17) (17) (18) (18) (19) (19) Here, β = “dn” for neutral dust and “dc” for charged dust.
We propose to study the weakly nonlinear bimodal fluctuation dynamics with a procedurally derived pair extended KdV-B equations. We apply a multi-scale analysis in a new stretched co-ordinate system defined by a pair of coordinates as [Citation29] (20) (20) where ϵ is a smallness (order) parameter, designating the strength of nonlinearity arising due to the fluid convection dynamics. μ is the normalized referral frame velocity in the new coordinate space. The dependent relevant physical variable (Z) is now nonlinearly perturbed and expanded in the various ϵ- powers around the hydrostatic homogeneous cloud equilibrium () as (21) (21) (22) (22) (23) (23) (24) (24) We apply the perturbation series expansion Equations (Equation21(21) (21) )–(Equation24(24) (24) ) in Equations (Equation11(11) (11) )–(Equation19(19) (19) ) for order by order analysis. Then we equate the like terms in various powers of ϵ and after simplification we get:
(25) (25) (26) (26)
with (27) (27) (28) (28) (29) (29) (30) (30) Equation (Equation17(17) (17) ) leads to: (31) (31) Equation (Equation18(18) (18) ) leads to: (32) (32) Finally we derive the pair extended KdV-B equations as: (33) (33) where, (34) (34) and (35) (35) and (36) (36) where
(37) (37) (37) (37) (39) (39) (40) (40) (41) (41) (41) (41) Also (43) (43) (44) (44) (45) (45) (46) (46) (47) (47) (48) (48)
The steady-state pulsational fluctuations can now be well realized by transforming Equations (Equation33(33) (33) ) and (Equation35(35) (35) ) into a time-stationary form of ordinary differential equations (ODEs) in a comoving Galilean frame as (49) (49) (50) (50) It is seen that Equations (Equation49(49) (49) ) and (Equation50(50) (50) ) are analytically non-integrable due to the involvement of non-linear complications. A numerical platform is therefore constructed for analysing the patterns of the fluctuations of interest.
3. Results and discussions
The presented work explores the nonlinear gravito-electrostatic (pulsational) properties of a complex astrophysical viscoelastic cloud with the help of derived extended pair KdV-B equations. The various multi-parametric inputs for our numerical analysis are taken from the different sources available in the literature [5, 30, 33, 34]. The different input initial values used are , , , , , . The results so obtained are focally indicative of diversified interstellar modified wave spectra which are graphically displayed in Figures .
In Figure , we show the profile of the normalized (a) electrostatic potential () varying with the normalized distance (ρ ), (b) potential gradient () with the normalized distance (ρ ), (c) potential gradient () over the potential (), and (d) potential curvature () with the normalized distance (ρ ). Various lines refer to (i) (blue curve), (ii) (red curve), (iii) (black curve). The other input parameters kept fixed are , , , , , , , , T = 1.0. It is found that the electrostatic fluctuations evolve as hybrid eigen-wave structures composed of compressive and rarefactive soliton pairs with amplitude increasing with referral frame velocity (Figure (a)). The corresponding amplitude gradients evolve as similar hybrid patterns with opposite phase but with wave amplitude varying negligibly with the referral frame velocity (Figure (b)). The nonlinear wave evolution, despite the Burger terms, evolve as conservative systems as evidenced confirmatorily with the help of phase space geometrical trajectories (Figure (c)). The wave amplitude curvature evolves as admixtures of compressive and rarefactive solitonic structures with wave amplitude decreasing with the referral frame velocity (Figure (d)), and so forth. The basic physical mechanism behind such fluctuation behaviour may be attributable to the complex non-local long-range gravito-electrostatic interplay on the astrophysical fluid scales of space and time as already reported in the literature [Citation5,Citation30,Citation33,Citation34].
In Figure , we plot the profile of the normalized (a) gravitational potential () varying with the normalized distance (ρ), (b) gravitational gradient () with the normalized distance (ρ ), (c) gravitational gradient () over the potential (), and (d) gravitational potential curvature () with the normalized distance (ρ). The reference of various lines as well as the different input parameters are the same as in Figure . From Figure , we find that gravitational fluctuations evolve as a typical hybrid structures composed of compressive and rarefactive soliton admixtures. The resultant wave amplitude is found to decrease with the referral frame velocity μ and vice versa (Figure (a)). The corresponding field fluctuations evolve accordingly but with an opposite phase (Figure (b)). This fluctuation evolution is still a conservative system indicated with the help of closed phase space (geometrical) trajectories (Figure (c)). The gravitational potential curvature also evolves as hybrid structures composed of compressive and rarefactive solitonic disturbances. The overall potential curvature amplitude is found to enhance with the referral frame velocity (Figure (d)).
In Figure , we depict the same as in Figure , but now for different values. Various lines refer to (i) (blue curve), (ii) (red curve), (iii) (black curve).The other input parameters kept fixed are , , , , , , , , T = 1.0. Here we find the same fluctuation scenarios as in Figure but a notable feature in this case is that the fluctuation dynamics are oblivious of variation in the equilibrium dust density. As a result, the expected distinct lines are no more distinguishable and so forth. The microphysical insight [Citation37] behind this feature lies in the gravity-induced electrostatic polarization effects on the astrophysical scales.
Figure is similar to Figure , but with different values. The reference of the various lines and the input parameters used are analogous to those of Figure . It is evident from this figure that the gravitational fluctuation amplitude decreases with . Furthermore the potential gradient also withnesses a downward slide with increase of equillibrium dust density.
Figure depicts the same as Figure , but for different values. Various lines refer to (i) (blue curve), (ii) (red curve), (iii) (black curve). The other input parameters kept fixed are , , , , , , , , T = 1.0. It is clear from this figure that the equilibrium grain charge has no role to play in modulating the wave amplitude of the electrostatic fluctuations under the action of gravito-electrostatic coupling. It might happen only because of the consideration of weak nonlinearity (negligible wave steepening with ) in our weakly nonlinear perturbative analysis (with very smallness order parameter, ).
Figure illustrates the same as Figure but with different values of . The reference of the different lines and the input parameters made use of here are similar to those of Figure . This figure reveals that the gravitational fluctuation amplitude decreases with .
Figure is analogous to Figure with different values. Various lines refer to (i) (blue curve), (ii) (red curve), (iii) (black curve). The other input parameters kept fixed are , , , , , , , , T = 1.0. The results obtained in this case are also akin to those of Figure . The fluctuation dynamics is found to be independent of variation in the dust mass. The physical mechanism responsible for this is attributable to gravity-induced electrostatic polarization effects as already discussed in the context of Figure .
Figure is similar to Figure but with various values of . The values of in the different coloured lines and the input parameters are as those in Figure . It is seen from the figure that gravitational fluctuation amplitude increases with .
The phase portraits presented here (Figures (c)–(c)) represent the geometrical trajectories of the weakly nonlinear local gravito-electrostatic fluctuations under consideration. This phase space is constructed with the defined phase variables -- potential and potential gradient. The closed nature of the phase trajectories physically ensures the fact that the fluctuation dynamics under consideration are quite conservative in nature in the considered plasma model configuration. It is quite in accordance with the phase plane analyses already extensively reported in the literature [Citation37].
It is worth repetition that the fluctuation dynamics investigated here have a unique feature that the net fluctuation effects are independent of equilibrium dust density, equilibrium dust charge and dust grain mass. The basic physical mechanism behind this is the gravity-induced electrostatic phenomena triggering non-local non-uniform polarization influences on the astrophysical scales of space and time as already discussed above. It is in qualitative accord with the literature [Citation38]. It is in addition to the rest of the figures presented above, although founded in the same physical setting.
4. Conclusions
We report a theoretical model development to speculate the nonlinear pulsational mode dynamics excitable in star-forming dust molecular clouds exhibiting strong collective cooperative effects. The basic governing equations are accordingly constructed in standard astronomical fluid scales of space and time. Application of multiple scaling techniques therein yields an extended pair KdV-B equation having gravito-electrostatic coupling. A numerical illustrative platform shows that the nonlinear gravito-electrostatic fluctuations evolve as a typical hybrid structures composed of compressive and rarefactive solitonic admixtures. The associated field fluctuations, geometrical trajectories, and the fluctuation curvatures are portrayed illustratively. We see that the electrostatic fluctuation amplitude increases with referral frame velocity (Figure ) but remains totally insensitive to equilibrium dust charge density (Figure ), equilibrium dust charge (Figure ) and dust mass (Figure ), whereas the gravitational amplitude decreases with referral frame velocity (Figure ), equilibrium dust charge density (Figure ), equilibrium dust charge (Figure ) and increases with dust mass (Figure ). The authenticity of our adopted calculation schemes is ensured in light of the similar gravito-electrostatic fluctuation results already reported in the literature [Citation5,Citation30,Citation33,Citation34] The analysis could be helpful in understanding astrophysical structure formation dynamics initiated via material transport processes trigerred by waves and instabilities of the non-local pulsational origin. In addition to the above, it can extensively be useful to understand flow-induced instability dynamics and saturation mechanism in different astrophysical situations in the presence of strong collective correlation effects [Citation39].
Disclosure statement
No potential conflict of interest was reported by the author(s).
References
- Draine BT, Salpeter EE. On the physics of dust grains in hot gas. Astrophys J. 1979;231:77–94.
- Pudritz RE. The stability of molecular clouds. Astrophys J. 1990;350:195–208.
- Spitzer Jr L. Physical processes in the interstellar medium. New York: Wiley; 1998.
- Murray DW, Chang P, Murray NW, et al. Collapse in self- gravitating turbulent fluids. Mon Not R Astron Soc. 2017;465:1316–1335.
- Karmakar PK. Nonlinear stability of pulsational mode of gravitational collapse in self gravitating hydrostatically bounded dust molecular cloud. Pramana J Phys. 2011;76:945–956.
- Pandey BP, Avinash K, Dwivedi CB. Jeans instability of a dusty plasma. Phys Rev E. 1994;49:5599–5606.
- Jeans JH. The stability of a spherical nebula. Philos Trans R Soc A. 1902;199:1–53.
- Yaroshenko VV, Verheest F, Morfill GE. Dust-acoustic waves in collisional dusty plasmas of planetary rings. Astron Astrophys. 2007;461:385–391.
- Haloi A, Karmakar PK. Nonlinear waves in nonuniform complex astrocloud. Contrib Plasma Phys. 2017;57:5–19. DOI: 10.1002/ctpp.201600042.
- Karmakar PK, Haloi A. Nonlinear waves in bipolar complex viscous astroclouds astrophys. Space Sci. 2017;362:94.(13pp). DOI: 10.1007/s10509-017-3067-2.
- Dwivedi CB, Sen AK, Bujarbarua S. Pulsational mode of gravitational collapse and its impact on the star formation. Astron Astrophys. 1999;345:1049–1053.
- Pandey BP, Vranjes J. The pulsational mode in the presence of dust charge fluctuations. Phys Scr. 2002;65:513–517.
- Karmakar PK, Borah B. Nonlinear pulsational eigen modes of a planar collisional dust molecular cloud with grain-charge fluctuation. Eur Phys J D. 2013;67:187. (14pp). DOI: 10.1140/epjd/e2013-40165-7.
- Borah B, Karmakar PK. Pulsational mode fluctuations and their basic conservation laws. Adv Space Res. 2015;55:416–427. DOI: 10.1016/j.asr.2014.09.009.
- Goertz CK. Dusty plasmas in the solar system. Rev Geophys. 1989;27:271–292.
- Barkan A, Merlino RL, Angelo ND. Laboratory observation of the dust-acoustic wave mode. Phys Plasmas. 1995;2:3563–3565. DOI: 10.1063/1.871121.
- Rao NN, Shukla PK, Yu MY. Dust-acoustic waves in dusty plasmas. Planet Space Sci. 1990;38:543–546. DOI: 10.1016/0032-0633(90)90147-I.
- Melandso F, Havnes O. Oscillations and resonances in electrostatically supported dust rings. J Geophys Res. 1991;96:5837–5845. DOI: 10.1029/90JA02360.
- Mendis DA, Rosenberg M. Cosmic dusty plasma. Ann Rev Astron Astrophys. 1994;32:419–463.
- Devi K, Nag A, Paul J, et al. Dynamics of sheath evolution in magnetized charge-fluctuating dusty plasmas. Chin J Phys. 2020;65:405–411. DOI: 10.1016/j.cjph.2020.02.028.
- Nag A, Devi K, Paul J, et al. Effect of levitated dust kinetics in the dynamic evolution of inhomogeneous plasma sheath. Chin J Phys. 2021;73:213–223. DOI: 10.1016/j.cjph.2021.05.016.
- Merlino RL, Heinrich JR, Kim SH, et al. Dusty plasmas: experiments on nonlinear dust acoustic waves, shocks and structures. Plasma Phys Control Fusion. 2012;54:124014(1)–124014(10). DOI:10.1088/0741-3335/54/12/124014.
- Hartquist TW, Pilipp W, Havnes O. Dusty plasmas in interstellar clouds and star forming regions. Astrophys Space Sci. 1997;246:243–289.
- Burman S, Chowdhury AR. Solitary waves in self gravitating dusty plasma. Chaos Sol Fract. 2002;13:973–979. DOI: 10.1016/S0960-0779(01)00044-3.
- Pandey BP, Holst BVD, Vranjes J, et al. Jeans instability of an inhomogeneous streaming dusty plasma. Pramana J Phys. 2003;61:109–120. DOI: 10.1007/BF02704515.
- Ida A, Sanuki H, Todoroki J. An extended K-dV equation for nonlinear magnetosonic wave in a multi-ion plasma. Phys Scr. 1996;53:85–88. DOI: 10.1088/0031-8949/53/1/016.
- Verheest F, Shukla PK. Nonlinear waves in multispecies self- gravitating dusty plasmas. Phys Scr. 1997;55:83–85. DOI: 10.1088/0031-8949/55/1/014.
- Bliokh P, Sinitsin V, Yaroshenko V. Dusty and self-gravitational plasmas in space. Netherlands: Springer; 1995.
- Borah B, Karmakar PK. Nonlinear waves in a self-gravitating charge-varying collisional dust molecular cloud in presence of diverse equilibrium inhomogeneities. Phys Scr. 2014;89:125602(1)–125602(26). DOI: 10.1088/0031-8949/89/12/125602.
- Karmakar PK, Haloi A, Roy S. Bimodal pair f-KdV dynamics in star-forming clouds. EPL. 2018;122:15001(1)–15001(8). DOI: 10.1209/0295-5075/122/15001.
- Gohain M, Karmakar PK. Gravito-electrostatic fluctuations of a polytropic charged dust cloud. Phys Scr. 2014;89:1256049(1)–256049(12).
- Karmakar PK, Das P. Instability analysis of cosmic viscoelastic gyro-gravitating clouds in the presence of dark matter. Astrophys Space Sci. 2017;362:142. (13pp). DOI: 10.1007/s10509-017-3125-9.
- Dutta P, Karmakar PK. Viscoelastic pulsational mode. Astrophys Space Sci. 2017;362:141. (12pp). DOI: 10.1007/s10509-017-3120-1.
- Chakraborty R, Karmakar PK. Solitary-train dynamics in viscoelastic nonthermal astroplasmas. EPL. 2018;123:65005(1)–65005(7). DOI: 10.1209/0295-5075/123/65005.
- El-Taibany WF, EL-Labany SK, El-Helbawy AS, et al. Dust-acoustic solitary, and periodic waves in magnetized self-gravito-electrostatic opposite polarity dusty plasmas. Eur Phys J Plus. 2022;137:261. DOI: 10.1140/epjp/s13360-022-02461-9.
- Yang W, Xiong Y, Chen H, et al. Jeans instability of dark-baryonic matter model in the context of Kaniadakis statistic distribution. J Taibah Univ Sci. 2022;16(1):337–343. DOI: 10.1080/16583655.2022.2058837.
- Borah B, Karmakar PK. A theoretical model for electromagnetic characterization of a spherical dust molecular cloud equilibrium structure. New Astron. 2015;40:49–63.
- Bally J, Harrison ER. The eletrically polarized universe. Astrophys J. 1978;220:743–744.
- Roy S, Bezbaruah T, Karmakar PK. Nonlinear streaming instability in viscoelastic quantum dusty plasmas. Waves Random Complex Media. 2023;1–23. DOI: 10.1080/17455030.2023.2177503.