Search

A star destroyed by a supermassive black hole (SMBH) in a tidal disruption event (TDE) is transformed into a filamentary structure known as a tidally disrupted stellar debris stream. We show that when ideal gas pressure dominates the thermodynamics of the stream, there is an exact solution to the hydrodynamics equations that describes the stream evolution and accounts for self-gravity, pressure, the dynamical expansion of the gas, and the transverse structure of the stream. We analyze the stability of this solution to cylindrically symmetric perturbations, and show that there is a critical stream density below which the stream is unstable and is not self-gravitating; this critical density is a factor of at least 40-50 smaller than the stream density in a TDE. Above this critical density the stream is overstable, self-gravity confines the stream, the oscillation period is exponentially long, and the growth rate of the overstability scales as $t^{1/6}$. The power-law growth and small power-law index of the overstability implies that the stream is effectively stable to cylindrically symmetric perturbations. We also use this solution to analyze the effects of hydrogen recombination, and suggest that even though recombination substantially increases the gas entropy, it is likely incapable of completely destroying the influence of self-gravity. We also show that the transient produced by recombination is far less luminous than previous estimates., Comment: 17 pages, 15 figures, MNRAS accepted

The dissociation and ionization of hydrogen, during the formation of giant planets via core accretion, reduce the effective adiabatic index γ of the gas and could trigger dynamical instability. We generalize the analysis of Chandrasekhar, who determined that the threshold for instability of a self-gravitating hydrostatic body lies at γ = 4/3, to account for the presence of a planetary core, which we model as an incompressible fluid. We show that the dominant effect of the core is to stabilize the envelope to radial perturbations, in some cases completely (i.e. for all γ > 1). When instability is possible, unstable planetary configurations occupy a strip of γ values whose upper boundary falls below γ = 4/3. Fiducial evolutionary tracks of giant planets forming through core accretion appear unlikely to cross the dynamical instability strip that we define.

Relativistic jets, or highly collimated and fast-moving outflows, are endemic to many astrophysical phenomena. The jets produced by gamma-ray bursts and tidal disruption events are accompanied by the accretion of material onto a black hole or neutron star, with the accretion rate exceeding the Eddington limit of the compact object by orders of magnitude. In such systems, radiation dominates the energy-momentum budget of the outflow, and the dynamical evolution of the jet is governed by the equations of radiation hydrodynamics. Here we show that there are analytic solutions to the equations of radiation hydrodynamics in the viscous (i.e., diffusive) regime that describe structured, relativistic jets, which consist of a fast-moving, highly relativistic core surrounded by a slower-moving, less relativistic sheath. In these solutions, the slower-moving, outer sheath contains most of the mass, and the jet structure is mediated by local anisotropies in the radiation field. We show that, depending on the pressure and density profile of the ambient medium, the angular profile of the jet Lorentz factor is Gaussian or falls off even more steeply with angle. These solutions have implications for the nature of jet production and evolution in hyperaccreting systems, and demonstrate that such jets -- and the corresponding jet structure -- can be sustained entirely by radiative processes. We discuss the implications of these findings in the context of jetted tidal disruption events and short and long gamma-ray bursts., Comment: MNRAS Accepted

The distribution of orbital energies imparted into stellar debris following the close encounter of a star with a supermassive black hole is the principal factor in determining the rate of return of debris to the black hole, and thus in determining the properties of the resulting lightcurves from such events. We present simulations of tidal disruption events for a range of $\beta\equiv r_{\rm t}/r_{\rm p}$ where $r_{\rm p}$ is the pericentre distance and $r_{\rm t}$ the tidal radius. We perform these simulations at different spatial resolutions to determine the numerical convergence of our models. We compare simulations in which the heating due to shocks is included or excluded from the dynamics. For $\beta \lesssim 8$ the simulation results are well-converged at sufficiently moderate-to-high spatial resolution, while for $\beta \gtrsim 8$ the breadth of the energy distribution can be grossly exaggerated by insufficient spatial resolution. We find that shock heating plays a non-negligible role only for $\beta \gtrsim 4$, and that typically the effect of shock heating is mild. We show that self-gravity can modify the energy distribution over time after the debris has receded to large distances for all $\beta$. Primarily, our results show that across a range of impact parameters, while the shape of the energy distribution varies with $\beta$, the width of the energy spread imparted to the bulk of the debris is closely matched to the canonical spread, $\Delta E = GM_\bullet R_\star/r_{\rm t}^2$, for the range of $\beta$ we have simulated., Comment: 17 pages, 7 figures, accepted for publication in ApJ

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

We present X-ray and radio observations of the Fast Blue Optical Transient (FBOT) CRTS-CSS161010 J045834-081803 (CSS161010 hereafter) at t=69-531 days. CSS161010 shows luminous X-ray ($L_x\sim5\times 10^{39}\,\rm{erg\,s^{-1}}$) and radio ($L_{\nu}\sim10^{29}\,\rm{erg\,s^{-1}Hz^{-1}}$) emission. The radio emission peaked at ~100 days post transient explosion and rapidly decayed. We interpret these observations in the context of synchrotron emission from an expanding blastwave. CSS161010 launched a mildly relativistic outflow with velocity $\Gamma\beta c\ge0.55c$ at ~100 days. This is faster than the non-relativistic AT2018cow ($\Gamma\beta c\sim0.1c$) and closer to ZTF18abvkwla ($\Gamma\beta c\ge0.3c$ at 63 days). The inferred initial kinetic energy of CSS161010 ($E_k\gtrsim10^{51}$ erg) is comparable to that of long Gamma Ray Bursts (GRBs), but the ejecta mass that is coupled to the mildly relativistic outflow is significantly larger ($\sim0.01-0.1\,\rm{M_{\odot}}$). This is consistent with the lack of observed gamma-rays. The luminous X-rays were produced by a different emission component to the synchrotron radio emission. CSS161010 is located at ~150 Mpc in a dwarf galaxy with stellar mass $M_{*}\sim10^{7}\,\rm{M_{\odot}}$ and specific star formation rate $sSFR\sim 0.3\,\rm{Gyr^{-1}}$. This mass is among the lowest inferred for host-galaxies of explosive transients from massive stars. Our observations of CSS161010 are consistent with an engine-driven aspherical explosion from a rare evolutionary path of a H-rich stellar progenitor, but we cannot rule out a stellar tidal disruption event on a centrally-located intermediate mass black hole. Regardless of the physical mechanism, CSS161010 establishes the existence of a new class of rare (rate $, Comment: Accepted to ApJL

An explosively generated shockwave with time-dependent radius $R(t)$ is characterized by a phase in which the shocked gas becomes radiative with an effective adiabatic index $\gamma \simeq 1$. Using the result that the post-shock gas is compressed into a shell of width $\Delta R/R \simeq \delta$, where $\delta = \gamma-1$, we show that a choice of self-similar variable that exploits this compressive behavior in the limit that $\gamma \rightarrow 1$ naturally leads to a series expansion of the post-shock fluid density, pressure, and velocity in the small quantity $\delta$. We demonstrate that the leading-order (in $\delta$) solutions, which are increasingly accurate as $\gamma \rightarrow 1$, can be written in simple, closed forms when the fluid is still approximated to be in the energy-conserving regime (i.e., the Sedov-Taylor limit), and that the density declines exponentially rapidly with distance behind the shock. We also analyze the solutions for the bubble surrounding a stellar or galactic wind that interacts with its surroundings, and derive expressions for the location of the contact discontinuity that separates the shocked ambient gas from the shocked wind. We discuss the implications of our findings in the context of the dynamical stability of nearly isothermal shocks., Comment: 5 pages, 2 figures, MNRASL Accepted

When a star passes close to a supermassive black hole (BH), the BH's tidal forces rip it apart into a thin stream, leading to a tidal disruption event (TDE). In this work, we study the post-disruption phase of TDEs in general relativistic hydrodynamics (GRHD) using our GPU-accelerated code H-AMR. We carry out the first grid-based simulation of a deep-penetration TDE ($\beta$=7) with realistic system parameters: a black-hole-to-star mass ratio of $10^6$, a parabolic stellar trajectory, and a nonzero BH spin. We also carry out a simulation of a tilted TDE whose stellar orbit is inclined relative to the BH midplane. We show that for our aligned TDE, an accretion disk forms due to the dissipation of orbital energy with $\sim$20 percent of the infalling material reaching the BH. The dissipation is initially dominated by violent self-intersections and later by stream-disk interactions near the pericenter. The self-intersections completely disrupt the incoming stream, resulting in five distinct self-intersection events separated by approximately 12 hours and a flaring in the accretion rate. We also find that the disk is eccentric with mean eccentricity e$\approx$0.88. For our tilted TDE, we find only partial self-intersections due to nodal precession near pericenter. Although these partial intersections eject gas out of the orbital plane, an accretion disk still forms with a similar accreted fraction of the material to the aligned case. These results have important implications for disk formation in realistic tidal disruptions. For instance, the periodicity in accretion rate induced by the complete stream disruption may explain the flaring events from Swift J1644+57., Comment: Accepted to MNRAS on November 23rd, 2021, 23 pages, 25 figures, uses mnras.cls. Comments welcome. Movies available at https://www.youtube.com/playlist?list=PL7YbfRC6zxzAqIPXWbJgXpr4Wq_nwHiDu

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

International audience; For the first time at Mars the statistical distribution of (1-D) electric field wave power in the magnetosphere is presented, along with the distribution of magnetic field wave power, as observed by the Mars Atmosphere and Volatile EvolutioN spacecraft from the first 14.5 months of the mission. Wave power in several different frequency bands was investigated, and the strongest wave powers were observed at the lowest frequencies. The presented statistical studies suggest that the full thermalization of ions within the magnetosheath does not appear to occur, as has been predicted by previous studies. Manual inspection of 140 periapsis passes on the dayside shows that Poynting fluxes (at 2-16 Hz) between ∼10-11 and 10-8 Wm-2 reach the upper ionosphere for all 140 cases. Wave power is not observed in the ionosphere for integrated electron densities greater than 1010.8 cm-2, corresponding to typical depths of 100-200 km. The observations presented support previous suggestions that energy from the Mars-solar wind interaction can propagate into the upper ionosphere and may provide an ionospheric heating source. Upstream of the shock, the orientation of the solar wind interplanetary magnetic field was shown to significantly affect the statistical distribution of wave power, based on whether the spacecraft was likely magnetically connected to the shock or not—something that is predicted but has not been quantitatively shown at Mars before. In flight performance and caveats of the Langmuir Probe and Waves electric field power spectra are also discussed.

Tidal disruption events, which occur when a star is shredded by the tidal field of a supermassive black hole, provide a means of fueling black hole accretion. Here we show, using a combination of three body orbit integrations and hydrodynamic simulations, that these events are also capable of generating circumbinary rings of gas around tight supermassive black hole binaries with small mass ratios. Depending on the thermodynamics, these rings can either fragment into clumps that orbit the binary, or evolve into a gaseous circumbinary disc. We argue that tidal disruptions provide a direct means of generating circumbinary discs around supermassive black hole binaries and, more generally, can replenish the reservoir of gas on very small scales in galactic nuclei., Comment: 5 pages, 2 figures, MNRAS Letters accepted

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

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

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

We analyze the outcome of the interaction between a stellar binary and a supermassive black hole binary (SMBHB) by performing a large number of gravitational scattering experiments. Most of the encounters result in either the ejection of an intact binary or the ejection of two individual stars following the tidal breakup of the binary. However, tidal disruption events (TDEs) and mergers constitute a few percent of the outcomes, and double, temporally distinct TDEs (i.e., separated by at least one orbit of the supermassive black hole binary) occur at the percent level. We also demonstrate that the properties of the ejected binaries are significantly altered through the interaction with the SMBHB, and their large eccentricities increase the merger rate and could lead to gravitational-wave inspirals far from the nucleus of the host galaxy. We discuss our results in the context of observed tidal disruption events, hypervelocity stars, and remote supernovae, such as calcium-rich transients., 9 pages, 4 figures, 1 table; ApJL Accepted

We present the first extensive radio to gamma-ray observations of a fast-rising blue optical transient (FBOT), AT2018cow, over its first ~100 days. AT2018cow rose over a few days to a peak luminosity $L_{pk}\sim4\times 10^{44}$ erg/s exceeding those of superluminous supernovae (SNe), before declining as $\propto t^{-2}$. Initial spectra at $\lesssim 15$ days were mostly featureless and indicated large expansion velocities v~0.1c and temperatures reaching 30000 K. Later spectra revealed a persistent optically-thick photosphere and the emergence of H and He emission features with v~sim 4000 km/s with no evidence for ejecta cooling. Our broad-band monitoring revealed a hard X-ray spectral component at $E\ge 10$ keV, in addition to luminous and highly variable soft X-rays, with properties unprecedented among astronomical transients. An abrupt change in the X-ray decay rate and variability appears to accompany the change in optical spectral properties. AT2018cow showed bright radio emission consistent with the interaction of a blastwave with $v_{sh}$~0.1c with a dense environment ($\dot M\sim10^{-3}-10^{-4}\,M_{\odot}yr^{-1}$ for $v_w=1000$ km\s). While these properties exclude Ni-powered transients, our multi-wavelength analysis instead indicates that AT2018cow harbored a "central engine", either a compact object (magnetar or black hole) or an embedded internal shock produced by interaction with a compact, dense circumstellar medium. The engine released $\sim10^{50}-10^{51.5}$ erg over $\sim10^3-10^5$ s and resides within low-mass fast-moving material with equatorial-polar density asymmetry ($M_{ej,fast}\lesssim0.3\,\rm{M_{\odot}}$). Successful SNe from low-mass H-rich stars (like electron-capture SNe) or failed explosions from blue supergiants satisfy these constraints. Intermediate-mass black-holes are disfavored by the large environmental density probed by the radio observations., Comment: 26 pages + tables. Submitted

Supermassive black holes (SMBHs) in galactic nuclei can eject hypervelocity stars (HVSs). Using restricted three-body integrations, we study the properties of stars ejected by circular, binary SMBHs as a function of their mass ratios $q = M_2 / M_1$ and separations $a$, focusing on the stellar velocity and angular distributions. We find that the ejection probability is an increasing function of $q$ and $a$, and that the mean ejected velocity scales with $q$ and $a$ similar to previous work but with modified scaling constants. Binary SMBHs tend to eject their fastest stars toward the binary orbital plane. We calculate the ejection rates as the binary SMBHs inspiral, and find that they eject stars with velocities $v_\infty > 1000$ km/s at rates of $\sim 4 \times 10^{-2} - 2 \times 10^{-1}$ yr$^{-1}$ for $q = 1$ ($\sim 10^{-4} - 10^{-3}$ yr$^{-1}$ for $q = 0.01$) over their lifetimes, and can emit a burst of HVSs with $v_\infty > 3000$ km/s as they coalesce. We integrate the stellar distributions over the binary SMBH inspiral and compare them to those produced by the "Hills mechanism" (in which a single SMBH ejects a star after tidally separating a binary star system), and find that $N \sim 100$ HVS velocity samples with $v_\infty \gtrsim 200$ km/s are needed to confidently distinguish between a binary and single SMBH origin., 19 pages, 16 figures

Stars approaching supermassive black holes (SMBHs) in the centers of galaxies can be torn apart by strong tidal forces. We study the physics of tidal disruption by a binary SMBH as a function of the binary mass ratio $q = M_2 / M_1$ and separation $a$, exploring a large set of points in the parameter range $q \in [0.01, 1]$ and $a/r_{t1} \in [10, 1000]$. We simulate encounters in which field stars approach the binary from the loss cone on parabolic, low angular momentum orbits. We present the rate of disruption and the orbital properties of the disrupted stars, and examine the fallback dynamics of the post-disruption debris in the "frozen-in" approximation. We conclude by calculating the time-dependent disruption rate over the lifetime of the binary. Throughout, we use a primary mass $M_1 = 10^6 M_\odot$ as our central example. We find that the tidal disruption rate is a factor of $\sim 2 - 7$ times larger than the rate for an isolated BH, and is independent of $q$ for $q \gtrsim 0.2$. In the "frozen-in" model, disruptions from close, nearly equal mass binaries can produce intense tidal fallbacks: for binaries with $q \gtrsim 0.2$ and $a/r_{t1} \sim 100$, roughly $\sim 18 - 40 \%$ of disruptions will have short rise times ($t_\textrm{rise} \sim 1 - 10$ d) and highly super-Eddington peak return rates ($\dot{M}_\textrm{peak} / \dot{M}_\textrm{Edd} \sim 2 \times 10^2 - 3 \times 10^3$)., Comment: 27 pages, 26 figures

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

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

The tidal forces close to massive black holes can rip apart stars that come too close to them. As the resulting stellar debris spirals toward the black hole, the debris heats up and emits x-rays. We report observations of a stable 131-second x-ray quasi-periodic oscillation from the tidal disruption event ASASSN-14li. Assuming the black hole mass indicated by host galaxy scaling relations, these observations imply that the periodicity originates from close to the event horizon and that the black hole is rapidly spinning. Our findings demonstrate that tidal disruption events can generate quasi-periodic oscillations that encode information about the physical properties of their black holes., United States. National Aeronautics and Space Administration (Grant PF6-170156), United States. National Aeronautics and Space Administration (Grant PF6-170150), United States. National Aeronautics and Space Administration (Grant PF5-160144), United States. National Aeronautics and Space Administration (Grant PF5-160145), United States. National Aeronautics and Space Administration (Award SV2-82023)

During the core collapse of massive stars, the formation of the protoneutron star is accompanied by the emission of a significant amount of mass-energy ($\sim 0.3 \, M_{\odot}$) in the form of neutrinos. This mass-energy loss generates an outward-propagating pressure wave that steepens into a shock near the stellar surface, potentially powering a weak transient associated with an otherwise-failed supernova. We analytically investigate this mass-loss-induced wave generation and propagation. Heuristic arguments provide an accurate estimate of the amount of energy contained in the outgoing sound pulse. We then develop a general formalism for analyzing the response of the star to centrally concentrated mass loss in linear perturbation theory. To build intuition, we apply this formalism to polytropic stellar models, finding qualitative and quantitative agreement with simulations and heuristic arguments. We also apply our results to realistic pre-collapse massive star progenitors (both giants and compact stars). Our analytic results for the sound pulse energy, excitation radius, and steepening in the stellar envelope are in good agreement with full time-dependent hydrodynamic simulations. We show that {prior} to the sound pulses arrival at the stellar photosphere, the photosphere has already reached velocities $\sim 20-100 \%$ of the local sound speed, thus likely modestly decreasing the stellar effective temperature prior to the star disappearing. Our results provide important constraints on the physical properties and observational appearance of failed supernovae., 15 pages, 10 figures, submitted to MNRAS

We study the ejection of mass during stellar core-collapse when the stalled shock does not revive and a black hole forms. Neutrino emission during the protoneutron star phase causes a decrease in the gravitational mass of the core, resulting in an outward going sound pulse that steepens into a shock as it travels out through the star. We explore the properties of this mass ejection mechanism over a range of stellar progenitors using spherically-symmetric, time-dependent hydrodynamic simulations that treat neutrino mass loss parametrically and follow the shock propagation over the entire star. We find that all types of stellar progenitor can eject mass through this mechanism. The ejected mass is a decreasing function of the surface gravity of the star, ranging from several $M_\odot$ for red supergiants to $\sim 0.1M_\odot$ for blue supergiants and $\sim 10^{-3} M_\odot $ for Wolf-Rayet stars. We find that the final shock energy at the surface is a decreasing function of the core-compactness, and is $\lesssim 10^{47}-10^{48}$ erg in all cases. In progenitors with a sufficiently large envelope, high core-compactness, or a combination of both, the sound pulse fails to unbind mass. Successful mass ejection is accompanied by significant fallback accretion that can last from hours to years. We predict the properties of shock breakout and thermal plateau emission produced by the ejection of the outer envelope of blue supergiant and Wolf-Rayet progenitors in otherwise failed supernovae., Comment: Accepted by MNRAS. Corrected errors in the evaluation of the analytical energy estimate (factor ~ 3) and in the derivation of the fallback accretion rate (factor ~ 2). Minor changes otherwise

Tidal disruption events (TDEs) observed in massive galaxies with inferred central black hole masses $M_h > 10^8 \ M_\odot$ are presumptive candidates for TDEs by lower mass secondaries in binary systems. We use hydrodynamic simulations to quantify the characteristics of such TDEs, focusing on extreme mass ratio binaries and mpc separations where the debris stream samples the binary potential. The simulations are initialised with disruption trajectories from 3-body integrations of stars with parabolic orbits with respect to the binary center of mass. The most common outcome is found to be the formation of an unbound debris stream, with either weak late-time accretion or no accretion at all. A substantial fraction of streams remain bound, however, and these commonly yield structured fallback rate curves that exhibit multiple peaks or sharp drops. We apply our results to the superluminous supernova candidate ASASSN-15lh and show that its features, including its anomalous rebrightening at $\sim 100$ days after detection, are consistent with the tidal disruption of a star by a supermassive black hole in a binary system., 10 pages, 3 figures; updated to include revisions from referee; fixed figures

Flares from tidal disruption events are unique tracers of quiescent black holes at the centre of galaxies. The appearance of these flares is very sensitive to whether the star is totally or partially disrupted, and in this paper we seek to identify the critical distance of the star from the black hole (r_d) that enables us to distinguish between these two outcomes. We perform here Mesh-free Finite Mass, traditional, and modern Smoothed Particle Hydrodynamical simulations of star-black hole close encounters, with the aim of checking if the value of r_d depends on the simulation technique. We find that the critical distance (or the so-called critical disruption parameter beta_d) depends only weakly on the adopted simulation method, being beta_d=0.92\pm 0.02 for a gamma=5/3 polytrope and beta_d=2.01\pm 0.01 for a gamma=4/3 polytrope., 8 pages, 8 figures, 2 tables; accepted to A&A

Coughlin et al. (2018) (Paper I) derived and analyzed a new regime of self-similarity that describes weak shocks (Mach number of order unity) in the gravitational field of a point mass. These solutions are relevant to low energy explosions, including failed supernovae. In this paper, we develop a formalism for analyzing the stability of shocks to radial perturbations, and we demonstrate that the self-similar solutions of Paper I are extremely weakly unstable to such radial perturbations. Specifically, we show that perturbations to the shock velocity and post-shock fluid quantities (the velocity, density, and pressure) grow with time as $t^{\alpha}$, where $\alpha \le 0.12$, implying that the ten-folding timescale of such perturbations is roughly ten orders of magnitude in time. We confirm these predictions by performing high-resolution, time-dependent numerical simulations. Using the same formalism, we also show that the Sedov-Taylor blastwave is trivially stable to radial perturbations provided that the self-similar, Sedov-Taylor solutions extend to the origin, and we derive simple expressions for the perturbations to the post-shock velocity, density, and pressure. Finally, we show that there is a third, self-similar solution (in addition to the the solutions in Paper I and the Sedov-Taylor solution) to the fluid equations that describes a rarefaction wave, i.e., an outward-propagating sound wave of infinitesimal amplitude. We interpret the stability of shock propagation in light of these three distinct self-similar solutions., Comment: ApJ Submitted

We present the exact solutions for the collapse of a spherically-symmetric, cold (i.e., pressureless) cloud under its own self-gravity, valid for arbitrary initial density profiles and not restricted to the realm of self-similarity. These solutions exhibit a number of remarkable features, including the self-consistent formation of and subsequent accretion onto a central point mass. A number of specific examples are provided, and we show that Penston's solution of pressureless, self-similar collapse is recovered for polytropic density profiles; importantly, however, we demonstrate that the time over which this solution holds is fleetingly narrow, implying that much of the collapse proceeds non-self-similarly. We show that our solutions can naturally incorporate turbulent pressure support, and we investigate the evolution of overdensities -- potentially generated by such turbulence -- as the collapse proceeds. Finally, we analyze the evolution of the angular velocity and magnetic fields in the limit that their dynamical influence is small, and we recover exact solutions for these quantities. Our results may provide important constraints on numerical models that attempt to elucidate the details of protostellar collapse when the initial conditions are far less idealized., Typos fixed in v2

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