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.

We use the astrometric and photometric data from Gaia Data Release 2 and line-of-sight velocities from various other surveys to study the 3d structure and kinematics of the Sagittarius dwarf galaxy. The combination of photometric and astrometric data makes it possible to obtain a very clean separation of Sgr member stars from the Milky Way foreground; our final catalogue contains ~2.6e5 candidate members with magnitudes G<18, more than half of them being red clump stars. We construct and analyze maps of the mean proper motion and its dispersion over the region ~30x12 degrees, which show a number of interesting features. The intrinsic 3d density distribution (orientation, thickness) is strongly constrained by kinematics; we find that the remnant is a prolate structure with the major axis pointing at 45deg from the orbital velocity and extending up to ~5 kpc, where it transitions into the stream. We perform a large suite of N-body simulations of a disrupting Sgr galaxy as it orbits the Milky Way over the past 2.5 Gyr, which are tailored to reproduce the observed properties of the remnant (not the stream). The richness of available constraints means that only a narrow range of parameters produce a final state consistent with observations. The total mass of the remnant is ~4e8 Msun, of which roughly a quarter resides in stars. The galaxy is significantly out of equilibrium, and even its central density is below the limit required to withstand tidal forces. We conclude that the Sgr galaxy is terminally wounded, and will be fully disrupted over the next Gyr.

Supermassive black holes (BHs) and their host galaxies are interlinked by virtue of feedbacks and are thought to be co-eval across the Hubble time. This relation is highlighted by an approximate proportionality between the BH mass $M_\bullet$ and the mass of a stellar bulge $M_\ast$ of the host galaxy. However, a large spread of the ratio $M_\bullet/M_\ast$ and a considerable excess of BH mass at redshifts $z\sim8$, indicate that the coevolution of central massive BHs and stellar populations in host galaxies may have experienced variations in its intensity. These issues require a robust determination of the relevant masses (BH, stars and gas), which is difficult in the case of distant high-redshift galaxies that are unresolved. In this paper, we seek to identify spectral diagnostics that may tell us about the relative masses of the BH, the gas mass and stellar mass. We consider general features of SEDs of galaxies that harbour growing massive BHs, forming stars and interstellar/circumgalactic gas. We focus on observational manifestations of possible predominances or intermittent variations in evolutionary episodes of growing massive BHs and forming stellar populations. We consider simplified scenarios for star formation and massive BHs growth, and simple models for chemical composition of gas, for dust free gas as well as for gas with dust mass fraction of $1/3$ of the metal content. We argue that wideband multi-frequency observations (X-ray to submillimeter) of the composite emission spectra of growing BH, stellar population and nebular emission of interstellar gas are sufficient to infer their masses.

Tidal disruption events occur rarely in any individual galaxy. Over the last decade, however, time-domain surveys have begun to accumulate statistical samples of these flares. What dynamical processes are responsible for feeding stars to supermassive black holes? At what rate are stars tidally disrupted in realistic galactic nuclei? What may we learn about supermassive black holes and broader astrophysical questions by estimating tidal disruption event rates from observational samples of flares? These are the questions we aim to address in this Chapter, which summarizes current theoretical knowledge about rates of stellar tidal disruption, and compares theoretical predictions to the current state of observations.

Agama (Action-based Galaxy Modelling Architecture) is a software library intended for a broad range of tasks within the field of stellar dynamics. As the name suggests, it is centered around the use of action/angle formalism to describe the structure of stellar systems, but this is only one of its many facets. The library contains a powerful framework for dealing with arbitrary density/potential profiles and distribution functions (analytic, extracted from N-body models, or fitted to the data), a vast collection of general-purpose mathematical routines, and covers many aspects of galaxy dynamics up to the very high-level interface for constructing self-consistent galaxy models. This document serves two purposes. First of all, this is a detailed reference for the library itself. Second, it describes various mathematical and numerical methods that could be applicable in a broader context. These include: (1) one- and multidimensional interpolation using cubic and quintic splines; (2) penalized spline fitting of noisy data and penalized spline density estimation; (3) adaptive rejection sampling from multidimensional probability distributions; (4) computation of gravitational potentials using spherical- and azimuthal-harmonic expansions; (5) Staeckel fudge method for computing action/angle variables; (6) a general discussion about good programming practices.

We present Forstand, a new code for constructing dynamical models of galaxies with the Schwarzschild orbit-superposition method. These models are constrained by line-of-sight kinematic observations and applicable to galaxies of all morphological types, including disks and triaxial rotating bars. Our implementation has several novel and improved features, is computationally efficient, and made publicly available. Using mock datasets taken from N-body simulations, we demonstrate that the pattern speed of a bar can be recovered with an accuracy of 10-20%, regardless of orientation, if the 3D shape of the galaxy is known or inferred correctly.

The presence of supermassive black holes (SMBHs) with masses up to $M_\bullet\sim10^9M_\odot$ at redshifts $z\simeq7.5$ suggests that their seeds may have started to grow long before the reionization in ambient medium with pristine chemical composition. During their latest 500Myr episode of growing from $z\geq10$ to $z\sim7$ the black holes shine as luminous as $10^{11}\hbox{--}10^{12}L_\odot$, with a cumulative spectrum consisting of the intrinsic continuum from hot accretion disk, nebular hydrogen and helium spectral lines and free-free continuum from gas of host halos. Here we address the question of whether such a plain spectrum would allow us to trace evolution of these growing SMBHs. In our calculations we assume that host galaxies have stellar populations with masses smaller than the mass of their central black holes -- the so-called obese black hole galaxies. Within this model we show that for a sufficiently high mass of gas in a host galaxy -- not smaller than the mass of a growing black hole, the cumulative spectrum in the far-infrared reveals a sharp transition from a quasi-blackbody Rayleigh-Jeans spectrum of the black hole $\propto\lambda^{-2}$ to a flat free-free nebular continuum $\lambda^{0.118}$ on longer wavelength limit. Once such a transition in the spectrum is resolved, the black hole mass can be inferred as a combination of the observed wavelength at the transition $\lambda_k$ and the corresponding spectral luminosity. Possible observability of this effect in spectra of growing high-$z$ SMBHs and determination of their mass with the upcoming JWST and the planned space project Spektr-M is briefly discussed.

We consider the possibility of the detection of $^3$HeII hyperfine line (rest frequency, $8.67 \, \rm GHz$) emission from ionized zones around accreting black holes (BHs) formed at high redshifts, $z=15\hbox{--}30$. We show that the brightness temperature in 8.67GHz line increases and reaches a peak value after the accretion onto the BH exhausts and HeIII recombines into HeII. This period of brightening last up to 40~million years. We find that during this period the maximum brightness temperature reaches $\simeq 0.2\hbox{--}0.5 \mu$K, depending on the epoch when such a black hole starts growing. The maximum angular size of the region emitting in the hyperfine line is around $0.5'$. The flux from such a region ($\simeq 0.3 \, \rm nJy$) is too small to be detected by SKA1-MID. The RMS of the collective flux from many emitting regions from a volume bounded by the synthesized beam and the band-width of SKA1-MID might reach 100~nJy, which is potentially detectable by SKA1-MID.

We review the Schwarzschild orbit-superposition approach and present a new implementation of this method, which can deal with a large class of systems, including rotating barred disk galaxies. We discuss two conceptuals problems in this field: the intrinsic degeneracy of determining the potential from line-of-sight kinematics, and the non-uniqueness of deprojection and related biases in potential inference, especially acute for triaxial bars. When applied to mock datasets with known 3d shape, our method correctly recovers the pattern speed and other potential parameters. However, more work is needed to systematically address these two problems for real observational datasets.

We review the implications of the Gaia Data Release 2 catalogue for studying the dynamics of Milky Way globular clusters, focusing on two separate topics. The first one is the analysis of the full 6-dimensional phase-space distribution of the entire population of Milky Way globular clusters: their mean proper motions (PM) can be measured with an exquisite precision (down to 0.05 mas/yr, including systematic errors). Using these data, and a suitable ansatz for the steady-state distribution function (DF) of the cluster population, we then determine simultaneously the best-fit parameters of this DF and the total Milky Way potential. We also discuss possible correlated structures in the space of integrals of motion. The second topic addresses the internal dynamics of a few dozen of the closest and richest globular clusters, again using the Gaia PM to measure the velocity dispersion and internal rotation, with a proper treatment of spatially correlated systematic errors. Clear rotation signatures are detected in 10 clusters, and a few more show weaker signatures at a level >=0.05$ mas/yr. PM dispersion profiles can be reliably measured down to 0.1 mas/yr, and agree well with the line-of-sight velocity dispersion profiles from the literature.

We use stellar proper motions (PM) from Gaia Data Release 2 for studying the internal kinematics of Milky Way globular clusters. In addition to statistical measurement errors, there are significant spatially correlated systematic errors, which cannot be ignored when studying the internal kinematics. We develop a mathematically consistent procedure for incorporating the spatial correlations in any model-fitting approach, and use it to determine rotation and velocity dispersion profiles of a few dozen clusters. We confirm detection of rotation in the sky plane for ~10 clusters reported in previous studies, and discover a few more clusters with rotation amplitudes exceeding ~0.05 mas/yr. However, in more than half of these cases the significance of this rotation signature is rather low when taking into account the systematic errors. We find that the PM dispersion is not sensitive to systematic errors in PM, however, it is quite sensitive to the selection criteria on the input sample, most importantly, in crowded central regions. When using the cleanest possible samples, PM dispersion can be reliably measured down to 0.1 mas/yr for ~60 clusters.

The Gaia Sausage is the major accretion event that built the stellar halo of the Milky Way galaxy. Here, we provide dynamical and chemical evidence for a second substantial accretion episode, distinct from the Gaia Sausage. The Sequoia Event provided the bulk of the high energy retrograde stars in the stellar halo, as well as the recently discovered globular cluster FSR 1758. There are up to 6 further globular clusters, including $\omega$~Centauri, as well as many of the retrograde substructures in Myeong et al. (2018), associated with the progenitor dwarf galaxy, named the Sequoia. The stellar mass in the Sequoia galaxy is $\sim 5 \times 10^{7} M_\odot$, whilst the total mass is $\sim 10^{10} M_\odot$, as judged from abundance matching or from the total sum of the globular cluster mass. Although clearly less massive than the Sausage, the Sequoia has a distinct chemo-dynamical signature. The strongly retrograde Sequoia stars have a typical eccentricity of $\sim0.6$, whereas the Sausage stars have no clear net rotation and move on predominantly radial orbits. On average, the Sequoia stars have lower metallicity by $\sim 0.3$ dex and higher abundance ratios as compared to the Sausage. We conjecture that the Sausage and the Sequoia galaxies may have been associated and accreted at a comparable epoch.

Binary supermassive black holes (SMBH) are expected to form naturally during galaxy mergers. After the dynamical friction phase, when the two SMBHs become gravitationally bound to each other, and a brief stage of initial rapid hardening, the orbit gradually continues to shrink due to three-body interactions with stars that enter the loss cone of the binary. Using the stellar-dynamical Monte Carlo code Raga, we explore the co-evolution of the binary SMBH and the nuclear star cluster in this slow stage, for various combinations of parameters (geometry of the star cluster, primary/secondary SMBH mass, initial eccentricity, inclusion of stellar captures/tidal disruptions). We compare the rates of stellar captures in galactic nuclei containing a binary to those of galaxies with a single SMBH. At early times, the rates are substantially higher in the case of a binary SMBH, but subsequently they drop to lower levels. Only in triaxial systems both the binary hardening rates and the capture rates remain sufficiently high during the entire evolution. We find that the hardening rate is not influenced by star captures, nor does it depend on eccentricity; however, it is higher when the difference between the black hole masses is greater. We confirm that the eccentricity of the binary tends to grow, which may significantly shorten the coalescence time due to earlier onset of gravitational-wave emission. We also explore the properties of the orbits entering the loss cone, and demonstrate that it remains partially full throughout the evolution in the triaxial case, but significantly depleted in the axisymmetric case. Finally, we study the distribution of ejected hypervelocity stars and the corresponding mass deficits in the central parts of the galaxies hosting a binary, and argue that the missing mass is difficult to quantify observationally.

We use Gaia Data Release 2 to determine the mean proper motions for 150 Milky Way globular clusters (almost the entire known population), with a typical uncertainty of 0.05 mas/yr limited mainly by systematic errors. Combining them with distance and line-of-sight velocity measurements from the literature, we analyze the distribution of globular clusters in the 6d phase space, using both position/velocity and action/angle coordinates. The population of clusters in the central 10 kpc has a mean rotational velocity reaching 50-80 km/s, and a nearly isotropic velocity dispersion 100-120 km/s, while in the outer galaxy the cluster orbits are strongly radially anisotropic. We confirm a concentration of clusters at high radial action in the outer region of the Galaxy. Finally, we explore a range of equilibrium distribution function-based models for the entire globular cluster system, and the information they provide about the potential of the Milky Way. The dynamics of clusters is best described by models with the circular velocity between 10 and 50 kpc staying in the range 210-240 km/s.

Galactic gas outflows are driven by stellar feedback with dominant contribution from supernovae (SN) explosions. The question of whether the energy deposited by SNe initiates a large scale outflow or gas circulation on smaller scales -- between discs and intermediate haloes, depends on SN rate and their distribution in space and time. We consider here gas circulation by disc-wide unclustered SNe with galactic star formation rate in the range from {$\simeq 6\times 10^{-4}$ to $\simeq 6\times 10^{-2}~M_\odot$~yr$^{-1}$~kpc$^{-2}$, corresponding to mid-to-high star formation observed in galaxies. We show that such disc-wide SN explosion regime can form circulation of warm ($T\sim 10^4$ K) and cold ($T<10^3$ K) phases within a few gas scale heights, and elevation of hot {($T>10^5$ K)} gas at higher ($z>1$ kpc) heights. We found that the threshold energy input rate for hot gas outflows with disc-wide supernovae explosions is estimated to be of the order $\sim 4\times 10^{-4}$~erg~s$^{-1}$~cm$^{-2}$. We discuss the observational manifestations of such phenomena in optical and X-ray bands. In particular, we found that for face-on galaxies with SF ($\Sigma_{_{\rm SF}}>0.02~M_\odot$~yr$^{-1}$~kpc$^{-2}$), the line profiles of ions typical for warm gas show a double-peak shape, corresponding to out-of-plane outflows. In the X-ray bands, galaxies with high SF rates ($\Sigma_{_{\rm SF}}>0.006~M_\odot$~yr$^{-1}$~kpc$^{-2}$) can be bright, with a smooth surface brightness in low-energy bands ($0.1\hbox{--}0.3$~keV) and patchy at higher energies ($1.6\hbox{--}8.3$~keV).}

Agama is a publicly available software library for a broad range of applications in the field of stellar dynamics. It provides methods for computing the gravitational potential of arbitrary analytic density profiles or N-body models; orbit integration and analysis; transformations between position/velocity and action/angle variables; distribution functions expressed in terms of actions and their moments; iterative construction of self-consistent multicomponent galaxy models. Applications include the inference about the structure of Milky Way or other galaxies from observations of stellar kinematics; preparation of equilibrium initial conditions for N-body simulations; analysis of snapshots from simulations. The library is written in C++, provides a Python interface, and can be coupled to other stellar-dynamical software: Amuse, Galpy and Nemo.

We use the proper motions (PM) of half a million red giant stars in Large Magellanic Cloud measured by Gaia to construct a 2d kinematic map of mean PM and its dispersion across the galaxy, out to 7 Kpc from its centre. We then explore a range of dynamical models and measure the rotation curve, mean azimuthal velocity, velocity dispersion profiles, and the orientation of the galaxy. We find that the circular velocity reaches 90 km/s at 5 Kpc, and that the velocity dispersion ranges from 30-40 km/s in the galaxy centre to 15-20 km/s at 7 Kpc.

We consider cosmological implications of the formation of first stellar size black holes (BHs) in the universe. Such BHs form and grow by accretion in minihaloes of masses $\simeq 10^5-10^7~M_\odot$, and emit non-thermal radiation which impact the ionization and thermal state of the IGM. We compute the implications of this process. We show that the influence regions for hydrogen increase to 10kpc (physical length) for non-growing BHs to more than 0.3--1Mpc for accreting BHs, the influence regions are ten times smaller for singly ionized helium. We consider three possible observables from the influence zones around accreting BHs during $8.5<z<25$: HI 21cm line, hyperfine line of $^3$HeII, and HI recombination lines. We show that the 21cm emitting region around a growing BH could produce brightness temperatures $\simeq 15$mK across an evolving structure of 1Mpc in size with hot, ionized gas closer to the BH and much cooler in outer regions. We show that the ongoing and upcoming radio interferometers such as LOFAR and SKA1-LOW might be able to detect these regions. $^3$HeII emission from regions surrounding the growing BH is weak: the corresponding brightness temperatures reaches tens of nano-Kelvin, which is below the range of upcoming SKA1-MED. We show that for growing BHs H$\alpha$ line could be detected by JWST with $S/N=10$ in $10^4$~seconds of integration. In light on the recent EDGES result, we show that with additional cooling of baryons owing to collision with dark matter the HI signal could be enhanced by more than an order of magnitude.

The ground-breaking detections of gravitational waves from black hole mergers by LIGO have rekindled interest in primordial black holes (PBHs) and the possibility of dark matter being composed of PBHs. It has been suggested that PBHs of tens of solar masses could serve as dark matter candidates. Recent analytical studies demonstrated that compact ultra-faint dwarf galaxies can serve as a sensitive test for the PBH dark matter hypothesis, since stars in such a halo-dominated system would be heated by the more massive PBHs, their present-day distribution can provide strong constraints on PBH mass. In this study, we further explore this scenario with more detailed calculations, using a combination of dynamical simulations and Bayesian inference methods. The joint evolution of stars and PBH dark matter is followed with a Fokker-Planck code PhaseFlow. We run a large suite of such simulations for different dark matter parameters, then use a Markov Chain Monte Carlo approach to constrain the PBH properties with observations of ultra-faint galaxies. We find that two-body relaxation between the stars and PBH drives up the stellar core size, and increases the central stellar velocity dispersion. Using the observed half-light radius and velocity dispersion of stars in the compact ultra-faint dwarf galaxies as joint constraints, we infer that these dwarfs may have a cored dark matter halo with the central density in the range of 1-2 $\rm{M_{\odot}/pc^3}$, and that the PBHs may have a mass range of 2-14 $\rm{M_{\odot}}$ if they constitute all or a substantial fraction of the dark matter.

We present an approach for simulating the collisional evolution of spherical isotropic stellar systems based on the one-dimensional Fokker-Planck equation. A novel aspect is that we use the phase volume as the argument of the distribution function, instead of the traditionally used energy, which facilitates the solution. The publicly available code, PhaseFlow, implements a high-accuracy finite-element method for the Fokker-Planck equation, and can handle multiple-component systems, optionally with the central black hole and taking into account loss-cone effects and star formation. We discuss the energy balance in the general setting, and in application to the Bahcall-Wolf cusp around a central black hole, for which we derive a perturbative solution. We stress that the cusp is not a steady-state structure, but rather evolves in amplitude while retaining an approximately $\rho\propto r^{-7/4}$ density profile. Finally, we apply the method to the nuclear star cluster of the Milky Way, and illustrate a possible evolutionary scenario in which a two-component system of lighter main-sequence stars and stellar-mass black holes develops a Bahcall-Wolf cusp in the heavier component and a weaker $\rho\propto r^{-3/2}$ cusp in the lighter, visible component, over the period of several Gyr. The present-day density profile is consistent with the recently detected mild cusp inside the central parsec, and is weakly sensitive to initial conditions.

Recent observations suggest that stellar tidal disruption events (TDE) are strongly overrepresented in rare, post-starburst galaxies. Several dynamical mechanisms have been proposed to elevate their TDE rates, ranging from central stellar overdensities to the presence of supermassive black hole (SMBH) binaries. Another such mechanism, introduced here, is a radial velocity anisotropy in the nuclear star cluster produced during the starburst. These, and other, dynamical hypotheses can be disentangled by comparing observations to theoretical predictions for the TDE delay time distribution (DTD). We show that SMBH binaries are a less plausible solution for the post-starburst preference, as they can only reproduce the observed DTD with extensive fine-tuning. The overdensity hypothesis produces a reasonable match to the observed DTD (based on the limited data currently available), provided that the initial stellar density profile created during the starburst, $\rho(r)$, is exceptional in both steepness and normalization. In particular, explaining the post-starburst preference requires $\rho \propto r^{-\gamma}$ with $\gamma \gtrsim 2.5$, i.e. much steeper than the classic Bahcall-Wolf equilibrium profile of $\gamma = 7/4$. For "ultrasteep" density cusps ($\gamma \ge 9/4$), we show that the TDE rate decays with time measured since the starburst as $\dot{N} \propto t^{-(4\gamma-9)/(2\gamma-3)} / \ln t$. Radial anisotropies also represent a promising explanation, provided that initial anisotropy parameters of $\beta_0 \approx 0.5$ are sustainable against the radial orbit instability. TDE rates in initially anisotropic cusps will decay roughly as $\dot{N} \propto t^{-\beta_0}$. As the sample of TDEs with well-studied host galaxies grows, the DTD will become a powerful tool for constraining the exceptional dynamical properties of post-starburst galactic nuclei.

We study the merging and evolution of isolated supernovae (SNe) remnants in a stellar cluster into a collective superbubble, with the help of 3-D hydrodynamic simulations. We particularly focus on the transition stage when the isolated SNe remnants gradually combine to form a superbubble. We find that when the SN rate is high ($\nu_{\rm sn}\sim 10^{-9}$ pc$^{-3}$ yr$^{-1}$), the merging phase lasts for $\sim 10^4$ yr, for $n=1\hbox{--}10$ cm$^{-3}$, and the merging phase lasts for a longer time ($\sim 0.1$ Myr or more) for lower SN rates ($\nu_{\rm sn}\le 10^{-10}$ pc$^{-3}$ yr$^{-1}$). During this transition phase, the growing superbubble is filled with dense and cool fragments of shells and most of the energy is radiated away during this merging process. After passing through the intermediate phase, the superbubble eventually settles on to a new power-law wind asymptote that is smaller than estimated in a continuous wind model. This results in a significant (more than {\it several times}) underestimation of the mechanical luminosity needed to feed the bubble. We determine the X-ray and H$\alpha$ surface brightnesses as functions of time for such merging SNe in a stellar cluster and find that clusters with high SN rate shine predominantly in soft X-rays and H$\alpha$. In particular, a low value of the volume averaged H$\alpha$ to H$\beta$ ratio and its large spread can be a good indicator of the transition phase of merging SNe.

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 explore the rates of tidal disruption events (TDEs) of stars by supermassive black holes (SBHs) in galactic nuclei formed in mergers followed by a formation and coalescence of a binary SBH. Such systems initially have a deficit of stars on low-angular-momentum orbits caused by the slingshot process during the binary SBH stage, which tends to reduce the flux of stars into the SBH compared to the steady-state value. On the other hand, a newly formed galactic nucleus has a non-spherical shape which enhances the mixing of stars in angular momentum and thus the TDE rate. In galaxies with relatively low SBH masses (<10^7 Msun), relaxation times are short enough to wash out the anisotropy in initial conditions, and for more massive SBH the enhancement of flux due to non-sphericity turns out to be more important than its suppression due to initial anisotropy. Therefore, the present-day TDE rates generally exceed conventional steady-state estimates based on a spherical isotropic approximation. We thus conjecture that the lower observationally inferred TDE rates compared to theoretical predictions cannot be attributed to the depletion of low-angular-momentum stars by SBH binaries.

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.

A binary supermassive black hole loses energy via ejection of stars in a galactic nucleus, until emission of gravitational waves becomes strong enough to induce rapid coalescence. Evolution via the gravitational slingshot requires that stars be continuously supplied to the binary, and it is known that in spherical galaxies the reservoir of such stars is quickly depleted, leading to stalling of the binary at parsec-scale separations. Recent N-body simulations of galaxy mergers and isolated nonspherical galaxies suggest that this stalling may not occur in less idealized systems. However, it remains unclear to what degree these conclusions are affected by collisional relaxation, which is much stronger in the numerical simulations than in real galaxies. In this study, we present a novel Monte Carlo method that can efficiently deal with both collisional and collisionless dynamics, and with galaxy models having arbitrary shapes. We show that without relaxation, the final-parsec problem may be overcome only in triaxial galaxies. Axisymmetry is not enough, but even a moderate departure from axisymmetry is sufficient to keep the binary shrinking. We find that the binary hardening rate is always substantially lower than the maximum possible, "full-loss-cone" rate, and that it decreases with time, but that stellar-dynamical interactions are nevertheless able to drive the binary to coalescence on a timescale <=1 Gyr in any triaxial galaxy.

The evolution of a supernova remnant in a cloudy medium as a function of the volume filling factor of the clouds is studied in a three-dimensional axially symmetrical model. The model includes the mixing of heavy elements (metals) ejected by the supernova and their contribution to radiative losses. The interaction of the supernova envelope with the cloudy phase of the interstellar medium leads to nonsimultaneous, and on average earlier, onsets of the radiative phase in different parts of the supernova envelope. Growth in the volume filling factor $f$ leads to a decrease in the time for the transition of the envelope to the radiative phase and a decrease in the envelope's mean radius, due to the increased energy losses by the envelope in the cloudy medium. When the development of hydrodynamical instabilities in the supernova envelope is efficient, the thermal energy falls as $E_t\sim t^{-2.3}$, for the propagation of the supernova remnant through either a homogeneous or a cloudy medium. When the volume filling factor is $f\simgt 0.1$, a layer with excess kinetic energy andmomentumforms far behind the global shock front from the supernova, which traps the hot gas of the cavity in the central part of the supernova remnant. Metals ejected by the supernova are also enclosed in the central region of the remnant, where the initial (high) metallicity is essentially preserved. Thus, the interaction of the supernova envelope with the cloudy interstellar medium appreciably changes the dynamics and structure of the distribution of the gas in the remnant. This affects the observational characteristics of the remnant, in particularly, leading to substantial fluctuations of the emission measure of the gas with $T>10^5$~K and the velocity dispersion of the ionized gas.

We compute the rates of capture of stars by supermassive black holes, using time-dependent Fokker--Planck equation with initial conditions that have a deficit of stars on low-angular-momentum orbits. One class of initial conditions has a gap in phase space created by a binary black hole, and the other has a globally tangentially-anisotropic velocity distribution. We find that for galactic nuclei that are younger than ~0.1 relaxation times, the flux of stars into the black hole is suppressed with respect to the steady-state value. This effect may substantially reduce the number of observable tidal disruption flares in galaxies with black hole masses M>10^7 Msun.

We present an implementation of the Schwarzschild orbit superposition method which can be used for constructing self-consistent equilibrium models of barred or non-barred disc galaxies, or of elliptical galaxies with figure rotation. This is a further development of the publicly available code SMILE; its main improvements include a new efficient representation of an arbitrary gravitational potential using two-dimensional spline interpolation of Fourier coefficients in the meridional plane, as well as the ability to deal with rotation of the density profile and with multicomponent mass models. We compare several published methods for constructing composite axisymmetric disc--bulge--halo models and demonstrate that our code produces the models that are closest to equilibrium. We also apply it to create models of triaxial elliptical galaxies with cuspy density profiles and figure rotation, and find that such models can be found and are stable over many dynamical times in a wide range of pattern speeds and angular momenta, covering both slow- and fast-rotator classes. We then attempt to create models of strongly barred disc galaxies, using an analytic three-component potential, and find that it is not possible to make a stable dynamically self-consistent model for this density profile. Finally, we take snapshots of two N-body simulations of barred disc galaxies embedded in nearly-spherical haloes, and construct equilibrium models using only information on the density profile of the snapshots. We demonstrate that such reconstructed models are in near-stationary state, in contrast with the original N-body simulations, one of which displayed significant secular evolution.

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.

I review the evolution of binary supermassive black holes and focus on the stellar-dynamical mechanisms that may help to overcome the final-parsec problem - the possible stalling of the binary at a separation much larger than is required for an efficient gravitational wave emission. Recent N-body simulations have suggested that a departure from spherical symmetry in the nucleus of the galaxy may keep the rate of interaction of stars with the binary at a high enough level so that the binary continues to shrink rather rapidly. However, a major problem of all these simulations is that they do not probe the regime where collisionless effects are dominant - in other words, the number of particles in the simulation is still not sufficient to reach the asymptotic behaviour of the system. I present a novel Monte Carlo method for simulating both collisional and collisionless evolution of non-spherical stellar systems, and apply it for the problem of binary supermassive black hole evolution. I show that in triaxial galaxies the final-parsec problem is largely non-existent, while in the axisymmetric case it seems to still exist in the limit of purely collisionless regime relevant for real galaxies, but disappears in the N-body simulations where the feasible values of N are still too low to get rid of collisional effects.

We have developed a novel Monte Carlo method for simulating the dynamical evolution of stellar systems in arbitrary geometry. The orbits of stars are followed in a smooth potential represented by a basis-set expansion and perturbed after each timestep using local velocity diffusion coefficients from the standard two-body relaxation theory. The potential and diffusion coefficients are updated after an interval of time that is a small fraction of the relaxation time, but may be longer than the dynamical time. Thus our approach is a bridge between the Spitzer's formulation of the Monte Carlo method and the temporally smoothed self-consistent field method. The primary advantages are the ability to follow the secular evolution of shape of the stellar system, and the possibility of scaling the amount of two-body relaxation to the necessary value, unrelated to the actual number of particles in the simulation. Possible future applications of this approach in galaxy dynamics include the problem of consumption of stars by a massive black hole in a non-spherical galactic nucleus, evolution of binary supermassive black holes, and the influence of chaos on the shape of galaxies, while for globular clusters it may be used for studying the influence of rotation.

We consider the problem of star consumption by supermassive black holes in non-spherical (axisymmetric, triaxial) galactic nuclei. We review the previous studies of the loss-cone problem and present a novel simulation method which allows to separate out the collisional (relaxation-related) and collisionless (related to non-conservation of angular momentum) processes and determine their relative importance for the capture rates in different geometries. We show that for black holes more massive than 10^7 Msun, the enhancement of capture rate in non-spherical galaxies is substantial, with even modest triaxiality being capable of keeping the capture rate at the level of few percent of black hole mass per Hubble time.

Using 3D numerical simulations we study the evolution of the H$\alpha$ intensity and velocity dispersion for single and multiple supenova (SN) explosions. We find that the $I_{\rm H\alpha}-\sigma$ diagram obtained for simulated gas flows is similar in shape to that observed in dwarf galaxies. We conclude that colliding SN shells with significant difference in age are resposible for high velocity dispersion that reaches values high as $\simgt 100$kms$^{-1}$. Such a high velocity dispersion could be hardly got for a single SN remnant. Peaks of velocity dispersion on the $I_{\rm H\alpha}-\sigma$ diagram may correspond to several stand-alone or merged SN remnants with moderately different ages. The procedure of the spatial resolution degrading in the H$\alpha$ intensity and velocity dispersion maps makes the simulated $I_{\rm H\alpha}-\sigma$ diagrams close to those observed in dwarf galaxies not only in shape, but also quantitatively.

We consider evolution of metal-enriched gas exposed to a superposition of time-dependent radiation field of a nearby starburst galaxy and nearly invariant (on timescales 100 Myr) extragalactic ionization background. Within nonequilibrium (time-dependent) photoionization models we determine ionization fraction of the OVI ion commonly observed in galactic circumference. We derive then conditions for OVI to appear in absorptions in extended galactic haloes depending on the galactic mass and star formation rate. We have found that the maximum OVI fraction can reach $\sim 0.4-0.9$ under combined action of the galactic and the extragalactic ionizing radiation fields. We conclude that soft X-ray emission with $E\simgt 113$~eV from the stellar population of central starforming galaxies is the main source of such a high fraction of OVI. This circumstance can explain high column densities ${\rm N(OVI)} \sim 10^{14.5 - 15.3}$~cm$^{-2}$ observed in the haloes of starforming galaxies at low redshifts (Tumlinson etal 2011) {\it even} for a relatively low ({ $\sim 0.01-0.1\zsun$}) metallicity. As a result, the requirements to the sources of oxygen in the extended haloes relax to a reasonably conservative level. We show that at $z\simlt 0.5$ ionization kinetics of oxygen in a relatively dense plasma $n\simgt 10^{-4}$ cm$^{-3}$ of outer halo exposed to a low extragalactic ionizing flux is dominated by nonequilibrium effects.

(Abridged) Heating of the interstellar medium by multiple supernovae (SNe) explosions is at the heart of producing galaxy-scale outflows. We use hydrodynamical simulations to study the efficiency of multiple SNe in heating the interstellar medium (ISM) and filling the volume with gas of high temperatures. We argue that it is important for SNe remnants to have a large filling factor {\it and} a large heating efficiency. For this, they have to be clustered in space and time, and keep exploding until the hot gas percolates through the whole region, in order to compensate for the radiative loss. In the case of a limited number of SNe, we find that although the filling factor can be large, the heating efficiency declines after reaching a large value. In the case of a continuous series of SNe, the hot gas ($T \ge 3 \times 10^6$ K) can percolate through the whole region after the total volume filling factor reaches a threshold of $\sim 0.3$. The efficiency of heating the gas to X-ray temperatures can be $\ge 0.1$ after this percolation epoch, which occurs after a period of $\approx 10$ Myr for a typical starburst SNe rate density of $\nu_{\rm SN} \approx 10^{-9}$ pc$^{-3}$ yr$^{-1}$ and gas density of $n\approx 10$ cm$^{-3}$ in starburst nuclei regions. This matches the recent observations of a time delay of similar order between the onset of star formation and galactic outflows. The efficiency to heat gas up to X-ray temperatures ($\ge 10^{6.5}$ K) roughly scales as $\nu_{\rm SN}^{0.2} n^{-0.6}$. For a typical SNe rate density and gas density in starburst nuclei, the heating efficiency is $\sim 0.15$, also consistent with previous interpretations from X-ray observations. We discuss the implications of our results with regard to observational diagnostics of ionic ratios and emission measures in starburst nuclei regions.

We consider the evolution of supermassive black hole binaries at the center of spherical, axisymmetric, and triaxial galaxies, using direct N-body integrations as well as analytic estimates. We find that the rates of binary hardening exhibit a significant N-dependence in all the models, at least for N in the investigated range of 10^5<=N<=10^6. Binary hardening rates are also substantially lower than would be expected if the binary loss cone remained full, as it would be if the orbits supplying stars to the binary were being efficiently replenished. The difference in binary hardening rates between the spherical and nonspherical models is less than a factor of two even in the simulations with the largest N. By studying the orbital populations of our models, we conclude that the rate of supply of stars to the binary via draining of centrophilic orbits is indeed expected to be much lower than the full-loss-cone rate, consistent with our simulations. We argue that the binary's evolution in the simulations is driven in roughly equal amounts by collisional and collisionless effects, even at the highest N-values currently accessible. While binary hardening rates would probably reach a limiting value for large N, our results suggest that we cannot approach that rate with currently available algorithms and computing hardware. The extrapolation of results from N-body simulations to real galaxies is therefore not straightforward, casting doubt on recent claims that triaxiality or axisymmetry alone are capable of solving the final-parsec problem in gas-free galaxies.

We study the impact of primordial magnetic fields on the HI absorption from the Epoch of Reionization. The presence of these fields result in two distinct effects: (a) the heating of the haloes from the decay of magnetic fields owing to ambipolar diffusion, and (b) an increase in the number of haloes owing to additional matter fluctuations induced by magnetic fields. We analyse both these effects and show that the latter is potentially observable because the number of haloes along of line of sight can increase by many orders of magnitude. While this effect is not strongly dependent on the magnetic field strength in the range $0.3\hbox{--}0.6$ nG, it is extremely sensitive to the magnetic field power spectral index for the near scale-free models. Therefore the detection of such absorption features could be a sensitive probe of the primordial magnetic field and its power spectrum. We discuss the detectability of these features with the ongoing and future radio interferometers. In particular, we show that LOFAR might be able to detect these absorption features at $z \simeq 10$ in less than 10 hrs of integration if the flux of the background source is 400 mJy.

The imprint of decaying dark matter (DM) particles in characteristics of the "21 cm fores" -- absorptions in 21 cm from minihalos in spectra of distant radio-loud sources -- is considered within a 1D self-consistent hydrodynamic description of minihalos from their turnaround point to virialization. The most pronounced influence of decaying DM on evolution of minihalos is found in the mass range $M=10^5-10^6\msun$, for which unstable DM with the current upper limit of the ionization rate $\xi_{L} = 0.59\times 10^{-25}$ s$^{-1}$ depresses 21 cm optical depth by an order of magnitude compared to the standard recombination scenario. Even rather a modest ionization $\xi \sim 0.3\xi_L$ practically "erases" absorption features and results in a considerable decrease (by factor of more than 2.5) of the number of strong ($W_\nu^{obs} \simgt 0.3$ kHz at $z\simeq 10$) absorptions. At such circumstances broad-band observations are to be more suitable for inferring physical conditions of the absorbing gas. X-ray photons from stellar activity of initial episodes of star formation can compete the contribution from decaying DM only at $z<10$. Therefore, when observed the 21 cm signal will allow to follow evolution of decaying DM particles in the redshift range $z=10-15$. Contrary, in case of non-detection of the 21 cm signal in the frequency range $\nu<140$ MHz a lower limit on the ionization rate from decaying dark matter can be established.

We consider the problem of consumption of stars by a supermassive black hole (SBH) at the center of an axisymmetric galaxy. Inside the SBH sphere of influence, motion of stars in the mean field is regular and can be described analytically in terms of three integrals of motion: the energy E, the z-component of angular momentum L_z, and the secular Hamiltonian H. There exist two classes of orbits, tubes and saucers; saucers occupy the low-angular-momentum parts of phase space and their fraction is proportional to the degree of flattening of the nucleus. Perturbations due to gravitational encounters lead to diffusion of stars in integral space, which can be described using the Fokker-Planck equation. We calculate the diffusion coefficients and solve this equation in the two-dimensional phase space (L_z, H), for various values of the capture radius and the degree of flattening. Capture rates are found to be modestly higher than in the spherical case, up to a factor of a few, and most captures take place from saucer orbits. We also carry out a set of collisional N-body simulations to confirm the predictions of the Fokker-Planck models. We discuss the implications of our results for rates of tidal disruption and capture in the Milky Way and external galaxies.

We review the methods used to study the orbital structure and chaotic properties of various galactic models and to construct self-consistent equilibrium solutions by Schwarzschild's orbit superposition technique. These methods are implemented in a new publicly available software tool, SMILE, which is intended to be a convenient and interactive instrument for studying a variety of 2D and 3D models, including arbitrary potentials represented by a basis-set expansion, a spherical-harmonic expansion with coefficients being smooth functions of radius (splines), or a set of fixed point masses. We also propose two new variants of Schwarzschild modelling, in which the density of each orbit is represented by the coefficients of the basis-set or spline spherical-harmonic expansion, and the orbit weights are assigned in such a way as to reproduce the coefficients of the underlying density model. We explore the accuracy of these general-purpose potential expansions and show that they may be efficiently used to approximate a wide range of analytic density models and serve as smooth representations of discrete particle sets (e.g. snapshots from an N-body simulation), for instance, for the purpose of orbit analysis of the snapshot. For the variants of Schwarzschild modelling, we use two test cases - a triaxial Dehnen model containing a central black hole, and a model re-created from an N-body snapshot obtained by a cold collapse. These tests demonstrate that all modelling approaches are capable of creating equilibrium models.

We present self-consistent calculations of non-equilibrium (time-dependent) cooling rates for a dust-free collisionally controlled gas in wide temperature ($10 K\le T\le 10^8 K$) and metallicity ($10^{-4} Z_\odot \le Z \le 2 Z_\odot$) ranges. We confirm that molecular hydrogen dominates cooling at $10^2 \simlt T\simlt 10^4$ K and $Z\simlt 10^{-3} Z_\odot$. We find that the contribution from H$_2$ into cooling rate around $T\sim (4-5)\times 10^3$ K stimulates thermal instability in the metallicity range $Z\simlt 10^{-2} Z_\odot$. Isobaric cooling rates are generally lower than isochoric ones, because the associated increase of gas density leads to both more efficient hydrogen recombination and equilibration of the fine-structure level populations. Isochoric cooling keeps the ionization fraction remains quite high at $T\simlt10^4$ K: up to $\sim0.01$ at $T\simeq 10^3$ K and $Z\simlt 0.1 Z_\odot$, and even higher at higher metallicity, contrary to isobaric cooling where it at least an order of magnitude lower. Despite this increase in ionization fraction the gas-phase formation rate of molecular hydrogen (via H$^-$) lowers with metallicity, because higher metallicity shorttens the evolution time. We implement our self-consistent cooling rates into the multi-dimensional parallel code ZEUS-MP in order to simulate evolution of a supernova remnant, and compare it with an analogous model with tabulated cooling rates published in previous works. We find significant differences between the two descriptions, which may appear, e.g., in mixing of the ejected metals in the circumstellar medium.

The paper considers the evolution of the supernova envelopes produced by Population III stars with masses of $M_*\sim 25-200 M_\odot$ located in non-rotating protogalaxies with masses of $M\sim 10^7 M_\odot$ at redshifts $z=12$, with dark-matter density profiles in the form of modified isothermal spheres. The supernova explosion occurs in the ionization zone formed by a single parent star. The properties of the distribution of heavy elements (metals) produced by the parent star are investigated, as well as the efficiency with which they are mixed with the primordial gas in the supernova envelope. In supernovae with high energies ($E\simgt 5\times 10^{52}$ erg), an appreciable fraction of the gas can be ejected from the protogalaxy, but nearly all the heavy elements remain in the protogalaxy. In explosions with lower energies ($E\simlt 3\times 10^{52}$ erg), essentially no gas and heavy elements are lost from the protogalaxy: during the first one to threemillion years, the gas and heavy elements are actively carried from the central region of the protogalaxy ($r\sim 0.1 r_{vir}$, where $r_{vir}$ is the virial radius of the protogalaxy), but an appreciable fraction of the mass of metals subsequently returns when the hot cavity cools and the envelope collapses. Supernovae with high energies ($E\simgt 5\times 10^{52}$ erg) are characterized by a very low efficiency of mixing of metals; their heavy elements are located in the small volume occupied by the disrupted envelope (in a volume comparable with that of the entire envelope), with most of the metals remaining inside the hot, rarified cavity of the envelope. (abridged)

The spin angular momentum S of a supermassive black hole (SBH) precesses due to torques from orbiting stars, and the stellar orbits precess due to dragging of inertial frames by the spinning hole. We solve the coupled post-Newtonian equations describing the joint evolution of S and the stellar angular momenta Lj, j = 1...N in spherical, rotating nuclear star clusters. In the absence of gravitational interactions between the stars, two evolutionary modes are found: (1) nearly uniform precession of S about the total angular momentum vector of the system; (2) damped precession, leading, in less than one precessional period, to alignment of S with the angular momentum of the rotating cluster. Beyond a certain distance from the SBH, the time scale for angular momentum changes due to gravitational encounters between the stars is shorter than spin-orbit precession times. We present a model, based on the Ornstein-Uhlenbeck equation, for the stochastic evolution of star clusters due to gravitational encounters and use it to evaluate the evolution of S in nuclei where changes in the Lj are due to frame dragging close to the SBH and to encounters farther out. Long-term evolution in this case is well described as uniform precession of the SBH about the cluster's rotational axis, with an increasingly important stochastic contribution when SBH masses are small. Spin precessional periods are predicted to be strongly dependent on nuclear properties, but typical values are 10-100 Myr for low-mass SBHs in dense nuclei, 100 Myr - 10 Gyr for intermediate mass SBHs, and > 10 Gyr for the most massive SBHs. We compare the evolution of SBH spins in stellar nuclei to the case of torquing by an inclined, gaseous accretion disk.

We consider the evolution of the ionization zone around Population III stars with $M_*\sim 25-200 M_\odot$ in protogalaxies with $M\sim 10^7 M_\odot$ at redshifts $z = 12$, assuming that the dark matter profile is a modified isothermal sphere. We study the conditions for the growth of instabilities in the ionization zones. The Rayleigh-Taylor and thermal instabilities develop efficiently in the ionization zones around 25-40 $M_\odot$ stars, while this efficiency is lower for stars with $\sim 120 M_\odot$. For more massive stars ($\sim 200 M_\odot$), the flux of ionizing photons is strong enough to considerably reduce the gas density in the ionization zone, and the typical lifetimes of stars ($\sim 2$ Myr) are insufficient for the growth of instabilities. The gas in a protogalaxy with $M\sim 10^7 M_\odot$ with a 200 $M_\odot$ central star is completely ionized by the end of the star's lifetime; in the case of a 120 $M_\odot$ central star, only one-third of the total mass of gas is ionized. Thus, ionizing photons from stars with $M_*\simlt 120 M_\odot$ cannot leave protogalaxies with $M\simgt 10^7 M_\odot$. If the masses of the central stars are 25 and 40 $M_\odot$, the gas in protogalaxies of this mass remains essentially neutral. We discuss the consequences of the evolution of the ionization zones for the propagation of the envelope after the supernova explosions of the stars and the efficiency of enrichment of the intergalactic medium in heavy elements.

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.

The absorption properties of the first low-mass protogalaxies (mini-halos) forming at high redshifts in the 21-cm line of atomic hydrogen are considered. The absorption properties of these protogalaxies are shown to depend strongly on both their mass and evolutionary status. The optical depths in the line reach $\sim$0.1-0.2 for small impact parameters of the line of sight. When a protogalaxy being compressed, the influence of gas accretion can be seen manifested in a non-monotonic frequency dependence of the optical depth. The absorption characteristics in the 21-cm line are determined by the thermal and dynamical evolution of the gas in protogalaxies. Since the theoretical line width in the observer's reference frame is 1-6 kHz and the expected separation between lines 8.4 kHz, the lines from low mass protogalaxies can be resolved using ongoing and future low frequency interferometers.

We consider the influence of decaying dark matter (DM) particles on the characteristics of 21 cm absorption in spectra of distant radio-loud sources - "21 cm forest" - from minihaloes with masses $M=10^5-10^7\msun$ virialized at $z_{vir} = 10$. We use 1D self-consistent hydrodynamic description to study evolution of minihaloes, and follow up their absorption characteristics from turnaround to virialization. We find that in the presence of decaying dark matter both thermal and dynamical evolution of minihaloes demonstrate significant deviation from those in the model without dark matter decay (standard recombination). We show that optical depth in 21 cm line is strongly suppressed in the presence of decaying particles: for $M=10^5-10^6\msun$ decaying dark matter with the energy rate deposited in baryonic gas $\xi_{L} = 0.59\times 10^{-25}$ s$^{-1}$ - the current upper limit of the energy deposit - decreases the optical depth and the equivalent width by an order of magnitude compared to the standard recombination. Thus additional ionization and heating from decaying DM particles almost "erases" absorption features from minihaloes with $M=10^5-10^6\msun$ for $\xi \simgt 0.3\xi_L$, which consequently considerably decreases the number of strong absorptions: for example, the number of absorptions with the equivalent width $W_\nu^{obs} \simgt 0.3$ kHz at $z\simeq10$ decreases more than 2.5 times for $\xi/\xi_{L} = 0.3$ and $\simgt$4.5 times for $\xi/\xi_{L} = 1$. We argue that "21 cm forest" absorptions might be a powerful probe of the presence of decaying dark matter in the early Universe.

We use both N-body simulations and integration in fixed potentials to explore the stability and the long-term secular evolution of self-consistent, equilibrium, non-rotating, triaxial spheroidal galactic models. More specifically, we consider Dehnen models built with the Schwarzschild method. We show that short-term stability depends on the degree of velocity anisotropy (radially anisotropic models are subject to rapid development of radial-orbit instability). Long-term stability, on the other hand, depends mainly on the properties of the potential, and in particular, on whether it admits a substantial fraction of strongly chaotic orbits. We show that in the case of a weak density cusp (gamma=1 Dehnen model) the N-body model is remarkably stable, while the strong-cusp (gamma=2) model exhibits substantial evolution of shape away from triaxiality, which we attribute to the effect of chaotic diffusion of orbits. The different behaviour of these two cases originates from the different phase space structure of the potential; in the weak-cusp case there exist numerous resonant orbit families that impede chaotic diffusion. We also find that it is hardly possible to affect the rate of this evolution by altering the fraction of chaotic orbits in the Schwarzschild model, which is explained by the fact that the chaotic properties of an orbit are not preserved by the N-body evolution. There are, however, parameters in Schwarzschild modelling that do affect the stability of an N-body model, so we discuss the recipes how to build a `good' Schwarzschild model.

We have presented the non-equilibrium (time-dependent) cooling rate and ionization state calculations for a gas behind shock waves with $v \sim 50-150$ km s$^{-1}$ ($T_s \sim 0.5 - 6\times 10^5$ K). Such shock waves do not lead to the radiative precursor formation, i.e. the thermal evolution of a gas behind the shock waves are controlled by collisions only. We have found that the cooling rate in a gas behind the shock waves with $v \sim 50-120$ km s$^{-1}$ ($T_s \sim 0.5 - 3\times 10^5$ K) differs considerably from the cooling rate for a gas cooled from $T = 10^8$ K. It is well-known that a gas cooled from $T = 10^8$ K is thermally unstable for isobaric and isochoric perturbations at $T \simgt 2\times 10^4$ K. We have studied the thermal instability in a collisionally controlled gas for shock waves with $v \sim 50-150$ km s$^{-1}$. We have found that the temperature range, where the postshock gas is thermally unstable, is significantly modified and depends on both gas metallicity and ionic composition of a gas before shock wave. For $Z \simgt 0.1Z_\odot$ the temperature range, where the thermal instability criterion for isochoric perturbations is not fulfilled, widens in comparison with that for a gas cooled from $T = 10^8$ K, while that for isobaric perturbations remains almost without a change. For $Z\sim Z_\odot$ a gas behind shock waves with $v \simlt 65$ km s$^{-1}$ ($T_s \simlt 10^5$ K) is thermally stable to isochoric perturbations during full its evolution. We have shown that the transition from isobaric to isochoric cooling for a gas with $Z \simgt 0.1Z_\odot$ behind shock waves with $T_s = 0.5 - 3\times 10^5$ K proceeds at lower column density layer behind a shock wave than that for a gas cooled from $T = 10^8$ K. (abridged)

Nonequilibrium (time-dependent) cooling rates and ionization state calculations are presented for low-density gas enriched with heavy elements (metals) and photoionized by external ultraviolet/X-ray radiation. We consider a wide range of gas densities and metallicities and also two types of external radiation field: a power-law and the extragalactic background spectra. We have found that both cooling efficiencies and ionic composition of enriched photoionized gas depend significantly on the gas metallicity and density, the flux amplitude, and the shape of ionizing radiation spectrum. The cooling rates and ionic composition of gas in nonequilibrium photoionization models differ strongly (by a factor of several) from those in photoequilibrium due to overionization of the ionic states in the nonequilibrium case. The difference is maximal at low values of the ionization parameter and similar in magnitude to that between the equlibrium and nonequilibrium cooling rates in the collisionally controlled gas. In general, the nonequilibrium effects are notable at $T\simlt 10^6$ K. In this temperature range, the mismatch of the ionic states and their ratios between the photoequilibrium and the photo-nonequilibrium models reach a factor of several. The net result is that the time-dependent energy losses due to each chemical element (i.e. the contributions to the total cooling rate) differ singificantly from the photoequilibrium ones. We advocate the use of nonequilibrium cooling rates and ionic states for gas with near-solar (and above) metallicity exposed to an arbitrary ionizing radiation flux. We provide a parameter space (in terms of temperature, density, metallicity and ionizing radiation flux), where the nonequilibrium cooling rates are to be used. (abridged)

We discuss the properties of orbits within the influence sphere of a supermassive black hole (BH), in the case that the surrounding star cluster is nonaxisymmetric. There are four major orbit families; one of these, the pyramid orbits, have the interesting property that they can approach arbitrarily closely to the BH. We derive the orbit-averaged equations of motion and show that in the limit of weak triaxiality, the pyramid orbits are integrable: the motion consists of a two-dimensional libration of the major axis of the orbit about the short axis of the triaxial figure, with eccentricity varying as a function of the two orientation angles, and reaching unity at the corners. Because pyramid orbits occupy the lowest angular momentum regions of phase space, they compete with collisional loss cone repopulation and with resonant relaxation in supplying matter to BHs. General relativistic advance of the periapse dominates the precession for sufficiently eccentric orbits, and we show that relativity imposes an upper limit to the eccentricity: roughly the value at which the relativistic precession time is equal to the time for torques to change the angular momentum. We argue that this upper limit to the eccentricity should apply also to evolution driven by resonant relaxation, with potentially important consequences for the rate of extreme-mass-ratio inspirals in low-luminosity galaxies. In giant galaxies, we show that capture of stars on pyramid orbits can dominate the feeding of BHs, at least until such a time as the pyramid orbits are depleted; however this time can be of order a Hubble time.

We constrain the spectrum of the cosmic ultraviolet background radiation by fitting the observed abundance ratios carbon ions at $z\sim 2\hbox{--}3$ with those expected from different models of the background radiation. We use the recently calculated modulation of the background radiation between 3 and 4 Ryd due to resonant line absorption by intergalactic HeII, and determine the ratios of CIII to CIV expected at these redshifts, as functions of metallicity, gas density and temperature. Our analysis of the observed ratios shows that 'delayed reionization' models, which assume a large fraction of HeII at $z\sim3$, is not favored by data. Our results suggest that HeII reionization was inhomogeneous, consistent with the predictions from recent simulations.

The evolution of dark matter in central areas of galaxies is considered (the Milky Way is taken as an example). It is driven by scattering off of dark matter particles by bulge stars, their absorption by the supermassive black hole and self-annihilation. This process is described by diffusion equation in the phase space of energy and angular momentum. The equation was integrated for several different models of initial dark matter distribution and using various assumptions about the dynamical factors. It turns out that because the Milky Way center is rather dynamically old (~5 relaxation times t_r), the difference in initial conditions almost vanishes. The density attains a nearly universal profile, and the gamma-ray flux from dark matter annihilation lies in rather narrow range, which enables more robust determination of the dark matter parameters. By present time the mass of dark matter inside the black hole sphere of influence (r<2 pc) has been reduced approximately twice, mostly because of heating by stars. It is shown that the dynamics of dark matter for t>t_r is determined mainly by stars outside the sphere of influence.

(Abridged) We explore numerically the chemical, thermal, and dynamical evolution of a shell formed by a high-energy supernova explosion ($10^{53}$ erg) in dwarf protogalaxies with total mass $10^7 M_\odot$ at a redshift $z=12$. We consider two initial configurations for the baryonic matter, one without rotation and the other having the ratio of rotational to gravitational energy $\beta=0.17$. The (non-rotating) dark matter halo is described by a quasi-isothermal sphere. We find that the dynamics of the shell is different in protogalaxies with and without rotation. For instance, the Rayleigh-Taylor instability in the shell develops faster in protogalaxies without rotation. The fraction of a blown-away baryonic mass is approximately twice as large in models with rotation ($\sim 20%$) than in models without rotation. On the other hand, the chemical evolution of gas in protogalaxies with and without rotation is found to be similar. The relative number densities of molecular hydrogen and HD molecules in the cold gas ($T \le 10^3$ K) saturate at typical values of $10^{-3}$ and $10^{-7}$, respectively. The clumps formed in the fragmented shell move with velocities that are at least twice as large as the escape velocity. The mass of the clumps is $\sim 0.1-10 \msun$, which is lower than the Jeans mass. We conclude that the clumps are pressure supported. A supernova explosion with energy $10^{53}$ ergs destructs our model protogalaxy. The clumps formed in the fragmented shell are pressure supported. We conclude that protogalaxies with total mass $\sim 10^{7} M_\odot$ are unlikely to form stars due to high-energy supernova explosions of the first stars.

Thermal regime of the baryons behind shock waves arising in the process of virialization of dark matter halos is governed at cetrain conditions by radiation of HD lines. A small fraction of the shocked gas can cool down to the temperature of the cosmic microwave background (CMB). We estimate an upper limit for this fraction: at $z=10$ it increases sharply from about $q_{_T}\sim 10^{-3}$ for dark halos of $M=5\times 10^7\msun$ to $\sim 0.1$ for halos with $M=10^8\msun$. Further increase of the halo mass does not lead however to a significant growth of $q_T$ -- the asymptotic value for $M\gg 10^8\msun$ is of 0.3. We estimate star formation rate associated with such shock waves, and show that they can provide a small but not negligible fraction of the star formation. We argue that extremely metal-poor low-mass stars in the Milky Way may have been formed from primordial gas behind such shocks.

Dark matter annihilation in so-called ``spikes'' near black holes is believed to be an important method of indirect dark matter detection. In the case of circular particle orbits, the density profile of dark matter has a plateau at small radii, the maximal density being limited by the annihilation cross-section. However, in the general case of arbitrary velocity anisotropy the situation is different. Particulary, for isotropic velocity distribution the density profile cannot be shallower than r^{-1/2} in the very centre. Indeed, a detailed study reveals that in many cases the term ``annihilation plateau'' is misleading, as the density actually continues to rise towards small radii and forms a weak cusp, rho ~ r^{-(beta+1/2)}, where beta is the anisotropy coefficient. The annihilation flux, however, does not change much in the latter case, if averaged over an area larger than the annihilation radius.

The influence of ultra-high energy cosmic rays (UHECRs) and decaying dark matter particles on the emission and absorption characteristics of neutral hydrogen in 21 cm at redshifts $z = 10-50$ is considered. In presence of UHECRs 21 cm can be seen in absorption with the brightness temperature $T_b=-(5-10)$ mK in the range $z=10-30$. Decayng particles can stimulate a 21 cm signal in emission with $T_b\sim 50-60$ mK at $z =50$, and $T_b \simeq 10$ mK at $z \sim 20$. Characteristics of the fluctuations of the brightness temperature, in particular, its power spectrum are also calculated. The maps of the power spectrum of the brightness temperature on the plane {\it wavenumber-redshift} are shown to be sensitive to the parameters of UHECRs and decaying dark matter. Observational possibilities to detect manifestations of UHECRs and/or decaying particles in 21 cm with the future radio telescopes (LOFAR, 21CMA and SKA), and to distinguish contributions from them are briefly discussed.

The presence of ultra-high-energy cosmic rays (UHECR) results in an increase in the degree of ionization in the post-recombination Universe, which stimulates the efficiency of the production of H$_2$ molecules and the formation of the first stellar objects. As a result, the onset of the formation of the first stars is shifted to higher redshifts, and the masses of the first stellar systems decrease. As a consequence, a sufficient increase in the ionizing radiation providing the reionization of the Universe can take place. We discuss possible observational manifestations of these effects and their dependence on the parameters of UHECR.

A simple analytical model for describing inner parts of dark matter halo is considered. It is assumed that dark matter density is power-law. The model deals with dark matter distribution function in phase space of adiabatic invariants (radial action and angular momentum). Two variants are considered for the angular part of the distribution function: narrow and broad distribution. The model allows to describe explicitly the process of adiabatic contraction of halo due to change of gravitational potential caused by condensation of baryonic matter in the centre. The modification of dark matter density in the centre is calculated, and is it shown that the standard algorithm of adiabatic contraction calculation overestimates the compressed halo density, especially in the case of strong radial anisotropy.

We investigate thermal regime of the baryons behind shock waves arising in the process of virialization of dark matter halos. We find a fraction of the shocked gas cooled by radiation of HD molecules down to the temperature of the cosmic microwave background (CMB): this fraction increases sharply from about $f_{\rm c}\sim 10^{-3}$ for dark halos of $M=5\times 10^7\msun$ to $\sim 0.1$ for halos with $M=10^8\msun$ at $z=10$. We show, however, that further increase of the mass does not lead to a significant growth of $f_{\rm c}$ -- the assymptotic value for $M\gg 10^8\msun$ is of 0.2. We estimate star formation rate associated with such shock waves, and show that it can be a small but not negligible fraction of the star formation connected with cooling by HI and H$_2$. We argue that extremely metal-poor low-mass stars in the Milky Way may have been formed from primordial gas behind such shocks.

HD molecules can be an important cooling agent of the primordial gas behind the shock waves originated through mergings of the dark matter haloes at epochs when first luminous objects were to form. We study the necessary conditions for the HD cooling to switch on in the low temperature range $T<200$ K. We show that these conditions are fulfiled in merging haloes with the total (dark matter and baryon) mass in excess of $M_{\rm cr}\sim 10^7[(1+z)/20]^{-2}\msun$. Haloes with masses $M>M_{\rm cr}$ may be the sites of low-mass star formation.

We study the effects of HD molecules on thermochemical evolution of the primordial gas behind shock waves, possibly arised in the process of galaxy formation. We find the critical shock velocity when deuterium transforms efficiently into HD molecules which then dominate gas cooling. Above this velocity the shocked gas is able to cool down to the temperature of the cosmic microwave background. Under these conditions the corresponding Jeans mass depends only on redshift and initial density of baryons $M_J \propto \delta_c^{-0.5} (1+z)^{0.5}$. At $z\simgt 45$ HD molecules heat shocked gas, and at larger redshift their contribution to thermal evolution becomes negligible.

The conditions for the fragmentation of the baryonic component during merging of dark matter halos in the early Universe are studied. We assume that the baryonic component undergoes a shock compression. The characteristic masses of protostellar molecular clouds and the minimum masses of protostars formed in these clouds decrease with increasing halo mass. This may indicate that the initial stellar mass function in more massive galaxies was shifted towards lower masses during the initial stages of their formation. This would result in an increase of the number of stars per unit halo mass, i.e., the efficiency of star formation.

The evolution of the dark matter distribution at the Galactic center is analyzed, which is caused by the combination of gravitational scattering on Galactic bulge stars and absorption by a supermassive black hole at the center of the bulge. Attention is focused on the boundary condition on the black hole. It is shown that its form depends on the energy of dark matter particles. The modified flux of dark matter particles onto the black hole is calculated. Estimates of the amount of dark matter absorbed show that the fraction of dark matter in the total mass of the black hole may be significant. The density of dark matter at the central part of the bulge is calculated. It is shown that recently observed gamma radiation from the Galactic center can be attributed to the annihilation of dark matter with this density.

The influence of dark matter (DM) on the growth of supermassive black holes (SMBHs) in galaxies is studied. It is shown that gravitational scattering of DM particles on bulge stars leads to diffusion of DM in phase space {m, m_z, I} ('m' denotes the angular momentum and 'I' is the radial action). Appropriate diffusion coefficients are calculated for different bulge models, and it is argued that the diffusion along 'm' axis is the most important effect. It is shown that this process leads to noticeable flow of DM into the black hole (BH), resulting in its power-law growth: M_{bh} ~ t^{9/16}. Comparison with observational data shows that, in principle, this effect may explain observed masses of SMBHs. Special attention is paid to the corrections related to the innermost region of BH gravitational influence and the diffusion along 'I' axis. Their influence on the BH growth law is shown to be negligible in most cases.

The presence of ultra-high energy cosmic rays in the universe in the pre-ionization epoch enhances significantly H$_2$ kinetics in virialized halos. It results in more efficient radiative cooling and decreases the lower mass limit of the first star-forming systems. Consequently, the fraction of baryons contained in the first luminous objects and their contribution to the reionization of the universe can significantly increase in comparison with the standard scenario.