Source author record

Jason Dexter

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

40works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

40 published item(s)

preprint2025arXiv

A Convolutional Neural Network for the Recovery of Transfer Functions From Velocity-Resolved Reverberation Mapping Data

One of the hallmarks of active galactic nuclei are that they are highly variable with time. In watching the spectra vary it has been observed that the emission-lines often appear to "reverberate" -- that is they vary in response to continuum variations assumed to originate close to the black hole. This critical observation underlies the reverberation mapping technique, an elegant physics experiment that has allowed us to characterize the environment around many supermassive black holes in nearby active galactic nuclei. Recent observations are of such quality that the response can be measured as a function of velocity across the emission-line, and in doing so we can construct velocity-delay maps that show the structure and physics of the gas in the broad-line region better than any other measurement to date. Unfortunately constructing such maps requires a deconvolution, and given that the data are often noisy and with gaps such deconvolutions are non-trivial. Here we present a novel deconvolution method for the recovery of velocity-delay maps using a custom convolutional neural network architecture, showcasing that such methods have great promise for the deconvolution of reverberation mapping data products. While we have designed this new method with the BLR in mind, in principle this technique could be applied to any reverberation deconvolution problem, including in the accretion disk and torus.

preprint2022arXiv

Millimeter light curves of Sagittarius A* observed during the 2017 Event Horizon Telescope campaign

The Event Horizon Telescope (EHT) observed the compact radio source, Sagittarius A* (Sgr A*), in the Galactic Center on 2017 April 5-11 in the 1.3 millimeter wavelength band. At the same time, interferometric array data from the Atacama Large Millimeter/submillimeter Array and the Submillimeter Array were collected, providing Sgr A* light curves simultaneous with the EHT observations. These data sets, complementing the EHT very-long-baseline interferometry, are characterized by a cadence and signal-to-noise ratio previously unattainable for Sgr A* at millimeter wavelengths, and they allow for the investigation of source variability on timescales as short as a minute. While most of the light curves correspond to a low variability state of Sgr A*, the April 11 observations follow an X-ray flare, and exhibit strongly enhanced variability. All of the light curves are consistent with a red noise process, with a power spectral density (PSD) slope measured to be between -2 and -3 on timescales between 1 min and several hours. Our results indicate a steepening of the PSD slope for timescales shorter than 0.3 h. The spectral energy distribution is flat at 220 GHz and there are no time-lags between the 213 and 229 GHz frequency bands, suggesting low optical depth for the event horizon scale source. We characterize Sgr A*'s variability, highlighting the different behavior observed just after the X-ray flare, and use Gaussian process modeling to extract a decorrelation timescale and a PSD slope. We also investigate the systematic calibration uncertainties by analyzing data from independent data reduction pipelines.

preprint2022arXiv

Nonthermal emission from the plunging region: a model for the high-energy tail of black hole X-ray binary soft states

X-ray binaries exhibit a soft spectral state comprising thermal blackbody emission at 1 keV and a power-law tail above 10 keV. Empirical models fit the high-energy power-law tail to radiation from a nonthermal electron distribution, but the physical location of the nonthermal electrons and the reason for their power-law index and high-energy cut-off are still largely unknown. Here, we propose that the nonthermal electrons originate from within the black hole's innermost stable circular orbit (the ''plunging region''). Using an analytic model for the plunging region dynamics and electron distribution function properties from particle-in-cell simulations, we outline a steady-state model that can reproduce the observed spectral features. In particular, our model reproduces photon indices of $Γ\gtrsim2$ and power-law luminosities on the order of a few percent of the disk luminosity for strong magnetic fields, consistent with observations of the soft state. Because the emission originates so close to the black hole, we predict that the power-law luminosity should strongly depend on the system inclination angle and black hole spin. This model could be extended to the power-law tails observed above 400 keV in the hard state of X-ray binaries.

preprint2022arXiv

Sgr A* X-ray flares from non-thermal particle acceleration in a magnetically arrested disc

Sgr A* exhibits flares in the near-infrared and X-ray bands, with the luminosity in these bands increasing by factors of $10-100$ for $\approx 60$ minutes. One of the models proposed to explain these flares is synchrotron emission of non-thermal particles accelerated by magnetic reconnection events in the accretion flow. We use the results from PIC simulations of magnetic reconnection to post-process 3D two-temperature GRMHD simulations of a magnetically arrested disc (MAD). We identify current sheets, retrieve their properties, estimate their potential to accelerate non-thermal particles and compute the expected non-thermal synchrotron emission. We find that the flux eruptions of MADs can provide suitable conditions for accelerating non-thermal particles to energies $γ_e \lesssim 10^6$ and producing simultaneous X-ray and near-infrared flares. For a suitable choice of current-sheet parameters and a simpified synchrotron cooling prescription, the model can simultaneously reproduce the quiescent and flaring X-ray luminosities as well as the X-ray spectral shape. While the near-infrared flares are mainly due to an increase in the temperature near the black hole during the MAD flux eruptions, the X-ray emission comes from narrow current sheets bordering highly magnetized, low-density regions near the black hole and equatorial current sheets where the flux on the black hole reconnects. As a result, not all infrared flares are accompanied by X-ray ones. The non-thermal flaring emission can extend to very hard ($\lesssim 100$ keV) X-ray energies.

preprint2022arXiv

What really makes an accretion disc MAD

Magnetically arrested accretion discs (MADs) around black holes (BH) have the potential to stimulate the production of powerful jets and account for recent ultra-high-resolution observations of BH environments. Their main properties are usually attributed to the accumulation of dynamically significant net magnetic (vertical) flux throughout the arrested region, which is then regulated by interchange instabilities. Here we propose instead that it is mainly a dynamically important {\it toroidal} field -- the result of dynamo action triggered by the significant but still relatively weak vertical field -- that defines and regulates the properties of MADs. We suggest that rapid convection-like instabilities, involving interchange of toroidal flux tubes and operating concurrently with the magnetorotatonal instability (MRI), can regulate the structure of the disc and the escape of net flux. We generalize the convective stability criteria and disc structure equations to include the effects of a strong toroidal field and show that convective flows could be driven towards two distinct marginally stable states, one of which we associate with MADs. We confirm the plausibility of our theoretical model by comparing its quantitative predictions to simulations of both MAD and SANE (strongly magnetized but not "arrested") discs, and suggest a set of criteria that could help to distinguish MADs from other accretion states. Contrary to previous claims in the literature, we argue that MRI is not suppressed in MADs and is probably responsible for the existence of the strong toroidal field.

preprint2020arXiv

Infrared interferometry to spatially and spectrally resolve jets in X-ray binaries

Infrared interferometry is a new frontier for precision ground based observing, with new instrumentation achieving milliarcsecond (mas) spatial resolutions for faint sources, along with astrometry on the order of 10 microarcseconds. This technique has already led to breakthroughs in the observations of the supermassive black hole at the Galactic centre and its orbiting stars, AGN, and exo-planets, and can be employed for studying X-ray binaries (XRBs), microquasars in particular. Beyond constraining the orbital parameters of the system using the centroid wobble and spatially resolving jet discrete ejections on mas scales, we also propose a novel method to discern between the various components contributing to the infrared bands: accretion disk, jets and companion star. We demonstrate that the GRAVITY instrument on the Very Large Telescope Interferometer (VLTI) should be able to detect a centroid shift in a number of sources, opening a new avenue of exploration for the myriad of transients expected to be discovered in the coming decade of radio all-sky surveys. We also present the first proof-of-concept GRAVITY observation of a low-mass X-ray binary transient, MAXI J1820+070, to search for extended jets on mas scales. We place the tightest constraints yet via direct imaging on the size of the infrared emitting region of the compact jet in a hard state XRB.

preprint2020arXiv

Magnetic flux inversion in a peculiar changing look AGN

We argue that the changing-look event in the active galactic nucleus 1ES 1927+654, followed by a dip of 3 orders of magnitude in the X-ray luminosity, is controlled by a change in the accretion rate and an inversion of magnetic flux in a magnetically arrested disk (MAD). Before the changing-look event, strong magnetic flux on the black hole powers X-ray emission via the Blandford-Znajek process while the UV emission is produced by a radiatively inefficient magnetized disk. An advection event, bringing flux of the opposite polarity, propagates inward leading, first, to a rise in the UV/optical luminosity and, then, to a dip in the X-ray luminosity when it reaches the black hole. We estimate the timescale for magnetic flux advection and find that the observed timescale between the beginning of the changing-look event and the minimum in the X-ray luminosity, $\approx200$ days, is in agreement with the time needed to cancel the magnetic flux in a MAD extending to $\approx180\:r_g$. Although flux inversion events might be rare due to the large ratio of flux-to-mass that is needed, we argue that AGN showing an unusually high ratio of X-ray to UV luminosity are prime candidates for such events. We also suggest that similar events may lead to jet interruptions in radio-loud objects.

preprint2020arXiv

The Effects of Tilt on the Images of Black Hole Accretion Flows

We analyze two 3D general-relativistic magnetohydrodynamic accretion simulations in the context of how they would manifest in Event Horizon Telescope (EHT) observations of supermassive black holes. The two simulations differ only in whether the initial angular momentum of the plasma is aligned with the rapid (a = 0.9) spin of the black hole. Both have low net magnetic flux. Ray tracing is employed to generate resolved images of the synchrotron emission. When using parameters appropriate for Sgr A* and assuming a viewing angle aligned with the black hole spin, we find the most prominent difference is that the central shadow in the image is noticeably eccentric in tilted models, with the ring of emission relatively unchanged. Applying this procedure to M87 with a viewing angle based on the large-scale jet, we find that adding tilt increases the angular size of the ring for fixed black hole mass and distance, while at the same time increasing the number of bright spots in the image. Our findings illustrate observable features that can distinguish tilted from aligned flows. They also show that tilted models can be viable for M87, and that not accounting for tilt can bias inferences of physical parameters. Future modeling of horizon-scale observations should account for potential angular momentum misalignment, which is likely generic at the low accretion rates appropriate for EHT targets.

preprint2019arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Accretion and Broad Emission Line Physics from a Hypervariable Quasar

We analyze extensive spectroscopic and photometric data of the hypervariable quasar SDSS J131424+530527 (RMID 017) at z=0.456, an optical "changing look" quasar from the Sloan Digital Sky Survey Reverberation Mapping project that increased in optical luminosity by a factor of 10 between 2014 and 2017. The observed broad emission lines all respond in luminosity and width to the changing optical continuum, as expected for photoionization in a stratified, virialized broad emission line region. The luminosity changes therefore result from intrinsic changes in accretion power rather than variable obscuration. The variability is continuous and apparently stochastic, disfavoring an origin as a discrete event such as a tidal disruption flare or microlensing event. It is coordinated on day timescales with blue leading red, consistent with reprocessing powering the entire optical SED. We show that this process cannot work in a standard thin disk geometry on energetic grounds, and would instead require a large covering factor reprocessor. Disk instability models could potentially also explain the data, provided that the instability sets in near the inner radius of a geometrically thick accretion disk.

preprint2016arXiv

A public code for general relativistic, polarised radiative transfer around spinning black holes

Ray tracing radiative transfer is a powerful method for comparing theoretical models of black hole accretion flows and jets with observations. We present a public code, grtrans, for carrying out such calculations in the Kerr metric, including the full treatment of polarised radiative transfer and parallel transport along geodesics. The code is written in Fortran 90 and efficiently parallelises with OpenMP, and the full code and several components have Python interfaces. We describe several tests which are used for verifiying the code, and we compare the results for polarised thin accretion disc and semi-analytic jet problems with those from the literature as examples of its use. Along the way, we provide accurate fitting functions for polarised synchrotron emission and transfer coefficients from thermal and power law distribution functions, and compare results from numerical integration and quadrature solutions of the polarised radiative transfer equations. We also show that all transfer coefficients can play an important role in predicted images and polarisation maps of the Galactic center black hole, Sgr A*, at submillimetre wavelengths.

preprint2016arXiv

Persistent Asymmetric Structure of Sagittarius A* on Event Horizon Scales

The Galactic Center black hole Sagittarius A* (Sgr A*) is a prime observing target for the Event Horizon Telescope (EHT), which can resolve the 1.3 mm emission from this source on angular scales comparable to that of the general relativistic shadow. Previous EHT observations have used visibility amplitudes to infer the morphology of the millimeter-wavelength emission. Potentially much richer source information is contained in the phases. We report on 1.3 mm phase information on Sgr A* obtained with the EHT on a total of 13 observing nights over 4 years. Closure phases, the sum of visibility phases along a closed triangle of interferometer baselines, are used because they are robust against phase corruptions introduced by instrumentation and the rapidly variable atmosphere. The median closure phase on a triangle including telescopes in California, Hawaii, and Arizona is nonzero. This result conclusively demonstrates that the millimeter emission is asymmetric on scales of a few Schwarzschild radii and can be used to break 180-degree rotational ambiguities inherent from amplitude data alone. The stability of the sign of the closure phase over most observing nights indicates persistent asymmetry in the image of Sgr A* that is not obscured by refraction due to interstellar electrons along the line of sight.

preprint2016arXiv

Young and Millisecond Pulsar GeV Gamma-ray Fluxes from the Galactic Center and Beyond

Gamma-ray observations have shown pulsars to be efficient converters of rotational energy into GeV photons and it is of wide-ranging interest to determine their contribution to the gamma-ray background. We arrive at flux predictions from both the young (<~ Myr) and millisecond (~Gyr) Galactic pulsar populations. We find that unresolved pulsars can yield both a significant fraction of the excess GeV gamma rays near the Galactic Center and an inverse Compton flux in the inner kpc similar to that inferred by Fermi. We compare models of the young pulsar population and millisecond pulsar population to constraints from gamma-ray and radio observations. Overall, we find that the young pulsars should outnumber millisecond pulsars as unassociated gamma-ray point sources in this region. The number of young radio pulsars discovered near the Galactic Center is in agreement with our model of the young pulsar population. Deeper radio observations at higher latitudes can constrain the total gamma-ray emission from both young and millisecond pulsars from the inner galaxy. While this is a step towards better understanding of pulsars, cosmic rays in the Milky Way, and searches for dark matter, we also discuss a few interesting puzzles that arise from the underlying physics of pulsar emission and evolution.

preprint2015arXiv

230 GHz VLBI observations of M87: event-horizon-scale structure at the enhanced very-high-energy $\rm γ$-ray state in 2012

We report on 230 GHz (1.3 mm) VLBI observations of M87 with the Event Horizon Telescope using antennas on Mauna Kea in Hawaii, Mt. Graham in Arizona and Cedar Flat in California. For the first time, we have acquired 230 GHz VLBI interferometric phase information on M87 through measurement of closure phase on the triangle of long baselines. Most of the measured closure phases are consistent with 0$^{\circ}$ as expected by physically-motivated models for 230 GHz structure such as jet models and accretion disk models. The brightness temperature of the event-horizon-scale structure is $\sim 1 \times 10^{10}$ K derived from the compact flux density of $\sim 1$ Jy and the angular size of $\sim 40 $ $\rm μ$as $\sim$ 5.5 $R_{\rm s}$, which is broadly consistent with the peak brightness of the radio cores at 1-86 GHz located within $\sim 10^2$ $R_{\rm s}$. Our observations occurred in the middle of an enhancement in very-high-energy (VHE) $\rm γ$-ray flux, presumably originating in the vicinity of the central black hole. Our measurements, combined with results of multi-wavelength observations, favor a scenario in which the VHE region has an extended size of $\sim$20-60 $R_{\rm s}$.

preprint2015arXiv

A Black Hole Mass-Variability Time Scale Correlation at Submillimeter Wavelengths

We analyze the light curves of 413 radio sources at submillimeter wavelengths using data from the Submillimeter Array calibrator database. The database includes more than 20,000 observations at 1.3 and 0.8 mm that span 13 years. We model the light curves as a damped random walk and determine a characteristic time scale $τ$ at which the variability amplitude saturates. For the vast majority of sources, primarily blazars and BL Lac objects, we find only lower limits on $τ$. For two nearby low luminosity active galactic nuclei, M81 and M87, however, we measure $τ=1.6^{+3.0}_{-0.9}$ days and $τ=45^{+61}_{-24}$ days, respectively ($2σ$ errors). Including the previously measured $τ=0.33\pm 0.16$ days for Sgr A*, we show an approximately linear correlation between $τ$ and black hole mass for these nearby LLAGN. Other LLAGN with spectra that peak in the submm are expected to follow this correlation. These characteristic time scales are comparable to the minimum time scale for emission processes close to an event horizon, and suggest that the underlying physics may be independent of black hole mass, accretion rate, and jet luminosity.

preprint2015arXiv

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

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

preprint2015arXiv

Resolved Magnetic-Field Structure and Variability Near the Event Horizon of Sagittarius A*

Near a black hole, differential rotation of a magnetized accretion disk is thought to produce an instability that amplifies weak magnetic fields, driving accretion and outflow. These magnetic fields would naturally give rise to the observed synchrotron emission in galaxy cores and to the formation of relativistic jets, but no observations to date have been able to resolve the expected horizon-scale magnetic-field structure. We report interferometric observations at 1.3-millimeter wavelength that spatially resolve the linearly polarized emission from the Galactic Center supermassive black hole, Sagittarius A*. We have found evidence for partially ordered fields near the event horizon, on scales of ~6 Schwarzschild radii, and we have detected and localized the intra-hour variability associated with these fields.

preprint2015arXiv

Statistics of X-ray flares of Sagittarius A*: evidence for solar-like self-organized criticality phenomenon

X-ray flares have routinely been observed from the supermassive black hole, Sagittarius A$^\star$ (Sgr A$^\star$), at our Galactic center. The nature of these flares remains largely unclear, despite of many theoretical models. In this paper, we study the statistical properties of the Sgr A$^\star$ X-ray flares, by fitting the count rate (CR) distribution and the structure function (SF) of the light curve with a Markov Chain Monte Carlo (MCMC) method. With the 3 million second \textit{Chandra} observations accumulated in the Sgr A$^\star$ X-ray Visionary Project, we construct the theoretical light curves through Monte Carlo simulations. We find that the $2-8$ keV X-ray light curve can be decomposed into a quiescent component with a constant count rate of $\sim6\times10^{-3}~$count s$^{-1}$ and a flare component with a power-law fluence distribution $dN/dE\propto E^{-α_{\rm E}}$ with $α_{\rm E}=1.65\pm0.17$. The duration-fluence correlation can also be modelled as a power-law $T\propto E^{α_{\rm ET}}$ with $α_{\rm ET} < 0.55$ ($95\%$ confidence). These statistical properties are consistent with the theoretical prediction of the self-organized criticality (SOC) system with the spatial dimension $S = 3$. We suggest that the X-ray flares represent plasmoid ejections driven by magnetic reconnection (similar to solar flares) in the accretion flow onto the black hole.

preprint2015arXiv

X-Ray Polarimetry with the Polarization Spectroscopic Telescope Array (PolSTAR)

This paper describes the Polarization Spectroscopic Telescope Array (PolSTAR), a mission proposed to NASA's 2014 Small Explorer (SMEX) announcement of opportunity. PolSTAR measures the linear polarization of 3-50 keV (requirement; goal: 2.5-70 keV) X-rays probing the behavior of matter, radiation and the very fabric of spacetime under the extreme conditions close to the event horizons of black holes, as well as in and around magnetars and neutron stars. The PolSTAR design is based on the technology developed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission launched in June 2012. In particular, it uses the same X-ray optics, extendable telescope boom, optical bench, and CdZnTe detectors as NuSTAR. The mission has the sensitivity to measure ~1% linear polarization fractions for X-ray sources with fluxes down to ~5 mCrab. This paper describes the PolSTAR design as well as the science drivers and the potential science return.

preprint2015arXiv

Young Pulsars and the Galactic Center GeV Gamma-ray Excess

Studies of Fermi data indicate an excess of GeV gamma rays around the Galactic center (GC), possibly due to dark matter. We show that young gamma-ray pulsars can yield a similar signal. First, a high concentration of GC supernovae naturally leads to a population of kicked pulsars symmetric about the GC. Second, while very-young pulsars with soft spectra reside near the Galactic plane, pulsars with spectra that have hardened with age accumulate at larger angles. This combination, including unresolved foreground pulsars, traces the morphology and spectrum of the Excess.

preprint2014arXiv

Evidence for Large Temperature Fluctuations in Quasar Accretion Disks From Spectral Variability

The well-known bluer-when-brighter trend observed in quasar variability is a signature of the complex processes in the accretion disk, and can be a probe of the quasar variability mechanism. Using a sample of 604 variable quasars with repeat spectra in SDSS-I/II, we construct difference spectra to investigate the physical causes of this bluer-when-brighter trend. The continuum of our composite difference spectrum is well-fit by a power-law, with a spectral index in excellent agreement with previous results. We measure the spectral variability relative to the underlying spectra of the quasars, which is independent of any extinction, and compare to model predictions. We show that our SDSS spectral variability results cannot be produced by global accretion rate fluctuations in a thin disk alone. However, we find that a simple model of a inhomogeneous disk with localized temperature fluctuations will produce power-law spectral variability over optical wavelengths. We show that the inhomogeneous disk will provide good fits to our observed spectral variability if the disk has large temperature fluctuations in many independently varying zones, in excellent agreement with independent constraints from quasar microlensing disk sizes, their strong UV spectral continuum, and single-band variability amplitudes. Our results provide an independent constraint on quasar variability models, and add to the mounting evidence that quasar accretion disks have large localized temperature fluctuations.

preprint2014arXiv

NuSTAR detection of high-energy X-ray emission and rapid variability from Sagittarius A* flares

Sagittarius A* harbors the supermassive black hole that lies at the dynamical center of our Galaxy. Sagittarius A* spends most of its time in a low luminosity emission state but flares frequently in the infrared and X-ray, increasing up to a few hundred fold in brightness for up to a few hours at a time. The physical processes giving rise to the X-ray flares are uncertain. Here we report the detection with the NuSTAR observatory in Summer and Fall 2012 of four low to medium amplitude X-ray flares to energies up to 79 keV. For the first time, we clearly see that the power-law spectrum of Sagittarius A* X-ray flares extends to high energy, with no evidence for a cut off. Although the photon index of the absorbed power-law fits are in agreement with past observations, we find a difference between the photon index of two of the flares (significant at the 95% confidence level). The spectra of the two brightest flares (~55 times quiescence in the 2-10 keV band) are compared to simple physical models in an attempt to identify the main X-ray emission mechanism, but the data do not allow us to significantly discriminate between them. However, we confirm the previous finding that the parameters obtained with synchrotron models are, for the X-ray emission, physically more reasonable than those obtained with inverse-Compton models. One flare exhibits large and rapid (< 100 s) variability, which, considering the total energy radiated, constrains the location of the flaring region to be within ~10 Schwarzschild radii of the black hole.

preprint2014arXiv

The Event Horizon Telescope: exploring strong gravity and accretion physics

The Event Horizon Telescope (EHT), a global sub-millimeter wavelength very long baseline interferometry array, is now resolving the innermost regions around the supermassive black holes Sgr A* and M87. Using black hole images from both simple geometric models and relativistic magnetohydrodynamical accretion flow simulations, we perform a variety of experiments to assess the promise of the EHT for studying strong gravity and accretion physics during the stages of its development. We find that (1) the addition of the LMT and ALMA along with upgraded instrumentation in the "Complete" stage of the EHT allow detection of the photon ring, a signature of Kerr strong gravity, for predicted values of its total flux; (2) the inclusion of coherently averaged closure phases in our analysis dramatically improves the precision of even the current array, allowing (3) significantly tighter constraints on plausible accretion models and (4) detections of structural variability at the levels predicted by the models. While observations at 345 GHz circumvent problems due to interstellar electron scattering in line-of-sight to the galactic center, short baselines provided by CARMA and/or the LMT could be required in order to constrain the overall shape of the accretion flow. Given the systematic uncertainties in the underlying models, using the full complement of two observing frequencies (230 and 345 GHz) and sources (Sgr A* and M87) may be critical for achieving transformative science with the EHT experiment.

preprint2014arXiv

The Peculiar Pulsar Population of the Central Parsec

Pulsars orbiting the Galactic center black hole, Sgr A*, would be potential probes of its mass, distance and spin, and may even be used to test general relativity. Despite predictions of large populations of both ordinary and millisecond pulsars in the Galactic center, none have been detected within 25 pc by deep radio surveys. One explanation has been that hyperstrong temporal scattering prevents pulsar detections, but the recent discovery of radio pulsations from a highly magnetized neutron star (magnetar) within 0.1 pc shows that the temporal scattering is much weaker than predicted. We argue that an intrinsic deficit in the ordinary pulsar population is the most likely reason for the lack of detections to date: a "missing pulsar problem" in the Galactic center. In contrast, we show that the discovery of a single magnetar implies efficient magnetar formation in the region. If the massive stars in the central parsec form magnetars rather than ordinary pulsars, their short lifetimes could explain the missing pulsars. Efficient magnetar formation could be caused by strongly magnetized progenitors, or could be further evidence of a top-heavy initial mass function. Furthermore, current high-frequency surveys should already be able to detect bright millisecond pulsars, given the measured degree of temporal scattering.

preprint2013arXiv

A geometric crescent model for black hole images

The Event Horizon Telescope (EHT), a global very long baseline interferometry array operating at millimetre wavelengths, is spatially resolving the immediate environments of black holes for the first time. The current observations of the Galactic center black hole, Sagittarius A* (Sgr A*), and M87 have been interpreted in terms of either geometric models (e.g., a symmetric Gaussian) or detailed calculations of the appearance of black hole accretion flows. The former are not physically motivated, while the latter are subject to large systematic uncertainties. Motivated by the dominant relativistic effects of Doppler beaming and gravitational lensing in many calculations, we propose a geometric crescent model for black hole images. We show that this simple model provides an excellent statistical description of the existing EHT data of Sgr A* and M87, superior to other geometric models for Sgr A*. It also qualitatively matches physically predicted models, bridging accretion theory and observation. Based on our results, we make predictions for the detectability of the black hole shadow, a signature of strong gravity, in future observations.

preprint2013arXiv

A model of the steep power law spectra and high-frequency quasi-periodic oscillations in luminous black hole X-ray binaries

We propose a new model of the steep power law state of luminous black hole X-ray binaries. The model uses the fact that at high luminosities, the inner radii of radiation pressure dominated accretion discs are expected to i) become effectively optically thin and ii) produce significant luminosities. The gas temperature therefore rises sharply inwards, producing local saturated Compton spectra with rapidly increasing peak energy. These spectra sum together to form a steep power law tail to the spectrum. A given photon energy on this tail corresponds to a narrow range in radius, so that local vertical oscillations of the disc naturally produce high quality high-frequency quasi-periodic oscillations (HFQPOs) in the hard X-ray band. The two lowest order modes have a robust frequency ratio of ~1.53. This model explains the appearance of steep power law spectra and HFQPOs at high luminosity, the 3:2 HFQPO frequency ratios, and their association with the power law spectral component. We predict an increase in QPO quality factor when the power spectrum is restricted to a narrower photon energy band, and an increase in HFQPO frequency at higher X-ray energies or lower luminosities. Future X-ray telescopes could detect additional HFQPOs from higher order modes. We demonstrate how this model could be used to measure black hole spin from HFQPOs, and qualitatively estimate the spin of GRO J1655-40 as a/M ~ 0.4 - 0.7.

preprint2013arXiv

An 8 hour characteristic time-scale in submillimetre light curves of Sagittarius A*

We compile and analyse long term (~10 year) submillimetre (1.3, 0.87, 0.43 mm, submm) wavelength light curves of the Galactic centre black hole, Sagittarius A*. The 0.87 and 0.43 mm data are taken from the literature, while the majority of the 1.3 mm light curve is from previously unpublished SMA and CARMA data. We use Monte Carlo simulations to show that on minute to few hour time-scales the variability is consistent with a red noise process with a 230 GHz power spectrum slope of 2.3+0.8-0.6 at 95% confidence. The light curve is de-correlated (white noise) on very long (month to year) times. In order to identify the transition time between red and white noise, we model the light curves as a stochastic damped random walk process. The models allow a quantitative estimate of this physical characteristic time-scale of 8-4+3 hours at 230 GHz at 95% confidence, with consistent results at 345 and 690 GHz. This corresponds to ~10 orbital times or ~1 inflow (viscous) time at R = 3 Rs, a typical radius producing the 230 GHz emission as measured by very long baseline interferometry and found in theoretical accretion flow and jet models. This time-scale is significantly shorter (longer) than those measured at radio (near-infrared, NIR) wavelengths, and is marginally inconsistent with the same variability mechanism operating in the submm and NIR for the expected t ~ R^3/2 scaling. It is crudely consistent with the analogous time-scale inferred in studies of quasar optical light curves after accounting for the difference in emission radius. We find evidence that the submm variability persists at least down to the ISCO, if not the event horizon. These results can be compared quantitatively with similar analyses at different wavebands to test for connections between the variability mechanisms, and with light curves from theoretical models of accreting black holes.

preprint2013arXiv

High-Angular-Resolution and High-Sensitivity Science Enabled by Beamformed ALMA

An international consortium is presently constructing a beamformer for the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile that will be available as a facility instrument. The beamformer will aggregate the entire collecting area of the array into a single, very large aperture. The extraordinary sensitivity of phased ALMA, combined with the extremely fine angular resolution available on baselines to the Northern Hemisphere, will enable transformational new very long baseline interferometry (VLBI) observations in Bands 6 and 7 (1.3 and 0.8 mm) and provide substantial improvements to existing VLBI arrays in Bands 1 and 3 (7 and 3 mm). The ALMA beamformer will have impact on a variety of scientific topics, including accretion and outflow processes around black holes in active galactic nuclei (AGN), tests of general relativity near black holes, jet launch and collimation from AGN and microquasars, pulsar and magnetar emission processes, the chemical history of the universe and the evolution of fundamental constants across cosmic time, maser science, and astrometry.

preprint2013arXiv

Transient jet formation and state transitions from large-scale magnetic reconnection in black hole accretion discs

Magnetically arrested accretion discs (MADs), where the magnetic pressure in the inner disc is dynamically important, provide an alternative mechanism for regulating accretion to what is commonly assumed in black hole systems. We show that a global magnetic field inversion in the MAD state can destroy the jet, significantly increase the accretion rate, and move the effective inner disc edge in to the marginally stable orbit. Reconnection of the MAD field in the inner radii launches a new type of transient outflow containing hot plasma generated by magnetic dissipation. This transient outflow can be as powerful as the steady magnetically-dominated Blandford-Znajek jet in the MAD state. The field inversion qualitatively describes many of the observational features associated with the high luminosity hard to soft state transition in black hole X-ray binaries: the jet line, the transient ballistic jet, and the drop in rms variability. These results demonstrate that the magnetic field configuration can influence the accretion state directly, and hence the magnetic field structure is an important second parameter in explaining observations of accreting black holes across the mass and luminosity scales.

preprint2012arXiv

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

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

preprint2012arXiv

Inhomogeneous accretion discs and the soft states of black hole X-ray binaries

Observations of black hole binaries (BHBs) have established a rich phenomenology of X-ray states. The soft states range from the low variability, accretion disc dominated thermal state (TD) to the higher variability, non-thermal steep power law state (SPL). The disc component in all states is typically modeled with standard thin disc accretion theory. However, this theory is inconsistent with optical/UV spectral, variability, and gravitational microlensing observations of active galactic nuclei (AGNs), the supermassive analogs of BHBs. An inhomogeneous disc (ID) model with large (~0.4 dex) temperature fluctuations in each radial annulus can qualitatively explain all of these AGN observations. The inhomogeneity may be a consequence of instabilities in radiation dominated discs, and therefore may be present in BHBs as well. We show that ID models can explain many features of the TD and SPL states of BHBs. The observed relationships between spectral hardness, disc fraction, and rms variability amplitude in BHBs are reproduced with temperature fluctuations similar to those inferred in AGNs, suggesting a unified picture of luminous accretion discs across orders of magnitude in black hole mass. This picture can be tested with spectral fitting of ID models, X-ray polarization observations, and radiation MHD simulations. If BHB accretion discs are indeed inhomogeneous, only the most disc dominated states (disc fraction > 0.95) can be used to robustly infer black hole spin using current continuum fitting methods.

preprint2012arXiv

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

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

preprint2012arXiv

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

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

preprint2012arXiv

Supernova Light Curves Powered by Fallback Accretion

Some fraction of the material ejected in a core collapse supernova explosion may remain bound to the compact remnant, and eventually turn around and fall back. We show that the late time (> days) power associated with the accretion of this "fallback" material may significantly affect the optical light curve, in some cases producing super-luminous or otherwise peculiar supernovae. We use spherically symmetric hydrodynamical models to estimate the accretion rate at late times for a range of progenitor masses and radii and explosion energies. The accretion rate onto the proto-neutron star or black hole decreases as Mdot ~ t^-5/3 at late times, but its normalization can be significantly enhanced at low explosion energies, in very massive stars, or if a strong reverse shock wave forms at the helium/hydrogen interface in the progenitor. If the resulting super-Eddington accretion drives an outflow which thermalizes in the outgoing ejecta, the supernova debris will be re-energized at a time when photons can diffuse out efficiently. The resulting light curves are different and more diverse than previous fallback supernova models which ignored the input of accretion power and produced short-lived, dim transients. The possible outcomes when fallback accretion power is significant include super-luminous (> 10^44 ergs / s) Type II events of both short and long durations, as well as luminous Type I events from compact stars that may have experienced significant mass loss. Accretion power may unbind the remaining infalling material, causing a sudden decrease in the brightness of some long duration Type II events. This scenario may be relevant for explaining some of the recently discovered classes of peculiar and rare supernovae.

preprint2012arXiv

The size of the jet launching region in M87

The supermassive black hole candidate at the center of M87 drives an ultra-relativistic jet visible on kiloparsec scales, and its large mass and relative proximity allow for event horizon scale imaging with very long baseline interferometry at millimeter wavelengths (mm-VLBI). Recently, relativistic magneto-hydrodynamic (MHD) simulations of black hole accretion flows have proven capable of launching magnetically-dominated jets. We construct time-dependent disc/jet models of the innermost portion of the M87 nucleus by performing relativistic radiative transfer calculations from one such simulation. We identify two types of models, jet-dominated or disc/jet, that can explain the spectral properties of M87, and use them to make predictions for current and future mm-VLBI observations. The Gaussian source size for the favored sky orientation and inclination from observations of the large-scale jet is 33-44 microarcseconds (~4-6 Schwarzschild radii) on current mm-VLBI telescopes, very similar to existing observations of Sgr A*. The black hole shadow, direct evidence of an event horizon, should be visible in future measurements using baselines between Hawaii and Mexico. Both models exhibit variability at millimeter wavelengths with factor of ~2 amplitudes on year timescales. For the low inclination of M87, the counter-jet dominates the event horizon scale millimeter wavelength emission from the jet-forming region.

preprint2012arXiv

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

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

preprint2011arXiv

Observational Signatures of Tilted Black Hole Accretion Disks from Simulations

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

preprint2010arXiv

Quasar Accretion Disks Are Strongly Inhomogeneous

Active galactic nuclei (AGN) have been observed to vary stochastically with 10-20 rms amplitudes over a range of optical wavelengths where the emission arises in an accretion disk. Since the accretion disk is unlikely to vary coherently, local fluctuations may be significantly larger than the global rms variability. We investigate toy models of quasar accretion disks consisting of a number of regions, n, whose temperatures vary independently with an amplitude of σ_T in dex. Models with large fluctuations (σ_T=0.35-0.50) in 100-1000 independently fluctuating zones for every factor of two in radius can explain the observed discrepancy between thin accretion disk sizes inferred from microlensing events and optical luminosity while matching the observed optical variability. For the same range of σ_T, inhomogeneous disk spectra provide excellent fits to the HST quasar composite without invoking global Compton scattering atmospheres to explain the high levels of observed UV emission. Simulated microlensing light curves for the Einstein cross from our time-varying toy models are well fit using a time-steady power-law temperature disk, and produce magnification light curves that are consistent with current microlensing observations. Deviations due to the inhomogeneous, time-dependent disk structure should occur above the 1% level in the light curves, detectable in future microlensing observations with millimag sensitivity.

preprint2010arXiv

The Submillimeter Bump in Sgr A* from Relativistic MHD Simulations

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

preprint2009arXiv

A Fast New Public Code for Computing Photon Orbits in a Kerr Spacetime

Relativistic radiative transfer problems require the calculation of photon trajectories in curved spacetime. We present a novel technique for rapid and accurate calculation of null geodesics in the Kerr metric. The equations of motion from the Hamilton-Jacobi equation are reduced directly to Carlson's elliptic integrals, simplifying algebraic manipulations and allowing all coordinates to be computed semi-analytically for the first time. We discuss the method, its implementation in a freely available FORTRAN code, and its application to toy problems from the literature.

preprint2009arXiv

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

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