Accretion Torus Research Articles

A Compact Binary Merger Model for the Short, Hard GRB 050509b

The first X-ray afterglow for a short ({approx}30ms), hard {gamma}-ray burst was detected by Swift on 9 May 2005 (GRB 050509b). No optical or radio counterpart was identified in follow-up observations. The tentative association of the GRB with a nearby giant elliptical galaxy at redshift z = 0.2248 would imply the progenitor had traveled several tens of kpc from its point of origin, in agreement with expectations linking these events to the final merger of compact binaries driven by gravitational wave emission. We model the dynamical merger of such a system and the time-dependent evolution of the accretion tori thus created. The resulting energetics, variability, and expected durations are consistent with GRB 050509b originating from the tidal disruption of a neutron star by a stellar mass black hole, or of the merger of two neutron stars followed by prompt gravitational collapse of the massive remnant. We discuss how the available {gamma}-ray and X-ray data provides a probe for the nature of the relativistic ejecta and the surrounding medium.

Gravitational waves from oscillating accretion tori: Comparison between different approaches

Quasiperiodic oscillations of high density thick accretion disks orbiting a Schwarzschild black hole have been recently addressed as interesting sources of gravitational waves. The aim of this paper is to compare the gravitational waveforms emitted from these sources when computed using (variations of) the standard quadrupole formula and gauge-invariant metric perturbation theory. To this goal we evolve representative disk models using an existing general relativistic hydrodynamics code which has been previously employed in investigations of such astrophysical systems. Two are the main results of this work: First, for stable and marginally stable disks, no excitation of the black hole quasinormal modes is found. Second, we provide a simple, relativistic modification of the Newtonian quadrupole formula which, in certain regimes, yields excellent agreement with the perturbative approach. This holds true as long as back-scattering of GWs is negligible. Otherwise, any functional form of the quadrupole formula yields systematic errors $\ensuremath{\sim}10%$.

Relativistic outflows from remnants of compact object mergers and their viability for short gamma-ray bursts

We present the first general relativistic hydrodynamic models of the launch and evolution of relativistic jets and winds, driven by thermal energy deposition, possibly due to neutrino-antineutrino annihilation, in the close vicinity of black hole-accretion torus systems. The latter are considered to be the remnants of compact object mergers. Our two-dimensional simulations establish the link between models of such mergers and future observations of short gamma-ray bursts by the SWIFT satellite. They show that ultrarelativistic outflow with maximum terminal Lorentz factors around 1000 develops for polar energy deposition rates above some erg s-1 per steradian, provided the merger environment has a sufficiently low baryon density. By the interaction with the dense accretion torus the ultrarelativistic outflow with Lorentz factors Γ above 100 is collimated into a sharp-edged cone that is embedded laterally by a wind with steeply declining Lorentz factor. The typical semi-opening angles of the cone are , corresponding to about of the hemisphere and apparent isotropized energies (kinetic plus internal) up to ≈erg although at most of the deposited energy is transferred to the outflow with . The viability of post-merger black hole-torus systems as engines of short, hard gamma-ray bursts is therefore confirmed. The annihilation of neutrino-antineutrino pairs radiated from the hot accretion torus appears as a suitable energy source for powerful axial outflow even if only ≈erg are deposited within a cone of half-opening angle around the system axis. Although the torus lifetimes are expected to be only between some 0.01 s and several 0.1 s, our models can explain the durations of all observed short gamma-ray bursts, because different propagation velocities of the front and rear ends will lead to a radial stretching of the ultrarelativistic fireball before transparency is reached. The ultrarelativistic flow reveals a highly non-uniform structure caused by the action of Kelvin-Helmholtz instabilities that originate at the fireball-torus interface. Large radial variations of the baryon density (up to several orders of magnitude) are uncorrelated with moderate variations of the Lorentz factor (factors of a few) and fluctuations of the gently declining radiation-dominated pressure. In the angular direction the Lorentz factor reveals a nearly flat plateau-like maximum with values of several hundreds, that can be located up to off the symmetry axis, and a steep decrease to less than 10 for polar angles larger than . Lateral expansion of the ultrarelativistic core of the flow is prevented by a subsonic velocity component of about towards the symmetry axis, whereas the moderately relativistic wings show a subsonic sideways inflation with less than (measured in the frame comoving with the radial flow).

Oscillations of thick accretion discs around black holes

Abstract We present a numerical study of the response of a thick accretion disc to a localized, external perturbation with the aim of exciting internal modes of oscillation. We find that the perturbations efficiently excite global modes recently identified as acoustic p-modes, and closely related to the epicyclic oscillations of test particles. The two strongest modes occur at eigenfrequencies which are in a 3:2 ratio. We have assumed a constant specific angular momentum distribution within the disc. Our models are in principle scale-free and can be used to simulate accretion tori around stellar or super-massive black holes.

Cosmic Gamma-Ray Bursts: The Big Picture

SummaryA “typical” GRB occurs in a star-forming region of a galaxy at a redshift z~1. In currently popular models, it is caused by the collapse of a massive star which has exhausted its nuclear fuel supply. The star collapses to a black hole threaded by a strong magnetic field, and possibly fed by an accretion torus. Through a variety of processes, electrons are accelerated and gamma-rays, X-rays, optical light, and radio emission ensue, with durations from seconds to years. In this talk, I will review the general observational properties of bursts, their afterglows and host galaxies, and some of the open questions about them.

Dynamics of oscillating relativistic tori around Kerr black holes

We present a comprehensive numerical study of the dynamics of relativistic axisymmetric accretion tori with a power-law distribution of specific angular momentum orbiting in the background space‐time of a Kerr black hole. By combining general relativistic hydrodynamics simulations with a linear perturbative approach we investigate the main dynamical properties of these objects over a large parameter space. The astrophysical implications of our results extend and improve two interesting results that have been recently reported in the literature. First, the induced quasi-periodic variation of the mass quadrupole moment makes relativistic tori of nuclear matter densities, as those formed during the last stages of binary neutron star mergers, promising sources of gravitational radiation, potentially detectable by interferometric instruments. Secondly, p-mode oscillations in relativistic tori of low rest-mass densities could be used to explain high-frequency quasi-periodic oscillations observed in X-ray binaries containing a black hole candidate under conditions more generic than those considered so far. Ke yw ords: accretion, accretion discs ‐ gravitational waves ‐ hydrodynamics ‐ relativity.

Disks, tori, and cocoons: emission and absorption diagnostics of AGN environments

One of the most important problems in the study of active galaxies is understanding the detailed geometry, physics, and evolution of the central engines and their environments. The leading models involve an accretion disk and torus structure around a central dense object, thought to be a supermassive black hole. Gas found in the environment of active galactic nuclei (AGN) is associated with different structures: molecular accretion disks, larger scale atomic tori, ionized and neutral “cocoons” in which the nuclear regions can be embedded. All of them can be studied at radio wavelengths by various means. Here, we summarize the work that has been done to date in the radio band to characterize these structures. Much has been learned about the central few parsecs of AGN in the last few decades with contemporary instruments but the picture remains incomplete. In order to be able to define a more accurate model of this region, significant advances in sensitivity, spectral and angular resolution, and bandpass stability are required. The necessary advances will only be provided by the Square Kilometer Array and we discuss the possibilities that these dramatic improvements will open for the study of the gas in the central region of AGN.

Accretion Signatures from Massive Young Stellar Objects

High resolution (lambda / Delta-lambda = 50,000) K-band spectra of massive, embedded, young stellar objects are presented. The present sample consists of four massive young stars located in nascent clusters powering Galactic giant H II regions. Emission in the 2.3 micron 2--0 vibrational--rotational bandhead of CO is observed. A range of velocity broadened profiles seen in three of the objects is consistent with the emission arising from a circumstellar disk seen at various inclination angles. Br gamma spectra of the same spectral and spatial resolution are also presented which support an accretion disk or torus model for massive stars. In the fourth object, Br emission suggesting a rotating torus is observed, but the CO profile is narrow, indicating that there may be different CO emission mechanisms in massive stars and this is consistent with earlier observations of the BN object and MWC 349. To--date, only young massive stars of late O or early B types have been identified with clear accretion disk signatures in such embedded clusters. Often such stars are found in the presence of other more massive stars which are revealed by their photospheric spectra but which exhibit no disk signatures. This suggests the timescale for dissipating their disks is much faster than the less massive OB stars or that the most massive stars do not form with accretion disks.

Radiation transfer of emission lines in curved space-time

We present a general formulation for ray-tracing calculations in curved space-time. The formulation takes full account of relativistic effects in the photon transport and the relative motions of the emitters and the light-of-sight absorbing material. We apply the formulation to calculate the emission from accretion disks and tori around rotating black holes.

Resonance in Forced Oscillations of an Accretion Disk and Kilohertz Quasi-periodic Oscillations

We have performed numerical simulations of a radially perturbed accretion torus around a black hole or neutron star and find that the torus performs radial and vertical motions at the appropriate epicyclic frequencies. We find clear evidence that vertical motions are excited in a nonlinear resonance when the applied perturbation is periodic in time. The strongest resonant response occurs when the frequency difference of the two oscillations is equal to one-half the forcing frequency, precisely as recently observed in the accreting pulsar SAX J1808.4-3658, where the observed kilohertz quasi-periodic oscillation peak separation is half the spin frequency of 401 Hz.

Global General Relativistic Magnetohydrodynamic Simulations of Accretion Tori

This paper presents an initial survey of the properties of accretion flows in the Kerr metric from three-dimensional, general relativistic magnetohydrodynamic simulations of accretion tori. We consider three fiducial models of tori around rotating, both prograde and retrograde, and nonrotating black holes; these three fiducial models are also contrasted with axisymmetric simulations and a pseudo-Newtonian simulation with equivalent initial conditions to delineate the limitations of these approximations.

A new simple model for high-frequency quasi-periodic oscillations in black hole candidates

Observations of X-ray emissions from binary systems have long since been considered important tools to test general relativity in strong-field regimes. The high-frequency quasi-periodic oscillations (HFQPOs) observed in binaries containing a black hole candidate, in particular, have been proposed as a means to measure more directly the properties of the black hole, such as its mass and spin. Numerous models have been suggested to explain the HFQPOs and the rich phenomenology accompanying them. Many of these models rest on a number of assumptions and are at times in conflict with the most recent observations. We here propose a new, simple model in which the HFQPOs result from basic p-mode oscillations of a small accretion torus orbiting close to the black hole. We show that within this model the key properties of the HFQPOs can be explained simply, given a single reasonable assumption. We also discuss observational tests that can refute the model.

Accretion of Low Angular Momentum Material onto Black Holes: Two‐dimensional Magnetohydrodynamic Case

We report on the second phase of our study of slightly rotating accretion flows onto black holes. We consider magnetohydrodynamical (MHD) accretion flows with a spherically symmetric density distribution at the outer boundary but with spherical symmetry broken by the introduction of a small, latitude-dependent angular momentum and a weak radial magnetic field. We study accretion flows by means of numerical two-dimensional, axisymmetric, MHD simulations with and without resistive heating. Our main result is that the properties of the accretion flow depend mostly on an equatorial accretion torus that is made of the material that has too much angular momentum to be accreted directly. The torus accretes, however, because of the transport of angular momentum due to the magnetorotational instability (MRI). Initially, accretion is dominated by the polar funnel, as in the hydrodynamic inviscid case, where material has zero or very low angular momentum. At the later phase of the evolution, the torus thickens toward the poles and develops a corona or an outflow or both. Consequently, the mass accretion through the funnel is stopped. The accretion of rotating gas through the torus is significantly reduced compared with the accretion of nonrotating gas (i.e., the Bondi rate). It is also much smaller than the accretion rate in the inviscid, weakly rotating case. Our results do not change if we switch on or off resistive heating. Overall our simulations are very similar to those presented by Stone, Pringle, Hawley, and Balbus despite different initial and outer boundary conditions. Thus, we confirm that MRI is very robust and controls the nature of radiatively inefficient accretion flows. Although the time-averaged properties of our models approach a steady state, we find that the instantaneous mass-accretion rate in the latter stages of our simulations is highly time-dependent, with the inner flow displaying three generic flow patterns.

What Happened to the NGC 6251 Counterjet?

We have used the Very Long Baseline Array to produce a high dynamic range image of the nucleus of NGC 6251 at 1.6 GHz and snapshot images at 5.0, 8.4, and 15.3 GHz to search for emission from a parsec-scale counterjet. Previous VLBI images at 1.6 GHz have set a lower limit for the jet/counterjet brightness ratio near the core at about 80:1, which is larger than expected, given the evidence that the radio axis is fairly close to the plane of the sky. A possible explanation is that the inner few parsecs of the counterjet are hidden by free-free absorption by ionized gas associated with an accretion disk or torus. This would be consistent with the nearly edge-on appearance of the arcsecond-scale dust disk seen in the center of NGC 6251 by Hubble Space Telescope (HST). We detect counterjet emission close to the core at 1.6 GHz, but not at the higher frequencies. Given that the optical depth of free-free absorption falls off more rapidly with increasing frequency than the optically thin synchrotron emission from a typical radio jet, this result implies that the absence of a detectable parsec-scale counterjet at high frequencies is not due to free-free absorption unless the density of ionized gas is extremely high and we have misidentified the core at 1.6 GHz. The most likely alternative is a large jet/counterjet brightness ratio caused by relativistic beaming, which in turn requires the inner radio axis to be closer to our line of sight than the orientation of the HST dust disk would suggest.

Can Neutrino‐cooled Accretion Disks Be an Origin of Gamma‐Ray Bursts?

It is often proposed that a massive torus with approximately solar mass surrounding a stellar-mass black hole could be a central engine of gamma-ray bursts. We study the properties of such massive accretion tori (or disks) based on the α viscosity model. For surface density exceeding about 1020 g cm-2, which is realized when ~1 M☉ of material is contained within a disk of size ~5 × 106 cm, we find that (1) the luminosity of photons is practically zero because of significant photon trapping, (2) neutrino cooling dominates over advective cooling, (3) the pressure of degenerate electrons dominates over the pressure of gas and photons, and (4) the magnetic field strength exceeds the critical value of about 4 × 1013 G, even if we take 0.1% of the equipartition value. The possible observable quantum electrodynamical (QED) effects arising from supercritical fields are discussed. Most interestingly, photon splitting may occur, producing a significant number of photons of energies below ~511 keV, thereby possibly suppressing e± pair creation.

Three‐dimensional Hydrodynamic Simulations of Accretion Tori in Kerr Spacetimes

This paper presents results of three-dimensional simulations of global hydrodynamic instabilities in black hole tori, extending earlier work by Hawley to Kerr spacetimes. This study probes a three-dimensional parameter space of torus angular momentum, torus size, and black hole angular momentum. We have observed the growth of the Papaloizou-Pringle instability for a range of torus configurations and the resultant formation of m=1 planets. We have also observed the quenching of this instability in the presence of early accretion flows; however, in one simulation both early accretion and planet formation occurred. Though most of the conclusions reached in Hawley's earlier work on Schwarzschild black holes carry over to Kerr spacetime, the presence of frame dragging in the Kerr geometry adds an element of complexity to the simulations; we have seen especially clear examples of this phenomenon in the accretion flows that arise from retrograde tori.

Models of rapidly rotating neutron stars: remnants of accretion-induced collapse

Equilibrium models of differentially rotating nascent neutron stars are constructed, which represent the result of the accretion induced collapse of rapidly rotating white dwarfs. The models are built in a two-step procedure: (1) a rapidly rotating pre-collapse white dwarf model is constructed; (2) a stationary axisymmetric neutron star having the same total mass and angular momentum distribution as the white dwarf is constructed. The resulting collapsed objects consist of a high density central core of size roughly 20 km, surrounded by a massive accretion torus extending over 1000 km from the rotation axis. The ratio of the rotational kinetic energy to the gravitational potential energy of these neutron stars ranges from 0.13 to 0.26, suggesting that some of these objects may have a non-axisymmetric dynamical instability that could emit a significant amount of gravitational radiation.

Newtonian hydrodynamics of the coalescence of black holes with neutron stars - IV. Irrotational binaries with a soft equation of state

We present the results of three-dimensional hydrodynamical simulations of the final stages of in-spiral in a black hole–neutron star binary, when the separation is comparable to the stellar radius. We use a Newtonian smooth particle hydrodynamics (SPH) code to model the evolution of the system, and take the neutron star to be a polytrope with a soft (adiabatic indices and equation of state and the black hole to be a Newtonian point mass. The only non-Newtonian effect we include is a gravitational radiation back reaction force, computed in the quadrupole approximation for point masses. We use irrotational binaries as initial conditions for our dynamical simulations, which are begun when the system is on the verge of initiating mass transfer and followed for approximately 23 ms. For all the cases studied we find that the star is disrupted on a dynamical time-scale, and forms a massive accretion torus around the spinning (Kerr) black hole. The rotation axis is clear of baryons (less than 10−5 M⊙ within 10°) to an extent that would not preclude the formation of a relativistic fireball capable of powering a cosmological gamma-ray burst. Some mass (the specific amount is sensitive to the stiffness of the equation of state) may be dynamically ejected from the system during the coalescence and could undergo r-process nucleosynthesis. We calculate the waveforms, luminosities and energy spectra of the gravitational radiation signal, and show how they reflect the global outcome of the coalescence process.

Black Hole-Neutron Star Mergers as Central Engines of Gamma-Ray Bursts.

Hydrodynamic simulations of the merger of stellar mass black hole-neutron star binaries are compared with mergers of binary neutron stars. The simulations are Newtonian but take into account the emission and back-reaction of gravitational waves. The use of a physical nuclear equation of state allows us to include the effects of neutrino emission. For low neutron star-to-black hole mass ratios, the neutron star transfers mass to the black hole during a few cycles of orbital decay and subsequent widening before finally being disrupted, whereas for ratios near unity the neutron star is destroyed during its first approach. A gas mass between approximately 0.3 and approximately 0.7 M middle dot in circle is left in an accretion torus around the black hole and radiates neutrinos at a luminosity of several times 1053 ergs s-1 during an estimated accretion timescale of about 0.1 s. The emitted neutrinos and antineutrinos annihilate into e+/- pairs with efficiencies of 1%-3% and rates of up to approximately 2x1052 ergs s-1, thus depositing an energy Enunu&d1; less, similar1051 ergs above the poles of the black hole in a region that contains less than 10-5 M middle dot in circle of baryonic matter. This could allow for relativistic expansion with Lorentz factors around 100 and is sufficient to explain apparent burst luminosities Lgamma approximately Enunu&d1;&solm0;&parl0;fOmegatgamma&parr0; up to several times 1053 ergs s-1 for burst durations tgamma approximately 0.1-1 s, if the gamma emission is collimated in two moderately focused jets in a fraction fOmega=2deltaOmega&solm0;&parl0;4pi&parr0; approximately 1&solm0;100-(1/10) of the sky.

Global Magnetohydrodynamical Simulations of Accretion Tori

Global time-dependent simulations provide a means to investigate time-dependent dynamic evolution in accretion disks. This paper seeks to extend previous local simulations by beginning a systematic effort to develop fully global three-dimensional simulations. The nonlinear development of the magnetorotational instability is investigated using a time-explicit finite difference code written in coordinates. The equations of ideal magnetohydrodynamics are solved with the assumption of an adiabatic equation of state. Both a Newtonian potential and a pseudo-Newtonian potential are used. Two simplifications are also explored: a gravitational potential (the cylindrical disk) and axisymmetry. The results from those simulations are compared with fully three-dimensional global simulations. The global simulations begin with equilibrium pressure-supported accretion tori. Two different initial field geometries are investigated: poloidal fields that are constant along initial equidensity surfaces and toroidal fields with a constant ratio of gas to magnetic pressure. In both cases the magnetorotational instability rapidly develops and the torus becomes turbulent. The resulting turbulence transports angular momentum, and the torus develops an angular momentum distribution that is near Keplerian. A comparison with axisymmetric simulations shows that in three dimensions the magnetorotational instability can act as a dynamo and regenerate poloidal field, thereby sustaining the turbulence. As previously observed in local simulations, the stress is dominated by the Maxwell component. The total stress in the interior of the disk is ≈0.1-0.2 times the thermal pressure. At late time the disks are characterized by relatively thick configurations, with rapid time dependence and tightly wrapped, low-m spiral structures.