Skip to content
BY 4.0 license Open Access Published by De Gruyter Open Access February 28, 2022

Small-scale sectorial perturbation modes against the background of a pulsating model of disk-like self-gravitating systems

  • Jakhongir Ganiev EMAIL logo and Salakhutdin Nuritdinov
From the journal Open Astronomy

Abstract

In this work, we consider small-scale sectorial perturbation modes in a disk-like model of a radially nonstationary spiral galaxy in order to study the gravitational instabilities of these modes. Calculations of horizontal sectorial small-scale perturbation modes, such as ( m ; N ) = ( 10 ; 10 ) , (11; 11), (12; 12), (13; 13), (14; 14), and (15; 15), against the background of a nonlinearly nonequilibrium anisotropic model of the self-gravitating disk have been carried out. For each of these perturbation modes, critical diagrams are plotted between the virial parameter and the degree of rotation. The growth rates of instability are calculated and compared for different values of the rotation parameter. The results of calculations and comparisons of instability regions show that with an increase in the degree of small scale, the instability region narrows as the wave numbers increase.

1 Introduction

Gravitational instabilities of disk-like models of self-gravitating systems are of great interest not only for galactic disks but also for accretion ones studied by Lodato (2007), Forgan and Rice (2011a, b), Lodato (2012), Rice (2016), Kratter and Lodato (2016), Paneque-Carreño et al. (2021), and Bethune and Latter (2021).

Both large-scale and various small-scale formations are observed in disk-like galaxies. Small-scale formations of disk-like galaxies are mainly open star clusters (OSC) and molecular clouds (MC). Formation and origin of these structural formations in disk-like self-gravitating systems occur due to gravitational instabilities (Mayer et al. 2016, Sharma 2016, Inoue and Yoshida 2018, Roshan and Rahvar 2019). They determine the global structure of the disk of galaxies. However, there is still no analysis of the problems of their origin and no one has studied in detail small-scale perturbations against the background of nonlinearly non-stationary models of disk-like self-gravitating systems, in particular, for our Galaxy. This raises a number of new questions. Can their global distribution be explained by the corresponding formation theory? It is not clear under what physical conditions these objects can form in disk-like systems, and what are the characteristic times of these phenomena? In addition, it is not clear which objects or structural formations of disk-like subsystems of galaxies are directly related to small-scale disturbances. These questions arise primarily due to the lack of knowledge of small-scale high-order perturbations in specific stationary and non-stationary models.

In the disks under study, self-gravity plays the main role in the studies of Rice et al. (2011), Trova et al. (2014), Meru (2015), and Young and Clarke (2015). The gravitational interaction between different parts of the system compresses matter. This process is called gravitational or Jeans instability, and it was studied by Fridman and Khoprskov (2013). It leads to a redistribution of mass, i.e., in one area of the system, the density increases, in another area, respectively, decreases. Clumps of MC are formed due to gravitational instabilities of gas disk-like galaxies in the studies of Forgan and Rice (2011a, b), Dipierro et al. (2014), and Tsukamoto et al. (2015).

Nowadays, many types of instabilities have already been identified for both equilibrium models (Binney and Tremaine 2008) and nonlinear nonstationary states of disk subsystems, which was considered by Nuritdinov (1993) and Mirtadjieva and Nuritdinov (2012), and the corresponding nonstationary dispersion equations (NDEs) were obtained for each perturbation mode.

2 Pulsating anisotropic model

This model is based on the well-known isotropic nonlinear – pulsating model of Nuritdinov (1993)

(1) Ψ ( r , v r , v , t ) = σ 0 2 π Π 1 Ω 2 1 Ω 2 Π 2 1 r 2 Π 2 ( v r v a ) 2 ( v v b ) 2 ] 1 2 χ ( R r ) .

Here σ 0 = σ ( 0 ; 0 ) , Ω is the dimensionless parameter characterizing the degree of rotation of the disk ( 0 Ω 1 ) . v r and v are radial and tangential particle velocities, respectively. The function Π ( t ) has the meaning of the stretching coefficient and has the following form:

Π ( t ) = ( 1 + λ cos ψ ) ( 1 λ 2 ) 1 , t = ψ + λ sin ψ ( 1 λ 2 ) 3 / 2 .

The model pulsates with an amplitude λ = 1 ( 2 T / U ) 0 , where ( 2 T / U ) 0 is the initial virial parameter.

Using Eq. (1), it is possible to make an anisotropic model by averaging the phase density of the rotation parameter Ω (see Bisnovatyj-Kogan and Zel’dovich 1970, Kalnajs 1972, Fridman and Polyachenko 1984).

(2) Ψ Aniz = 1 + 1 Ψ Izot ρ ( Ω ) d Ω 1 + 1 ρ ( Ω ) d Ω ,

where ρ ( Ω ) is the weight function of Ω . We have constructed an anisotropic model, taking the weight function, for example, in the form ρ ( Ω ) = 2 / π ( 1 Ω 2 ) 1 / 2 . After integration and some transformations, we obtain a nonlinearly pulsating model with the anisotropic velocity diagram as follows:

(3) Ψ Aniz = σ 0 π [ 1 + Ω ( x v y y v x ) ] χ ( ( 1 r 2 / Π 2 ) ( 1 Π 2 v 2 ) Π 2 ( v r v a ) 2 ) .

Next, we can find the NDE of model (3). It can be in the following form:

(4) a Aniz ( ψ ) = 4 N ( N 2 1 ) Π 3 γ m N ( N 1 ) ( N 2 + N m 2 ) D N ¯ + i Ω m ( N 2 + N m 2 ) d N ,

where γ m N = ( N + m 1 ) !! ( N m 1 ) !! ( N + m ) !! ( N m ) !! ,

D N ¯ = ψ E W N 1 d P N d ( cosh ) d ψ 1 , W = 1 + λ cos ψ 1 1 + λ cos ψ , d N = ψ E I N d ψ 1 , E = Π 3 ( ψ 1 ) S ( ψ , ψ 1 ) a ( ψ 1 ) , I N = W N 1 1 N + 2 d 2 P N + 1 d ( cosh ) 2 P N d ( cosh ) t g h , cosh = ( cos ψ + λ ) ( cos ψ 1 + λ ) + ( 1 λ 2 ) sin ψ sin ψ 1 ( 1 + λ cos ψ ) ( 1 + λ cos ψ 1 ) ,

where P N ( cosh ) is the Legendre polynomial, N , m are the radial and the azimuthal wave numbers, respectively.

3 Calculation results

In this section, we explore the results of calculations for six values of the sectorial small-scale oscillation modes taking m = N = 10 , 11 , 12 , 13 , 14 , 15 . The gravitational instabilities of these perturbation modes have been calculated using the NDE expression (4).

For example, NDE for the case m = N = 10 has the following form:

(5) a ( ψ ) 10 ; 10 = 1 ( 1 + λ cos ψ ) 21 K 10 ; 10 ( l τ ) . τ = 0 ; 9 ¯ .

Here

l τ ( ψ ) = ψ Π 3 ( ψ 1 ) S ( ψ , ψ 1 ) D ( ψ 1 ) ( λ + cos ψ 1 ) m 1 τ sin τ ψ 1 d ψ 1 ,

where λ is the amplitude of the disk pulsation. The expression for K 10 ; 10 is given in the Appendix.

Numerical calculation of the NDE shows that for the oscillation mode m = N = 10 , the instability starts from the virial parameter’s ( 2 T / U ) 0 = 0.056 at Ω = 0 and reaches 0.325 when Ω = 1 (Figure 1).

Figure 1 
               Critical dependence of the virial ratio on the rotation parameter for 
                     
                        
                        
                           m
                           =
                           10
                        
                        m=10
                     
                  ; 
                     
                        
                        
                           N
                           =
                           10
                        
                        N=10
                     
                  .
Figure 1

Critical dependence of the virial ratio on the rotation parameter for m = 10 ; N = 10 .

For the perturbation mode m = N = 11 , the calculation of the NDE shows that the instability starts from the virial parameter’s value ( 2 T / U ) 0 = 0.049 and reaches 0.282 (Figure 2).

Figure 2 
               Critical dependence of the virial ratio on the rotation parameter for 
                     
                        
                        
                           m
                           =
                           11
                        
                        m=11
                     
                  ; 
                     
                        
                        
                           N
                           =
                           11
                        
                        N=11
                     
                  .
Figure 2

Critical dependence of the virial ratio on the rotation parameter for m = 11 ; N = 11 .

It is shown in Figure 3 for the case m = N = 12 , the instability starts from the virial parameter’s ( 2 T / U ) 0 = 0.044 and reaches 0.248.

Figure 3 
               Critical dependence of the virial ratio on the rotation parameter for 
                     
                        
                        
                           m
                           =
                           12
                        
                        m=12
                     
                  ; 
                     
                        
                        
                           N
                           =
                           12
                        
                        N=12
                     
                  .
Figure 3

Critical dependence of the virial ratio on the rotation parameter for m = 12 ; N = 12 .

Instability in the fourth case m = N = 13 (Figure 4) starts from the value of the virial parameter’s ( 2 T / U ) 0 = 0.039 and reaches 0.221, respectively.

Figure 4 
               Critical dependence of the virial ratio on the rotation parameter for 
                     
                        
                        
                           m
                           =
                           13
                        
                        m=13
                     
                  ; 
                     
                        
                        
                           N
                           =
                           13
                        
                        N=13
                     
                  .
Figure 4

Critical dependence of the virial ratio on the rotation parameter for m = 13 ; N = 13 .

For the oscillation mode m = N = 14 , the calculation of the NDE shows that the instability starts from the value of the virial parameter’s ( 2 T / U ) 0 = 0.035 and reaches 0.199 (Figure 5).

Figure 5 
               Critical dependence of the virial ratio on the rotation parameter for 
                     
                        
                        
                           m
                           =
                           14
                        
                        m=14
                     
                  ; 
                     
                        
                        
                           N
                           =
                           14
                        
                        N=14
                     
                  .
Figure 5

Critical dependence of the virial ratio on the rotation parameter for m = 14 ; N = 14 .

The instability in the sixth case, m = N = 15 , starts from the virial parameter’s value ( 2 T / U ) 0 = 0.032 at Ω = 0 and reaches 0.180 when Ω = 1 .

The shaded area defines the instability regions for the indicated disturbance modes in the figures (Figures 16). It can be seen from the figures that as the speed of rotation of the model grows, the region of instabilities always grows, and with increasing values of m and N , the region of instability, correspondingly, narrows.

Figure 6 
               Critical dependence of the virial ratio on the rotation parameter for 
                     
                        
                        
                           m
                           =
                           15
                        
                        m=15
                     
                  ; 
                     
                        
                        
                           N
                           =
                           15
                        
                        N=15
                     
                  .
Figure 6

Critical dependence of the virial ratio on the rotation parameter for m = 15 ; N = 15 .

We also calculated and compared the increments of instabilities of small-scale oscillation modes for different values of the parameters (Figures 7 and 8). In Figures 7 and 8 show the increments of instabilities for all small-scale oscillation modes at different values of the virial ratio and degree of rotation. It can be seen from the figures that the curves of the growth rates of instabilities of the studied perturbation modes are arranged in descending order and it should be noted that with an increase in the value of the wavenumbers, the increments of instabilities also increase. The curved lines intersect each other, merging at some values of the increment of instabilities, except for the value of the rotation parameter Ω = 0.2 , without changing their position as the wave numbers increase.

Figure 7 
               (a to c) Comparison of the increments of instabilities of all small-scale oscillation modes for different values of the rotation parameter.
Figure 7

(a to c) Comparison of the increments of instabilities of all small-scale oscillation modes for different values of the rotation parameter.

Figure 8 
               (a to c) Comparison of the increments of instabilities of all small-scale oscillation modes for different values of the rotation parameter.
Figure 8

(a to c) Comparison of the increments of instabilities of all small-scale oscillation modes for different values of the rotation parameter.

4 Conclusion

The gravitational instabilities are deeply studied for oscillation modes m = N = 10 , 11 , 12 , 13 , 14 , 15 in this work. Diagrams of the critical dependence of the virial ratio on the rotation parameter are plotted and the increments of instabilities of all small-scale oscillation modes for various values of the model parameters are compared. Figures 16 show that with an increase in the value of the model rotation speed, an increase in the region of instabilities is observed.

From a comparison of the graphs (Figures 7 and 8), we can see that the intersection of the curves at some values of the instability increments except for the case Ω = 0.2 , while the location of the curves does not change. Note that an increase in the value of wave numbers leads to an increase in the value of the increments of instabilities.

Acknowledgements

This work was partly supported by the Ministry of Innovative Development of the Republic of Uzbekistan.

  1. Funding information: The authors state no funding involved.

  2. Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

  3. Conflict of interest: The authors state no conflict of interest.

Appendix

Substituting all the expressions of the functions in (4) and calculating each of them step by step, we find that

K 10 ; 10 = 4199 65536 55 c 9 495 c 7 e 2 b 2 + 3465 4 c 5 e 4 b 4 5775 16 c 3 e 6 b 6 + 3465 128 c e 8 b 8 l 0 + 1485 c 8 e 2 b 6930 c 6 e 4 b 3 + 51975 8 e 6 b 5 c 4 10395 8 e 8 b 7 c 2 + 3465 128 e 10 b 9 l 1 + 28215 2 e 4 b 2 c 7 530145 16 e 6 b 4 c 5 495 c 9 e 2 + 509355 32 e 8 b 6 c 3 10395 8 e 10 b 8 c l 2 + 234465 4 e 6 b 3 c 6 557865 8 e 8 b 5 c 4 6930 e 4 b c 8 + 509355 32 e 10 b 7 c 2 5775 16 e 12 b 9 l 3 + 7494795 64 e 8 b 4 c 5 557865 8 e 10 b 6 c 3 530145 16 e 6 b 2 c 7 + 51975 8 e 12 b 8 c + 3465 4 c 9 e 4 l 4 + 7494795 64 e 10 b 5 c 4 530145 16 e 12 b 7 c 2 557865 8 e 8 b 3 c 6 + 3465 4 e 14 b 9 + 51975 8 e 6 b c 8 l 5 + 234465 4 e 12 b 6 c 3 6930 e 14 b 8 c 557865 8 e 10 b 4 c 5 + 509355 32 e 8 b 2 c 7 5775 16 c 9 e 6 l 6 + 28215 2 e 14 b 7 c 2 495 e 16 b 9 530145 16 e 12 b 5 c 4 + 509355 32 e 10 b 3 c 6 10395 8 e 8 b c 8 l 7 + 1485 e 16 b 8 c 6930 e 14 b 6 c 3 + 51975 8 e 12 b 4 c 5 10395 8 e 10 b 2 c 7 + 3465 128 c 9 e 8 l 8 + 55 e 18 b 9 495 e 16 b 7 c 2 + 3465 4 e 14 b 5 c 4 5775 16 e 12 b 3 c 6 + 3465 128 e 10 b c 8 l 9

+ 4 i Ω e 495 4 b c 8 + 1155 2 e 2 b 3 c 6 17325 32 e 4 b 5 c 4 + 3465 32 e 6 b 7 c 2 1155 512 e 8 b 9 l 0 + 5445 2 c 7 e 2 b 2 + 495 4 c 9 + 197505 32 c 5 e 4 b 4 93555 32 c 3 e 6 b 6 + 121275 512 c e 8 b 8 l 1 + 5445 2 c 7 e 2 b 2 + 495 4 c 9 + 197505 32 c 5 e 4 b 4 93555 32 c 3 e 6 b 6 + 121275 512 c e 8 b 8 l 2 + 1867635 32 e 6 b 4 c 5 + 308385 16 e 4 b 2 c 7 + 4163775 128 e 8 b 6 c 3 1155 2 c 9 e 2 93555 32 e 10 b 8 c l 3 + 21417165 256 e 8 b 5 c 4 + 1867635 32 e 6 b 3 c 6 + 696465 32 e 10 b 7 c 2 197505 32 e 4 b c 8 17325 32 e 12 b 9 l 4 + 1867635 32 e 10 b 6 c 3 + 21417165 256 e 8 b 4 c 5 + 197505 32 e 12 b 8 c 696465 32 e 6 b 2 c 7 + 17325 32 c 9 e 4 l 5 + 308385 16 e 12 b 7 c 2 + 186763532 e 10 b 5 c 4 + 1155 2 e 14 b 9 4163775 128 e 8 b 3 c 6 + 93555 32 e 6 b c 8 l 6 + 5445 2 e 14 b 8 c + 308385 16 e 12 b 6 c 3 696465 32 e 10 b 4 c 5 + 620235 128 e 8 b 2 c 7 3465 32 c 9 e 6 l 7 + 495 4 e 16 b 9 + 5445 2 e 14 b 7 c 2 19750532 e 12 b 5 c 4 + 93555 32 e 10 b 3 c 6 121275 512 e 8 b c 8 l 8 + 495 4 e 16 b 8 c 1155 2 e 14 b 6 c 3 + 17325 32 e 12 b 4 c 5 3465 32 e 10 b 2 c 7 + 1155 512 c 9 e 8 l 9 ,

where c = λ + cos ψ , b = sin ψ , e = 1 λ 2 .

References

Bethune W, Latter H. 2021. Spiral structures in gravito-turbulent gaseous disks. A&A. 650:20. 10.1051/0004-6361/202040094Search in Google Scholar

Binney J, Tremaine S. 2008. Galactic Dynamics. Second edition ISBN 978-0-691-13026-2 (HB). USA: Princeton University Press. p. 383. 10.1515/9781400828722Search in Google Scholar

Bisnovatyj-Kogan GS, Zel’dovich YaB. 1970. Models of point mass clusters with quadratic gravitational potential. Astrofizika. 6:387–396. 10.1007/BF01002656Search in Google Scholar

Dipierro G, Lodato G, Testi L, de Gregorio Monsalvo I. 2014. How to detect the signatures of self-gravitating circumstellar discs with the Atacama Large Millimeter/sub-millimeter Array. MNRAS. 444:1919–1929. 10.1093/mnras/stu1584Search in Google Scholar

Forgan D, Rice K. 2011a. The nature of angular momentum transport in radiative self-gravitating protostellar discs. MNRAS. 410:994–1006. 10.1111/j.1365-2966.2010.17500.xSearch in Google Scholar

Forgan D, Rice K. 2011b. The Jeans mass as a fundamental measure of self-gravitating disc fragmentation and initial fragment mass. MNRAS. 417:1928–1937. 10.1111/j.1365-2966.2011.19380.xSearch in Google Scholar

Fridman AM, Polyachenko VL. 1984. Physics of gravitating systems II. New-York: Springer-Verlag. 10.1007/978-3-642-87833-6Search in Google Scholar

Fridman AM, Khoprskov AV. 2013. Physics of Galactic Disks. UK: Cambridge International Science Publishing. p. 754. Search in Google Scholar

Inoue Sh, Yoshida N. 2018. Spiral-arm instability: giant clump formation via fragmentation of a galactic spiral arm. MNRAS. 474:3466–3487. 10.1093/mnras/stx2978Search in Google Scholar

Kalnajs AJ. 1972. On the chemical composition and the origin of the globular clusters of our Galaxy. ApJ. 175:63. 10.1086/151538Search in Google Scholar

Kratter K, Lodato G. 2016. Gravitational instabilities in circumstellar disks. Annual Rev Astron Astrophys. 54:271–311. 10.1146/annurev-astro-081915-023307Search in Google Scholar

Lodato G. 2007. Self-gravitating accretion discs. Nuovo Cimento Rivista Serie. 30:293. Search in Google Scholar

Lodato G. 2012. The role of gravitational instabilities in the feeding of supermassive black holes. Adv Astron. 2012:846875. 10.1155/2012/846875Search in Google Scholar

Mayer L, Peters T, Pineda JE, Wadsley J, Rogers P. 2016. Direct detection of precursors of gas giants formed by gravitational instability with the atacama large millimeter/submillimeter array. ApJ. 823:L36. 10.3847/2041-8205/823/2/L36Search in Google Scholar

Meru F. 2015. Triggered fragmentation in self-gravitating discs: forming fragments at small radii. MNARS. 454:2529–2538. 10.1093/mnras/stv2128Search in Google Scholar

Mirtadjieva KT, Nuritdinov SN. 2012. Instabilities in a nonstationary model of self-gravitating disks. IV. Generalization of the models and comparison of results. Astrophys. 55:551–564. 10.1007/s10511-012-9260-7Search in Google Scholar

Nuritdinov S. 1993. Author’s abstract of dissertation for thesis. uch. doctoras degree nat. - mat. sciences. C - Petersburg. Search in Google Scholar

Paneque-Carreño T, Pérez LM, Benisty M, Hall C, Veronesi B, Lodato G, et al. 2021. Spiral arms and a massive dust disk with non-Keplerian kinematics: possible evidence for gravitational instability in the disk of Elias 2–27. ApJ. 914:88. 10.3847/1538-4357/abf243Search in Google Scholar

Rice WK, Armitage PJ, Mamatsashvili GR, Lodato G, Clarke CJ. 2011. Stability of self-gravitating discs under irradiation. MNRAS. 418:1356–1362. 10.1111/j.1365-2966.2011.19586.xSearch in Google Scholar

Rice K. 2016. The evolution of self-gravitating accretion discs. PASA. 33:e012. 10.1017/pasa.2016.12Search in Google Scholar

Roshan M, Rahvar S. 2019. Evolution of spiral galaxies in nonlocal gravity. ApJ. 872:6. 10.3847/1538-4357/aafc2aSearch in Google Scholar

Sharma P. 2016. Self-gravitational instability of dusty plasma with dissipative effects. Astrophys Space Sci. 361:114. 10.1007/s10509-016-2700-9Search in Google Scholar

Trova A, Huré JM, Hersant F. 2014. Self-gravity in thin discs and edge effects: an extension of Paczynski’s approximation. A&A. 563:A132. 10.1051/0004-6361/201322125Search in Google Scholar

Tsukamoto Y, Takahashi SZ, Machida MN, Inutsuka SI. 2015. Effects of radiative transfer on the structure of self-gravitating discs, their fragmentation and the evolution of the fragments. MNRAS. 446:1175–1190. 10.1093/mnras/stu2160Search in Google Scholar

Young MD, Clarke CJ. 2015. Dependence of fragmentation in self-gravitating accretion discs on small-scale structure. MNRAS. 451:3987–3994. 10.1093/mnras/stv1266Search in Google Scholar

Received: 2021-10-31
Accepted: 2022-01-30
Published Online: 2022-02-28

© 2022 Jakhongir Ganiev and Salakhutdin Nuritdinov, published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Downloaded on 28.3.2024 from https://www.degruyter.com/document/doi/10.1515/astro-2022-0011/html
Scroll to top button