Search

While the subclass of interacting supernovae with narrow hydrogen emission lines (SNe IIn) consists of some of the longest-lasting and brightest SNe ever discovered, their progenitors are still not well understood. Investigating SNe IIn as they emit across the electromagnetic spectrum is the most robust way to understand the progenitor evolution before the explosion. This work presents X-Ray, optical, infrared, and radio observations of the strongly interacting Type IIn SN 2020ywx covering a period $>1200$ days after discovery. Through multiwavelength modeling, we find that the progenitor of 2020ywx was losing mass at $\sim10^{-2}$--$10^{-3} \mathrm{\,M_{\odot}\,yr^{-1}}$ for at least 100 yr pre-explosion using the circumstellar medium (CSM) speed of $120$ km/s measured from our optical and NIR spectra. Despite the similar magnitude of mass-loss measured in different wavelength ranges, we find discrepancies between the X-ray and optical/radio-derived mass-loss evolution, which suggest asymmetries in the CSM. Furthermore, we find evidence for dust formation due to the combination of a growing blueshift in optical emission lines and near-infrared continuum emission which we fit with blackbodies at $\sim$ 1000 K. Based on the observed elevated mass loss over more than 100 years and the configuration of the CSM inferred from the multiwavelength observations, we invoke binary interaction as the most plausible mechanism to explain the overall mass-loss evolution. SN 2020ywx is thus a case that may support the growing observational consensus that SNe IIn mass loss is explained by binary interaction., Comment: Submitted to ApJ, 31 pages, 19 figures

The Legacy Survey of Space and Time (LSST) at Vera C. Rubin Observatory is planned to begin in the Fall of 2025. The LSST survey cadence has been designed via a community-driven process regulated by the Survey Cadence Optimization Committee (SCOC), which recommended up to 3% of the observing time to carry out Target of Opportunity (ToO) observations. Experts from the scientific community, Rubin Observatory personnel, and members of the SCOC were brought together to deliver a recommendation for the implementation of the ToO program during a workshop held in March 2024. Four main science cases were identified: gravitational wave multi-messenger astronomy, high energy neutrinos, Galactic supernovae, and small potentially hazardous asteroids possible impactors. Additional science cases were identified and briefly addressed in the documents, including lensed or poorly localized gamma-ray bursts and twilight discoveries. Trigger prioritization, automated response, and detailed strategies were discussed for each science case. This document represents the outcome of the Rubin ToO 2024 workshop, with additional contributions from members of the Rubin Science Collaborations. The implementation of the selection criteria and strategies presented in this document has been endorsed in the SCOC Phase 3 Recommendations document (PSTN-056). Although the ToO program is still to be finalized, this document serves as a baseline plan for ToO observations with the Rubin Observatory.

We present a comprehensive multi-epoch photometric and spectroscopic study of SN 2024bch, a nearby (19.9 Mpc) Type II supernova (SN) with prominent early high ionization emission lines. Optical spectra from 2.9 days after the estimated explosion reveal narrow lines of H I, He II, C IV, and N IV that disappear by day 6. High cadence photometry from the ground and TESS show that the SN brightened quickly and reached a peak M$_V \sim$ $-$17.8 mag within a week of explosion, and late-time photometry suggests a $^{56}$Ni mass of 0.050 M$_{\odot}$. High-resolution spectra from day 8 and 43 trace the unshocked circumstellar medium (CSM) and indicate a wind velocity of 30--40 km s$^{-1}$, a value consistent with a red supergiant (RSG) progenitor. Comparisons between models and the early spectra suggest a pre-SN mass-loss rate of $\dot{M} \sim 10^{-3}-10^{-2}\ M_\odot\ \mathrm{yr}^{-1}$, which is too high to be explained by quiescent mass loss from RSGs, but is consistent with some recent measurements of similar SNe. Persistent blueshifted H I and [O I] emission lines seen in the optical and NIR spectra could be produced by asymmetries in the SN ejecta, while the multi-component H$\alpha$ may indicate continued interaction with an asymmetric CSM well into the nebular phase. SN 2024bch provides another clue to the complex environments and mass-loss histories around massive stars., Comment: Accepted in ApJ 2024 Dec 30

We present supernova (SN) 2023ufx, a unique Type IIP SN with the shortest known plateau duration ($t_\mathrm{PT}$ $\sim$47 days), a luminous V-band peak ($M_{V}$ = $-$18.42 $\pm$ 0.08 mag), and a rapid early decline rate ($s1$ = 3.47 $\pm$ 0.09 mag (50 days)$^{-1}$). By comparing observed photometry to a hydrodynamic MESA+STELLA model grid, we constrain the progenitor to be a massive red supergiant with M$_\mathrm{ZAMS}$ $\simeq$19 - 25 M$_{\odot}$. Independent comparisons with nebular spectral models also suggest an initial He-core mass of $\sim$6 M$_{\odot}$, and thus a massive progenitor. For a Type IIP, SN 2023ufx produced an unusually high amount of nickel ($^{56}$Ni) $\sim$0.14 $\pm$ 0.02 M$_{\odot}$, during the explosion. We find that the short plateau duration in SN 2023ufx can be explained with the presence of a small hydrogen envelope (M$_\mathrm{H_\mathrm{env}}$ $\simeq$1.2 M$_{\odot}$), suggesting partial stripping of the progenitor. About $\simeq$0.09 M$_{\odot}$ of CSM through mass loss from late-time stellar evolution of the progenitor is needed to fit the early time ($\lesssim$10 days) pseudo-bolometric light curve. Nebular line diagnostics of broad and multi-peak components of [O I] $\lambda\lambda$6300, 6364, H$\alpha$, and [Ca II] $\lambda \lambda$7291, 7323 suggest that the explosion of SN 2023ufx could be inherently asymmetric, preferentially ejecting material along our line-of-sight., Comment: Submitted to ApJ, 30 pages, 19 figures

We present multi-epoch optical spectropolarimetric and imaging polarimetric observations of the nearby Type II supernova (SN) 2023ixf discovered in M101 at a distance of 6.85 Mpc. The first imaging polarimetric observations were taken +2.33 days (60085.08 MJD) after the explosion, while the last imaging polarimetric data points (+73.19 and +76.19 days) were acquired after the fall from the light curve plateau. At +2.33 days there is strong evidence of circumstellar material (CSM) interaction in the spectra and the light curve. A significant level of polarization $P_r = 0.88\pm 0.06 \% $ seen during this phase indicates that this CSM is aspherical. We find that the polarization evolves with time toward the interstellar polarization level ($0.35\%$) during the photospheric phase, which suggests that the recombination photosphere is spherically symmetric. There is a jump in polarization ($P_r =0.65 \pm 0.08 \% $) at +73.19 days when the light curve falls from the plateau. This is a phase where polarimetric data is sensitive to non-spherical inner ejecta or a decrease in optical depth into the single scattering regime. We also present spectropolarimetric data that reveal line (de)polarization during most of the observed epochs. In addition, at +14.50 days we see an "inverse P Cygn" profile in the H and He line polarization, which clearly indicates the presence of asymmetrically distributed material overlying the photosphere. The overall temporal evolution of polarization is typical for Type II SNe, but the high level of polarization during the rising phase has only been observed in SN 2023ixf., Comment: 14 pages, 7 figures, submitted to ApJL, comments welcome

We present and analyze the extensive optical broadband photometry of the Type II SN 2023ixf up to one year after explosion. We find that, when compared to two pre-existing model grids, the pseudo-bolometric light curve is consistent with drastically different combinations of progenitor and explosion properties. This may be an effect of known degeneracies in Type IIP light-curve models. We independently compute a large grid of ${\tt MESA+STELLA}$ single-star progenitor and light-curve models with various zero-age main-sequence masses, mass-loss efficiencies, and convective efficiencies. Using the observed progenitor variability as an additional constraint, we select stellar models consistent with the pulsation period and explode them according to previously established scaling laws to match plateau properties. Our hydrodynamic modeling indicates that SN 2023ixf is most consistent with a moderate-energy ($E_{\rm exp}\approx7\times10^{50}$ erg) explosion of an initially high-mass red supergiant progenitor ($\gtrsim 17\ M_{\odot}$) that lost a significant amount of mass in its prior evolution, leaving a low-mass hydrogen envelope ($\lesssim 3\ M_{\odot}$) at the time of explosion, with a radius $\gtrsim 950\ R_{\odot}$ and a synthesized $^{56}$Ni mass of $0.07\ M_{\odot}$. We posit that previous mass transfer in a binary system may have stripped the envelope of SN 2023ixf's progenitor. The analysis method with pulsation period presented in this work offers a way to break degeneracies in light-curve modeling in the future, particularly with the upcoming Vera C.~Rubin Observatory Legacy Survey of Space and Time, when a record of progenitor variability will be more common., Comment: 18 pages, 7 figures, submitted to ApJ. Comments welcome

SN 2023ixf was discovered in M101 within a day of explosion and rapidly classified as a Type II supernova with flash features. Here we present ultraviolet (UV) spectra obtained with the Hubble Space Telescope 14, 19, 24, and 66 days after explosion. Interaction between the supernova ejecta and circumstellar material (CSM) is seen in the UV throughout our observations in the flux of the first three epochs and asymmetric Mg II emission on day 66. We compare our observations to CMFGEN supernova models that include CSM interaction ($\dot{M}<10^{-3}$ M$_{\odot}$ yr$^{-1}$) and find that the power from CSM interaction is decreasing with time, from $L_{\rm sh}\approx5\times10^{42}$ erg s$^{-1}$ to $L_{\rm sh}\approx1\times10^{40}$ erg s$^{-1}$ between days 14 and 66. We examine the contribution of individual atomic species to the spectra on days 14 and 19, showing that the majority of the features are dominated by iron, nickel, magnesium, and chromium absorption in the ejecta. The UV spectral energy distribution of SN 2023ixf sits between that of supernovae which show no definitive signs of CSM interaction and those with persistent signatures assuming the same progenitor radius and metallicity. Finally, we show that the evolution and asymmetric shape of the Mg II $\lambda\lambda$ 2796, 2802 emission are not unique to SN 2023ixf. These observations add to the early measurements of dense, confined CSM interaction, tracing the mass-loss history of SN 2023ixf to $\sim33$ yr prior to the explosion and the density profile to a radius of $\sim5.7\times10^{15}$ cm. They show the relatively short evolution from a quiescent red supergiant wind to high mass loss., Comment: Accepted ApJL

We present multiwavelength observations of the Swift short $\gamma$-ray burst GRB 231117A, localized to an underlying galaxy at redshift $z = 0.257$ at a small projected offset ($\sim 2~$kpc). We uncover long-lived X-ray (Chandra) and radio/millimeter (VLA, MeerKAT, and ALMA) afterglow emission, detected to $\sim 37~$days and $\sim 20~$days (rest frame), respectively. We measure a wide jet ($\sim 10.4^\circ$) and relatively high circumburst density ($\sim 0.07~{\rm cm}^{-3}$) compared to the short GRB population. Our data cannot be easily fit with a standard forward shock model, but they are generally well fit with the incorporation of a refreshed forward shock and a reverse shock at $< 1~$day. We incorporate GRB 231117A into a larger sample of 132 X-ray detected events, 71 of which were radio-observed (17 cm-band detections), for a systematic study of the distributions of redshifts, jet and afterglow properties, galactocentric offsets, and local environments of events with and without detected radio afterglows. Compared to the entire short GRB population, the majority of radio-detected GRBs are at relatively low redshifts ($z < 0.6$) and have high circumburst densities ($> 10^{-2}~{\rm cm}^{-3}$), consistent with their smaller ($< 8~$kpc) projected galactocentric offsets. We additionally find that 70% of short GRBs with opening angle measurements were radio-detected, indicating the importance of radio afterglows in jet measurements, especially in the cases of wide ($> 10^\circ$) jets where observational evidence of collimation may only be detectable at radio wavelengths. Owing to improved observing strategies and the emergence of sensitive radio facilities, the number of radio-detected short GRBs has quadrupled in the past decade., Comment: 30 pages, 11 figures, submitted to ApJ

Variability can be the pathway to understanding the physical processes in astrophysical jets, however, the high-cadence observations required to test particle acceleration models are still missing. Here we report on the first attempt to produce continuous, >24 hour polarization light curves of blazars using telescopes distributed across the globe and the rotation of the Earth to avoid the rising Sun. Our campaign involved 16 telescopes in Asia, Europe, and North America. We observed BL Lacertae and CGRaBS J0211+1051 for a combined 685 telescope hours. We find large variations in the polarization degree and angle for both sources in sub-hour timescales as well as a ~180 degree rotation of the polarization angle in CGRaBS J0211+1051 in less than two days. We compared our high-cadence observations to Particle-In-Cell magnetic reconnection and turbulent plasma simulations. We find that although the state of the art simulation frameworks can produce a large fraction of the polarization properties, they do not account for the entirety of the observed polarization behavior in blazar jets., Comment: 20 pages, 15 figures, 2 tables, accepted for publication in A&A. The data used in the paper are available here: https://doi.org/10.7910/DVN/IETSXS

We present high-cadence photometric and spectroscopic observations of supernova (SN) 2024ggi, a Type II SN with flash spectroscopy features which exploded in the nearby galaxy NGC 3621 at $\sim$7 Mpc. The light-curve evolution over the first 30 hours can be fit by two power law indices with a break after 22 hours, rising from $M_V \approx -12.95$ mag at +0.66 days to $M_V \approx -17.91$ mag after 7 days. In addition, the densely sampled color curve shows a strong blueward evolution over the first few days and then behaves as a normal SN II with a redward evolution as the ejecta cool. Such deviations could be due to interaction with circumstellar material (CSM). Early high- and low-resolution spectra clearly show high-ionization flash features from the first spectrum to +3.42 days after the explosion. From the high-resolution spectra, we calculate the CSM velocity to be 37 $\pm~4~\mathrm{km\,s^{-1}} $. We also see the line strength evolve rapidly from 1.22 to 1.49 days in the earliest high-resolution spectra. Comparison of the low-resolution spectra with CMFGEN models suggests that the pre-explosion mass-loss rate of SN 2024ggi falls in a range of $10^{-3}$ to $10^{-2}$ M$_{\odot}$ yr$^{-1}$, which is similar to that derived for SN 2023ixf. However, the rapid temporal evolution of the narrow lines in the spectra of SN 2024ggi ($R_\mathrm{CSM} \sim 2.7 \times 10^{14} \mathrm{cm}$) could indicate a smaller spatial extent of the CSM than in SN 2023ixf ($R_\mathrm{CSM} \sim 5.4 \times 10^{14} \mathrm{cm}$) which in turn implies lower total CSM mass for SN 2024ggi., Comment: 23 pages, 15 figures, 4 tables, accepted for publication in ApJL

We present photometric and spectroscopic observations of SN 2023fyq, a type Ibn supernova in the nearby galaxy NGC 4388 (D$\simeq$18~Mpc). In addition, we trace long-standing precursor emission at the position of SN 2023fyq using data from DLT40, ATLAS, ZTF, ASAS-SN, Swift, and amateur astronomer Koichi Itagaki. Precursor activity is observed up to nearly three years before the supernova explosion, with a relatively rapid rise in the final 100 days. The double-peaked post-explosion light curve reaches a luminosity of $\sim10^{43}~\rm erg\,s^{-1}$. The strong intermediate-width He lines observed in the nebular spectrum of SN 2023fyq imply the interaction is still active at late phases. We found that the precursor activity in SN 2023fyq is best explained by the mass transfer in a binary system involving a low-mass He star and a compact companion. An equatorial disk is likely formed in this process ($\sim$0.6$\rm M_{\odot}$), and the interaction of SN ejecta with this disk powers the main peak of the supernova. The early SN light curve reveals the presence of dense extended material ($\sim$0.3$\rm M_{\odot}$) at $\sim$3000$\rm R_{\odot}$ ejected weeks before the SN explosion, likely due to final-stage core silicon burning or runaway mass transfer resulting from binary orbital shrinking, leading to rapid rising precursor emission within $\sim$30 days prior to explosion. The final explosion could be triggered either by the core-collapse of the He star or by the merger of the He star with a compact object. SN 2023fyq, along with SN 2018gjx and SN 2015G, forms a unique class of Type Ibn SNe which originate in binary systems and are likely to exhibit detectable long-lasting pre-explosion outbursts with magnitudes ranging from $-$10 to $-$13., Comment: submitted to ApJ

We present comprehensive optical observations of SN~2021gmj, a Type II supernova (SN~II) discovered within a day of explosion by the Distance Less Than 40~Mpc (DLT40) survey. Follow-up observations show that SN~2021gmj is a low-luminosity SN~II (LL~SN~II), with a peak magnitude $M_V = -15.45$ and Fe~II velocity of $\sim 1800 \ \mathrm{km} \ \mathrm{s}^{-1}$ at 50 days past explosion. Using the expanding photosphere method, we derive a distance of $17.8^{+0.6}_{-0.4}$~Mpc. From the tail of the light curve we obtain a radioactive nickel mass of $0.014 \pm 0.001$ M$_{\odot}$. The presence of circumstellar material (CSM) is suggested by the early-time light curve, early spectra, and high-velocity H$\alpha$ in absorption. Analytical shock-cooling models of the light curve cannot reproduce the fast rise, supporting the idea that the early-time emission is partially powered by the interaction of the SN ejecta and CSM. The inferred low CSM mass of 0.025 M$_{\odot}$ in our hydrodynamic-modeling light curve analysis is also consistent with our spectroscopy. We observe a broad feature near 4600 \AA, which may be high-ionization lines of C, N, or/and He~II. This feature is reproduced by radiation-hydrodynamic simulations of red supergiants with extended atmospheres. Several LL~SNe~II show similar spectral features, implying that high-density material around the progenitor may be common among them., Comment: Accepted version at ApJ

We present high cadence optical and ultraviolet observations of the Type II supernova (SN), SN 2022jox which exhibits early spectroscopic high ionization flash features of \ion{H}{1}, \ion{He}{2}, \ion{C}{4}, and \ion{N}{4} that disappear within the first few days after explosion. SN 2022jox was discovered by the Distance Less than 40 Mpc (DLT40) survey $\sim$0.75 days after explosion with followup spectra and UV photometry obtained within minutes of discovery. The SN reached a peak brightness of M$_V \sim$ $-$17.3 mag, and has an estimated $^{56}$Ni mass of 0.04 M$_{\odot}$, typical values for normal Type II SNe. The modeling of the early lightcurve and the strong flash signatures present in the optical spectra indicate interaction with circumstellar material (CSM) created from a progenitor with a mass loss rate of $\dot{M} \sim 10^{-3}-10^{-2}\ M_\odot\ \mathrm{yr}^{-1}$. There may also be some indication of late-time CSM interaction in the form of an emission line blueward of H$\alpha$ seen in spectra around 200 days. The mass-loss rate is much higher than the values typically associated with quiescent mass loss from red supergiants, the known progenitors of Type II SNe, but is comparable to inferred values from similar core collapse SNe with flash features, suggesting an eruptive event or a superwind in the progenitor in the months or years before explosion., Comment: ApJ, accepted 2024 Feb 14

We present upgraded infrastructure for Searches after Gravitational Waves Using ARizona Observatories (SAGUARO) during LIGO, Virgo, and KAGRA's fourth gravitational-wave (GW) observing run (O4). These upgrades implement many of the lessons we learned after a comprehensive analysis of potential electromagnetic counterparts to the GWs discovered during the previous observing run. We have developed a new web-based target and observation manager (TOM) that allows us to coordinate sky surveys, vet potential counterparts, and trigger follow-up observations from one centralized portal. The TOM includes software that aggregates all publicly available information on the light curves and possible host galaxies of targets, allowing us to rule out potential contaminants like active galactic nuclei, variable stars, solar-system objects, and preexisting supernovae, as well as to assess the viability of any plausible counterparts. We have also upgraded our image-subtraction pipeline by assembling deeper reference images and training a new neural network-based real-bogus classifier. These infrastructure upgrades will aid coordination by enabling the prompt reporting of observations, discoveries, and analysis to the GW follow-up community, and put SAGUARO in an advantageous position to discover kilonovae in the remainder of O4 and beyond. Many elements of our open-source software stack have broad utility beyond multimessenger astronomy, and will be particularly relevant in the "big data" era of transient discoveries by the Vera C. Rubin Observatory., Comment: updated to match accepted version

We present high-cadence photometric and spectroscopic observations of SN~2023axu, a classical Type II supernova with an absolute $V$-band peak magnitude of $-16.5 \pm 0.1$ mag. SN~2023axu was discovered by the Distance Less Than 40 Mpc (DLT40) survey within 1 day of the last non-detection in the nearby galaxy NGC 2283 at 13.7 Mpc. We modeled the early light curve using a recently updated shock cooling model that includes the effects of line blanketing and found the explosion epoch to be MJD 59971.48 $\pm$ 0.03 and the probable progenitor to be a red supergiant with a radius of 417 $\pm$ 28 $R_\odot$. The shock cooling model cannot match the rise of observed data in the $r$ and $i$ bands and underpredicts the overall UV data which points to possible interaction with circumstellar material. This interpretation is further supported by spectral behavior. We see a ledge feature around 4600 \AA\ in the very early spectra (+1.1 and +1.5 days after the explosion) which can be a sign of circumstellar interaction. The signs of circumstellar material are further bolstered by the presence of absorption features blueward of H$\alpha$ and H$\beta$ at day $>$40 which is also generally attributed to circumstellar interaction. Our analysis shows the need for high-cadence early photometric and spectroscopic data to decipher the mass-loss history of the progenitor., Comment: 18 pages, 12 figures, to be submitted to the AAS Journals

We present optical, infrared, ultraviolet, and radio observations of SN 2022xkq, an underluminous fast-declining type Ia supernova (SN Ia) in NGC 1784 ($\mathrm{D}\approx31$ Mpc), from $<1$ to 180 days after explosion. The high-cadence observations of SN 2022xkq, a photometrically transitional and spectroscopically 91bg-like SN Ia, cover the first days and weeks following explosion which are critical to distinguishing between explosion scenarios. The early light curve of SN 2022xkq has a red early color and exhibits a flux excess which is more prominent in redder bands; this is the first time such a feature has been seen in a transitional/91bg-like SN Ia. We also present 92 optical and 19 near-infrared (NIR) spectra, beginning 0.4 days after explosion in the optical and 2.6 days after explosion in the NIR. SN 2022xkq exhibits a long-lived C I 1.0693 $\mu$m feature which persists until 5 days post-maximum. We also detect C II $\lambda$6580 in the pre-maximum optical spectra. These lines are evidence for unburnt carbon that is difficult to reconcile with the double detonation of a sub-Chandrasekhar mass white dwarf. No existing explosion model can fully explain the photometric and spectroscopic dataset of SN 2022xkq, but the considerable breadth of the observations is ideal for furthering our understanding of the processes which produce faint SNe Ia., Comment: 38 pages, 16 figures, accepted for publication in ApJ, the figure 15 input models and synthetic spectra are now available at https://zenodo.org/record/8379254

We present optical and near-infrared observations of SN~2022crv, a stripped envelope supernova in NGC~3054, discovered within 12 hrs of explosion by the Distance Less Than 40 Mpc Survey. We suggest SN~2022crv is a transitional object on the continuum between SNe Ib and SNe IIb. A high-velocity hydrogen feature ($\sim$$-$20,000 -- $-$16,000 $\rm km\,s^{-1}$) was conspicuous in SN~2022crv at early phases, and then quickly disappeared around maximum light. By comparing with hydrodynamic modeling, we find that a hydrogen envelope of $\sim 10^{-3}$ \msun{} can reproduce the behaviour of the hydrogen feature observed in SN~2022crv. The early light curve of SN~2022crv did not show envelope cooling emission, implying that SN~2022crv had a compact progenitor with extremely low amount of hydrogen. The analysis of the nebular spectra shows that SN~2022crv is consistent with the explosion of a He star with a final mass of $\sim$4.5 -- 5.6 \msun{} that has evolved from a $\sim$16 -- 22 \msun{} zero-age main sequence star in a binary system with about 1.0 -- 1.7 \msun{} of oxygen finally synthesized in the core. The high metallicity at the supernova site indicates that the progenitor experienced a strong stellar wind mass loss. In order to retain a small amount of residual hydrogen at such a high metallicity, the initial orbital separation of the binary system is likely larger than $\sim$1000~$\rm R_{\odot}$. The near-infrared spectra of SN~2022crv show a unique absorption feature on the blue side of He I line at $\sim$1.005~$\mu$m. This is the first time that such a feature has been observed in a Type Ib/IIb, and could be due to \ion{Sr}{2}. Further detailed modelling on SN~2022crv can shed light on the progenitor and the origin of the mysterious absorption feature in the near infrared., Comment: accepted for publication in ApJ

Nebular-phase observations of peculiar Type Ia supernovae (SNe Ia) provide important constraints on progenitor scenarios and explosion dynamics for both these rare SNe and the more common, cosmologically useful SNe Ia. We present observations from an extensive ground-based and space-based follow-up campaign to characterize SN 2022pul, a "super-Chandrasekhar" mass SN Ia (alternatively "03fg-like" SN), from before peak brightness to well into the nebular phase across optical to mid-infrared (MIR) wavelengths. The early rise of the light curve is atypical, exhibiting two distinct components, consistent with SN Ia ejecta interacting with dense carbon-oxygen rich circumstellar material (CSM). In the optical, SN 2022pul is most similar to SN 2012dn, having a low estimated peak luminosity ($M_{B}=-18.9$ mag) and high photospheric velocity relative to other 03fg-like SNe. In the nebular phase, SN 2022pul adds to the increasing diversity of the 03fg-like subclass. From 168 to 336 days after peak $B$-band brightness, SN 2022pul exhibits asymmetric and narrow emission from [O I] $\lambda\lambda 6300,\ 6364$ (${\rm FWHM} \approx 2{,}000$ km s$^{-1}$), strong, broad emission from [Ca II] $\lambda\lambda 7291,\ 7323$ (${\rm FWHM} \approx 7{,}300$ km s$^{-1}$), and a rapid Fe III to Fe II ionization change. Finally, we present the first-ever optical-to-mid-infrared (MIR) nebular spectrum of an 03fg-like SN Ia using data from JWST. In the MIR, strong lines of neon and argon, weak emission from stable nickel, and strong thermal dust emission (with $T \approx 500$ K), combined with prominent [O I] in the optical, suggest that SN 2022pul was produced by a white dwarf merger within carbon/oxygen-rich CSM., Comment: 23 pages, 11 figures, submitted to ApJ

We present an analysis of ground-based and JWST observations of SN~2022pul, a peculiar "03fg-like" (or "super-Chandrasekhar") Type Ia supernova (SN Ia), in the nebular phase at 338d post explosion. Our combined spectrum continuously covers 0.4--14 $\mu$m and includes the first mid-infrared spectrum of an 03fg-like SN Ia. Compared to normal SN Ia 2021aefx, SN 2022pul exhibits a lower mean ionization state, asymmetric emission-line profiles, stronger emission from the intermediate-mass elements (IMEs) argon and calcium, weaker emission from iron-group elements (IGEs), and the first unambiguous detection of neon in a SN Ia. Strong, broad, centrally peaked [Ne II] line at 12.81 $\mu$m was previously predicted as a hallmark of "violent merger'' SN Ia models, where dynamical interaction between two sub-$M_{ch}$ white dwarfs (WDs) causes disruption of the lower mass WD and detonation of the other. The violent merger scenario was already a leading hypothesis for 03fg-like SNe Ia; in SN 2022pul it can explain the large-scale ejecta asymmetries seen between the IMEs and IGEs and the central location of narrow oxygen and broad neon. We modify extant models to add clumping of the ejecta to better reproduce the optical iron emission, and add mass in the innermost region ($< 2000$ km s$^{-1}$) to account for the observed narrow [O I]~$\lambda\lambda6300$, 6364 emission. A violent WD-WD merger explains many of the observations of SN 2022pul, and our results favor this model interpretation for the subclass of 03fg-like SN Ia., Comment: 20 pages, 10 figures, published in ApJ

We perform a comprehensive search for optical precursor emission at the position of SN~2023ixf using data from the DLT40, ZTF and ATLAS surveys. By comparing the current data set with precursor outburst hydrodynamical model light curves, we find that the probability of a significant outburst within five years of explosion is low, and the circumstellar material (CSM) ejected during any possible precursor outburst is likely smaller than $\sim$0.015\msun. By comparing to a set of toy models, we find that, if there was a precursor outburst, the duration must have been shorter than $\sim$100 days for a typical brightness of $M_{r}\simeq-9$ mag or shorter than 200 days for $M_{r}\simeq-8$ mag; brighter, longer outbursts would have been discovered. Precursor activity like that observed in the normal type II SN~2020tlf ($M_{r}\simeq-11.5$) can be excluded in SN~2023ixf. If the dense CSM inferred by early flash spectroscopy and other studies is related to one or more precursor outbursts, then our observations indicate that any such outburst would have to be faint and only last for days to months, or it occurred more than five years prior to the explosion. Alternatively, any dense, confined CSM may not be due to eruptive mass loss from a single red supergiant (RSG) progenitor. Taken together, the results of SN~2023ixf and SN~2020tlf indicate that there may be more than one physical mechanism behind the dense CSM inferred around some normal type II SNe., Comment: 10 pages, 4 figures

The semi-regular variable star V Canum Venaticorum (V CVn) is well-known for its unusual linear polarization position angle (PA). Decades of observing V CVn reveal a nearly constant PA spanning hundreds of pulsation cycles. This phenomenon has persisted through variability that has ranged by 2 magnitudes in optical brightness and through variability in the polarization amplitude over 0.3% and 6.9%. Additionally, the polarization fraction of V CVn varies inversely with brightness. This paper presents polarization measurements obtained over three pulsation cycles. We find that the polarization maximum does not always occur precisely at the same time as the brightness minimum. Instead, we observe a small lead or lag in relation to the brightness minimum, spanning a period of a few days up to three weeks. Furthermore, the PA sometimes exhibits a non-negligible rotation, especially at lower polarization levels. To elucidate the unusual optical behavior of V CVn, we present a list of literature sources that also exhibit polarization variability with a roughly fixed PA. We find this correlation occurs in stars with high tangential space velocities, i.e., "runaway" stars, suggesting that the long-term constant PA is related to how the circumstellar gas is shaped by the star's high-speed motion through the interstellar medium., Comment: 9 pages + appendices, 7 figures, accepted for publication in A&A

We present the optical spectroscopic evolution of SN~2023ixf seen in sub-night cadence spectra from 1.18 to 14 days after explosion. We identify high-ionization emission features, signatures of interaction with material surrounding the progenitor star, that fade over the first 7 days, with rapid evolution between spectra observed within the same night. We compare the emission lines present and their relative strength to those of other supernovae with early interaction, finding a close match to SN~2020pni and SN~2017ahn in the first spectrum and SN~2014G at later epochs. To physically interpret our observations we compare them to CMFGEN models with confined, dense circumstellar material around a red supergiant progenitor from the literature. We find that very few models reproduce the blended \NC{} emission lines observed in the first few spectra and their rapid disappearance thereafter, making this a unique diagnostic. From the best models, we find a mass-loss rate of $10^{-3}-10^{-2}$ \mlunit{}, which far exceeds the mass-loss rate for any steady wind, especially for a red supergiant in the initial mass range of the detected progenitor. These mass-loss rates are, however, similar to rates inferred for other supernovae with early circumstellar interaction. Using the phase when the narrow emission features disappear, we calculate an outer dense radius of circumstellar material $R_\mathrm{CSM, out}\sim5\times10^{14}~\mathrm{cm}$ and a mean circumstellar material density of $\rho=5.6\times10^{-14}~\mathrm{g\,cm^{-3}}$. This is consistent with the lower limit on the outer radius of the circumstellar material we calculate from the peak \Halpha{} emission flux, $R_\text{CSM, out}\gtrsim9\times10^{13}~\mathrm{cm}$., Comment: Published in ApJL

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 $\Delta 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., Comment: 13 pages, 5 figures, published in ApJL, replacement with revisions to match published version

We present a series of high-resolution echelle spectra of SN~2023ixf in M101, obtained nightly during the first week or so after discovery using PEPSI on the LBT. NaID absorption in these spectra indicates a host reddening of $E(B-V)$=0.031~mag and a systemic velocity of $+$7~km~s$^{-1}$ relative to the average redshift of M101. Dramatic changes are seen in in the strength and shape of strong emission lines emitted by CSM, including HeII4686, CIV5801,5811, H$\alpha$, and NIV7109,7123. In general, these narrow lines broaden to become intermediate-width lines before disappearing from the spectrum within a few days, indicating a limited extent to the dense CSM of around 20-30 AU (or $\la$10$^{14.7}$ cm). H$\alpha$ persists in the spectrum for about a week as an intermediate-width emission line with P~Cyg absorption at 700-1300 km s$^{-1}$ arising in the post-shock shell of swept-up CSM. Early narrow emission lines are blueshifted and indicate an expansion speed in the pre-shock CSM of about 115 km s$^{-1}$, but with even broader emission in higher ionization lines. This is faster than the normal winds of red supergiants, suggesting some mode of eruptive mass loss from the progenitor or radiative acceleration of the CSM. A lack of narrow blueshifted absorption suggests that most of the CSM is not along our line of sight. This and several other clues indicate that the CSM of SN~2023ixf is significantly aspherical. We find that CSM lines disappear after a few days because the asymmetric CSM is engulfed by the SN photosphere., Comment: Accepted verion, 8/14

We present the densely sampled early light curve of the Type II supernova (SN) 2023ixf, first observed within hours of explosion in the nearby Pinwheel Galaxy (Messier 101; 6.7 Mpc). Comparing these data to recently updated models of shock-cooling emission, we find that the progenitor likely had a radius of $410 \pm 10\ R_\odot$. Our estimate is model dependent but consistent with a red supergiant. These models provide a good fit to the data starting about 1 day after the explosion, despite the fact that the classification spectrum shows signatures of circumstellar material around SN 2023ixf during that time. Photometry during the first day after the explosion, provided almost entirely by amateur astronomers, does not agree with the shock-cooling models or a simple power-law rise fit to data after 1 day. We consider the possible causes of this discrepancy, including precursor activity from the progenitor star, circumstellar interaction, and emission from the shock before or after it breaks out of the stellar surface. The very low luminosity ($-11\mathrm{\ mag} > M > -14\mathrm{\ mag}$) and short duration of the initial excess lead us to prefer a scenario related to prolonged emission from the SN shock traveling through the progenitor system., Comment: updated to match accepted version

We present very early photometric and spectroscopic observations of the Type Ia supernova (SN Ia) 2023bee, starting about 8 hr after the explosion, which reveal a strong excess in the optical and nearest UV (U and UVW1) bands during the first several days of explosion. This data set allows us to probe the nature of the binary companion of the exploding white dwarf and the conditions leading to its ignition. We find a good match to the Kasen model in which a main-sequence companion star stings the ejecta with a shock as they buzz past. Models of double detonations, shells of radioactive nickel near the surface, interaction with circumstellar material, and pulsational delayed detonations do not provide good matches to our light curves. We also observe signatures of unburned material, in the form of carbon absorption, in our earliest spectra. Our radio nondetections place a limit on the mass-loss rate from the putative companion that rules out a red giant but allows a main-sequence star. We discuss our results in the context of other similar SNe Ia in the literature., Comment: updated to match accepted version

We present five far- and near-ultraviolet spectra of the Type II plateau supernova, SN 2022acko, obtained 5, 6, 7, 19, and 21 days after explosion, all observed with the Hubble Space Telescope/Space Telescope Imaging Spectrograph. The first three epochs are earlier than any Type II plateau supernova has been observed in the far-ultraviolet revealing unprecedented characteristics. These three spectra are dominated by strong lines, primarily from metals, which contrasts with the relatively featureless early optical spectra. The flux decreases over the initial time series as the ejecta cools and line-blanketing takes effect. We model this unique dataset with the non-local thermodynamic equilibrium radiation transport code CMFGEN, finding a good match to the explosion of a low mass red supergiant with energy Ekin = 6 x 10^50 erg. With these models we identify, for the first time, the ions that dominate the early UV spectra. We also present optical photometry and spectroscopy, showing that SN 2022acko has a peak absolute magnitude of V = -15.4 mag and plateau length of ~115d. The spectra closely resemble those of SN 2005cs and SN 2012A. Using the combined optical and UV spectra, we report the fraction of flux redwards of the uvw2, U, B, and V filters on days 5, 7, and 19. We also create a spectral time-series of Type II supernovae in the ultraviolet, demonstrating the rapid decline of UV flux over the first few weeks of evolution. Future observations of Type II supernovae will continue to explore the diversity seen in the limited set of high-quality UV spectra., Comment: Published in ApJL

A merger of binary neutron stars creates heavy unstable elements whose radioactive decay produces a thermal emission known as a kilonova. In this paper, we predict the photometric and polarimetric behaviour of this emission by performing 3-D Monte Carlo radiative transfer simulations. In particular, we choose three hydrodynamical models for merger ejecta, two including jets with different luminosities and one without a jet structure, to help decipher the impact of jets on the light curve and polarimetric behaviour. In terms of photometry, we find distinct color evolutions across the three models. Models without a jet show the highest variation in light curves for different viewing angles. In contrast, to previous studies, we find models with a jet to produce fainter kilonovae when viewed from orientations close to the jet axis, compared to a model without a jet. In terms of polarimetry, we predict relatively low levels (<~0.3-0.4%) at all orientations that, however, remain non-negligible until a few days after the merger and longer than previously found. Despite the low levels, we find that the presence of a jet enhances the degree of polarization at wavelengths ranging from 0.25 to 2.5\micron, an effect that is found to increase with the jet luminosity. Thus, future photometric and polarimetric campaigns should observe kilonovae in blue and red filters for a few days after the merger to help constrain the properties of the ejecta (e.g. composition) and jet., Comment: 11 pages, 7 figures, accepted for publication and in press

We present photometric and spectroscopic observations of the extraordinary gamma-ray burst (GRB) 221009A in search of an associated supernova. Some past GRBs have shown bumps in the optical light curve that coincide with the emergence of supernova spectral features, but we do not detect any significant light curve features in GRB~221009A, nor do we detect any clear sign of supernova spectral features. Using two well-studied GRB-associated supernovae (SN~2013dx, $M_{r,max} = -19.54$; SN~2016jca, $M_{r,max} = -19.04$) at a similar redshift as GRB~221009A ($z=0.151$), we modeled how the emergence of a supernova would affect the light curve. If we assume the GRB afterglow to decay at the same rate as the X-ray data, the combination of afterglow and a supernova component is fainter than the observed GRB brightness. For the case where we assume the best-fit power law to the optical data as the GRB afterglow component, a supernova contribution should have created a clear bump in the light curve, assuming only extinction from the Milky Way. If we assume a higher extinction of $E(B-V)$=$1.74$ mag (as has been suggested elsewhere), the supernova contribution would have been hard to detect, with a limit on the associated supernova of $M_{r,max} \approx-$19.54. We do not observe any clear supernova features in our spectra, which were taken around the time of expected maximum light. The lack of a bright supernova associated with GRB~221009A may indicate that the energy from the explosion is mostly concentrated in the jet, leaving a lower energy budget available for the supernova., Comment: 19 pages, 8 figures, 4 tables, accepted for publication in ApJL

We report on analysis using the James Webb Space Telescope (JWST) to identify a candidate progenitor star of the Type II-plateau supernova SN 2022acko in the nearby, barred spiral galaxy NGC 1300. To our knowledge, our discovery represents the first time JWST has been used to localize a progenitor system in pre-explosion archival Hubble Space Telescope (HST) images. We astrometrically registered a JWST NIRCam image from 2023 January, in which the SN was serendipitously captured, to pre-SN HST F160W and F814W images from 2017 and 2004, respectively. An object corresponding precisely to the SN position has been isolated with reasonable confidence. That object has a spectral energy distribution (SED) and overall luminosity consistent with a single-star model having an initial mass possibly somewhat less than the canonical 8 Msun theoretical threshold for core collapse (although masses as high as 9 Msun for the star are also possible); however, the star's SED and luminosity are inconsistent with that of a super-asymptotic giant branch star which might be a forerunner of an electron-capture SN. The properties of the progenitor alone imply that SN 2022acko is a relatively normal SN II-P, albeit most likely a low-luminosity one. The progenitor candidate should be confirmed with follow-up HST imaging at late times, when the SN has sufficiently faded. This potential use of JWST opens a new era of identifying SN progenitor candidates at high spatial resolution., Comment: 8 pages, substantial changes from v1, to appear in MNRAS

The sunlight reflected from the Moon during a total lunar eclipse has been transmitted through the Earth's atmosphere on the way to the Moon. The combination of multiple scattering and inhomogeneous atmospheric characteristics during that transmission can potentially polarize that light. A similar (although much smaller) effect should also be observable from the atmosphere of a transiting exoplanet. We present the results of polarization observations during the first 15 minutes of totality of the lunar eclipse of 2022 May 16. We find degrees of polarization of 2.1 +/- 0.4 per cent in B, 1.2 +/- 0.3 per cent in V, 0.5 +/- 0.2 per cent in R and 0.2 +/- 0.2 per cent in I. Our polarization values lie in the middle of the range of those reported for previous eclipses, providing further evidence that the induced polarization can change from event to event. We found no significant polarization difference (<0.02 per cent) between a region of dark Mare and nearby bright uplands or between the lunar limb and regions closer to the disk centre due to the different angle of incidence. This further strengthens the interpretation of the polarization's origin being due to scattering in the Earth's atmosphere rather than by the lunar regolith., Comment: Accepted for publication in MNRAS

Massive star binaries are critical laboratories for measuring masses and stellar wind mass-loss rates. A major challenge is inferring viewing inclination and extracting information about the colliding wind interaction (CWI) region. Polarimetric variability from electron scattering in the highly ionized winds provides important diagnostic information about system geometry. We combine for the first time the well-known generalized treatment of \citet{brown_polarisation_1978} for variable polarization from binaries with the semi-analytic solution for the geometry and surface density CWI shock interface between the winds based on Canto et al 1996. Our calculations include some simplifications in the form of inverse square-law wind densities and the assumption of axisymmetry, but in so doing arrive at several robust conclusions. One is that when the winds are nearly equal (e.g., O\,+\,O binaries), the polarization has a relatively mild decline with binary separation. Another is that despite Thomson scattering being a gray opacity, the continuum polarization can show chromatic effects at ultraviolet wavelengths but will be mostly constant at longer wavelengths. Finally, when one wind dominates the other, as for example in WR+OB binaries, the polarization is expected to be larger at wavelengths where the OB component is more luminous, and generally smaller at wavelengths where the WR component is more luminous. This behavior arises because from the perspective of the WR star, the distortion of the scattering envelope from spherical is a minor perturbation situated far from the WR star. By contrast, the polarization contribution from the OB star is dominated by the geometry of the CWI shock., Comment: accepted for publication in the Astrophysical Journal. arXiv admin note: text overlap with arXiv:2111.11552

Understanding the evolution of massive binary stars requires accurate estimates of their masses. This understanding is critically important because massive star evolution can potentially lead to gravitational wave sources such as binary black holes or neutron stars. For Wolf-Rayet stars with optically thick stellar winds, their masses can only be determined with accurate inclination angle estimates from binary systems which have spectroscopic $M \sin i$ measurements. Orbitally-phased polarization signals can encode the inclination angle of binary systems, where the Wolf-Rayet winds act as scattering regions. We investigated four Wolf-Rayet + O star binary systems, WR 42, WR 79, WR 127, and WR 153, with publicly available phased polarization data to estimate their masses. To avoid the biases present in analytic models of polarization while retaining computational expediency, we used a Monte Carlo radiative transfer model accurately emulated by a neural network. We used the emulated model to investigate the posterior distribution of parameters of our four systems. Our mass estimates calculated from the estimated inclination angles put strong constraints on existing mass estimates for three of the systems, and disagrees with the existing mass estimates for WR 153. We recommend a concerted effort to obtain polarization observations that can be used to estimate the masses of Wolf-Rayet binary systems and increase our understanding of their evolutionary paths., Comment: 24 pages, 13 figures, accepted to ApJ

The observational signatures of black holes in x-ray binary systems depend on their masses, spins, accretion rate and the misalignment angle between the black hole spin and the orbital angular momentum. We present optical polarimetric observations of the black hole x-ray binary MAXI J1820+070, from which we constrain the position angle of the binary orbital axis. Combining this with previous determinations of the relativistic jet orientation axis, which traces the black hole spin, and the inclination of the orbit, we determine a lower limit of 40 deg on the spin-orbit misalignment angle. The misalignment has to originate from either the binary or black hole formation stage. If other x-ray binaries have similarly large misalignments, these would bias measurements of black hole masses and spins from x-ray observations., Comment: 31 pages, 9 figures

We study the polarization produced by scattering from dust in a bow shock-shaped region of enhanced density surrounding a stellar source, using the Monte Carlo radiative transfer code SLIP. Bow shocks are structures formed by the interaction of the winds of fast-moving stars with the interstellar medium. Our previous study focused on the polarization produced in these structures by electron scattering; we showed that polarization is highly dependent on inclination angle and that multiple scattering changes the shape and degree of polarization. In contrast to electron scattering, dust scattering is wavelength-dependent, which changes the polarization behaviour. Here we explore different dust particle sizes and compositions and generate polarized spectral energy distributions for each case. We find that the polarization SED behaviour depends on the dust composition and grain size. Including dust emission leads to polarization changes with temperature at higher optical depth in ways that are sensitive to the orientation of the bow shock. In various scenarios and under certain assumptions, our simulations can constrain the optical depth and dust properties of resolved and unresolved bow shock-shaped scattering regions.Constraints on optical depth can provide estimates of local ISM density for observed bow shocks. We also study the impact of dust grains filling the region between the star and bow shock. We see that as the density of dust between the star and bow shock increases, the resulting polarization is suppressed for all the optical depth regimes., Comment: 21 pages, accepted for publication in MNRAS

Polarization plays an important role in various time-domain astrophysics to understand the magnetic fields, geometry, and environments of spatially unresolved variable sources. In this paper we present the results of laboratory and on-sky testing of a novel dual-beam, dual-camera optical imaging polarimeter (MOPTOP) exploiting high sensitivity, low-noise CMOS technology and designed to monitor variable and transient sources with low systematic errors and high sensitivity. We present a data reduction algorithm that corrects for sensitivity variations between the cameras on a source-by-source basis. Using our data reduction algorithm, we show that our dual-beam, dual-camera technique delivers the benefits of low and stable instrumental polarization ($<0.05$\% for lab data and $<0.25$\% for on sky data) and high throughput while avoiding the additional sky brightness and image overlap problems associated with dual-beam, single-camera polarimeters., Comment: 12 pages, 11 figures, Accepted for publication in MNRAS

Bow shocks and related density enhancements produced by the winds of massive stars moving through the interstellar medium provide important information regarding the motions of the stars, the properties of their stellar winds, and the characteristics of the local medium. Since bow shocks are aspherical structures, light scattering within them produces a net polarization signal even if the region is spatially unresolved. Scattering opacity arising from free electrons and dust leads to a distribution of polarized intensity across the bow shock structure. That polarization encodes information about the shape, composition, opacity, density, and ionisation state of the material within the structure. In this paper we use the Monte Carlo radiative transfer code SLIP to investigate the polarization created when photons scatter in a bow shock-shaped region of enhanced density surrounding a stellar source. We present results assuming electron scattering, and investigate the polarization behaviour as a function of optical depth, temperature, and source of photons for two different cases: pure scattering and scattering with absorption. In both regimes we consider resolved and unresolved cases. We discuss the implication of these results as well as their possible use along with observational data to constrain the properties of observed bow shock systems. In different situations and under certain assumptions, our simulations can constrain viewing angle, optical depth and temperature of the scattering region, and the relative luminosities of the star and shock., Comment: 20 pages, 13 figures, 1 table; Accepted for publication in MNRAS