Damping of the previously discovered resonant drag instability (RDI) of dust streaming in protoplanetary disc is studied using the local approach to dynamics of gas-dust perturbations in the limit of the small dust fraction. Turbulence in a disc is represented by the effective viscosity and diffusivity in equations of motion for gas and dust, respectively. In the standard case of the Schmidt number (ratio of the effective viscosity to diffusivity) Sc = 1, the reduced description of RDI in terms of the inertial wave (IW) and the streaming dust wave (SDW) falling in resonance with each other reveals that damping solution differs from the inviscid solution simply by adding the characteristic damping frequency to its growth rate. RDI is fully suppressed at the threshold viscosity, which is estimated analytically, first, for radial drift, next, for vertical settling of dust, and at last, in the case of settling combined with radial drift of the dust. In the last case, RDI survives up to the highest threshold viscosity, with a greater excess for smaller solids. Once Sc \neq 1, a new instability specific for dissipative perturbations on the dust settling background emerges. This instability of the quasi-resonant nature is referred to as settling viscous instability (SVI). The mode akin to SDW (IW) becomes growing in a region of long waves provided that Sc > 1 (Sc < 1). SVI leads to an additional increase of the threshold viscosity.

The ionospheric scattering of pulses emitted by PSR B0950+08 is measured using the 10-m RadioAstron Space Radio Telescope, the 300-m Arecibo Radio Telescope and the 14x25-m Westerbork Synthesis Radio Telescope (WSRT) at a frequency band between 316 and 332 MHz. We analyse this phenomenon based on a simulated model of the phase difference obtained between antennas that are widely separated by nearly 25 Earth diameters. We present a technique for processing and analysing the ionospheric total electron content (TEC) at the ground stations of the ground-space interferometer. This technique allows us to derive almost synchronous half-hour structures of the TEC in the ionosphere at an intercontinental distance between the Arecibo and WSRT stations. We find that the amplitude values of the detected structures are approximately twice as large as the values for the TEC derived in the International Reference Ionosphere (IRI) project. Furthermore, the values of the TEC outside these structures are almost the same as the corresponding values found by the IRI. According to a preliminary analysis, the detected structures were observed during a geomagnetic storm with a minimum Dst index of ~75 nT generated by interplanetary disturbances, and may be due to the influence of interplanetary and magnetospheric phenomena on ionospheric disturbances. We show that the Space Very Long Baseline Interferometry provides us with new opportunities to study the TEC, and we demonstrate the capabilities of this instrument to research the ionosphere.

The recently discovered resonant drag instability (RDI) of dust streaming in protoplanetary disc is considered as the mode coupling of subsonic gas-dust mixture perturbations. This mode coupling is coalescence of two modes with nearly equal phase velocities: inertial wave (IW) having positive energy and a streaming dust wave (SDW) having negative energy as measured in the frame of gas environment being at rest in vertical hydrostatic equilibrium. SDW is a trivial mode produced by the bulk streaming of dust, which transports perturbations of dust density. In this way, settling combined with radial drift of the dust makes possible coupling of SDW with IW and the onset of the instability. In accordance with the concept of the mode coupling, RDI growth rate is proportional to the square root of the coupling term of the dispersion equation, which itself is proportional to mass fraction of dust, $f\ll 1$. This clarifies why RDI growth rate $\propto f^{1/2}$. When SDW has positive energy, its resonance with IW provides an avoided crossing instead of the mode coupling. In the high wavenumber limit RDI with unbounded growth rate $\propto f^{1/3}$ is explained by the triple mode coupling, which is coupling of SDW with two IW. It coexists with a new quasi-resonant instability accompanied by bonding of two oppositely propagating low-frequency IW. The mode coupling does not exist for dust streaming only radially in a disc. In this case RDI is provided by the obscured mechanism associated with the inertia of solids.

The recently discovered resonant drag instability of dust settling in protoplanetary disc is considered as the mode coupling of subsonic gas-dust mixture perturbations. This mode coupling is coalescence of two modes with nearly equal phase velocities: the first mode is inertial wave having positive energy, while the second mode is a settling dust wave (SDW) having negative energy as measured in the frame of gas environment being at rest in vertical hydrostatic equilibrium. SDW is a trivial mode produced by the bulk settling of dust, which transports perturbations of dust density. The phase velocity of SDW is equal to the bulk settling velocity times the cosine of the angle formed by the wave vector and the rotation axis. In this way, the bulk settling of dust makes possible the coupling of SDW with the inertial wave and the onset of the instability. In accordance with the concept of the mode coupling, the instability growth rate is proportional to the square root of the dispersion equation coupling term, which itself contains the small mass fraction of dust in gas-dust mixture, the squared radial wavenumber of the modes, and the squared bulk settling velocity. Thus, the higher is the bulk settling velocity, the heavier clumps of dust can be aggregated by the instability of the same rate.

Planet migration originally refers to protoplanetary disks, which are more massive and dense than typical accretion disks in binary systems. We study planet migration in an accretion disk in a binary system consisting of a solar-like star hosting a planet and a red giant donor star. The accretion disk is fed by a stellar wind. %, disk self-gravity is neglected. We use the $\alpha$-disk model and consider that the stellar wind is time-dependent. Assuming the disk is quasi-stationary we calculate its temperature and surface density profiles. In addition to the standard disk model, when matter is captured by the disk at its outer edge, we study the situation when the stellar wind delivers matter on the whole disc surface inside the accretion radius with the rate depending on distance from the central star. Implying that a planet experiences classical type I/II migration we calculate migration time for a planet on a circular orbit coplanar with the disk. Potentially, rapid inward planet migration can result in a planet-star merger which can be accompanied by an optical or/and UV/X-ray transient. We calculate timescale of migration for different parameters of planets and binaries. Our results demonstrate that planets can fall on their host stars within the lifetime of the late-type donor for realistic sets of parameters.

We consider misaligned accretion discs formed after tidal disruption events occurring when a star encounters a supermassive rotating black hole. We use the linear theory of warped accretion discs to find the disc shape when the stream produced by the disrupted star provides a source of mass and angular momentum that is misaligned with the black hole. The evolution of the surface density and aspect ratio is found from a one dimensional vertically averaged model. We extend previous work which assumed a quasi-stationary disc to allow unrestricted dynamical propagation of disc tilt and twist through time dependent backgrounds. We consider a smaller value of the viscosity parameter, $\alpha =0.01,$ finding the dynamics varies significantly. At early times the disc inclination is found to be nearly uniform at small radii where the aspect ratio is large. However, since torques arise from the Lense-Thirring effect and the stream there is non uniform precession. We propose a simple model for this requiring only the background surface density and aspect ratio. At these times the $\alpha \sim 0.01$ disc exhibits a new feature. An inclined hot inner region joins an outer low inclination cool region via a thin transition front propagating outwards with a speed exceeding that of bending waves in the cool region. These waves accumulate where the propagation speeds match producing an inclination spike separating inner and outer discs. At late times a sequence of quasi-stationary configurations approximates disc shapes at small radii. We discuss observational implications of our results.

Enhanced angular momentum transfer through the boundary layer near the surface of weakly magnetised accreting star is required in order to explain the observed accretion timescales in low-mass X-ray binaries, cataclysmic variables or young stars with massive protoplanetary discs. Accretion disc boundary layer is locally represented by incompressible homogeneous and boundless flow of the cyclonic type, which is linearly stable. Its non-linear instability at the shear rates of order of the rotational frequency remains an issue. We argue that hydrodynamical subcritical turbulence in such a flow is sustained by the non-linear feedback from essentially three-dimensional vortices, which are generated by quasi-two-dimensional trailing shearing spirals grown to high amplitude via the swing amplification. We refer to those three-dimensional vortices as cross-rolls, since they are aligned in the shearwise direction in contrast to streamwise rolls generated by the anti-lift-up mechanism in rotating shear flow on the Rayleigh line. Transient growth of cross-rolls is studied analytically and further confronted with direct numerical simulations (DNS) of dynamics of non-linear erturbations. DNS performed in a tall box show that transition Reynolds number $R_T$ as function of shear rate accords with the line of constant maximum transient growth of cross-rolls. The transition in the tall box has been observed until the shear rate three times higher than the rotational frequency, when $R_T\sim 50000$. Assuming that the cross-rolls are also responsible for turbulence in the Keplerian flow, we estimate $R_T\lesssim 10^8$ in this case. The most favourable shear rate to test the existence of turbulence in the quasi-Keplerian regime may be sub-Keplerian and equal approximately to $1/2$.

Turbulent state of spectrally stable shear flows may be developed and sustained according to the bypass scenario of transition. If it works in non-magnetised boundless and homogeneous quasi-Keplerian flow, transiently growing shearing vortices should supply turbulence with energy. Employing the large shearing box approximation, as well as a set of global disc models, we study the optimal growth of the shearing vortices in such a flow in the whole range of azimuthal length-scales, $\lambda_y$, as compared to the flow scale-height, $H$. It is shown that with the account of the viscosity the highest possible amplification of shearing vortices, $G_{max}$, attains maximum at $\lambda_y\lesssim H$ and declines towards both the large scales $\lambda_y\gg H$ and the small scales $\lambda_y\ll H$ in a good agreement with analytical estimations based on balanced solutions. We pay main attention to the large-scale vortices $\lambda_y\gg H$, which produce $G_{max}\propto (\Omega/\kappa)^4$, where $\Omega$ and $\kappa$ denote local rotational and epicyclic frequencies, respectively. It is demonstrated that the large-scale vortices acquire high density perturbation as they approach the instant of swing. At the same time, their growth is not affected by bulk viscosity. We check that $G_{max}$ obtained globally is comparable to its local counterpart and the shape and localisation of global optimal vortices can be explained in terms of the local approach. The obtained results allow us to suggest that the critical Reynolds number of subcritical transition to turbulence in quasi-Keplerian flow, as well as the corresponding turbulent effective azimuthal stress should substantially depend on shear rate.

The paper is a brief overview of the works by Iosif S. Shklovsky (1916--1985), carried out over almost 30 years (1955--1985), on the nature of activity (primarily in the radio frequency range) in nuclei of some galaxies. Worthy of note is Shklovsky's pioneering work of 1962, in which he made an attempt to consider possible evolutionary tracks of extragalactic radio sources by constructing an analog of the Herzsprung--Russel diagram for stars (radio luminosity at 160 MHz was taken instead of optical luminosity; total radio size at the same frequency, as the other parameter). Later works by other authors are also discussed, where similar diagrams were plotted using a larger observational material. Special attention is paid to the evolution of Shklovsky's views regarding the possible ways of gas getting into radio galaxies' central regions, followed by high-velocity ejections of magnetized plasmons from their nuclei. In his assumptions, Shklovsky was mainly based on the observational data for the properties of the closest radio galaxy, NGC 4486 (Virgo A, M87), which he believed to be the same reference standard for extragalactic radio astronomy as the Crab Nebula for galactic radio astronomy. Shklovsky's approach to the recurrence of the activity phenomenon in galactic nuclei and the one-sided character of radio ejections from them is discussed. Modern views on these issues are also briefly considered.

We discovered fine-scale structure within the scattering disk of PSR B0329+54 in observations with the RadioAstron ground-space radio interferometer. Here, we describe this phenomenon, characterize it with averages and correlation functions, and interpret it as the result of decorrelation of the impulse-response function of interstellar scattering between the widely-separated antennas. This instrument included the 10-m Space Radio Telescope, the 110-m Green Bank Telescope, the 14x25-m Westerbork Synthesis Radio Telescope, and the 64-m Kalyazin Radio Telescope. The observations were performed at 324 MHz, on baselines of up to 235,000 km in November 2012 and January 2014. In the delay domain, on long baselines the interferometric visibility consists of many discrete spikes within a limited range of delays. On short baselines it consists of a sharp spike surrounded by lower spikes. The average envelope of correlations of the visibility function show two exponential scales, with characteristic delays of $\tau_1=4.1\pm 0.3\ \mu{\rm s}$ and $\tau_2=23\pm 3\ \mu{\rm s}$, indicating the presence of two scales of scattering in the interstellar medium. These two scales are present in the pulse-broadening function. The longer scale contains 0.38 times the scattered power of the shorter one. We suggest that the longer tail arises from highly-scattered paths, possibly from anisotropic scattering or from substructure at large angles.

This paper reviews some aspects of one of the major unsolved problems in understanding astrophysical (in particular, accretion) disks: whether the disk interiors may be effectively viscous in spite of the absence of marnetorotational instability? In this case a rotational homogeneous inviscid flow with a Keplerian angular velocity profile is spectrally stable, making the transient growth of perturbations a candidate mechanism for energy transfer from the regular motion to perturbations. Transient perturbations differ qualitatively from perturbation modes and can grow substantially in shear flows due to the nonnormality of their dynamical evolution operator. Since the eigenvectors of this operator, alias perturbation modes, are mutually nonorthogonal, they can mutually interfere, resulting in the transient growth of their linear combinations. Physically, a growing transient perturbation is a leading spiral whose branches are shrunk as a result of the differential rotation of the flow. This paper discusses in detail the transient growth of vortex shear harmonics in the spatially local limit as well as methods for identifying the optimal (fastest growth) perturbations. Special attention is given to obtaining such solutions variationally, by integrating the direct and adjoint equations forward and backward in time, respectively. The material is presented in a newcomer-friendly style.

We present not a literature review but a description, as detailed and consistent as possible, of two analytic models of disk accretion onto a rotating black hole: a standard relativistic disk and a twisted relativistic disk. Although one of these models is much older than the other, both are of topical current interest for black hole studies. The way the exposition is presented, the reader with only a limited knowledge of general relativity and relativistic hydrodynamics can --- with little or no use of additional sources -- gain good insight into many technical details lacking in the original papers.

We present the description of the FAKERAT software developed for planning Very Long Baseline Interferometry observations from space (space-VLBI). The results of the planned observations using the FAKERAT package during the first two years after launch of the space radio telescopy (SRT) in the space-ground interferometer modes are reported.

In this paper we introduce the first results that use data extracted directly from numerical simulations as inputs to the analytic twisted disk model of Zhuravlev & Ivanov. In both numerical and analytic approaches, fully relativistic models of tilted and twisted disks having a moderate effective viscosity around a slowly rotating Kerr black hole are considered. Qualitatively, the analytic model demonstrates the same dynamics as the simulations, although with some quantitative offset. Namely, the GRMHD simulations generally give smaller variations of tilt and twist across the disk. When the black hole and the disk rotate in the same sense, the simulated tilted disk and analytic model show no sign of Bardeen-Petterson alignment, even in the innermost parts of the disk, where the characteristic time for relaxation to a quasi-stationary configuration is of the same order as the computation time. In the opposite case, when the direction of the disk's rotation is opposite to that of black hole, a partial alignment is observed, in agreement with previous theoretical estimates. Thus, both fully numerical and analytic schemes demonstrate that the Bardeen-Petterson effect may not be possible for the case of prograde rotation provided that disk's effective viscosity is sufficiently small. This may have implications in modeling of different astrophysical phenomena, such as disk spectra and jet orientation.

This paper presents our latest numerical simulations of accretion disks that are misaligned with respect to the rotation axis of a Kerr black hole. In this work we use a new, fully conservative version of the Cosmos++ general relativistic magnetohydrodynamics (GRMHD) code, coupled with an ad hoc cooling function designed to control the thickness of the disk. Together these allow us to simulate the thinnest tilted accretion disks ever using a GRMHD code. In this way, we are able to probe the regime where the dimensionless stress and scale height of the disk become comparable. We present results for both prograde and retrograde cases. The simulated prograde tilted disk shows no sign of Bardeen-Petterson alignment even in the innermost parts of the disk. The simulated retrograde tilted disk, however, does show modest alignment. The implication of these results is that the parameter space associated with Bardeen-Petterson alignment for prograde disks may be rather small, only including very thin disks. Unlike our previous work, we find no evidence for standing shocks in our simulated tilted disks. We ascribe this to the combination of small black hole spin, small tilt angle, and small disk scale height in these simulations. We also add to the growing body of literature pointing out that the turbulence driven by the magnetorotational instability in global simulations of accretion disks is not isotropic. Finally, we provide a comparison between our moderately thin, untilted reference simulation and other numerical simulations of thin disks in the literature.

We study linear transient dynamics in a thin Keplerian disc employing a method based on variational formulation of optimisation problem. It is shown that in a shearing sheet approximation due to a prominent excitation of density waves by vortices the most rapidly growing shearing harmonic has azimuthal wavelength, $\lambda_y$, of order of the disc thickness, $H$, and its initial shape is always nearly identical to a vortex having the same potential vorticity. Also, in the limit $\lambda_y\gg H$ the optimal growth $G\propto (\Omega/\kappa)^4$, where $\Omega$ and $\kappa$ stand for local rotational and epicyclic frequencies, respectively, what suggests that transient growth of large scale vortices can be much stronger in areas with non-Keplerian rotation, e.g. in the inner parts of relativistic discs around the black holes. We estimate that if disc is already in a turbulent state with effective viscosity given by the Shakura parameter $\alpha<1$, the considered large scale vortices with wavelengths $H/\alpha>\lambda_y>H$ have the most favourable conditions to be transiently amplified before they are damped. At the same time, turbulence is a natural source of the potential vorticity for this transient activity. We extend our study to a global spatial scale showing that global perturbations with azimuthal wavelengths more than an order of magnitude greater than the disc thickness still are able to attain the growth of dozens of times in a few Keplerian periods at the inner boundary of disc.

RadioAstron space-ground VLBI observations of the pulsar B0950+08, conducted with the 10-m space radio telescope in conjunction with the Arecibo 300-m telescope and Westerbork Synthesis Radio Telescope at a frequency of 324 MHz, were analyzed in order to investigate plasma inhomogeneities in the direction of this nearby pulsar. The observations were conducted at a spacecraft distance of 330,000 km, resulting in a projected baseline of 220,000 km, providing the greatest angular resolution ever achieved at meter wavelengths. Our analysis is based on fundamental behavior of structure and coherence functions. We find that the pulsar shows scintillation on two frequency scales, both much less than the observing frequency; but modulation is less than 100%. We infer that the scattering is weak, but a refracting wedge disperses the scintillation pattern. The refraction angle of this "cosmic prism" is measured as theta_0=1.1 - 4.4 mas, with the refraction direction being approximately perpendicular to the observer velocity. We show that the observed parameters of scintillation effects indicate that two plasma layers lie along the line of sight to the pulsar, at distances of 4.4 - 16.4 pc and 26 - 170 pc, and traveling in different directions relative to the line of sight. Spectra of turbulence for the two layers are found to follow a power law with the indices gamma_1 = gamma_2 = 3.00 +/- 0.08, significantly different from the index expected for a Kolmogorov spectrum of turbulence, gamma=11/3.

We have studied the statistical and polarization properties of giant pulses (GPs) emitted by the millisecond pulsar B1937+21, with high sensitivity and time resolution. The observations were made in June 2005 with the 100-m Robert C. Byrd Green Bank Telescope at S-band (2052-2116 MHz) using the Mk5A VLBI recording system, with formal time resolution of 16 ns. The total observing time was about 4.5 hours; the rate of detection of GPs was about 130 per hour at the average longitudes of the main pulse (MPGPs) and 60 per hour at the interpulse (IPGPs). While the average profile shows well-defined polarization behavior, with regular evolution of the linear polarization position angle (PA), GPs exhibit random properties, occasionally having high linear or circular polarization. Neither MPGPs nor IPGPs show a preferred PA. The cumulative probability distribution (CPD) of GP pulse energy was constructed down to the level where GPs merge with regular pulses and noise. For both MPGPs and IPGPs, the CPD follows a power law with a break, the power index changing from -2.4 at high energy to -1.6 for low energy. Pulse smearing due to scattering masks the intrinsic shape and duration of the detected GPs. The smearing time varied during the observing session within a range of a few hundred nanoseconds. The measured polarization and statistical properties of GPs impose strong constraints on physical models of GPs. Some of these properties support a model in which GPs are generated by the electric discharge caused by magnetic reconnection of field lines connecting the opposite magnetic poles of a neutron star.

Results of analysis of long time variations of solar activities based on monthly Wolf number series (1749-2011 years) are represented. Shown the presence of stable oscillations in autocorrelation function for the Wolf numbers series with a period near 42.5 years. The statistical model of variations of autocorrelation function based on assumption of nonstationarity process are constructed. The question of using the model to predict the solar activity is discussed.

A thin gaseous disc with an almost keplerian angular velocity profile, bounded by a free surface and rotating around point-mass gravitating object is nearly spectrally stable. Despite that the substantial transient growth of linear perturbations measured by the evolution of their acoustic energy is possible. This fact is demonstrated for the simple model of a non-viscous polytropic thin disc of a finite radial size where the small adiabatic perturbations are considered as a linear combination of neutral modes with a corotational radius located beyond the outer boundary of the flow.

(abbreviated) In this paper we derive equations describing dynamics and stationary configurations of a twisted fully relativistic thin accretion disc around a slowly rotating black hole. We find that the disc dynamics and stationary shapes are determined by a pair of equations for two complex variables describing orientation of the disc rings and velocity perturbations in the disc. We analyse shapes of stationary twisted configurations. It is shown that the stationary configurations depend on two parameters - the $\alpha $parameter and $\tilde \delta = \delta_{*}/\sqrt a$, where $\delta_{*}\sim h/r$ is the disc opening angle (h is the disc halfthickness) and $a$ is the black hole rotational parameter. When $a > 0$ and $\tilde \delta \ll 1$ the shapes depend drastically on value of $\alpha$. When $\alpha $ is small the disc inclination angle oscillates with radius with amplitude and radial frequency of the oscillations dramatically increasing towards the last stable orbit. For moderate values of $\alpha $ the oscillations do not take place but the disc does not align with the equatorial plane at small radii. Its inclination angle is either increasing towards $R_{ms}$ or exhibits a non-monotonic dependence on the radial coordinate. Finally, when $\alpha $ is sufficiently large the disc aligns with the equatorial plane at small radii. When $a < 0$ the disc aligns with the equatorial plane for all values of $\alpha $. The results reported here may have implications for determining structure and variability of accretion discs close to $R_{ms}$ as well as for modelling of emission spectra coming from different sources, which are supposed to contain black holes.

The dynamics of small global perturbations in the form of linear combination of a finite number of non-axisymmetric eigenmodes is studied in two-dimensional approximation. The background flow is assumed to be an axisymmetric perfect fluid with the adiabatic index $\gamma=5/3$ rotating with power law angular velocity distribution $\Omega \propto r^{-q}$, $1.5<q<2.0$, confined by free boundaries in the radial direction. The substantial transient growth of acoustic energy of optimized perturbations is discovered. An optimal energy growth $G$ is calculated numerically for a variety of parameters. Its value depends essentially on the perturbation azimuthal wavenumber $m$ and increases for higher values of $m$. The closer the rotation profile to the Keplerian law, the larger growth factors can be obtained but over a longer time. The highest acoustic energy increase found numerically is of order $\sim 10^2$ over $\sim 6$ typical Keplerian periods. Slow neutral eigenmodes with corotation radius beyond the outer boundary mostly contribute to the transient growth. The revealed linear temporal behaviour of perturbations may play an important role in angular momentum transfer in toroidal flows near compact relativistic objects.

The instability of non-homoentropic axisymmetric flow of perfect fluid with respect to non-axisymmetric infinitesimal perturbations was investigated by numerical integration of hydrodynamical differential equations in two-dimensional approximation. The non-trivial influence of entropy gradient on unstable sound and surface gravity waves was revealed. In particular, both decrease and growth of entropy against the direction of effective gravitational acceleration $g_{eff}$ give rise to growing surface gravity modes which are stable with the same parameters in the case of homoentropic flow. At the same time increment of sound modes either grows monotonically while the rate of entropy decrease against $g_{eff}$ gets higher or vanishes at some values of positive and negative entropy gradient in the basic flow. The calculations have showed also that growing internal gravity modes appear only in the flow unstable to axisymmetric perturbations. At last, the analysis of boundary problem with free boundaries uncovered that's incorrect to set the entropy distribution according to polytropic law with polytropic index different from adiabatic value, since in this case perturbations don't satisfy the free boundary conditions.

We investigate the problem of self-gravitating dust disk dynamics in a static state taking into account nonlinear effects. For this purpose Schrodinger-type equation including the mass conservation law is used for the whole description of hydrodynamic flows of self-gravitating dust. We have shown a purely hydrodynamic mechanism of ring formation in the radial direction by taking into account nonlinearity in the lowest order of expansion parameter, which determines an order of magnitude of flow.