sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Guillochon, James

Normalized to: Guillochon, J.

72 article(s) in total. 779 co-authors, from 1 to 32 common article(s). Median position in authors list is 3,0.

[1]  [pdf] - 1981225
The Tidal Disruption of Sun-like Stars by Massive Black Holes
Comments: 9 pages, 5 figures. Accepted for publication in ApJ Letters
Submitted: 2019-07-10, last modified: 2019-08-06
We present the first simulations of the tidal disruption of stars with realistic structures and compositions by massive black holes (BHs). We build stars in the stellar evolution code MESA and simulate their disruption in the 3D adaptive-mesh hydrodynamics code FLASH, using an extended Helmholtz equation of state and tracking 49 elements. We study the disruption of a 1$M_\odot$ star and 3$M_\odot$ star at zero-age main sequence (ZAMS), middle-age, and terminal-age main sequence (TAMS). The maximum BH mass for tidal disruption increases by a factor of $\sim$2 from stellar radius changes due to MS evolution; this is equivalent to varying BH spin from 0 to 0.75. The shape of the mass fallback rate curves is different from the results for polytropes of Guillochon & Ramirez-Ruiz (2013). The peak timescale $t_{\rm peak}$ increases with stellar age, while the peak fallback rate $\dot M_{\rm peak}$ decreases with age, and these effects diminish with increasing impact parameter $\beta$. For a $\beta=1$ disruption of a 1$M_\odot$ star by a $10^6 M_\odot$ BH, from ZAMS to TAMS, $t_{\rm peak}$ increases from 30 to 54 days, while $\dot M_{\rm peak}$ decreases from 0.66 to 0.14 $M_\odot$/yr. Compositional anomalies in nitrogen, helium, and carbon can occur before the peak timescale for disruptions of MS stars, which is in contrast to predictions from the "frozen-in" model. More massive stars can show stronger anomalies at earlier times, meaning that compositional constraints can be key in determining the mass of the disrupted star. The abundance anomalies predicted by these simulations provide a natural explanation for the spectral features and varying line strengths observed in tidal disruption events.
[2]  [pdf] - 1927448
The tidal disruption event AT2017eqx: spectroscopic evolution from hydrogen rich to poor suggests an atmosphere and outflow
Comments: Updated to match published version
Submitted: 2019-04-23, last modified: 2019-08-02
We present and analyse a new tidal disruption event (TDE), AT2017eqx at redshift z=0.1089, discovered by Pan-STARRS and ATLAS. The position of the transient is consistent with the nucleus of its host galaxy; it peaks at a luminosity of $L \approx 10^{44}$ erg s$^{-1}$; and the spectrum shows a persistent blackbody temperature $T \gtrsim 20,000$ K with broad H I and He II emission. The lines are initially centered at zero velocity, but by 100 days the H I lines disappear while the He II develops a blueshift of $\gtrsim 5,000$ km s$^{-1}$. Both the early- and late-time morphologies have been seen in other TDEs, but the complete transition between them is unprecedented. The evolution can be explained by combining an extended atmosphere, undergoing slow contraction, with a wind in the polar direction becoming visible at late times. Our observations confirm that a lack of hydrogen a TDE spectrum does not indicate a stripped star, while the proposed model implies that much of the diversity in TDEs may be due to the observer viewing angle. Modelling the light curve suggests AT2017eqx resulted from the complete disruption of a solar-mass star by a black hole of $\sim 10^{6.3} M_\odot$. The host is another quiescent, Balmer-strong galaxy, though fainter and less centrally concentrated than most TDE hosts. Radio limits rule out a relativistic jet, while X-ray limits at 500 days are among the deepest for a TDE at this phase.
[3]  [pdf] - 1924999
Models and Simulations for the Photometric LSST Astronomical Time Series Classification Challenge (PLAsTiCC)
Submitted: 2019-03-27, last modified: 2019-07-10
We describe the simulated data sample for the "Photometric LSST Astronomical Time Series Classification Challenge" (PLAsTiCC), a publicly available challenge to classify transient and variable events that will be observed by the Large Synoptic Survey Telescope (LSST), a new facility expected to start in the early 2020s. The challenge was hosted by Kaggle, ran from 2018 September 28 to 2018 December 17, and included 1,094 teams competing for prizes. Here we provide details of the 18 transient and variable source models, which were not revealed until after the challenge, and release the model libraries at We describe the LSST Operations Simulator used to predict realistic observing conditions, and we describe the publicly available SNANA simulation code used to transform the models into observed fluxes and uncertainties in the LSST passbands (ugrizy). Although PLAsTiCC has finished, the publicly available models and simulation tools are being used within the astronomy community to further improve classification, and to study contamination in photometrically identified samples of type Ia supernova used to measure properties of dark energy. Our simulation framework will continue serving as a platform to improve the PLAsTiCC models, and to develop new models.
[4]  [pdf] - 1846780
Constraining the Stellar Mass Function from the Deficiency of Tidal Disruption Flares in the Nuclei of Massive Galaxies
Comments: Accepted for publication in MNRAS
Submitted: 2018-06-29, last modified: 2019-02-28
The rate of tidal disruption flares (TDFs) per mass of the disrupting black hole encodes information on the present-day mass function (PDMF) of stars in the clusters surrounding super massive black holes. We explore how the shape of the TDF rate with black hole mass can constrain the PDMF, with only weak dependence on black hole spin. We show that existing data can marginally constrain the minimum and maximum masses of stars in the cluster, and the high-mass end of the PDMF slope, as well as the overall TDF rate. With of order 100 TDFs expected to be identified with the Zwicky Transient Facility, the overall rate can be highly constrained, but still with only marginal constraints on the PDMF. However, if less than of order 10% of the TDFs expected to be found by LSST over a decade (of order 1000 TDFs) are identified, then precise and accurate estimates can be made for the minimum stellar mass (within a factor of two) and the average slope of the high-mass PDMF (to within ~10%) in nuclear star clusters. This technique could be adapted in the future to probe, in addition to the PDMF, the local black hole mass function and possibly the massive black hole binary population.
[5]  [pdf] - 1828935
Catching Element Formation In The Act
Fryer, Chris L.; Timmes, Frank; Hungerford, Aimee L.; Couture, Aaron; Adams, Fred; Aoki, Wako; Arcones, Almudena; Arnett, David; Auchettl, Katie; Avila, Melina; Badenes, Carles; Baron, Eddie; Bauswein, Andreas; Beacom, John; Blackmon, Jeff; Blondin, Stephane; Bloser, Peter; Boggs, Steve; Boss, Alan; Brandt, Terri; Bravo, Eduardo; Brown, Ed; Brown, Peter; Budtz-Jorgensen, Steve Bruenn. Carl; Burns, Eric; Calder, Alan; Caputo, Regina; Champagne, Art; Chevalier, Roger; Chieffi, Alessandro; Chipps, Kelly; Cinabro, David; Clarkson, Ondrea; Clayton, Don; Coc, Alain; Connolly, Devin; Conroy, Charlie; Cote, Benoit; Couch, Sean; Dauphas, Nicolas; deBoer, Richard James; Deibel, Catherine; Denisenkov, Pavel; Desch, Steve; Dessart, Luc; Diehl, Roland; Doherty, Carolyn; Dominguez, Inma; Dong, Subo; Dwarkadas, Vikram; Fan, Doreen; Fields, Brian; Fields, Carl; Filippenko, Alex; Fisher, Robert; Foucart, Francois; Fransson, Claes; Frohlich, Carla; Fuller, George; Gibson, Brad; Giryanskaya, Viktoriya; Gorres, Joachim; Goriely, Stephane; Grebenev, Sergei; Grefenstette, Brian; Grohs, Evan; Guillochon, James; Harpole, Alice; Harris, Chelsea; Harris, J. Austin; Harrison, Fiona; Hartmann, Dieter; Hashimoto, Masa-aki; Heger, Alexander; Hernanz, Margarita; Herwig, Falk; Hirschi, Raphael; Hix, Raphael William; Hoflich, Peter; Hoffman, Robert; Holcomb, Cole; Hsiao, Eric; Iliadis, Christian; Janiuk, Agnieszka; Janka, Thomas; Jerkstrand, Anders; Johns, Lucas; Jones, Samuel; Jose, Jordi; Kajino, Toshitaka; Karakas, Amanda; Karpov, Platon; Kasen, Dan; Kierans, Carolyn; Kippen, Marc; Korobkin, Oleg; Kobayashi, Chiaki; Kozma, Cecilia; Krot, Saha; Kumar, Pawan; Kuvvetli, Irfan; Laird, Alison; Laming, Martin; Larsson, Josefin; Lattanzio, John; Lattimer, James; Leising, Mark; Lennarz, Annika; Lentz, Eric; Limongi, Marco; Lippuner, Jonas; Livne, Eli; Lloyd-Ronning, Nicole; Longland, Richard; Lopez, Laura A.; Lugaro, Maria; Lutovinov, Alexander; Madsen, Kristin; Malone, Chris; Matteucci, Francesca; McEnery, Julie; Meisel, Zach; Messer, Bronson; Metzger, Brian; Meyer, Bradley; Meynet, Georges; Mezzacappa, Anthony; Miller, Jonah; Miller, Richard; Milne, Peter; Misch, Wendell; Mitchell, Lee; Mosta, Philipp; Motizuki, Yuko; Muller, Bernhard; Mumpower, Matthew; Murphy, Jeremiah; Nagataki, Shigehiro; Nakar, Ehud; Nomoto, Ken'ichi; Nugent, Peter; Nunes, Filomena; O'Shea, Brian; Oberlack, Uwe; Pain, Steven; Parker, Lucas; Perego, Albino; Pignatari, Marco; Pinedo, Gabriel Martinez; Plewa, Tomasz; Poznanski, Dovi; Priedhorsky, William; Pritychenko, Boris; Radice, David; Ramirez-Ruiz, Enrico; Rauscher, Thomas; Reddy, Sanjay; Rehm, Ernst; Reifarth, Rene; Richman, Debra; Ricker, Paul; Rijal, Nabin; Roberts, Luke; Ropke, Friedrich; Rosswog, Stephan; Ruiter, Ashley J.; Ruiz, Chris; Savin, Daniel Wolf; Schatz, Hendrik; Schneider, Dieter; Schwab, Josiah; Seitenzahl, Ivo; Shen, Ken; Siegert, Thomas; Sim, Stuart; Smith, David; Smith, Karl; Smith, Michael; Sollerman, Jesper; Sprouse, Trevor; Spyrou, Artemis; Starrfield, Sumner; Steiner, Andrew; Strong, Andrew W.; Sukhbold, Tuguldur; Suntzeff, Nick; Surman, Rebecca; Tanimori, Toru; The, Lih-Sin; Thielemann, Friedrich-Karl; Tolstov, Alexey; Tominaga, Nozomu; Tomsick, John; Townsley, Dean; Tsintari, Pelagia; Tsygankov, Sergey; Vartanyan, David; Venters, Tonia; Vestrand, Tom; Vink, Jacco; Waldman, Roni; Wang, Lifang; Wang, Xilu; Warren, MacKenzie; West, Christopher; Wheeler, J. Craig; Wiescher, Michael; Winkler, Christoph; Winter, Lisa; Wolf, Bill; Woolf, Richard; Woosley, Stan; Wu, Jin; Wrede, Chris; Yamada, Shoichi; Young, Patrick; Zegers, Remco; Zingale, Michael; Zwart, Simon Portegies
Comments: 14 pages including 3 figures
Submitted: 2019-02-07
Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.
[6]  [pdf] - 1838119
Weighing Black Holes using Tidal Disruption Events
Comments: accepted to ApJ; added table 3, table 4, figure 5
Submitted: 2018-01-24, last modified: 2019-02-05
While once rare, observations of stars being tidally disrupted by supermassive black holes are quickly becoming commonplace. To continue to learn from these events it is necessary to robustly and systematically compare our growing number of observations with theory. We present a tidal disruption module for the Modular Open Source Fitter for Transients (MOSFiT) and the results from fitting 14 tidal disruption events (TDEs). Our model uses FLASH simulations of TDEs to generate bolometric luminosities and passes these luminosities through viscosity and reprocessing transformation functions to create multi-wavelength light curves. It then uses an MCMC fitting routine to compare these theoretical light curves with observations. We find that none of the events show evidence for viscous delays exceeding a few days, supporting the theory that our current observing strategies in the optical/UV are missing a significant number of viscously delayed flares. We find that the events have black hole masses of $10^6 - 10^8 M_{\odot}$, and that the masses we predict are as reliable as those based on bulk galaxy properties. We also find that there is a preference for stars with mass $< 1 M_{\odot}$, as expected when low-mass stars greatly outnumber high-mass stars.
[7]  [pdf] - 1830648
astroquery: An Astronomical Web-Querying Package in Python
Comments: Accepted to AJ. The documentation is at and the repository is at
Submitted: 2019-01-14
astroquery is a collection of tools for requesting data from databases hosted on remote servers with interfaces exposed on the internet, including those with web pages but without formal application program interfaces (APIs). These tools are built on the Python requests package, which is used to make HTTP requests, and astropy, which provides most of the data parsing functionality. astroquery modules generally attempt to replicate the web page interface provided by a given service as closely as possible, making the transition from browser-based to command-line interaction easy. astroquery has received significant contributions from throughout the astronomical community, including several significant contributions from telescope archives. astroquery enables the creation of fully reproducible workflows from data acquisition through publication. This paper describes the philosophy, basic structure, and development model of the astroquery package. The complete documentation for astroquery can be found at
[8]  [pdf] - 1823078
Discovery and Follow-up of the Unusual Nuclear Transient OGLE17aaj
Comments: 5 pages, 3 figures, accepted for publication in Astronomy & Astrophysics as Letter to Editor
Submitted: 2019-01-11
We report on the discovery and follow-up of a peculiar transient, OGLE17aaj, which occurred in the nucleus of a weakly active galaxy. We investigate whether it can be interpreted as a new candidate for a tidal disruption event (TDE). We present the OGLE-IV light curve that covers the slow 60-day-long rise to maximum along with photometric, spectroscopic, and X-ray follow-up during the first year. OGLE17aaj is a nuclear transient exhibiting some properties similar to previously found TDEs, including a long rise time, lack of colour-temperature evolution, and high black-body temperature. On the other hand, its narrow emission lines and slow post-peak evolution are different from previously observed TDEs. Its spectrum and light-curve evolution is similar to F01004-2237 and AT 2017bgt. Signatures of historical low-level nuclear variability suggest that OGLE17aaj may instead be related to a new type of accretion event in active super-massive black holes.
[9]  [pdf] - 1833029
Systematic investigation of the fallback accretion powered model for hydrogen-poor superluminous supernovae
Comments: 12 pages, 8 figures, 3 tables, accepted by The Astrophysical Journal
Submitted: 2018-05-31, last modified: 2018-09-28
The energy liberated by fallback accretion has been suggested as a possible engine to power hydrogen-poor superluminous supernovae. We systematically investigate this model using the Bayesian light-curve fitting code MOSFiT (Modular Open Source Fitter for Transients), fitting the light curves of 37 hydrogen-poor superluminous supernovae assuming a fallback accretion central engine. We find that this model can yield good fits to their light curves, with a fit quality that rivals the popular magnetar engine models. Examining our derived parameters for the fallback model, we find the total energy requirements from the accretion disk are estimated to be 0.002 - 0.7 Msun c^2. If we adopt a typical conversion efficiency ~ 1e-3, the required mass to accrete is thus 2 - 700 Msun. Many superluminous supernovae, therefore, require an unrealistic accretion mass, and so only a fraction of these events could be powered by fallback accretion unless the true efficiency is much greater than our fiducial value. The superluminous supernovae that require the smallest amounts of fallback mass still remain to be the fallback accretion powered supernova candidates, but they are difficult to be distinguished solely by their light curve properties.
[10]  [pdf] - 1771675
Shock Breakouts from Tidal Disruption Events
Submitted: 2018-08-30, last modified: 2018-09-06
Tidal disruption events of stars by supermassive black holes have so far been discovered months to years after the fact. In this paper we explore the short, faint and hard burst of radiation is emitted at maximum compression, as a result of shock breakout. The detection of this burst can be used to capture tidal disruption events in real time. We verify that shock breakout from main sequence stars produces radiation in the X-ray range, but find that it is difficult to detect using all sky X-ray surveying telescopes. In the case of shock breakout from red giants, most of the radiation is emitted in the UV and visible range, which is significantly easier to detect. A similar burst of UV/optical radiation will also be emitted by stars puffed by tidal heating from a previous passage close to the central black hole. This radiation can be detected by surveys like ZTF and LSST. We calculate detection rates for different types of galactic nuclei. For the case of a very full loss cone we predict a detection rate of once per month with LSST, whereas for the case of a very empty loss cone we predict a rate of once per year with LSST. Evidence from a recent tidal disruption event, ASASSN-14li, seems to favour a very full loss cone, in which case LSST is expected to detect one such event every month.
[11]  [pdf] - 1743955
Serendipitous Discovery of a 14-year-old Supernova at 16 Mpc
Comments: 4 pages, 1 figure. Accepted to RNAAS
Submitted: 2018-09-01
In this Research Note we present a serendipitous discovery of the transient CGS2004A (AT2004iu) in an image collected by the Carnegie-Irvine Galaxy Survey of the Scd: galaxy NGC 1892, which we determine is most likely to have been a type IIP supernova.
[12]  [pdf] - 1755904
Three Hypervelocity White Dwarfs in Gaia DR2: Evidence for Dynamically Driven Double-Degenerate Double-Detonation Type Ia Supernovae
Comments: Accepted for publication in ApJ. Minor corrections for clarity. D6 spectra are available as ancillary data files
Submitted: 2018-04-30, last modified: 2018-07-18
Double detonations in double white dwarf (WD) binaries undergoing unstable mass transfer have emerged in recent years as one of the most promising Type Ia supernova (SN Ia) progenitor scenarios. One potential outcome of this "dynamically driven double-degenerate double-detonation" (D^6) scenario is that the companion WD survives the explosion and is flung away with a velocity equal to its > 1000 km/s pre-SN orbital velocity. We perform a search for these hypervelocity runaway WDs using Gaia's second data release. In this paper, we discuss seven candidates followed up with ground-based instruments. Three sources are likely to be some of the fastest known stars in the Milky Way, with total Galactocentric velocities between 1000 and 3000 km/s, and are consistent with having previously been companion WDs in pre-SN Ia systems. However, although the radial velocity of one of the stars is > 1000 km/s, the radial velocities of the other two stars are puzzlingly consistent with 0. The combined five-parameter astrometric solutions from Gaia and radial velocities from follow-up spectra yield tentative 6D confirmation of the D^6 scenario. The past position of one of these stars places it within a faint, old SN remnant, further strengthening the interpretation of these candidates as hypervelocity runaways from binary systems that underwent SNe Ia.
[13]  [pdf] - 1715387
Multi-Messenger Astrophysics: Harnessing the Data Revolution
Comments: workshop summary, 20 pages, by the participants of the Cyberinfrastructure for Multi-Messenger Astrophysics NSF-funded Workshop Held on 23-24 May 2018 at the University of Maryland
Submitted: 2018-07-12
The past year has witnessed discovery of the first identified counterparts to a gravitational wave transient (GW 170817A) and a very high-energy neutrino (IceCube-170922A). These source identifications, and ensuing detailed studies, have realized longstanding dreams of astronomers and physicists to routinely carry out observations of cosmic sources by other than electromagnetic means, and inaugurated the era of "multi-messenger" astronomy. While this new era promises extraordinary physical insights into the universe, it brings with it new challenges, including: highly heterogeneous, high-volume, high-velocity datasets; globe-spanning cross-disciplinary teams of researchers, regularly brought together into transient collaborations; an extraordinary breadth and depth of domain-specific knowledge and computing resources required to anticipate, model, and interpret observations; and the routine need for adaptive, distributed, rapid-response observing campaigns to fully exploit the scientific potential of each source. We argue, therefore, that the time is ripe for the community to conceive and propose an Institute for Multi-Messenger Astrophysics that would coordinate its resources in a sustained and strategic fashion to efficiently address these challenges, while simultaneously serving as a center for education and key supporting activities. In this fashion, we can prepare now to realize the bright future that we see, beyond, through these newly opened windows onto the universe.
[14]  [pdf] - 1709438
Revisiting hypervelocity stars after Gaia DR2
Comments: 9 pages, 6 figures, 1 table. Accepted to MNRAS. The Open Fast Stars Catalog is available at
Submitted: 2018-04-26, last modified: 2018-06-05
Hypervelocity stars are intriguing rare objects traveling at speeds large enough to be unbound from the Milky Way. Several mechanisms have been proposed for producing them, including the interaction of the Galaxy's super-massive black hole (SMBH) with a binary; rapid mass-loss from a companion to a star in a short-period binary; the tidal disruption of an infalling galaxy and finally ejection from the Large Magellanic Cloud. While previously discovered high-velocity early-type stars are thought to be the result of an interaction with the SMBH, the origin of high-velocity late type stars is ambiguous. The second data release of Gaia (DR2) enables a unique opportunity to resolve this ambiguity and determine whether any late-type candidates are truly unbound from the Milky Way. In this paper, we utilize the new proper motion and velocity information available from DR2 to re-evaluate a collection of historical data compiled on the newly-created Open Fast Stars Catalog. We find that almost all previously-known high-velocity late-type stars are most likely bound to the Milky Way. Only one late-type object (LAMOST J115209.12+120258.0) is unbound from the Galaxy. Performing integrations of orbital histories, we find that this object cannot have been ejected from the Galactic centre and thus may be either debris from the disruption of a satellite galaxy or a disc runaway.
[15]  [pdf] - 1728995
Spitzer Space Telescope Infrared Observations of the Binary Neutron Star Merger GW170817
Comments: 6 pages, 2 figures, submitted to ApJL
Submitted: 2018-05-21
We present Spitzer Space Telescope 3.6 and 4.5 micron observations of the binary neutron star merger GW170817 at 43, 74, and 264 days post-merger. Using the final observation as a template, we uncover a source at the position of GW170817 at 4.5 micron with a brightness of 22.9+/-0.3 AB mag at 43 days and 23.8+/-0.3 AB mag at 74 days (the uncertainty is dominated by systematics from the image subtraction); no obvious source is detected at 3.6 micron to a 3-sigma limit of >23.3 AB mag in both epochs. The measured brightness is dimmer by a factor of about 2-3 times compared to our previously published kilonova model, which is based on UV, optical, and near-IR data at <30 days. However, the observed fading rate and color (m_{3.6}-m_{4.5}> 0 AB mag) are consistent with our model. We suggest that the discrepancy is likely due to a transition to the nebular phase, or a reduced thermalization efficiency at such late time. Using the Spitzer data as a guide, we briefly discuss the prospects of observing future binary neutron star mergers with Spitzer (in LIGO/Virgo Observing Run 3) and the James Webb Space Telescope (in LIGO/Virgo Observing Run 4 and beyond).
[16]  [pdf] - 1673995
Open Astronomy Catalogs API
Comments: Published in RNAAS. A living description of the API is available at
Submitted: 2018-04-28
We announce the public release of the application program interface (API) for the Open Astronomy Catalogs (OACs), the OACAPI. The OACs serve near-complete collections of supernova, tidal disruption, kilonova, and fast stars data (including photometry, spectra, radio, and X-ray observations) via a user-friendly web interface that displays the data interactively and offers full data downloads. The OACAPI, by contrast, enables users to specifically download particular pieces of the OAC dataset via a flexible programmatic syntax, either via URL GET requests, or via a module within the astroquery Python package.
[17]  [pdf] - 1672460
Comparison of the X-ray emission from Tidal Disruption Events with those of Active Galactic Nuclei
Comments: 11 pages, 6 figures. Accepted to ApJ
Submitted: 2017-03-17, last modified: 2018-04-25
One of the main challenges of current tidal disruption events (TDEs) studies is that emission arising from AGN activity may potentially mimic the expected X-ray emission of a TDE. Here we compare the X-ray properties of TDEs and AGN to determine a set of characteristics which would allow us to discriminate between flares arising from these two objects. We find that at peak, TDEs are brighter than AGN found at similar redshifts. However, compared to preflare upperlimits, highly variable AGN can produce flares of a similar order of magnitude as those seen from X-ray TDEs. Nevertheless, TDEs decay significantly more monotonically, and their emission exhibits little variation in spectral hardness as a function of time. We also find that X-ray TDEs are less absorbed, and their emission is much softer than the emission detected from AGN found at similar redshifts. We derive the X-ray luminosity function (LF) for X-ray TDEs using the events from Auchettl et al. (2017). Interestingly, our X-ray LF matches closely the theoretically derived LF by Milosavljevi\'c et al. (2006). which assumes a higher TDE rate currently estimated from observations. Using our results and the results of Stone & Metzger (2016), we estimate a TDE rate of $(0.7-4.7)\times10^{-4}$ yr$^{-1}$ per galaxy, higher than current observational estimates. We find that TDEs can contribute significantly to the LF of AGN for $z\lesssim0.4$, while there is no evidence that TDEs influence the growth of $10^{6-7}M_{\odot}$ BHs. However, BHs $<10^{6}M_{\odot}$ can grow from TDEs arising from super-Eddington accretion without contributing significantly to the observed AGN LF at $z=0$.
[18]  [pdf] - 1693944
Observations of SN 2015F suggest a correlation between the intrinsic luminosity of Type Ia supernovae and the shape of their light curves >900 days after explosion
Comments: Accepted for publication in ApJ (14 pages, 6 figures, 3 tables)
Submitted: 2017-11-03, last modified: 2018-04-12
The late-time light curves of Type Ia supernovae (SNe Ia), observed $>900$ days after explosion, present the possibility of a new diagnostic for SN Ia progenitor and explosion models. First, however, we must discover what physical process (or combination of processes) leads to the slow-down of the late-time light curve relative to a pure $^{56}$Co decay, as observed in SNe 2011fe, 2012cg, and 2014J. We present Hubble Space Telescope observations of SN 2015F, taken $\approx 600-1040$ days past maximum light. Unlike those of the three other SNe Ia, the light curve of SN 2015F remains consistent with being powered solely by the radioactive decay of $^{56}$Co. We fit the light curves of these four SNe Ia in a consistent manner and measure possible correlations between the light curve stretch - a proxy for the intrinsic luminosity of the SN - and the parameters of the physical model used in the fit (e.g., the mass ratio of $^{56}$Co and $^{57}$Co produced in the explosion, or the time at which freeze-out sets in). We propose a new, late-time Phillips-like correlation between the stretch of the SNe and the shape of their late-time light curves, which we parametrize as the difference between their pseudo-bolometric luminosities at 600 and 900 days: $\Delta L_{900} = {\rm log}(L_{600}/L_{900})$. This model-independent correlation provides a new way to test which physical process lies behind the slow-down of SN Ia light curves $>900$ days after explosion, and, ultimately, fresh constraints on the various SN Ia progenitor and explosion models.
[19]  [pdf] - 1686642
MOSFiT: Modular Open-Source Fitter for Transients
Comments: 18 pages, 6 figures, accepted to ApJS. Revision cleans up descriptions and adds additional code examples. Documentation available at
Submitted: 2017-10-05, last modified: 2018-02-20
Much of the progress made in time-domain astronomy is accomplished by relating observational multi-wavelength time series data to models derived from our understanding of physical laws. This goal is typically accomplished by dividing the task in two: collecting data (observing), and constructing models to represent that data (theorizing). Owing to the natural tendency for specialization, a disconnect can develop between the best available theories and the best available data, potentially delaying advances in our understanding new classes of transients. We introduce MOSFiT: the Modular Open-Source Fitter for Transients, a Python-based package that downloads transient datasets from open online catalogs (e.g., the Open Supernova Catalog), generates Monte Carlo ensembles of semi-analytical light curve fits to those datasets and their associated Bayesian parameter posteriors, and optionally delivers the fitting results back to those same catalogs to make them available to the rest of the community. MOSFiT is designed to help bridge the gap between observations and theory in time-domain astronomy; in addition to making the application of existing models and creation of new models as simple as possible, MOSFiT yields statistically robust predictions for transient characteristics, with a standard output format that includes all the setup information necessary to reproduce a given result. As large-scale surveys such as LSST discover entirely new classes of transients, tools such as MOSFiT will be critical for enabling rapid comparison of models against data in statistically consistent, reproducible, and scientifically beneficial ways.
[20]  [pdf] - 1620649
The magnetar model for Type I superluminous supernovae I: Bayesian analysis of the full multicolour light curve sample with MOSFiT
Submitted: 2017-06-02, last modified: 2018-01-16
We use the new Modular Open Source Fitter for Transients (MOSFiT) to model 38 hydrogen-poor superluminous supernovae (SLSNe). We fit their multicolour light curves with a magnetar spin-down model and present the posterior distributions of magnetar and ejecta parameters. The colour evolution can be well matched with a simple absorbed blackbody. We find the following medians (1$\sigma$ ranges): spin period 2.4 ms (1.2-4 ms); magnetic field $0.8\times 10^{14}$ G (0.2-1.8 $\times 10^{14}$ G); ejecta mass 4.8 Msun (2.2-12.9 Msun); kinetic energy $3.9\times 10^{51}$ erg (1.9-9.8 $\times 10^{51}$ erg). This significantly narrows the parameter space compared to our priors, showing that although the model is flexible, the parameter space relevant to SLSNe is well constrained by existing data. The requirement that the instantaneous engine power is $\sim 10^{44}$ erg at the light curve peak necessitates either a large rotational energy (P<2 ms), or more commonly that the spin-down and diffusion timescales be well-matched. We find no evidence for separate populations of fast- and slow-declining SLSNe, which instead form a continuum both in light curve widths and inferred parameters. Variations in the spectra are well explained through differences in spin-down power and photospheric radii at maximum-light. We find no correlations between any model parameters and the properties of SLSN host galaxies. Comparing our posteriors to stellar evolution models, we show that SLSNe require rapidly rotating (fastest 10%) massive stars (> 20 Msun), and that this is consistent with the observed SLSN rate. High mass, low metallicity, and likely binary interaction all serve to maintain rapid rotation essential for magnetar formation. By reproducing the full set of SLSN light curves, our posteriors can be used to inform photometric searches for SLSNe in future survey data.
[21]  [pdf] - 1630156
A dependence of the tidal disruption event rate on global stellar surface mass density and stellar velocity dispersion
Comments: Accepted for publication in ApJ
Submitted: 2017-07-10, last modified: 2017-12-21
The rate of tidal disruption events (TDEs), $R_\text{TDE}$, is predicted to depend on stellar conditions near the super-massive black hole (SMBH), which are on difficult-to-measure sub-parsec scales. We test whether $R_\text{TDE}$ depends on kpc-scale global galaxy properties, which are observable. We concentrate on stellar surface mass density, $\Sigma_{M_\star}$, and velocity dispersion, $\sigma_v$, which correlate with the stellar density and velocity dispersion of the stars around the SMBH. We consider 35 TDE candidates, with and without known X-ray emission. The hosts range from star-forming to quiescent to quiescent with strong Balmer absorption lines. The last (often with post-starburst spectra) are overrepresented in our sample by a factor of $35^{+21}_{-17}$ or $18^{+8}_{-7}$, depending on the strength of the H$\delta$ absorption line. For a subsample of hosts with homogeneous measurements, $\Sigma_{M_\star}=10^9$-$10^{10}~{\rm M_\odot / kpc^2}$, higher on average than for a volume-weighted control sample of Sloan Digital Sky Survey galaxies with similar redshifts and stellar masses. This is because: (1) most of the TDE hosts here are quiescent galaxies, which tend to have higher $\Sigma_{M_\star}$ than the star-forming galaxies that dominate the control, and (2) the star-forming hosts have higher average $\Sigma_{M_\star}$ than the star-forming control. There is also a weak suggestion that TDE hosts have lower $\sigma_v$ than for the quiescent control. Assuming that $R_{\rm TDE}\propto \Sigma_{M_\star}^\alpha \times \sigma_v^\beta$, and applying a statistical model to the TDE hosts and control sample, we estimate $\hat{\alpha}=0.9 \pm 0.2$ and $\hat{\beta}=-1.0 \pm 0.6$. This is broadly consistent with $R_\text{TDE}$ being tied to the dynamical relaxation of stars surrounding the SMBH.
[22]  [pdf] - 1605053
The Combined Ultraviolet, Optical, and Near-Infrared Light Curves of the Kilonova Associated with the Binary Neutron Star Merger GW170817: Unified Data Set, Analytic Models, and Physical Implications
Comments: 35 pages; 5 figures; 3 tables. Significant changes made to Table 3. Accepted to ApJL
Submitted: 2017-10-31, last modified: 2017-11-28
We present the first effort to aggregate, homogenize, and uniformly model the combined ultraviolet, optical, and near-infrared dataset for the electromagnetic counterpart of the binary neutron star merger GW170817. By assembling all of the available data from 18 different papers and 46 different instruments, we are able to identify and mitigate systematic offsets between individual datasets, and to identify clear outlying measurements, with the resulting pruned and adjusted dataset offering an opportunity to expand the study of the kilonova. The unified dataset includes 647 individual flux measurements, spanning 0.45 to 29.4 days post-merger, and thus has greater constraining power for physical models than any single dataset. We test a number of semi-analytical models and find that the data are well modeled with a three-component kilonova model: a "blue" lanthanide-poor component with Mej~0.020 Msol and vej~0.27c; an intermediate opacity "purple" component with Mej~0.047 Msol and vej~0.15c; and a "red" lanthanide-rich component with Mej~0.011 Msol and vej~0.14c. We further explore the possibility of ejecta asymmetry and its impact on the estimated parameters. From the inferred parameters we draw conclusions about the physical mechanisms responsible for the various ejecta components, the properties of the neutron stars, and, combined with an up-to-date merger rate, the implications for r-process enrichment via this channel. To facilitate future studies of this keystone event we make the unified dataset and our modeling code public.
[23]  [pdf] - 1589828
The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. II. UV, Optical, and Near-IR Light Curves and Comparison to Kilonova Models
Cowperthwaite, P. S.; Berger, E.; Villar, V. A.; Metzger, B. D.; Nicholl, M.; Chornock, R.; Blanchard, P. K.; Fong, W.; Margutti, R.; Soares-Santos, M.; Alexander, K. D.; Allam, S.; Annis, J.; Brout, D.; Brown, D. A.; Butler, R. E.; Chen, H. -Y.; Diehl, H. T.; Doctor, Z.; Drout, M. R.; Eftekhari, T.; Farr, B.; Finley, D. A.; Foley, R. J.; Frieman, J. A.; Fryer, C. L.; García-Bellido, J.; Gill, M. S. S.; Guillochon, J.; Herner, K.; Holz, D. E.; Kasen, D.; Kessler, R.; Marriner, J.; Matheson, T.; Neilsen,, E. H.; Quataert, E.; Palmese, A.; Rest, A.; Sako, M.; Scolnic, D. M.; Smith, N.; Tucker, D. L.; Williams, P. K. G.; Balbinot, E.; Carlin, J. L.; Cook, E. R.; Durret, F.; Li, T. S.; Lopes, P. A. A.; Lourenço, A. C. C.; Marshall, J. L.; Medina, G. E.; Muir, J.; Muñoz, R. R.; Sauseda, M.; Schlegel, D. J.; Secco, L. F.; Vivas, A. K.; Wester, W.; Zenteno, A.; Zhang, Y.; Abbott, T. M. C.; Banerji, M.; Bechtol, K.; Benoit-Lévy, A.; Bertin, E.; Buckley-Geer, E.; Burke, D. L.; Capozzi, D.; Rosell, A. Carnero; Kind, M. Carrasco; Castander, F. J.; Crocce, M.; Cunha, C. E.; D'Andrea, C. B.; da Costa, L. N.; Davis, C.; DePoy, D. L.; Desai, S.; Dietrich, J. P.; Drlica-Wagner, A.; Eifler, T. F.; Evrard, A. E.; Fernandez, E.; Flaugher, B.; Fosalba, P.; Gaztanaga, E.; Gerdes, D. W.; Giannantonio, T.; Goldstein, D. A.; Gruen, D.; Gruendl, R. A.; Gutierrez, G.; Honscheid, K.; Jain, B.; James, D. J.; Jeltema, T.; Johnson, M. W. G.; Johnson, M. D.; Kent, S.; Krause, E.; Kron, R.; Kuehn, K.; Kuropatkin, N.; Lahav, O.; Lima, M.; Lin, H.; Maia, M. A. G.; March, M.; Martini, P.; McMahon, R. G.; Menanteau, F.; Miller, C. J.; Miquel, R.; Mohr, J. J.; Neilsen, E.; Nichol, R. C.; Ogando, R. L. C.; Plazas, A. A.; Roe, N.; Romer, A. K.; Roodman, A.; Rykoff, E. S.; Sanchez, E.; Scarpine, V.; Schindler, R.; Schubnell, M.; Sevilla-Noarbe, I.; Smith, M.; Smith, R. C.; Sobreira, F.; Suchyta, E.; Swanson, M. E. C.; Tarle, G.; Thomas, D.; Thomas, R. C.; Troxel, M. A.; Vikram, V.; Walker, A. R.; Wechsler, R. H.; Weller, J.; Yanny, B.; Zuntz, J.
Comments: 13 Pages, 3 Figures, 2 Tables. ApJL, In Press. Keywords: GW170817, LVC
Submitted: 2017-10-16
We present UV, optical, and NIR photometry of the first electromagnetic counterpart to a gravitational wave source from Advanced LIGO/Virgo, the binary neutron star merger GW170817. Our data set extends from the discovery of the optical counterpart at $0.47$ days to $18.5$ days post-merger, and includes observations with the Dark Energy Camera (DECam), Gemini-South/FLAMINGOS-2 (GS/F2), and the {\it Hubble Space Telescope} ({\it HST}). The spectral energy distribution (SED) inferred from this photometry at $0.6$ days is well described by a blackbody model with $T\approx 8300$ K, a radius of $R\approx 4.5\times 10^{14}$ cm (corresponding to an expansion velocity of $v\approx 0.3c$), and a bolometric luminosity of $L_{\rm bol}\approx 5\times10^{41}$ erg s$^{-1}$. At $1.5$ days we find a multi-component SED across the optical and NIR, and subsequently we observe rapid fading in the UV and blue optical bands and significant reddening of the optical/NIR colors. Modeling the entire data set we find that models with heating from radioactive decay of $^{56}$Ni, or those with only a single component of opacity from $r$-process elements, fail to capture the rapid optical decline and red optical/NIR colors. Instead, models with two components consistent with lanthanide-poor and lanthanide-rich ejecta provide a good fit to the data, the resulting "blue" component has $M_\mathrm{ej}^\mathrm{blue}\approx 0.01$ M$_\odot$ and $v_\mathrm{ej}^\mathrm{blue}\approx 0.3$c, and the "red" component has $M_\mathrm{ej}^\mathrm{red}\approx 0.04$ M$_\odot$ and $v_\mathrm{ej}^\mathrm{red}\approx 0.1$c. These ejecta masses are broadly consistent with the estimated $r$-process production rate required to explain the Milky Way $r$-process abundances, providing the first evidence that BNS mergers can be a dominant site of $r$-process enrichment.
[24]  [pdf] - 1584320
iPTF15eqv: Multi-wavelength Expos\'e of a Peculiar Calcium-rich Transient
Comments: 24 pages, 16 figures. Closely matches version published in The Astrophysical Journal
Submitted: 2017-06-06, last modified: 2017-10-12
The progenitor systems of the class of "Ca-rich transients" is a key open issue in time domain astrophysics. These intriguing objects exhibit unusually strong calcium line emissions months after explosion, fall within an intermediate luminosity range, are often found at large projected distances from their host galaxies, and may play a vital role in enriching galaxies and the intergalactic medium. Here we present multi-wavelength observations of iPTF15eqv in NGC 3430, which exhibits a unique combination of properties that bridge those observed in Ca-rich transients and Type Ib/c supernovae. iPTF15eqv has among the highest [Ca II]/[O I] emission line ratios observed to date, yet is more luminous and decays more slowly than other Ca-rich transients. Optical and near-infrared photometry and spectroscopy reveal signatures consistent with the supernova explosion of a < 10 solar mass star that was stripped of its H-rich envelope via binary interaction. Distinct chemical abundances and ejecta kinematics suggest that the core collapse occurred through electron capture processes. Deep limits on possible radio emission made with the Jansky Very Large Array imply a clean environment ($n <$ 0.1 cm$^{-3}$) within a radius of $\sim 10^{17}$ cm. Chandra X-ray Observatory observations rule out alternative scenarios involving tidal disruption of a white dwarf by a black hole, for masses > 100 solar masses). Our results challenge the notion that spectroscopically classified Ca-rich transients only originate from white dwarf progenitor systems, complicate the view that they are all associated with large ejection velocities, and indicate that their chemical abundances may vary widely between events.
[25]  [pdf] - 1585151
The superluminous supernova SN 2017egm in the nearby galaxy NGC 3191: a metal-rich environment can support a typical SLSN evolution
Comments: Accepted for publication in ApJL
Submitted: 2017-06-27, last modified: 2017-07-27
At redshift z=0.03, the recently-discovered SN 2017egm is the nearest Type I superluminous supernova (SLSN) to date, and first near the center of a massive spiral galaxy (NGC 3191). Using SDSS spectra of NGC 3191, we find a metallicity ~2 Z$_\odot$ at the nucleus and ~1.3 Z$_\odot$ for a star forming region at a radial offset similar to SN 2017egm. Archival radio-to-UV photometry reveals a star-formation rate ~15 M$_\odot$ yr$^{-1}$ (with ~70% dust-obscured), which can account for a Swift X-ray detection, and stellar mass ~$10^{10.7}$ M$_\odot$. We model the early UV-optical light curves with a magnetar central-engine model, using the Bayesian light curve fitting tool MOSFiT. The fits indicate ejecta mass 2-4 M$_\odot$, spin period 4-6 ms, magnetic field (0.7-1.7)$\times 10^{14}$G, and kinetic energy 1-2 $\times10^{51}$ erg. These parameters are consistent with the overall distributions for SLSNe, modeled by Nicholl et al (2017), although the derived mass and spin are towards the low end, possibly indicating enhanced loss of mass and angular momentum before explosion. This has two implications: (i) SLSNe can occur at solar metallicity, although with a low fraction ~10%; and (ii) metallicity has at most a modest effect on their properties. Both conclusions are in line with results for long gamma-ray bursts. Assuming a monotonic rise gives an explosion date MJD $57889\pm1$. However, a short-lived excess in the data relative to the best-fitting models may indicate an early-time `bump'. If confirmed, SN 2017egm would be the first SLSN with a spectrum during the bump-phase; this shows the same O II lines seen at maximum light, which may be an important clue for explaining these bumps.
[26]  [pdf] - 1586380
Theoretical Models of Optical Transients. I. A Broad Exploration of the Duration-Luminosity Phase Space
Comments: 30 pages, 20 figures, submitted to ApJ; interactive figure:
Submitted: 2017-07-25
The duration-luminosity phase space of optical transients is used, mostly heuristically, to compare various classes of transient events, to explore the origin of new transients, and to influence optical survey observing strategies. For example, several observational searches have been guided by intriguing voids and gaps in this phase space. However we should ask: Do we expect to find transients in these voids given our understanding of the various heating sources operating in astrophysical transients? In this work, we explore a broad range of theoretical models and empirical relations to generate optical light curves and to populate the duration-luminosity phase space (DLPS). We explore transients powered by adiabatic expansion, radioactive decay, magnetar spin-down, and circumstellar interaction. For each heating source, we provide a concise summary of the basic physical processes, a physically motivated choice of model parameter ranges, an overall summary of the resulting light curves and their the occupied range in the DLPS, and how the various model input parameters affect the light curves. We specifically explore the key voids discussed in the literature: the intermediate luminosity gap between classical novae and supernovae, and short-duration transients (<10 days). We find that few physical models lead to transients that occupy these voids. Moreover, we find that only relativistic expansion can produce fast and luminous transients, while for all other heating sources, events with durations <10 days are dim (> -15 mag). Finally, we explore the detection potential of optical surveys (e.g., LSST) in the DLPS and quantify the notion that short-duration and dim transients are exponentially more difficult to discover in untargeted surveys.
[27]  [pdf] - 1585810
The Galactic Distribution of Fragments Formed from Tidally Disrupted Stars
Comments: 6 pages, 5 figures. Submitted to ApJL
Submitted: 2017-07-11
Approximately once every $10^4$ years, a star passes close enough to the supermassive black hole Sgr A* at the center of the Milky Way to be pulled apart by the black hole's tidal forces. The star is then "spaghettified" into a long stream of matter, with approximately one half being bound to Sgr A* and the other half unbound. Within this stream, the local self-gravity dominates the tidal field of Sgr A*, which at minimum restricts the stream to a small finite width. As the stream cools from adiabatic expansion and begins to recombine, the residual self-gravity allows for planetary-mass fragments to form along the length of the stream; these fragments are then shot out into the galaxy at range of velocities, with the fastest moving at ~10% c. We determine the phase space distributions of these fragments for a realistic ensemble of stellar disruptions, along with the local density of fragments in the solar neighborhood. We find that ~$10^7$ fragments produced by Sgr A* accumulate within the Milky Way over its lifetime, that there are ~$10^7$ fragments that lie within 1 Mpc of the Milky Way originating from other galaxies, and that the nearest fragment to our Sun is on average 500 pc distant.
[28]  [pdf] - 1579685
Optical Thermonuclear Transients From Tidal Compression of White Dwarfs as Tracers of the Low End of the Massive Black Hole Mass Function
Comments: This version corrects a normalization error in the optical/radio afterglow lightcurves in Fig 12
Submitted: 2015-08-10, last modified: 2017-06-14
In this paper, we model the observable signatures of tidal disruptions of white dwarf (WD) stars by massive black holes (MBHs) of moderate mass, $\approx 10^3 - 10^5 M_\odot$. When the WD passes deep enough within the MBH's tidal field, these signatures include thermonuclear transients from burning during maximum compression. We combine a hydrodynamic simulation that includes nuclear burning of the disruption of a $0.6 M_\odot$ C/O WD with a Monte Carlo radiative transfer calculation to synthesize the properties of a representative transient. The transient's emission emerges in the optical, with lightcurves and spectra reminiscent of type I SNe. The properties are strongly viewing-angle dependent, and key spectral signatures are $\approx 10,000$ km s$^{-1}$ Doppler shifts due to the orbital motion of the unbound ejecta. Disruptions of He WDs likely produce large quantities of intermediate-mass elements, offering a possible production mechanism for Ca-rich transients. Accompanying multiwavelength transients are fueled by accretion and arise from the nascent accretion disk and relativistic jet. If MBHs of moderate mass exist with number densities similar to those of supermassive MBHs, both high energy wide-field monitors and upcoming optical surveys should detect tens to hundreds of WD tidal disruptions per year. The current best strategy for their detection may therefore be deep optical follow up of high-energy transients of unusually-long duration. The detection rate or the non-detection of these transients by current and upcoming surveys can thus be used to place meaningful constraints on the extrapolation of the MBH mass function to moderate masses.
[29]  [pdf] - 1581304
Low-mass White Dwarfs with Hydrogen Envelopes as a Missing Link in the Tidal Disruption Menu
Comments: 19 pages, 13 figures. Published in ApJ
Submitted: 2017-01-27, last modified: 2017-06-05
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.
[30]  [pdf] - 1582082
PS16dtm: A Tidal Disruption Event in a Narrow-line Seyfert 1 Galaxy
Comments: 19 pages, 17 figures, Submitted to ApJ
Submitted: 2017-03-22
[Abridged] We present observations of PS16dtm, a luminous transient that occurred at the nucleus of a known Narrow-line Seyfert 1 galaxy hosting a 10$^6$ M$_\odot$ black hole. The transient was previously claimed to be a Type IIn SLSN due to its luminosity and hydrogen emission lines. The light curve shows that PS16dtm brightened by about two magnitudes in ~50 days relative to the archival host brightness and then exhibited a plateau phase for about 100 days followed by the onset of fading in the UV. During the plateau PS16dtm showed no color evolution, maintained a blackbody temperature of 1.7 x 10$^4$ K, and radiated at approximately $L_{Edd}$ of the SMBH. The spectra exhibit multi-component hydrogen emission lines and strong FeII emission, show little evolution with time, and closely resemble the spectra of NLS1s while being distinct from those of Type IIn SNe. Moreover, PS16dtm is undetected in the X-rays to a limit an order of magnitude below an archival X-ray detection of its host galaxy. These observations strongly link PS16dtm to activity associated with the SMBH and are difficult to reconcile with a SN origin or any known form of AGN variability, and therefore we argue that it is a TDE in which the accretion of the stellar debris powers the rise in the continuum and excitation of the pre-existing broad line region, while providing material that obscures the X-ray emitting region of the pre-existing AGN accretion disk. A detailed TDE model fit to the light curve indicates that PS16dtm will remain bright for several years; we further predict that the X-ray emission will reappear on a similar timescale as the accretion rate declines. Finally, we place PS16dtm in the context of other TDEs and find that TDEs in AGN galaxies are an order of magnitude more efficient and reach Eddington luminosities, likely due to interaction of the stellar debris with the pre-existing accretion disk.
[31]  [pdf] - 1554091
Optical/UV-to-X-Ray Echoes from the Tidal Disruption Flare ASASSN-14li
Comments: Publish in ApJ Letters
Submitted: 2017-03-20
We carried out the first multi-wavelength (optical/UV and X-ray) photometric reverberation mapping of a tidal disruption flare (TDF) ASASSN-14li. We find that its X-ray variations are correlated with and lag the optical/UV fluctuations by 32$\pm$4 days. Based on the direction and the magnitude of the X-ray time lag, we rule out X-ray reprocessing and direct emission from a standard circular thin disk as the dominant source of its optical/UV emission. The lag magnitude also rules out an AGN disk-driven instability as the origin of ASASSN-14li and thus strongly supports the tidal disruption picture for this event and similar objects. We suggest that the majority of the optical/UV emission likely originates from debris stream self-interactions. Perturbations at the self-interaction sites produce optical/UV variability and travel down to the black hole where they modulate the X-rays. The time lag between the optical/UV and the X-rays variations thus correspond to the time taken by these fluctuations to travel from the self-interaction site to close to the black hole. We further discuss these time lags within the context of the three variants of the self-interaction model. High-cadence monitoring observations of future TDFs will be sensitive enough to detect these echoes and would allow us to establish the origin of optical/UV emission in TDFs in general.
[32]  [pdf] - 1547792
Enabling New ALMA Science with Improved Support for Time-Domain Observations
Comments: 9 pages; whitepaper submitted to the ALMA Science Advisory Council; corresponding author P. K. G. Williams (
Submitted: 2017-03-14
While the Atacama Large Millimeter/submillimeter Array (ALMA) is a uniquely powerful telescope, its impact in certain fields of astrophysics has been limited by observatory policies rather than the telescope's innate technical capabilities. In particular, several observatory policies present challenges for observations of variable, mobile, and/or transient sources --- collectively referred to here as "time-domain" observations. In this whitepaper we identify some of these policies, describe the scientific applications they impair, and suggest changes that would increase ALMA's science impact in Cycle 6 and beyond. Parties interested in time-domain science with ALMA are encouraged to join the ALMA Time-domain Special Interest Group (ATSIG) by signing up for the ATSIG mailing list at .
[33]  [pdf] - 1564078
The fine line between total and partial tidal disruption events
Comments: 8 pages, 8 figures, 2 tables; accepted to A&A
Submitted: 2017-02-24, last modified: 2017-03-10
Flares from tidal disruption events are unique tracers of quiescent black holes at the centre of galaxies. The appearance of these flares is very sensitive to whether the star is totally or partially disrupted, and in this paper we seek to identify the critical distance of the star from the black hole (r_d) that enables us to distinguish between these two outcomes. We perform here Mesh-free Finite Mass, traditional, and modern Smoothed Particle Hydrodynamical simulations of star-black hole close encounters, with the aim of checking if the value of r_d depends on the simulation technique. We find that the critical distance (or the so-called critical disruption parameter beta_d) depends only weakly on the adopted simulation method, being beta_d=0.92\pm 0.02 for a gamma=5/3 polytrope and beta_d=2.01\pm 0.01 for a gamma=4/3 polytrope.
[34]  [pdf] - 1535355
A likely decade-long sustained tidal disruption event
Comments: Published in Nature Astronomy on Feb 6. Including supplementary materials
Submitted: 2017-02-02
Multiwavelength flares from tidal disruption and accretion of stars can be used to find and study otherwise dormant massive black holes in galactic nuclei. Previous well-monitored candidate flares are short-lived, with most emission confined to within ~1 year. Here we report the discovery of a well observed super-long (>11 years) luminous soft X-ray flare from the nuclear region of a dwarf starburst galaxy. After an apparently fast rise within ~4 months a decade ago, the X-ray luminosity, though showing a weak trend of decay, has been persistently high at around the Eddington limit (when the radiation pressure balances the gravitational force). The X-ray spectra are generally soft (steeply declining towards higher energies) and can be described with Comptonized emission from an optically thick low-temperature corona, a super-Eddington accretion signature often observed in accreting stellar-mass black holes. Dramatic spectral softening was also caught in one recent observation, implying either a temporary transition from the super-Eddington accretion state to the standard thermal state or the presence of a transient highly blueshifted (~0.36c) warm absorber. All these properties in concert suggest a tidal disruption event of an unusually long super-Eddington accretion phase that has never been observed before.
[35]  [pdf] - 1561040
New physical insights about Tidal Disruption Events from a comprehensive observational inventory at X-ray wavelengths
Comments: Submitted to ApJ. Text, figures and tables updated. 76 pages, 24 Figures, and 22 Tables. We make our results publicly available at the open TDE catalog:
Submitted: 2016-11-07, last modified: 2017-01-10
We perform a comprehensive study of the X-ray emission from 70 transient sources which have been classified as a tidal disruption event (TDE) in the literature. We explore the properties of these candidates using nearly three decades of X-ray observations to quantify the properties and characteristics of X-ray TDEs observationally. We find that the emission from X-ray TDEs increase by two to four orders of magnitude compared to pre-flare constraints, which evolves significantly with time and decays with powerlaw indices that are typically shallower than the canonical $t^{-5/3}$ decay law, implying that X-ray TDEs are viscously delayed. These events exhibit enhanced column densities relative to Galactic and are quite soft in nature, with no strong correlation between the amount of detected soft and hard emission. At peak, jetted events have an X-ray to optical ratio $\gg$1, while non-jetted events have a ratio $\sim$1, which suggests that these events undergo reprocessing at different rates. X-ray TDEs have long $T_{90}$ values consistent with that expected from a viscously-driven accretion disk formed by the disruption of a main-sequence star by a black hole with a mass $<$$10^{7}M_{\odot}$. The isotropic luminosities of X-ray TDEs is bimodal such that jetted and non-jetted events are separated by a "reprocessing valley" which we suggest is naturally populated by optical/UV TDEs that most likely produce X-rays, but due to reprocessing this emission is "veiled" from observations. Our results suggest that non-jetted X-ray TDEs likely originate from partial disruptions and/or disruptions of low mass stars.
[36]  [pdf] - 1531743
Simulations of Magnetic Fields in Tidally-Disrupted Stars
Comments: Accepted to ApJL. 7 pages, 4 figures
Submitted: 2016-09-26, last modified: 2016-12-22
We perform the first magnetohydrodynamical simulations of tidal disruptions of stars by supermassive black holes. We consider stars with both tangled and ordered magnetic fields, for both grazing and deeply disruptive encounters. When the star survives disruption, we find its magnetic field amplifies by a factor of up to twenty, but see no evidence for the a self-sustaining dynamo that would yield arbitrary field growth. For stars that do not survive, and within the tidal debris streams produced in partial disruptions, we find that the component of the magnetic field parallel to the direction of stretching along the debris stream only decreases slightly with time, eventually resulting in a stream where the magnetic pressure is in equipartition with the gas. Our results suggest that the returning gas in most (if not all) stellar tidal disruptions is already highly magnetized by the time it returns to the black hole.
[37]  [pdf] - 1530658
An Open Catalog for Supernova Data
Comments: 26 pages, 6 figures, 1 table. Updated in response to referee comments, added Luke Zoltan Kelley to author list. Catalog available at
Submitted: 2016-05-03, last modified: 2016-11-07
We present the Open Supernova Catalog, an online collection of observations and metadata for presently 36,000+ supernovae and related candidates. The catalog is freely available on the web (, with its main interface having been designed to be a user-friendly, rapidly-searchable table accessible on desktop and mobile devices. In addition to the primary catalog table containing supernova metadata, an individual page is generated for each supernova which displays its available metadata, light curves, and spectra spanning X-ray to radio frequencies. The data presented in the catalog is automatically rebuilt on a daily basis and is constructed by parsing several dozen sources, including the data presented in the supernova literature and from secondary sources such as other web-based catalogs. Individual supernova data is stored in the hierarchical, human- and machine-readable JSON format, with the entirety of each supernova's data being contained within a single JSON file bearing its name. The setup we present here, which is based upon open source software maintained via git repositories hosted on github, enables anyone to download the entirety of the supernova dataset to their home computer in minutes, and to make contributions of their own data back to the catalog via git. As the supernova dataset continues to grow, especially in the upcoming era of all-sky synoptic telescopes which will increase the total number of events by orders of magnitude, we hope that the catalog we have designed will be a valuable tool for the community to analyze both historical and contemporary supernovae.
[38]  [pdf] - 1530814
OGLE16aaa - a Signature of a Hungry Super Massive Black Hole
Comments: Accepted in MNRAS Letters
Submitted: 2016-06-09, last modified: 2016-10-12
We present the discovery and first three months of follow-up observations of a currently on-going unusual transient detected by the OGLE-IV survey, located in the centre of a galaxy at redshift z=0.1655. The long rise to absolute magnitude of -20.5 mag, slow decline, very broad He and H spectral features make OGLE16aaa similar to other optical/UV Tidal Disruption Events (TDEs). Weak narrow emission lines in the spectrum and archival photometric observations suggest the host galaxy is a weak-line Active Galactic Nucleus (AGN), which has been accreting at higher rate in the past. OGLE16aaa, along with SDSS J0748, seems to form a sub-class of TDEs by weakly or recently active super-massive black holes (SMBHs). This class might bridge the TDEs by quiescent SMBHs and flares observed as "changing-look QSOs", if we interpret the latter as TDEs. If this picture is true, the previously applied requirement for identifying a flare as a TDE that it had to come from an inactive nucleus, could be leading to observational bias in TDE selection, thus affecting TDE-rate estimations.
[39]  [pdf] - 1494915
The Astropy Problem
Muna, Demitri; Alexander, Michael; Allen, Alice; Ashley, Richard; Asmus, Daniel; Azzollini, Ruyman; Bannister, Michele; Beaton, Rachael; Benson, Andrew; Berriman, G. Bruce; Bilicki, Maciej; Boyce, Peter; Bridge, Joanna; Cami, Jan; Cangi, Eryn; Chen, Xian; Christiny, Nicholas; Clark, Christopher; Collins, Michelle; Comparat, Johan; Cook, Neil; Croton, Darren; Davids, Isak Delberth; Depagne, Éric; Donor, John; Santos, Leonardo A. dos; Douglas, Stephanie; Du, Alan; Durbin, Meredith; Erb, Dawn; Faes, Daniel; Fernández-Trincado, J. G.; Foley, Anthony; Fotopoulou, Sotiria; Frimann, Søren; Frinchaboy, Peter; Garcia-Dias, Rafael; Gawryszczak, Artur; George, Elizabeth; Gonzalez, Sebastian; Gordon, Karl; Gorgone, Nicholas; Gosmeyer, Catherine; Grasha, Katie; Greenfield, Perry; Grellmann, Rebekka; Guillochon, James; Gurwell, Mark; Haas, Marcel; Hagen, Alex; Haggard, Daryl; Haines, Tim; Hall, Patrick; Hellwing, Wojciech; Herenz, Edmund Christian; Hinton, Samuel; Hlozek, Renee; Hoffman, John; Holman, Derek; Holwerda, Benne Willem; Horton, Anthony; Hummels, Cameron; Jacobs, Daniel; Jensen, Jens Juel; Jones, David; Karick, Arna; Kelley, Luke; Kenworthy, Matthew; Kitchener, Ben; Klaes, Dominik; Kohn, Saul; Konorski, Piotr; Krawczyk, Coleman; Kuehn, Kyler; Kuutma, Teet; Lam, Michael T.; Lane, Richard; Liske, Jochen; Lopez-Camara, Diego; Mack, Katherine; Mangham, Sam; Mao, Qingqing; Marsh, David J. E.; Mateu, Cecilia; Maurin, Loïc; McCormac, James; Momcheva, Ivelina; Monteiro, Hektor; Mueller, Michael; Munoz, Roberto; Naidu, Rohan; Nelson, Nicholas; Nitschelm, Christian; North, Chris; Nunez-Iglesias, Juan; Ogaz, Sara; Owen, Russell; Parejko, John; Patrício, Vera; Pepper, Joshua; Perrin, Marshall; Pickering, Timothy; Piscionere, Jennifer; Pogge, Richard; Poleski, Radek; Pourtsidou, Alkistis; Price-Whelan, Adrian M.; Rawls, Meredith L.; Read, Shaun; Rees, Glen; Rein, Hanno; Rice, Thomas; Riemer-Sørensen, Signe; Rusomarov, Naum; Sanchez, Sebastian F.; Santander-García, Miguel; Sarid, Gal; Schoenell, William; Scholz, Aleks; Schuhmann, Robert L.; Schuster, William; Scicluna, Peter; Seidel, Marja; Shao, Lijing; Sharma, Pranav; Shulevski, Aleksandar; Shupe, David; Sifón, Cristóbal; Simmons, Brooke; Sinha, Manodeep; Skillen, Ian; Soergel, Bjoern; Spriggs, Thomas; Srinivasan, Sundar; Stevens, Abigail; Streicher, Ole; Suchyta, Eric; Tan, Joshua; Telford, O. Grace; Thomas, Romain; Tonini, Chiara; Tremblay, Grant; Tuttle, Sarah; Urrutia, Tanya; Vaughan, Sam; Verdugo, Miguel; Wagner, Alexander; Walawender, Josh; Wetzel, Andrew; Willett, Kyle; Williams, Peter K. G.; Yang, Guang; Zhu, Guangtun; Zonca, Andrea
Submitted: 2016-10-10
The Astropy Project ( 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.
[40]  [pdf] - 1457122
The X-ray through Optical Fluxes and Line Strengths of Tidal Disruption Events
Comments: Updated to match version accepted to ApJ
Submitted: 2015-10-28, last modified: 2016-08-10
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.
[41]  [pdf] - 1501917
Prompt Radiation and Mass Outflows from the Stream-Stream Collisions of Tidal Disruption Events
Comments: 14 pages, 12 figures, 2 tables, movies available at
Submitted: 2016-03-24
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.
[42]  [pdf] - 1374120
Discovery of an outflow from radio observations of the tidal disruption event ASASSN-14li
Comments: 10 pages, 6 figures, 2 tables; published in ApJ Letters. Fixes an error in our density calculation and expands discussion, major conclusions unchanged from previous version
Submitted: 2015-10-05, last modified: 2016-03-07
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, $\rho(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.
[43]  [pdf] - 1370336
An Ultraviolet Spectrum of the Tidal Disruption Flare ASASSN-14li
Comments: Accepted for publication in ApJ Letters. 8 pages, 4 figures
Submitted: 2016-01-13, last modified: 2016-02-02
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 $\Delta 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.
[44]  [pdf] - 1370299
The Gamma-ray Afterglows of Tidal Disruption Events
Comments: 12 pages, 3 figures, submitted. Comments welcome!
Submitted: 2015-12-18, last modified: 2015-12-23
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 $\gamma$-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 $\gamma$-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.
[45]  [pdf] - 1304221
Flows of X-ray gas reveal the disruption of a star by a massive black hole
Comments: Published in the Oct 22 2015 issue of Nature
Submitted: 2015-10-21
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.
[46]  [pdf] - 1290246
Observational Cosmology With Semi-Relativistic Stars
Comments: 5 pages, 5 figures. Accepted for publication in Annals of Mathematical Sciences and Applications
Submitted: 2014-11-18, last modified: 2015-10-08
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.
[47]  [pdf] - 1304052
Thermonuclear detonations ensuing white dwarf mergers
Comments: Accepted for publication in MNRAS
Submitted: 2015-08-10, last modified: 2015-10-01
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.
[48]  [pdf] - 1411171
Unbound Debris Streams and Remnants Resulting From the Tidal Disruptions of Stars by Supermassive Black Holes
Comments: 16 pages, 9 figures. Submitted to ApJ
Submitted: 2015-09-29
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.
[49]  [pdf] - 1284930
SETI via Leakage from Light Sails in Exoplanetary Systems
Comments: 6 pages, 5 figures. Accepted to ApJL. Revisions correct minor algebra errors and make bandwidth assumptions explicit
Submitted: 2015-08-12, last modified: 2015-09-03
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).
[50]  [pdf] - 1273058
A Dark Year for Tidal Disruption Events
Comments: Accepted to ApJ. Updated version adds new figure (Figure 9). 11 pages, 9 figures. Movie of stream intersections available at
Submitted: 2015-01-21, last modified: 2015-08-06
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]
[51]  [pdf] - 1205699
The Fastest Unbound Stars in the Universe
Comments: 20 pages, 18 figures. Accepted by ApJ
Submitted: 2014-11-18, last modified: 2015-06-16
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.
[52]  [pdf] - 1124034
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
Comments: 16 pages, 13 figures (2 new figures in revised version). Published in ApJ. Latest version incorporates erratum that fixes issue with fitting formulae not including enough significant digits
Submitted: 2012-06-11, last modified: 2015-04-05
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]
[53]  [pdf] - 968972
Relativistic effects on tidal disruption kicks of solitary stars
Comments: 10 pages, 1 table, 11 figures; accepted for publication in MNRAS
Submitted: 2015-02-06, last modified: 2015-02-17
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 $\gamma$ = 5/3 and $\gamma$ = 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 $\beta$ alone, reaching at most the escape velocity of the original star. For a given $\beta$, 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-$\beta$ encounters are statistically more likely than high-$\beta$ 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.
[54]  [pdf] - 928751
Death by Dynamics: Planetoid-Induced Explosions on White Dwarfs
Comments: 15 pages, 1 figure, 1 table
Submitted: 2015-01-30
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.
[55]  [pdf] - 1222723
A Luminous, Fast Rising UV-Transient Discovered by ROTSE: a Tidal Disruption Event?
Comments: 13 pages, accepted for publication in ApJ; some references added
Submitted: 2014-10-22, last modified: 2014-10-31
We present follow-up observations of an optical transient (OT) discovered by ROTSE on Jan. 21, 2009. Photometric monitoring was carried out with ROTSE-IIIb in the optical and Swift in the UV up to +70 days after discovery. The light curve showed a fast rise time of ~10 days followed by a steep decline over the next 60 days, which was much faster than that implied by 56Ni - 56Co radioactive decay. The SDSS DR10 database contains a faint, red object at the position of the OT, which appears slightly extended. This and other lines of evidence suggest that the OT is of extragalactic origin, and this faint object is likely the host galaxy. A sequence of optical spectra obtained with the 9.2-m Hobby-Eberly Telescope (HET) between +8 and +45 days after discovery revealed a hot, blue continuum with no visible spectral features. A few weak features that appeared after +30 days probably originated from the underlying host. Fitting synthetic templates to the observed spectrum of the host galaxy revealed a redshift of z = 0.19. At this redshift the peak magnitude of the OT is close to -22.5, similar to the brightest super-luminous supernovae; however, the lack of identifiable spectral features makes the massive stellar death hypothesis less likely. A more plausible explanation appears to be the tidal disruption of a sun-like star by the central super-massive black hole. We argue that this transient likely belongs to a class of super-Eddington tidal disruption events.
[56]  [pdf] - 1209390
Illuminating Massive Black Holes With White Dwarfs: Orbital Dynamics and High Energy Transients from Tidal Interactions
Comments: Revised following response from referee
Submitted: 2014-05-06, last modified: 2014-07-06
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.
[57]  [pdf] - 1202718
Possible Origin of the G2 Cloud from the Tidal Disruption of a Known Giant Star by Sgr A*
Comments: 5 pages, 3 figures, 1 table. Revised version includes list of candidate disrupted giants and clarification on the numerical method. Accepted to ApJL
Submitted: 2014-01-13, last modified: 2014-04-03
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-$\gamma$ 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.
[58]  [pdf] - 1171472
Turbovelocity Stars: Kicks Resulting From the Tidal Disruption of Solitary Stars
Comments: 6 pages, 5 figures. Published in ApJL, includes erratum which corrects algebra errors in Section 3.2
Submitted: 2013-05-20, last modified: 2014-01-16
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.
[59]  [pdf] - 1166176
PS1-10jh: The Disruption of a Main-Sequence Star of Near-Solar Composition
Comments: 21 pages, 11 figures, 1 table. Accepted to ApJ, revised version addresses referee's comments and includes 1 new figure
Submitted: 2013-04-23, last modified: 2014-01-16
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.
[60]  [pdf] - 1172636
Spoon-Feeding Giant Stars to Supermassive Black Holes: Episodic Mass Transfer from Evolving Stars and Their Contribution to the Quiescent Activity of Galactic Nuclei
Comments: 12 pages, 5 figures
Submitted: 2013-07-10, last modified: 2013-08-29
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.
[61]  [pdf] - 679279
Conditions For Successful Helium Detonations In Astrophysical Environments
Comments: 9 pages, 8 figures. Accepted for publication in the Astrophysical Journal
Submitted: 2013-02-25, last modified: 2013-05-10
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 ...
[62]  [pdf] - 1164581
Circumstellar absorption in double detonation Type Ia supernovae
Comments: Accepted for publication in The Astrophysical Journal Letters
Submitted: 2013-02-12, last modified: 2013-05-06
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.
[63]  [pdf] - 1157633
On the Survivability and Metamorphism of Tidally Disrupted Giant Planets: the Role of Dense Cores
Comments: Accepted for publication in ApJ. 15 pages, 9 figures, 3 tables. Two movies at and
Submitted: 2012-11-08
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.
[64]  [pdf] - 1124100
The Tidal Disruption of Giant Stars and Their Contribution to the Flaring Supermassive Black Hole Population
Comments: 18 pages, 18 figures, submitted to ApJ
Submitted: 2012-06-13
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.
[65]  [pdf] - 1123207
The dynamics, appearance and demographics of relativistic jets triggered by tidal disruption of stars in quiescent supermassive black holes
Comments: 16 pages, 16 figures, submitted to ApJ
Submitted: 2012-05-07
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]
[66]  [pdf] - 1092885
How the merger of two white dwarfs depends on their mass ratio: orbital stability and detonations at contact
Comments: 13 pages, 10 figures, submitted to MNRAS
Submitted: 2012-01-11
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.
[67]  [pdf] - 1053425
A relativistic jetted outburst from a massive black hole fed by a tidally disrupted star
Comments: Submitted, 32 pages including supplemental online material
Submitted: 2011-04-16
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.
[68]  [pdf] - 1051705
Prelude to a double degenerate merger: the onset of mass transfer and its impact on gravitational waves and surface detonations
Comments: 20 pages, 25 figures, submitted to ApJ. Movies available at
Submitted: 2011-01-26
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.
[69]  [pdf] - 1042498
Consequences of the Ejection and Disruption of Giant Planets
Comments: 17 pages, 14 figures, submitted to ApJ. Movies and high-resolution snapshots available at
Submitted: 2010-12-10
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.
[70]  [pdf] - 1018187
Surface Detonations in Double Degenerate Binary Systems Triggered by Accretion Stream Instabilities
Comments: 6 pages, 4 figures, 1 table. Submitted to ApJL. For a high-resolution version, movies, and other supporting material see
Submitted: 2009-11-03
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-$\alpha$ 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.
[71]  [pdf] - 1002959
Collisions of white dwarfs as a new progenitor channel for type Ia supernovae
Comments: 3 figures, minor changes according to referee comments, ApJL in press
Submitted: 2009-07-18, last modified: 2009-10-01
We present the results of a systematic numerical study of an alternative progenitor scenario to produce type Ia supernova explosions, which is not restricted to the ignition of a CO white dwarf near the Chandrasekhar mass. In this scenario, a shock-triggered thermonuclear explosion ensues from the collision of two white dwarfs. Consistent modeling of the gas dynamics together with nuclear reactions using both a smoothed particle and a grid-based hydrodynamics code are performed to study the viability of this alternative progenitor channel. We find that shock-triggered ignition and the synthesis of Ni are in fact a natural outcome for moderately massive white dwarf pairs colliding close to head-on. We use a multi-dimensional radiative transfer code to calculate the emergent broadband light curves and spectral time-series of these events. The synthetic spectra and lightcurves compare well with those of normal type Ia supernovae over a similar B-band decline rate and are broadly consistent with the Phillips relation, although a mild dependence on viewing angle is observed due to the asymmetry of the ejected debris. While event rates within galactic centers and globular clusters are found to be much too low to explain the bulk of the type Ia supernovae population, they may be frequent enough to make as much as a one percent contribution to the overall rate. Although these rate estimates are still subject to substantial uncertainties, they do suggest that dense stellar systems should provide upcoming supernova surveys with hundreds of such collision-induced thermonuclear explosions per year.
[72]  [pdf] - 1001060
Three-Dimensional Simulations of Tidally Disrupted Solar-Type Stars and the Observational Signatures of Shock Breakout
Comments: 10 pages, 7 figures, Accepted to ApJ. For movies see
Submitted: 2008-11-10, last modified: 2009-09-17
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.