We present a crossmatch between a combined catalog of X-ray sources and the Vera C. Rubin Observatory Data Preview 1 (DP1) to identify optical counterparts. The six fields targeted as part of DP1 include the Extended Chandra Deep Field South (E-CDF-S), the Euclid Deep Field South, the Fornax Dwarf Spheroidal Galaxy (Fornax dSph), 47 Tucanae (47 Tuc), and science validation fields with low Galactic and ecliptic latitude (SV_95_-25 and SV_38_7, respectively). We find matches to 2314 of 3830 X-ray sources. We also compare our crossmatch to DP1 in the E-CDF-S field to previous efforts to identify optical counterparts. The probability of a chance coincidence match varies across each DP1 field, with overall high reliability in the E-CDF-S field, and a lower proportion of high-reliability matches in the other fields. The majority of previously known sources that we detect are, unsurprisingly, active galaxies. We plot the X-ray-to-optical flux ratio against optical magnitude and color in an effort to identify Galactic accreting compact objects using a Gaia color threshold transformed to LSST g - i, but do not find any strong candidates in these primarily extragalactic counterparts. The DP1 dataset contains high-cadence photometry collected over a number of nights. We calculate the Stetson J variability index for each object under the hypothesis that X-ray counterparts tend to exhibit higher optical variability; however, the evidence is inconclusive whether our sample is more variable over DP1 timescales when compared to field objects.

The Zwicky Transient Facility Census of the Local Universe survey yielded a sample of 330 Type IIP supernovae (SNe) with well-constrained peak luminosities. In paper I, we measured their luminosity function and volumetric rate. Here (paper II), we present the largest systematic study of lightcurve properties for Type IIP SNe from a volume-limited survey, analyzing a selected subset of 129 events, including 16 low-luminosity Type IIP (LLIIP) SNe with Mr,peak >= -16 mag. We find that plateau slope correlates with peak brightness, with many LLIIP SNe showing positive slopes-suggesting smaller progenitor radii and distinct density profiles compared to brighter Type IIP SNe. The plateau duration shows only a weak dependence on peak brightness, likely suggesting binary interaction. One SN exhibits a plateau-to-tail drop of >3.5 mag, consistent with an electron-capture or failed SN with very low or zero nickel mass. We derive explosion and progenitor parameters of the entire Type IIP SN sample using semi-analytical and radiation-hydrodynamical models. Based on radiation-hydrodynamical model fitting, LLIIP SNe are characterized by low nickel masses (0.001-0.025 M-circle dot), low explosion energies (0.1-0.28 & times; 10(51 )erg), low ejecta masses ( 8.1(-1.7)(+0.8) M-circle dot), and ZAMS masses below 11 M-circle dot. In comparison, the full Type IIP SN sample spans a wider range with nickel masses (0.001-0.222 M-circle dot), explosion energies (0.10-4.43 & times; 10(51 )erg), ejecta masses (5.4-24.8 M-circle dot), and ZAMS masses (9.3-16.7 M-circle dot). We find strong correlations between peak brightness, explosion energy, and nickel mass that extend to the low-luminosity end. We conclude that LLIIP SNe represent the faint, low-energy end of the Type IIP population and originate from the lowest-mass core-collapse progenitors.

Time-domain and multimessenger astronomy (MMA/TDA) targets demand rapid-response follow-up observations. In many cases, it is the only way to make discoveries and advance our understanding of the astrophysical phenomena, for example, kilonovae accompanying gravitational waves from compact object mergers, shock breakout in supernovae, prompt emission from GRBs, etc. Presently the MMA/TDA follow-up workflow requires wrangling disparate software packages and user interfaces. We present an end-to-end software tool for the community, the Gemini Observation and Analysis of Targets System (GOATS), which unifies and simplifies the workflow, particularly for Gemini follow-up observations. GOATS achieves this by integrating services from Gemini Observatory and its parent organization, NSF NOIRLab. From a single platform, GOATS enables enhanced target selection via NOIRLab's ANTARES alert broker, triggering of Gemini (and other facilities within the Astronomical Event Observatory Network), automated data retrieval from the Gemini Observatory Archive, and interactive data reduction and analysis through Gemini's DRAGONS software and NOIRLab's Astro Data Lab science platform. GOATS was successfully deployed in an end-to-end demonstration of real-time follow-up of Rubin/LSST alerts with NOIRLab facilities. As part of this demonstration, we selected targets from the Rubin alert stream and triggered follow-up observations within minutes of the Rubin detections. We obtained spectra for several targets and classified them as supernova of various types (Ia, IIP, Ib/c) with redshifts ranging from 0.05 to 0.35. By eliminating the need to manually connect tools and automating repetitive tasks, GOATS lowers the entry barrier and allows users to focus on the scientific interpretation of the observation results.

We present the searches conducted with the Zwicky Transient Facility (ZTF) in response to S250206dm, a bona fide event with an online false alarm rate of one in 25 yr, detected by the International Gravitational Wave Network. Although the event is significant, the nature of the compact objects involved remains unclear, with at least one likely neutron star. ZTF covered 68% of the last refined Bilby localization region, though we did not identify any likely optical counterpart. We describe the ZTF strategy, potential candidates, and the observations that helped rule out candidates, including sources circulated by other collaborations. Similar to Ahumada et al., we perform a frequentist analysis, using simsurvey, as well as Bayesian analysis, using nimbus, to quantify the efficiency of our searches. We find that, given the nominal up-to-date distance to this event of 373 +/- 104 Mpc, our efficiencies are above 10% for KNe brighter than -17.5 absolute magnitude. Assuming the optical counterpart known as kilonova (KN) lies within the ZTF footprint, our limits constrain the brightest end of the KN parameter space. Through dedicated radiative transfer simulations of KNe from binary neutron star (BNS) and black hole-neutron star mergers, we exclude parts of the BNS KN parameter space. Up to 35% of the models with high wind ejecta mass (Mwind approximate to 0.13 M circle dot) are ruled out when viewed face-on ( cos theta obs=1.0 ). Finally, we present a joint analysis using the combined coverage from ZTF and the Gravitational Wave Multimessenger Dark Energy Camera Survey. The joint observations cover 73% of the Bilby localization region, and the combined efficiency has a stronger impact on rising and slowly fading models, allowing us to rule out 55% of the high-mass KN models viewed face-on.

The Taurid Resonant Swarm (TRS) within the Taurid Complex hosts dynamically-concentrated debris in a 7:2 mean-motion resonance with Jupiter. Fireball observations have confirmed that the TRS is rich in sub-meter-sized particles, but whether this enhancement extends to larger, asteroid-sized objects remains unclear. Here we reanalyze the data obtained by a Zwicky Transient Facility (ZTF) campaign during the 2022 TRS encounter, and find that the TRS may host up to similar to 102 Tunguska-sized objects and up to similar to 103 Chelyabinsk-sized objects, the latter of which agrees the estimate derived from bolide records. This translates to an impact frequency of less than once every 4 million years. However, we caution that these numbers are based on the unverified assumption that the orbital distribution of the TRS asteroids follows that of fireball-sized meteoroids. Future wide-field facilities, such as the Vera C. Rubin Observatory, could take advantage of TRS's close approaches in the 2020-30s and validate the constraints of the asteroid-sized objects in the TRS.

Many white dwarfs are observed in compact double white dwarf binaries, and through the emission of gravitational waves, a large fraction are destined to merge. The merger remnants that do not explode in a Type Ia supernova are expected to initially be rapidly rotating and highly magnetized. In this work, we present our discovery of the variable white dwarf ZTF J200832.79+444939.67, hereafter ZTF J2008+4449, as a likely merger remnant showing signs of circumstellar material without a stellar or substellar companion. The nature of ZTF J2008+4449 as a merger remnant is supported by its physical properties: it is hot (35 500 ± 300 K) and massive (1.12 ± 0.03 M⊙), rapidly rotating with a period of ≈6.6 minutes, and likely possesses exceptionally strong magnetic fields (∼400−600 MG) at its surface. Remarkably, we detect a significant period derivative of (1.80 ± 0.09)×10−12 s/s, indicating that the white dwarf is spinning down, and a soft X-ray emission that is inconsistent with photospheric emission. As the presence of a mass-transferring stellar or brown dwarf companion is excluded by infrared photometry, the detected spin-down and X-ray emission could be tell-tale signs of a magnetically driven wind or of interaction with circumstellar material, possibly originating from the fallback of gravitationally bound merger ejecta or from the tidal disruption of a planetary object. We also detect Balmer emission, which requires the presence of ionized hydrogen in the vicinity of the white dwarf, showing Doppler shifts as high as ≈2000 km s−1. The unusual variability of the Balmer emission on the spin period of the white dwarf is consistent with the trapping of a half ring of ionized gas in the magnetosphere of the white dwarf.

We report on the observation and measurement of astrometry, photometry, morphology, and activity of the interstellar object 3I/ATLAS, also designated C/2025 N1 (ATLAS) with the NSF-DOE Vera C. Rubin Observatory. Comet 3I/ATLAS, the third known interstellar object, was discovered on UT 2025 July 1. Rubin Observatory had coincidentally collected images of the object's region of the sky during routine commissioning. Facilitated by Rubin's high resolution and large aperture, we successfully recovered object detections from Rubin observations spanning UT 2025 June 21 (10 days before discovery, when 3I/ATLAS was 4.5 au from the Sun) through the date of discovery, and we acquired additional images through UT 2025 July 20 as part of commissioning. We measure on-sky locations of 3I/ATLAS in Rubin ugrizy bands, with a typical precision of similar to 70 mas, and briefly describe the reason this is coarser than our measured static source astrometric precision of similar to 3 mas in Rubin images. We measure grizy magnitudes of 3I/ATLAS photometry at similar to 0.01 mag precision, detecting no short-term photometric variability above 0.01 mag. We derive an estimated near-nucleus dust-to-nucleus scattering cross-sectional ratio of eta greater than or similar to 13 on UT 2025 July 2 based on Rubin photometry and an upper limit nucleus size computed from Hubble Space Telescope observations. We find Rubin colors of g - r = (0.657 +/- 0.013) mag, r - i = (0.235 +/- 0.018) mag, i - z = (0.147 +/- 0.042) mag, and z - y = (0.047 +/- 0.052) mag. These data represent the earliest observations of this object by a large (greater than or similar to 8 m class) telescope and illustrate the type of measurements (and discoveries) Rubin's Legacy Survey of Space and Time will provide after it begins in early 2026.

Machine learning methods are well established in the classification of quasars (QSOs). However, the advent of light-curve observations adds a great amount of complexity to the problem. Our goal is to use the Zwicky Transient Facility (ZTF) to create a catalog of QSOs. We process the ZTF DR20 light curves with a transformer artificial neural network and combine different surveys with extreme gradient boosting. Based on ZTF g -band and Wide-field Infrared Survey Explorer (WISE) observations, we find 4,849,574 objects classified as QSOs with confidence higher than 90% (QZO). We robustly classify objects fainter than the 5 σ signal-to-noise ratio (SNR) limit at g = 20.8 by requiring g < n _obs /80 + 20.375. For 33% of QZO objects, with available WISE data, we publish redshifts with estimated error Δ z /(1 + z ) = 0.14. We find that ZTF classification is superior to the Pan-STARRS static bands, and on par with WISE and Gaia measurements, but the light curves provide the most important features for QSO classification in the ZTF data set. Using ZTF g -band data with at least 100 observational epochs per light curve, we obtain a 97% F1 score for QSOs. We find that with 3 day median cadence, a survey time span of at least 900 days is required to achieve a 90% QSO F1 score. However, one can obtain the same score with a survey time span of 1800 days and the median cadence prolonged to 12 days.

We present new rotational period estimates for 216 Jupiter Trojans using photometric data from the Zwicky Transient Facility (ZTF), including 80 Trojans with previously unknown periods. Our analysis reveals rotation periods ranging from 4.6 hours to 447.8 hours. These results support the existence of a spin barrier for Trojans larger than 10 km, with periods clustering between 4 and 4.8 hours. This spin barrier is roughly twice as long as that observed for main-belt asteroids, suggesting that Jupiter Trojans have significantly lower bulk densities, likely due to a higher fraction of ices and volatile materials in their composition. We identify three new Trojans with reliable rotation periods near the spin barrier, doubling the number of known Trojans in this critical period range. Using these results, we estimate a mean density of approximately 0.52 g/cm^3 for rubble-pile Trojans. Our findings support the growing evidence that many Trojans are rubble-pile bodies with distinct physical properties compared to main-belt asteroids. Looking forward, we anticipate that data from the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will provide rotational period estimates for several hundred thousand Trojans, down to objects as small as 1 km, enabling a more detailed investigation of their rotational properties and internal structure.

Tidal disruption events (TDEs) occur when a star enters the tidal radius of a supermassive black hole (SMBH). If the star only grazes the tidal radius, a fraction of the stellar mass will be accreted in a partial TDE (pTDE). The remainder can continue orbiting and may re-disrupted at pericenter, causing a repeating pTDE. pTDEs may be as or more common than full TDEs (fTDEs), yet few are known. In this work, we present the discovery of the first repeating pTDE from a systematically-selected sample, AT\,2020vdq. AT\,2020vdq was originally identified as an optically- and radio-flaring TDE. Around $3$ years after its discovery, it rebrightened dramatically and rapidly in the optical. The optical flare was remarkably fast and luminous compared to previous TDEs. It was accompanied by extremely broad (${\sim}0.1c$) optical/UV spectral features and faint X-ray emission ($L_X \sim 3\times10^{41}$\,erg\,s$^{-1}$), but no new radio-emitting component. Based on the transient optical/UV spectral features and the broadband light curve, we show that AT\,2020vdq is a repeating pTDE. We then use it to constrain TDE models; in particular, we favor a star originally in a very tight binary system that is tidally broken apart by the Hills mechanism. We also constrain the repeating pTDE rate to be $10^{-6}$ to $10^{-5}$ yr$^{-1}$ galaxy$^{-1}$, with uncertainties dominated by the unknown distribution of pTDE repeat timescales. In the Hills framework, this means the binary fraction in the galactic nucleus is of the order few percent.

Context. Gravitationally lensed type Ia supernovae (glSNe Ia) are unique astronomical tools that can be used to study cosmological parameters, distributions of dark matter, the astrophysics of the supernovae, and the intervening lensing galaxies themselves. A small number of highly magnified glSNe Ia have been discovered by ground-based telescopes such as the Zwicky Transient Facility (ZTF), but simulations predict that a fainter, undetected population may also exist. Aims. We present a systematic search for glSNe Ia in the ZTF archive of alerts distributed from June 1 2019 to September 1 2022. Methods. Using the AMPEL platform, we developed a pipeline that distinguishes candidate glSNe Ia from other variable sources. Initial cuts were applied to the ZTF alert photometry (with constraints on the peak absolute magnitude and the distance to a catalogue-matched galaxy, as examples) before forced photometry was obtained for the remaining candidates. Additional cuts were applied to refine the candidates based on their light curve colours, lens galaxy colours, and the resulting parameters from fits to the SALT2 SN Ia template. The candidates were also cross-matched with the DESI spectroscopic catalogue. Results. Seven transients were identified that passed all the cuts and had an associated galaxy DESI redshift, which we present as glSN Ia candidates. Although superluminous supernovae (SLSNe) cannot be fully rejected as contaminants, two events, ZTF19abpjicm and ZTF22aahmovu, are significantly different from typical SLSNe and their light curves can be modelled as two-image glSN Ia systems. From this two-image modelling, we estimate time delays of 22 +/- 3 and 34 +/- 1 days for the two events, respectively, which suggests that we have uncovered a population of glSNe Ia with longer time delays. Conclusions. The pipeline is efficient and sensitive enough to parse full alert streams. It is currently being applied to the live ZTF alert stream to identify and follow-up future candidates while active. This pipeline could be the foundation for glSNe Ia searches in future surveys, such as the Rubin Observatory Legacy Survey of Space and Time.

Tidal disruption events (TDEs) that are spatially offset from the nuclei of their host galaxies offer a new probe of massive black hole (MBH) wanderers, binaries, triples, and recoiling MBHs. Here we present AT2024tvd, the first off-nuclear TDE identified through optical sky surveys. High-resolution imaging with the Hubble Space Telescope shows that AT2024tvd is 0 . ″ 914 ± 0 . ″ 010 offset from the apparent center of its host galaxy, corresponding to a projected distance of 0.808 ± 0.009 kpc at z = 0.045. Chandra and Very Large Array observations support the same conclusion for the TDE’s X-ray and radio emission. AT2024tvd exhibits typical properties of nuclear TDEs, including a persistent hot UV/optical component that peaks at L bb ∼ 6 × 10 43 erg s −1 , broad hydrogen lines in its optical spectra, and delayed brightening of luminous ( L X,peak ∼ 3 × 10 43 erg s −1 ), highly variable soft X-ray emission. The MBH mass of AT2024tvd is 10 6±1 M ⊙ , at least 10 times lower than its host galaxy’s central black hole mass (≳10 8 M ⊙ ). The MBH in AT2024tvd has two possible origins: a wandering MBH from the lower-mass galaxy in a minor merger during the dynamical friction phase or a recoiling MBH ejected by triple interactions. Combining AT2024tvd with two previously known off-nuclear TDEs discovered in X-rays (3XMM J2150 and EP240222a), which likely involve intermediate-mass black holes in satellite galaxies, we find that the parent galaxies of all three events are very massive (∼10 10.9 M ⊙ ). This result aligns with expectations from cosmological simulations that the number of offset MBHs scales linearly with the host halo mass.

We present multiwavelength analysis of ZTF23abelseb (AT 2023sva), an optically discovered fast-fading (Delta m(r)=2.2 mag in Delta t=0.74 d), luminous (M-r similar to-30.0 mag), and red (g-r=0.50 mag) transient at z=2.28 with accompanying luminous radio emission. AT 2023sva does not possess a gamma-ray burst (GRB) counterpart to an isotropic equivalent energy limit of E-gamma,(iso)<1.6x10(52) erg, determined through searching gamma-ray satellite archives between the last non-detection and first detection, making it the sixth example of an optically discovered afterglow with a redshift measurement and no detected GRB counterpart. We analyse AT 2023sva's optical, radio, and X-ray observations to characterize the source. From radio analyses, we find the clear presence of strong interstellar scintillation (ISS) 72 d after the initial explosion, allowing us to place constraints on the source's angular size and bulk Lorentz factor. When comparing the source sizes derived from ISS of orphan events to those of the classical GRB population, we find orphan events have statistically smaller source sizes. We also utilize Bayesian techniques to model the multiwavelength afterglow. Within this framework, we find evidence that AT 2023sva possesses a shallow power-law structured jet viewed slightly off-axis (theta(v)=0.07 +/- 0.02) just outside of the jet's core opening angle (theta(c)=0.06 +/- 0.02). We determine this is likely the reason for the lack of a detected GRB counterpart, but also investigate other scenarios. AT 2023sva's evidence for possessing a structured jet stresses the importance of broadening orphan afterglow search strategies to a diverse range of GRB jet angular energy profiles, to maximize the return of future optical surveys.

Although all Type II supernovae (SNe) originate from massive stars possessing a hydrogen-rich envelope, their light-curve morphology is diverse, reflecting poorly characterized heterogeneity in the physical properties of their progenitor systems. Here, we present a detailed light-curve analysis of a magnitude-limited sample of 639 Type II SNe from the Zwicky Transient Facility Bright Transient Survey. Using Gaussian processes, we systematically measure empirical light-curve features (e.g. rise times, peak colours, and luminosities) in a robust sampling-independent manner. We focus on rise times as they are highly sensitive to pre-explosion progenitor properties, especially the presence of a dense circumstellar medium (CSM) shed by the progenitor in the years immediately pre-explosion. By correlating our feature measurements with physical parameters from an extensive grid of stella hydrodynamical models with varying progenitor properties (CSM structure, M-center dot, R-CSM, and M-ZAMS), we quantify the proportion of events with sufficient pre-explosion mass loss to significantly alter the initial light curve (roughly MCSM >= 10(-2.5)M(circle dot)) in a highly complete sample of 377 spectroscopically classified Type II SNe. We find that 67 +/- 6 per cent of observed SNe in our magnitude-limited sample show evidence for substantial CSM (M-CSM >= 10(-2.5)M(circle dot)) close to the progenitor (R-CSM < 10(15) cm) at the time of explosion. After applying a volumetric-correction, we find 36(-7)(+5) per cent of all Type II SN progenitors possess substantial CSM within 1015 cm at the time of explosion. This high fraction of progenitors with dense CSM, supported by photometric and spectroscopic evidence of previous SNe, reveals mass-loss rates significantly exceeding those measured in local group red supergiants or predicted by current theoretical models.

Type Ia supernova (SN Ia) cosmology relies on the estimation of lightcurve parameters to derive precision distances that leads to the estimation of cosmological parameters. The empirical SALT2 lightcurve modeling that relies on only two parameters, a stretch x1, and a color c, has been used by the community for almost two decades. In this paper we study the ability of the SALT2 model to fit the nearly 3000 cosmology-grade SN Ia lightcurves from the second release of the Zwicky Transient Facility (ZTF) cosmology science working group. While the ZTF data was not used to train SALT2, the algorithm is modeling the ZTF SN Ia optical lightcurves remarkably well, except for lightcurve points prior to -10 d from maximum, where the training critically lacks statistics. We find that the lightcurve fitting is robust against the considered choice of phase-range, but we show the [-10; +40] d range to be optimal in terms of statistics and accuracy. We do not detect any significant features in the lightcurve fit residuals that could be connected to the host environment. Potential systematic population differences related to the SN Ia host properties might thus not be accountable for by the addition of extra lightcurve parameters. However, a small but significant inconsistency between residuals of blue- and red-SN Ia strongly suggests the existence of a phase-dependent color term, with potential implications for the use of SNe Ia in precision cosmology. We thus encourage modellers to explore this avenue and we emphasize the importance that SN Ia cosmology must include a SALT2 retraining to accurately model the lightcurves and avoid biasing the derivation of cosmological parameters.

We present the discovery of PSR J1947–1120, a new huntsman millisecond pulsar with a red giant companion star in a 10.3 day orbit. This pulsar was found via optical, X-ray, and radio follow-up of the previously unassociated γ -ray source 4FGL J1947.6–1121. PSR J1947–1120 is the second confirmed pulsar in the huntsman class and establishes this as a bona fide subclass of millisecond pulsars. We use MESA models to show that huntsman pulsars can be naturally explained as neutron star binaries whose secondaries are currently in the “red bump” region of the red giant branch, temporarily underfilling their Roche lobes and hence halting mass transfer. Huntsman pulsars offer a new view of the formation of typical millisecond pulsars, allowing novel constraints on the efficiency of mass transfer and recycling at an intermediate stage in the process.

Main-sequence dipper stars, characterized by irregular and often aperiodic luminosity dimming events, offer a unique opportunity to explore the variability of circumstellar material and its potential links to planet formation, debris disks, and broadly star–planet interactions. The advent of all-sky time-domain surveys has enabled the rapid discovery of these unique systems. We present the results of a large systematic search for main-sequence dipper stars, conducted across a sample of 63 million FGK main-sequence stars using data from Gaia eDR3 and the Zwicky Transient Facility survey. Using a novel light-curve scoring algorithm and a scalable workflow tailored for analyzing millions of light curves, we have identified 81 new dipper star candidates. Our sample reveals a diverse phenomenology of light-curve dimming shapes, such as skewed and symmetric dimmings with timescales spanning days to years, some of which closely resemble exaggerated versions of KIC 8462852. Our sample reveals no clear periodicity patterns in many of these dippers and no IR excess or variability. Using archival data collated for this study, we thoroughly investigate several classification scenarios and hypothesize that the mechanisms of such dimming events are either driven by circumstellar clumps or occultations by stellar/substellar companions with disks. Our study marks a significant step forward in understanding main-sequence dipper stars.

We analyze imaging from Data Preview 1 of the Vera C. Rubin Observatory to explore the performance of early LSST pipelines in the 47 Tucanae field. The coadd- object catalog demonstrates the depth and precision possible with Rubin, recovering well-defined color–magnitude diagrams for 47 Tuc Small Magellanic Cloud. Unfortunately, the existing pipelines fail to recover sources within ∼28 pc of the cluster center, due to the extreme source density. Using Rubin’s forced photometry on stars identified via Difference Imaging, we can recover sources down to ∼14 pc from the cluster center, and find 14,744 potential cluster members with this extended data set. While this forced photometry has significant systematics, our analysis showcases the potential for detailed structural studies of crowded fields with the Rubin Observatory.

We present the optical discovery and multiwavelength follow-up observations of AT 2024kmq, a likely tidal disruption event (TDE) associated with a supermassive ( M _BH ∼ 10 ^8 M _⊙ ) black hole in a massive galaxy at z = 0.192. The optical light curve of AT 2024kmq exhibits two distinct peaks: an early fast (timescale 1 day) and luminous ( M ≈ −20 mag) red peak, then a slower (timescale 1 month) blue peak with a higher optical luminosity ( M ≈ −22 mag) and featureless optical spectra. The second component is similar to the spectroscopic class of “featureless TDEs” in the literature, and during this second component we detect highly variable, luminous ( L _X ≈ 10 ^44 erg s ^−1 ), and hard ( f _ν ∝ ν ^−1.5 ) X-ray emission. Luminous (10 ^29 erg s ^−1 Hz ^−1 at 10 GHz) but unchanging radio emission likely arises from an underlying active galactic nucleus. The luminosity, timescale, and color of the early red optical peak can be explained by synchrotron emission, or alternatively by thermal emission from material at a large radius ( R ≈ a few × 10 ^15 cm). Possible physical origins for this early red component include an off-axis relativistic jet, and shocks from self-intersecting debris leading to the formation of the accretion disk. Late-time radio observations will help distinguish between the two possibilities.

Astrophysical high-energy (TeV-PeV) neutrinos were first discovered in 2013, but their origin remains largely unknown. Here we present SN 2023uqf, a supernova found in coincidence with high-energy neutrino IC231004A, as part of a systematic optical follow-up program with the Zwicky Transient Facility. SN 2023uqf had a luminous and rapidly-evolving lightcurve, and spectroscopic observations indicated that the source was a Type Ibn supernova. Spectroscopic signatures confirm ongoing interaction between the supernova ejecta and a dense circumstellar medium, as expected for high-energy neutrino production in a core-collapse supernova. Given the rare nature of Type Ibn supernovae, SN 2023uqf is unlikely to have been discovered by chance over the course of our program (p=0.3