-
Repeating nuclear transients from repeating partial tidal disruption events: reproducing ASASSN-14ko and AT2020vdq
Authors:
Ananya Bandopadhyay,
Eric R. Coughlin,
C. J. Nixon,
Dheeraj R. Pasham
Abstract:
Some electromagnetic outbursts from the nuclei of distant galaxies have been found to repeat on months-to-years timescales, and each of these sources can putatively arise from the accretion flares generated through the repeated tidal stripping of a star on a bound orbit about a supermassive black hole (SMBH), i.e., a repeating partial tidal disruption event (rpTDE). Here we test the rpTDE model th…
▽ More
Some electromagnetic outbursts from the nuclei of distant galaxies have been found to repeat on months-to-years timescales, and each of these sources can putatively arise from the accretion flares generated through the repeated tidal stripping of a star on a bound orbit about a supermassive black hole (SMBH), i.e., a repeating partial tidal disruption event (rpTDE). Here we test the rpTDE model through analytical estimates and hydrodynamical simulations of the interaction between a range of stars, which differ from one another in mass and age, and an SMBH. We show that higher-mass ($\gtrsim 1 M_{\odot}$), evolved stars can survive many ($\gtrsim 10-100$) encounters with an SMBH while simultaneously losing $few \times 0.01 M_{\odot}$, resulting in accretion flares that are approximately evenly spaced in time with nearly the same amplitude, quantitatively reproducing ASASSN-14ko. We also show that the energy imparted to the star via tides can lead to a change in its orbital period that is comparable to the observed decay in the recurrence time of ASASSN-14ko's flares, $\dot{P}\simeq-0.0026$. Contrarily, lower-mass and less-evolved stars lose progressively more mass and produce brighter accretion flares on subsequent encounters for the same pericenter distances, leading to the rapid destruction of the star and cessation of flares. Such systems cannot reproduce ASASSN-14ko-like transients, but are promising candidates for recreating events such as AT2020vdq, which displayed a second and much brighter outburst compared to the first. Our results imply that the lightcurves of repeating transients are tightly coupled with stellar type.
△ Less
Submitted 28 June, 2024; v1 submitted 5 June, 2024;
originally announced June 2024.
-
Quasi-stars as a Means of Rapid Black Hole Growth in the Early Universe
Authors:
Eric R. Coughlin,
Mitchell C. Begelman
Abstract:
JWST observations demonstrate that supermassive black holes (SMBHs) exist by redshifts $z \gtrsim 10$, providing further evidence for "direct collapse" black hole (BH) formation, whereby massive ($\sim 10^{3-5} M_{\odot}$) SMBH seeds are generated within a few Myr as a byproduct of the rapid inflow of gas into the centers of protogalaxies. Here we analyze the intermediate "quasi-star" phase that a…
▽ More
JWST observations demonstrate that supermassive black holes (SMBHs) exist by redshifts $z \gtrsim 10$, providing further evidence for "direct collapse" black hole (BH) formation, whereby massive ($\sim 10^{3-5} M_{\odot}$) SMBH seeds are generated within a few Myr as a byproduct of the rapid inflow of gas into the centers of protogalaxies. Here we analyze the intermediate "quasi-star" phase that accompanies some direct collapse models, during which a natal BH accretes mass from and energetically sustains (through accretion) an overlying gaseous envelope. We argue that previous estimates of the maximum BH mass that can be reached during this stage, $\sim 1\%$ of the total quasi-star mass, are unphysical, and arise from underestimating the efficiency with which energy can be transported outward from regions close to the BH. We construct new quasi-star models that consist of an inner, "saturated-convection" region (which conforms to a convection-dominated accretion flow near the BH) matched to an outer, adiabatic envelope. These solutions exist up to a BH mass of $\sim 60\%$ the total quasi-star mass, at which point the adiabatic envelope contains only 2\% of the mass (with the remaining $\sim 38\%$ in the saturated-convection region), and this upper limit is reached within a time of $20-40$ Myr. We conclude that quasi-stars remain a viable route for producing SMBHs at large redshifts, consistent with recent JWST observations.
△ Less
Submitted 10 June, 2024; v1 submitted 30 April, 2024;
originally announced May 2024.
-
Ultra-deep cover: an exotic and jetted tidal disruption event candidate disguised as a gamma-ray burst
Authors:
R. A. J. Eyles-Ferris,
C. J. Nixon,
E. R. Coughlin,
P. T. O'Brien
Abstract:
Gamma-ray bursts (GRBs) are traditionally classified as either short GRBs with durations $\lesssim 2$ s that are powered by compact object mergers, or long GRBs with durations $\gtrsim 2$ s that powered by the deaths of massive stars. Recent results, however, have challenged this dichotomy and suggest that there exists a population of merger-driven long bursts. One such example, GRB 191019A, has a…
▽ More
Gamma-ray bursts (GRBs) are traditionally classified as either short GRBs with durations $\lesssim 2$ s that are powered by compact object mergers, or long GRBs with durations $\gtrsim 2$ s that powered by the deaths of massive stars. Recent results, however, have challenged this dichotomy and suggest that there exists a population of merger-driven long bursts. One such example, GRB 191019A, has a $t_{90} \approx 64$ s but many of its other properties -- including its host galaxy, afterglow luminosity and lack of associated supernova -- are more consistent with a short GRB. Here we propose an alternative interpretation: that GRB 191019A (which is located in the nucleus of its host) is an atypical jetted tidal disruption event (TDE). In particular, we suggest the short timescale and rapid decline, not expected for standard TDEs, are the result of an "ultra-deep" encounter, in which the star came well within the tidal radius of the black hole and promptly self-intersected, circularised, accreted, and launched a relativistic outflow. This model reproduces the timescale and luminosity through a prompt super-Eddington accretion phase and accounts for the lack of late optical emission. This would make GRB 191019A only the fifth jetted TDE and the first discovered ultra-deep TDE. The ultra-deep TDE model can be distinguished from merger-driven long GRBs via the soft X-ray flash that results from prompt self-intersection of the debris stream; the detection of this flash will be possible with wide-field and soft-X-ray satellites such as $\textit{Einstein Probe}$ or $\textit{SVOM}$.
△ Less
Submitted 5 April, 2024;
originally announced April 2024.
-
Alive but Barely Kicking: News from 3+ years of Swift and XMM-Newton X-ray Monitoring of Quasi-Periodic Eruptions from eRO-QPE1
Authors:
Dheeraj R. Pasham,
Eric R. Coughlin,
Michal Zajacek,
Itai Linial,
Petra Sukova,
Christopher J. Nixon,
Agnieszka Janiuk,
Marzena Sniegowska,
Vojtech Witzany,
Vladimir Karas,
M. Krumpe,
Diego Altamirano,
Thomas Wevers,
Riccardo Arcodia
Abstract:
Quasi-periodic Eruptions (QPEs) represent a novel class of extragalactic X-ray transients that are known to repeat at roughly regular intervals of a few hours to days. Their underlying physical mechanism is a topic of heated debate, with most models proposing that they originate either from instabilities within the inner accretion flow or from orbiting objects. At present, our knowledge of how QPE…
▽ More
Quasi-periodic Eruptions (QPEs) represent a novel class of extragalactic X-ray transients that are known to repeat at roughly regular intervals of a few hours to days. Their underlying physical mechanism is a topic of heated debate, with most models proposing that they originate either from instabilities within the inner accretion flow or from orbiting objects. At present, our knowledge of how QPEs evolve over an extended timescale of multiple years is limited, except for the unique QPE source GSN 069. In this study, we present results from strategically designed Swift observing programs spanning the past three years, aimed at tracking eruptions from eRO-QPE1. Our main results are: 1) the recurrence time of eruptions can vary between 0.6 and 1.2 days, 2) there is no detectable secular trend in evolution of the recurrence times, 3) consistent with prior studies, their eruption profiles can have complex shapes, and 4) the peak flux of the eruptions has been declining over the past 3 years with the eruptions barely detected in the most recent Swift dataset taken in June of 2023. This trend of weakening eruptions has been reported recently in GSN 069. However, because the background luminosity of eRO-QPE1 is below our detection limit, we cannot verify if the weakening is correlated with the background luminosity (as is claimed to be the case for GSN 069). We discuss these findings within the context of various proposed QPE models.
△ Less
Submitted 14 February, 2024;
originally announced February 2024.
-
Lense-Thirring Precession after a Supermassive Black Hole Disrupts a Star
Authors:
Dheeraj R. Pasham,
Michal Zajacek,
C. J. Nixon,
Eric R. Coughlin,
Marzena Sniegowska,
Agnieszka Janiuk,
Bozena Czerny,
Thomas Wevers,
Muryel Guolo,
Yukta Ajay,
Michael Loewenstein
Abstract:
An accretion disk formed around a supermassive black hole (SMBH) after it disrupts a star is expected to be initially misaligned with respect to the black hole's equatorial plane. This misalignment induces relativistic torques (the Lense-Thirring effect) on the disk, causing the disk to precess at early times, while at late times the disk aligns with the black hole and precession terminates. Here,…
▽ More
An accretion disk formed around a supermassive black hole (SMBH) after it disrupts a star is expected to be initially misaligned with respect to the black hole's equatorial plane. This misalignment induces relativistic torques (the Lense-Thirring effect) on the disk, causing the disk to precess at early times, while at late times the disk aligns with the black hole and precession terminates. Here, using high-cadence X-ray monitoring observations of a TDE, we report the discovery of strong, quasi-periodic X-ray flux and temperature modulations from a TDE. These X-ray modulations are separated by 17.0$^{+1.2}_{-2.4}$ days and persist for roughly 130 days during the early phase of the TDE. Lense-Thirring precession of the accretion flow can produce this X-ray variability, but other physical mechanisms, such as the radiation-pressure instability, cannot be ruled out. Assuming typical TDE parameters, i.e., a solar-like star with the resulting disk extending at-most to so-called circularization radius, and that the disk precesses as a rigid body, we constrain the disrupting black hole's dimensionless spin parameter to be 0.05<|a|<0.5.
△ Less
Submitted 14 February, 2024;
originally announced February 2024.
-
Episodic X-ray Outflows from the Tidal Disruption Event ASASSN-14li
Authors:
Yukta Ajay,
Dheeraj R. Pasham,
Thomas Wevers,
Eric R. Coughlin,
Francesco Tombesi,
Muryel Guolo,
James F. Steiner
Abstract:
ASASSN-14li is a low-redshift ($z= 0.0206$) tidal disruption event (TDE) that has been studied extensively across the entire electromagnetic spectrum, and has provided one of the most sensitive measurements of a TDE to-date. Its X-ray spectrum is soft and thermal (kT$\sim$ 0.05 keV) and shows a residual broad absorption-like feature between 0.6-0.8 keV, which can be associated with a blue-shifted…
▽ More
ASASSN-14li is a low-redshift ($z= 0.0206$) tidal disruption event (TDE) that has been studied extensively across the entire electromagnetic spectrum, and has provided one of the most sensitive measurements of a TDE to-date. Its X-ray spectrum is soft and thermal (kT$\sim$ 0.05 keV) and shows a residual broad absorption-like feature between 0.6-0.8 keV, which can be associated with a blue-shifted O VII (rest-frame energy 0.57 keV) resulting from an ultrafast outflow (UFO) at early times (within 40 days of optical discovery). By carefully accounting for pile-up and using precise XSTAR photo-ionization table models, we analyze the entire archival X-ray data from XMM-Newton and track the evolution of this absorption feature for $\sim$4.5 years post disruption. Our main finding is that, contrary to the previous literature, the absorption feature is transient and intermittent. Assuming the same underlying physical basis (i.e. outflows) for the recurring absorption feature in ASASSN-14li, the outflow is seen to disappear and reappear multiple times during the first $\sim$2 years of its evolution. No observable spectral imprint is detected thereafter. While theoretical studies suggest the launch of outflows in the early phases of the outburst during the super-Eddington regime, the outflow's intermittent behavior for multiple years after disruption is unusual. We discuss this peculiar behavior within the context of varying inner disk truncation, radiation pressure, and magnetically-driven outflow scenarios.
△ Less
Submitted 23 January, 2024;
originally announced January 2024.
-
Delayed X-ray brightening accompanied by variable ionized absorption following a tidal disruption event
Authors:
T. Wevers,
M. Guolo,
D. R. Pasham,
E. R. Coughlin,
F. Tombesi,
Y. Yao,
S. Gezari
Abstract:
Supermassive black holes can experience super-Eddington peak mass fallback rates following the tidal disruption of a star. The theoretical expectation is that part of the infalling material is expelled by means of an accretion disk wind, whose observational signature includes blueshifted absorption lines of highly ionized species in X-ray spectra. To date, however, only one such ultra-fast outflow…
▽ More
Supermassive black holes can experience super-Eddington peak mass fallback rates following the tidal disruption of a star. The theoretical expectation is that part of the infalling material is expelled by means of an accretion disk wind, whose observational signature includes blueshifted absorption lines of highly ionized species in X-ray spectra. To date, however, only one such ultra-fast outflow (UFO) has been reported in the tidal disruption event (TDE) ASASSN-14li. Here we report on the discovery of transient absorption-like signatures in X-ray spectra of the TDE AT2020ksf/Gaia20cjk (at a redshift of $z$=0.092), following an X-ray brightening $\sim 230$ days after UV/optical peak. We find that while no statistically significant absorption features are present initially, they appear on a timescale of several days, and remain detected up to 770 days after peak. Simple thermal continuum models, combined with a power-law or neutral absorber, do not describe these features well. Adding a partial covering, low velocity ionized absorber improves the fit at early times, but fails at late times. A high velocity (v$_w$ $\sim$ 42000 km s$^{-1}$, or -0.15c), ionized absorber (ultra-fast outflow) provides a good fit to all data. The few day timescale of variability is consistent with expectations for a clumpy wind. We discuss several scenarios that could explain the X-ray delay, as well as the potential for larger scale wind feedback. The serendipitous nature of the discovery could suggest a high incidence of UFOs in TDEs, alleviating some of the tension with theoretical expectations.
△ Less
Submitted 8 January, 2024; v1 submitted 15 November, 2023;
originally announced November 2023.
-
The Peak of the Fallback Rate from Tidal Disruption Events: Dependence on Stellar Type
Authors:
Ananya Bandopadhyay,
Julia Fancher,
Aluel Athian,
Valentino Indelicato,
Sarah Kapalanga,
Angela Kumah,
Daniel A. Paradiso,
Matthew Todd,
Eric R. Coughlin,
C. J. Nixon
Abstract:
A star completely destroyed in a tidal disruption event (TDE) ignites a luminous flare that is powered by the fallback of tidally stripped debris to a supermassive black hole (SMBH) of mass $M_{\bullet}$. We analyze two estimates for the peak fallback rate in a TDE, one being the "frozen-in" model, which predicts a strong dependence of the time to peak fallback rate, $t_{\rm peak}$, on both stella…
▽ More
A star completely destroyed in a tidal disruption event (TDE) ignites a luminous flare that is powered by the fallback of tidally stripped debris to a supermassive black hole (SMBH) of mass $M_{\bullet}$. We analyze two estimates for the peak fallback rate in a TDE, one being the "frozen-in" model, which predicts a strong dependence of the time to peak fallback rate, $t_{\rm peak}$, on both stellar mass and age, with $15\textrm{ days} \lesssim t_{\rm peak} \lesssim 10$ yr for main sequence stars with masses $0.2\le M_{\star}/M_{\odot} \le 5$ and $M_{\bullet} = 10^6M_{\odot}$. The second estimate, which postulates that the star is completely destroyed when tides dominate the maximum stellar self-gravity, predicts that $t_{\rm peak}$ is very weakly dependent on stellar type, with $t_{\rm peak} = \left(23.2\pm4.0\textrm{ days}\right)\left(M_{\bullet}/10^6M_{\odot}\right)^{1/2}$ for $0.2\le M_{\star}/M_{\odot} \le 5$, while $t_{\rm peak} = \left(29.8\pm3.6\textrm{ days}\right)\left(M_{\bullet}/10^6M_{\odot}\right)^{1/2}$ for a Kroupa initial mass function truncated at $1.5 M_{\odot}$. This second estimate also agrees closely with hydrodynamical simulations, while the frozen-in model is discrepant by orders of magnitude. We conclude that (1) the time to peak luminosity in complete TDEs is almost exclusively determined by SMBH mass, and (2) massive-star TDEs power the largest accretion luminosities. Consequently, (a) decades-long extra-galactic outbursts cannot be powered by complete TDEs, including massive-star disruptions, and (b) the most highly super-Eddington TDEs are powered by the complete disruption of massive stars, which -- if responsible for producing jetted TDEs -- would explain the rarity of jetted TDEs and their preference for young and star-forming host galaxies.
△ Less
Submitted 17 October, 2023;
originally announced October 2023.
-
Low-energy Explosions in a Gravitational Field: Implications for Sub-energetic Supernovae and Fast X-ray Transients
Authors:
Daniel A. Paradiso,
Eric R. Coughlin,
Jonathan Zrake,
Dheeraj R. Pasham
Abstract:
Observations and theory suggest that core-collapse supernovae can span a range of explosion energies, and when sub-energetic, the shockwave initiating the explosion can decelerate to speeds comparable to the escape speed of the progenitor. In these cases, gravity will complicate the explosion hydrodynamics and conceivably cause the shock to stall at large radii within the progenitor star. To under…
▽ More
Observations and theory suggest that core-collapse supernovae can span a range of explosion energies, and when sub-energetic, the shockwave initiating the explosion can decelerate to speeds comparable to the escape speed of the progenitor. In these cases, gravity will complicate the explosion hydrodynamics and conceivably cause the shock to stall at large radii within the progenitor star. To understand these unique properties of weak explosions, we develop a perturbative approach for modeling the propagation of an initially strong shock into a time-steady, infalling medium in the gravitational field of a compact object. This method writes the shock position and the post-shock velocity, density, and pressure as series solutions in the (time-dependent) ratio of the freefall speed to the shock speed, and predicts that the shock stalls within the progenitor if the explosion energy is below a critical value. We show that our model agrees very well with hydrodynamic simulations, and accurately predicts (e.g.) the time-dependent shock position and velocity and the radius at which the shock stalls. Our results have implications for black hole formation and the newly detected class of fast X-ray transients (FXTs). In particular, we propose that a ``phantom shock breakout'' -- where the outer edge of the star falls through a stalled shock -- can yield a burst of X-rays without a subsequent optical/UV signature, similar to FXTs. This model predicts the rise time of the X-ray burst, $t_{\rm d}$, and the mean photon energy, $kT$, are anti-correlated, approximately as $T \propto t_{\rm d}^{-5/8}$.
△ Less
Submitted 4 December, 2023; v1 submitted 16 October, 2023;
originally announced October 2023.
-
On the relative importance of shocks and self-gravity in modifying tidal disruption event debris streams
Authors:
Julia Fancher,
Eric R. Coughlin,
Chris Nixon
Abstract:
In a tidal disruption event (TDE), a star is destroyed by the gravitational field of a supermassive black hole (SMBH) to produce a stream of debris, some of which accretes onto the SMBH and creates a luminous flare. The distribution of mass along the stream has a direct impact on the accretion rate, and thus modeling the time-dependent evolution of this distribution provides insight into the relev…
▽ More
In a tidal disruption event (TDE), a star is destroyed by the gravitational field of a supermassive black hole (SMBH) to produce a stream of debris, some of which accretes onto the SMBH and creates a luminous flare. The distribution of mass along the stream has a direct impact on the accretion rate, and thus modeling the time-dependent evolution of this distribution provides insight into the relevant physical processes that drive the observable properties of TDEs. Analytic models that only account for the ballistic evolution of the debris do not capture salient and time-dependent features of the mass distribution, suggesting that fluid dynamical effects significantly modify the debris dynamics. Previous investigations have claimed that shocks are primarily responsible for these modifications, but here we show -- with high-resolution hydrodynamical simulations -- that self-gravity is the dominant physical mechanism responsible for the anomalous (i.e., not predicted by ballistic models) debris stream features and its time dependence. These high-resolution simulations also show that there is a specific length scale on which self-gravity modifies the debris mass distribution, and as such there is enhanced power in specific Fourier modes. Our results have implications for the stability of the debris stream under the influence of self-gravity, particularly at late times and the corresponding observational signatures of TDEs.
△ Less
Submitted 13 September, 2023;
originally announced September 2023.
-
X-ray eruptions every 22 days from the nucleus of a nearby galaxy
Authors:
Muryel Guolo,
Dheeraj R. Pasham,
Michal Zajaček,
Eric R. Coughlin,
Suvi Gezari,
Petra Suková,
Thomas Wevers,
Vojtěch Witzany,
Francesco Tombesi,
Sjoert van Velzen,
Kate D. Alexander,
Yuhan Yao,
Riccardo Arcodia,
Vladimır Karas,
James Miller-Jones,
Ronald Remillard,
Keith Gendreau,
Elizabeth C. Ferrara
Abstract:
Galactic nuclei showing recurrent phases of activity and quiescence have recently been discovered, with recurrence times as short as a few hours to a day -- known as quasi-periodic X-ray eruption (QPE) sources -- to as long as hundreds to a thousand days for repeating nuclear transients (RNTs). Here we present a multi-wavelength overview of Swift J023017.0+283603 (hereafter Swift J0230+28), a sour…
▽ More
Galactic nuclei showing recurrent phases of activity and quiescence have recently been discovered, with recurrence times as short as a few hours to a day -- known as quasi-periodic X-ray eruption (QPE) sources -- to as long as hundreds to a thousand days for repeating nuclear transients (RNTs). Here we present a multi-wavelength overview of Swift J023017.0+283603 (hereafter Swift J0230+28), a source that exhibits repeating and quasi-periodic X-ray flares from the nucleus of a previously unremarkable galaxy at $\sim$ 165 Mpc, with a recurrence time of approximately 22 days, an intermediary timescale between known RNTs and QPE sources. The source also shows transient radio emission, likely associated with the X-ray emission. Such recurrent soft X-ray eruptions, with no accompanying UV/optical emission, are strikingly similar to QPE sources. However, in addition to having a recurrence time that is $\sim 25$ times longer than the longest-known QPE source, Swift J0230+28's eruptions exhibit somewhat distinct shapes and temperature evolution than the known QPE sources. Scenarios involving extreme mass ratio inspirals are favored over disk instability models. The source reveals an unexplored timescale for repeating extragalactic transients and highlights the need for a wide-field, time-domain X-ray mission to explore the parameter space of recurring X-ray transients.
△ Less
Submitted 15 January, 2024; v1 submitted 6 September, 2023;
originally announced September 2023.
-
The division between weak and strong explosions from failed supernovae
Authors:
Eric R. Coughlin
Abstract:
Some massive stars likely fail to produce core-collapse supernovae, but these failed supernovae (FSNe) can generate an electromagnetic outburst prior to the disappearance of the star, as the mass lost to neutrinos during the stellar core-collapse results in the formation and breakout of a second shock. We show that when the mass lost to neutrinos is sufficiently small, there are two self-similar s…
▽ More
Some massive stars likely fail to produce core-collapse supernovae, but these failed supernovae (FSNe) can generate an electromagnetic outburst prior to the disappearance of the star, as the mass lost to neutrinos during the stellar core-collapse results in the formation and breakout of a second shock. We show that when the mass lost to neutrinos is sufficiently small, there are two self-similar solutions that describe the propagation of a weak shock into a hydrodynamically expanding envelope that simultaneously yield accretion onto the black hole. The larger-Mach number solution is unstable and yields the minimum Mach number that a shock must have to strengthen into the energy-conserving regime. Above a critical mass loss there are no weak-shock solutions, implying that there are only strong explosions if the neutrino mass loss is above a critical value, and this value is a few percent of the mass of the star (and is physically achievable) for typical parameters. Our results imply that the fate of the explosion from a FSN -- weak with little to no mass ejection or strong with the expulsion of the majority of the envelope -- is a sensitive function of the stellar properties and the neutrino mass loss. We also show that there is a second type of self-similar solution for the shock that results in the ``settling'' of the gas near the compact object, which may be applicable to non-terminal stellar eruptions and the response of a gaseous disc to gravitational-wave induced mass loss from a binary black hole merger.
△ Less
Submitted 8 August, 2023;
originally announced August 2023.
-
The dynamics of debris streams from tidal disruption events: exact solutions, critical stream density, and hydrogen recombination
Authors:
Eric R. Coughlin
Abstract:
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, t…
▽ More
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.
△ Less
Submitted 2 May, 2023;
originally announced May 2023.
-
The luminosity function of TDEs from fallback-powered emission: implications for the black hole mass function
Authors:
Eric R. Coughlin,
Matt Nicholl
Abstract:
Tidal disruption events (TDEs), in which a star is destroyed by the gravitational field of a supermassive black hole (SMBH), are being observed at a high rate owing to the advanced state of survey science. One of the properties of TDEs that is measured with increasing statistical reliability is the TDE luminosity function, $d\dot{N}_{\rm TDE}/dL$, which is the TDE rate per luminosity (i.e., how ma…
▽ More
Tidal disruption events (TDEs), in which a star is destroyed by the gravitational field of a supermassive black hole (SMBH), are being observed at a high rate owing to the advanced state of survey science. One of the properties of TDEs that is measured with increasing statistical reliability is the TDE luminosity function, $d\dot{N}_{\rm TDE}/dL$, which is the TDE rate per luminosity (i.e., how many TDEs are within a given luminosity range). Here we show that if the luminous emission from a TDE is directly coupled to the rate of return of tidally destroyed debris to the SMBH, then the TDE luminosity function is in good agreement with observations and scales as $\propto L^{-2.5}$ for high luminosities, provided that the SMBH mass function $dN_{\bullet}/dM_{\bullet}$ -- the number of SMBHs ($N_{\bullet}$) per SMBH mass ($M_{\bullet}$) -- is approximately flat in the mass range over which we observe TDEs. We also show that there is a cutoff in the luminosity function at low luminosities that is a result of direct captures, and this cutoff has been tentatively observed. If $dN_{\bullet}/dM_{\bullet}$ is flat, which is in agreement with some observational campaigns, these results suggest that the fallback rate feeds the accretion rate in TDEs. Contrarily, if $dN_{\bullet}/d\log M_{\bullet}$ is flat, which has been found theoretically and is suggested by other observational investigations, then the emission from TDEs is likely powered by another mechanism. Future observations and more TDE statistics, provided by the Rubin Observatory/LSST, will provide additional evidence as to the reality of this tension.
△ Less
Submitted 27 April, 2023;
originally announced April 2023.
-
Measuring the properties of $f-$mode oscillations of a protoneutron star by third generation gravitational-wave detectors
Authors:
Chaitanya Afle,
Suman Kumar Kundu,
Jenna Cammerino,
Eric R Coughlin,
Duncan A. Brown,
David Vartanyan,
Adam Burrows
Abstract:
Core-collapse supernovae are among the astrophysical sources of gravitational waves that could be detected by third-generation gravitational-wave detectors. Here, we analyze the gravitational-wave strain signals from two- and three-dimensional simulations of core-collapse supernovae generated using the code F{\sc{ornax}}. A subset of the two-dimensional simulations has non-zero core rotation at th…
▽ More
Core-collapse supernovae are among the astrophysical sources of gravitational waves that could be detected by third-generation gravitational-wave detectors. Here, we analyze the gravitational-wave strain signals from two- and three-dimensional simulations of core-collapse supernovae generated using the code F{\sc{ornax}}. A subset of the two-dimensional simulations has non-zero core rotation at the core bounce. A dominant source of time changing quadrupole moment is the $l=2$ fundamental mode ($f-$ mode) oscillation of the proto-neutron star. From the time-frequency spectrogram of the gravitational-wave strain we see that, starting $\sim 400$ ms after the core bounce, most of the power lies within a narrow track that represents the frequency evolution of the $f-$mode oscillations. The $f-$mode frequencies obtained from linear perturbation analysis of the angle-averaged profile of the protoneutron star corroborate what we observe in the spectrograms of the gravitational-wave signal. We explore the measurability of the $f-$mode frequency evolution of protoneutron star for a supernova signal observed in the third-generation gravitational-wave detectors. Measurement of the frequency evolution can reveal information about the masses, radii, and densities of the proto-neutron stars. We find that if the third generation detectors observe a supernova within 10 kpc, we can measure these frequencies to within $\sim$90\% accuracy. We can also measure the energy emitted in the fundamental $f-$mode using the spectrogram data of the strain signal. We find that the energy in the $f-$mode can be measured to within 20\% error for signals observed by Cosmic Explorer using simulations with successful explosion, assuming source distances within 10 kpc.
△ Less
Submitted 9 April, 2023;
originally announced April 2023.
-
Tidal capture of stars by supermassive black holes: implications for periodic nuclear transients and quasi-periodic eruptions
Authors:
M. Cufari,
C. J. Nixon,
Eric R. Coughlin
Abstract:
Stars that plunge into the center of a galaxy are tidally perturbed by a supermassive black hole (SMBH), with closer encounters resulting in larger perturbations. Exciting these tides comes at the expense of the star's orbital energy, which leads to the naive conclusion that a smaller pericenter (i.e., a closer encounter between the star and SMBH) always yields a more tightly bound star to the SMB…
▽ More
Stars that plunge into the center of a galaxy are tidally perturbed by a supermassive black hole (SMBH), with closer encounters resulting in larger perturbations. Exciting these tides comes at the expense of the star's orbital energy, which leads to the naive conclusion that a smaller pericenter (i.e., a closer encounter between the star and SMBH) always yields a more tightly bound star to the SMBH. However, once the pericenter distance is small enough that the star is partially disrupted, morphological asymmetries in the mass lost by the star can yield an \emph{increase} in the orbital energy of the surviving core, resulting in its ejection -- not capture -- by the SMBH. Using smoothed-particle hydrodynamics simulations, we show that the combination of these two effects -- tidal excitation and asymmetric mass loss -- result in a maximum amount of energy lost through tides of $\sim 2.5\%$ of the binding energy of the star, which is significantly smaller than the theoretical maximum of the total stellar binding energy. This result implies that stars that are repeatedly partially disrupted by SMBHs many ($\gtrsim 10$) times on short-period orbits ($\lesssim$ few years), as has been invoked to explain the periodic nuclear transient ASASSN-14ko and quasi-periodic eruptions, must be bound to the SMBH through a mechanism other than tidal capture, such as a dynamical exchange (i.e., Hills capture).
△ Less
Submitted 3 January, 2023;
originally announced January 2023.
-
The Birth of a Relativistic Jet Following the Disruption of a Star by a Cosmological Black Hole
Authors:
Dheeraj R. Pasham,
Matteo Lucchini,
Tanmoy Laskar,
Benjamin P. Gompertz,
Shubham Srivastav,
Matt Nicholl,
Stephen J. Smartt,
James C. A. Miller-Jones,
Kate D. Alexander,
Rob Fender,
Graham P. Smith,
Michael D. Fulton,
Gulab Dewangan,
Keith Gendreau,
Eric R. Coughlin,
Lauren Rhodes,
Assaf Horesh,
Sjoert van Velzen,
Itai Sfaradi,
Muryel Guolo,
N. Castro Segura,
Aysha Aamer,
Joseph P. Anderson,
Iair Arcavi,
Sean J. Brennan
, et al. (41 additional authors not shown)
Abstract:
A black hole can launch a powerful relativistic jet after it tidally disrupts a star. If this jet fortuitously aligns with our line of sight, the overall brightness is Doppler boosted by several orders of magnitude. Consequently, such on-axis relativistic tidal disruption events (TDEs) have the potential to unveil cosmological (redshift $z>$1) quiescent black holes and are ideal test beds to under…
▽ More
A black hole can launch a powerful relativistic jet after it tidally disrupts a star. If this jet fortuitously aligns with our line of sight, the overall brightness is Doppler boosted by several orders of magnitude. Consequently, such on-axis relativistic tidal disruption events (TDEs) have the potential to unveil cosmological (redshift $z>$1) quiescent black holes and are ideal test beds to understand the radiative mechanisms operating in super-Eddington jets. Here, we present multi-wavelength (X-ray, UV, optical, and radio) observations of the optically discovered transient \target at $z=1.193$. Its unusual X-ray properties, including a peak observed luminosity of $\gtrsim$10$^{48}$ erg s$^{-1}$, systematic variability on timescales as short as 1000 seconds, and overall duration lasting more than 30 days in the rest-frame are traits associated with relativistic TDEs. The X-ray to radio spectral energy distributions spanning 5-50 days after discovery can be explained as synchrotron emission from a relativistic jet (radio), synchrotron self-Compton (X-rays), and thermal emission similar to that seen in low-redshift TDEs (UV/optical). Our modeling implies a beamed, highly relativistic jet akin to blazars but requires extreme matter-domination, i.e, high ratio of electron-to-magnetic field energy densities in the jet, and challenges our theoretical understanding of jets.
△ Less
Submitted 29 November, 2022;
originally announced November 2022.
-
Simulated optical light curves of super-Eddington tidal disruption events with ZEBRA flows
Authors:
R. A. J. Eyles-Ferris,
R. L. C. Starling,
P. T. O'Brien,
C. J. Nixon,
Eric R. Coughlin
Abstract:
We present simulated optical light curves of super-Eddington tidal disruption events (TDEs) using the zero-Bernoulli accretion (ZEBRA) flow model, which proposes that during the super-Eddington phase, the disc is quasi-spherical, radiation-pressure dominated, and accompanied by the production of strong jets. We construct light curves for both on- and off-axis (with respect to the jet) observers to…
▽ More
We present simulated optical light curves of super-Eddington tidal disruption events (TDEs) using the zero-Bernoulli accretion (ZEBRA) flow model, which proposes that during the super-Eddington phase, the disc is quasi-spherical, radiation-pressure dominated, and accompanied by the production of strong jets. We construct light curves for both on- and off-axis (with respect to the jet) observers to account for the anisotropic nature of the jetted emission. We find that at optical wavelengths, emission from the accretion flow is orders of magnitude brighter than that produced by the jet, even with boosting from synchrotron self-Compton. Comparing to the observed jetted TDE Swift J2058.4+0516, we find that the ZEBRA model accurately captures the timescale for which accretion remains super-Eddington and reproduces the luminosity of the transient. However, we find the shape of the light curves deviate at early times and the radius and temperature of our modelled ZEBRA are $\sim2.7 - 4.1$ times smaller and $\sim1.4 - 2.3$ times larger, respectively, than observed. We suggest that this indicates the ZEBRA inflates more, and more rapidly, than currently predicted by the model, and we discuss possible extensions to the model to account for this. Such refinements, coupled with valuable new data from upcoming large scale surveys, could help to resolve the nature of super-Eddington TDEs and how they are powered.
△ Less
Submitted 21 October, 2022;
originally announced October 2022.
-
Stars Crushed by Black Holes. III. Mild Compression of Radiative Stars by Supermassive Black Holes
Authors:
Suman Kumar Kundu,
Eric R. Coughlin,
C. J. Nixon
Abstract:
A tidal disruption event (TDE) occurs when the gravitational field of a supermassive black hole (SMBH) destroys a star. For TDEs in which the star enters deep within the tidal radius, such that the ratio of the tidal radius to the pericenter distance $β$ satisfies $β\gg 1$, the star is tidally compressed and heated. It was predicted that the maximum density and temperature attained during deep TDE…
▽ More
A tidal disruption event (TDE) occurs when the gravitational field of a supermassive black hole (SMBH) destroys a star. For TDEs in which the star enters deep within the tidal radius, such that the ratio of the tidal radius to the pericenter distance $β$ satisfies $β\gg 1$, the star is tidally compressed and heated. It was predicted that the maximum density and temperature attained during deep TDEs scale as $\propto β^3$ and $\propto β^2$, respectively, and nuclear detonation triggered by $β\gtrsim 5$, but these predictions have been debated over the last four decades. We perform Newtonian smoothed-particle hydrodynamics (SPH) simulations of deep TDEs between a Sun-like star and a $10^6 M_\odot$ SMBH for $2 \le β\le 10$. We find that neither the maximum density nor temperature follow the $\propto β^3$ and $\propto β^2$ scalings or, for that matter, any power-law dependence, and that the maximum-achieved density and temperature are reduced by $\sim$ an order of magnitude compared to past predictions. We also perform simulations in the Schwarzschild metric, and find that relativistic effects modestly increase the maximum density (by a factor of $\lesssim 1.5$) and induce a time lag relative to the Newtonian simulations, which is induced by time dilation. We also confirm that the time the star spends at high density and temperature is a very small fraction of its dynamical time. We therefore predict that the amount of nuclear burning achieved by radiative stars during deep TDEs is minimal.
△ Less
Submitted 4 October, 2022;
originally announced October 2022.
-
Live to die another day: the rebrightening of AT2018fyk as a repeating partial tidal disruption event
Authors:
T. Wevers,
E. R. Coughlin,
D. R. Pasham,
M. Guolo,
Y. Sun,
S. Wen,
P. G. Jonker,
A. Zabludoff,
A. Malyali,
R. Arcodia,
Z. Liu,
A. Merloni,
A. Rau,
I. Grotova,
P. Short,
Z. Cao
Abstract:
Stars that interact with supermassive black holes (SMBHs) can either be completely or partially destroyed by tides. In a partial tidal disruption event (TDE) the high-density core of the star remains intact, and the low-density, outer envelope of the star is stripped and feeds a luminous accretion episode. The TDE AT2018fyk, with an inferred black hole mass of $10^{7.7\pm0.4}$ M$_{\odot}$, experie…
▽ More
Stars that interact with supermassive black holes (SMBHs) can either be completely or partially destroyed by tides. In a partial tidal disruption event (TDE) the high-density core of the star remains intact, and the low-density, outer envelope of the star is stripped and feeds a luminous accretion episode. The TDE AT2018fyk, with an inferred black hole mass of $10^{7.7\pm0.4}$ M$_{\odot}$, experienced an extreme dimming event at X-ray (factor of $>$6000) and UV (factor $\sim$15) wavelengths $\sim$500--600 days after discovery. Here we report on the re-emergence of these emission components roughly 1200 days after discovery. We find that the source properties are similar to those of the pre-dimming accretion state, suggesting that the accretion flow was rejuvenated to a similar state. We propose that a repeating partial TDE, where the partially disrupted star is on a $\sim 1200$ day orbit about the SMBH and is periodically stripped of mass during each pericenter passage, powers its unique lightcurve. This scenario provides a plausible explanation for AT2018fyk's overall properties, including the rapid dimming event and the rebrightening at late times. We also provide testable predictions for the behavior of the accretion flow in the future: if the second encounter was also a partial disruption then we predict another strong dimming event around day 1800 (August 2023), and a subsequent rebrightening around day 2400 (March 2025). This source provides strong evidence of the partial disruption of a star by a SMBH.
△ Less
Submitted 7 November, 2022; v1 submitted 15 September, 2022;
originally announced September 2022.
-
A simple and accurate prescription for the tidal disruption radius of a star and the peak accretion rate in tidal disruption events
Authors:
Eric R. Coughlin,
Chris Nixon
Abstract:
A star destroyed by a supermassive black hole (SMBH) in a tidal disruption event (TDE) enables the study of SMBHs. We propose that the distance within which a star is completely destroyed by a SMBH, defined $r_{\rm t, c}$, is accurately estimated by equating the SMBH tidal field (including numerical factors) to the maximum gravitational field in the star. We demonstrate that this definition accura…
▽ More
A star destroyed by a supermassive black hole (SMBH) in a tidal disruption event (TDE) enables the study of SMBHs. We propose that the distance within which a star is completely destroyed by a SMBH, defined $r_{\rm t, c}$, is accurately estimated by equating the SMBH tidal field (including numerical factors) to the maximum gravitational field in the star. We demonstrate that this definition accurately reproduces the critical $β_{\rm c} = r_{\rm t}/r_{\rm t, c}$, where $r_{\rm t} = R_{\star}\left(M_{\bullet}/M_{\star}\right)^{1/3}$ is the standard tidal radius with $R_{\star}$ and $M_{\star}$ the stellar radius and mass and $M_{\bullet}$ the SMBH mass, for multiple stellar progenitors at various ages, and can be reasonably approximated by $β_{\rm c} \simeq \left[ρ_{\rm c}/(4ρ_{\star})\right]^{1/3}$, where $ρ_{\rm c}$ ($ρ_{\star}$) is the central (average) stellar density. We also calculate the peak fallback rate and time at which the fallback rate peaks, finding excellent agreement with hydrodynamical simulations, and also suggest that the partial disruption radius -- the distance at which any mass is successfully liberated from the star -- is $β_{\rm partial} \simeq 4^{-1/3} \simeq 0.6$. For given stellar and SMBH populations, this model yields, e.g., the fraction of partial TDEs, the peak luminosity distribution of TDEs, and the number of directly captured stars.
△ Less
Submitted 8 September, 2022;
originally announced September 2022.
-
Spherically Symmetric Accretion onto a Compact Object through a Standing Shock: The Effects of General Relativity in the Schwarzschild Geometry
Authors:
Suman Kumar Kundu,
Eric R. Coughlin
Abstract:
A core-collapse supernova is generated by the passage of a shockwave through the envelope of a massive star, where the shock wave is initially launched from the ``bounce'' of the neutron star formed during the collapse of the stellar core. Instead of successfully exploding the star, however, numerical investigations of core-collapse supernovae find that this shock tends to ``stall'' at small radii…
▽ More
A core-collapse supernova is generated by the passage of a shockwave through the envelope of a massive star, where the shock wave is initially launched from the ``bounce'' of the neutron star formed during the collapse of the stellar core. Instead of successfully exploding the star, however, numerical investigations of core-collapse supernovae find that this shock tends to ``stall'' at small radii ($\lesssim$ 10 neutron star radii), with stellar material accreting onto the central object through the standing shock. Here, we present time-steady, adiabatic solutions for the density, pressure, and velocity of the shocked fluid that accretes onto the compact object through the stalled shock, and we include the effects of general relativity in the Schwarzschild metric. Similar to previous works that were carried out in the Newtonian limit, we find that the gas ``settles'' interior to the stalled shock; in the relativistic regime analyzed here, the velocity asymptotically approaches zero near the Schwarzschild radius. These solutions can represent accretion onto a material surface if the radius of the compact object is outside of its event horizon, such as a neutron star; we also discuss the possibility that these solutions can approximately represent the accretion of gas onto a newly formed black hole following a core-collapse event. Our findings and solutions are particularly relevant in weak and failed supernovae, where the shock is pushed to small radii and relativistic effects are large.
△ Less
Submitted 31 August, 2022;
originally announced September 2022.
-
On the impact of relativistic gravity on the rate of tidal disruption events
Authors:
Eric R. Coughlin,
Chris Nixon
Abstract:
The tidal disruption of stars by supermassive black holes (SMBHs) probes relativistic gravity. In the coming decade, the number of observed tidal disruption events (TDEs) will grow by several orders of magnitude, allowing statistical inferences of the properties of the SMBH and stellar populations. Here we analyse the probability distribution functions of the pericentre distances of stars that enc…
▽ More
The tidal disruption of stars by supermassive black holes (SMBHs) probes relativistic gravity. In the coming decade, the number of observed tidal disruption events (TDEs) will grow by several orders of magnitude, allowing statistical inferences of the properties of the SMBH and stellar populations. Here we analyse the probability distribution functions of the pericentre distances of stars that encounter an SMBH in the Schwarzschild geometry, where the results are completely analytic, and the Kerr metric. From this analysis we calculate the number of observable TDEs, defined to be those that come within the tidal radius $r_{\rm t}$ but outside the direct capture radius (which is, in general, larger than the horizon radius). We find that relativistic effects result in a steep decline in the number of stars that have pericenter distances $r_{\rm p} \lesssim 10\,r_{\rm g}$, where $r_{\rm g} = GM/c^2$, and that for maximally spinning SMBHs the distribution function of $r_{\rm p}$ at such distances scales as $f_{\rm r_{\rm p}}\propto r_{\rm p}^{4/3}$, or in terms of $β\equiv r_{\rm t}/r_{\rm p}$ scales as $f_β \propto β^{-10/3}$. We find that spin has little effect on the TDE fraction until the very high-mass end, where instead of being identically zero the rate is small ($\lesssim 1\%$ of the expected rate in the absence of relativistic effects). Effectively independent of spin, if the progenitors of TDEs reflect the predominantly low-mass stellar population and thus have masses $\lesssim 1M_{\odot}$, we expect a substantial reduction in the rate of TDEs above $10^{7}M_{\odot}$.
△ Less
Submitted 28 July, 2022;
originally announced July 2022.
-
Using the Hills Mechanism to Generate Repeating Partial Tidal Disruption Events and ASASSN-14ko
Authors:
M. Cufari,
Eric R. Coughlin,
C. J. Nixon
Abstract:
Periodic nuclear transients have been detected with increasing frequency, with one such system -- ASASSN-14ko -- exhibiting highly regular outbursts on a timescale of $114 \pm 1$ days. It has been postulated that the outbursts from this source are generated by the repeated partial disruption of a star, but how the star was placed onto such a tightly bound orbit about the supermassive black hole re…
▽ More
Periodic nuclear transients have been detected with increasing frequency, with one such system -- ASASSN-14ko -- exhibiting highly regular outbursts on a timescale of $114 \pm 1$ days. It has been postulated that the outbursts from this source are generated by the repeated partial disruption of a star, but how the star was placed onto such a tightly bound orbit about the supermassive black hole remains unclear. Here we use analytic arguments and three-body integrations to demonstrate that the Hills mechanism, where a binary system is destroyed by the tides of the black hole, can lead to the capture of a star on a $\sim 114$ day orbit and with a pericenter distance that is comparable to the tidal radius of one of the stars within the binary. Thus, Hills capture can produce stars on tightly bound orbits that undergo repeated partial disruption, leading to a viable mechanism for generating not only the outbursts detected from ASASSN-14ko, but for periodic nuclear transients in general. We also show that the rate of change of the period of the captured star due to gravitational-wave emission is likely too small to produce the observed value for ASASSN-14ko, indicating that in this system there must be additional effects that contribute to the decay of the orbit. In general, however, gravitational-wave emission can be important for limiting the lifetimes of these systems, and could produce observable period decay rates in future events.
△ Less
Submitted 15 March, 2022;
originally announced March 2022.
-
Stellar Revival and Repeated Flares in Deeply Plunging Tidal Disruption Events
Authors:
Chris Nixon,
Eric R. Coughlin
Abstract:
Tidal disruption events with tidal radius $r_{\rm t}$ and pericenter distance $r_{\rm p}$ are characterized by the quantity $β= r_{\rm t}/r_{\rm p}$, and "deep encounters" have $β\gg 1$. It has been assumed that there is a critical $β\equiv β_{\rm c} \sim 1$ that differentiates between partial and full disruption: for $β< β_{\rm c}$ a fraction of the star survives the tidal interaction with the bl…
▽ More
Tidal disruption events with tidal radius $r_{\rm t}$ and pericenter distance $r_{\rm p}$ are characterized by the quantity $β= r_{\rm t}/r_{\rm p}$, and "deep encounters" have $β\gg 1$. It has been assumed that there is a critical $β\equiv β_{\rm c} \sim 1$ that differentiates between partial and full disruption: for $β< β_{\rm c}$ a fraction of the star survives the tidal interaction with the black hole, while for $β> β_{\rm c}$ the star is completely destroyed, and hence all deep encounters should be full. Here we show that this assumption is incorrect by providing an example of a $β= 16$ encounter between a $γ= 5/3$, solar-like polytrope and a $10^6 M_{\odot}$ black hole -- for which previous investigations have found $β_{\rm c} \simeq 0.9$ -- that results in the reformation of a stellar core post-disruption that comprises approximately 25\% of the original stellar mass. We propose that the core reforms under self-gravity, which remains important because of the compression of the gas both near pericenter, where the compression occurs out of the orbital plane, and substantially after pericenter, where compression is within the plane. We find that the core forms on a bound orbit about the black hole, and we discuss the corresponding implications of our findings in the context of recently observed, repeating nuclear transients.
△ Less
Submitted 31 January, 2022;
originally announced February 2022.
-
A Physical Model of Delayed Rebrightenings in Shock-interacting Supernovae without Narrow Line Emission
Authors:
Eric R. Coughlin,
Jonathan Zrake
Abstract:
Core-collapse supernovae can display evidence of interaction with pre-existing, circumstellar shells of material by rebrightening and forming spectral lines, and can even change types as Hydrogen appears in previously Hydrogen-poor spectra. However, a recently observed core-collapse supernova -- SN 2019tsf -- was found to brighten after roughly 100 days after it was first observed, suggesting that…
▽ More
Core-collapse supernovae can display evidence of interaction with pre-existing, circumstellar shells of material by rebrightening and forming spectral lines, and can even change types as Hydrogen appears in previously Hydrogen-poor spectra. However, a recently observed core-collapse supernova -- SN 2019tsf -- was found to brighten after roughly 100 days after it was first observed, suggesting that the supernova ejecta was interacting with surrounding material, but it lacked any observable emission lines and thereby challenged the standard supernova-interaction picture. We show through linear perturbation theory that delayed rebrightenings without the formation of spectral lines are generated as a consequence of the finite sound crossing time of the post-shock gas left in the wake of a supernova explosion. In particular, we demonstrate that sound waves -- generated in the post-shock flow as a consequence of the interaction between a shock and a density enhancement -- traverse the shocked ejecta and impinge upon the shock from behind in a finite time, generating sudden changes in the shock properties in the absence of ambient density enhancements. We also show that a blastwave dominated by gas pressure and propagating in a wind-fed medium is unstable from the standpoint that small perturbations lead to the formation of reverse shocks within the post-shock flow, implying that the gas within the inner regions of these blastwaves should be highly turbulent.
△ Less
Submitted 29 November, 2021;
originally announced November 2021.
-
Stars Crushed by Black Holes. II. A Physical Model of Adiabatic Compression and Shock Formation in Tidal Disruption Events
Authors:
Eric R. Coughlin,
Chris Nixon
Abstract:
We develop a Newtonian model of a deep tidal disruption event (TDE), for which the pericenter distance of the star, $r_{\rm p}$, is well within the tidal radius of the black hole, $r_{\rm t}$, i.e., when $β\equiv r_{\rm t}/r_{\rm p} \gg 1$. We find that shocks form for $β\gtrsim 3$, but they are weak (with Mach numbers $\sim 1$) for all $β$, and that they reach the center of the star prior to the…
▽ More
We develop a Newtonian model of a deep tidal disruption event (TDE), for which the pericenter distance of the star, $r_{\rm p}$, is well within the tidal radius of the black hole, $r_{\rm t}$, i.e., when $β\equiv r_{\rm t}/r_{\rm p} \gg 1$. We find that shocks form for $β\gtrsim 3$, but they are weak (with Mach numbers $\sim 1$) for all $β$, and that they reach the center of the star prior to the time of maximum adiabatic compression for $β\gtrsim 10$. The maximum density and temperature reached during the TDE follow much shallower relations with $β$ than the previously predicted $ρ_{\rm max} \propto β^3$ and $T_{\rm max} \propto β^2$ scalings. Below $β\simeq 10$, this shallower dependence occurs because the pressure gradient is dynamically significant before the pressure is comparable to the ram pressure of the freefalling gas, while above $β\simeq 10$ we find that shocks prematurely halt the compression and yield the scalings $ρ_{\rm max} \propto β^{1.62}$ and $T_{\rm max} \propto β^{1.12}$. We find excellent agreement between our results and high-resolution simulations. Our results demonstrate that, in the Newtonian limit, the compression experienced by the star is completely independent of the mass of the black hole. We discuss our results in the context of existing (affine) models, polytropic vs.~non-polytropic stars, and general relativistic effects, which become important when the pericenter of the star nears the direct capture radius, at $β\sim 12.5$ (2.7) for a solar-like star disrupted by a $10^6M_{\odot}$ ($10^{7}M_{\odot}$) supermassive black hole.
△ Less
Submitted 24 November, 2021;
originally announced November 2021.
-
Stars Crushed by Black Holes. I. On the Energy Distribution of Stellar Debris in Tidal Disruption Events
Authors:
Sarah Norman,
Chris Nixon,
Eric R. Coughlin
Abstract:
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…
▽ More
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 $β\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 $β\lesssim 8$ the simulation results are well-converged at sufficiently moderate-to-high spatial resolution, while for $β\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 $β\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 $β$. Primarily, our results show that across a range of impact parameters, while the shape of the energy distribution varies with $β$, the width of the energy spread imparted to the bulk of the debris is closely matched to the canonical spread, $ΔE = GM_\bullet R_\star/r_{\rm t}^2$, for the range of $β$ we have simulated.
△ Less
Submitted 24 November, 2021;
originally announced November 2021.
-
The Eccentric Nature of Eccentric Tidal Disruption Events
Authors:
M. Cufari,
Eric. R. Coughlin,
C. J. Nixon
Abstract:
Upon entering the tidal sphere of a supermassive black hole, a star is ripped apart by tides and transformed into a stream of debris. The ultimate fate of that debris, and the properties of the bright flare that is produced and observed, depends on a number of parameters, including the energy of the center of mass of the original star. Here we present the results of a set of smoothed particle hydr…
▽ More
Upon entering the tidal sphere of a supermassive black hole, a star is ripped apart by tides and transformed into a stream of debris. The ultimate fate of that debris, and the properties of the bright flare that is produced and observed, depends on a number of parameters, including the energy of the center of mass of the original star. Here we present the results of a set of smoothed particle hydrodynamics simulations in which a $1~M_\odot $, $γ= 5/3$ polytrope is disrupted by a $10^6 ~M_\odot$ supermassive black hole. Each simulation has a pericenter distance of $r_{\rm p} = r_{\rm t}$ (i.e., $β\equiv r_{\rm t}/r_{\rm p} = 1$ with $r_{\rm t}$ the tidal radius), and we vary the eccentricity $e$ of the stellar orbit from $e = 0.8$ up to $e = 1.20$ and study the nature of the fallback of debris onto the black hole and the long-term fate of the unbound material. For simulations with eccentricities $e \lesssim 0.98$, the fallback curve has a distinct, three-peak structure that is induced by self-gravity. For simulations with eccentricities $e \gtrsim 1.06$, the core of the disrupted star reforms following its initial disruption. Our results have implications for, e.g., tidal disruption events produced by supermassive black hole binaries.
△ Less
Submitted 21 October, 2021;
originally announced October 2021.
-
Dynamical stability of giant planets: the critical adiabatic index in the presence of a solid core
Authors:
Suman Kumar Kundu,
Eric R. Coughlin,
Andrew N. Youdin,
Philip J. Armitage
Abstract:
The dissociation and ionization of hydrogen, during the formation of giant planets via core accretion, reduces 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, whi…
▽ More
The dissociation and ionization of hydrogen, during the formation of giant planets via core accretion, reduces 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.
△ Less
Submitted 6 August, 2021;
originally announced August 2021.
-
Nonthermal filaments from the tidal destruction of clouds in the Galactic center
Authors:
Eric R. Coughlin,
C. J. Nixon,
Adam Ginsburg
Abstract:
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…
▽ More
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.
△ Less
Submitted 23 November, 2020; v1 submitted 26 October, 2020;
originally announced October 2020.
-
Structured, relativistic jets driven by radiation
Authors:
Eric R. Coughlin,
Mitchell C. Begelman
Abstract:
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 en…
▽ More
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.
△ Less
Submitted 28 September, 2020; v1 submitted 8 September, 2020;
originally announced September 2020.
-
The persistence of pancakes and the revival of self-gravity in tidal disruption events
Authors:
Eric R. Coughlin,
C. J. Nixon,
Patrick R. Miles
Abstract:
The destruction of a star by the tides of a supermassive black hole (SMBH) powers a bright accretion flare, and the theoretical modeling of such tidal disruption events (TDEs) can provide a direct means of inferring SMBH properties from observations. Previously it has been shown that TDEs with $β= r_{\rm t}/r_{\rm p} = 1$, where $r_{\rm t}$ is the tidal disruption radius and $r_{\rm p}$ is the per…
▽ More
The destruction of a star by the tides of a supermassive black hole (SMBH) powers a bright accretion flare, and the theoretical modeling of such tidal disruption events (TDEs) can provide a direct means of inferring SMBH properties from observations. Previously it has been shown that TDEs with $β= r_{\rm t}/r_{\rm p} = 1$, where $r_{\rm t}$ is the tidal disruption radius and $r_{\rm p}$ is the pericenter distance of the star, form an in-plane caustic, or ``pancake,'' where the tidally disrupted debris is compressed into a one-dimensional line within the orbital plane of the star. Here we show that this result applies generally to all TDEs for which the star is fully disrupted, i.e., that satisfy $β\gtrsim 1$. We show that the location of this caustic is always outside of the tidal disruption radius of the star and the compression of the gas near the caustic is at most mildly supersonic, which results in an adiabatic increase in the gas density above the tidal density of the black hole. As such, this in-plane pancake revitalizes the influence of self-gravity even for large $β$, in agreement with recent simulations. This finding suggests that for all TDEs in which the star is fully disrupted, self-gravity is revived post-pericenter, keeps the stream of debris narrowly confined in its transverse directions, and renders the debris prone to gravitational instability.
△ Less
Submitted 27 August, 2020;
originally announced August 2020.
-
Short gamma-ray bursts and the decompression of neutron star matter in tidal streams
Authors:
C. J. Nixon,
E. R. Coughlin,
J. E. Pringle
Abstract:
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 ma…
▽ More
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.
△ Less
Submitted 14 August, 2020;
originally announced August 2020.
-
Tidal Disruption Disks Formed and Fed by Stream-Stream and Stream-Disk Interactions in Global GRHD Simulations
Authors:
Zachary L. Andalman,
Matthew T. P. Liska,
Alexander Tchekhovskoy,
Eric R. Coughlin,
Nicholas Stone
Abstract:
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 ($β$=7) with realistic system p…
▽ More
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 ($β$=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.
△ Less
Submitted 23 November, 2021; v1 submitted 11 August, 2020;
originally announced August 2020.
-
Fallback Rates from Partial Tidal Disruption Events
Authors:
Patrick R. Miles,
Eric R. Coughlin,
C. J. Nixon
Abstract:
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. Anal…
▽ More
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.
△ Less
Submitted 16 June, 2020;
originally announced June 2020.
-
Emission from elliptical streams of dusty debris around white dwarfs
Authors:
C. J. Nixon,
J. E. Pringle,
E. R. Coughlin,
A. Swan,
J. Farihi
Abstract:
White dwarfs are routinely observed to have polluted atmospheres, and sometimes significant infrared excesses, that indicate ongoing accretion of circumstellar dust and rocky debris. Typically this debris is assumed to be in the form of a (circular) disc, and to originate from asteroids that passed close enough to the white dwarf to be pulled apart by tides. However, theoretical considerations sug…
▽ More
White dwarfs are routinely observed to have polluted atmospheres, and sometimes significant infrared excesses, that indicate ongoing accretion of circumstellar dust and rocky debris. Typically this debris is assumed to be in the form of a (circular) disc, and to originate from asteroids that passed close enough to the white dwarf to be pulled apart by tides. However, theoretical considerations suggest that the circularisation of the debris, which initially occupies highly eccentric orbits, is very slow. We therefore hypothesise that the observations may be readily explained by the debris remaining on highly eccentric orbits, and we explore the properties of such debris. For the generic case of an asteroid originating at several au from the white dwarf, we find that all of the tidal debris is always bound to the white dwarf and that the orbital energy distribution of the debris is narrow enough that it executes similar elliptical orbits with only a narrow spread. Assuming that the tidal field of the white dwarf is sufficient to minimise the effects of self-gravity and collisions within the debris, we estimate the time over which the debris spreads into a single elliptical ring, and we generate toy spectra and lightcurves from the initial disruption to late times when the debris distribution is essentially time steady. Finally we speculate on the connection between these simple considerations and the observed properties of these systems, and on additional physical processes that may change this simple picture.
△ Less
Submitted 13 June, 2020;
originally announced June 2020.
-
Variability in Short Gamma-ray Bursts: Gravitationally Unstable Tidal Tails
Authors:
Eric R. Coughlin,
C. J. Nixon,
Jennifer Barnes,
Brian D. Metzger,
R. Margutti
Abstract:
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…
▽ More
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 $γ\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 $γ$ equal), $γ=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.
△ Less
Submitted 23 June, 2020; v1 submitted 4 June, 2020;
originally announced June 2020.
-
The structure of nearly isothermal, adiabatic shockwaves
Authors:
Eric R. Coughlin
Abstract:
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 $γ\simeq 1$. Using the result that the post-shock gas is compressed into a shell of width $ΔR/R \simeq δ$, where $δ= γ-1$, we show that a choice of self-similar variable that exploits this compressive behavior in the limit tha…
▽ More
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 $γ\simeq 1$. Using the result that the post-shock gas is compressed into a shell of width $ΔR/R \simeq δ$, where $δ= γ-1$, we show that a choice of self-similar variable that exploits this compressive behavior in the limit that $γ\rightarrow 1$ naturally leads to a series expansion of the post-shock fluid density, pressure, and velocity in the small quantity $δ$. We demonstrate that the leading-order (in $δ$) solutions, which are increasingly accurate as $γ\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.
△ Less
Submitted 7 May, 2020;
originally announced May 2020.
-
A mildly relativistic outflow from the energetic, fast-rising blue optical transient CSS161010 in a dwarf galaxy
Authors:
D. L. Coppejans,
R. Margutti,
G. Terreran,
A. J. Nayana,
E. R. Coughlin,
T. Laskar,
K. D. Alexander,
M. Bietenholz,
D. Caprioli,
P. Chandra,
M. Drout,
D. Frederiks,
C. Frohmaier,
K. Hurley,
C. S. Kochanek,
M. MacLeod,
A. Meisner,
P. E. Nugent,
A. Ridnaia,
D. J. Sand,
D. Svinkin,
C. Ward,
S. Yang,
A. Baldeschi,
I. V. Chilingarian
, et al. (13 additional authors not shown)
Abstract:
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_ν\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 obse…
▽ More
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_ν\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 $Γβc\ge0.55c$ at ~100 days. This is faster than the non-relativistic AT2018cow ($Γβc\sim0.1c$) and closer to ZTF18abvkwla ($Γβ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 $<0.4\%$ of the local core-collapse supernova rate) H-rich transients that can launch mildly relativistic outflows.
△ Less
Submitted 23 May, 2020; v1 submitted 23 March, 2020;
originally announced March 2020.
-
The Gravitational Instability of Adiabatic Filaments
Authors:
Eric R. Coughlin,
C. J. Nixon
Abstract:
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…
▽ More
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^{σ_{\rm u}τ}$, where $τ$ is time relative to the sound crossing time across the diameter of the cylinder, and we derive the growth rate $σ_{\rm u}$ as a function of wavenumber. We find that there is a maximum growth rate $σ_{\rm max} \sim 1$ that occurs at a specific wavenumber $k_{\rm max} \sim 1$, and we derive the growth rate $σ_{\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.
△ Less
Submitted 17 February, 2020;
originally announced February 2020.
-
The Influence of Black Hole Binarity on Tidal Disruption Events
Authors:
Eric R. Coughlin,
Philip J. Armitage,
Giuseppe Lodato,
C. J. Nixon
Abstract:
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 companio…
▽ More
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.
△ Less
Submitted 27 September, 2019;
originally announced September 2019.
-
On the Diversity of Fallback Rates from Tidal Disruption Events with Accurate Stellar Structure
Authors:
Elen C. A. Golightly,
C. J. Nixon,
Eric R. Coughlin
Abstract:
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 outcom…
▽ More
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 $γ= 4/3$ and $γ= 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 $γ= 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.
△ Less
Submitted 5 August, 2019; v1 submitted 12 July, 2019;
originally announced July 2019.
-
Partial Stellar Disruption by a Supermassive Black Hole: Is the Lightcurve Really Proportional to $t^{-9/4}$?
Authors:
Eric R. Coughlin,
C. J. Nixon
Abstract:
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…
▽ More
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.
△ Less
Submitted 4 September, 2019; v1 submitted 5 July, 2019;
originally announced July 2019.
-
Energy-conserving, Relativistic Corrections to Strong Shock Propagation
Authors:
Eric R. Coughlin
Abstract:
Astrophysical explosions are accompanied by the propagation of a shock wave through an ambient medium. Depending on the mass and energy involved in the explosion, the shock velocity $V$ can be non-relativistic ($V \ll c$, where $c$ is the speed of light), ultra-relativistic ($V \simeq c$), or moderately relativistic ($V \sim few\times 0.1c$). While self-similar, energy-conserving solutions to the…
▽ More
Astrophysical explosions are accompanied by the propagation of a shock wave through an ambient medium. Depending on the mass and energy involved in the explosion, the shock velocity $V$ can be non-relativistic ($V \ll c$, where $c$ is the speed of light), ultra-relativistic ($V \simeq c$), or moderately relativistic ($V \sim few\times 0.1c$). While self-similar, energy-conserving solutions to the fluid equations that describe the shock propagation are known in the non-relativistic (the Sedov-Taylor blastwave) and ultra-relativistic (the Blandford-McKee blastwave) regimes, the finite speed of light violates scale invariance and self-similarity when the flow is only mildly relativistic. By treating relativistic terms as perturbations to the fluid equations, here we derive the $\mathcal{O}(V^2/c^2)$, energy-conserving corrections to the non-relativistic, Sedov-Taylor solution for the propagation of a strong shock. We show that relativistic terms modify the post-shock fluid velocity, density, pressure, and the shock speed itself, the latter being constrained by global energy conservation. We derive these corrections for a range of post-shock adiabatic indices $γ$ (which we set as a fixed number for the post-shock gas) and ambient power-law indices $n$, where the density of the ambient medium $ρ_{\rm a}$ into which the shock advances declines with spherical radius $r$ as $ρ_{\rm a} \propto r^{-n}$. For Sedov-Taylor blastwaves that terminate in a contact discontinuity with diverging density, we find that there is no relativistic correction to the Sedov-Taylor solution that simultaneously satisfies the fluid equations and conserves energy. These solutions have implications for relativistic supernovae, the transition from ultra- to sub-relativistic velocities in gamma-ray bursts, and other high-energy phenomena.
△ Less
Submitted 28 May, 2019;
originally announced May 2019.
-
Ultra-deep tidal disruption events: prompt self-intersections and observables
Authors:
Siva Darbha,
Eric R. Coughlin,
Daniel Kasen,
Chris Nixon
Abstract:
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 $β$ vs SMBH masses $M$ that produce…
▽ More
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 $β$ 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 $β\sim 50 - 127$ for $M/M_\odot = 10^4$, down to $β\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.
△ Less
Submitted 16 May, 2019;
originally announced May 2019.
-
Weak Shock Propagation with Accretion. III. A Numerical Study On Shock Propagation & Stability
Authors:
Stephen Ro,
Eric R. Coughlin,
Eliot Quataert
Abstract:
Core-collapse supernovae span a wide range of energies, from much less than to much greater than the binding energy of the progenitor star. As a result, the shock wave generated from a supernova explosion can have a wide range of Mach numbers. In this paper, we investigate the propagation of shocks with arbitrary initial strengths in polytropic stellar envelopes using a suite of spherically symmet…
▽ More
Core-collapse supernovae span a wide range of energies, from much less than to much greater than the binding energy of the progenitor star. As a result, the shock wave generated from a supernova explosion can have a wide range of Mach numbers. In this paper, we investigate the propagation of shocks with arbitrary initial strengths in polytropic stellar envelopes using a suite of spherically symmetric hydrodynamic simulations. We interpret these results using the three known self-similar solutions for this problem: the Sedov-Taylor blastwave describes an infinitely strong shock and the self-similar solutions from Coughlin et al. (2018b) (Paper I) and Coughlin et al. (2019) (Paper II) describe a weak and infinitely weak shock (the latter being a rarefaction wave). We find that shocks, no matter their initial strengths, evolve toward either the infinitely strong or infinitely weak self-similar solutions at sufficiently late times. For a given density profile, a single function characterizes the long-term evolution of a shock's radius and strength. However, shocks with strengths near the self-similar solution for a weak shock (from Paper I) evolve extremely slowly with time. Therefore, the self-similar solutions for infinitely strong and infinitely weak shocks are not likely to be realized in low-energy stellar explosions, which will instead retain memory of the shock strength initiated in the stellar interior.
△ Less
Submitted 29 April, 2019; v1 submitted 2 April, 2019;
originally announced April 2019.
-
Thawing the frozen-in approximation: implications for self-gravity in deeply plunging tidal disruption events
Authors:
Elad Steinberg,
Eric R. Coughlin,
Nicholas C. Stone,
Brian D. Metzger
Abstract:
The tidal destruction of a star by a massive black hole, known as a tidal disruption event (TDE), is commonly modeled using the "frozen-in" approximation. Under this approximation, the star maintains exact hydrostatic balance prior to entering the tidal sphere (radius $r_{\rm t}$), after which point its internal pressure and self-gravity become instantaneously negligible and the debris undergoes b…
▽ More
The tidal destruction of a star by a massive black hole, known as a tidal disruption event (TDE), is commonly modeled using the "frozen-in" approximation. Under this approximation, the star maintains exact hydrostatic balance prior to entering the tidal sphere (radius $r_{\rm t}$), after which point its internal pressure and self-gravity become instantaneously negligible and the debris undergoes ballistic free fall. We present a suite of hydrodynamical simulations of TDEs with high penetration factors $β\equiv r_{\rm t}/r_{\rm p} = 5-7$, where $r_{\rm p}$ is the pericenter of the stellar center of mass, calculated using a Voronoi-based moving-mesh technique. We show that basic assumptions of the frozen-in model, such as the neglect of self-gravity inside $r_{\rm t}$, are violated. Indeed, roughly equal fractions of the final energy spread accumulate exiting and entering the tidal sphere, though the frozen-in prediction is correct at the order-of-magnitude level. We also show that an $\mathcal{O}(1)$ fraction of the debris mass remains transversely confined by self-gravity even for large $β$ which has implications for the radio emission from the unbound debris and, potentially, for the circularization efficiency of the bound streams.
△ Less
Submitted 2 April, 2019; v1 submitted 9 March, 2019;
originally announced March 2019.
-
Weak Shock Propagation with Accretion. II. Stability of Self-similar Solutions to Radial Perturbations
Authors:
Eric R. Coughlin,
Stephen Ro,
Eliot Quataert
Abstract:
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 s…
▽ More
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^α$, where $α\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.
△ Less
Submitted 14 January, 2019;
originally announced January 2019.
-
Black Hole Accretion Discs and Luminous Transients in Failed Supernovae from Non-Rotating Supergiants
Authors:
Eliot Quataert,
Daniel Lecoanet,
Eric R. Coughlin
Abstract:
We show that for supergiants, net angular momentum is not a necessary condition for forming accretion discs during core-collapse. Even absent net rotation, convective motions in the outer parts of supergiants generate mean horizontal flows at a given radius with velocities of order 1 km/s; the direction of the mean flow will vary as a function of height through the convection zone. We confirm thes…
▽ More
We show that for supergiants, net angular momentum is not a necessary condition for forming accretion discs during core-collapse. Even absent net rotation, convective motions in the outer parts of supergiants generate mean horizontal flows at a given radius with velocities of order 1 km/s; the direction of the mean flow will vary as a function of height through the convection zone. We confirm these analytic estimates using Cartesian Boussinesq convection simulations. These mean horizontal flows lead to a random angular momentum in supergiant convection zones that exceeds that of the last stable circular orbit of a black hole by a factor of order 10. As a result, failed explosions of supergiants - in which the accretion shock onto the neutron star does not revive, leading to black hole formation - may often produce accretion discs that can power day-week (blue supergiants) or week-year (yellow and red supergiants) non-thermal and thermal transients through winds and jets. These transients will be especially time variable because the angular momentum of the accreting material will vary substantially in time. Observed sources such as Swift J1644+57, iPTF14hls, and SN 2018cow, as well as energetic Type II supernovae (OGLE-2014-SN-073) may be produced by this mechanism.
△ Less
Submitted 29 November, 2018;
originally announced November 2018.