skip to main content


Search for: All records

Creators/Authors contains: "Lau, Ryan M."

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. 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 ≈1000 days and an amplitude of Δm≈ 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-IRJandKsbands 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 (MKs=10.7mag) and color (JKs= 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 (Teff=35001400+800K) and luminosity (logL/L=5.1±0.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 ofMinit= 17 ± 4M. We estimate the pre-supernova mass-loss rate of the star between 3 and 19 yr before explosion from the SED modeling atṀ3×105to 3 × 10−4Myr−1for an assumed wind velocity ofvw= 10 km s−1, perhaps pointing to enhanced mass loss in a pulsation-driven wind.

     
    more » « less
  2. We present optical and infrared (IR) light curves of the enshrouded massive binary NaSt1 (WR 122) with observations from Palomar Gattini-IR (PGIR), the Zwicky Transient Facility (ZTF), the Katzman Automatic Imaging Telescope (KAIT), and the All-Sky Automated Survey for Supernovae (ASAS-SN). The optical and IR light curves span between 2014 July and 2020 Oct., revealing periodic, sinusoidal variability from NaSt1 with a P=305.2±1.0 d period. We also present historical IR light curves taken between 1983 July and 1989 May that also indicate NaSt1 exhibits long-term IR variability on timescales of ∼decades. Fixed-period sinusoidal fits to the recent optical and IR light curves show that amplitude of NaSt1's variability is different at different wavelengths and also reveal significant phase offsets of ∼18 d between the ZTF r and PGIR J light curves.We interpret the ∼300 d period of the observed variability as the orbital period of a binary system in NaSt1. Assuming a circular orbit and adopting a range of combined stellar mass values in the range 20-100 M⊙ in NaSt1, we estimate orbital separations of ∼2-4 au. We suggest that the sinusoidal 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 provide an interpretation on the nature of NaSt1 and speculate that the mass-transfer process may have been triggered by Roche-lobe overflow (RLOF) during an eruptive phase of a Ofpe/WN9 star. Lastly, we claim that NaSt1 ceased RLOF mass transfer ≲3400 yr ago. 
    more » « less
  3. 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 to be in the range of ≈10–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 (≈5–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 the J band (≈1.2 μ m), which is significantly less affected by dust extinction compared to optical bands, the extinction distribution of the PGIR sample is highly skewed to a large extinction values (>50% of events obscured by A V ≳ 5 mag). Using recent estimates for the distribution of Galactic mass and dust, we show that the extinction distribution of the PGIR sample is commensurate with 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 ≈ 43.7 − 8.7 + 19.5 yr −1 . Our results suggest that all-sky near-infrared time-domain surveys are well poised to uncover the Galactic nova population. 
    more » « less
  4. null (Ed.)
  5. Abstract

    We report our Spitzer Space Telescope observations and detections of the binary neutron star merger GW170817. At 4.5 μm, GW170817 is detected at 21.9  mag AB at +43 days and 23.9 mag AB at +74 days after merger. At 3.6 μm, GW170817 is not detected to a limit of 23.2 mag AB at +43 days and 23.1 mag AB at  +74 days. Our detections constitute the latest and reddest constraints on the kilonova/macronova emission and composition of heavy elements. The 4.5 μm luminosity at this late phase cannot be explained by elements exclusively from the first abundance peak of the r-process. Moreover, the steep decline in the Spitzer band, with a power-law index of 3.4 ± 0.2, can be explained by a few of the heaviest isotopes with half-life around 14  d dominating the luminosity (e.g. 140Ba, 143Pr, 147Nd, 156Eu, 191Os, 223Ra, 225Ra, 233Pa, 234Th) or a model with lower deposition efficiency. This data offers evidence that the heaviest elements in the second and third r-process abundance peak were indeed synthesized. Our conclusion is verified by both analytics and network simulations and robust despite intricacies and uncertainties in the nuclear physics. Future observations with Spitzer and James Webb Space Telescope will further illuminate the relative abundance of the synthesized heavy elements.

     
    more » « less
  6. null (Ed.)