Researcher profile

Eric R. Coughlin

Eric R. Coughlin contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2026arXiv

The Maximum Gravity Model for partial Tidal Disruption Events: Mass Loss, Peak Fallback Rate and Dependence on Stellar Properties

A star entering the tidal sphere of a supermassive black hole (SMBH) can be partially stripped of mass, resulting in a partial tidal disruption event (TDE). Here we develop an analytical model for properties of these events, including the peak fallback rate, $\dot{M}_{\rm peak}$, the time at which the peak occurs, $t_{\rm peak}$, and the amount of mass removed from the star, $ΔM$, for any star and any pericenter distance associated with the stellar orbit about the black hole. We compare the model predictions to 1276 hydrodynamical simulations of partial TDEs of main-sequence stars by a $10^6 M_\odot$ SMBH. The model yields $t_{\rm peak}$ predictions that are in good agreement (to within tens of percent) with the numerical simulations for any stellar mass and age. The agreement for $\dot{M}_{\rm peak}$ is weaker due to the influence of self-gravity on the debris stream dynamics, which remains dynamically important for partial TDEs; the agreement for $\dot{M}_{\rm peak}$ is, however, to within a factor of $\sim 2-3$ in the majority of cases considered, with larger differences for low-mass stars ($M_\star \lesssim 0.5 M_\odot$) on grazing orbits with small mass loss. We show that partial TDE lightcurves for disruptions caused by $\sim 10^6M_\odot$ SMBHs can span $\sim 20-100$ day peak timescales, whereas grazing encounters of high-mass stars with high-mass SMBHs can yield longer peak timescales ($t\gtrsim 1000$ days), associated with some observed transients. Our model provides a significant step toward an analytical prescription for TDE lightcurves and luminosity functions.

preprint2023arXiv

Tidal capture of stars by supermassive black holes: implications for periodic nuclear transients and quasi-periodic eruptions

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).

preprint2022arXiv

On the impact of relativistic gravity on the rate of tidal disruption events

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}$.

preprint2022arXiv

Spherically Symmetric Accretion onto a Compact Object through a Standing Shock: The Effects of General Relativity in the Schwarzschild Geometry

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.

preprint2022arXiv

Stellar Revival and Repeated Flares in Deeply Plunging Tidal Disruption Events

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 &#34;deep encounters&#34; 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.

preprint2022arXiv

Using the Hills Mechanism to Generate Repeating Partial Tidal Disruption Events and ASASSN-14ko

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.

preprint2021arXiv

A Physical Model of Delayed Rebrightenings in Shock-interacting Supernovae without Narrow Line Emission

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.

preprint2021arXiv

Stars Crushed by Black Holes. I. On the Energy Distribution of Stellar Debris in Tidal Disruption Events

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.

preprint2021arXiv

Stars Crushed by Black Holes. II. A Physical Model of Adiabatic Compression and Shock Formation in Tidal Disruption Events

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.

preprint2020arXiv

Fallback Rates from Partial Tidal Disruption Events

A tidal disruption event (TDE) occurs when a star plunges through a supermassive black hole&#39;s tidal radius, at which point the star&#39;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&#39;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.

preprint2020arXiv

The Gravitational Instability of Adiabatic Filaments

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.

preprint2020arXiv

The persistence of pancakes and the revival of self-gravity in tidal disruption events

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,&#39;&#39; 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.

preprint2020arXiv

The structure of nearly isothermal, adiabatic shockwaves

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.

preprint2020arXiv

Variability in Short Gamma-ray Bursts: Gravitationally Unstable Tidal Tails

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 &#34;tails&#34; 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.