Search

Among the potential milliHz gravitational wave (GW) sources for the upcoming space-based interferometer LISA are extreme- or intermediate-mass ratio inspirals (EMRI/IMRIs). These events involve the coalescence of supermassive black holes in the mass range $10^5 M_{\odot} \lesssim M \lesssim 10^7 M_{\odot}$ with companion BHs of much lower masses. A subset of E/IMRIs are expected to occur in the accretion discs of active galactic nuclei (AGN), where torques exerted by the disc can interfere with the inspiral and cause a phase shift in the GW waveform. Here we use a suite of two-dimensional hydrodynamical simulations with the moving-mesh code DISCO to present a systematic study of disc torques. We measure torques on an inspiraling BH and compute the corresponding waveform deviations as a function of the binary mass ratio $q\equiv M_2/M_1$, the disc viscosity ($\alpha$), and gas temperature (or equivalently Mach number; $\mathcal{M}$). We find that the absolute value of the gas torques is within an order of magnitude of previously determined planetary migration torques, but their precise value and sign depends non-trivially on the combination of these parameters. The gas imprint is detectable by LISA for binaries embedded in AGN discs with surface densities above $\Sigma_0\ge10^{4-6} \rm \, g cm^{-2}$, depending on $q$, $\alpha$ and $\mathcal{M}$. Deviations are most pronounced in discs with higher viscosities, and for E/IMRIs detected at frequencies where LISA is most sensitive. Torques in colder discs exhibit a noticeable dependence on the GW-driven inspiral rate as well as strong fluctuations at late stages of the inspiral. Our results further suggest that LISA may be able to place constraints on AGN disc parameters and the physics of disc-satellite interaction.

Variable active galactic nuclei showing periodic light curves have been proposed as massive black hole binary (MBHB) candidates. In such scenarios the periodicity can be due to relativistic Doppler-boosting of the emitted light. This hypothesis can be tested through the timing of scattered polarized light. Following the results of polarization studies in type I nuclei and of dynamical studies of MBHBs with circumbinary discs, we assume a coplanar equatorial scattering ring, whose elements contribute differently to the total polarized flux, due to different scattering angles, levels of Doppler boost, and line-of-sight time delays. We find that in the presence of a MBHB, both the degree of polarization and the polarization angle have periodic modulations. The minimum of the polarization degree approximately coincides with the peak of the light curve, regardless of the scattering ring size. The polarization angle oscillates around the semi-minor axis of the projected MBHB orbital ellipse, with a frequency equal either to the binary's orbital frequency (for large scattering screen radii), or twice this value (for smaller scattering structures). These distinctive features can be used to probe the nature of periodic MBHB candidates and to compile catalogs of the most promising sub-pc MBHBs. The identification of such polarization features in gravitational-wave detected MBHBs would enormously increase the amount of physical information about the sources, allowing the measurement of the individual masses of the binary components, and the orientation of the line of nodes on the sky, even for monochromatic gravitational wave signals., 11 pages, 6 figures, submitted to MNRAS

We demonstrate how pulsar timing arrays (PTAs) can, in principle, yield a purely gravitational wave (GW) measurement of the luminosity distance and comoving distance to a supermassive black hole binary source, hence providing an estimate of the source redshift and the Hubble constant. The luminosity distance is derived through standard measurement of the chirp mass, which for the slowly evolving binary sources in the PTA band can be found by comparing the frequency of GW-timing residuals at the Earth compared to those at distant pulsars in the array. The comoving distance can be measured from GW-timing parallax caused by the curvature of the GW wavefronts. This can be detected for single sources at the high-frequency end of the PTA band out to Gpc distances with a future PTA containing well-timed pulsars out to $\mathcal{O}(10)$ kpc, when the pulsar distance is constrained to less than a GW wavelength. Such a future PTA, with $\gtrsim 30$ pulsars with precise distance measurements between 1 and 20~kpc, could measure the Hubble constant at the tens of percent level for a single source at $0.1 \lesssim z \lesssim 1.5$. At $z \lesssim 0.1$, the luminosity and comoving distances are too similar to disentangle, unless the fractional error in the luminosity distance measurement is decreased below 10%. At $z\gtrsim 1.5$, this measurement will likely be restricted by a signal-to-noise ratio threshold. Generally, clarification of the different types of cosmological distances that can be probed by PTAs, and their relation to pulsar distance measurements is important for ongoing PTA experiments aimed at detecting and characterizing GWs., Submitted to Physical Review D. Comments welcome

The past year has seen numerous publications underlining the importance of a space mission to the ice giants in the upcoming decade. Proposed mission plans involve a $\sim$10 year cruise time to the ice giants. This cruise time can be utilized to search for low-frequency gravitational waves (GWs) by observing the Doppler shift caused by them in the Earth-spacecraft radio link. We calculate the sensitivity of prospective ice giant missions to GWs. Then, adopting a steady-state black hole binary population, we derive a conservative estimate for the detection rate of extreme mass ratio inspirals (EMRIs), supermassive- (SMBH) and stellar mass binary black hole (sBBH) mergers. We link the SMBH population to the fraction of quasars $f_\rm{bin}$ resulting from galaxy mergers that pair SMBHs to a binary. For a total of ten 40-day observations during the cruise of a single spacecraft, $\mathcal{O}(f_\rm{bin})\sim0.5$ detections of SMBH mergers are likely, if Allan deviation of Cassini-era noise is improved by $\sim 10^2$ in the $10^{-5}-10^{-3}$ Hz range. For EMRIs the number of detections lies between $\mathcal{O}(0.1) - \mathcal{O}(100)$. Furthermore, ice giant missions combined with the Laser Interferometer Space Antenna (LISA) would improve the localisation by an order of magnitude compared to LISA by itself., Accepted for publication in MNRAS: Letters

The community may be on the verge of detecting low-frequency gravitational waves from massive black hole binaries (MBHBs), but no examples of binary active galactic nuclei (AGN) have been confirmed. Because MBHBs are intrinsically rare, the most promising detection methods utilize photometric data from all-sky surveys. Recently D'Orazio & Di Stefano 2018 (arXiv:1707.02335) suggested gravitational self-lensing as a method of detecting AGN in close separation binaries. In this study we calculate the detectability of lensing signatures in realistic populations of simulated MBHBs. Within our model assumptions, we find that VRO's LSST should be able to detect 10s to 100s of self-lensing binaries, with the rate uncertainty depending primarily on the orientation of AGN disks relative to their binary orbits. Roughly a quarter of lensing detectable systems should also show detectable Doppler boosting signatures. If AGN disks tend to be aligned with the orbit, lensing signatures are very nearly achromatic, while in misaligned configurations the bluer optical bands are lensed more than redder ones. Whether substantial obscuring material (e.g.~a dusty torus) will be present in close binaries remains uncertain, but our estimates suggest that a substantial fraction of systems would still be observable in this case., 10+3 pages, 4+2 figures, comments welcome and appreciated

We examine the light curves of two quasars, motivated by recent suggestions that a supermassive black hole binary (SMBHB) can exhibit sharp lensing spikes. We model the variability of each light curve as due to a combination of two relativistic effects: the orbital relativistic Doppler boost and gravitational binary self-lensing. In order to model each system we extend previous Doppler plus self-lensing models to include eccentricity. The first quasar is identified in optical data as a binary candidate with a 20-yr period (Ark 120), and shows a prominent spike. For this source, we rule out the lensing hypothesis and disfavor the Doppler-boost hypothesis due to discrepancies in the measured vs. recovered values of the binary mass and optical spectral slope. The second source, which we nickname Spikey, is the rare case of an active galactic nucleus (AGN) identified in Kepler's high-quality, high-cadence photometric data. For this source, we find a model, consisting of a combination of a Doppler modulation and a narrow symmetric lensing spike, consistent with an eccentric SMBHB with a total mass of approximately 30 million solar masses, rest-frame orbital period T=418 days, eccentricity e=0.5, and seen at an inclination of 8 degrees from edge-on. This interpretation can be tested by monitoring Spikey for periodic behavior and recurring flares in the next few years. In preparation for such monitoring we present the first X-ray observations of this object taken by the Neil Gehrels Swift observatory., Comment: Accepted to MNRAS

Using a Milky Way double neutron star (DNS) merger rate of 210 Myr$^{-1}$, as derived by the Laser Interferometer Gravitational-Wave Observatory (LIGO), we demonstrate that the Laser Interferometer Space Antenna (LISA) will detect on average 240 (330) DNSs within the Milky Way for a 4-year (8-year) mission with a signal-to-noise ratio greater than 7. Even adopting a more pessimistic rate of 42 Myr$^{-1}$, as derived by the population of Galactic DNSs, we find a significant detection of 46 (65) Milky Way DNSs. These DNSs can be leveraged to constrain formation scenarios. In particular, traditional NS-discovery methods using radio telescopes are unable to detect DNSs with $P_{\rm orb}$ $\lesssim$1 hour (merger times $\lesssim$10 Myr). If a fast-merging channel exists that forms DNSs at these short orbital periods, LISA affords, perhaps, the only opportunity to observationally characterize these systems; we show that toy models for possible formation scenarios leave unique imprints on DNS orbital eccentricities, which may be measured by LISA for values as small as $\sim$10$^{-2}$., 10 pages, 4 figures, submitted to ApJL

We show that nearly half of all binary black hole (BBH) mergers dynamically assembled in globular clusters have measurable eccentricities ($e>0.01$) in the LISA band ($10^{-2}$ Hz), when General Relativistic corrections are properly included in the $N$-body evolution. If only Newtonian gravity is included, the derived fraction of eccentric LISA sources is significantly lower, which explains why recent studies all have greatly underestimated this fraction. Our findings have major implications for how to observationally distinguish between BBH formation channels using eccentricity with LISA, which is one of the key science goals of the mission. We illustrate that the relatively large population of eccentric LISA sources reported here originates from BBHs that merge between hardening binary-single interactions inside their globular cluster. These results indicate a bright future for using LISA to probe the origin of BBH mergers., accepted version. MNRAS

Using numerical hydrodynamics calculations and a novel method for densely sampling parameter space, we measure the accretion and torque on a binary system from a circumbinary disk. In agreement with previous studies, we find that the net torque on the binary is positive for mass ratios close to unity, and that accretion always drives the binary towards equal mass. Accretion variability depends sensitively on the numerical sink prescription, but the torque and relative accretion onto each component do not depend on the sink timescale. Positive torque and highly variable accretion occurs only for mass ratios greater than around $0.05$. This means that for mass ratios below $0.05$, the binary would migrate inward until the secondary accreted sufficient mass, after which it would execute a U-turn and migrate outward. We explore a range of viscosities, from $\alpha = 0.03$ to $\alpha = 0.15$, and find that this outward torque is proportional to the viscous torque, simply proportional to viscosity in this range. Dependence of accretion and torque on mass ratio is explored in detail, densely sampling mass ratios between $0.01$ and unity. For mass ratio $q > 0.6$, accretion variability is found to exhibit a distinct sawtooth pattern, typically with a five-orbit cycle that provides a "smoking gun" prediction for variable quasars observed over long periods, as a potential means to confirm the presence of a binary., Comment: ApJ Accepted

The optical and UV emission from sub-parsec massive black hole binaries (MBHBs) in active galactic nuclei (AGN) is believed to vary periodically, on timescales comparable to the binary's orbital time. If driven by accretion rate fluctuations, the variability could be isotropic. If dominated by relativistic Doppler modulation, the variability should instead be anisotropic, resembling a rotating forward-beamed lighthouse. We consider the infrared (IR) reverberation of either type of periodic emission by pc-scale circumbinary dust tori. We predict the phase and amplitude of IR variability as a function of the ratio of dust light crossing time to the source variability period, and of the torus inclination and opening angle. We enumerate several differences between the isotropic and anisotropic cases. Interestingly, for a nearly face-on binary with an inclined dust torus, the Doppler boost can produce IR variability without any observable optical/UV variability. Such orphan-IR variability would have been missed in optical searches for periodic AGN. We apply our models to time-domain WISE IR data from the MBHB candidate PG 1302-102 and find consistency with dust reverberation by both isotropically emitting and Doppler-boosted sources in the shorter wavelength W1-W2 (2.8-5.3 micrometer) bands. We constrain the dust torus to be thin (aspect ratio ~0.1), with an inner radius at 1-5 pc. More generally, our dust echo models will aid in identifying new MBHB candidates, determining their nature, and constraining the physical properties of MBHBs and their dust tori., We dedicate this work to the memory of Arlin Crotts, an expert on supernova light echoes who passed away on November 19, 2015

We present and analyse a new tidal disruption event (TDE), AT2017eqx at redshift z=0.1089, discovered by Pan-STARRS and ATLAS. The position of the transient is consistent with the nucleus of its host galaxy; it peaks at a luminosity of $L \approx 10^{44}$ erg s$^{-1}$; and the spectrum shows a persistent blackbody temperature $T \gtrsim 20,000$ K with broad H I and He II emission. The lines are initially centered at zero velocity, but by 100 days the H I lines disappear while the He II develops a blueshift of $\gtrsim 5,000$ km s$^{-1}$. Both the early- and late-time morphologies have been seen in other TDEs, but the complete transition between them is unprecedented. The evolution can be explained by combining an extended atmosphere, undergoing slow contraction, with a wind in the polar direction becoming visible at late times. Our observations confirm that a lack of hydrogen a TDE spectrum does not indicate a stripped star, while the proposed model implies that much of the diversity in TDEs may be due to the observer viewing angle. Modelling the light curve suggests AT2017eqx resulted from the complete disruption of a solar-mass star by a black hole of $\sim 10^{6.3} M_\odot$. The host is another quiescent, Balmer-strong galaxy, though fainter and less centrally concentrated than most TDE hosts. Radio limits rule out a relativistic jet, while X-ray limits at 500 days are among the deepest for a TDE at this phase., Comment: Updated to match published version

We explore a unique electromagnetic signature of stellar-mass compact-object binaries long before they are detectable in gravitational waves. We show that gravitational lensing of light emitting components of a compact-object binary, by the other binary component, could be detectable in the nearby universe. This periodic lensing signature could be detected from present and future X-ray observations, identifying the progenitors of binaries that merge in the LIGO band, and also unveiling populations that do not merge, thus providing a tracer of the compact-object binary population in an enigmatic portion of its life. We argue that periodically repeating lensing flares could be observed for ~100 ksec orbital-period binaries with the future Lynx X-ray mission, possibly concurrent with gravitational wave emission in the LISA band. Binaries with longer orbital periods could be more common and be detectable as single lensing flares, though with reliance on a model for the flare that can be tested by observations of succeeding flares. Non-detection of such events, even with existing X-ray observations, will help to constrain the population of EM bright compact-object binaries., Published in MNRAS

The bright quasar PG1302-102 has been identified as a candidate supermassive black hole binary from its near-sinusoidal optical variability. While the significance of its optical periodicity has been debated due to the stochastic variability of quasars, its multi-wavelength variability in the ultraviolet (UV) and optical bands is consistent with relativistic Doppler boost caused by the orbital motion in a binary. However, this conclusion was based previously on sparse UV data which were not taken simultaneously with the optical data. Here we report simultaneous follow-up observations of PG1302-102 with the Ultraviolet Optical Telescope on the Neil Gehrels Swift Observatory in six optical + UV bands. The additional nine Swift observations produce light curves roughly consistent with the trend under the Doppler boost hypothesis, which predicts that UV variability should track the optical, but with a ~2.2 times higher amplitude. We perform a statistical analysis to quantitatively test this hypothesis. We find that the data are consistent with the Doppler boost hypothesis when we compare the the amplitudes in optical B-band and UV light curves. However, the ratio of UV to V-band variability is larger than expected and is consistent with the Doppler model, only if either the UV/optical spectral slopes vary, the stochastic variability makes a large contribution in the UV, or the sparse new optical data underestimate the true optical variability. We have evidence for the latter from comparison with the optical light curve from ASAS-SN. Additionally, the simultaneous analysis of all four bands strongly disfavors the Doppler boost model whenever Swift V-band is involved. Additional, simultaneous optical + UV observations tracing out another cycle of the 5.2-year proposed periodicity should lead to a definitive conclusion., Comment: accepted in MNRAS

We study the formation rate of binary black hole mergers formed through gravitational-wave emission between unbound, single black holes in globular clusters. While the formation of these binaries in very dense systems such as galactic nuclei has been well studied, we show here that this process can operate in lower-density stellar systems as well, forming binaries at a rate similar to other proposed pathways for creating eccentric mergers. Recent advances in post-Newtonian cluster dynamics indicate that a large fraction of dynamically-assembled binary black holes merge inside their host clusters during weak and strong binary-single and binary-binary interactions, and that these systems may retain measurable eccentricities as they travel through the LIGO and LISA sensitivity bands. Using an analytic approach to modeling binary black holes from globular clusters, we show that the formation of merging binaries from previously unbound black holes can operate at a similar rate to mergers forming during strong binary encounters, and that these binaries inhabit a unique region of the gravitational-wave frequency space which can be identified by proposed deci-Hertz space-based detectors., Comment: 16 pages, 5 figures. Accepted version (PRD)

We consider binary black holes (BBHs) in a hierarchical triple system where a more compact, less-massive binary is emitting detectable gravitational waves (GWs), and the tertiary is a supermassive BH at the center of a nuclear star cluster. As previous works have shown, the orbital motion of the outer binary can generate a detectable relativistic Doppler boost of the GWs emitted by the orbiting inner binary. We show here that for outer-binary orbits with a period of order one year, there can be a non-negligible probability for repeated gravitational lensing of the GWs emitted by the inner binary. Repeating gravitational lensing events could be detected by the LISA observatory as periodic GW amplitude spikes before the BBH enters the LIGO band. Such a detection would confirm the origin of some BBH mergers in nuclear star clusters. GW lensing also offers new testing grounds for strong gravity., Comment: Accepted to Physical Review D. Comments welcome

We present analysis of Chandra X-ray observations of seven quasars that were identified as candidate subparsec binary supermassive black hole (SMBH) systems in the Catalina Real-Time Transient Survey based on the apparent periodicity in their optical light curves. Simulations predict that close-separation accreting SMBH binaries will have different X-ray spectra than single accreting SMBHs, including harder or softer X-ray spectra, ripple-like profiles in the Fe K-α line, and distinct peaks in the spectrum due to the separation of the accretion disk into a circumbinary disk and mini disks around each SMBH. We obtained Chandra observations to test these models and assess whether these quasars could contain binary SMBHs. We instead find that the quasar spectra are all well fit by simple absorbed power-law models, with the rest-frame 2–10 keV photon indices, Γ, and the X-ray-to-optical power slopes, α_(OX), indistinguishable from those of the larger quasar population. This may indicate that these seven quasars are not truly subparsec binary SMBH systems, or it may simply reflect that our sample size was too small to robustly detect any differences. Alternatively, the X-ray spectral changes might only be evident at energies higher than probed by Chandra. Given the available models and current data, no firm conclusions are drawn. These observations will help motivate and direct further work on theoretical models of binary SMBH systems, such as modeling systems with thinner accretion disks and larger binary separations.

Among the growing list of confirmed exoplanets, the number of planets identified in dense star clusters remains sparse. Previous analyses have suggested this may be due in part to dynamical interactions that unbind large fractions of planets from their host stars, limiting the survival of planetary systems in clusters. Thus, alternative detection strategies may be necessary to study planets in clusters that may no longer be bound to a host star. Here, we use the cluster Monte Carlo code CMC to explore the evolution of planetary systems in dense star clusters. Depending on a number of initial conditions, we show that $10-50\%$ of primordial planetary systems are broken through dynamical encounters over a cluster's full lifetime, populating clusters with "free-floating" planets. Furthermore, a large number ($30-80\%$) of planets are ejected from their host cluster through strong dynamical encounters and/or tidal loss. Additionally, we show that planets naturally mix with stellar-mass black holes (BHs) in the central regions of their host cluster. As a consequence, up to a few hundreds of planets will be tidally disrupted through close passages of BHs. We show these BH-planet tidal disruption events (TDEs) occur in clusters at a rate of up to $10^{-5}\,\rm{yr}^{-1}$ in a Milky Way-type galaxy. We predict BH-planet TDEs should be detected by upcoming transient surveys such as LSST at a rate of roughly a few events per year. The observed rate of BH-planet TDEs would place new constraints upon the formation and survival of both planetary systems and BHs in dense star clusters. Additionally, depending on various assumptions including the initial number of planets and their orbital properties, we predict that typical globular clusters may contain a few dynamically-formed NS-planet systems at present as well as up to roughly 100 dynamically-formed WD-planet systems., 18 pages, 8 figures. Submitted to ApJ. Comments welcome

We derive the observable gravitational wave (GW) peak frequency ($f$) distribution of binary black holes (BBHs) that currently reside inside their globular clusters (GCs), with and without 2.5 Post-Newtonian (2.5PN) effects included in the dynamical evolution of the BBHs. Recent Newtonian studies have reported that a notable number of nearby non-merging BBHs, i.e. those BBHs that are expected to undergo further dynamical interactions before merger, in GCs are likely to be observable by LISA. However, our 2.5PN calculations show that the distribution of $\log f$ for the non-merging BBH population above $\sim 10^{-3.5}$ Hz scales as $f^{-34/9}$ instead of the $f^{-2/3}$ scaling found in the Newtonian case. This leads to an approximately two-orders-of-magnitude reduction in the expected number of GW sources at $\sim 10^{-3}$ Hz, which lead us to conclude that observing nearby BBHs with LISA is not as likely as has been claimed in the recent literature. In fact, our results suggest that it might be more likely that LISA detects the population of BBHs that will merge before undergoing further interactions. This interestingly suggests that the BBH merger rate derived from LIGO can be used to forecast the number of nearby LISA sources, as well as providing insight into the fraction of BBH mergers forming in GCs., Comment: 9 pages. 1 figure

We show that a 10 year Gaia mission could astrometrically detect the orbital motion of ~1 sub-parsec separation supermassive black hole binary in the heart of nearby, bright active galactic nuclei (AGN). Candidate AGN lie out to a redshift of z=0.02 and in the V-band magnitude range $10 \lesssim m_V \lesssim 13$. The distribution of detectable binary masses peaks at a few times $10^7$ solar masses and is truncated above a few times $10^8$ solar masses., Comment: Published in Phys. Rev. D

Supermassive black hole binaries (SMBHBs) should be common in galactic nuclei as a result of frequent galaxy mergers. Recently, a large sample of sub-parsec SMBHB candidates was identified as bright periodically variable quasars in optical surveys. If the observed periodicity corresponds to the redshifted binary orbital period, the inferred orbital velocities are relativistic (v/c~0.1). The optical and UV luminosities are expected to arise from gas bound to the individual BHs, and would be modulated by the relativistic Doppler effect. The optical and UV light curves should vary in tandem with relative amplitudes which depend on the respective spectral slopes. We constructed a control sample of 42 quasars with aperiodic variability, to test whether this Doppler colour signature can be distinguished from intrinsic chromatic variability. We found that the Doppler signature can arise by chance in ~20% (~37%) of quasars in the nUV (fUV) band. These probabilities reflect the limited quality of the control sample and represent upper limits on how frequently quasars mimic the Doppler brightness+colour variations. We performed separate tests on the periodic quasar candidates, and found that for the majority, the Doppler boost hypothesis requires an unusually steep UV spectrum or an unexpectedly large BH mass and orbital velocity. We conclude that at most ~1/3rd of these periodic candidates can harbor Doppler-modulated SMBHBs., Comment: submitted in MNRAS

In anticipation of a LIGO detection of a black hole/neutron star merger, we expand on the intriguing possibility of an electromagnetic counterpart. Black hole/Neutron star mergers could be disappointingly dark since most black holes will be large enough to swallow a neutron star whole, without tidal disruption and without the subsequent fireworks. Encouragingly, we previously found a promising source of luminosity since the black hole and the highly-magnetized neutron star establish an electronic circuit -- a black hole battery. In this paper, arguing against common lore, we consider the electric charge of the black hole as an overlooked source of electromagnetic radiation. Relying on the well known Wald mechanism by which a spinning black hole immersed in an external magnetic field acquires a stable net charge, we show that a strongly-magnetized neutron star in such a binary system will give rise to a large enough charge in the black hole to allow for potentially observable effects. Although the maximum charge is stable, we show there is a continuous flux of charges contributing to the luminosity. Most interestingly, the spinning charged black hole then creates its own magnetic dipole to power a black hole pulsar., Comment: Published in Physical Review D

In paper I of this series we showed that a large percentage of the binary black hole (BBH) mergers that form through dynamical interactions in globular clusters will have significant eccentricity in the ~10^{-3}-10^{-1} Hz LISA band. In this work we quantify the evolution of these highly eccentric binaries through the LISA and LIGO bands, and compute the stochastic gravitational wave background from the merging, eccentric population. We find that the population of BBHs that merge in-between three-body encounters inside their cluster (~50% of all cluster-formed BBH mergers) will have measurable eccentricity for their entire lifetime in the LISA band. The population of BBHs that merge during three-body encounters (~5% of all cluster-formed BBH mergers), will be detectable by LIGO with eccentricities of e~0.1. The gravitational wave background from dynamically assembled BBHs encodes a characteristic bump due to the high initial eccentricities of these systems. The location and amplitude of this bump depends on globular cluster properties., Comment: Published in MNRAS. Companion paper to arxiv:1804.06519

The rate of tidal disruption flares (TDFs) per mass of the disrupting black hole encodes information on the present-day mass function (PDMF) of stars in the clusters surrounding super massive black holes. We explore how the shape of the TDF rate with black hole mass can constrain the PDMF, with only weak dependence on black hole spin. We show that existing data can marginally constrain the minimum and maximum masses of stars in the cluster, and the high-mass end of the PDMF slope, as well as the overall TDF rate. With of order 100 TDFs expected to be identified with the Zwicky Transient Facility, the overall rate can be highly constrained, but still with only marginal constraints on the PDMF. However, if less than of order 10% of the TDFs expected to be found by LSST over a decade (of order 1000 TDFs) are identified, then precise and accurate estimates can be made for the minimum stellar mass (within a factor of two) and the average slope of the high-mass PDMF (to within ~10%) in nuclear star clusters. This technique could be adapted in the future to probe, in addition to the PDMF, the local black hole mass function and possibly the massive black hole binary population., Comment: Accepted for publication in MNRAS

The coalescence of a compact object with a $10^{4}-10^{7} {\rm M_\odot}$ supermassive black hole (SMBH) produces mHz gravitational waves (GWs) detectable by the future Laser Interferometer Space Antenna (LISA). If such an inspiral occurs in the accretion disc of an active galactic nucleus (AGN), the gas torques imprint a small deviation in the GW waveform. Here we present two-dimensional hydrodynamical simulations with the moving-mesh code DISCO of a BH inspiraling at the GW rate in a binary system with a mass ratio $q\!=\!M_2/M_1\!=\!10^{-3}$, embedded in an accretion disc. We assume a locally isothermal equation of state for the gas (with Mach number $\mathcal{M}=20$) and implement a standard $\alpha$-prescription for its viscosity (with $\alpha = 0.03$). We find disc torques on the binary that are weaker than in previous semi-analytic toy models, and are in the opposite direction: the gas disc slows down, rather than speeds up the inspiral. We compute the resulting deviations in the GW waveform, which scale linearly with the mass of the disc. The SNR of these deviations accumulates mostly at high frequencies, and becomes detectable in a 5-year LISA observation if the total phase shift exceeds a few radians. We find that this occurs if the disc surface density exceeds $\Sigma_0 \gtrsim 10^{2-3}\rm g\,cm^{-2}$, as may be the case in thin discs with near-Eddington accretion rates. Since the characteristic imprint on the GW signal is strongly dependent on disc parameters, a LISA detection of an intermediate mass ratio inspiral would probe the physics of AGN discs and migration., Comment: Incorporates published erratum - revised figures in Section 5, results unchanged

Very long baseline interferometry (VLBI) at millimeter (mm) wavelengths is being employed to resolve event-horizon scale structure of the environment surrounding the Milky-Way black hole, at an angular resolution of a few tens of micro-arcseconds. The same approach could also resolve the orbital separation of a population of massive black hole binaries (MBHBs). Modeling the inspiral of binaries due to gravitational wave emission and gas and requiring binary orbital periods of less than 10 years, we estimate that there may exist ~100 resolvable MBHBs that are bright enough to be observed by mm-wavelength VLBI instruments over the entire sky, at redshifts z, Published in The Astrophysical Journal

Nearly 150 massive black hole binary (MBHB) candidates at sub-pc orbital separations have been reported in recent literature. Nevertheless, the definitive detection of even a single such object remains elusive. If at least one of the black holes is accreting, the light emitted from its accretion disc will be lensed by the other black hole for binary orbital inclinations near to the line of sight. This binary self-lensing could provide a unique signature of compact MBHB systems. We show that, for MBHBs with masses in the range 10^6 - 10^10 solar masses and with orbital periods less than ~10 yr, strong lensing events should occur in one to 10s of percent of MBHB systems that are monitored for an entire orbit. Lensing events will last from days for the less massive, shorter period MBHBs up to a year for the most massive ~10 year orbital period MBHBs. At small inclinations of the binary orbit to the line of sight, lensing must occur and will be accompanied by periodicity due to the relativistic Doppler boost. Flares at the same phase as the otherwise average flux of the Doppler modulation would be a smoking gun signature of self-lensing and can be used to constrain binary parameters. For MBHBs with separation >100 Schwarzschild radii, we show that finite-sized source effects could serve as a probe of MBH accretion disc structure. Finally, we stress that our lensing probability estimate implies that ~10 of the known MBHB candidates identified through quasar periodicity should exhibit strong lensing flares., Accepted to MNRAS

Because most large galaxies contain a central black hole, and galaxies often merge, black-hole binaries are expected to be common in galactic nuclei. Although they cannot be imaged, periodicities in the light curves of quasars have been interpreted as evidence for binaries, most recently in PG~1302-102, with a short rest-frame optical period of 4 years. If the orbital period matches this value, then for the range of estimated black hole masses the components would be separated by 0.007-0.017 pc, implying relativistic orbital speeds. There has been much debate over whether black hole orbits could be smaller than 1 pc. Here we show that the amplitude and the sinusoid-like shape of the variability of PG~1302-102 can be fit by relativistic Doppler boosting of emission from a compact, steadily accreting, unequal-mass binary. We predict that brightness variations in the ultraviolet light curve track those in the optical, but with a 2-3 times larger amplitude. This prediction is relatively insensitive to the details of the emission process, and is consistent with archival UV data. Follow-up UV and optical observations in the next few years can test this prediction and confirm the existence of a binary black hole in the relativistic regime., To appear as a Letter in the September 17, 2015 issue of Nature. Press embargo until 18:00 London time / 13:00 US Eastern Time on 16 September 2015

Direct detection of black hole-neutron star pairs is anticipated with the advent of aLIGO. Electromagnetic counterparts may be crucial for a confident gravitational-wave detection as well as for extraction of astronomical information. Yet black hole-neutron star pairs are notoriously dark and so inaccessible to telescopes. Contrary to this expectation, a bright electromagnetic transient can occur in the final moments before merger as long as the neutron star is highly magnetized. The orbital motion of the neutron star magnet creates a Faraday flux and corresponding power available for luminosity. A spectrum of curvature radiation ramps up until the rapid injection of energy ignites a fireball, which would appear as an energetic blackbody peaking in the x ray to $\gamma$ rays for neutron star field strengths ranging from $10^{12}$G to $10^{16}$G respectively and a $10M_{\odot} $ black hole. The fireball event may last from a few milliseconds to a few seconds depending on the neutron star magnetic-field strength, and may be observable with the Fermi Gamma-Ray Burst Monitor with a rate up to a few per year for neutron star field strengths $\gtrsim 10^{14}$G. We also discuss a possible decaying post-merger event which could accompany this signal. As an electromagnetic counterpart to these otherwise dark pairs, the black-hole battery should be of great value to the development of multi-messenger astronomy in the era of aLIGO., Comment: Published in Physical Review D

We study circumbinary accretion discs in the framework of the restricted three-body problem (R3Bp) and via numerically solving the height-integrated equations of viscous hydrodynamics. Varying the mass ratio of the binary, we find a pronounced change in the behaviour of the disc near mass ratio $q \equiv M_s/M_p \sim 0.04$. For mass ratios above $q=0.04$, solutions for the hydrodynamic flow transition from steady, to strongly-fluctuating; a narrow annular gap in the surface density around the secondary's orbit changes to a hollow central cavity; and a spatial symmetry is lost, resulting in a lopsided disc. This phase transition is coincident with the mass ratio above which stable orbits do not exist around the L4 and L5 equilibrium points of the R3B problem. Using the DISCO code, we find that for thin discs, for which a gap or cavity can remain open, the mass ratio of the transition is relatively insensitive to disc viscosity and pressure. The $q=0.04$ transition has relevance for the evolution of massive black hole binary+disc systems at the centers of galactic nuclei, as well as for young stellar binaries and possibly planets around brown dwarfs., published in MNRAS

It is well known that the classical gravitational two body problem can be transformed into a spherical harmonic oscillator by regularization. We find that a modification of the regularization transformation has a similar result to leading order in general relativity. In the resulting harmonic oscillator, the leading-order relativistic perturbation is formally a negative centrifugal force. The net centrifugal force changes sign at three Schwarzschild radii, which interestingly mimics the innermost stable circular orbit (ISCO) of the full Schwarzschild problem. Transforming the harmonic-oscillator solution back to spatial coordinates yields, for both timelike and null weak-field Schwarzschild geodesics, a solution for $t,r,\phi$ in terms of elementary functions of a variable that can be interpreted as a generalized eccentric anomaly. The textbook expressions for relativistic precession and light deflection are easily recovered. We suggest how this solution could be combined with additional perturbations into numerical methods suitable for applications such as relativistic accretion or dynamics of the Galactic-centre stars.

We describe an undergraduate laboratory that combines an accurate measurement of the speed of light, a fundamental investigation of a basic laser system, and a nontrivial use of statistical analysis. Students grapple with the existence of longitudinal modes in a laser cavity as they change the cavity length of an adjustable-cavity HeNe laser and tune the cavity to produce lasing in the TEM$_{00}$ mode. For appropriate laser cavity lengths, the laser gain curve of a HeNe laser allows simultaneous operation of multiple longitudinal modes. The difference frequency between the modes is measured using a self-heterodyne detection with a diode photodetector and a radio frequency spectrum analyzer. Asymmetric effects due to frequency pushing and frequency pulling, as well as transverse modes, are minimized by simultaneously monitoring and adjusting the mode structure as viewed with a Fabry-Perot interferometer. The frequency spacing of longitudinal modes is proportional to the inverse of the cavity length with a proportionality constant equal to half the speed of light. By changing the length of the cavity, without changing the path length within the HeNe gas, the speed of light in air can be measured to be ($2.9972 \pm0.0002) \times 10^{8}$ m/s, which is to high enough precision to distinguish between the speed of light in air and that in a vacuum., Comment: 6 pages; 6 figures; Published in the American Journal of Physics

Graham et al. (2015) have detected a 5.2 year periodic optical variability of the quasar PG 1302-102 at redshift $z=0.3$, which they interpret as the redshifted orbital period $(1+z)t_{\rm bin}$ of a putative supermassive black hole binary (SMBHB). Here we consider the implications of a $3-8$ times shorter orbital period, suggested by hydrodynamical simulations of circumbinary discs (CBDs) with nearly equal--mass SMBHBs ($q\equiv M_2/M_1\gtrsim 0.3$). With the corresponding $2-4$ times tighter binary separation, PG 1302 would be undergoing gravitational wave dominated inspiral, and serve as a proof that the BHs can be fueled and produce bright emission even in this late stage of the merger. The expected fraction of binaries with the shorter $t_{\rm bin}$, among bright quasars, would be reduced by 1-2 orders of magnitude, compared to the 5.2 year period, in better agreement with the rarity of candidates reported by Graham et al. (2015). Finally, shorter periods would imply higher binary speeds, possibly imprinting periodicity on the light curves from relativistic beaming, as well as measurable relativistic effects on the Fe K$\alpha$ line. The CBD model predicts additional periodic variability on timescales of $t_{\rm bin}$ and $\approx 0.5 t_{\rm bin}$, as well as periodic variation of broad line widths and offsets relative to the narrow lines, which are consistent with the observations. Future observations will be able to test these predictions and hence the binary+CBD hypothesis for PG 1302., Comment: Published in MNRAS

Most standard descriptions of Type II migration state that massive, gap-opening planets must migrate at the viscous drift rate. This is based on the idea that the disk is separated into an inner and outer region and gas is considered unable to cross the gap. In fact, gas easily crosses the gap on horseshoe orbits, nullifying this necessary premise which would set the migration rate. In this work, it is demonstrated using highly accurate numerical calculations that the actual migration rate is dependent on disk and planet parameters, and can be significantly larger or smaller than the viscous drift rate. In the limiting case of a disk much more massive than the secondary, the migration rate saturates to a constant which is sensitive to disk parameters and is not necessarily of order viscous rate. In the opposite limit of a low-mass disk, the migration rate decreases linearly with disk mass. Steady-state solutions in the low disk mass limit show no pile-up outside the secondary's orbit, and no corresponding drainage of the inner disk., Comment: ApJ Accepted

As a black hole and neutron star approach during inspiral, the field lines of a magnetized neutron star eventually thread the black hole event horizon and a short-lived electromagnetic circuit is established. The black hole acts as a battery that provides power to the circuit, thereby lighting up the pair just before merger. Although originally suggested as a promising electromagnetic counterpart to gravitational-wave detection, the luminous signals are promising more generally as potentially detectable phenomena, such as short gamma-ray bursts. To aid in the theoretical understanding, we present analytic solutions for the electromagnetic fields of a magnetic dipole in the presence of an event horizon. In the limit that the neutron star is very close to a Schwarzschild horizon, the Rindler limit, we can solve Maxwell's equations exactly for a magnetic dipole on an arbitrary worldline. We present these solutions here and investigate a proxy for a small segment of the neutron star orbit around a big black hole. We find that the voltage the black hole battery can provide is in the range ~10^16 statvolts with a projected luminosity of 10^42 ergs/s for an M=10M_sun black hole, a neutron star with a B-field of 10^12 G, and an orbital velocity ~0.5c at a distance of 3M from the horizon. Larger black holes provide less power for binary separations at a fixed number of gravitational radii. The black hole/neutron star system therefore has a significant power supply to light up various elements in the circuit possibly powering jets, beamed radiation, or even a hot spot on the neutron star crust., Published in Physical Review D: http://link.aps.org/doi/10.1103/PhysRevD.88.064059