Showing 1–49 of 49 results for author: Hankins, M

An infrared census of R Coronae Borealis Stars II -- Spectroscopic classifications and implications for the rate of low-mass white dwarf mergers

Authors: Viraj R. Karambelkar, Mansi M. Kasliwal, Patrick Tisserand, Shreya Anand, Michael C. B. Ashley, Lars Bildsten, Geoffrey C. Clayton, Courtney C. Crawford, Kishalay De, Nicholas Earley, Matthew J. Hankins, Xander Hall, Astrid Lamberts, Ryan M. Lau, Dan McKenna, Anna Moore, Eran O. Ofek, Roger M. Smith, Roberto Soria, Jamie Soon, Tony Travouillon

Abstract: We present results from a systematic infrared (IR) census of R Coronae Borealis (RCB) stars in the Milky Way, using data from the Palomar Gattini IR (PGIR) survey. R Coronae Borealis stars are dusty, erratic variable stars presumably formed from the merger of a He-core and a CO-core white dwarf (WD). PGIR is a 30 cm $J$-band telescope with a 25 deg$^{2}$ camera that surveys 18000 deg$^{2}$ of the… ▽ More We present results from a systematic infrared (IR) census of R Coronae Borealis (RCB) stars in the Milky Way, using data from the Palomar Gattini IR (PGIR) survey. R Coronae Borealis stars are dusty, erratic variable stars presumably formed from the merger of a He-core and a CO-core white dwarf (WD). PGIR is a 30 cm $J$-band telescope with a 25 deg$^{2}$ camera that surveys 18000 deg$^{2}$ of the northern sky ($δ>-28^{o}$) at a cadence of 2 days. Using PGIR J-band lightcurves for $\sim$60 million stars together with mid-IR colors from WISE, we selected a sample of 530 candidate RCB stars. We obtained near-IR spectra for these candidates and identified 53 RCB stars in our sample. Accounting for our selection criteria, we find that there are a total of $\approx350^{+150}_{-100}$ RCB stars in the Milky Way. Assuming typical RCB lifetimes, this corresponds to an RCB formation rate of 0.8 - 5 $\times$ 10$^{-3}$ yr$^{-1}$, consistent with observational and theoretical estimates of the He-CO WD merger rate. We searched for quasi-periodic pulsations in the PGIR lightcurves of RCB stars and present pulsation periods for 16 RCB stars. We also examined high-cadenced TESS lightcurves for RCB and the chemically similar, but dustless hydrogen-deficient carbon (dLHdC) stars. We find that dLHdC stars show variations on timescales shorter than RCB stars, suggesting that they may have lower masses than RCB stars. Finally, we identified 3 new spectroscopically confirmed and 12 candidate Galactic DY Per type stars - believed to be colder cousins of RCB stars - doubling the sample of Galactic DY Per type stars. △ Less

Submitted 11 July, 2024; originally announced July 2024.

Comments: accepted for publication in PASP

SOFIA/FORCAST Galactic Center Source Catalog

Authors: Angela S. Cotera, Matthew J. Hankins, John Bally, Ashley T. Barnes, Cara D. Battersby, H Perry Hatchfield, Terry L. Herter, Ryan M. Lau, Steven N. Longmore, Elisabeth A. C. Mills, Mark R. Morris, James T. Radomski, Janet P. Simpson, Zachary Stephens, Daniel L. Walker

Abstract: The central regions of the Milky Way constitute a unique laboratory for a wide swath of astrophysical studies, consequently the inner $\sim$400 pc has been the target of numerous large surveys at all accessible wavelengths. In this paper we present a catalog of sources at 25 and 37 $μ$m located within all of the regions observed with the SOFIA/FORCAST instrument in the inner $\sim$200 pc of the Ga… ▽ More The central regions of the Milky Way constitute a unique laboratory for a wide swath of astrophysical studies, consequently the inner $\sim$400 pc has been the target of numerous large surveys at all accessible wavelengths. In this paper we present a catalog of sources at 25 and 37 $μ$m located within all of the regions observed with the SOFIA/FORCAST instrument in the inner $\sim$200 pc of the Galaxy. The majority of the observations were obtained as part of the SOFIA Cycle 7 Galactic Center Legacy program survey, which was designed to complement the Spitzer/MIPS 24 $μ$m catalog in regions saturated in the MIPS observations. Due to the wide variety of source types captured by our observations at 25 and 37 $μ$m, we do not limit the FORCAST source catalog to unresolved point sources, or treat all sources as if they are point-like sources. The catalog includes all detectable sources in the regions, resulting in a catalog of 950 sources, including point sources, compact sources, and extended sources. We also provide the user with metrics to discriminate between the source types. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: 29 pages, 13 figures, Accepted to ApJ

An Automated Catalog of Long Period Variables using Infrared Lightcurves from Palomar Gattini-IR

Authors: Aswin Suresh, Viraj Karambelkar, Mansi M. Kasliwal, Michael C. B. Ashley, Kishalay De, Matthew J. Hankins, Anna M. Moore, Jamie Soon, Roberto Soria, Tony Travouillon, Kayton K. Truong

Abstract: Stars in the Asymptotic Giant Branch (AGB) phase, dominated by low to intermediate-mass stars in the late stage of evolution, undergo periodic pulsations, with periods of several hundred days, earning them the name Long Period Variables (LPVs). These stars gradually shed their mass through stellar winds and mass ejections, enveloping themselves in dust. Infrared (IR) surveys can probe these dust-e… ▽ More Stars in the Asymptotic Giant Branch (AGB) phase, dominated by low to intermediate-mass stars in the late stage of evolution, undergo periodic pulsations, with periods of several hundred days, earning them the name Long Period Variables (LPVs). These stars gradually shed their mass through stellar winds and mass ejections, enveloping themselves in dust. Infrared (IR) surveys can probe these dust-enshrouded phases and uncover populations of LPV stars in the Milky Way. In this paper, we present a catalog of 159,696 Long Period Variables using near-IR lightcurves from the Palomar Gattini - IR (PGIR) survey. PGIR has been surveying the entire accessible northern sky ($δ> -28^{\circ}$) in the J-band at a cadence of 2-3 days since September 2018, and has produced J-band lightcurves for more than 60 million sources. We used a gradient-boosted decision tree classifier trained on a comprehensive feature set extracted from PGIR lightcurves to search for LPVs in this dataset. We developed a parallelized and optimized code to extract features at a rate of ~0.1 seconds per lightcurve. Our model can successfully distinguish LPVs from other stars with a true positive rate and weighted g-mean of 0.95. 73,346 (~46%) of the sources in our catalog are new, previously unknown LPVs. △ Less

Submitted 12 February, 2024; originally announced February 2024.

Comments: 15 pages, 15 figures

SOFIA/HAWC+ Far-Infrared Polarimetric Large Area CMZ Exploration (FIREPLACE) II: Detection of a Magnetized Dust Ring in the Galactic Center

Authors: Natalie O. Butterfield, Jordan A. Guerra, David T. Chuss, Mark R. Morris, Dylan Pare, Edward J. Wollack, Allison H. Costa, Matthew J. Hankins, Johannes Staguhn, Ellen Zweibel

Abstract: We present the detection of a magnetized dust ring (M0.8-0.2) in the Central Molecular Zone (CMZ) of the Galactic Center. The results presented in this paper utilize the first data release (DR1) of the Far-Infrared Polarimetric Large Area CMZ Exploration (FIREPLACE) survey (i.e., FIREPLACE I; Butterfield et al. 2023). The FIREPLACE survey is a 214 $μ$m polarimetic survey of the Galactic Center usi… ▽ More We present the detection of a magnetized dust ring (M0.8-0.2) in the Central Molecular Zone (CMZ) of the Galactic Center. The results presented in this paper utilize the first data release (DR1) of the Far-Infrared Polarimetric Large Area CMZ Exploration (FIREPLACE) survey (i.e., FIREPLACE I; Butterfield et al. 2023). The FIREPLACE survey is a 214 $μ$m polarimetic survey of the Galactic Center using the SOFIA/HAWC+ telescope. The M0.8-0.2 ring is a region of gas and dust that has a circular morphology with a central depression. The dust polarization in the M0.8-0.2 ring implies a curved magnetic field that traces the ring-like structure of the cloud. We posit an interpretation in which an expanding shell compresses and concentrates the ambient gas and magnetic field. We argue that this compression results in the strengthening of the magnetic field, as we infer from the observations toward the interior of the ring. △ Less

Submitted 29 April, 2024; v1 submitted 3 January, 2024; originally announced January 2024.

Comments: Accepted for publication in ApJ

A First Look with JWST Aperture Masking Interferometry (AMI): Resolving Circumstellar Dust around the Wolf-Rayet Binary WR 137 beyond the Rayleigh Limit

Authors: Ryan M. Lau, Matthew J. Hankins, Joel Sanchez-Bermudez, Deepashri Thatte, Anthony Soulain, Rachel A. Cooper, Anand Sivaramakrishnan, Michael F. Corcoran, Alexandra Z. Greenbaum, Theodore R. Gull, Yinuo Han, Olivia C. Jones, Thomas Madura, Anthony F. J. Moffat, Mark R. Morris, Takashi Onaka, Christopher M. P. Russell, Noel D. Richardson, Nathan Smith, Peter Tuthill, Kevin Volk, Gerd Weigelt, Peredur M. Williams

Abstract: We present infrared aperture masking interferometry (AMI) observations of newly formed dust from the colliding winds of the massive binary system Wolf-Rayet (WR) 137 with JWST using the Near Infrared Imager and Slitless Spectrograph (NIRISS). NIRISS AMI observations of WR 137 and a point-spread-function calibrator star, HD~228337, were taken using the F380M and F480M filters in 2022 July and Augus… ▽ More We present infrared aperture masking interferometry (AMI) observations of newly formed dust from the colliding winds of the massive binary system Wolf-Rayet (WR) 137 with JWST using the Near Infrared Imager and Slitless Spectrograph (NIRISS). NIRISS AMI observations of WR 137 and a point-spread-function calibrator star, HD~228337, were taken using the F380M and F480M filters in 2022 July and August as part of the Director's Discretionary Early Release Science (DD-ERS) program 1349. Interferometric observables (squared visibilities and closure phases) from the WR 137 "interferogram" were extracted and calibrated using three independent software tools: ImPlaneIA, AMICAL, and SAMpip. The analysis of the calibrated observables yielded consistent values except for slightly discrepant closure phases measured by ImPlaneIA. Based on all three sets of calibrated observables, images were reconstructed using three independent software tools: BSMEM, IRBis, and SQUEEZE. All reconstructed image combinations generated consistent images in both F380M and F480M filters. The reconstructed images of WR 137 reveal a bright central core with a $\sim300$ mas linear filament extending to the northwest. A geometric colliding-wind model with dust production constrained to the orbital plane of the binary system and enhanced as the system approaches periapsis provided a general agreement with the interferometric observables and reconstructed images. Based on a colliding-wind dust condensation analysis, we suggest that dust formation within the orbital plane of WR 137 is induced by enhanced equatorial mass-loss from the rapidly rotating O9 companion star, whose axis of rotation is aligned with that of the orbit. △ Less

Submitted 22 December, 2023; v1 submitted 27 November, 2023; originally announced November 2023.

Comments: 18 pages, 8 figures, Accepted for publication in ApJ. Updated plotting error in Fig. 2

The JWST Galactic Center Survey -- A White Paper

Authors: Rainer Schoedel, Steve Longmore, Jonny Henshaw, Adam Ginsburg, John Bally, Anja Feldmeier, Matt Hosek, Francisco Nogueras Lara, Anna Ciurlo, Mélanie Chevance, J. M. Diederik Kruijssen, Ralf Klessen, Gabriele Ponti, Pau Amaro-Seoane, Konstantina Anastasopoulou, Jay Anderson, Maria Arias, Ashley T. Barnes, Cara Battersby, Giuseppe Bono, Lucía Bravo Ferres, Aaron Bryant, Miguel Cano Gonzáalez, Santi Cassisi, Leonardo Chaves-Velasquez , et al. (85 additional authors not shown)

Abstract: The inner hundred parsecs of the Milky Way hosts the nearest supermassive black hole, largest reservoir of dense gas, greatest stellar density, hundreds of massive main and post main sequence stars, and the highest volume density of supernovae in the Galaxy. As the nearest environment in which it is possible to simultaneously observe many of the extreme processes shaping the Universe, it is one of… ▽ More The inner hundred parsecs of the Milky Way hosts the nearest supermassive black hole, largest reservoir of dense gas, greatest stellar density, hundreds of massive main and post main sequence stars, and the highest volume density of supernovae in the Galaxy. As the nearest environment in which it is possible to simultaneously observe many of the extreme processes shaping the Universe, it is one of the most well-studied regions in astrophysics. Due to its proximity, we can study the center of our Galaxy on scales down to a few hundred AU, a hundred times better than in similar Local Group galaxies and thousands of times better than in the nearest active galaxies. The Galactic Center (GC) is therefore of outstanding astrophysical interest. However, in spite of intense observational work over the past decades, there are still fundamental things unknown about the GC. JWST has the unique capability to provide us with the necessary, game-changing data. In this White Paper, we advocate for a JWST NIRCam survey that aims at solving central questions, that we have identified as a community: i) the 3D structure and kinematics of gas and stars; ii) ancient star formation and its relation with the overall history of the Milky Way, as well as recent star formation and its implications for the overall energetics of our galaxy's nucleus; and iii) the (non-)universality of star formation and the stellar initial mass function. We advocate for a large-area, multi-epoch, multi-wavelength NIRCam survey of the inner 100\,pc of the Galaxy in the form of a Treasury GO JWST Large Program that is open to the community. We describe how this survey will derive the physical and kinematic properties of ~10,000,000 stars, how this will solve the key unknowns and provide a valuable resource for the community with long-lasting legacy value. △ Less

Submitted 14 March, 2024; v1 submitted 18 October, 2023; originally announced October 2023.

Comments: This White Paper will be updated when required (e.g. new authors joining, editing of content). Most recent update: 24 Oct 2023

Modeling the Multiwavelength Evolution of the V960 Mon System

Authors: Adolfo S. Carvalho, Lynne A. Hillenbrand, Franz-Josef Hambsch, Shawn Dvorak, Michael Sitko, Ray W. Russell, Victoria Hammond, Michael Connelley, Michael C. B. Ashley, Matthew J. Hankins

Abstract: We study the evolution of the FU Ori object V960 Mon since its outburst, using available multi-wavelength photometric time series over 8 years, complemented by several epochs of moderate-dispersion spectrophotometry. We find that the source fading can be well-described by a decrease in the temperature of the inner disk, which results from a combination of decreasing accretion rate and increasing i… ▽ More We study the evolution of the FU Ori object V960 Mon since its outburst, using available multi-wavelength photometric time series over 8 years, complemented by several epochs of moderate-dispersion spectrophotometry. We find that the source fading can be well-described by a decrease in the temperature of the inner disk, which results from a combination of decreasing accretion rate and increasing inner disk radius. We model the system with a disk atmosphere model that produces the observed variations in multi-band photometry (this paper) and high resolution spectral lines (a companion paper). △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: 15 pages, 13 figures, 2 tables, Accepted to ApJ

SOFIA/HAWC+ Far-InfraRed Polarimetric Large Area CMZ Exploration (FIREPLACE) Survey I: General Results from the Pilot Program

Authors: Natalie O. Butterfield, David T. Chuss, Jordan A. Guerra, Mark R. Morris, Dylan Pare, Edward J. Wollack, C. Darren Dowell, Matthew J. Hankins, Kaitlyn Karpovich, Javad Siah, Johannes Staguhn, Ellen Zweibel

Abstract: We present the first data release (DR1) of the Far-Infrared Polarimetric Large Area CMZ Exploration (FIREPLACE) survey. The survey was taken using the 214-micron band of the HAWC+ instrument with the SOFIA telescope (19.6$'$ resolution; 0.7 pc). In this first data release we present dust polarization observations covering a ~0.5$°$ region of the Galactic Center's Central Molecular Zone (CMZ), appr… ▽ More We present the first data release (DR1) of the Far-Infrared Polarimetric Large Area CMZ Exploration (FIREPLACE) survey. The survey was taken using the 214-micron band of the HAWC+ instrument with the SOFIA telescope (19.6$'$ resolution; 0.7 pc). In this first data release we present dust polarization observations covering a ~0.5$°$ region of the Galactic Center's Central Molecular Zone (CMZ), approximately centered on the Sgr B2 complex. We detect ~25,000 Nyquist-sampled polarization pseudovectors, after applying the standard SOFIA cuts for minimum signal-to-noise in fractional polarization and total intensity of 3 and 200, respectively. Analysis of the magnetic field orientation suggests a bimodal distribution in the field direction. This bimodal distribution shows enhancements in the distribution of field directions for orientations parallel and perpendicular to the Galactic plane, which is suggestive of a CMZ magnetic field configuration with polodial and torodial components. Furthermore, a detailed analysis of individual clouds included in our survey (i.e., Sgr B2, Sgr B2-NW, Sgr B2-Halo, Sgr B1, and Clouds-E/F) shows these clouds have fractional polarization values of 1--10% at 214-micron, with most of the emission having values $<$5%. A few of these clouds (i.e., Sgr B2, Clouds-E/F) show relatively low fractional polarization values toward the cores of the cloud, with higher fractional polarization values toward the less dense periphery. We also observe higher fractional polarization towards compact HII regions which could indicate an enhancement in the grain alignment in the dust surrounding these sources. △ Less

Submitted 4 December, 2023; v1 submitted 2 June, 2023; originally announced June 2023.

Comments: Accepted for publication in ApJ

From Dust to Nanodust: Resolving Circumstellar Dust from the Colliding-Wind Binary Wolf-Rayet (WR) 140

Authors: Ryan M. Lau, Jason Wang, Matthew J. Hankins, Thayne Currie, Vincent Deo, Izumi Endo, Olivier Guyon, Yinuo Han, Anthony P. Jones, Nemanja Jovanovic, Julien Lozi, Anthony F. J. Moffat, Takashi Onaka, Garreth Ruane, Andreas A. C. Sander, Samaporn Tinyanont, Peter G. Tuthill, Gerd Weigelt, Peredur M. Williams, Sebastien Vievard

Abstract: Wolf-Rayet (WR) 140 is the archetypal periodic dust-forming colliding-wind binary that hosts a carbon-rich WR (WC) star and an O-star companion with an orbital period of 7.93 years and an orbital eccentricity of 0.9. Throughout the past several decades, multiple dust-formation episodes from WR 140 have been observed that are linked to the binary orbit and occur near the time of periastron passage.… ▽ More Wolf-Rayet (WR) 140 is the archetypal periodic dust-forming colliding-wind binary that hosts a carbon-rich WR (WC) star and an O-star companion with an orbital period of 7.93 years and an orbital eccentricity of 0.9. Throughout the past several decades, multiple dust-formation episodes from WR 140 have been observed that are linked to the binary orbit and occur near the time of periastron passage. Given its predictable dust-formation episodes, WR 140 presents an ideal astrophysical laboratory for investigating the formation and evolution of dust in the hostile environment around a massive binary system. In this paper, we present near- and mid-infrared (IR) spectroscopic and imaging observations of WR 140 with Subaru/SCExAO+CHARIS, Keck/NIRC2+PyWFS, and Subaru/COMICS taken between 2020 June and Sept that resolve the circumstellar dust emission linked to its most recent dust-formation episode in 2016 Dec. Our spectral energy distribution (SED) analysis of WR 140's resolved circumstellar dust emission reveals the presence of a hot ($T_\mathrm{d}\sim1000$ K) near-IR dust component that is co-spatial with the previously known and cooler ($T_\mathrm{d}\sim500$ K) mid-IR dust component composed of $300-500$ Å-sized dust grains. We attribute the hot near-IR dust emission to the presence of nano-sized ("nanodust") grains and suggest they were formed from grain-grain collisions or the rotational disruption of the larger grain size population by radiative torques in the strong radiation field from the central binary. Lastly, we speculate on the astrophysical implications of nanodust formation around colliding-wind WC binaries, which may present an early source of carbonaceous nanodust in the interstellar medium. △ Less

Submitted 23 May, 2023; originally announced May 2023.

Comments: 21 pages, 8 figures, Accepted for publication in ApJ

A survey for variable young stars with small telescopes: VI -- Analysis of the outbursting Be stars NSW284, Gaia19eyy, and VES263

Authors: Dirk Froebrich, Lynne A. Hillenbrand, Carys Herbert, Kishalay De, Jochen Eislöffel, Justyn Campbell-White, Ruhee Kahar, Franz-Josef Hambsch, Thomas Urtly, Adam Popowicz, Krzysztof Bernacki, Andrzej Malcher, Slawomir Lasota, Jerzy Fiolka, Piotr Jozwik-Wabik, Franky Dubois, Ludwig Logie, Steve Rau, Mark Phillips, George Fleming, Rafael Gonzalez Farfán, Francisco C. Soldán Alfaro, Tim Nelson, Stephen R. L. Futcher, Samantha M. Rolfe , et al. (22 additional authors not shown)

Abstract: This paper is one in a series reporting results from small telescope observations of variable young stars. Here, we study the repeating outbursts of three likely Be stars based on long-term optical, near-infrared, and mid-infrared photometry for all three objects, along with follow-up spectra for two of the three. The sources are characterised as rare, truly regularly outbursting Be stars. We inte… ▽ More This paper is one in a series reporting results from small telescope observations of variable young stars. Here, we study the repeating outbursts of three likely Be stars based on long-term optical, near-infrared, and mid-infrared photometry for all three objects, along with follow-up spectra for two of the three. The sources are characterised as rare, truly regularly outbursting Be stars. We interpret the photometric data within a framework for modelling light curve morphology, and find that the models correctly predict the burst shapes, including their larger amplitudes and later peaks towards longer wavelengths. We are thus able to infer the start and end times of mass loading into the circumstellar disks of these stars. The disk sizes are typically 3-6 times the areas of the central star. The disk temperatures are ~40%, and the disk luminosities are ~10% of those of the central Be star, respectively. The available spectroscopy is consistent with inside-out evolution of the disk. Higher excitation lines have larger velocity widths in their double-horned shaped emission profiles. Our observations and analysis support the decretion disk model for outbursting Be stars. △ Less

Submitted 6 February, 2023; originally announced February 2023.

Comments: Accepted for publication by MNRAS. 20 pages, 11 figures, 2 tables

A Systematic Study of Ia-CSM Supernovae from the ZTF Bright Transient Survey

Authors: Yashvi Sharma, Jesper Sollerman, Christoffer Fremling, Shrinivas R. Kulkarni, Kishalay De, Ido Irani, Steve Schulze, Nora Linn Strotjohann, Avishay Gal-Yam, Kate Maguire, Daniel A. Perley, Eric C. Bellm, Erik C. Kool, Thomas Brink, Rachel Bruch, Maxime Deckers, Richard Dekany, Alison Dugas, Samantha Goldwasser, Matthew J. Graham, Melissa L. Graham, Steven L. Groom, Matt Hankins, Jacob Jencson, Joel P. Johansson , et al. (13 additional authors not shown)

Abstract: Among the supernovae (SNe) that show strong interaction with the circumstellar medium, there is a rare subclass of Type Ia supernovae, SNe Ia-CSM, that show strong narrow hydrogen emission lines much like SNe IIn but on top of a diluted over-luminous Type Ia spectrum. In the only previous systematic study of this class (Silverman et al. 2013), 16 objects were identified, 8 historic and 8 from the… ▽ More Among the supernovae (SNe) that show strong interaction with the circumstellar medium, there is a rare subclass of Type Ia supernovae, SNe Ia-CSM, that show strong narrow hydrogen emission lines much like SNe IIn but on top of a diluted over-luminous Type Ia spectrum. In the only previous systematic study of this class (Silverman et al. 2013), 16 objects were identified, 8 historic and 8 from the Palomar Transient Factory (PTF). Now using the successor survey to PTF, the Zwicky Transient Facility (ZTF), we have classified 12 additional objects of this type through the systematic Bright Transient Survey (BTS). In this study, we present and analyze the optical and mid-IR light curves, optical spectra, and host galaxy properties of this sample. Consistent with previous studies, we find the objects to have slowly evolving light curves compared to normal SNe Ia with peak absolute magnitudes between -19.1 and -21, spectra having weak H$β$, large Balmer decrements of ~7 and strong Ca NIR emission. Out of 10 SNe from our sample observed by NEOWISE, 9 have $3σ$ detections, along with some showing a clear reduction in red-wing of H$α$, indicative of newly formed dust. We do not find our SN Ia-CSM sample to have a significantly different distribution of equivalent width of He I $\lambda5876$ than SNe IIn as observed in Silverman et al. 2013. The hosts tend to be late-type galaxies with recent star formation. We also derive a rate estimate of 29$^{+27}_{-21}$ Gpc$^{-3}$ yr$^{-1}$ for SNe Ia-CSM which is ~0.02--0.2 % of the SN Ia rate. This work nearly doubles the sample of well-studied Ia-CSM objects in Silverman et al. 2013, increasing the total number to 28. △ Less

Submitted 11 January, 2023; originally announced January 2023.

Comments: 26 pages, 14 figures, 5 tables. Submitted to ApJ

A search for relativistic ejecta in a sample of ZTF broad-lined Type Ic supernovae

Authors: Alessandra Corsi, Anna Y. Q. Ho, S. Bradley Cenko, Shrinivas R. Kulkarni, Shreya Anand, Sheng Yang, Jesper Sollerman, Gokul P. Srinivasaragavan, Conor M. B. Omand, Arvind Balasubramanian, Dale A. Frail, Christoffer Fremling, Daniel A. Perley, Yuhan Yao, Aishwarya S. Dahiwale, Kishalay De, Alison Dugas, Matthew Hankins, Jacob Jencson, Mansi M. Kasliwal, Anastasios Tzanidakis, Eric C. Bellm, Russ R. Laher, Frank J. Masci, Josiah N. Purdum , et al. (1 additional authors not shown)

Abstract: The dividing line between gamma-ray bursts (GRBs) and ordinary stripped-envelope core-collapse supernovae (SNe) is yet to be fully understood. Observationally mapping the variety of ejecta outcomes (ultra-relativistic, mildly-relativistic or non-relativistic) in SNe of Type Ic with broad lines (Ic-BL) can provide a key test to stellar explosion models. However, this requires large samples of the r… ▽ More The dividing line between gamma-ray bursts (GRBs) and ordinary stripped-envelope core-collapse supernovae (SNe) is yet to be fully understood. Observationally mapping the variety of ejecta outcomes (ultra-relativistic, mildly-relativistic or non-relativistic) in SNe of Type Ic with broad lines (Ic-BL) can provide a key test to stellar explosion models. However, this requires large samples of the rare Ic-BL events with follow-up observations in the radio, where fast ejecta can be probed largely free of geometry and viewing angle effects. Here, we present the results of a radio (and X-ray) follow-up campaign of 16 SNe Ic-BL detected by the Zwicky Transient Facility (ZTF). Our radio campaign resulted in 4 counterpart detections and 12 deep upper limits. None of the events in our sample is as relativistic as SN 1998bw and we constrain the fraction of SN 1998bw-like explosions to $< 19\%$ (3$σ$ Gaussian equivalent), a factor of $\approx 2$ smaller than previously established. We exclude relativistic ejecta with radio luminosity densities in between $\approx 5\times10^{27}$ erg s$^{-1}$ Hz$^{-1}$ and $\approx 10^{29}$ erg s$^{-1}$ Hz$^{-1}$ at $t\gtrsim 20$ d since explosion for $\approx 60\%$ of the events in our sample. This shows that SNe Ic-BL similar to the GRB-associated SN 1998bw, SN 2003lw, SN 2010dh, or to the relativistic SN 2009bb and iPTF17cw, are rare. Our results also exclude an association of the SNe Ic-BL in our sample with largely off-axis GRBs with energies $E\gtrsim 10^{50}$ erg. The parameter space of SN2006aj-like events (faint and fast-peaking radio emission) is, on the other hand, left largely unconstrained and systematically exploring it represents a promising line of future research. △ Less

Submitted 17 October, 2022; originally announced October 2022.

Comments: 28 pages, 11 figures, submitted to AAS journal

Nested Dust Shells around the Wolf-Rayet Binary WR 140 observed with JWST

Authors: Ryan M. Lau, Matthew J. Hankins, Yinuo Han, Ioannis Argyriou, Michael F. Corcoran, Jan J. Eldridge, Izumi Endo, Ori D. Fox, Macarena Garcia Marin, Theodore R. Gull, Olivia C. Jones, Kenji Hamaguchi, Astrid Lamberts, David R. Law, Thomas Madura, Sergey V. Marchenko, Hideo Matsuhara, Anthony F. J. Moffat, Mark R. Morris, Patrick W. Morris, Takashi Onaka, Michael E. Ressler, Noel D. Richardson, Christopher M. P. Russell, Joel Sanchez-Bermudez , et al. (7 additional authors not shown)

Abstract: Massive colliding-wind binaries that host a Wolf-Rayet (WR) star present a potentially important source of dust and chemical enrichment in the interstellar medium (ISM). However, the chemical composition and survival of dust formed from such systems is not well understood. The carbon-rich WR (WC) binary WR~140 presents an ideal astrophysical laboratory for investigating these questions given its w… ▽ More Massive colliding-wind binaries that host a Wolf-Rayet (WR) star present a potentially important source of dust and chemical enrichment in the interstellar medium (ISM). However, the chemical composition and survival of dust formed from such systems is not well understood. The carbon-rich WR (WC) binary WR~140 presents an ideal astrophysical laboratory for investigating these questions given its well-defined orbital period and predictable dust-formation episodes every 7.93 years around periastron passage. We present observations from our Early Release Science program (ERS1349) with the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) Medium-Resolution Spectrometer (MRS) and Imager that reveal the spectral and spatial signatures of nested circumstellar dust shells around WR~140. MIRI MRS spectroscopy of the second dust shell and Imager detections of over 17 shells formed throughout the past $\gtrsim130$ years confirm the survival of carbonaceous dust grains from WR~140 that are likely carriers of "unidentified infrared" (UIR)-band features at 6.4 and 7.7 $μ$m. The observations indicate that dust-forming WC binaries can enrich the ISM with organic compounds and carbonaceous dust. △ Less

Submitted 12 October, 2022; originally announced October 2022.

Comments: Published in Nature Astronomy on Oct 12, 2022; 21 pages, 5 figures, 2 tables

Journal ref: Lau, R.M., Hankins, M.J., Han, Y. et al. Nested dust shells around the Wolf-Rayet binary WR 140 observed with JWST. Nat Astron (2022)

SRGA J181414.6-225604: A new Galactic symbiotic X-ray binary outburst triggered by an intense mass loss episode of a heavily obscured Mira variable

Authors: Kishalay De, Ilya Mereminskiy, Roberto Soria, Charlie Conroy, Erin Kara, Shreya Anand, Michael C. B. Ashley, Martha L. Boyer, Deepto Chakrabarty, Brian Grefenstette, Matthew J. Hankins, Lynne A. Hillenbrand, Jacob E. Jencson, Viraj Karambelkar, Mansi M. Kasliwal, Ryan M. Lau, Alexander Lutovinov, Anna M. Moore, Mason Ng, Christos Panagiotou, Dheeraj R. Pasham, Andrey Semena, Robert Simcoe, Jamie Soon, Gokul P. Srinivasaragavan , et al. (2 additional authors not shown)

Abstract: We present the discovery and multi-wavelength characterization of SRGA J181414.6-225604, a Galactic hard X-ray transient discovered during the ongoing SRG/ART-XC sky survey. Using data from the Palomar Gattini-IR survey, we identify a spatially and temporally coincident variable infrared (IR) source, IRAS 18111-2257, and classify it as a very late-type (M7-M8), long period ($1502 \pm 24$ days) and… ▽ More We present the discovery and multi-wavelength characterization of SRGA J181414.6-225604, a Galactic hard X-ray transient discovered during the ongoing SRG/ART-XC sky survey. Using data from the Palomar Gattini-IR survey, we identify a spatially and temporally coincident variable infrared (IR) source, IRAS 18111-2257, and classify it as a very late-type (M7-M8), long period ($1502 \pm 24$ days) and luminous ($M_K\approx -9.9 \pm 0.2$) O-rich Mira donor star located at a distance of $\approx 14.6^{+2.9}_{-2.3}$ kpc. Combining multi-color photometric data over the last $\approx 25$ years, we show that the IR counterpart underwent a recent (starting $\approx 800$ days before the X-ray flare) enhanced mass loss (reaching $\approx 2.1 \times 10^{-5}$ M$_\odot$ yr$^{-1}$) episode resulting in an expanding dust shell obscuring the underlying star. Multi-epoch follow-up from Swift, NICER and NuSTAR reveal a $\approx 200$ day long X-ray outburst reaching a peak luminosity of $L_X \approx 2.5 \times 10^{36}$ erg s$^{-1}$, characterized by a heavily absorbed ($N_{\rm H} \approx 6\times 10^{22}$ cm$^{-2}$) X-ray spectrum consistent with an optically thick Comptonized plasma. The X-ray spectral and timing behavior suggest the presence of clumpy wind accretion together with a dense ionized nebula overabundant in silicate material surrounding the compact object. Together, we show that SRGA J181414.6-225604 is a new symbiotic X-ray binary in outburst, triggered by an intense dust formation episode of a highly evolved donor. Our results offer the first direct confirmation for the speculated connection between enhanced late-stage donor mass loss and active lifetimes of the symbiotic X-ray binaries. △ Less

Submitted 18 May, 2022; originally announced May 2022.

Comments: 32 pages, 15 figures, Accepted for publication in ApJ

The Hydrogen-Poor Superluminous Supernovae from the Zwicky Transient Facility Phase-I Survey: I. Light Curves and Measurements

Authors: Z. H. Chen, Lin Yan, T. Kangas, R. Lunnan, S. Schulze, J. Sollerman, D. A. Perley, T. -W. Chen, K. Taggart, K. R. Hinds, A. Gal-Yam, X. F. Wang, I. Andreoni, E. Bellm, J. S. Bloom, K. Burdge, A. Burgos, D. Cook, A. Dahiwale, K. De, R. Dekany, A. Dugas, S. Frederik, C. Fremling, M. Graham , et al. (18 additional authors not shown)

Abstract: During the Zwicky Transient Facility (ZTF) Phase-I operation, 78 hydrogen-poor superluminous supernovae (SLSNe-I) were discovered in less than three years, making up the largest sample from a single survey. This paper (Paper I) presents the data, including the optical/ultraviolet light curves and classification spectra, while Paper II in this series will focus on the detailed analysis of the light… ▽ More During the Zwicky Transient Facility (ZTF) Phase-I operation, 78 hydrogen-poor superluminous supernovae (SLSNe-I) were discovered in less than three years, making up the largest sample from a single survey. This paper (Paper I) presents the data, including the optical/ultraviolet light curves and classification spectra, while Paper II in this series will focus on the detailed analysis of the light curves and modeling. Our photometry is primarily taken by the ZTF in the $g,r,i$ bands, and with additional data from other ground-based facilities and Swift. The events of our sample cover a redshift range of $z = 0.06 - 0.67$, with a median and $1σ$ error (16\% and 84\% percentiles) $z_{\rm med} = 0.265^{+0.143}_{-0.135}$. The peak luminosity covers $-22.8\,{\rm mag} \leq M_{g,\rm peak} \leq -19.8$\,mag, with a median value of $-21.48^{+1.13}_{-0.61}$\,mag. Their light curves evolve slowly with the mean rest-frame rise time of $t_{\rm rise} = 41.9\pm17.8$\,days. The luminosity and time scale distributions suggest that low luminosity SLSNe-I with peak luminosity $\sim -20$\,mag or extremely fast rising events ($<10$\,days) exist but are rare. We confirm previous findings that slowly rising SLSNe-I also tend to fade slowly. The rest-frame color and temperature evolution show large scatters, suggesting that the SLSN-I population may have diverse spectral energy distributions. The peak rest-frame color shows a moderate correlation with the peak absolute magnitude, i.e. brighter SLSNe-I tend to have bluer colors. With optical and ultraviolet photometry, we construct bolometric luminosity and derive a bolometric correction relation generally applicable for converting $g,r$-band photometry to bolometric luminosity for SLSNe-I. △ Less

Submitted 3 November, 2022; v1 submitted 4 February, 2022; originally announced February 2022.

Comments: 38 pages, 25 figures, Accepted by APJ

A massive AGB donor in Scutum X-1: Identification of the first Mira variable in an X-ray binary

Authors: Kishalay De, Deepto Chakrabarty, Roberto Soria, Michael C. B. Ashley, Charlie Conroy, Matthew J. Hankins, Mansi M. Kasliwal, Ryan M. Lau, Anna M. Moore, Robert Simcoe, Jamie Soon, Tony Travouillon

Abstract: The symbiotic X-ray binary Sct X-1 was suggested as the first known neutron star accreting from a red supergiant companion. Although known for nearly 50 years, detailed characterization of the donor remains lacking, particularly due to the extremely high reddening towards the source ($A_V\gtrsim25$ mag). Here, we present i) improved localization of the counterpart using Gaia and Chandra observatio… ▽ More The symbiotic X-ray binary Sct X-1 was suggested as the first known neutron star accreting from a red supergiant companion. Although known for nearly 50 years, detailed characterization of the donor remains lacking, particularly due to the extremely high reddening towards the source ($A_V\gtrsim25$ mag). Here, we present i) improved localization of the counterpart using Gaia and Chandra observations, ii) the first broadband infrared spectrum ($\approx1-5\,μ$m; $R\approx 2000$) obtained with SpeX on the NASA Infrared Telescope Facility and iii) $J$-band light curve from the Palomar Gattini-IR survey. The infrared spectrum is characterized by i) deep water absorption features (H$_2$O index $\approx 40$%), ii) strong TiO, VO and CO features, and iii) weak/absent CN lines. We show that these features are inconsistent with known red supergiants, but suggest a M8-9 III type O-rich Mira donor star. We report the discovery of large amplitude ($ΔJ\approx3.5$ mag) periodic photometric variability suggesting a pulsation period of $621\pm36\,{\rm(systematic)}\pm8\,{\rm(statistical)}$ days, which we use to constrain the donor to be a relatively luminous Mira ($M_K=-8.6\pm0.3$ mag) at a distance of $3.6^{+0.8}_{-0.7}$ kpc. Comparing these characteristics to recent models, we find the donor to be consistent with a $\approx 3-5$ M$_\odot$ star at an age of $\approx 0.1-0.3$ Gyr. Together, we show that Sct X-1 was previously mis-classified as an evolved High Mass X-ray Binary; instead it is an intermediate mass system with the first confirmed Mira donor in an X-ray binary. We discuss the implications of Mira donors in symbiotic X-ray binaries, and highlight the potential of wide-field infrared time domain surveys and broadband infrared spectroscopy to unveil their demographics. △ Less

Submitted 2 March, 2022; v1 submitted 24 January, 2022; originally announced January 2022.

Comments: 10 pages, 4 figures. Accepted for publication in ApJ Letters

Supernova Siblings and their Parent Galaxies in the Zwicky Transient Facility Bright Transient Surve

Authors: M. L. Graham, C. Fremling, D. A. Perley, R. Biswas, C. A. Phillips, J. Sollerman, P. E. Nugent, S. Nance, S. Dhawan, J. Nordin, A. Goobar, A. Miller, J. D. Neill, X. J. Hall, M. J. Hankins, D. A. Duev, M. M. Kasliwal, M. Rigault, E. C. Bellm, D. Hale, P. Mróz, S. R. Kulkarni

Abstract: Supernova (SN) siblings -- two or more SNe in the same parent galaxy -- are useful tools for exploring progenitor stellar populations as well as properties of the host galaxies such as distance, star formation rate, dust extinction, and metallicity. Since the average SN rate for a Milky Way-type galaxy is just one per century, a large imaging survey is required to discover an appreciable sample of… ▽ More Supernova (SN) siblings -- two or more SNe in the same parent galaxy -- are useful tools for exploring progenitor stellar populations as well as properties of the host galaxies such as distance, star formation rate, dust extinction, and metallicity. Since the average SN rate for a Milky Way-type galaxy is just one per century, a large imaging survey is required to discover an appreciable sample of SN siblings. From the wide-field Zwicky Transient Facility (ZTF) Bright Transient Survey (BTS; which aims for spectroscopic completeness for all transients which peak brighter than $r{<}$18.5 mag) we present 10 SN siblings in 5 parent galaxies. For each of these families we analyze the SN's location within the host and its underlying stellar population, finding agreement with expectations that SNe from more massive progenitors are found nearer to their host core and in regions of more active star formation. We also present an analysis of the relative rates of core collapse and thermonuclear SN siblings, finding a significantly lower ratio than past SN sibling samples due to the unbiased nature of the ZTF. △ Less

Submitted 29 December, 2021; originally announced December 2021.

Comments: 14 pages, 4 tables, 7 figures, accepted to MNRAS

Second Timescale Photometry of the Very Fast Nova V1674 Her with Palomar Gattini-IR

Authors: Kylie Y. Hansen, Kishalay De, Michael C. B. Ashley, Mansi M. Kasliwal, Alexander Delacroix, Tim Greffe, David Hale, Matthew J. Hankins, Ryan Lau, Chengkui Li, Daniel McKenna, Anna M. Moore, Eran O. Ofek, Roger M. Smith, Jamie Soon, Roberto Soria, Gokul P. Srinivasaragavan, Tony Travouillon

Abstract: We report second-timescale infrared photometry of the nova V1674 Her using Palomar Gattini-IR. These observations constitute the first infrared and highest temporal resolution data (resolution of ~ 0.84 s) of the nova reported to date. PGIR observed in this fast readout mode for more than an hour on three nights between 3 and 6 days after discovery. We searched for periodic variability using a Lom… ▽ More We report second-timescale infrared photometry of the nova V1674 Her using Palomar Gattini-IR. These observations constitute the first infrared and highest temporal resolution data (resolution of ~ 0.84 s) of the nova reported to date. PGIR observed in this fast readout mode for more than an hour on three nights between 3 and 6 days after discovery. We searched for periodic variability using a Lomb-Scargle periodogram and did not detect anything down to a three sigma upper limit of 0.093 mag. This suggests that the periodic variability detected in the nova by Patterson et al. 2021 was lower by at least a factor of about 1.65 in the first week of the eruption. △ Less

Submitted 25 October, 2021; originally announced October 2021.

Comments: 3 pages, 1 figure

Progenitor and Close-In Circumstellar Medium of Type II Supernova 2020fqv from High-Cadence Photometry and Ultra-Rapid UV Spectroscopy

Authors: Samaporn Tinyanont, Ryan Ridden-Harper, Ryan Foley, Viktoriya Morozova, Charles Kilpatrick, Georgios Dimitriadis, Lindsay DeMarchi, Alexander Gagliano, Wynn V. Jacobson-Galán, Alexander Messick, Justin Pierel, Anthony Piro, Enrico Ramirez-Ruiz, Matthew Siebert, Kenneth Chambers, Karoli Clever, David Coulter, Kishalay De, Matthew Hankins, Tiara Hung, Saurabh Jha, Camilo Jiménez-Ángel, David Jones, Mansi Kasliwal, Chien-Cheng Lin , et al. (13 additional authors not shown)

Abstract: We present observations of SN 2020fqv, a Virgo-cluster Type II core-collapse supernova (CCSN) with a high temporal resolution light curve from the Transiting Exoplanet Survey Satellite (TESS) covering the time of explosion; ultraviolet (UV) spectroscopy from the Hubble Space Telescope (HST) starting 3.3 days post-explosion; ground-based spectroscopic observations starting 1.1~days post-explosion;… ▽ More We present observations of SN 2020fqv, a Virgo-cluster Type II core-collapse supernova (CCSN) with a high temporal resolution light curve from the Transiting Exoplanet Survey Satellite (TESS) covering the time of explosion; ultraviolet (UV) spectroscopy from the Hubble Space Telescope (HST) starting 3.3 days post-explosion; ground-based spectroscopic observations starting 1.1~days post-explosion; along with extensive photometric observations. Massive stars have complicated mass-loss histories leading up to their death as CCSNe, creating circumstellar medium (CSM) with which the SNe interact. Observations during the first few days post-explosion can provide important information about the mass-loss rate during the late stages of stellar evolution. Model fits to the quasi-bolometric light curve of SN 2020fqv reveal ~0.23 $M_{\odot}$ of CSM confined within ~1450 $R_{\odot}$ ($10^{14}$ cm) from its progenitor star. Early spectra (<4 days post-explosion), both from HST and ground-based observatories, show emission features from high-ionization metal species from the outer, optically thin part of this CSM. We find that the CSM is consistent with an eruption caused by the injection of $\sim$$5\times 10^{46}$ erg into the stellar envelope $\sim$300 days pre-explosion, potentially from a nuclear burning instability at the onset of oxygen burning. Light-curve fitting, nebular spectroscopy, and pre-explosion \textit{HST} imaging consistently point to a red supergiant (RSG) progenitor with $M_{\rm ZAMS}$$\approx$$13.5$--$15 \, M_{\odot}$, typical for SN~II progenitor stars. This finding demonstrates that a typical RSG, like the progenitor of SN 2020fqv, has a complicated mass-loss history immediately before core collapse. △ Less

Submitted 20 October, 2021; originally announced October 2021.

Comments: Published in MNRAS

PGIR 20eid (SN2020qmp): A Type IIP Supernova at 15.6 Mpc discovered by the Palomar Gattini-IR survey

Authors: G. P. Srinivasaragavan, I. Sfaradi, J. Jencson, K. De, A. Horesh, M. M. Kasliwal, S. Tinyanont, M. Hankins, S. Schulze, M. C. B. Ashley, M. J. Graham, V. Karambelkar, R. Lau, A. A. Mahabal, A. M. Moore, E. O. Ofek, Y. Sharma, J. Sollerman, J. Soon, R. Soria, T. Travouillon, R. Walters

Abstract: We present a detailed analysis of SN 2020qmp, a nearby type IIP core-collapse supernova (CCSN), discovered by the Palomar Gattini-IR (PGIR) survey in the galaxy UGC07125. We illustrate how the multiwavelength study of this event helps our general understanding of stellar progenitors and circumstellar medium (CSM) interactions in CCSNe. We also highlight the importance of near-infrared (NIR) survey… ▽ More We present a detailed analysis of SN 2020qmp, a nearby type IIP core-collapse supernova (CCSN), discovered by the Palomar Gattini-IR (PGIR) survey in the galaxy UGC07125. We illustrate how the multiwavelength study of this event helps our general understanding of stellar progenitors and circumstellar medium (CSM) interactions in CCSNe. We also highlight the importance of near-infrared (NIR) surveys for early detections of SNe in dusty environments. SN 2020qmp displays characteristic hydrogen lines in its optical spectra, as well as a plateau in its optical LC, hallmarks of a type IIP SN. We do not detect linear polarization during the plateau phase, with a 3 sigma upper limit of 0.78%. Through hydrodynamical LC modeling and an analysis of its nebular spectra, we estimate a progenitor mass of around 11 solar masses, and an explosion energy of around 0.8e51 erg. We find that the spectral energy distribution cannot be explained by a simple CSM interaction model, assuming a constant shock velocity and steady mass-loss rate, and the excess X-ray luminosity compared with the synchrotron radio luminosity suggests deviations from equipartition. Finally, we demonstrate the advantages of NIR surveys over optical surveys for the detection of dust-obscured CCSNe in the local universe. Specifically, our simulations show that the Wide-Field Infrared Transient Explorer will detect about 14 more CCSNe out of 75 expected in its footprint within 40 Mpc, over five years than an optical survey equivalent to the Zwicky Transient Facility would detect. We have determined or constrained the main properties of SN 2020qmp and of its progenitor, highlighting the value of multiwavelength follow-up observations of nearby CCSNe. We have also shown that forthcoming NIR surveys will finally enable us to do a nearly complete census of CCSNe in the local universe. △ Less

Submitted 2 March, 2022; v1 submitted 5 September, 2021; originally announced September 2021.

Comments: 15 pages, 14 figures, Accepted in Astronomy & Astrophysics; Revised version after referee comments and language edits from journal

Journal ref: A&A 660, A138 (2022)

Discovery of a 310-day Period from the Enshrouded Massive System NaSt1 (WR 122)

Authors: Ryan M. Lau, Samaporn Tinyanont, Matthew J. Hankins, Michael C. B. Ashley, Kishalay De, Alexei V. Filippenko, Lynne A. Hillenbrand, Mansi M. Kasliwal, Jon C. Mauerhan, Anthony F. J. Moffat, Anna M. Moore, Nathan Smith, Jamie Soon, Roberto Soria, Tony Travouillon, Karel A. van der Hucht, Peredur M. Williams, WeiKang Zheng

Abstract: We present optical and infrared (IR) light curves of NaSt1, also known as Wolf-Rayet (WR) 122, with observations from Palomar Gattini-IR (PGIR), the Zwicky Transient Facility (ZTF), the Katzman Automatic Imaging Telescope (KAIT), the Asteroid Terrestrial-impact Last Alert System (ATLAS), and the All-Sky Automated Survey for Supernovae (ASAS-SN). We identify a $P=309.7\pm0.7$ d photometric period f… ▽ More We present optical and infrared (IR) light curves of NaSt1, also known as Wolf-Rayet (WR) 122, with observations from Palomar Gattini-IR (PGIR), the Zwicky Transient Facility (ZTF), the Katzman Automatic Imaging Telescope (KAIT), the Asteroid Terrestrial-impact Last Alert System (ATLAS), and the All-Sky Automated Survey for Supernovae (ASAS-SN). We identify a $P=309.7\pm0.7$ d photometric period from the optical and IR light curves that reveal periodic, sinusoidal variability between 2014 July and 2021 July. We also present historical IR light curves taken between 1983 July and 1989 May, which show variability consistent with the period of the present-day light curves. In the past, NaSt1 was brighter in the $J$ band with larger variability amplitudes than the present-day PGIR values, suggesting that NaSt1 exhibits variability on longer ($\gtrsim$ decade) timescales. Sinusoidal fits to the recent optical and IR light curves show that the amplitude of NaSt1's variability differs at various wavelengths and also reveal significant phase offsets of $17.0\pm2.5$ d between the ZTF $r$ and PGIR $J$ light curves. We interpret the $310$ d photometric period from NaSt1 as the orbital period of an enshrouded massive binary. We suggest that the photometric variability of NaSt1 may arise from variations in the line-of-sight optical depth toward circumstellar optical/IR emitting regions throughout its orbit due to colliding-wind dust formation. We speculate that past mass transfer in NaSt1 may have been triggered by Roche-lobe overflow (RLOF) during an eruptive phase of an Ofpe/WN9 star. Lastly, we argue that NaSt1 is no longer undergoing RLOF mass transfer. △ Less

Submitted 29 August, 2021; v1 submitted 15 March, 2021; originally announced March 2021.

Comments: 20 Pages, 7 Figures, 5 Tables, 1 Table Attachment, Accepted to ApJ

Infrared Spectropolarimetric Detection of Intrinsic Polarization from a Core-Collapse Supernova

Authors: Samaporn Tinyanont, Maxwell Millar-Blanchaer, Mansi Kasliwal, Dimitri Mawet, Douglas C Leonard, Mattia Bulla, Kishalay De, Nemanja Jovanovic, Matthew Hankins, Gautam Vasisht, Eugene Serabyn

Abstract: Massive stars die an explosive death as a core-collapse supernova (CCSN). The exact physical processes that cause the collapsing star to rebound into an explosion are not well-understood, and the key in resolving this issue may lie in the measurement of the shape of CCSNe ejecta. Spectropolarimetry is the only way to perform this measurement for CCSNe outside of the Milky Way and Magellanic Clouds… ▽ More Massive stars die an explosive death as a core-collapse supernova (CCSN). The exact physical processes that cause the collapsing star to rebound into an explosion are not well-understood, and the key in resolving this issue may lie in the measurement of the shape of CCSNe ejecta. Spectropolarimetry is the only way to perform this measurement for CCSNe outside of the Milky Way and Magellanic Clouds. We present an infrared (IR) spectropolarimetric detection of a CCSN, enabled by the new highly sensitive WIRC+Pol instrument at Palomar Observatory, that can observe CCSNe (M = -17 mags) out to 20 Mpc to ~0.1% polarimetric precision. IR spectropolarimetry is less affected than optical by dust scattering in the circumstellar and interstellar media, thereby providing a more unbiased probe of the intrinsic geometry of the SN ejecta. SN 2018hna, a SN 1987A-like explosion, shows 2.0+-0.3% continuum polarization in the J band oriented at ~160 degree on-sky at 182 d after the explosion. Assuming prolate geometry like in SN 1987A, we infer an ejecta axis ratio of <0.48 with the axis of symmetry pointing at 70 degree position angle. The axis ratio is similar to that of SN 1987A suggesting that they may share intrinsic geometry and inclination angle. Our data do not rule out oblate ejecta. We also observe one other core-collapse and two thermonuclear SNe in the J band. SN 2020oi, a stripped-envelope Type Ic SN in Messier 100 has p = 0.37+-0.09% at peak light, indicative of either a 10% asymmetry or host interstellar polarization. The SNe Ia, 2019ein and 2020ue have p < 0.33% and < 1.08% near peak light, indicative of asymmetries of less than 10% and 20%, respectively. △ Less

Submitted 3 February, 2021; originally announced February 2021.

Outbursting Young Stellar Object PGIR 20dci in the Perseus Arm

Authors: Lynne A. Hillenbrand, Kishalay De, Matthew Hankins, Mansi M. Kasliwal, Luisa M. Rebull, Ryan M. Lau, Roc M. Cutri, Michael C. B. Ashley, Viraj R. Karambelkar, Anna M. Moore, T. Travouillon, A. K. Mainzer

Abstract: We report the discovery of a likely outbursting Class I young stellar object, associated with the star-forming region NGC 281-W (distance $\sim 2.8$ kpc). The source is currently seen only at infrared wavelengths, appearing in both the Palomar Gattini InfraRed ($1.2~μ$m) and the Near Earth Object Widefield Infrared Survey Explorer ($3.4$ and $4.6~μ$m) photometric time-domain surveys. Recent near-i… ▽ More We report the discovery of a likely outbursting Class I young stellar object, associated with the star-forming region NGC 281-W (distance $\sim 2.8$ kpc). The source is currently seen only at infrared wavelengths, appearing in both the Palomar Gattini InfraRed ($1.2~μ$m) and the Near Earth Object Widefield Infrared Survey Explorer ($3.4$ and $4.6~μ$m) photometric time-domain surveys. Recent near-infrared imaging reveals a new, extended scattered light nebula. Recent near-infrared spectroscopy confirms the similarity of PGIR 20dci to FU Ori type sources, based on strong molecular absorption in CO, H$_2$O, and OH, weak absorption in several atomic lines, and a warm wind/outflow as indicated by a P Cygni profile in the HeI 10830 A line. This is a rare case of an FU Ori star with a well-measured long term photometric rise before a sharper outburst, and the second instance of an FU Ori star with a documented two-step brightening in the mid-infrared. △ Less

Submitted 9 February, 2021; v1 submitted 11 January, 2021; originally announced January 2021.

Comments: accepted to AAS Journals on 2021, February 5

A population of heavily reddened, optically missed novae from Palomar Gattini-IR: Constraints on the Galactic nova rate

Authors: Kishalay De, Mansi M. Kasliwal, Matthew J. Hankins, Jennifer L. Sokoloski, Scott M. Adams, Michael C. B. Ashley, Aliya-Nur Babul, Ashot Bagdasaryan, Alexandre Delacroix, Richard Dekany, Timothee Greffe, David Hale, Jacob E. Jencson, Viraj R. Karambelkar, Ryan M. Lau, Ashish Mahabal, Daniel McKenna, Anna M. Moore, Eran O. Ofek, Manasi Sharma, Roger M. Smith, Jamie Soon, Roberto Soria, Gokul Srinivasaragavan, Samaporn Tinyanont , et al. (3 additional authors not shown)

Abstract: The nova rate in the Milky Way remains largely uncertain, despite its vital importance in constraining models of Galactic chemical evolution as well as understanding progenitor channels for Type Ia supernovae. The rate has been previously estimated in the range of $\approx10-300$ yr$^{-1}$, either based on extrapolations from a handful of very bright optical novae or the nova rates in nearby galax… ▽ More The nova rate in the Milky Way remains largely uncertain, despite its vital importance in constraining models of Galactic chemical evolution as well as understanding progenitor channels for Type Ia supernovae. The rate has been previously estimated in the range of $\approx10-300$ yr$^{-1}$, either based on extrapolations from a handful of very bright optical novae or the nova rates in nearby galaxies; both methods are subject to debatable assumptions. The total discovery rate of optical novae remains much smaller ($\approx5-10$ yr$^{-1}$) than these estimates, even with the advent of all-sky optical time domain surveys. Here, we present a systematic sample of 12 spectroscopically confirmed Galactic novae detected in the first 17 months of Palomar Gattini-IR (PGIR), a wide-field near-infrared time domain survey. Operating in $J$-band ($\approx1.2$ $μ$m) that is relatively immune to dust extinction, the extinction distribution of the PGIR sample is highly skewed to large extinction values ($> 50$% of events obscured by $A_V\gtrsim5$ mag). Using recent estimates for the distribution of mass and dust in the Galaxy, we show that the observed extinction distribution of the PGIR sample is commensurate with that expected from dust models. The PGIR extinction distribution is inconsistent with that reported in previous optical searches (null hypothesis probability $< 0.01$%), suggesting that a large population of highly obscured novae have been systematically missed in previous optical searches. We perform the first quantitative simulation of a $3π$ time domain survey to estimate the Galactic nova rate using PGIR, and derive a rate of $\approx 46.0^{+12.5}_{-12.4}$ yr$^{-1}$. Our results suggest that all-sky near-infrared time-domain surveys are well poised to uncover the Galactic nova population. △ Less

Submitted 24 January, 2021; v1 submitted 11 January, 2021; originally announced January 2021.

Comments: 25 pages, 12 figures. Submitted to ApJ

Census of R Coronae Borealis stars I: Infrared light curves from Palomar Gattini IR

Authors: Viraj R. Karambelkar, Mansi M. Kasliwal, Patrick Tisserand, Kishalay De, Shreya Anand, Michael C. B. Ashley, Alex Delacroix, Matthew Hankins, Jacob E. Jencson, Ryan M. Lau, Dan McKenna, Anna Moore, Eran O. Ofek, Roger M. Smith, Roberto Soria, Jamie Soon, Samaporn Tinyanont, Tony Travouillon, Yuhan Yao

Abstract: We are undertaking the first systematic infrared (IR) census of R Coronae Borealis (RCB) stars in the Milky Way, beginning with IR light curves from the Palomar Gattini IR (PGIR) survey. PGIR is a 30 cm $J$-band telescope with a 25 deg$^{2}$ camera that is surveying 18000 deg$^{2}$ of the northern sky ($δ>-28^{o}$) at a cadence of 2 days. We present PGIR light curves for 922 RCB candidates selecte… ▽ More We are undertaking the first systematic infrared (IR) census of R Coronae Borealis (RCB) stars in the Milky Way, beginning with IR light curves from the Palomar Gattini IR (PGIR) survey. PGIR is a 30 cm $J$-band telescope with a 25 deg$^{2}$ camera that is surveying 18000 deg$^{2}$ of the northern sky ($δ>-28^{o}$) at a cadence of 2 days. We present PGIR light curves for 922 RCB candidates selected from a mid-IR color-based catalog (Tisserand et al. 2020). Of these 922, 149 are promising RCB candidates as they show pulsations or declines similar to RCB stars. Majority of the candidates that are not RCB stars are either long period variables (LPVs) or RV-Tauri stars. We identify IR color-based criteria to better distinguish between RCB stars and LPVs. As part of a pilot spectroscopic run, we obtained NIR spectra for 26 out of the 149 promising candidates and spectroscopically confirm 11 new RCB stars. We detect strong He I $λ10830$ features in spectra of all RCB stars, likely originating within high velocity (200-400 km-s$^{-1}$) winds in their atmospheres. 9 of these RCB stars show $^{12}$C$^{16}$O and $^{12}$C$^{18}$O molecular absorption features, suggesting that they are formed through a white dwarf merger. We detect quasi-periodic pulsations in the light curves of 5 RCB stars. The periods range between 30-125 days and likely originate from the strange-mode instability in these stars. Our pilot run results motivate a dedicated IR spectroscopic campaign to classify all RCB candidates. △ Less

Submitted 21 December, 2020; originally announced December 2020.

Comments: 20 pages, 14 figures, Submitted to ApJ

Revealing Efficient Dust-Formation at Low Metallicity in Extragalactic Carbon-Rich Wolf-Rayet Binaries

Authors: Ryan M. Lau, Matthew J. Hankins, Mansi M. Kasliwal, Howard E. Bond, Kishalay De, Jacob E. Jencson, Anthony F. J. Moffat, Nathan Smith, Peredur M. Williams

Abstract: We present Spitzer/IRAC observations of dust formation from six extragalactic carbon-rich Wolf-Rayet (WC) binary candidates in low-metallicity (Z $\lesssim0.65$ Z$_\odot$) environments using multi-epoch mid-infrared (IR) imaging data from the SPitzer InfraRed Intensive Transients Survey (SPIRITS). Optical follow-up spectroscopy of SPIRITS~16ln, 19q, 16df, 18hb, and 14apu reveals emission features… ▽ More We present Spitzer/IRAC observations of dust formation from six extragalactic carbon-rich Wolf-Rayet (WC) binary candidates in low-metallicity (Z $\lesssim0.65$ Z$_\odot$) environments using multi-epoch mid-infrared (IR) imaging data from the SPitzer InfraRed Intensive Transients Survey (SPIRITS). Optical follow-up spectroscopy of SPIRITS~16ln, 19q, 16df, 18hb, and 14apu reveals emission features from C IV $\lambda5801\text{-}12$~and/or the C III-IV $\lambda4650$ He II $\lambda4686$~blend that are consistent with early-type WC stars. We identify SPIRITS~16ln as the variable mid-IR counterpart of the recently discovered colliding-wind WC4+O binary candidate, N604-WRXc, located in the sub-solar metallicity NGC 604 H II~region in M33. We interpret the mid-IR variability from SPIRITS~16ln as a dust-formation episode in an eccentric colliding-wind WC binary. SPIRITS~19q, 16df, 14apu, and 18hb exhibit absolute [3.6] magnitudes exceeding one of most IR-luminous dust-forming WC systems known, WR~104 (M$_\mathrm{[3.6]}\lesssim-12.3$). An analysis of dust formation in the mid-IR outburst from SPIRITS~19q reveals a high dust production rate of $\dot{M}_d\gtrsim2\times10^{-6}$ M$_\odot$ yr$^{-1}$, which may therefore exceed that of the most efficient dust-forming WC systems known. We demonstrate that efficient dust-formation is feasible from early-type WC binaries in the theoretical framework of colliding-wind binary dust formation if the systems host an O-type companion with a high mass-loss rate ($\dot{M}\gtrsim1.6\times10^{-6}$ M$_\odot$ yr$^{-1}$). This efficient dust-formation from early-type WC binaries highlights their potential role as significant sources of dust in low-metallicity environments. △ Less

Submitted 6 January, 2021; v1 submitted 19 November, 2020; originally announced November 2020.

Comments: 28 pages, 10 Figures, and 5 Tables. Accepted to ApJ on 2021 Jan 4 (Replacement of version submitted to arXiv on 2020 Nov 19)

Bright, months-long stellar outbursts announce the explosion of interaction-powered supernovae

Authors: Nora L. Strotjohann, Eran O. Ofek, Avishay Gal-Yam, Rachel Bruch, Steve Schulze, Nir Shaviv, Jesper Sollerman, Alexei V. Filippenko, Ofer Yaron, Christoffer Fremling, Jakob Nordin, Erik C. Kool, Dan A. Perley, Anna Y. Q. Ho, Yi Yang, Yuhan Yao, Maayane T. Soumagnac, Melissa L. Graham, Cristina Barbarino, Leonardo Tartaglia, Kishalay De, Daniel A. Goldstein, David O. Cook, Thomas G. Brink, Kirsty Taggart , et al. (31 additional authors not shown)

Abstract: Interaction-powered supernovae (SNe) explode within an optically-thick circumstellar medium (CSM) that could be ejected during eruptive events. To identify and characterize such pre-explosion outbursts we produce forced-photometry light curves for 196 interacting SNe, mostly of Type IIn, detected by the Zwicky Transient Facility between early 2018 and June 2020. Extensive tests demonstrate that we… ▽ More Interaction-powered supernovae (SNe) explode within an optically-thick circumstellar medium (CSM) that could be ejected during eruptive events. To identify and characterize such pre-explosion outbursts we produce forced-photometry light curves for 196 interacting SNe, mostly of Type IIn, detected by the Zwicky Transient Facility between early 2018 and June 2020. Extensive tests demonstrate that we only expect a few false detections among the 70,000 analyzed pre-explosion images after applying quality cuts and bias corrections. We detect precursor eruptions prior to 18 Type IIn SNe and prior to the Type Ibn SN2019uo. Precursors become brighter and more frequent in the last months before the SN and month-long outbursts brighter than magnitude -13 occur prior to 25% (5 - 69%, 95% confidence range) of all Type IIn SNe within the final three months before the explosion. With radiative energies of up to $10^{49}\,\text{erg}$, precursors could eject $\sim1\,\text{M}_\odot$ of material. Nevertheless, SNe with detected precursors are not significantly more luminous than other SNe IIn and the characteristic narrow hydrogen lines in their spectra typically originate from earlier, undetected mass-loss events. The long precursor durations require ongoing energy injection and they could, for example, be powered by interaction or by a continuum-driven wind. Instabilities during the neon and oxygen burning phases are predicted to launch precursors in the final years to months before the explosion; however, the brightest precursor is 100 times more energetic than anticipated. △ Less

Submitted 12 March, 2021; v1 submitted 21 October, 2020; originally announced October 2020.

Journal ref: ApJ 907 99 (2021)

Optical follow-up of the neutron star-black hole mergers S200105ae and S200115j

Authors: Shreya Anand, Michael W. Coughlin, Mansi M. Kasliwal, Mattia Bulla, Tomás Ahumada, Ana Sagués Carracedo, Mouza Almualla, Igor Andreoni, Robert Stein, Francois Foucart, Leo P. Singer, Jesper Sollerman, Eric C. Bellm, Bryce Bolin, M. D. Caballero-García, Alberto J. Castro-Tirado, S. Bradley Cenko, Kishalay De, Richard G. Dekany, Dmitry A. Duev, Michael Feeney, Christoffer Fremling, Daniel A. Goldstein, V. Zach Golkhou, Matthew J. Graham , et al. (24 additional authors not shown)

Abstract: LIGO and Virgo's third observing run (O3) revealed the first neutron star-black hole (NSBH) merger candidates in gravitational waves. These events are predicted to synthesize r-process elements creating optical/near-IR "kilonova" (KN) emission. The joint gravitational-wave (GW) and electromagnetic detection of an NSBH merger could be used to constrain the equation of state of dense nuclear matter,… ▽ More LIGO and Virgo's third observing run (O3) revealed the first neutron star-black hole (NSBH) merger candidates in gravitational waves. These events are predicted to synthesize r-process elements creating optical/near-IR "kilonova" (KN) emission. The joint gravitational-wave (GW) and electromagnetic detection of an NSBH merger could be used to constrain the equation of state of dense nuclear matter, and independently measure the local expansion rate of the universe. Here, we present the optical follow-up and analysis of two of the only three high-significance NSBH merger candidates detected to date, S200105ae and S200115j, with the Zwicky Transient Facility (ZTF). ZTF observed $\sim$\,48\% of S200105ae and $\sim$\,22\% of S200115j's localization probabilities, with observations sensitive to KNe brighter than $-$17.5\,mag fading at 0.5\,mag/day in g- and r-bands; extensive searches and systematic follow-up of candidates did not yield a viable counterpart. We present state-of-the-art KN models tailored to NSBH systems that place constraints on the ejecta properties of these NSBH mergers. We show that with depths of $\rm m_{\rm AB}\approx 22$ mag, attainable in meter-class, wide field-of-view survey instruments, strong constraints on ejecta mass are possible, with the potential to rule out low mass ratios, high BH spins, and large neutron star radii. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Characterization of Temporarily-Captured Minimoon 2020 CD$_3$ by Keck Time-resolved Spectrophotometry

Authors: Bryce T. Bolin, Christoffer Fremling, Timothy R. Holt, Matthew J. Hankins, Tomás Ahumada, Shreya Anand, Varun Bhalerao, Kevin B. Burdge, Chris M. Copperwheat, Michael Coughlin, Kunal P. Deshmukh, Kishalay De, Mansi M. Kasliwal, Alessandro Morbidelli, Josiah N. Purdum, Robert Quimby, Dennis Bodewits, Chan-Kao Chang, Wing-Huen Ip, Chen-Yen Hsu, Russ R. Laher, Zhong-Yi Lin, Carey M. Lisse, Frank J. Masci, Chow-Choong Ngeow , et al. (20 additional authors not shown)

Abstract: We present time-resolved visible spectrophotometry of minimoon 2020 CD$_3$, the second asteroid known to become temporarily captured by the Earth-Moon system's gravitational field. The spectrophotometry was taken with Keck I/LRIS between wavelengths 434 nm and 912 nm in $B$, $g$, $V$, $R$, $I$ and RG850 filters as it was leaving the Earth-Moon system on 2020 March 23 UTC. The spectrophotometry of… ▽ More We present time-resolved visible spectrophotometry of minimoon 2020 CD$_3$, the second asteroid known to become temporarily captured by the Earth-Moon system's gravitational field. The spectrophotometry was taken with Keck I/LRIS between wavelengths 434 nm and 912 nm in $B$, $g$, $V$, $R$, $I$ and RG850 filters as it was leaving the Earth-Moon system on 2020 March 23 UTC. The spectrophotometry of 2020 CD$_3$ most closely resembles the spectra of V-type asteroids and some Lunar rock samples with a reddish slope of ~18$\%$/100 nm between 434 nm and 761 nm corresponding to colors of $g$-$r$ = 0.62$\pm$0.08, $r$-$i$ = 0.21 $\pm$ 0.06 and an absorption band at ~900 nm corresponding to $i$-$z$ = -0.54$\pm$0.10. Combining our measured 31.9$\pm$0.1 absolute magnitude with an albedo of 0.35 typical for V-type asteroids, we determine 2020 CD$_3$'s diameter to be ~0.9$\pm$0.1 m making it the first minimoon and one of the smallest asteroids to be spectrally studied. We use our time-series photometry to detect periodic lightcurve variations with a $<$10$^{-4}$ false alarm probability corresponding to a lightcurve period of ~573 s and a lightcurve amplitude of ~1 mag implying 2020 CD$_3$ possesses a $b/a$ axial ratio of ~2.5. In addition, we extend the observational arc of 2020 CD$_3$ to 37 days between 2020 February 15 UTC and 2020 March 23 UTC. From the improved orbital solution for 2020 CD$_3$, we estimate its likely duration of its capture to be ~2 y, and we measure the non-gravitation perturbation on its orbit due to radiation pressure with an area-to-mass ratio of 6.9$\pm$2.4$\times$10$^{-4}$ m$^2$/kg implying a density of 2.3$\pm$0.8 g/cm$^3$, broadly compatible with the densities of other meter-scale asteroids and Lunar rock. We searched for pre-discovery detections of 2020 CD$_3$ in the ZTF archive as far back as 2018 October, but were unable to locate any positive detections. △ Less

Submitted 12 August, 2020; originally announced August 2020.

Comments: 26 pages, 8 figures, 3 Tables, accepted for publication in ApJL

Resolving Decades of Periodic Spirals from the Wolf-Rayet Dust Factory WR 112

Authors: Ryan M. Lau, Matthew J. Hankins, Yinuo Han, Izumi Endo, Anthony F. J. Moffat, Michael E. Ressler, Itsuki Sakon, Joel Sanchez-Bermudez, Anthony Soulain, Ian R. Stevens, Peter G. Tuthill, Peredur M. Williams

Abstract: WR 112 is a dust-forming carbon-rich Wolf-Rayet (WC) binary with a dusty circumstellar nebula that exhibits a complex asymmetric morphology, which traces the orbital motion and dust formation in the colliding winds of the central binary. Unraveling the complicated circumstellar dust emission around WR 112 therefore provides an opportunity to understand the dust formation process in colliding-wind… ▽ More WR 112 is a dust-forming carbon-rich Wolf-Rayet (WC) binary with a dusty circumstellar nebula that exhibits a complex asymmetric morphology, which traces the orbital motion and dust formation in the colliding winds of the central binary. Unraveling the complicated circumstellar dust emission around WR 112 therefore provides an opportunity to understand the dust formation process in colliding-wind WC binaries. In this work, we present a multi-epoch analysis of the circumstellar dust around WR 112 using seven high spatial resolution (FWHM $\sim0.3-0.4''$) N-band ($λ\sim12$ $μ$m) imaging observations spanning almost 20 years and includes newly obtained images from Subaru/COMICS in Oct 2019. In contrast to previous interpretations of a face-on spiral morphology, we observe clear evidence of proper motion of the circumstellar dust around WR 112 consistent with a nearly edge-on spiral with a $θ_s=55^\circ$ half-opening angle and a $\sim20$-yr period. The revised near edge-on geometry of WR 112 reconciles previous observations of highly variable non-thermal radio emission that was inconsistent with a face-on geometry. We estimate a revised distance to WR 112 of $d = 3.39^{+0.89}_{-0.84}$ kpc based on the observed dust expansion rate and a spectroscopically derived WC terminal wind velocity of $v_\infty= 1230\pm260$ km s$^{-1}$. With the newly derived WR 112 parameters we fit optically-thin dust spectral energy distribution models and determine a dust production rate of $\dot{M}_d=2.7^{+1.0}_{-1.3}\times10^{-6}$ M$_\odot$ yr$^{-1}$, which demonstrates that WR 112 is one of the most prolific dust-making WC systems known. △ Less

Submitted 3 August, 2020; originally announced August 2020.

Comments: 17 pages, 9 figures, 1 animated gif, accepted for publication in ApJ

Constraining the X-ray - Infrared spectral index of second-timescale flares from SGR1935+2154 with Palomar Gattini-IR

Authors: Kishalay De, Michael C. B. Ashley, Igor Andreoni, Mansi M. Kasliwal, Roberto Soria, Gokul P. Srinivasaragavan, Ce Cai, Alexander Delacroix, Tim Greffe, David Hale, Matthew J. Hankins, Chengkui Li, Daniel McKenna, Anna M. Moore, Eran O. Ofek, Roger M. Smith, Jamie Soon, Tony Travouillon, Shuangnan Zhang

Abstract: The Galactic magnetar SGR1935+2154 has been reported to produce the first known example of a bright millisecond duration radio burst (FRB 200428) similar to the cosmological population of fast radio bursts (FRBs), bolstering the association of FRBs to active magnetars. The detection of a coincident bright X-ray burst has revealed the first observed multi-wavelength counterpart of a FRB. However, t… ▽ More The Galactic magnetar SGR1935+2154 has been reported to produce the first known example of a bright millisecond duration radio burst (FRB 200428) similar to the cosmological population of fast radio bursts (FRBs), bolstering the association of FRBs to active magnetars. The detection of a coincident bright X-ray burst has revealed the first observed multi-wavelength counterpart of a FRB. However, the search for similar emission at optical wavelengths has been hampered by the high inferred extinction on the line of sight. Here, we present results from the first search for second-timescale emission from the source at near-infrared wavelengths using the Palomar Gattini-IR observing system in J-band, made possible by a recently implemented detector read-out mode that allowed for short exposure times of 0.84 s with 99.9% observing efficiency. With a total observing time of 12 hours (47728 images) on source, we place median $3\,σ$ limits on the second-timescale emission of $< 20$ mJy (13.1 AB mag). We present non-detection limits from epochs of four simultaneous X-ray bursts detected by the Insight-{\it HXMT} and {\it NuSTAR} telescopes during our observing campaign. The limits translate to an extinction corrected fluence limit of $< 125$ Jy ms for an estimated extinction of $A_J = 2.0$ mag. These limits provide the most stringent constraints to date on the fluence of flares at frequencies of $\sim 10^{14}$ Hz, and constrain the ratio of the near-infrared (NIR) fluence to that of coincident X-ray bursts to $R_{\rm NIR} < 2.5 \times 10^{-2}$. Our observations were sensitive enough to easily detect a near-infrared counterpart of FRB 200428 if the NIR emission falls on the same power law as that observed across its radio to X-ray spectrum. The non-detection of NIR emission around the coincident X-ray bursts constrains the fluence index of the brightest burst to be steeper than $0.35$. △ Less

Submitted 6 July, 2020; originally announced July 2020.

Comments: 10 pages, 4 figures, submitted to ApJL. Comments welcome

Helium-rich Superluminous Supernovae From the Zwicky Transient Facility

Authors: Lin Yan, D. Perley, S. Schulze, R. Lunnan, J. Sollerman, K. De, Z. Chen, C. Fremling, A. Gal-Yam, K. Taggart, T. W. Chen, I. Andreoni, E. C. Bellm, V. Cunningham, R. Dekany, D. Duev, C. Fransson, R. Laher, M. Hankins, A. Ho, J. Jencson, S. Kaye, S. Kulkarni, M. Kasliwal, V. Golkhou , et al. (17 additional authors not shown)

Abstract: Helium is expected to be present in the ejecta of some hydrogen-poor superluminous supernovae (SLSN-I). However, so far only one event, PTF10hgi has been identified with He features in its photospheric spectra (Quimby et al. 2018). We present the discovery of a new He-rich SLSN-I, ZTF19aawfbtg (SN2019hge) at $z=0.0866$. This event has more than 10 optical spectra at phases from $-41$ to $+103$\,da… ▽ More Helium is expected to be present in the ejecta of some hydrogen-poor superluminous supernovae (SLSN-I). However, so far only one event, PTF10hgi has been identified with He features in its photospheric spectra (Quimby et al. 2018). We present the discovery of a new He-rich SLSN-I, ZTF19aawfbtg (SN2019hge) at $z=0.0866$. This event has more than 10 optical spectra at phases from $-41$ to $+103$\,days relative to the peak, most of which match well with that of PTF10hgi. Confirmation comes from a near-IR spectrum taken at $+34$ days, revealing He I features with P-Cygni profiles at 1.083 and 2.058$μ$m. Using the optical spectra of PTF10hgi and SN2019hge as templates, we examine 70 SLSN-I discovered by ZTF in the first two years of operation and found additional five SLSN-I with distinct He-features. The excitation of He\,I atoms in normal core collapse supernovae requires non-thermal radiation, as proposed by previous studies. These He-rich events can not be explained by traditional $^{56}$Ni mixing model because of their blue spectra, high peak luminosity and slow rising time scales. Magnetar models offer a possible solution since pulsar winds naturally generate high energy particles as sources of non-thermal excitation. An alternative model is ejecta interaction with H-poor CSM which may be supported by the observed light curve undulations. These six SLSN-Ib appear to have relatively low-peak luminosities (rest-frame $M_g = -20.06\pm0.16$ mag). △ Less

Submitted 24 August, 2020; v1 submitted 24 June, 2020; originally announced June 2020.

Comments: 12 pages, 6 figures, Accepted for publication in ApJL. Matched with the published version

Kilonova Luminosity Function Constraints based on Zwicky Transient Facility Searches for 13 Neutron Star Mergers

Authors: Mansi M. Kasliwal, Shreya Anand, Tomas Ahumada, Robert Stein, Ana Sagues Carracedo, Igor Andreoni, Michael W. Coughlin, Leo P. Singer, Erik C. Kool, Kishalay De, Harsh Kumar, Mouza AlMualla, Yuhan Yao, Mattia Bulla, Dougal Dobie, Simeon Reusch, Daniel A. Perley, S. Bradley Cenko, Varun Bhalerao, David L. Kaplan, Jesper Sollerman, Ariel Goobar, Christopher M. Copperwheat, Eric C. Bellm, G. C. Anupama , et al. (78 additional authors not shown)

Abstract: We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GR… ▽ More We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged) △ Less

Submitted 19 June, 2020; originally announced June 2020.

Comments: Submitted to ApJ

Revisiting the Impact of Dust Production from Carbon-Rich Wolf-Rayet Binaries

Authors: Ryan M. Lau, J. J. Eldridge, Matthew J. Hankins, Astrid Lamberts, Itsuki Sakon, Peredur M. Williams

Abstract: We present a dust spectral energy distribution (SED) and binary stellar population analysis revisiting the dust production rates (DPRs) in the winds of carbon-rich Wolf-Rayet (WC) binaries and their impact on galactic dust budgets. DustEM SED models of 19 Galactic WC ``dustars" reveal DPRs of $\dot{M}_d\sim10^{-10}-10^{-6}$ M$_\odot$ yr$^{-1}$ and carbon dust condensation fractions, $χ_C$, between… ▽ More We present a dust spectral energy distribution (SED) and binary stellar population analysis revisiting the dust production rates (DPRs) in the winds of carbon-rich Wolf-Rayet (WC) binaries and their impact on galactic dust budgets. DustEM SED models of 19 Galactic WC ``dustars" reveal DPRs of $\dot{M}_d\sim10^{-10}-10^{-6}$ M$_\odot$ yr$^{-1}$ and carbon dust condensation fractions, $χ_C$, between $0.002 - 40\%$. A large ($0.1 - 1.0$ $μ$m) dust grain size composition is favored for efficient dustars where $χ_C\gtrsim1\%$. Results for dustars with known orbital periods verify a power-law relation between $χ_C$, orbital period, WC mass-loss rate, and wind velocity consistent with predictions from theoretical models of dust formation in colliding-wind binaries. We incorporated dust production into Binary Population and Spectral Synthesis (BPASS) models to analyze dust production rates from WC dustars, asymptotic giant branch stars (AGBs), red supergiants (RSGs), and core-collapse supernovae (SNe). BPASS models assuming constant star formation (SF) and a co-eval $10^6$ M$_\odot$ stellar population were performed at low, Large Magellanic Cloud (LMC)-like, and solar metallicities (Z = 0.001, 0.008, and 0.020). Both constant SF and co-eval models show that SNe are net dust destroyers at all metallicities. Constant SF models at LMC-like metallicities show that AGB stars slightly outproduce WC binaries and RSGs by factors of $2-3$, whereas at solar metallicites WC binaries are the dominant source of dust for $\sim60$ Myr until the onset of AGBs, which match the dust input of WC binaries. Co-eval population models show that for "bursty" SF, AGB stars dominate dust production at late times ($t\gtrsim 70$ Myr). △ Less

Submitted 15 June, 2020; originally announced June 2020.

Comments: 34 pages, 12 figures, 7 tables, Accepted for publication in ApJ

The Zwicky Transient Facility Census of the Local Universe I: Systematic search for Calcium rich gap transients reveal three related spectroscopic sub-classes

Authors: Kishalay De, Mansi M. Kasliwal, Anastasios Tzanidakis, U. Christoffer Fremling, Scott Adams, Igor Andreoni, Ashot Bagdasaryan, Eric C. Bellm, Lars Bildsten, Christopher Cannella, David O. Cook, Alexandre Delacroix, Andrew Drake, Dmitry Duev, Alison Dugas, Sara Frederick, Avishay Gal-Yam, Daniel Goldstein, V. Zach Golkhou, Matthew J. Graham, David Hale, Matthew Hankins, George Helou, Anna Y. Q. Ho, Ido Irani , et al. (25 additional authors not shown)

Abstract: (Abridged) Using the Zwicky Transient Facility alert stream, we are conducting a large campaign to spectroscopically classify all transients occurring in galaxies in the Census of the Local Universe (CLU) catalog. The aim of the experiment is to construct a spectroscopically complete, volume-limited sample of transients coincident within 100" of CLU galaxies out to 200 Mpc, and to a depth of 20 ma… ▽ More (Abridged) Using the Zwicky Transient Facility alert stream, we are conducting a large campaign to spectroscopically classify all transients occurring in galaxies in the Census of the Local Universe (CLU) catalog. The aim of the experiment is to construct a spectroscopically complete, volume-limited sample of transients coincident within 100" of CLU galaxies out to 200 Mpc, and to a depth of 20 mag. We describe the survey design and spectroscopic completeness from the first 16 months of operations. We present results from a systematic search for Calcium rich gap transients in the sample of 22 low luminosity (peak absolute magnitude $M > -17$), hydrogen poor events found in the experiment (out of 754 spectroscopically classified SNe). We report the detection of eight Calcium rich gap transients, and constrain their volumetric rate to be at least $\approx 15\pm5$% of the SN Ia rate. Combining this sample with ten events from the literature, we find a likely continuum of spectroscopic properties ranging from events with SN Ia-like features (Ca-Ia objects) to SN Ib/c-like features (Ca-Ib/c objects) at peak light. Within the Ca-Ib/c events, we find two populations of events distinguished by their red ($g - r \approx 1.5$ mag) or green ($g - r \approx 0.5$ mag) spectral colors at $r$-band peak, wherein redder events show strong line blanketing signatures, slower light curves, weaker He lines and lower [Ca II]/[O I] in the nebular phase. Together, we find that the spectroscopic continuum, volumetric rates and striking old environments are consistent with the explosive burning of He shells on low mass white dwarfs. We posit that Ca-Ia and red Ca-Ib/c objects are consistent with the double detonation of He shells with high He burning efficiency, while green Ca-Ib/c objects could arise from less efficient He burning scenarios such as detonations in low density He shells or He shell deflagrations. △ Less

Submitted 19 April, 2020; originally announced April 2020.

Comments: 30 figures, 6 tables. Submitted to ApJ

Early Ultra-Violet observations of type IIn supernovae constrain the asphericity of their circumstellar material

Authors: Maayane T. Soumagnac, Eran O. Ofek, Jingyi Liang, Avishay Gal-yam, Peter Nugent, Yi Yang, S. Bradley Cenko, Jesper Sollerman, Daniel A. Perley, Igor Andreoni, Cristina Barbarino, Kevin B. Burdge, Rachel J. Bruch, Kishalay De, Alison Dugas, Christoffer Fremling, Melissa L. Graham, Matthew J. Hankins, Nora Linn Strotjohann, Shane Moran, James D. Neill, Steve Schulze, David L. Shupe, Brigitta M. Sipocz, Kirsty Taggart , et al. (19 additional authors not shown)

Abstract: We present a survey of the early evolution of 12 Type IIn supernovae (SNe IIn) in the Ultra-Violet (UV) and visible light. We use this survey to constrain the geometry of the circumstellar material (CSM) surrounding SN IIn explosions, which may shed light on their progenitor diversity. In order to distinguish between aspherical and spherical circumstellar material (CSM), we estimate the blackbody… ▽ More We present a survey of the early evolution of 12 Type IIn supernovae (SNe IIn) in the Ultra-Violet (UV) and visible light. We use this survey to constrain the geometry of the circumstellar material (CSM) surrounding SN IIn explosions, which may shed light on their progenitor diversity. In order to distinguish between aspherical and spherical circumstellar material (CSM), we estimate the blackbody radius temporal evolution of the SNe IIn of our sample, following the method introduced by Soumagnac et al. We find that higher luminosity objects tend to show evidence for aspherical CSM. Depending on whether this correlation is due to physical reasons or to some selection bias, we derive a lower limit between 35% and 66% on the fraction of SNe IIn showing evidence for aspherical CSM. This result suggests that asphericity of the CSM surrounding SNe IIn is common - consistent with data from resolved images of stars undergoing considerable mass loss. It should be taken into account for more realistic modelling of these events. △ Less

Submitted 15 January, 2020; originally announced January 2020.

Comments: Submitted to ApJ. Comments are welcome

SOFIA/FORCAST Galactic Center Legacy Survey: Overview

Authors: Matthew J. Hankins, Ryan M. Lau, James T. Radomski, Angela S. Cotera, Mark R. Morris, Elisabeth A. C. Mills, Daniel L. Walker, Ashley T. Barnes, Janet P. Simpson, Terry L. Herter, Steven N. Longmore, John Bally, Mansi M. Kasliwal, Nadeen B. Sabha, Macarena Garcia-Marin

Abstract: The Galactic Center contains some of the most extreme conditions for star formation in our Galaxy as well as many other phenomena that are unique to this region. Given our relative proximity to the Galactic Center, we are able to study details of physical processes to a level that is simply not yet possible for more distant galaxies, yielding an otherwise inaccessible view of the nuclear region of… ▽ More The Galactic Center contains some of the most extreme conditions for star formation in our Galaxy as well as many other phenomena that are unique to this region. Given our relative proximity to the Galactic Center, we are able to study details of physical processes to a level that is simply not yet possible for more distant galaxies, yielding an otherwise inaccessible view of the nuclear region of a galaxy. We recently carried out a targeted imaging survey of mid-infrared bright portions of the Galactic Center at 25 and 37 $μ$m using the FORCAST instrument on SOFIA. This survey was one of the inaugural Legacy Programs from SOFIA cycle 7, observing a total area of 403 arcmin$^2$ (2180 pc$^2$), including the Sgr A, B, and C complexes. Here we present an overview of the survey strategy, observations, and data reduction as an accompaniment to the initial public release of the survey data. We discuss interesting regions and features within the data including extended features near the circumnuclear disk, structures in the Arched Filaments and Sickle H II regions, and signs of embedded star formation in Sgr B2 and Sgr C. We also feature a handful of less well studied mid-infrared sources located between Sgr A and Sgr C that could be sites of relatively isolated star formation activity. Last, we discuss plans for subsequent publications and future data releases from the survey. △ Less

Submitted 14 April, 2020; v1 submitted 15 January, 2020; originally announced January 2020.

Comments: 21 pages, 11 figures, Accepted to ApJ

GROWTH on S190814bv: Deep Synoptic Limits on the Optical/Near-Infrared Counterpart to a Neutron Star-Black Hole Merger

Authors: Igor Andreoni, Daniel A. Goldstein, Mansi M. Kasliwal, Peter E. Nugent, Rongpu Zhou, Jeffrey A. Newman, Mattia Bulla, Francois Foucart, Kenta Hotokezaka, Ehud Nakar, Samaya Nissanke, Geert Raaijmakers, Joshua S. Bloom, Kishalay De, Jacob E. Jencson, Charlotte Ward, Tomás Ahumada, Shreya Anand, David A. H. Buckley, Maria D. Caballero-García, Alberto J. Castro-Tirado, Christopher M. Copperwheat, Michael W. Coughlin, S. Bradley Cenko, Mariusz Gromadzki , et al. (27 additional authors not shown)

Abstract: On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community wi… ▽ More On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to a NSBH merger. To search for potential counterparts, the GROWTH collaboration performed real-time image subtraction on 6 nights of public Dark Energy Camera (DECam) images acquired in the three weeks following the merger, covering $>$98\% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and we placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be $M_\mathrm{ej} < 0.04$~$M_{\odot}$ at polar viewing angles, or $M_\mathrm{ej} < 0.03$~$M_{\odot}$ if the opacity is $κ< 2$~cm$^2$g$^{-1}$. Assuming a tidal deformability for the neutron star at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be $ χ< 0.7$ for mass ratios $Q < 6$, with weaker constraints for more compact neutron stars. We publicly release the photometry from this campaign at http://www.astro.caltech.edu/~danny/static/s190814bv. △ Less

Submitted 31 December, 2019; v1 submitted 29 October, 2019; originally announced October 2019.

Comments: Accepted to ApJ with no major changes

Palomar Gattini-IR: Survey overview, data processing system, on-sky performance and first results

Authors: Kishalay De, Matthew J. Hankins, Mansi M. Kasliwal, Anna M. Moore, Eran O. Ofek, Scott M. Adams, Michael C. B. Ashley, Aliya-Nur Babul, Ashot Bagdasaryan, Kevin B. Burdge, Jill Burnham, Richard G. Dekany, Alexander Declacroix, Antony Galla, Tim Greffe, David Hale, Jacob E. Jencson, Ryan M. Lau, Ashish Mahabal, Daniel McKenna, Manasi Sharma, Patrick L. Shopbell, Roger M. Smith, Jamie Soon, Jennifer Sokoloski , et al. (2 additional authors not shown)

Abstract: (Abridged) Palomar Gattini-IR is a new wide-field, near-infrared robotic time domain survey operating at Palomar Observatory. Using a 30 cm telescope mounted with a H2RG detector, Gattini-IR achieves a field of view of 25 sq. deg. with a pixel scale of 8.7" in J-band. Here, we describe the system design, survey operations, data processing system and on-sky performance of Palomar Gattini-IR. As a p… ▽ More (Abridged) Palomar Gattini-IR is a new wide-field, near-infrared robotic time domain survey operating at Palomar Observatory. Using a 30 cm telescope mounted with a H2RG detector, Gattini-IR achieves a field of view of 25 sq. deg. with a pixel scale of 8.7" in J-band. Here, we describe the system design, survey operations, data processing system and on-sky performance of Palomar Gattini-IR. As a part of the nominal survey, Gattini-IR scans $\approx 7500$ square degrees of the sky every night to a median 5$σ$ depth of $15.7$ AB mag outside the Galactic plane. The survey covers $\approx 15000$ square degrees of the sky visible from Palomar with a median cadence of 2 days. A real-time data processing system produces stacked science images from dithered raw images taken on sky, together with PSF-fit source catalogs and transient candidates identified from subtractions within a median delay of $\approx 4$ hours from the time of observation. The calibrated data products achieve an astrometric accuracy (RMS) of $\approx 0.7$" with respect to Gaia DR2 for sources with S/N $> 10$, and better than $\approx 0.35$" for sources brighter than $\approx 12$ Vega mag. The photometric accuracy (RMS) achieved in the PSF-fit source catalogs is better than $\approx 3$% for sources brighter than $\approx 12$ Vega mag, as calibrated against the 2MASS catalog. With a field of view $\approx 40\times$ larger than any other existing near infrared imaging instrument, Gattini-IR is probing the reddest and dustiest transients in the local universe such as dust obscured supernovae in nearby galaxies, novae behind large columns of extinction within the galaxy, reddened micro-lensing events in the Galactic plane and variability from cool and dust obscured stars. We present results from transients and variables identified since the start of the commissioning period. △ Less

Submitted 29 October, 2019; originally announced October 2019.

Comments: 35 pages, 29 figures. Submitted to PASP

GROWTH on S190425z: Searching thousands of square degrees to identify an optical or infrared counterpart to a binary neutron star merger with the Zwicky Transient Facility and Palomar Gattini IR

Authors: Michael W. Coughlin, Tomás Ahumada, Shreya Anand, Kishalay De, Matthew J. Hankins, Mansi M. Kasliwal, Leo P. Singer, Eric C. Bellm, Igor Andreoni, S. Bradley Cenko, Jeff Cooke, Christopher M. Copperwheat, Alison M. Dugas, Jacob E. Jencson, Daniel A. Perley, Po-Chieh Yu, Varun Bhalerao, Harsh Kumar, Joshua S. Bloom, G. C. Anupama, Michael C. B. Ashley, Ashot Bagdasaryan, Rahul Biswas, David A. H. Buckley, Kevin B. Burdge , et al. (54 additional authors not shown)

Abstract: The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most… ▽ More The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg$^2$ of the initial skymap over the next two nights, corresponding to 46\% integrated probability, ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% integrated probability based on the initial sky map. The revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient "alerts" over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical lightcurves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled-out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible. △ Less

Submitted 4 October, 2019; v1 submitted 29 July, 2019; originally announced July 2019.

Discovery of an intermediate-luminosity red transient in M51 and its likely dust-obscured, infrared-variable progenitor

Authors: Jacob E. Jencson, Scott M. Adams, Howard E. Bond, Schuyler D. van Dyk, Mansi M. Kasliwal, John Bally, Nadejda Blagorodnova, Kishalay De, Christoffer Fremling, Yuhan Yao, Andrew Fruchter, David Rubin, Cristina Barbarino, Jesper Sollerman, Adam A. Miller, Erin K. S. Hicks, Matthew A. Malkan, Igor Andreoni, Eric C. Bellm, Robert Buchheim, Richard Dekany, Michael Feeney, Sara Frederick, Avishay Gal-Yam, Robert D. Gehrz , et al. (27 additional authors not shown)

Abstract: We present the discovery of an optical transient (OT) in Messier 51, designated M51 OT2019-1 (also ZTF19aadyppr, AT 2019abn, ATLAS19bzl), by the Zwicky Transient Facility (ZTF). The OT rose over 15 days to an observed luminosity of $M_r=-13$ ($νL_ν=9\times10^6~L_{\odot}$), in the luminosity gap between novae and typical supernovae (SNe). Spectra during the outburst show a red continuum, Balmer emi… ▽ More We present the discovery of an optical transient (OT) in Messier 51, designated M51 OT2019-1 (also ZTF19aadyppr, AT 2019abn, ATLAS19bzl), by the Zwicky Transient Facility (ZTF). The OT rose over 15 days to an observed luminosity of $M_r=-13$ ($νL_ν=9\times10^6~L_{\odot}$), in the luminosity gap between novae and typical supernovae (SNe). Spectra during the outburst show a red continuum, Balmer emission with a velocity width of $\approx400$ km s$^{-1}$, Ca II and [Ca II] emission, and absorption features characteristic of an F-type supergiant. The spectra and multiband light curves are similar to the so-called "SN impostors" and intermediate-luminosity red transients (ILRTs). We directly identify the likely progenitor in archival Spitzer Space Telescope imaging with a $4.5~μ$m luminosity of $M_{[4.5]}\approx-12.2$ and a $[3.6]-[4.5]$ color redder than 0.74 mag, similar to those of the prototype ILRTs SN 2008S and NGC 300 OT2008-1. Intensive monitoring of M51 with Spitzer further reveals evidence for variability of the progenitor candidate at [4.5] in the years before the OT. The progenitor is not detected in pre-outburst Hubble Space Telescope optical and near-IR images. The optical colors during outburst combined with spectroscopic temperature constraints imply a higher reddening of $E(B-V)\approx0.7$ mag and higher intrinsic luminosity of $M_r\approx-14.9$ ($νL_ν=5.3\times10^7~L_{\odot}$) near peak than seen in previous ILRT candidates. Moreover, the extinction estimate is higher on the rise than on the plateau, suggestive of an extended phase of circumstellar dust destruction. These results, enabled by the early discovery of M51 OT2019-1 and extensive pre-outburst archival coverage, offer new clues about the debated origins of ILRTs and may challenge the hypothesis that they arise from the electron-capture induced collapse of extreme asymptotic giant branch stars. △ Less

Submitted 29 July, 2019; v1 submitted 15 April, 2019; originally announced April 2019.

Comments: 21 pages, 5 figures, published in ApJL

Journal ref: The Astrophysical Journal Letters, 880 (2019) L20

The Dynamic Infrared Sky

Authors: Mansi M. Kasliwal, Scott Adams, Igor Andreoni, Michael Ashley, Nadia Blagorodnova, Kishalay De, Danielle Frostig, Gabor Furesz, Jacob Jencson, Matthew Hankins, George Helou, Ryan Lau, Anna Moore, Eran Ofek, Rob Simcoe, Jennifer Sokoloski, Jamie Soon, Samaporn Tinyanont, Tony Travouillon

Abstract: Opening up the dynamic infrared sky for systematic time-domain exploration would yield many scientific advances. Multi-messenger pursuits such as localizing gravitational waves from neutron star mergers and quantifying the nucleosynthetic yields require the infrared. Another multi-messenger endeavor that needs infrared surveyors is the study of the much-awaited supernova in our own Milky Way. Unde… ▽ More Opening up the dynamic infrared sky for systematic time-domain exploration would yield many scientific advances. Multi-messenger pursuits such as localizing gravitational waves from neutron star mergers and quantifying the nucleosynthetic yields require the infrared. Another multi-messenger endeavor that needs infrared surveyors is the study of the much-awaited supernova in our own Milky Way. Understanding shocks in novae, true rates of supernovae and stellar mergers are some other examples of stellar evolution and high energy physics wherein the answers are buried in the infrared. We discuss some of the challenges in the infrared and pathfinders to overcome them. We conclude with recommendations on both infrared discovery engines and infrared follow-up machines that would enable this field to flourish in the next decade. △ Less

Submitted 19 March, 2019; originally announced March 2019.

Comments: Astro2020 Science White Paper for Decadal Survey

Dust Production Rates in the Fomalhaut Debris Disk from SOFIA/FORCAST Mid-infrared Imaging

Authors: J. D. Adams, T. L. Herter, R. M. Lau, C. Trinh, M. Hankins

Abstract: We present the first spatially resolved mid-infrared (37.1 $μ$m) image of the Fomalhaut debris disk. We use PSF fitting and subtraction to distinctly measure the flux from the unresolved component and the debris disk. We measure an infrared excess in the point source of $0.9 \pm 0.2$ Jy, consistent with emission from warm dust in an inner disk structure (Su et al. 2016), and inconsistent with a st… ▽ More We present the first spatially resolved mid-infrared (37.1 $μ$m) image of the Fomalhaut debris disk. We use PSF fitting and subtraction to distinctly measure the flux from the unresolved component and the debris disk. We measure an infrared excess in the point source of $0.9 \pm 0.2$ Jy, consistent with emission from warm dust in an inner disk structure (Su et al. 2016), and inconsistent with a stellar wind origin. We cannot confirm or rule out the presence of a pileup ring (Su et al. 2016) near the star. In the cold region, the 37 $μ$m imaging is sensitive to emission from small, blowout grains, which is an excellent probe of the dust production rate from planetesimal collisions. Under the assumptions that the dust grains are icy aggregates and the debris disk is in steady state, this result is consistent with the dust production rates predicted by Kenyon & Bromley (2008) from theoretical models of icy planet formation. We find a dust luminosity of $(7.9 \pm 0.8) \times 10^{-4}$ L$_\odot$ and a dust mass of 8 -- 16 lunar masses, depending on grain porosity, with $\sim 1$ lunar mass in grains with radius 1 $μ$m -- 1 mm. If the grains are icy and highly porous, meter-sized objects must be invoked to explain the far-IR, submm, and mm emission. If the grains are composed of astronomical silicates, there is a dearth of blowout grains (Pawellek et al. 2014) and the mass loss rate is well below the predicted dust production values. △ Less

Submitted 25 June, 2018; originally announced June 2018.

Comments: 8 pages, 8 figures, 2 tables, accepted for publication in ApJ

An Infrared Study of the Circumstellar Material Associated with the Carbon Star R Sculptoris

Authors: Matthew J. Hankins, Terry L. Herter, Matthias Maercker, Ryan M. Lau, Gregory C. Sloan

Abstract: The asymptotic giant branch (AGB) star R Sculptoris (R Scl) is one of the most extensively studied stars on the AGB. R Scl is a carbon star with a massive circumstellar shell ($M_{shell}\sim 7.3\times10^{-3}~M_{\odot}$) which is thought to have been produced during a thermal pulse event $\sim2200$ years ago. To study the thermal dust emission associated with its circumstellar material, observation… ▽ More The asymptotic giant branch (AGB) star R Sculptoris (R Scl) is one of the most extensively studied stars on the AGB. R Scl is a carbon star with a massive circumstellar shell ($M_{shell}\sim 7.3\times10^{-3}~M_{\odot}$) which is thought to have been produced during a thermal pulse event $\sim2200$ years ago. To study the thermal dust emission associated with its circumstellar material, observations were taken with the Faint Object InfraRed CAMera for the SOFIA Telescope (FORCAST) at 19.7, 25.2, 31.5, 34.8, and 37.1 $μ$m. Maps of the infrared emission at these wavelengths were used to study the morphology and temperature structure of the spatially extended dust emission. Using the radiative transfer code DUSTY and fitting the spatial profile of the emission, we find that a geometrically thin dust shell cannot reproduce the observed spatially resolved emission. Instead, a second dust component in addition to the shell is needed to reproduce the observed emission. This component, which lies interior to the dust shell, traces the circumstellar envelope of R Scl. It is best fit by a density profile with $n \propto r^α$ where $α=0.75^{+0.45}_{-0.25}$ and dust mass of $M_d=9.0^{+2.3}_{-4.1}\times10^{-6}~M_{\odot}$. The strong departure from an $r^{-2}$ law indicates that the mass-loss rate of R Scl has not been constant. This result is consistent with a slow decline in the post-pulse mass-loss which has been inferred from observations of the molecular gas. △ Less

Submitted 6 December, 2017; originally announced December 2017.

Comments: 10 pages, 10 figures, accepted to ApJ

An Infrared Study of the Dust Properties and Geometry of the Arched Filaments HII Region with SOFIA/FORCAST

Authors: Matthew J. Hankins, Ryan M. Lau, Mark R. Morris, Terry L. Herter

Abstract: Massive stellar clusters provide radiation ($\mathrm{\sim 10^7-10^8~L_{\odot}}$) and winds ($\mathrm{\sim 1000~km/s}$) that act to heat dust and shape their surrounding environment. In this paper, the Arched Filaments in the Galactic center were studied to better understand the influence of the Arches cluster on its nearby interstellar medium (ISM). The Arched Filaments were observed with the Fain… ▽ More Massive stellar clusters provide radiation ($\mathrm{\sim 10^7-10^8~L_{\odot}}$) and winds ($\mathrm{\sim 1000~km/s}$) that act to heat dust and shape their surrounding environment. In this paper, the Arched Filaments in the Galactic center were studied to better understand the influence of the Arches cluster on its nearby interstellar medium (ISM). The Arched Filaments were observed with the Faint Object InfraRed CAMera for the SOFIA Telescope (FORCAST) at 19.7, 25.2, 31.5, and 37.1 $μ$m. Color-temperature maps of the region created with the 25.2 and 37.1 $μ$m data reveal relatively uniform dust temperatures (70-100 K) over the extent of the filaments ($\sim 25$ pc). Distances between the cluster and the filaments were calculated assuming equilibrium heating of standard size ISM dust grains ($\sim$0.1 $μ$m). The distances inferred by this method are in conflict with the projected distance between the filaments and the cluster, although this inconsistency can be explained if the characteristic grain size in the filaments is smaller ($\sim$0.01 $μ$m) than typical values. DustEM models of selected locations within the filaments show evidence of depleted abundances of polycyclic aromatic hydrocarbons (PAHs) by factors of $\sim$1.6-10 by mass compared to the diffuse ISM. The evidence for both PAH depletion and a smaller characteristic grain size points to processing of the ISM within the filaments. We argue that the eroding of dust grains within the filaments is not likely attributable to the radiation or winds from the Arches cluster, but may be related to the physical conditions in the Galactic center. △ Less

Submitted 14 February, 2017; originally announced February 2017.

Comments: 38 pages, 10 figures. Accepted to ApJ

Stagnant Shells in the Vicinity of the Dusty Wolf-Rayet-OB Binary WR 112

Authors: R. M. Lau, M. J. Hankins, R. Schödel, J. Sanchez-Bermudez, A. F. J. Moffat, M. E. Ressler

Abstract: We present high spatial resolution mid-infrared images of the nebula around the late-type carbon-rich Wolf-Rayet (WC)-OB binary system WR~112 taken by the recently upgraded VLT spectrometer and imager for the mid-infrared (VISIR) with the PAH1, NeII\_2, and Q3 filters. The observations reveal a morphology resembling a series of arc-like filaments and broken shells. Dust temperatures and masses are… ▽ More We present high spatial resolution mid-infrared images of the nebula around the late-type carbon-rich Wolf-Rayet (WC)-OB binary system WR~112 taken by the recently upgraded VLT spectrometer and imager for the mid-infrared (VISIR) with the PAH1, NeII\_2, and Q3 filters. The observations reveal a morphology resembling a series of arc-like filaments and broken shells. Dust temperatures and masses are derived for each of the identified filamentary structures, which exhibit temperatures ranging from $179_{-6}^{+8}$ K at the exterior W2 filament to $355_{-25}^{+37}$ K in the central 3". The total dust mass summed over the features is $2.6\pm0.4\times10^{-5}$ $\mathrm{M}_\odot$. A multi-epoch analysis of mid-IR photometry of WR~112 over the past $\sim20$ yr reveals no significant variability in the observed dust temperature and mass. The morphology of the mid-IR dust emission from WR~112 also exhibits no significant expansion from imaging data taken in 2001, 2007, and 2016, which disputes the current interpretation of the nebula as a high expansion velocity ($\sim1200$ km s$^{-1}$) "pinwheel"-shaped outflow driven by the central WC-OB colliding-wind binary. An upper limit of $\lesssim120$ km s$^{-1}$ is derived for the expansion velocity assuming a distance of $4.15$ kpc. The upper limit on the average total mass-loss rate from the central 3" of WR~112 is estimated to be $\lesssim8\times10^{-6}$ $\mathrm{M}_\odot$ yr$^{-1}$. We leave its true nature as an open question, but propose that the WR~112 nebula may have formed in the outflow during a previous red or yellow supergiant phase of the central Wolf-Rayet star. △ Less

Submitted 16 December, 2016; originally announced December 2016.

Comments: 7 pages, 3 figures, 2 tables, Accepted for publication in ApJL

Infrared Observations of the Quintuplet Proper Members using SOFIA/FORCAST and Gemini/TReCS

Authors: Matthew J. Hankins, Ryan M. Lau, Mark R. Morris, Joel Sanchez-Bermudez, Jörg-Uwe Pott, Joseph D. Adams, Terry L. Herter

Abstract: Since their discovery, the Quintuplet proper members (QPMs) have been somewhat mysterious in nature. Originally dubbed the "cocoon stars" due to their cool featureless spectra, high-resolution near-infrared imaging observations have shown that at least two of the objects exhibit "pinwheel" nebulae consistent with binary systems with a carbon-rich Wolf-Rayet star and O/B companion. In this paper, w… ▽ More Since their discovery, the Quintuplet proper members (QPMs) have been somewhat mysterious in nature. Originally dubbed the "cocoon stars" due to their cool featureless spectra, high-resolution near-infrared imaging observations have shown that at least two of the objects exhibit "pinwheel" nebulae consistent with binary systems with a carbon-rich Wolf-Rayet star and O/B companion. In this paper, we present 19.7, 25.2, 31.5, and 37.1 $μ$m observations of the QPMs (with an angular resolution of 3.2-3.8") taken with the Faint Object Infrared Camera for the SOFIA Telescope (FORCAST) in conjunction with high-resolution ($\sim$0.1-0.2") images at 8.8 and 11.7 $μ$m from the Thermal-Region Camera Spectrograph (TReCS). DUSTY models of the thermal dust emission of two of the four detected QPMs, Q2 and Q3, are fitted by radial density profiles which are consistent with constant mass loss rates ($ρ_d \propto r^{-2}$). For the two remaining sources, Q1 and Q9, extended structures ($\sim$ 1") are detected around these objects in high-resolution imaging data. Based on the fitted dust masses, Q9 has an unusually large dust reservoir ($\mathrm{M_d}=1.3^{+0.8}_{-0.4}\times 10^{-3} \mathrm{M_{\odot}}$) compared to typical dusty Wolf-Rayet stars which suggests that it may have recently undergone an episode of enhanced mass loss. △ Less

Submitted 10 July, 2016; originally announced July 2016.

Comments: 28 pages, 10 figures. Accepted to ApJ

An Apparent Precessing Helical Outflow from a Massive Evolved Star: Evidence for Binary Interaction

Authors: Ryan M. Lau, Matthew J. Hankins, Terry L. Herter, Mark R. Morris, Elisabeth A. C. Mills, Michael E. Ressler

Abstract: Massive, evolved stars play a crucial role in the metal-enrichment, dust budget, and energetics of the interstellar medium, however, the details of their evolution are uncertain because of their rarity and short lifetimes before exploding as supernovae. Discrepancies between theoretical predictions from single-star evolutionary models and observations of massive stars have evoked a shifting paradi… ▽ More Massive, evolved stars play a crucial role in the metal-enrichment, dust budget, and energetics of the interstellar medium, however, the details of their evolution are uncertain because of their rarity and short lifetimes before exploding as supernovae. Discrepancies between theoretical predictions from single-star evolutionary models and observations of massive stars have evoked a shifting paradigm that implicates the importance of binary interaction. We present mid- to far-infrared observations from the Stratospheric Observatory for Infrared Astronomy (SOFIA) of a conical "helix" of warm dust ($\sim180$ K) that appears to extend from the Wolf-Rayet star WR102c. Our interpretation of the helix is a precessing, collimated outflow that emerged from WR102c during a previous evolutionary phase as a rapidly rotating luminous blue variable. We attribute the precession of WR102c to gravitational interactions with an unseen compact binary companion whose orbital period can be constrained to $800\,\mathrm{d}<P<1400$ d from the inferred precession period, $τ_p\sim1.4\times10^4$ yr, and limits imposed on the stellar and orbital parameters of the system. Our results concur with the range of orbital periods ($P\lesssim1500$ d) where spin-up via mass exchange is expected to occur for massive binary systems. △ Less

Submitted 16 January, 2016; v1 submitted 23 December, 2015; originally announced December 2015.

Comments: 13 pages, 6 figures, Accepted to ApJ

SOFIA/FORCAST Observations of Warm Dust in S106: A Fragmented Environment

Authors: J. D. Adams, T. L. Herter, J. L. Hora, N. Schneider, R. M. Lau, J. G. Staughn, R. Simon, N. Smith, R. D. Gehrz, L. E. Allen, S. Bontemps, S. J. Carey, G. G. Fazio, R. A. Gutermuth, A. Guzman Fernandez, M. Hankins, T. Hill, E. Keto, X. P. Koenig, K. E. Kraemer, S. T. Megeath, D. R. Mizuno, F. Motte, P. C. Myers, H. A. Smith

Abstract: We present mid-IR (19 - 37 microns) imaging observations of S106 from SOFIA/FORCAST, complemented with IR observations from Spitzer/IRAC (3.6 - 8.0 microns), IRTF/MIRLIN (11.3 and 12.5 microns), and Herschel/PACS (70 and 160 microns). We use these observations, observations in the literature, and radiation transfer modeling to study the heating and composition of the warm (~ 100 K) dust in the reg… ▽ More We present mid-IR (19 - 37 microns) imaging observations of S106 from SOFIA/FORCAST, complemented with IR observations from Spitzer/IRAC (3.6 - 8.0 microns), IRTF/MIRLIN (11.3 and 12.5 microns), and Herschel/PACS (70 and 160 microns). We use these observations, observations in the literature, and radiation transfer modeling to study the heating and composition of the warm (~ 100 K) dust in the region. The dust is heated radiatively by the source S106 IR, with little contributions from grain-electron collisions and Ly-alpha radiation. The dust luminosity is >~ (9.02 +/- 1.01) x 10^4 L_sun, consistent with heating by a mid- to late-type O star. We find a temperature gradient (~ 75 - 107 K) in the lobes, which is consistent with a dusty equatorial geometry around S106 IR. Furthermore, the SOFIA observations resolve several cool (~ 65 - 70 K) lanes and pockets of warmer (~ 75 - 90 K) dust in the ionization shadow, indicating that the environment is fragmented. We model the dust mass as a composition of amorphous silicates, amorphous carbon, big grains, very small grains, and PAHs. We present the relative abundances of each grain component for several locations in S106. △ Less

Submitted 14 October, 2015; originally announced October 2015.

Comments: 18 pages, 11 figures, 4 tables, accepted by ApJ