Search for: All records

Abstract Data derived from general relativistic magnetohydrodynamic simulations of accretion onto black holes can be used as input to a postprocessing scheme that predicts the radiated spectrum. Combining a relativistic Compton scattering radiation transfer solution in the corona with detailed local atmosphere solutions incorporating local ionization and thermal balance within the disk photosphere, it is possible to study both spectral formation and intrinsic spectral variability in the radiation from relativistic accretion disks. With this method, we find that radiatively efficient systems with black holes of 10M_⊙accreting at ≈0.01 in Eddington units produce spectra very similar to those observed in the hard states of X-ray binaries. The spectral shape above 10 keV is well described by a power law with an exponential cutoff. Intrinsic turbulent variations lead to order-unity changes in bolometric luminosity, variations in the logarithmic spectral slope ∼0.1, and factor of 2 alterations in the cutoff energy on timescales ∼50 (M_BH/10M_⊙) ms. Within the corona, the range of gas temperature spans more than 1 order of magnitude. The wide distribution of temperatures is central to defining the spectrum: the logarithmic spectral slope is harder by ∼0.3 and the cutoff energy larger by a factor ∼10–30 than if the coronal temperature everywhere were its mass-weighted mean. 

Abstract Black hole–neutron star binaries are of interest in many ways: they are intrinsically transient, radiate gravitational waves detectable by LIGO, and may produceγ-ray bursts. Although it has long been assumed that their late-stage orbital evolution is driven entirely by gravitational wave emission, we show here that in certain circumstances, mass transfer from the neutron star onto the black hole can both alter the binary's orbital evolution and significantly reduce the neutron star's mass: when the fraction of its mass transferred per orbit is ≳10⁻², the neutron star's mass diminishes by order unity, leading to mergers in which the neutron star mass is exceptionally small. The mass transfer creates a gas disk around the black holebeforemerger that can be comparable in mass to the debris remaining after merger, i.e., ~0.1M_⊙. These processes are most important when the initial neutron star–black hole mass ratioqis in the range ≈0.2–0.8, the orbital semimajor axis is 40 ≲ a₀/r_g ≲ 300 (r_g ≡ GM_BH/c²), and the eccentricity is large ate₀ ≳ 0.8. Systems of this sort may be generated through the dynamical evolution of a triple system, as well as by other means. 

Abstract The detection of GW170817/AT2017gfo inaugurated an era of multimessenger astrophysics, in which gravitational-wave and multiwavelength photon observations complement one another to provide unique insight into astrophysical systems. A broad theoretical consensus exists, in which the photon phenomenology of neutron star mergers largely rests upon the evolution of the small amount of matter left on bound orbits around the black hole or massive neutron star remaining after the merger. Because this accretion disk is far from inflow equilibrium, its subsequent evolution depends very strongly on its initial state, yet very little is known about how this state is determined. Using both snapshot and tracer particle data from a numerical relativity/MHD simulation of an equal-mass neutron star merger that collapses to a black hole, we show how gravitational forces arising in a nonaxisymmetric, dynamical spacetime supplement hydrodynamical effects in shaping the initial structure of the bound debris disk. The work done by hydrodynamical forces is ∼10 times greater than that due to time-dependent gravity. Although gravitational torques prior to remnant relaxation are an order of magnitude larger than hydrodynamical torques, their intrinsic sign symmetry leads to strong cancellation; as a result, hydrodynamical and gravitational torques have a comparable effect. We also show that the debris disk’s initial specific angular momentum distribution is sharply peaked at roughly the specific angular momentum of the merged neutron star’s outer layers, a fewr_gc, and identify the regulating mechanism. 

Abstract Many studies have found that neutron star mergers leave a fraction of the stars’ mass in bound orbits surrounding the resulting massive neutron star or black hole. This mass is a site ofr-process nucleosynthesis and can generate a wind that contributes to a kilonova. However, comparatively little is known about the dynamics determining its mass or initial structure. Here we begin to investigate these questions, starting with the origin of the disk mass. Using tracer particle as well as discretized fluid data from numerical simulations, we identify where in the neutron stars the debris came from, the paths it takes in order to escape from the neutron stars’ interiors, and the times and locations at which its orbital properties diverge from those of neighboring fluid elements that end up remaining in the merged neutron star. 

Abstract Quasiperiodic erupters are a remarkable class of objects exhibiting very-large-amplitude quasiperiodic X-ray flares. Although numerous dynamical models have been proposed to explain them, relatively little attention has been given to using the properties of their radiation to constrain their dynamics. Here we show that the observed luminosity, spectrum, repetition period, duty cycle, and fluctuations in the latter two quantities point toward a model in which a main-sequence star on a moderately eccentric orbit around a supermassive black hole periodically transfers mass to the Roche lobe of the black hole; orbital dynamics lead to mildly relativistic shocks near the black hole; and thermal X-rays at the observed temperature are emitted by the gas as it flows away from the shock. Strong X-ray irradiation of the star by the flare itself augments the mass transfer, creates fluctuations in flare timing, and stirs turbulence in the stellar atmosphere that amplifies the magnetic field to a level at which magnetic stresses can accelerate infall of the transferred mass toward the black hole. 

Abstract The magnetorotational instability (MRI) has been extensively studied in circular magnetized disks, and its ability to drive accretion has been demonstrated in a multitude of scenarios. There are reasons to expect eccentric magnetized disks to also exist, but the behavior of the MRI in these disks remains largely uncharted territory. Here we present the first simulations that follow the nonlinear development of the MRI in eccentric disks. We find that the MRI in eccentric disks resembles circular disks in two ways, in the overall level of saturation and in the dependence of the detailed saturated state on magnetic topology. However, in contrast with circular disks, the Maxwell stress in eccentric disks can be negative in some disk sectors, even though the integrated stress is always positive. The angular momentum flux raises the eccentricity of the inner parts of the disk and diminishes the same of the outer parts. Because material accreting onto a black hole from an eccentric orbit possesses more energy than material tracing the innermost stable circular orbit, the radiative efficiency of eccentric disks may be significantly lower than circular disks. This may resolve the “inverse energy problem” seen in many tidal disruption events. 

Abstract Accreting supermassive binary black holes (SMBBHs) are potential multimessenger sources because they emit both gravitational-wave and electromagnetic (EM) radiation. Past work has shown that their EM output may be periodically modulated by an asymmetric density distribution in the circumbinary disk, often called an “overdensity” or “lump;” this modulation could possibly be used to identify a source as a binary. We explore the sensitivity of the overdensity to SMBBH mass ratio and magnetic flux through the accretion disk. We find that the relative amplitude of the overdensity and its associated EM periodic signal both degrade with diminishing mass ratio, vanishing altogether somewhere between 1:2 and 1:5. Greater magnetization also weakens the lump and any modulation of the light output. We develop a model to describe how lump formation results from internal stress degrading faster in the lump region than it can be rejuvenated through accretion inflow, and predicts a threshold value in specific internal stress below which lump formation should occur and which all our lump-forming simulations satisfy. Thus, detection of such a modulation would provide a constraint on both mass ratio and magnetic flux piercing the accretion flow. 

Abstract We perform a full 3D general relativistic magnetohydrodynamical (GRMHD) simulation of an equal-mass, spinning, binary black hole approaching merger, surrounded by a circumbinary disk and with a minidisk around each black hole. For this purpose, we evolve the ideal GRMHD equations on top of an approximated spacetime for the binary that is valid in every position of space, including the black hole horizons, during the inspiral regime. We use relaxed initial data for the circumbinary disk from a previous long-term simulation, where the accretion is dominated by am= 1 overdensity called the lump. We compare our new spinning simulation with a previous non-spinning run, studying how spin influences the minidisk properties. We analyze the accretion from the inner edge of the lump to the black hole, focusing on the angular momentum budget of the fluid around the minidisks. We find that minidisks in the spinning case have more mass over a cycle than the non-spinning case. However, in both cases we find that most of the mass received by the black holes is delivered by the direct plunging of material from the lump. We also analyze the morphology and variability of the electromagnetic fluxes, and we find they share the same periodicities of the accretion rate. In the spinning case, we find that the outflows are stronger than the non-spinning case. Our results will be useful to understand and produce realistic synthetic light curves and spectra, which can be used in future observations. 

Abstract We present a survey of how the spectral features of black hole X-ray binary systems depend on spin, accretion rate, viewing angle, and Fe abundance when predicted on the basis of first-principles physical calculations. The power-law component hardens with increasing spin. The thermal component strengthens with increasing accretion rate. The Compton bump is enhanced by higher accretion rate and lower spin. The Fe K α equivalent width grows sublinearly with Fe abundance. Strikingly, the K α profile is more sensitive to accretion rate than to spin because its radial surface brightness profile is relatively flat, and higher accretion rate extends the production region to smaller radii. The overall radiative efficiency is at least 30%–100% greater than as predicted by the Novikov–Thorne model.