Search

Context. Recent observational data show that the Milky Way (MW) galaxy contains about 170 globular clusters (GCs). A fraction of them is likely formed in dwarf galaxies accreted onto the MW in the past, while the remaining of clusters are formed in-situ. Therefore, different parameters, including orbits, of the globular clusters is a valuable tool for studying the Milky Way evolution. However, since the evolution of the 3D mass distribution of the MW is poorly constrained, the orbits of the clusters are usually calculated in static potentials. Aims. In this work, we study the evolution of the GCs in several external potentials, where we aim to quantify the effects of the evolving galaxy potential on the orbits of the GCs. Methods. For the orbits calculation we used five MW-like potentials from IllustrisTNG-100 simulation. The orbits of 159 GCs were integrated using a high-order N-body parallel dynamic code phi-GPU, with initial conditions obtained from recent Gaia DR3 catalogues. Results. We provide a classification of the GCs orbits according to their 3D shapes and association with different components of the MW (disk, halo, bulge). We also found that the globular clusters in the external potentials have roughly similar energy-angular momentum distributions at the present time. However, both total energy and total angular momentum of the GCs are not conserved due to time-varying nature of the potentials. In some extreme cases, the total energy can change up to 40% (18 objects) over the last 5 Gyr of evolution. We found that the in-situ formed GCs are less affected by the evolution of the TNG potentials as compared to the clusters which are likely formed ex-situ. Therefore, our results suggest that time-varying potentials significantly affect the orbits of the GC, thus making it vital for understanding the formation of the MW., Comment: 19 pages, 17 figures, 2 tables, accepted for publication in Astronomy and Astrophysics

Using the data from Gaia (ESA) Data Release 2 we performed the orbital calculations of globular clusters (GCs) of the Milky Way. To explore possible close encounters (or collisions) between the GCs, using our own developed high-order phi-GRAPE code, we integrated backward and forward the orbits of 119 objects with reliable positions and proper motions. In calculations, we adopted a realistic axisymmetric Galactic potential (bulge + disk + halo). Using different impact conditions, we found four pairs of the six GCs that may have experienced an encounter within twice the sum of the half-mass radii (collisions) over the last 5Gyr: Terzan3 - NGC 6553, Terzan 3 - NGC 6218, Liller 1 - NGC 6522 and Djorg 2 - NGC 6553., Comment: Accepted for publication in Advances in Astronomy and Space Physics (http://aasp.kiev.ua). arXiv admin note: substantial text overlap with arXiv:2201.06891

Understanding radial migration is a crucial point to build relevant chemical and dynamical evolution models of the Milky Way disk. In this paper, we analyze a high-resolution N-body simulation of a Milky Way-type galaxy to study the role that the slowing down of a stellar bar has is generating migration from the inner to the outer disk. Stellar particles are trapped by the main resonances (corotation and Outer Lindblad resonance) which then propagate outwards across the disk due to the bar slowing down. Once the bar strength reaches its maximal amplitude, some of the stars, delivered to the outer disk, escape the resonances and some of them settle on nearly circular orbits. The number of the escaped stars gradually increases also due to the decrease of the bar strength when the boxy/peanut bulge forms. We show that this mechanism is not limited only to stars on nearly circular orbits: also stars initially on more eccentric orbits can be transferred outwards (out to the OLR location) and can end up on nearly circular orbits. Therefore, the propagation of the bar resonances outwards can induce the circularization of the orbits of some of the migrating stars. The mechanism investigated in this paper can explain the presence of metal-rich stars at the solar vicinity and more generally in the outer galactic disk. Our dynamical model predicts that up to 3% of stars in between of corotation and the OLR can be formed in the innermost region of the Milky Way. The epoch of the Milky Way bar formation can be potentially constrained by analyzing the age distribution of the most metal-rich stars at the solar vicinity., Comment: 9 pages, 8 figures, A&A in press

In this paper we introduce a new method for analysing Milky Way phase-space which allows us to reveal the imprint left by the Milky Way bar and spiral arms on the stars with full phase-space data in Gaia Data Release 2. The unprecedented quality and extended spatial coverage of these data enable us to discover six prominent stellar density structures in the disc to a distance of 5 kpc from the Sun. Four of these structures correspond to the spiral arms detected previously in the gas and young stars (Scutum-Centaurus, Sagittarius, Local and Perseus). The remaining two are associated with the main resonances of the Milky Way bar where corotation is placed at around 6.2 kpc and the outer Lindblad resonance beyond the Solar radius, at around 9 kpc. For the first time we provide evidence of the imprint left by spiral arms and resonances in the stellar densities not relying on a specific tracer, through enhancing the signatures left by these asymmetries. Our method offers new avenues for studying how the stellar populations in our Galaxy are shaped., Comment: 10 pages, 6 figures, accepted for publication in A&A Letters

We present a scenario of the chemical enrichment of the solar neighborhood that solves the G-dwarf problem by taking into account constraints on a larger scale. We argue that the Milky Way disk within 10 kpc has been enriched to solar metallicity by a massive stellar population: the thick disk, which itself formed from a massive turbulent gaseous disk. The key new consideration is that the pre-enrichment provided by the thick disk is not related to the mass fraction of this stellar population at the solar radius, as is classically assumed in inside-out scenarios, but is actually related to the formation of the entire massive thick disk, due to the vigorous gas phase mixing that occurred during its formation. Hence, the fact that this population represents only 15-25% of the local stellar surface density today is irrelevant for `solving' the G-dwarf problem. The only condition for this scenario to work is that the thick disk was formed from a turbulent gaseous disk that permitted a homogeneous -- not radially dependent -- distribution of metals, allowing the solar ring to be enriched to solar metallicity. At the solar radius, the gas flowing from the outer disk combined with the solar metallicity gas left over from thick disk formation, providing the fuel necessary to form the thin disk at the correct metallicity to solve the G-dwarf problem. Chemical evolution at R$>$6 kpc, and in particular beyond the solar radius, can be reproduced with the same scheme. These results imply that the local metallicity distribution is not connected to the gas accretion history of the Milky Way. Finally, we argue that the Sun is the result of the evolution typical of stars in the disk beyond $\sim$6 kpc (i.e., also undergoing dilution), and has none of the characteristics of inner disk stars. [Abridged], Comment: 15 pages, 7 figures. Accepted for publication in A&A, published version

We explore the origin of phase-space substructures revealed by the second Gaia data release in the disc of the Milky Way, such as the ridges in the $V_{\phi}$-$r$ plane, the undulations in the $V_{\phi}$-$r$-$V_r$ space and the streams in the $V_{\phi}$-$V_r$ plane. We use a collisionless N-body simulation with co-spatial thin and thick discs, along with orbit integration, to study the orbital structure close to the Outer Lindblad Resonance (OLR) of the bar. We find that a prominent, long-lived ridge is formed in the $V_{\phi}$-$r$ plane due to the OLR which translates to streams in the $V_{\phi}$-$V_r$ plane and examine which closed periodic and trapped librating orbits are responsible for these features. We find that orbits which carry out small librations around the $x_1(1)$ family are preferentially found at negative $V_r$, giving rise to a `horn'-like feature, while orbits with larger libration amplitudes, trapped around the $x_1(2)$ and $x_1(1)$ families, constitute the positive $V_r$ substructure, i.e. the Hercules-like feature. This changing libration amplitude of orbits will translate to a changing ratio of thin/thick disc stars, which could have implications on the metallicity distribution in this plane. We find that a scenario in which the Sun is placed close to the OLR gives rise to a strong asymmetry in $V_r$ in the $V_{\phi}$-$V_r$ plane (i.e. Hercules vs. `the horn') and subsequently to undulations in the $V_{\phi}$-$r$-$V_r$ space. We also explore a scenario in which the Sun is placed closer to the bar corotation and find that the bar perturbation $alone$ cannot give rise to the these features., Comment: 17 pages, accepted for publication in MNRAS

In the MW bulge, metal-rich stars form a strong bar and are more peanut-shaped than metal-poor stars. It has been recently claimed that this behavior is driven by the initial in-plane radial velocity dispersion of these populations, rather than by their initial vertical random motions. This has led to the suggestion that a thick disc is not necessary to explain the characteristics of the MW bulge. We rediscuss this issue by analyzing two dissipationless N-body simulations of boxy/peanut (b/p)-shaped bulges formed from composite stellar discs, made of kinematically cold and hot stellar populations, and we conclude that initial vertical random motions are as important as in-plane random motions in determining the relative contribution of cold and hot disc populations with height above the plane, the metallicity and age trends. Previous statements emphasizing the dominant role of in-plane motions in determining these trends are not confirmed. However, differences exist in the morphology and strength of the resulting b/p-shaped bulges: a model where disc populations have initially only different in-plane random motions, but similar thickness, results into a b/p bulge where all populations have a similar peanut shape, independently on their initial kinematics, or metallicity. We discuss the reasons behind these differences, and also predict the signatures that these two extreme initial conditions would leave on the vertical age and metallicity gradients of disc stars, outside the bulge region. We conclude that a metal-poor, kinematically (radial and vertical) hot component, that is a thick disc, is necessary in the MW before bar formation, supporting the scenario traced in previous works. [abridged], Comment: 17 pages, 16 figures, A&A, in press

Previous studies based on the analysis of Gaia DR2 data have revealed that accreted stars, possibly originating from a single progenitor satellite, are a significant component of the halo of our Galaxy, potentially constituting most of the halo stars at $\rm [Fe/H] < -1$ within a few kpc from the Sun and beyond. In this paper, we couple astrometric data from Gaia DR2 with elemental abundances from APOGEE DR14 to characterize the kinematics and chemistry of in-situ and accreted populations up to $\rm [Fe/H] \sim -2$. Accreted stars appear to significantly impact the Galactic chemo-kinematic relations, not only at $\rm [Fe/H] < -1$, but also at metallicities typical of the thick and metal-poor thin discs. They constitute about 60% of all stars at $\rm [Fe/H] < -1$, the remaining 40% being made of (metal-weak) thick disc stars. We find that the stellar kinematic fossil record shows the imprint left by this accretion event which heated the old Galactic disc. We are able to age-date this kinematic imprint, showing that the accretion occurred between 9 and 11 Gyr ago, and that it led to the last significant heating of the Galactic disc. An important fraction of stars with abundances typical of the (metal-rich) thick disc, and heated by this interaction, is now found in the Galactic halo. Indeed about half of the kinematically defined halo at few kpc from the Sun is composed of metal-rich thick disc stars. Moreover, we suggest that this metal-rich thick disc component dominates the stellar halo of the inner Galaxy. The new picture that emerges from this study is one where the standard non-rotating in-situ halo population, the collapsed halo, seems to be more elusive than ever., Comment: 25 pages, 24 figures, accepted for publication on A&A

This paper focuses on the parallel implementation of a direct $N$-body method~(particle-particle algorithm) and the application of multiple GPUs for galactic dynamics simulations. Application of a hybrid OpenMP-CUDA technology is considered for models with a number of particles $N \sim 10^5 \div 10^7$. By means of $N$-body simulations of gravitationally unstable stellar galactic we have investigated the algorithms parallelization efficiency for various Nvidia Tesla graphics processors~(K20, K40, K80). Particular attention was paid to the parallel performance of simulations and accuracy of the numerical solution by comparing single and double floating-point precisions~(SP and DP). We showed that the double-precision simulations are slower by a factor of~$1.7$ than the single-precision runs performed on Nvidia Tesla K-Series processors. We also claim that application of the single-precision operations leads to incorrect result in the evolution of the non-axisymmetric gravitating $N$-body systems. In particular, it leads to significant quantitative and even qualitative distortions in the galactic disk evolution. For instance, after $10^4$ integration time steps for the single-precision numbers the total energy, momentum, and angular momentum of a system with $N = 2^{20}$ conserve with accuracy of $10^{-3}$, $10^{-2}$ and $10^{-3}$ respectively, in comparison to the double-precision simulations these values are $10^{-5}$, $10^{-15}$ and $10^{-13}$, respectively. Our estimations evidence in favour of usage of the second-order accuracy schemes with double-precision numbers since it is more efficient than in the fourth-order schemes with single-precision numbers., Comment: 12 pages, 7 figures

There is a long-standing debate on the origin of the metal-poor stellar populations of the Milky Way (MW) bulge, with the two leading scenarios being that these populations are either i) part of a classical metal-poor spheroid or ii) the same population as the chemically defined thick disc seen at the Solar neighbourhood. Here we test whether the latter scenario can reproduce the observed chemical properties of the MW bulge. To do so we compare an N-body simulation of a composite (thin+thick) stellar disc -- which evolves secularly to form a bar and a boxy/peanut (b/p) bulge -- to data from APOGEE DR13. This model, in which the thick disc is massive and centrally concentrated, can reproduce the morphology of the metal-rich and metal-poor stellar populations in the bulge, as well as the mean metallicity and [$\alpha$/Fe] maps as obtained from the APOGEE data. It also reproduces the trends, in both longitude and latitude, of the bulge metallicity distribution function (MDF). Additionally, we show that the model predicts small but measurable azimuthal metallicity variations in the inner disc due to the differential mapping of the thin and thick disc in the bar. We therefore see that the chemo-morphological relations of stellar populations in the MW bulge are naturally reproduced by mapping the thin and thick discs of the inner MW into a b/p., Comment: 15 pages, 14 figures. Accepted for publication in A&A

This letter studies the formation of azimuthal metallicity variations in the disks of spiral galaxies in the absence of initial radial metallicity gradients. Using high-resolution $N$-body simulations, we model composite stellar discs, made of kinematically cold and hot stellar populations, and study their response to spiral arm perturbations. We find that, as expected, disk populations with different kinematics respond differently to a spiral perturbation, with the tendency for dynamically cooler populations to show a larger fractional contribution to spiral arms than dynamically hotter populations. By assuming a relation between kinematics and metallicity, namely the hotter the population, the more metal-poor it is, this differential response to the spiral arm perturbations naturally leads to azimuthal variations in the mean metallicity of stars in the simulated disk. Thus, azimuthal variations in the mean metallicity of stars across a spiral galaxy are not necessarily a consequence of the reshaping, by radial migration, of an initial radial metallicity gradient. They indeed arise naturally also in stellar disks which have initially only a negative vertical metallicity gradient., Comment: 5 pages, 6 figures, Accepted for publication in A&A Letters

Dark matter in galaxies, its abundance, and its distribution remain a subject of long-standing discussion, especially in view of the fact that neither dark matter particles nor dark matter bodies have yet been found. Experts' opinions range from a very large number of completely dark galaxies exist to nonbaryonic dark matter does not exist at all in any significant amounts. We discuss astronomical evidence for the existence of dark matter and its connection with visible matter and examine attempts to estimate its mass and distribution in galaxies from photometry, dynamics, gravitational lensing, and other observations (the cosmological aspects of the existence of dark matter are not considered in this review). In our view, the presence of dark matter in and around galaxies is a well-established fact. We conclude with an overview of mechanisms by which a dark halo can influence intragalactic processes., Comment: 82 pages, 35 figures

We examine the metallicity trends in the Milky Way (MW) bulge - using APOGEE DR13 data - and explore their origin by comparing two N-body models of isolated galaxies which develop a bar and a boxy/peanut (b/p) bulge. Both models have been proposed as scenarios for reconciling a disc origin of the MW bulge with a negative vertical metallicity gradient. The first model is a superposition of co-spatial, i.e. overlapping, disc populations with different scaleheights, kinematics and metallicities. In this model the thick, metal-poor, and centrally concentrated disc populations contribute significantly to the stellar mass budget in the inner galaxy. The second model is a single disc with an initial steep radial metallicity gradient, which is mapped by the bar into the b/p bulge in such a way that the vertical metallicity gradient of the MW bulge is reproduced -- as shown already in previous works in the literature. However, as we show here, the latter model does not reproduce the positive longitudinal metallicity gradient of the inner disc, nor the metal-poor innermost regions of the Bulge seen in the data. On the other hand, the model with co-spatial thin and thick disc populations reproduces all the aforementioned trends. We therefore see that it is possible to reconcile a (primarily) disc origin for the MW bulge with the observed trends in metallicity by mapping the inner thin and thick discs of the MW into a b/p. For this scenario to reproduce the observations, the $\alpha$-enhanced, metal-poor, thick disc populations must have a significant mass contribution in the inner regions -- as has been suggested for the Milky Way., Comment: 6 pages, accepted for publication in A&A Letters

We study numerically large-scale magnetic field evolution and its enhancement in gaseous disks of spiral galaxies. We consider a set of models with the various spiral pattern parameters and the initial magnetic field strength with taking into account gas self-gravity and cooling/heating. In agreement with previous studies, we find out that galactic magnetic field is mostly aligned with gaseous structures, however small-scale gaseous structures (spurs and clumps) are more chaotic than the magnetic field structure. In spiral arms magnetic field strongly coexists with the gas distribution, in the inter-arm region we see filamentary magnetic field structure. Simulations reveal the presence of the small-scale irregularities of the magnetic field as well as the reversal of magnetic field at the outer edge of the large-scale spurs. We provide evidences that the magnetic field in the spiral arms has a stronger mean-field component, and there is a clear inverse correlation between gas density and plasma-beta parameter, compared to the rest of the disk with a more turbulent component of the field and an absence of correlation between gas density and plasma-beta. We show the mean field growth up to 3-10$\mu G$ in the cold gas during several rotation periods (500-800 Myr), whereas ratio between azimuthal and radial field is equal to 4/1. Mean field strength increases by a factor of 1.5-2.5 for models with various spiral pattern parameters. Random magnetic field component can reach up to 25 % from the total strength. By making an analysis of the time-depended evolution of radial Poynting flux we point out that the magnetic field strength is enhanced stronger at the galactic outskirts which is due to the radial transfer of magnetic energy by the spiral arms pushing the magnetic field outward. Our results also support the presence of sufficient conditions for development of MRI at distances >11 kpc., Comment: 15 pages, 15 figures, accepted for publication in A&A

Galaxy surveys have suggested that rapid and sustained decrease in the star-formation rate, "quenching", in massive disk galaxies is frequently related to the presence of a bar. Optical and near-IR observations reveal that nearly 60% of disk galaxies in the local universe are barred, thus it is important to understand the relationship between bars and star formation in disk galaxies. Recent observational results imply that the Milky Way quenched about 9-10 Gyr ago, at the transition between the cessation of the growth of the kinematically hot, old, metal-poor thick disk and the kinematically colder, younger, and more metal-rich thin disk. Although perhaps coincidental, the quenching episode could also be related to the formation of the bar. To explore the relation between bar formation and star formation in gas rich galaxies quantitatively, we simulated gas-rich disk isolated galaxies. Our simulations include prescriptions for star formation, stellar feedback, and for regulating the multi-phase interstellar medium. We find that the action of stellar bar efficiently quenches star formation, reducing the star-formation rate by a factor of 10 in less than 1 Gyr. Analytical and self-consistent galaxy simulations with bars suggest that the action of the stellar bar increases the gas random motions within the co-rotation radius of the bar. Indeed, we detect an increase in the gas velocity dispersion up to 20-35 km/s at the end of the bar formation phase. The star-formation efficiency decreases rapidly, and in all of our models, the bar quenches the star formation in the galaxy. The star-formation efficiency is much lower in simulated barred compared to unbarred galaxies and more rapid bar formation implies more rapid quenching., Comment: 17 pages, 17 figure, accepted for publication in A&A

We performed long-slit spectral observations of two SB-type galaxies: NGC 5347, UGC 1344. They were previously suspected as the galaxies with unusually low mass-to-light ratios (on the ground of mass estimates from the Hi linewidths), which are in conflict with their observed colours. The observations were conducted at the Russian 6-m telescope. The aim of the study was to clarify the kinematics and structure, as well as the properties of stellar populations of the galaxies. The results of observations disproved the peculiarly low mass-to-light ratios of both galaxies. The most probable reasons of underestimation of their masses are discussed. We tried to reproduce the main observed features of kinematical profiles of the galaxies in the N-body simulations of barred galaxies. We found that both galaxies possess central components of different structure. Indeed, the age and velocity dispersion of stellar population in NGC 5347 are low in its innermost part in comparison to that of the bulge or a bar, which agrees with the presence of nuclear kinematically decoupled disc. It probably was formed due to the bar which supplied the inner region with gas. The kinematical profiles of the second galaxy UGC 1344 give evidences in favour of the central peanut-shaped bulge. In spite of the different luminosities of the two galaxies, they possess nearly equal (close to solar) central stellar abundance and the flattening of the stellar metallicity gradient in the bar regions. However, in the less luminous NGC 5347 the mean stellar age is younger than that in UGC 1344., Comment: 14 pages, 11 figures, 3 tables, accepted for publication in MNRAS

We present evolutionary calculations for the size and aromatization degree distributions of interstellar dust grains, driven by their destruction by radiation, collisions with gas particles, and shattering due to grain-grain collisions. Based on these calculations we model dust emission spectra. The initial grain size distribution play the important role in the evolution of an ensemble of dust particles. Radiation in the considered intensity range mostly aromatizes grains. The smallest grains are mainly destructed via sputtering by collisions with gas particles. There are no grains smaller than 20~\AA\ in the medium at relative gas-dust velocities more than 50~km/s, which is typical for shocks in supernova remnants. The IR emission spectrum changes significantly due to the dust evolution depending on the adopted grain properties, in particular, on the energy of the C--C bonds ($E_0$). Aromatic bands in the near-IR (2--15~$\mu$m) are absent, if $E_0$ is low, even when the medium properties correspond to the average interstellar medium in our Galaxy. As in reality these bands are observed, high $E_0$ values are more preferable. We consider dependence of the emission intensity ratios for various photometric bands on the medium properties. The aromatization degree of small grains shows up most strongly in the $I_{3.4}/I_{11.3}$ intensity ratio while the fraction of aromatic grains in the total dust mass influences the $I_{3.3}/I_{70+160}$ ratio., Comment: Published in Astronomy Reports

A model of key processes influencing the evolution of a hydrocarbon grain of an arbitrary size under astrophysical conditions corresponding to ionized hydrogen regions (HII regions) and supernova remnants is presented. The considered processes include aromatization and photodestruction, sputtering by electrons and ions, and shattering due to collisions between grains. The model can be used to simulate the grain size distribution and the aromatization degree during the evolution of HII regions and supernova remnants for a specified radiation field, relative velocity of gas and dust, etc. The contribution of various processes to the evolution of hydrocarbon dust grains for parameters typical for the interstellar medium of our Galaxy is presented. Small grains (less than 50 carbon atoms) should be fully aromatized in the general interstellar medium. If larger grains initially have an aliphatic structure, it is preserved to a substantial extent. Variations in the size distribution of the grains due to their mutual collisions depend appreciably on the adopted initial size distribution. For the MRN initial distribution a significant redistribution of grain sizes is obtained, which increases the mass fraction of smaller grains. Characteristic for an initial distribution from the work of Jones et al. (2013), with high initial fraction of small grains, is a general decrease in the number of grains of all sizes., Comment: Published in Astronomy Reports

By means of 3D hydrodynamical simulations, in a separate paper we have discussed the properties of non-axisymmetric density wave trains in the outermost regions of galaxy disks, based on the picture that self-excited global spiral modes in the bright optical stellar disk are accompanied by low-amplitude short trailing wave signals outside corotation; in the gas, such wave trains can penetrate through the outer Lindblad resonance and propagate outwards, forming prominent spiral patterns. In this paper we present the synthetic 21~cm velocity maps expected from simulated models of the outer gaseous disk, focusing on the case when the disk is dominated by a two-armed spiral pattern, but considering also other more complex situations. We discuss some aspects of the spiral pattern in the gaseous periphery of galaxy disks noted in our simulations that might be interesting to compare with specific observed cases., Comment: 14 pages, 7 figures, accepted for publication in Journal of Plasma Physics

We investigate physical properties of molecular clouds in disc galaxies with different morphology: a galaxy without prominent structure, a spiral barred galaxy and a galaxy with flocculent structure. Our $N$-body/hydrodynamical simulations take into account non-equilibrium H$_2$ and CO chemical kinetics, self-gravity, star formation and feedback processes. For the simulated galaxies the scaling relations of giant molecular clouds or so called Larson's relations are studied for two types of a cloud definition (or extraction methods): the first one is based on total column density position-position (PP) datasets and the second one is indicated by the CO~(1-0) line emission used position-position-velocity (PPV) data. We find that the cloud populations obtained by using both cloud extraction methods generally have similar physical parameters. Except that for the CO data the mass spectrum of clouds has a tail with low-massive objects $M\sim 10^3-10^4$~\Msun. Varying column density threshold the power-law indices in the scaling relations are significantly changed. In contrast, the relations are invariant to CO brightness temperature threshold. Finally, we find that the mass spectra of clouds for the PPV data are almost insensitive to the galactic morphology, whereas the spectra for the PP data demonstrate significant variations., Comment: Accepted for publication in MNRAS, 16 pages, 9 figures

Most of giant molecular clouds (GMCs) in M 33 are connected with spiral-like gaseous arms (filaments) with the exception of the inner 2 kpc region where the link between the arms and GMCs disappears (see Tosaki et al. 2011). We check whether it may be caused by the dynamic friction retarding the clouds. Using semi-analytical model for this galaxy we calculate the dynamics of GMCs of different masses situated at different initial galactocentric distances in the disk plane. We demonstrate that the dynamical friction may really change the orbits of GMCs in the central 2 kpc-size region. However in this case the typical lifetimes of GMCs should be close to or greater than $10^8$~yr, which is larger than the usually accepted values., Comment: 5 pages, 5 figures, Accepted for publication in MNRAS

Deep HI observations of the outer parts of disc galaxies demonstrate the frequent presence of extended, well-developed spiral arms far beyond the optical radius. To understand the nature and the origin of such outer spiral structure, we investigate the propagation in the outer gaseous disc of large-scale spiral waves excited in the bright optical disc. Using hydrodynamical simulations, we show that non-axisymmetric density waves, penetrating in the gas through the outer Lindblad resonance, can exhibit relatively regular spiral structures outside the bright optical stellar disc. For low-amplitude structures, the results of numerical simulations match the predictions of a simple WKB linear theory. The amplitude of spiral structure increases rapidly with radius. Beyond $\approx 2$ optical radii, spirals become nonlinear (the linear theory becomes quantitatively and qualitatively inadequate) and unstable to Kelvin-Helmholtz instability. In numerical simulations, in models for which gas is available very far out, spiral arms can extend out to 25 disc scale-lengths. A comparison between the properties of the models we have investigated and the observed properties of individual galaxies may shed light into the problem of the amount and distribution of dark matter in the outer halo., Comment: 12 pages, 13 figures (accepted for publication in MNRAS)

We present a numerical code for multi-component simulation of the galactic evolution. Our code includes the following parts: $N$-body is used to evolve dark matter, stellar dynamics and dust grains, gas dynamics is based on TVD-MUSCL scheme with the extra modules for thermal processes, star formation, magnetic fields, chemical kinetics and multi-species advection. We describe our code in brief, but we give more details for the magneto-gas dynamics. We present several tests for our code and show that our code have passed the tests with a reasonable accuracy. Our code is parallelized using the MPI library. We apply our code to study the large scale dynamics of galactic discs., Comment: 13 pages, 11 figures, 25th IUPAP Conference on Computational Physics (CCP2013) 20-24 August 2013, Moscow, Russia

The polar ring galaxies (PRGs) represent an interesting type of peculiar systems in which the outer matter is rotating in the plane which is roughly perpendicular to the disk of the main galaxy. Despite the long lasting study of the PRGs there is a lack of the detailed enough observational data, there are still many open questions. Among the most interesting issues there are the estimate of the flattening of the dark matter halos in these systems and the verification of the assumption that the most massive polar structures were formed by the accretion of the matter from the intergalactic filaments. The new catalog recently collected by our team using the SDSS images increased by several times the number of known PRGs. The current paper gives the overview of our latest results on the study of morphological and photometric structure of the PRGs. Using the stellar and ionized gas kinematics data based on spectral observations with the Russian 6-m telescope we estimate the dark matter halo shape in the individual galaxies., Comment: 8 pages, to appear in Baltic Astronomy, based on the talk presented on "Modern stellar astronomy-2014" (Rostov-on-Don)

With the aim to determine the spatial distribution of the dark matter halo, we investigate two polar ring galaxies NGC 4262 and SPRC-7. For both galaxies the stellar kinematics data for the central galaxy were obtained from optical spectroscopy at the 6-m telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences. The information about polar gaseous components was taken from the optical 3D-spectroscopic observations of ionized gas (for SPRC-7) and H {\sc i} radio observations (for NGC 4262). SPRC-7 is the system with a relative angle $\delta=73^{\circ}$ towards the central galaxy and a quite massive stellar-gaseous polar component. Meanwhile NGC 4262 is the classic polar case with $\delta=88^{\circ}$ where the polar ring mainly consists of neutral gas with a negligible stellar contribution to the mass. We are hence dealing with two different systems and the results are quite diverse too. The observed properties of both galaxies were compared with the results of self consistent simulations of velocity fields of the polar component along with the rotation curve of the central lenticular galaxy. For SPRC-7 we have found a slightly flattened halo towards the polar plane with the axis ratio $c/a \simeq 1.7 \pm 0.2$ for the isothermal halo model and $c/a \simeq 1.5 \pm 0.2$ for the NFW model. The case of NGC 4262 is more unusual, the shape of the dark matter distribution varies strongly with radius. Namely, the dark matter halo is fattened in the vicinity of the galactic disc ($c/a \approx 0.4 \pm 0.1$), however it is prolate far beyond the central galaxy ($c/a \approx 1.7$ for the isothermal halo and $c/a \approx 2.3$ for NFW)., Comment: 15 pages, 11 figures, Accepted for publication in MNRAS

We develop a 3D dynamical model of dust outflows from galactic discs. The outflows are initiated by multiple SN explosions in a magnetized interstellar medium (ISM) with a gravitationally stratified density distribution. Dust grains are treated as particles in cells interacting collisionally with gas, and forced by stellar radiation of the disc and Lorenz force. We show that magnetic field plays a crucial role in accelerating the charged dust grains and expelling them out of the disc: in 10--20~Myr they can be elevated at distances up to 10~kpc above the galactic plane. The dust-to-gas ratio in the outflowing medium varies in the range $5 \cdot 10^{-4} - 5 \cdot 10^{-2}$ along the vertical stream. Overall the dust mass loss rate depends on the parameters of ISM and may reach up to $3\times 10^{-2}$~\Msun~yr$^{-1}$, Comment: 5 pages, 3 figures. to appear in Proceedings of the conference "The Life Cycle of Dust in the Universe", November 2013, Taipei, Taiwan

With the aim to determine the spatial distribution of dark matter, we investigate the polar ring galaxy NGC 4262. We used the stellar kinematics data for the central galaxy obtained from optical spectroscopy together with information about the kinematics of the neutral hydrogen for polar component. It was shown that NGC 4262 is the classic polar ring galaxy case with the relative angle of $88^{\circ}$ between components. From simulations of the central galaxy and ring kinematics we found that the shape of the dark matter distribution varies strongly with the radius. Namely, the dark matter halo is flattened towards the galactic disk plane $c/a=0.4$, however it is prolate to the orthogonal (polar) plane far beyond the central galaxy $c/a=1.7$. Also, the simulations of the ring evolution let us to confirm the stability of the ring and the formation of quasi-spiral structures within it., Comment: 4 pages, 2 figures, proceedings of "Multi-spin galaxies", Napoli, October 2013

Using $N$-body simulations ($N\sim 10^6 - 10^7$), we examine how a non-axisymmetric dark halo affects the dynamical evolution of the structure in collisionless (stellar) discs. We demonstrate how the model parameters such as mass of the halo, initial conditions in the disc and the halo axes ratio affect morphology and kinematics of the stellar discs. We show that a non-axisymmetric halo can generate a large-scale spiral density pattern in the embedded stellar disc. The pattern is observed in the disc for many periods of its revolution, even if the disc is gravitationally over-stable. The growth of the spiral arms is not accompanied by significant dynamical heating of the disc, irrelevant to its initial parameters. We also investigate transformation of the dark halo's shape driven by the long-lived spiral pattern in the disc . We show that the analysis of the velocity field in the stellar disc and in the spiral pattern gives us a possibility to figure out the spatial orientation of the triaxial-shaped dark halo and to measure the triaxiality., Comment: 11 pages, 13 figures, accepted to MNRAS

Spectroscopic observations at the Russian 6-m telescope are used to study the two polar ring galaxies (PRGs) from the catalogue by Moiseev et al.: SPRC-7 and SPRC-260. We have analyzed the kinematics of the stellar component of the central galaxies as well as the ionized gas kinematics in the external ring structures. The disc-halo decomposition of rotation curves in two perpendicular directions are considered. The observed 2D velocity fields are compared with the model predictions for different dark halo shapes. Based on these data, we constrain that for potential of DM halo semiaxis ratios is $s=0.8$, $q=1$ for SPRC-7 and $s=0.95$, $q=1.1$ for SPRC-260. Using 3D hydrodynamic simulations we also study the dynamics and evolution of the polar component in the potential of the galactic disc and dark halo for these two galaxies. We show that the polar component is dynamically quasi-stable on the scale of $\sim10$ dynamical times (about a few Gyr). This is demonstrate the possibility for the growth of a spiral structure, which then steadily transforms to a lopsided gaseous system in the polar pane., Comment: 4 pages, 4 figures, to pe published in EWASS-SPS4, Roma, P. Di Matteo and C. Jog, eds

Using observational data, we build numerical N-body, hydrodynamical and combined equilibrium models for the spiral galaxy NGC 5247. The models turn out to be unstable towards spiral structure formation. We simulate scenarios of spiral structure formation for different sets of equilibrium rotation curves, radial velocity dispersion profiles and disk thickness and demonstrate that in all cases a simulated spiral pattern qualitatively agrees with the observed morphology of NGC 5247. We also demonstrate that an admixture of a gaseous component with mass of about a few percent of the total mass of the disk increases a lifetime of a spiral pattern by approximately 30%. The simulated spiral pattern in this case lasts for about 3 Gyr from the beginning of the growth of perturbations., Comment: 14 pages, 21 figures, accepted for publication in MNRAS

Using 3D hydrodynamic calculations we simulate formation of molecular clouds in the Galaxy. The simulations take into account molecular hydrogen chemical kinetics, cooling and heating processes. Comprehensive gravitational potential accounts for contributions from the stellar bulge, two and four armed spiral structure, stellar disk, dark halo and takes into account self-gravitation of the gaseous component. Gas clouds in our model form in the spiral arms due to shear and wiggle instabilities and turn into molecular clouds after $t\simgt 100$ Myr. At the times $t\sim 100 - 300$ Myr the clouds form hierarchical structures and agglomerations with the sizes of 100 pc and greater. We analyze physical properties of the simulated clouds and find that synthetic statistical distributions like mass spectrum, "mass-size" relation and velocity dispersion are close to those observed in the Galaxy. The synthetic $l-v$ (galactic longitude - radial velocity) diagram of the simulated molecular gas distribution resembles observed one and displays a structure with appearance similar to Molecular Ring of the Galaxy. Existence of this structure in our modelling can be explained by superposition of emission from the galactic bar and the spiral arms at $\sim$3-4 kpc., Comment: 10 pages, 8 figures

In this paper we examine the effect of the formation and evolution of the disk galaxy on the distribution of dark halo matter. We have made simulations of isolated dark matter (DM) halo and two component (DM + baryons). N-body technique was used for stellar and DM particles and TVD MUSCL scheme for gas-dynamic simulations. The simulations include the processes of star formation, stellar feedback, heating and cooling of the interstellar medium. The results of numerical experiments with high spatial resolution let us to conclude in two main findings. First, accounting of star formation and supernova feedback resolves the so-called problem of cusp in distribution of dark matter predicted by cosmological simulations. Second, the interaction of dark matter with dynamic substructures of stellar and gaseous galactic disk (e.g., spiral waves, bar) has an impact on the shape of the dark halo. In particular, the in-plane distribution of dark matter is more symmetric in runs, where the baryonic component was taken into account., Comment: 7 pages, 6 figures

The dynamics of a galactic disk in a non-axisymmetric (triaxial) dark halo is studied in detail using high-resolution, numerical, hydrodynamical models. A long-lived, two-armed spiral pattern is generated for a wide range of parameters. The spiral structure is global, and the number of turns can be two or three, depending on the model parameters. The morphology and kinematics of the spiral pattern are studied as functions of the halo and disk parameters. The spiral structure rotates slowly, and its angular velocity varies quasi-periodically. Models with differing relative halo masses, halo semi-axis ratios, distributions of matter in the disk, Mach numbers in the gaseous component, and angular rotational velocities of their halos are considered., Comment: 22 pages, 11 figures

The results of numerical simulations of a gaseous disk in the potential of a stellar spiral density wave are presented. The conditions under which straightened spiral arm segments (rows) form in the gas component are studied. These features of the spiral structure were identified in a series of works by A.D. Chernin with coauthors. Gas-dynamic simulations have been performed for a wide range of model parameters: the pitch angle of the spiral pattern, the amplitude of the stellar spiral density wave, the disk rotation speed, and the temperature of the gas component. The results of 2D- and 3D-disk simulations are compared. The rows in the numerical simulations are shown to be an essentially nonstationary phenomenon. A statistical analysis of the distribution of geometric parameters for spiral patterns with rows in the observed galaxies and the constructed hydrodynamic models shows good agreement. In particular, the numerical simulations and observations of galaxies give $<{\alpha}>\simeq 120^\circ$ for the average angles between straight segments., Comment: 22 pages, 10 figures