Harnessing the sheer discovery potential of gravitational wave astronomy will require bold, deliberate, and sustained efforts to train and develop the requisite workforce. The next decade requires a strategic plan to build -- from the ground up -- a robust, open, and well-connected gravitational wave astronomy community with deep participation from traditional astronomers, physicists, data scientists, and instrumentalists. This basic infrastructure is sorely needed as an enabling foundation for research. We outline a set of recommendations for funding agencies, universities, and professional societies to help build a thriving, diverse, and inclusive new field.

One of the central goals of LISA is the detection of gravitational waves from the merger of supermassive black holes. Contrary to stellar-mass black hole mergers, such events are expected to be rich X-ray sources due to the accretion of material from the circumbinary disks onto the black holes. The orbital dynamics before merger is also expected to modulate the resulting X-ray emission via Doppler boosting in a quasi-periodic way, and in a simple phase relation with the gravitational wave from the inspiral of the black holes. Detecting the X-ray source would enable a precise and early localization of the binary, thus allowing many telescopes to observe the very moment of the merger. Although identifying the correct X-ray source in the relatively large LISA sky localization will be challenging due to the presence of many confounding point sources, the quasi-periodic modulation may aid in the identification. We explore the practical feasibility of such idea. We simulate populations of merging supermassive black holes, their detection with LISA and their X-ray lightcurves using a simple model. Taking the parameters of the X-ray Telescope on the proposed NASA Transient Astrophysics Probe, we then design and simulate an observation campaign that searches for the modulated X-ray source while LISA is still observing the inspiral of the black holes. Assuming a fiducial LISA detection rate of $10$ mergers per year at redshift closer than $3.5$, we expect a few detections of modulated X-ray counterparts over the nominal duration of the LISA mission.

A natural consequence of the galaxy formation paradigm is the existence of supermassive black hole (SMBH) binaries. Gravitational perturbations from a far-away SMBH companion can induce high orbital eccentricities on dark matter particles orbiting the primary SMBH via the eccentric Kozai-Lidov mechanism. This process yields an influx of dark matter particles into the primary SMBH ergosphere, where test particles linger for long timescales. This influx results in high self-gravitating densities, forming a dark matter clump that is extremely close to the SMBH. In such a situation, the gravitational wave emission between the dark matter clump and the SMBH is potentially detectable by LISA. If dark matter self-annihilates, the high densities of the clump will result in a unique co-detection of gravitational wave emission and high energy electromagnetic signatures.

Shortly after the discovery of the Kerr metric in 1963, it was realized that a region existed outside of the black hole's event horizon where no time-like observer could remain stationary. In 1969, Roger Penrose showed that particles within this ergosphere region could possess negative energy, as measured by an observer at infinity. When captured by the horizon, these negative energy particles essentially extract mass and angular momentum from the black hole. While the decay of a single particle within the ergosphere is not a particularly efficient means of energy extraction, the "collision" of multiple particles can reach arbitrarily high center-of-mass energy in the limit of extremal black hole spin. The resulting particles can escape with high efficiency, potentially serving as a probe of high-energy particle physics as well as general relativity. In this paper, we briefly review the history of the field and highlight a specific astrophysical application of the collisional Penrose process: the potential to enhance annihilation of dark matter particles in the vicinity of a supermassive black hole.

In the science fiction film $Interstellar$, a band of intrepid astronauts sets out to explore a system of planets orbiting a supermassive black hole, searching for a world that may be conducive to hosting human life. While the film legitimately boasts a relatively high level of scientific accuracy, it is still restricted by Hollywood sensitivities and limitations. In this paper, we discuss a number of additional astrophysical effects that may be important in determining the (un)inhabitable environment of a planet orbiting close to a giant, accreting black hole. Foremost among these effects is the blueshift and beaming of incident radiation on the planet, due to the time dilation of an observer orbiting very close to the black hole. This results in high-energy flux incoming from surrounding stars and background radiation, with significant implications for habitability.

The first terrestrial gravitational wave interferometers have dramatically underscored the scientific value of observing the Universe through an entirely different window, and of folding this new channel of information with traditional astronomical data for a multimessenger view. The Laser Interferometer Space Antenna (LISA) will broaden the reach of gravitational wave astronomy by conducting the first survey of the millihertz gravitational wave sky, detecting tens of thousands of individual astrophysical sources ranging from white-dwarf binaries in our own galaxy to mergers of massive black holes at redshifts extending beyond the epoch of reionization. These observations will inform - and transform - our understanding of the end state of stellar evolution, massive black hole birth, and the co-evolution of galaxies and black holes through cosmic time. LISA also has the potential to detect gravitational wave emission from elusive astrophysical sources such as intermediate-mass black holes as well as exotic cosmological sources such as inflationary fields and cosmic string cusps.

The Laser Interferometer Space Antenna (LISA) will open three decades of gravitational wave (GW) spectrum between 0.1 and 100 mHz, the mHz band. This band is expected to be the richest part of the GW spectrum, in types of sources, numbers of sources, signal-to-noise ratios and discovery potential. When LISA opens the low-frequency window of the gravitational wave spectrum, around 2034, the surge of gravitational-wave astronomy will strongly compel a subsequent mission to further explore the frequency bands of the GW spectrum that can only be accessed from space. The 2020s is the time to start developing technology and studying mission concepts for a large-scale mission to be launched in the 2040s. The mission concept would then be proposed to Astro2030. Only space based missions can access the GW spectrum between 10 nHz and 1 Hz because of the Earths seismic noise. This white paper surveys the science in this band and mission concepts that could accomplish that science. The proposed small scale activity is a technology development program that would support a range of concepts and a mission concept study to choose a specific mission concept for Astro2030. In this white paper, we will refer to a generic GW mission beyond LISA as bLISA.

Experience suggests that structural issues in how institutional Astrophysics approaches data-driven science and the development of discovery technology may be hampering the community's ability to respond effectively to a rapidly changing environment in which increasingly complex, heterogeneous datasets are challenging our existing information infrastructure and traditional approaches to analysis. We stand at the confluence of a new epoch of multimessenger science, remote co-location of data and processing power and new observing strategies based on miniaturized spacecraft. Significant effort will be required by the community to adapt to this rapidly evolving range of possible discovery moduses. In the suggested creation of a new Astrophysics element, Advanced Astrophysics Discovery Technology, we offer an affirmative solution that places the visibility of discovery technologies at a level that we suggest is fully commensurate with their importance to the future of the field.

This white paper highlights compact object and fundamental physics science opportunities afforded by high-throughput broadband (0.1-60 keV) X-ray polarization observations. X-ray polarimetry gives new observables with geometric information about stellar remnants which are many orders of magnitude too small for direct imaging. The X-ray polarimetric data also reveal details about the emission mechanisms and the structure of the magnetic fields in and around the most extreme objects in the Universe. Whereas the Imaging X-ray Polarimetry Explorer (IXPE) to be launched in 2021 will obtain first results for bright objects, a follow-up mission could be one order of magnitude more sensitive and would be able to use a broader bandpass to perform physics type experiments for representative samples of sources.

Coalescing, massive black-hole (MBH) binaries are the most powerful sources of gravitational waves (GWs) in the Universe, which makes MBH science a prime focus for ongoing and upcoming GW observatories. The Laser Interferometer Space Antenna (LISA) -- a gigameter scale space-based GW observatory -- will grant us access to an immense cosmological volume, revealing MBHs merging when the first cosmic structures assembled in the Dark Ages. LISA will unveil the yet unknown origin of the first quasars, and detect the teeming population of MBHs of $10^4 - 10^7$ solar masses. forming within protogalactic halos. The Pulsar Timing Array, a galactic-scale GW survey, can access the largest MBHs the Universe, detecting the cosmic GW foreground from inspiraling MBH binaries of about 10^9 solar masses. LISA can measure MBH spins and masses with precision far exceeding that from electromagnetic (EM) probes, and together, both GW observatories will provide the first full census of binary MBHs, and their orbital dynamics, across cosmic time. Detecting the loud gravitational signal of these MBH binaries will also trigger alerts for EM counterpart searches, from decades (PTAs) to hours (LISA) prior to the final merger. By witnessing both the GW and EM signals of MBH mergers, precious information will be gathered about the rich and complex environment in the aftermath of a galaxy collision. The unique GW characterization of MBHs will shed light on the deep link between MBHs of $10^4-10^{10}$ solar masses and the grand design of galaxy assembly, as well as on the complex dynamics that drive MBHs to coalescence.

The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditions favorable for the origin, development and sustainment of life as we know it. This requires the understanding of global (astrospheric) and local (atmospheric, surface and internal) environments of exoplanets in the framework of the physical processes of the interaction between evolving planet-hosting stars along with exoplanetary evolution over geological timescales, and the resulting impact on climate and habitability of exoplanets. Feedbacks between astrophysical, physico-chemical atmospheric and geological processes can only be understood through interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary, Earth sciences, astrobiology, and the origin of life communities. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets and potential exomoons around them may significantly modify the extent and the location of the habitable zone and provide new directions for searching for signatures of life. Thus, characterization of stellar ionizing outputs becomes an important task for further understanding the extent of habitability in the universe. The goal of this white paper is to identify and describe promising key research goals to aid the theoretical characterization and observational detection of ionizing radiation from quiescent and flaring upper atmospheres of planet hosts as well as properties of stellar coronal mass ejections and stellar energetic particle events.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

LISA will open the mHz band of gravitational waves (GWs) to the astronomy community. The strong gravity which powers the variety of GW sources in this band is also crucial in a number of important astrophysical processes at the current frontiers of astronomy. These range from the beginning of structure formation in the early universe, through the origin and cosmic evolution of massive black holes in concert with their galactic environments, to the evolution of stellar remnant binaries in the Milky Way and in nearby galaxies. These processes and their associated populations also drive current and future observations across the electromagnetic (EM) spectrum. We review opportunities for science breakthroughs, involving either direct coincident EM+GW observations, or indirect multimessenger studies. We argue that for the US community to fully capitalize on the opportunities from the LISA mission, the US efforts should be accompanied by a coordinated and sustained program of multi-disciplinary science investment, following the GW data through to its impact on broad areas of astrophysics. Support for LISA-related multimessenger observers and theorists should be sized appropriately for a flagship observatory and may be coordinated through a dedicated mHz GW research center.

The LIGO/Virgo gravitational-wave (GW) interferometers have to-date detected ten merging black hole (BH) binaries, some with masses considerably larger than had been anticipated. Stellar-mass BH binaries at the high end of the observed mass range (with "chirp mass" ${\cal M} \gtrsim 25 M_{\odot}$) should be detectable by a space-based GW observatory years before those binaries become visible to ground-based GW detectors. This white paper discusses some of the synergies that result when the same binaries are observed by instruments in space and on the ground. We consider intermediate-mass black hole binaries (with total mass $M \sim 10^2 -10^4 M_{\odot}$) as well as stellar-mass black hole binaries. We illustrate how combining space-based and ground-based data sets can break degeneracies and thereby improve our understanding of the binary's physical parameters. While early work focused on how space-based observatories can forecast precisely when some mergers will be observed on the ground, the reverse is also important: ground-based detections will allow us to "dig deeper" into archived, space-based data to confidently identify black hole inspirals whose signal-to-noise ratios were originally sub-threshold, increasing the number of binaries observed in both bands by a factor of $\sim 4 - 7$.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

Low-frequency gravitational-wave astronomy can perform precision tests of general relativity and probe fundamental physics in a regime previously inaccessible. A space-based detector will be a formidable tool to explore gravity's role in the cosmos, potentially telling us if and where Einstein's theory fails and providing clues about some of the greatest mysteries in physics and astronomy, such as dark matter and the origin of the Universe.

We describe results from a new technique for the prediction of complete, self-consistent X-ray spectra from three-dimensional General Relativistic magnetohydrodynamic (GRMHD) simulations of black hole accretion flows. Density and cooling rate data from a HARM3D GRMHD simulation are processed by both an improved version of the Monte Carlo radiation transport code PANDURATA (in the corona) and the Feautrier solver PTRANSX (in the disk), with XSTAR subroutines. The codes are run in a sequential but iterative fashion to achieve globally energy-conserving and self-consistent radiation fields, temperature maps, and photoionization equilibria. The output is the X-ray spectrum as seen by a distant observer. For the example cases we consider here---a non-rotating $10 M_\odot$ black hole with solar abundances, accreting at 0.01, 0.03, 0.1, or 0.3 Eddington---we find spectra resembling actual observations of stellar-mass black holes in the soft or steep power-law state: broad thermal peaks (at 1-3 keV), steep power-laws extending to high energy ($\Gamma$ = 2.7-4.5), and prominent, asymmetric Fe K$\alpha$ emission lines with equivalent widths in the range 40-400 eV (larger EW at lower accretion rates). By starting with simulation data, we obviate the need for parameterized descriptions of the accretion flow geometry---no a priori specification of the corona's shape or flux, or the disk temperature or density, etc., are needed. Instead, we apply the relevant physical principles to simulation output using appropriate numerical techniques; this procedure allows us to calculate inclination-dependent spectra after choosing only a small number of physically meaningful parameters: black hole mass and spin, accretion rate, and elemental abundances.

It is currently unknown how common life is on exoplanets, or how long planets can remain viable for life. To date, we have a superficial notion of habitability, a necessary first step, but so far lacking an understanding of the detailed interaction between stars and planets over geological timescales, dynamical evolution of planetary systems, and atmospheric evolution on planets in other systems. A planet mass, net insolation, and atmospheric composition alone are insufficient to determine the probability that life on a planet could arise or be detected. The latter set of planetary considerations, among others, underpin the concept of the habitable zone (HZ), defined as the circumstellar region where standing bodies of liquid water could be supported on the surface of a rocky planet. However, stars within the same spectral class are often treated in the same way in HZ studies, without any regard for variations in activity among individual stars. Such formulations ignore differences in how nonthermal emission and magnetic energy of transient events in different stars affect the ability of an exoplanet to retain its atmosphere.In the last few years there has been a growing appreciation that the atmospheric chemistry, and even retention of an atmosphere in many cases, depends critically on the high-energy radiation and particle environments around these stars. Indeed, recent studies have shown stellar activity and the extreme space weather, such as that created by the frequent flares and coronal mass ejections (CMEs) from the active stars and young Sun, may have profoundly affected the chemistry and climate and thus habitability of the early Earth and terrestrial type exoplanets. The goal of this white paper is to identify and describe promising key research goals to aid the field of the exoplanetary habitability for the next 20 years.

Binary black hole (BBH) mergers provide a prime source for current and future interferometric GW observatories. Massive BBH mergers may often take place in plasma-rich environments, leading to the exciting possibility of a concurrent electromagnetic (EM) signal observable by traditional astronomical facilities. However, many critical questions about the generation of such counterparts remain unanswered. We explore mechanisms that may drive EM counterparts with magnetohydrodynamic simulations treating a range of scenarios involving equal-mass black-hole binaries immersed in an initially homogeneous fluid with uniform, orbitally aligned magnetic fields. We find that the time development of Poynting luminosity, which may drive jet-like emissions, is relatively insensitive to aspects of the initial configuration. In particular, over a significant range of initial values, the central magnetic field strength is effectively regulated by the gas flow to yield a Poynting luminosity of $10^{45}-10^{46} \rho_{-13} M_8^2 \, {\rm erg}\,{\rm s}^{-1}$, with BBH mass scaled to $M_8 \equiv M/(10^8 M_{\odot})$ and ambient density $\rho_{-13} \equiv \rho/(10^{-13} \, {\rm g} \, {\rm cm}^{-3})$. We also calculate the direct plasma synchrotron emissions processed through geodesic ray-tracing. Despite lensing effects and dynamics, we find the observed synchrotron flux varies little leading up to merger.

We calculate the electromagnetic signal of a gamma-ray flare coming from the surface of a neutron star shortly before merger with a black hole companion. Using a new version of the Monte Carlo radiation transport code Pandurata that incorporates dynamic spacetimes, we integrate photon geodesics from the neutron star surface until they reach a distant observer or are captured by the black hole. The gamma-ray light curve is modulated by a number of relativistic effects, including Doppler beaming and gravitational lensing. Because the photons originate from the inspiraling neutron star, the light curve closely resembles the corresponding gravitational waveform: a chirp signal characterized by a steadily increasing frequency and amplitude. We propose to search for these electromagnetic chirps using matched filtering algorithms similar to those used in LIGO data analysis.

We present first results from a new technique for the prediction of Fe K$\alpha$ profiles directly from general relativistic magnetohydrodynamic (GRMHD) simulations. Data from a GRMHD simulation are processed by a Monte Carlo global radiation transport code, which determines the X-ray flux irradiating the disk surface and the coronal electron temperature self-consistently. With that irradiating flux and the disk's density structure drawn from the simulation, we determine the reprocessed Fe K$\alpha$ emission from photoionization equilibrium and solution of the radiation transfer equation. We produce maps of the surface brightness of Fe K$\alpha$ emission over the disk surface, which---for our example of a $10 M_\odot$, Schwarzschild black hole accreting at $1\%$ the Eddington value---rises steeply one gravitational radius outside the radius of the innermost stable circular orbit and then falls $\propto r^{-2}$ at larger radii. We explain these features of the Fe K$\alpha$ radial surface brightness profile as consequences of the disk's ionization structure and an extended coronal geometry, respectively. We also present the corresponding Fe K$\alpha$ line profiles as would be seen by distant observers at several inclinations. Both the shapes of the line profiles and the equivalent widths of our predicted K$\alpha$ lines are qualitatively similar to those typically observed from accreting black holes. Most importantly, this work represents a direct link between theory and observation: in a fully self-consistent way, we produce observable results---iron fluorescence line profiles---from the theory of black hole accretion with almost no phenomenological assumptions.

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

We present a new upper limit on the energy that may be extracted from a Kerr black hole by means of particle collisions in the ergosphere (i.e., the "collisional Penrose process"). Earlier work on this subject has focused largely on particles with critical values of angular momentum falling into an extremal Kerr black hole from infinity and colliding just outside the horizon. While these collisions are able to reach arbitrarily high center-of-mass energies, it is very difficult for the reaction products to escape back to infinity, effectively limiting the peak efficiency of such a process to roughly $130\%$. When we allow one of the initial particles to have impact parameter $b > 2M$, and thus not get captured by the horizon, it is able to collide along outgoing trajectories, greatly increasing the chance that the products can escape. For equal-mass particles annihilating to photons, we find a greatly increased peak energy of $E_{\rm out} \approx 6\times E_{\rm in}$. For Compton scattering, the efficiency can go even higher, with $E_{\rm out} \approx 14\times E_{\rm in}$, and for repeated scattering events, photons can both be produced {\it and} escape to infinity with Planck-scale energies.

We use a Monte Carlo code to calculate the geodesic orbits of test particles around Kerr black holes, generating a distribution function of both bound and unbound populations of dark matter particles. From this distribution function, we calculate annihilation rates and observable gamma-ray spectra for a few simple dark matter models. The features of these spectra are sensitive to the black hole spin, observer inclination, and detailed properties of the dark matter annihilation cross section and density profile. Confirming earlier analytic work, we find that for rapidly spinning black holes, the collisional Penrose process can reach efficiencies exceeding $600\%$, leading to a high-energy tail in the annihilation spectrum. The high particle density and large proper volume of the region immediately surrounding the horizon ensures that the observed flux from these extreme events is non-negligible.

We consider the evolution of a supermassive black hole binary (SMBHB) surrounded by a retrograde accretion disk. Assuming the disk is exactly in the binary plane and transfers energy and angular momentum to the binary via direct gas accretion, we calculate the time evolution of the binary's semi-major axis $a$ and eccentricity $e$. Because the gas is predominantly transferred when the binary is at apocenter, we find the eccentricity grows rapidly while maintaining constant $a(1+e)$. After accreting only a fraction of the secondary's mass, the eccentricity grows to nearly unity; from then on, gravitational wave emission dominates the evolution, preserving constant $a(1-e)$. The high-eccentricity waveforms redistribute the peak gravitational wave power from the nHz to $\mu$Hz bands, substantially affecting the signal that might be detected with pulsar timing arrays. We also estimate the torque coupling binaries of arbitrary eccentricity with obliquely aligned circumbinary disks. If the outer edge of the disk is not an extremely large multiple of the binary separation, retrograde accretion can drive the binary into the gravitational wave-dominated state before these torques align the binary with the angular momentum of the mass supply.

Supermassive black hole binaries are the strongest gravitational wave sources in the universe. The systems most likely to be observed with pulsar timing arrays (PTAs) will have particularly high masses ($\gtrsim 10^9 M_\odot$), long periods ($T_{\rm orb} \gtrsim 1$ yr), and be in the local universe ($z \lesssim 1$). These features are also the most favorable for bright electromagnetic counterparts, which should be easily observable with existing ground- and space-based telescopes. Wide-field X-ray observatories such as ISS-Lobster will provide independent candidates that can be used to lower the threshold for PTA detections of resolvable binary sources. The primary challenge lies in correctly identifying and characterizing binary sources with long orbital periods, as opposed to "normal" active galactic nuclei (AGN) hosting single black holes. Here too ISS-Lobster will provide valuable new understanding into the wide range of behaviors seen in AGN by vastly expanding our sample of X-ray light curves from accreting supermassive black holes.

We present here an overview of recent work in the subject of astrophysical manifestations of super-massive black hole (SMBH) mergers. This is a field that has been traditionally driven by theoretical work, but in recent years has also generated a great deal of interest and excitement in the observational astronomy community. In particular, the electromagnetic (EM) counterparts to SMBH mergers provide the means to detect and characterize these highly energetic events at cosmological distances, even in the absence of a space-based gravitational-wave observatory. In addition to providing a mechanism for observing SMBH mergers, EM counterparts also give important information about the environments in which these remarkable events take place, thus teaching us about the mechanisms through which galaxies form and evolve symbiotically with their central black holes.

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

In this document, we describe the scientific potential of blazar observations with a X-ray polarimetry mission like GEMS (Gravity and Extreme Magnetism SMEX). We describe five blazar science investigations that such a mission would enable: (i) the structure and the role of magnetic fields in AGN jets, (ii) analysis of the polarization of the synchrotron X-ray emission from AGN jets, (iii) discrimination between synchrotron self-Compton and external Compton models for blazars with inverse Compton emission in the X-ray band, (iv) a precision study of the polarization properties of the X-ray emission from Cen-A, (v) tests of Lorentz Invariance based on X-ray polarimetric observations of blazars. We conclude with a discussion of a straw man observation program and recommended accompanying multiwavelength observations.

We present a new code for radiation transport around Kerr black holes, including arbitrary emission and absorption terms, as well as electron scattering and polarization. The code is particularly useful for analyzing accretion flows made up of optically thick disks and optically thin coronae. We give a detailed description of the methods employed in the code, and also present results from a number of numerical tests to assess its accuracy and convergence.

We examine the expected X-ray polarization properties of neutron-star X-ray sources of various types, e.g., accretion and rotation powered pulsars, magnetars, and low-mass X-ray binaries. We summarize the model calculations leading to these expected properties. We describe how a comparison of these with their observed properties, as inferred from GEMS data, will probe the essential dynamical, electromagnetic, plasma, and emission processes in neutron-star binaries, discriminate between models of these processes, and constrain model parameters. An exciting goal is the first observational demonstration in this context of the existence of vacuum resonance, a fundamental quantum electrodynamical phenomenon first described in the 1930s.

We present here a summary of the scientific goals behind the Gravity and Extreme Magnetism SMEX (GEMS) X-ray polarimetry mission's black hole (BH) observing program. The primary targets can be divided into two classes: stellar-mass galactic BHs in accreting binaries, and super-massive BHs in the centers of active galactic nuclei (AGN). The stellar-mass BHs can in turn be divided into various X-ray spectral states: thermal-dominant (disk), hard (radio jet), and steep power-law (hot corona). These different spectral states are thought to be generated by different accretion geometries and emission mechanisms. X-ray polarization is an ideal tool for probing the geometry around these BHs and revealing the specific properties of the accreting gas.

In this white paper, we discuss the concept of a next-generation X-ray mission called BEST (Black hole Evolution and Space Time). The mission concept uses a 3000 square centimeter effective area mirror (at 6 keV) to achieve unprecedented sensitivities for hard X-ray imaging spectrometry (5-70 keV) and for broadband X-ray polarimetry (2-70 keV). BEST can make substantial contributions to our understanding of the inner workings of accreting black holes, our knowledge about the fabric of extremely curved spacetime, and the evolution of supermassive black holes. BEST will allow for time resolved studies of accretion disks. With a more than seven times larger mirror area and a seven times wider bandpass than GEMS, BEST will take X-ray polarimetry to a new level: it will probe the time variability of the X-ray polarization from stellar mass and supermassive black holes, and it will measure the polarization properties in 30 independent energy bins. These capabilities will allow BEST to conduct tests of accretion disk models and the underlying spacetimes. With three times larger mirror area and ten times better angular resolution than NuSTAR, BEST will be able to make deep field observations with a more than 15 times better sensitivity than NuSTAR. The mission will be able to trace the evolution of obscured and unobscured black holes in the redshift range from zero to six, covering the most important epoch of supermassive black hole growth. The hard X-ray sensitivity of BEST will enable a deep census of non-thermal particle populations. BEST will give us insights into AGN feedback by measuring the particle luminosity injected by AGNs into the interstellar medium (ISM) of their hosts, and will map the emission from particles accelerated at large scale structure shocks. Finally, BEST has the potential to constrain the equation of state of neutron stars (NS).

X-ray polarimetry promises to give new information about high-energy astrophysical sources, such as binary black hole systems, micro-quasars, active galactic nuclei, and gamma-ray bursts. We designed, built and tested a hard X-ray polarimeter X-Calibur to be used in the focal plane of the InFOCuS grazing incidence hard X-ray telescope. X-Calibur combines a low-Z Compton scatterer with a CZT detector assembly to measure the polarization of 10-80 keV X-rays making use of the fact that polarized photons Compton scatter preferentially perpendicular to the electric field orientation. X-Calibur achieves a high detection efficiency of order unity.

X-ray polarimetry promises to give qualitatively new information about high-energy astrophysical sources, such as binary black hole systems, micro-quasars, active galactic nuclei, and gamma-ray bursts. We designed, built and tested a hard X-ray polarimeter X-Calibur to be used in the focal plane of the InFOCuS grazing incidence hard X-ray telescope. X-Calibur combines a low-Z Compton scatterer with a CZT detector assembly to measure the polarization of 10-80 keV X-rays making use of the fact that polarized photons Compton scatter preferentially perpendicular to the electric field orientation. X-Calibur achieves a high detection efficiency of order unity.

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

During the final moments of a binary black hole (BH) merger, the gravitational wave (GW) luminosity of the system is greater than the combined electromagnetic output of the entire observable universe. However, the extremely weak coupling between GWs and ordinary matter makes these waves very difficult to detect directly. Fortunately, the inspiraling BH system will interact strongly--on a purely Newtonian level--with any surrounding material in the host galaxy, and this matter can in turn produce unique electromagnetic (EM) signals detectable at Earth. By identifying EM counterparts to GW sources, we will be able to study the host environments of the merging BHs, in turn greatly expanding the scientific yield of a mission like LISA.

We calculate the location and stability of the L_4 and L_5 Lagrange equilibrium points in the circular restricted three-body problem as the binary system evolves via gravitational radiation losses. Relative to the purely Newtonian case, we find that the L_4 equilibrium point moves towards the secondary mass and becomes slightly less stable, while the L_5 point moves away from the secondary and gains in stability. We discuss a number of astrophysical applications of these results, in particular as a mechanism for producing electromagnetic counterparts to gravitational-wave signals.

We present new calculations of X-ray polarization from accreting black holes (BHs), using a Monte-Carlo ray-tracing code in full general relativity. In our model, an optically thick disk in the BH equatorial plane produces thermal seed photons with polarization oriented parallel to the disk surface. These seed photons are then inverse-Compton scattered through a hot (but thermal) corona, producing a hard X-ray power-law spectrum. We consider three different models for the corona geometry: a wedge ``sandwich'' with aspect ratio H/R and vertically-integrated optical depth tau_0 constant throughout the disk; an inhomogeneous ``clumpy'' corona with a finite number of hot clouds distributed randomly above the disk within a wedge geometry; and a spherical corona of uniform density, centered on the BH and surrounded by a truncated thermal disk with inner radius R_edge. In all cases we find a characteristic transition from horizontal polarization at low energies to vertical polarization above the thermal peak; the vertical direction is defined as the projection of the BH spin axis on the plane of the sky. We show how the details of the spectropolarization signal can be used to distinguish between these models and infer various properties of the corona and BH. Although the bulk of this paper focuses on stellar-mass BHs, we also consider the effects of coronal scattering on the X-ray polarization signal from supermassive BHs in active galactic nuclei.

We present new calculations of X-ray polarization from black hole accretion disks in the thermally-dominated state, using a Monte-Carlo ray-tracing code in full general relativity. In contrast to many previously published studies, our approach allows us to include returning radiation that is deflected by the strong-field gravity of the BH and scatters off of the disk before reaching a distant observer. Although carrying a relatively small fraction of the total observed flux, the scattered radiation tends to be highly polarized and in a direction perpendicular to the direct radiation. We show how these new features of the polarization spectra may be developed into a powerful tool for measuring black hole spin and probing the gas flow in the innermost disk.

The association of an electromagnetic signal with the merger of a pair of supermassive black holes would have many important implications. For example, it would provide new information about gas and magnetic field interactions in dynamical spacetimes as well as a combination of redshift and luminosity distance that would enable precise cosmological tests. A proposal first made by Bode & Phinney (2007) is that because radiation of gravitational waves during the final inspiral and merger of the holes is abrupt and decreases the mass of the central object by a few percent, there will be waves in the disk that can steepen into shocks and thus increase the disk luminosity in a characteristic way. We evaluate this process analytically and numerically. We find that shocks only occur when the fractional mass loss exceeds the half-thickness (h/r) of the disk, hence significant energy release only occurs for geometrically thin disks which are thus at low Eddington ratios. This strongly limits the effective energy release, and in fact our simulations show that the natural variations in disk luminosity are likely to obscure this effect entirely. However, we demonstrate that the reduction of luminosity caused by the retreat of the inner edge of the disk following mass loss is potentially detectable. This decrease occurs even if the disk is geometrically thick, and lasts for a duration on the order of the viscous time of the modified disk. Observationally, the best prospect for detection would be a sensitive future X-ray instrument with a field of view of on the order of a square degree, or possibly a wide-field radio array such as the Square Kilometer Array, if the disk changes produce or interrupt radio emission from a jet.

A supermassive black hole ejected from the center of a galaxy by gravitational wave recoil carries a retinue of bound stars - a "hypercompact stellar system" (HCSS). The numbers and properties of HCSSs contain information about the merger histories of galaxies, the late evolution of binary black holes, and the distribution of gravitational-wave kicks. We relate the structural properties of HCSSs to the properties of their host galaxies, in two regimes: collisional, i.e. short nuclear relaxation times; and collisionless, i.e. long nuclear relaxtion times. HCSSs are expected to be similar in size and luminosity to globular clustersbut in extreme cases their stellar mass can approach that of UCDs. However they differ from all other classes of compact stellar system in having very high internal velocities. We show that the kick velocity is encoded in the velocity dispersion of the bound stars. Given a large enough sample of HCSSs, the distribution of gravitational-wave kicks can therefore be empirically determined. We combine a hierarchical merger algorithm with stellar population models to compute the rate of production of HCSSs over time and the probability of observing HCSSs in the local universe as a function of their apparent magnitude, color, size and velocity dispersion, under two assumptions about the star formation history prior to the kick. We predict that roughly 100 should be detectable within 2 Mpc of the center of the Virgo cluster and that many of these should be bright enough that their high internal velocity dispersions could be measured with reasonable exposure times.

X-ray timing observations of accreting stellar mass black holes have shown that they can produce signals with such short time scales that we must be probing very close to the innermost stable circular orbit that is predicted by the theory of General Relativity (GR). These signals are quasi-periodic oscillations (QPOs), and both the high-frequency variety (HFQPOs, which have frequencies in the 40-450 Hz range) as well as the 0.1-10 Hz low-frequency type have the potential to provide tests of GR in the strong field limit. An important step on the path to GR tests is to constrain the physical black hole properties, and the straightforward frequency measurements that are possible with X-ray timing may provide one of the cleanest measurements of black hole spins. While current X-ray satellites have uncovered these phenomenona, the HFQPOs are weak signals, and future X-ray timing missions with larger effective area are required for testing the candidate theoretical QPO mechanisms. Another main goal in the study of accreting black holes is to understand the production of relativistic jets. Here, we have also made progress during the past decade by finding clear connections between the radio emission that traces the strength of the jet and the properties of the X-ray emission. With new radio capabilities just coming on-line, continuing detailed X-ray studies of accreting black holes is crucial for continuing to make progress.

We present new calculations of X-ray polarization from black hole (BH) accretion disks in the thermally-dominated state, using a Monte-Carlo ray-tracing code in full general relativity. In contrast to many previously published studies, our approach allows us to include returning radiation that is deflected by the strong-field gravity of the BH and scatters off of the disk before reaching a distant observer. Although carrying a relatively small fraction of the total observed flux, the scattered radiation tends to be highly polarized and in a direction perpendicular to the direct radiation. For moderately large spin parameters (a/M >~ 0.9), this scattered returning radiation dominates the polarization signal at energies above the thermal peak, giving a net rotation in the polarization angle of 90 deg. We show how these new features of the polarization spectra from BHs in the thermal state may be developed into a powerful tool for measuring BH spin and probing the gas flow in the innermost disk. In addition to determining the emission profile, polarization observations can be used to constrain other properties of the system such as BH mass, inclination, and distance. New instruments currently under development should be able to exploit this tool in the near future.

It is widely expected that the coming decade will witness the first direct detection of gravitational waves (GWs). The ground-based LIGO and Virgo GW observatories are being upgraded to advanced sensitivity, and are expected to observe a significant binary merger rate. The launch of The Laser Interferometer Space Antenna (LISA) would extend the GW window to low frequencies, opening new vistas on dynamical processes involving massive (M >~ 10^5 M_Sun) black holes. GW events are likely to be accompanied by electromagnetic (EM) counterparts and, since information carried electromagnetically is complementary to that carried gravitationally, a great deal can be learned about an event and its environment if it becomes possible to measure both forms of radiation in concert. Measurements of this kind will mark the dawn of trans-spectral astrophysics, bridging two distinct spectral bands of information. The aim of this whitepaper is to articulate future directions in both theory and observation that are likely to impact broad astrophysical inquiries of general interest. What will EM observations reflect on the nature and diversity of GW sources? Can GW sources be exploited as complementary probes of cosmology? What cross-facility coordination will expand the science returns of gravitational and electromagnetic observations?

SDSSJ092712.65+294344.0 was identified by the SDSS as a quasar, but has the unusual property of having two emission line systems offset by 2650 km/s. One of these contains the usual combination of broad and narrow lines, the other only narrow lines. In the first paper commenting on this system (Komossa et al. 2008), it was interpreted as a galaxy in which a pair of black holes had merged, imparting a several thousand km/s recoil to the new, larger black hole. In two other papers (Bogdanovic, Eracleous & Sigurdsson 2008; Dotti et al. 2008), it was interpreted as a small-separation binary black hole. We propose a new interpretation: that this system is a more distant analog of NGC1275, a large and a small galaxy interacting near the center of a rich cluster.

We model the spectra and light curves of circumbinary accretion disks during the time after the central black holes merge. The most immediate effect of this merger is the dissipation of energy in the outer regions of the disk due to the gravitational wave energy and linear momentum flux released at merger. This has the effect of perturbing the gas in the disk, which then radiates the dissipated energy over a cooling timescale, giving a characteristic infrared signal for tens of thousands of years when the total black hole mass is M~10^8 M_sun. On the basis of a simple cosmological merger model in which a typical supermassive black hole undergoes a few major mergers during its lifetime, we predict that ~10^4-10^5 of these IR sources should be observable today and discuss the possibility of identifying them with multi-wavelength surveys such as SWIRE/XMM-LSS/XBootes and COSMOS.

We present a multipolar analysis of the gravitational recoil computed in recent numerical simulations of binary black hole (BH) coalescence, for both unequal masses and non-zero, non-precessing spins. We show that multipole moments up to and including l=4 are sufficient to accurately reproduce the final recoil velocity (within ~2%) and that only a few dominant modes contribute significantly to it (within ~5%). We describe how the relative amplitudes, and more importantly, the relative phases, of these few modes control the way in which the recoil builds up throughout the inspiral, merger, and ringdown phases. We also find that the numerical results can be reproduced by an ``effective Newtonian'' formula for the multipole moments obtained by replacing the radial separation in the Newtonian formulae with an effective radius computed from the numerical data. Beyond the merger, the numerical results are reproduced by a superposition of three Kerr quasi-normal modes (QNMs). Analytic formulae, obtained by expressing the multipole moments in terms of the fundamental QNMs of a Kerr BH, are able to explain the onset and amount of ``anti-kick'' for each of the simulations. Lastly, we apply this multipolar analysis to help explain the remarkable difference between the amplitudes of planar and non-planar kicks for equal-mass spinning black holes.

Recent numerical simulations of binary black hole mergers show the possibility of producing very large recoil velocities (> 3000 km/s). Kicks of this magnitude should be sufficient to eject the final black hole from virtually any galactic potential. This result has been seen as a potential contradiction with observations of supermassive black holes residing in the centers of most galaxies in the local universe. Using an extremely simplified merger tree model, we show that, even in the limit of very large ejection probability, after a small number of merger generations there should still be an appreciable fraction (>50%) of galaxies with supermassive black holes today. We go on to argue that the inclusion of more realistic physics ingredients in the merger model should systematically increase this retention fraction, helping to resolve a potential conflict between theory and observation. Lastly, we develop a more realistic Monte Carlo model to confirm the qualitative arguments and estimate occupation fractions as a function of the central galactic velocity dispersion.

We calculate the linear momentum flux from merging black holes (BHs) with arbitrary masses and spin orientations, using the effective-one-body (EOB) model. This model includes an analytic description of the inspiral phase, a short merger, and a superposition of exponentially damped quasi-normal ringdown modes of a Kerr BH. By varying the matching point between inspiral and ringdown, we can estimate the systematic errors generated with this method. Within these confidence limits, we find close agreement with previously reported results from numerical relativity. Using a Monte Carlo implementation of the EOB model, we are able to sample a large volume of BH parameter space and estimate the distribution of recoil velocities. For a range of mass ratios 1 <= m_1/m_2 <= 10, spin magnitudes of a_{1,2}=0.9, and uniform random spin orientations, we find that a fraction f_{500}=0.12^{+0.06}_{-0.05} of binaries have recoil velocities greater than 500 km/s and f_{1000}=0.027^{+0.021}_{-0.014} have kicks greater than 1000 km/s. These velocities likely are capable of ejecting the final BH from its host galaxy. Limiting the sample to comparable-mass binaries with m_1/m_2 <= 4, the typical kicks are even larger, with f_{500}=0.31_{-0.12}^{+0.13} and f_{1000}=0.079^{+0.062}_{-0.042}.

We consider a model in which massive stars form in a self-gravitating accretion disk around an active galactic nucleus (AGN). These stars may evolve and collapse to form compact objects on a time scale shorter than the accretion time, thus producing an important family of sources for LISA. Assuming the compact object formation/inspiral rate is proportional to the steady-state gas accretion rate, we use the intrinsic hard X-ray AGN luminosity function to estimate expected event rates and signal strengths. We find that these sources will produce a continuous low-frequency (<~ mHz) background detectable by LISA if more than 1% of the accreted matter is in the form of compact objects. For compact objects with masses >~ 10 solar masses the last stages of the inspiral events should be resolvable above a few mHz, with rates as high as a few hundred per year.

We consider a model in which massive stars form in a self-gravitating accretion disk around an active galactic nucleus. These stars may evolve and collapse to form compact objects on a time scale shorter than the accretion time, thus producing an important family of sources for LISA. Assuming the compact object formation/inspiral rate is proportional to the steady-state gas accretion rate, we use the observed extra-galactic X-ray luminosity function to estimate expected event rates and signal strengths. We find that these sources will produce a continuous low-frequency background detectable by LISA if more than >~ 1% of the accreted matter is in the form of compact objects. For compact objects with m >~ 10 Msun, the last stages of the inspiral events should be resolvable above a few mHz, at a rate of ~10-100 per year.

We use a relativistic ray-tracing code to calculate the light curves observed from a global general relativistic magneto-hydrodynamic simulation of an accretion flow onto a Schwarzschild black hole. We apply three basic emission models to sample different properties of the time-dependent accretion disk. With one of these models, which assumes thermal blackbody emission and free-free absorption, we can predict qualitative features of the high-frequency power spectrum from stellar-mass black holes in the "Thermal Dominant" state. The simulated power spectrum is characterized by a power law of index Gamma ~ 3 and total rms fractional variance of <~ 2% above 10 Hz. For each emission model, we find that the variability amplitude should increase with increasing inclination angle. On the basis of a newly-developed formalism for quantifying the significance of quasi-periodic oscillations (QPOs) in simulation data, we find that these simulations are able to identify any such features with (rms/mean) amplitudes >~ 1 % near the orbital frequency at the inner-most stable orbit. Initial results indicate the existence of transient QPO peaks with frequency ratios of nearly 2:3 at a 99.9% confidence limit, but they are not generic features because at any given time they are seen only from certain observer directions. Additionally, we present detailed analysis of the azimuthal structure of the accretion disk and the evolution of density perturbations in the inner disk. These "hot spot" structures appear to be roughly self-similar over a range of disk radii, with a single characteristic size \delta\phi=25 deg and \delta r/r=0.3, and typical lifetimes T_l ~ 0.3 T_orb.

This Thesis describes the basic framework and applications of a relativistic ray-tracing code for analyzing accretion processes around Kerr black holes. We begin in Chapter 1 with a brief historical summary of the major advances in black hole astrophysics over the past few decades. In Chapter 2 we present a detailed description of the ray-tracing code, which is used to calculate the transfer function between the accretion disk and the detector. In Chapter 3, we employ a simple ``hot spot'' model to explain the frequencies and amplitudes of quasi-periodic oscillations (QPOs). In Chapter 4, we introduce additional features to the hot spot model to explain the broadening of the QPO peaks as well as the damping of higher-frequency harmonics in the power spectrum. In Chapter 5 we present a description of the structure of a relativistic alpha-disk around a Kerr black hole, and the observed spectrum from such a disk. The features of this modified thermal spectrum may be used to infer the physical properties of the accretion disk and the central black hole. In Chapter 6 we develop a Monte Carlo code to calculate the detailed propagation of photons from a hot spot emitter scattering through a corona surrounding the black hole. The coronal scattering has two major observable effects: the inverse-Compton process alters the photon spectrum by adding a high energy power-law tail, and the random scattering of each photon effectively damps out the highest frequency modulations in the X-ray light curve.

We develop a simple physical model to describe the most common type of low-frequency quasi-periodic oscillations (QPOs) seen in a number of accreting black hole systems, as well as the shape of the relativistically broadened iron emission lines that often appear simultaneously in such sources. The model is based on an inclined ring of hot gas that orbits the black hole along geodesic trajectories. For spinning black holes, this ring will precess around the spin axis of the black hole at the Lense-Thirring (``frame-dragging'') frequency. Using a relativistic ray-tracing code, we calculate X-ray light curves and observed energy spectra as a function of the radius and tilt angle of the ring, the spin magnitude, and the inclination of the black hole. The model predicts higher-amplitude QPOs for systems with high inclinations, as seen in a growing number of black hole binary systems. We find that the Rossi X-ray Timing Explorer observations of low-frequency QPOs in GRS 1915+105 are consistent with a ring of radius R ~ 10M orbiting a black hole with spin a/M ~0.5 and inclination angle of i ~ 70 deg. Finally, we describe how future X-ray missions may be able to use simultaneous timing and spectroscopic observations to measure the black hole spin and probe the inner-most regions of the accretion disk.

We use a relativistic ray-tracing code to analyze the X-ray emission from a pressure-supported oscillating relativistic torus around a black hole. We show that a strong correlation exists between the {\it intrinsic} frequencies of the torus normal modes and the {\it extrinsic} frequencies seen in the observed light curve power spectrum. This correlation demonstrates the feasibility of the oscillating-torus model to explain the multiple peaks seen in black hole high-frequency quasi-periodic oscillations. Using an optically thin, monochromatic emission model, we also determine how a relativistically broadened emission line and the amplitude of the X-ray modulations are dependent on the observer's inclination angle and on the torus oscillation amplitudes. Observations of these features can provide important information about the torus as well as the black hole.

Using a Monte Carlo ray-tracing code in full general relativity, we calculate the transport of photons from a geodesic hot spot emitter through a corona of hot electrons surrounding a black hole. Each photon is followed until it is either captured by the black hole or is detected by a distant observer. The source is assumed to be a low-energy thermal emitter (T_em ~ 1 keV), isotropic in the rest frame of a massive geodesic test particle. The coronal scattering has two major observable effects: the Comptonization of the photon spectrum due to the high energy electrons, and the convolution of the time-dependent light curve as each photon is effectively scattered into a different time bin. Both of these effects are clearly present in the Rossi X-Ray Timing Explorer (RXTE) observations of high frequency quasi-periodic oscillations (QPOs) seen in black hole binaries. These QPOs tend to occur when the system is in the Steep Power Law spectral state and also show no evidence for significant power at higher harmonic frequencies, consistent with the smoothing out of the light curve by multiple random time delays. We present simulated photon spectra and light curves and compare with RXTE data, allowing us to infer the properties of the corona as well as the hot spot emitter.

We discuss the use of the bicoherence - a measure of the phase coupling of oscillations at different frequencies - as a diagnostic between different models for high frequency quasi-periodic oscillations from Galactic black hole candidates. We show that this statistic is capable of finding qualitative distinctions between different hot spot models which produce nearly identical Fourier power density spectra. Finally, we show that proposed new timing missions should detect enough counts to make real use of this statistic.

Starting with a post-Newtonian description of compact binary systems, we derive a set of equations that describes the evolution of the orbital angular momentum and both spin vectors during inspiral. We find regions of phase space that exhibit resonance behavior, characterized by small librations of the spin vectors around a fixed orientation. Due to the loss of energy and orbital angular momentum through radiation reaction, systems can eventually be captured into these resonance orientations. By investigating the long-term evolution of compact binaries with a variety of initial conditions, we find that the distribution in parameter space can be strongly affected by resonance captures. This has the effect of significantly reducing the size of search space for gravitational wave sources, in turn improving the chances of detecting such sources through methods of template matching. Furthermore, by calculating the expected spin distribution at the end of the inspiral phase, we can predict what are the most likely initial conditions for the plunge phase, a result of great interest for numerical relativity calculations.

In the context of a relativistic hot spot model, we investigate different physical mechanisms to explain the behavior of quasi-periodic oscillations (QPOs) from accreting black holes. The locations and amplitudes of the QPO peaks are determined by the ray-tracing calculations presented in Schnittman & Bertschinger (2004a): the black hole mass and angular momentum give the geodesic coordinate frequencies, while the disk inclination and the hot spot size, shape, and overbrightness give the amplitudes of the different peaks. In this paper additional features are added to the existing model to explain the broadening of the QPO peaks as well as the damping of higher frequency harmonics in the power spectrum. We present a number of analytic results that closely agree with more detailed numerical calculations. Four primary pieces are developed: the addition of multiple hot spots with random phases, a finite width in the distribution of geodesic orbits, Poisson sampling of the detected photons, and the scattering of photons from the hot spot through a corona of hot electrons around the black hole. Finally, the complete model is used to fit the observed power spectra of both type A and type B QPOs seen in XTE J1550-564, giving confidence limits on each of the model parameters.

In at least two black hole binary systems, the Rossi X-Ray Timing Explorer has detected high frequency quasi-periodic oscillations (HFQPOs) with a 2:3 frequency commensurability. We propose a simple hot spot model to explain the positions, amplitudes, and widths of the HFQPO peaks. Using the exact geodesic equations for the Kerr metric, we calculate the trajectories of massive test particles, which are treated as isotropic, monochromatic emitters in their rest frames. By varying the hot spot parameters, we are able to explain the different features observed in ``Type A'' and ``Type B'' QPOs from XTE J1550-564. In the context of this model, the observed power spectra allow us to infer values for the black hole mass and angular momentum, and also constrain the parameters of the model.

Observations from the Rossi X-Ray Timing Explorer have shown the existence of high frequency quasi-periodic oscillations (HFQPOs) in the X-ray flux from accreting black hole binary systems. In at least two systems, these HFQPOs come in pairs with a 2:3 frequency commensurability. We propose a simple ``hot spot'' model to explain the position and amplitude of the HFQPO peaks. Using the exact geodesic equations for the Kerr metric, we calculate the trajectories of massive test particles, which are treated as isotropic, monochromatic emitters in their rest frames. Photons are traced from the accretion disk to a distant observer to produce time- and frequency-dependent images of the orbiting hot spot and background disk. The power spectrum of the X-ray light curve consists of multiple peaks at integral combinations of the black hole coordinate frequencies. In particular, if the radial frequency is one-third of the azimuthal frequency (as is the case near the inner-most stable circular orbit), beat frequencies appear in the power spectrum at two-thirds and four-thirds of the fundamental azimuthal orbital frequency, in agreement with observations. We also model the effects of shearing the hot spot in the disk, producing an arc of emission that also follows a geodesic orbit, as well as the effects of non-planar orbits that experience Lens-Thirring precession around the black hole axis. By varying the arc length, we are able to explain the different QPO features observed in ``Type A'' and ``Type B'' X-ray outbursts from XTE J1550-564. In the context of this model, the observed power spectra allow us to infer values for the black hole mass and angular momentum, and also constrain the parameters of the model, such as the hot spot size and luminosity.

We investigate the orbits of compact binary systems during the final inspiral period before coalescence by integrating numerically the second-order post-Newtonian equations of motion. We include spin-orbit and spin-spin coupling terms, which, according to a recent study by Levin [J. Levin, Phys. Rev. Lett. 84, 3515 (2000)], may cause the orbits to become chaotic. To examine this claim, we study the divergence of initially nearby phase-space trajectories and attempt to measure the Lyapunov exponent gamma. Even for systems with maximally spinning objects and large spin-orbit misalignment angles, we find no chaotic behavior. For all the systems we consider, we can place a strict lower limit on the divergence time t_L=1/gamma that is many times greater than the typical inspiral time, suggesting that chaos should not adversely affect the detection of inspiral events by upcoming gravitational-wave detectors.