Keplerian Accretion Disks Research Articles

Может ли нейтронная звезда перейти из состояния эжектора в состояние аккретора, минуя состояние пропеллера? I. Приближение немагнитного аккреционного потока

We discuss evolution of a neutron star in the massive binary system in terms of evolutionary chain of its states: ejector, propeller and accretor. We show that a direct transition between the states of ejector and accretor (without the propeller state) can take place only in the systems with relatively low-velocity stellar wind of the massive component in case of formation of the Keplerian accretion disk of low viscosity.

A probable Keplerian disk feeding an optically revealed massive young star.

The canonical picture of star formation involves disk-mediated accretion, with Keplerian accretion disks and associated bipolar jets primarily observed in nearby, low-mass young stellar objects (YSOs). Recently, rotating gaseous structures and Keplerian disks have been detected around several massive (M > 8 M⊙) YSOs (MYSOs)1-4, including several disk-jet systems5-7. All the known MYSO systems are in the Milky Way, and all are embedded in their natal material. Here we report the detection of a rotating gaseous structure around an extragalactic MYSO in the Large Magellanic Cloud. The gas motion indicates that there is a radial flow of material falling from larger scales onto a central disk-like structure. The latter exhibits signs of Keplerian rotation, so that there is a rotating toroid feeding an accretion disk and thus the growth of the central star. The system is in almost all aspects comparable to Milky Way high-mass YSOs accreting gas from a Keplerian disk. The key difference between this source and its Galactic counterparts is that it is optically revealed rather than being deeply embedded in its natal material as is expected of such a massive young star. We suggest that this is the consequence of the star having formed in a low-metallicity and low-dust content environment. Thus, these results provide important constraints for models of the formation and evolution of massive stars and their circumstellar disks.

Spin Equilibrium of Rapidly Spinning Neutron Stars via Transient Accretion

The concept of spin equilibrium due to an interaction between the stellar magnetosphere and a thin, Keplerian accretion disk, and a well-known formula of the corresponding equilibrium spin frequency, provide a key understanding of spin evolution and the distribution of rapidly spinning neutron stars, viz., millisecond pulsars. However, this concept and formula are for stable accretion, but the mass transfer to most accreting millisecond pulsars is transient and the accretion rate evolves by orders of magnitude during an outburst. In this short and focussed review, we briefly discuss a relatively new concept of the spin equilibrium condition and a new formula for the equilibrium spin frequency for transiently accreting millisecond pulsars. We also review a new method to estimate this equilibrium spin frequency for observed transiently accreting millisecond pulsars, even when a pulsar has not yet attained the spin equilibrium. These will be crucial to probe the spin evolution and distribution of millisecond pulsars, and should also be applicable to all magnetic stars transiently accreting via a thin, Keplerian accretion disk.

Магниторотационная неустойчивость в кеплеровских дисках: нелокальный подход

Магниторотационная неустойчивость в кеплеровских дисках: нелокальный подход

Магниторотационная неустойчивость в кеплеровских дисках: нелокальный подход

Магниторотационная неустойчивость в кеплеровских дисках: нелокальный подход

A discontinuity at the base of the transition layer located between the Keplerian accretion disk and the compact object

Context.We study the geometry of the transition layer (TL) between the classical Keplerian accretion disk (the TL outer boundary) and the compact object at the TL inner boundary.Aims.Our goal is to use the hydrodynamical formalism to demonstrate that the TL is created, together with a shock, in response to a discontinuity and to adjust the Keplerian disk motion to a central object (CO).Methods.We apply hydrodynamical equations to describe a plasma motion near a CO in the TL.Results.We point out that, before matter accretes to a CO, the TL cloud is formed between an adjustment radius and the TL inner boundary, which is probably a site where the emergent Compton spectrum originates. Using a generalization of the Randkine–Hugoniot relation and a solution of the azimutal force balance equation, we are able to reproduce the geometric characteristics of the TL.

Торможение аккреционного диска магнитным полем аккретора в приближении коротации

We analyze the process of interaction of the viscous keplerian accretion disk with the dipole magnetic field of the accretor within the corotation approximation. We estimate the radial dependence of the gaseous pressure in such non-magnetized accretion disk. We calculate the equilibrium radius at which the gaseous pressure in the disk is balanced by the pressure of the dipole magnetic field of the accreting star. We show that this radius depends on the disk viscosity and is significantly smaller than the canonical Alfvén radius realized in the process of spherical accretion.

Charged particle dynamics in parabolic magnetosphere around Schwarzschild black hole

The study of charged particle dynamics in the combined gravitational and magnetic field can provide important theoretical insight into astrophysical processes around black holes. In this paper, we explore the charged particle dynamics in parabolic magnetic field configuration around Schwarzschild black hole, since the paraboloidal shapes of magnetic field lines around black holes are well motivated by the numerical simulations and supported by observations of relativistic jets. Analysing the stability of bounded orbits and using the effective potential approach, we show the possibility of existence of stable circular off-equatorial orbits around the symmetry axis. We also show the influence of radiation reaction force on the dynamics of charged particles, in particular on the chaoticity of the motion and Poincaré sections, oscillatory frequencies, and emitted electromagnetic spectrum. Applied to Keplerian accretion disks, we show that in parabolic magnetic field configuration, the thin accretion configurations can be either destroyed or transformed into a thick toroidal structure given the radiation reaction and electromagnetic-disk interactions included. Calculating the Fourier spectra for radiating charged particle trajectories, we find that the radiation reaction force does not affect the main frequency peaks, however, it lowers the higher harmonics making the spectrum more flat and diluted in high frequency range.

Numerical simulations of MHD jets from Keplerian accretion disks

Context. The most successful scenario for the origin of astrophysical jets requires a large-scale magnetic field anchored in a rotating object (black hole or star) and/or its surrounding accretion disk. Platform jet simulations, where the mass load onto the magnetic field is not computed by solving the vertical equilibrium of the disk but is imposed as a boundary condition, are very useful for probing the jet acceleration and collimation mechanisms. The drawback of such simulations is the very large parameter space: despite many previous attempts, it is very difficult to determine the generic results that can be derived from them. Aims. We wish to establish a firm link between jet simulations and analytical studies of magnetically driven steady-state jets from Keplerian accretion disks. In particular, the latter have predicted the existence of recollimation shocks – due to the dominant hoop stress –, which have so far never been observed in platform simulations. Methods. We performed a set of axisymmetric magnetohydrodynamics (MHD) simulations of nonrelativistic jets using the PLUTO code. The simulations are designed to reproduce the boundary conditions generally expected in analytical studies. We vary two parameters: the magnetic flux radial exponent α and the jet mass load κ. In order to reach the huge unprecedented spatial scales implied by the analytical solutions, we used a new method allowing us to boost the temporal evolution. Results. We confirm the existence of standing recollimation shocks at large distances. As in self-similar studies, their altitude evolves with the mass load κ. The shocks are weak and correspond to oblique shocks in a moderately high, fast magnetosonic flow. The jet emitted from the disk is focused toward the inner axial spine, which is the outflow connected to the central object. The presence of this spine is shown to have a strong influence on jet asymptotics. We also argue that steady-state solutions with α ≥ 1 are numerically out of range. Conclusions. Internal recollimation shocks may produce observable features such as standing knots of enhanced emission and a decrease in the flow rotation rate. However, more realistic simulations (e.g. fully three-dimensional) must be carried out in order to investigate nonaxisymmetric instabilities and with ejection only from a finite zone in the disk, so as to to verify whether these MHD recollimation shocks and their properties are maintained.

Geodesic Model of HF QPOs Tested for Black Holes in Spacetimes Reflecting the Effect of Surrounding Dark Matter

Using the simple but robust model of a shell of dark matter (DM) around a Schwarzschild black hole (BH), represented by the mass ratio of the shell and BH ΔM/M, the shell extension Δr s and its inner radius r s, we study the influence of DM on the spacetime structure and geodesic motion, and provide a classification of the BH+DM shell spacetimes according to the properties of the stable circular geodesics governing Keplerian disks. We focus our attention on the epicyclic motion around circular geodesics that can be related to observational phenomena in X-ray radiation from Keplerian accretion disks, assumed to be influenced by the DM shell only gravitationally. We give the frequencies of the orbital and epicyclic motions and discuss their properties in terms of the parameters governing the DM shell. Using the frequencies in relevant variants of the standard geodesic model of high-frequency quasiperiodic oscillations (HF QPOs), we test the role of DM by fitting the HF QPO data from some microquasars and active galactic nuclei with supermassive BHs where no variant of the geodesic model applied in the standard vacuum BH background is able to explain the data. We thus provide a robust review of the applicability of the geodesic model of HF QPOs, and also provide limits on the amount of DM around a BH. We demonstrate that the geodesic model could be well applied to most observations of active galactic nuclei, with strong restrictions on the amount of invisible matter around BHs.

Flux and spectral variability of Mrk 421 during its moderate activity state usingNuSTAR: Possible accretion disc contribution?

Context.The X-ray emission in BL Lac objects is believed to be dominated by synchrotron emission from their relativistic jets. However, when the jet emission is not strong, one could expect signatures of X-ray emission from inverse Compton scattering of accretion disc photons by hot and energetic electrons in the corona. Moreover, the observed X-ray variability can also originate in the disc, and get propagated and amplified by the jet.Aims.Here, we present results on the BL Lac object Mrk 421 using the Nuclear Spectroscopic Telescope Array data acquired during 2017 when the source was in a moderate X-ray brightness state. For comparison with high jet activity state, we also considered one epoch of data in April 2013 when the source was in a very high X-ray brightness state. Our aim is to explore the possibility of the signature of accretion disc emission in the overall X-ray emission from Mrk 421 and also examine changes in accretion parameters considering their contribution to spectral variations.Methods.We divided each epoch of data into different segments in order to find small-scale variability. Data for all segments were fitted using a simple power-law model. We also fitted the full epoch data using the two component advective flow (TCAF) model to extract the accretion flow parameters. Furthermore, we estimated the X-ray flux coming from the different components of the flow using the lowest normalisation method and analysed the relations between them. For consistency, we performed the spectral analysis using models available in the literature.Results.The simple power-law function does not fit the spectra well, and a cutoff needs to be added. The spectral fitting of the data using the TCAF model shows that the data can be explained with a model where (a) the size of the dynamic corona at the base of the jet is from ∼28 to 10 rs, (b) the disc mass accretion rate is from 0.021 to 0.051ṀEdd, (c) the halo mass accretion rate is from 0.22 to 0.35ṀEdd, and (d) the viscosity parameter of the Keplerian accretion disc from 0.18 to 0.25. In the assumed model, the total flux, disc and jet flux correlate with the radio flux observed during these epochs.Conclusions.From the spectral analysis, we conclude that the spectra of all the epochs of Mrk 421 in 2017 are well described by the accretion-disc-based TCAF model. The estimated disc and jet flux relations with radio flux show that accretion disc can contribute to the observed X-ray emission, when X-ray data (that cover a small portion of the broad band spectral energy distribution of Mrk 421) are considered in isolation. However, the present disc-based models are disfavoured with respect to the relativistic jet models when considering the X-ray data in conjunction with data at other wavelengths.

Evolution of braneworld Kerr-Newman naked singularities

We study evolution of the braneworld Kerr--Newman (K-N) naked singularities, namely their mass $M$ , spin $a$, and tidal charge $b$ characterizing the role of the bulk space, due to matter in-falling from Keplerian accretion disk. We construct the evolution in two limiting cases applied to the tidal charge. In the first case we assume $b$ = const during the evolution, in the second one we assume that the dimensionless tidal charge $\beta \equiv b/M^2$ = const. For positive values of the tidal charge the evolution is equivalent to the case of the standard K-N naked singularity under accretion of electrically neutral matter. We demonstrate that counter-rotating accretion always converts a K-N naked singularity into an extreme K-N black hole and that the corotating accretion leads to variety of outcomes. The conversion to an extreme K-N black hole is possible for naked singularity with dimensionless tidal charge $\beta < 0.25$, and $\beta \in (0.25, 1)$ with sufficiently low spin. In other cases the accretion ends in a transcendental state. For $0.25 < \beta < 1$ this is a mining unstable K-N naked singularity enabling formally unlimited energy extraction from the naked singularity. In the case of $\beta > 1$, the corotating accretion creates unlimited torodial structure of mater orbiting the naked singularity. Both non-standard outcomes of the corotating accretion imply a transcendence of such naked singularity due to nonlinear gravitational effects.

Spectroscopic observations of V455 Andromedae superoutburst in 2007: The most exotic spectral features in dwarf nova outbursts

Abstract We present our spectroscopic observations of V455 Andromedae during the 2007 superoutburst. Our observations cover this superoutburst from around the optical peak of the outburst to the post-superoutburst stage. During the early superhump phase, the emission lines of the Balmer series, He i, He ii, Bowen blend, and C iv/N iv blend were detected. The He ii 4686 line exhibited a double-peaked emission profile, where Balmer emission lines were single-peaked, which is unexpected from its high inclination. In the ordinary superhump phase, the Balmer series transitioned to double-peaked emission profiles, and high-ionization lines were significantly weakened. These transitions of the line profiles should be related to the structural transformation of the accretion disk, as expected between the early and ordinary superhump transition in the thermal–tidal instability model. The Doppler map of Hα during the early superhump phase exhibits a compact blob centered at the primary white dwarf. In analogy to SW Sex-type cataclysmic variables, this feature could emerge from the disk wind and/or the mass accretion column on to the magnetized white dwarf. The Doppler map of He ii 4686 Å is dominated by the ring-like structure and imposed two flaring regions with the velocity of ∼300 km s−1, which is too slow for a Keplerian accretion disk. The phase of the flaring regions was coincident with the inner spiral arm structure identified during the early superhump phase. Our disk wind model with the enhanced emission from the wind component launched from the inner arm structure successfully reproduced the observed properties of He ii 4686 Å. Therefore, V455 And is the first case in dwarf nova outbursts where the presence of the disk wind is inferred from an optical spectrum.

LkHα 225 (V1318 Cyg) South in Outburst

Abstract Magakian et al. called attention to the current bright state of LkHα 225 South, a well-known highly embedded, intermediate-mass young stellar object that over the past two decades has brightened visually from &gt;20 m to &lt;13 m . We present recent optical photometric monitoring showing colorless, nonsinusoidal, periodic brightness oscillations occurring every 43 days with amplitude ∼0.7 mag. We also present new flux-calibrated optical and near-infrared spectroscopy, which we model in terms of a Keplerian accretion disk, and high-dispersion spectra that demonstrate similarity to some categories of “mixed-temperature” accretion-outburst objects. At blue wavelengths, LkHα 225 South has a pure absorption spectrum and is a good spectral match to the FU Ori stars V1515 Cyg and V1057 Cyg. At red optical and infrared wavelengths, however, the spectrum is more similar to Gaia 19ajj, showing emission in TiO, CO, and metals. Sr ii absorption indicates a low-surface-gravity atmosphere. There are also signatures of a strong wind/outflow. LkHα 225 South was moderately bright in the early 1950s as well as in the late 1980s, with evidence for deep fades during intervening epochs. The body of evidence suggests that LkHα 225 South is another case of a source with episodically enhanced accretion that causes brightening by orders of magnitude, and development of a hot absorption spectrum and warm wind. It is similar to Gaia 19ajj, but also reminiscent in its long brightening time and brightness oscillation near peak to the embedded sources L1634 IRS7 and ESO Ha 99.

OSIRIS: a new code for ray tracing around compact objects

The radiation observed in quasars and active galactic nuclei is mainly produced by a relativistic plasma orbiting close to the black hole event horizon, where strong gravitational effects are relevant. The observational data of such systems can be compared with theoretical models to infer the black hole and plasma properties. In the comparison process, ray-tracing algorithms are essential to computing the trajectories followed by the photons from the source to our telescopes. In this paper, we present OSIRIS: a new stable FORTRAN code capable of efficiently computing null-geodesics around compact objects, including general relativistic effects such as gravitational lensing, redshift, and relativistic boosting. The algorithm is based on the Hamiltonian formulation and uses different integration schemes to evolve null-geodesics while tracking the error in the Hamiltonian constrain to ensure physical results. We found from an error analysis that the integration schemes are all stable, and the best one maintains an error below 10^{-11}. Particularly, to test the robustness and ability of the code to evolve geodesics in curved space-time, we compute the shadow and Einstein rings of a Kerr black hole with different rotation parameters and obtain the image of a thin Keplerian accretion disk around a Schwarzschild black hole. Although OSIRIS is parallelized neither with MPI nor with CUDA, the computation times are of the same order as those reported by other codes with these types of parallel computing platforms.

Forced Linear Shear Flows with Rotation: Rotating Couette–Poiseuille Flow, Its Stability, and Astrophysical Implications

Abstract We explore the effect of forcing on the linear shear flow or plane Couette flow, which is also the background flow in the very small region of the Keplerian accretion disk. We show that depending on the strength of forcing and boundary conditions suitable for the systems under consideration, the background plane shear flow, and hence the accretion disk velocity profile, is modified into parabolic flow, which is a plane Poiseuille flow or Couette–Poiseuille flow, depending on the frame of reference. In the presence of rotation, the plane Poiseuille flow becomes unstable at a smaller Reynolds number under pure vertical as well as three-dimensional perturbations. Hence, while rotation stabilizes the plane Couette flow, the same destabilizes the plane Poiseuille flow faster and hence the forced local accretion disk. Depending on the various factors, when the local linear shear flow becomes a Poiseuille flow in the shearing box due to the presence of extra force, the flow becomes unstable even for Keplerian rotation, and hence turbulence will ensue. This helps to resolve the long-standing problem of subcritical transition to turbulence in hydrodynamic accretion disks and the laboratory plane Couette flow.

A slim disc approach to external photoevaporation of discs

ABSTRACT The photoevaporation of protoplanetary discs by nearby massive stars present in their birth cluster plays a vital role in their evolution. Previous modelling assumes that the disc behaves like a classical Keplerian accretion disc out to a radius where the photoevaporative outflow is launched. There is then an abrupt change in the angular velocity profile, and the outflow is modelled by forcing the fluid parcels to conserve their specific angular momenta. Instead, we model externally photoevaporating discs using the slim disc formalism. The slim disc approach self-consistently includes the advection of radial and angular momentum as well as angular momentum redistribution by internal viscous torques. Our resulting models produce a smooth transition from a rotationally supported Keplerian disc to a photoevaporative driven outflow, where this transition typically occurs over ∼4–5 scale heights. The penetration of ultraviolet photons predominately sets the radius of the transition and the viscosity’s strength plays a minor role. By studying the entrainment of dust particles in the outflow, we find a rapid change in the dust size and surface density distribution in the transition region due to the steep gas density gradients present. This rapid change in the dust properties leaves a potentially observable signature in the continuum spectral index of the disc at mm wavelengths. Using the slim disc formalism in future evolutionary calculations will reveal how both the gas and dust evolve in their outer regions and the observable imprints of the external photoevaporation process.

Origin of hydrodynamic instability from noise: From laboratory flow to accretion disk

We attempt to address the old problem of plane shear flows: the origin of turbulence and hence transport of angular momentum in accretion flows as well as laboratory flows, such as plane Couette flow. We undertake the problem by introducing an extra force in Orr-Sommerfeld and Squire equations along with the Coriolis force mimicking the local region of the accretion disk. For plane Couette flow, the Coriolis term drops. Subsequently we solve the equations by WKB approximation method. We investigate the dispersion relation for the Keplerian flow and plane Couette flow for all possible combinations of wave vectors. Due to the very presence of extra force, we show that both the flows are unstable for a certain range of wave vectors. However, the nature of instability between the flows is different. We also study the Argand diagrams of the perturbation eigenmodes. It helps us to compare the different time scales corresponding to the perturbations as well as accretion. We ultimately conclude with this formalism that fluid gets enough time to be unstable and hence plausibly turbulent particularly in the local regime of the Keplerian accretion disks. Repetition of the analysis throughout the disk explains the transport of angular momentum and matter along outward and inward direction respectively.

Hydrodynamical instability with noise in the Keplerian accretion discs: modified Landau equation

ABSTRACT Origin of hydrodynamical instability and turbulence in the Keplerian accretion disc as well as similar laboratory shear flows, e.g. plane Couette flow, is a long-standing puzzle. These flows are linearly stable. Here we explore the evolution of perturbation in such flows in the presence of an additional force. Such a force, which is expected to be stochastic in nature hence behaving as noise, could be result of thermal fluctuations (however small be), Brownian ratchet, grain–fluid interactions, feedback from outflows in astrophysical discs, etc. We essentially establish the evolution of nonlinear perturbation in the presence of Coriolis and external forces, which is modified Landau equation. We show that even in the linear regime, under suitable forcing and Reynolds number, the otherwise least stable perturbation evolves to a very large saturated amplitude, leading to nonlinearity and plausible turbulence. Hence, forcing essentially leads a linear stable mode to unstable. We further show that nonlinear perturbation diverges at a shorter time-scale in the presence of force, leading to a fast transition to turbulence. Interestingly, emergence of nonlinearity depends only on the force but not on the initial amplitude of perturbation, unlike original Landau equation based solution.

The Rise and Fall of the King: The Correlation between FO Aquarii’s Low States and the White Dwarf’s Spin-down

Abstract The intermediate polar FO Aquarii (FO Aqr) experienced its first reported low-accretion states in 2016, 2017, and 2018, and, using newly available photographic plates, we identify prediscovery low states in 1965, 1966, and 1974. The primary focus of our analysis, however, is an extensive set of time-series photometry obtained between 2002 and 2018, with particularly intensive coverage of the 2016–2018 low states. After computing an updated spin ephemeris for the white dwarf (WD), we show that its spin period began to increase in 2014 after having spent 27 yr decreasing; no other intermediate polar has experienced a sign change of its period derivative, but FO Aqr has now done so twice. Our central finding is that the recent low states all occurred shortly after the WD began to spin down, even though no low states were reported in the preceding quarter-century, when it was spinning up. Additionally, the system’s mode of accretion is extremely sensitive to the mass-transfer rate, with accretion being almost exclusively disk-fed when FO Aqr is brighter than V ∼ 14 and substantially stream-fed when it is not. Even in the low states, a grazing eclipse remains detectable, confirming the presence of a disklike structure (but not necessarily a Keplerian accretion disk). We relate these various observations to theoretical predictions that during the low state, the system’s accretion disk dissipates into a non-Keplerian ring of diamagnetic blobs. Finally, a new XMM-Newton observation from a high state in 2017 reveals an anomalously soft X-ray spectrum and diminished X-ray luminosity compared to pre-2016 observations.