Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userQuasi-periodic oscillations as global hydrodynamic modes in the boundary layers of viscous accretion disks
82
0
0.0
(
0
)
Authors
M. Hakan Erkut
Ask ChatGPT about the research
No Arabic abstract
The observational characteristics of quasi-periodic oscillations (QPOs) from accreting neutron stars strongly indicate the oscillatory modes in the innermost regions of accretion disks as a likely source of the QPOs. The inner regions of accretion disks around neutron stars can harbor very high frequency modes related to the radial epicyclic frequency $kappa $. The degeneracy of $kappa $ with the orbital frequency $Omega $ is removed in a non-Keplerian boundary or transition zone near the magnetopause between the disk and the compact object. We show, by analyzing the global hydrodynamic modes of long wavelength in the boundary layers of viscous accretion disks, that the fastest growing mode frequencies are associated with frequency bands around $kappa $ and $kappa pm Omega $. The maximum growth rates are achieved near the radius where the orbital frequency $Omega $ is maximum. The global hydrodynamic parameters such as the surface density profile and the radial drift velocity determine which modes of free oscillations will grow at a given particular radius in the boundary layer. In accordance with the peak separation between kHz QPOs observed in neutron-star sources, the difference frequency between two consecutive bands of the fastest growing modes is always related to the spin frequency of the neutron star. This is a natural outcome of the boundary condition imposed by the rotating magnetosphere on the boundary region of the inner disk.
rate research
Read More
Disk accretion onto weakly magnetized objects possessing a material surface must proceed via the so-called boundary layer (BL) -- a region at the inner edge of the disk, in which the velocity of accreting material abruptly decreases from its Keplerian value. Supersonic shear arising in the BL is known to be conducive to excitation of acoustic waves that propagate into the accretor and the disk, enabling angular momentum and mass transport across the BL. We carry out a numerical exploration of different wave modes that operate near the BL, focusing on morphological characteristics of the modes in the innermost parts of accretion disk. Using a large suite of simulations covering a broad range of Mach numbers (of the supersonic shear flow in the BL), we provide accurate characterization of the different types of modes, verifying their properties against analytical results, when available. We discover new types of modes, in particular, global spiral density waves launched by vortices forming in the disk near the BL as a result of the Rossby wave instability; this instability is triggered by the vortensity production in that region caused by the nonlinear damping of acoustic waves. Azimuthal wavenumbers of the dominant modes that we observe appear to increase monotonically with the Mach number of the runs, but a particular mix of modes found in a simulation is mildly stochastic. Our results provide a basis for better understanding of the angular momentum and mass transport across the BL as well as the emission variability in accreting objects.
We use Smoothed Particle Hydrodynamics to study viscous accretion flows around a weakly magnetic neutron star. We show the formation of multiple ``boundary layers in presence of both cooling and viscosity. We find that with the introduction of a small viscosity in a sub-Keplerian flow, much like the wind accretion in HMXBs such as Cir X-1, only a single Normal Boundary Layer (NBOL) forms to adjust the rotational velocity component. With the increase of viscosity, the region extends radially and beyond some critical value, a RAdiative KEplerian Disk/layer (RAKED) forms between the sub-Keplerian flow and the NBOL. When viscosity is increased further only NBOL and RAKED remain. In all such cases, the CENtrifugal pressure dominated BOundary Layer (CENBOL) is formed, away from the star, as in the case of black holes. This is the first self-consistent study where such a transition from sub-Keplerian flows has been reported for neutron stars. We also identify the connection between accretion and ejection of matter, following the Two-Component Advective Flow for black holes, for neutron stars. The results are crucial in the understanding of the formation of disks, boundary layers and outflows in wind dominated neutron star systems.
When the accretion disc around a weakly magnetised neutron star (NS) meets the stellar surface, it should brake down to match the rotation of the NS, forming a boundary layer. As the mechanisms potentially responsible for this braking are apparently inefficient, it is reasonable to consider this layer as a spreading layer (SL) with negligible radial extent and structure. We perform hydrodynamical 2D spectral simulations of an SL, considering the disc as a source of matter and angular momentum. Interaction of new, rapidly rotating matter with the pre-existing, relatively slow material co-rotating with the star leads to instabilities capable of transferring angular momentum and creating variability on dynamical timescales. For small accretion rates, we find that the SL is unstable for heating instability that disrupts the initial latitudinal symmetry and produces large deviations between the two hemispheres. This instability also results in breaking of the axial symmetry as coherent flow structures are formed and escape from the SL intermittently. At enhanced accretion rates, the SL is prone to shearing instability and acts as a source of oblique waves that propagate towards the poles, leading to patterns that again break the axial symmetry. We compute artificial light curves of an SL viewed at different inclination angles. Most of the simulated light curves show oscillations at frequencies close to 1kHz. We interpret these oscillations as inertial modes excited by shear instabilities near the boundary of the SL. Their frequencies, dependence on flux, and amplitude variations can explain the high-frequency pair quasi-periodic oscillations observed in many low-mass X-ray binaries.
Recently, the vertical shear instability (VSI) has become an attractive purely hydrodynamic candidate for the anomalous angular momentum transport required for weakly ionized accretion disks. In direct three-dimensional numerical simulations of VSI turbulence in disks, a meridional circulation pattern was observed that is opposite to the usual viscous flow behavior. Here, we investigate whether this feature can possibly be explained by an anisotropy of the VSI turbulence. Using three-dimensional hydrodynamical simulations, we calculate the turbulent Reynolds stresses relevant for angular momentum transport for a representative section of a disk. We find that the vertical stress is significantly stronger than the radial stress. Using our results in viscous disk simulations with different viscosity coefficients for the radial and vertical direction, we find good agreement with the VSI turbulence for the stresses and meridional flow; this provides additional evidence for the anisotropy. The results are important with respect to the transport of small embedded particles in disks.
Axisymmetric magnetorotational instability (MRI) in viscous accretion disks is investigated by linear analysis and two-dimensional nonlinear simulations. The linear growth of the viscous MRI is characterized by the Reynolds number defined as $R_{rm MRI} equiv v_A^2/ uOmega $, where $v_A$ is the Alfv{e}n velocity, $ u$ is the kinematic viscosity, and $Omega$ is the angular velocity of the disk. Although the linear growth rate is suppressed considerably as the Reynolds number decreases, the nonlinear behavior is found to be almost independent of $R_{rm MRI}$. At the nonlinear evolutionary stage, a two-channel flow continues growing and the Maxwell stress increases until the end of calculations even though the Reynolds number is much smaller than unity. A large portion of the injected energy to the system is converted to the magnetic energy. The gain rate of the thermal energy, on the other hand, is found to be much larger than the viscous heating rate. Nonlinear behavior of the MRI in the viscous regime and its difference from that in the highly resistive regime can be explained schematically by using the characteristics of the linear dispersion relation. Applying our results to the case with both the viscosity and resistivity, it is anticipated that the critical value of the Lundquist number $S_{rm MRI} equiv v_A^2/etaOmega$ for active turbulence depends on the magnetic Prandtl number $S_{{rm MRI},c} propto Pm^{1/2}$ in the regime of $Pm gg 1$ and remains constant when $Pm ll 1$, where $Pm equiv S_{rm MRI}/R_{rm MRI} = u/eta$ and $eta$ is the magnetic diffusivity.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
