Showing 1–50 of 88 results for author: Fragile, P C

Tilted Accretion Disks

Authors: P. Chris Fragile, Matthew Liska

Abstract: In this chapter, we review some of the interesting consequences that tilt between the spin axis of the black hole and angular momentum axis of the accretion disk can have on the dynamics, thermodynamics, and observational appearance of accreting systems, from precessing coronae and jets to standing nozzle shocks and quasi-periodic oscillations. We begin the chapter by examining some of the reasons… ▽ More In this chapter, we review some of the interesting consequences that tilt between the spin axis of the black hole and angular momentum axis of the accretion disk can have on the dynamics, thermodynamics, and observational appearance of accreting systems, from precessing coronae and jets to standing nozzle shocks and quasi-periodic oscillations. We begin the chapter by examining some of the reasons tilted disks are interesting as well as present arguments for how ubiquitous they may be. We then review the existing simulation results in the literature, broadly dividing them into sections on thick disks, thin disks, and magnetically arrested disks (MADs). We finish by highlighting some of the phenomenology that is unique to tilted disk simulations and discuss how this may apply to observations. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: 29 pages, 15 figures, chapter for book on GRMHD Simulations

Truncated, Tilted Discs as a Possible Source of Quasi-Periodic Oscillations

Authors: A. Deepika Bollimpalli, P. Chris Fragile, W. Janosz Dewberry, Włodek Kluźniak

Abstract: Many accreting black holes and neutron stars exhibit rapid variability in their X-ray light curves, termed quasi-periodic oscillations (QPOs). The most commonly observed type is the low-frequency ($\lesssim 10$ Hz), type-C QPO, while only a handful of sources exhibit high-frequency QPOs ($\gtrsim 60$ Hz). The leading model for the type-C QPO is Lense-Thirring precession of a hot, geometrically thi… ▽ More Many accreting black holes and neutron stars exhibit rapid variability in their X-ray light curves, termed quasi-periodic oscillations (QPOs). The most commonly observed type is the low-frequency ($\lesssim 10$ Hz), type-C QPO, while only a handful of sources exhibit high-frequency QPOs ($\gtrsim 60$ Hz). The leading model for the type-C QPO is Lense-Thirring precession of a hot, geometrically thick accretion flow that is misaligned with the black hole's spin axis. However, existing versions of this model have not taken into account the effects of a surrounding, geometrically thin disc on the precessing, inner, geometrically thick flow. In Bollimpalli et. al 2023, using a set of GRMHD simulations of tilted, truncated accretion discs, we confirmed that the outer thin disc slows down the precession rate of the precessing torus, which has direct observational implications for type-C QPOs. In this paper, we provide a detailed analysis of those simulations and compare them with an aligned truncated disc simulation. We find that the misalignment of the disc excites additional variability in the inner hot flow, which is absent in the comparable aligned-disc simulations. This suggests that the misalignment may be a crucial requirement for producing QPOs. We attribute this variability to global vertical oscillations of the inner torus at epicyclic frequencies corresponding to the transition radius. This explanation is consistent with current observations of higher frequency QPOs in black hole X-ray binary systems. △ Less

Submitted 22 December, 2023; originally announced December 2023.

Comments: 17 pages, 15 figures, 1 table; accepted for publication in MNRAS

Multi-Frequency General Relativistic Radiation-Magnetohydrodynamic Simulations of Thin Disks

Authors: P. Chris Fragile, Peter Anninos, Nathaniel Roth, Bhupendra Mishra

Abstract: We present a set of six general relativistic, multi-frequency, radiation magnetohydrodynamic simulations of thin accretion disks with different target mass accretion rates around black holes with spins ranging from non-rotating to rapidly spinning. The simulations use the $\mathbf{M}_1$ closure scheme with twelve, independent frequency (or energy) bins ranging logarithmically from… ▽ More We present a set of six general relativistic, multi-frequency, radiation magnetohydrodynamic simulations of thin accretion disks with different target mass accretion rates around black holes with spins ranging from non-rotating to rapidly spinning. The simulations use the $\mathbf{M}_1$ closure scheme with twelve, independent frequency (or energy) bins ranging logarithmically from $5\times 10^{-3}$ to $5\times 10^3$ keV. The multi-frequency capability allows us to generate crude spectra and energy-dependent light curves directly from the simulations without a need for special post-processing. While we generally find roughly thermal spectra with peaks around 1 to 4 keV, our high-spin cases showed harder than expected tails for the soft or thermally dominant state. This leads to radiative efficiencies that are up to five times higher than expected for a Novikov-Thorne disk at the same spin. We attribute these high efficiencies to the high-energy, coronal emission. These coronae mostly occupy the effectively optically thin regions near the inner edges of the disks and also cover or sandwich the inner $\sim 15 GM/c^2$ of the disks. △ Less

Submitted 31 October, 2023; originally announced November 2023.

Comments: 15 pages, 8 figures

Impact of neutron star spin on Poynting-Robertson drag during a Type I X-ray burst

Authors: J. Speicher, P. C. Fragile, D. R. Ballantyne

Abstract: External irradiation of a neutron star (NS) accretion disc induces Poynting-Robertson (PR) drag, removing angular momentum and increasing the mass accretion rate. Recent simulations show PR drag significantly enhancing the mass accretion rate during Type I X-ray bursts, which could explain X-ray spectral features such as an increase in the persistent emission and a soft excess. However, prograde s… ▽ More External irradiation of a neutron star (NS) accretion disc induces Poynting-Robertson (PR) drag, removing angular momentum and increasing the mass accretion rate. Recent simulations show PR drag significantly enhancing the mass accretion rate during Type I X-ray bursts, which could explain X-ray spectral features such as an increase in the persistent emission and a soft excess. However, prograde spin of the NS is expected to weaken PR drag, challenging its importance during bursts. Here, we study the effect of spin on PR drag during X-ray bursts. We run four simulations, with two assuming a non-spinning NS and two using a spin parameter of $a_*=0.2$, corresponding to a rotation frequency of 500 Hz. For each scenario, we simulate the disc evolution subject to an X-ray burst and compare it to the evolution found with no burst. PR drag drains the inner disc region during a burst, moving the inner disc radius outward by $\approx1.6$ km in the $a_*=0$ and by $\approx2.2$ km in the $a_*=0.2$ simulation. The burst enhances the mass accretion rate across the innermost stable circular orbit $\approx7.9$ times when the NS is not spinning and $\approx11.2$ times when it is spinning. The explanation for this seemingly contradictory result is that the disc is closer to the NS when $a_*=0.2$, and the resulting stronger irradiating flux offsets the weakening effect of spin on the PR drag. Hence, PR drag remains a viable explanation for the increased persistent emission and soft excess observed during X-ray bursts in spinning NS systems. △ Less

Submitted 31 August, 2023; originally announced September 2023.

Comments: 9 pages, 8 figures, accepted for publication in MNRAS

The Luminous, Hard State Can't Be MAD

Authors: P. Chris Fragile, Koushik Chatterjee, Adam Ingram, Matthew Middleton

Abstract: We present a straightforward argument for why the luminous, hard state of black hole X-ray binaries (BHXRBs) cannot always be associated with a magnetically arrested accretion disc (MAD). It relies on three core premises: 1) that the type-C quasi-periodic oscillation (QPO) is best explained by Lense-Thirring (LT) precession of a tilted, inner, hot flow; 2) that observed optical and infrared (IR) Q… ▽ More We present a straightforward argument for why the luminous, hard state of black hole X-ray binaries (BHXRBs) cannot always be associated with a magnetically arrested accretion disc (MAD). It relies on three core premises: 1) that the type-C quasi-periodic oscillation (QPO) is best explained by Lense-Thirring (LT) precession of a tilted, inner, hot flow; 2) that observed optical and infrared (IR) QPOs with the same or lower frequency as the type-C QPO suggest the jet, too, must precess in these systems; and 3) that numerical simulations of MADs show that their strong magnetic fields promote alignment of the disc with the black hole and, thereby, suppress LT precession. If all three premises hold true, then, at least whenever the optical and IR QPOs are observed alongside the type-C QPO, these systems cannot be in the MAD state. Extending the argument further, if the type-C QPO is always associated with LT precession, then it would rule out MADs anytime this timing feature is seen, which covers nearly all BHXRBs when they are in the luminous, hard and hard-intermediate states. △ Less

Submitted 17 July, 2023; originally announced July 2023.

Comments: 5 pages, 2 figures, accepted to MNRAS

Effect of geometrically thin discs on precessing, thick flows: Relevance to type-C QPOs

Authors: D. A. Bollimpalli, P. C. Fragile, W. Kluźniak

Abstract: Type-C quasi-periodic oscillations (QPOs) are the low-frequency QPOs most commonly observed during the hard spectral state of X-ray binary systems. The leading model for these QPOs is the Lense-Thirring precession of a hot, geometrically thick accretion flow that is misaligned with respect to the black hole spin axis. However, none of the work done to date has accounted for the effects of a surrou… ▽ More Type-C quasi-periodic oscillations (QPOs) are the low-frequency QPOs most commonly observed during the hard spectral state of X-ray binary systems. The leading model for these QPOs is the Lense-Thirring precession of a hot, geometrically thick accretion flow that is misaligned with respect to the black hole spin axis. However, none of the work done to date has accounted for the effects of a surrounding, geometrically thin disc on this precession, as would be the case in the truncated disc picture of the hard state. To address this, we perform a set of GRMHD simulations of truncated discs misaligned with the spin axes of their central black holes. Our results confirm that the inner-hot flow still undergoes precession, though at a rate that is only 5 percent of what is predicted for an isolated, precessing torus. We find that the exchange of angular momentum between the outer, thin and the inner, thick disc causes this slow-down in the precession rate and discuss its relevance to type-C QPOs. △ Less

Submitted 6 October, 2022; originally announced October 2022.

Comments: Submitted to MNRAS Letters

The Role of Strong Magnetic Fields in Stabilizing Highly Luminous, Thin Disks

Authors: Bhupendra Mishra, P. Chris Fragile, Jessica Anderson, Aidan Blankenship, Hui Li, Krzysztof Nalewajko

Abstract: We present a set of three-dimensional, global, general relativistic radiation magnetohydrodynamic simulations of thin, radiation-pressure-dominated accretion disks surrounding a non-rotating, stellar-mass black hole. The simulations are initialized using the Shakura-Sunyaev model with a mass accretion rate of $\dot{M} = 3 L_\mathrm{Edd}/c^2$ (corresponding to $L=0.17 L_\mathrm{Edd}$). Our previous… ▽ More We present a set of three-dimensional, global, general relativistic radiation magnetohydrodynamic simulations of thin, radiation-pressure-dominated accretion disks surrounding a non-rotating, stellar-mass black hole. The simulations are initialized using the Shakura-Sunyaev model with a mass accretion rate of $\dot{M} = 3 L_\mathrm{Edd}/c^2$ (corresponding to $L=0.17 L_\mathrm{Edd}$). Our previous work demonstrated that such disks are thermally unstable when accretion is driven by an $α$-viscosity. In the present work, we test the hypothesis that strong magnetic fields can both drive accretion through the magneto-rotational instability and restore stability to such disks. We test four initial magnetic field configurations: 1) a zero-net-flux case with a single, radially extended set of magnetic field loops (dipole); 2) a zero-net-flux case with two radially extended sets of magnetic field loops of opposite polarity stacked vertically (quadrupole); 3) a zero-net-flux case with multiple radially concentric rings of alternating polarity (multi-loop); and 4) a net-flux, vertical magnetic field configuration (vertical). In all cases, the fields are initially weak, with the gas-to-magnetic pressure ratio $\gtrsim 100$. Based on the results of these simulations, we find that the dipole and multi-loop configurations remain thermally unstable like their $α$-viscosity counterpart, in our case collapsing vertically on the local thermal timescale and never fully recovering. The vertical case, on the other hand, stabilizes and remains so for the duration of our tests (many thermal timescales). The quadrupole case is intermediate, showing signs of both stability and instability. The key stabilizing criteria is, $P_\mathrm{mag} \gtrsim 0.5P_\mathrm{tot}$ with strong toroidal fields near the disk midplane. We also report a comparison of our models to the standard Shakura-Sunyaev disk. △ Less

Submitted 20 September, 2022; v1 submitted 7 September, 2022; originally announced September 2022.

Comments: 19 pages, 10 figures

Observations and Simulations of Radio Emission and Magnetic Fields in Minkowski's Object

Authors: C. Nolting, M. Lacy, S. Croft, P. C. Fragile, S. T. Linden, K. Nyland, P. Patil

Abstract: We combine new data from the Karl G. Jansky Very Large Array with previous radio observations to create a more complete picture of the ongoing interactions between the radio jet from galaxy NGC 541 and the star-forming system known as Minkowski's Object (MO). We then compare those observations with synthetic radio data generated from a new set of magnetohydrodynamic simulations of a jet-cloud inte… ▽ More We combine new data from the Karl G. Jansky Very Large Array with previous radio observations to create a more complete picture of the ongoing interactions between the radio jet from galaxy NGC 541 and the star-forming system known as Minkowski's Object (MO). We then compare those observations with synthetic radio data generated from a new set of magnetohydrodynamic simulations of a jet-cloud interaction specifically tailored to the parameters of MO. The combination of radio intensity, polarization, and spectral index measurements all convincingly support the interaction scenario and provide additional constraints on the local dynamical state of the intracluster medium and the time since the jet-cloud interaction first began. In particular, we show that only a simulation with a bent radio jet can reproduce the observations. △ Less

Submitted 23 August, 2022; v1 submitted 9 June, 2022; originally announced June 2022.

Comments: 17 pages, 6 figures, links to simulation movies in figure captions, accepted for publication in ApJ

Evolution of accretion disc reflection spectra due to a Type I X-ray burst

Authors: J. Speicher, D. R. Ballantyne, P. C. Fragile

Abstract: Irradiation of the accretion disc causes reflection signatures in the observed X-ray spectrum, encoding important information about the disc structure and density. A Type I X-ray burst will strongly irradiate the accretion disc and alter its properties. Previous numerical simulations predicted the evolution of the accretion disc due to an X-ray burst. Here, we process time-averaged simulation data… ▽ More Irradiation of the accretion disc causes reflection signatures in the observed X-ray spectrum, encoding important information about the disc structure and density. A Type I X-ray burst will strongly irradiate the accretion disc and alter its properties. Previous numerical simulations predicted the evolution of the accretion disc due to an X-ray burst. Here, we process time-averaged simulation data of six time intervals to track changes in the reflection spectrum from the burst onset to just past its peak. We divide the reflecting region of the disc within $r\lesssim50$ km into 6-7 radial zones for every time interval and compute the reflection spectra for each zone. We integrate these reflection spectra to obtain a total reflection spectrum per time interval. The burst ionizes and heats the disc, which gradually weakens all emission lines. Compton scattering and bremsstrahlung rates increase in the disc during the burst rise, and the soft excess at $<$3 keV rises from $\approx4$% to $\approx38$% of the total emission at the burst peak. A soft excess is expected to be ubiquitous in the reflection spectra of X-ray bursts. Structural disc changes such as inflation because of heating or drainage of the inner disc due to Poynting-Robertson drag affect the strength of the soft excess. Further studies on the dependence of the reflection spectrum characteristics to changes in the accretion disc during an X-ray burst may lead to probes of the disc geometry. △ Less

Submitted 22 October, 2021; originally announced October 2021.

Comments: 9 pages, 10 figures, accepted for publication in MNRAS

Magneto-rotational instability in magnetically polarized discs

Authors: Oscar M. Pimentel, P. Chris Fragile, F. D. Lora-Clavijo, Bridget Ierace, Deepika Bollimpalli

Abstract: The magneto-rotational instability (MRI) is the most likely mechanism for transportation of angular momentum and dissipation of energy within hot, ionized accretion discs. This instability is produced through the interactions of a differentially rotating plasma with an embedded magnetic field. Like all substances in nature, the plasma in an accretion disc has the potential to become magnetically p… ▽ More The magneto-rotational instability (MRI) is the most likely mechanism for transportation of angular momentum and dissipation of energy within hot, ionized accretion discs. This instability is produced through the interactions of a differentially rotating plasma with an embedded magnetic field. Like all substances in nature, the plasma in an accretion disc has the potential to become magnetically polarized when it interacts with the magnetic field. In this paper, we study the effect of this magnetic susceptibility, parameterized by $χ_m$, on the MRI, specifically within the context of black hole accretion. We find from a linear analysis within the Newtonian limit that the minimum wavelength of the first unstable mode and the wavelength of the fastest growing mode are shorter in paramagnetic ($χ_m>0$) than in diamagnetic ($χ_m<0$) discs, all other parameters being equal. Furthermore, the magnetization parameter (ratio of gas to magnetic pressure) in the saturated state should be smaller when the magnetic susceptibility is positive than when it is negative. We confirm this latter prediction through a set of numerical simulations of magnetically polarized black hole accretion discs. We additionally find that the vertically integrated stress and mass accretion rate are somewhat larger when the disc is paramagnetic than when it is diamagnetic. If astrophysical discs are able to become magnetically polarized to any significant degree, then our results would be relevant to properly interpreting observations. △ Less

Submitted 24 May, 2021; originally announced May 2021.

Comments: 12 pages, 8 figures, accepted for publication in MNRAS

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

Authors: Peter Anninos, P. Chris Fragile

Abstract: 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… ▽ More 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. △ Less

Submitted 23 July, 2020; originally announced July 2020.

Comments: 22 pages, 11 figures, submitted to Astrophysical Journal

Looking for the underlying cause of black hole X-ray variability in GRMHD simulations

Authors: D. A. Bollimpalli, R. Mahmoud, C. Done, P. C. Fragile, W. Kluźniak, R. Narayan, C. J. White

Abstract: Long-term observations have shown that black hole X-ray binaries exhibit strong, aperiodic variability on time-scales of a few milliseconds to seconds. The observed light curves display various characteristic features like a log-normal distribution of flux and a linear rms-flux relation, which indicate that the underlying variability process is stochastic in nature. It is also thought to be intrin… ▽ More Long-term observations have shown that black hole X-ray binaries exhibit strong, aperiodic variability on time-scales of a few milliseconds to seconds. The observed light curves display various characteristic features like a log-normal distribution of flux and a linear rms-flux relation, which indicate that the underlying variability process is stochastic in nature. It is also thought to be intrinsic to accretion. This variability has been modelled as inward propagating fluctuations of mass accretion rate, although the physical process driving the fluctuations remains puzzling. In this work, we analyse five exceptionally long duration general relativistic magnetohydrodynamic (GRMHD) simulations of optically thin, geometrically thick, black hole accretion flows to look for hints of propagating fluctuations in the simulation data. We find that the accretion profiles from these simulations do show evidence for inward propagating fluctuations below the viscous frequency by featuring strong radial coherence and positive time lags when comparing smaller to larger radii, although these time lags are generally shorter than the viscous time-scale and frequency independent. Our simulations also support the notion that the fluctuations in $\dot{M}$ build up in a multiplicative manner, as the simulations exhibit linear rms-mass flux relations, as well as log-normal distributions of their mass fluxes. When combining the mass fluxes from the simulations with an assumed emissivity profile, we additionally find broad agreement with observed power spectra and time lags, including a recovery of the frequency dependency of the time lags. △ Less

Submitted 27 July, 2020; v1 submitted 1 June, 2020; originally announced June 2020.

Comments: 22 pages, 20 figures, 4 tables

Journal ref: MNRAS 496 (2020), 3808-3828

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

Authors: P. Chris Fragile, David R. Ballantyne, Aidan Blankenship

Abstract: 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 outflow… ▽ More 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. △ Less

Submitted 3 January, 2020; originally announced January 2020.

Comments: 26 pages, 10 figures

Neutron Star QPOs from Oscillating, Precessing Hot, Thick Flow

Authors: P. Chris Fragile

Abstract: 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… ▽ More 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. △ Less

Submitted 21 November, 2019; originally announced November 2019.

Comments: 6 pages, 3 figures, accepted for publication in MNRAS

Nuclear Ignition of White Dwarf Stars by Relativistic Encounters with Rotating Intermediate Mass Black Holes

Authors: Peter Anninos, Robert D. Hoffman, Manvir Grewal, Michael J. Lavell, P. Chris Fragile

Abstract: We present results from general relativistic calculations of nuclear ignition in white dwarf stars triggered by near encounters with rotating intermediate mass black holes with different spin and alignment parameters. These encounters create thermonuclear environments characteristic of Type Ia supernovae capable of producing both calcium and iron group elements in arbitrary ratios, depending prima… ▽ More We present results from general relativistic calculations of nuclear ignition in white dwarf stars triggered by near encounters with rotating intermediate mass black holes with different spin and alignment parameters. These encounters create thermonuclear environments characteristic of Type Ia supernovae capable of producing both calcium and iron group elements in arbitrary ratios, depending primarily on the proximity of the interaction which acts as a strong moderator of nucleosynthesis. We explore the effects of black hole spin and spin-orbital alignment on burn product synthesis to determine whether they might also be capable of moderating reactive flows. When normalized to equivalent impact penetration, accounting for frame dragging corrections, the influence of spin is weak, no more than 25% as measured by nuclear energy release and mass of burn products, even for near maximally rotating black holes. Stars on prograde trajectories approach closer to the black hole and produce significantly more unbound debris and iron group elements than is possible by encounters with nonrotating black holes or by retrograde orbits, at more than 50% mass conversion efficiency. The debris contains several radioisotopes, most notably Ni56, made in amounts that produce sub-luminous (but still observable) light curves compared to branch-normal SNe Ia. △ Less

Submitted 2 October, 2019; originally announced October 2019.

Comments: 20 pages, 7 figures, submitted to ApJ

Journal ref: The Astrophysical Journal, Volume 885, Issue 2, article id. 136, 14 pp. (2019)

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

Authors: Bhupendra Mishra, Wlodek Kluźniak, P. Chris Fragile

Abstract: 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… ▽ More 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. △ Less

Submitted 13 September, 2019; originally announced September 2019.

Comments: submitted MNRAS, comments are welcome

Chandra Spectral and Timing Analysis of Sgr A*'s Brightest X-ray Flares

Authors: Daryl Haggard, Melania Nynka, Brayden Mon, Noelia de la Cruz Hernandez, Michael Nowak, Craig Heinke, Joseph Neilsen, Jason Dexter, P. Chris Fragile, Fred Baganoff, Geoffrey C. Bower, Lia R. Corrales, Francesco Coti Zelati, Nathalie Degenaar, Sera Markoff, Mark R. Morris, Gabriele Ponti, Nanda Rea, Joern Wilms, Farhad Yusef-Zadeh

Abstract: We analyze the two brightest Chandra X-ray flares detected from Sagittarius A*, with peak luminosities more than 600 x and 245 x greater than the quiescent X-ray emission. The brightest flare has a distinctive double-peaked morphology --- it lasts 5.7 ksec ($\sim 2$ hours), with a rapid rise time of 1500 sec and a decay time of 2500 sec. The second flare lasts 3.4 ksec, with rise and decay times o… ▽ More We analyze the two brightest Chandra X-ray flares detected from Sagittarius A*, with peak luminosities more than 600 x and 245 x greater than the quiescent X-ray emission. The brightest flare has a distinctive double-peaked morphology --- it lasts 5.7 ksec ($\sim 2$ hours), with a rapid rise time of 1500 sec and a decay time of 2500 sec. The second flare lasts 3.4 ksec, with rise and decay times of 1700 sec and 1400 sec. These luminous flares are significantly harder than quiescence: the first has a power law spectral index $Γ= 2.06\pm 0.14$ and the second has $Γ= 2.03\pm 0.27$, compared to $Γ= 3.0\pm0.2$ for the quiescent accretion flow. These spectral indices (as well as the flare hardness ratios) are consistent with previously-detected Sgr A* flares, suggesting that bright and faint flares arise from similar physical processes. Leveraging the brightest flare's long duration and high signal-to-noise, we search for intraflare variability and detect excess X-ray power at a frequency of $ν\approx 3$ mHz, but show that it is an instrumental artifact and not of astrophysical origin. We find no other evidence (at the 95% confidence level) for periodic or quasi-periodic variability in either flares' time series. We also search for non-periodic excess power but do not find compelling evidence in the power spectrum. Bright flares like these remain our most promising avenue for identifying Sgr A*'s short timescale variability in the X-ray, which may probe the characteristic size scale for the X-ray emission region. △ Less

Submitted 4 December, 2019; v1 submitted 5 August, 2019; originally announced August 2019.

Comments: Updated to match published version; 19 pages, 7 figures, 3 tables

Journal ref: The Astrophysical Journal, Volume 886, Issue 2, article id. 96, 14 pp. (2019)

The Lense-Thirring timing-accretion plane for ULXs

Authors: M. Middleton, P. C. Fragile, A. Ingram, T. P. Roberts

Abstract: Identifying the compact object in ultraluminous X-ray sources (ULXs) has to-date required detection of pulsations or a cyclotron resonance scattering feature (CRSF), indicating a magnetised neutron star. However, pulsations are observed to be transient and it is plausible that accretion onto the neutron star may have suppressed the surface magnetic field such that pulsations and CRSFs will be enti… ▽ More Identifying the compact object in ultraluminous X-ray sources (ULXs) has to-date required detection of pulsations or a cyclotron resonance scattering feature (CRSF), indicating a magnetised neutron star. However, pulsations are observed to be transient and it is plausible that accretion onto the neutron star may have suppressed the surface magnetic field such that pulsations and CRSFs will be entirely absent. We may therefore lack direct means to identify neutron star systems whilst we presently lack an effective means by which to identify black hole ULXs. Here we present a possible method for separating the ULX population by assuming the X-ray, mHz quasi-periodic oscillations (QPOs) and day timescale periods/QPOs are associated with Lense-Thirring precession of the inflow and outflowing wind respectively. The precession timescales combined with the temperature of the soft X-ray component produce planes where the accretor mass enters as a free parameter. Depending on the properties of the wind, use of these planes may be robust to a range in the angular momentum (spin) and, for high accretion rates, essentially independent of the neutron star's surface dipole field strength. Our model also predicts the mHz QPO frequency and magnitude of the phase-lag imprinted due to propagation through the optically thick wind; in the case of NGC 5408 X-1 we subsequently infer a black hole mass and high spin. Finally, we note that observing secular QPO evolution over sufficient baselines may indicate a neutron star, as the precession responds to spin-up which is not readily observable for black hole primaries. △ Less

Submitted 9 September, 2019; v1 submitted 7 May, 2019; originally announced May 2019.

Comments: 17 pages, accepted to MNRAS (updated, pre-epoof)

The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project

Authors: Oliver Porth, Koushik Chatterjee, Ramesh Narayan, Charles F. Gammie, Yosuke Mizuno, Peter Anninos, John G. Baker, Matteo Bugli, Chi-kwan Chan, Jordy Davelaar, Luca Del Zanna, Zachariah B. Etienne, P. Chris Fragile, Bernard J. Kelly, Matthew Liska, Sera Markoff, Jonathan C. McKinney, Bhupendra Mishra, Scott C. Noble, Héctor Olivares, Ben Prather, Luciano Rezzolla, Benjamin R. Ryan, James M. Stone, Niccolò Tomei , et al. (3 additional authors not shown)

Abstract: Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic… ▽ More Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL. Agreement between the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems. △ Less

Submitted 5 August, 2019; v1 submitted 9 April, 2019; originally announced April 2019.

Comments: Accepted version for publication in ApJS (May 28th, 2019)

A Remarkably Loud Quasi-Periodicity after a Star is Disrupted by a Massive Black Hole

Authors: Dheeraj R. Pasham, Ronald A. Remillard, P. Chris Fragile, Alessia Franchini, Nicholas C. Stone, Giuseppe Lodato, Jeroen Homan, Deepto Chakrabarty, Frederick K. Baganoff, James F. Steiner, Eric R. Coughlin, Nishanth R. Pasham

Abstract: The immense tidal forces of massive black holes can rip apart stars that come too close to them. As the resulting stellar debris spirals inwards, it heats up and emits x-rays when near the black hole. Here, we report the discovery of an exceptionally stable 131-second x-ray quasi-periodicity from a black hole after it disrupted a star. Using a black hole mass indicated from host galaxy scaling rel… ▽ More The immense tidal forces of massive black holes can rip apart stars that come too close to them. As the resulting stellar debris spirals inwards, it heats up and emits x-rays when near the black hole. Here, we report the discovery of an exceptionally stable 131-second x-ray quasi-periodicity from a black hole after it disrupted a star. Using a black hole mass indicated from host galaxy scaling relations implies that, (1) this periodicity originates from very close to the black hole's event horizon, and (2) the black hole is rapidly spinning. Our findings suggest that other disruption events with similar highly sensitive observations likely also exhibit quasi-periodicities that encode information about the fundamental properties of their black holes. △ Less

Submitted 25 October, 2018; originally announced October 2018.

Comments: Re-submitted after minor comments from referees

Journal ref: Science 01 Feb 2019

Simulating the Collapse of a Thick Accretion Disk due to a Type I X-ray Burst from a Neutron Star

Authors: P. Chris Fragile, David R. Ballantyne, Thomas J. Maccarone, Jason W. L. Witry

Abstract: We use two-dimensional, general relativistic, viscous, radiation hydrodynamic simulations to study the impact of a Type I X-ray burst on a hot and geometrically thick accretion disk surrounding an unmagnetized, non-rotating neutron star. The disk is initially consistent with a system in its low/hard spectral state, and is subject to a burst which rises to a peak luminosity of $10^{38}$ erg s… ▽ More We use two-dimensional, general relativistic, viscous, radiation hydrodynamic simulations to study the impact of a Type I X-ray burst on a hot and geometrically thick accretion disk surrounding an unmagnetized, non-rotating neutron star. The disk is initially consistent with a system in its low/hard spectral state, and is subject to a burst which rises to a peak luminosity of $10^{38}$ erg s$^{-1}$ in $2.05$ s. At the peak of the burst, the temperature of the disk has dropped by more than three orders of magnitude and its scale height has gone down by more than one order of magnitude. The simulations show that these effects predominantly happen due to Compton cooling of the hot plasma, and clearly illustrate the potential cooling effects of bursts on accretion disk coronae. In addition, we demonstrate the presence of Poynting-Robertson drag, though it only enhances the mass accretion rate onto the neutron star by a factor of $\sim 3$-$4$ compared to a simulation with no burst. Simulations such as these are important for building a general understanding of the response of an accretion disk to an intense X-ray impulse, which, in turn, will be crucial for deciphering burst spectra. Detailed analysis of such spectra offers the potential to measure neutron star radii, and hence constrain the neutron star equation of state, but only if the contributions coming from the impacted disk and its associated corona can be understood. △ Less

Submitted 18 October, 2018; originally announced October 2018.

Comments: 7 pages, 3 figures, to appear in ApJ Letters

Breathing Oscillations in a Global Simulation of a Thin Accretion Disk

Authors: Bhupendra Mishra, Wlodek Kluzniak, P. Chris Fragile

Abstract: We study the oscillations of an axisymmetric, viscous, radiative, general relativistic hydrodynamical simulation of a geometrically thin disk around a non-rotating, $6.62\,M_\odot$ black hole. The numerical setup is initialized with a Novikov-Thorne, gas-pressure-dominated accretion disk, with an initial mass accretion rate of $\dot{m} = 0.01\,L_\mathrm{Edd}/c^2$ (where $L_\mathrm{Edd}$ is the Edd… ▽ More We study the oscillations of an axisymmetric, viscous, radiative, general relativistic hydrodynamical simulation of a geometrically thin disk around a non-rotating, $6.62\,M_\odot$ black hole. The numerical setup is initialized with a Novikov-Thorne, gas-pressure-dominated accretion disk, with an initial mass accretion rate of $\dot{m} = 0.01\,L_\mathrm{Edd}/c^2$ (where $L_\mathrm{Edd}$ is the Eddington luminosity and $c$ is the speed of light). Viscosity is treated with the $α$-prescription. The simulation was evolved for about $1000$ Keplerian orbital periods at three Schwarzschild radii (ISCO radius). Power density spectra of the radial and vertical fluid velocity components, the total (gas $+$ radiation) midplane pressure, and the vertical component of radiative flux from the photosphere, all reveal strong power at the local breathing oscillation frequency. The first, second and third harmonics of the breathing oscillation are also clearly seen in the data. We quantify the properties of these oscillations by extracting eigenfunctions of the radial and vertical velocity components and total pressure. This confirms that these oscillations are associated with breathing motion. △ Less

Submitted 12 October, 2018; originally announced October 2018.

Comments: 10 pages, 6figures, submitted to MNRAS

Relativistic Tidal Disruption and Nuclear Ignition of White Dwarf Stars by Intermediate Mass Black Holes

Authors: Peter Anninos, P. Chris Fragile, Samuel S. Olivier, Robert Hoffman, Bhupendra Mishra, Karen Camarda

Abstract: We present results from general relativistic calculations of the tidal disruption of white dwarf stars from near encounters with intermediate mass black holes. We follow the evolution of 0.2 and $0.6 M_\odot$ stars on parabolic trajectories that approach $10^3$ - $10^4 M_\odot$ black holes as close as a few Schwarzschild radii at periapsis, paying particular attention to the effect tidal disruptio… ▽ More We present results from general relativistic calculations of the tidal disruption of white dwarf stars from near encounters with intermediate mass black holes. We follow the evolution of 0.2 and $0.6 M_\odot$ stars on parabolic trajectories that approach $10^3$ - $10^4 M_\odot$ black holes as close as a few Schwarzschild radii at periapsis, paying particular attention to the effect tidal disruption has on thermonuclear reactions and the synthesis of intermediate to heavy ion elements. These encounters create diverse thermonuclear environments characteristic of Type I supernovae and capable of producing both intermediate and heavy mass elements in arbitrary ratios, depending on the strength (or proximity) of the interaction. Nuclear ignition is triggered in all of our calculations, even at weak tidal strengths $β\sim 2.6$ and large periapsis radius $R_P \sim 28$ Schwarzschild radii. A strong inverse correlation exists between the mass ratio of calcium to iron group elements and tidal strength, with $β\lesssim 5$ producing predominately calcium-rich debris. At these moderate to weak interactions, nucleosynthesis is not especially efficient, limiting the total mass and outflows of calcium group elements to $< 15$\% of available nuclear fuel. Iron group elements however continue to be produced in greater quantity and ratio with increasing tidal strength, peaking at $\sim 60$\% mass conversion efficiency in our closest encounter cases. These events generate short bursts of gravitational waves with characteristic frequencies 0.1-0.7 Hz and strain amplitudes $0.5\times10^{-22}$ - $3.5\times10^{-22}$ at 10 Mpc source distance. △ Less

Submitted 16 August, 2018; originally announced August 2018.

Comments: 24 pages, 17 figures, accepted for publication in The Astrophysical Journal

Relativistic Viscous Radiation Hydrodynamic Simulations of Geometrically Thin Disks: I. Thermal and Other Instabilities

Authors: P. Chris Fragile, Sarina M. Etheridge, Peter Anninos, Bhupendra Mishra, Wlodek Kluzniak

Abstract: We present results from two-dimensional, general relativistic, viscous, radiation hydrodynamic numerical simulations of Shakura-Sunyaev thin disks accreting onto stellar mass Schwarzschild black holes. We consider cases on both the gas- and 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$ t… ▽ More We present results from two-dimensional, general relativistic, viscous, radiation hydrodynamic numerical simulations of Shakura-Sunyaev thin disks accreting onto stellar mass Schwarzschild black holes. We consider cases on both the gas- and 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$. The simulations directly test the stability of this standard disk model on the different branches. We find clear evidence of thermal instability for all radiation-pressure-dominated disks, resulting universally in the vertical collapse of the disks, which in some cases then settle onto the stable, gas-pressure-dominated branch. Although these results are consistent with decades-old theoretical predictions, they appear to be in conflict with available observational data from black hole X-ray binaries. We also find evidence for a radiation-pressure-driven instability that breaks the unstable disks up into alternating rings of high and low surface density on a timescale comparable to the thermal collapse. Since radiation is included self-consistently in the simulations, we are able to calculate lightcurves and power density spectra (PDS). For the most part, we measure radiative efficiencies (ratio of luminosity to mass accretion rate) close to 6%, as expected for a non-rotating black hole. The PDS appear as broken power laws, with a break typically around 100 Hz. There is no evidence of significant excess power at any frequencies, i.e. no quasi-periodic oscillations are observed. △ Less

Submitted 16 March, 2018; originally announced March 2018.

Comments: 16 pages, 15 figures, accepted for publication in ApJ

Lense-Thirring precession in ULXs as a possible means to constrain the neutron star equation-of-state

Authors: M. J. Middleton, P. C. Fragile, M. Bachetti, M. Brightman, Y-F. Jiang, W. C. G. Ho, T. P. Roberts, A. R. Ingram, T. Dauser, C. Pinto, D. J. Walton, F. Fuerst, A. C. Fabian, N. Gehrels

Abstract: The presence of neutron stars in at least three ultraluminous X-ray sources is now firmly established and offers an unambiguous view of super-critical accretion. All three systems show long-timescale periods (60-80 days) in the X-rays and/or optical, two of which are known to be super-orbital in nature. Should the flow be classically super critical, i.e. the Eddington limit is reached locally in t… ▽ More The presence of neutron stars in at least three ultraluminous X-ray sources is now firmly established and offers an unambiguous view of super-critical accretion. All three systems show long-timescale periods (60-80 days) in the X-rays and/or optical, two of which are known to be super-orbital in nature. Should the flow be classically super critical, i.e. the Eddington limit is reached locally in the disc (implying surface dipole fields that are sub-magnetar in strength), then the large scale-height flow can precess through the Lense-Thirring effect which could provide an explanation for the observed super-orbital periods. By connecting the details of the Lense-Thirring effect with the observed pulsar spin period, we are able to infer the moment-of-inertia and therefore equation-of-state of the neutron star without relying on the inclination of, or distance to the system. We apply our technique to the case of NGC 7793 P13 and demonstrate that stronger magnetic fields imply stiffer equations of state. We discuss the caveats and uncertainties, many of which can be addressed through forthcoming radiative magnetohydrodynamic (RMHD) simulations and their connection to observation. △ Less

Submitted 20 November, 2017; originally announced November 2017.

Comments: 14 pages, 9 figures, accepted to MNRAS

Simultaneous Monitoring of X-ray and Radio Variability in Sagittarius A*

Authors: Daniel M. Capellupo, Daryl Haggard, Nicolas Choux, Fred Baganoff, Geoffrey C. Bower, Bill Cotton, Nathalie Degenaar, Jason Dexter, Heino Falcke, P. Chris Fragile, Craig O. Heinke, Casey J. Law, Sera Markoff, Joey Neilsen, Gabriele Ponti, Nanda Rea, Farhad Yusef-Zadeh

Abstract: Monitoring of Sagittarius A* from X-ray to radio wavelengths has revealed structured variability --- including X-ray flares --- but it is challenging to establish correlations between them. Most studies have focused on variability in the X-ray and infrared, where variations are often simultaneous, and because long time series at sub-millimeter and radio wavelengths are limited. Previous work on su… ▽ More Monitoring of Sagittarius A* from X-ray to radio wavelengths has revealed structured variability --- including X-ray flares --- but it is challenging to establish correlations between them. Most studies have focused on variability in the X-ray and infrared, where variations are often simultaneous, and because long time series at sub-millimeter and radio wavelengths are limited. Previous work on sub-mm and radio variability hints at a lag between X-ray flares and their candidate sub-millimeter or radio counterparts, with the long wavelength data lagging the X-ray. However, there is only one published time lag between an X-ray flare and a possible radio counterpart. Here we report 9 contemporaneous X-ray and radio observations of Sgr A*. We detect significant radio variability peaking $\gtrsim$176 minutes after the brightest X-ray flare ever detected from Sgr A*. We also report other potentially associated X-ray and radio variability, with the radio peaks appearing $\lesssim$80 minutes after these weaker X-ray flares. Taken at face value, these results suggest that stronger X-ray flares lead to longer time lags in the radio. However, we also test the possibility that the variability at X-ray and radio wavelengths is not temporally correlated. We cross-correlate data from mismatched X-ray and radio epochs and obtain comparable correlations to the matched data. Hence, we find no overall statistical evidence that X-ray flares and radio variability are correlated, underscoring a need for more simultaneous, long duration X-ray--radio monitoring of Sgr A*. △ Less

Submitted 6 July, 2017; originally announced July 2017.

Comments: Accepted for publication in ApJ

CosmosDG: An hp-Adaptive Discontinuous Galerkin Code for Hyper-Resolved Relativistic MHD

Authors: Peter Anninos, Colton Bryant, P. Chris Fragile, A. Miguel Holgado, Cheuk Lau, Daniel Nemergut

Abstract: We have extended Cosmos++, a multi-dimensional unstructured adaptive mesh code for solving the covariant Newtonian and general relativistic radiation magnetohydrodynamic (MHD) equations, to accommodate both discrete finite volume and arbitrarily high order finite element structures. The new finite element implementation, called CosmosDG, is based on a discontinuous Galerkin (DG) formulation, using… ▽ More We have extended Cosmos++, a multi-dimensional unstructured adaptive mesh code for solving the covariant Newtonian and general relativistic radiation magnetohydrodynamic (MHD) equations, to accommodate both discrete finite volume and arbitrarily high order finite element structures. The new finite element implementation, called CosmosDG, is based on a discontinuous Galerkin (DG) formulation, using both entropy-based artificial viscosity and slope limiting procedures for regularization of shocks. High order multi-stage forward Euler and strong stability preserving Runge-Kutta time integration options complement high order spatial discretization. We have also added flexibility in the code infrastructure allowing for both adaptive mesh and adaptive basis order refinement to be performed separately or simultaneously in a local (cell-by-cell) manner. We discuss in this report the DG formulation and present tests demonstrating the robustness, accuracy, and convergence of our numerical methods applied to special and general relativistic MHD, though we note an equivalent capability currently also exists in CosmosDG for Newtonian systems. △ Less

Submitted 29 June, 2017; originally announced June 2017.

Comments: 21 pages, 14 figures, submitted to The Astrophysical Journal

On the Decay of Strong Magnetization in Global Disc Simulations with Toroidal Fields

Authors: P. Chris Fragile, Aleksander Sadowski

Abstract: Strong magnetization in accretion discs could resolve a number of outstanding issues related to stability and state transitions in low-mass X-ray binaries. However, it is unclear how real discs become strongly magnetized and, even if they do, whether they can remain in such a state. In this paper, we address the latter issue through a pair of global disc simulations. Here, we only consider cases o… ▽ More Strong magnetization in accretion discs could resolve a number of outstanding issues related to stability and state transitions in low-mass X-ray binaries. However, it is unclear how real discs become strongly magnetized and, even if they do, whether they can remain in such a state. In this paper, we address the latter issue through a pair of global disc simulations. Here, we only consider cases of initially purely toroidal magnetic fields contained entirely within a compact torus. We find that, over only a few tens of orbital periods, the magnetization of an initially strongly magnetized disc, $P_\mathrm{mag}/P_\mathrm{gas} \ge 10$, drops to $\lesssim 0.1$, similar to the steady-state value reached in initially weakly magnetized discs. This is consistent with recent shearing box simulations with initially strong toroidal fields, the robust conclusion being that strongly magnetized toroidal fields can not be locally self-sustaining. These results appear to leave net poloidal flux or extended radial fields as the only avenues for establishing strongly magnetized discs, ruling out the thermal collapse scenario. △ Less

Submitted 27 January, 2017; v1 submitted 4 January, 2017; originally announced January 2017.

Comments: 6 pages, 7 figures, accepted to MNRAS

Numerical Simulations of a Jet-Cloud Collision and Starburst: Application to Minkowski's Object

Authors: P. Chris Fragile, Peter Anninos, Steve Croft, Mark Lacy, Jason W. L. Witry

Abstract: We present results of three-dimensional, multi-physics simulations of an AGN jet colliding with an intergalactic cloud. The purpose of these simulations is to assess the degree of "positive feedback," i.e. jet-induced star formation, that results. We have specifically tailored our simulation parameters to facilitate comparison with recent observations of Minkowski's Object (M.O.), a stellar nurser… ▽ More We present results of three-dimensional, multi-physics simulations of an AGN jet colliding with an intergalactic cloud. The purpose of these simulations is to assess the degree of "positive feedback," i.e. jet-induced star formation, that results. We have specifically tailored our simulation parameters to facilitate comparison with recent observations of Minkowski's Object (M.O.), a stellar nursery located at the termination point of a radio jet coming from galaxy NGC 541. As shown in our simulations, such a collision triggers shocks which propagate around and through the cloud. These shocks condense the gas and under the right circumstances may trigger cooling instabilities, creating runaway increases in density, to the point that individual clumps can become Jeans unstable. Our simulations provide information about the expected star formation rate, total mass converted to \ion{H}{1}, H$_2$, and stars, and the relative velocity of the stars and gas. Our results confirm the possibility of jet-induced star formation, and agree well with the observations of M.O. △ Less

Submitted 26 October, 2017; v1 submitted 30 December, 2016; originally announced January 2017.

Comments: 11 pages, 6 figures, ApJ accepted version

Three-dimensional, global, radiative GRMHD simulations of a thermally unstable disc

Authors: B. Mishra, P. C. Fragile, L. C. Johnson, W. Kluźniak

Abstract: We present results of a set of three-dimensional, general relativistic radiation magnetohydrodynamics simulations of thin accretion discs around a non-rotating black hole to test their thermal stability. We consider two cases, one that is initially radiation pressure dominated and expected to be thermally unstable and another that is initially gas-pressure dominated and expected to remain stable.… ▽ More We present results of a set of three-dimensional, general relativistic radiation magnetohydrodynamics simulations of thin accretion discs around a non-rotating black hole to test their thermal stability. We consider two cases, one that is initially radiation pressure dominated and expected to be thermally unstable and another that is initially gas-pressure dominated and expected to remain stable. Indeed, we find that cooling dominates over heating in the radiation pressure dominated model, causing the disc to collapse vertically on roughly the local cooling timescale. We also find that heating and cooling within the disc have a different dependence on the mid-plane pressure, a prerequisite of thermal instability. Comparison of our data with the relevant thin-disc thermal equilibrium curve suggests that our disc may be headed for the thermally stable, gas-pressure-dominated branch. However, because the disc collapses to the point that we are no longer able to resolve it, we had to terminate the simulation. On the other hand, the gas pressure dominated model, which was run for twice as long as the radiation pressure dominated one, remains stable, with heating and cooling roughly in balance. Finally, the radiation pressure dominated simulation shows some evidence of viscous instability. The strongest evidence is in plots of surface density, which show the disc breaking up into rings. △ Less

Submitted 7 September, 2016; v1 submitted 13 March, 2016; originally announced March 2016.

Comments: 12 pages, 20 figures, accepted in MNRAS

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

Authors: P. Chris Fragile, Odele Straub, Omer Blaes

Abstract: 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 occu… ▽ More 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. △ Less

Submitted 17 June, 2016; v1 submitted 25 February, 2016; originally announced February 2016.

Comments: 7 pages, 4 figures, accepted to MNRAS

Quasi-periodic oscillations from relativistic hydrodynamical slender tori

Authors: B. Mishra, F. H. Vincent, A. Manousakis, P. C. Fragile, T. Paumard, W. Kluźniak

Abstract: We simulate an oscillating purely hydrodynamical torus with constant specific angular mo- mentum around a Schwarzschild black hole. The goal is to search for quasi-periodic oscil- lations (QPOs) in the light curve of the torus. The initial torus setup is subjected to radial, vertical and diagonal (combination of radial and vertical) velocity perturbations. The hydro- dynamical simulations are perf… ▽ More We simulate an oscillating purely hydrodynamical torus with constant specific angular mo- mentum around a Schwarzschild black hole. The goal is to search for quasi-periodic oscil- lations (QPOs) in the light curve of the torus. The initial torus setup is subjected to radial, vertical and diagonal (combination of radial and vertical) velocity perturbations. The hydro- dynamical simulations are performed using the general relativistic magnetohydrodynamics code Cosmos++ and ray-traced using the GYOTO code. We found that a horizontal velocity perturbation triggers the radial and plus modes, while a vertical velocity perturbation trig- gers the vertical and X modes. The diagonal perturbation gives a combination of the modes triggered in the radial and vertical perturbations. △ Less

Submitted 26 October, 2015; originally announced October 2015.

Comments: 9 pages, 6 figures, submitted in MNRAS

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

Authors: Geoffrey C. Bower, Sera Markoff, Jason Dexter, Mark A. Gurwell, James M. Moran, Andreas Brunthaler, Heino Falcke, P. Chris Fragile, Dipankar Maitra, Dan Marrone, Alison Peck, Anthony Rushton, Melvyn C. H. Wright

Abstract: 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… ▽ More 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. △ Less

Submitted 23 February, 2015; originally announced February 2015.

Comments: Accepted for publication in ApJ

Local Stability of Strongly Magnetized Black Hole Tori

Authors: Maciek Wielgus, P. Chris Fragile, Ziming Wang, Julia Wilson

Abstract: 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, regard… ▽ More 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. △ Less

Submitted 15 December, 2014; originally announced December 2014.

Comments: 9 pages, 11 figures, accepted in MNRAS

Journal ref: MNRAS 447 (4): 3593-3601, 2015

The X-ray Flux Distribution of Sagittarius A* as Seen by Chandra

Authors: J. Neilsen, S. Markoff, M. A. Nowak, J. Dexter, G. Witzel, N. Barrière, Y. Li, F. K. Baganoff, N. Degenaar, P. C. Fragile, C. Gammie, A. Goldwurm, N. Grosso, D. Haggard

Abstract: We present a statistical analysis of the X-ray flux distribution of Sgr A* from the Chandra X-ray Observatory's 3 Ms Sgr A* X-ray Visionary Project (XVP) in 2012. Our analysis indicates that the observed X-ray flux distribution can be decomposed into a steady quiescent component, represented by a Poisson process with rate $Q=(5.24\pm0.08)\times10^{-3}$ cts s$^{-1},$ and a variable component, repre… ▽ More We present a statistical analysis of the X-ray flux distribution of Sgr A* from the Chandra X-ray Observatory's 3 Ms Sgr A* X-ray Visionary Project (XVP) in 2012. Our analysis indicates that the observed X-ray flux distribution can be decomposed into a steady quiescent component, represented by a Poisson process with rate $Q=(5.24\pm0.08)\times10^{-3}$ cts s$^{-1},$ and a variable component, represented by a power law process ($dN/dF\propto F^{-ξ},$ $ξ=1.92_{-0.02}^{+0.03}$). This slope matches our recently-reported distribution of flare luminosities. The variability may also be described by a log-normal process with a median unabsorbed 2-8 keV flux of $1.8^{+0.9}_{-0.6}\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$ and a shape parameter $σ=2.4\pm0.2,$ but the power law provides a superior description of the data. In this decomposition of the flux distribution, all of the intrinsic X-ray variability of Sgr A* (spanning at least three orders of magnitude in flux) can be attributed to flaring activity, likely in the inner accretion flow. We confirm that at the faint end, the variable component contributes ~10% of the apparent quiescent flux, as previously indicated by our statistical analysis of X-ray flares in these Chandra observations. Our flux distribution provides a new and important observational constraint on theoretical models of Sgr A*, and we use simple radiation models to explore the extent to which a statistical comparison of the X-ray and infrared can provide insights into the physics of the X-ray emission mechanism. △ Less

Submitted 9 December, 2014; originally announced December 2014.

Comments: 13 pages, 7 figures, accepted for publication in ApJ. Comments welcome

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

Authors: P. Chris Fragile, Ally Olejar, Peter Anninos

Abstract: 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… ▽ More 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. △ Less

Submitted 19 August, 2014; originally announced August 2014.

Comments: 18 pages, 9 figures

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

Authors: Viacheslav V. Zhuravlev, Pavel B. Ivanov, P. Chris Fragile, Danilo Morales Teixeira

Abstract: 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 d… ▽ More 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. △ Less

Submitted 9 September, 2014; v1 submitted 20 June, 2014; originally announced June 2014.

Comments: 15 pages, 11 figures, accepted to Astrophysical Journal

Conservative GRMHD Simulations of Moderately Thin, Tilted Accretion Disks

Authors: Danilo Morales Teixeira, P. Chris Fragile, Viacheslav V. Zhuravlev, Pavel B. Ivanov

Abstract: 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… ▽ More 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. △ Less

Submitted 4 September, 2014; v1 submitted 20 June, 2014; originally announced June 2014.

Comments: 13 pages, 14 figures, accepted to Astrophysical Journal

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

Authors: P. Chris Fragile, Peter Anninos, Stephen D. Murray

Abstract: 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,… ▽ More 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. △ Less

Submitted 2 January, 2014; originally announced January 2014.

Comments: Proceedings of IAU Symposium No. 303, The Galactic Center: Feeding and Feedback in a Normal Galactic Nucleus

The 3 Ms Chandra Campaign on Sgr A*: A Census of X-ray Flaring Activity from the Galactic Center

Authors: J. Neilsen, M. A. Nowak, C. Gammie, J. Dexter, S. Markoff, D. Haggard, S. Nayakshin, Q. D. Wang, N. Grosso, D. Porquet, J. A. Tomsick, N. Degenaar, P. C. Fragile, J. C. Houck, R. Wijnands, J. M. Miller, F. K. Baganoff

Abstract: Over the last decade, X-ray observations of Sgr A* have revealed a black hole in a deep sleep, punctuated roughly once per day by brief flares. The extreme X-ray faintness of this supermassive black hole has been a long-standing puzzle in black hole accretion. To study the accretion processes in the Galactic Center, Chandra (in concert with numerous ground- and space-based observatories) undertook… ▽ More Over the last decade, X-ray observations of Sgr A* have revealed a black hole in a deep sleep, punctuated roughly once per day by brief flares. The extreme X-ray faintness of this supermassive black hole has been a long-standing puzzle in black hole accretion. To study the accretion processes in the Galactic Center, Chandra (in concert with numerous ground- and space-based observatories) undertook a 3 Ms campaign on Sgr A* in 2012. With its excellent observing cadence, sensitivity, and spectral resolution, this Chandra X-ray Visionary Project (XVP) provides an unprecedented opportunity to study the behavior of the closest supermassive black hole. We present a progress report from our ongoing study of X-ray flares, including the brightest flare ever seen from Sgr A*. Focusing on the statistics of the flares and the quiescent emission, we discuss the physical implications of X-ray variability in the Galactic Center. △ Less

Submitted 26 November, 2013; originally announced November 2013.

Comments: Proceedings of IAU Symposium No. 303, The Galactic Center: Feeding and Feedback in a Normal Galactic Nucleus

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

Authors: Aleksey Generozov, Omer Blaes, P. Chris Fragile, Ken B. Henisey

Abstract: 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… ▽ More 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. △ Less

Submitted 21 November, 2013; originally announced November 2013.

Comments: Accepted to The Astrophysical Journal

A Chandra/HETGS Census of X-ray Variability From Sgr A* During 2012

Authors: J. Neilsen, M. A. Nowak, C. Gammie, J. Dexter, S. Markoff, D. Haggard, S. Nayakshin, Q. D. Wang, N. Grosso, D. Porquet, J. A. Tomsick, N. Degenaar, P. C. Fragile, R. Wijnands, J. M. Miller, F. K. Baganoff

Abstract: We present the first systematic analysis of the X-ray variability of Sgr A* during the Chandra X-ray Observatory's 2012 Sgr A* X-ray Visionary Project (XVP). With 38 High Energy Transmission Grating Spectrometer (HETGS) observations spaced an average of 7 days apart, this unprecedented campaign enables detailed study of the X-ray emission from this supermassive black hole at high spatial, spectral… ▽ More We present the first systematic analysis of the X-ray variability of Sgr A* during the Chandra X-ray Observatory's 2012 Sgr A* X-ray Visionary Project (XVP). With 38 High Energy Transmission Grating Spectrometer (HETGS) observations spaced an average of 7 days apart, this unprecedented campaign enables detailed study of the X-ray emission from this supermassive black hole at high spatial, spectral and timing resolution. In 3 Ms of observations, we detect 39 X-ray flares from Sgr A*, lasting from a few hundred seconds to approximately 8 ks, and ranging in 2-10 keV luminosity from ~1e34 erg/s to 2e35 erg/s. Despite tentative evidence for a gap in the distribution of flare peak count rates, there is no evidence for X-ray color differences between faint and bright flares. Our preliminary X-ray flare luminosity distribution dN/dL is consistent with a power law with index -1.9 (+0.3 -0.4); this is similar to some estimates of Sgr A*'s NIR flux distribution. The observed flares contribute one-third of the total X-ray output of Sgr A* during the campaign, and as much as 10% of the quiescent X-ray emission could be comprised of weak, undetected flares, which may also contribute high-frequency variability. We argue that flares may be the only source of X-ray emission from the inner accretion flow. △ Less

Submitted 22 July, 2013; originally announced July 2013.

Comments: Version resubmitted to the Astrophysical Journal. 15 pages, 5 figures. Comments welcome

Current Status of Simulations

Authors: P. Chris Fragile

Abstract: 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 so… ▽ More 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. △ Less

Submitted 19 April, 2013; originally announced April 2013.

Comments: 15 pages, 14 figures, to appear in the proceedings of the ISSI-Bern workshop on "The Physics of Accretion onto Black Holes" (8-12 October 2012)

Variability from Nonaxisymmetric Fluctuations Interacting with Standing Shocks in Tilted Black Hole Accretion Disks

Authors: Ken B. Henisey, Omer M. Blaes, P. Chris Fragile

Abstract: We study the spatial and temporal behavior of fluid in fully three-dimensional, general relativistic, magnetohydrodynamical simulations of both tilted and untilted black hole accretion flows. We uncover characteristically greater variability in tilted simulations at frequencies similar to those predicted by the formalism of trapped modes, but ultimately conclude that its spatial structure is incon… ▽ More We study the spatial and temporal behavior of fluid in fully three-dimensional, general relativistic, magnetohydrodynamical simulations of both tilted and untilted black hole accretion flows. We uncover characteristically greater variability in tilted simulations at frequencies similar to those predicted by the formalism of trapped modes, but ultimately conclude that its spatial structure is inconsistent with a modal interpretation. We find instead that previously identified, transient, over-dense clumps orbiting on roughly Keplerian trajectories appear generically in our global simulations, independent of tilt. Associated with these fluctuations are acoustic spiral waves interior to the orbits of the clumps. We show that the two nonaxisymmetric standing shock structures that exist in the inner regions of these tilted flows effectively amplify the variability caused by these spiral waves to markedly higher levels than in untilted flows, which lack standing shocks. Our identification of clumps, spirals, and spiral-shock interactions in these fully general relativistic, magnetohydrodynamical simulations suggests that these features may be important dynamical elements in models which incorporate tilt as a way to explain the observed variability in black hole accretion flows. △ Less

Submitted 14 November, 2012; v1 submitted 9 November, 2012; originally announced November 2012.

Comments: accepted in ApJ October 1, 2012

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

Authors: Samia Drappeau, Salomé Dibi, Jason Dexter, Sera Markoff, P. Chris Fragile

Abstract: 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 p… ▽ More 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. △ Less

Submitted 20 September, 2012; originally announced September 2012.

Comments: Submitted to MNRAS. 13 pages, 15 figures

Report number: NSF-KITP-12-179

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

Authors: Peter Anninos, P. Chris Fragile, Julia Wilson, Stephen D. Murray

Abstract: 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… ▽ More 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. △ Less

Submitted 2 January, 2014; v1 submitted 7 September, 2012; originally announced September 2012.

Comments: 14 pages, 8 figures, final accepted version

Report number: NSF-KITP-12-058

Journal ref: ApJ, 759, 132 (2012)

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

Authors: Salomé Dibi, Samia Drappeau, P. Chris Fragile, Sera Markoff, Jason Dexter

Abstract: 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 ignori… ▽ More 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. △ Less

Submitted 18 June, 2012; originally announced June 2012.

Comments: Resubmitted to MNRAS, including modifications in response to referee report. 13 pages, 15 figures

Report number: NSF-KITP-12-053

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

Authors: P. Chris Fragile, Julia Wilson, Marco Rodriguez

Abstract: 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 genera… ▽ More 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. △ Less

Submitted 10 May, 2012; v1 submitted 1 May, 2012; originally announced May 2012.

Comments: 9 pages, 7 figures, accepted for publication in MNRAS

Report number: NSF-KITP-12-052

Journal ref: MNRAS, 424, 524 (2012)

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

Authors: P. Chris Fragile, Anna Gillespie, Timothy Monahan, Marco Rodriguez, Peter Anninos

Abstract: 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 s… ▽ More 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. △ Less

Submitted 8 May, 2012; v1 submitted 24 April, 2012; originally announced April 2012.

Comments: 15 pages, 6 figures, ApJS accepted version

Report number: NSF-KITP-12-054

Journal ref: ApJS, 201, 9 (2012)

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

Authors: Jason Dexter, P. Chris Fragile

Abstract: 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 a… ▽ More 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). △ Less

Submitted 19 April, 2012; originally announced April 2012.

Comments: 17 pages, 18 figures, submitted to MNRAS; for movies and version with high-res figures see http://astro.berkeley.edu/~jdexter/tiltedsgra

Report number: NSF-KITP-12-057