X-structures are often observed in galaxies hosting the so-called B/PS (boxy/peanuts) bulges and are visible from the edge-on view. They are the most notable features of B/PS bulges and appear as four rays protruding from the disk of the host galaxy and distinguishable against the B/PS bulge background. In some works their origin is thought to be connected with the so-called banana-shaped orbits with a vertical resonance 2:1. A star in such an orbit performs two oscillations in the vertical direction per one revolution in the bar frame. Several recent studies that analyzed ensembles of orbits arising in different $N$-body models do not confirm the dominance of the resonant 2:1 orbits in X-structures. In our work we analyze two $N$-body models and show how the X-structure in our models is gradually assembled from the center to the periphery from orbits with less than 2:1 frequency ratio. The most number of such orbits is enclosed in a 'farfalle'-shape (Italian pasta) form and turns out to be non-periodic. We conclude that the X-structure is only the envelope of regions of high density caused by the crossing or folding of different types of orbits at their highest points, and does not have a "backbone" similar to that of the in-plane bar. Comparing the orbital structure of two different numerical models, we show that the dominance of one or another family of orbits with a certain ratio of the vertical oscillations frequency to the in-plane frequency depends on the parameters of the underlying galaxy and ultimately determines the morphology of the X-structure and the opening angle of its rays.

By means of self-consistent numerical simulations we investigated the dynamical impact of classical bulges on the growth of the secondary buckling of a bar. Overall we considered 14 models with different disc and bulge parameters. We obtained that a bulge with a quite modest mass $B/D=0.1$ leads to completely symmetrical evolution of the bar almost independently of the initial stellar disc parameters and even can damp the first bending. At the same time, the bars in all our bulgeless models suffer from the short primary and prolonged secondary buckling. Given the smallness of the mass suppressing secondary buckling, we conclude that a classical bulge along with the gas central concentration may be the main culprits for the rarity of bars with ongoing buckling in the local Universe.

Recent observational studies of X-shaped structures revealed that values of their opening angles lie in a narrow range: from 20 deg to 43 deg with smaller X-shaped structures appearing to follow a characteristic opening angle ~ 27 deg - 31 deg. We use self-consistent numerical simulations to uncover what parameters of host galaxies govern the opening angle spread. We constructed a series of equilibrium models of galaxies with high spatial resolution, varying the dark halo and bulge contribution in the overall gravitational potential, the initial disc thickness of models and the Toomre parameter Q and followed their evolution for almost 8 Gyrs. Each model demonstrated the formation of clear X-structures with different flatness. We have obtained that opening angles lie in the range from 25 deg to 42 deg throughout the entire evolution. These values are roughly consistent with observational data. The greatest variation in the opening angles is obtained by varying the mass of the dark halo. The initial thickness of the disc and the Toomre parameter are responsible for smaller variations of the angle and shorter X-structures. An increase of both parameters changes the morphology of in-plane bars and X-structures. In some cases we observed even double X-structures. The main effect of the bulge is to prevent buckling at intermediate and late stages of the disc evolution. Comparison of models with different halo masses indicates that the smallest observable values of opening angles can be associated with the presence of a heavy dark halo (more than 3 masses of the disc within the optical radius).

We applied a criterion of gravitational instability, valid for two-component and infinitesimally thin discs, to observational data along the major axis for 7 spiral galaxies of early types. Unlike most papers, the dispersion equation corresponding to the criterion was solved directly without using any approximation. The velocity dispersion of stars in the radial direction $\sigma_R$ was limited by the range of possible values instead of a fixed value. For all galaxies, the outer regions of the disc were analyzed up to $R < 130^{"}$. The maximal and sub-maximal disc models were used to translate surface brightness into surface density. The largest destabilizing disturbance stars can exert on a gaseous disc was estimated. It was shown that the two-component criterion differs a little from the one-fluid criterion for galaxies with a large surface gas density, but it allows to explain large-scale star formation in those regions where the gaseous disc is stable. In the galaxy NGC 1167 star formation is entirely driven by the self-gravity of the stars. A comparison is made with the conventional approximations which also include the thickness effect and with models for different sound speed $c_\mathrm{g}$. It is shown that values of the effective Toomre parameter correspond to the instability criterion of a two-component disc $Q_\mathrm{eff}<1.5-2.5$. This result is consistent with previous theoretical and observational studies.

We present a detailed photometric study of a sample of 22 edge-on galaxies with clearly visible X-shaped structures. We propose a novel method to derive geometrical parameters of these features, along with the parameters of their host galaxies based on the multi-component photometric decomposition of galactic images. To include the X-shaped structure into our photometric model, we use the IMFIT package, in which we implement a new component describing the X-shaped structure. This method is applied for a sample of galaxies with available SDSS and Spitzer IRAC 3.6 $\mu$m observations. In order to explain our results, we perform realistic $N$-body simulations of a Milky Way-type galaxy and compare the observed and the model X-shaped structures. Our main conclusions are as follows: (1) galaxies with strong X-shaped structures reside in approximately the same local environments as field galaxies; (2) the characteristic size of the X-shaped structures is about 2/3 of the bar size; (3) there is a correlation between the X-shaped structure size and its observed flatness: the larger structures are more flattened; (4) our $N$-body simulations qualitatively confirm the observational results and support the bar-driven scenario for the X-shaped structure formation.

We present a modification of the method for reconstructing the stellar velocity ellipsoid (SVE) in disc galaxies. Our version does not need any parametrization of the velocity dispersion profiles and uses only one assumption that the ratio $\sigma_z/\sigma_R$ remains constant along the profile or along several pieces of the profile. The method was tested on two galaxies from the sample of other authors and for the first time was applied to three lenticular galaxies NGC~1167, NGC~3245 and NGC~4150 as well as to one Sab galaxy NGC~338. We found that for galaxies with a high inclination ($i >55-60^\circ$) it is difficult or rather impossible to extract the information about SVE while for galaxies at an intermediate inclination the procedure of extracting is successful. For NGC~1167 we managed to reconstruct SVE, provided that the value of $\sigma_z/\sigma_R$ is piecewise constant. We found $\sigma_z/\sigma_R=0.7$ for the inner parts of the disc and $\sigma_z/\sigma_R=0.3$ for the outskirts. We also obtained a rigid constrain on the value of the radial velocity dispersion $\sigma_R$ for highly inclined galaxies and tested the result using the asymmetric drift equation, provided that the gas rotation curve is available.

We study the properties of galaxies with very thin discs using a sample of 85 objects whose stellar disc radial-to-vertical scale ratio determined from photometric decomposition, exceeds nine. We present evidences of similarities between the very thin disc galaxies (VTD galaxies) and low surface brightness (LSB) disc galaxies, and conclude that both small and giant LSB galaxies may reveal themselves as VTD, edge-on galaxies. Our VTD galaxies are mostly bulgeless, and those with large radial scale length tend to have redder colors. We performed spectral observations of 22 VTD galaxies with the Dual Imaging Spectrograph on the 3.5m telescope at the Apache Point Observatory. The spectra with good resolution (R ~ 5000) allow us to determine the distance and the ionized gas rotation curve maximum for the galaxies. Our VTD galaxies have low dust content, in contrast to regular disc galaxies. Apparently, VTD galaxies reside in specific cosmological low-density environments and tend to have less connection with filaments. Comparing a toy model that assumes marginally low star formation in galactic discs with obtained gas kinematics data, we conclude that there is a threshold central surface density of about 88 Mo/pc**2, which we observe in the case of very thin, rotationally supported galactic discs.

Based on SDSS data, we have selected a sample of nine edge-on spiral galaxies with bulges whose major axes show a high inclination to the disk plane. Such objects are called polar-bulge galaxies. They are similar in their morphology to polar-ring galaxies, but the central objects in them have small size and low luminosity. We have performed a photometric analysis of the galaxies in the g and r bands and determined the main characteristics of their bulges and disks. We show that the disks of such galaxies are typical for the disks of spiral galaxies of late morphological types. The integrated characteristics of their bulges are similar to the parameters of normal bulges. The stellar disks of polar-bulge galaxies often show large-scale warps, which can be explained by their interaction with neighboring galaxies or external accretion from outside.

We analyze the dependence of the stellar disc flatness on the galaxy morphological type using 2D decomposition of galaxies from the reliable subsample of the Edge-on Galaxies in SDSS (EGIS) catalogue. Combining these data with the retrieved models of the edge-on galaxies from the Two Micron All Sky Survey (2MASS) and the Spitzer Survey of Stellar Structure in Galaxies (S$^4$G) catalogue, we make the following conclusions: (1) The disc relative thickness $z_0/h$ in the near- and mid-infrared passbands correlates weakly with morphological type and does not correlate with the bulge-to-total luminosity ratio $B/T$ in all studied bands. (2) Applying an 1D photometric profile analysis overestimates the disc thickness in galaxies with large bulges making an illusion of the relationship between the disc flattening and the ratio $B/T$. (3) In our sample the early-type disc galaxies (S0/a) have both flat and "puffed" discs. The early spirals and intermediate-type galaxies have a large scatter of the disc flatness, which can be caused by the presence of a bar: barred galaxies have thicker stellar discs, on average. On the other hand, the late-type spirals are mostly thin galaxies, whereas irregular galaxies have puffed stellar discs.

We present a catalog of true edge-on disk galaxies automatically selected from the Seventh Data Release (DR7) of the Sloan Digital Sky Survey. A visual inspection of the $g$, $r$ and $i$ images of about 15000 galaxies allowed us to split the initial sample of edge-on galaxy candidates into 4768 (31.8% of the initial sample) genuine edge-on galaxies, 8350 (55.7%) non-edge-ons, and 1865 (12.5%) edge-on galaxies not suitable for simple automatic analysis because these objects show signs of interaction, warps, or nearby bright stars project on it. We added more candidate galaxies from RFGC, EFIGI, RC3, and Galaxy Zoo catalogs found in the SDSS footprints. Our final sample consists of 5747 genuine edge-on galaxies. We estimate the structural parameters of the stellar disks (the stellar disk thickness, radial scale length, and central surface brightness) in the galaxies by analyzing photometric profiles in each of the g, r, and i images. We also perform simplified 3-D modeling of the light distribution in the stellar disks of edge-on galaxies from our sample. Our large sample is intended to be used for studying scaling relations in the stellar disks and bulges and for estimating parameters of the thick disks in different types of galaxies via the image stacking. In this paper we present the sample selection procedure and general description of the sample.

We analyzed several basic correlations between structural parameters of galaxies. The data were taken from various samples in different passbands which are available in the literature. We discuss disc scaling relations as well as some debatable issues concerning the so-called Photometric Plane for bulges and elliptical galaxies in different forms and various versions of the famous Kormendy relation. We show that some of the correlations under discussion are artificial (self-correlations), while others truly reveal some new essential details of the structural properties of galaxies. Our main results are as follows: (1) At present, we can not conclude that faint stellar discs are, on average, more thin than discs in high surface brightness galaxies. The ``central surface brightness -- thickness'' correlation appears only as a consequence of the transparent exponential disc model to describe real galaxy discs. (2) The Photometric Plane appears to have no independent physical sense. Various forms of this plane are merely sophisticated versions of the Kormendy relation or of the self-relation involving the central surface brightness of a bulge/elliptical galaxy and the Sersic index n. (3) The Kormendy relation is a physical correlation presumably reflecting the difference in the origin of bright and faint ellipticals and bulges. We present arguments that involve creating artificial samples to prove our main idea.

We demonstrate that in N-body simulations of isolated disc galaxies there is numerical vertical heating which slowly increases the vertical velocity dispersion and the disc thickness. Even for models with over a million particles in a disc, this heating can be significant. Such an effect is just the same as in numerical experiments by Sellwood (2013). We also show that in a stellar disc, outside a boxy/peanut bulge, if it presents, the saturation level of the bending instability is rather close to the value predicted by the linear theory. We pay attention to the fact that the bending instability develops and decays very fast, so it couldn't play any role in secular vertical heating. However the bending instability defines the minimal value of the ratio between the vertical and radial velocity dispersions $\sigma_z / \sigma_R \approx 0.3$ (so indirectly the minimal thickness) which could have stellar discs in real galaxies. We demonstrate that observations confirm last statement.

The model of an young star with a protoplanetary disk and a low-mass companion, which is moving on a circular orbit, inclined to the disk plane, is considered. The hydrodynamic models of such a system were calculated by SPH method. It was shown the distortions in the disk, caused by the orbital motion of the companion, lead to the strong dependence of illumination conditions of the disk on the azimuth (because of extinction between the star and the disk surface) and, therefore, it leads to the appearance of large-scale asymmetry in images of disks. The calculations showed the dependence of the illumination on the azimuth was stronger in the central part of the disk than on the periphery. Bright and dark domains are located not symmetry with respect to the line of nodes. The sizes and locations of the domains are depended on the model parameters as well as on the phase of the orbital period. The bright and dark domains do not follow the companion, but they make small amplitude oscillations with respect to some direction. The properties of the model, which were written above, open new opportunities of searching low-mass companions in the vicinity of young stars. The stars with protoplanetary disks, which are observed face-on or under the small inclination angle i, are the best ones for this purpose.

The cyclic activity model of a young star with the low-mass secondary component (q = M2/M1 <= 0.1) accreting a matter from circumbinary disk is considered. It is assumed that the orbit is circular and the disk and orbital planes are non-coplanar. Sets of hydrodynamics models of such a system have been calculated by the SPH method and then the variations of the circumstellar extinction and phase light curves were determined. The calculations showed that depending on the model parameters and orientation of the system in regards to an observer the different in shape and amplitude light curves can be observed. An important property of the considered models is also the dependence of the mass accretion rate onto the components on the phase of the orbital period. The results of the calculation can be used for analysis of the cyclic activity of UX Ori stars and young stars with the long-lasting eclipses.

We present the results of our analysis of the geometrical tidal tail characteristics for nearby and distant interacting galaxies. The sample includes more than two hundred nearby galaxies and about seven hundred distant ones. The distant galaxies have been selected in several deep fields of the Hubble Space Telescope (HDF-N, HDF-S, HUDF, GOODS, GEMS) and they are at mean redshift z=0.65. We analyze the distributions of lengths and thicknesses for the tidal structures and show that the tails in distant galaxies are shorter than those in nearby ones. This effect can be partly attributed to observational selection effects, but, on the other hand, it may result from the general evolution of the sizes of spiral galaxies wih z. The location of interacting galaxies on the galaxy luminosity ($L$) -- tidal tail length ($l$) plane are shown to be explained by a simple geometrical model, with the upper envelope of the observed distribution being $l \propto \sqrt{L}$. We have solved the problem on the relationship between the observed distribution of tail flattening and the tail length in angular measure by assuming the tidal tails to be arcs of circumferences visible at arbitrary angles to the line of sight. We conclude that the angular length of the tidal tails visually distinguished in nearby and distant galaxies, on average, exceeds 180 degrees.

We consider a model for the cyclic brightness variations of a young star with a low-mass companion that accretes matter from the remnants of a protostellar cloud. At small inclinations of the binary orbit to the line of sight, the streams of matter and the density waves excited in the circumbinary disk can screen the primary component of the binary from the observer. To study these phenomena, we have computed grids of hydrodynamic models for binary systems by the SPH method based on which we have constructed the phase light curves as a function of the rotation angle of the apsidal line relative to the observer. The model parameters were varied within the following ranges: the component mass ratio q = 0.01-0.1 and the eccentricity e = 0-0.5. We adopted optical grain characteristics typical of circumstellar dust. Our computations have shown that the brightness oscillations with orbital phase can have a complex structure. The amplitudes and shapes of the light curves depend strongly on the inclination of the binary orbit and its orientation relative to the observer and on the accretion rate. The results of our computations are used to analyze the cyclic activity of UX Ori stars.

We consider a model for the cyclic activity of young binary stars that accrete matter from the remnants of a protostellar cloud. If the orbit of such a binary system is inclined at a small angle to the line of sight, then the streams of matter and the density waves excited in the circumbinary disk can screen the primary component of the binary from the observer. To study these phenomena by the SPH (smoothed particle hydrodynamics) method, we have computed grids of hydrodynamic models for binary systems based on which we have constructed the light curves as a function of the orbital phase.The main emphasis is on investigating the properties of the brightness oscillations. Therefore, the model parameters were varied within the following ranges: the component mass ratio q = M2 : M1 = 0.2 - 0.5 and the eccentricity e = 0 - 0.7. The parameter that defined the binary viscosity was also varied. We adopted optical grain characteristics typical of circumstellar dust. Our computations have shown that bimodal oscillations are excited in binaries with eccentric orbits, provided that the binary components do not differ too much in mass. In this case, the ratios of the periods and amplitudes of the bimodal oscillations and their shape depend strongly on the inclination of the binary plane and its orientation relative to the observer. Our analysis shows that the computed light curves can be used in interpreting the cyclic activity of UX Ori stars.

NGC7217 is an unbarred early-type spiral galaxy having a multi-segment exponential light profile and a system of starforming rings of the unknown origin; it also possesses a circumnuclear gaseous polar disc. We analysed new long slit spectroscopic data for NGC7217 and derived the radial distributions of its stellar population parameters and stellar and gaseous kinematics up to the radius of r~100 arcsec (~8 kpc). We performed the dynamical analysis of the galaxy by recovering its velocity ellipsoid at different radii, and estimated the scaleheights of its two exponential discs. The inner exponential stellar disc of NGC7217 appears to be thin and harbours intermediate age stars (t(SSP) ~ 5 Gyr). The outer stellar disc seen between the radii of 4 and 7 kpc is very thick (z0 = 1...3 kpc), metal-poor, [Fe/H]<-0.4 dex, and has predominantly young stars, t(SSP) = 2 Gyr. The remnants of minor mergers of gas-rich satellites with an early-type giant disc galaxy available in the GalMer database well resemble different structural components of NGC7217, suggesting two minor merger events in the past responsible for the formation of the inner polar gaseous disc and large outer starforming ring. Another possibility to form the outer ring is the re-accretion of the tidal streams created by the first minor merger.

We present a sample of edge-on spiral galaxies both of early and late types.The sample consists of 175 galaxies in the Ks-filter, 169 galaxies in the H-filter and 165 galaxies in the J-filter. Bulge and disc decompositions of each galaxy image, taken from the Two Micron All Sky Survey (2MASS), were performed. We discuss several scaling relations for bulges and discs which indicate a tight link between their formation and evolution. We show that galaxies with bulges fitted by the Sersic index n<2 (pseudobulges) have quite different distributions of their structural parameters than galaxies with n>=2 bulges (classical bulges). First of all, the distribution of the apparent bulge axis ratio q_b for the subsample with n<2 can be attributed to triaxial, nearly prolate bulges, while n>=2 bulges seem to be oblate spheroids with moderate flattening. Secondly, the Photometric Plane of the sample bulges is not flat and has a prominent curvature towards small values of n. Thirdly, despite of the existence of a clear relation between the flattening of stellar discs h/z_0 and the relative mass of a spherical component, the distributions over both parameters are quite different for galaxies possesing bulges and pseudobulges.

Malin 1 is a unique, extraordinarily large low surface brightness galaxy. The structure and the origins of the galaxy are poorly understood. The reason for such a situation is an absence of detailed observational data, especially, of high-resolution kinematics. In this Letter we study the stellar kinematics of the inner part (r < 15 kpc) of Malin 1. We present spectroscopic arguments in favour of a small galaxy - Malin 1B - being a companion probably interacting with the main galaxy - Malin 1. This object is clearly seen in many published images of Malin 1 but is not mentioned in any astronomical databases. Malin 1B is located at the projected distance of 14 kpc from the Malin 1's nucleus and has small - 65$\pm$16 km/s - relative velocity, which we determined for the first time. We suggest that ongoing interaction with Malin 1B can explain main morphological features of the Malin 1's central region - two-armed spiral structure, a bar, and an external one-armed spiral pattern. We also investigated the large scale environment of Malin 1 and postulate that the galaxy SDSS J123708.91+142253.2 might be responsible for the formation of extended low-surface brightness envelope by means of head-on collision with Malin 1 (in the framework of collision scenario proposed by Mapelli et al. 2008). To test the collisional origins of Malin 1 global structure, more observational data and new numerical models are needed.

A sample of edge-on spiral galaxies aimed at a study of the main structural and photometric parameters of edge-on galaxies both of early and late types is presented. The data were taken from the 2MASS in the J, H and K_s filters. The sample consists of 175 galaxies in the K_s-filter, 169 galaxies in the H-filter and 165 galaxies in the J-filter. We present bulge and disc decompositions of each galaxy image. All galaxies have been modelled with a Sersic bulge and exponential disc with the BUDDA v2.1 package. The main conclusions of our general statistical analysis of the sample are: (1) The distribution of the apparent bulge axis ratio q_b for the subsample with n < 2 can be attributed to triaxial, nearly prolate bulges that are seen from different projections, while n > 2 bulges seem to be oblate spheroids with moderate flattening. (2) For the sample galaxies, the effective radius of the bulge r_{e,b}, the disc scalelength h and the disc scaleheight z_0 are well correlated. However, there is a clear trend for the ratio r_{e,b}/h to increase with n. (3) There is a hint that the fundamental planes of discs, which links only disc parameters and the maximum rotational velocity of gas, are different for galaxies with different bulges. (4) The investigation of the Photometric Plane of sample bulges shows that the plane is not flat and has a prominent curvature towards small values of n. For bulges this fact was not noticed earlier. (5) The clear relation between the flattening of stellar discs h/z_0 and the relative mass of a spherical component, including a dark halo, is confirmed not for bulgeless galaxies but for galaxies with massive bulges. (Abridged)

We present a new method for constructing equilibrium phase models for stellar systems, which we call the iterative method. It relies on constrained, or guided evolution, so that the equilibrium solution has a number of desired parameters and/or constraints. This method is very powerful, to a large extent due to its simplicity. It can be used for mass distributions with an arbitrary geometry and a large variety of kinematical constraints. We present several examples illustrating it. Applications of this method include the creation of initial conditions for N-body simulations and the modelling of galaxies from their photometric and kinematic observations.

Hydrodynamic models of a young binary system accreting matter from the remnants of a protostellar cloud have been calculated by the SPH method. It is shown that periodic variations in column density in projection onto the primary component take place at low inclinations of the binary plane to the line of sight. They can result in periodic extinction variations. Three periodic components can exist in general case. The first component has a period equal to the orbital one and is attributable to the streams of matter penetrating into the inner regions of the binary. The second component has a period that is a factor of 5-8 longer than the orbital one and is related to the density waves generated in a circumbinary (CB) disk. The third, longest period is attributable to the precession of the inner asymmetric region of CB disk. The relationship between the amplitudes of these cycles depends on the model parameters as well as on the inclination and orientation of the binary in space. We show that at a dust-to-gas ratio of 1:100 and and a mass extinction coefficient of 250 cm$^2$ g$^{-1}$, the amplitude of the brightness variations of the primary component in the V-band can reach $1^m$ at a mass accretion rate onto the binary components of $10^{-8} M_{\odot}$ yr$^{-1}$ and a $10^{\rm o}$ inclination of the binary plane to the line of sight. We discuss possible applications of the model to pre-main-sequence stars.

An iterative approach is used to construct spherically symmetric equilibrium models with an anisotropic velocity distribution. The potentialities of the method have been tested on models with known distribution functions, the Osipkov-Merritt models. It is shown that models that differ significantly from the Osipkov-Merritt models can be constructed. An N-body model of a dark halo with a density distribution that approximates the results of cosmological simulations (the Navarro--Frenk--White model) has been constructed. The anisotropy profile has been taken to be similar to that yielded by cosmological simulations. The constructed models can serve as direct input data for investigating the dynamics and stability of such systems in N-body simulations.

One widely used technique for the construction of equilibrium models of stellar disks is based on the Jeans equations and the moments of velocity distribution functions derived using these equations. Stellar disks constructed using this technique are shown to be "not entirely" in equilibrium. Our attempt to abandon the epicyclic approximation and the approximation of infinite isothermal layers, which are commonly adopted in this technique, failed to improve the situation substantially. We conclude that the main drawback of techniques based on the Jeans equations is that the system of equations employed is not closed, and therefore requires adopting an essentially ad hoc additional closure condition. A new iterative approach to constructing equilibrium N-body models with a given density distribution is proposed. The main idea behind this approach is that a model is first constructed using some approximation method, and is then allowed to adjust to an equilibrium state with the specified density and the required parameters of the velocity distribution remaining fixed in the process. This iterative approach was used to construct isotropic, spherically symmetric models and models of stellar disks embedded in an external potential. The numerical models constructed prove to be close to equilibrium. It is shown that the commonly adopted assumption that the profile of the radial velocity dispersion is exponential may be wrong. The technique proposed can be applied to a wide range of problems involving the construction of models of stellar systems with various geometries.

We analyze the relationship between the mass of a spherical component and the minimum possible thickness of stable stellar disks. This relationship for real galaxies allows the lower limit on the dark halo mass to be estimated (the thinner the stable stellar disk is, the more massive the dark halo must be). In our analysis, we use both theoretical relations and numerical N-body simulations of the dynamical evolution of thin disks in the presence of spherical components with different density profiles and different masses. We conclude that the theoretical relationship between the thickness of disk galaxies and the mass of their spherical components is a lower envelope for the model data points. We recommend using this theoretical relationship to estimate the lower limit for the dark halo mass in galaxies. The estimate obtained turns out to be weak. Even for the thinnest galaxies, the dark halo mass within four exponential disk scale lengths must be more than one stellar disk mass.

A criterion for the choice of optimal softening length $\epsilon$ and time-step $dt$ for $N$-body simulations of a collisionless stellar system is analyzed. Plummer and Hernquist spheres are used as models to follow how changes in various parameters of an initially equilibrium stable model depend on $\epsilon$ and $dt$. These dependences are used to derive a criterion for choosing $\epsilon$ and $dt$. The resulting criterion is compared to Merritt's criterion for choosing the softening length, which is based on minimizing the mean irregular force acting on a particle with unit mass. Our criterion for choosing $\epsilon$ and $dt$ indicate that $\epsilon$ must be a factor of 1.5-2 smaller than the mean distance between particles in the densest regions to be resolved. The time-step must always be adjusted to the chosen $\epsilon$ (the particle must, on average, travel a distance smaller than $0.5\epsilon$ during one time-step). An algorithm for solving N-body problems with adaptive variations of the softening length is discussed in connection with the task of choosing $\epsilon$, but is found not to be promising.

The saturation conditions for bending modes in inhomogeneous thin stellar disks that follow from an analysis of the dispersion relation are compared with those derived from $N$-body simulations. In the central regions of inhomogeneous disks, the reserve of disk strength against the growth of bending instability is smaller than that for a homogeneous layer. The spheroidal component (a dark halo, a bulge) is shown to have a stabilizing effect. The latter turns out to depend not only on the total mass of the spherical component, but also on the degree of mass concentration toward the center. We conclude that the presence of a compact (not necessarily massive) bulge in spiral galaxies may prove to be enough to suppress the bending perturbations that increase the disk thickness. This conclusion is corroborated by our $N$-body simulations in which we simulated the evolution of almost equilibrium, but unstable finite-thickness disks in the presence of spheroidal components. The final disk thickness at the same total mass of the spherical component (dark halo + bulge) has been found to be much smaller than that in the simulations where a concentrated bulge is present.

We present results of numerical modeling of the cyclic brightness modulation in the young binary systems with the eccentric orbits and low-mass secondary components. Brightness variations of the primary is due to the periodical extinction variations on the line-of-sight caused by the disk wind of the secondary and a common envelope it produces. A matter distribution in the envelope has been calculated in the ballistic approach. Calculations showed that for the young binaries with the elliptic orbits parameters of the photometric minima (their depth, duration and the shape of light curves) depend not only on the disk wind parameters and an inclination of the binary orbit to the line-of-sight but also on the longitude of the periastron. A modulation of the scattered radiation of the common envelope with a phase of the orbital period has been investigated in the single scattering approach. It is shown that an amplitude of the modulation is maximal when the system is seen edge-on and has also a non-zero value in the binaries observed pole-on. Possible applications of the theory to the young stellar objects are discussed. In particular, an attention is payed to a resemblance of the light curves in some models with light curves of the objects suspected as candidates to FUORs.

We investigate the nonlinear growth stages of bending instability in stellar disks with exponential radial density profiles.We found that the unstable modes are global (the wavelengths are larger than the disk scale lengths) and that the instability saturation level is much higher than that following from a linear criterion. The instability saturation time scales are of the order of one billion years or more. For this reason, the bending instability can play an important role in the secular heating of a stellar disk in the $z$ direction. In an extensive series of numerical $N$-body simulations with a high spatial resolution, we were able to scan in detail the space of key parameters (the initial disk thickness $z_0$, the Toomre parameter $Q$, and the ratio of dark halo mass to disk mass $M_{\rm h} / M_{\rm d}$). We revealed three distinct mechanisms of disk heating in the $z$ direction: bending instability of the entire disk, bending instability of the bar, and heating on vertical inhomogeneities in the distribution of stellar matter.

We present an analysis of the location of interacting E/S0 galaxies on the Fundamental Plane (FP). Using the NEMO package, we performed N-body simulations of close encounters and mergers of two self-gravitating spherical galaxies. Two models for encounting galaxies - Plummer's and Hernquist's - were used. The number of particles ranged from 20000 to 50000 per galaxy. The changes of central density, half-mass radius and central velocity dispersion were analysed. It was found that close encounters between galaxies alter noticeably the above parameters within a very short time interval (10^7 - 10^8 years) just before the final merger. The amplitudes of parameter changes strongly depend on the initial mass concentration of a model. In some experiments we considered an encounter of two spherical galaxies with dark matter components. The effect of parameter changes was less pronounced than for the experiments without a dark halo. The results of the simulations were used to discuss the FP for interacting early-type galaxies.

We review recent results on the tidal structures of spiral galaxies. Topics included are general characteristics of tails; kinematics of tidal structures and dark haloes of host galaxies; frequency of tidal distortions at z~1.

Numerical simulations of the disruption of a dwarf companion moving in the polar plane of a massive galaxy are presented. The constructed model is compared with observational data on the recently discovered low-surface-brightness stellar ring around galaxy NGC 5907. Constraints on the ring lifetime (<=1.5 gyr after the first approach of the galaxies), on the structure of companion -- the ring presecursor, and on mass of dark halo of the main galaxy in whose gravitational field the companion moves are provided. The dark halo mass within 50 kpc of NGC 5907 center cannot exceed 3 or 4 "visible" masses.

We present an analysis of interacting E/S0 galaxies location on the Fundamental Plane. Using the NEMO package, we performed N-body simulations of close encounters and mergers between two spherical galaxies. We followed how structural and dynamical parameters (central density, half-mass radius and velocity dispersion)of galaxies are changed during the encounter. We analysed the dependence of these changes on initial mass concentration and presence of dark halo. The results of our simulations are used to discuss the Fundamental Plane for interacting early-type galaxies.

We present an analysis of structural features of all known galaxies with optical polar rings. We find a clear dichotomy for objects of this peculiar class. Bulge-dominated S0 galaxies possess only short narrow rings, while disk-dominated objects always have wide extended polar rings. We try by gas dynamical simulations to explain such a segregation by dependence of the ring-forming process on different galaxy potentials. It is found that the total mass captured into the ring during an encounter of a host-ring system with a gas-rich spiral galaxy of comparable mass exceeds 10^9 M_{\odot} (or about 10% of all gas in the donor galaxy), which is of the order of that found by observation. The process of gas to gather into a steady-state ring takes approximately (7 - 9) x 10^8 years. This time is somewhat shorter for rings forming around bulge-dominated galaxies. We also present observational arguments for S0 galaxies with extended rings to be similar to late-type spirals by their photometric properties, while numerical modelling of the extended ring formation suggests that these galaxies must possess massive dark halos as well. In this case, the sizes of the modelled rings turn out large enough (up to 30 kpc in diameter), and the time scale for ring formation is prolonged up to several Gyrs.