Source author record

James Guillochon

James Guillochon appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

34works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

34 published item(s)

preprint2020arXiv

Stellar Tidal Disruption Events with Abundances and Realistic Structures (STARS): Library of Fallback Rates

We present the STARS library, a grid of tidal disruption event (TDE) simulations interpolated to provide the mass fallback rate ($dM/dt$) to the black hole for a main-sequence star of any stellar mass, stellar age, and impact parameter. We use a one-dimensional stellar evolution code to construct stars with accurate stellar structures and chemical abundances, then perform tidal disruption simulations in a three-dimensional adaptive-mesh hydrodynamics code with a Helmholtz equation of state, in unprecedented resolution: from 131 to 524 cells across the diameter of the star. The interpolated library of fallback rates is available on GitHub (https://github.com/jamielaw-smith/STARS_library) and version 1.0.0 is archived on Zenodo; one can query the library for any stellar mass, stellar age, and impact parameter. We provide new fitting formulae for important disruption quantities ($β_{\rm crit}, ΔM, \dot M_{\rm peak}, t_{\rm peak}, n_\infty$) as a function of stellar mass, stellar age, and impact parameter. Each of these quantities vary significantly with stellar mass and stellar age, but we are able to reduce all of our simulations to a single relationship that depends only on stellar structure, characterized by a single parameter $ρ_c/\barρ$, and impact parameter $β$. We also find that, in general, more centrally concentrated stars have steeper $dM/dt$ rise slopes and shallower decay slopes. For the same $ΔM$, the $dM/dt$ shape varies significantly with stellar mass, promising the potential determination of stellar properties from the TDE light curve alone. The $dM/dt$ shape depends strongly on stellar structure and to a certain extent stellar mass, meaning that fitting TDEs using this library offers a better opportunity to determine the nature of the disrupted star and the black hole.

preprint2017arXiv

Low-mass White Dwarfs with Hydrogen Envelopes as a Missing Link in the Tidal Disruption Menu

We construct a menu of objects that can give rise to bright flares when disrupted by massive black holes (BHs), ranging from planets to evolved stars. Through their tidal disruption, main sequence and evolved stars can effectively probe the existence of otherwise quiescent supermassive BHs and white dwarfs can probe intermediate mass BHs. Many low-mass white dwarfs possess extended hydrogen envelopes, which allow for the production of prompt flares in disruptive encounters with moderately massive BHs of $10^5$ to $10^7~M_\odot$--masses that may constitute the majority of massive BHs by number. These objects are a missing link in two ways: (1) for probing moderately massive BHs and (2) for understanding the hydrodynamics of the disruption of objects with tenuous envelopes. A flare arising from the tidal disruption of a $0.17~M_\odot$ white dwarf by a $10^5~M_\odot$ BH reaches a maximum between 0.6 and 11 days, with a peak fallback rate that is usually super-Eddington and results in a flare that is likely brighter than a typical tidal disruption event. Encounters stripping only the envelope can provide hydrogen-only fallback, while encounters disrupting the core evolve from H- to He-rich fallback. While most tidal disruption candidates observed thus far are consistent with the disruptions of main sequence stars, the rapid timescales of nuclear transients such as Dougie and PTF10iya are naturally explained by the disruption of low-mass white dwarfs. As the number of observed flares continues to increase, the menu presented here will be essential for characterizing nuclear BHs and their environments through tidal disruptions.

preprint2016arXiv

An Ultraviolet Spectrum of the Tidal Disruption Flare ASASSN-14li

We present a Hubble Space Telescope STIS spectrum of ASASSN-14li, the first rest-frame UV spectrum of a tidal disruption flare (TDF). The underlying continuum is well fit by a blackbody with $T_{\mathrm{UV}} = 3.5 \times 10^{4}$ K, an order of magnitude smaller than the temperature inferred from X-ray spectra (and significantly more precise than previous efforts based on optical and near-UV photometry). Super-imposed on this blue continuum, we detect three classes of features: narrow absorption from the Milky Way (probably a high-velocity cloud), and narrow absorption and broad (FWHM $\approx 2000$-8000 km s$^{-1}$) emission lines at/near the systemic host velocity. The absorption lines are blueshifted with respect to the emission lines by $Δv = -(250$-400) km s$^{-1}$. Due both to this velocity offset and the lack of common low-ionization features (Mg II, Fe II), we argue these arise from the same absorbing material responsible for the low-velocity outflow discovered at X-ray wavelengths. The broad nuclear emission lines display a remarkable abundance pattern: N III], N IV], He II are quite prominent, while the common quasar emission lines of C III] and Mg II are weak or entirely absent. Detailed modeling of this spectrum will help elucidate fundamental questions regarding the nature of the emission process(es) at work in TDFs, while future UV spectroscopy of ASASSN-14li would help to confirm (or refute) the previously proposed connection between TDFs and "N-rich" quasars.

preprint2016arXiv

Discovery of an outflow from radio observations of the tidal disruption event ASASSN-14li

We report the discovery of transient radio emission from the nearby optically-discovered TDE ASASSN-14li (distance of 90 Mpc), making it the first typical TDE detected in the radio, and unambiguously pointing to the formation of a non-relativistic outflow with a kinetic energy of $\approx (4-10)\times10^{47}$ erg, a velocity of $\approx 12,000-36,000$ km s$^{-1}$, and a mass of $\approx 3\times10^{-5}-7\times10^{-4}$ M$_{\odot}$. We show that the outflow was ejected on 2014 August 11--25, in agreement with an independent estimate of the timing of super-Eddington accretion based on the optical, UV, and X-ray observations, and that the ejected mass corresponds to about $1-10\%$ of the mass accreted in the super-Eddington phase. The temporal evolution of the radio emission also uncovers the circumnuclear density profile, $ρ(R)\propto R^{-2.5}$ on a scale of about 0.01 pc, a scale that cannot be probed via direct measurements even in the nearest SMBHs. Our discovery of radio emission from the nearest well-studied TDE to date, with a radio luminosity lower than all previous limits, indicates that non-relativistic outflows are ubiquitous in TDEs, and that future, more sensitive, radio surveys will uncover similar events.

preprint2016arXiv

Prompt Radiation and Mass Outflows from the Stream-Stream Collisions of Tidal Disruption Events

Stream-stream collisions play an important role for the circularization of highly eccentric streams resulting from tidal disruption events (TDEs). We perform three dimensional radiation hydrodynamic simulations to show that stream collisions can contribute significant optical and ultraviolet light to the flares produced by TDEs, and can sometimes explain the majority of the observed emission. Our simulations focus on the region near the radiation pressure dominated shock produced by a collision and track how the kinetic energy of the stream is dissipated by the associated shock. When the mass flow rate of the stream $\dot{M}$ is a significant fraction of the Eddington accretion rate, $\gtrsim2\%$ of the initial kinetic energy is converted to radiation directly as a result of the collision. In this regime, the collision redistributes the specific kinetic energy into the downstream gas and more than $16\%$ of the mass can become unbound. The fraction of unbound gas decreases rapidly as $\dot{M}$ drops significantly below the Eddington limit, with no unbound gas being produced when $\dot{M}$ drops to $1\%$ of Eddington; we find however that the radiative efficiency increases slightly to $\lesssim 8\%$ in these low $\dot{M}$ cases. The effective radiation temperature and size of the photosphere is determined by the stream velocity and $\dot{M}$, which we find to be a few times $10^4$~K and $10^{14}$~cm in our calculations, comparable to the inferred values of some TDE candidates. The photosphere size is directly proportional to $\dot{M}$, which can explain the rapidly changing photosphere sizes seen in TDE candidates such as PS1-10jh.

preprint2016arXiv

The Astropy Problem

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

preprint2016arXiv

The X-ray through Optical Fluxes and Line Strengths of Tidal Disruption Events

Observations of luminous flares resulting from the possible tidal disruption of stars by supermassive black holes have raised a number of puzzles. Outstanding questions include the origin of the optical and ultraviolet (UV) flux, the weakness of hydrogen lines in the spectrum, and the occasional simultaneous observation of x-rays. Here we study the emission from tidal disruption events (TDEs) produced as radiation from black hole accretion propagates through an extended, optically thick envelope formed from stellar debris. We analytically describe key physics controlling spectrum formation, and present detailed radiative transfer calculations that model the spectral energy distribution (SED) and optical line strengths of TDEs near peak brightness. The steady-state transfer is coupled to a non local thermodynamic equilibrium treatment of the excitation and ionization states of hydrogen, helium and oxygen (as a representative metal). Our calculations show how an extended envelope can reprocess a fraction of soft x-rays and produce the observed optical fluxes of order 10^43 ergs per second. Variations in the mass or size of the envelope may help explain how the optical flux changes over time with roughly constant color. For high enough accretion luminosities, x-rays can highly ionize the reprocessing region and escape to be observed simultaneously with the optical flux, producing an SED not described by a single blackbody. Due to optical depth effects, hydrogen Balmer line emission is often strongly suppressed relative to helium line emission (with HeII-to-H line ratios of at least 5:1 in some cases) even in the disruption of a solar-composition star. We discuss the implications of our results to understanding the type of stars destroyed in TDEs and the physical processes responsible for producing the observed flares.

preprint2015arXiv

A Dark Year for Tidal Disruption Events

The disruption of a main-sequence star by a supermassive black hole results in the initial production of an extended debris stream that winds repeatedly around the black hole, producing a complex three-dimensional figure that may self-intersect. Both analytical work and simulations have shown that typical encounters generate streams that are extremely thin. In this paper we show that this implies that even small relativistic precessions attributed to black hole spin can induce deflections that prevent the stream from self-intersecting even after many windings. Additionally, hydrodynamical simulations have demonstrated that energy is deposited very slowly via hydrodynamic processes alone, resulting in the liberation of very little gravitational binding energy in the absence of stream-stream collisions. This naturally leads to a "dark period" in which the flare is not observable for some time, persisting for up to a dozen orbital periods of the most bound material, which translates to years for disruptions around black holes with mass $\sim 10^{7} M_{\odot}$. We find that more-massive black holes tend to have more violent stream self-intersections, resulting in short viscous times that lead to prompt accretion onto the black hole. For these tidal disruption events (TDEs), the accretion rate onto the black hole should still closely follow the original fallback rate after a fixed delay time $t_{\rm delay}$. For lower black hole masses ($M_{\rm h} \lesssim 10^{6}$), we find that flares are typically slowed down by about an order of magnitude, and because the accretion rates for TDEs about higher-mass black hole are already sub-Eddington, this results in the majority of TDEs being sub-Eddington at peak. This also implies that current searches for TDEs are biased towards prompt flares, with slowed flares likely having been unidentified. [abridged]

preprint2015arXiv

Death by Dynamics: Planetoid-Induced Explosions on White Dwarfs

At intervals as short as ten thousand years, each white dwarf (WD) passes within a solar radius of a planetoid, i.e., a comet, asteroid, or planet. Gravitational tidal forces tear the planetoid apart; its metal-rich debris falls onto the WD, enriching the atmosphere. A third of WDs exhibit atmospheric "pollution". For roughly every hundred planetoid disruptions, a planetoid collides with a WD. We simulate a small number of collisions, in which "death-by-dynamics" refers to the fate of the planetoid. We also compute the energies and likely durations of a broad sample of collision events, and identify detection strategies at optical and X-ray wavelengths. Collisions with the most massive planetoids can be detected in external galaxies. Some may trigger nuclear burning. If one in $\sim 10^7-10^8$ of WD-planetoid collisions creates the conditions needed for a Type Ia supernova (SN~Ia), "death-by-dynamics" would also refer to the fate of the WD, and could provide a novel channel for the production of SN~Ia. We consider the circumstances under which the rate of SNe~Ia can be increased by interactions with planetoids.

preprint2015arXiv

Flows of X-ray gas reveal the disruption of a star by a massive black hole

Tidal forces close to massive black holes can violently disrupt stars that make a close approach. These extreme events are discovered via bright X-ray and optical/UV flares in galactic centers. Prior studies based on modeling decaying flux trends have been able to estimate broad properties, such as the mass accretion rate. Here we report the detection of flows of highly ionized X-ray gas in high-resolution X-ray spectra of a nearby tidal disruption event. Variability within the absorption-dominated spectra indicates that the gas is relatively close to the black hole. Narrow line widths indicate that the gas does not stretch over a large range of radii, giving a low volume filling factor. Modest outflow speeds of a few hundred kilometers per second are observed, significantly below the escape speed from the radius set by variability. The gas flow is consistent with a rotating wind from the inner, super-Eddington region of a nascent accretion disk, or with a filament of disrupted stellar gas near to the apocenter of an elliptical orbit. Flows of this sort are predicted by fundamental analytical theory and more recent numerical simulations.

preprint2015arXiv

Hydrodynamical Simulations to Determine the Feeding Rate of Black Holes by the Tidal Disruption of Stars: The Importance of the Impact Parameter and Stellar Structure

The disruption of stars by supermassive black holes has been linked to more than a dozen flares in the cores of galaxies out to redshift $z \sim 0.4$. Modeling these flares properly requires a prediction of the rate of mass return to the black hole after a disruption. Through hydrodynamical simulation, we show that aside from the full disruption of a solar mass star at the exact limit where the star is destroyed, the common assumptions used to estimate $\dot{M}(t)$, the rate of mass return to the black hole, are largely invalid. While the analytical approximation to tidal disruption predicts that the least-centrally concentrated stars and the deepest encounters should have more quickly-peaked flares, we find that the most-centrally concentrated stars have the quickest-peaking flares, and the trend between the time of peak and the impact parameter for deeply-penetrating encounters reverses beyond the critical distance at which the star is completely destroyed. We also show that the most-centrally concentrated stars produced a characteristic drop in $\dot{M}(t)$ shortly after peak when a star is only partially disrupted, with the power law index $n$ being as extreme as -4 in the months immediately following the peak of a flare. Additionally, we find that $n$ asymptotes to $\simeq -2.2$ for both low- and high-mass stars for approximately half of all stellar disruptions. Both of these results are significantly steeper than the typically assumed $n = -5/3$. As these precipitous decay rates are only seen for events in which a stellar core survives the disruption, they can be used to determine if an observed tidal disruption flare produced a surviving remnant. These results should be taken into consideration when flares arising from tidal disruptions are modeled. [abridged]

preprint2015arXiv

Observational Cosmology With Semi-Relativistic Stars

Galaxy mergers lead to the formation of massive black hole binaries which can accelerate background stars close to the speed of light. We estimate the comoving density of ejected stars with a peculiar velocity in excess of $0.1c$ or $0.5c$ to be $\sim 10^{10}$ and $10^5$ Gpc$^{-3}$ respectively, in the present-day Universe. Semi-relativistic giant stars will be detectable with forthcoming telescopes out to a distance of a few Mpc, where their proper motion, radial velocity, and age, can be spectroscopically measured. In difference from traditional cosmological messengers, such as photons, neutrinos, or cosmic-rays, these stars shine and so their trajectories need not be directed at the observer for them to be detected. Tracing the stars to their parent galaxies as a function of speed and age will provide a novel test of the equivalence principle and the standard cosmological parameters. Semi-relativistic stars could also flag black hole binaries as gravitational wave sources for the future eLISA observatory.

preprint2015arXiv

Relativistic effects on tidal disruption kicks of solitary stars

Solitary stars that wander too close to their galactic centres can become tidally disrupted, if the tidal forces due to the supermassive black hole (SMBH) residing there overcome the self-gravity of the star. If the star is only partially disrupted, so that a fraction survives as a self-bound object, this remaining core will experience a net gain in specific orbital energy, which translates into a velocity "kick" of up to $\sim 10^3$ km/s. In this paper, we present the result of smoothed particle hydrodynamics (SPH) simulations of such partial disruptions, and analyse the velocity kick imparted on the surviving core. We compare $γ$ = 5/3 and $γ$ = 4/3 polytropes disrupted in both a Newtonian potential, and a generalized potential that reproduces most relativistic effects around a Schwarzschild black hole either exactly or to excellent precision. For the Newtonian case, we confirm the results of previous studies that the kick velocity of the surviving core is virtually independent of the ratio of the black hole to stellar mass, and is a function of the impact parameter $β$ alone, reaching at most the escape velocity of the original star. For a given $β$, relativistic effects become increasingly important for larger black hole masses. In particular, we find that the kick velocity increases with the black hole mass, making larger kicks more common than in the Newtonian case, as low-$β$ encounters are statistically more likely than high-$β$ encounters. The analysis of the tidal tensor for the generalized potential shows that our results are robust lower limits on the true relativistic kick velocities, and are generally in very good agreement with the exact results.

preprint2015arXiv

SETI via Leakage from Light Sails in Exoplanetary Systems

The primary challenge of rocket propulsion is the burden of needing to accelerate the spacecraft's own fuel, resulting in only a logarithmic gain in maximum speed as propellant is added to the spacecraft. Light sails offer an attractive alternative in which fuel is not carried by the spacecraft, with acceleration being provided by an external source of light. By artificially illuminating the spacecraft with beamed radiation, speeds are only limited by the area of the sail, heat resistance of its material, and power use of the accelerating apparatus. In this paper, we show that leakage from a light sail propulsion apparatus in operation around a solar system analogue would be detectable. To demonstrate this, we model the launch and arrival of a microwave beam-driven light sail constructed for transit between planets in orbit around a single star, and find an optimal beam frequency on the order of tens of GHz. Leakage from these beams yields transients with flux densities of Jy and durations of tens of seconds at 100 pc. Because most travel within a planetary system would be conducted between the habitable worlds within that system, multiply-transiting exoplanetary systems offer the greatest chance of detection, especially when the planets are in projected conjunction as viewed from Earth. If interplanetary travel via beam-driven light sails is commonly employed in our galaxy, this activity could be revealed by radio follow-up of nearby transiting exoplanetary systems. The expected signal properties define a new strategy in the search for extraterrestrial intelligence (SETI).

preprint2015arXiv

The Fastest Unbound Stars in the Universe

The discovery of hypervelocity stars (HVS) leaving our galaxy with speeds of nearly $10^{3}$ km s$^{-1}$ has provided strong evidence towards the existence of a massive compact object at the galaxy's center. HVS ejected via the disruption of stellar binaries can occasionally yield a star with $v_{\infty} \lesssim 10^4$ km s$^{-1}$, here we show that this mechanism can be extended to massive black hole (MBH) mergers, where the secondary star is replaced by a MBH with mass $M_2 \gtrsim 10^5 M_{\odot}$. We find that stars that are originally bound to the secondary MBH are frequently ejected with $v_{\infty} > 10^4$ km s$^{-1}$, and occasionally with velocities $\sim 10^5$ km s$^{-1}$ (one third the speed of light), for this reason we refer to stars ejected from these systems as "semi-relativistic" hypervelocity stars (SHS). Bound to no galaxy, the velocities of these stars are so great that they can cross a significant fraction of the observable universe in the time since their ejection (several Gpc). We demonstrate that if a significant fraction of MBH mergers undergo a phase in which their orbital eccentricity is $\gtrsim 0.5$ and their periapse distance is tens of the primary's Schwarzschild radius, the space density of fast-moving ($v_{\infty} > 10^{4}$ km s$^{-1}$) SHS may be as large as $10^{3}$ Mpc$^{-3}$. Hundreds of the SHS will be giant stars that could be detected by future all-sky infrared surveys such as WFIRST or Euclid and proper motion surveys such as LSST, with spectroscopic follow-up being possible with JWST.

preprint2015arXiv

The Gamma-ray Afterglows of Tidal Disruption Events

A star wandering too close to a supermassive black hole (SMBH) will be tidally disrupted. Previous studies of such "tidal disruption event" (TDE) mostly focus on the stellar debris that are bound to the system, because they give rise to luminous flares. On the other hand, half of the stellar debris in principle are unbound and can stream to a great distance, but so far there is no clear evidence that this "unbound debris stream" (UDS) exists. Motivated by the fact that the circum-nuclear region around SMBHs is usually filled with dense molecular clouds (MCs), here we investigate the observational signatures resulting from the collision between an UDS and a MC, which is likely to happen hundreds of years after a TDE. We focus on $γ$-ray emission ($0.1-10^5$ GeV), which comes from the encounter of shock-accelerated cosmic rays with background protons and, more importantly, is not subject to extinction. We show that because of the high proton density inside a MC, the peak $γ$-ray luminosity, about $10^{39}~{\rm erg~s^{-1}}$, is at least $100$ times greater than that in the case without a MC (only with a smooth interstellar medium). The luminosity decays on a time-scale of decades, depending on the distance of the MC, and about a dozen of these "TDE afterglows" could be detected within a distance of about $16$ Mpc by the future Cherenkov Telescope Array. Without careful discrimination, these sources potentially could contaminate the searches for starburt galaxies, galactic nuclei containing millisecond pulsars or dark-matter annihilation signals.

preprint2015arXiv

Thermonuclear detonations ensuing white dwarf mergers

The merger of two white dwarfs (WDs) has for many years not been considered as the favoured model for the progenitor system of type Ia supernovae (SNe Ia). But recent years have seen a change of opinion as a number of studies, both observational and theoretical, have concluded that they should contribute significantly to the observed type Ia supernova rate. In this paper, we study the ignition and propagation of detonation through post-merger remnants and we follow the resulting nucleosynthesis up to the point where a homologous expansion is reached. In our study we cover the entire range of WD masses and compositions. For the emergence of a detonation we study three different setups. The first two are guided by the merger remnants from our earlier simulations (Dan et al. 2014}, while for the third one the ignitions were set by placing hotspots with properties determined by spatially resolved calculations taken from the literature. There are some caveats to our approach which we investigate. We carefully compare the nucleosynthetic yields of successful explosions with SN Ia observations. Only three of our models are consistent with all the imposed constraints and potentially lead to a standard type Ia event. The first one, a $0.45\,M_\odot\ {\rm helium\ (He)} + 0.9 \,M_\odot$ carbon-oxygen (CO) WD system produces a sub-luminous, SN 1991bg-like event while the other two, a $0.45\,M_\odot\ {\rm He} + 1.1\,M_\odot$ oxygen-neon (ONe) WD system and a $1.05+1.05\,M_\odot$ system with two CO WDs, are good candidates for common SNe Ia.

preprint2015arXiv

Unbound Debris Streams and Remnants Resulting From the Tidal Disruptions of Stars by Supermassive Black Holes

The kinetic energy of a star in orbit about a supermassive black hole is a significant fraction of its rest mass energy when its periapse is comparable to its tidal radius. Upon its destruction, a fraction of this energy is extracted and injected into the stellar debris, half of which becomes unbound from the black hole, with the fastest material moving at $\sim 0.03 c$. In this paper, we present a formalism for determining the fate of these unbound debris streams (UDSs) as they depart from the black hole and interact with the surrounding gas. As the density and velocity varies along the length of a UDS, we find that hydrodynamical drag quickly shapes UDSs into loop-like structures, with the densest portions of the streams leading portions of lower density. As UDSs travel outwards, their drag against the ISM increases quadratically with distance, which causes UDSs to deposit their momentum and energy into the ambient medium before the surrounding shocked ISM has a chance to cool. This sudden injection of $\sim 10^{50}$ erg into the ambient medium generates a Sedov-like unbound debris remnant (UDR) that mimics supernova remnants (SNRs) in energetics and appearance, accelerates particles which will produce cosmic rays and synchrotron emission, and provides momentum feedback into the molecular clouds surrounding a black hole. We estimate that a few of these UDRs might be present within a couple degrees of the Galactic Center masquerading as SNRs, and that the UDR scenario is a plausible explanation for Sgr A East.

preprint2014arXiv

Illuminating Massive Black Holes With White Dwarfs: Orbital Dynamics and High Energy Transients from Tidal Interactions

White dwarfs (WDs) can be tidally disrupted only by massive black holes (MBHs) with masses less than $\sim10^5 M_\odot$. These tidal interactions feed material to the MBH well above its Eddington limit, with the potential to launch a relativistic jet. The corresponding beamed emission is a promising signpost to an otherwise quiescent MBH of relatively low mass. We show that the mass transfer history, and thus the lightcurve, are quite different when the disruptive orbit is parabolic, eccentric, or circular. The mass lost each orbit exponentiates in the eccentric-orbit case leading to the destruction of the WD after several tens of orbits. We examine the stellar dynamics of clusters surrounding MBHs to show that single-passage WD disruptions are substantially more common than repeating encounters. The $10^{49}$ erg s$^{-1}$ peak luminosity of these events makes them visible to cosmological distances. They may be detectible at rates of as many as tens per year by instruments like Swift. In fact, WD-disruption transients significantly outshine their main-sequence star counterparts, and are the most likely tidal interaction to be detected arising from MBHs with masses less than $10^5 M_\odot$. The detection or non-detection of such WD-disruption transients by Swift is, therefore, a powerful tool to constrain lower end of the MBH mass function. The emerging class of ultra-long gamma ray bursts all have peak luminosities and durations reminiscent of WD disruptions, offering a hint that WD-disruption transients may already be present in existing datasets.

preprint2014arXiv

Possible Origin of the G2 Cloud from the Tidal Disruption of a Known Giant Star by Sgr A*

The discovery of the gas cloud G2 on a near-radial orbit about Sgr A* has prompted much speculation on its origin. In this Letter, we propose that G2 formed out of the debris stream produced by the removal of mass from the outer envelope of a nearby giant star. We perform hydrodynamical simulations of the returning tidal debris stream with cooling, and find that the stream condenses into clumps that fall periodically onto Sgr A*. We propose that one of these clumps is the observed G2 cloud, with the rest of the stream being detectable at lower Br-$γ$ emissivity along a trajectory that would trace from G2 to the star that was partially disrupted. By simultaneously fitting the orbits of S2, G2, and $\sim$ 2,000 candidate stars, and by fixing the orbital plane of each candidate star to G2 (as is expected for a tidal disruption), we find that several stars have orbits that are compatible with the notion that one of them was tidally disrupted to produce G2. If one of these stars were indeed disrupted, it last encountered Sgr A* hundreds of years ago, and has likely encountered Sgr A* repeatedly. However, while these stars are compatible with the giant disruption scenario given their measured positions and proper motions, their radial velocities are currently unknown. If one of these stars' radial velocity is measured to be compatible with a disruptive orbit, it would strongly suggest its disruption produced G2.

preprint2014arXiv

PS1-10jh: The Disruption of a Main-Sequence Star of Near-Solar Composition

When a star comes within a critical distance to a supermassive black hole (SMBH), immense tidal forces disrupt the star, resulting in a stream of debris that falls back onto the SMBH and powers a luminous flare. In this paper, we perform hydrodynamical simulations of the disruption of a main-sequence star by a SMBH to characterize the evolution of the debris stream after a tidal disruption. We demonstrate that this debris stream is confined by self-gravity in the two directions perpendicular to the original direction of the star's travel, and as a consequence has a negligible surface area and makes almost no contribution to either the continuum or line emission. We therefore propose that any observed emission lines are not the result of photoionization in this unbound debris, but are produced in the region above and below the forming elliptical accretion disk, analogous to the broad-line region (BLR) in steadily-accreting active galactic nuclei. As each line within a BLR is observationally linked to a particular location in the accretion disk, we suggest that the absence of a line indicates that the accretion disk does not yet extend to the distance required to produce that line. This model can be used to understand the spectral properties of the tidal disruption event (TDE) PS1-10jh, for which He II lines are observed, but the Balmer series and He I are not. Using a maximum likelihood analysis, we show that the disruption of a main-sequence star of near-solar composition can reproduce this event.

preprint2014arXiv

Turbovelocity Stars: Kicks Resulting From the Tidal Disruption of Solitary Stars

The centers of most known galaxies host supermassive black holes (SMBHs). In orbit around these black holes are a centrally-concentrated distribution of stars, both in single and in binary systems. Occasionally, these stars are perturbed onto orbits that bring them close to the SMBH. If the star is in a binary system, the three-body interaction with the SMBH can lead to large changes in orbital energy, depositing one of the two stars on a tightly-bound orbit, and its companion into a hyperbolic orbit that may escape the galaxy. In this Letter, we show that the disruption of solitary stars can also lead to large positive increases in orbital energy. The kick velocity depends on the amount of mass the star loses at pericenter, but not on the ratio of black hole to stellar mass, and are at most the star's own escape velocity. We find that these kicks are usually too small to result in the ejection of stars from the Milky Way, but can eject the stars from the black hole's sphere of influence, reducing their probability of being disrupted again. We estimate that ~10^5 stars, ~1% of all stars within 10 pc of the galactic center, are likely to have had mass removed by the central black hole through tidal interaction, and speculate that these "turbovelocity" stars will at first be redder, but eventually bluer, and always brighter than their unharrassed peers.

preprint2013arXiv

Circumstellar absorption in double detonation Type Ia supernovae

Upon formation, degenerate He core white dwarfs are surrounded by a radiative H-rich layer primarily supported by ideal gas pressure. In this Letter, we examine the effect of this H-rich layer on mass transfer in He+C/O double white dwarf binaries that will eventually merge and possibly yield a Type Ia supernova (SN Ia) in the double detonation scenario. Because its thermal profile and equation of state differ from the underlying He core, the H-rich layer is transferred stably onto the C/O white dwarf prior to the He core's tidal disruption. We find that this material is ejected from the binary system and sweeps up the surrounding interstellar medium hundreds to thousands of years before the SN Ia. The close match between the resulting circumstellar medium profiles and values inferred from recent observations of circumstellar absorption in SNe Ia gives further credence to the resurgent double detonation scenario.

preprint2013arXiv

Conditions For Successful Helium Detonations In Astrophysical Environments

Several models for type Ia-like supernovae events rely on the production of a self-sustained detonation powered by nuclear reactions. In the absence of hydrogen, the fuel that powers these detonations typically consists of either pure helium (He) or a mixture of carbon and oxygen (C/O). Studies that systematically determine the conditions required to initiate detonations in C/O material exist, but until now no analogous investigation of He matter has been conducted. We perform one-dimensional reactive hydrodynamical simulations at a variety of initial density and temperature combinations and find critical length scales for the initiation of He detonations that range between 1 - $10^{10}$ cm. A simple estimate of the length scales over which the total consumption of fuel will occur for steady-state detonations is provided by the Chapman-Jouguet (CJ) formalism. Our initiation lengths are consistently smaller than the corresponding CJ length scales by a factor of $\Sim 100$, providing opportunities for thermonuclear explosions in a wider range of low-mass white dwarfs (WDs) than previously thought possible. We find that virialized WDs with as little mass as 0.24 $M_\odot$ can be detonated, and that even less massive WDs can be detonated if a sizable fraction of their mass is raised to a higher adiabat. That the initiation length is exceeded by the CJ length implies that certain systems may not reach nuclear statistical equilibrium within the time it takes a detonation to traverse the object. In support of this hypothesis, we demonstrate that incomplete burning will occur in the majority of He WD detonations and that $^{40}$Ca, $^{44}$Ti, or $^{48}$Cr, rather than $^{56}$Ni, is the predominant burning product for many of these events. We anticipate that ...

preprint2013arXiv

Spoon-Feeding Giant Stars to Supermassive Black Holes: Episodic Mass Transfer from Evolving Stars and Their Contribution to the Quiescent Activity of Galactic Nuclei

Stars may be tidally disrupted if, in a single orbit, they are scattered too close to a supermassive black hole (SMBH). Tidal disruption events are thought to power luminous but short-lived accretion episodes that can light up otherwise quiescent SMBHs in transient flares. Here we explore a more gradual process of tidal stripping where stars approach the tidal disruption radius by stellar evolution while in an eccentric orbit. After the onset of mass transfer, these stars episodically transfer mass to the SMBH every pericenter passage giving rise to low-level flares that repeat on the orbital timescale. Giant stars, in particular, will exhibit a runaway response to mass loss and "spoon-feed" material to the black hole for tens to hundreds of orbital periods. In contrast to full tidal disruption events, the duty cycle of this feeding mode is of order unity for black holes with mass greater than approximately 10 million solar masses. This mode of quasi-steady SMBH feeding is competitive with indirect SMBH feeding through stellar winds, and spoon-fed giant stars may play a significant role in determining the quiescent luminosity of local SMBHs.

preprint2012arXiv

How the merger of two white dwarfs depends on their mass ratio: orbital stability and detonations at contact

Despite their unique astrophysical relevance, the outcome of white dwarf binary mergers has so far only been studied for a very restricted number of systems. Here we present the results of a survey with more than two hundred simulations systematically scanning the white dwarf binary parameter space. We consider white dwarf masses ranging from 0.2 to 1.2 $M_\odot$ and account for their different chemical compositions. We find excellent agreement with the orbital evolution predicted by mass transfer stability analysis. Much of our effort in this paper is dedicated to determining which binary systems are prone to a thermonuclear explosion just prior to merger or at surface contact. We find that a large fraction of He-accreting binary systems explode: all dynamically unstable systems with accretor masses below 1.1 $M_\odot$ and donor masses above $\sim$ 0.4 $M_\odot$ are found to trigger a helium detonation at surface contact. A substantial fraction of these systems could explode at earlier times via detonations induced by instabilities in the accretion stream, as we have demonstrated in our previous work. We do not find definitive evidence for an explosion prior to merger or at surface contact in any of the studied double carbon-oxygen systems. Although we cannot exclude their occurrence if some helium is present, the available parameter space for a successful detonation in a white dwarf binary of pure carbon-oxygen composition is small. We demonstrate that a wide variety of dynamically unstable systems are viable type Ia candidates. The next decade thus holds enormous promise for the study of these events, in particular with the advent of wide-field synoptic surveys allowing a detailed characterization of their explosive properties.

preprint2012arXiv

On the Survivability and Metamorphism of Tidally Disrupted Giant Planets: the Role of Dense Cores

A large population of planetary candidates in short-period orbits have been found through transit searches. Radial velocity surveys have also revealed several Jupiter-mass planets with highly eccentric orbits. Measurements of the Rossiter-McLaughlin effect indicate some misaligned planetary systems. This diversity could be induced by post-formation dynamical processes such as planet-planet scattering, the Kozai effect, or secular chaos which brings planets to the vicinity of their host stars. In this work, we propose a novel mechanism to form close-in super-Earths and Neptune-like planets through the tidal disruption of giant planets as a consequence of these dynamical processes. We model the core-envelope structure of giant planets with composite polytropes. Using three-dimensional hydrodynamical simulations of close encounters between planets and their host stars, we find that the presence of a core with a mass more than ten Earth masses can significantly increase the fraction of envelope which remains bound to it. After the encounter, planets with cores are more likely to be retained by their host stars in contrast with previous studies which suggested that coreless planets are often ejected. As a substantial fraction of their gaseous envelopes is preferentially lost while the dense incompressible cores retain most of their original mass, the resulting metallicity of the surviving planets is increased. Our results suggest that some gas giant planets can be effectively transformed into either super-Earths or Neptune-like planets after multiple close stellar passages. Finally, we analyze the orbits and structure of known planets and Kepler candidates and find that our model is capable producing some of the shortest-period objects.

preprint2012arXiv

The dynamics, appearance and demographics of relativistic jets triggered by tidal disruption of stars in quiescent supermassive black holes

We examine the consequences of a model in which relativistic jets can be triggered in quiescent massive black holes when a geometrically thick and hot accretion disk forms as a result of the tidal disruption of a star. To estimate the power, thrust and lifetime of the jet, we use the mass accretion history onto the black hole as calculated by detailed hydrodynamic simulations of the tidal disruption of stars. We go on to determine the states of the interstellar medium in various types of quiescent galactic nuclei, and describe how this external matter can affect jets propagating through it. We use this information, together with a two-dimensional hydrodynamic model of the structure of the relativistic flow, to study the dynamics of the jet, the propagation of which is regulated by the density stratification of the environment and by its injection history. The breaking of symmetry involved in transitioning from one to two dimensions is crucial and leads to qualitatively new phenomena. Many of the observed properties of the Swift 1644+57/GRB 110328A event can be understood as resulting from accretion onto and jets driven by a $10^6 M_\odot$ central mass black hole following the disruption of sun-like star. With the inclusion of a stochastic contribution to the luminosity due to variations in the feeding rate driven by instabilities near the tidal radius, we find that our model can explain the X-ray light curve without invoking a rarely-occurring deep encounter. In conjunction with the number density of black holes in the local universe, we hypothesize that the conditions required to produce the Swift event are not anomalous, but are in fact representative of the jet-driven flare population arising from tidal disruptions. [abridged]

preprint2012arXiv

The Tidal Disruption of Giant Stars and Their Contribution to the Flaring Supermassive Black Hole Population

Sun-like stars are thought to be regularly disrupted by supermassive black holes (SMBHs) within galactic nuclei. Yet, as stars evolve off the main sequence their vulnerability to tidal disruption increases drastically as they develop a bifurcated structure consisting of a dense core and a tenuous envelope. Here we present the first hydrodynamic simulations of the tidal disruption of giant stars and show that the core has a substantial influence on the star's ability to survive the encounter. Stars with more massive cores retain large fractions of their envelope mass, even in deep encounters. Accretion flares resulting from the disruption of giant stars should last for tens to hundreds of years. Their characteristic signature in transient searches would not be the $t^{-5/3}$ decay typically associated with tidal disruption events, but a correlated rise over many orders of magnitude in brightness on months to years timescales. We calculate the relative disruption rates of stars of varying evolutionary stages in typical galactic centers, then use our results to produce Monte Carlo realizations of the expected flaring event populations. We find that the demographics of tidal disruption flares are strongly dependent on both stellar and black hole mass, especially near the limiting SMBH mass scale of $\sim 10^8 M_\odot$. At this black hole mass, we predict a sharp transition in the SMBH flaring diet beyond which all observable disruptions arise from evolved stars, accompanied by a dramatic cutoff in the overall tidal disruption flaring rate. Black holes less massive than this limiting mass scale will show observable flares from both main sequence and evolved stars, with giants contributing up to 10% of the event rate. The relative fractions of stars disrupted at different evolutionary states can constrain the properties and distributions of stars in galactic nuclei other than our own.

preprint2011arXiv

A relativistic jetted outburst from a massive black hole fed by a tidally disrupted star

While gas accretion onto some massive black holes (MBHs) at the centers of galaxies actively powers luminous emission, the vast majority of MBHs are considered dormant. Occasionally, a star passing too near a MBH is torn apart by gravitational forces, leading to a bright panchromatic tidal disruption flare (TDF). While the high-energy transient Swift J164449.3+573451 ("Sw 1644+57") initially displayed none of the theoretically anticipated (nor previously observed) TDF characteristics, we show that the observations (Levan et al. 2011) suggest a sudden accretion event onto a central MBH of mass ~10^6-10^7 solar masses. We find evidence for a mildly relativistic outflow, jet collimation, and a spectrum characterized by synchrotron and inverse Compton processes; this leads to a natural analogy of Sw 1644+57 with a smaller-scale blazar. The phenomenologically novel Sw 1644+57 thus connects the study of TDFs and active galaxies, opening a new vista on disk-jet interactions in BHs and magnetic field generation and transport in accretion systems.

preprint2011arXiv

Prelude to a double degenerate merger: the onset of mass transfer and its impact on gravitational waves and surface detonations

We present the results of a systematic numerical study of the onset of mass transfer in double degenerate binary systems and its impact on the subsequent evolution. All investigated systems belong to the regime of direct impact, unstable mass transfer. In all of the investigated cases, even those considered unstable by conventional stability analysis, we find a long-lived mass transfer phase continuing for as many as several dozen orbital periods. This settles a recent debate sparked by a discrepancy between earlier SPH calculations that showed disruptions after a few orbital periods and newer grid-based studies in which mass transfer continued for tens of orbits. As we show that these binaries can survive at small separation for hundreds of orbital periods, their associated gravitational wave signal should be included when calculating the gravitational wave foreground (although expected to below LISA's sensitivity at these high frequencies). We also show that the inclusion of the entropy increase associated with shock-heating of the accreted material reduces the number of orbits a binary survives given the same initial conditions, although the effect is not as pronounced when using the appropriate initial conditions. The use of accurate initial conditions and a correct treatment of shock heating allows for a reliable time evolution of the temperature, density, and angular momentum, which are important when considering thermonuclear events that may occur during the mass transfer phase and/or after merger. Our treatment allows us to accurately identify when surface detonations may occur in the lead-up to the merger, as well as the properties of final merger products.

preprint2010arXiv

Consequences of the Ejection and Disruption of Giant Planets

The discovery of Jupiter-mass planets in close orbits about their parent stars has challenged models of planet formation. Recent observations have shown that a number of these planets have highly inclined, sometimes retrograde orbits about their parent stars, prompting much speculation as to their origin. It is known that migration alone cannot account for the observed population of these misaligned hot Jupiters, which suggests that dynamical processes after the gas disc dissipates play a substantial role in yielding the observed inclination and eccentricity distributions. One particularly promising candidate is planet-planet scattering, which is not very well understood in the non-linear regime of tides. Through three-dimensional hydrodynamical simulations of multi-orbit encounters, we show that planets that are scattered into an orbit about their parent stars with closest approach distance being less than approximately three times the tidal radius are either destroyed or completely ejected from the system. We find that as few as 5 and as many as 18 of the currently known hot Jupiters have a maximum initial apastron for scattering that lies well within the ice line, implying that these planets must have migrated either before or after the scattering event that brought them to their current positions. If stellar tides are unimportant $(Q_\ast \gtrsim 10^7)$, disk migration is required to explain the existence of the hot Jupiters present in these systems. Additionally, we find that the disruption and/or ejection of Jupiter-mass planets deposits a Sun's worth of angular momentum onto the host star. For systems in which planet-planet scattering is common, we predict that planetary hosts have up to a 35% chance of possessing an obliquity relative to the invariable plane of greater than 90 degrees.

preprint2009arXiv

Surface Detonations in Double Degenerate Binary Systems Triggered by Accretion Stream Instabilities

We present three-dimensional simulations on a new mechanism for the detonation of a sub-Chandrasekhar CO white dwarf in a dynamically unstable system where the secondary is either a pure He white dwarf or a He/CO hybrid. For dynamically unstable systems where the accretion stream directly impacts the surface of the primary, the final tens of orbits can have mass accretion rates that range from $10^{-5}$ to $10^{-3} M_{\odot}$ s$^{-1}$, leading to the rapid accumulation of helium on the surface of the primary. After $\sim 10^{-2}$ $M_{\odot}$ of helium has been accreted, the ram pressure of the hot helium torus can deflect the accretion stream such that the stream no longer directly impacts the surface. The velocity difference between the stream and the torus produces shearing which seeds large-scale Kelvin-Helmholtz instabilities along the interface between the two regions. These instabilities eventually grow into dense knots of material that periodically strike the surface of the primary, adiabatically compressing the underlying helium torus. If the temperature of the compressed material is raised above a critical temperature, the timescale for triple-$α$ reactions becomes comparable to the dynamical timescale, leading to the detonation of the primary's helium envelope. This detonation drives shockwaves into the primary which tend to concentrate at one or more focal points within the primary's CO core. If a relatively small amount of mass is raised above a critical temperature and density at these focal points, the CO core may itself be detonated.

preprint2009arXiv

Three-Dimensional Simulations of Tidally Disrupted Solar-Type Stars and the Observational Signatures of Shock Breakout

We describe a three-dimensional simulation of a $1 M_{\odot}$ solar-type star approaching a $10^{6} M_{\odot}$ black hole on a parabolic orbit with a pericenter distance well within the tidal radius. While falling towards the black hole, the star is not only stretched along the orbital direction but even more severely compressed at right angles to the orbit. The overbearing degree of compression achieved shortly after pericenter leads to the production of strong shocks which largely homogenize the temperature profile of the star, resulting in surface temperatures comparable to the initial temperature of the star's core. This phenomenon, which precedes the fallback accretion phase, gives rise to a unique double-peaked X-ray signature that, if detected, may be one of the few observable diagnostics of how stars behave under the influence of strong gravitational fields. If $\sim 10^{6} M_{\odot}$ black holes were prevalent in small or even dwarf galaxies, the nearest of such flares may be detectable by EXIST from no further away than the Virgo Cluster.