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.

Using observational data on the kinematical properties and density distributions of the subsystems of the Milky Way galaxy, we construct a set of multi-component equilibrium models of its disk. The dynamics of the disk is studied numerically using collisionless-gaseous numerical simulations. After approximately one Gyr, a prominent central bar is formed with a semi-axis of about three kiloparsecs. Outside the central regions, a multi-armed spiral pattern develops, which can be characterized by the superposition of m=2, 3 and 4-armed spiral patterns. The spiral structure and the bar exist for at least 3 Gyr in our simulations. The presence of the bar in the disk of the MilkyWay galaxy imposes rather strict limitations on the density distributions in the subsystems of the MilkyWay galaxy.We find that bar does not form if the radial scale length of the surface density distribution of the disk is more than 2.6 kpc. Analogously, the formation of bar is suppressed in the Milky Way disk in models with a massive and compact stellar bulge. If future GAIA observations confirm the existence of the three-four armed spiral pattern in the disk of the MilkyWay, this will prove the long-term existence of spiral patterns in galactic disks.

Parallel implementation features of self-gravitating gas dynamics modeling on multiple GPUs are considered applying the GPU-Direct technology. The parallel algorithm for solving of the self-gravitating gas dynamics problem based on hybrid OpenMP-CUDA parallel programming model has been described in detail. The gas-dynamic forces are calculated by the modified SPH-method (Smoothed Particle Hydrodynamics) while the N-body problem gravitational interaction is obtained by the direct method (so-called Particle-Particle algorithm). The key factor in the SPH-method performance is creation of the neighbor lists of the particles which contribute into the gas-dynamic forces calculation. Our implementation is based on hierarchical grid sorting method using a cascading algorithm for parallel computations of partial sums at CUDA block. The parallelization efficiency of the algorithm for various GPUs of the Nvidia Tesla line (K20, K40, K80) is studied in the framework of galactic' gaseous halos collisions models by the SPH-method.

Galaxies with 'rows' in Vorontsov-Velyaminov's terminology stand out among the variety of spiral galactic patterns. A characteristic feature of such objects is the sequence of straight-line segments that forms the spiral arm. In 2001 A. Chernin and co-authors published a catalog of such galaxies which includes 204 objects from the Palomar Atlas. In this paper, we supplement the catalog with 276 objects based on an analysis of all the galaxies from the New General Catalogue and Index Catalogue. The total number of NGC and IC galaxies with rows is 406, including the objects of Chernin et al. (2001). The use of more recent galaxy images allowed us to detect more 'rows' on average, compared with the catalog of Chernin et al. When comparing the principal galaxy properties we found no significant differences between galaxies with rows and all S-type NGC/IC galaxies.We discuss twomechanisms for the formation of polygonal structures based on numerical gas-dynamic and collisionless N-body calculations, which demonstrate that a spiral pattern with rows is a transient stage in the evolution of galaxies and a system with a powerful spiral structure can pass through this stage. The hypothesis of A. Chernin et al. (2001) that the occurrence frequency of interacting galaxies is twice higher among galaxies with rows is not confirmed for the combined set of 480 galaxies. The presence of a central stellar bar appears to be a favorable factor for the formation of a system of 'rows'.

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.

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

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

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.

We present the results of spectroscopic observations of three S0-Sa galaxies: NGC338, NGC 3245, and NGC 5440 at the SAO RAS 6-m BTA telescope. The radial distributions of the line-of-sight velocities and radial velocity dispersions of stars and ionized gas were obtained, and rotation curves of galaxies were computed. We construct the numerical dynamic N-body galaxy models with N > 10^6 point masses. The models include three components: a "live" bulge, a collisionless disk, dynamically evolving to the marginally stable state, and a pseudo-isothermal dark halo. The estimates of radial velocities and velocity dispersions of stars obtained from observations are compared with model estimates, projected onto the line of sight. We show that the disks of NGC 5440 and the outer regions of NGC 338 are dynamically overheated. Taking into account the previously obtained observations, we conclude that the dynamic heating of the disk is present in a large number of early-type disk galaxies, and it seems to ensue from the external effects. The estimates of the disk mass and relative mass of the dark halo are given for seven galaxies, observed at the BTA.

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.

We use the marginal stability condition for galactic disks and the stellar velocity dispersion data published by different authors to place upper limits on the disk local surface density at two radial scalelengths $R=2h$. Extrapolating these estimates, we constrain the total mass of the disks and compare these estimates to those based on the photometry and color of stellar populations. The comparison reveals that the stellar disks of most of spiral galaxies in our sample cannot be substantially overheated and are therefore unlikely to have experienced a significant merging event in their history. The same conclusion applies to some, but not all of the S0 galaxies we consider. However, a substantial part of the early type galaxies do show the stellar velocity dispersion well in excess of the gravitational stability threshold suggesting a major merger event in the past. We find dynamically overheated disks among both seemingly isolated galaxies and those forming pairs. The ratio of the marginal stability disk mass estimate to the total galaxy mass within four radial scalelengths remains within a range of 0.4---0.8. We see no evidence for a noticeable running of this ratio with either the morphological type or color index.

The non-linear dynamics of bending instability and vertical structure of a galactic stellar disc embedded into a spherical halo are studied with N-body numerical modelling. Development of the bending instability in stellar galactic disc is considered as the main factor that increases the disc thickness. Correlation between the disc vertical scale height and the halo-to-disc mass ratio is predicted from the simulations. The method of assessment of the spherical-to-disc mass ratio for edge-on spiral galaxies with small bulges is considered. Modelling of eight edge-on galaxies: NGC 891, NGC 4738, NGC 5170, UGC 6080, UGC 7321, UGC 8286, UGC 9422 and UGC 9556 is performed. Parameters of stellar discs, dark haloes and bulges are estimated. The lower limit of the dark-to-luminous mass ratio in our galaxies is of the order of 1 within the limits of their stellar discs. The dark haloes dominate by mass in the galaxies with very thin stellar discs (NGC 5170, UGC 7321 and UGC 8286).

We present results of numerical modeling made for the galactic stellar and stellar-gas disk embedded in the spherical halo and bulge. The stellar disk is simulated by N-body system, the equations of hydrodynamics are solved by TVD-method. We used TREEcode-algorithm for calculation of a self-gravity in stellar and gaseous components. The possibility of bars birth in a hot stellar disk because of gravitational instability of a cold gas component is investigated. The conditions of occurrence lopsided-galaxies from a axisymmetric disk as a result of gravitational instability are explored. The self-consistent models of double bars are constructed and the dynamical stability of these structures is discussed.

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.

The static pressure of the hot gas that fills clusters and groups of galaxies can affect significantly the volume density and thickness of the gas disks in galaxies. In combination with the dynamic pressure, the static pressure allows several observed peculiarities of spiral galaxies surrounded by a hot medium to be explained.

Spectroscopic observations of three lenticular (S0) galaxies (NGC 1167, NGC 4150, and NGC 6340) and one SBa galaxy (NGC 2273) have been taken with the 6-m telescope of the Special AstrophysicalObservatory of the Russian Academy of Sciences aimed to study the structure and kinematic properties of early-type disk galaxies. The radial profiles of the stellar radial velocities and the velocity dispersion are measured. N-body simulations are used to construct dynamical models of galaxies containing a stellar disk, bulge, and halo. The masses of individual components are estimated formaximum-mass disk models. A comparison of models with estimated rotational velocities and the stellar velocity dispersion suggests that the stellar disks in lenticular galaxies are "overheated"; i.e., there is a significant excess velocity dispersion over the minimum level required to maintain the stability of the disk. This supports the hypothesis that the stellar disks of S0 galaxies were subject to strong gravitational perturbations. The relative thickness of the stellar disks in the S0 galaxies considered substantially exceed the typical disk thickness of spiral galaxies.

We present N-body simulations of unstable spiral modes in a dynamically cool collisionless disc. We show that spiral modes grow in a thin collisionless disk in accordance with the analytical perturbation theory. We use the particle-mesh code SUPERBOX with nested grids to follow the evolution of unstable spirals that emerge from an unstable equilibrium state. We use a large number of particles (up to 40 million particles) and high-resolution spatial grids in our simulations (128^3 cells). These allow us to trace the dynamics of the unstable spiral modes until their wave amplitudes are saturated due to nonlinear effects. In general, the results of our simulations are in agreement with the analytical predictions. The growth rate and the pattern speed of the most unstable bar-mode measured in N-body simulations agree with the linear analysis. However the parameters of secondary unstable modes are in lesser agreement because of the still limited resolution of our simulations.

We propose that the lower bound of the stellar radial velocity dispersion c_r of an equilibrium stellar disk is determined by the gravitational stability condition. We compared the estimates of stellar velocity dispersion at radii r>(1.5 - 2)R_0 (where R_0 is the photometric radial scalelength of a disk), found in the literature, with the minimal values of c_r necessary for the disk to be in a stable state, using the results of numerical simulations of 3D collisional disks. This approach enables to estimate an upper limit of the local surface density and (if R_0 is known) a total masses of a disk and a dark halo. We argue that the old stellar disks of spiral galaxies with active star formation usually have the velocity dispersions which are close to the expected marginal values. A rough values of disk-to-total mass ratios (within the fixed radius) are found for about twenty spiral galaxies. Unlike spirals, the disks of the ``red'' Sa--S0 galaxies are evidently ``overheated'': their radial dispersion of velocities at r ~ 2 R_0 exceeds significantly the marginal values for gravitational stability.

We present results of numerical N-body simulations of a galactic stellar disk embedded into a spherical dark halo. The non-linear dynamics of bending instabilities developing in the disk is studied. The bending modes, axisymmetric and not, are considered as main factors increasing the disk thickness. The relation between the disk vertical scale height and the halo+bulge-to-disk mass ratio is inferred. The method of estimation of the spherical-to-disk mass ratio for edge-on spiral galaxies based on this relation is studied and applied to constrain the spherical subsystem mass and the mass of dark halos (M_h/M_d) in seven edge-on galaxies. The values of M_h/M_d are of order 1 for our galaxies.

N-body dynamical simulations are used to analyze the conditions for the gravitational stability of a three-dimensional stellar disk in the gravitational field of two rigid spherical components--a bulge and a halo whose central concentrations and relative masses vary over wide ranges. The number of point masses N in the simulations varies from 40 to 500 thousands and the evolution of the simulated models is followed over 10--20 rotation periods of the outer edge of the disk. The initially unstable disks are heated and, as a rule, reach a quasi-stationary equilibrium with a steady-state radial-velocity dispersion $c_r$ over five to eight periods of rotation. The radial behavior of the Toomre stability parameter $Q_T (r)$ for the final state of the disk is estimated. Numerical models are used to analyze the dependence of the gravitational stability of the disk on the relative masses of the spherical components, disk thickness, degree of differential rotation, and initial state of the disk. Formal application of existing, analytical, local criteria for marginal stability of the disk can lead to errors in radial velocity dispersion $c_r$ of more than a factor of 1.5. It is suggested that the approximate constancy of $Q_T \simeq 1.2 -- 1.5$ for $r\simeq (1\div 2)\times L$ (where L is the radial scale of disk surface density), valid for a wide range of models, can be used to estimate upper limits for the mass and density of a disk based on the observed distributions of the rotational velocity of the gaseous component and of the stellar velocity dispersion.

We present results of numerical modeling made for the galactic stellar disk embedded in the spherical halo. The non-linear dynamics of bending instabilities developed in the disk is studied. The axisymmetrical bending mode is considered as a main factor increasing the disk thickness. The model dependence between the disk thickness and the halo to disk mass ratio $M_h/M_d$ is inferred. The modeling results enable us to estimate the dark halo mass for several edge-on galaxies. The obtained values of $M_h/M_d$ for the galaxies considered are of order 1 - 2.

In this paper we present calculations of period distribution for old accreting isolated neutron stars (INSs). At the age about a few billions years low velocity INSs come to the stage of accretion. At that stage their period evolution is governed by magnetic breaking and accreted angular momentum. Due to turbulence of the interstellar medium (ISM) accreted momentum can both accelerate and decelerate rotation of an INS and spin evolution has chaotic character. Calculations show that for constant magnetic fields INSs have relatively long spin periods, $\ge 10^4$--$10^5$ s, depending on parameters of INSs and an ISM. Due to long periods INSs have high spin up/spin down rates, which should fluctuate on a time scale about a few years.

In this paper we present calculations of period distribution for old accreting isolated neutron stars (INSs). After few billion years of evolution low velocity INSs come to the stage of accretion. At this stage INS's period evolution is governed by magnetic braking and angular momentum accreted. Since the interstellar medium is turbulized accreted momentum can either accelerate or decelerate spin of an INS, therefore the evolution of period has chaotic character. Our calculations show that in the case of constant magnetic field accreting INSs have relatively long spin periods (some hours and more, depending on INS's spatial velocity, its magnetic field and density of the surrounding medium). Such periods are much longer than the values measured by {\emph ROSAT} for 3 radio-silent isolated neutron stars. Due to long periods INSs should have high spin up/down rates, $\dot p$, which should fluctuate on a time scale of $\sim 1$ yr.

Different ways of the determination of masses of galactic disks, based on the kinematic data, are briefly discussed. The analysis of the rotation curves which reach maximum inside of a disk, and N-body modeling, reproducing the rotation curves and stellar velocity dispersion of real galaxies, enable to conclude that the mass of a disk is usually significantly less than the total mass of spherical components (bulge + dark halo) inside of optical borders of a galaxy, although the exceptions also exist.

Dynamics of linear perturbations in a differentially rotating accretion disk with non-homogeneous vertical structure is investigated. It has been found that turbulent viscosity results in instability of both pinching oscillations, and bending modes. Not only the low-frequency fundamental modes, but also the high-frequency reflective harmonics appear to be unstable. The question of the limits of applicability of the thin disk model (MTD) is also investigated. The insignificant distinctions in the dispersion properties of MTD and three-dimensional model appear for wave numbers k<(1-3)/h (h is the half-thickness of a disk). In the long-wavelenght limit, the relative difference between eigenfrequencies of the unstable acoustic mode in the 3D-model and the MTD is smaller than 5%.

The model of a geometrically thin gaseous disk in the external gravitational potential is considered. The dinamics of small nonaxisymmetric perturbations in the plane of the accretion disk with dissipative effects is investigated. It is showed, that conditions of development and parameters of unstable oscillation modes in the opticaly thick accretion disk are strongly depended on the models of viscosity and opacity.