Theoretical astrophysics and plasma physics at RPC

Abstract:

We derive an analytical model for the so-called phenomenon of resonant dynamical friction, where a disc of stars around a supermassive black hole interacts with a massive perturber, so as to align its inclination with the disc’s orientation. We show that it stems from a singular behaviour of the orbit-averaged equations of motion, which leads to a rapid alignment of the argument of the ascending node of each of the disc stars, with that of the perturber, p, with a phase difference of 90◦. This phenomenon occurs for all stars whose maximum possible ˙ (maximized over all values of for all the disc stars) is greater than ˙ p; this corresponds approximately to all stars whose semi-major axes are less than twice that of the perturber. The rate at which the perturber’s inclination decreases with time is proportional to its mass and is shown to be much faster than Chandrasekhar’s dynamical friction. We find that the total alignment time is inversely proportional to the root of the perturber’s mass. This persists until the perturber enters the disc. The predictions of this model agree with a suite of numerical N-body simulations, which we perform to explore this phenomenon, for a wide range of initial conditions, masses, etc., and are an instance of a general phenomenon. Similar effects could occur in the context of planetary systems, too.

Abstract:

The equations of electrostatic drift kinetics are observed to possess a symmetry associated with their intrinsic scale invariance. Under the assumptions of spatial periodicity, stationarity, and locality, this symmetry implies a particular scaling of the turbulent heat flux with the system's parallel size, from which its scaling with the equilibrium temperature gradient can be deduced under some additional assumptions. This macroscopic transport prediction is then confirmed numerically for a reduced model of electron-temperature-gradient-driven turbulence in slab geometry. The system realises this scaling through a turbulent cascade from large to small perpendicular spatial scales. The route of this cascade through wavenumber space (i.e. the relationship between parallel and perpendicular scales in the inertial range) is shown to be determined by a balance between nonlinear-decorrelation and parallel-dissipation timescales. This type of ‘critically balanced’ cascade, which maintains a constant energy flux despite the presence of parallel dissipation throughout the inertial range (as well as order-unity dissipative losses at the outer scale) is expected to be a generic feature of plasma turbulence. The outer scale of the turbulence, on which the turbulent heat flux depends, is determined by the breaking of drift-kinetic scale invariance due to the existence of large-scale parallel inhomogeneity (the parallel system size).

Abstract:

This paper examines the energetics of a convective flow subject to an oscillation with a period $t_{\rm osc}$ much smaller than the convective time-scale $t_{\rm conv}$, allowing for compressibility and uniform rotation. We show that the energy of the oscillation is exchanged with the kinetic energy of the convective flow at a rate $D_R$ that couples the Reynolds stress of the oscillation with the convective velocity gradient. For the equilibrium tide and inertial waves, this is the only energy exchange term, whereas for p modes there are also exchanges with the potential and internal energy of the convective flow. Locally, $\left| D_R \right| \sim u^{\prime 2} / t_{\rm conv}$, where $u^{\prime}$ is the oscillating velocity. If $t_{\rm conv} \ll t_{\rm osc}$ and assuming mixing length theory, $\left| D_R \right|$ is $\left( \lambda_{\rm conv} / \lambda_{\rm osc} \right)^2$ smaller, where $\lambda_{\rm conv}$ and $\lambda_{\rm osc}$ are the characteristic scales of convection and the oscillation. Assuming local dissipation, we show that the equilibrium tide lags behind the tidal potential by a phase $\delta(r) \sim r \omega_{\rm osc} / \left( g(r) t_{\rm conv}(r) \right)$, where g is the gravitational acceleration. The equilibrium tide can be described locally as a harmonic oscillator with natural frequency $\left( g/r \right)^{1/2}$ and subject to a damping force $-u^{\prime}/t_{\rm conv}$. Although $\delta(r)$ varies by orders of magnitude through the flow, it is possible to define an average phase shift $\overline{\delta }$ which is in good agreement with observations for Jupiter and some of the moons of Saturn. Finally, $1 / \overline{\delta }$ is shown to be equal to the standard tidal dissipation factor.

Abstract:

Motivated by the increasing number of detections of merging black holes by LIGO-VIRGO-KAGRA, black hole (BH) binary mergers in the discs of active galactic nuclei (AGNs) is investigated as a possible merger channel. In this pathway, BH encounters in the gas disc form mutually bound BH binary systems through interaction with the gas in the disc and subsequently inspiral through gravitational torques induced by the local gas. To determine the feasibility of this merger pathway, we present the first three-dimensional global hydrodynamic simulations of the formation and evolution of a stellar-mass BH binaries AGN discs with three different AGN disc masses and five different initial radial separations. These 15 simulations show binary capture of prograde and retrograde binaries can be successful in a range of disc densities including cases well below that of a standard radiatively efficient alpha disc, identifying that the majority of these captured binaries are then subsequently hardened by the surrounding gas. The eccentricity evolution depends strongly on the orbital rotation where prograde binaries are governed by gravitational torques form their circumbinary mini disc, with eccentricities being damped, while for retrograde binaries the eccentricities are excited to >∼ 0.9 by accretion torques. In two cases, retrograde binaries ultimately undergo a close periapsis passage which results in a merger via gravitational waves after only a few thousand binary orbits. Thus, the merger time-scale can be far shorter than the AGN disc lifetime. These simulations support an efficient AGN disc merger pathway for BHs.