We present a set of self-consistent chemo-dynamical simulations of MW-type galaxies formation to study the origin of the bimodality of $\alpha$-elements in stellar populations. We explore how the bimodality is related to the geometrically and kinematically defined stellar discs, gas accretion and radial migration. We find that the two $\alpha$-sequences are formed in quite different physical environments. The high-$\alpha$ sequence is formed early from a burst of star formation (SF) in a turbulent, compact gaseous disc which forms a thick disc. The low-$\alpha$ stellar populations is the result of quiescent SF supported by the slow accretion of enriched gas onto a radially extended thin disc. Stellar feedback-driven outflows during the formation of the thick disc are responsible for the enrichment of the surrounding gaseous halo, which subsequently feeds the disc on a longer time-scale. During the thin disc phase, chemical evolution reaches an equilibrium metallicity and abundance, where the stars pile-up. This equilibrium metallicity decreases towards the outer disc, generating the ridgeline that forms the low-$\alpha$ sequence. We identify a second mechanism capable of creating a low-$\alpha$ sequence in one of our simulations. Rapid shutdown of the SF, provoked by the feedback at the end of the thick disc phase, suppresses the chemical enrichment of the halo gas, which, once accreted onto the star-forming disc, dilutes the ISM at the beginning of the thin disc formation. Both mechanisms can operate in a galaxy, but the former is expected to occur when SF efficiency ceases to be dominated by the formation of the thick disc, while the latter can occur in the inner regions. Being the result of the presence of low and high gas density environments, the bimodality is independent of any particular merger history, suggesting that it could be much more widespread than has been claimed.

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.

Boxy, peanut- or X-shaped "bulges" are observed in a large fraction of barred galaxies viewed in, or close to, edge-on projection, as well as in the Milky Way. They are the product of dynamical instabilities occurring in stellar bars, which cause the latter to buckle and thicken vertically. Recent studies have found nearby galaxies that harbour two such features arising at different radial scales, in a nested configuration. In this paper we explore the formation of such double peanuts, using a collisionless N-body simulation of a pure disc evolving in isolation within a live dark matter halo, which we analyse in a completely analogous way to observations of real galaxies. In the simulation we find a stable double configuration consisting of two X/peanut structures associated to the same galactic bar - rotating with the same pattern speed - but with different morphology, formation time, and evolution. The inner, conventional peanut-shaped structure forms early via the buckling of the bar, and experiences little evolution once it stabilises. This feature is consistent in terms of size, strength and morphology, with peanut structures observed in nearby galaxies. The outer structure, however, displays a strong X, or "bow-tie", morphology. It forms just after the inner peanut, and gradually extends in time (within 1 to 1.5 Gyr) to almost the end of the bar, a radial scale where ansae occur. We conclude that, although both structures form, and are dynamically coupled to, the same bar, they are supported by inherently different mechanisms.

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.

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.

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.

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]

Using a single N-body simulation ($N=0.14\times 10^9$) we explore the formation, evolution and spatial variation of the phase-space spirals similar to those recently discovered by Antoja et al. in the Milky Way disk, with Gaia DR2. For the first time in the literature, we use a self-consistent N-body simulation of an isolated Milky Way-type galaxy to show that the phase-space spirals develop naturally from vertical oscillations driven by the buckling of the stellar bar. We claim that the physical mechanism standing behind the observed incomplete phase-space mixing process can be internal and not necessarily due to the perturbation induced by a massive satellite. In our model, the bending oscillations propagate outwards and produce axisymmetric variations of the mean vertical coordinate and of the vertical velocity component. As a consequence, the phase-space wrapping results in the formation of patterns with various morphology across the disk, depending on the bar orientation, distance to the galactic center and time elapsed since the bar buckling. Once bending waves appear, they are supported for a long time via disk self-gravity. The underlying physical mechanism implies the link between in-plane and vertical motion that leads directly to phase-space structures whose amplitude and shape are in remarkable agreement with those of the phase-space spirals observed in the Milky Way disk. In our isolated galaxy simulation, phase-space spirals are still distinguishable, at the solar neighbourhood, 3 Gyr after the buckling phase. The long-lived character of the phase-space spirals generated by the bar buckling instability cast doubts on the timing argument used so far to get back at the time of the onset of the perturbation: phase-space spirals may have been caused by perturbations originated several Gyrs ago, and not as recent as suggested so far.

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]

Globular clusters (GCs), the oldest stellar systems observed in the Milky Way, have for long been considered single stellar populations. As such, they provided an ideal laboratory to understand stellar dynamics and primordial star formation processes. However, during the last two decades, observations unveiled their real, complex nature. Beside their pristine stars, GCs host one or more helium enriched and possibly younger stellar populations whose formation mechanism is still unknown. Even more puzzling is the existence of GCs showing star by star iron spreads. Using detailed N-body simulations we explore the hypothesis that these anomalies in metallicity could be the result of mutual stripping and mergers between a primordial population of disc GCs. In the first paper of this series we proved, both with analytical arguments and short-term N-body simulations, that disc GCs have larger fly-by and close encounter rates with respect to halo clusters. These interactions lead to mass exchange and even mergers that form new GCs, possibly showing metallicity spreads. Here, by means of long-term direct N-body simulations, we provide predictions on the dynamical properties of GCs that underwent these processes. The comparison of our predictions with available and future observational data could provide insights on the origin of GCs and on the Milky Way build-up history as a whole.

The formation scenario for giant low surface brightness (gLSB) galaxies with discs as large as 100 kpc still remains unclear. These stellar systems are rare and very hard to observe, therefore a detailed insight on every additional object helps to understand their nature. Here we present a detailed observational study of the gLSB UGC 1922 performed using deep optical imaging and spectroscopic observations combined with archival ultraviolet data. We derived spatially resolved properties of stellar and ionized gas kinematics and characteristics of stellar populations and interstellar medium. We reveal the presence of a kinematically decoupled central component, which counter rotates with respect to the main disc of UGC 1922. The radial metallicity gradient of the ionised gas is in agreement with that found for moderate-size LSB galaxies. At the same time, a slowly rotating and dynamically hot central region of the galaxy hosts a large number of old metal-rich stars, which creates an appearance of a giant elliptical galaxy, that grew an enormous star forming disc. We reproduce most of the observed features of UGC 1922 in N-body/hydrodynamical simulations of an in-plane merger of giant Sa and Sd galaxies. We also discuss alternative formation scenarios of this unusual system.

We present the results of a self-consistent $N$-body simulation following the evolution of a primordial population of thick disc globular clusters (GCs). We study how the internal properties of such clusters evolve under the action of mutual interactions, while they orbit a Milky Way-like galaxy. For the first time, through analytical and numerical considerations, we find that physical encounters between disc GCs are a crucial factor that contributed to the shape of the current properties of the Galactic GC system. Close passages or motion on similar orbits may indeed have a significant impact on the internal structure of clusters, producing multiple gravitationally bound sub-populations through the exchange of mass and even mergers. Our model produces two major mergers and a few small mass exchanges between pairs of GCs. Two of our GCs accrete stars from two companions, ending up with three internal sub-populations. We propose these early interactions and mergers between thick disc GCs with slightly different initial chemical compositions as a possible explanation for the presence of the spreads in metallicity observed in some of the massive Milky Ways GCs.

We investigate the nature of the double color-magnitude sequence observed in the Gaia DR2 HR diagram of stars with high transverse velocities. The stars in the reddest-color sequence are likely dominated by the dynamically-hot tail of the thick disk population. Information from Nissen & Schuster (2010) and from the APOGEE survey suggests that stars in the blue-color sequence have elemental abundance patterns that can be explained by this population having a relatively low star-formation efficiency during its formation. In dynamical and orbital spaces, such as the `Toomre diagram', the two sequences show a significant overlap, but with a tendency for stars on the blue-color sequence to dominate regions with no or retrograde rotation and high total orbital energy. In the plane defined by the maximal vertical excursion of the orbits versus their apocenters, stars of both sequences redistribute into discrete wedges. We conclude that stars which are typically assigned to the halo in the solar vicinity are actually both accreted stars lying along the blue sequence in the HR diagram, and the low velocity tail of the old Galactic disk, possibly dynamically heated by past accretion events. Our results imply that a halo population formed in situ and responsible for the early chemical enrichment prior to the formation of the thick disk is yet to be robustly identified, and that what has been defined as the stars of the in situ stellar halo of the Galaxy may be in fact fossil records of its last significant merger.

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.

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.

We show that the bulge and the disk of the Milky Way (MW) at R$\lesssim$7~kpc are well described by a unique chemical evolution and a two-phase star-formation history (SFH). We argue that the populations within this inner disk, not the entire disk, are the same, and that the outer Lindblad resonance (OLR) of the bar plays a key role in explaining this uniformity. In our model of a two-phase star formation history, the metallicity, [$\alpha$/Fe] and [$\alpha$/H] distributions, and age-metallicity relation are all compatible with the observations of both the inner disk and bulge. The dip at [Fe/H]$\sim$0 dex seen in the metallicity distributions of the bulge and inner disk reflects the quenching episode in the SFH of the inner MW at age $\sim$8 Gyr, and the common evolution of the bulge and inner disk stars. We show that at z$\le$1.5, when the MW was starting to quench, transitioning between the end of the $\alpha$-enhanced thick disk formation to the start of the thin disk, and yet was still gas rich, the gas accretion rate could not have been significant. The [$\alpha$/Fe] abundance ratio before and after this quenching phase would be different, which is not observed. The present analysis suggests that Milky Way history, and in particular at the transition from the thick to the thin disk -- the epoch of the quenching -- must have been driven by a decrease of the star formation efficiency. We argue that the decline in the intensity of gas accretion, the formation of the bar, and the quenching of the SFR at the same epoch may be causally connected thus explaining their temporal coincidence (abridged).

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.

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.

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.

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.

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.

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.

Using N-body/gasdynamic simulations of a Milky Way-like galaxy we analyse a Kennicutt-Schmidt relation, $\Sigma_{SFR} \propto \Sigma_{gas}^N$, at different spatial scales. We simulate synthetic observations in CO lines and UV band. We adopt the star formation rate defined in two ways: based on free fall collapse of a molecular cloud - $\Sigma_{SFR, cl}$, and calculated by using a UV flux calibration - $\Sigma_{SFR, UV}$. We study a KS relation for spatially smoothed maps with effective spatial resolution from molecular cloud scales to several hundred parsecs. We find that for spatially and kinematically resolved molecular clouds the $\Sigma_{SFR, cl} \propto \Sigma_{\rm gas}^N$ relation follows the power-law with index $N \approx 1.4$. Using UV flux as SFR calibrator we confirm a systematic offset between the $\Sigma_{\rm UV}$ and $\Sigma_{\rm gas}$ distributions on scales compared to molecular cloud sizes. Degrading resolution of our simulated maps for surface densities of gas and star formation rates we establish that there is no relation $\Sigma_{\rm SFR, UV} - \Sigma_{\rm gas}$ below the resolution $\sim 50$ pc. We find a transition range around scales $\sim 50-120$ pc, where the power-law index $N$ increases from 0 to 1-1.8 and saturates for scales larger $\sim 120$ pc. A value of the index saturated depends on a surface gas density threshold and it becomes steeper for higher $\Sigma_{gas}$ threshold. Averaging over scales with size of $>150$ pc the power-law index $N$ equals 1.3-1.4 for surface gas density threshold $\sim 5 M_\odot$pc$^{-2}$. At scales $>120$ pc surface SFR densities determined by using CO data and UV flux, $\Sigma_{\rm SFR, UV}/\Sigma_{\rm SFR, cl}$, demonstrate a discrepancy about a factor of 3. We argue that this may be originated from overestimating (constant) values of conversion factor, star formation efficiency or UV calibration used in our analysis.

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.

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.

With the help of high-resolution long-slit and integral-field spectroscopy observations, the number of confirmed cases of galaxies with counterrotation is increasing rapidly. The evolution of such counterrotating galaxies remains far from being well understood. In this paper we study the dynamics of counterrotating collisionless stellar disks by means of $N$-body simulations. We show that, in the presence of counterrotation, an otherwise gravitationally stable disk can naturally generate bending waves accompanied by strong disk heating across the disk plane, that is in the vertical direction. Such conclusion is found to hold even for dynamically warm systems with typical values of the initial vertical-to-radial velocity dispersion ratio $\sigma_{\rm z}/\sigma_{\rm R} \approx 0.5$, for which the role of pressure anisotropy should be unimportant. We note that, during evolution, the $\sigma_{\rm z}/\sigma_{\rm R}$ ratio tends to rise up to values close to unity in the case of locally Jeans-stable disks, whereas in disks that are initially Jeans-unstable it may reach even higher values, especially in the innermost regions. This unusual behavior of the $\sigma_{\rm z}/\sigma_{\rm R}$ ratio in galaxies with counterrotation appears not to have been noticed earlier. Our investigations of systems made of two counterrotating components with different mass-ratios suggest that even apparently normal disk galaxies (i.e., with a minor counterrotating component so as to escape detection in current observations) might be subject to significant disk heating especially in the vertical direction.

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.

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.

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.

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.

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.

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.

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).

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}$

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.

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.

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.

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.

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.

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.

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.

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.

The results of numerical simulations of the dynamics of galactic disks, which are submerged into nonaxisymmetric dark massive halo are discussed. Galaxy disks dynamics in the nonaxisymmetric (triaxial) dark halo were investigated in detail by the high resolution numerical hydrodynamical methods (TVD \& SPH) and N-body models. The long-lived two arms spiral structure generates for a wide range of parameters. The spiral structure is global and number of turns can be 2-3 in depends of model parameters. Morphology and kinematics of spiral structures were investigated in depends of the halo and the disk parameters. The spiral structure rotates slowly and the angular velocity varies is quasiperiodic.