Search

We present the results of a multi-wavelength follow up campaign for the luminous nuclear transient Gaia16aax, which was first identified in January 2016. The transient is spatially consistent with the nucleus of an active galaxy at z=0.25, hosting a black hole of mass $\rm \sim6\times10^8M_\odot$. The nucleus brightened by more than 1 magnitude in the Gaia G-band over a timescale of less than one year, before fading back to its pre-outburst state over the following three years. The optical spectra of the source show broad Balmer lines similar to the ones present in a pre-outburst spectrum. During the outburst, the $\rm H\alpha$ and $\rm H\beta$ emission lines develop a secondary peak. We also report on the discovery of two transients with similar light curve evolution and spectra: Gaia16aka and Gaia16ajq. We consider possible scenarios to explain the observed outbursts. We exclude that the transient event could be caused by a microlensing event, variable dust absorption or a tidal encounter between a neutron star and a stellar mass black hole in the accretion disk. We consider variability in the accretion flow in the inner part of the disk, or a tidal disruption event of a star $\geq 1 M_{\odot}$ by a rapidly spinning supermassive black hole as the most plausible scenarios. We note that the similarity between the light curves of the three Gaia transients may be a function of the Gaia alerts selection criteria., Comment: 21 pages, 17 figure - accepted for publication in MNRAS main journal

We use well-established observational evidence to draw conclusions about the fundamental nature of the viscosity in accretion discs. To do this, we first summarise the observational evidence for the value of the dimensionless accretion disc viscosity parameter $\alpha$, defined by Shakura & Sunyaev (1973, 1976). We find that, for fully ionized discs, the value of $\alpha$ is readily amenable to reliable estimation and that the observations are consistent with the hypothesis that $\alpha \sim 0.2 - 0.3$. In contrast in discs that are not fully ionized, estimates of the value of $\alpha$ are generally less direct and the values obtained are generally $ < 0.01$ and often $ \ll 0.01$. We conclude that this gives us crucial information about the nature of viscosity in accretion discs. First, in fully ionized discs the strength of the turbulence is always limited by being at most trans-sonic. This implies that it is necessary that credible models of the turbulence reflect this fact. Second, the smaller values of $\alpha$ found for less ionized, and therefore less strongly conducting, discs imply that magnetism plays a dominant role. This provides important observational support for the concept of magneto-rotational instability (MRI) driven hydromagnetic turbulence., Comment: Accepted for publication in New Astronomy

We present long-duration numerical simulations of the tidal disruption of stars modelled with accurate stellar structures and spanning a range of pericentre distances, corresponding to cases where the stars are partially and completely disrupted. We substantiate the prediction that the late-time power-law index of the fallback rate $n_{\infty} \simeq -5/3$ for full disruptions, while for partial disruptions---in which the central part of the star survives the tidal encounter intact---we show that $n_{\infty} \simeq -9/4$. For the subset of simulations where the pericenter distance is close to that which delineates full from partial disruption, we find that a stellar core can reform after the star has been completely destroyed; for these events the energy of the zombie core is slightly positive, which results in late-time evolution from $n \simeq -9/4$ to $n \simeq -5/3$. We find that self-gravity can generate an $n(t)$ that deviates from $n_{\infty}$ by a small but significant amount for several years post-disruption. In one specific case with the stellar pericenter near the critical value, we find self-gravity also drives the re-collapse of the central regions of the debris stream into a collection of several cores while the rest of the stream remains relatively smooth. We also show that it is possible for the surviving stellar core in a partial disruption to acquire a circumstellar disc that is shed from the rapidly rotating core. Finally, we provide a novel analytical fitting function for the fallback rates that may also be useful in a range of contexts beyond TDEs., Comment: 21 pages, 12 figures, accepted for publication in ApJ

Synchrotron-emitting, nonthermal filaments (NTFs) have been observed near the Galactic center for nearly four decades, yet their physical origin remains unclear. Here we investigate the possibility that NTFs are produced by the destruction of molecular clouds by the gravitational potential of the Galactic center. We show that this model predicts the formation of a filamentary structure with length on the order of tens to hundreds of pc, a highly ordered magnetic field along the axis of the filament, and conditions conducive to magnetic reconnection that result in particle acceleration. This model therefore yields the observed magnetic properties of NTFs and a population of relativistic electrons, without the need to appeal to a dipolar, $\sim$ mG, Galactic magnetic field. As the clouds can be both completely or partially disrupted, this model provides a means of establishing the connection between filamentary structures and molecular clouds that is observed in some, but not all, cases., Figure added, includes referee suggestions

In order to provide an explanation for the unexpected radial brightness distribution of the steady accretion discs seen in nova-like variables, Nixon & Pringle (2019) proposed that the accretion energy is redistributed outwards by means of strong, magnetically driven, surface flows. In this paper we note that the "powerful, rotating disc winds" observed in the soft states of black hole X-ray binaries, and also in the disc around a magnetized neutron star in Her X-1, have the properties of the outflows postulated by Nixon & Pringle to exist in the nova-like variable accretion discs around white dwarfs. The relevant properties are that the flows are not winds, but are, instead, bound flows (traveling at less than the escape velocity) and that the mass fluxes in the flows are a substantial fraction of the accretion rate in the disc., 4 pages, 2 figures, accepted for publication in Astronomy & Astrophysics

Many accretion discs are thought to be warped. Recent hydrodynamical simulations show that (i) discs can break into distinct planes when the amplitude of an imposed warp is sufficiently high and the viscosity sufficiently low, and that (ii) discs can tear up into discrete rings when an initially planar disc is subject to a forced precession. Previously, we investigated the local stability of isolated, Keplerian, warped discs in order to understand the physics causing an accretion disc to break into distinct planes, finding that antidiffusion of the warp amplitude is the underlying cause. Here, we explore the behaviour of this instability in disc regions where the rotation profile deviates from Keplerian. We find that at small warp amplitudes non-Keplerian rotation can stabilize the disc by increasing the critical warp amplitude for instability, while at large warp amplitudes non-Keplerian rotation can lead to an increased growth rate for discs that are unstable. Tidal effects on discs in binary systems are typically weak enough such that the disc remains close to Keplerian rotation. However, the inner regions of discs around black holes are strongly affected, with the smallest radius at which the disc can break into discrete planes being a function of the black hole spin. We suggest that interpreting observed frequencies in the power spectra of light curves from accreting compact objects as nodal and apsidal precession of discrete orbits requires an instability that can break the disc into discrete rings such as the one explored here. © 2020 The Author(s)., Science and Technology Facilities Council, STFC: ST/M005917/1 Türkiye Bilimsel ve Teknolojik Araştirma Kurumu, TÜBITAK: TÜB?TAK – 117F280, We thank the referee for useful comments. We thank Gordon Ogilvie for providing the code to calculate the torque coefficients. We thank Jim Pringle and Eric Coughlin for useful discussions. SD acknowledges the warm hospitality of the Department of Physics and Astronomy at University of Leicester during her visit. SD is supported by The Scientific and Technological Research Council of Turkey (TÜB?TAK – 117F280). CJN is supported by the Science and Technology Facilities Council (grant number ST/M005917/1).

Be stars are rapidly rotating B stars with Balmer emission lines that indicate the presence of a Keplerian, rotationally supported, circumstellar gas disc. Current disc models, referred to as "decretion discs", make use of the zero torque inner boundary condition typically applied to accretion discs, with the 'decretion' modelled by adding mass to the disc at a radius of about two per cent larger than the inner disc boundary. We point out that, in this model, the rates at which mass and energy need to be added to the disc are implausibly large. What is required is that the disc has not only a source of mass but also a continuing source of angular momentum. We argue that the disc evolution may be more physically modelled by application of the non-zero torque inner boundary condition of Nixon & Pringle (2020), which determines the torque applied at the boundary as a fraction of the advected angular momentum flux there and approaches the accretion and decretion disc cases in the appropriate limits. We provide supporting arguments for the suggestion that the origin of the disc material is small-scale magnetic flaring events on the stellar surface, which, when combined with rapid rotation, can provide sufficient mass to form, and sufficient angular momentum to maintain, a Keplerian Be star disc. We discuss the origin of such small-scale magnetic fields in radiative stars with differential rotation. We conclude that small-scale magnetic fields on the stellar surface, may be able to provide the necessary mass flux and the necessary time-dependent torque on the disc inner regions to drive the observed disc evolution., 6 pages, no figures, accepted for publication in The Astrophysical Journal Letters

Short gamma-ray bursts are thought to result from the mergers of two neutron stars or a neutron star and stellar mass black hole. The final stages of the merger are generally accompanied by the production of one or more tidal "tails" of ejecta, which fall back onto the remnant-disc system at late times. Using the results of a linear stability analysis, we show that if the material comprising these tails is modeled as adiabatic and the effective adiabatic index satisfies $\gamma \ge 5/3$, then the tails are gravitationally unstable and collapse to form small-scale knots. We analytically estimate the properties of these knots, including their spacing along the tidal tail and the total number produced, and their effect on the mass return rate to the merger remnant. We perform hydrodynamical simulations of the disruption of a polytropic (with the polytropic and adiabatic indices $\gamma$ equal), $\gamma =2$ neutron star by a black hole, and find agreement between the predictions of the linear stability analysis and the distribution of knots that collapse out of the instability. The return of these knots to the black hole induces variability in the fallback rate, which can manifest as variability in the lightcurve of the GRB and -- depending on how rapidly the instability operates -- the prompt emission. The late-time variability induced by the return of these knots is also consistent with the extended emission observed in some GRBs., Comment: Small corrections, additional references included to reflect ApJL published version

A tidal disruption event (TDE) occurs when a star plunges through a supermassive black hole's tidal radius, at which point the star's self-gravity is overwhelmed by the tidal gravity of the black hole. In a partial TDE, where the star does not reach the full disruption radius, only a fraction of the star's mass is tidally stripped while the rest remains intact in the form of a surviving core. Analytical arguments have recently suggested that the temporal scaling of the fallback rate of debris to the black hole asymptotes to $t^{-9/4}$ for partial disruptions, effectively independently of the mass of the intact core. We present hydrodynamical simulations that verify the existence of this predicted, $t^{-9/4}$ scaling. We also define a break timescale -- the time at which the fallback rate transitions from a $t^{-5/3}$ scaling to the characteristic $t^{-9/4}$ scaling -- and measure this break timescale as a function of the impact parameter and the surviving core mass. These results deepen our understanding of the properties and breadth of possible fallback curves expected from TDEs and will therefore facilitate more accurate interpretation of data from wide-field surveys., Comment: Accepted for publication in The Astrophysical Journal, 13th June 2020. 11 pages, 8 figures

Blueshifted absorption lines in the X-ray spectra of AGN show that ultra-fast outflows with typical velocities $v \sim 0.1c$ are a common feature of these luminous objects. Such powerful AGN winds offer an explanation of the observed M-$\sigma$ relation linking the mass of the supermassive black hole and the velocity dispersion in the galaxy's stellar bulge. An extended XMM-Newton study of the luminous Seyfert galaxy PG1211+143 recently revealed a variable multi-velocity wind. Here we report the detection of a short-lived, ultra-fast inflow during the same observation. Previous reports of inflows used single absorption lines with uncertain identifications, but this new result identifies an array of resonance absorption lines of highly ionised Fe, Ca, Ar, S and Si, sharing a common redshift when compared with a grid of realistic photoionization spectra. The redshifted absorption arises in a column of highly ionized matter close to the black hole, with a line-of-sight velocity, $v \sim 0.3c$, inconsistent with the standard picture of a plane circular accretion disc. This may represent the first direct evidence for chaotic accretion in AGN, where accretion discs are generally misaligned to the black hole spin. For sufficient inclinations, the Lense-Thirring effect can break the discs into discrete rings, which then precess, collide and shock, causing near free-fall of gas towards the black hole. The observed accretion rate for the reported infall is comparable to the hard X-ray luminosity in PG1211+143, suggesting that direct infall may be a significant contributor to inner disc accretion., Comment: 7 pages, 3 figures, accepted for publication in MNRAS

A circumbinary disc around a pair of merging stellar-mass black holes may be shocked and heated during the recoil of the merged hole, causing a near-simultaneous electromagnetic counterpart to the gravitational wave event. The shocks occur around the recoil radius, where the disc orbital velocity is equal to the recoil velocity. The amount of mass present near this radius at the time of the merger is critical in determining how much radiation is released. We explore the evolution of a circumbinary disc in two limits. First, we consider an accretion disc that feels no torque from the binary. The disc does not survive until the merger unless there is a dead zone, a region of low turbulence. Even with the dead zone, the surface density in this case may be small. Second, we consider a disc that feels a strong binary torque that prevents accretion on to the binary. In this case there is significantly more mass in regions of interest at the time of the merger. A dead zone in this disc increases the mass close to the recoil radius. For typical binary-disc parameters we expect accretion to be significantly slowed by the resonant torque from the binary, and for a dead zone to be present. We conclude that provided significant mass orbits the binary after the formation of the black hole binary and that the radiation produced in recoil shocks can escape the flow efficiently, there is likely to be an observable electromagnetic signal from black hole binary mergers., Comment: Accepted for publication in MNRAS

We investigate the effects of gradual heating on the evolution of turbulent molecular clouds of mass $2\times 10^6$ M$_\odot$ and virial parameters ranging between $0.7-1.2$. This gradual heating represents the energy output from processes such as winds from massive stars or feedback from High Mass X-ray binaries (HMXBs), contrasting the impulsive energy injection from supernovae (SNe). For stars with a mass high enough that their lifetime is shorter than the life of the cloud, we include a SN feedback prescription. Including both effects, we investigate the interplay between slow and fast forms of feedback and their effectiveness at triggering/suppressing star formation. We find that SN feedback can carve low density chimneys in the gas, offering a path of least resistance for the energy to escape. Once this occurs the more stable, but less energetic, gradual feedback is able to keep the chimneys open. By funneling the hot destructive gas away from the centre of the cloud, chimneys can have a positive effect on both the efficiency and duration of star formation. Moreover, the critical factor is the number of high mass stars and SNe (and any subsequent HMXBs) active within the free-fall time of each cloud. This can vary from cloud to cloud due to the stochasticity of SN delay times and in HMXB formation. However, the defining factor in our simulations is the efficiency of the cooling, which can alter the Jeans mass required for sink particle formation, along with the number of massive stars in the cloud., 35 pages, 46 figures, accepted for publication in MNRAS

We investigate the formation of discs within binary systems where at least one component has left the main sequence. In particular we calculate the occurrence rates of systems which can host long-lived, massive discs that may be able to support the formation of planets. We synthesize a population based on Milky Way properties, using both theoretical and observational inputs to constrain key properties such as the shape of the initial mass function, binary fraction, and mass transfer physics. We predict 0.26\% of binary systems will host Second generation discs (where the primary has evolved), and 0.13\% of systems will host Third generation discs (where the secondary also evolves). For the Milky Way, this translates into 130 million and 90 million Second and Third generation systems respectively from an estimated total of 50 billion binary systems. Of these systems that form discs, we estimate approximately 20\% of Second and 3.8\% of Third generation discs have enough mass to form a planetary system. We speculate on how the process of planet formation differs in these systems from conventional planet formation in protostellar discs., 14 pages, 11 figures

The Aquistore CO2 Storage Site is located in southeastern Saskatchewan, Canada. CO2 is injected into a brine-filled sandstone formation at ∼3200 m depth immediately above the Precambrian basement. Sustained injection rates of 400-600 tonnes/day were achieved at the site starting in the fall of 2015 with a total of 88 ktonnes having been injected by the end of September, 2016. Seismic monitoring methods have been employed to track the subsurface CO2 plume and to record any injection-induced seismicity. Passive seismic monitoring is being conducted using two orthogonal arrays of short-period geophones, 3 broadband seismographs, and an array of downhole geophones. No significant injection-related seismicity (Mw > -1) has been detected during the first 17 months of CO2 injection. The first post-injection time-lapse 3D seismic surveys (surface and VSP) were conducted at the site in February, 2016. The VSP data were acquired with a distributed acoustic sensing system using a 2750 m casing-conveyed optical fibre cable in the observation well. 3D seismic modelling of fluid flow simulations in conjunction with seismic repeatability estimates obtained from field data indicate that the time-lapse VSP should be capable of imaging the CO2 plume after a total injection of ∼30 ktonnes. In addition, this first monitor survey tests the ability of surface seismic data acquired with a sparse permanent array to detect or image the CO2 plume after limited injection. Time-lapse logging is being conducted on a regular basis to provide in situ measurement of the change in seismic velocity associated with changes in CO2 saturation.

Fifteen years of XMM-Newton observations have established that ultrafast highly ionized winds (UFOs) are common in radio-quiet active galactic nucleus (AGN). A simple theory of Eddington-limited accretion correctly predicts the typical velocity (∼0.1c) and high ionization of such winds, with observed flow energy capable of ejecting star-forming gas. An extended XMM-Newton observation of the archetypal UFO PG 1211+143 recently found a more complex flow pattern, suggesting that intensive XMM-Newton observations offer exciting potential for probing the inner accretion disk structure and super-massive black hole (SMBH) growth.

Mergers are fundamental to the standard paradigm of galaxy evolution, and provide a natural formation mechanism for supermassive black hole binaries. The formation process of such a binary can have a direct impact on the rate at which stars are tidally disrupted by one or the other black hole, and the luminous signature of the tidal disruption itself can have distinct imprints of a binary companion. In this chapter we review our current understanding of the influence of black hole binarity on the properties of tidal disruption events. We discuss the rates of tidal disruption by supermassive black hole binaries, the impact of a second black hole on the fallback of debris and the formation of an accretion flow, and the prospects for detection of tidal disruption events by supermassive black hole binaries., 22 pages, 5 figures. Accepted chapter for Springer Space Science Reviews book on tidal disruption events

In a standard, steady, thin accretion disc, the radial distribution of the dissipation of the accretion energy is determined simply by energy considerations. Here we draw attention to the fact that while the (quasi-)steady discs in dwarf novae in outburst are in agreement with the expected emission distribution, the steady discs in the nova-like variables are not. We note that essentially the only difference between these two sets of discs is the time for which they have been in the high viscosity, high accretion rate state. In such discs, the major process by which angular momentum is transported outwards is MHD turbulence. We speculate that such turbulence gives rise to corona-like structures (here called magnetically controlled zones, or MCZs) which are also able to provide non-negligible angular momentum transport, the magnitude of which depends on the spatial scale $L$ of the magnetic field structures in such zones. For short-lived, high accretion rate discs (such as those in dwarf novae) we expect $L \sim H$ and the MCZ to have little effect. But, with time (such as in the nova-like variables) an inverse cascade in the MHD turbulence enables $L$, and the net effect of the MCZ, to grow. We present a simple toy model which demonstrates that such ideas can provide an explanation for the difference between the dwarf novae and the nova-like variable discs., 8 pages, 4 figures, accepted for publication in Astronomy and Astrophysics

The tidal disruption of stars by supermassive black holes (SMBHs) can be used to probe the SMBH mass function, the properties of individual stars, and stellar dynamics in galactic nuclei. Upcoming missions will detect thousands of TDEs, and accurate theoretical modeling is required to interpret the data with precision. Here we analyze the influence of more realistic stellar structure on the outcome of TDEs; in particular, we compare the fallback rates -- being the rate at which tidally-disrupted debris returns to the black hole -- from progenitors generated with the stellar evolution code {\sc mesa} to $\gamma = 4/3$ and $\gamma = 5/3$ polytropes. We find that {\sc mesa}-generated density profiles yield qualitatively-different fallback rates as compared to polytropic approximations, and that only the fallback curves from low-mass ($1M_{\odot}$ or less), zero-age main-sequence stars are well fit by either a $\gamma = 4/3$ or $5/3$ polytrope. Stellar age has a strong affect on the shape of the fallback curve, and can produce characteristic timescales (e.g., the time to the peak of the fallback rate) that greatly differ from the polytropic values. We use these differences to assess the degree to which the inferred black hole mass from the observed lightcurve can deviate from the true value, and find that the discrepancy can be at the order of magnitude level. Accurate stellar structure also leads to a substantial variation in the critical impact parameter at which the star is fully disrupted, and can increase the susceptibility of the debris stream to fragmentation under its own self-gravity. These results suggest that detailed modeling is required to accurately interpret observed lightcurves of TDEs., Comment: 14 pages, 9 figures, 2 tables. Updated version following referee report. Accepted for publication in ApJ Letters

The tidal disruption of a star by a supermassive black hole, and the subsequent accretion of the disrupted debris by that black hole, offers a direct means to study the inner regions of otherwise-quiescent galaxies. These tidal disruption events (TDEs) are being discovered at an ever-increasing rate. We present a model for the evolution of the tidally-disrupted debris from a partial TDE, in which a stellar core survives the initial tidal encounter and continues to exert a gravitational influence on the expanding stream of tidally-stripped debris. We use this model to show that the asymptotic fallback rate of material to the black hole in a partial TDE scales as $\propto t^{-2.26\pm0.01}$, and is effectively independent of the mass of the core that survives the encounter; we also estimate the rate at which TDEs approach this asymptotic scaling as a function of the core mass. These findings suggest that the late-time accretion rate onto a black hole from a TDE either declines as $t^{-5/3}$ if the star is completely disrupted or $t^{-9/4}$ if a core is left behind. We emphasize that previous investigations have not recovered this result due to the assumption of a Keplerian energy-period relationship for the debris orbits, which is no longer valid when a surviving core generates a time-dependent, gravitational potential. This dichotomy of fallback rates has important implications for the characteristic signatures of TDEs in the current era of wide-field surveys., ApJL Accepted

The tidal force from a supermassive black hole can rip apart a star that passes close enough in what is known as a Tidal Disruption Event. Typically half of the destroyed star remains bound to the black hole and falls back on highly eccentric orbits, forming an accretion flow which powers a luminous flare. In this paper we use analytical and numerical calculations to explore the effect of stellar rotation on the fallback rate of material. We find that slowly spinning stars ($\Omega_* \lesssim 0.01 \Omega_{\rm{breakup}}$) provide only a small perturbation to fallback rates found in the non-spinning case. However when the star spins faster, there can be significant effects. If the star is spinning retrograde with respect to its orbit the tidal force from the black hole has to spin down the star first before disrupting it, causing delayed and sometimes only partial disruption events. However, if the star is spinning prograde this works with the tidal force and the material falls back sooner and with a higher peak rate. We examine the power-law index of the fallback curves, finding that in all cases the fallback rate overshoots the canonical $t^{-5/3}$ rate briefly after the peak, with the depth of the overshoot dependent on the stellar spin. We also find that in general the late time evolution is slightly flatter than the canonical $t^{-5/3}$ rate. We therefore conclude that considering the spin of the star may be important in modelling observed TDE lightcurves., Comment: 14 pages, 11 figures, accepted for publication in ApJ

We investigate the impact of time-resolved `gradual' stellar feedback processes in high redshift dwarf spheroidal galaxies. Here `gradual' feedback refers to individual stellar feedback events which deposit energy over a period of time. We conduct high resolution hydrodynamical simulations of dwarf spheroidal galaxies with halo masses of 10$^7$ M$_{\odot}$ - 10$^8$ M$_{\odot}$, based on z = 6 progenitors of the Milky Way's dwarf spheroidal galaxies. We also include a novel feedback prescription for individual massive stars, which includes stellar winds and a HMXB (High Mass X-ray Binary) phase, on top of supernovae. We find the mass of gas unbound across a 1 Gyr starburst is uniformly lowered if gradual feedback mechanisms are included across the range of metallicities, halo concentration parameters and galaxy masses studied here. Furthermore, we find including gradual feedback in the smallest galaxies delays the unbinding of the majority of the gas and facilitates the production of `chimneys' in the dense shell surrounding the feedback generated hot, pressurised `superbubble'. These `chimneys' vent hot gas from the galaxy interior, lowering the temperature of the central 10 kpc of the gaseous halo. Additionally, we find radiative cooling has little effect on the energetics of simulations that include a short, violent starburst compared with those that have a longer, less concentrated starburst. Finally, we investigate the relative impact of HMXB feedback and stellar winds on our results, finding the ubiquity of stellar winds throughout each starburst makes them a defining factor in the final state of the interstellar medium., Comment: accepted for publication in MNRAS, 18 pages, 35 figures in main body of text (bibliography and appendices extra)

A star approaching a supermassive black hole (SMBH) can be torn apart in a tidal disruption event (TDE). We examine ultra-deep TDEs, a new regime in which the disrupted debris approaches close to the black hole's Schwarzschild radius, and the leading part intersects the trailing part at the first pericenter passage. We calculate the range of penetration factors $\beta$ vs SMBH masses $M$ that produce these prompt self-intersections using a Newtonian analytic estimate and a general relativistic (GR) geodesic model. We find that significant self-intersection of Solar-type stars requires $\beta \sim 50 - 127$ for $M/M_\odot = 10^4$, down to $\beta \sim 5.6 - 5.9$ for $M/M_\odot = 10^6$. We run smoothed-particle hydrodynamic (SPH) simulations to corroborate our calculations and find close agreement, with a slightly shallower dependence on $M$. We predict that the shock from the collision emits an X-ray flare lasting $t \sim 2$ s with $L \sim 10^{47}$ ergs/s at $E \sim 2$ keV, and the debris has a prompt accretion episode lasting $t \sim$ several min. The events are rare and occur with a rate $\dot{N} \lesssim 10^{-7}$ Mpc$^{-3}$ yr$^{-1}$. Ultra-deep TDEs can probe the strong gravity and demographics of low-mass SMBHs., Comment: 12 pages, 5 figures

Short gamma-ray bursts (sGRBs) are generally thought to result from the merger of two neutron stars or the merger of a neutron star with a black hole. It is becoming standard practise to model these mergers with hydrodynamical simulations that employ equations of state that are derived, for example, for determining the behaviour of matter in core-collapse supernovae (CCSNe), and which therefore make use of the assumption that the matter is hot and in nuclear statistical equilibrium (NSE). In this Letter we draw attention to the fact that the hydrodynamical timescale (roughly the gravitational timescale of the neutron star) may be several orders of magnitude shorter than the timescale on which such equilibrium can be re-established in the tidal debris ejected during a sGRB, and that on the initial decompression timescales the unshocked tidal ejecta may remain sufficiently cool that the employed equations of state are not appropriate for modelling the dynamics of this part of the flow. On timescales short compared with the timescale on which NSE can be (re)established, the equation of state can remain relatively stiff and thus the stream of tidal debris can remain narrow and vulnerable to gravitational instability, as has recently been suggested. These findings suggest that estimates of the type and abundances of heavy elements formed in short gamma-ray bursts need to be revisited. We suggest that the most direct method of testing the physical and dynamical properties of tidal ejecta in sGRBs will come from modelling of their light curves, which provides the cleanest source of information on the system dynamics., Comment: 8 pages. Accepted for publication in The Astrophysical Journal Letters

Filamentary structures, or long and narrow streams of material, arise in many areas of astronomy. Here we investigate the stability of such filaments by performing an eigenmode analysis of adiabatic and polytropic fluid cylinders, which are the cylindrical analog of spherical polytropes. We show that these cylinders are gravitationally unstable to perturbations along the axis of the cylinder below a critical wavenumber $k_{\rm crit} \simeq few$, where $k_{\rm crit}$ is measured relative to the radius of the cylinder. Below this critical wavenumber perturbations grow as $\propto e^{\sigma_{\rm u}\tau}$, where $\tau$ is time relative to the sound crossing time across the diameter of the cylinder, and we derive the growth rate $\sigma_{\rm u}$ as a function of wavenumber. We find that there is a maximum growth rate $\sigma_{\rm max} \sim 1$ that occurs at a specific wavenumber $k_{\rm max} \sim 1$, and we derive the growth rate $\sigma_{\rm max}$ and the wavenumbers $k_{\rm max}$ and $k_{\rm crit}$ for a range of adiabatic indices. To the extent that filamentary structures can be approximated as adiabatic and fluid-like, our results imply that these filaments are unstable without the need to appeal to magnetic fields or external media. Further, the objects that condense out of the instability of such filaments are separated by a preferred length scale, form over a preferred timescale, and possess a preferred mass scale., Comment: ApJS Accepted

Recent observations of several protoplanetary discs have found evidence of departures from flat, circular motion in the inner regions of the disc. One possible explanation for these observations is a disc warp, which could be induced by a planet on a misaligned orbit. We present three-dimensional numerical simulations of the tidal interaction between a protoplanetary disc and a misaligned planet. For low planet masses we show that our simulations accurately model the evolution of inclined planet orbit (up to moderate inclinations). For a planet massive enough to carve a gap, the disc is separated into two components and the gas interior and exterior to the planet orbit evolve separately, forming an inner and outer disc. Due to the inclination of the planet, a warp develops across the planet orbit such that there is a relative tilt and twist between these discs. We show that when other parameters are held constant, the relative inclination that develops between the inner and outer disc depends on the outer radius of the total disc modelled. For a given disc mass, our results suggest that the observational relevance of the warp depends more strongly on the mass of the planet rather than the inclination of the orbit., Accepted for publication in MNRAS. 17 pages, 12 Figures

Fluid injection is known to induce seismic events if the injection causes fracturing of the surrounding rock or if resulting pressure changes reactivate pre-existing faults and fractures. Carbon dioxide (CO2) storage projects where CO2 is injected into deep geological formations for permanent containment are one sector where induced seismicity has been observed. The Aquistore storage project in Saskatchewan, Canada began CO2 injection into the basal Cambrian sandstone at ∼3.2 km deep in April 2015 and the site has been extensively monitored for seismicity. Passive seismic monitoring instrumentation includes a small network of broadband seismometers, a continuously recording array of near-surface geophones and temporary deployments of downhole geophones at depths from 2950 m to 3010 m in an observation well. To date no injection-related induced seismicity has been observed. The seismic arrays are functioning as expected and local mine blasts, orientation shots and perforation shots have been detected using standard detection algorithms. Data stacking algorithms have also been tested on short-periods of data. Using synthetic data added to noise models, the estimated minimum detectable event local magnitude is −0.8 for the broadband stations and between −1.6 and −0.6 for the near-surface geophones. Thus far, small volumes of CO2 have been injected at Aquistore (∼140 kt) and injection has generally occurred below the fracture pressure. As a result, predicted pore pressure changes are small and periods without injection have allowed relaxation of the pressure plume. Geomechanical modelling suggests insignificant effective stress changes at an identified fault near the Aquistore injection well. It is therefore not surprising that no induced seismicity has been detected. With further injection, continued seismic monitoring is essential to provide warning of any fault reactivation and thus any potential increase in seismic risk or CO2 leakage risk.

Joint International Ocean Discovery Program and International Continental Scientific Drilling Program Expedition 364 drilled into the peak ring of the Chicxulub impact crater. We present P-wave velocity, density, and porosity measurements from Hole M0077A that reveal unusual physical properties of the peak-ring rocks. Across the boundary between post-impact sedimentary rock and suevite (impact melt-bearing breccia) we measure a sharp decrease in velocity and density, and an increase in porosity. Velocity, density, and porosity values for the suevite are 2900–3700 m/s, 2.06–2.37 g/cm3, and 20–35%, respectively. The thin (25 m) impact melt rock unit below the suevite has velocity measurements of 3650–4350 m/s, density measurements of 2.26–2.37 g/cm3, and porosity measurements of 19–22%. We associate the low velocity, low density, and high porosity of suevite and impact melt rock with rapid emplacement, hydrothermal alteration products, and observations of pore space, vugs, and vesicles. The uplifted granitic peak ring materials have values of 4000–4200 m/s, 2.39–2.44 g/cm3, and 8–13% for velocity, density, and porosity, respectively; these values differ significantly from typical unaltered granite which has higher velocity and density, and lower porosity. The majority of Hole M0077A peak-ring velocity, density, and porosity measurements indicate considerable rock damage, and are consistent with numerical model predictions for peak-ring formation where the lithologies present within the peak ring represent some of the most shocked and damaged rocks in an impact basin. We integrate our results with previous seismic datasets to map the suevite near the borehole. We map suevite below the Paleogene sedimentary rock in the annular trough, on the peak ring, and in the central basin, implying that, post impact, suevite covered the entire floor of the impact basin. Suevite thickness is 100–165 m on the top of the peak ring but 200 m in the central basin, suggesting that suevite flowed downslope from the collapsing central uplift during and after peak-ring formation, accumulating preferentially within the central basin.

Disc accretion on to a black hole is often misaligned from its spin axis. If the disc maintains a significant magnetic field normal to its local plane, we show that dipole radiation from Lense-Thirring precessing disc annuli can extract a significant fraction of the accretion energy, sharply peaked towards small disc radii $R$ (as $R^{-17/2}$ for fields with constant equipartition ratio). This low-frequency emission is immediately absorbed by surrounding matter or refracted towards the regions of lowest density. The resultant mechanical pressure, dipole angular pattern, and much lower matter density towards the rotational poles create a strong tendency to drive jets along the black hole spin axis, similar to the spin-axis jets of radio pulsars, also strong dipole emitters. The coherent primary emission may explain the high brightness temperatures seen in jets. The intrinsic disc emission is modulated at Lense-Thirring frequencies near the inner edge, providing a physical mechanism for low-frequency QPOs. Dipole emission requires nonzero hole spin, but uses only disc accretion energy. No spin energy is extracted, unlike the Blandford-Znajek process. MHD/GRMHD formulations do not directly give radiation fields, but can be checked post-process for dipole emission and so self-consistency, given sufficient resolution. Jets driven by dipole radiation should be more common in AGN than in X-ray binaries, and in low accretion rate states than high, agreeing with observation. In non-black-hole accretion, misaligned disc annuli precess because of the accretor's mass quadrupole moment, similarly producing jets and QPOs., 6 pages, no figures, accepted for publication in ApJL

We present Phantom, a fast, parallel, modular and low-memory smoothed particle hydrodynamics and magnetohydrodynamics code developed over the last decade for astrophysical applications in three dimensions. The code has been developed with a focus on stellar, galactic, planetary and high energy astrophysics and has already been used widely for studies of accretion discs and turbulence, from the birth of planets to how black holes accrete. Here we describe and test the core algorithms as well as modules for magnetohydrodynamics, self-gravity, sink particles, H_2 chemistry, dust-gas mixtures, physical viscosity, external forces including numerous galactic potentials as well as implementations of Lense-Thirring precession, Poynting-Robertson drag and stochastic turbulent driving. Phantom is hereby made publicly available., Comment: 88 pages, 60 figures, accepted to PASA. Code available from https://phantomsph.bitbucket.io/

Current planet formation theories provide successful frameworks with which to interpret the array of new observational data in this field. However, each of the two main theories (core accretion, gravitational instability) is unable to explain some key aspects. In many planet formation calculations, it is usual to treat the initial properties of the planet forming disc (mass, radius, etc.) as free parameters. In this paper, we stress the importance of setting the formation of planet forming discs within the context of the formation of the central stars. By exploring the early stages of disc formation, we introduce the concept of the Maximum Mass Solar Nebula (MMSN), as opposed to the oft-used Minimum Mass Solar Nebula (here mmsn). It is evident that almost all protoplanetary discs start their evolution in a strongly self-gravitating state. In agreement with almost all previous work in this area, we conclude that on the scales relevant to planet formation these discs are not gravitationally unstable to gas fragmentation, but instead form strong, transient spiral arms. These spiral arms can act as efficient dust traps allowing the accumulation and subsequent fragmentation of the dust (but not the gas). This phase is likely to populate the disc with relatively large planetesimals on short timescales while the disc is still veiled by a dusty-gas envelope. Crucially, the early formation of large planetesimals overcomes the main barriers remaining within the core accretion model. A prediction of this picture is that essentially all observable protoplanetary discs are already planet hosting., 6 pages, no figures, accepted for publication in MNRAS

Accretion discs are generally warped. If a warp in a disc is too large, the disc can `break' apart into two or more distinct planes, with only tenuous connections between them. Further if an initially planar disc is subject to a strong differential precession, then it can be torn apart into discrete annuli that precess effectively independently. In previous investigations, torque-balance formulae have been used to predict where and when the disc breaks into distinct parts. In this work, focusing on discs with Keplerian rotation and where the shearing motions driving the radial communication of the warp are damped locally by turbulence (the `diffusive' regime), we investigate the stability of warped discs to determine the precise criterion for an isolated warped disc to break. We find and solve the dispersion relation, which in general yields three roots. We provide a comprehensive analysis of this viscous-warp instability and the emergent growth rates and their dependence on disc parameters. The physics of the instability can be understood as a combination of (1) a term which would generally encapsulate the classical Lightman-Eardley instability in planar discs (given by $\partial(\nu\Sigma)/\partial\Sigma < 0$) but is here modified by the warp to include $\partial(\nu_1|\psi|)/\partial|\psi| < 0$ and (2) a similar condition acting on the diffusion of the warp amplitude given in simplified form by $\partial(\nu_2|\psi|)/\partial|\psi| < 0$. We discuss our findings in the context of discs with an imposed precession, and comment on the implications for different astrophysical systems., Comment: 13 pages, 7 figures. Accepted for publication in MNRAS

After the tidal disruption of a star by a massive black hole, disrupted stellar debris can fall back to the hole at a rate significantly exceeding its Eddington limit. To understand how black hole mass affects the duration of super-Eddington accretion in tidal disruption events, we first run a suite of simulations of the disruption of a Solar-like star by a supermassive black hole of varying mass to directly measure the fallback rate onto the hole, and we compare these fallback rates to the analytic predictions of the "frozen-in" model. Then, adopting a Zero-Bernoulli Accretion flow as an analytic prescription for the accretion flow around the hole, we investigate how the accretion rate onto the black hole evolves with the more accurate fallback rates calculated from the simulations. We find that numerically-simulated fallback rates yield accretion rates onto the hole that can, depending on the black hole mass, be nearly an order of magnitude larger than those predicted by the frozen-in approximation. Our results place new limits on the maximum black hole mass for which super-Eddington accretion occurs in tidal disruption events., Comment: 10 pages, 10 figures, MNRAS Accepted

Magnetic fields are observed in star forming regions. However simulations of the late stages of star formation that do not include magnetic fields provide a good fit to the properties of young stars including the initial mass function (IMF) and the multiplicity. We argue here that the simulations that do include magnetic fields are unable to capture the correct physics, in particular the high value of the magnetic Prandtl number, and the low value of the magnetic diffusivity. The artificially high (numerical and uncontrolled) magnetic diffusivity leads to a large magnetic flux pervading the star forming region. We argue further that in reality the dynamics of high magnetic Prandtl number turbulence may lead to local regions of magnetic energy dissipation through reconnection, meaning that the regions of molecular clouds which are forming stars might be essentially free of magnetic fields. Thus the simulations that ignore magnetic fields on the scales on which the properties of stellar masses, stellar multiplicities and planet-forming discs are determined, may be closer to reality than those which include magnetic fields, but can only do so in an unrealistic parameter regime., Comment: 28 pages, 0 figures, accepted for publication in New Astronomy

A tidal disruption event, which occurs when a star is destroyed by the gravitational field of a supermassive black hole, produces a stream of debris, the evolution of which ultimately determines the observational properties of the event. Here we show that a post-periapsis caustic -- a location where the locus of gas parcels comprising the stream would collapse into a two-dimensional surface if they evolved solely in the gravitational field of the hole -- occurs when the pericenter distance of the star is on the order of the tidal radius of the hole. It is demonstrated that this "pancake" induces significant density perturbations in the debris stream, and, for stiffer equations of state (adiabatic index $\gamma \gtrsim 5/3$), these fluctuations are sufficient to gravitationally destabilize the stream, resulting in its fragmentation into bound clumps. The results of our findings are discussed in the context of the observational properties of tidal disruption events., Comment: 18 pages, 14 figures. Accepted for publication in MNRAS

Accretion discs are common in binary systems, and they are often found to be misaligned with respect to the binary orbit. The gravitational torque from a companion induces nodal precession in misaligned disc orbits. We calculate whether this precession is strong enough to overcome the internal disc torques communicating angular momentum. For typical parameters precession wins: the disc breaks into distinct planes that precess effectively independently. We run hydrodynamical simulations to check these results, and confirm that disc breaking is widespread and generally enhances accretion on to the central object. This applies in many cases of astrophysical accretion, e.g. supermassive black hole binaries and X--ray binaries., 8 pages, 6 figures, accepted for publication in MNRAS

We analyse the interaction of an eccentric binary with a circular coplanar circumbinary disc that rotates in a retrograde sense with respect to the binary. In the circular binary case, no Lindblad resonances lie within the disc and no Lindblad resonant torques are produced, as was previously known. By analytic means, we show that when the binary orbit is eccentric, there exist components of the gravitational potential of the binary which rotate in a retrograde sense to the binary orbit and so rotate progradely with respect to this disc, allowing a resonant interaction to occur between the binary and the disc. The resulting resonant torques distinctly alter the disc response from the circular binary case. We describe results of three-dimensional hydrodynamic simulations to explore this effect and categorise the response of the disc in terms of modes whose strengths vary as a function of binary mass ratio and eccentricity. These mode strengths are weak compared to the largest mode strengths expected in the prograde case where the binary and disc rotate in the same sense. However, for sufficiently high binary eccentricity, resonant torques open a gap in a retrograde circumbinary disc, while permitting gas inflow on to the binary via gas streams. The inflow results in a time varying accretion rate on to the binary that is modulated over the binary orbital period, as was previously found to occur in the prograde case., 17 pages, 15 figures. Accepted to MNRAS

We investigate the effect of black hole spin on warped or misaligned accretion discs - in particular i) whether or not the inner disc edge aligns with the black hole spin and ii) whether the disc can maintain a smooth transition between an aligned inner disc and a misaligned outer disc, known as the Bardeen-Petterson effect. We employ high resolution 3D smoothed particle hydrodynamics simulations of $\alpha$-discs subject to Lense-Thirring precession, focussing on the bending wave regime where the disc viscosity is smaller than the aspect ratio $\alpha \lesssim H/R$. We first address the controversy in the literature regarding possible steady-state oscillations of the tilt close to the black hole. We successfully recover such oscillations in 3D at both small and moderate inclinations ($\lesssim 15^{\circ}$), provided both Lense-Thirring and Einstein precession are present, sufficient resolution is employed, and provided the disc is not so thick so as to simply accrete misaligned. Second, we find that discs inclined by more than a few degrees in general steepen and break rather than maintain a smooth transition, again in contrast to previous findings, but only once the disc scale height is adequately resolved. Finally, we find that when the disc plane is misaligned to the black hole spin by a large angle, the disc 'tears' into discrete rings which precess effectively independently and cause rapid accretion, consistent with previous findings in the diffusive regime ($\alpha \gtrsim H/R$). Thus misalignment between the disc and the spin axis of the black hole provides a robust mechanism for growing black holes quickly, regardless of whether the disc is thick or thin., Comment: 15 pages, 18 figures, movies available at http://users.monash.edu.au/~rnealon/ or YouTube

A hundred days out from the election a number of issues are buzzing: housing, immigration, water, climate change, electricity bills, and the perennials, economic growth, incomes and taxes and law and order. Based on previous contests, some of these will become the raw material of the political debates while others will fade to the background.

In some galaxies, the stars orbiting the supermassive black hole take the form of an eccentric nuclear disk, in which every star is on a coherent, apsidally-aligned orbit. The most famous example of an eccentric nuclear disk is the double nucleus of Andromeda, and there is strong evidence for many more in the local universe. Despite their apparent ubiquity however, a dynamical explanation for their longevity has remained a mystery: differential precession should wipe out large-scale apsidal-alignment on a short timescale. Here we identify a new dynamical mechanism which stabilizes eccentric nuclear disks, and explain for first time the negative eccentricity gradient seen in the Andromeda nucleus. The stabilizing mechanism drives oscillations of the eccentricity vectors of individual orbits, both in direction (about the mean body of the disk) and in magnitude. Combined with the negative eccentricity gradient, the eccentricity oscillations push some stars near the inner edge of the disk extremely close to the black hole, potentially leading to tidal disruption events. Order of magnitude calculations predict extremely high rates in recently-formed eccentric nuclear disks ($\sim0.1 - 1$ ${\rm yr}^{-1} {\rm gal}^{-1}$). Unless the stellar disks are replenished, these rates should decrease with time as the disk depletes in mass. If eccentric nuclear disks form during gas-rich major mergers, this may explain the preferential occurrence of tidal disruption events in recently-merged and post-merger (E+A/K+A) galaxies., Comment: Accepted to ApJ

Microlensing observations indicate that quasar accretion discs have half-light radii larger than expected from standard theoretical predictions based on quasar fluxes or black hole masses. Blackburne and colleagues have also found a very weak wavelength dependence of these half-light radii. We consider disc temperature profile models that might match these observations. Nixon and colleagues have suggested that misaligned accretion discs around spinning black holes will be disrupted at radii small enough for the Lense-Thirring torque to overcome the disc's viscous torque. Gas in precessing annuli torn off a disc will spread radially and intersect with the remaining disc, heating the disc at potentially large radii. However, if the intersection occurs at an angle of more than a degree or so, highly supersonic collisions will shock-heat the gas to a Compton temperature of T~10^7 K, and the spectral energy distributions (SEDs) of discs with such shock-heated regions are poor fits to observations of quasar SEDs. Torn discs where heating occurs in intermittent weak shocks that occur whenever the intersection angle reaches a tenth of a degree pose less of a conflict with observations, but do not have significantly larger half-light radii than standard discs. We also study two phenomenological disc temperature profile models. We find that discs with a temperature spike at relatively large radii and lowered temperatures at radii inside the spike yield improved and acceptable fits to microlensing sizes in most cases. Such temperature profiles could in principle occur in sub-Keplerian discs partially supported by magnetic pressure. However, such discs overpredict the fluxes from quasars studied with microlensing except in the limit of negligible continuum emission from radii inside the temperature spike.

The manner in which warps in accretion disks evolve depends on the magnitude of the viscosity. For small viscosity $(\alpha < H/R)$ the warp evolves in a wave-like manner; for large viscosity $H/R, Comment: Accepted for publication in ApJ

We present an accretion cycle which can explain state transitions and other observed phenomena in black hole X-ray binaries. This model is based on the process of disc tearing, where individual rings of gas break off the disc and precess effectively independently. This occurs when the Lense-Thirring effect is stronger than the local disc viscosity. We discuss implications of this model for quasi-periodic oscillations and the disc-jet-corona coupling. We also speculate on applying this model to active galactic nuclei and other accreting systems., Comment: 6 pages, 3 figures. Accepted to MNRAS

We investigate systematically whether accreting black hole systems are likely to reach global alignment of the black hole spin and its accretion disc with the binary plane. In low-mass X-ray binaries (LMXBs) there is only a modest tendency to reach such global alignment, and it is difficult to achieve fully: except for special initial conditions we expect misalignment of the spin and orbital planes by ~1 radian for most of the LMXB lifetime. The same is expected in high-mass X-ray binaries (HMXBs). A fairly close approach to global alignment is likely in most stellar-mass ultraluminous X-ray binary systems (ULXs) where the companion star fills its Roche lobe and transfers on a thermal timescale to a black hole of lower mass. These systems are unlikely to show orbital eclipses, as their emission cones are close to the hole's spin axis. This offers a potential observational test, as models for ULXs invoking intermediate-mass black holes do predict eclipses for ensembles of > ~10 systems. Recent observational work shows that eclipses are either absent or extremely rare in ULXs, supporting the picture that most ULXs are stellar-mass binaries with companion stars more massive than the accretor., Accepted for publication in MNRAS in original form. 4 pages, 1 figure

With hydrodynamical simulations we determine the conditions under which an initially coplanar planet-disc system that orbits a member of a misaligned binary star evolves to form a planet that undergoes Kozai-Lidov (KL) oscillations once the disc disperses. These oscillations may explain the large orbital eccentricities, as well as the large misalignments with respect to the spin of the central star, observed for some exoplanets. The planet is assumed to be massive enough to open a gap in the disc. The planet's tilt relative to the binary orbital plane is subject to two types of oscillations. The first type, present at even small inclination angles relative to the binary orbital plane, is due to the interaction of the planet with the disc and binary companion and is amplified by a secular resonance. The second type of oscillation is the KL oscillation that operates on both the planet and disc at larger binary inclination angles. We find that for a sufficiently massive disc, even a relatively low inclination planet-disc system can force a planet to an inclination above the critical KL angle, as a consequence of the first type of tilt oscillation, allowing it to undergo the second type of oscillation. We conclude that the hydrodynamical evolution of a sufficiently massive and inclined disc in a binary system broadens the range of systems that form eccentric and misaligned giant planets to include a wide range of initial misalignment angles (20 to 160 degrees)., Accepted for publication in MNRAS

We investigate the pre-disruption gravitational dynamics and post-disruption hydrodynamics of the tidal disruption of stars by supermassive black hole (SMBH) binaries. We focus on binaries with relatively low mass primaries ($10^6M_{\odot}$), moderate mass ratios, and separations with reasonably long gravitational wave inspiral times (tens of Myr). First, we generate a large ensemble (between 1 and 10 million) of restricted three-body integrations to quantify the statistical properties of tidal disruptions by circular SMBH binaries of initially-unbound stars. Compared to the reference case of a disruption by a single SMBH, the binary potential induces significant variance into the specific energy and angular momentum of the star at the point of disruption. Second, we use Newtonian numerical hydrodynamics to study the detailed evolution of the fallback debris from 120 disruptions randomly selected from the three-body ensemble (excluding only the most deeply penetrating encounters). We find that the overall morphology of the debris is greatly altered by the presence of the second black hole, and the accretion rate histories display a wide range of behaviors, including order of magnitude dips and excesses relative to control simulations that include only one black hole. Complex evolution persists, in some cases, for many orbital periods of the binary. We find evidence for power in the accretion curves on timescales related to the binary orbital period, though there is no exact periodicity. We discuss our results in the context of future wide-field surveys, and comment on the prospects of identifying and characterizing the subset of events occurring in nuclei with binary SMBHs., Comment: 31 pages, 25 figures, MNRAS Accepted; updated to incorporate referee comments. Movies of simulations can be found here: http://w.astro.berkeley.edu/~eric_coughlin/movies.html