Source author record

Benjamin J. Shappee

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

34works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

34 published item(s)

preprint2026arXiv

JCMT Constraints on the Early-Time HCN and CO Emission and HCN Temporal Evolution of 3I/ATLAS

Interstellar objects (ISOs), particularly those with cometary activity, provide unique insight into the primordial physical and chemical conditions present during the formation of the planetary system in which they originated. Observations in the sub-mm regime allow for direct measurements of several parent molecules released from the comet nucleus into the coma. Here we present observations of the third ISO, 3I/ATLAS, with the `Ū`ū heterodyne receiver on the James Clerk Maxwell Telescope (JCMT), which targeted emission from HCN($J = 3 - 2$) and CO($J = 2 - 1$). Our observations, taken between 16 July 2025 and 21 July 2025 (UT), when 3I/ATLAS was at a heliocentric distance between 4.01 and 3.84 au, provide the earliest sub-mm constraints on its activity. We do not detect HCN or CO in these epochs, with 3$σ$ upper-limits on the production rates of $Q(HCN) < 1.7 \times 10^{24}$ s$^{-1}$ at $r_h = 4.01 - 3.97$ au and $Q(CO) < 1.1 \times 10^{27}$ s$^{-1}$ at $r_h = 3.94 - 3.84$ au, respectively. We combine this HCN limit with later JCMT observations of HCN to constrain its temporal evolution. Fitting the HCN detections with a $Q(HCN) \propto r_h^{-n}$ model and accounting for the upper-limits yields $n = 12.7^{+6.9}_{-2.5}$. This slope is steeper than those of typical Solar System comets, but consistent with the production rate slopes measured for other species in the coma of 3I/ATLAS.

preprint2022arXiv

ASAS-SN follow-up of IceCube high-energy neutrino alerts

We report on the search for optical counterparts to IceCube neutrino alerts released between April 2016 and August 2021 with the All-Sky Automated Survey for SuperNovae (ASAS-SN). Despite the discovery of a diffuse astrophysical high-energy neutrino flux in 2013, the source of those neutrinos remains largely unknown. Since 2016, IceCube has published likely-astrophysical neutrinos as public realtime alerts. Through a combination of normal survey and triggered target-of-opportunity observations, ASAS-SN obtained images within 1 hour of the neutrino detection for 20% (11) of all observable IceCube alerts and within one day for another 57% (32). For all observable alerts, we obtained images within at least two weeks from the neutrino alert. ASAS-SN provides the only optical follow-up for about 17% of IceCube's neutrino alerts. We recover the two previously claimed counterparts to neutrino alerts, the flaring-blazar TXS 0506+056 and the tidal disruption event AT2019dsg. We investigate the light curves of previously-detected transients in the alert footprints, but do not identify any further candidate neutrino sources. We also analysed the optical light curves of Fermi 4FGL sources coincident with high-energy neutrino alerts, but do not identify any contemporaneous flaring activity. Finally, we derive constraints on the luminosity functions of neutrino sources for a range of assumed evolution models.

preprint2022arXiv

COol Companions ON Ultrawide orbiTS (COCONUTS). III. A Very Red L6 Benchmark Brown Dwarf around a Young M5 Dwarf

We present the third discovery from the COol Companions ON Ultrawide orbiTS (COCONUTS) program, the COCONUTS-3 system, composed of a young M5 primary star UCAC4 374-046899 and a very red L6 dwarf WISEA J081322.19$-$152203.2. These two objects have a projected separation of 61$''$ (1891 au) and are physically associated given their common proper motions and estimated distances. The primary star, COCONUTS-3A, has a mass of $0.123\pm0.006$ M$_{\odot}$ and we estimate its age as 100 Myr to 1 Gyr based on its kinematics and spectrophotometric properties. We derive its metallicity as $0.21 \pm 0.07$ dex using empirical calibrations established by older higher-gravity M dwarfs and find this [Fe/H] could be slightly underestimated according to PHOENIX models given COCONUTS-3A's younger age. The companion, COCONUTS-3B, has a near-infrared spectral type of L6$\pm$1 INT-G, and we infer physical properties of $T_{\rm eff} = 1362^{+48}_{-73}$ K, $\log{(g)}= 4.96^{+0.15}_{-0.34}$ dex, $R = 1.03^{+0.12}_{-0.06}$ R$_{\rm Jup}$, and $M = 39^{+11}_{-18}$ M$_{\rm Jup}$, using its bolometric luminosity, its host star's age, and hot-start evolution models. We construct cloudy atmospheric model spectra at the evolution-based physical parameters and compare them to COCONUTS-3B's spectrophotometry. We find this companion possesses ample condensate clouds in its photosphere with the data-model discrepancies likely due to the models using an older version of the opacity database. Compared to field-age L6 dwarfs, COCONUTS-3B has fainter absolute magnitudes and a 120 K cooler $T_{\rm eff}$. Also, the J-K color of this companion is among the reddest for ultracool benchmarks with ages older than a few 100 Myr. COCONUTS-3 likely formed in the same fashion as stellar binaries given the companion-to-host mass ratio of 0.3 and represents a valuable benchmark to quantify the systematics of substellar model atmospheres.

preprint2022arXiv

Exploring metallicity-dependent rates of Type Ia supernovae and their impact on galaxy formation

Type Ia supernovae play a critical role in stellar feedback and elemental enrichment in galaxies. Recent transient surveys like the All-Sky Automated Survey for Supernova (ASAS-SN) and the Dark Energy Survey (DES) find that the specific Ia rate at z ~ 0 may be ~ 15-50 times higher in lower-mass galaxies than at Milky Way-mass. Independently, Milky Way observations show that the close-binary fraction of solar-type stars is higher at lower metallicity. Motivated by these observations, we use the FIRE-2 cosmological zoom-in simulations to explore the impact of varying Ia rate models, including metallicity dependence, on galaxies across a range of stellar masses: 10^7 Msun - 10^{11} Msun. First, we benchmark our simulated star-formation histories (SFHs) against observations. We show that assumed SFHs and stellar mass functions play a major role in determining the degree of tension between observations and metallicity-independent Ia rate models, and potentially cause ASAS-SN and DES observations to be much more consistent with each other than might naively appear. Models in which the Ia rate increases with decreasing metallicity (as ~ Z^{-0.5} to Z^{-1}) provide significantly better agreement with observations. Encouragingly, these increases in Ia rate (> 10 times in low-mass galaxies) do not significantly impact galaxy stellar masses and morphologies: effective radii, axis ratios, and v/sigma remain largely unaffected except for our most extreme rate models. We explore implications for both [Fe/H] and [alpha/Fe] enrichment: metallicity-dependent Ia rate models can improve agreement with observed stellar mass-metallicity relations in low-mass galaxies. Our results demonstrate that a wide range of metallicity-dependent Ia models are viable for galaxy formation and motivate future work in this area.

preprint2022arXiv

Investigating the Nature of the Luminous Ambiguous Nuclear Transient ASASSN-17jz

We present observations of the extremely luminous but ambiguous nuclear transient (ANT) ASASSN-17jz, spanning roughly 1200 days of the object's evolution. ASASSN-17jz was discovered by the All-Sky Automated Survey for Supernovae (ASAS-SN) in the galaxy SDSS J171955.84+414049.4 on UT 2017 July 27 at a redshift of $z=0.1641$. The transient peaked at an absolute $B$-band magnitude of $M_{B,{\rm peak}}=-22.81$, corresponding to a bolometric luminosity of $L_{\rm bol,peak}=8.3\times10^{44}$~erg~s$^{-1}$, and exhibited late-time ultraviolet emission that was still ongoing in our latest observations. Integrating the full light curve gives a total emitted energy of $E_{\rm tot}=(1.36\pm0.08)\times10^{52}$~erg, with $(0.80\pm0.02)\times10^{52}$~erg of this emitted within 200 days of peak light. This late-time ultraviolet emission is accompanied by increasing X-ray emission that becomes softer as it brightens. ASASSN-17jz exhibited a large number of spectral emission lines most commonly seen in active galactic nuclei (AGNs) with little evidence of evolution. It also showed transient Balmer features which became fainter and broader over time, and are still being detected $>1000$ days after peak brightness. We consider various physical scenarios for the origin of the transient, including supernovae (SNe), tidal disruption events (TDEs), AGN outbursts, and ANTs. We find that the most likely explanation is that ASASSN-17jz was an SN~IIn occurring in or near the disk of an existing AGN, and that the late-time emission is caused by the AGN transitioning to a more active state.

preprint2022arXiv

Variability Selected Active Galactic Nuclei from ASAS-SN Survey: Constraining the Low Luminosity AGN Population

Low luminosity active galactic nuclei (LLAGN) probe accretion physics in the low Eddington regime and can provide additional clues about galaxy evolution. AGN variability is ubiquitous and thus provides a reliable tool for finding AGN. We analyze the All-Sky Automated Survey for SuperNovae light curves of 1218 galaxies with $g<14$ mag and Sloan Digital Sky Survey spectra in search of AGN. We find 37 objects that are both variable and have AGN-like structure functions, which is about 3% of the sample. The majority of the variability selected AGN are LLAGN with Eddington ratios ranging from $10^{-4}$ to $10^{-2}$. We thus estimate the fraction of LLAGN in the population of galaxies as 2% down to a median Eddington ratio of $2\times10^{-3}$. Combining the BPT line ratio diagnostics and the broad-line AGN, up to $\sim$60% of the AGN candidates are confirmed spectroscopically. The BPT diagnostics also classified 10-30% of the candidates as star forming galaxies rather than AGN.

preprint2021arXiv

Discovery of a Highly Eccentric, Chromospherically Active Binary: ASASSN-V J192114.84+624950.8

As part of an All-Sky Automated Survey for SuperNovae (ASAS-SN) search for sources with large flux decrements, we discovered a transient where the quiescent, stellar source, ASASSN-V J192114.84+624950.8, rapidly decreased in flux by $\sim55\%$ ($\sim0.9$ mag) in the g-band. The \textit{TESS} light curve revealed that the source is a highly eccentric, eclipsing binary. Fits to the light curve using \textsc{phoebe} find the binary orbit to have $e=0.79$, $P_{\rm orb}=18.462~\text{days}$, and $i=88.6^{\circ}$ and the ratios of the stellar radii and temperatures to be $R_2/R_1 = 0.71$ and $T_{e,2}/T_{e,1} = 0.82$. Both stars are chromospherically active, allowing us to determine their rotational periods of $P_1=1.52$ days and $P_2=1.79$ days, respectively. A LBT/MODS spectrum shows that the primary is a late-G or early-K type dwarf. Fits to the SED show that the luminosities and temperatures of the two stars are $L_1 = 0.48~L_{\sun}$, $T_1= 5050~K$, $L_2 = 0.12~L_{\sun}$, and $T_{2} = 4190~K$. We conclude that ASASSN-V J192114.84+624950.8 consists of two chromospherically active, rotational variable stars in a highly elliptical eclipsing orbit.

preprint2021arXiv

Late-Onset Circumstellar Medium Interactions are Rare: An Unbiased GALEX View of Type Ia Supernovae

Using ultraviolet (UV) light curves we constrain the circumstellar environments of 1080 Type Ia supernovae (SNe Ia) within $z<0.5$ from archival Galaxy Evolution Explorer (GALEX) observations. All SNe Ia are required to have pre- and post-explosion GALEX observations to ensure adequate subtraction of the host-galaxy flux. Using the late-time GALEX observations we look for the UV excess expected from any interaction between the SN ejecta and circumstellar material (CSM). Four SNe Ia are detected near maximum light and we compare the GALEX photometry to archival data, but we find none of our targets show convincing evidence of CSM interaction. A recent Hubble Space Telescope (HST) survey estimates that $\sim6\%$ of SNe Ia may interact with distant CSM, but statistical inferences are complicated by the small sample size and selection effects. By injecting model light curves into our data and then recovering them, we constrain a broad range of CSM interactions based on the CSM interaction start time and the maximum luminosity. Combining our GALEX non-detections with the HST results, we constrain occurrence of late-onset CSM interaction among SNe Ia with moderate CSM interaction, similar to that observed in PTF11kx, to $f_\text{CSM}\lesssim5.1\%$ between $0-500$ days after discovery and $\lesssim2.7\%$ between $500-1000$ days after discovery at $90\%$ confidence. For weaker CSM interactions similar to SN 2015cp, we obtain limits of $\lesssim16\%$ and $\lesssim4.8\%$, respectively, for the same time ranges.

preprint2021arXiv

The Rapid X-ray and UV Evolution of ASASSN-14ko

ASASSN-14ko is a recently discovered periodically flaring transient at the center of the AGN ESO 253-G003 with a slowly decreasing period. Here we show that the flares originate from the northern, brighter nucleus in this dual-AGN, post-merger system. The light curves for the two flares that occurred in May 2020 and September 2020 are nearly identical over all wavelengths. For both events, Swift observations showed that the UV and optical wavelengths brightened in unison. The effective temperature of the UV/optical emission rises and falls with the increase and subsequent decline in the luminosity. The X-ray flux, in contrast, first rapidly drops over $\sim$2.6 days, rises for $\sim$5.8 days, drops again over $\sim$4.3 days and then recovers. The X-ray spectral evolution of the two flares differ, however. During the May 2020 peak the spectrum softened with increases in the X-ray luminosity, while we observed the reverse for the September 2020 peak.

preprint2020arXiv

COol Companions ON Ultrawide orbiTS (COCONUTS). I. A High-Gravity T4 Benchmark around an Old White Dwarf and A Re-Examination of the Surface-Gravity Dependence of the L/T Transition

We present the first discovery from the COol Companions ON Ultrawide orbiTS (COCONUTS) program, a large-scale survey for wide-orbit planetary and substellar companions. We have discovered a co-moving system COCONUTS-1, composed of a hydrogen-dominated white dwarf (PSO J058.9855+45.4184; $d=31.5$ pc) and a T4 companion (PSO J058.9869+45.4296) at a $40.6''$ (1280 au) projected separation. We derive physical properties for COCONUTS-1B from (1) its near-infrared spectrum using cloudless Sonora atmospheric models, and (2) its luminosity and the white dwarf's age ($7.3_{-1.6}^{+2.8}$ Gyr) using Sonora evolutionary models. The two methods give consistent temperatures and radii, but atmospheric models infer a lower surface gravity and therefore an unphysically young age. Assuming evolutionary model parameters ($T_{\rm eff}=1255^{+6}_{-8}$ K, $\log{g}=5.44^{+0.02}_{-0.03}$ dex, $R=0.789^{+0.011}_{-0.005}$ R$_{\rm Jup}$), we find cloudless model atmospheres have brighter Y- and J-band fluxes than the data, suggesting condensate clouds have not fully dispersed around 1300 K. The W2 flux (4.6 $μ$m) of COCONUTS-1B is fainter than models, suggesting non-equilibrium mixing of CO. To investigate the gravity dependence of the L/T transition, we compile all 60 known L6-T6 benchmarks and derive a homogeneous set of temperatures, surface gravities, and masses. As is well-known, young, low-gravity late-L dwarfs have significantly fainter, redder near-infrared photometry and $\approx200-300$ K cooler temperatures than old, high-gravity objects. Our sample now reveals such gravity dependence becomes weaker for T dwarfs, with young objects having comparable near-infrared photometry and $\approx100$ K cooler temperatures compared to old objects. Finally, we find that young objects have a larger amplitude J-band brightening than old objects, and also brighten at H band as they cross the L/T transition.

preprint2020arXiv

Direct evidence for shock-powered optical emission in a nova

Classical novae are thermonuclear explosions that occur on the surfaces of white dwarf stars in interacting binary systems (Bode & Evans 2008). It has long been thought that the luminosity of classical novae is powered by continued nuclear burning on the surface of the white dwarf after the initial runaway (Gallaher & Starrfield 1978). However, recent observations of GeV $γ$-rays from classical novae have hinted that shocks internal to the nova ejecta may dominate the nova emission. Shocks have also been suggested to power the luminosity of events as diverse as stellar mergers (Metzger & Pejcha 2017), supernovae (Moriya et al. 2018), and tidal disruption events (Roth et al. 2016), but observational confirmation has been lacking. Here we report simultaneous space-based optical and $γ$-ray observations of the 2018 nova V906 Carinae (ASASSN-18fv), revealing a remarkable series of distinct correlated flares in both bands. The optical and $γ$-ray flares occur simultaneously, implying a common origin in shocks. During the flares, the nova luminosity doubles, implying that the bulk of the luminosity is shock-powered. Furthermore, we detect concurrent but weak X-ray emission from deeply embedded shocks, confirming that the shock power does not appear in the X-ray band and supporting its emergence at longer wavelengths. Our data, spanning the spectrum from radio to $γ$-ray, provide direct evidence that shocks can power substantial luminosity in classical novae and other optical transients.

preprint2020arXiv

Does Gravity Fall Down? Evidence for Gravitational Wave Deflection Along the Line of Sight to GW 170817

We present a novel test of general relativity (GR): measuring the geometric component of the time delay due to gravitational lensing. GR predicts that photons and gravitational waves follow the same geodesic paths and thus experience the same geometric time delay. We show that for typical systems, the time delays are tens of seconds, and thus can dominate over astrophysical delays in the timing of photon emission. For the case of GW 170817, we use a multi-plane lensing code to evaluate the time delay due to four massive halos along the line of sight. From literature mass and distance measurements of these halos, we establish at high confidence (significantly greater than 5 sigma) that the gravitational waves of GW 170817 underwent gravitational deflection to arrive within 1.7 seconds of the photons.

preprint2020arXiv

H$α$ Luminosity of ATLAS18qtd Does Not Plateau in the Nebular Phase

We present new spectroscopic and photometric observations of ATLAS18qtd/SN 2018cqj, a fast-declining Type Ia supernova with variable H$α$ emission in previously-published nebular phase spectra. ATLAS18qtd is undetected in both spectroscopic and photometric observations which occurred at $\sim 540~\rm d$ after maximum light and $\sim 230~\rm d$ after the last H$α$ detection. With these new non-detections, we place an upper limit on the H$α$ luminosity of $\lesssim 1.1\times 10^{36}~\rm{erg}~\rm s^{-1}$ indicating the H$α$ flux decreased by a factor of $\gtrsim 4$ since the previous detection. This upper limit excludes H$α$ emission that plateaus or increases since the previous detection but cannot confirm that the H$α$ emission decay rate is equivalent to the supernova decay rate.

preprint2020arXiv

The First Habitable Zone Earth-sized Planet from TESS. I: Validation of the TOI-700 System

We present the discovery and validation of a three-planet system orbiting the nearby (31.1 pc) M2 dwarf star TOI-700 (TIC 150428135). TOI-700 lies in the TESS continuous viewing zone in the Southern Ecliptic Hemisphere; observations spanning 11 sectors reveal three planets with radii ranging from 1 R$_\oplus$ to 2.6 R$_\oplus$ and orbital periods ranging from 9.98 to 37.43 days. Ground-based follow-up combined with diagnostic vetting and validation tests enable us to rule out common astrophysical false-positive scenarios and validate the system of planets. The outermost planet, TOI-700 d, has a radius of $1.19\pm0.11$ R$_\oplus$ and resides in the conservative habitable zone of its host star, where it receives a flux from its star that is approximately 86% of the Earth's insolation. In contrast to some other low-mass stars that host Earth-sized planets in their habitable zones, TOI-700 exhibits low levels of stellar activity, presenting a valuable opportunity to study potentially-rocky planets over a wide range of conditions affecting atmospheric escape. While atmospheric characterization of TOI-700 d with the James Webb Space Telescope (JWST) will be challenging, the larger sub-Neptune, TOI-700 c (R = 2.63 R$_\oplus$), will be an excellent target for JWST and beyond. TESS is scheduled to return to the Southern Hemisphere and observe TOI-700 for an additional 11 sectors in its extended mission, which should provide further constraints on the known planet parameters and searches for additional planets and transit timing variations in the system.

preprint2020arXiv

The Outburst of the Young Star Gaia19bey

We report photometry and spectroscopy of the outburst of the young stellar object Gaia19bey. We have established the outburst light curve with archival Gaia G, ATLAS Orange, ZTF r-band and Pan-STARRS rizy-filter photometry, showing an outburst of approximately 4 years duration, longer than typical EXors but shorter than FUors. Its pre-outburst SED shows a flat far-infrared spectrum, confirming the early evolutionary state of Gaia19bey and its similarity to other deeply embedded young stars experiencing outbursts. A lower limit to the peak outburst luminosity is approximately 182 L_sun at an assumed distance of 1.4 kpc, the minimum plausible distance. Infrared and optical spectroscopy near maximum light showed an emission line spectrum, including HI lines, strong red CaII emission, other metal emission lines, infrared CO bandhead emission, and a strong infrared continuum. Towards the end of the outburst, the emission lines have all but disappeared and the spectrum has changed into an almost pure continuum spectrum. This indicates a cessation of magnetospheric accretion activity. The near-infrared colors have become redder as Gaia19bey has faded, indicating a cooling of the continuum component. Near the end of the outburst, the only remaining strong emission lines are forbidden shock-excited emission lines. Adaptive optics integral field spectroscopy shows the H_2 1--0 S(1) emission with the morphology of an outflow cavity and the extended emission in the [FeII] line at 1644 nm with the morphology of an edge-on disk. However, we do not detect any large-scale jet from Gaia19bey.

preprint2020arXiv

The Rise and Fall of ASASSN-18pg: Following a TDE from Early To Late Times

We present nearly 500 days of observations of the tidal disruption event ASASSN-18pg, spanning from 54 days before peak light to 441 days after peak light. Our dataset includes X-ray, UV, and optical photometry, optical spectroscopy, radio observations, and the first published spectropolarimetric observations of a TDE. ASASSN-18pg was discovered on 2018 July 11 by the All-Sky Automated Survey for Supernovae (ASAS-SN) at a distance of $d=78.6$ Mpc, and with a peak UV magnitude of $m\simeq14$ it is both one of the nearest and brightest TDEs discovered to-date. The photometric data allow us to track both the rise to peak and the long-term evolution of the TDE. ASASSN-18pg peaked at a luminosity of $L\simeq2.2\times10^{44}$ erg s$^{-1}$, and its late-time evolution is shallower than a flux $\propto t^{-5/3}$ power-law model, similar to what has been seen in other TDEs. ASASSN-18pg exhibited Balmer lines and spectroscopic features consistent with Bowen fluorescence prior to peak which remained detectable for roughly 225 days after peak. Analysis of the two-component H$α$ profile indicates that, if they are the result of reprocessing of emission from the accretion disk, the different spectroscopic lines may be coming from regions between $\sim10$ and $\sim60$ light-days from the black hole. No X-ray emission is detected from the TDE and there is no evidence of a jet or strong outflow detected in the radio. Our spectropolarimetric observations give no strong evidence for significant asphericity in the emission region, with the emission region having an axis ratio of at least $\sim0.65$.

preprint2019arXiv

A Catalog of M-dwarf Flares with ASAS-SN

We analyzed the light curves of 1376 early-to-late, nearby M dwarfs to search for white-light flares using photometry from the All-Sky Automated Survey for Supernovae (ASAS-SN). We identified 480 M dwarfs with at least one potential flare employing a simple statistical algorithm that searches for sudden increases in $V$-band flux. After more detailed evaluation, we identified 62 individual flares on 62 stars. The event amplitudes range from $0.12 <ΔV < 2.04$ mag. Using classical-flare models, we place lower limits on the flare energies and obtain $V$-band energies spanning $2.0\times10^{30} \lesssim E_{V} \lesssim 6.9\times10^{35}$ erg. The fraction of flaring stars increases with spectral type, and most flaring stars show moderate to strong H$α$ emission. Additionally, we find that 14 of the 62 flaring stars are rotational variables, and they have shorter rotation periods and stronger H$α$ emission than non-flaring rotational variable M dwarfs.

preprint2019arXiv

The New EXor Outburst of ESO-H-alpha~99 observed by Gaia ATLAS and TESS

We report photometry and spectroscopy of the outburst of the young stellar object ESO-Halpha 99. The outburst was first noticed in Gaia alert Gaia18dvc and later by ATLAS. We have established the outburst light curve with archival ATLAS ``Orange'' filter photometry, Gaia data, new V-band photometry, and J, H, and K_s photometry from IRIS and UKIRT. The brightness has fluctuated several times near the light curve maximum. The TESS satellite observed ESO-Halpha 99 with high cadence during one of these minor minima and found brightness fluctuations on timescales of days and hours. Imaging with UKIRT shows the outline of an outflow cavity, and we find one knot of H_2 1-0 S(1) emission, now named MHO 1520, on the symmetry axis of this nebula, indicating recent collimated outflow activity from ESO-Halpha 99. Its pre-outburst SED shows a flat FIR spectrum, confirming its early evolutionary state and its similarity to other deeply embedded objects in the broader EXor class. The pre-outburst luminosity is 34 L_sun, a much higher luminosity than typical EXors, indicating that ESO-Halpha 99 may be a star of intermediate mass. Infrared and optical spectroscopy show a rich emission line spectrum, including HI lines, strong red CaII emission, as well as infrared CO bandhead emission, all characteristic EXors in the broadest sense. Comparison of the present spectra with an optical spectrum obtained in 1993, presumably in the quiescent state of the object, shows that during the present outburst the continuum component of the spectrum has increased notably more than the emission lines. The Halpha equivalent width during the outburst is down to one half of its 1993 level, and shock-excited emission lines are much less prominent.

preprint2017arXiv

Where Is the Flux Going? The Long-Term Photometric Variability of Boyajian's Star

We present ~800 days of photometric monitoring of Boyajian's Star (KIC 8462852) from the All-Sky Automated Survey for Supernovae (ASAS-SN) and ~4000 days of monitoring from the All Sky Automated Survey (ASAS). We show that from 2015 to the present the brightness of Boyajian's Star has steadily decreased at a rate of 6.3 +/- 1.4 mmag yr^-1, such that the star is now 1.5% fainter than it was in February 2015. Moreover, the longer time baseline afforded by ASAS suggests that Boyajian's Star has also undergone two brightening episodes in the past 11 years, rather than only exhibiting a monotonic decline. We analyze a sample of ~1000 comparison stars of similar brightness located in the same ASAS-SN field and demonstrate that the recent fading is significant at >99.4% confidence. The 2015-2017 dimming rate is consistent with that measured with Kepler data for the time period from 2009 to 2013. This long-term variability is difficult to explain with any of the physical models for the star's behavior proposed to date.

preprint2016arXiv

An Extreme Analogue of $ε$ Aurigae: An M-giant Eclipsed Every 69 Years by a Large Opaque Disk Surrounding a Small Hot Source

We present TYC 2505-672-1 as a newly discovered and remarkable eclipsing system comprising an M-type red giant that undergoes a ~3.45 year long, near-total eclipse (depth of ~4.5 mag) with a very long period of ~69.1 yr. TYC 2505-672-1 is now the longest-period eclipsing binary system yet discovered, more than twice as long as that of the currently longest-period system, $ε$ Aurigae. We show from analysis of the light curve including both our own data and historical data spanning more than 120 yr and from modeling of the spectral energy distribution, both before and during eclipse, that the red giant primary is orbited by a moderately hot source (T$_{eff}$~8000 K) that is itself surrounded by an extended, opaque circumstellar disk. From the measured ratio of luminosities, the radius of the hot companion must be in the range 0.1-0.5 Rsun (depending on the assumed radius of the red giant primary), which is an order of magnitude smaller than that for a main sequence A star and 1-2 orders of magnitude larger than that for a white dwarf. The companion is therefore most likely a "stripped red giant" subdwarf-B type star destined to become a He white dwarf. It is however somewhat cooler than most sdB stars, implying a very low mass for this "pre-He-WD" star. The opaque disk surrounding this hot source may be a remnant of the stripping of its former hydrogen envelope. However, it is puzzling how this object became stripped, given that it is at present so distant (orbital semi-major axis of ~24 AU) from the current red giant primary star. Extrapolating from our calculated ephemeris, the next eclipse should begin in early UT 2080 April and end in mid UT 2083 September (eclipse center UT 2081 December 24). This system is poised to become an exemplar of a very rare class of systems, even more extreme in several respects than the well studied archetype $ε$ Aurigae.

preprint2016arXiv

ASASSN-16ae: A Powerful White-Light Flare on an Early-L Dwarf

We report the discovery and classification of SDSS~J053341.43+001434.1 (SDSS0533), an early-L dwarf first discovered during a powerful $ΔV < -11$ magnitude flare observed as part of the ASAS-SN survey. Optical and infrared spectroscopy indicate a spectral type of L0 with strong H$α$ emission and a blue NIR spectral slope. Combining the photometric distance, proper motion, and radial velocity of SDSS0533 yields three-dimensional velocities of $(U,V,W)=(14\pm13,-35\pm14,-94\pm22)$~km~s$^{-1}$, indicating that it is most likely part of the thick disk population and probably old. The three detections of SDSS0533 obtained during the flare are consistent with a total $V$-band flare energy of at least $4.9\times10^{33}$~ergs (corresponding to a total thermal energy of at least $E_{\rm tot}>3.7\times10^{34}$~erg), placing it among the strongest detected M dwarf flares. The presence of this powerful flare on an old L0 dwarf may indicate that stellar-type magnetic activity persists down to the end of the main sequence and on older ML transition dwarfs.

preprint2016arXiv

DM Ori: A Young Star Occulted by a Disturbance in its Protoplanetary Disk

In some planet formation theories, protoplanets grow gravitationally within a young star's protoplanetary disk, a signature of which may be a localized disturbance in the disk's radial and/or vertical structure. Using time-series photometric observations by the Kilodegree Extremely Little Telescope South (KELT-South) project and the All-Sky Automated Survey for SuperNovae (ASAS-SN), combined with archival observations, we present the discovery of two extended dimming events of the young star, DM Ori. This young system faded by $\sim$1.5 mag from 2000 March to 2002 August and then again in 2013 January until 2014 September (depth $\sim$1.7 mag). We constrain the duration of the 2000-2002 dimming to be $<$ 860 days, and the event in 2013-2014 to be $<$ 585 days, separated by $\sim$12.5 years. A model of the spectral energy distribution (SED) indicates a large infrared excess consistent with an extensive circumstellar disk. Using basic kinematic arguments, we propose that DM Ori is likely being periodically occulted by a feature (possibly a warp or perturbation) in its circumstellar disk. In this scenario, the occulting feature is located $>$6 AU from the host star, moving at $\sim$14.6 km/s, and is $\sim$4.9 AU in width. This localized structure may indicate a disturbance such as may be caused by a protoplanet early in its formation.

preprint2016arXiv

Hello Darkness My Old Friend: The Fading of the Nearby TDE ASASSN-14ae

We present late-time optical spectroscopy taken with the Large Binocular Telescope's Multi-Object Double Spectrograph, an improved ASAS-SN pre-discovery non-detection, and late-time SWIFT observations of the nearby ($d=193$ Mpc, $z=0.0436$) tidal disruption event (TDE) ASASSN-14ae. Our observations span from $\sim$20 days before to $\sim$750 days after discovery. The proximity of ASASSN-14ae allows us to study the optical evolution of the flare and the transition to a host dominated state with exceptionally high precision. We measure very weak H$α$ emission 300 days after discovery ($L_{\rm Hα} \simeq 4\times 10^{39}$ ergs s$^{-1}$) and the most stringent upper limit to date on the H$α$ luminosity $\sim$750 days after discovery ($L_{\rm Hα} \lesssim 10^{39}$ ergs s$^{-1}$), suggesting that the optical emission arising from a TDE can vanish on a timescale as short as 1 year. Our results have important implications for both spectroscopic detection of TDE candidates at late times, as well as the nature of TDE host galaxies themselves.

preprint2016arXiv

The Afterglow and Early-Type Host Galaxy of the Short GRB 150101B at z=0.1343

We present the discovery of the X-ray and optical afterglows of the short-duration GRB 150101B, pinpointing the event to an early-type host galaxy at z=0.1343 +/- 0.0030. This makes GRB 150101B the most nearby short GRB with an early-type host galaxy discovered to date. Fitting the spectral energy distribution of the host galaxy results in an inferred stellar mass of ~7x10^10 M_sol, stellar population age of ~2-2.5 Gyr, and star formation rate of <0.4 M_sol yr^-1. The host of GRB 150101B is one of the largest and most luminous short GRB host galaxies, with a B-band luminosity of ~4.3L* and half-light radius of ~8 kpc. GRB 150101B is located at a projected distance of 7.35 +/- 0.07 kpc from its host center, and lies on a faint region of its host rest-frame optical light. Its location, combined with the lack of associated supernova, is consistent with a NS-NS/NS-BH merger progenitor. From modeling the evolution of the broad-band afterglow, we calculate isotropic-equivalent gamma-ray and kinetic energies of ~1.3x10^49 erg and ~(6-14)x10^51 erg, respectively, a circumburst density of ~(0.8-4)x10^-5 cm^-3, and a jet opening angle of >9 deg. Using observations extending to ~30 days, we place upper limits of <(2-4)x10^41 erg s^-1 on associated kilonova emission. We compare searches following previous short GRBs to existing kilonova models, and demonstrate the difficulty of performing effective kilonova searches from cosmological short GRBs using current ground-based facilities. We show that at the Advanced LIGO/VIRGO horizon distance of 200 Mpc, searches reaching depths of ~23-24 AB mag are necessary to probe a meaningful range of kilonova models.

preprint2016arXiv

The Eruption of the Candidate Young Star ASASSN-15qi

Outbursts on young stars are usually interpreted as accretion bursts caused by instabilities in the disk or the star-disk connection. However, some protostellar outbursts may not fit into this framework. In this paper, we analyze optical and near-infrared spectra and photometry to characterize the 2015 outburst of the probable young star ASASSN-15qi. The $\sim 3.5$ mag brightening in the $V$ band was sudden, with an unresolved rise time of less than one day. The outburst decayed exponentially by 1 mag for 6 days and then gradually back to the pre-outburst level after 200 days. The outburst is dominated by emission from $\sim10,000$ K gas. An explosive release of energy accelerated matter from the star in all directions, seen in a spectacular cool, spherical wind with a maximum velocity of 1000 km/s. The wind and hot gas both disappeared as the outburst faded and the source returned to its quiescent F-star spectrum. Nebulosity near the star brightened with a delay of 10-20 days. Fluorescent excitation of H$_2$ is detected in emission from vibrational levels as high as $v=11$, also with a possible time delay in flux increase. The mid-infrared spectral energy distribution does not indicate the presence of warm dust emission, although the optical photospheric absorption and CO overtone emission could be related to a gaseous disk. Archival photometry reveals a prior outburst in 1976. Although we speculate about possible causes for this outburst, none of the explanations are compelling.

preprint2014arXiv

Rapid Eccentricity Oscillations and the Mergers of Compact Objects in Hierarchical Triples

Kozai-Lidov (KL) oscillations can accelerate compact object mergers via gravitational wave (GW) radiation by driving the inner binaries of hierarchical triples to high eccentricities. We perform direct three-body integrations of high mass ratio compact object triple systems using Fewbody including post-Newtonian terms. We find that the inner binary undergoes rapid eccentricity oscillations (REOs) on the timescale of the outer orbital period which drive it to higher eccentricities than secular theory would otherwise predict, resulting in substantially reduced merger times. For a uniform distribution of tertiary eccentricity ($e_2$), ~40% of systems merge within ~1-2 eccentric KL timescales whereas secular theory predicts that only ~20% of such systems merge that rapidly. This discrepancy becomes especially pronounced at low $e_2$, with secular theory overpredicting the merger time by many orders of magnitude. We show that a non-negligible fraction of systems have eccentricity > 0.8 when they merge, in contrast to predictions from secular theory. Our results are applicable to high mass ratio triple systems containing black holes or neutron stars. In objects in which tidal effects are important, such as white dwarfs, stars, and planets, REOs can reduce the tidal circularization timescale by an order of magnitude and bring the components of the inner binary into closer orbits than would be possible in the secular approximation.

preprint2013arXiv

Characterizing a Dramatic $ΔV\sim-9$ Flare on an Ultracool Dwarf Found by the ASAS-SN Survey

We analyze a $ΔV\sim-9$ magnitude flare on the newly identified M8 dwarf SDSS J022116.84+194020.4 (hereafter SDSSJ0221) detected as part of the All-Sky Automated Survey for Supernovae (ASAS-SN). Using infrared and optical spectra, we confirm that SDSSJ0221 is a relatively nearby (d$\sim$76 pc) M8 dwarf with strong quiescent H$α$ emission. Based on kinematics and the absence of features consistent with low-gravity (young) ultracool dwarfs, we place a lower limit of 200 Myr on the age of SDSSJ0221. When modeled with a simple, classical flare light-curve, this flare is consistent with a total $U$-band flare energy $E_U\sim$ 10$^{34}$ erg, confirming that the most dramatic flares are not limited to warmer, more massive stars. Scaled to include a rough estimate of the emission line contribution to the $V$ band, we estimate a blackbody filling factor of $\sim$$10-30\%$ during the flare peak and $\sim$$0.5-1.6\%$ during the flare decay phase. These filling factors correspond to flare areas that are an order of magnitude larger than those measured for most mid-M dwarf flares.

preprint2013arXiv

Greatly Enhanced Eccentricity Oscillations in Quadruple Systems Composed of Two Binaries: Implications for Stars, Planets, and Transients

We study the orbital evolution of hierarchical quadruple systems composed of two binaries on a long mutual orbit, where each binary acts as a Kozai-Lidov (KL) perturber on the other. We find that the coupling between the two binaries qualitatively changes the behavior of their KL cycles. The binaries can experience coherent eccentricity oscillations as well as excursions to very high eccentricity that occur over a much larger fraction of the parameter space than in triple systems. For a ratio of outer to inner semi-major axes of 10 to 20, about 30 to 50% of equal-mass quadruples reach eccentricity 1-e < 10^{-3} in one of the binaries. This is about 4 to 12 times more than for triples with equivalent parameters. Orbital "flips" and collisions without previous tidal interaction are similarly enhanced in quadruples relative to triples. We argue that the frequency of evolutionary paths influenced by KL cycles is comparable in the triple and quadruple populations even though field quadruples are a factor of ~5 less frequent than triples. Additionally, quadruples might be a non-negligible source of triples and provide fundamentally new evolutionary outcomes involving close binaries, mergers, collisions, and associated transients, which occur without any fine tuning of parameters. Finally, we study the perturbations to a planetary orbit due to a distant binary and we find that the fraction of orbital flips is a factor of 3 to 4 higher than for the corresponding triple system given our fiducial parameters with implications for hot Jupiters and star-planet collisions.

preprint2013arXiv

No Stripped Hydrogen in the Nebular Spectra of Nearby Type Ia Supernova 2011fe

A generic prediction of the single-degenerate model for Type Ia supernovae (SNe Ia) is that a significant amount of material will be stripped from the donor star (~0.5 Msun for a giant donor and ~0.15 Msun for a main sequence donor) by the supernova ejecta. This material, excited by gamma-rays from radioactive decay, would then produce relatively narrow (1000 km s-1) emission features observable once the supernova enters the nebular phase. Such emission has never been detected, which already provides strong constraints on Type Ia progenitor models. In this Letter we report the deepest limit yet on the presence of H-alpha emission originating from the stripped hydrogen in the nebular spectrum of a Type Ia supernova obtained using a high signal-to-noise spectrum of the nearby normal SN Ia 2011fe 274 days after B-band maximum light with the Large Binocular Telescope's Multi-Object Double Spectrograph. We put a conservative upper limit on the H-alpha flux of 3.14 x 10^-17 erg s-1 cm-2, which corresponds to a luminosity of 1.57 x 10^35 erg s-1. By scaling models from the literature, our flux limit translates into an upper limit of 0.001 Msun of stripped material. This is an order of magnitude stronger than previous limits. SN 2011fe was a typical Type Ia supernova, special only in its proximity, and we argue that lack of hydrogen emission in its nebular spectrum adds yet another strong constraint on the single degenerate class of models for SNe Ia.

preprint2013arXiv

The Mass-Loss Induced Eccentric Kozai Mechanism: A New Channel for the Production of Close Compact Object-Stellar Binaries

Over a broad range of initial inclinations and eccentricities an appreciable fraction of hierarchical triple star systems with similar masses are essentially unaffected by the Kozai-Lidov mechanism (KM) until the primary in the central binary evolves into a compact object. Once it does, it may be much less massive than the other components in the ternary, enabling the "eccentric Kozai mechanism (EKM):" the mutual inclination between the inner and outer binary can flip signs driving the inner binary to very high eccentricity, leading to a close binary or collision. We demonstrate this "Mass-loss Induced Eccentric Kozai" (MIEK) mechanism by considering an example system and defining an ad-hoc minimal separation between the inner two members at which tidal affects become important. For fixed initial masses and semi-major axes, but uniform distributions of eccentricity and cosine of the mutual inclination, ~10% of systems interact tidally or collide while the primary is on the MS due to the KM or EKM. Those affected by the EKM are not captured by earlier quadrupole-order secular calculations. We show that fully ~30% of systems interact tidally or collide for the first time as the primary swells to AU scales, mostly as a result of the KM. Finally, ~2% of systems interact tidally or collide for the first time after the primary sheds most of its mass and becomes a WD, mostly as a result of the MIEK mechanism. These findings motivate a more detailed study of mass-loss in triple systems and the formation of close NS/WD-MS and NS/WD-NS/WD binaries without an initial common envelope phase.

preprint2013arXiv

The Mid-Infrared and Optical Decay of SN 2011fe

We measure the decay rate of the mid-IR luminosity from type Ia supernova 2011fe between six months and one year after explosion using Spitzer/IRAC observations. The fading in the 3.6 micron channel is 1.48+/-0.02 mag/100d, which is similar to that seen in blue optical bands. The supernova brightness fades at 0.78+/-0.02 mag/100d in the 4.5 micron channel which is close to that observed in the near-IR. We argue that the difference is a result of doubly ionized iron-peak elements dominating the bluer IRAC band while singly ionized species are controlling the longer wavelength channel. To test this, we use Large Binocular Telescope spectra taken during the same phases to show that doubly ionized emission lines do fade more slowly than their singly ionized cousins. We also find that [Co III] emission fades at more than twice the radioactive decay rate due to the combination of decreasing excitation in the nebula, recombination and cobalt decaying to iron. The nebular emission velocities of [Fe III] and [Co III] lines show a smaller blue-shift than emission from singly ionized atoms. The Si II velocity gradient near maximum light combined with our nebular velocity measurements suggest SN 2011fe was a typical member of the `low velocity gradient' class of type Ia. Analyzing IRAC photometry from other supernovae we find that mid-IR color of type Ia events is correlated with the early light curve width and can be used as an indicator of the radioactive nickel yield.

preprint2013arXiv

Type Ia Single Degenerate Survivors Must Be Overluminous

In the single-degenerate (SD) channel of a Type Ia supernovae (SN Ia) explosion, a main-sequence (MS) donor star survives the explosion but it is stripped of mass and shock heated. An essentially unavoidable consequence of mass loss during the explosion is that the companion must have an overextended envelope after the explosion. While this has been noted previously, it has not been strongly emphasized as an inevitable consequence. We calculate the future evolution of the companion by injecting 2-6 10^47 ergs into the stellar evolution model of a 1 Msun donor star based on the post-explosion progenitors seen in simulations. We find that, due to the Kelvin-Helmholtz collapse of the envelope, the companion must become significantly more luminous (10 - 10^3 Lsun) for a long period of time (10^3 - 10^4 years). The lack of such a luminous "leftover" star in the LMC supernova remnant SNR 0609-67.5 provides another piece of evidence against the SD scenario. We also show that none of the stars proposed as the survivors of the Tycho supernova, including Tycho G, could plausibly be the donor star. Additionally, luminous donors closer than ~10 Mpc should be observable with the Hubble Space Telescope starting ~2 years post-peak. Such systems include SN 1937C, SN 1972E, SN 1986G, and SN 2011fe. Thus, the SD channel is already ruled out for at least two nearby SNe Ia and can be easily tested for a number of additional ones. We also discuss similar implications for the companions of core-collapse SNe.

preprint2012arXiv

Very Low-mass Stellar and Substellar Companions to Solar-like Stars from MARVELS II: A Short-period Companion Orbiting an F Star with Evidence of a Stellar Tertiary And Significant Mutual Inclination

We report the discovery via radial velocity of a short-period (P = 2.430420 \pm 0.000006 days) companion to the F-type main sequence star TYC 2930-00872-1. A long-term trend in the radial velocities indicates the presence of a tertiary stellar companion with $P > 2000$ days. High-resolution spectroscopy of the host star yields T_eff = 6427 +/- 33 K, log(g) = 4.52 +/- 0.14, and [Fe/H]=-0.04 +/- 0.05. These parameters, combined with the broad-band spectral energy distribution and parallax, allow us to infer a mass and radius of the host star of M_1=1.21 +/- 0.08 M_\odot and R_1=1.09_{-0.13}^{+0.15} R_\odot. We are able to exclude transits of the inner companion with high confidence. The host star's spectrum exhibits clear Ca H and K core emission indicating stellar activity, but a lack of photometric variability and small v*sin(I) suggest the primary's spin axis is oriented in a pole-on configuration. The rotational period of the primary from an activity-rotation relation matches the orbital period of the inner companion to within 1.5 σ, suggesting they are tidally locked. If the inner companion's orbital angular momentum vector is aligned with the stellar spin axis, as expected through tidal evolution, then it has a stellar mass of M_2 ~ 0.3-0.4 M_\odot. Direct imaging limits the existence of stellar companions to projected separations < 30 AU. No set of spectral lines and no significant flux contribution to the spectral energy distribution from either companion are detected, which places individual upper mass limits of M < 1.0 M_\odot, provided they are not stellar remnants. If the tertiary is not a stellar remnant, then it likely has a mass of ~0.5-0.6 M_\odot, and its orbit is likely significantly inclined from that of the secondary, suggesting that the Kozai-Lidov mechanism may have driven the dynamical evolution of this system.

preprint2010arXiv

A New Cepheid Distance to the Giant Spiral M101 Based On Image Subtraction of HST/ACS Observations

We accurately determine a new Cepheid distance to M101 (NGC 5457) using archival HST/ACS V and I time series photometry of two fields within the galaxy. We make a slight modification to the ISIS image subtraction package to obtain optimal differential light curves from HST data. We discovered 827 Cepheids with periods between 3 and 80 days, the largest extragalactic sample of Cepheids observed with HST by a factor of 2. With this large Cepheid sample we find that the relative distance of M101 from the LMC is delta μLMC = 10.63 \pm 0.04 (random) \pm 0.07 (systematic) mag. If we use the geometrically determined maser distance to NGC 4258 as our distance anchor, the distance modulus of M101 is μ0 = 29.04 \pm 0.05 (random) \pm 0.18 (systematic) mag or D = 6.4 \pm 0.2 (random) \pm 0.5 (systematic) Mpc. The uncertainty is dominated by the maser distance estimate (\pm 0.15 mag), which should improve over the next few years. We determine a steep metallicity dependence, gamma, for our Cepheid sample through two methods, yielding gamma = -0.84 \pm 0.22 (random) \pm 0.07 (systematic) mag dex-1 and gamma = -0.72+0.20 -0.22 (random) \pm 0.06 (systematic) mag dex-1. We see marginal evidence for variations in the Wesenheit P-L relation slope as a function of deprojected galactocentric radius. We also use the TRGB method to independently determine the distance modulus to M101 of μ0 = 29.05 \pm 0.06 (random) \pm 0.12 (systematic) mag.