Source author record

David R. Ballantyne

David R. Ballantyne 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

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

16 published item(s)

preprint2022arXiv

Burst-Disk Interaction in 4U 1636-536 as observed by NICER

We present the detection of 51 thermonuclear X-ray bursts observed from 4U 1636-536 by the Neutron Star Interior Composition Explorer (NICER) over the course of a three year monitoring campaign. We performed time resolved spectroscopy for 40 of these bursts and showed the existence of a strong soft excess in all the burst spectra. The excess emission can be characterized by the use of a scaling factor (f_a method) to the persistent emission of the source, which is attributed to the increased mass accretion rate on to the neutron star due to Poynting-Robertson drag. The soft excess emission can also be characterized by the use of a model taking into account the reflection of the burst emission off of the accretion disk. We also present time resolved spectral analysis of 5 X-ray bursts simultaneously observed by NICER and AstroSat, which confirm the main results with even greater precision. Finally, we present evidence for Compton cooling using 7 X-ray bursts observed contemporaneously with \nustar, by means of a correlated decrease in the hard X-ray lightcurve of 4U 1636-536 as the bursts start.

preprint2022arXiv

Massive Black Hole Binaries from the TNG50-3 Simulation: I. Coalescence and LISA Detection Rates

We evaluate the cosmological coalescence and detection rates for massive black hole (MBH) binaries targeted by the gravitational wave observatory Laser Interferometer Space Antenna (LISA). Our calculation starts with a population of gravitationally unbound MBH pairs, drawn from the TNG50-3 cosmological simulation, and follows their orbital evolution from kpc scales all the way to coalescence using a semi-analytic model developed in our previous work. We find that for a majority of MBH pairs that coalesce within a Hubble time dynamical friction is the most important mechanism that determines their coalescence rate. Our model predicts a MBH coalescence rate < 0.45/ yr and a LISA detection rate < 0.34/ yr. Most LISA detections should originate from 10^6 - 10^6.8 solar masses MBHs in gas-rich galaxies at redshifts 1.6 < z < 2.4, and have a characteristic signal to noise ratio SNR ~ 100. We however find a dramatic reduction in the coalescence and detection rates, as well as the average SNR, if the effects of radiative feedback from accreting MBHs are taken into account. In this case, the MBH coalescence rate is reduced by 78% (to < 0.1/ yr), and the LISA detection rate is reduced by 94% (to 0.02/ yr), whereas the average SNR is ~ 10. We emphasize that our model provides a lower limit on the LISA detection rate, consistent with other works in the literature that draw their MBH pairs from cosmological simulations.

preprint2022arXiv

Massive Black Hole Binaries from the TNG50-3 Simulation: II. Using Dual AGNs to Predict the Rate of Black Hole Mergers

Dual active galaxy nuclei (dAGNs) trace the population of post-merger galaxies and are the precursors to massive black hole (MBH) mergers, an important source of gravitational waves that may be observed by LISA. In Paper I of this series, we used the population of nearly 2000 galaxy mergers predicted by the TNG50-3 simulation to seed semi-analytic models of the orbital evolution and coalescence of MBH pairs with initial separations of about 1 kpc. Here, we calculate the dAGN luminosities and separation of these pairs as they evolve in post-merger galaxies, and show how the coalescence fraction of dAGNs changes with redshift. We find that because of the several Gyr long dynamical friction timescale for orbital evolution, the fraction of dAGNs that eventually end in a MBH merger grows with redshift and does not pass 50% until a redshift of 1. However, dAGNs in galaxies with bulge masses >10^10 solar masses, or comprised of near-equal mass MBHs, evolve more quickly and have higher than average coalescence fractions. At any redshift, dAGNs observed with small separations (> 0.7 kpc) have a higher probability of merging in a Hubble time than more widely separated systems. As found in Paper I, radiation feedback effects can significantly reduce the number of MBH mergers, and this could be manifested as a larger than expected number of widely separated dAGNs. We present a method to estimate the MBH coalescence rate as well as the potential LISA detection rate given a survey of dAGNs. Comparing these rates to the eventual LISA measurements will help determine the efficiency of dynamical friction in post-merger galaxies.

preprint2020arXiv

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

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

preprint2016arXiv

The Geometry of the Infrared and X-ray Obscurer in a Dusty Hyperluminous Quasar

We study the geometry of the AGN obscurer in IRAS 09104+4109, an IR-luminous, radio-intermediate FR-I source at $z=0.442$, using infrared data from Spitzer and Herschel, X-ray data from Nustar, Swift, Suzaku, and Chandra, and an optical spectrum from Palomar. The infrared data imply a total rest-frame 1-1000$μ$m luminosity of $5.5\times10^{46}$erg s$^{-1}$ and require both an AGN torus and starburst model. The AGN torus has an anisotropy-corrected IR luminosity of $4.9\times10^{46}$erg s$^{-1}$, and a viewing angle and half opening angle both of approximately $36$ degrees from pole-on. The starburst has a star formation rate of $(110\pm34)$M$_{\odot}$ yr$^{-1}$ and an age of $<50$Myr. These results are consistent with two epochs of luminous activity in 09104: one approximately $150$Myr ago, and one ongoing. The X-ray data suggest a photon index of $Γ\simeq 1.8$ and a line-of-sight column of $N_{\rm H} \simeq 5\times10^{23}$cm$^{-2}$. This argues against a reflection-dominated hard X-ray spectrum, which would have implied a much higher $N_{\rm H}$ and luminosity. The X-ray and infrared data are consistent with a bolometric AGN luminosity of $L_{\rm bol}\sim(0.5-2.5)\times10^{47}$erg s$^{-1}$. The X-ray and infrared data are further consistent with coaligned AGN obscurers in which the line of sight "skims" the torus. This is also consistent with the optical spectra, which show both coronal iron lines and broad lines in polarized but not direct light. Combining constraints from the X-ray, optical, and infrared data suggests that the AGN obscurer is within a vertical height of $20$pc, and a radius of $125$pc, of the nucleus.

preprint2015arXiv

NuSTAR observations of the powerful radio-galaxy Cygnus A

We present NuSTAR observations of the powerful radio galaxy Cygnus A, focusing on the central absorbed active galactic nucleus (AGN). Cygnus A is embedded in a cool-core galaxy cluster, and hence we also examine archival XMM-Newton data to facilitate the decomposition of the spectrum into the AGN and intracluster medium (ICM) components. NuSTAR gives a source-dominated spectrum of the AGN out to >70keV. In gross terms, the NuSTAR spectrum of the AGN has the form of a power law (Gamma~1.6-1.7) absorbed by a neutral column density of N_H~1.6x10^23 cm^-2. However, we also detect curvature in the hard (>10keV) spectrum resulting from reflection by Compton-thick matter out of our line-of-sight to the X-ray source. Compton reflection, possibly from the outer accretion disk or obscuring torus, is required even permitting a high-energy cutoff in the continuum source; the limit on the cutoff energy is E_cut>111keV (90% confidence). Interestingly, the absorbed power-law plus reflection model leaves residuals suggesting the absorption/emission from a fast (15,000-26,000km/s), high column-density (N_W>3x10^23 cm^-2), highly ionized (xi~2,500 erg cm/s) wind. A second, even faster ionized wind component is also suggested by these data. We show that the ionized wind likely carries a significant mass and momentum flux, and may carry sufficient kinetic energy to exercise feedback on the host galaxy. If confirmed, the simultaneous presence of a strong wind and powerful jets in Cygnus A demonstrates that feedback from radio-jets and sub-relativistic winds are not mutually exclusive phases of AGN activity but can occur simultaneously.

preprint2015arXiv

On the equivalent width of the Fe K$α$ line produced by a dusty absorber in active galactic nuclei

Obscured AGNs provide an opportunity to study the material surrounding the central engine. Geometric and physical constraints on the absorber can be deduced from the reprocessed AGN emission. In particular, the obscuring gas may reprocess the nuclear X-ray emission producing a narrow Fe K$α$ line and a Compton reflection hump. In recent years, models of the X-ray reflection from an obscuring torus have been computed; however, although the reflecting gas may be dusty, the models do not yet take into account the effects of dust on the predicted spectrum. We study this problem by analyzing two sets of models, with and without the presence of dust, using the one dimensional photo-ionization code Cloudy. The calculations are performed for a range of column densities ($22 <{\rm log}[N_H(\rm cm^{-2})]< 24.5$ ) and hydrogen densities ( $6 <{\rm log}[n_H(\rm cm^{-3})]< 8$). The calculations show the presence of dust can enhance the Fe K$α$ equivalent width (EW) in the reflected spectrum by factors up to $\approx$ 8 for Compton thick (CT) gas and a typical ISM grain size distribution. The enhancement in EW with respect to the reflection continuum is due to the reduction in the reflected continuum intensity caused by the anisotropic scattering behaviour of dust grains. This effect will be most relevant for reflection from distant, predominately neutral gas, and is a possible explanation for AGNs which show a strong Fe K$α$ EW and a relatively weak reflection continuum. Our results show it is an important to take into account dust while modeling the X-ray reflection spectrum, and that inferring a CT column density from an observed Fe K$α$ EW may not always be valid. Multi-dimensional models are needed to fully explore the magnitude of the effect.

preprint2013arXiv

On the contribution of active galactic nuclei to reionization (Research Note)

The electron scattering optical depth constraints on reionization suggest that there may be other sources that contribute to the ionization of hydrogen aside from observable star forming galaxies. Often the calculated value of the electron scattering optical depth, τ_{es}, falls below the measurements derived from observations of the CMB or an assumption about non-observable sources must be made in order to reach agreement. Here, we calculate the hydrogen ionization fraction as a function of redshift and the electron scattering optical depth from both galaxies and active galactic nuclei (AGN) factoring in the secondary collisional ionizations from the AGN X-ray emission. In this paper we use the most current determination of the evolving hard X-ray luminosity function and extrapolate its evolution beyond z = 6. The AGN spectral energy distributions (SEDs) include both UV and X-ray ionizing photons. To search for the largest possible effect, all AGN are assumed to have λ_{Edd} = 1.0 and be completely unobscured. The results show that AGNs produce a perturbative effect on the reionization of hydrogen and remains in agreement with current constraints. Our calculations find the epoch of reionization still ends at z \approx 6 and only increases the electron scattering optical depth by ~2.3% under the most optimal conditions. This can only be moderately increased by assuming a constant black hole mass of M_{BH} = 10^(5) solar masses. As a result, we conclude that there is a need for other sources beyond observable galaxies and AGNs that contribute to the reionization of hydrogen at z > 6.

preprint2013arXiv

The Nuclear Spectroscopic Telescope Array (NuSTAR) Mission

The Nuclear Spectroscopic Telescope Array (NuSTAR) mission, launched on 13 June 2012, is the first focusing high-energy X-ray telescope in orbit. NuSTAR operates in the band from 3 -- 79 keV, extending the sensitivity of focusing far beyond the ~10 keV high-energy cutoff achieved by all previous X-ray satellites. The inherently low-background associated with concentrating the X-ray light enables NuSTAR to probe the hard X-ray sky with a more than one-hundred-fold improvement in sensitivity over the collimated or coded-mask instruments that have operated in this bandpass. Using its unprecedented combination of sensitivity, spatial and spectral resolution, NuSTAR will pursue five primary scientific objectives, and will also undertake a broad program of targeted observations. The observatory consists of two co-aligned grazing-incidence X-ray telescopes pointed at celestial targets by a three-axis stabilized spacecraft. Deployed into a 600 km, near-circular, 6degree inclination orbit, the Observatory has now completed commissioning, and is performing consistent with pre-launch expectations. NuSTAR is now executing its primary science mission, and with an expected orbit lifetime of ten years, we anticipate proposing a guest investigator program, to begin in Fall 2014.

preprint2012arXiv

A Tale of Two Populations: The Contribution of Merger and Secular Processes to the Evolution of Active Galactic Nuclei

Due to the co-evolution of supermassive black holes and their host galaxies, understanding the mechanisms that trigger active galactic nuclei (AGN) are imperative to understanding galaxy evolution and the formation of massive galaxies. It is observationally difficult to determine the trigger of a given AGN due to the difference between the AGN lifetime and triggering timescales. Here, we utilize AGN population synthesis modeling to determine the importance of different AGN triggering mechanisms. An AGN population model is computed by combining an observationally motivated AGN triggering rate and a theoretical AGN light curve. The free parameters of the AGN light curve are constrained by minimizing a χsquared test with respect to the observed AGN hard X-ray luminosity function. The observed black hole space density, AGN number counts, and X-ray background spectrum are also considered as observational constraints. It is found that major mergers are not able to account for the entire AGN population. Therefore, non-merger processes, such as secular mechanisms, must also trigger AGN. Indeed, non-merger processes are the dominant AGN triggering mechanism at z \lesssim 1--1.5. Furthermore, the shape and evolution of the black hole mass function of AGN triggered by major mergers is intrinsically different from the shape and evolution of the black hole mass function of AGN triggered by secular processes.

preprint2012arXiv

CO Line Emission from Compact Nuclear Starburst Disks Around Active Galactic Nuclei

There is substantial evidence for a connection between star formation in the nuclear region of a galaxy and growth of the central supermassive black hole. Furthermore, starburst activity in the region around an active galactic nucleus (AGN) may provide the obscuration required by the unified model of AGN. Molecular line emission is one of the best observational avenues to detect and characterize dense, star-forming gas in galactic nuclei over a range of redshift. This paper presents predictions for the carbon monoxide (CO) line features from models of nuclear starburst disks around AGN. These small scale ($\la 100$ pc), dense and hot starbursts have CO luminosities similar to scaled-down ultra-luminous infrared galaxies and quasar host galaxies. Nuclear starburst disks that exhibit a pc-scale starburst and could potentially act as the obscuring torus show more efficient CO excitation and higher brightness temperature ratios than those without such a compact starburst. In addition, the compact starburst models predict strong absorption when $J_{\mathrm{Upper}} \ga 10$, a unique observational signature of these objects. These findings allow for the possibility that CO SLEDs could be used to determine if starburst disks are responsible for the obscuration in $z \la 1$ AGN. Directly isolating the nuclear CO line emission of such compact regions around AGN from galactic-scale emission will require high resolution imaging or selecting AGN host galaxies with weak galactic-scale star formation. Stacking individual CO SLEDs will also be useful in detecting the predicted high-$J$ features.

preprint2012arXiv

The Merger-Triggered Active Galactic Nuclei Contribution to the Ultraluminous Infrared Galaxy Population

It has long been thought that there is a connection between ultraluminous infrared galaxies (ULIRGs), quasars, and major mergers. Indeed, simulations show that major mergers are capable of triggering massive starbursts and quasars. However, observations by the Herschel Space Observatory suggest that, at least at high redshift, there may not always be a simple causal connection between ULIRGs and mergers. Here, we combine an evolving merger-triggered AGN luminosity function with a merger-triggered starburst model to calculate the maximum contribution of major mergers to the ULIRG population. We find that major mergers can account for the entire local population of ULIRGs hosting AGN and ~25% of the total local ULIRG luminosity density. By z ~ 1, major mergers can no longer account for the luminosity density of ULIRGs hosting AGN and contribute \lesssim 12% of the total ULIRG luminosity density. This drop is likely due to high redshift galaxies being more gas rich and therefore able to achieve high star formation rates through secular evolution. Additionally, we find that major mergers can account for the local population of warm ULIRGs. This suggests that selecting high redshift warm ULIRGs will allow for the identification of high redshift merger-triggered ULIRGs. As major mergers are likely to trigger very highly obscured AGN, a significant fraction of the high redshift warm ULIRG population may host Compton thick AGN.

preprint2011arXiv

Are Active Galactic Nuclei the Solution to the Excess Cosmic Radio Background at 1.4 GHz?

Recently the ARCADE 2 experiment measured the cosmic radio background (CRB) and found the brightness temperature of the CRB at 1.4 GHz to be ~480 mK. Integrating the flux density from the observed 1.4 GHz radio source count produces a brightness temperature of ~100 mK---less than a quarter of the observed CRB at 1.4 GHz. Radio quiet AGN are a large fraction of the 1.4 GHz uJy sources and typically host significant star formation. Thus, it is possible that AGN and host star formation could be responsible for some fraction of the excess CRB at 1.4 GHz. Here, an X-ray background population synthesis model is used in conjunction with empirical radio to X-ray luminosity ratios to calculate the AGN contribution to the CRB at 1.4 GHz including the emission from host star formation. It is found that AGN and host star formation contribute <~9% of the CRB at 1.4 GHz. When all known 1.4 GHz radio source classes are considered, <~60% of the CRB at 1.4 GHz is accounted for; therefore, it is necessary that either known radio sources evolve significantly at flux densities below current survey sensitivity limits or a new population of low flux density radio sources exist.

preprint2011arXiv

Properties and Expected Number Counts of Active Galactic Nuclei and their Hosts in the Far Infrared

Telescopes like Herschel and the Atacama Large Millimeter/submillimeter Array (ALMA) are creating new opportunities to study sources in the far infrared (FIR), a wavelength region dominated by cold dust emission. Probing cold dust in active galaxies allows for study of the star formation history of active galactic nuclei (AGN) hosts. The FIR is also an important spectral region for observing AGN which are heavily enshrouded by dust, such as Compton thick AGN. By using information from deep X-ray surveys and cosmic X-ray background synthesis models, we compute Cloudy photoionization simulations which are used to predict the spectral energy distribution (SED) of AGN in the FIR. Expected differential number counts of AGN and their host galaxies are calculated in the Herschel bands. The expected contribution of AGN and their hosts to the cosmic infrared background (CIRB) and the infrared luminosity density are also computed. Multiple star formation scenarios are investigated using a modified blackbody star formation SED. It is found that FIR observations at ~500 um are an excellent tool in determining the star formation history of AGN hosts. Additionally, the AGN contribution to the CIRB can be used to determine whether star formation in AGN hosts evolves differently than in normal galaxies. The contribution of Compton thick AGN to the bright end differential number counts and to the bright source infrared luminosity density is a good test of AGN evolution models where quasars are triggered by major mergers.

preprint2011arXiv

The Young, the Old, and the Dusty: Stellar Populations of AGN Hosts

Studying the average properties of active galactic nuclei (AGN) host stellar populations is an important step in understanding the role of AGN in galaxy evolution and the processes which trigger and fuel AGN activity. Here we calculate model spectral energy distributions (SEDs) that include emission from the AGN, the host galaxy stellar population, and dust enshrouded star formation. Using the framework of cosmic X-ray background population synthesis modeling, the model AGN hosts are constrained using optical (B band) and near infrared (J band, 3.6 um, 5.7 um, 8.0 um, and 24 um) luminosity functions and number counts. It is found that at z < 1, type 1 and type 2 AGN hosts have similar stellar populations, in agreement with the orientation based unified model and indicative of secular evolution. At z > 1, type 2 AGN hosts are intrinsically different from type 1 AGN hosts, suggesting that the simple orientation based unified model does not hold at z > 1. Also, it is found that if Compton thick (CT) AGN evolve like less obscured type 2 AGN, then, on average, CT AGN hosts are similar to type 2 AGN hosts; however, if CT obscuration is connected to an evolutionary stage of black hole growth, then CT AGN hosts will also be in specific evolutionary stages. Multi-wavelength selection criteria of CT AGN are discussed.

preprint1999arXiv

Emission-line Helium Abundances in Highly Obscured Nebulae

This paper outlines a way to determine the ICF using only infrared data. We identify four line pairs, [NeIII] 36$\micron$/[NeII] 12.8$\micron$, [NeIII]~15.6$\micron$ /[NeII] 12.8$\micron$, [ArIII] 9$\micron$/[ArII] 6.9$\micron$, and [ArIII] 21$\micron$/[ArII] 6.9$\micron$, that are sensitive to the He ICF. This happens because the ions cover a wide range of ionization, the line pairs are not sensitive to electron temperature, they have similar critical densities, and are formed within the He$^+$/H$^+$ region of the nebula. We compute a very wide range of photoionization models appropriate for galactic HII regions. The models cover a wide range of densities, ionization parameters, stellar temperatures, and use continua from four very different stellar atmospheres. The results show that each line pair has a critical intensity ratio above which the He ICF is always small. Below these values the ICF depends very strongly on details of the models for three of the ratios, and so other information would be needed to determine the helium abundance. The [Ar III] 9$\micron$/[ArII] 6.9$\micron$ ratio can indicate the ICF directly due to the near exact match in the critical densities of the two lines. Finally, continua predicted by the latest generation of stellar atmospheres are sufficiently hard that they routinely produce significantly negative ICFs.