Unstable Accretion Disks Research Articles

Constraining the Evolution of the Unstable Accretion Disk in SMC X-1 with NICER

Neutron star high-mass X-ray binaries with superorbital modulations in luminosity host warped inner accretion disks that occult the neutron star during precession. In SMC X-1, the instability in the warped disk geometry causes superorbital period “excursions”: times of instability when the superorbital period decreases from its typical value of 55 to ∼40 days. Disk instability makes SMC X-1 an ideal system in which to investigate the effects of variable disk geometry on the inner accretion flow. Using the high-resolution spectral and timing capabilities of the Neutron Star Interior Composition Explorer, we examined the high state of four different superorbital cycles of SMC X-1 to search for changes in spectral shape and connections to the unstable disk geometry. We performed pulse phase-averaged and phase-resolved spectroscopy to closely compare the changes in spectral shape and any cycle-to-cycle variations. While some parameters, including the photon index and absorbing column density, show slight variations with superorbital phase, these changes are most evident during the intermediate state of the superorbital cycle. Few spectral changes are observed within the high state of the superorbital cycle, possibly indicating the disk instability does not significantly change SMC X-1's accretion process.

An Outburst by AM CVn Binary SDSS J113732.32+405458.3

Abstract We report the discovery of a one magnitude increase in the optical brightness of the 59.63 minutes orbital period AM CVn binary SDSS J113732.32+405458.3. Public g, r, and i band data from the Zwicky Transient Facility exhibit a subsequent decline over a 300 days period, while a few data points from commissioning show that the peak was likely seen. Such an outburst is likely due to a change in the state of the accretion disk, making this the longest period AM CVn binary to reveal an unstable accretion disk. The object is now back to its previously observed (by SDSS and PS-1) quiescent brightness that is likely set by the accreting white dwarf. Prior observations of this object also imply that the recurrence times for such outbursts are likely more than 12 yr.

Gravitational wave emission from unstable accretion discs in tidal disruption events

Abstract Gravitational waves can be emitted by accretion discs if they undergo instabilities that generate a time varying mass quadrupole. In this work we investigate the gravitational signal generated by a thick accretion disc of 1 M⊙ around a static supermassive black hole of 106 M⊙, assumed to be formed after the tidal disruption of a solar type star. This torus has been shown to be unstable to a global non-axisymmetric hydrodynamic instability, the Papaloizou–Pringle instability, in the case where it is not already accreting and has a weak magnetic field. We start by deriving analytical estimates of the maximum amplitude of the gravitational wave signal, with the aim to establish its detectability by the Laser Interferometer Space Antenna (LISA). Then, we compare these estimates with those obtained through a numerical simulation of the torus, made with a 3D smoothed particle hydrodynamics code. Our numerical analysis shows that the measured strain is two orders of magnitude lower than the maximum value obtained analytically. However, accretion discs affected by the Papaloizou–Pringle instability may still be interesting sources for LISA, if we consider discs generated after deeply penetrating tidal disruptions of main-sequence stars of higher mass.

The First Bird’s-eye View of a Gravitationally Unstable Accretion Disk in High-mass Star Formation

Abstract We report on the first bird’s-eye view of the innermost accretion disk around the high-mass protostellar object G353.273+0.641, taken by Atacama Large Millimeter/submillimeter Array long baselines. The disk traced by dust continuum emission has a radius of 250 au, surrounded by the infalling rotating envelope traced by thermal CH3OH lines. This disk radius is consistent with the centrifugal radius estimated from the specific angular momentum in the envelope. The lower-limit envelope mass is ∼5–7 M ☉ and accretion rate onto the stellar surface is 3 × 10−3 M ☉ yr−1 or higher. The expected stellar age is well younger than 104 yr, indicating that the host object is one of the youngest high-mass objects at present. The disk mass is 2–7 M ☉, depending on the dust opacity index. The estimated Toomre’s Q parameter is typically 1–2 and can reach 0.4 at the minimum. These Q values clearly satisfy the classical criteria for gravitational instability, and are consistent with recent numerical studies. Observed asymmetric and clumpy structures could trace a spiral arm and/or disk fragmentation. We found that 70% of the angular momentum in the accretion flow could be removed via the gravitational torque in the disk. Our study has indicated that the dynamical nature of a self-gravitating disk could dominate the early phase of high-mass star formation. This is remarkably consistent with the early evolutionary scenario of a low-mass protostar.

Gravitational stability and fragmentation condition for discs around accreting supermassive stars

Supermassive stars (SMSs) with mass $\sim10^{5}~\rm{M}_{\odot}$ are promising candidates for the origin of supermassive black holes observed at redshift $\gtrsim6$. They are supposed to form as a result of rapid accretion of primordial gas, although it can be obstructed by the time variation caused by circum-stellar disc fragmentation due to gravitational instability. To assess the occurrence of fragmentation, we study the structure of marginally gravitationally unstable accretion discs, by using a steady one-dimensional thin disc model with detailed treatment of chemical and thermal processes. Motivated by two SMS formation scenarios, i.e., those with strong ultraviolet radiation background or with large velocity difference between the baryon and the dark matter, we consider two types of flows, i.e., atomic and molecular flows, respectively, for a wide range of the central stellar mass $10-10^5~\rm{M}_{\odot}$ and the accretion rate $10^{-3}-1~\rm{M}_{\odot}~\rm{yr}^{-1}$. In the case of a mostly atomic gas flowing to the disc outer boundary, the fragmentation condition is expressed as the accretion rate being higher than the critical value of $10^{-1}~\rm{M}_{\odot}~\rm{yr}^{-1}$ regardless of the central stellar mass. On the other hand, in the case of molecular flows, there is a critical disc radius outside of which the disc becomes unstable. Those conditions appears to be marginally satisfied according to numerical simulations, suggesting that disc fragmentation can be common during SMS formation.

VISCOSITY PRESCRIPTION FOR GRAVITATIONALLY UNSTABLE ACCRETION DISKS

Gravitationally unstable accretion disks emerge in a variety of astrophysical contexts - giant planet formation, FU Orioni outbursts, feeding of AGNs, and the origin of Pop III stars. When a gravitationally unstable disk is unable to cool rapidly it settles into a quasi-stationary, fluctuating gravitoturbulent state, in which its Toomre Q remains close to a constant value Q_0~1. Here we develop an analytical formalism describing the evolution of such a disk, which is based on the assumptions of Q=Q_0 and local thermal equilibrium. Our approach works in the presence of additional sources of angular momentum transport (e.g. MRI), as well as external irradiation. Thermal balance dictates a unique value of the gravitoturbulent stress \alpha_{gt} driving disk evolution, which is a function of the local surface density and angular frequency. We compare this approach with other commonly used gravitoturbulent viscosity prescriptions, which specify the explicit dependence of stress \alpha_{gt} on Toomre Q in an ad hoc fashion, and identify the ones that provide consistent results. We nevertheless argue that our Q=Q_0 approach is more flexible, robust, and straightforward, and should be given preference in applications. We illustrate this with a couple of analytical calculations - locations of the snow line and of the outer edge of the dead zone in a gravitoturbulent protoplanetary disk - which clearly show the simplicity and versatility of the Q=Q_0 approach.

P-mode oscillation on slim discs

AbstractWe numerically investigate the thermally unstable accretion discs around spinning black holes with different spins. We adopted an additional evolutionary viscosity equation to replace the standard alpha-prescription based on the results of two MHD simulations. We find an interesting oscillation when accretion switches to slim disc mode. The oscillation arises from the sonic point of accretion flow and propagates outwards. We mimic the bolometric light-curve and find a series of harmonics on its power spectrum. The frequency ratio of those harmonics is a regular integer series. The lowest frequency of the harmonics is identical to the prediction of trapped p-mode in QPO theory.

ORBITAL DECAY AND EVIDENCE OF DISK FORMATION IN THE X-RAY BINARY PULSAR OAO 1657–415

OAO 1657-415 is an eclipsing X-ray binary wind-fed pulsar that has exhibited smooth spin-up/spin-down episodes and has undergone several torque reversals throughout its long history of observation. We present a frequency history spanning nearly 19 years of observations from the Burst and Transient Source Experiment (CGRO/BATSE) and from the Gamma-Ray Burst Monitor (Fermi/GBM). The analysis suggests two modes of accretion: one resulting in steady spin-up during which we believe a stable accretion disk is present and one that results in what appears to be a random walk in spin frequency where an unstable accretion disk forms alternating in direction ("flip flop"). Orbital elements of the pulsar system are determined at several intervals throughout this history. With these ephemerides, statistically significant orbital decay ($\dot{P}/P =(-3.40 \pm0.15)\times10^{-6} \mathrm{yr^{-1}}$) is established suggesting a transition between wind-fed and disk-mediated accretion.

PERIODIC STRUCTURE IN THE MEGAPARSEC-SCALE JET OF PKS 0637–752

We present 18 GHz Australia Telescope Compact Array imaging of the Mpc-scale quasar jet PKS 0637-752 with angular resolution ~0.58 arcseconds. We draw attention to a spectacular train of quasi-periodic knots along the inner 11 arcseconds of the jet, with average separation ~1.1 arcsec (7.6 kpc projected). We consider two classes of model to explain the periodic knots: those that involve a static pattern through which the jet plasma travels (e.g. stationary shocks); and those that involve modulation of the jet engine. Interpreting the knots as re-confinement shocks implies the jet kinetic power Q ~ 10^{46} erg/s, but the constant knot separation along the jet is not expected in a realistic external density profile. For models involving modulation of the jet engine, we find that the required modulation period is 2 x 10^3 yr < \tau < 3 x 10^5 yr. The lower end of this range is applicable if the jet remains highly relativistic on kpc-scales, as implied by the IC/CMB model of jet X-ray emission. We suggest that the quasi-periodic jet structure in PKS 0637-752 may be analogous to the quasi-periodic jet modulation seen in the microquasar GRS 1915+105, believed to result from limit cycle behaviour in an unstable accretion disk. If variations in the accretion rate are driven by a binary black hole, the predicted orbital radius is 0.7 < a < 30 pc, which corresponds to a maximum angular separation of ~0.1 - 5 mas.

Inside-out outbursts in the thermally unstable accretion disk of the neutron star microquasar Aquila X-1

The observed features of typical, symmetric outbursts in the neutron star X-ray transient Aquila X-1 such as the rise, the plateau in the maximum, and the decay, are reproduced with a timedependent disk instability model. No external effect such as irradiation is necessary to account for the observation. In the model, a truncated inner disk radius that is more than ten times larger than the magnetically truncated one and that is plausibly appropriate to the scale proposed by hot flow models must be adopted. The model shows that such symmetric outbursts in Aql X-1 may be inside-out outbursts, while outside-in outbursts have been conjectured based on observations. The implication of the presented model, with observations, is discussed together.

Effects of Black Hole Spin on the Limit-Cycle Behaviour of Accretion Disks

We present a spatially 1.5-dimensional, time-dependent numerical study of accretion disks around Kerr black holes. Our study focuses on the limit-cycle behavior of thermally unstable accretion disks. We find that maximal luminosity may be a more appropriate probe of black hole spin than the cycle duration and influence radius.

NUMERICAL TREATMENT OF THIN ACCRETION DISK DYNAMICS AROUND ROTATING BLACK HOLES

In the present study, we perform the numerical simulation of a relativistic thin accretion disk around the nonrotating and rapidly rotating black holes using the general relativistic hydrodynamic code with Kerr in Kerr–Schild coordinate that describes the central rotating black hole. Since the high energy X-rays are produced close to the event horizon resulting the black hole–disk interaction, this interaction should be modeled in the relativistic region. We have set up two different initial conditions depending on the values of thermodynamical variables around the black hole. In the first setup, the computational domain is filled with constant parameters without injecting gas from the outer boundary. In the second, the computational domain is filled with the matter which is then injected from the outer boundary. The matter is assumed to be at rest far from the black hole. Both cases are modeled over a wide range of initial parameters such as the black hole angular momentum, adiabatic index, Mach number and asymptotic velocity of the fluid. It has been found that initial values and setups play an important role in determining the types of the shock cone and in designating the events on the accretion disk. The continuing injection from the outer boundary presents a tail shock to the steady state accretion disk. The opening angle of shock cone grows as long as the rotation parameter becomes larger. A more compressible fluid (bigger adiabatic index) also presents a bigger opening angle, a spherical shock around the rotating black hole, and less accumulated gas in the computational domain. While results from [J. A. Font, J. M. A. Ibanez and P. Papadopoulos, Mon. Not. R. Astron. Soc.305 (1999) 920] indicate that the tail shock is warped around for the rotating hole, our study shows that it is the case not only for the warped tail shock but also for the spherical and elliptical shocks around the rotating black hole. The warping around the rotating black hole in our case is much smaller than the one by [J. A. Font, J. M. A. Ibanez and P. Papadopoulos, Mon. Not. R. Astron. Soc.305 (1999) 920], due to the representation of results at the different coordinates. Contrary to the nonrotating black hole, the tail shock is slightly warped around the rotating black hole. The filled computational domain without any injection leads to an unstable accretion disk. However much of it reaches a steady state for a short period of time and presents quasi-periodic oscillation (QPO). Furthermore, the disk tends to loose mass during the whole dynamical evolution. The time-variability of these types of accretion flowing close to the black hole may clarify the light curves in Sgr A*.

RS Ophiuchi: thermonuclear explosion or disc instability?

Abstract Sokoloski, Rupen &amp; Mioduszewski have recently reported evidence that the recurrent nova RS Ophiuchi produced a pair of highly collimated radio jets within days of its 2006 outburst. This suggests that an accretion disc must be present during the outburst. However, in the standard picture of recurrent novae as thermonuclear events, any such disc must be expelled from the white dwarf vicinity, as the nuclear energy yield greatly exceeds its binding energy. We suggest instead that the outbursts of RS Oph are thermal-viscous instabilities in a disc irradiated by the central accreting white dwarf. The distinctive feature of RS Oph is the very large size of its accretion disc. Given this, it fits naturally into a consistent picture of systems with unstable accretion discs. This picture explains the presence and speed of the jets, the brightness and duration of the outburst, and its rise time and linear decay, as well as the faintness of the quiescence. By contrast, the hitherto standard picture of recurrent thermonuclear explosions has a number of severe difficulties. These include the presence of jets, the faintness of quiescence, and the fact the accretion disc must be unstable unless it is far smaller than any reasonable estimate.

Effects of the energy equation in studies of limit-cycle behaviors of black hole accretion disks

The inconsistency of the energy equation used in the literature is pointed out and a new consistent energy equation is given. With this new energy equation, calculations are made for the limit-cycle behaviors of thermally unstable accretion disks around black holes. From the comparison of our numerical results with those obtained using the inconsistent energy equation, it is found that the inconsistent energy equation undervalues the temperature and overvalues the effective optical depth when the accreted gas becomes effectively optically thin. Thus, it is dangerous if the inconsistent energy equation is used in the studies of very hot and optically thin accretion flows such as advection-dominated accretion flows (ADAFs), and our new energy equation is likely to be a better alternative.

Self‐Gravitating Fragmentation of Eccentric Accretion Disks

We consider the effects of eccentricity on the fragmentation of gravitationally unstable accretion disks, using numerical hydrodynamics. We find that eccentricity does not affect the overall stability of the disk against fragmentation, but significantly alters the manner in which such fragments accrete gas. Variable tidal forces around an eccentric orbit slow the accretion process, and suppress the formation of weakly-bound clumps. The "stellar" mass function resulting from the fragmentation of an eccentric disk is found to have a significantly higher characteristic mass than that from a corresponding circular disk. We discuss our results in terms of the disk(s) of massive stars at ~0.1pc from the Galactic Center, and find that the fragmentation of an eccentric accretion disk, due to gravitational instability, is a viable mechanism for the formation of these systems.

Studies of Thermally Unstable Accretion Disks around Black Holes with Adaptive Pseudospectral Domain Decomposition. I. Limit‐Cycle Behavior in the Case of Moderate Viscosity

We present a numerical method for spatially 1.5-dimensional, time-dependent studies of accretion disks around black holes. The method originates from a combination of the standard pseudospectral and adaptive domain decomposition methods found in the literature, but with a number of improvements in both the numerical and physical senses. In particular, we introduce a new treatment for the connection at the interfaces of decomposed subdomains, construct an adaptive function for the mapping between the Chebyshev-Gauss-Lobatto collocation points and the physical collocation points in each subdomain, and modify the oversimplified one-dimensional basic equations of accretion flows to account for the effects of viscous stresses in both the azimuthal and radial directions. Our method is verified by reproducing the best results obtained previously by Szuszkiewicz & Miller on the limit-cycle behavior of thermally unstable accretion disks with moderate viscosity. A new finding is that, according to our computations, the Bernoulli function of the matter in such disks is always and everywhere negative, so that outflows are unlikely to originate from these disks. We are encouraged to study the more difficult case of thermally unstable accretion disks with strong viscosity in ongoing work.

On Be star candidates and possible blue pre-main sequence objects in the Small Magellanic Cloud

Recently the OGLE experiment has provided accurate light curves and colours for about 2 millions stars in the Small Magellanic Cloud. We have examined this database for its content of Be stars, applying some selection criteria, and we have found a sample of ~1000 candidates. Some of these stars show beautiful light curves with amazing variations never observed in any Galactic variable. We find outbursts in 13% of the sample (type-1 stars), high and low states in 15%, periodic variations in 7%, and the usual variations seen in Galactic Be stars in 65% of the cases. The Galactic counterparts of type-1 objects could be the outbursting Be stars found by Hubert & Floquet ([CITE]) after the analysis of Hipparcos photometry. We discuss the possibility that type-1 stars could correspond to Be stars with accreting white dwarf companions or alternatively, blue pre-main sequence stars surrounded by thermally unstable accretion disks. We provide coordinates and basic photometric information for these stars and some examples of light curves.

The Structure and Evolution of Circumbinary Disks in Cataclysmic Variable Systems

We investigate the structure and evolution of a geometrically thin, viscous, Keplerian circumbinary (CB) disk, using detailed models of their radiative/convective vertical structure. We use a simplified description for the evolution of the cataclysmic binary and focus on cases where the circumbinary disk causes accelerated mass transfer (10-8 M☉ yr-1). The inner edge of the disk is assumed to be determined by the tidal truncation radius, and the mass input rate into the disk is assumed to be a small fraction (10-5 to 0.01) of the mass transfer rate. Under the action of the viscous stresses in the disk the matter drifts outward with the optically thick region extending to several AU. The inner part of the disk is cool with maximum effective temperatures 3000 K, while the outermost parts of the disk are 30 K and optically thin. We calculate the effects of thermal instability on a sufficiently massive CB disk. It leads to outbursts reminiscent of those in thermally unstable accretion disks, with the instability remaining confined to the inner regions of the CB disk. However, for most of the evolutionary sequences, the surface densities required to trigger instability are not reached. The spectral energy distributions from circumbinary disks are calculated, and the prospects for the detection of such disks in the infrared and submillimeter wavelength regions are discussed.

An RXTE observation of the intermediate polar XY arietis

RXTE's observation of XY Ari covered more eclipses of this close binary at a higher count rate than ever before. The eclipses located the accretion regions on the white dwarf and showed that they covered < 0.002 of the white dwarf surface. Additionally we recorded the first outburst of XY Ari seen, allowing us to watch as an unstable accretion disk overwhelmed the magnetic field of the white dwarf and pushed inwards, cutting off our line-of-sight to the lower accretion pole. We also find limits on the mass of the white dwarf.

Deriving the Quasar Luminosity Function from Accretion-Disk Instabilities

We have derived the quasar luminosity function assuming that the quasar activity is driven by a thermal-viscous unstable accretion disk around a supermassive black hole. The instabilities produce large amplitude, long-term variability of a single source. We take a light curve of a single source and calculate the luminosity function, from the function of time it spends at each luminosity. Convolving this with an assumed mass distribution we fit well the observed optical luminosity function of quasars at four redshifts. As a result we obtain the evolution of the mass distribution between redshifts 2.5 and 0.5. The maximum of the active black hole mass function moves towards lower mass by a factor ~10 at the low redshift. The number of high mass sources declines rapidly, and so low mass sources become dominant at lower redshift. The main conclusions are following: 1) The quasar long-term variability due to the disk thermal-viscous instabilities provides a natural explanation for the observed quasar luminosity function. 2) The peak of the mass function evolves towards lower black hole masses at lower redshifts by a factor ~10. 3) High mass sources die subsequently when redshift gets smaller. 4) The number of high mass sources declines rapidly, and so low mass sources become dominant at lower redshift. 5) The periodic outbursts of activity appear as long as the matter is supplied to the accretion disk. 6) Since the time-averaged accretion rate is low, the remnant sources (or sources in the low activity phase) do not grow to very massive black holes. 7) A continuous fuel supply at a relatively low accretion rate (~0.01 - 0.1 \dot M_{Edd}$) for each single source is required over the lifetime of the entire quasar population.