Showing 1–50 of 79 results for author: Lau, R

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

Overview Results of JWST Observations of Star-Forming Clusters in the Extreme Outer Galaxy

Authors: Natsuko Izumi, Michael E. Ressler, Ryan M. Lau, Patrick M. Koch, Masao Saito, Naoto Kobayashi, Chikako Yasui

Abstract: The extreme outer Galaxy (EOG), which we define as the region of the Milky Way with a galactocentric radius of more than 18 kpc, provides an excellent opportunity to study star formation in an environment significantly different from that in the solar neighborhood because of its lower metallicity and lower gas density. We carried out near- and mid-infrared (NIR and MIR) imaging observations toward… ▽ More The extreme outer Galaxy (EOG), which we define as the region of the Milky Way with a galactocentric radius of more than 18 kpc, provides an excellent opportunity to study star formation in an environment significantly different from that in the solar neighborhood because of its lower metallicity and lower gas density. We carried out near- and mid-infrared (NIR and MIR) imaging observations toward two star-forming clusters located in the EOG using JWST NIRCam and MIRI with nine filters: F115W, F150W, F200W, F350W, F405N, F444W, F770W, F1280W, and F2100W. In this paper, we present an overview of the observations, data reduction, and initial results. The NIR sensitivity is approximately 10--80 times better than our previous observation with the Subaru 8.2 m telescope. Accordingly, the mass detection limit reaches to about 0.01--0.05 $M_\odot$, which is about 10 times better than the previous observations. At MIR wavelengths, the high sensitivity and resolution data enable us to resolve individual young stellar objects in such a distant region for the first time. The mass detection limit at MIR F770W filter reaches about 0.1--0.3 $M_\odot$. With these new observations, we have identified components of the clusters that previous surveys did not detect, including class 0 candidates, outflow/jet components, and distinctive nebular structures. These data will enable us to investigate the properties of star formation in the EOG at the same depth of detail as previous observations of star formation in the solar neighborhood. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: 39 pages, 21 figures, and 5 tables. Accepted for publication in Astronomical Journal. The figure with the full resolution is shown in the publication in the Astronomical Journal

Towards quantum-enhanced long-baseline optical/near-IR interferometry

Authors: Jayadev K. Rajagopal, Ryan M. Lau, Isack Padilla, Stephen T. Ridgway, Chaohan Cui, Brittany McClinton, Aqil Sajjad, Stuartt Corder, Mark Rawlings, Fredrik Rantakyro, J. Gabriel Richardson, Amit Ashok, Saikat Guha

Abstract: Microarcsecond resolutions afforded by an optical-NIR array with kilometer-baselines would enable breakthrough science. However significant technology barriers exist in transporting weakly coherent photon states over these distances: primarily photon loss and phase errors. Quantum telescopy, using entangled states to link spatially separated apertures, offers a possible solution to the loss of pho… ▽ More Microarcsecond resolutions afforded by an optical-NIR array with kilometer-baselines would enable breakthrough science. However significant technology barriers exist in transporting weakly coherent photon states over these distances: primarily photon loss and phase errors. Quantum telescopy, using entangled states to link spatially separated apertures, offers a possible solution to the loss of photons. We report on an initiative launched by NSF NOIRLab in collaboration with the Center for Quantum Networks and Arizona Quantum Initiative at the University of Arizona, Tucson, to explore these concepts further. A brief description of the quantum concepts and a possible technology roadmap towards a quantum-enhanced very long baseline optical-NIR interferometric array is presented. An on-sky demonstration of measuring spatial coherence of photons with apertures linked through the simplest Gottesman protocol over short baselines and with limited phase fluctuations is envisaged as the first step. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: Proceeding of SPIE Conference "Astronomical Telescopes + Instrumentation" (June 2024)

Report number: Paper No. 13095-58

The first Palomar Gattini-IR catalog of J-band light curves: construction and public data release

Authors: Shion Murakawa, Kishalay De, Michael C. B. Ashley, Nicholas Earley, Lynne A. Hillenbrand, Mansi M. Kasliwal, Ryan M. Lau, Anna M. Moore, J. L. Sokoloski, Roberto Soria

Abstract: Palomar Gattini-IR is a wide-field, synoptic infrared time domain survey covering $\approx 15000$ sq. deg. of the sky at $\approx 1-3$ night cadence to a depth of $J\approx 13.0$ and $\approx 14.9$ Vega mag in and outside the Galactic plane, respectively. Here, we present the first data release of $J$-band light curves of 2MASS sources within the survey footprint covering approximately the first f… ▽ More Palomar Gattini-IR is a wide-field, synoptic infrared time domain survey covering $\approx 15000$ sq. deg. of the sky at $\approx 1-3$ night cadence to a depth of $J\approx 13.0$ and $\approx 14.9$ Vega mag in and outside the Galactic plane, respectively. Here, we present the first data release of $J$-band light curves of 2MASS sources within the survey footprint covering approximately the first four years of operations. We describe the construction of the source catalog based on the 2MASS point sources, followed by exposure filtering criteria and forced PSF photometry. The catalog contains light curves of $\approx 286$ million unique sources with 2MASS magnitudes of $J < 15.5$ mag, with a total of $\approx 50$ billion photometric measurements and $\approx 20$ billion individual source detections at signal-to-noise-ratio $> 3$. We demonstrate the photometric fidelity of the catalog by i) quantifying the magnitude-dependent accuracy and uncertainty of the photometry with respect to 2MASS and ii) comparing against forced PGIR aperture photometry for known variable sources. We present simple filtering criteria for selecting reliable photometric measurements as well as example Python notebooks for users. This catalog is the largest compilation of nightly cadence, synoptic infrared light curves to date, comparable to those in the largest optical surveys, providing a stepping stone to upcoming infrared surveys in the coming decade. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: 10 Pages, 5 figures, submitted to PASP. Full catalog will be available upon publication, individual requests welcome!

WTP19aalnxx: Discovery of a bright mid-infrared transient in the emerging class of low luminosity supernovae revealed by delayed circumstellar interaction

Authors: Charlotte Myers, Kishalay De, Lin Yan, Jacob E. Jencson, Nicholas Earley, Christoffer Fremling, Daichi Hiramatsu, Mansi M. Kasliwal, Ryan M. Lau, Morgan MacLeod, Megan Masterson, Christos Panagiotou, Robert Simcoe, Samaporn Tinyanont

Abstract: While core-collapse supernovae (SNe) often show early and consistent signs of circumstellar (CSM) interaction, some exhibit delayed signatures due to interaction with distant material around the progenitor star. Here we present the discovery in NEOWISE data of WTP19aalnxx, a luminous mid-infrared (IR) transient in the outskirts of the galaxy KUG 0022-007 at $\approx 190$ Mpc. First detected in 201… ▽ More While core-collapse supernovae (SNe) often show early and consistent signs of circumstellar (CSM) interaction, some exhibit delayed signatures due to interaction with distant material around the progenitor star. Here we present the discovery in NEOWISE data of WTP19aalnxx, a luminous mid-infrared (IR) transient in the outskirts of the galaxy KUG 0022-007 at $\approx 190$ Mpc. First detected in 2018, WTP19aalnxx reaches a peak absolute (Vega) magnitude of $\approx-22$ at $4.6 \, μ$m in $\approx3$ yr, comparable to the most luminous interacting SNe. Archival data reveal a $\gtrsim 5\times$ fainter optical counterpart detected since 2015, while follow-up near-IR observations in 2022 reveal an extremely red ($Ks-W2 \approx 3.7$ mag) active transient. Deep optical spectroscopy confirm strong CSM interaction signatures via intermediate-width Balmer emission lines and coronal metal lines. Modeling the broadband spectral energy distribution, we estimate the presence of $\gtrsim 10^{-2}$ M$_\odot$ of warm dust, likely formed in the shock interaction region. Together with the lack of nebular Fe emission, we suggest that WTP19aalnxx is a missed, low (optical) luminosity SN in an emerging family of core-collapse SNe distinguished by their CSM-interaction-powered mid-IR emission that outshines the optical bands. Investigating the Zwicky Transient Facility sample of SNe in NEOWISE data, we find $17$ core-collapse SNe ($\gtrsim 3$% in a volume-limited sample) without early signs of CSM interaction that exhibit delayed IR brightening, suggestive of dense CSM shells at $\lesssim 10^{17}$cm. We suggest that synoptic IR surveys offer a new route to revealing late-time CSM interaction and the prevalence of intense terminal mass loss in massive stars. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 15 pages, 5 figures, submitted to ApJ

The long-period spectroscopic orbit and dust creation in the Wolf-Rayet binary system WR 125

Authors: Noel D. Richardson, Andrea R. Daly, Peredur M. Williams, Grant M. Hill, Victor I. Shenavrin, Izumi Endo, André-Nicolas Chené, Nicole Karnath, Ryan M. Lau, Anthony F. J. Moffat, Gerd Weigelt

Abstract: Several long-period binaries with a carbon-rich Wolf-Rayet star and an O star produce dust in their wind collisions. In eccentric binaries, this is seen most strongly near periastron passage. The exact conditions leading to dust creation require orbital properties to be determined, which is difficult owing to their long periods. Recently, the binary system WR 125 (WC7+O9III) began a dust creation… ▽ More Several long-period binaries with a carbon-rich Wolf-Rayet star and an O star produce dust in their wind collisions. In eccentric binaries, this is seen most strongly near periastron passage. The exact conditions leading to dust creation require orbital properties to be determined, which is difficult owing to their long periods. Recently, the binary system WR 125 (WC7+O9III) began a dust creation episode seen through an infrared outburst first detected by NEOWISE-R, which was the first outburst detected since 1991. We present new near- and mid-infrared photometry, which we use to show consistency between the two outbursts and derive an orbital period of 28.12$^{+0.10}_{-0.05}$ yr. We use a long time-series of optical spectra to place the first constraints on its orbital elements, on the assumption that this system will produce dust near periastron. The orbit has a mild eccentricity of 0.29$\pm$0.12 and is only derived for the Wolf-Rayet component, as the O star's radial velocities have noise that is likely larger than the expected semi-amplitude of the orbit. We also present SOFIA/FORCAST grism spectroscopy to examine the infrared spectral energy distribution (SED) of the dust during this outburst, comparing its properties to other WCd binaries, deriving a dust temperature of 580 K in 2021. This collection of observations will allow us to plan future observations of this system and place the system in the context of dust-creating Wolf-Rayet binaries. △ Less

Submitted 16 May, 2024; originally announced May 2024.

Comments: accepted to ApJ

Emission lines due to ionizing radiation from a compact object in the remnant of Supernova 1987A

Authors: C. Fransson, M. J. Barlow, P. J. Kavanagh, J. Larsson, O. C. Jones, B. Sargent, M. Meixner, P. Bouchet, T. Temim, G. S. Wright, J. A. D. L. Blommaert, N. Habel, A. S. Hirschauer, J. Hjorth, L. Lenkić, T. Tikkanen, R. Wesson, A. Coulais, O. D. Fox, R. Gastaud, A. Glasse, J. Jaspers, O. Krause, R. M. Lau, O. Nayak , et al. (9 additional authors not shown)

Abstract: The nearby Supernova 1987A was accompanied by a burst of neutrino emission, which indicates that a compact object (a neutron star or black hole) was formed in the explosion. There has been no direct observation of this compact object. In this work, we observe the supernova remnant with JWST spectroscopy finding narrow infrared emission lines of argon and sulphur. The line emission is spatially unr… ▽ More The nearby Supernova 1987A was accompanied by a burst of neutrino emission, which indicates that a compact object (a neutron star or black hole) was formed in the explosion. There has been no direct observation of this compact object. In this work, we observe the supernova remnant with JWST spectroscopy finding narrow infrared emission lines of argon and sulphur. The line emission is spatially unresolved and blueshifted in velocity relative to the supernova rest frame. We interpret the lines as gas illuminated by a source of ionizing photons located close to the center of the expanding ejecta. Photoionization models show that the line ratios are consistent with ionization by a cooling neutron star or pulsar wind nebula. The velocity shift could be evidence for a neutron star natal kick. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: Authors version of manuscript published in Science on 22 Feb 2024

Journal ref: SCIENCE 22 Feb 2024 Vol 383, Issue 6685 pp. 898-903

JWST MIRI Imager Observations of Supernova SN 1987A

Authors: P. Bouchet, R. Gastaud, A. Coulais, M. J. Barlow, C. Fransson, P. J. Kavanagh, J. Larsson, T. Temim, O. C. Jones, A. S. Hirschauer, T. Tikkanen, J. A. D. L. Blommaert, O. D. Fox, A. Glasse, N. Habel, J. Hjorth, J. Jaspers, O. Krause, R. M. Lau, L. Lenkić, M. Meixner, O. Nayak, A. Rest, B. Sargent, R. Wesson , et al. (9 additional authors not shown)

Abstract: There exist very few mid-infrared (IR) observations of supernovae (SNe) in general. Therefore, SN 1987A, the closest visible SN in 400 years, gives us the opportunity to explore the mid-IR properties of SNe, the dust in their ejecta and surrounding medium, and to witness the birth of a SN remnant (SNR). The James Webb Space Telescope (JWST), with its high spatial resolution and extreme sensitivity… ▽ More There exist very few mid-infrared (IR) observations of supernovae (SNe) in general. Therefore, SN 1987A, the closest visible SN in 400 years, gives us the opportunity to explore the mid-IR properties of SNe, the dust in their ejecta and surrounding medium, and to witness the birth of a SN remnant (SNR). The James Webb Space Telescope (JWST), with its high spatial resolution and extreme sensitivity, gives a new view on these issues. We report on the first imaging observations obtained with the Mid-InfraRed Instrument (MIRI). We build temperature maps and discuss the morphology of the nascent SNR. Our results show that the temperatures in the equatorial ring (ER) are quite non-uniform. This could be due to dust destruction in some parts of the ring, as had been assumed in some previous works. We show that the IR emission extends beyond the ER, illustrating the fact that the shock wave has now passed through this ring to affect the circumstellar medium on a larger scale. Finally, while sub-mm Atacama Large Millimeter Array (ALMA) observations have hinted at the location of the compact remnant of SN 1987A, we note that our MIRI data have found no such evidence. △ Less

Submitted 21 February, 2024; originally announced February 2024.

Comments: 19 pages, 19 figures, 2 tables; Accepted for publication in the Astrophysical Journal (February 2, 2024)

Windows on the Universe: Establishing the Infrastructure for a Collaborative Multi-messenger Ecosystem

Authors: The 2023 Windows on the Universe Workshop White Paper Working Group, T. Ahumada, J. E. Andrews, S. Antier, E. Blaufuss, P. R. Brady, A. M. Brazier, E. Burns, S. B. Cenko, P. Chandra, D. Chatterjee, A. Corsi, M. W. Coughlin, D. A. Coulter, S. Fu, A. Goldstein, L. P. Guy, E. J. Hooper, S. B. Howell, T. B. Humensky, J. A. Kennea, S. M. Jarrett, R. M. Lau, T. R. Lewis, L. Lu , et al. (21 additional authors not shown)

Abstract: In this White Paper, we present recommendations for the scientific community and funding agencies to foster the infrastructure for a collaborative multi-messenger and time-domain astronomy (MMA/TDA) ecosystem. MMA/TDA is poised for breakthrough discoveries in the coming decade. In much the same way that expanding beyond the optical bandpass revealed entirely new and unexpected discoveries, cosmic… ▽ More In this White Paper, we present recommendations for the scientific community and funding agencies to foster the infrastructure for a collaborative multi-messenger and time-domain astronomy (MMA/TDA) ecosystem. MMA/TDA is poised for breakthrough discoveries in the coming decade. In much the same way that expanding beyond the optical bandpass revealed entirely new and unexpected discoveries, cosmic messengers beyond light (i.e., gravitational waves, neutrinos, and cosmic rays) open entirely new windows to answer some of the most fundamental questions in (astro)physics: heavy element synthesis, equation of state of dense matter, particle acceleration, etc. This field was prioritized as a frontier scientific pursuit in the 2020 Decadal Survey on Astronomy and Astrophysics via its "New Windows on the Dynamic Universe" theme. MMA/TDA science presents technical challenges distinct from those experienced in other disciplines. Successful observations require coordination across myriad boundaries -- different cosmic messengers, ground vs. space, international borders, etc. -- all for sources that may not be well localized, and whose brightness may be changing rapidly with time. Add that all of this work is undertaken by real human beings, with distinct backgrounds, experiences, cultures, and expectations, that often conflict. To address these challenges and help MMA/TDA realize its full scientific potential in the coming decade (and beyond), the second in a series of community workshops sponsored by the U.S. National Science Foundation (NSF) and NASA titled "Windows on the Universe: Establishing the Infrastructure for a Collaborative Multi-Messenger Ecosystem" was held on October 16-18, 2023 in Tucson, AZ. Here we present the primary recommendations from this workshop focused on three key topics -- hardware, software, and people and policy. [abridged] △ Less

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

Comments: Workshop white paper

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

Serendipitous detection of the dusty Type IIL SN 1980K with JWST/MIRI

Authors: Szanna Zsíros, Tamás Szalai, Ilse De Looze, Arkaprabha Sarangi, Melissa Shahbandeh, Ori D. Fox, Tea Temim, Dan Milisavljevic, Schuyler D. Van Dyk, Nathan Smith, Alexei V. Filippenko, Thomas G. Brink, WeiKang Zheng, Luc Dessart, Jacob Jencson, Joel Johansson, Justin Pierel, Armin Rest, Samaporn Tinyanont, Maria Niculescu-Duvaz, M. J. Barlow, Roger Wesson, Jennifer Andrews, Geoff Clayton, Kishalay De , et al. (17 additional authors not shown)

Abstract: We present mid-infrared (mid-IR) imaging of the Type IIL supernova (SN) 1980K with the James Webb Space Telescope (JWST) more than 40 yr post-explosion. SN 1980K, located in the nearby ($D\approx7$ Mpc) "SN factory" galaxy NGC 6946, was serendipitously captured in JWST/MIRI images taken of the field of SN 2004et in the same galaxy. SN 1980K serves as a promising candidate for studying the transiti… ▽ More We present mid-infrared (mid-IR) imaging of the Type IIL supernova (SN) 1980K with the James Webb Space Telescope (JWST) more than 40 yr post-explosion. SN 1980K, located in the nearby ($D\approx7$ Mpc) "SN factory" galaxy NGC 6946, was serendipitously captured in JWST/MIRI images taken of the field of SN 2004et in the same galaxy. SN 1980K serves as a promising candidate for studying the transitional phase between young SNe and older SN remnants and also provides a great opportunity to investigate its the close environment. SN 1980K can be identified as a clear and bright point source in all eight MIRI filters from F560W up to F2550W. We fit analytical dust models to the mid-IR spectral energy distribution that reveal a large amount ($M_d \approx 0.002 {M}_{\odot}$) of Si-dominated dust at $T_{dust}\approx 150$ K (accompanied by a hotter dust/gas component), and also computed numerical SED dust models. Radiative transfer modeling of a late-time optical spectrum obtained recently with Keck discloses that an even larger ($\sim 0.24-0.58~{M}_{\odot}$) amount of dust is needed in order for selective extinction to explain the asymmetric line profile shapes observed in SN 1980K. As a conclusion, with JWST, we may see i) pre-existing circumstellar dust heated collisionally (or, partly radiatively), analogous to the equatorial ring of SN 1987A, or ii) the mid-IR component of the presumed newly-formed dust, accompanied by much more colder dust present in the ejecta (as suggested by the late-time the optical spectra). △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 14 pages, 9 figures

FORCASTing the spectroscopic dust properties of the WC+O binary WR137 with SOFIA

Authors: Megan J. Peatt, Noel D. Richardson, Peredur M. Williams, Nicole Karnath, Victor I. Shenavrin, Ryan M. Lau, Anthony F. J. Moffat, Gerd Weigelt

Abstract: WR 137 (HD 192641) is a binary system consisting of a carbon-rich Wolf-Rayet star and an Oe companion star in a 13-year orbit. Near periastron, the winds of the two stars collide and form carbonaceous dust. We obtained three mid-infrared grism spectra of the system with SOFIA and FORCAST during the last year of SOFIA's operations in July 2021, February 2021, and May 2022 (Cycle 9). Within these sp… ▽ More WR 137 (HD 192641) is a binary system consisting of a carbon-rich Wolf-Rayet star and an Oe companion star in a 13-year orbit. Near periastron, the winds of the two stars collide and form carbonaceous dust. We obtained three mid-infrared grism spectra of the system with SOFIA and FORCAST during the last year of SOFIA's operations in July 2021, February 2021, and May 2022 (Cycle 9). Within these spectra, we have identified several wind lines from He I, He II, C III, and C IV that are emitted from the Wolf-Rayet wind as well as a weak emission feature around 6.3-6.4 $μ$m that may have shifted its peak flux from 6.29 to 6.41$μ$m through this time period. The weak feature grew as the continuum dust emission grew while the WR emission appeared to decline due to lower contrast with the continuum. Furthermore, we observe that the peak of the feature shifts to redder wavelengths during the observations. We compare this feature to the UIR feature and other emission lines identified in dusty WC binaries. For WR 137, we speculate that mixing of the winds in the system with the Oe star's disk is important for starting the dust formation and that it is less important as dust formation continues. Previous infrared photometry shows "mini-eruptions" of dust production which could then be explained with variations of the Oe star disk. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 13 pages, accepted to ApJ

Ejecta, Rings, and Dust in SN 1987A with JWST MIRI/MRS

Authors: O. C. Jones, P. J. Kavanagh, M. J. Barlow, T. Temim, C. Fransson, J. Larsson, J. A. D. L. Blommaert, M. Meixner, R. M. Lau, B. Sargent, P. Bouchet, J. Hjorth, G. S. Wright, A. Coulais, O. D. Fox, R. Gastaud, A. Glasse, N. Habel, A. S. Hirschauer, J. Jaspers, O. Krause, Lenkić, O. Nayak, A. Rest, T. Tikkanen , et al. (9 additional authors not shown)

Abstract: Supernova (SN) 1987A is the nearest supernova in $\sim$400 years. Using the {\em JWST} MIRI Medium Resolution Spectrograph, we spatially resolved the ejecta, equatorial ring (ER) and outer rings in the mid-infrared 12,927 days after the explosion. The spectra are rich in line and dust continuum emission, both in the ejecta and the ring. Broad emission lines (280-380~km~s$^{-1}$ FWHM) seen from all… ▽ More Supernova (SN) 1987A is the nearest supernova in $\sim$400 years. Using the {\em JWST} MIRI Medium Resolution Spectrograph, we spatially resolved the ejecta, equatorial ring (ER) and outer rings in the mid-infrared 12,927 days after the explosion. The spectra are rich in line and dust continuum emission, both in the ejecta and the ring. Broad emission lines (280-380~km~s$^{-1}$ FWHM) seen from all singly-ionized species originate from the expanding ER, with properties consistent with dense post-shock cooling gas. Narrower emission lines (100-170~km~s$^{-1}$ FWHM) are seen from species originating from a more extended lower-density component whose high ionization may have been produced by shocks progressing through the ER, or by the UV radiation pulse associated with the original supernova event. The asymmetric east-west dust emission in the ER has continued to fade, with constant temperature, signifying a reduction in dust mass. Small grains in the ER are preferentially destroyed, with larger grains from the progenitor surviving the transition from SN into SNR. The ER is fit with a single set of optical constants, eliminating the need for a secondary featureless hot dust component. We find several broad ejecta emission lines from [Ne~{\sc ii}], [Ar~{\sc ii}], [Fe~{\sc ii}], and [Ni~{\sc ii}]. With the exception of [Fe~{\sc ii}]~25.99$μ$m, these all originate from the ejecta close to the ring and are likely being excited by X-rays from the interaction. The [Fe~{\sc ii}]~5.34$μ$m to 25.99$μ$m line ratio indicates a temperature of only a few hundred K in the inner core, consistent with being powered by ${}^{44}$Ti decay. △ Less

Submitted 29 February, 2024; v1 submitted 13 July, 2023; originally announced July 2023.

Comments: 27 pages, 16 figures, 4 tables. Accepted ApJ

Impact of Pycnonuclear Fusion Uncertainties on the Cooling of Accreting Neutron Star Crusts

Authors: R. Jain, E. F. Brown, H. Schatz, A. V. Afanasjev, M. Beard, L. R. Gasques, S. S. Gupta, G. W. Hitt, W. R. Hix, R. Lau, P. Moller, W. J. Ong, M. Wiescher, Y. Xu

Abstract: The observation of X-rays during quiescence from transiently accreting neutron stars provides unique clues about the nature of dense matter. This, however, requires extensive modeling of the crusts and matching the results to observations. The pycnonuclear fusion reaction rates implemented in these models are theoretically calculated by extending phenomenological expressions and have large uncerta… ▽ More The observation of X-rays during quiescence from transiently accreting neutron stars provides unique clues about the nature of dense matter. This, however, requires extensive modeling of the crusts and matching the results to observations. The pycnonuclear fusion reaction rates implemented in these models are theoretically calculated by extending phenomenological expressions and have large uncertainties spanning many orders of magnitude. We present the first sensitivity studies of these pycnonuclear fusion reactions in realistic network calculations. We also couple the reaction network with the thermal evolution code dStar to further study their impact on the neutron star cooling curves in quiescence. Varying the pycnonuclear fusion reaction rates alters the depth at which nuclear heat is deposited although the total heating remains constant. The enhancement of the pycnonuclear fusion reaction rates leads to an overall shallower deposition of nuclear heat. The impurity factors are also altered depending on the type of ashes deposited on the crust. These total changes correspond to a variation of up to 9 eV in the modeled cooling curves. While this is not sufficient to explain the shallow heat source, it is comparable to the observational uncertainties and can still be important for modeling the neutron star crust. △ Less

Submitted 1 October, 2023; v1 submitted 20 June, 2023; originally announced June 2023.

Comments: AASTeX63, 11 pages with 9 figures

Journal ref: ApJ 955 51 (2023)

A Luminous Red Supergiant and Dusty Long-period Variable Progenitor for SN 2023ixf

Authors: Jacob E. Jencson, Jeniveve Pearson, Emma R. Beasor, Ryan M. Lau, Jennifer E. Andrews, K. Azalee Bostroem, Yize Dong, Michael Engesser, Sebastian Gomez, Muryel Guolo, Emily Hoang, Griffin Hosseinzadeh, Saurabh W. Jha, Viraj Karambelkar, Mansi M. Kasliwal, Michael Lundquist, Nicolas E. Meza Retamal, Armin Rest, David J. Sand, Melissa Shahbandeh, Manisha Shrestha, Nathan Smith, Jay Strader, Stefano Valenti, Qinan Wang , et al. (1 additional authors not shown)

Abstract: We analyze pre-explosion near- and mid-infrared (IR) imaging of the site of SN 2023ixf in the nearby spiral galaxy M101 and characterize the candidate progenitor star. The star displays compelling evidence of variability with a possible period of $\approx$1000 days and an amplitude of $Δm \approx 0.6$ mag in extensive monitoring with the Spitzer Space Telescope since 2004, likely indicative of rad… ▽ More We analyze pre-explosion near- and mid-infrared (IR) imaging of the site of SN 2023ixf in the nearby spiral galaxy M101 and characterize the candidate progenitor star. The star displays compelling evidence of variability with a possible period of $\approx$1000 days and an amplitude of $Δm \approx 0.6$ mag in extensive monitoring with the Spitzer Space Telescope since 2004, likely indicative of radial pulsations. Variability consistent with this period is also seen in the near-IR $J$ and $K_{s}$ bands between 2010 and 2023, up to just 10 days before the explosion. Beyond the periodic variability, we do not find evidence for any IR-bright pre-supernova outbursts in this time period. The IR brightness ($M_{K_s} = -10.7$ mag) and color ($J-K_{s} = 1.6$ mag) of the star suggest a luminous and dusty red supergiant. Modeling of the phase-averaged spectral energy distribution (SED) yields constraints on the stellar temperature ($T_{\mathrm{eff}} = 3500_{-1400}^{+800}$ K) and luminosity ($\log L/L_{\odot} = 5.1\pm0.2$). This places the candidate among the most luminous Type II supernova progenitors with direct imaging constraints, with the caveat that many of these rely only on optical measurements. Comparison with stellar evolution models gives an initial mass of $M_{\mathrm{init}} = 17\pm4 M_{\odot}$. We estimate the pre-supernova mass-loss rate of the star between 3 and 19 yr before explosion from the SED modeling at $\dot M \approx 3\times10^{-5}$ to $3\times10^{-4} M_{\odot}$ yr$^{-1}$ for an assumed wind velocity of $v_w = 10$ km s$^{-1}$, perhaps pointing to enhanced mass loss in a pulsation-driven wind. △ Less

Submitted 1 August, 2023; v1 submitted 14 June, 2023; originally announced June 2023.

Comments: 13 pages, 5 figures, published in ApJL, replacement with revisions to match published version

Journal ref: ApJL 952 (2023) L30

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

JWST NIRSpec observations of Supernova 1987A -- from the inner ejecta to the reverse shock

Authors: J. Larsson, C. Fransson, B. Sargent, O. C. Jones, M. J. Barlow, P. Bouchet, M. Meixner, J. A. D. L. Blommaert, A. Coulais, O. D. Fox, R. Gastaud, A. Glasse, N. Habel, A. S. Hirschauer, J. Hjorth, J. Jaspers, P. J. Kavanagh, O. Krause, R. M. Lau, L. Lenkic, O. Nayak, A. Rest, T. Temim, T. Tikkanen, R. Wesson , et al. (1 additional authors not shown)

Abstract: We present initial results from JWST NIRSpec integral field unit observations of the nearby Supernova (SN) 1987A. The observations provide the first spatially-resolved spectroscopy of the ejecta and equatorial ring (ER) over the 1-5 μm range. We construct 3D emissivity maps of the [Fe I] 1.443 μm line from the inner ejecta and the He I 1.083 μm line from the reverse shock (RS), where the former pr… ▽ More We present initial results from JWST NIRSpec integral field unit observations of the nearby Supernova (SN) 1987A. The observations provide the first spatially-resolved spectroscopy of the ejecta and equatorial ring (ER) over the 1-5 μm range. We construct 3D emissivity maps of the [Fe I] 1.443 μm line from the inner ejecta and the He I 1.083 μm line from the reverse shock (RS), where the former probes the explosion geometry and the latter traces the structure of the circumstellar medium. We also present a model for the integrated spectrum of the ejecta. The [Fe I] 3D map reveals a highly-asymmetric morphology resembling a broken dipole, dominated by two large clumps with velocities of ~2300 km/s. We also find evidence that the Fe-rich inner ejecta have started to interact with the RS. The RS surface traced by the He I line extends from just inside the ER to higher latitudes on both sides of the ER with a half-opening angle ~45 degrees, forming a bubble-like structure. The spectral model for the ejecta allows us to identify the many emission lines, including numerous H_2 lines. We find that the H_2 is most likely excited by far-UV emission, while the metal lines ratios are consistent with a combination of collisional excitation and recombination in the low-temperature ejecta. We also find several high-ionization coronal lines from the ER, requiring a temperature > 2 \times 10^6 K. △ Less

Submitted 16 May, 2023; v1 submitted 7 February, 2023; originally announced February 2023.

Comments: Accepted for publication in ApJL

JWST Discovery of Dust Reservoirs in Nearby Type IIP Supernovae 2004et and 2017eaw

Authors: Melissa Shahbandeh, Arkaprabha Sarangi, Tea Temim, Tamas Szalai, Ori D. Fox, Samaporn Tinyanont, Eli Dwek, Luc Dessart, Alexei V. Filippenko, Thomas G. Brink, Ryan J. Foley, Jacob Jencson, Justin Pierel, Szanna Zsiros, Armin Rest, WeiKang Zheng, Jennifer Andrews, Geoffrey C. Clayton, Kishalay De, Michael Engesser, Suvi Gezari, Sebastian Gomez, Shireen Gonzaga, Joel Johansson, Mansi Kasliwal , et al. (14 additional authors not shown)

Abstract: Supernova (SN) explosions have been sought for decades as a possible source of dust in the Universe, providing the seeds of galaxies, stars, and planetary systems. SN 1987A offers one of the most promising examples of significant SN dust formation, but until the James Webb Space Telescope (JWST), instruments have traditionally lacked the sensitivity at both late times (>1 yr post-explosion) and lo… ▽ More Supernova (SN) explosions have been sought for decades as a possible source of dust in the Universe, providing the seeds of galaxies, stars, and planetary systems. SN 1987A offers one of the most promising examples of significant SN dust formation, but until the James Webb Space Telescope (JWST), instruments have traditionally lacked the sensitivity at both late times (>1 yr post-explosion) and longer wavelengths (i.e., >10 um) to detect analogous dust reservoirs. Here we present JWST/MIRI observations of two historic Type IIP SNe, 2004et and SN 2017eaw, at nearly 18 and 5 yr post-explosion, respectively. We fit the spectral energy distributions as functions of dust mass and temperature, from which we are able to constrain the dust geometry, origin, and heating mechanism. We place a 90% confidence lower limit on the dust masses for SNe 2004et and 2017eaw of >0.014 and >4e-4 M_sun, respectively. More dust may exist at even colder temperatures or may be obscured by high optical depths. We conclude dust formation in the ejecta to be the most plausible and consistent scenario. The observed dust is radiatively heated to ~100-150 K by ongoing shock interaction with the circumstellar medium. Regardless of the best fit or heating mechanism adopted, the inferred dust mass for SN 2004et is the second highest (next to SN 1987A) inferred dust mass in extragalactic SNe thus far, promoting the prospect of SNe as potential significant sources of dust in the Universe. △ Less

Submitted 25 January, 2023; originally announced January 2023.

Comments: 12 pages, 7 figures, submitting to MNRAS

Recurring outbursts of the supernova impostor AT 2016blu in NGC 4559

Authors: Mojgan Aghakhanloo, Nathan Smith, Peter Milne, Jennifer E. Andrews, Schuyler D. Van Dyk, Alexei V. Filippenko, Jacob E. Jencson, Ryan M. Lau, David J. Sand, Samuel Wyatt, WeiKang Zheng

Abstract: We present the first photometric analysis of the supernova (SN) impostor AT 2016blu in NGC 4559. This transient was discovered by the Lick Observatory Supernova Search in 2012 and has continued its outbursts since then. Optical and infrared photometry of AT 2016blu reveals at least 19 outbursts in 2012-2022. Similar photometry from 1999-2009 shows no outbursts, indicating that the star was relativ… ▽ More We present the first photometric analysis of the supernova (SN) impostor AT 2016blu in NGC 4559. This transient was discovered by the Lick Observatory Supernova Search in 2012 and has continued its outbursts since then. Optical and infrared photometry of AT 2016blu reveals at least 19 outbursts in 2012-2022. Similar photometry from 1999-2009 shows no outbursts, indicating that the star was relatively stable in the decade before discovery. Archival {\it Hubble Space Telescope} observations suggest that the progenitor had a minimum initial mass of $M >= 33$ M$_{\odot}$ and a luminosity of $L >= 10^{5.7}$ L$_{\odot}$. AT 2016blu's outbursts show irregular variability with multiple closely spaced peaks having typical amplitudes of 1-2 mag and durations of 1-4 weeks. While individual outbursts have irregular light curves, concentrations of these peaks recur with a period of $\sim 113 \pm 2$ d. Based on this period, we predict times for upcoming outbursts in 2023 and 2024. AT 2016blu shares similarities with SN 2000ch in NGC 3432, where outbursts may arise from periastron encounters in an eccentric binary containing a luminous blue variable (LBV). We propose that AT 2016blu's outbursts are also driven by interactions that intensify around periastron in an eccentric system. Intrinsic variability of the LBV-like primary star may cause different intensity and duration of binary interaction at each periastron passage. AT 2016blu also resembles the periastron encounters of $η$ Carinae prior to its Great Eruption and the erratic pre-SN eruptions of SN 2009ip. This similarity and the onset of eruptions in the past decade hint that AT 2016blu may also be headed for a catastrophe, making it a target of great interest. △ Less

Submitted 20 September, 2023; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: 18 pages, 14 figures, 6 tables, MNRAS Accepted

Radiation-driven acceleration in the expanding WR140 dust shell

Authors: Yinuo Han, Peter G. Tuthill, Ryan M. Lau, Anthony Soulain

Abstract: The Wolf-Rayet (WR) binary system WR140 is a close (0.9-16.7 mas) binary star consisting of an O5 primary and WC7 companion and is known as the archetype of episodic dust-producing WRs. Dust in WR binaries is known to form in a confined stream originating from the collision of the two stellar winds, with orbital motion of the binary sculpting the large-scale dust structure into arcs as dust is swe… ▽ More The Wolf-Rayet (WR) binary system WR140 is a close (0.9-16.7 mas) binary star consisting of an O5 primary and WC7 companion and is known as the archetype of episodic dust-producing WRs. Dust in WR binaries is known to form in a confined stream originating from the collision of the two stellar winds, with orbital motion of the binary sculpting the large-scale dust structure into arcs as dust is swept radially outwards. It is understood that sensitive conditions required for dust production in WR140 are only met around periastron when the two stars are sufficiently close. Here we present multiepoch imagery of the circumstellar dust shell of WR140. We constructed geometric models that closely trace the expansion of the intricately structured dust plume, showing that complex effects induced by orbital modulation may result in a 'Goldilocks zone' for dust production. We find that the expansion of the dust plume cannot be reproduced under the assumption of a simple uniform-speed outflow, finding instead the dust to be accelerating. This constitutes a direct kinematic record of dust motion under acceleration by radiation pressure and further highlights the complexity of the physical conditions in colliding-wind binaries. △ Less

Submitted 12 October, 2022; originally announced October 2022.

Comments: Published in Nature

Journal ref: Nature 610, 269-272 (2022)

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)

Smoke on the wind: dust nucleation in archetype colliding wind pinwheel WR104

Authors: A. Soulain, A. Lamberts, F. Millour, P. Tuthill, R. M. Lau

Abstract: A handful of binary Wolf-Rayet stars are known to harbour spectacular spiral structures spanning a few hundred AU. These systems host some of the highest dust production rates in the Universe and are therefore interesting candidates to address the origin of the enigmatic dust excess observed across galactic evolution. The substantial interaction between the winds of the Wolf-Rayet star and its com… ▽ More A handful of binary Wolf-Rayet stars are known to harbour spectacular spiral structures spanning a few hundred AU. These systems host some of the highest dust production rates in the Universe and are therefore interesting candidates to address the origin of the enigmatic dust excess observed across galactic evolution. The substantial interaction between the winds of the Wolf-Rayet star and its companion constitutes a unique laboratory to study the mechanisms of dust nucleation in a hostile environment. Using the grid-based $\texttt{RAMSES}$ code, we investigate this problem by performing a 3D hydrodynamic simulation of the inner region of the prototypical spiral nebula around WR104. We then process the $\texttt{RAMSES}$ results using the radiative transfer code $\texttt{RADMC3d}$ to generate a candidate observable scene. This allows us to estimate the geometrical parameters of the shocked region. We link those quantities to the specific chemical pathway for dust nucleation, where the hydrogen-rich companion's wind catalyses dust formation. The scaling laws we derive constitute a unique tool that can be directly compared to observations. Depending on the dust nucleation locus, the velocity field reveals a differential wind speed. Thus, the initial dust speed could be more balanced between the speeds of the two stellar winds ($\sim$1600 km/s). With $\texttt{RADMC3d}$, we provide constraints on the dust nucleation radius for different combinations of dust-to-gas ratio, hydrogen enrichment and dust grain properties. Finally, our models reveal that dust may escape beyond the boundaries of the spiral due to hydrodynamical instabilities in the wind collision zone. △ Less

Submitted 5 September, 2022; originally announced September 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

Detection of a broad 8 $μ$m UIR feature in the mid-infrared spectrum of WR 125 observed with Subaru/COMICS

Authors: Izumi Endo, Ryan M. Lau, Itsuki Sakon, Takashi Onaka, Peredur M. Williams, Victor I. Shenavrin

Abstract: We present the detection of a broad 8 $μ$m feature in newly formed dust around the carbon-rich Wolf-Rayet (WC) binary WR 125 from N-band low-resolution (NL; R$\sim$250) spectroscopy between 7.3-13.6 $μ$m and N-band (11.7 $μ$m) and Q-band (18.8 $μ$m) imaging with Subaru/COMICS in 2019 October. WR 125 is a colliding wind binary (${\rm WC7+O9}$) that exhibited renewed dust formation starting in 2018,… ▽ More We present the detection of a broad 8 $μ$m feature in newly formed dust around the carbon-rich Wolf-Rayet (WC) binary WR 125 from N-band low-resolution (NL; R$\sim$250) spectroscopy between 7.3-13.6 $μ$m and N-band (11.7 $μ$m) and Q-band (18.8 $μ$m) imaging with Subaru/COMICS in 2019 October. WR 125 is a colliding wind binary (${\rm WC7+O9}$) that exhibited renewed dust formation starting in 2018, $\sim$28 years after its first dust formation episode had been observed. We also compare our infrared photometry with historical observations and revise the dust-formation period of WR 125 to 28.1 years. Archival infrared spectra of five dusty WC stars, WR 48a, WR 98a, WR 104, WR 112 and WR 118, obtained with ISO/SWS are reanalyzed and compared with the WR 125 spectrum to search for a similar feature. We analyze the dusty WC spectra using two different extinction curves to investigate the impact of interstellar extinction correction on the presence and/or properties of the 8 $μ$m feature. All of the dusty WC spectra dereddened with the two different extinction curves show a broad feature around 8 $μ$m (FWHM$\sim$1-2 $μ$m). We suggest that these 8 $μ$m features seen in the dusty WC spectra are related to the Class C unidentified infrared (UIR) features. △ Less

Submitted 8 May, 2022; v1 submitted 27 April, 2022; originally announced April 2022.

Comments: 17 pages, 7 figures, 4 tables, accepted to 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

The Impact of Neutron Transfer Reactions on Heating and Cooling of Accreted Neutron Star Crusts

Authors: H. Schatz, Z. Meisel, E. F. Brown, S. S. Gupta, G. W. Hitt, W. R. Hix, R. Jain, R. Lau, P. Möller, W. -J. Ong, P. S. Shternin, Y. Xu, M. Wiescher

Abstract: Nuclear reactions heat and cool the crust of accreting neutron stars and need to be understood to interpret observations of X-ray bursts and of long-term cooling in transiently accreting systems. It was recently suggested that previously neglected neutron transfer reactions may play a significant role in the nuclear processes. We present results from full nuclear network calculations that now incl… ▽ More Nuclear reactions heat and cool the crust of accreting neutron stars and need to be understood to interpret observations of X-ray bursts and of long-term cooling in transiently accreting systems. It was recently suggested that previously neglected neutron transfer reactions may play a significant role in the nuclear processes. We present results from full nuclear network calculations that now include these reactions and determine their impact on crust composition, crust impurity, heating, and cooling. We find that a large number of neutron transfer reactions indeed occur and impact crust models. In particular, we identify a new type of reaction cycle that brings a pair of nuclei across the nuclear chart into equilibrium via alternating neutron capture and neutron release, interspersed with a neutron transfer. While neutron transfer reactions lead to changes in crust model predictions, and need to be considered in future studies, previous conclusions concerning heating, cooling, and compositional evolution are remarkably robust. △ Less

Submitted 21 December, 2021; originally announced December 2021.

Comments: 8 pages, 10 figures, accepted for Astrophysical Journal

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

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

AT 2019qyl in NGC 300: Internal Collisions in the Early Outflow from a Very Fast Nova in a Symbiotic Binary

Authors: Jacob E. Jencson, Jennifer E. Andrews, Howard E. Bond, Viraj Karambelkar, David J. Sand, Schuyler D. van Dyk, Nadejda Blagorodnova, Martha L. Boyer, Mansi M. Kasliwal, Ryan M. Lau, Shazrene Mohamed, Robert Williams, Patricia A. Whitelock, Rachael C. Amaro, K. Azalee Bostroem, Yize Dong, Michael J. Lundquist, Stefano Valenti, Samuel D. Wyatt, Jamie Burke, Kishalay De, Saurabh W. Jha, Joel Johansson, César Rojas-Bravo, David A. Coulter , et al. (17 additional authors not shown)

Abstract: Nova eruptions, thermonuclear explosions on the surfaces of white dwarfs (WDs), are now recognized to be among the most common shock-powered astrophysical transients. We present the early discovery and rapid ultraviolet (UV), optical, and infrared (IR) temporal development of AT 2019qyl, a recent nova in the nearby Sculptor Group galaxy NGC 300. The light curve shows a rapid rise lasting… ▽ More Nova eruptions, thermonuclear explosions on the surfaces of white dwarfs (WDs), are now recognized to be among the most common shock-powered astrophysical transients. We present the early discovery and rapid ultraviolet (UV), optical, and infrared (IR) temporal development of AT 2019qyl, a recent nova in the nearby Sculptor Group galaxy NGC 300. The light curve shows a rapid rise lasting $\lesssim 1$ day, reaching a peak absolute magnitude of $M_V = -9.2$ mag, and a very fast decline, fading by 2 mag over 3.5 days. A steep dropoff in the light curves after 71 days and the rapid decline timescale suggest a low-mass ejection from a massive WD with $M_{\rm WD} \gtrsim 1.2~M_{\odot}$. We present an unprecedented view of the early spectroscopic evolution of such an event. Three spectra prior to the peak reveal a complex, multicomponent outflow giving rise to internal collisions and shocks in the ejecta of an He/N-class nova. We identify a coincident IR-variable counterpart in the extensive preeruption coverage of the transient location and infer the presence of a symbiotic progenitor system with an O-rich asymptotic-giant-branch donor star, as well as evidence for an earlier UV-bright outburst in 2014. We suggest that AT 2019qyl is analogous to the subset of Galactic recurrent novae with red-giant companions such as RS Oph and other embedded nova systems like V407 Cyg. Our observations provide new evidence that internal shocks between multiple, distinct outflow components likely contribute to the generation of the shock-powered emission from such systems. △ Less

Submitted 5 September, 2021; v1 submitted 22 February, 2021; originally announced February 2021.

Comments: 32 pages, 10 figures, accepted in ApJ

The luminous red nova AT 2018bwo in NGC 45 and its binary yellow supergiant progenitor

Authors: Nadejda Blagorodnova, Jakub Klencki, Ondrej Pejcha, Paul M. Vreeswijk, Howard E. Bond, Kevin B. Burdge, Kishalay De, Christoffer Fremling, Robert D. Gehrz, Jacob E. Jencson, Mansi M. Kasliwal, Thomas Kupfer, Ryan M. Lau, Frank J. Masci, R. Michael Rich

Abstract: Luminous Red Novae (LRNe) are astrophysical transients associated with the partial ejection of a binary system's common envelope (CE) shortly before its merger. Here we present the results of our photometric and spectroscopic follow-up campaign of AT2018bwo (DLT18x), a LRN discovered in NGC45, and investigate its progenitor system using binary stellar-evolution models. The transient reached a peak… ▽ More Luminous Red Novae (LRNe) are astrophysical transients associated with the partial ejection of a binary system's common envelope (CE) shortly before its merger. Here we present the results of our photometric and spectroscopic follow-up campaign of AT2018bwo (DLT18x), a LRN discovered in NGC45, and investigate its progenitor system using binary stellar-evolution models. The transient reached a peak magnitude of $M_r=-10.97\pm0.11$ and maintained this brightness during its optical plateau of $t_p = 41\pm5$days. During this phase, it showed a rather stable photospheric temperature of ~3300K and a luminosity of ~$10^{40}$erg/s. The photosphere of AT2018bwo at early times appeared larger and cooler than other similar LRNe, likely due to an extended mass-loss episode before the merger. Towards the end of the plateau, optical spectra showed a reddened continuum with strong molecular absorption bands. The reprocessed emission by the cooling dust was also detected in the mid-infrared bands ~1.5 years after the outburst. Archival Spitzer and Hubble Space Telescope data taken 10-14 years before the transient event suggest a progenitor star with $T_{prog}\sim 6500$K, $R_{prog}\sim 100R_{\odot}$ and $L_{prog}\sim 2\times10^4L_{\odot}$, and an upper limit for optically thin warm (1000 K) dust mass of $M_d<10^{-6}M_{\odot}$. Using stellar binary-evolution models, we determined the properties of binary systems consistent with the progenitor parameter space. For AT2018bwo, we infer a primary mass of 12-16 $M_{\odot}$, which is 9-45% larger than the ~11$M_{\odot}$ obtained using single-star evolution models. The system, consistent with a yellow-supergiant primary, was likely in a stable mass-transfer regime with -2.4<log ($\dot{M}/M_{\odot}$/yr)<-1.2 a decade before the main instability occurred. During the dynamical merger, the system would have ejected 0.15-0.5$M_{\odot}$ with a velocity of ~500 km/s. △ Less

Submitted 27 January, 2022; v1 submitted 10 February, 2021; originally announced February 2021.

Comments: 21 pages, 16 figures, accepted version for A&A

Journal ref: A&A 653, A134 (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)

The extreme colliding-wind system Apep: resolved imagery of the central binary and dust plume in the infrared

Authors: Y. Han, P. G. Tuthill, R. M. Lau, A. Soulain, J. R. Callingham, P. M. Williams, P. A. Crowther, B. J. S. Pope, B. Marcote

Abstract: The recent discovery of a spectacular dust plume in the system 2XMM J160050.7-514245 (referred to as "Apep") suggested a physical origin in a colliding-wind binary by way of the "Pinwheel" mechanism. Observational data pointed to a hierarchical triple-star system, however several extreme and unexpected physical properties seem to defy the established physics of such objects. Most notably, a stark… ▽ More The recent discovery of a spectacular dust plume in the system 2XMM J160050.7-514245 (referred to as "Apep") suggested a physical origin in a colliding-wind binary by way of the "Pinwheel" mechanism. Observational data pointed to a hierarchical triple-star system, however several extreme and unexpected physical properties seem to defy the established physics of such objects. Most notably, a stark discrepancy was found in the observed outflow speed of the gas as measured spectroscopically in the line-of-sight direction compared to the proper motion expansion of the dust in the sky plane. This enigmatic behaviour arises at the wind base within the central Wolf-Rayet binary: a system that has so far remained spatially unresolved. Here we present an updated proper motion study deriving the expansion speed of Apep's dust plume over a two-year baseline that is four times slower than the spectroscopic wind speed, confirming and strengthening the previous finding. We also present the results from high-angular-resolution near-infrared imaging studies of the heart of the system, revealing a close binary with properties matching a Wolf-Rayet colliding-wind system. Based on these new observational constraints, an improved geometric model is presented yielding a close match to the data, constraining the orbital parameters of the Wolf-Rayet binary and lending further support to the anisotropic wind model. △ Less

Submitted 13 August, 2020; originally announced August 2020.

Comments: This article has been accepted for publication in the Monthly Notices of the Royal Astronomical Society. 17 pages, 9 figures, 6 tables

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

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

From SN2010da to NGC 300 ULX-1: Ten Years of Observations of an Unusual High Mass X-ray Binary in NGC 300

Authors: Breanna A. Binder, Stefania Carpano, Marianne Heida, Ryan Lau

Abstract: In 2010 May, an intermediate luminosity optical transient was discovered in the nearby galaxy NGC 300 by a South African amateur astronomer. In the decade since its discovery, multi-wavelength observations of the misnamed ``SN 2010da'' have continually re-shaped our understanding of this high mass X-ray binary system. In this review, we present an overview of the multi-wavelength observations and… ▽ More In 2010 May, an intermediate luminosity optical transient was discovered in the nearby galaxy NGC 300 by a South African amateur astronomer. In the decade since its discovery, multi-wavelength observations of the misnamed ``SN 2010da'' have continually re-shaped our understanding of this high mass X-ray binary system. In this review, we present an overview of the multi-wavelength observations and attempts to understand the 2010 transient event and, later, the re-classification of this system as NGC~300 ULX-1: a red supergiant + neutron star ultraluminous X-ray source. △ Less

Submitted 10 February, 2020; originally announced February 2020.

Comments: Accepted to special edition of Galaxies: High-Mass X-ray Binaries and Ultraluminous X-ray Sources Hosting Neutron Stars

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

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

Supernova 2014C: ongoing interaction with extended circumstellar material with silicate dust

Authors: Samaporn Tinyanont, Ryan M Lau, Mansi M Kasliwal, Keiichi Maeda, Nathan Smith, Ori D Fox, Robert D Gehrz, Kishalay De, Jacob Jencson, John Bally, Frank Masci

Abstract: Supernova (SN) 2014C is a unique explosion where a seemingly typical hydrogen-poor stripped envelope SN started to interact with a dense, hydrogen-rich circumstellar medium (CSM) a few months after the explosion. The delayed interaction suggests a detached CSM shell, unlike in a typical SN IIn where the CSM is much closer and the interaction commences earlier post-explosion; indicating a different… ▽ More Supernova (SN) 2014C is a unique explosion where a seemingly typical hydrogen-poor stripped envelope SN started to interact with a dense, hydrogen-rich circumstellar medium (CSM) a few months after the explosion. The delayed interaction suggests a detached CSM shell, unlike in a typical SN IIn where the CSM is much closer and the interaction commences earlier post-explosion; indicating a different mass loss history. We present near- to mid-infrared observations of SN 2014C from 1-5 years after the explosion, including uncommon 9.7 $μ$m imaging with COMICS on the Subaru telescope. Spectroscopy shows that the interaction is still ongoing, with the intermediate-width He I 1.083 $μ$m emission present out to our latest epoch 1639 days post-explosion. The last Spitzer/IRAC photometry at 1920 days post-explosion further confirms ongoing CSM interaction. The 1-10 $μ$m spectral energy distributions (SEDs) can be explained by a dust model with a mixture of 69% carbonaceous and 31% silicate dust, pointing to a chemically inhomogeneous CSM. The inference of silicate dust is the first among interacting SNe. An SED model with purely carbonaceous CSM dust is possible, but would require more than 0.22 $M_{\odot}$ of dust, which is an order of magnitude larger than what observed in any other SNe, measured in the same way, at this epoch. The light curve beyond 500 days is well fit by an interaction model with a wind-driven CSM and a mass loss rate of $\sim 10^{-3} \, M_{\odot}\,\rm yr^{-1}$, which presents an additional CSM density component exterior to the constant density shell reported previously in the literature. SN 2014C could originate in a binary system, similar to RY Scuti, which would explain the observed chemical and density profile inhomogeneity in the CSM. △ Less

Submitted 13 September, 2019; originally announced September 2019.

Comments: Submitted to ApJ

Discovery of a red supergiant donor star in SN2010da/NGC 300 ULX-1

Authors: M. Heida, R. M. Lau, B. Davies, M. Brightman, F. Fürst, B. W. Grefenstette, J. A. Kennea, F. Tramper, D. J. Walton, F. A. Harrison

Abstract: SN2010da/NGC 300 ULX-1 was first detected as a supernova impostor in May 2010 and was recently discovered to be a pulsating ultraluminous X-ray source. In this letter, we present VLT/X-shooter spectra of this source obtained in October 2018, covering the wavelength range 350-2300 nm. The $J$- and $H$-bands clearly show the presence of a red supergiant donor star that is best matched by a MARCS ste… ▽ More SN2010da/NGC 300 ULX-1 was first detected as a supernova impostor in May 2010 and was recently discovered to be a pulsating ultraluminous X-ray source. In this letter, we present VLT/X-shooter spectra of this source obtained in October 2018, covering the wavelength range 350-2300 nm. The $J$- and $H$-bands clearly show the presence of a red supergiant donor star that is best matched by a MARCS stellar atmosphere with $T_{\rm eff} = 3650 - 3900$ K and $\log(L_{\rm bol}/L_{\odot}) = 4.25\pm0.10$, which yields a stellar radius $R = 310 \pm 70 R_{\odot}$. To fit the full spectrum, two additional components are required: a blue excess that can be fitted either by a hot blackbody (T $\gtrsim 20,000$ K) or a power law (spectral index $α\approx 4$) and is likely due to X-ray emission reprocessed in the outer accretion disk or the donor star; and a red excess that is well fitted by a blackbody with a temperature of $\sim 1100$ K, and is likely due to warm dust in the vicinity of SN2010da. The presence of a red supergiant in this system implies an orbital period of at least 0.8-2.1 years, assuming Roche lobe overflow. Given the large donor-to-compact object mass ratio, orbital modulations of the radial velocity of the red supergiant are likely undetectable. However, the radial velocity amplitude of the neutron star is large enough (up to 40-60 km s$^{-1}$) to potentially be measured in the future, unless the system is viewed at a very unfavorable inclination. △ Less

Submitted 4 September, 2019; originally announced September 2019.

Comments: 10 pages, 6 figures, accepted for publication in ApJ Letters

Uncovering Red and Dusty Ultraluminous X-ray Sources with Spitzer

Authors: Ryan M. Lau, Marianne Heida, Dominic J. Walton, Mansi M. Kasliwal, Scott M. Adams, Ann Marie Cody, Kishalay De, Robert D. Gehrz, Felix Furst, Jacob E. Jencson, jamie A. Kennea, Frank Masci

Abstract: We present a mid-infrared (IR) sample study of nearby ultraluminous X-ray sources (ULXs) using multi-epoch observations with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. Spitzer/IRAC observations taken after 2014 were obtained as part of the Spitzer Infrared Intensive Transients Survey (SPIRITS). Our sample includes 96 ULXs located within 10 Mpc. Of the 96~ULXs, 12 have candida… ▽ More We present a mid-infrared (IR) sample study of nearby ultraluminous X-ray sources (ULXs) using multi-epoch observations with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. Spitzer/IRAC observations taken after 2014 were obtained as part of the Spitzer Infrared Intensive Transients Survey (SPIRITS). Our sample includes 96 ULXs located within 10 Mpc. Of the 96~ULXs, 12 have candidate counterparts consistent with absolute mid-IR magnitudes of supergiants, and 16 counterparts exceeded the mid-IR brightness of single supergiants and are thus more consistent with star clusters or non-ULX background active galactic nuclei (AGN). The supergiant candidate counterparts exhibit a bi-modal color distribution in a Spitzer/IRAC color-magnitude diagram, where "red" and "blue" ULXs fall in IRAC colors $[3.6] - [4.5]\sim0.7$ and $[3.6] - [4.5]\sim0.0$, respectively. The mid-IR colors and absolute magnitudes of 4 "red" and 5 "blue" ULXs are consistent with that of supergiant B[e] (sgB[e]) and red supergiant (RSG) stars, respectively. While "blue", RSG-like mid-IR ULX counterparts likely host RSG mass donors, we propose the "red" counterparts are ULXs exhibiting the "B[e] phenomenon'' rather than hosts of sgB[e] mass donors. We show that the mid-IR excess from the "red" ULXs is likely due to thermal emission from circumstellar or circumbinary dust. Using dust as a probe for total mass, we estimate mass-loss rates of $\dot{M}\sim1\times10^{-4}$ M$_\odot$ yr$^{-1}$ in dust-forming outflows of red ULXs. Based on the transient mid-IR behavior and its relatively flat spectral index, $α=-0.19\pm0.1$, we suggest that the mid-IR emission from Holmberg IX X-1 originates from a variable jet. △ Less

Submitted 18 April, 2019; originally announced April 2019.

Comments: 40 pages, 11 figures, accepted for publication in ApJ

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

SPIRITS Catalog of Infrared Variables: Identification of Extremely Luminous Long Period Variables

Authors: V. R. Karambelkar, S. M. Adams, P. A. Whitelock, M. M. Kasliwal, J. E. Jencson, M. L. Boyer, S. R. Goldman, F. Masci, A. M. Cody, J. Bally, H. E. Bond, R. D. Gehrz, M. Parthasarathy, R. M. Lau

Abstract: We present a catalog of 417 luminous infrared variable stars with periods exceeding 250 days. These were identified in 20 nearby galaxies by the ongoing SPIRITS survey with the Spitzer Space Telescope. Of these, 359 variables have $M_{[4.5]}$ (phase-weighted mean magnitudes) fainter than $-12$ and periods and luminosities consistent with previously reported variables in the Large Magellanic Cloud.… ▽ More We present a catalog of 417 luminous infrared variable stars with periods exceeding 250 days. These were identified in 20 nearby galaxies by the ongoing SPIRITS survey with the Spitzer Space Telescope. Of these, 359 variables have $M_{[4.5]}$ (phase-weighted mean magnitudes) fainter than $-12$ and periods and luminosities consistent with previously reported variables in the Large Magellanic Cloud. However, 58 variables are more luminous than $M_{[4.5]} = -12$, including 11 that are brighter than $M_{[4.5]} = -13$ with the brightest having $M_{[4.5]} = -15.51$. Most of these bright variable sources have quasi-periods longer than 1000 days, including four over 2000 days. We suggest that the fundamental period-luminosity relationship, previously measured for the Large Magellanic Cloud, extends to much higher luminosities and longer periods in this large galaxy sample. We posit that these variables include massive AGB stars (possibly super-AGB stars), red supergiants experiencing exceptionally high mass-loss rates, and interacting binaries. We also present 3.6, 4.5, 5.8 and 8.0 $μ$m photometric catalogs for all sources in these 20 galaxies. △ Less

Submitted 6 May, 2019; v1 submitted 22 January, 2019; originally announced January 2019.

Comments: 18 pages, 25 figures

Supernova 2017eaw: molecule and dust formation from infrared observations

Authors: Samaporn Tinyanont, Mansi M Kasliwal, Kelsie Krafton, Ryan Lau, Jeonghee Rho, Douglas C Leonard, Kishalay De, Jacob Jencson, Dimitri Mawet, Maxwell Millar-Blanchaer, Ricky Nilsson, Lin Yan, Robert D Gehrz, George Helou, Schuyler D Van Dyk, Eugene Serabyn, Ori D Fox, Geoffrey Clayton

Abstract: We present infrared (IR) photometry and spectroscopy of the Type II-P SN 2017eaw and its progenitor in the nearby galaxy NGC 6946. Progenitor observations in the Ks band in 4 epochs from 1 year to 1 day before the explosion reveal no significant variability in the progenitor star greater than 6% that last longer than 200 days. SN 2017eaw is a typical SN II-P with near-IR and mid-IR photometric evo… ▽ More We present infrared (IR) photometry and spectroscopy of the Type II-P SN 2017eaw and its progenitor in the nearby galaxy NGC 6946. Progenitor observations in the Ks band in 4 epochs from 1 year to 1 day before the explosion reveal no significant variability in the progenitor star greater than 6% that last longer than 200 days. SN 2017eaw is a typical SN II-P with near-IR and mid-IR photometric evolution similar to those of SNe 2002hh and 2004et, other normal SNe II-P in the same galaxy. Spectroscopic monitoring between 389 and 480 days post explosion reveals strong CO first overtone emission at 389 d, with a line profile matching that of SN 1987A from the same epoch, indicating $\sim 10^{-3} \, M_{\odot}$ of CO at 1,800 K. From the 389 d epoch until the most recent observation at 566 d, the first overtone feature fades while the 4.5 $μ$m excess, likely from the CO fundamental band, remains. This behavior indicates that the CO has not been destroyed, but that the gas has cooled enough that the levels responsible for first overtone emissions are no longer populated. Finally, the evolution of Spitzer 3.6 $μ$m photometry shows evidence for dust formation in SN 2017eaw, with a dust mass of $10^{-6}$ or $10^{-4}\,M_{\odot}$ assuming carbonaceous or silicate grains respectively. △ Less

Submitted 7 January, 2019; originally announced January 2019.