Next Article in Journal
Features of Secular Changes in the Flux Density of the Cas a Supernova Remnant, from Centimeter to Decameter Wavelengths
Next Article in Special Issue
Naturally-Coupled Dark Sectors
Previous Article in Journal
Eclipsing Systems with Pulsating Components (Types β Cep, δ Sct, γ Dor or Red Giant) in the Era of High-Accuracy Space Data
Previous Article in Special Issue
History of Astronomy in Australia: Big-Impact Astronomy from World War II until the Lunar Landing (1945–1969)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling of Spiral Structure in a Multi-Component Milky Way-Like Galaxy

by
Sergey Khrapov
1,†,
Alexander Khoperskov
1,† and
Vladimir Korchagin
2,*
1
Volgograd State University, Universitetsky pr., 100, 400062 Volgograd, Russia
2
Institute of Physics, Southern Federal University, Stachki Street 194, 344090 Rostov-on-Don, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Galaxies 2021, 9(2), 29; https://doi.org/10.3390/galaxies9020029
Submission received: 2 March 2021 / Revised: 14 April 2021 / Accepted: 25 April 2021 / Published: 28 April 2021

Abstract

:
Using recent observational data, we construct a set of multi-component equilibrium models of the disk of a Milky Way-like galaxy. The disk dynamics are studied using collisionless-gaseous numerical simulations, based on the joined integration of the equations of motion for the collision-less particles using direct integration of gravitational interaction and the gaseous SPH-particles. We find that after approximately one Gyr, a prominent central bar is formed having a semi-axis length of about three kpc, together with a multi-armed spiral pattern represented by a superposition of m = 2-, 3-, and 4-armed spirals. The spiral structure and the bar exist for at least 3 Gyr in our simulations. The existence of the Milky Way bar imposes limitations on the density distributions in the subsystems of the Milky Way galaxy. We find that a bar does not form if the radial scale length of the density distribution in the disk exceeds 2.6 kpc. As expected, the bar formation is also suppressed by a compact massive stellar bulge. We also demonstrate that the maximum value in the rotation curve of the disk of the Milky Way galaxy, as found in its central regions, is explained by non-circular motion due to the presence of a bar and its orientation relative to an observer.

1. Introduction

Attempts to understand the phenomenon of spiral structure in galaxies have a long history [1,2,3,4,5]. There is a general consensus that the spiral structure is a manifestation of density perturbations propagating in a multi-component stellar-gaseous disk. However, there is not agreement as to the mechanism for spiral generation, one that can successfully explain the rich morphological variety observed in galaxies. To describe the observed patterns, a few mechanisms have been suggested. Some researchers treat the spiral structure as the long-lived global modes that last in galactic disks for tens of galactic rotations [1,3,6]. Others consider the spiral structure to be a transient phenomenon, so that the spiral pattern changes many times during a galaxy’s lifetime [7,8,9,10,11]. The tidal encounter of a small mass companion with a disk galaxy can also be a mechanism for the formation of spiral structure under certain conditions [12]. Khoperskov et al. [13] showed that a nonaxisymmetric halo can generate a large-scale spiral density wave in an embedded stellar disk. These authors showed that the pattern is observed during many revolutions, even if the disk is gravitationally stable.
Recent observations have provided a significant step in understanding the nature of the spiral structure in galaxies. Using stellar-cluster catalogs for three nearby spiral galaxies (NGC 1566, M51a, and NGC 628) from the Legacy ExtraGalactic UV Survey, Shabani et al. [14] measured the gradients in stellar ages across the spiral arms within these galaxies. The authors find that the grand design spiral galaxy NGC 1566 shows a significant age gradient across the spiral arms, consistent “with the prediction of a stationary density wave theory”. The authors attribute the presence of a global mode in this galaxy to a “strong bar”. While an age gradient of the star clusters across the spiral arms in the other two galaxies (M51a and NGC 628) was not confirmed by [14], the authors do suggest that the spiral structure in M51a is a result of tidal interaction. Peterken et al. published, almost simultaneously, a paper [15] entitled “A direct test of density wave theory in a grand-design spiral galaxy”. They used the SDSS MaNGA survey to study an age gradient in the galaxy UGC 3825. They not only determined an age gradient across the arms, but also measured the pattern speed of the spiral arms, finding that it varies little with radius. The authors concluded that their observations are consistent with the presence of a quasi-stationary density wave in UGC 3825. Recently, Bialopetravicius & Narbutis [16] confirmed the existence of an age gradient across the spiral arms in the galaxy M83. Abdeen et al. [17] looked for an age gradient across the spiral arms in fifteen nearby spiral galaxies and confirmed the existence of a gradient in all of them. Remarkably, Abdeen et al. [17] found an age gradient in the galaxy NGC 5194 (M51a). This galaxy is a non-barred galaxy of Sb-type and, thus, the galaxy shows the presence of a global density perturbation not associated with a bar.
Our knowledge of the spiral structure of the Milky Way galaxy is more limited. The morphology of the Milky-Way spiral structure (the number of arms, pitch angle, position(s) of corotation resonance(s)) remains obscured because of our vantage point, i.e., we view the Galactic disk and its spiral structure edge-on [18,19,20,21,22]. Difficulties in observing the far side of the Galactic disk, primarily due to extinction, prevent us from obtaining reliable information regarding the morphology and kinematics of nearly half the Milky Way disk. The use of radio interferometry has allowed some initial progress in the determination of the properties of the spiral structure on the opposite side of the Galactic disk [23].
Masers are the most important source of information for the determination of the morphology of the spiral pattern in the Milky Way galaxy [24]. Complementing the maser data, the GAIA astrometric mission [25,26,27,28,29] has recently provided positions and kinematics for more than a billion stars in the Milky Way. This will enable us to more reliably infer the morphology of the spiral pattern in our own galaxy. It is important to note that the GAIA data on Milky Way spiral structure in the solar neighborhood agrees with the results that are based on VLBI maser observations [26]. The third GAIA release [29] will give new impetus in understanding the Milky Way spiral structure.
The structure of the gas flow in the stellar bar area is complex [30,31,32,33,34]. For this reason, the ability to reproduce the kinematics of features in this region is a useful touchstone for verifying theoretical models. The asymmetry of the central gas layer can be caused either by a lopsided nuclear bar or by infalling gas flow [30]. The observed galactic HI and CO l-v diagrams demonstrate a clumpy morphology, due to the multiphase interstellar medium, gravitational instability at different scales, and star-formation and supernovae feedback [31]. The bulge of our galaxy shows some significant differences with a classical bulge, which formed as a result of major mergers of pre-existing disks [33]. A series of simulations with a huge number (eight billion equal-mass particles) of particles provides the initial conditions for the dynamic models that are required to reproduce the observed characteristics of the Milky Way [34].
Korchagin et al. [35] modeled the dynamics of a two-component stellar-gaseous Milky Way disk using 2D-simulations. These authors found that a three-armed spiral pattern is generated in the Milky Way disk, and its spiral structure is sustained during at least 3 Gyr.
Here, we present the results of three-dimensional simulations of the dynamics of a Milky Way stellar-gaseous disk. We model the dynamics of three-dimensional multicomponent disks using a set of equilibrium models constrained by observational errors and demonstrate that, in most of the models, a multi-armed spiral pattern (a superposition of spirals with different azimuthal wavenumbers) is generated in the disk.

2. Observational Data

2.1. Density Distribution of Stars

We use the commonly accepted exponential density distribution to approximate the surface density distribution of stars in the disk of the Milky Way
σ * ( r ) = σ 0 exp { r / r d } ,
where σ 0 is the central surface density and r d is the radial exponential scale of the stellar disk.
Using COBE/DIRBE observational data Drimmel & Spergel [36] found that the distribution of stars in the Milky Way disk can be approximated by an exponential function with radial scale length 0.28 R , where R is the distance of the Sun from the center of the Galaxy. Adopting the distance from the Sun to the center of galaxy to be 8 kpc leads to a galactic disk radial scale length of 2.24 kpc. Recent estimates of the radial scale length for Milky Way old stellar populations within 4 < r < 15 (kpc) give a value of 2.2 ± 0.2 kpc [37]. From a study of the density distribution of red clump stars outside the solar circle [38], find that the radial scale length of the old stellar disk has a value of r d = 2.37 ± 0.02 kpc. We assume the density distribution of the Milky Way old stellar disk varies exponentially with radius, having a scale length of 2.25 kpc. However, it is necessary to stress that estimates of the radial scale length of the Milky Way disk vary over a rather large range r d = (2–4) kpc [36,39,40], and the references in these papers. The value of the central surface density of the disk, σ 0 , is very sensitive to uncertainty in the radial scale length r d (1). This, in turn, imposes limitations on the density distributions of the central bulge and of the halo of the Milky Way disk.
The two key parameters that determine the disk stellar density distribution are the radial scale length of the old stellar disk and the stellar surface density in the solar neighborhood [41]. Estimates that were undertaken by different authors give values varying over a rather wide range. We adopt for the surface density of stars in the solar neighborhood a value of 33.4 ± 3 M / pc 2 from the recent paper by McKee et al. [42]. With the these parameters, the total mass of the old stellar disk of the Milky Way is approxiamtely (4–5) × 10 10 M .
Observations of edge-on galaxies show that the vertical density distribution of stellar disks can be approximated by functions ρ * ( z ) exp ( z / h e ) , sech ( z / h 1 ) , sech 2 ( z / h * ) or their combinations [43,44,45], where h e , h 1 , and h * are the vertical scale heights. The vertical scale height of a stellar disk often depends on radial coordinate, which complicates the situation. In a simplest model of a one-component isothermal disk, the solution is [46,47]
ρ ( r , z ) = ρ 0 ( r ) sech 2 ( z / h ) ,
where h is the vertical scale height of the disk and ϱ 0 is the central volume density.
Deviations from (2) are determined by large number of factors, such as an inhomogeneity of disk velocity dispersion in vertical direction, multi-component nature of the disk that reflects its chemical and dynamical evolution, etc. An additional factor causing deviation from (2) is a dark matter halo, as seen from the equilibrium Jeans equation in vertical direction [46,47]:
d 2 ρ d z 2 + 2 d ln c z d z d ρ d z 1 ρ ( d ρ d z ) 2 + 4 π G c z 2 σ 2 z 0 ρ [ ρ ρ 0 + 2 z 0 σ { ρ h 1 4 π G r V c 2 r } ] = 0 ,
where G is the gravitational constant, ρ ( z ) is the volume mass density, ρ 0 = ρ ( z = 0 ) , σ is the surface density, z 0 = 0 ( ρ ( z ) / ρ 0 ) d z , V c is the circular velocity in the plane z = 0 , ρ h is the local halo density, and c z is the velocity dispersion of stars in the vertical direction. The solution (2) corresponds to the case c z = const , so that the expression in curly brackets is equal to zero. The initial equilibrium of the disk in our simulations was built using the iteration procedure that is described in [47], which was modified for presence of a thin gaseous component.
Juric et al. [48] find that the three-dimensional distribution of old stellar population of Milky Way can be represented by two subsystems: the thin disk with a radial scale length about 2.5 kpc and a vertical scale height of 300 pc and the thick disk with a scale length and height of 3.6 kpc and 900 pc, respectively.
The vertical scale height of the density distribution of the Milky Way disk is another important parameter influencing the gravitational stability of the system and, thus, in determining the morphology of the spiral pattern. To explore this, in our numerical simulations we vary the vertical scale height of the collisionless disk from 200 to 380 parsec.
The observational data show that the Milky Way scale height increases with radius [49]. Nonetheless, we build our equilibrium disk models with constant vertical scale heights.

2.2. Rotation Curves

Figure 1 shows the observed rotational curve of the Milky Way disk based on various observations. There is fairly good agreement between the different determinations of the rotation curve in the range for r (6–14) kpc. However, the curves differ considerably in the inner regions of the Galaxy.
New kinematical data that are based on maser observations provide independent and promising data on the kinematical properties of the cold component of the Milky Way disk [53,55,56]. However, despite the fact that maser observations directly determine the parallaxes of sources, one has to keep in mind that the reconstruction of the Galactic rotation curve involves a number of assumptions, such as, for example, the model of the Galactic potential, which is necessary to reconstruct the rotation curve [57].
The Galactic rotation curve has a local maximum in the inner disk region, as observed in atomic and molecular hydrogen lines [58,59,60]. The maximum of the rotational velocity is probably associated with the non-circular motion of gas, as caused by the bar, in the central regions of the Galactic disk.
The angular velocity of the stellar disk is usually lower than that of the gaseous disk. Estimates of the rotational velocity of stars in the solar neighborhood vary from 200 km s 1 [61] to 237 km s 1 [62]. A recent estimate by [63] gives the value of the rotational velocity of the old stellar disk (G- and K- dwarfs) as V = 205 ± 20 km s 1 .
Outside the stellar disk, the rotation curve of the Milky Way slowly decreases with a radius reaching about 160 km s 1 at a distance of 100 kpc from the center of the Galaxy [64], which apparently points to the presence of the extended and massive halo of the Galaxy.

2.3. Velocity Dispersion of Stars

Figure 2 shows a summary of data on the observational determinations of the radial dependence of the velocity dispersion of stars c r in the Milky Way.
In general, the observational data can be approximated by the function c r = 91 [km s 1 ] / exp ( r / 7.4 [ kpc ] ) . The velocity dispersion profile that was obtained by Tiede & Terndrup [73] gives good agreement with the data in the disk’s central regions of the Galaxy, but in the solar neighborhood with R = 8 kpc the estimate by Tiede & Terndrup [73] gives the value 17 km s 1 , which is too low, and contradicts observations. The velocity dispersion in the central regions of the Milky Way disk is about c r (110–130) km s 1 [65,70,72] and is likely associated with the bulge stellar populations.
Using the Gaia-ESO spectroscopic survey, [71] determined the velocity dispersions ( c r , c φ , c z ) for the thin and thick disk separately, within 6.5 kpc r 9 kpc. They found, for the thin disk, a velocity dispersion of c r = 33 km s 1 (See points 10 in Figure 2). Piffl et al. [40] used a sample of 200,000 red giant stars from the Gaia-ESO spectroscopic survey and estimated the radial velocity dispersion for the thin disk within 1.5 kpc of the plane to be c r = 34 km s 1 [40]. It is known that the presence of young stellar populations decreases the effective velocity dispersion in a stellar disk. A detailed study of this effect was conducted by Rafikov [74], who discussed the axisymmetric stability criterion of a disk consisting of a few stellar components and a gaseous component. We take this into account by using a multiplicative factor, c r c r / 1.1 , based on a statistically averaged estimate of the velocity dispersion for the multi-component stellar disk (giants, main sequence stars, white and brown dwarfs):
c r ( a v ) = l σ l c l 2 / l σ l 1 / 2 ,
where σ l and c l are the surface density and dispersion velocities for the l-th component, respectively. Using data of [74] for nine types of stellar populations, we obtain the above mentioned decrease in the radial velocity dispersion of the stars.
It is known that, for a collisionless disk in equilibrium, the condition
c φ = c r ϰ 2 Ω
is satisfied ( ϰ is the epicyclic frequency, Ω is the angular speed). We use this condition to determine the initial dependence of the azimuthal velocity dispersion along the radial direction. Our numerical simulations show that this condition is satisfied during the disk evolution.
Different authors provide values of 0.45 0.6 for the ratio c z / c r in the solar neighborhood. A recent estimate by [63] of the velocity dispersions c z = 16 km s 1 , c r = 35 km s 1 , gives a value for c z / c r of 0.46, which is what we used to build our equilibrium models. The ratio c z / c r = 16 km s 1 /35 km s 1 is equal to 0.46, based on SCUSS and SDSS observational data for G and K dwarfs of the thin disk [63]. Estimates by Katz et al. [25] and Adibekyan et al. [75] give c z / c r = 0.5 for this ratio. We choose a value of 0.5 in our simulations for the ratio c z / c r . We note that the kinematical properties of the stellar disk change with time during disk evolution, due to disk heating Hayden et al. [76].

2.4. Stellar Bulge and Bar

A number of observational studies demonstrate that the Milky Way has a boxy/peanutshaped bulge. Recent evidence points to the presence of an additional spheroidal component [77]. The structural properties and origin of the bulge are the subject of intense study. In this paper, we neglect the possible evolution of the bulge and consider it to be a fixed structure of the Galaxy contributing to the total gravitational potential. Estimates of the mass of the Milky Way bulge vary over a rather wide range, (0.9–2.0) × 10 10 M [78,79], which, in part, is a result of different assumptions regarding the stellar density distribution of the bulge. Estimates of the Galactic bulge’s mass are based on photometric and kinematical data. A recent estimate for the bulge, based on photometric data [78], gives the total mass of stars and remnants in the bulge as 2.0 ± 0.3 × 10 10 M , which is consistent with the estimate by Portail et al. [79], who derive a total dynamical mass for the bulge equal of 1.85 ± 0.05 × 10 10 M . The kinematical estimates of the bulge’s mass are based on the assumption that the maximum of the rotational velocity of gas close to the galactic center at 0.5 kpc is caused only by the galactic bulge. This gives a total bulge mass of approximately M b = 0.99 × 10 10 M [80]. We vary the total mass of the bulge in our simulations within (0.7–1.6) × 10 10 M .
The presence of a prominent stellar bar with a large semi-axis of about (3.1–3.5) kpc and with an axial ratio of about 1:0.4 is another peculiarity of the Milky Way galaxy [77]. The major semi-axis of the bar deviates by approximately 23 from the direction to the Galactic center [81]. A successful dynamical model of the Milky Way must agree with the available observational data and reproduce the observed parameters of the bar.

2.5. Gas Distribution

Gas plays a significant, possibly crucial, role in forming and sustaining the spiral structure in galaxies. Goldreich & Lynden-Bell [82] pointed to the importance of gas in the formation of spiral structure: S0 galaxies are topographically similar to normal spirals, but they have no gas, no dust, and no spiral arms. This suggests that stellar dynamics is not solely responsible for arm formation.
The total amount of gas within 30 kpc of the Milky Way disk is approoximately 8 × 10 9 M [83]. In our simulations, we vary the total amount of the disk’s gaseous component from 3.5 × 10 9 M to 6.5 × 10 9 M within the optical radius of the disk, R o p t = 4 r d = 9 kpc, s listed in Table 1. As a result, the observationally limited ratio of masses of gas and stars in the Milky Way disk varies within μ g = M g / M d = (0.13–0.17). In our simulations, the gas surface density in the solar neighborhood was kept within the observed estimates of ( H 2 , H I and H I I ) σ g = 13.7 ± 1.6 M pc 2 [42]. The observational data for distant galaxies indicate a significant role of gas accretion, minor and major mergers, which makes the stellar-gas disk a non-conservative system, exerting a strong influence on the kinematics and morphology of galaxies [84].
The current computational possibilities allow following the dynamics of Milky Way-like galaxy for more than 10 Gyr. Recent cosmological simulations undertaken by a few groups for following the formation of galaxies in isolated Milky Way-mass are impressive haloes [85,86,87,88]. These simulations follow many aspects of galaxy formation, including black hole accretion and its feedback, feedback from massive stars, stellar and chemical evolution, metallicity-dependent cooling, star formation, and an influence of magnetic fields, e.g., mock images of the Milky Way-like galaxy in FIRE-2 simulations demonstrate, at z = 0, a developed multi-armed spiral pattern. Our simulations use a more simplified model when gas is treated as a quasi-isothermal one, but we base the simulations of disk dynamics on the current observational knowledge of Milky Way stellar and gaseous rotation curves, stellar, and gaseous velocity dispersions, the density distributions of gaseous and stellar components in the Milky Way disk, as well as on the knowledge of size and density distribution in the stellar bulge of the Milky Way galaxy.

2.6. Dark Matter Halo

The equilibrium of the galactic disks depends on the dark mass distribution in halos [89,90]. The rotation of the Galactic disk is determined by the common gravity of the disk itself, the bulge, and of the massive dark matter halo. Once the rotation curve is specified, as well as the mass and density distributions of the disk and bulge, one can determine the mass of the dark-matter spherical halo. Within the optical radius of the Milky Way disk, the halo mass exceeds that of the disk (see Table 1). The Galactic halo reaches a mass of approximately (5.5–10.3) × 10 10 M within R o p t . In numerical models, with the value of R o p t differing depending on the adopted radial scale length of the disk r d (See Table 1). Based on kinematical data for 200,000 giant stars, [40] estimated the mass of the dark-matter halo within the solar radius to be M h ( r R ) = ( 6 ± 0.9 ) × 10 10 M . However, the galactic halo extends to much larger distances than the optical radius of the stellar disk (for a review see [89,91] and the references there).

3. Numerical Code and Stability Criteria

3.1. Equations and Numerical Algorithm

We treat the dynamics of the 3D stellar-gaseous galactic disk self-consistently. The gaseous subsystem is modeled by N g Smoothed Particle Hydrodynamics (SPH) particles. The dynamics of the collisionless (stellar) disk is modeled using N particles and a direct particle-particle integration scheme. We set the numbers of stellar and gaseous particles equal, N g = N = N / 2 , where N = N g + N is the total number of particles used in the numerical simulations.
Because the star formation is not taken into account in our model, the relative number of the gaseous and stellar particles does not influence the results of simulations. Therefore, based on the convenience of modeling, we chose N = N g .
The equations of motion of gravitationally interacting SPH and collisionless disk particles in the external gravitational fields of an axisymmetric isothermal halo f h and bulge f b are as follows:
d 2 r i d t 2 = p i ϱ i + f i h + f i b + j = 1 , j i N f i j , 0 i < N g ; f i h + f i b + j = 1 , j i N f i j , N g i < N ; ,
where t is the time, ∇ is the Hamilton operator, the radius-vector r i ( t ) specifies the position of the i-th particle in space, while ϱ i and p i are the mass density and gas pressure of the i-th particle, respectively. We use the analytical halo model described in the works [92,93]. The gravitational interaction between the i-th and j-th particles is given by the equation:
f i j = G m j ( r i r j ) | ( r i r j ) 2 + δ 2 | 3 / 2 ,
where m j and r j are the mass and radius-vector of the j-th particle, respectively, and δ is the gravitational softening length, preventing unrealistic accelerations during close encounters of the particles [94,95]. In our simulations, the number of particles is N = ( 2 19 2 21 ) , and the softening length is taken to be δ 40 pc.
With the maximum resolution N = N g = 2 20 in simulations listed in Table 1, masses of the stellar and the gaseous particles are within m 2 3 × 10 4 M , m g 2.4 6.5 × 10 3 M .
In order to describe the dynamics of the gaseous disk, one must supplement the equations of motion with the equation of specific internal energy conservation, e i , and with the equation of state of the gas e = e ( p , ϱ ) :
d e i d t = p i ϱ i · v i ,
e i = p i ( γ 1 ) ϱ i ,
where v i is the velocity vector of the i-th particle and γ is the adiabatic index. We choose an equation of state for the gas close to an isothermal one with γ = 1.05 to decrease the effect of gas heating by shock fronts during the evolution of the gaseous disk. Such an approach is employed as the simplest cooling model in various studies of astrophysical gas dynamics [96,97,98,99]. A more realistic consideration of thermal processes requires a cooling function that depends on temperature, density, and gas metallicity [100,101,102,103]. Acorrect description of gas cooling is possible for multicomponent gas models that take chemical reactions and gas metallicity into account, as is done in the modeling of the molecular clouds in the disks of spiral galaxies [104,105,106]. An accurate accounting of gas cooling should also include radiation transfer (see, e.g., [101]), and we leave this for further consideration.
In accordance with the SPH-approach [107], the density of gas that is associated with the i-th gas particle, the equation of motion (6), and the energy conservation Equation (8) can be written in the following form (details can be found in [108]):
ϱ i = ϱ ( r i ) = j = 1 N g m j W ( | r i r j | , h i j ) .
d v i d t = j = 1 , j i N g m j Π i j W p | r i r j | , h i j + f i h + f i b + j = 1 , j i N f i j ,
d e i d t = 1 2 j = 1 , j i N g m j Π i j ( v i v j ) · W p | r i r j | , h i j .
Here, W is the Monaghan smoothing kernel [107] and W p is the smoothing kernel that is used for the approximation of pressure forces [108,109], and h i j = 0.5 ( h i + h j ) is the effective smoothing length, where the smoothing length for each particle depends on its mass and density as h i = 1.3 m i / ϱ i 1 / 3 [108,110]. Π i j = p i ϱ i 2 + p j ϱ j 2 + ν i j a is the symmetric approximation of the pressure forces, and ν i j a is the artificial viscosity [108]. We have, for the artificial viscosity, ν i j a :
ν i j a = μ i j ( β μ i j α c i j ) ϱ i j , where μ i j = h i j Δ r i j · Δ v i j | Δ r i j | 2 + η h i j 2 , Δ r i j · Δ v i j < 0 0 , Δ r i j · Δ v i j 0 ,
Δ r i j = r i r j , Δ v i j = v i v j , ϱ i j = 0.5 ϱ i + ϱ j , c i j = 0.5 γ p i / ϱ i + γ p j / ϱ j . The empirical constants α , β , and η determine the value of the artificial viscosity. In our simulations ee used the values α = 0.5 , β = 1 , and η = 0.1 .
If the smoothing kernel W (Monaghan cubic spline [107]) is used to calculate the pressure gradient, a nonphysical numerical clustering of particles will occur in highpressure regions [111]. This is caused by the weakening of the interaction between the particles in the vicinity of 0 < ξ < 2 3 and lim ξ 0 W ξ = 0 , where ξ = | r i r j | / h i j is the distance between the i-th and j-th particles. To eliminate this and improve the stability of the numerical algorithm, we use the smoothing kernel W p taken from [109]:
W p ( ξ , h ) = 15 64 π h 3 ( 2 ξ ) 3 , 0 ξ 2 ; 0 , ξ 2 .
From Equation (13), it follows that lim ξ 0 W p ξ = 45 64 π h 4 0 .
For the numerical integration of the differential Equations (11) and (12), we use the predictor–corrector scheme of second-order accuracy (the so-called leapfrog method). We use the direct particle-particle algorithm to calculate the gravitational forces. The leapfrog method allows for us to simulate the dynamics of the disk systems, even in the cartesian coordinate system. Khrapov et al. [112] have shown that use of the leapfrog method in double-precision simulations conserves the total energy, momentum, and angular momentum of the equilibrium system with N = 2 20 particles with an accuracy of 10 5 , 10 15 and 10 13 , respectively. For single-precision simulations, the accuracy of conservation of the above-mentioned quantities is equal to 10 3 , 10 2 and 10 3 , respectively see [112].
Details of the realization of the predictor-corrector method are described by [108,112]. Here, we outline the coordination procedure for simulations of gaseous and stellar disks. For the gaseous disk, the integration time step Δ t g ( t ) is limited by the stability condition of the SPH-algorithm (see e.g., [108]), while, for the stellar disk, the time step is limited by the condition of applicability of Newtonian gravity, namely the time step in the collisionless simulations should be greater than the time of light propagation in the region of the simulations Δ t Δ t c r i t . Here, Δ t c r i t is the propagation time of light within the region. If we choose Δ t = Δ t g , the angular momentum conservation is satisfied with an accuracy of 10 13 . However, in this case, the condition of applicability of Newtonian gravity fails, and the total integration time of the problem increases by factors of tens to hundreds. Therefore, we choose the value Δ t = Δ t c r i t 2 × 10 5 years in our simulations. The calculation of the gravitational interaction between all of the particles was carried out once for a time interval ( t , t + Δ t ) using the expression (7). For the gaseous disk, the number of time steps is large ( n g 1 ) for each time interval ( t , t + Δ t ) , in order that the gravitational force vector is constant during each time interval, which leads to an accuracy of conservation of angular momentum of 10 2 . The following correction procedure for the velocities of gaseous particles at the last time step Δ t g ( t n g ) allows us to increase the accuracy of angular momentum conservation to 10 8 :
v x ( t + Δ t ) = v x ( t n g ) + x Δ v R y Δ v φ R ,
v y ( t + Δ t ) = v y ( t n g ) + x Δ v φ + y Δ v R R ,
v z ( t + Δ t ) = v z ( t n g ) + τ [ F z ( t + Δ t ) F z ( t ) ] ,
where v x , v y , and v z are the components of the velocity vector of gaseous particles in the Cartesian coordinate system ( x , y , z ) , R = x 2 + y 2 , Δ v R = τ [ F R ( t + Δ t ) F R ( t ) ] , Δ v φ = τ [ F φ ( t + Δ t ) F φ ( t ) ] , τ = 0.5 [ Δ t Δ t g ( t n g ) ] . Here, ( F R , F φ , F z ) are the components of the total gravitational force in the cylindrical system of coordinates ( R , φ , z ) .
For multiple GPUs, the details of a parallel OpenMP-CUDA implementation of SPH and N-body numerical algorithms are presented in the following papers [108,112].

3.2. Stability Criteria

Toomre [113] derived a stability criterion that determines the growth of small-scale spiral perturbations. The criterion is:
Q T = c r c T , c T = 3.36 G σ ϰ ,
where σ is the surface density of stars, ϰ is the epicyclic frequency of the disk, and c r is the velocity dispersion of the stellar disk in the direction along the disk radius. For a gaseous disk, the criterion is:
Q T g = c s c T g , c T g = π G σ g ϰ ,
where c s is the sound speed of gas and σ g is the surface density of the gaseous disk. Disks of real galaxies are multi-component, and the formulation of a local stability criterion for such systems is not an easy task. A few different local stability criteria have been suggested. We use the two-component criterion, as suggested by [114]:
Q T = Q T Q T g Q T g + W g Q T , W g = 2 c r c s c r 2 + c s 2 .
We note that gravitating disks are unstable with respect to large-scale spiral perturbations (global modes), even when the local Q-stability parameter is greater than unity, which, in particular, is the case in our numerical simulations.
The nonlinear stage of the large-scale instability that is discussed in Section 4.2 and Section 4.3 re-distributes the surface density and velocity dispersion of the disk.

4. Simulations

It is convenient to choose the following units for the numerical modeling:
m = k m × 10 10 M , r = k r × 10 kpc , v = G m r k m k r × 65.75 km s 1 , t k r 3 / 2 k m 1 / 2 × 142.7 Myr .
With such a definition, the dimensionless gravitational constant and dimensionless mass of the stellar disk for any model will be equal to one for arbitrary values of the parameters k m and k r that were taken from the Table 1. We assume that the optical radius of the disk is equal to four radial scale lengths R o p t = 4 r d = r , where r d is the radial scale length of the thin stellar disk, meaning that the dimensionless value of the optical radius of the disk is equal to the one as well ( R ¯ o p t = 1 ). For example, for models 706 and 707, we have k m = 3.72 , k r = 0.9 which gives m = 3.72 × 10 10 M , r = 9 kpc, v 133.7 km s 1 and t 63.2 Myr.
Henceforth, we will use the dimensionless units unless the dimension of the quantity is explicitly indicated.
  • Rotation curve of the cold gaseous component V g ( r ) .
  • Rotation curve of the stellar disk V ( r ) .
  • Velocity dispersion of the stars in the solar neighborhood c r together with its radial dependence c r ( r ) .
  • Surface density of the thin stellar disk in the solar neighborhood σ .
  • Radial scale length of the density distribution of the thin stellar disk r d .
  • Surface density of the gaseous disk in the solar neighborhood σ g together with the surface density dependence on radius σ g ( r ) .
  • Vertical scale height of the thin stellar disk h.
  • Size of the stellar bar r b a r .
  • Radial scale length of the bulge density distribution b and its mass M b .
We have constructed over one hundred numerical models of the Milky Way galaxy using different sets of equilibrium parameters. Table 1 lists the models satisfying the observational properties of the Milky Way. To discuss the results of the numerical simulations, we choose model 706 (see Table 1), which agrees with the available observational data for the Milky Way galaxy.
When building equilibrium models of the Milky Way galaxy, it is important that the model parameters are consistent with observational data. Below are listed the specific parameters and equilibrium profiles that must be matched with the available observational data for the Milky Way.
Equilibrium particle distributions were constructed using the force balance equations. For the collisionless (stellar) particles, the Jeans equations were solved along the radial and vertical coordinates. For the gas particles, the hydrostatic-balance equations using a polytropic equation of state were solved in the z and r -directions. Because of the assumptions used to constrain the equilibrium state, and due to the finite number of particles modeling both the collisionless and the gaseous disks, the equilibrium solution is approximate and it contains so-called numerical noise. Thus, the evolution of a stellar-gas system begins from a quasistationary equilibrium that includes practically a full spectrum of initial perturbations. This instability leads to the growth of the most unstable modes that form the observed spiral structure.

4.1. Spiral Structure

For the models considered, Figure 3 shows the typical radial dependence of the Toomre Q-parameter ( Q T ( r ) ). There is a moderate increase of the Q-parameter in the region r d < r < 4 r d with a broad minimum at the value Q T = (1–1.5). As is known from previous studies, in this case we have conditions for a spiral global instability and the formation of spiral structure in the disk, including a central bar.
Figure 4 shows the typical temporal evolution of a stellar gaseous disk with equilibrium parameters in agreement with observational data. The formation of an open spiral/bar in the central regions of the disk occurs rather quickly after approximately 1.5 disk revolutions and, after approximately two disk revolutions, the bar is formed. Outside the stellar bar, a complex spiral pattern grows. Even in the stellar disk, with its more regular spirals, a superposition of two- and three- armed spiral patterns with different amplitudes grows and reaches the saturation stage. The gaseous disk exhibits an even more complicated spiral structure due to nonlinear interactions of the unstable modes, which leads to the branching of spirals, the appearance of rows, and that of a multi-level spiral pattern at the periphery of the disk.
Fourier analysis of the surface density perturbations provides a means for visualizing the growing spiral and bar-like structures. It is convenient to characterize the growth of perturbations by the time dependence of the global Fourier amplitudes for different azimuthal numbers m [94,115]:
A m ( t ) = 1 N k j = 1 N k exp I m φ j ( t ) ,
where I = 1 , φ j is the azimuthal angle (radians) of the j-particle, and the summation is taken for the particles inside a radial layer r k r j r k + 1 , N k is the total number of particles in this layer.
Figure 5 and Figure 6 show the time dependence of the amplitudes of the m = (2–4) harmonics for different regions of the stellar and gaseous disks. In the central regions of the stellar disk r < 2 r d (see the left and middle panels in Figure 5), the dominant mode is m = 2 harmonics, related to the bar. In addition to the m = 2 mode, the disk is susceptible to the instability of m = 3 and m = 4 spirals, which are more noticeable at the periphery of the disk. The growth rates of perturbations in the central regions of the disk are noticeably larger when compared to that in the periphery of the disk. The temporal behavior of perturbations in the gaseous disk is qualitatively similar to that of the stellar disk. As one can see from Figure 6, in the central regions of the gaseous disk, the mode m = 2 also prevails, while, outside the disk’s central regions, the perturbations are represented by a combination of m = 2 , m = 3 , and m = 4 spiral modes.
The number of spiral arms in the Milky Way galaxy, as well as their pitch angle, has been the subject of studies for a long time without consensus among the researchers as to the results. From observations of the stellar and the gaseous disks, estimates of the pitch angle of the Milky Way spiral pattern range from 5 to 25 [116,117,118]. Unsuccessful attempts to understand the properties of the global spiral pattern in the Milky Way disk have led to an approach where segments of the spiral arms are treated separately [26,118,119]. The disagreement between various determinations of the properties of the Milky Way spiral arms that are based on observational data (e.g., masers, ionized. neutral and molecular hydrogen gas, and 2MASS sources) that we suggest is a manifestation of the complex and non-stationary spiral structure of the Milky Way disk. We believe that this structure is a superposition of nonlinear spiral patterns with different azimuthal wavenumbers, angular speeds of the patterns, and amplitudes.
Here we have paid special attention to the reconstruction of the maximum of the Milky Way disk’s rotational velocity in its central region r 1 kpc. Figure 7 illustrates the typical spatial distributions of the gas and of stars in models 401, 405, 455, and 474, where the rotational velocity in the disk has an internal maximum in the region r 1 kpc due to the bulge-disk mass distribution. We find that, due to a massive and centrally concentrated bulge, which is necessary for reconstructing the internal peak on the velocity rotation curve, the central bar either does not form or has a short major axis and is short-lived. Even if the bar appears at the initial stages of the disk evolution in a relatively cold stellar disk, it is destroyed at later stages. A strong gradient in the gravitational potential in the central regions of the disk has a well-known [120,121] destructive role on bar formation that is caused by the presence of a compact massive bulge. An additional factor suppressing bar instability is the increase of gas density in the central regions of the disk in the process of disk evolution. These arguments lead us to choose models in which the central peak of the disk’s rotation curve is not explained by the presence in the disk of a massive and compact bulge. Instead, the central peak of the disk’s rotation curve is formed by the process of nonlinear evolution of spiral/bar perturbations. Examples of these models are 600, 602, 610, 701 shown in Figure 7 and 706 in Figure 4. These models develop a pronounced stellar bar with gas kinematics that are similar to that observed for the Milky Way bar.
Outside the central regions of the stellar disk, models with a bar demonstrate a mixture of two-, three-, and four-armed spiral patterns of lower amplitude. The gaseous disks show a more complicated behavior, with the appearance of a non-stationary spiral pattern.
We note that models with a relatively low stellar velocity dispersion (e.g., model 701 in Table 1) develop unrealistically high-amplitude spiral perturbations. Actually, the development of a strong spiral perturbations leads to a heating of the stellar disk, and the velocity dispersion of the stars increases to approximately 40 km s 1 in a short period of time. Thus, the observed amplitudes of spiral perturbations in the disk of the Milky Way impose additional restrictions on the velocity dispersion of stars in the solar neighborhood. Models with a velocity dispersion c r = (29–33) km s 1 support a long-lived stellar bar with a semi-axis of ∼3 kpc, together with a spiral pattern that has observable amplitudes. Edge-on images of the numerical models demonstrate features that are similar to those observed in the central regions of the Milky Way, such as the so-called X-structure (see the edge-on images of models 401, 455, 602, 610, and 701 in Figure 7.
The radial scale length of the stellar disk is of critical importance in the development of a bar in the disk’s central regions. The models with a radial scale length of the stellar disk larger than r d = 2.6 kpc, such as model 730 with r d = 2.6 kpc or model 720 with r d = 3 kpc, lead to a relatively low surface density disk in the central region σ 0 700 M pc 2 . That, in turn, requires the presence of a massive dark matter halo to explain the observed disk rotation curve. Such models do not allow the formation of a bar-mode in the disk (see, e.g., model 730 in Figure 7). Besides that, for model 730, the values of the parameters Q T and Q T g are (15–20)% larger within the optical radius when compared to models 706 and 707, which also leads to the stabilisation of the bar-mode. Attempts to explain the observed central maximum of the rotational velocity of the disk in these modes (curves 1, 2 in Figure 1) were unsuccessful. We then considered a set of models with a less concentrated bulge, which explained the appearance of the maximum of the disk’s rotational velocity as a result of bar appearance and the nonlinear evolution of perturbations.

4.2. Density Profile and Rotation Curve

Here, we discuss the behavior of model 707 (see Table 1), which satisfies all of the limitations imposed by the observational data for the Milky Way galaxy, as discussed above. Figure 8 shows the decomposition of the rotation curve for this model.
As one can see from Figure 8 and Figure 9, the nonlinear stage of instability leads to a steepening of the gas’s rotation curve in the disk’s central regions. However, note that the amplitude of the central peak in the gaseous rotation curve is smaller than that observed for the velocity peak in the rotation curve of the Milky Way gaseous disk.
Figure 10 shows the azimuthal density profiles of the stellar σ ( φ ) and of the gaseous disks σ g ( φ ) at three fixed radii r = 0.15 (inside the bar), r = 0.51 (at two exponential scale lengths of the disk), and at the periphery of the stellar disk close to the Sun’s orbit r = 0.89 . Dashed black lines correspond to the averaged values of the density over the azimuthal angle φ . Azimuthal density distributions demonstrate the features of the spiral disturbances. In the central zone, the stellar density distribution shows a pronounced two-arm mode. As for the gaseous component, in addition to the main mode, it shows the presence of two other modes at lower amplitudes. At two radial scale lengths of the disk ( r 2 r d ), the stellar system demonstrates a four-armed spiral structure. At the periphery of the stellar disk, the amplitudes of the density perturbations decrease considerably and remain within 10–20 percent. Gas shows more complicated behavior, which can be described as a superposition of two-, three-, and four-armed spiral patterns (see Figure 6).
The radial scale length of the stellar density distribution is one of the key parameters that determines the mass of the stellar disk and its stability properties. The estimated radial scale length of the Milky Way’s stellar density distribution is about 2.25 kpc [36] and we use that value in our simulations. In the set of models with r d = 3 kpc and r d = 2.6 kpc, the surface density of the stellar disk decreases in the central regions to σ 0 = 481 M / pc 2 and σ 0 = 725 M / pc 2 (corresponding to models 720 and 730). In order to satisfy the observed rotation curves in these models, one has to increase the mass of the dark matter halo relative to the mass of the disk to μ 5 . This, as a result, suppresses the formation of the central bar in the Milky Way galaxy. In other words, with the observed surface density of the solar neighborhood, the dynamical models do not allow us to reproduce the stellar bar in the central regions of the disk if the exponential scale length of the surface density distribution is more than 2.5 kpc.
We find that the dynamical models that conform not only to the observed rotation curve of the galactic disk, but also to the surface densities of the stellar and gaseous components in the solar neighborhood, together with their observed velocity dispersions, give a relatively small total mass of the disk M d = M + M g 4.5 × 10 10 M within 2 R o p t = 18 kpc. Therefore, dark matter dominates within 2 R o p t = 18 kpc, composing about 74 percent of the total mass of the Milky Way galaxy.

4.3. Stellar Velocity Dispersion

As an initial condition, we choose a radial dependence of the velocity dispersion that corresponds to the observational data for the Milky Way (curve 6 on Figure 2). At the galactocentric radius of the Sun, this choice gives the value c r (30–32) km s 1 . The tangential component of the velocity dispersion in the collisionless disk c φ is prescribed according to Equation (5).
During the evolution of the perturbations, the velocity dispersion of the stars evolves. Figure 11 illustrates the temporal evolution of the velocity dispersion for a number of models. As one can see, the typical behavior of the velocity dispersion is its slow evolution during approximately one billion years. At roughly this time, the formation of spiral structure and of a bar takes place, and it is accompanied by essentially no heating of the collisionless disk. At later times, some increase of the velocity dispersion of the disk is observed.
During the development of perturbations, the velocity dispersions deviate from those of the axisymmetric case. Figure 12 shows the dependence of the velocity dispersion c r as a function of the azimuthal angle φ at three different radii of the disk for model 707 at time t = 9 . As one sees, large amplitudes of the density perturbations in the central regions of the disk are accompanied by strong variations in the velocity dispersion c r . The variations of the velocity dispersion are less prominent at the periphery of the disk.

4.4. Disk Kinematics in Central Regions

Because of presence of the bar, the kinematical properties of the stellar and gaseous disks inside r 3 kpc demonstrate some specific features. Figure 13 shows the rotation curves of the disk projected along lines representing different viewing orientations with respect to the semimajor axis of the bar. As one can see, the rotation curve of the disk along the semi-minor-axis of the bar has a local maximum inside 1 kpc from the center of galaxy, similar to the observed kinematical properties of the Milky Way disk within the central kiloparsec (see Figure 1).
The estimates of the angular velocity and mass of the Milky Way bar [81] provide values of Ω b a r < 50 km   s 1 1 kpc 1 and M b a r 2 × 10 9 M , respectively.
The conditions for the formation of a bar are rather sensitive to the mass of the dark matter halo μ and that of the bulge μ b , relative to the mass of the thin disk. An important parameter is the radial scale length of the disk density distribution together with the scale length of the mass distribution of the bulge, as mentioned above. In the models with μ 2 , the bar mode does not develop due to the suppression of bar instability by the massive halo. The classical stability criterion of the global bar mode [122] requires that the ratio of the kinetic energy of disk rotation E r o t to the total gravitational energy | E Φ | does not exceed some critical value ( E r o t / | E Φ | ) c r i t 0.14 . Thus, galaxies with a dark halo are more stable with respect to the formation of the bar mode. The authors [122] made a fundamental conclusion, that stabilization of the bar mode requires the value of the parameter μ within (1–2). Criteria of this kind have been tested using better N-body models and they show that the exact value of μ depends on a number of parameters of the galactic subsystems that determine the properties of the bulge, the gaseous disk, the radial profile of the dark matter, etc. [22,27,81,89,123].
Let us estimate the influence of the relative mass of the halo on the velocity dispersion of stars that are close to the boundary of stability of the system. The rotational velocity of the the disk is V R Φ R where R is some radius, and Φ is the gravitational potential. Distracting from the particular mass distribution, we can write: Φ / R G ( M d + M h ) / R 2 . On the boundary of stability, the velocity dispersion of stars is inversely proportional to the epicyclic frequency: c r 1 / ϰ 1 / M d + M h . As a result c r / V 1 / ( M d + M h ) 1 / ( 1 + μ ) for the fixed mass of the disk M d , so the disk at the boundary of stability is colder in the models with a more massive halo.
Similarly, in the models with a moderate ratio of masses of the halo and of the disk μ 1.6 , but with a massive and compact bulge b / R o p t 0.014 , the bar is not reproduced in the numerical simulations due to scattering of the orbits of the collisionless particles in the central regions where the gravitational potential has a large gradient. Thus, the presence of a bar in the Milky Way disk puts rather strict limitations on its equilibrium properties.

5. Discussion

Athanassola [124] showed that a bar grows faster in galactic models with a live halo as compared to models with a rigid one. Athanassoula attributed this effect to the destabilizing influence of the resonant halo particles. The total number of particles in Athanassoula’s numerical simulations was of order of 10 6 with the mass of individual particle being about 2.5 × 10 5 M . This leads, apart from the resonant effects, to the generation of high amplitude noise in the disk component. Polyachenko et al. [123] returned to the problem of bar formation using one to two orders of magnitudes larger numbers of particles. They found that, while a live halo decreases the time of growth for the bar instability from 500 to 250 Myr, the pattern speed and other parameters of the bar remain approximately the same. Because we are interested in the late stages of disk evolution, the possible influence of a live halo on the growth rate of the bar instability is insignificant for our purpose.
A possible influence of the “live” halo on the dynamics of the disks can be obscured in the numerical simulations by the particle resolution when massive dark matter “particles” generate a high amplitude noise in the disk.
Long et al. [125] and Collier et al. [126] showed that spinning dark live halos suppresses the growth of stellar bars. Recent studies based on GAIA DR2 observational data [127] demonstrate that the Milky Way halo has close to zero rotational velocity with average value of 1 ± 4 km   s 1 .
Chemin discussed the dependence of the observed rotation curve on bar orientation [128]. He used the results of high-resolution hydrodynamical simulations by Renaud et al. [129] and demonstrated that the rotation curve of the disk depends on bar orientation relative to an observer. We confirm this result in a more realistic model using combined, direct N-body-gas dynamical simulations.
Deg et al. [130] presented a new version of the GlactICS code, and tested it simulating the dynamics of a Milky Way-like galaxy. The authors used a two-component collisionlessgaseous model, and included bulge and dark matter halo in their simulations. Similar to our results, they find that, after approximately 1 Gyr, a bar formation with semi-axis of about 3 kpc occurs in the disk, and a multi-armed spiral pattern, which is more prominent in the gaseous component, is formed outside the bar region. Strictly speaking, to model the dynamics of a particular galaxy, one must guess its initial equilibrium, one that fits the final stage of observed equilibrium properties for the galactic disk. However, the disk evolution does not noticeably alter the initial density distribution of the massive stellar disk, and up to 0.5 Gyr does not affect considerably the density distribution of the gaseous disk, thus justifying our approach, as seen in Figure 9. Therefore, the set of unstable spiral modes is well-determined by the current state of a slowly evolving disk. The succcessful modeling of spiral structure in nearby spiral galaxies NGC 1566 [131] and NGC 5247 [93], which was done using the observed equilibrium properties of these systems, further justifies the validity of our approach.

Swing Amplification vs. Global Modes

In the review paper by Dobbs and Baba [7] the authors state “We argue that, with the possible exception of barred galaxies, spiral arms are transient, recurrent, and initiated by swing amplified instabilities in the disc”. The swing amplification mechanism of spiral formation was suggested by A. Toomre (see, e.g., [132]) and it remains in many studies the governing mechanism for explaining the origin of spiral structure. However, there are a number of unanswered questions related to this approach. The observed amplitudes of the spirals is one of the oldest of issues with the swing amplification theory. Swing amplification can amplify the initial perturbations by about a hundred times [132], which is not sufficient for explaining the observed amplitudes of the spiral arms in galaxies.
Similar to other hydrodynamical or plasma systems, galactic disks oscillate, and they have their own global modes that can be unstable. Why are these global modes not considered in the explanation of spiral structure in galaxies? What happens to these global modes in an explanation involving the swing amplification mechanism? Why are global modes observed in barred galaxies described as the possible exceptionato the recurrent and swing amplified instabilities? Differential rotation and gravity also work in barred galaxies, so the swing amplified noise should be seen in these galaxies too.
It is assumed that the origin of the spiral pattern recently discovered in the galaxy M51a is caused by its satellite galaxy, NGC 5195. However, the closest approach between these two systems was about 900 Myr ago, according to Wahde and Donner [133]. Is the observed two-armed global spiral structure attributed to an interaction that occurred about one Gyr ago? Why would swing amplification, operating on much shorted time scales, be suppressed in this galaxy?
Sellwood & Carlberg concluded in their recent paper [134] “We argue that apparently shearing transient spirals in simulations result from the superposition of two or more steadily rotating patterns, each of which is best accounted for as a normal mode of the non-smooth disc.” We completely agree with this statement, and will also cite here Vadim Antonov, the author of classical works in stellar dynamics (Antonov’s theorems in Binney and Tremaine’s Galactic Dynamics), who said: “The galaxies are similar to copper trampets”. We can add here that some galaxies, like NGC 1566, ”play a pure spiral tone”, while others “play an accord” of spiral patterns.
In light of this, a number of questions arise:
What is the mechanism of spiral saturation at the nonlinear stage of instability? Laughlin et al. [135] concluded that nonlinear self-interaction of a growing mode is responsible for mode saturation, but further study is needed.
Why does the presence of gas in the stellar-gaseous gravitating disk extend the lifetime of the spiral pattern? Goldreich & Lynden-Bell [82] pointed out, back in the sixties, the importance of gas to the formation of spiral structure: S0 galaxies are topographically similar to normal spirals but they have no gas, no dust, and no spiral arms. The importance of gas in sustaining the spiral structure was also confirmed in numerical simulations [93]. In light of new observational data, the spiral structure in galaxies continues to pose questions to theory.

6. Summary

To model the dynamics of a Milky Way-like collisionless-gaseous disk, we use, for the first time in such a study, direct integration in calculating the gravitational potential, contrary to the approximate tree-code or particle-mesh codes that were used before. It is not clear whether an approximate treatment of the gravity can correctly reproduce the dynamics of the multi-component disks over cosmological time intervals. An agreement of simulation results in the hydrodynamical, collisionless, and the collisionless–hydrodynamical approaches [93] provides validity to our model. As for particular results, we demonstrate that:
An axisymmetric two-component gravitating disk with parameters close to those of the azimuthally averaged Milky Way galaxy—in terms of observed rotation curve, velocity dispersion profile, and masses of the bulge and of the stellar and gaseous disk components—is unstable towards near-exponentially growing spirals having numbers of arms m = (2–4).
At the nonlinear stage of instability, the spirals saturate at few tens of percent in the central regions of the disk, to a few percent at the disk’s periphery.
At the nonlinear stage, a prominent bar is formed in the central regions of the disk with a large semi-axis of about 3 kpc.
Outside the bar region, a complex spiral structure, represented by a superposition of two-, three-, and four-armed spiral patterns, rotating with different angular velocities, is formed. The spiral structure of the Milky Way galaxy is interpreted in some papers as single spiral pattern with a fixed number of arms and spiral pitch angle, and one rotation resonance with fixed position. We show that this is not the case for the Milky Way-like disk.
We demonstrate that the peak in the rotation curve of the disk of the Milky Way galaxy, which is located in its central regions, is a result of the non-circular motions caused by the bar which develops in the disk.
We confirm that the presence of a massive and centrally concentrated bulge prevents the formation of a bar, in agreement with [122].
We also show that the presence of gas of about ten percent the disk’s mass extends the lifetime of the spiral structure to a few Gyr as compared to what is found in purely collisionless models.
This conclusion though should be validated in further studies with a multi-phase gas model using proper feedback, which is beyond the modelling capabiities of this study.

Author Contributions

Conceptualization, methodology: V.K. and A.K.; Writing—original draft preparation, V.K. and A.K.; Writing—review & editing, all authors; software and validation, S.K.; numerical simulations, S.K. and A.K.; visualization, S.K. and A.K.; supervision and project administration, V.K. All authors have read and agreed to the published version of the manuscript.

Funding

The development of software for modeling galaxies and numerical experiments were funded by the Ministry of Science and Higher Education of the Russian Federation (the government task no. 0633-2020-0003, S.K. and A.K.).

Institutional Review Board Statement

The approval was waived because research is fundamental and is not included into research plans of the universities the authors are affiliated with.

Data Availability Statement

All data used in this paper are taken from the open sources and the references are given.

Acknowledgments

Authors thank T. Girard and the referees for valuable comments. Numerical simulations were carried by using the equipment of the shared research facilities of “Supercomputer Center of Volgograd State University”. V.K. acknowledges financial support by Southern Federal University, 2020 (Ministry of Science and Higher Education of Russian Federation).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
COBE/DIRBECosmic Background Explorer/Diffuse Infrared Background Experiment
SPHSmoothed-particle hydrodynamics
OpenMP-CUDAOpen Multi-Processing - Compute Unified Device Architecture
VLBIVery Long Baseline Interferometry

References

  1. Bertin, G.; Lin, C.C. Spiral Structure in Galaxies a Density Wave Theory; MA MIT Press: Cambridge, MA, USA, 1996; 271p. [Google Scholar]
  2. Shu, F.H. Six Decades of Spiral Density Wave Theory. Annu. Rev. Astron. Astrophys. 2016, 54, 667–724. [Google Scholar] [CrossRef]
  3. Bertin, G.; Lin, C.C.; Lowe, S.A.; Thurstans, R.P. Modal approach to the morphology of spiral galaxies. I—Basic structure and astrophysical viability. Astrophys. J. 1989, 338, 78–103. [Google Scholar] [CrossRef]
  4. Gerola, H.; Seiden, P.E. Stochastic star formation and spiral structure of galaxies. Astrophys. J. 1978, 223, 129–139. [Google Scholar] [CrossRef]
  5. Sellwood, J.A.; Trick, W.H.; Carlberg, R.G.; Coronado, J.; Rix, H.-W. Discriminating among theories of spiral structure using Gaia DR2. Mon. Not. R. Astron. Soc. 2019, 484, 3154–3167. [Google Scholar] [CrossRef] [Green Version]
  6. Bertin, G.; Lin, C.C.; Lowe, S.A.; Thurstans, R.P. Modal Approach to the Morphology of Spiral Galaxies. II. Dynamical Mechanisms. Astrophys. J. 1989, 338, 104–120. [Google Scholar] [CrossRef]
  7. Dobbs, C.; Baba, J. Dawes Review 4: Spiral Structures in Disc Galaxies. Publ. Astron. Soc. Aust. 2014, 31, e035. [Google Scholar] [CrossRef] [Green Version]
  8. D’Onghia, E. Disk-stability Constraints on the Number of Arms in Spiral Galaxies. Astrophys. J. Lett. 2015, 808, L8. [Google Scholar] [CrossRef]
  9. Fujii, M.S.; Bédorf, J.; Baba, J.; Portegies Zwart, S. The dynamics of stellar discs in live dark-matter haloes. Mon. Not. R. Astron. Soc. 2018, 477, 1451–1471. [Google Scholar] [CrossRef] [Green Version]
  10. Michikoshi, S.; Kokubo, E. Galactic Spiral Arms by Swing Amplification. Astrophys. J. 2016, 821, 35. [Google Scholar] [CrossRef] [Green Version]
  11. Sellwood, J.A.; Carlberg, R.G. Spiral instabilities: Mechanism for recurrence. Mon. Not. R. Astron. Soc. 2019, 489, 116–131. [Google Scholar] [CrossRef] [Green Version]
  12. Pettitt, A.R.; Tasker, E.J.; Wadsley, J.W.; Keller, B.W.; Benincasa, S.M. Star formation and ISM morphology in tidally induced spiral structures. Mon. Not. R. Astron. Soc. 2017, 468, 4189–4204. [Google Scholar] [CrossRef] [Green Version]
  13. Khoperskov, A.V.; Khoperskov, S.A.; Zasov, A.V.; Bizyaev, D.V.; Khrapov, S.S. Interaction between collisionless galactic discs and non-axisymmetric dark matter haloes. Mon. Not. R. Astron. Soc. 2013, 431, 1230–1239. [Google Scholar] [CrossRef] [Green Version]
  14. Shabani, F.; Grebel, E.K.; Pasquali, A.; D’Onghia, E.D.; Gallagher, J.S., III; Adamo, A.; Messa, M.; Elmegreen, B.G.; Dobbs, C.; Gouliermis, D.A.; et al. Search for star cluster age gradients across spiral arms of three LEGUS disc galaxies. Mon. Not. R. Astron. Soc. 2018, 478, 3590–3604. [Google Scholar] [CrossRef] [Green Version]
  15. Peterken, T.G.; Merrifield, M.R.; Aragon-Salamanca, A.; Coleman, N.D.; Krawczyk, M.; Masters, K.L.; Wijmans, A.-M.; Westfall, K.B. A direct test of density wave theory in a grand-design spiral galaxy. Nat. Astron. 2019, 3, 178–182. [Google Scholar] [CrossRef]
  16. Bialopetravicius, J.; Narbutis, D. Study of star clusters in the M83 galaxy with a convolutional neural network. Astron. J. 2020, 160, 264. [Google Scholar] [CrossRef]
  17. Abdeen, S.; Kennefick, D.; Kennefick, J.; Eufrasio, R.; Davies, B.L.; Miller, R.; Shields, D.W.; Monson, E.B.; Basset, C.; O’Mara, H. Evidence in favour of density wave theory through star formation history maps and spatially resolved stellar clusters. arXiv 2020, arXiv:2010.14540. [Google Scholar]
  18. Hou, L.G.; Han, J.L.; Shi, W.B. The spiral structure of our Milky Way Galaxy. Astron. Astrophys. 2009, 499, 473–482. [Google Scholar] [CrossRef] [Green Version]
  19. Hou, L.G.; Han, J.L. The observed spiral structure of the Milky Way. Astron. Astrophys. 2014, 569, A125. [Google Scholar] [CrossRef] [Green Version]
  20. Xu, Y.; Reid, M.; Dame, T.; Menten, K.; Sakai, N.; Li, J.; Brunthaler, A.; Moscadelli, L.; Zhang, B.; Zheng, X. The local spiral structure of the Milky Way. Sci. Adv. 2016, 2, e1600878. [Google Scholar] [CrossRef] [Green Version]
  21. Bland-Hawthorn, J.; Gerhard, O. The Galaxy in Context: Structural, Kinematic, and Integrated Properties. Annu. Rev. Astron. Astrophys. 2016, 54, 529–596. [Google Scholar] [CrossRef] [Green Version]
  22. Shen, J.; Zheng, X.-W. The bar and spiral arms in the Milky Way: Structure and kinematics. Res. Astron. Astrophys. 2020, 20, 159. [Google Scholar] [CrossRef]
  23. Sanna, A.; Reid, M.J.; Dame, T.M.; Menten, K.M.; Brunthaler, A. Mapping spiral structure on the far side of the Milky Way. Science 2017, 358, 227–230. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Zhang, B.; Reid, M.J.; Zhang, L.; Wu, Y.; Hu, B.; Sakai, N.; Menten, K.M.; Zheng, X.; Brunthaler, A.; Dame, T.M.; et al. Parallaxes for Star-forming Regions in the Inner Perseus Spiral Arm. Astron. J. 2019, 157, 200. [Google Scholar] [CrossRef] [Green Version]
  25. Katz, D.; Antoja, T.; Romero-Gómez, M.; Drimmel, R.; Reylé, C.; Seabroke, G.M.; Soubiran, C.; Babusiaux, C.; Di Matteo, P.; et al. Gaia Data Release 2. Mapping the Milky Way disc kinematics. Astron. Astrophys. 2018, 616, A11. [Google Scholar] [CrossRef] [Green Version]
  26. Xu, Y.; Bian, S.B.; Reid, M.J.; Li, J.J.; Zhang, B.; Yan, Q.Z.; Dame, T.M.; Menten, K.M.; He, Z.H.; Liao, S.L.; et al. A comparison of the local spiral structure from Gaia DR2 and VLBI maser parallaxes. Astron. Astrophys. 2018, 616, L15. [Google Scholar] [CrossRef] [Green Version]
  27. Khoperskov, S.; Gerhard, O.; Di Matteo, P.; Haywood, M.; Katz, D.; Khrapov, S.; Khoperskov, A.; Arnaboldi, M. Hic sunt dracones: Cartography of the Milky Way spiral arms and bar resonances with Gaia Data Release 2. Astron. Astrophys. 2020, 4, L8. [Google Scholar] [CrossRef] [Green Version]
  28. Griv, E.; Gedalin, M.; Shih, I.-C.; Hou, L.-G.; Jiang, I.-G. The nearby spiral density-wave structure of the Galaxy: Line-of-sight velocities of the Gaia DR2 main-sequence A, F, G, and K stars. Mon. Not. R. Astron. Soc. 2020, 493, 2111–2126. [Google Scholar] [CrossRef]
  29. Brown, A.G.A.; Vallenari, A.; Prusti, T.; de Bruijne, J.H.J.; Babusiaux, C.; Biermann, M. Gaia Early Data Release 3: Summary of the contents and survey properties. arXiv 2020, arXiv:2012.01533. [Google Scholar]
  30. Rodriguez-Fernandez, N.J.; Combes, F. Gas flow models in the Milky Way embedded bars. Astron. Astrophys. 2008, 489, 115–133. [Google Scholar] [CrossRef] [Green Version]
  31. Baba, J.; Saitoh, T.R.; Wada, K. On the Interpretation of the l-v Features in the Milky Way Galaxy. Publ. Astron. Soc. Jpn. 2010, 62, 1413–1422. [Google Scholar] [CrossRef] [Green Version]
  32. Sormani, M.C.; Binney, J.; Magorrian, J. Gas flow in barred potentials—III. Effects of varying the quadrupole. Mon. Not. R. Astron. Soc. 2015, 454, 1818–1839. [Google Scholar] [CrossRef] [Green Version]
  33. Shen, J.; Rich, R.M.; Kormendy, J.; Howard, C.D.; De Propris, R.; Kunder, A. Our Milky Way as a Pure-disk Galaxy—A Challenge for Galaxy Formation. Astrophys. J. Lett. 2010, 720, L72–L76. [Google Scholar] [CrossRef] [Green Version]
  34. Fujii, M.S.; Bédorf, J.; Baba, J.; Portegies Zwart, S. Modelling the Milky Way as a dry Galaxy. Mon. Not. R. Astron. Soc. 2019, 482, 1983–2015. [Google Scholar] [CrossRef]
  35. Korchagin, V.I.; Khoperskov, S.A.; Khoperskov, A.V. Role of gaseous disk in the formation of the spiral structure of the Milky Way galaxy. Balt. Astron. 2016, 25, 356–361. [Google Scholar] [CrossRef] [Green Version]
  36. Drimmel, R.; Spergel, D.N. Three-dimensional Structure of the Milky Way Disk: The Distribution of Stars and Dust beyond 0.35 Rsolar. Astrophys. J. 2001, 556, 181–202. [Google Scholar] [CrossRef] [Green Version]
  37. Bovy, J.; Rix, H.-W.; Schlafly, E.F.; Nidever, D.L.; Holtzman, J.A.; Shetrone, M.; Beers, T.C. The Stellar Population Structure of the Galactic Disk. Astrophys. J. 2016, 823, 30. [Google Scholar] [CrossRef] [Green Version]
  38. Liu, C.; Xu, Y.; Wang, H.; Wan, J. Rediscovering the Galactic outer disk with LAMOST data. Proc. Int. Astron. Union 2017, 13, 109–115. [Google Scholar] [CrossRef] [Green Version]
  39. Gerhard, O.E. Structure and Mass Distribution of the Milky Way Bulge and Disk. Astron. Soc. Pac. Conf. Ser. 2001, 230, 21–30. [Google Scholar]
  40. Piffl, T.; Binney, J.; McMillan, P.J.; Steinmetz, M.; Helmi, A.; Wyse, R.F.G.; Bienaym, O.; Bland-Hawthorn, J.; Freeman, K.; Gibson, B.; et al. Constraining the Galaxy’s dark halo with RAVE stars. Mon. Not. R. Astron. Soc. 2014, 445, 3133–3151. [Google Scholar] [CrossRef] [Green Version]
  41. Khoperskov, A.V.; Tyurina, N.V. A Dynamical Model of the Galaxy. Astron. Rep. 2003, 47, 443–457. [Google Scholar] [CrossRef]
  42. McKee, C.F.; Parravano, A.; Hollenbach, D.J. Stars, Gas, and Dark Matter in the Solar Neighborhood. Astrophys. J. 2015, 814, 13. [Google Scholar] [CrossRef] [Green Version]
  43. Mosenkov, A.V.; Sotnikova, N.Y.; Reshetnikov, V.P.; Bizyaev, D.V.; Kautsch, S.J. Does the stellar disc flattening depend on the galaxy type? Mon. Not. R. Astron. Soc. 2015, 451, 2376–2389. [Google Scholar] [CrossRef] [Green Version]
  44. Comerón, S.; Salo, H.; Knapen, J.H. The reports of thick discs’ deaths are greatly exaggerated. Thick discs are NOT artefacts caused by diffuse scattered light. Astron. Astrophys. 2018, 610, A5. [Google Scholar] [CrossRef] [Green Version]
  45. Sarkar, S.; Jog, C.J. Vertical stellar density distribution in a non-isothermal galactic disc. Mon. Not. R. Astron. Soc. 2020, 499, 2523–2533. [Google Scholar] [CrossRef]
  46. Bahcall, J.N. The distribution of stars perpendicular to galactic disk. Astrophys. J. 1984, 276, 156–168. [Google Scholar] [CrossRef]
  47. Khoperskov, A.; Bizyaev, D.; Tiurina, N.; Butenko, M. Numerical modelling of the vertical structure and dark halo parameters in disc galaxies. Astron. Nachrichten 2010, 331, 731–745. [Google Scholar] [CrossRef] [Green Version]
  48. Juric, M.; Ivezic, Z.; Brooks, A.; Lupton, R.H.; Schlegel, D. The Milky Way Tomography with SDSS. I. Stellar Number Density Distribution. Astrophys. J. 2008, 673, 864–914. [Google Scholar] [CrossRef]
  49. Wan, J.-C.; Liu, C.; Deng, L.-C. Red Clump Stars from LAMOST II: The outer disc of the Milky Way. Res. Astron. Astrophys. 2017, 17, 79. [Google Scholar] [CrossRef] [Green Version]
  50. Eilers, A.-C.; Hogg, D.W.; Rix, H.-W.; Ness, M.K. The Circular Velocity Curve of the Milky Way from 5 to 25 kpc. Astrophys. J. 2019, 871, 120. [Google Scholar] [CrossRef] [Green Version]
  51. Sofue, Y. Rotation Curve of the Milky Way and the Dark Matter Density. Galaxies 2020, 8, 37. [Google Scholar] [CrossRef]
  52. Chrobakova, Z.; Lopez-Corredoira, M.; Sylos Labini, F.; Wang, H.-F.; Nagy, R. Gaia-DR2 extended kinematical maps. III. Rotation curves analysis, dark matter, and MOND tests. Astron. Astrophys. 2020, 642, A95. [Google Scholar] [CrossRef]
  53. Rastorguev, A.S.; Utkin, N.D.; Zabolotskikh, M.V.; Dambis, A.K.; Bajkova, A.T.; Bobylev, V.V. Galactic masers: Kinematics, spiral structure and the disk dynamic state. Astrophys. Bull. 2017, 72, 122–140. [Google Scholar] [CrossRef] [Green Version]
  54. Marasco, A.; Fraternali, F.; van der Hulst, J.M.; Oosterloo, T. Distribution and kinematics of atomic and molecular gas inside the solar circle. Astron. Astrophys. 2017, 607, A106. [Google Scholar] [CrossRef] [Green Version]
  55. Reid, M.J.; Menten, K.M.; Brunthaler, A.; Zheng, X.W.; Dame, T.M.; Xu, Y.; Wu, Y.; Zhang, B.; Sanna, A.; Sato, M.; et al. Trigonometric Parallaxes of High Mass Star Forming Regions The Structure and Kinematics of the Milky Way. Astrophys. J. 2014, 783, 130. [Google Scholar] [CrossRef] [Green Version]
  56. Quiroga-Nunez, L.H.; van Langevelde, H.J.; Reid, M.J.; Green, J.A. Simulated Galactic methanol maser distribution to constrain Milky Way parameters. Astron. Astrophys. 2017, 604, A72. [Google Scholar] [CrossRef] [Green Version]
  57. Bobylev, V.V.; Bajkova, A.T. Galactic rotation curve and spiral density wave parameters from 73 masers. Astron. Lett. 2013, 39, 809–818. [Google Scholar] [CrossRef] [Green Version]
  58. Clemens, D.P. Massachusetts-Stony Brook Galactic plane CO survey: The galactic disk rotation curve. Astrophys. J. 1985, 295, 422–436. [Google Scholar] [CrossRef]
  59. Sofue, Y. Rotation Curve and Mass Distribution in the Galactic Center—From Black Hole to Entire Galaxy. Publ. Astron. Soc. Jpn. 2013, 65, 118. [Google Scholar] [CrossRef]
  60. Sofue, Y. Rotation and mass in the Milky Way and spiral galaxies. Publ. Astron. Soc. Jpn. 2017, 69, R1. [Google Scholar] [CrossRef] [Green Version]
  61. Dambis, A.K.; Mel’Nik, A.M.; Rastorguev, A.S. Rotation curve of the system of classical Cepheids and the distance to the galactic center. Astron. Lett. 1995, 21, 291–307. [Google Scholar]
  62. Lopez-Corredoira, M. Milky Way rotation curve from proper motions of red clump giants. Astron. Astrophys. 2014, 563, A128. [Google Scholar] [CrossRef]
  63. Peng, X.; Wu, Z.; Qi, Z.; Du, C.; Ma, J.; Zhou, X.; Jia, Y.; Wang, S. The Study of Galactic Disk Kinematics with SCUSS and SDSS Data. Publ. Astron. Soc. Pac. 2018, 130, 074102. [Google Scholar] [CrossRef]
  64. Bajkova, A.T.; Bobylev, V.V. Rotation curve and mass distribution in the Galaxy from the velocities of objects at distances up to 200 kpc. Astron. Lett. 2016, 42, 567–582. [Google Scholar] [CrossRef] [Green Version]
  65. Blum, R.D.; Carr, J.S.; Depoy, D.L.; Sellgren, K.; Terndrup, D.M. Radial Velocities of M Giants at 300 Parsec Projected Radius from the Galactic Center. Astrophys. J. 1994, 422, 111–117. [Google Scholar] [CrossRef]
  66. Rich, R.M. Kinematics and Abundances of K Giants in the Nuclear Bulge of the Galaxy. Astrophys. J. 1990, 362, 604–619. [Google Scholar] [CrossRef]
  67. Minniti, D.; White, S.D.M.; Olszewski, E.W.; Hill, J.M. Rotation of the Galactic Bulge. Astrophys. J. Lett. 1992, 393, 47–50. [Google Scholar] [CrossRef]
  68. Lewis, J.R.; Freeman, K.C. Kinematics and Chemical Properties of the Old Disk in the Galaxy. Astron. J. 1989, 97, 139–162. [Google Scholar] [CrossRef]
  69. Rich, R.M.; Reitzel, D.B.; Howard, C.D.; Zhao, H. The Bulge Radial Velocity Assay: Techniques and a Rotation Curve. Astrophys. J. 2007, 658, L29–L32. [Google Scholar] [CrossRef]
  70. Tremaine, S.; Gebhardt, K.; Bender, R.; Bower, G.; Dressler, A.; Faber, S.M.; Filippenko, A.V.; Green, R.; Grillmair, C.; Ho, L.C.; et al. The Slope of the Black Hole Mass versus Velocity Dispersion Correlation. Astrophys. J. 2002, 574, 740–753. [Google Scholar] [CrossRef]
  71. Guiglion, G.; Recio-Blanco, A.; de Laverny, P.; Kordopatis, G. The Gaia-ESO Survey: New constraints on the Galactic disc velocity dispersion and its chemical dependencies. Astron. Astrophys. 2015, 583, A91. [Google Scholar] [CrossRef] [Green Version]
  72. Kunder, A.; Koch, A.; Rich, R.M.; de Propris, R.; Howard, C.D.; Stubbs, S.A.; Johnson, C.I.; Shen, J.; Wang, Y.; Robin, A.C.; et al. The Bulge Radial Velocity Assay (BRAVA). II. Complete Sample and Data Release. Astron. J. 2012, 143, 57. [Google Scholar] [CrossRef] [Green Version]
  73. Tiede, G.P.; Terndrup, D.M. Kinematics, Metallicities, and Stellar Distributions in the Inner Disk and Bulge of the Milky Way. Astron. J. 1999, 118, 895–910. [Google Scholar] [CrossRef]
  74. Rafikov, R.R. The local axisymmetric instability criterion in a thin, rotating, multicomponent disc. Mon. Not. R. Astron. Soc. 2001, 323, 445–452. [Google Scholar] [CrossRef] [Green Version]
  75. Adibekyan, V.Z.; Figueira, P.; Santos, N.C.; Hakobyan, A.A.; Sousa, S.G.; Pace, G.; Delgado, M.E.; Robin, A.C.; Israelian, G.; Gonzalez Hernandez, J.I. Kinematics and chemical properties of the Galactic stellar populations. The HARPS FGK dwarfs sample. Astron. Astrophys. 2013, 554, A44. [Google Scholar] [CrossRef] [Green Version]
  76. Hayden, M.R.; Recio-Blanco, A.; de Laverny, P.; Mikolaitis, S.; Worley, C.C. The AMBRE project: The thick thin disk and thin thick disk of the Milky Way. Astron. Astrophys. 2017, 608, L1. [Google Scholar] [CrossRef] [Green Version]
  77. Gonzalez, O.A.; Gadotti, D. The Milky Way Bulge: Observed Properties and a Comparison to External Galaxies. Astrophys. Space Sci. Libr. 2016, 418, 199–232. [Google Scholar] [CrossRef]
  78. Valenti, E.; Zoccali, M.; Gonzalez, O.A.; Minniti, D.; Alonso-García, J.; Marchetti, E.; Hempel, M.; Renzini, A.; Rejkuba, M. Stellar density profile and mass of the Milky Way bulge from VVV data. Astron. Astrophys. 2016, 587, L6. [Google Scholar] [CrossRef]
  79. Portail, M.; Gerhard, O.; Wegg, C.; Ness, M. Dynamical modelling of the galactic bulge and bar: The Milky Way’s pattern speed, stellar and dark matter mass distribution. Mon. Not. R. Astron. Soc. 2017, 465, 1621–1644. [Google Scholar] [CrossRef] [Green Version]
  80. Bobylev, V.V.; Bajkova, A.T. Galactic orbits of globular clusters in the region of the Galactic bulge. Astron. Rep. 2017, 61, 551–559. [Google Scholar] [CrossRef] [Green Version]
  81. Michtchenko, T.A.; Lépine, J.R.D.; Barros, D.A.; Vieira, R.S.S. Combined dynamical effects of the bar and spiral arms in a Galaxy model. Application to the solar neighbourhood. Astron. Astrophys. 2018, 615, A10. [Google Scholar] [CrossRef] [Green Version]
  82. Goldreich, P.; Lynden-Bell, D., II. Spiral arms as sheared gravitational instabilities. Mon. Not. R. Astron. Soc. 1965, 130, 125–158. [Google Scholar] [CrossRef]
  83. Nakanishi, H.; Sofue, Y. Three-dimensional distribution of the ISM in the Milky Way galaxy. III. The total neutral gas disk. Publ. Astron. Soc. Jpn. 2016, 68, 5. [Google Scholar] [CrossRef] [Green Version]
  84. Khoperskov, S.; Zinchenko, I.; Avramov, B.; Khrapov, S.; Berczik, P.; Saburova, A.; Ishchenko, M.; Khoperskov, A.; Pulsoni, C.; Venichenko, Y.; et al. Extreme kinematic misalignment in IllustrisTNG galaxies: The origin, structure, and internal dynamics of galaxies with a large-scale counterrotation. Mon. Not. R. Astron. Soc. 2021, 500, 3870–3888. [Google Scholar] [CrossRef]
  85. Graziani, L.; Salvadori, S.; Schneider, R.; Kawata, D.; de Bennassuti, M.; Maselli, A. Galaxy formation with radiative and chemical feedback. Mon. Not. R. Astron. Soc. 2015, 449, 3137–3148. [Google Scholar] [CrossRef] [Green Version]
  86. Graziani, L.; de Bennassuti, M.; Schneider, R.; Kawata, D.; Salvadori, S. The history of the dark and luminous side of Milky Way-like progenitors. Mon. Not. R. Astron. Soc. 2017, 469, 1101–1116. [Google Scholar] [CrossRef] [Green Version]
  87. Grand, R.J.J.; Gomez, F.A.; Marinacci, F.; Pakmor, R.; Springel, V.; Campbell, D.J.R.; Frenk, C.S.; Jenkins, A.; White, S.D.M. The Auriga Project: The properties and formation mechanisms of discgalaxies across cosmic time. Mon. Not. R. Astron. Soc. 2017, 467, 179–207. [Google Scholar] [CrossRef] [Green Version]
  88. Hopkins, P.F.; Wetzel, A.; Keres, D.; Faucher-Giguere, C.-A.; Quataert, E.; Boylan-Kolchin, M.; Murray, N.; Hayward, C.C.; Garrison-Kimmel, S.; Feldmann, R.; et al. FIRE-2 simulations: Physics versus numerics in galaxy formation. Mon. Not. R. Astron. Soc. 2018, 480, 800–863. [Google Scholar] [CrossRef] [Green Version]
  89. Zasov, A.V.; Saburova, A.S.; Khoperskov, A.V.; Khoperskov, S.A. Dark matter in galaxies. Physics-Uspekhi 2017, 60, 3–39. [Google Scholar] [CrossRef] [Green Version]
  90. De Lucia, G. Lighting Up Dark Matter Haloes. Galaxies 2019, 7, 56. [Google Scholar] [CrossRef] [Green Version]
  91. Deason, A.J.; Erkal, D.; Belokurov, V.; Fattahi, A.; Gómez, F.A.; Grand, R.J.J.; Pakmor, R.; Xue, X.-X.; Liu, C.; Yang, C.; et al. The mass of the Milky Way out to 100 kpc using halo stars. Mon. Not. R. Astron. Soc. 2021, 501, 5964–5972. [Google Scholar] [CrossRef]
  92. Khoperskov, A.V.; Zasov, A.V.; Tyurina, N.V. Minimum Velocity Dispersion in Stable Stellar Disks. Numerical Simulations. Astron. Rep. 2003, 47, 357–376. [Google Scholar] [CrossRef]
  93. Khoperskov, S.A.; Khoperskov, A.V.; Khrykin, I.S.; Korchagin, V.I.; Casetti-Dinescu, D.I.; Girard, T.; van Altena, W.; Maitra, D. Global gravitationally organized spiral waves and the structure of NGC 5247. Mon. Not. R. Astron. Soc. 2012, 427, 1983–1993. [Google Scholar] [CrossRef] [Green Version]
  94. Smirnov, A.A.; Sotnikova, N.Y.; Koshkin, A.A. Simulations of slow bars in anisotropic disk systems. Astron. Lett. 2017, 43, 61–74. [Google Scholar] [CrossRef]
  95. Maureira-Fredes, C.; Amaro-Seoane, P. GRAVIDY, a GPU modular, parallel direct-summation N-body integrator: Dynamics with softening. Mon. Not. R. Astron. Soc. 2018, 473, 3113–3127. [Google Scholar] [CrossRef]
  96. Sawada, K.; Matsuda, T.; Hachisu, I. Spiral shocks on a Roche lobe overflow in a semi-detached binary system. Mon. Not. R. Astron. Soc. 1986, 219, 75–88. [Google Scholar] [CrossRef] [Green Version]
  97. Bisikalo, D.V.; Boyarchuk, A.A.; Kuznetsov, O.A. The effect of viscosity on the flow morphology in semidetached binary systems. Results of 3D simulations. II. Astron. Rep. 2000, 44, 26–35. [Google Scholar] [CrossRef] [Green Version]
  98. Khoperskov, S.A.; Khoperskov, A.V.; Eremin, M.A.; Butenko, M.A. Polygonal structures in a gaseous disk: Numerical simulations. Astron. Lett. 2011, 37, 563–575. [Google Scholar] [CrossRef] [Green Version]
  99. Dobbs, C.L.; Bonnell, I.A.; Clark, P.C. Centrally condensed turbulent cores: Massive stars or fragmentation? Mon. Not. R. Astron. Soc. 2005, 360, 2–8. [Google Scholar] [CrossRef] [Green Version]
  100. Vasiliev, E.O. Thermal instability in a collisionally cooled gas. Mon. Not. R. Astron. Soc. 2012, 419, 3641–3648. [Google Scholar] [CrossRef] [Green Version]
  101. Khoperskov, S.A.; Vasiliev, E.O.; Sobolev, A.M.; Khoperskov, A.V. The simulation of molecular clouds formation in the Milky Way. Mon. Not. R. Astron. Soc. 2013, 428, 2311–2320. [Google Scholar] [CrossRef] [Green Version]
  102. Vorobyov, E.I.; Recchi, S.; Hensler, G. Stellar hydrodynamical modeling of dwarf galaxies: Simulation methodology, tests, and first results. Astron. Astrophys. Suppl. Ser. 2015, 579, A9. [Google Scholar] [CrossRef] [Green Version]
  103. Vasiliev, E.O.; Shchekinov, Y.A. Supernova Remnants in the Hα and Hβ Lines. Astrophysics 2017, 60, 1–18. [Google Scholar] [CrossRef]
  104. Dobbs, C.L.; Pringle, J.E. The exciting lives of giant molecular clouds. Mon. Not. R. Astron. Soc. 2013, 432, 653–667. [Google Scholar] [CrossRef]
  105. Khoperskov, S.A.; Vasiliev, E.O.; Ladeyschikov, D.A.; Sobolev, A.M.; Khoperskov, A.V. Giant molecular cloud scaling relations: The role of the cloud definition. Mon. Not. R. Astron. Soc. 2016, 455, 1782–1795. [Google Scholar] [CrossRef] [Green Version]
  106. Duarte-Cabral, A.; Dobbs, C.L. What can simulated molecular clouds tell us about real molecular clouds? Mon. Not. R. Astron. Soc. 2016, 458, 3667–3683. [Google Scholar] [CrossRef] [Green Version]
  107. Monaghan, J.J. Smoothed particle hydrodynamics. Annu. Rev. Astron. Astrophys. 1992, 30, 543–574. [Google Scholar] [CrossRef]
  108. Khrapov, S.; Khoperskov, A. Smoothed-particle hydrodynamics models: Implementation features on GPUs. Commun. Comput. Inf. Sci. 2017, 793, 266–277. [Google Scholar] [CrossRef] [Green Version]
  109. Muller, M.; Charypar, D.; Gross, M. Particle-based fluid simulation for interactive applications. In Proceedings of the 2003 ACM SIGGRAPH/Eurographics Symposium on Computer Animation, San Diego, CA, USA, 26–27 July 2003; pp. 154–159. [Google Scholar]
  110. Monaghan, J.J. Smoothed particle hydrodynamics. Rep. Prog. Phys. 2005, 68, 1703–1759. [Google Scholar] [CrossRef]
  111. Desbrun, M.; Gascuel, M.P. Smoothed particles: A new paradigm for animating highly deformable bodies. In Computer Animation and Simulation; Eurographics Workshop on Computer Animation and Simulation (EGCAS): Poitiers, France, 1996; pp. 61–76. [Google Scholar]
  112. Khrapov, S.S.; Khoperskov, S.A.; Khoperskov, A.V. New features of parallel implementation of N-body problems on GPU. Bull. South Ural State Univ. Ser. Math. Model. Program. Comput. Softw. 2018, 11, 124–136. [Google Scholar] [CrossRef]
  113. Toomre, A. On the gravitational stability of a disk of stars. Astrophys. J. 1964, 139, 1217–1238. [Google Scholar] [CrossRef]
  114. Romeo, A.B.; Falstad, N. A simple and accurate approximation for the Q stability parameter in multicomponent and realistically thick discs. Mon. Not. R. Astron. Soc. 2013, 433, 1389–1397. [Google Scholar] [CrossRef] [Green Version]
  115. Sellwood, J.A.; Athanassoula, E. Unstable modes from galaxy simulations. Mon. Not. R. Astron. Soc. 1986, 221, 195–212. [Google Scholar] [CrossRef] [Green Version]
  116. Vallée, J.P. New Velocimetry and Revised Cartography of the Spiral Arms in the Milky Way—A Consistent Symbiosis. Astron. J. 2008, 135, 1301–1310. [Google Scholar] [CrossRef]
  117. Vallée, J.P. Constraining the pitch angle of the galactic spiral arms in the Milky Way. New Astron. Rev. 2017, 79, 49–58. [Google Scholar] [CrossRef] [Green Version]
  118. Nikiforov, I.I.; Veselova, A.V. Geometric Aspects and Testing of the Galactic Center Distance Determination from Spiral Arm Segments. Astron. Lett. 2018, 44, 81–102. [Google Scholar] [CrossRef] [Green Version]
  119. Dame, T.M.; Elmegreen, B.G.; Cohen, R.S.; Thaddeus, P. The largest molecular cloud complexes in the first galactic quadrant. Astrophys. J. 1986, 305, 892–908. [Google Scholar] [CrossRef]
  120. Norman, C.A.; Sellwood, J.A.; Hasan, H. Bar Dissolution and Bulge Formation: An Example of Secular Dynamical Evolution in Galaxies. Astrophys. J. 1996, 462, 114–124. [Google Scholar] [CrossRef]
  121. Kataria, S.K.; Das, M. A study of the effect of bulges on bar formation in disc galaxies. Mon. Not. R. Astron. Soc. 2018, 475, 1653–1664. [Google Scholar] [CrossRef]
  122. Ostriker, J.P.; Peebles, P.J.E. A Numerical Study of the Stability of Flattened Galaxies: Or, can Cold Galaxies Survive? Astrophys. J. 1973, 186, 467–480. [Google Scholar] [CrossRef]
  123. Polyacheko, E.V.; Berczik, P.; Just, A. On the bar formation mechanism in galaxies with cuspy bulges. Mon. Not. R. Astron. Soc. 2016, 462, 3727–3738. [Google Scholar] [CrossRef]
  124. Athanassoula, E. Bar-Halo Interaction and Bar Growth. Astrophys. J. 2002, 569, L83–L86. [Google Scholar] [CrossRef]
  125. Long, S.; Shlosman, I.; Heller, C. Secular damping of stellar bars in spinning dark matter Halos. Astrophys. J. Lett. 2014, 783, L18. [Google Scholar] [CrossRef] [Green Version]
  126. Collier, A.; Shlosman, I.; Heller, C. Stellar bars in counter-rotating dark matter haloes: The role of halo orbit reversals. Mon. Not. R. Astron. Soc. 2019, 489, 3102–3115. [Google Scholar] [CrossRef] [Green Version]
  127. Bajkova, A.T.; Carraro, G.; Korchagin, V.I.; Budanova, N.O.; Bobylev, V.V. Milky Way Subsystems from Globular Cluster Kinematics Using Gaia DR2 and HST Data. Astrophys. J. 2020, 895, 69–87. [Google Scholar] [CrossRef]
  128. Chemin, L.; Renaud, F.; Soubiran, C. Incorrect rotation curve of the Milky Way. Astron. Astrophys. 2015, 578, A14. [Google Scholar] [CrossRef]
  129. Renaud, F.; Bournaud, F.; Emsellem, E.; Elmegreen, B.; Teyssier, R.; Alves, J.; Chapon, D.; Combes, F.; Dekel, A.; Gabor, J.; et al. A sub-parsec resolution simulation of the Milky Way: Global structure of the interstellar medium and properties of molecular clouds. Mon. Not. R. Astron. Soc. 2013, 436, 1836–1851. [Google Scholar] [CrossRef] [Green Version]
  130. Deg, N.; Widrow, L.M.; Randriam Ampandry, T.; Carignan, C. GALACTICS with gas. Mon. Not. R. Astron. Soc. 2019, 486, 5391–5399. [Google Scholar] [CrossRef]
  131. Korchagin, V.; Kikuchi, N.; Miyama, S.M.; Orlova, N.; Peterson, B.A. Global spiral modes in NGC 1566: Observations and theory. Astrophys. J. 2000, 541, 565–578. [Google Scholar] [CrossRef] [Green Version]
  132. Toomre, A. Theories of spiral structure. Annu. Rev. Astron. Astrophys. 1977, 15, 437–478. [Google Scholar] [CrossRef]
  133. Wahde, M.; Donner, K.J. Determination of the orbital parameters of the M 51 system using a genetic algorithm. Astron. Astrophys. 2001, 379, 115–124. [Google Scholar] [CrossRef] [Green Version]
  134. Sellwood, J.A.; Calberg, R.G. Spiral instabilities: Linear and nonlinear effects. Monthly Notices of the Royal Astronomical Society. arXiv 2020, arXiv:2011.03041. [Google Scholar]
  135. Laughlin, G.; Korchagin, V.; Adams, F. Spiral Mode Saturation in Self-Gravitating Disks. Astrophys. J. 1997, 477, 410–423. [Google Scholar] [CrossRef]
Figure 1. The rotation curve of the Milky Way galaxy based on different observational data. The circular velocity V c and rotation velocity of the gas disk V g are shown: curve 1 is the circular velocity from WISE, 2MASS, and Gaia [50]; 2 is the result of compiling a large amount of data from various spectroscopic and trigonometric measurements from radio to optical wavelengths [51]; 3 is the reconstructed circular velocity using Gaia Collaboration data and simulations that are based on the Jeans equation [52]; 4—maser observations from [53]; gaseous rotation curves inside the solar circle [54]: atomic hydrogen (curve 5), molecular hydrogen (curve 6); the maser rotation curve from Reid et al. [55] is shown by line 7.
Figure 1. The rotation curve of the Milky Way galaxy based on different observational data. The circular velocity V c and rotation velocity of the gas disk V g are shown: curve 1 is the circular velocity from WISE, 2MASS, and Gaia [50]; 2 is the result of compiling a large amount of data from various spectroscopic and trigonometric measurements from radio to optical wavelengths [51]; 3 is the reconstructed circular velocity using Gaia Collaboration data and simulations that are based on the Jeans equation [52]; 4—maser observations from [53]; gaseous rotation curves inside the solar circle [54]: atomic hydrogen (curve 5), molecular hydrogen (curve 6); the maser rotation curve from Reid et al. [55] is shown by line 7.
Galaxies 09 00029 g001
Figure 2. Radial component of the velocity dispersion of stars as a function of radius. 1—velocity dispersion of late type giants within 300 pc radius [65], 2—[66], 3—[67], 4—old disk K-giants [68], 5—decrease by 10 percents from [68] due to taking into account younger stellar populations, 6—approximation by c r = 91 [km s 1 ] exp ( r / 7.4 [kpc]), 7—[69], 8—outside the solar circle [38] taking into account a decrease by 16% for RGB-stars, 9—velocity dispersion within bulge region [70], 10—as based on Gaia-ESO Survey Guiglion et al. [71], 11 and 12 for giant stars from Gaia in layer 200 pc z 200 pc for negative and positive azimuths, respectively [71], 13 is the central velocity dispersion of bulge [72], 14 is the velocity dispersion from Gaia data analysis [50].
Figure 2. Radial component of the velocity dispersion of stars as a function of radius. 1—velocity dispersion of late type giants within 300 pc radius [65], 2—[66], 3—[67], 4—old disk K-giants [68], 5—decrease by 10 percents from [68] due to taking into account younger stellar populations, 6—approximation by c r = 91 [km s 1 ] exp ( r / 7.4 [kpc]), 7—[69], 8—outside the solar circle [38] taking into account a decrease by 16% for RGB-stars, 9—velocity dispersion within bulge region [70], 10—as based on Gaia-ESO Survey Guiglion et al. [71], 11 and 12 for giant stars from Gaia in layer 200 pc z 200 pc for negative and positive azimuths, respectively [71], 13 is the central velocity dispersion of bulge [72], 14 is the velocity dispersion from Gaia data analysis [50].
Galaxies 09 00029 g002
Figure 3. Radial dependence of the Q T parameters on the radius in two models. Curve 1— Q T parameter of the stellar disk (the same for both 706 and 707 models), curve 2— Q T g parameter for the gaseous component in model 706, curve 3—two-fluid Q T parameter for model 706, 4— Q T g parameter of gaseous disk in model 707, 5—two-fluid Q T parameter in model 707, curve 6— Q T parameter of the stellar disk in model 730, curve 7— Q T g parameter for the gaseous component in model 730, curve 8—two-fluid Q T parameter for model 730.
Figure 3. Radial dependence of the Q T parameters on the radius in two models. Curve 1— Q T parameter of the stellar disk (the same for both 706 and 707 models), curve 2— Q T g parameter for the gaseous component in model 706, curve 3—two-fluid Q T parameter for model 706, 4— Q T g parameter of gaseous disk in model 707, 5—two-fluid Q T parameter in model 707, curve 6— Q T parameter of the stellar disk in model 730, curve 7— Q T g parameter for the gaseous component in model 730, curve 8—two-fluid Q T parameter for model 730.
Galaxies 09 00029 g003
Figure 4. Dynamics of gaseous (orange frames) and stellar (blue frames) disks at different times for model 706 shown sequentially, from left to right and from top to bottom, taken at times t = { 3 , 7 , 11 , 15 , 19 , 23 } .
Figure 4. Dynamics of gaseous (orange frames) and stellar (blue frames) disks at different times for model 706 shown sequentially, from left to right and from top to bottom, taken at times t = { 3 , 7 , 11 , 15 , 19 , 23 } .
Galaxies 09 00029 g004
Figure 5. Time dependence of global amplitudes of m = 2 , 3 , 4 Fourier harmonics in model 706 for the three regions of the stellar disk: (1) r 1.8 kpc, (2) 3.6 r 5.4 (kpc), (3) 7.2 r 9 (kpc).
Figure 5. Time dependence of global amplitudes of m = 2 , 3 , 4 Fourier harmonics in model 706 for the three regions of the stellar disk: (1) r 1.8 kpc, (2) 3.6 r 5.4 (kpc), (3) 7.2 r 9 (kpc).
Galaxies 09 00029 g005
Figure 6. The same as in Figure 5, but for the gaseous disk.
Figure 6. The same as in Figure 5, but for the gaseous disk.
Galaxies 09 00029 g006
Figure 7. Morphology of the spiral structure of gaseous (orange frames) and stellar (blue frames) disks in different models (404, 405, 455, 474, 600, 602, 610, 701, and 730) at time t = 15 .
Figure 7. Morphology of the spiral structure of gaseous (orange frames) and stellar (blue frames) disks in different models (404, 405, 455, 474, 600, 602, 610, 701, and 730) at time t = 15 .
Galaxies 09 00029 g007
Figure 8. Azimuthally averaged rotation curve of the gas in numerical model 706 at different moments of time compared to observational data: 1—HI, 2—masers, 3—numerical model at t = 0 , 4— t = 5 , 5— t = 10 , 6— t = 15 , 7— t = 20 ).
Figure 8. Azimuthally averaged rotation curve of the gas in numerical model 706 at different moments of time compared to observational data: 1—HI, 2—masers, 3—numerical model at t = 0 , 4— t = 5 , 5— t = 10 , 6— t = 15 , 7— t = 20 ).
Galaxies 09 00029 g008
Figure 9. Azimuthally averaged radial dependencies of the surface densities of gas σ g and stars σ at different times. The times shown are the same as in Figure 8.
Figure 9. Azimuthally averaged radial dependencies of the surface densities of gas σ g and stars σ at different times. The times shown are the same as in Figure 8.
Galaxies 09 00029 g009
Figure 10. Azimuthal dependences of the logarithm of the density distributions for stellar (upper panel) and for gaseous (lower panel) components taken at three different radii at time t = 10 (experiment 706). The azimuthal angle φ is in radians.
Figure 10. Azimuthal dependences of the logarithm of the density distributions for stellar (upper panel) and for gaseous (lower panel) components taken at three different radii at time t = 10 (experiment 706). The azimuthal angle φ is in radians.
Galaxies 09 00029 g010
Figure 11. Dependence on time of the velocity dispersion c r ( t ) for different models.
Figure 11. Dependence on time of the velocity dispersion c r ( t ) for different models.
Galaxies 09 00029 g011
Figure 12. Azimuthal dependencies of the velocity dispersion c r ( φ ) at different radii for model 707. The azimuthat angle φ is measured in radians.
Figure 12. Azimuthal dependencies of the velocity dispersion c r ( φ ) at different radii for model 707. The azimuthat angle φ is measured in radians.
Galaxies 09 00029 g012
Figure 13. The typical profiles of the azimuthal velocities of the stellar (left frame) and gaseous disks (right frame) along the minor (curves 3 and 4) and major (curves 5 and 6) axis of the bar, taken at t = 15 for model 706. We show two rotation curves along major and minor semi-axes of the bar. For comparison, the azimuthally averaged velocity profiles are also shown. Curve 1 shows the azimuthal velocity at t = 0 , while 2 is the azimuthally averaged rotational curve at time t = 15 . The vertical arrow shows the approximate size of the bar.
Figure 13. The typical profiles of the azimuthal velocities of the stellar (left frame) and gaseous disks (right frame) along the minor (curves 3 and 4) and major (curves 5 and 6) axis of the bar, taken at t = 15 for model 706. We show two rotation curves along major and minor semi-axes of the bar. For comparison, the azimuthally averaged velocity profiles are also shown. Curve 1 shows the azimuthal velocity at t = 0 , while 2 is the azimuthally averaged rotational curve at time t = 15 . The vertical arrow shows the approximate size of the bar.
Galaxies 09 00029 g013
Table 1. Model parameters.
Table 1. Model parameters.
Experiment M d × 10 10 M r d kpc σ 0 M / pc 2 σ M / pc 2 h pc M h × 10 10 M a kpc μ × 10 10 M M b b kpc r b max kpc M g × 10 10 M μ g M / pc 2 σ g min Q T min Q Tg min Q T
123456789101112131415161718
4554.802.25150842.73005.533.61.151.320.21.20.240.056.721.17/1.055.061.01
4703.612.251136323006.283.51.740.9560.160.70.610.0858.571.39/1.213.41.16
4733.612.251136323006.283.51.740.9560.160.70.610.0858.571.25/1.123.391.07
4743.612.251136321926.283.51.740.9560.160.70.610.0858.571.12/1.023.360.98
4753.612.251136322436.283.51.740.9560.160.70.610.0858.571.26/1.123.371.08
4814.872.25153143.93005.513.51.130.9560.160.70.6330.0659.061.16/1.033.811
4824.872.25153143.92435.513.51.130.9560.160.70.6330.0659.061.17/1.033.71
4905.012.251576443005.232.251.040.6720.160.70.6510.0659.071.30/1.093.631.06
6003.802.25119433.73006.033.81.590.9640.301.40.6420.0859.021.36/1.223.31.13
6013.802.25119433.72436.033.81.590.9640.301.40.6420.0859.021.23/1.133.281.05
6023.802.25119433.73006.033.81.590.9640.301.40.6420.0859.021.49/1.273.311.21
6104.102.25129036.83005.972.811.461.1580.301.50.5330.1307.581.40/1.244.11.19
7013.722.25117033.43155.903.331.591.1650.301.40.4840.13013.71.27/1.142.181.05
7023.722.25117033.43155.903.331.591.1650.301.40.4840.13013.70.98/0.952.180.88
7033.722.25117033.43885.903.331.591.1650.301.40.4840.13013.71.32/1.192.191.09
7053.722.25117033.43885.903.01.591.1590.311.60.4840.13013.71.34/1.212.181.11
7063.722.25117033.43886.023.001.621.010.391.80.4840.13013.71.33/1.22.191.10
7073.722.25117033.43886.023.001.621.010.391.80.3910.10511.11.33/1.22.721.12
7202.723.0048133.430010.323.013.791.640.212.00.4120.15113.71.94/1.792.941.58
7303.082.6072533.430010.323.013.791.640.212.00.5050.16413.71.59/1.442.581.32
1—number of experiment; 2—the disk mass inside the optical radius ( r R o p t ); 3—the exponential scale of the stellar disk; 4—the central surface density of the stellar disk (See (1)); 5—the surface density of the stellar disk near the solar neighborhood; 6—the vertical scale height (See (2)); 7—the halo mass inside the optical radius; 8—the radial dark halo scale in the isothermal model; 9—the relative mass of the dark halo in units of the stellar disk within the optical radius; 10—the bulge mass; 11—the bulge scale length; 12—the bulge size; 13—the gaseous disk mass within r 2 R o p t ; 14—the relative mass of the gaseous disk in units of the stellar disk mass; 15— the surface density of the gaseous disk near the solar neighborhood; 16—the minimum value of the Toomre parameter for the stellar disk; 17—the minimum value of the Toomre parameter for the gaseous disk; 18—the minimum value of the Toomre parameter for the two-component model.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Khrapov, S.; Khoperskov, A.; Korchagin, V. Modeling of Spiral Structure in a Multi-Component Milky Way-Like Galaxy. Galaxies 2021, 9, 29. https://doi.org/10.3390/galaxies9020029

AMA Style

Khrapov S, Khoperskov A, Korchagin V. Modeling of Spiral Structure in a Multi-Component Milky Way-Like Galaxy. Galaxies. 2021; 9(2):29. https://doi.org/10.3390/galaxies9020029

Chicago/Turabian Style

Khrapov, Sergey, Alexander Khoperskov, and Vladimir Korchagin. 2021. "Modeling of Spiral Structure in a Multi-Component Milky Way-Like Galaxy" Galaxies 9, no. 2: 29. https://doi.org/10.3390/galaxies9020029

APA Style

Khrapov, S., Khoperskov, A., & Korchagin, V. (2021). Modeling of Spiral Structure in a Multi-Component Milky Way-Like Galaxy. Galaxies, 9(2), 29. https://doi.org/10.3390/galaxies9020029

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop