BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Scott C. Noble

Scott C. Noble appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEgr-qcastro-ph.COastro-ph.GAastro-ph.SRastro-phastro-ph.IMhep-phphysics.comp-ph

Catalog footprint

What is connected

22works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2026arXiv

The Advanced X-ray Imaging Satellite (AXIS) Community Science Book

The AXIS Community Science Book represents the collective effort of 592 scientists worldwide to define the transformative science enabled by the Advanced X-ray Imaging Satellite (AXIS), a next-generation X-ray mission selected by NASA's Astrophysics Probe Program for Phase A study. AXIS will advance the legacy of high-angular-resolution X-ray astronomy with ~1.5'' imaging over a wide 24' field of view and an order of magnitude greater collecting area than Chandra in the 0.3-12 keV band. Combining sharp imaging, high throughput, and rapid response capabilities, AXIS will open new windows on virtually every aspect of modern astrophysics, exploring the birth and growth of supermassive black holes, the feedback processes that shape galaxies, the life cycles of stars and exoplanet environments, and the nature of compact stellar remnants, supernova remnants, and explosive transients. This book compiles 138 community-contributed science cases developed by five Science Working Groups focused on AGN and supermassive black holes, galaxy evolution and feedback, compact objects and supernova remnants, stellar physics and exoplanets, and time-domain and multi-messenger astrophysics. Together, these studies establish the scientific foundation for next-generation X-ray exploration in the 2030s and highlight strong synergies with facilities of the 2030s, such as JWST, Roman, Rubin/LSST, SKA, ALMA, ngVLA, and next-generation gravitational-wave and neutrino networks.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Advancing the Landscape of Multimessenger Science in the Next Decade

The last decade has brought about a profound transformation in multimessenger science. Ten years ago, facilities had been built or were under construction that would eventually discover the nature of objects in our universe could be detected through multiple messengers. Nonetheless, multimessenger science was hardly more than a dream. The rewards for our foresight were finally realized through IceCube's discovery of the diffuse astrophysical neutrino flux, the first observation of gravitational waves by LIGO, and the first joint detections in gravitational waves and photons and in neutrinos and photons. Today we live in the dawn of the multimessenger era. The successes of the multimessenger campaigns of the last decade have pushed multimessenger science to the forefront of priority science areas in both the particle physics and the astrophysics communities. Multimessenger science provides new methods of testing fundamental theories about the nature of matter and energy, particularly in conditions that are not reproducible on Earth. This white paper will present the science and facilities that will provide opportunities for the particle physics community renew its commitment and maintain its leadership in multimessenger science.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Electromagnetic Signatures from Supermassive Binary Black Holes Approaching Merger

We present fully relativistic predictions for the electromagnetic emission produced by accretion disks surrounding spinning and nonspinning supermassive binary black holes on the verge of merging. We use the code Bothros to post-process data from 3D general relativistic magnetohydrodynamic simulations via ray-tracing calculations. These simulations model the dynamics of a circumbinary disk and the mini-disks that form around two equal-mass black holes orbiting each other at an initial separation of 20 gravitational radii, and evolve the system for more than 10 orbits in the inspiral regime. We model the emission as the sum of thermal blackbody radiation emitted by an optically thick accretion disk and a power-law spectrum extending to hard X-rays emitted by a hot optically thin corona. We generate time-dependent spectra, images, and light curves at various frequencies to investigate intrinsic periodic signals in the emission, as well as the effects of the black hole spin. We find that prograde black hole spin makes mini-disks brighter since the smaller innermost stable circular orbit angular momentum demands more dissipation before matter plunges to the horizon. However, compared to mini-disks in larger separation binaries with spinning black holes, our mini-disks are less luminous: unlike those systems, their mass accretion rate is lower than in the circumbinary disk, and they radiate with lower efficiency because their inflow times are shorter. Compared to a single black hole system matched in mass and accretion rate, these binaries have spectra noticeably weaker and softer in the UV. Finally, we discuss the implications of our findings for the potential observability of these systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

HARM3D+NUC: A new method for simulating the post-merger phase of binary neutron star mergers with GRMHD, tabulated EOS and neutrino leakage

The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mini-disk accretion onto spinning black hole binaries: quasi-periodicities and outflows

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 mini-disks 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 an $m=1$ overdensity called the lump. We compare our new spinning simulation with a previous non-spinning run, studying how spin influences the mini-disk 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 mini-disks. We find that mini-disks in the spinning case have more mass over a cycle than the non-spinning case. However, in both cases, we find 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 $8$ times 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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

PATOKA: Simulating Electromagnetic Observables of Black Hole Accretion

The Event Horizon Telescope (EHT) has released analyses of reconstructed images of horizon-scale millimeter emission near the supermassive black hole at the center of the M87 galaxy. Parts of the analyses made use of a large library of synthetic black hole images and spectra, which were produced using numerical general relativistic magnetohydrodynamics fluid simulations and polarized ray tracing. In this article, we describe the PATOKA pipeline, which was used to generate the Illinois contribution to the EHT simulation library. We begin by describing the relevant accretion systems and radiative processes. We then describe the details of the three numerical codes we use, iharm, ipole, and igrmonty, paying particular attention to differences between the current generation of the codes and the originally published versions. Finally, we provide a brief overview of simulated data as produced by PATOKA and conclude with a discussion of limitations and future directions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

PatchworkWave: A Multipatch Infrastructure for Multiphysics/Multiscale/Multiframe/Multimethod Simulations at Arbitrary Order

We present an extension of the PatchworkMHD code [1], itself an MHD-capable extension of the Patchwork code [2], for which several algorithms presented here were co-developed. Its purpose is to create a multipatch scheme compatible with numerical simulations of arbitrary equations of motion at any discretization order in space and time. In the Patchwork framework, the global simulation is comprised of an arbitrary number of moving, local meshes, or patches, which are free to employ their own resolution, coordinate system/topology, physics equations, reference frame, and in our new approach, numerical method. Each local patch exchanges boundary data with a single global patch on which all other patches reside through a client-router-server parallelization model. In generalizing Patchwork to be compatible with arbitrary order time integration, PatchworkMHD and PatchworkWave have significantly improved the interpatch interpolation accuracy by removing an interpolation of interpolated data feedback present in the original Patchwork code. Furthermore, we extend Patchwork to be multimethod by allowing multiple state vectors to be updated simultaneously, with each state vector providing its own interpatch interpolation and transformation procedures. As such, our scheme is compatible with nearly any set of hyperbolic partial differential equations. We demonstrate our changes through the implementation of a scalar wave toy-model that is evolved on arbitrary, time dependent patch configurations at 4th order accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tidal radii of main sequence stars -- III. Partial disruptions

In this paper, the third in this series, we continue our study of tidal disruption events of main-sequence stars by a non-spinning $10^{6}~\rm{M}_\odot$ supermassive black hole. Here we focus on the outcomes of partial disruptions. As the encounter becomes weaker, the debris mass is increasingly concentrated near the outer edges of the energy distribution. As a result, the mass fallback rate can deviate substantially from a $t^{-5/3}$ power-law, becoming more like a single peak with a tail declining as $t^{-p}$ with $p\simeq2-5$. Surviving remnants are spun-up in the prograde direction and are hotter than MS stars of the same mass. Their specific orbital energy is $\simeq10^{-3}\times$ that of the debris (but of either sign with respect to the black hole potential) while their specific angular momentum is close to that of the original star. Even for strong encounters, remnants have speeds at infinity relative to the black hole potential $\lesssim 300$ km s$^{-1}$, so they are unable to travel far out into the galactic bulge. Remnants bound to the black hole can possibly go through a second tidal disruption event.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Disk Emission from Magneto-hydrodynamic Simulations of Spinning Black Holes

We present the results of a new series of global 3D relativistic magneto-hydrodynamic (MHD) simulations of thin accretion disks around spinning black holes. The disks have aspect ratios of $H/R\sim 0.05$ and spin parameters $a/M=0, 0.5, 0.9$, and $0.99$. Using the ray-tracing code Pandurata, we generate broad-band thermal spectra and polarization signatures from the MHD simulations. We find that the simulated spectra can be well fit with a simple, universal emissivity profile that better reproduces the behavior of the emission from the inner disk, compared to traditional analyses carried out using a Novikov-Thorne thin disk model. Lastly, we show how spectropolarization observations can be used to convincingly break the spin-inclination degeneracy well-known to the continuum fitting method of measuring black hole spin.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Driven neutron star collapse: Type~I critical phenomena and the initial black hole mass distribution

We study the general relativistic collapse of neutron star (NS) models in spherical symmetry. Our initially stable models are driven to collapse by the addition of one of two things: an initially in-going velocity profile, or a shell of minimally coupled, massless scalar field that falls onto the star. Tolman-Oppenheimer-Volkoff (TOV) solutions with an initially isentropic, gamma-law equation of state serve as our NS models. The initial values of the velocity profile's amplitude and the star's central density span a parameter space which we have surveyed extensively and which we find provides a rich picture of the possible end states of NS collapse. This parameter space survey elucidates the boundary between Type I and Type II critical behavior in perfect fluids which coincides, on the subcritical side, with the boundary between dispersed and bound end states. For our particular model, initial velocity amplitudes greater than 0.3c are needed to probe the regime where arbitrarily small black holes can form. In addition, we investigate Type I behavior in our system by varying the initial amplitude of the initially imploding scalar field. In this case we find that the Type I critical solutions resemble TOV solutions on the 1-mode unstable branch of equilibrium solutions, and that the critical solutions' frequencies agree well with the fundamental mode frequencies of the unstable equilibria. Additionally, the critical solution's scaling exponent is shown to be well approximated by a linear function of the initial star's central density.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

General Relativistic Hydrodynamic Simulation of Accretion Flow from a Stellar Tidal Disruption

We study how the matter dispersed when a supermassive black hole tidally disrupts a star joins an accretion flow. Combining a relativistic hydrodynamic simulation of the stellar disruption with a relativistic hydrodynamics simulation of the tidal debris motion, we track such a system until ~80% of the stellar mass bound to the black hole has settled into an accretion flow. Shocks near the stellar pericenter and also near the apocenter of the most tightly-bound debris dissipate orbital energy, but only enough to make the characteristic radius comparable to the semi-major axis of the most-bound material, not the tidal radius as previously thought. The outer shocks are caused by post-Newtonian effects, both on the stellar orbit during its disruption and on the tidal forces. Accumulation of mass into the accretion flow is non-monotonic and slow, requiring ~3--10x the orbital period of the most tightly-bound tidal streams, while the inflow time for most of the mass may be comparable to or longer than the mass accumulation time. Deflection by shocks does, however, remove enough angular momentum and energy from some mass for it to move inward even before most of the mass is accumulated into the accretion flow. Although the accretion rate rises sharply and then decays roughly as a power-law, its maximum is ~0.1x the previous expectation, and the duration of the peak is ~5x longer than previously predicted. The geometric mean of the black hole mass and stellar mass inferred from a measured event timescale is therefore ~0.2x the value given by classical theory.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Resolving the relative influence of strong field spacetime dynamics and MHD on circumbinary disk physics

In this paper we evolve magnetized and unmagnetized circumbinary accretion disks around supermassive black hole binaries in the relativistic regime. We use a post-Newtonian expansion to construct an analytical spacetime and determine how the order of the post-Newtonian (PN) expansion affects the dynamics of the gas. We find very small differences in the late-time bulk dynamics of non-magnetized hydrodynamic evolutions between the two spacetimes down to separations of approximately $40GM/c^2$ where $M$ is the total mass of the binary. For smaller separations, the differences due to PN-order become comparable to differences caused by using initial data further from equilibrium. For magnetized gas, MHD stresses, which drives the accretion dynamics, tends to mask all higher order PN effects even at separations of $20GM/c^2$, leading to essentially the same observed electromagnetic luminosity. This implies that our calculations of the EM signal may be robust down to small binary separations. Our investigation is the first to demonstrate how the level of PN accuracy affects a circumbinary disk's evolution and informs us of the range in separation within which to trust the PN approximation for this kind of study. We also address the influence the initial conditions and binary separation have on simulation predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Approximate black hole binary spacetime via asymptotic matching

We construct a fully analytic, general relativistic, nonspinning black hole binary spacetime that approximately solves the vacuum Einstein equations everywhere in space and time for black holes sufficiently well separated. The metric is constructed by asymptotically matching perturbed Schwarzschild metrics near each black hole to a two-body post-Newtonian metric far from them, and a two-body post-Minkowskian metric farther still. Asymptotic matching is done without linearizing about a particular time slice, and thus it is valid dynamically and for all times, provided the binary is sufficiently well separated. This approximate global metric can be used for long dynamical evolutions of relativistic magnetohydrodynamical, circumbinary disks around inspiraling supermassive black holes to study a variety of phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamic Fisheye Grids for Binary Black Holes Simulations

We present a new warped gridding scheme adapted to simulating gas dynamics in binary black hole spacetimes. The grid concentrates grid points in the vicinity of each black hole to resolve the smaller scale structures there, and rarefies grid points away from each black hole to keep the overall problem size at a practical level. In this respect, our system can be thought of as a "double" version of the fisheye coordinate system, used before in numerical relativity codes for evolving binary black holes. The gridding scheme is constructed as a mapping between a uniform coordinate system--in which the equations of motion are solved--to the distorted system representing the spatial locations of our grid points. Since we are motivated to eventually use this system for circumbinary disc calculations, we demonstrate how the distorted system can be constructed to asymptote to the typical spherical polar coordinate system, amenable to efficiently simulating orbiting gas flows about central objects with little numerical diffusion. We discuss its implementation in the Harm3d code, tailored to evolve the magnetohydrodynamics equations in curved spacetimes. We evaluate the performance of the system's implementation in Harm3d with a series of tests, such as the advected magnetic field loop test, magnetized Bondi accretion, and evolutions of hydrodynamic discs about a single black hole and about a binary black hole. Like we have done with Harm3d, this gridding scheme can be implemented in other unigrid codes as a (possibly) simpler alternative to adaptive mesh refinement.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

GRHydro: A new open source general-relativistic magnetohydrodynamics code for the Einstein Toolkit

We present the new general-relativistic magnetohydrodynamics (GRMHD) capabilities of the Einstein Toolkit, an open-source community-driven numerical relativity and computational relativistic astrophysics code. The GRMHD extension of the Toolkit builds upon previous releases and implements the evolution of relativistic magnetised fluids in the ideal MHD limit in fully dynamical spacetimes using the same shock-capturing techniques previously applied to hydrodynamical evolution. In order to maintain the divergence-free character of the magnetic field, the code implements both hyperbolic divergence cleaning and constrained transport schemes. We present test results for a number of MHD tests in Minkowski and curved spacetimes. Minkowski tests include aligned and oblique planar shocks, cylindrical explosions, magnetic rotors, Alfvén waves and advected loops, as well as a set of tests designed to study the response of the divergence cleaning scheme to numerically generated monopoles. We study the code's performance in curved spacetimes with spherical accretion onto a black hole on a fixed background spacetime and in fully dynamical spacetimes by evolutions of a magnetised polytropic neutron star and of the collapse of a magnetised stellar core. Our results agree well with exact solutions where these are available and we demonstrate convergence. All code and input files used to generate the results are available on http://einsteintoolkit.org. This makes our work fully reproducible and provides new users with an introduction to applications of the code.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

X-ray Spectra from MHD Simulations of Accreting Black Holes

We present the results of a new global radiation transport code coupled to a general relativistic magneto-hydrodynamic simulation of an accreting, non-rotating black hole. For the first time, we are able to explain from first principles in a self-consistent way all the components seen in the X-ray spectra of stellar-mass black holes, including a thermal peak and all the features associated with strong hard X-ray emission: a power-law extending to high energies, a Compton reflection hump, and a broad iron line. Varying only the mass accretion rate, we are able to reproduce a wide range of X-ray states seen in most galactic black hole sources. The temperature in the corona is T_e ~ 10 keV in a boundary layer near the disk and rises smoothly to T_e >~ 100 keV in low-density regions far above the disk. Even as the disk's reflection edge varies from the horizon out to ~ 6M as the accretion rate decreases, we find that the shape of the Fe Kαline is remarkably constant. This is because photons emitted from the plunging region are strongly beamed into the horizon and never reach the observer. We have also carried out a basic timing analysis of the spectra and find that the fractional variability increases with photon energy and viewer inclination angle, consistent with the coronal hot spot model for X-ray fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Circumbinary MHD Accretion into Inspiraling Binary Black Holes

As 2 black holes bound to each other in a close binary approach merger their inspiral time becomes shorter than the characteristic inflow time of surrounding orbiting matter. Using an innovative technique in which we represent the changing spacetime in the region occupied by the orbiting matter with a 2.5PN approximation and the binary orbital evolution with 3.5PN, we have simulated the MHD evolution of a circumbinary disk surrounding an equal-mass non-spinning binary. Prior to the beginning of the inspiral, the structure of the circumbinary disk is predicted well by extrapolation from Newtonian results. The binary opens a low-density gap whose radius is roughly two binary separations, and matter piles up at the outer edge of this gap as inflow is retarded by torques exerted by the binary; nonetheless, the accretion rate is diminished relative to its value at larger radius by only about a factor of 2. During inspiral, the inner edge of the disk at first moves inward in coordination with the shrinking binary, but as the orbital evolution accelerates, the rate at which the inner edge moves toward smaller radii falls behind the rate of binary compression. In this stage, the rate of angular momentum transfer from the binary to the disk slows substantially, but the net accretion rate decreases by only 10-20%. When the binary separation is tens of gravitational radii, the rest-mass efficiency of disk radiation is a few percent, suggesting that supermassive binary black holes in galactic nuclei could be very luminous at this stage of their evolution. If the luminosity were optically thin, it would be modulated at a frequency that is a beat between the orbital frequency of the disk's surface density maximum and the binary orbital frequency. However, a disk with sufficient surface density to be luminous should also be optically thick; as a result, the periodic modulation may be suppressed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Near-Infrared and X-ray Quasi-Periodic Oscillations in Numerical Models of Sgr A*

We report transient quasi-periodic oscillations (QPOs) on minute timescales in relativistic, radiative models of the galactic center source Sgr A*. The QPOs result from nonaxisymmetric $m=1$ structure in the accretion flow excited by MHD turbulence. Near-infrared (NIR) and X-ray power spectra show significant peaks at frequencies comparable to the orbital frequency at the innermost stable circular orbit (ISCO) $f_o$. The excess power is associated with inward propagating magnetic filaments inside the ISCO. The amplitudes of the QPOs are sensitive to the electron distribution function. We argue that transient QPOs appear at a range of frequencies in the neighborhood of $f_o$ and that the power spectra, averaged over long times, likely show a broad bump near $f_o$ rather than distinct, narrow QPO features.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Global GRMHD Simulations of Black Hole Accretion Flows: a Convergence Study

Global, general relativistic magnetohydrodynamic (GRMHD) simulations of nonradiative, magnetized disks are widely used to model accreting black holes. We have performed a convergence study of GRMHD models computed with HARM3D. The models span a factor of 4 in linear resolution, from 96x96x64 to 384x384x256. We consider three diagnostics of convergence: (1) dimensionless shell-averaged quantities such as plasma β; (2) the azimuthal correlation length of fluid variables; and (3) synthetic spectra of the source including synchrotron emission, absorption, and Compton scattering. Shell-averaged temperature is, except for the lowest resolution run, nearly independent of resolution; shell-averaged plasma β decreases steadily with resolution but shows signs of convergence. The azimuthal correlation lengths of density, internal energy, and temperature decrease steadily with resolution but show signs of convergence. In contrast, the azimuthal correlation length of magnetic field decreases nearly linearly with grid size. We argue by analogy with local models, however, that convergence should be achieved with another factor of 2 in resolution. Synthetic spectra are, except for the lowest resolution run, nearly independent of resolution. The convergence behavior is consistent with that of higher physical resolution local model (shearing box) calculations and with the recent nonrelativistic global convergence studies of Hawley et al. (2011).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Radiative efficiency and thermal spectrum of accretion onto Schwarzschild black holes

Recent general relativistic magneto-hydrodynamic (MHD) simulations of accretion onto black holes have shown that, contrary to the basic assumptions of the Novikov-Thorne model, there can be substantial magnetic stress throughout the plunging region. Additional dissipation and radiation can therefore be expected. We use data from a particularly well-resolved simulation of accretion onto a non-spinning black hole to compute both the radiative efficiency of such a flow and its spectrum if all emitted light is radiated with a thermal spectrum whose temperature matches the local effective temperature. This disk is geometrically thin enough (H/r ~= 0.06) that little heat is retained in the flow. In terms of light reaching infinity (i.e., after allowance for all relativistic effects and for photon capture by the black hole), we find that the radiative efficiency is at least ~=6-10% greater than predicted by the Novikov-Thorne model (complete radiation of all heat might yield another ~6%). We also find that the spectrum more closely resembles the Novikov-Thorne prediction for a/M ~= 0.2--0.3 than for the correct value, a/M=0. As a result, if the spin of a non-spinning black hole is inferred by model-fitting to a Novikov-Thorne model with known black hole mass, distance, and inclination, the inferred a/M is too large by ~= 0.2--0.3.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dependence of inner accretion disk stress on parameters: the Schwarzschild case

We explore the parameter dependence of inner disk stress in black hole accretion by contrasting the results of a number of simulations, all employing 3-d general relativistic MHD in a Schwarzschild spacetime. Five of these simulations were performed with the intrinsically conservative code HARM3D, which allows careful regulation of the disk aspect ratio, H/R; our simulations span a range in H/R from 0.06 to 0.17. We contrast these simulations with two previously reported simulations in a Schwarzschild spacetime in order to investigate possible dependence of the inner disk stress on magnetic topology. In all cases, much care was devoted to technical issues: ensuring adequate resolution and azimuthal extent, and averaging only over those time-periods when the accretion flow is in approximate inflow equilibrium. We find that the time-averaged radial-dependence of fluid-frame electromagnetic stress is almost completely independent of both disk thickness and poloidal magnetic topology. It rises smoothly inward at all radii (exhibiting no feature associated with the ISCO) until just outside the event horizon, where the stress plummets to zero. Reynolds stress can also be significant near the ISCO and in the plunging region; the magnitude of this stress, however, depends on both disk thickness and magnetic topology. The two stresses combine to make the net angular momentum accreted per unit rest-mass 7-15% less than the angular momentum of the ISCO.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Direct calculation of the radiative efficiency of an accretion disk around a black hole

Numerical simulation of magnetohydrodynamic (MHD) turbulence makes it possible to study accretion dynamics in detail. However, special effort is required to connect inflow dynamics (dependent largely on angular momentum transport) to radiation (dependent largely on thermodynamics and photon diffusion). To this end we extend the flux-conservative, general relativistic MHD code HARM from axisymmetry to full 3D. The use of an energy conserving algorithm allows the energy dissipated in the course of relativistic accretion to be captured as heat. The inclusion of a simple optically thin cooling function permits explicit control of the simulated disk's geometric thickness as well as a direct calculation of both the amplitude and location of the radiative cooling associated with the accretion stresses. Fully relativistic ray-tracing is used to compute the luminosity received by distant observers. For a disk with aspect ratio H/r ~ 0.1 accreting onto a black hole with spin parameter a/M = 0.9, we find that there is significant dissipation beyond that predicted by the classical Novikov-Thorne model. However, much of it occurs deep in the potential, where photon capture and gravitational redshifting can strongly limit the net photon energy escaping to infinity. In addition, with these parameters and this radiation model, significant thermal and magnetic energy remains with the gas and is accreted by the black hole. In our model, the net luminosi ty reaching infinity is 6% greater than the Novikov-Thorne prediction. If the accreted thermal energy were wholly radiated, the total luminosity of the accretion flow would be ~20% greater than the Novikov-Thorne value.

0 citations0 reviewsNo code linkedOpen access