Search

It has been recently discovered that a few super-Eddington sources undergoing black hole super-Eddington accretion exhibit X-ray reflection signatures. In such new systems, one expects that the coronal X-ray emissions are mainly reflected by optically thick super-Eddington winds instead of thin disks. In this paper, we conduct a series of general relativistic ray-tracing and Monte Carlo radiative transfer simulations to model the X-ray reflection signatures, especially the characteristic Fe K$\alpha$ line, produced from super-Eddington accretion flows. In particular, we allow the photons emitted by a lamppost corona to be reflected multiple times in a cone-like funnel surrounded by fast winds. We find that the Fe K$\alpha$ line profile most sensitively depends on the wind kinematics, while its exact shape also depends on the funnel open angle and corona height. Furthermore, very interestingly, we find that the Fe K$\alpha$ line can have a prominent double-peak profile in certain parameter spaces even with a face-on orientation. Moreover, we compare the Fe K$\alpha$ line profiles produced from super-Eddington and thin disks and show that such lines can provide important insights into the understanding of black hole systems undergoing super-Eddington accretion., Comment: 23 pages, 21 figures, 2 tables. Comments are welcome

In the search for quantum advantage with near-term quantum devices, navigating the optimization landscape is significantly hampered by the barren plateaus phenomenon. This study presents a strategy to overcome this obstacle without changing the quantum circuit architecture. We propose incorporating auxiliary control qubits to shift the circuit from a unitary $2$-design to a unitary $1$-design, mitigating the prevalence of barren plateaus. We then remove these auxiliary qubits to return to the original circuit structure while preserving the unitary $1$-design properties. Our experiment suggests that the proposed structure effectively mitigates the barren plateaus phenomenon. A significant experimental finding is that the gradient of $\theta_{1,1}$, the first parameter in the quantum circuit, displays a broader distribution as the number of qubits and layers increases. This suggests a higher probability of obtaining effective gradients. This stability is critical for the efficient training of quantum circuits, especially for larger and more complex systems. The results of this study represent a significant advance in the optimization of quantum circuits and offer a promising avenue for the scalable and practical implementation of quantum computing technologies. This approach opens up new opportunities in quantum learning and other applications that require robust quantum computing power., Comment: 15 pages, 9 figures

We present the tidal disruption event (TDE) AT2022lri, hosted in a nearby ($\approx\!144$ Mpc) quiescent galaxy with a low-mass massive black hole ($10^4\,M_\odot < M_{\rm BH} < 10^6\,M_\odot$). AT2022lri belongs to the TDE-H+He subtype. More than 1 Ms of X-ray data were collected with NICER, Swift, and XMM-Newton from 187 d to 672 d after peak. The X-ray luminosity gradually declined from $1.5\times 10^{44}\,{\rm erg\,s^{-1}}$ to $1.5\times 10^{43}\,{\rm erg\,s^{-1}}$ and remains much above the UV and optical luminosity, consistent with a super-Eddington accretion flow viewed face-on. Sporadic strong X-ray dips atop a long-term decline are observed, with variability timescale of $\approx\!0.5$ hr--1 d and amplitude of $\approx\!2$--8. When fitted with simple continuum models, the X-ray spectrum is dominated by a thermal disk component with inner temperature going from $\sim\! 146$ eV to $\sim\! 86$ eV. However, there are residual features that peak around 1 keV, which, in some cases, cannot be reproduced by a single broad emission line. We analyzed a subset of time-resolved spectra with two physically motivated models describing either a scenario where ionized absorbers contribute extra absorption and emission lines or where disk reflection plays an important role. Both models provide good and statistically comparable fits, show that the X-ray dips are correlated with drops in the inner disk temperature, and require the existence of sub-relativistic (0.1--0.3$c$) ionized outflows. We propose that the disk temperature fluctuation stems from episodic drops of the mass accretion rate triggered by magnetic instabilities or/and wobbling of the inner accretion disk along the black hole's spin axis., Comment: 35 pages, 20 figures, submitted

Multi-peaked supernovae with precursors, dramatic light-curve rebrightenings, and spectral transformation are rare, but are being discovered in increasing numbers by modern night-sky transient surveys like the Zwicky Transient Facility (ZTF). Here, we present the observations and analysis of SN 2023aew, which showed a dramatic increase in brightness following an initial luminous (-17.4 mag) and long (~100 days) unusual first peak (possibly precursor). SN 2023aew was classified as a Type IIb supernova during the first peak but changed its type to resemble a stripped-envelope supernova (SESN) after the marked rebrightening. We present comparisons of SN 2023aew's spectral evolution with SESN subtypes and argue that it is similar to SNe Ibc during its main peak. P-Cygni Balmer lines are present during the first peak, but vanish during the second peak's photospheric phase, before H$\alpha$ resurfaces again during the nebular phase. The nebular lines ([O I], [Ca II], Mg I], H$\alpha$) exhibit a double-peaked structure which hints towards a clumpy or non-spherical ejecta. We analyze the second peak in the light curve of SN 2023aew and find it to be broader than normal SESNe as well as requiring a very high $^{56}$Ni mass to power the peak luminosity. We discuss the possible origins of SN 2023aew including an eruption scenario where a part of the envelope is ejected during the first peak which also powers the second peak of the light curve through SN-CSM interaction., Comment: 22 pages, 11 figures, 5 tables

Classical gamma-ray bursts (GRBs) have two distinct emission episodes: prompt emission from ultra-relativistic ejecta and afterglow from shocked circumstellar material. While both components are extremely luminous in known GRBs, a variety of scenarios predict the existence of luminous afterglow emission with little or no associated high-energy prompt emission. We present AT 2019pim, the first secure example of this phenomenon to be identified. Serendipitously discovered during follow-up observations of a gravitational-wave trigger and located in a contemporaneous TESS sector, it is hallmarked by a fast-rising (t ~ 2 hr), luminous (M_UV,peak ~ -24.4 mag) optical transient with accompanying luminous X-ray and radio emission. No gamma-ray emission consistent with the time and location of the transient was detected by Fermi-GBM or by Konus, placing strong limits on an accompanying GRB. We investigate several independent observational aspects of the afterglow in the context of constraints on relativistic motion and find all of them are consistent with an initial Lorentz factor of Gamma_0 ~ 30-50, significantly lower than in any well-observed GRB and consistent with the theoretically-predicted "dirty fireball" scenario in which the high-energy prompt emission is stifled by pair production. However, we cannot rule out a structured jet model in which only the line-of-sight material was ejected at low-Gamma, off-axis from a classical high-Gamma jet core. This event represents a milestone in orphan afterglow searches, demonstrating that luminous afterglows with weak or no detectable gamma-ray radiation exist in nature and can be discovered by high-cadence optical surveys., Comment: Submitted to MNRAS

When a star passes through the tidal disruption radius of a massive black hole (BH), it can be torn apart by the tidal force of the BH, known as the Tidal Disruption Event (TDE). Since the UV/optical emitting region inferred from the blackbody radius is significantly larger than the circularization radius predicted by the classical TDE theory, two competing models, stream collision and envelope reprocessing, were proposed to explain the unexpectedly large UV/optical emitting size. Here, we investigate the variability behaviors (cross-correlation and time delay) of three representative TDEs with continuum reverberation mapping. Our results demonstrate that TDE behavior is clearly inconsistent with the envelope reprocessing scenario. In contrast, the picture of the stream collision, together with the late-time formed accretion disk, can explain heterogeneous observations. This provides compelling evidence that the UV/optical emission originates from stream collisions during the early-stage of TDE evolution and gradually transitions to being dominated by accretion disk with detectable X-ray emission in a late stage. After fading back to a quiescent state, recurrent flares may be observed in some occasions, such as partial TDEs., Comment: 56 pages, 21 figures, 5 tables. Submitted, comments welcome!

Optical surveys have become increasingly adept at identifying candidate Tidal Disruption Events (TDEs) in large numbers, but classifying these generally requires extensive spectroscopic resources. Here we present $\texttt{tdescore}$, a simple binary photometric classifier that is trained using a systematic census of $\sim$3000 nuclear transients from the Zwicky Transient Facility (ZTF). The sample is highly imbalanced, with TDEs representing $\sim$2% of the total. $\texttt{tdescore}$ is nonetheless able to reject non-TDEs with 99.6% accuracy, yielding a sample of probable TDEs with recall of 77.5% for a precision of 80.2%. $\texttt{tdescore}$ is thus substantially better than any available TDE photometric classifier scheme in the literature, with performance not far from spectroscopy as a method for classifying ZTF nuclear transients, despite relying solely on ZTF data and multi-wavelength catalogue cross-matching. In a novel extension, we use `SHapley Additive exPlanations' ($\texttt{SHAP}$) to provide a human-readable justification for each individual $\texttt{tdescore}$ classification, enabling users to understand and form opinions about the underlying classifier reasoning. $\texttt{tdescore}$ can serve as a model for photometric identification of TDEs with time-domain surveys, such as the upcoming Rubin observatory., Comment: Accepted for publication, 16 pages, 7 figures, 2 tables

In recent years, certain luminous extragalactic optical transients have been observed to last only a few days. Their short observed duration implies a different powering mechanism from the most common luminous extragalactic transients (supernovae) whose timescale is weeks. Some short-duration transients, most notably AT2018cow, display blue optical colours and bright radio and X-ray emission. Several AT2018cow-like transients have shown hints of a long-lived embedded energy source, such as X-ray variability, prolonged ultraviolet emission, a tentative X-ray quasiperiodic oscillation, and large energies coupled to fast (but subrelativistic) radio-emitting ejecta. Here we report observations of minutes-duration optical flares in the aftermath of an AT2018cow-like transient, AT2022tsd (the "Tasmanian Devil"). The flares occur over a period of months, are highly energetic, and are likely nonthermal, implying that they arise from a near-relativistic outflow or jet. Our observations confirm that in some AT2018cow-like transients the embedded energy source is a compact object, either a magnetar or an accreting black hole., Comment: 79 pages, 3 figures (main text) + 7 figures (extended data) + 2 figures (supplementary information). Published online in Nature on 15 November 2023

The Advanced X-ray Imaging Satellite (AXIS) promises revolutionary science in the X-ray and multi-messenger time domain. AXIS will leverage excellent spatial resolution (<1.5 arcsec), sensitivity (80x that of Swift), and a large collecting area (5-10x that of Chandra) across a 24-arcmin diameter field of view to discover and characterize a wide range of X-ray transients from supernova-shock breakouts to tidal disruption events to highly variable supermassive black holes. The observatory's ability to localize and monitor faint X-ray sources opens up new opportunities to hunt for counterparts to distant binary neutron star mergers, fast radio bursts, and exotic phenomena like fast X-ray transients. AXIS will offer a response time of <2 hours to community alerts, enabling studies of gravitational wave sources, high-energy neutrino emitters, X-ray binaries, magnetars, and other targets of opportunity. This white paper highlights some of the discovery science that will be driven by AXIS in this burgeoning field of time domain and multi-messenger astrophysics., Comment: This White Paper is part of a series commissioned for the AXIS Probe Concept Mission; additional AXIS White Papers can be found at http://axis.astro.umd.edu

Neutron stars and stellar-mass black holes are the remnants of massive star explosions. Most massive stars reside in close binary systems, and the interplay between the companion star and the newly formed compact object has been theoretically explored, but signatures for binarity or evidence for the formation of a compact object during a supernova explosion are still lacking. Here we report a stripped-envelope supernova, SN 2022jli, which shows 12.4-day periodic undulations during the declining light curve. Narrow H$\alpha$ emission is detected in late-time spectra with concordant periodic velocity shifts, likely arising from hydrogen gas stripped from a companion and accreted onto the compact remnant. A new Fermi/LAT $\gamma$-ray source is temporally and positionally consistent with SN 2022jli. The observed properties of SN 2022jli, including periodic undulations in the optical light curve, coherent H$\alpha$ emission shifting, and evidence for association with a $\gamma$-ray source, point to the explosion of a massive star in a binary system leaving behind a bound compact remnant. Mass accretion from the companion star onto the compact object powers the light curve of the supernova and generates the $\gamma$-ray emission., Comment: Accepted for publication in Nature

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., Comment: 24 pages, 13 figures, submitted to ApJ

In this work, we use the Jansky VLA Sky Survey (VLASS) to compile the first sample of six radio-selected tidal disruption events (TDEs) with transient optical counterparts. While we still lack the statistics to do detailed population studies of radio-selected TDEs, we use these events to suggest trends in host galaxy and optical light curve properties that may correlate with the presence of radio emission, and hence can inform optically-selected TDE radio follow-up campaigns. We find that radio-selected TDEs tend to have faint and cool optical flares, as well as host galaxies with low SMBH masses. Our radio-selected TDEs also tend to have more energetic, larger radio emitting regions than radio-detected, optically-selected TDEs. We consider possible explanations for these trends, including by invoking super-Eddington accretion and enhanced circumnuclear media. Finally, we constrain the radio-emitting TDE rate to be $\gtrsim 10$ Gpc$^{-3}$ yr$^{-1}$., Comment: 26 pages, 5 tables, 11 figures, submitted to ApJ

Galactic nuclei showing recurrent phases of activity and quiescence have recently been discovered, with recurrence times as short as a few hours to a day -- known as quasi-periodic X-ray eruption (QPE) sources -- to as long as hundreds to a thousand days for repeating nuclear transients (RNTs). Here we present a multi-wavelength overview of Swift J023017.0+283603 (hereafter Swift J0230+28), a source that exhibits repeating and quasi-periodic X-ray flares from the nucleus of a previously unremarkable galaxy at $\sim$ 165 Mpc, with a recurrence time of approximately 22 days, an intermediary timescale between known RNTs and QPE sources. The source also shows transient radio emission, likely associated with the X-ray emission. Such recurrent soft X-ray eruptions, with no accompanying UV/optical emission, are strikingly similar to QPE sources. However, in addition to having a recurrence time that is $\sim 25$ times longer than the longest-known QPE source, Swift J0230+28's eruptions exhibit somewhat distinct shapes and temperature evolution than the known QPE sources. Scenarios involving extreme mass ratio inspirals are favored over disk instability models. The source reveals an unexplored timescale for repeating extragalactic transients and highlights the need for a wide-field, time-domain X-ray mission to explore the parameter space of recurring X-ray transients., Comment: Final version, appeared on Nature Astronomy on 12 January 2024

We present a systematic analysis of the X-ray emission of a sample of 17 optically selected, X-ray-detected tidal disruption events (TDEs) discovered between 2014 and 2021. The X-ray light curves show a diverse range of temporal behaviors, with most sources not following the expected power-law decline. The X-ray spectra are mostly extremely soft and consistent with thermal emission from the innermost region of an accretion disk, which cools as the accretion rate decreases. Three sources show formation of a hard X-ray corona, at late-times. The spectral energy distribution shape, probed by the ratio ($L_{\rm\,BB}/L_{\rm\,X}$) between the UV/optical and X-ray, shows a wide range $L_{\rm BB}/L_{\rm X}\,\in\,(0.5,\,3000)$ at early-times, and converges to disk-like values $L_{\rm\,BB}/L_{\rm\,X}\,\in\,(0.5,\,10)$ at late-times. We estimate the fraction of optically discovered TDEs with $L_{\rm\,X}\,\geq 10^{42}~\rm{erg}~\rm{s}^{-1} $ to be at least $40\%$, and show that X-ray loudness is independent of black hole mass. We argue that distinct disk formation time scales are unlikely to fully explain the diverse range of X-ray evolutions. We combine our sample with X-ray discovered ones to construct an X-ray luminosity function, best fitted by a broken power-law, with a break at $L_{\rm\,X} \approx 10^{44}~\rm{erg}~\rm{s}^{-1} $. We show that there is no dichotomy between optically and X-ray selected TDEs, instead there is a continuum of early time $L_{\rm\,BB}/L_{\rm\,X}$, at least as wide as $L_{\rm\,BB}/L_{\rm\,X}\,\in\,(0.1,\,3000)$, with optical/X-ray surveys selecting preferentially, but not exclusively, from the higher/lower end of the distribution. Our findings are consistent with unification models for the overall TDE population., Comment: Accepted (ApJ), moderate modifications after peer review

AT2022cmc was recently reported as the first on-axis jetted tidal disruption event (TDE) discovered in the last decade, and the fourth on-axis jetted TDE candidate known so far. In this work, we present NuSTAR hard X-ray (3--30 keV) observations of AT2022cmc, as well as soft X-ray (0.3--6 keV) observations obtained by NICER, Swift, and XMM-Newton. Our analysis reveals that the broadband X-ray spectra can be well described by a broken power-law with $f_\nu \propto \nu^{-0.5}$ ($f_\nu \propto \nu^{-1}$) below (above) the rest-frame break energy of $E_{\rm bk}\sim 10$ keV at observer-frame $t_{\rm obs}=7.8$ and 17.6 days since discovery. At $t_{\rm obs} = 36.2$ days, the X-ray spectrum is consistent with either a single power-law or a broken power-law. By modeling the spectral energy distribution evolution from radio to hard X-ray across the three NuSTAR observing epochs, we find that the sub-millimeter/radio emission originates from external shocks at large distances $\gtrsim\! 10^{17}$ cm from the black hole, the UV/optical light comes from a thermal envelope with radius $\sim\!10^{15}$ cm, and the X-ray emission is consistent with synchrotron radiation powered by energy dissipation at intermediate radii within the (likely magnetically dominated) jet. We constrain the bulk Lorentz factor of the jet to be of the order 10--100. Our interpretation differs from the model proposed by Pasham et al. (2023) where both the radio and X-rays come from the same emitting zone in a matter-dominated jet. Our model for the jet X-ray emission has broad implications on the nature of relativistic jets in other sources such as gamma-ray bursts., Comment: 17 pages, 11 figures, ApJ accepted

The host galaxies of tidal disruption events (TDEs) have been shown to possess peculiar properties, including high central light concentrations, unusual star-formation histories, and ``green'' colors. The ubiquity of these large-scale galaxy characteristics among TDE host populations suggests they may serve to boost the TDE rate in such galaxies by influencing the nuclear stellar dynamics. We present the first population study of integral field spectroscopy for thirteen TDE host galaxies across all spectral classes and X-ray brightnesses with the purpose of investigating their large-scale properties. We derive the black hole masses via stellar kinematics (i.e., the $M-\sigma$ relation) and find masses in the range $5.0 \lesssim \log(M_{\rm BH}/M_\odot) \lesssim 8.0$, with a distribution dominated by black holes with $M_{\rm BH} \sim 10^6 M_\odot$. We find one object with $M_{\rm BH} \gtrsim 10^8 M_\odot$, above the ``Hills mass'', which if the disrupted star was of solar type, allows a lower limit of $a \gtrsim 0.16$ to be placed on its spin, lending further support to the proposed connection between featureless TDEs and jetted TDEs. We also explore the level of rotational support in the TDE hosts, quantified by $(V/\sigma)_e$, a parameter which has been shown to correlate with stellar age and may explain the peculiar host galaxy preferences of TDEs. We find that the TDE hosts exhibit a broad range in $(V/\sigma)_e$ following a similar distribution as E+A galaxies, which have been shown to be overrepresented among TDE host populations., Comment: 20 pages, 13 figures, 4 tables; accepted to ApJ

Eruptive mass loss of massive stars prior to supernova (SN) explosion is key to understanding their evolution and end fate. An observational signature of pre-SN mass loss is the detection of an early, short-lived peak prior to the radioactive-powered peak in the lightcurve of the SN. This is usually attributed to the SN shock passing through an extended envelope or circumstellar medium (CSM). Such an early peak is common for double-peaked Type IIb SNe with an extended Hydrogen envelope but is uncommon for normal Type Ibc SNe with very compact progenitors. In this paper, we systematically study a sample of 14 double-peaked Type Ibc SNe out of 475 Type Ibc SNe detected by the Zwicky Transient Facility. The rate of these events is ~ 3-9 % of Type Ibc SNe. A strong correlation is seen between the peak brightness of the first and the second peak. We perform a holistic analysis of this sample's photometric and spectroscopic properties. We find that six SNe have ejecta mass less than 1.5 Msun. Based on the nebular spectra and lightcurve properties, we estimate that the progenitor masses for these are less than ~ 12 Msun. The rest have an ejecta mass > 2.4 Msun and a higher progenitor mass. This sample suggests that the SNe with low progenitor masses undergo late-time binary mass transfer. Meanwhile, the SNe with higher progenitor masses are consistent with wave-driven mass loss or pulsation-pair instability-driven mass loss simulations., Comment: Accepted for publication in ApJ

SN 1987A was an unusual hydrogen-rich core-collapse supernova originating from a blue supergiant star. Similar blue supergiant explosions remain a small family of events, and are broadly characterized by their long rises to peak. The Zwicky Transient Facility (ZTF) Census of the Local Universe (CLU) experiment aims to construct a spectroscopically complete sample of transients occurring in galaxies from the CLU galaxy catalog. We identify 13 long-rising (>40 days) Type II supernovae from the volume-limited CLU experiment during a 3.5 year period from June 2018 to December 2021, approximately doubling the previously known number of these events. We present photometric and spectroscopic data of these 13 events, finding peak r-band absolute magnitudes ranging from -15.6 to -17.5 mag and the tentative detection of Ba II lines in 9 events. Using our CLU sample of events, we derive a long-rising Type II supernova rate of $1.37^{+0.26}_{-0.30}\times10^{-6}$ Mpc$^{-3}$ yr$^{-1}$, $\approx$1.4% of the total core-collapse supernova rate. This is the first volumetric rate of these events estimated from a large, systematic, volume-limited experiment., Comment: 33 pages, 18 figures, 5 tables. Updated to ApJ accepted version

Abridged - Stars with ZAMS masses between 140 and $260 M_\odot$ are thought to explode as pair-instability supernovae (PISNe). During their thermonuclear runaway, PISNe can produce up to several tens of solar masses of radioactive nickel, resulting in luminous transients similar to some superluminous supernovae (SLSNe). Yet, no unambiguous PISN has been discovered so far. SN2018ibb is a H-poor SLSN at $z=0.166$ that evolves extremely slowly compared to the hundreds of known SLSNe. Between mid 2018 and early 2022, we monitored its photometric and spectroscopic evolution from the UV to the NIR with 2-10m class telescopes. SN2018ibb radiated $>3\times10^{51} \rm erg$ during its evolution, and its bolometric light curve reached $>2\times10^{44} \rm erg\,s^{-1}$ at peak. The long-lasting rise of $>93$ rest-frame days implies a long diffusion time, which requires a very high total ejected mass. The PISN mechanism naturally provides both the energy source ($^{56}$Ni) and the long diffusion time. Theoretical models of PISNe make clear predictions for their photometric and spectroscopic properties. SN2018ibb complies with most tests on the light curves, nebular spectra and host galaxy, potentially all tests with the interpretation we propose. Both the light curve and the spectra require 25-44 $M_\odot$ of freshly nucleosynthesised $^{56}$Ni, pointing to the explosion of a metal-poor star with a He-core mass of 120-130 $M_\odot$ at the time of death. This interpretation is also supported by the tentative detection of [Co II]$\lambda$1.025$\mu$m, which has never been observed in any other PISN candidate or SLSN before. Powering by a central engine, such as a magnetar or a black hole, can be excluded with high confidence. This makes SN2018ibb by far the best candidate for being a PISN, to date., Comment: Accepted in A&A, the revised version includes a PISN rate estimate and an additional test with PISN models. 47 pages, main text 41 pages, 38 figures, 16 Tables

We conduct a systematic tidal disruption event (TDE) demographics analysis using the largest sample of optically selected TDEs. A flux-limited, spectroscopically complete sample of 33 TDEs is constructed using the Zwicky Transient Facility over three years (from October 2018 to September 2021). We infer the black hole (BH) mass ($M_{\rm BH}$) with host galaxy scaling relations, showing that the sample $M_{\rm BH}$ ranges from $10^{5.1}\,M_\odot$ to $10^{8.2}\,M_\odot$. We developed a survey efficiency corrected maximum volume method to infer the rates. The rest-frame $g$-band luminosity function (LF) can be well described by a broken power-law of $\phi (L_g)\propto [(L_g / L_{\rm bk})^{0.3} + (L_g / L_{\rm bk})^{2.6}]^{-1}$, with $L_{\rm bk}=10^{43.1}\,{\rm erg\,s^{-1}}$. In the BH mass regime of $10^{5.3}\lesssim (M_{\rm BH}/M_\odot) \lesssim 10^{7.3}$, the TDE mass function follows $\phi(M_{\rm BH})\propto M_{\rm BH}^{-0.25}$, which favors a flat local BH mass function ($dn_{\rm BH}/d{\rm log}M_{\rm BH}\approx{\rm constant}$). We confirm the significant rate suppression at the high-mass end ($M_{\rm BH}\gtrsim 10^{7.5}\,M_\odot$), which is consistent with theoretical predictions considering direct capture of hydrogen-burning stars by the event horizon. At a host galaxy mass of $M_{\rm gal}\sim 10^{10}\,M_\odot$, the average optical TDE rate is $\approx 3.2\times 10^{-5}\,{\rm galaxy^{-1}\,yr^{-1}}$. We constrain the optical TDE rate to be [3.7, 7.4, and 1.6$]\times 10^{-5}\,{\rm galaxy^{-1}\,yr^{-1}}$ in galaxies with red, green, and blue colors., Comment: Replaced following peer-review process. 38 pages, 23 figures. Accepted for publication in ApJL

One of the open questions following the discovery of GW170817 is whether neutron star mergers are the only astrophysical sites capable of producing $r$-process elements. Simulations have shown that 0.01-0.1M$_\odot$ of $r$-process material could be generated in the outflows originating from the accretion disk surrounding the rapidly rotating black hole that forms as a remnant to both neutron star mergers and collapsing massive stars associated with long-duration gamma-ray bursts (collapsars). The hallmark signature of $r$-process nucleosynthesis in the binary neutron star merger GW170817 was its long-lasting near-infrared emission, thus motivating a systematic photometric study of the light curves of broadlined stripped-envelope (Ic-BL) supernovae (SNe) associated with collapsars. We present the first systematic study of 25 SNe Ic-BL -- including 18 observed with the Zwicky Transient Facility and 7 from the literature -- in the optical/near-infrared bands to determine what quantity of $r$-process material, if any, is synthesized in these explosions. Using semi-analytic models designed to account for $r$-process production in SNe Ic-BL, we perform light curve fitting to derive constraints on the $r$-process mass for these SNe. We also perform independent light curve fits to models without $r$-process. We find that the $r$-process-free models are a better fit to the light curves of the objects in our sample. Thus we find no compelling evidence of $r$-process enrichment in any of our objects. Further high-cadence infrared photometric studies and nebular spectroscopic analysis would be sensitive to smaller quantities of $r$-process ejecta mass or indicate whether all collapsars are completely devoid of $r$-process nucleosynthesis., Comment: 37 pages, 8 figures, 4 tables

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

Tidal disruption events (TDEs) are bursts of electromagnetic energy released when supermassive black holes (SMBHs) at the centers of galaxies violently disrupt a star that passes too close. TDEs provide a new window to study accretion onto SMBHs; in some rare cases, this accretion leads to launching of a relativistic jet, but the necessary conditions are not fully understood. The best studied jetted TDE to date is Swift J1644+57, which was discovered in gamma-rays, but was too obscured by dust to be seen at optical wavelengths. Here we report the optical discovery of AT2022cmc, a rapidly fading source at cosmological distance (redshift z=1.19325) whose unique lightcurve transitioned into a luminous plateau within days. Observations of a bright counterpart at other wavelengths, including X-rays, sub-millimeter, and radio, supports the interpretation of AT2022cmc as a jetted TDE containing a synchrotron "afterglow", likely launched by a SMBH with spin $a \gtrsim 0.3$. Using 4 years of Zwicky Transient Facility (ZTF) survey data, we calculate a rate of $0.02 ^{+ 0.04 }_{- 0.01 }$ Gpc$^{-3}$ yr$^{-1}$ for on-axis jetted TDEs based on the luminous, fast-fading red component, thus providing a measurement complementary to the rates derived from X-ray and radio observations. Correcting for the beaming angle effects, this rate confirms that about 1% of TDEs have relativistic jets. Optical surveys can use AT2022cmc as a prototype to unveil a population of jetted TDEs., Comment: Accepted for publication in Nature

Luminous red novae (LRNe) are transients characterized by low luminosities and expansion velocities, and are associated with mergers or common envelope ejections in stellar binaries. Intermediate-luminosity red transients (ILRTs) are an observationally similar class with unknown origins, but generally believed to either be electron capture supernovae (ECSN) in super-AGB stars, or outbursts in dusty luminous blue variables (LBVs). In this paper, we present a systematic sample of 8 LRNe and 8 ILRTs detected as part of the Census of the Local Universe (CLU) experiment on the Zwicky Transient Facility (ZTF). The CLU experiment spectroscopically classifies ZTF transients associated with nearby ($<150$ Mpc) galaxies, achieving 80% completeness for m$_{r}<20$\,mag. Using the ZTF-CLU sample, we derive the first systematic LRNe volumetric-rate of 7.8$^{+6.5}_{-3.7}\times10^{-5}$ Mpc$^{-3}$ yr$^{-1}$ in the luminosity range $-16\leq$M$_{\rm{r}}$$\leq -11$ mag. We find that in this luminosity range, the LRN rate scales as dN/dL $\propto L^{-2.5\pm0.3}$ - significantly steeper than the previously derived scaling of $L^{-1.4\pm0.3}$ for lower luminosity LRNe (M$_{V}\geq-10$). The steeper power law for LRNe at high luminosities is consistent with the massive merger rates predicted by binary population synthesis models. We find that the rates of the brightest LRNe (M$_{r}\leq-13$ mag) are consistent with a significant fraction of them being progenitors of double compact objects (DCOs) that merge within a Hubble time. For ILRTs, we derive a volumetric rate of $2.6^{+1.8}_{-1.4}\times10^{-6}$ Mpc$^{-3}$yr$^{-1}$ for M$_{\rm{r}}\leq-13.5$, that scales as dN/dL $\propto L^{-2.5\pm0.5}$. This rate is $\approx1-5\%$ of the local core-collapse supernova rate, and is consistent with theoretical ECSN rate estimates., Comment: 32 pages, 16 figures. Submitted to ApJ, comments welcome

A rare class of supernovae (SNe) is characterized by strong interaction between the ejecta and several solar masses of circumstellar matter (CSM) as evidenced by strong Balmer-line emission. Within the first few weeks after the explosion, they may display spectral features similar to overluminous Type Ia SNe, while at later phase their observation properties exhibit remarkable similarities with some extreme case of Type IIn SNe that show strong Balmer lines years after the explosion. We present polarimetric observations of SN2018evt obtained by the ESO Very Large Telescope from 172 to 219 days after the estimated time of peak luminosity to study the geometry of the CSM. The nonzero continuum polarization decreases over time, suggesting that the mass loss of the progenitor star is aspherical. The prominent H$\alpha$ emission can be decomposed into a broad, time-evolving component and an intermediate-width, static component. The former shows polarized signals, and it is likely to arise from a cold dense shell (CDS) within the region between the forward and reverse shocks. The latter is significantly unpolarized, and it is likely to arise from shocked, fragmented gas clouds in the H-rich CSM. We infer that SN2018evt exploded inside a massive and aspherical circumstellar cloud. The symmetry axes of the CSM and the SN appear to be similar. SN\,2018evt shows observational properties common to events that display strong interaction between the ejecta and CSM, implying that they share similar circumstellar configurations. Our preliminary estimate also suggests that the circumstellar environment of SN2018evt has been significantly enriched at a rate of $\sim0.1$ M$_\odot$ yr$^{-1}$ over a period of $>100$ yr., Comment: 29 pages, 20 figures, 10 tables, accepted for publication in MNRAS

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

The fate of stars in the zero-age main-sequence (ZAMS) range $\approx 8-12$ Msun is unclear. They could evolve to form white dwarfs or explode as electron-capture supernovae (SNe) or iron core-collapse SNe (CCSNe). Even though the initial mass function indicates that this mass range should account for over 40% of all CCSNe progenitors, few have been observationally confirmed, likely owing to the faintness and rapid evolution of these transients. In this paper, we present a sample of nine Ca-rich/O-poor Type IIb SNe detected by the Zwicky Transient Facility with progenitors likely in this mass range. We perform a holistic analysis of the spectroscopic and photometric properties of the sample. These sources have a flux ratio of [Ca II] $\lambda \lambda$7291, 7324 to [O I] $\lambda \lambda$6300, 6364 of $\gtrsim$ 2 in their nebular spectra. Comparing the measured [O I] luminosity ($\lesssim 10^{39} \mathrm{erg\ s^{-1}}$) and derived oxygen mass ($\lesssim 0.1$ Msun) with theoretical models, we infer that the progenitor ZAMS mass for these explosions is less than 12 Msun. These correspond to He-stars with core masses less than around 3 Msun. We find that the ejecta properties (Mej $\lesssim 1$ Msun) are also consistent with those expected for such low mass He-stars. The low ejecta mass of these sources indicates a class of strongly-stripped SNe that is a transition between the regular stripped-envelope SNe and ultra-stripped SNe. The progenitor could be stripped by a main sequence companion and result in the formation of a neutron star $-$ main sequence binary. Such binaries have been suggested to be progenitors of neutron star $-$ white dwarf systems that could merge within a Hubble time, and be detectable with LISA., Comment: Accepted for publication in ApJ

We present a sample of 14 hydrogen-rich superluminous supernovae (SLSNe II) from the Zwicky Transient Facility (ZTF) between 2018 and 2020. We include all classified SLSNe with peaks $M_{g}<-20$ mag and with observed \emph{broad} but not narrow Balmer emission, corresponding to roughly 20 per cent of all hydrogen-rich SLSNe in ZTF phase I. We examine the light curves and spectra of SLSNe II and attempt to constrain their power source using light-curve models. The brightest events are photometrically and spectroscopically similar to the prototypical SN 2008es, while others are found spectroscopically more reminiscent of non-superluminous SNe II, especially SNe II-L. $^{56}$Ni decay as the primary power source is ruled out. Light-curve models generally cannot distinguish between circumstellar interaction (CSI) and a magnetar central engine, but an excess of ultraviolet (UV) emission signifying CSI is seen in most of the SNe with UV data, at a wide range of photometric properties. Simultaneously, the broad H$\alpha$ profiles of the brightest SLSNe II can be explained through electron scattering in a symmetric circumstellar medium (CSM). In other SLSNe II without narrow lines, the CSM may be confined and wholly overrun by the ejecta. CSI, possibly involving mass lost in recent eruptions, is implied to be the dominant power source in most SLSNe II, and the diversity in properties is likely the result of different mass loss histories. Based on their radiated energy, an additional power source may be required for the brightest SLSNe II, however -- possibly a central engine combined with CSI., Comment: 27 pages including appendices; 14 figures. Revised version, resubmitted to MNRAS after the referee's comments

We present X-ray, UV, optical, and radio observations of the nearby ($\approx78$ Mpc) tidal disruption event (TDE) AT2021ehb/ZTF21aanxhjv during its first 430 days of evolution. AT2021ehb occurs in the nucleus of a galaxy hosting a $\approx 10^{7}\,M_\odot$ black hole ($M_{\rm BH}$ inferred from host galaxy scaling relations). High-cadence Swift and NICER monitoring reveals a delayed X-ray brightening. The spectrum first undergoes a gradual ${\rm soft }\rightarrow{\rm hard}$ transition and then suddenly turns soft again within 3 days at $\delta t\approx 272$ days during which the X-ray flux drops by a factor of ten. In the joint NICER+NuSTAR observation ($\delta t =264$ days, harder state), we observe a prominent non-thermal component up to 30 keV and an extremely broad emission line in the iron K band. The bolometric luminosity of AT2021ehb reaches a maximum of $6.0^{+10.4}_{-3.8}\% L_{\rm Edd}$ when the X-ray spectrum is the hardest. During the dramatic X-ray evolution, no radio emission is detected, the UV/optical luminosity stays relatively constant, and the optical spectra are featureless. We propose the following interpretations: (i) the ${\rm soft }\rightarrow{\rm hard}$ transition may be caused by the gradual formation of a magnetically dominated corona; (ii) hard X-ray photons escape from the system along solid angles with low scattering optical depth ($\sim\,$a few) whereas the UV/optical emission is likely generated by reprocessing materials with much larger column density -- the system is highly aspherical; (iii) the abrupt X-ray flux drop may be triggered by the thermal-viscous instability in the inner accretion flow leading to a much thinner disk., Comment: 35 pages, 21 figures, accepted by ApJ

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

We present results from extensive broadband follow-up of GRB 210204A over the period of thirty days. We detect optical flares in the afterglow at 7.6 x 10^5 s and 1.1 x 10^6 s after the burst: the most delayed flaring ever detected in a GRB afterglow. At the source redshift of 0.876, the rest-frame delay is 5.8 x 10^5 s (6.71 d). We investigate possible causes for this flaring and conclude that the most likely cause is a refreshed shock in the jet. The prompt emission of the GRB is within the range of typical long bursts: it shows three disjoint emission episodes, which all follow the typical GRB correlations. This suggests that GRB 210204A might not have any special properties that caused late-time flaring, and the lack of such detections for other afterglows might be resulting from the paucity of late-time observations. Systematic late-time follow-up of a larger sample of GRBs can shed more light on such afterglow behaviour. Further analysis of the GRB 210204A shows that the late time bump in the light curve is highly unlikely due to underlying SNe at redshift (z) = 0.876 and is more likely due to the late time flaring activity. The cause of this variability is not clearly quantifiable due to the lack of multi-band data at late time constraints by the bad weather conditions. The flare of GRB 210204A is the latest flare detected to date., Comment: 18 pages, 12 figures

The Fermi Gamma-ray Burst Monitor (GBM) triggers on-board in response to $\sim$ 40 short gamma-ray bursts (SGRBs) per year; however, their large localization regions have made the search for optical counterparts a challenging endeavour. We have developed and executed an extensive program with the wide field of view of the Zwicky Transient Facility (ZTF) camera, mounted on the Palomar 48 inch Oschin telescope (P48), to perform target-of-opportunity (ToO) observations on 10 Fermi-GBM SGRBs during 2018 and 2020-2021. Bridging the large sky areas with small field of view optical telescopes in order to track the evolution of potential candidates, we look for the elusive SGRB afterglows and kilonovae (KNe) associated with these high-energy events. No counterpart has yet been found, even though more than 10 ground based telescopes, part of the Global Relay of Observatories Watching Transients Happen (GROWTH) network, have taken part in these efforts. The candidate selection procedure and the follow-up strategy have shown that ZTF is an efficient instrument for searching for poorly localized SGRBs, retrieving a reasonable number of candidates to follow-up and showing promising capabilities as the community approaches the multi-messenger era. Based on the median limiting magnitude of ZTF, our searches would have been able to retrieve a GW170817-like event up to $\sim$ 200 Mpc and SGRB afterglows to z = 0.16 or 0.4, depending on the assumed underlying energy model. Future ToOs will expand the horizon to z = 0.2 and 0.7 respectively.

Tumer et al. 2022, ATel #15171, have recently reported the discovery of an X-ray source, NuSTAR J053449+2126.0, during a calibration observation which took place on 25 April 2020. We scan the Zwicky Transient Facility (ZTF) alerts and archival photometry to determine the nature of the source. Palomar Gattini-IR is searched as well. We identify no obvious counterpart candidate. Follow-up X-ray and optical studies are needed to determine the true counterpart., Comment: Accepted and Published in RNAAS

Tidal disruption events (TDEs) offer a unique way to study dormant black holes. While the number of observed TDEs has grown thanks to the emergence of wide-field surveys in the past few decades, questions regarding the nature of the observed optical, UV, and X-ray emission remain. We present a uniformly selected sample of 30 spectroscopically classified TDEs from the Zwicky Transient Facility Phase I survey operations with follow-up \textit{Swift} UV and X-ray observations. Through our investigation into correlations between light curve properties, we recover a shallow positive correlation between the peak bolometric luminosity and decay timescales. We introduce a new spectroscopic class of TDE, TDE-featureless, which are characterized by featureless optical spectra. The new TDE-featureless class shows larger peak bolometric luminosities, peak blackbody temperatures, and peak blackbody radii. We examine the differences between the X-ray bright and X-ray faint populations of TDEs in this sample, finding that X-ray bright TDEs show higher peak blackbody luminosities than the X-ray faint sub-sample. This sample of optically selected TDEs is the largest sample of TDEs from a single survey yet, and the systematic discovery, classification, and follow-up of this sample allows for robust characterization of TDE properties, an important stepping stone looking forward toward the Rubin era., Comment: 44 pages, 22 figures, 9 tables, accepted to ApJ

We demonstrate a fully-integrated multipurpose microwave frequency identification system on silicon-on-insulator platform. Thanks to its multipurpose features, the chip is able to identify different types of microwave signals, including single-frequency, multiple-frequency, chirped and frequency-hopping microwave signals, as well as discriminate instantaneous frequency variation among the frequency-modulated signals. This demonstration exhibits fully integrated solution and fully functional microwave frequency identification, which can meet the requirements in reduction of size, weight and power for future advanced microwave photonic processor., Comment: 23 pages,6 figures

Dirty fireballs are a hypothesized class of relativistic massive-star explosions with an initial Lorentz factor $\Gamma_\mathrm{init}$ below the $\Gamma_\mathrm{init}\sim100$ required to produce a long-duration gamma-ray burst (LGRB), but which could still produce optical emission resembling LGRB afterglows. Here we present the results of a search for on-axis optical afterglows using the Zwicky Transient Facility (ZTF). Our search yielded seven optical transients that resemble on-axis LGRB afterglows in terms of their red colors ($g-r>0$ mag), faint host galaxy ($r>23$ mag), and rapid fading ($dr/dt>1$ mag/day). Spectroscopy of the transient emission within a few days of discovery established cosmological distances ($z=0.876$ to $z=2.9$) for six events, tripling the number of afterglows with redshift measurements discovered by optical surveys without a $\gamma$-ray trigger. Upon a retrospective search, four events (ZTF20abbiixp/AT2020kym, ZTF21aagwbjr/AT2021buv, ZTF21aakruew/AT2021cwd, ZTF21abfmpwn/AT2021qbd) turned out to have a likely associated LGRB (GRB200524A, GRB210204A, GRB210212B, GRB210610B), while three did not (ZTF20aajnksq/AT2020blt, ZTF21aaeyldq/AT2021any, ZTF21aayokph/AT2021lfa). Our search revealed no definitive new class of events: the simplest explanation for the apparently "orphan" events is that they were regular LGRBs missed by high-energy satellites due to detector sensitivity and duty cycle, although it is possible that they were intrinsically faint in $\gamma$-rays or viewed slightly off-axis. We rule out a scenario in which dirty fireballs have a similar energy per solid angle to LGRBs and are an order of magnitude more common. In addition, we set the first direct constraint on the ratio of the opening angles of the material producing $\gamma$-rays and the material producing early optical afterglow emission, finding that they must be comparable., Comment: 37 pages, 12 figures. Accepted for publication in ApJ

We present AT2020mrf (SRGe J154754.2$+$443907), an extra-galactic ($z=0.1353$) fast blue optical transient (FBOT) with a rise time of $t_{g,\rm rise}=3.7$ days and a peak luminosity of $M_{g,\rm peak}=-20.0$. Its optical spectrum around peak shows a broad ($v\sim0.1c$) emission feature on a blue continuum ($T\sim2\times10^4$ K), which bears a striking resemblance to AT2018cow. Its bright radio emission ($\nu L_\nu = 1.2\times 10^{39}\,{\rm erg\,s^{-1}}$; $\nu_{\rm rest}= 7.4$ GHz; 261 days) is similar to four other AT2018cow-like events, and can be explained by synchrotron radiation from the interaction between a sub-relativistic ($\gtrsim0.07$-$0.08c$) forward shock and a dense environment ($\dot M \lesssim 10^{-3}\,M_\odot \,{\rm yr^{-1}}$ for $v_{\rm w}=10^3\,{\rm km\,s^{-1}}$). AT2020mrf occurs in a galaxy with $M_\ast \sim 10^8\,M_\odot$ and specific star formation rate $\sim 10^{-10}\, {\rm yr^{-1}}$, supporting the idea that AT2018cow-like events are preferentially hosted by dwarf galaxies. The X-ray luminosity of AT2020mrf is the highest among FBOTs. At 35-37 days, SRG/eROSITA detected luminous ($L_{\rm X}\sim 2\times 10^{43}\,{\rm erg\,s^{-1}}$; 0.3-10 keV) X-ray emission. The X-ray spectral shape ($f_\nu \propto \nu^{-0.8}$) and erratic intraday variability are reminiscent of AT2018cow, but the luminosity is a factor of $\sim20$ greater than AT2018cow. At 328 days, Chandra detected it at $L_{\rm X}\sim10^{42}\,{\rm erg\,s^{-1}}$, which is $>200$ times more luminous than AT2018cow and CSS161010. At the same time, the X-ray emission remains variable on the timescale of $\sim1$ day. We show that a central engine, probably a millisecond magnetar or an accreting black hole, is required to power the explosion. We predict the rates at which events like AT2018cow and AT2020mrf will be detected by SRG and Einstein Probe., Comment: 26 pages, 18 figures, submitted to ApJ

UVEX is a proposed medium class Explorer mission designed to provide crucial missing capabilities that will address objectives central to a broad range of modern astrophysics. The UVEX design has two co-aligned wide-field imagers operating in the FUV and NUV and a powerful broadband medium resolution spectrometer. In its two-year baseline mission, UVEX will perform a multi-cadence synoptic all-sky survey 50/100 times deeper than GALEX in the NUV/FUV, cadenced surveys of the Large and Small Magellanic Clouds, rapid target of opportunity followup, as well as spectroscopic followup of samples of stars and galaxies. The science program is built around three pillars. First, UVEX will explore the low-mass, low-metallicity galaxy frontier through imaging and spectroscopic surveys that will probe key aspects of the evolution of galaxies by understanding how star formation and stellar evolution at low metallicities affect the growth and evolution of low-metallicity, low-mass galaxies in the local universe. Such galaxies contain half the mass in the local universe, and are analogs for the first galaxies, but observed at distances that make them accessible to detailed study. Second, UVEX will explore the dynamic universe through time-domain surveys and prompt spectroscopic followup capability will probe the environments, energetics, and emission processes in the early aftermaths of gravitational wave-discovered compact object mergers, discover hot, fast UV transients, and diagnose the early stages of stellar explosions. Finally, UVEX will become a key community resource by leaving a large all-sky legacy data set, enabling a wide range of scientific studies and filling a gap in the new generation of wide-field, sensitive optical and infrared surveys provided by the Rubin, Euclid, and Roman observatories. This paper discusses the scientific potential of UVEX, and the broad scientific program., Comment: 69 pages, 43 figures

We present observations of SN 2021csp, the second example of a newly-identified type of supernova (Type Icn) hallmarked by strong, narrow, P Cygni carbon features at early times. The SN appears as a fast and luminous blue transient at early times, reaching a peak absolute magnitude of -20 within 3 days due to strong interaction between fast SN ejecta (v ~ 30000 km/s) and a massive, dense, fast-moving C/O wind shed by the WC-like progenitor months before explosion. The narrow line features disappear from the spectrum 10-20 days after explosion and are replaced by a blue continuum dominated by broad Fe features, reminiscent of Type Ibn and IIn supernovae and indicative of weaker interaction with more extended H/He-poor material. The transient then abruptly fades ~60 days post-explosion when interaction ceases. Deep limits at later phases suggest minimal heavy-element nucleosynthesis, a low ejecta mass, or both, and imply an origin distinct from that of classical Type Ic supernovae. We place SN 2021csp in context with other fast-evolving interacting transients, and discuss various progenitor scenarios: an ultrastripped progenitor star, a pulsational pair-instability eruption, or a jet-driven fallback supernova from a Wolf-Rayet star. The fallback scenario would naturally explain the similarity between these events and radio-loud fast transients, and suggests a picture in which most stars massive enough to undergo a WR phase collapse directly to black holes at the end of their lives., Comment: Accepted to ApJ 29th December 2021. Originally submitted to ApJ on 6th August 2021

The origin of cosmic high-energy neutrinos remains largely unexplained. For high-energy neutrino alerts from IceCube, a coincidence with time-variable emission has been seen for three different types of accreting black holes: (1) a gamma-ray flare from a blazar (TXS 0506+056), (2) an optical transient following a stellar tidal disruption event (TDE; AT2019dsg), and (3) an optical outburst from an active galactic nucleus (AGN; AT2019fdr). For the latter two sources, infrared follow-up observations revealed a powerful reverberation signal due to dust heated by the flare. This discovery motivates a systematic study of neutrino emission from all supermassive black hole with similar dust echoes. Because dust reprocessing is agnostic to the origin of the outburst, our work unifies TDEs and high-amplitude flares from AGN into a population that we dub accretion flares. Besides the two known events, we uncover a third flare that is coincident with a PeV-scale neutrino (AT2019aalc). Based solely on the optical and infrared properties, we estimate a significance of 3.6$\sigma$ for this association of high-energy neutrinos with three accretion flares. Our results imply that at least ~10% of the IceCube high-energy neutrino alerts could be due to accretion flares. This is surprising because the sum of the fluence of these flares is at least three orders of magnitude lower compared to the total fluence of normal AGN. It thus appears that the efficiency of high-energy neutrino production in accretion flares is increased compared to non-flaring AGN. We speculate that this can be explained by the high Eddington ratio of the flares., Comment: Accepted for publication in MNRAS

Observations of the extragalactic ($z=0.0141$) transient AT2018cow established a new class of energetic explosions shocking a dense medium, which produce luminous emission at millimeter and sub-millimeter wavelengths. Here we present detailed millimeter- through centimeter-wave observations of a similar transient, ZTF20acigmel (AT2020xnd) at $z=0.2433$. Using observations from the NOrthern Extended Millimeter Array and the Very Large Array, we model the unusual millimeter and radio emission from AT2020xnd under several different assumptions, and ultimately favor synchrotron radiation from a thermal electron population (relativistic Maxwellian). The thermal-electron model implies a fast but sub-relativistic ($v\approx0.3c$) shock and a high ambient density ($n_e\approx4\times10^{3}$cm$^{-3}$ at $\Delta t\approx40$ days). The X-ray luminosity of $L_X\approx10^{43}$ erg sec$^{-1}$ exceeds simple predictions from the radio and UVOIR luminosity and likely has a separate physical origin, such as a central engine. Using the fact that month-long luminous ($L_\nu\approx 2\times10^{30}$ erg sec$^{-1}$ Hz$^{-1}$ at 100 GHz) millimeter emission appears to be a generic feature of transients with fast ($t_{1/2}\approx3$ days) and luminous ($M_\mathrm{peak}\approx -21 $mag) optical light curves, we estimate the rate at which transients like AT2018cow and AT2020xnd will be detected by future wide-field millimeter transient surveys such as CMB-S4, and conclude that energetic explosions in dense environments may represent a significant population of extragalactic transients in the 100 GHz sky., Comment: 27 pages, 13 figures. Accepted for publication in ApJ on 23 Jan 2022

We present the discovery of ZTF 21aaoryiz/SN 2021fcg -- an extremely low-luminosity Type Iax supernova. SN 2021fcg was discovered by the Zwicky Transient Facility in the star-forming galaxy IC0512 at a distance of $\approx$ 27 Mpc. It reached a peak absolute magnitude of $M_{r} =$ $-12.66\pm0.20$ mag, making it the least luminous thermonuclear supernova discovered to date. The E(B-V) contribution from the underlying host galaxy is unconstrained. However, even if it were as large as 0.5 mag, the peak absolute magnitude would be $M_{r} = -13.78\pm0.20$ mag -- still consistent with being the lowest luminosity SN. Optical spectra of SN 2021fcg taken at 37 and 65 days post maximum show strong [Ca II], Ca II and Na I D emission and several weak [Fe II] emission lines. The [Ca II] emission in the two spectra has extremely low velocities of $\approx 1300$ and $1000$ km s$^{-1}$ respectively. The spectra very closely resemble those of the very low luminosity Type Iax supernovae SN 2008ha, SN 2010ae and SN 2019gsc taken at similar phases. The peak bolometric luminosity of SN 2021fcg is $\approx$ $2.5^{+1.5}_{-0.3}\times10^{40}$ erg s$^{-1}$ which is a factor of three lower than that for SN 2008ha. The bolometric lightcurve of SN 2021fcg is consistent with a very low ejected nickel mass (M$_{\rm{Ni}} \approx 0.8^{+0.4}_{-0.5}\times10^{-3}$ M$_{\odot}$). The low luminosity and nickel mass of SN 2021fcg pose a challenge to the picture that low luminosity SNe Iax originate from deflagrations of near M$_{\rm{ch}}$ hybrid carbon-oxygen-neon white dwarfs. Instead, the merger of a carbon-oxygen and oxygen-neon white dwarf is a promising model to explain SN 2021fcg., Comment: accepted for publication in ApJL. 10 pages, 5 figures

We present a search for extragalactic fast blue optical transients (FBOTs) during Phase I of the Zwicky Transient Facility (ZTF). We identify 38 candidates with durations above half-maximum light 1 d < t1/2 < 12 d, of which 28 have blue (g-r<-0.2 mag) colors at peak light. Of the 38 transients (28 FBOTs), 19 (13) can be spectroscopically classified as core-collapse supernovae (SNe): 11 (8) H- or He-rich (Type II/IIb/Ib) SNe, 6 (4) interacting (Type IIn/Ibn) SNe, and 2 (1) H&He-poor (Type Ic/Ic-BL) SNe. Two FBOTs (published previously) had high-S/N predominantly featureless spectra and luminous radio emission: AT2018lug and AT2020xnd. Seven (five) did not have a definitive classification: AT 2020bdh showed tentative broad H$\alpha$ in emission, and AT 2020bot showed unidentified broad features and was 10 kpc offset from the center of an early-type galaxy. Ten (six) have no spectroscopic observations or redshift measurements. We present multiwavelength (radio, millimeter, and/or X-ray) observations for five FBOTs (three Type Ibn, one Type IIn/Ibn, one Type IIb). Additionally, we search radio-survey (VLA and ASKAP) data to set limits on the presence of radio emission for 22 of the transients. All X-ray and radio observations resulted in non-detections; we rule out AT2018cow-like X-ray and radio behavior for five FBOTs and more luminous emission (such as that seen in the Camel) for four additional FBOTs. We conclude that exotic transients similar to AT2018cow, the Koala, and the Camel represent a rare subset of FBOTs, and use ZTF's SN classification experiments to measure the rate to be at most 0.1% of the local core-collapse SN rate., Comment: Replaced following peer-review process. 46 pages, 20 figures. Accepted for publication in ApJ

While optical surveys regularly discover slow transients like supernovae on their own, the most common way to discover extragalactic fast transients, fading away in a few nights, is via follow-up observations of gamma-ray burst and gravitational-wave triggers. However, wide-field surveys have the potential to also identify rapidly fading transients independently of such external triggers. The volumetric survey speed of the Zwicky Transient Facility (ZTF) makes it sensitive to faint and fast-fading objects as kilonovae, the optical counterparts to binary neutron stars and neutron star-black hole mergers, out to almost 200Mpc. We introduce an open-source software infrastructure, the ZTF REaltime Search and Triggering, ZTFReST, designed to identify kilonovae and fast optical transients in ZTF data. Using the ZTF alert stream combined with forced photometry, we have implemented automated candidate ranking based on their photometric evolution and fitting to kilonova models. Automated triggering of follow-up systems, such as Las Cumbres Observatory, has also been implemented. In 13 months of science validation, we found several extragalactic fast transients independent of any external trigger (though some counterparts were identified later), including at least one supernova with post-shock cooling emission, two known afterglows with an associated gamma-ray burst, two known afterglows without any known gamma-ray counterpart, and three new fast-declining sources (ZTF20abtxwfx, ZTF20acozryr, and ZTF21aagwbjr) that are likely associated with GRB200817A, GRB201103B, and GRB210204A. However, we have not found any objects which appear to be kilonovae; therefore, we constrain the rate of GW170817-like kilonovae to $R < 900$Gpc$^{-3}$yr$^{-1}$. A framework such as ZTFReST could become a prime tool for kilonova and fast transient discovery with the Vera C. Rubin Observatory., Comment: Submitted for publication

We study a class of McKean-Vlasov type stochastic differential equations (SDEs) which arise from the random vortex dynamics and other physics models. By introducing a new approach we resolve the existence and uniqueness of both the weak and strong solutions for the McKean-Vlasov stochastic differential equations whose coefficients are defined in terms of singular integral kernels such as the Biot-Savart kernel. These SDEs which involve the distributions of solutions are in general not Lipschitz continuous with respect to the usual distances on the space of distributions such as the Wasserstein distance. Therefore there is an obstacle in adapting the ordinary SDE method for the study of this class of SDEs, and the conventional methods seem not appropriate for dealing with such distributional SDEs which appear in applications such as fluid mechanics., Comment: 24 pages

The many unusual properties of the enigmatic AT2018cow suggested that at least some subset of the empirical class of fast blue optical transients (FBOTs) represents a genuinely new astrophysical phenomenon. Unfortunately, the intrinsic rarity and fleeting nature of these events have made it difficult to identify additional examples early enough to acquire the observations necessary to constrain theoretical models. We present here the Zwicky Transient Facility discovery of AT2020xnd (ZTF20acigmel, the "Camel") at z=0.243, the first unambiguous AT2018cow analog to be found and confirmed in real time. AT2018cow and AT2020xnd share all key observational properties: a fast optical rise, sustained high photospheric temperature, absence of a second peak attributable to ejection of a radioactively-heated stellar envelope, extremely luminous radio, millimetre, and X-ray emission, and a dwarf-galaxy host. This supports the argument that AT2018cow-like events represent a distinct phenomenon from slower-evolving radio-quiet supernovae, likely requiring a different progenitor or a different central engine. The sample properties of the four known members of this class to date disfavour tidal disruption models but are consistent with the alternative model of an accretion powered jet following the direct collapse of a massive star to a black hole. Contextual filtering of alert streams combined with rapid photometric verification using multi-band imaging provides an efficient way to identify future members of this class, even at high redshift., Comment: Accepted to MNRAS

We observed Episodically Active Asteroid (6478) Gault in 2020 with multiple telescopes in Asia and North America and have found that it is no longer active after its recent outbursts at the end of 2018 and start of 2019. The inactivity during this apparation allowed us to measure the absolute magnitude of Gault of H_r = 14.63 +/- 0.02, G_r = 0.21 +/- 0.02 from our secular phasecurve observations. In addition, we were able to constrain Gault's rotation period using time-series photometric lightcurves taken over 17 hours on multiple days in 2020 August, September and October. The photometric lightcurves have a repeating $\lesssim$0.05 magnitude feature suggesting that (6478) Gault has a rotation period of ~2.5 hours and may have a semi-spherical or top-like shape, much like Near-Earth Asteroids Ryugu and Bennu. The rotation period of ~2.5 hours is near to the expected critical rotation period for an asteroid with the physical properties of (6478) Gault suggesting that its activity observed over multiple epochs is due to surface mass shedding from its fast rotation spun up by the Yarkovsky-O'Keefe-Radzievskii-Paddack effect., Comment: 23 pages, 18 figures; Accepted by ApJ Letters

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

AT2019wey is a Galactic low mass X-ray binary with a candidate black hole accretor first discovered as an optical transient by ATLAS in December 2019. It was then associated with an X-ray source discovered by SRG in March 2020. After observing a brightening in X-rays in August 2020, VLA observations of the source revealed an optically thin spectrum that subsequently shifted to optically thick, as the source continued to brighten in radio. This motivated observations of the source with the VLBA. We found a resolved source that we interpret to be a steady compact jet, a feature associated with black hole X-ray binary systems in the hard X-ray spectral state. The jet power is comparable to the accretion-disk X-ray luminosity. Here, we summarize the results from these observations., Comment: 7 pages, 5 figures, accepted to ApJL

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