Supermassive black holes (SMBHs) are thought to provide energy that prevents catastrophic cooling in the centers of massive galaxies and galaxy clusters. However, it remains unclear how this "feedback" process operates. We use high-resolution optical data to study the kinematics of multi-phase filamentary structures by measuring the velocity structure function (VSF) of the filaments over a wide range of scales in the centers of three nearby galaxy clusters: Perseus, Abell 2597 and Virgo. We find that the motions of the filaments are turbulent in all three clusters studied. There is a clear correlation between features of the VSFs and the sizes of bubbles inflated by SMBH driven jets. Our study demonstrates that SMBHs are the main driver of turbulent gas motions in the centers of galaxy clusters and suggests that this turbulence is an important channel for coupling feedback to the environment. Our measured amplitude of turbulence is in good agreement with Hitomi Doppler line broadening measurement and X-ray surface brightness fluctuation analysis, suggesting that the motion of the cold filaments is well-coupled to that of the hot gas. The smallest scales we probe are comparable to the mean free path in the intracluster medium (ICM). Our direct detection of turbulence on these scales provides the clearest evidence to date that isotropic viscosity is suppressed in the weakly-collisional, magnetized intracluster plasma.

We present three-dimensional, hybrid-kinetic numerical simulations of driven Alfv\'{e}n-wave turbulence of relevance to the collisionless near-Earth solar wind. Special attention is paid to the spectral transition that occurs near the ion-Larmor scale and to the origins of preferential perpendicular ion heating and of non-thermal wings in the parallel distribution function. Several novel diagnostics are used to show that the ion heating rate increases as the kinetic-Alfv\'{e}n-wave fluctuations, which comprise the majority of the sub-ion-Larmor turbulent cascade, attain near-ion-cyclotron frequencies. We find that ${\approx}75$-$80\%$ of the cascade energy goes into heating the ions, broadly consistent with the near-Earth solar wind. This heating is accompanied by clear velocity-space signatures in the particle energization rates and the distribution functions, including a flattened core in the perpendicular-velocity distribution and non-Maxwellian wings in the parallel-velocity distribution. The latter are attributed to transit-time damping and the pitch-angle scattering of perpendicularly heated particles into the parallel direction. Accompanying these features is a steepening of the spectral index of sub-ion-Larmor magnetic-field fluctuations beyond the canonical $-2.8$, as field energy is transferred to thermal energy. These predictions may be tested by measurements in the near-Earth solar wind.

Spiral density waves are known to exist in many astrophysical disks, potentially affecting disk structure and evolution. We conduct a numerical study of the effects produced by a density wave, evolving into a shock, on the characteristics of the underlying disk. We measure the deposition of angular momentum in the disk by spiral shocks of different strength and verify the analytical prediction of Rafikov (2016) for the behavior of this quantity, using shock amplitude (which is potentially observable) as the input variable. Good agreement between the theory and numerics is found as we vary shock amplitude (including highly nonlinear shocks), disk aspect ratio, equation of state, radial profiles of the background density and temperature, and pattern speed of the wave. We show that high numerical resolution is required to properly capture shock-driven transport, especially at low wave amplitudes. We also demonstrate that relating local mass accretion rate to shock dissipation in rapidly evolving disks requires accounting for the time-dependent contribution to the angular momentum budget, caused by the time dependence of the radial pressure support. We provide a simple analytical prescription for the behavior of this contribution and demonstrate its excellent agreement with the simulation results. Using these findings we formulate a theoretical framework for studying one-dimensional (in radius) evolution of the shock-mediated accretion disks, which can be applied to a variety of astrophysical systems.

The interaction of a satellite with a gaseous disk results in the excitation of spiral density waves which remove angular momentum from the orbit. In addition, if the orbit is not coplanar with the disk, three-dimensional effects will excite bending and eccentricity waves. We perform three-dimensional hydrodynamic simulations to study nonlinear disk-satellite interaction in inviscid protoplanetary disks for a variety of orbital inclinations from $0^\circ$ to $180^\circ$. It is well known that three-dimensional effects are important even for zero inclination. In this work we (1) show that for planets with small inclinations (as in the Solar system), effects such as the total torque and migration rate strongly depend on the inclination and are significantly different (about 2.5 times smaller) from the two-dimensional case, (2) give formulae for the migration rate, inclination damping, and precession rate of planets with different inclination angles in disk with different scale heights, and (3) present the observational signatures of a planet on an inclined orbit with respect to the protoplanetary disk. For misaligned planets we find good agreement with linear theory in the limit of small inclinations, and with dynamical friction estimates for intermediate inclinations. We find that in the latter case, the dynamical friction force is not parallel to the relative planetary velocity. Overall, the derived formulae will be important for studying exoplanets with obliquity.

We investigate the internal structure of the current sheet in the pulsar wind within force-free and two-fluid MHD approximations. Within the force-free approximation we obtain general asymptotic solution of the Grad-Shafranov equation for quasi-spherical pulsar wind up to the second order in small parameter $\varepsilon = (\Omega r/c)^{-1}$. The solution allows an arbitrary latitudinal structure of the radial magnetic field, including that obtained in the numerical simulations of oblique rotators. It is also shown that the shape of the current sheet does not depend on the latitudinal structure. For the internal region of the current sheet outside the fast magnetosonic surface where the force-free approximation is not valid we use two-fluid MHD approximation. Carrying out calculations in the comoving reference frame we succeed in determining intrinsic electric and magnetic fields of a sheet. It allows us to analyze time-dependent effects which were not investigated up to now. In particular, we estimate the efficiency of the particle acceleration inside the sheet. Finally, after investigating the motion of individual particles in the time-dependent current sheet, we find self-consistently the width of the sheet and its time evolution.

At present, there are theoretical models of radio pulsar evolution which predict both the alignment, i.e., evolution of inclination angle $\chi$ between magnetic and rotational axes to $0^{\circ}$, and its counter-alignment, i.e., evolution to $90^{\circ}$. At the same time, both models well describe the pulsar distribution on $P$-$\dot P$ diagram. For this reason, up to now it was impossible to determine the braking mechanisms since it was rather difficult to estimate inclination angle evolution on the basis of observation. In this paper we demonstrate that statistics of interpulse pulsars can give us the key to solve alignment/counter-alignment problem as the number of interpulse pulsars (both, having $\chi \sim 0^{\circ}$ and $\chi \sim 90^{\circ}$) drastically depends on evolution of inclination angle.

We investigate the motion of the primary beam outside the light cylinder in the pulsar wind. Inside the light cylinder both primary and secondary plasma move along dipole magnetic field lines where their energies can be arbitrary. But at larger distances the theory predicts quasi-radial motion with the velocity exactly corresponding to the drift velocity which cannot be the same for primary and secondary plasma. Hence, the deceleration of the primary beam is to take place simultaneously resulting in the acceleration of the secondary plasma. We investigate this process in the three-fluid MHD approximation and demonstrate that for most pulsars the energy of the beam remains practically unchanged. Only for young radio pulsars (Crab, Vela) essential deceleration up to the energy of the secondary plasma takes place outside the fast magnetosonic surface $r_{\rm F} \sim (10$-$100) R_{\rm L}$, the energy of secondary plasma itself increasing insufficiently.

Pulsars are famous for their rotational stability. Most of them steadily spin down and display a highly repetitive pulse shape. But some pulsars experience timing irregularities such as nulling, intermittency, mode changing and timing noise. As changes in the pulse shape are often correlated with timing irregularities, precession is a possible cause of these phenomena. Whereas pulsar magnetospheres are filled with plasma, most pulsar precession studies were carried out within the vacuum approximation and neglected the effects of magnetospheric currents and charges. Recent numerical simulations of plasma-filled pulsar magnetospheres provide us with a detailed quantitative description of magnetospheric torques exerted on the pulsar surface. In this paper, we present the study of neutron star evolution using these new torque expressions. We show that they lead to (1) much slower long-term evolution of pulsar parameters and (2) much less extreme solutions for these parameters than the vacuum magnetosphere models. To facilitate the interpretation of observed pulsar timing residuals, we derive an analytic model that (1) describes the time evolution of non-spherical pulsars and (2) translates the observed pulsar timing residuals into the geometrical parameters of the pulsar. We apply this model to two pulsars with very different temporal behaviours. For the pulsar B1828-11, we demonstrate that the timing residual curves allow two pulsar geometries: one with stellar deformation pointing along the magnetic axis and one along the rotational axis. For the Crab pulsar, we use the model show that the recent observation of its magnetic and rotational axes moving away from each other can be explained by precession.

Previously unexplored accretion regimes associated with the rotation of accreting matter, namely the perturbations of a quasi-spherical subsonic settling flow and Bondi-Hoyle accretion in the presence of axial rotation, are considered within the framework of ideal hydrodynamics. For subsonic settling accretion, the perturbations are shown to grow rapidly as the gravitating center is approached, so that the flow in the inner regions can no longer be considered as quasi-spherical. For Bondi-Hoyle accretion, a vacuum cylindrical cavity is shown to be formed at large distances from the gravitating center near the flow axis, with the flow velocity outside this cavity being virtually independent of the distance to the rotation axis.