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

P. Chris Fragile

P. Chris Fragile appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEastro-ph.GAastro-phastro-ph.IMastro-ph.SRgr-qc

Catalog footprint

What is connected

30works
6topics
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

30 published item(s)

preprint2020arXiv

Interactions of Type I X-ray Bursts with Thin Accretion Discs

We perform a set of numerical experiments studying the interaction of Type I X-ray bursts with thin, Shakura-Sunyaev type accretion discs. Careful observations of X-ray spectra during such bursts have hinted at changes occurring in the inner regions of the disc. We now clearly demonstrate a number of key effects that take place simultaneously, including: evidence for weak, radiation-driven outflows along the surface of the disc; significant levels of Poynting-Robertson (PR) drag, leading to enhanced accretion; and prominent heating in the disc, which increases the height, while lowering the density and optical depth. The PR drag causes the inner edge of the disc to retreat from the neutron star surface toward larger radii and then recover on the timescale of the burst. We conclude that the rich interaction of an X-ray burst with the surrounding disc provides a novel way to study the physics of accretion onto compact objects.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multi-Frequency General Relativistic Radiation-Hydrodynamics with $\bf{M}_1$ Closure

We report on recent upgrades to our general relativistic radiation-magnetohydrodynamics code, Cosmos++, which expands the two-moment, $\bf{M}_1$, radiation treatment from grey to multi-frequency transport, including Doppler and gravitational frequency shifts. The solver accommodates either photon (Bose-Einstein) or neutrino (Fermi-Dirac) statistical distribution functions with absorption, emission, and elastic scattering processes. An implicit scheme is implemented to simultaneously solve the primitive inversion problem together with the radiation-matter coupling source terms, providing stability over a broad range of opacities and optical depths where the interactions terms can be stiff. We discuss our formulations and numerical methods, and validate our methods against a wide variety of test problems spanning optically thin to thick regimes in flat, weakly curved, and strongly curved spacetimes.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Neutron Star QPOs from Oscillating, Precessing Hot, Thick Flow

Across black hole (BH) and neutron star (NS) low-mass X-ray binaries (LMXBs), there appears to be some correlation between certain high- and low-frequency quasi-periodic oscillations (QPOs). In a previous paper, we showed that for BH LMXBs, this could be explained by the simultaneous oscillation and precession of a hot, thick, torus-like corona. In the current work, we extend this idea to NS LMXBs by associating the horizontal branch oscillations (HBO) with precession and the upper-kiloHertz (ukHz) QPO with vertical epicyclic motion. For the Atoll source 4U 1608-52, the model can match many distinct, simultaneous observations of the HBO and ukHz QPO by varying the inner and outer radius of the torus, while maintaining fixed values for the mass (M_{NS}) and spin (a_*) of the neutron star. The best fit values are M_{NS} = 1.38 \pm 0.03 M_\odot and a_* = 0.325 \pm 0.005. By combining these constraints with the measured spin frequency, we are able to obtain an estimate for the moment of inertia of I_{NS} = 1.40 \pm 0.02 \times 10^{45} g cm^2, which places constraints on the equation of state. The model is unable to fit the lower-kHz QPO, but evidence suggests that QPO may be associated with the boundary layer between the accretion flow and the neutron star surface, which is not treated in this work.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Relativistic, Viscous, Radiation Hydrodynamic Simulations of Geometrically Thin Disks. II. Disk Variability

We perform detailed variability analysis of two-dimensional viscous, radiation hydrodynamic numerical simulations of Shakura-Sunyaev thin disks around a stellar mass black hole. Disk models are initialized on both the gas-, as well as radiation-, pressure-dominated branches of the thermal equilibrium curve, with mass accretion rates spanning the range from $\dot{M} = 0.01 L_\mathrm{Edd}/c^2$ to $10 L_\mathrm{Edd}/c^2$. An analysis of temporal variations of the numerically simulated disk reveals multiple robust, coherent oscillations. Considering the local mass flux variability, we find an oscillation occurring at the maximum radial epicyclic frequency, $3.5\times 10^{-3}\,t_\mathrm{g}^{-1}$, a possible signature of a trapped fundamental ${\it g}$-mode. Although present in each of our simulated models, the trapped ${\it g}$-mode feature is most prominent in the gas-pressure-dominated case. The total pressure fluctuations in the disk suggest strong evidence for standing-wave ${\it p}$-modes, some trapped in the inner disk close to the ISCO, others present in the middle/outer parts of the disk. Knowing that the trapped ${\it g}$-mode frequency and maximum radial epicyclic frequency differ by only $0.01\%$ in the case of a non-rotating black hole, we simulated an additional initially gas-pressure-dominated disk with a dimensionless black hole spin parameter $a_* = 0.5$. The oscillation frequency in the spinning black hole case confirms that this oscillation is indeed a trapped ${\it g}$-mode. All the numerical models we report here also show a set of high frequency oscillations at the vertical epicyclic and breathing mode frequencies. The vertical oscillations show a 3:2 frequency ratio with oscillations occurring approximately at the radial epicyclic frequency, which could be of astrophysical importance in observed twin peak, high-frequency quasi-periodic oscillations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

High-Frequency and Type-C QPOs from Oscillating, Precessing Hot, Thick Flow

Motivated by recent studies showing an apparent correlation between the high-frequency quasi-periodic oscillations (QPOs) and the low-frequency, type-C QPO in low-mass, black hole X-ray binaries (LMXBs), we explore a model that explains all three QPOs in terms of an oscillating, precessing hot flow in the truncated-disk geometry. Our model favors attributing the two high-frequency QPOs, often occurring in a near 3:2 frequency ratio, to the breathing and vertical epicyclic frequency modes of the hot, thick flow, although we can not rule out the Keplerian and m=-1 radial epicyclic modes. In either case, the type-C QPO is attributed to precession. The correlation of the QPOs comes from the fact that all three frequencies are associated with the same geometrical structure. While the exact QPO frequencies are sensitive to the black hole mass and spin, their evolution over the course of an outburst is mainly tied to the truncation radius between the geometrically thin, optically thick disk and the inner, hot flow. We show that, in the case of the LMXB GRO J1655-40, this model can explain the one simultaneous observation of all three QPOs and that an extrapolation of the model appears to match lower frequency observations where only two of the three components are seen. Thus, this model may be able to unify multiple QPO observations using the properties of a single, simple, geometrical model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Radio and Millimeter Monitoring of Sgr A*: Spectrum, Variability, and Constraints on the G2 Encounter

We report new observations with the Very Large Array, Atacama Large Millimeter Array, and Submillimeter Array at frequencies from 1.0 to 355 GHz of the Galactic Center black hole, Sagittarius A*. These observations were conducted between October 2012 and November 2014. While we see variability over the whole spectrum with an amplitude as large as a factor of 2 at millimeter wavelengths, we find no evidence for a change in the mean flux density or spectrum of Sgr A* that can be attributed to interaction with the G2 source. The absence of a bow shock at low frequencies is consistent with a cross-sectional area for G2 that is less than $2 \times 10^{29}$ cm$^2$. This result fits with several model predictions including a magnetically arrested cloud, a pressure-confined stellar wind, and a stellar photosphere of a binary merger. There is no evidence for enhanced accretion onto the black hole driving greater jet and/or accretion flow emission. Finally, we measure the millimeter wavelength spectral index of Sgr A* to be flat; combined with previous measurements, this suggests that there is no spectral break between 230 and 690 GHz. The emission region is thus likely in a transition between optically thick and thin at these frequencies and requires a mix of lepton distributions with varying temperatures consistent with stratification.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Conservative GRMHD Simulations of Moderately Thin, Tilted Accretion Disks

This paper presents our latest numerical simulations of accretion disks that are misaligned with respect to the rotation axis of a Kerr black hole. In this work we use a new, fully conservative version of the Cosmos++ general relativistic magnetohydrodynamics (GRMHD) code, coupled with an ad hoc cooling function designed to control the thickness of the disk. Together these allow us to simulate the thinnest tilted accretion disks ever using a GRMHD code. In this way, we are able to probe the regime where the dimensionless stress and scale height of the disk become comparable. We present results for both prograde and retrograde cases. The simulated prograde tilted disk shows no sign of Bardeen-Petterson alignment even in the innermost parts of the disk. The simulated retrograde tilted disk, however, does show modest alignment. The implication of these results is that the parameter space associated with Bardeen-Petterson alignment for prograde disks may be rather small, only including very thin disks. Unlike our previous work, we find no evidence for standing shocks in our simulated tilted disks. We ascribe this to the combination of small black hole spin, small tilt angle, and small disk scale height in these simulations. We also add to the growing body of literature pointing out that the turbulence driven by the magnetorotational instability in global simulations of accretion disks is not isotropic. Finally, we provide a comparison between our moderately thin, untilted reference simulation and other numerical simulations of thin disks in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Local Stability of Strongly Magnetized Black Hole Tori

We investigate the local stability of strongly magnetized relativistic tori orbiting Kerr black holes, for the case of a purely toroidal magnetic field topology. Our approach encompasses both solving the full dispersion relation numerically and solving a simplified analytic treatment. Both approaches indicate that such tori are subject to an unstable non-axisymmetric magnetorotational mode, regardless of such background properties as the gas-to-magnetic pressure ratio and specific angular momentum distribution. We demonstrate that our modal analysis matches well with what is seen in global, three-dimensional, general relativistic magnetohydrodynamic simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

No Evidence for Bardeen-Petterson Alignment in GRMHD Simulations and Semi-Analytic Models of Moderately Thin, Prograde, Tilted Accretion Disks

In this paper we introduce the first results that use data extracted directly from numerical simulations as inputs to the analytic twisted disk model of Zhuravlev & Ivanov. In both numerical and analytic approaches, fully relativistic models of tilted and twisted disks having a moderate effective viscosity around a slowly rotating Kerr black hole are considered. Qualitatively, the analytic model demonstrates the same dynamics as the simulations, although with some quantitative offset. Namely, the GRMHD simulations generally give smaller variations of tilt and twist across the disk. When the black hole and the disk rotate in the same sense, the simulated tilted disk and analytic model show no sign of Bardeen-Petterson alignment, even in the innermost parts of the disk, where the characteristic time for relaxation to a quasi-stationary configuration is of the same order as the computation time. In the opposite case, when the direction of the disk's rotation is opposite to that of black hole, a partial alignment is observed, in agreement with previous theoretical estimates. Thus, both fully numerical and analytic schemes demonstrate that the Bardeen-Petterson effect may not be possible for the case of prograde rotation provided that disk's effective viscosity is sufficiently small. This may have implications in modeling of different astrophysical phenomena, such as disk spectra and jet orientation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical Simulations of Optically Thick Accretion onto a Black Hole - II. Rotating Flow

In this paper we report on recent upgrades to our general relativistic radiation magnetohydrodynamics code, Cosmos++, including the development of a new primitive inversion scheme and a hybrid implicit-explicit solver with a more general closure relation for the radiation equations. The new hybrid solver helps stabilize the treatment of the radiation source terms, while the new closure allows for a much broader range of optical depths to be considered. These changes allow us to expand by orders of magnitude the range of temperatures, opacities, and mass accretion rates, and move a step closer toward our goal of performing global simulations of radiation-pressure-dominated black hole accretion disks. In this work we test and validate the new method against an array of problems. We also demonstrate its ability to handle super-Eddington, quasi-spherical accretion. Even with just a single proof-of-principle simulation, we already see tantalizing hints of the interesting phenomenology associated with the coupling of radiation and gas in super-Eddington accretion flows.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Three-Dimensional Moving-Mesh Simulations of Galactic Center Cloud G2

Using three-dimensional, moving-mesh simulations, we investigate the future evolution of the recently discovered gas cloud G2 traveling through the galactic center. We consider the case of a spherical cloud initially in pressure equilibrium with the background. Our suite of simulations explores the following parameters: the equation of state, radial profiles of the background gas, and start times for the evolution. Our primary focus is on how the fate of this cloud will affect the future activity of Sgr A*. From our simulations we expect an average feeding rate in the range of $5-19 \times 10^{-8}$ solar masses per year beginning in 2013 and lasting for at least 7 years (our simulations stop in year 2020). The accretion varies by less than a factor of three on timescales <1 month, and shows no more than a factor of 10 difference between the maximum and minimum. These rates are comparable to the current estimated accretion rate in the immediate vicinity of Sgr A*, although they represent only a small (<5%) increase over the current expected feeding rate at the effective inner boundary of our simulations. Therefore, the break up of cloud G2 may have only a minimal effect on the brightness and variability of Sgr A* over the next decade. This is because current models of the galactic center predict that most of the gas will be caught up in outflows. However, if the accreted G2 material can remain cold, it may not mix well with the hot, diffuse background gas, and instead accrete efficiently onto Sgr A*. Further observations of G2 will give us an unprecedented opportunity to test this idea. The break up of the cloud itself may also be observable. By tracking the amount of cloud energy that is dissipated during our simulations, we are able to get a rough estimate of the luminosity associated with its tidal disruption; we find values of a few $10^{36}$ erg/s.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Update on 3D Moving Mesh Simulations of Galactic Center Cloud G2

Using three-dimensional, moving-mesh simulations, we investigate the future evolution of the recently discovered gas cloud G2 traveling through the galactic center. From our simulations we expect an average feeding rate onto Sgr A* in the range of $(5-19) \times 10^{-8} M_\odot\mathrm{~yr}^{-1}$ beginning in 2014. This accretion varies by less than a factor of three on timescales of about 1 month, and shows no more than a factor of 10 difference between the maximum and minimum observed rates within any given model. These rates are comparable to the current estimated accretion rate in the immediate vicinity of Sgr A*, although they represent only a small (<10%) increase over the current expected feeding rate at the effective inner boundary of our simulations $(r_\mathrm{acc} = 750 R_S \sim 10^{15} \mathrm{cm})$. We also explore multiple possible equations of state to describe the gas. In examining the Br-$γ$ light curves produced from our simulations, we find that all of our isothermal models predict significant (factor of 10) enhancements in the luminosity of G2 as it approaches pericenter, in conflict with observations. Models that instead allow the cloud to heat as it is compressed do better at matching observations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Current Status of Simulations

As the title suggests, the purpose of this chapter is to review the current status of numerical simulations of black hole accretion disks. This chapter focuses exclusively on global simulations of the accretion process within a few tens of gravitational radii of the black hole. Most of the simulations discussed are performed using general relativistic magnetohydrodynamic (MHD) schemes, although some mention is made of Newtonian radiation MHD simulations and smoothed particle hydrodynamics. The goal is to convey some of the exciting work that has been going on in the past few years and provide some speculation on future directions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Foundations of Black Hole Accretion Disk Theory

This review covers the main aspects of black hole accretion disk theory. We begin with the view that one of the main goals of the theory is to better understand the nature of black holes themselves. In this light we discuss how accretion disks might reveal some of the unique signatures of strong gravity: the event horizon, the innermost stable circular orbit, and the ergosphere. We then review, from a first-principles perspective, the physical processes at play in accretion disks. This leads us to the four primary accretion disk models that we review: Polish doughnuts (thick disks), Shakura-Sunyaev (thin) disks, slim disks, and advection-dominated accretion flows (ADAFs). After presenting the models we discuss issues of stability, oscillations, and jets. Following our review of the analytic work, we take a parallel approach in reviewing numerical studies of black hole accretion disks. We finish with a few select applications that highlight particular astrophysical applications: measurements of black hole mass and spin, black hole vs. neutron star accretion disks, black hole accretion disk spectral states, and quasi-periodic oscillations (QPOs).

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Physical Properties of the Inner Shocks in Hot, Tilted Black Hole Accretion Flows

Simulations of hot, pressure supported, tilted black hole accretion flows, in which the angular momentum of the flow is misaligned with the black hole spin axis, can exhibit two non-axisymmetric shock structures in the inner regions of the flow. We analyze the strength and significance of these shock structures in simulations with tilt angles of 10 and 15 degrees. By integrating fluid trajectories in the simulations through the shocks, and tracking the variations of fluid quantities along these trajectories, we show that these shocks are strong, with substantial compression ratios, in contrast to earlier claims. However, they are only moderately relativistic. We also show that the two density enhancements resembling flow streams in their shape are in fact merely post-shock compressions, as fluid trajectories cut across, rather than flow along, them. The dissipation associated with the shocks is a substantial fraction ($\simeq3-12$ percent) of the rest mass energy advected into the hole, and therefore comparable to the dissipation expected from turbulence. The shocks should therefore make order unity changes in the observed properties of black hole accretion flows that are tilted.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

GRMHD simulations of accretion onto Sgr A*: How important are radiative losses?

We present general relativistic magnetohydrodynamic (GRMHD) numerical simulations of the accretion flow around the supermassive black hole in the Galactic centre, Sagittarius A* (Sgr A*). The simulations include for the first time radiative cooling processes (synchrotron, bremsstrahlung, and inverse Compton) self-consistently in the dynamics, allowing us to test the common simplification of ignoring all cooling losses in the modeling of Sgr A*. We confirm that for Sgr A*, neglecting the cooling losses is a reasonable approximation if the Galactic centre is accreting below ~10^{-8} Msun/yr i.e. Mdot < 10^{-7} Mdot_Edd. But above this limit, we show that radiative losses should be taken into account as significant differences appear in the dynamics and the resulting spectra when comparing simulations with and without cooling. This limit implies that most nearby low-luminosity active galactic nuclei are in the regime where cooling should be taken into account. We further make a parameter study of axisymmetric gas accretion around the supermassive black hole at the Galactic centre. This approach allows us to investigate the physics of gas accretion in general, while confronting our results with the well studied and observed source, Sgr A*, as a test case. We confirm that the nature of the accretion flow and outflow is strongly dependent on the initial geometry of the magnetic field. For example, we find it difficult, even with very high spins, to generate powerful outflows from discs threaded with multiple, separate poloidal field loops.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

No Correlation Between Disc Scale-Height and Jet Power in GRMHD Simulations

It is now well established that changes in the X-ray spectral state of black hole low-mass X-ray binaries are correlated with changes in the radio properties of those systems. Assuming radio power is a proxy for jet power, we can say that the jet is continuously present in the hard state and undetectable (and therefore weaker) in the soft state. Since the different accretion states are also generally assumed to be associated with different disc geometries -- the hard state with a hot, thick flow, and the soft state with a cold, thin disc -- we investigate the possibility that these two phenomena are linked; i.e., that the difference in disc geometry is the cause of the difference in observed jet power. We do this by comparing various measures of jet power in numerical simulations of accretion discs of differing temperatures and thicknesses. We perform these simulations using the general relativistic magnetohydrodynamic code Cosmos++ and a newly added cooling function, which allows us to regulate the disc scale height H/r at different radii. We find no apparent correlation between the disc scale height and jet power whenever we normalize the latter by the mass accretion history of each simulation. We attribute this result to the role that the "corona" plays in confining and accelerating the jet (our corona may also be considered a failed MHD "wind"). The properties of the corona do not vary significantly from one simulation to another, even though the scale heights of the discs vary by up to a factor of four. If this holds true in nature, then it suggests that the correlation between spectral state and jet power must be attributable to some other property, possibly the topology of the magnetic field. Alternatively, it could be that the corona disappears altogether in the soft state, which would be consistent with observations, but has so far not been seen in simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Numerical Simulations of Optically Thick Accretion onto a Black Hole - I. Spherical Case

Modeling the radiation generated by accreting matter is an important step towards realistic simulations of black hole accretion disks, especially at high accretion rates. To this end, we have recently added radiation transport to the existing general relativistic magnetohydrodynamic code, Cosmos++. However, before attempting to model radiative accretion disks, we have tested the new code using a series of shock tube and Bondi (spherical inflow) problems. The four radiative shock tube tests, first presented by Farris et al. (2008), have known analytic solutions, allowing us to calculate errors and convergence rates for our code. The Bondi problem only has an analytic solution when radiative processes are ignored, but is pertinent because it is closer to the physics we ultimately want to study. In our simulations, we include Thomson scattering and thermal bremsstrahlung in the opacity, focusing exclusively on the super-Eddington regime. Unlike accretion onto bodies with solid surfaces, super-Eddington accretion onto black holes does not produce super-Eddington luminosity. In our examples, despite accreting at up to 300 times the Eddington rate, our measured luminosity is always several orders of magnitude below Eddington.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Radiative Models of Sagittarius A* and M87 from Relativistic MHD Simulations

Ongoing millimeter VLBI observations with the Event Horizon Telescope allow unprecedented study of the innermost portion of black hole accretion flows. Interpreting the observations requires relativistic, time-dependent physical modeling. We discuss the comparison of radiative transfer calculations from general relativistic MHD simulations of Sagittarius A* and M87 with current and future mm-VLBI observations. This comparison allows estimates of the viewing geometry and physical conditions of the Sgr A* accretion flow. The viewing geometry for M87 is already constrained from observations of its large-scale jet, but, unlike Sgr A*, there is no consensus for its millimeter emission geometry or electron population. Despite this uncertainty, as long as the emission region is compact, robust predictions for the size of its jet launching region can be made. For both sources, the black hole shadow may be detected with future observations including ALMA and/or the LMT, which would constitute the first direct evidence for a black hole event horizon.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Self-consistent spectra from radiative GRMHD simulations of accretion onto Sgr A*

We present the first spectral energy distributions produced self-consistently by 2.5D general relativistic magneto-hydrodynamical (GRMHD) numerical simulations, where radiative cooling is included in the dynamical calculation. As a case study, we focus on the accretion flow around the supermassive black hole in the Galactic Centre, Sagittarius A* (Sgr A*), which has the best constrained physical parameters. We compare the simulated spectra to the observational data of Sgr A* and explore the parameter space of our model to determine the effect of changing the initial magnetic field configuration, ion to electron temperature ratio T_i/T_e and the target accretion rate. We find the best description of the data for a mass accretion rate of ~ 1e-9 Msun/yr, and rapid spin (0.7 < a_* < 0.9). The submillimeter peak flux seems largely independent of initial conditions, while the higher energies can be very sensitive to the initial magnetic field configuration. Finally, we also discuss flaring features observed in some simulations, that may be due to artifacts of the 2D configuration.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Tilted black hole accretion disc models of Sagittarius A*: time-variable millimetre to near-infrared emission

High-resolution, multi-wavelength, and time-domain observations of the Galactic centre black hole candidate, Sgr A*, allow for a direct test of contemporary accretion theory. To date, all models have assumed alignment between the accretion disc and black hole angular momentum axes, but this is unjustified for geometrically thick accretion flows like that onto Sgr A*. Instead, we calculate images and spectra from a set of simulations of accretion flows misaligned ('tilted') by 15 degrees from the black hole spin axis and compare them with millimetre (mm) to near-infrared (NIR) observations. Non-axisymmetric standing shocks from eccentric fluid orbits dominate the emission, leading to a wide range of possible image morphologies. These effects invalidate previous parameter estimates from model fitting, including estimates of the dimensionless black hole spin, except possibly at low values of spin or tilt. At 1.3mm, the images have crescent morphologies, and the black hole shadow may still be accessible to future mm-VLBI observations. Shock heating leads to high energy electrons (T > 10^12 K), which can naturally produce the observed NIR flux, spectral index, and rapid variability ('flaring'). This NIR emission is uncorrelated with that in the mm, which also agrees with observations. These are the first models to self-consistently explain the time-variable mm to NIR emission of Sgr A*. Predictions of the model include significant structural changes observable with mm-VLBI on both the dynamical (hour) and Lense-Thirring precession (day-year) timescales; and ~30-50 microarcsecond changes in centroid position from extreme gravitational lensing events during NIR flares, detectable with the future VLT instrument GRAVITY. If the observed NIR emission is caused by shock heating in a tilted accretion disc, then the Galactic centre black hole has a positive, non-zero spin parameter (a > 0).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Observational Signatures of Tilted Black Hole Accretion Disks from Simulations

Geometrically thick accretion flows may be present in black hole X-ray binaries observed in the low/hard state and in low-luminosity active galactic nuclei. Unlike in geometrically thin disks, the angular momentum axis in these sources is not expected to align with the black hole spin axis. We compute images from three-dimensional general relativistic magnetohydrodynamic simulations of misaligned (tilted) accretion flows using relativistic radiative transfer, and compare the estimated locations of the radiation edge with expectations from their aligned (untilted) counterparts. The radiation edge in the tilted simulations is independent of black hole spin for a tilt of 15 degrees, in stark contrast to the results for untilted simulations, which agree with the monotonic dependence on spin expected from thin accretion disk theory. Synthetic emission line profiles from the tilted simulations depend strongly on the observer's azimuth, and exhibit unique features such as broad "blue wings." Coupled with precession, the azimuthal variation could generate time fluctuations in observed emission lines, which would be a clear "signature" of a tilted accretion flow. Finally, we evaluate the possibility that the observed low- and high-frequency quasi-periodic oscillations (QPOs) from black hole binaries could be produced by misaligned accretion flows. Although low-frequency QPOs from precessing, tilted disks remains a viable option, we find little evidence for significant power in our light curves in the frequency range of high-frequency QPOs.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Submillimeter Bump in Sgr A* from Relativistic MHD Simulations

Recent high resolution observations of the Galactic center black hole allow for direct comparison with accretion disk simulations. We compare two-temperature synchrotron emission models from three dimensional, general relativistic magnetohydrodynamic simulations to millimeter observations of Sgr A*. Fits to very long baseline interferometry and spectral index measurements disfavor the monochromatic face-on black hole shadow models from our previous work. Inclination angles \le 20 degrees are ruled out to 3 σ. We estimate the inclination and position angles of the black hole, as well as the electron temperature of the accretion flow and the accretion rate, to be i=50+35-15 degrees, ξ=-23+97-22 degrees, T_e=(5.4 +/- 3.0)x10^10 K and Mdot=(5+15-2)x10^-9 M_sun / yr respectively, with 90% confidence. The black hole shadow is unobscured in all best fit models, and may be detected by observations on baselines between Chile and California, Arizona or Mexico at 1.3mm or .87mm either through direct sampling of the visibility amplitude or using closure phase information. Millimeter flaring behavior consistent with the observations is present in all viable models, and is caused by magnetic turbulence in the inner radii of the accretion flow. The variability at optically thin frequencies is strongly correlated with that in the accretion rate. The simulations provide a universal picture of the 1.3mm emission region as a small region near the midplane in the inner radii of the accretion flow, which is roughly isothermal and has ν/ν_c ~ 1-20, where ν_c is the critical frequency for thermal synchrotron emission.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Dynamical Bar-Mode Instability in Differentially Rotating Magnetized Neutron Stars

This paper presents a numerical study over a wide parameter space of the likelihood of the dynamical bar-mode instability in differentially rotating magnetized neutron stars. The innovative aspect of this study is the incorporation of magnetic fields in such a context, which have thus far been neglected in the purely hydrodynamical simulations available in the literature. The investigation uses the Cosmos++ code which allows us to perform three dimensional simulations on a cylindrical grid at high resolution. A sample of Newtonian magneto-hydrodynamical simulations starting from a set of models previously analyzed by other authors without magnetic fields has been performed, providing estimates of the effects of magnetic fields on the dynamical bar-mode deformation of rotating neutron stars. Overall, our results suggest that the effect of magnetic fields are not likely to be very significant in realistic configurations. Only in the most extreme cases are the magnetic fields able to suppress growth of the bar mode.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Effective Inner Radius of Tilted Black Hole Accretion Disks

One of the primary means of determining the spin of an astrophysical black hole is by actually measuring the inner radius of a surrounding accretion disk and using that to infer the spin. By comparing a number of different estimates of the inner radius from simulations of tilted accretion disks with differing black-hole spins, we show that such a procedure can give quite wrong answers. Over the range 0 <= a/M <= 0.9, we find that, for moderately thick disks (H/r ~ 0.2) with modest tilt (15 degrees), the inner radius is nearly independent of spin. This result is likely dependent on tilt, such that for larger tilts, it may even be that the inner radius would increase with increasing spin. In the opposite limit, we confirm through numerical simulations of untilted disks that, in the limit of zero tilt, the inner radius recovers approximately the expected dependence on spin.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Excitation of Trapped Waves in Simulations of Tilted Black Hole Accretion Disks with Magnetorotational Turbulence

We analyze the time dependence of fluid variables in general relativistic, magnetohydrodynamic simulations of accretion flows onto a black hole with dimensionless spin parameter a/M=0.9. We consider both the case where the angular momentum of the accretion material is aligned with the black hole spin axis (an untilted flow) and where it is misaligned by 15 degrees (a tilted flow). In comparison to the untilted simulation, the tilted simulation exhibits a clear excess of inertial variability, that is, variability at frequencies below the local radial epicyclic frequency. We further study the radial structure of this inertial-like power by focusing on a radially extended band at 118 (M/10Msol)^-1 Hz found in each of the three analyzed fluid variables. The three dimensional density structure at this frequency suggests that the power is a composite oscillation whose dominant components are an over dense clump corotating with the background flow, a low order inertial wave, and a low order inertial-acoustic wave. Our results provide preliminary confirmation of earlier suggestions that disk tilt can be an important excitation mechanism for inertial waves.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Millimeter Flares and VLBI Visibilities from Relativistic Simulations of Magnetized Accretion onto the Galactic Center Black Hole

The recent VLBI observation of the Galactic center black hole candidate Sgr A* at 1.3mm shows source structure on event-horizon scales. This detection enables a direct comparison of the emission region with models of the accretion flow onto the black hole. We present the first results from time-dependent radiative transfer of general relativistic MHD simulation data, and compare simulated synchrotron images at black hole spin a=0.9 with the VLBI measurements. After tuning the accretion rate to match the millimeter flux, we find excellent agreement between predicted and observed visibilities, even when viewed face-on (i < 30 degrees). VLBI measurements on 2000-3000km baselines should constrain the inclination. The data constrain the accretion rate to be (1.0-2.3)x10^-9 M_sun / yr with 99% confidence, consistent with but independent of prior estimates derived from spectroscopic and polarimetric measurements. Finally, we compute light curves, which show that magnetic turbulence can directly produce flaring events with .5 hour rise times, 2-3.5 hour durations and 40-50% flux modulation, in agreement with observations of Sgr A* at millimeter wavelengths.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Application of the Cubed-Sphere Grid to Tilted Black-Hole Accretion Disks

In recent work we presented the first results of global general relativistic magnetohydrodynamic (GRMHD) simulations of tilted (or misaligned) accretion disks around rotating black holes. The simulated tilted disks showed dramatic differences from comparable untilted disks, such as asymmetrical accretion onto the hole through opposing "plunging streams" and global precession of the disk powered by a torque provided by the black hole. However, those simulations used a traditional spherical-polar grid that was purposefully underresolved along the pole, which prevented us from assessing the behavior of any jets that may have been associated with the tilted disks. To address this shortcoming we have added a block-structured "cubed-sphere" grid option to the Cosmos++ GRMHD code, which will allow us to simultaneously resolve the disk and polar regions. Here we present our implementation of this grid and the results of a small suite of validation tests intended to demonstrate that the new grid performs as expected. The most important test in this work is a comparison of identical tilted disks, one evolved using our spherical-polar grid and the other with the cubed-sphere grid. We also demonstrate an interesting dependence of the early-time evolution of our disks on their orientation with respect to the grid alignment. This dependence arises from the differing treatment of current sheets within the disks, especially whether they are aligned with symmetry planes of the grid or not.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

General Relativistic Magnetohydrodynamic Simulations of the Hard State as a Magnetically-Dominated Accretion Flow

(Abridged) We present one of the first physically-motivated two-dimensional general relativistic magnetohydrodynamic (GRMHD) numerical simulations of a radiatively-cooled black-hole accretion disk. The fiducial simulation combines a total-energy-conserving formulation with a radiative cooling function, which includes bremsstrahlung, synchrotron, and Compton effects. By comparison with other simulations we show that in optically thin advection-dominated accretion flows, radiative cooling can significantly affect the structure, without necessarily leading to an optically thick, geometrically thin accretion disk. We further compare the results of our radiatively-cooled simulation to the predictions of a previously developed analytic model for such flows. For the very low stress parameter and accretion rate found in our simulated disk, we closely match a state called the "transition" solution between an outer advection-dominated accretion flow and what would be a magnetically-dominated accretion flow (MDAF) in the interior. The qualitative and quantitative agreement between the numerical and analytic models is quite good, with only a few well-understood exceptions. According to the analytic model then, at significantly higher stress or accretion, we would expect a full MDAF to form. The collection of simulations in this work also provide important data for interpreting other numerical results in the literature, as they span the most common treatments of thermodynamics, including simulations evolving: 1) the internal energy only; 2) the internal energy plus an explicit cooling function; 3) the total energy without cooling; and 4) total energy including cooling. We find that the total energy formulation is a necessary prerequisite for proper treatment of radiative cooling in MRI accretion flows.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

3D Relativistic MHD Simulation of a Tilted Accretion Disk Around a Rapidly Rotating Black Hole

We posit that accreting compact objects, including stellar mass black holes and neutron stars as well as supermassive black holes, may undergo extended periods of accretion during which the angular momentum of the disk at large scales is misaligned with that of the compact object. In such a scenario, Lense-Thirring precession caused by the rotating compact object can dramatically affect the disk. In this presentation we describe results from a three-dimensional relativistic magnetohydrodynamic simulation of an MRI turbulent disk accreting onto a tilted rapidly rotating black hole. For this case, the disk does not achieve the commonly described Bardeen-Petterson configuration; rather, it remains nearly planar, undergoing a slow global precession. Accretion from the disk onto the hole occurs predominantly through two opposing plunging streams that start from high latitudes with respect to both the black-hole and disk midplanes. This is a consequence of the non-sphericity of the gravitational spacetime of the black hole.

0 citations0 reviewsNo code linkedOpen access