Researcher profile

Eric C. Bellm

Eric C. Bellm contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
33works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

33 published item(s)

preprint2022arXiv

A 62-minute orbital period black widow binary in a wide hierarchical triple

Over a dozen millisecond pulsars are ablating low-mass companions in close binary systems. In the original "black widow", the 8-hour orbital period eclipsing pulsar PSR J1959+2048 (PSR B1957+20), high energy emission originating from the pulsar is irradiating and may eventually destroy a low-mass companion. These systems are not only physical laboratories that reveal the dramatic result of exposing a close companion star to the relativistic energy output of a pulsar, but are also believed to harbour some of the most massive neutron stars, allowing for robust tests of the neutron star equation of state. Here, we report observations of ZTF J1406+1222, a wide hierarchical triple hosting a 62-minute orbital period black widow candidate whose optical flux varies by a factor of more than 10. ZTF J1406+1222 pushes the boundaries of evolutionary models, falling below the 80 minute minimum orbital period of hydrogen-rich systems. The wide tertiary companion is a rare low metallicity cool subdwarf star, and the system has a Galactic halo orbit consistent with passing near the Galactic center, making it a probe of formation channels, neutron star kick physics, and binary evolution.

preprint2022arXiv

A Uniform Type Ia Supernova Distance Ladder with the Zwicky Transient Facility: Absolute Calibration Based on the Tip of the Red Giant Branch (TRGB) Method

The current Cepheid-calibrated distance ladder measurement of $H_0$ is reported to be in tension with the values inferred from the cosmic microwave background (CMB), assuming standard cosmology. However, some tip of the red giant branch (TRGB) estimates report $H_0$ in better agreement with the CMB. Hence, it is critical to reduce systematic uncertainties in local measurements to understand the Hubble tension. In this paper, we propose a uniform distance ladder between the second and third rungs, combining SNe~Ia observed by the Zwicky Transient Facility (ZTF) with a TRGB calibration of their absolute luminosity. A large, volume-limited sample of both calibrator and Hubble flow SNe~Ia from the \emph{same} survey minimizes two of the largest sources of systematics: host-galaxy bias and non-uniform photometric calibration. We present results from a pilot study using existing TRGB distance to the host galaxy of ZTF SN~Ia SN 2021rhu (aka ZTF21abiuvdk) in NGC7814. Combining the ZTF calibrator with a volume-limited sample from the first data release of ZTF Hubble flow SNe~Ia, we infer $H_0 = 76.94 \pm 6.4\, {\rm km}\,{\rm s^{-1}}\,{\rm Mpc^{-1}}$, an $8.3 \%$ measurement. The error budget is dominated by the single object calibrating the SN~Ia luminosity in this pilot study. However, the ZTF sample includes already five other SNe~Ia within $\sim$ 20 Mpc for which TRGB distances can be obtained with HST. Finally, we present the prospects of building this distance ladder out to 80 Mpc with JWST observations of more than one hundred ZTF SNe~Ia.

preprint2022arXiv

Candidate Tidal Disruption Event AT2019fdr Coincident with a High-Energy Neutrino

The origins of the high-energy cosmic neutrino flux remain largely unknown. Recently, one high-energy neutrino was associated with a tidal disruption event (TDE). Here we present AT2019fdr, an exceptionally luminous TDE candidate, coincident with another high-energy neutrino. Our observations, including a bright dust echo and soft late-time X-ray emission, further support a TDE origin of this flare. The probability of finding two such bright events by chance is just 0.034%. We evaluate several models for neutrino production and show that AT2019fdr is capable of producing the observed high-energy neutrino, reinforcing the case for TDEs as neutrino sources.

preprint2022arXiv

Discovery of a double detonation thermonuclear supernova progenitor

We present the discovery of a new double detonation progenitor system consisting of a hot subdwarf B (sdB) binary with a white dwarf companion with an P=76.34179(2) min orbital period. Spectroscopic observations are consistent with an sdB star during helium core burning residing on the extreme horizontal branch. Chimera light curves are dominated by ellipsoidal deformation of the sdB star and a weak eclipse of the companion white dwarf. Combining spectroscopic and light curve fits we find a low mass sdB star, $M_{\rm sdB}=0.383\pm0.028$ M$_\odot$ with a massive white dwarf companion, $M_{\rm WD}=0.725\pm0.026$ M$_\odot$. From the eclipses we find a blackbody temperature for the white dwarf of 26,800 K resulting in a cooling age of $\approx$25 Myrs whereas our MESA model predicts an sdB age of $\approx$170 Myrs. We conclude that the sdB formed first through stable mass transfer followed by a common envelope which led to the formation of the white dwarf companion $\approx$25 Myrs ago. Using the MESA stellar evolutionary code we find that the sdB star will start mass transfer in $\approx$6 Myrs and in $\approx$60 Myrs the white dwarf will reach a total mass of $0.92$ M$_\odot$ with a thick helium layer of $0.17$ M$_\odot$. This will lead to a detonation that will likely destroy the white dwarf in a peculiar thermonuclear supernova. PTF1 2238+7430 is only the second confirmed candidate for a double detonation thermonuclear supernova. Using both systems we estimate that at least $\approx$1% of white dwarf thermonuclear supernovae originate from sdB+WD binaries with thick helium layers, consistent with the small number of observed peculiar thermonuclear explosions.

preprint2022arXiv

In search of short gamma-ray burst optical counterpart with the Zwicky Transient Facility

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.

preprint2022arXiv

The Final Season Reimagined: 30 Tidal Disruption Events from the ZTF-I Survey

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.

preprint2022arXiv

The Type Icn SN 2021csp: Implications for the Origins of the Fastest Supernovae and the Fates of Wolf-Rayet Stars

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.

preprint2022arXiv

Variability-selected intermediate mass black hole candidates in dwarf galaxies from ZTF and WISE

While it is difficult to observe the first black hole seeds in the early Universe, we can study intermediate mass black holes (IMBHs) in local dwarf galaxies for clues about their origins. In this paper we present a sample of variability--selected AGN in dwarf galaxies using optical photometry from the Zwicky Transient Facility (ZTF) and forward--modeled mid--IR photometry of time--resolved Wide--field Infrared Survey Explorer ({\it WISE}) coadded images. We found that 44 out of 25,714 dwarf galaxies had optically variable AGN candidates, and 148 out of 79,879 dwarf galaxies had mid--IR variable AGN candidates, corresponding to active fractions of $0.17\pm0.03$\% and $0.19\pm0.02$\% respectively. We found that spectroscopic approaches to AGN identification would have missed 81\% of our ZTF IMBH candidates and 69\% of our {\it WISE} IMBH candidates. Only $9$ candidates have been detected previously in radio, X-ray, and variability searches for dwarf galaxy AGN. The ZTF and {\it WISE} dwarf galaxy AGN with broad Balmer lines have virial masses down to $10^{5.5}M_\odot$ and for the rest of the sample, BH masses predicted from host galaxy mass range between $10^{5.2}M_\odot<M_{\text{BH}}<10^{7.3}M_\odot$. We found that only 5 of 152 previously reported variability--selected AGN candidates from the Palomar Transient Factory in common with our parent sample were variable in ZTF. We also determined a nuclear supernova fraction of $0.05\pm0.01$\% year$^{-1}$ for dwarf galaxies in ZTF. Our ZTF and {\it WISE} IMBH candidates show the promise of variability searches for the discovery of otherwise hidden low mass AGN.

preprint2021arXiv

Blazar variability with the Vera C. Rubin Legacy Survey of Space and Time (LSST)

With their emission mainly coming from a relativistic jet pointing towards us, blazars are fundamental sources to study extragalactic jets and their central engines, consisting of supermassive black holes (SMBHs) fed by accretion discs. They are also candidate sources of high-energy neutrinos and cosmic rays. Because of the jet orientation, the non-thermal blazar emission is Doppler beamed; its variability is unpredictable and occurs on time-scales from less than one hour to years. The comprehension of the diverse mechanisms producing the flux and spectral changes requires well-sampled multiband light curves on long time periods. In particular, outbursts are the best test bench to shed light on the underlying physics, especially when studied in a multiwavelength context. The Vera C. Rubin Legacy Survey of Space and Time (Rubin-LSST) will monitor the southern sky for ten years in six photometric bands, offering a formidable tool to study blazar variability features in a statistical way. The alert system will allow us to trigger follow-up observations of outstanding events, especially at high (keV-to-GeV) and very high (TeV) energies. We here examine the simulated Rubin-LSST survey strategies with the aim of understanding which cadences are more suitable for the blazar variability science. Our metrics include light curve and colour sampling. We also investigate the problem of saturation, which will affect the brightest and many flaring sources, and will have a detrimental impact on follow-up observations.

preprint2021arXiv

Cataclysmic Variables in the Second Year of the Zwicky Transient Facility

Using a filter in the GROWTH Marshal based on color and the amplitude and the timescale of variability, we have identified 372 objects as known or candidate cataclysmic variables (CVs) during the second year of operation of the Zwicky Transient Facility (ZTF). From the available difference imaging data, we found that 93 are previously confirmed CVs, and 279 are strong candidates. Spectra of four of the candidates confirm them as CVs by the presence of Balmer emission lines, while one of the four has prominent HeII lines indicative of containing a magnetic white dwarf. Gaia EDR3 parallaxes are available for 154 of these systems, resulting in distances from 108-2096 pc and absolute magnitudes in the range of 7.5-15.0, with the largest number of candidates between 10.5-12.5. The total numbers are 21% higher than from the previous year of the survey with a greater number of distances available but a smaller percentage of systems close to the Galactic plane. Comparison of these findings with a machine learning method of searching all the light curves reveals large differences in each dataset related to the parameters involved in the search process.

preprint2021arXiv

Give Me a Few Hours: Exploring Short Timescales in Rubin Observatory Cadence Simulations

The limiting temporal resolution of a time-domain survey in detecting transient behavior is set by the time between observations of the same sky area. We analyze the distribution of visit separations for a range of Vera C. Rubin Observatory cadence simulations. Current simulations are strongly peaked at the 22 minute visit pair separation and provide effectively no constraint on temporal evolution within the night. This choice will necessarily prevent Rubin from discovering a wide range of astrophysical phenomena in time to trigger rapid followup. We present a science-agnostic metric to supplement detailed simulations of fast-evolving transients and variables and suggest potential approaches for improving the range of timescales explored.

preprint2021arXiv

Initial Characterization of Active Transitioning Centaur, P/2019 LD2 (ATLAS), using Hubble, Spitzer, ZTF, Keck, APO and GROWTH Visible & Infrared Imaging and Spectroscopy

We present visible and mid-infrared imagery and photometry of temporary Jovian co-orbital comet P/2019 LD$_2$ taken with HST/WFC3, Spitzer/IRAC, the GROWTH telescope network, visible spectroscopy from Keck/LRIS and archival ZTF observations taken between 2019 April and 2020 August. Our observations indicate that the nucleus of LD$_2$ has a radius between 0.2-1.8 km assuming a 0.08 albedo and a coma dominated by $\sim$100$μ$ m-scale dust ejected at $\sim$1 m/s speeds with a $\sim$1&#39;&#39; jet pointing in the SW direction. LD$_2$ experienced a total dust mass loss of $\sim$10$^8$ kg at a loss rate of $\sim$6 kg/s with Af$ρ$/cross-section varying between $\sim$85 cm/125 km$^2$ and $\sim$200 cm/310 km$^2$ from 2019 April 9 to 2019 Nov 8. If the increase in Af$ρ$/cross-section remained constant, it implies LD$_2$&#39;s activity began $\sim$2018 November when within 4.8 au of the Sun, implying the onset of H$_2$O sublimation. We measure CO/CO$_2$ gas production of $\lesssim$10$^{27}$ mol/s /$\lesssim$10$^{26}$ mol/s from our 4.5 $μ$m Spitzer observations, $g$-$r$ = 0.59$\pm$0.03, $r$-$i$ = 0.18$\pm$0.05, $i$-$z$ = 0.01$\pm$0.07 from GROWTH observations, H$_2$O gas production of $\lesssim$80 kg/s scaling from our estimated $C_2$ production of $Q_{C_2}\lesssim$7.5$\times10^{24}$ mol/s from Keck/LRIS spectroscopy. We determine that the long-term orbit of LD$_2$ is similar to Jupiter family comets having close encounters with Jupiter within $\sim$0.5 Hill radius in the last $\sim$3 y, within 0.8 Hill radius in $\sim$9 y. Additionally, 78.8$\%$ of our orbital clones are ejected from the Solar System within $1 \times 10^{6}$ years having a dynamical half-life of 3.4 $\times 10^5$ years.

preprint2021arXiv

Microlensing Events in the Galactic Plane Using the Zwicky Transient Facility

Microlensing is a powerful technique to study the Galactic population of &#34;dark&#34; objects such as exoplanets both bound and unbound, brown dwarfs, low-luminosity stars, old white dwarfs, neutron stars, and almost the only way to study isolated stellar-mass black holes. The majority of previous efforts to search for gravitational microlensing events have concentrated towards high-density fields such as the Galactic bulge. Microlensing events in the Galactic plane have the advantage of closer proximity and better constrained relative proper motions, leading to better constrained lens mass estimates at the expense of a lower optical depth compared to events towards the Galactic bulge. We use the Zwicky Transient Facility (ZTF) Data Release 5 (DR5) compiled from 2018--2021 to survey the Galactic plane in the region of $|b| < 20^\circ$. We find a total of 60 candidate microlensing events including three that show a strong microlensing parallax effect. The rate of events traces Galactic structure, decreasing exponentially as a function Galactic longitude with scale length $\ell_0 \sim 37^\circ$. On average, we find Einstein timescales of our microlensing events to be about three times as long ($\sim60$ days) compared to those towards the Galactic bulge ($\sim20$ days). This pilot project demonstrates that microlensing towards the Galactic plane shows strong promise for characterization of dark objects within the Galatic disk.

preprint2021arXiv

Optimization of the Observing Cadence for the Rubin Observatory Legacy Survey of Space and Time: a pioneering process of community-focused experimental design

Vera C. Rubin Observatory is a ground-based astronomical facility under construction, a joint project of the National Science Foundation and the U.S. Department of Energy, designed to conduct a multi-purpose 10-year optical survey of the southern hemisphere sky: the Legacy Survey of Space and Time. Significant flexibility in survey strategy remains within the constraints imposed by the core science goals of probing dark energy and dark matter, cataloging the Solar System, exploring the transient optical sky, and mapping the Milky Way. The survey&#39;s massive data throughput will be transformational for many other astrophysics domains and Rubin&#39;s data access policy sets the stage for a huge potential users&#39; community. To ensure that the survey science potential is maximized while serving as broad a community as possible, Rubin Observatory has involved the scientific community at large in the process of setting and refining the details of the observing strategy. The motivation, history, and decision-making process of this strategy optimization are detailed in this paper, giving context to the science-driven proposals and recommendations for the survey strategy included in this Focus Issue.

preprint2021arXiv

Optimizing Cadences with Realistic Light Curve Filtering for Serendipitous Kilonova Discovery with Vera Rubin Observatory

Current and future optical and near-infrared wide-field surveys have the potential of finding kilonovae, the optical and infrared counterparts to neutron star mergers, independently of gravitational-wave or high-energy gamma-ray burst triggers. The ability to discover fast and faint transients such as kilonovae largely depends on the area observed, the depth of those observations, the number of re-visits per field in a given time frame, and the filters adopted by the survey; it also depends on the ability to perform rapid follow-up observations to confirm the nature of the transients. In this work, we assess kilonova detectability in existing simulations of the LSST strategy for the Vera C. Rubin Wide Fast Deep survey, with focus on comparing rolling to baseline cadences. Although currently available cadences can enable the detection of more than 300 kilonovae out to 1400 Mpc over the ten-year survey, we can expect only 3-32 kilonovae similar to GW170817 to be recognizable as fast-evolving transients. We also explore the detectability of kilonovae over the plausible parameter space, focusing on viewing angle and ejecta masses. We find that observations in redder izy bands are crucial for identification of nearby (within 300 Mpc) kilonovae that could be spectroscopically classified more easily than more distant sources. Rubin&#39;s potential for serendipitous kilonova discovery could be increased by gain of efficiency with the employment of individual 30s exposures (as opposed to 2x15s snap pairs), with the addition of red-band observations coupled with same-night observations in g- or r-bands, and possibly with further development of a new rolling-cadence strategy.

preprint2020arXiv

A new class of Roche lobe-filling hot subdwarf binaries

We present the discovery of the second binary with a Roche lobe-filling hot subdwarf transferring mass to a white dwarf (WD) companion. This 56 minute binary was discovered using data from the Zwicky Transient Facility. Spectroscopic observations reveal an He-sdOB star with an effective temperature of $T_{\rm eff}=33,700\pm1000$ K and a surface gravity of $log(g)=5.54\pm0.11$. The GTC+HiPERCAM light curve is dominated by the ellipsoidal deformation of the He-sdOB star and shows an eclipse of the He-sdOB by an accretion disk as well as a weak eclipse of the WD. We infer a He-sdOB mass of $M_{\rm sdOB}=0.41\pm0.04$ M$_\odot$ and a WD mass of $M_{\rm WD}=0.68\pm0.05$ M$_\odot$. The weak eclipses imply a WD black-body temperature of $63,000\pm10,000$ K and a radius $R_{\rm WD}=0.0148\pm0.0020$ M$_\odot$ as expected for a WD of such high temperature. The He-sdOB star is likely undergoing hydrogen shell burning and will continue transferring mass for $\approx1$ Myrs at a rate of $10^{-9} M_\odot {\rm yr}^{-1}$ which is consistent with the high WD temperature. The hot subdwarf will then turn into a WD and the system will merge in $\approx30$ Myrs. We suggest that Galactic reddening could bias discoveries towards preferentially finding Roche lobe-filling systems during the short-lived shell burning phase. Studies using reddening corrected samples should reveal a large population of helium core-burning hot subdwarfs with $T_{\rm eff}\approx25,000$ K in binaries of 60-90 minutes with WDs. Though not yet in contact, these binaries would eventually come into contact through gravitational wave emission and explode as a sub-luminous thermonuclear supernova or evolve into a massive single WD.

preprint2020arXiv

Cataclysmic Variables in the First Year of the Zwicky Transient Facility

Using selection criteria based on amplitude, time and color, we have identified 329 objects as known or candidate cataclysmic variable (CVs) during the first year of testing and operation of the Zwicky Transient Facility (ZTF). Of these, 90 are previously confirmed CVs, 218 are strong candidates based on the shape and color of their light curves obtained during 3-562 days of observations, and the remaining 21 are possible CVs but with too few data points to be listed as good candidates. Almost half the strong candidates are within 10 deg of the galactic plane, in contrast to most other large surveys which have avoided crowded fields. The available Gaia parallaxes are consistent with sampling the low mass transfer CVs, as predicted by population models. Our followup spectra have confirmed Balmer/helium emission lines in 27 objects, with four showing high excitation HeII emission, including candidates for an AM CVn, a polar and an intermediate polar. Our results demonstrate that a complete survey of the galactic plane is needed to accomplish an accurate determination of the number of CVs existing in the Milky Way.

preprint2020arXiv

Characterization of the Nucleus, Morphology and Activity of Interstellar Comet 2I/Borisov by Optical and Near-Infrared GROWTH, Apache Point, IRTF, ZTF and Keck Observations

We present visible and near-infrared photometric and spectroscopic observations of interstellar object 2I/Borisov taken from 2019 September 10 to 2019 November 29 using the GROWTH, the APO ARC 3.5 m and the NASA/IRTF 3.0 m combined with post and pre-discovery observations of 2I obtained by ZTF from 2019 March 17 to 2019 May 5. Comparison with imaging of distant Solar System comets shows an object very similar to mildly active Solar System comets with an out-gassing rate of $\sim$10$^{27}$ mol/sec. The photometry, taken in filters spanning the visible and NIR range shows a gradual brightening trend of $\sim0.03$ mags/day since 2019 September 10 UTC for a reddish object becoming neutral in the NIR. The lightcurve from recent and pre-discovery data reveals a brightness trend suggesting the recent onset of significant H$_2$O sublimation with the comet being active with super volatiles such as CO at heliocentric distances $>$6 au consistent with its extended morphology. Using the advanced capability to significantly reduce the scattered light from the coma enabled by high-resolution NIR images from Keck adaptive optics taken on 2019 October 04, we estimate a diameter of 2I&#39;s nucleus of $\lesssim$1.4 km. We use the size estimates of 1I/&#39;Oumuamua and 2I/Borisov to roughly estimate the slope of the ISO size-distribution resulting in a slope of $\sim$3.4$\pm$1.2, similar to Solar System comets and bodies produced from collisional equilibrium.

preprint2020arXiv

Early Ultra-Violet observations of type IIn supernovae constrain the asphericity of their circumstellar material

We present a survey of the early evolution of 12 Type IIn supernovae (SNe IIn) in the Ultra-Violet (UV) and visible light. We use this survey to constrain the geometry of the circumstellar material (CSM) surrounding SN IIn explosions, which may shed light on their progenitor diversity. In order to distinguish between aspherical and spherical circumstellar material (CSM), we estimate the blackbody radius temporal evolution of the SNe IIn of our sample, following the method introduced by Soumagnac et al. We find that higher luminosity objects tend to show evidence for aspherical CSM. Depending on whether this correlation is due to physical reasons or to some selection bias, we derive a lower limit between 35% and 66% on the fraction of SNe IIn showing evidence for aspherical CSM. This result suggests that asphericity of the CSM surrounding SNe IIn is common - consistent with data from resolved images of stars undergoing considerable mass loss. It should be taken into account for more realistic modelling of these events.

preprint2020arXiv

Kilonova Luminosity Function Constraints based on Zwicky Transient Facility Searches for 13 Neutron Star Mergers

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo&#39;s third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged)

preprint2020arXiv

Optical follow-up of the neutron star-black hole mergers S200105ae and S200115j

LIGO and Virgo&#39;s third observing run (O3) revealed the first neutron star-black hole (NSBH) merger candidates in gravitational waves. These events are predicted to synthesize r-process elements creating optical/near-IR &#34;kilonova&#34; (KN) emission. The joint gravitational-wave (GW) and electromagnetic detection of an NSBH merger could be used to constrain the equation of state of dense nuclear matter, and independently measure the local expansion rate of the universe. Here, we present the optical follow-up and analysis of two of the only three high-significance NSBH merger candidates detected to date, S200105ae and S200115j, with the Zwicky Transient Facility (ZTF). ZTF observed $\sim$\,48\% of S200105ae and $\sim$\,22\% of S200115j&#39;s localization probabilities, with observations sensitive to KNe brighter than $-$17.5\,mag fading at 0.5\,mag/day in g- and r-bands; extensive searches and systematic follow-up of candidates did not yield a viable counterpart. We present state-of-the-art KN models tailored to NSBH systems that place constraints on the ejecta properties of these NSBH mergers. We show that with depths of $\rm m_{\rm AB}\approx 22$ mag, attainable in meter-class, wide field-of-view survey instruments, strong constraints on ejecta mass are possible, with the potential to rule out low mass ratios, high BH spins, and large neutron star radii.

preprint2020arXiv

Seventeen Tidal Disruption Events from the First Half of ZTF Survey Observations: Entering a New Era of Population Studies

While tidal disruption events (TDEs) have long been heralded as laboratories for the study of quiescent black holes, the small number of known TDEs and uncertainties in their emission mechanism have hindered progress towards this promise. Here present 17 new TDEs that have been detected recently by the Zwicky Transient Facility along with Swift UV and X-ray follow-up observations. Our homogeneous analysis of the optical/UV light curves, including 22 previously known TDEs from the literature, reveals a clean separation of light curve properties with spectroscopic class. The TDEs with Bowen fluorescence features in their optical spectra have smaller blackbody radii, as well as longer rise times and higher disruption rates compared to the rest of the sample. The Bowen fluorescence mechanism requires a high density which can be reached at smaller radii, which in turn yields longer diffusion timescales. Thus, the difference in rise times suggests the pre-peak TDE light curves are governed not by the fallback timescale, but instead by the diffusion of photons through the tidal debris. The small subset of TDEs that show only helium emission lines in their spectra have the longest rise times, the highest luminosities and the lowest rates. We also report, for the first time, the detection of soft X-ray flares from a TDE on day timescales. Based on the fact the flares peak at a luminosity similar to the optical/UV blackbody luminosity, we attribute them to brief glimpses through a reprocessing layer that otherwise obscures the inner accretion flow.

preprint2020arXiv

SN2019dge: a Helium-rich Ultra-Stripped Envelope Supernova

We present observations of ZTF18abfcmjw (SN2019dge), a helium-rich supernova with a fast-evolving light curve indicating an extremely low ejecta mass ($\approx 0.3\,M_\odot$) and low kinetic energy ($\approx 1.2\times 10^{50}\,{\rm erg}$). Early-time (<4 d after explosion) photometry reveal evidence of shock cooling from an extended helium-rich envelope of $\sim0.1\,M_\odot$ located at $\sim 3\times 10^{12}\,{\rm cm}$ from the progenitor. Early-time He II line emission and subsequent spectra show signatures of interaction with helium-rich circumstellar material, which extends from $\gtrsim 5\times 10^{13}\,{\rm cm}$ to $\gtrsim 2\times 10^{16}\,{\rm cm}$. We interpret SN2019dge as a helium-rich supernova from an ultra-stripped progenitor, which originates from a close binary system consisting of a mass-losing helium star and a low-mass main sequence star or a compact object (i.e., a white dwarf, a neutron star, or a black hole). We infer that the local volumetric birth rate of 19dge-like ultra-stripped SNe is in the range of 1400--8200$\,{\rm Gpc^{-3}\, yr^{-1}}$ (i.e., 2--12% of core-collapse supernova rate). This can be compared to the observed coalescence rate of compact neutron star binaries that are not formed by dynamical capture.

preprint2020arXiv

TDE Hosts are Green and Centrally Concentrated: Signatures of a Post-Merger System

We study the properties of the galaxies hosting the first 19 tidal disruption events (TDEs) detected with the Zwicky Transient Facility (ZTF) within the context of a carefully constructed, representative host galaxy sample. We find that the ZTF sample of TDE hosts is dominated by compact &#34;green valley&#34; galaxies. After we restrict the comparison sample to galaxies with a similar concentration, as measured by Sersic index, we find this green valley over representation is even larger. That is, concentrated red sequence galaxies are not producing TDEs at elevated levels. We present host galaxy spectra which show that E+A galaxies are overrepresented in the ZTF sample by a factor of $\approx$22, which is lower than previous TDE host galaxy studies have found. We find that this overrepresentation can be fully accounted for when taking into account the masses, colors, and Sérsic indices of the ZTF TDE hosts. The combination of both green colors and high Sérsic index of the typical TDE host galaxy could be explained if the TDE rate is temporarily enhanced following a merger that leads to a higher central concentration of stars.

preprint2020arXiv

The first ultracompact Roche lobe-filling hot subdwarf binary

We report the discovery of the first short period binary in which a hot subdwarf star (sdOB) fills its Roche lobe and started mass transfer to its companion. The object was discovered as part of a dedicated high-cadence survey of the Galactic Plane named the Zwicky Transient Facility and exhibits a period of $P_{\rm orb}=39.3401(1)$ min, making it the most compact hot subdwarf binary currently known. Spectroscopic observations are consistent with an intermediate He-sdOB star with an effective temperature of $T_{\rm eff}=42,400\pm300$ K and a surface gravity of $\log(g)=5.77\pm0.05$. A high-signal-to noise GTC+HiPERCAM light curve is dominated by the ellipsoidal deformation of the sdOB star and an eclipse of the sdOB by an accretion disk. We infer a low-mass hot subdwarf donor with a mass $M_{\rm sdOB}=0.337\pm0.015$ M$_\odot$ and a white dwarf accretor with a mass $M_{\rm WD}=0.545\pm0.020$ M$_\odot$. Theoretical binary modeling indicates the hot subdwarf formed during a common envelope phase when a $2.5-2.8$ M$_\odot$ star lost its envelope when crossing the Hertzsprung Gap. To match its current $P_{\rm orb}$, $T_{\rm eff}$, $\log(g)$, and masses, we estimate a post-common envelope period of $P_{\rm orb}\approx150$ min, and find the sdOB star is currently undergoing hydrogen shell burning. We estimate that the hot subdwarf will become a white dwarf with a thick helium layer of $\approx0.1$ M$_\odot$ and will merge with its carbon/oxygen white dwarf companion after $\approx17$ Myr and presumably explode as a thermonuclear supernova or form an R CrB star.

preprint2020arXiv

ZTF J1901+5309: A 40.6-Minute Orbital Period Eclipsing Double White Dwarf System

The Zwicky Transient Facility has begun to discover binary systems with orbital periods that are less than 1\,hr. Combined with dedicated follow-up systems, which allow for high-cadence photometry of these sources, systematic confirmation and characterization of these sources are now possible. Here, we report the discovery of ZTF J190125.42+530929.5, a 40.6\,min orbital period, eclipsing double white-dwarf binary. Both photometric and spectroscopic modeling confirm its nature, yielding an estimated inclination of $i = 86.2^{+0.6}_{-0.2}\,\rm degrees$ and primary and secondary effective temperatures of $\textrm{T}_\textrm{eff} = 28,000^{+500}_{-500}\,K$ and $\textrm{T}_\textrm{eff} = 17,600^{+400}_{-400}\,K$ respectively. This system adds to a growing list of sources for future gravitational-wave detectors and contributes to the demographic analysis of double degenerates.

preprint2020arXiv

Zwicky Transient Facility constraints on the optical emission from the nearby repeating FRB 180916.J0158+65

The discovery rate of fast radio bursts (FRBs) is increasing dramatically thanks to new radio facilities. Meanwhile, wide-field instruments such as the 47 deg$^2$ Zwicky Transient Facility (ZTF) survey the optical sky to study transient and variable sources. We present serendipitous ZTF observations of the CHIME repeating source FRB 180916.J0158+65, that was localized to a spiral galaxy 149 Mpc away and is the first FRB suggesting periodic modulation in its activity. While 147 ZTF exposures corresponded to expected high-activity periods of this FRB, no single ZTF exposure was at the same time as a CHIME detection. No $>3σ$ optical source was found at the FRB location in 683 ZTF exposures, totalling 5.69 hours of integration time. We combined ZTF upper limits and expected repetitions from FRB 180916.J0158+65 in a statistical framework using a Weibull distribution, agnostic of periodic modulation priors. The analysis yielded a constraint on the ratio between the optical and radio fluences of $η\lesssim 200$, corresponding to an optical energy $E_{\rm opt} \lesssim 3 \times 10^{46}$ erg for a fiducial 10 Jy ms FRB (90% confidence). A deeper (but less statistically robust) constraint of $η\lesssim 3$ can be placed assuming a rate of $r(>5$ Jy ms)= hr$^{-1}$ and $1.2\pm 1.1$ FRB occurring during exposures taken in high-activity windows. The constraint can be improved with shorter per-image exposures and longer integration time, or observing FRBs at higher Galactic latitudes. This work demonstrated how current surveys can statistically constrain multi-wavelength counterparts to FRBs even without deliberately scheduled simultaneous radio observation.

preprint2019arXiv

A Twilight Search for Atiras, Vatiras and Co-orbital Asteroids: Preliminary Results

Near-Earth Objects (NEOs) that orbit the Sun on or within Earth&#39;s orbit are tricky to detect for Earth-based observers due to their proximity to the Sun in the sky. These small bodies hold clues to the dynamical history of the inner solar system as well as the physical evolution of planetesimals in extreme environments. Populations in this region include the Atira and Vatira asteroids, as well as Venus and Earth co-orbital asteroids. Here we present a twilight search for these small bodies, conducted using the 1.2-m Oschin Schmidt and the Zwicky Transient Facility (ZTF) camera at Palomar Observatory. The ZTF twilight survey operates at solar elongations down to $35^\circ$ with limiting magnitude of $r=19.5$. During a total of 40 evening sessions and 62 morning sessions conducted between 2018 November 15 and 2019 June 23, we detected 6 Atiras, including 2 new discoveries 2019 AQ$_3$ and 2019 LF$_6$, but no Vatiras or Earth/Venus co-orbital asteroids. NEO population models show that these new discoveries are likely only the tip of the iceberg, with the bulk of the population yet to be found. The population models also suggest that we have only detected 5--$7\%$ of the $H<20$ Atira population over the 7-month survey. Co-orbital asteroids are smaller in diameters and require deeper surveys. A systematic and efficient survey of the near-Sun region will require deeper searches and/or facilities that can operate at small solar elongations.

preprint2019arXiv

Comet 240P/NEAT is Stirring

Comets are primitive objects that formed in the protoplanetary disk, and have been largely preserved over the history of the Solar System. However, they are not pristine, and surfaces of cometary nuclei do evolve. In order to understand the extent of their primitive nature, we must define the mechanisms that affect their surfaces and comae. We examine the lightcurve of comet 240P/NEAT over three consecutive orbits, and investigate three events of significant brightening ($Δm \sim -2$ mag). Unlike typical cometary outbursts, each of the three events are long-lived, with enhanced activity for at least 3 to 6 months. The third event, observed by the Zwicky Transient Facility, occurred in at least two stages. The anomalous behavior appears to have started after the comet was perturbed by Jupiter in 2007, reducing its perihelion distance from 2.53 to 2.12 au. We suggest that the brightening events are temporary transitions to a higher baseline activity level, brought on by the increased insolation, which has warmed previously insulated sub-surface layers. The new activity is isolated to one or two locations on the nucleus, indicating that the surface or immediate sub-surface is heterogeneous. Further study of this phenomenon may provide insight into cometary outbursts, the structure of the near-surface nucleus, and cometary nucleus mantling.

preprint2019arXiv

Evidence for Late-stage Eruptive Mass-loss in the Progenitor to SN2018gep, a Broad-lined Ic Supernova: Pre-explosion Emission and a Rapidly Rising Luminous Transient

We present detailed observations of ZTF18abukavn (SN2018gep), discovered in high-cadence data from the Zwicky Transient Facility as a rapidly rising ($1.4\pm0.1$ mag/hr) and luminous ($M_{g,\mathrm{peak}}=-20$ mag) transient. It is spectroscopically classified as a broad-lined stripped-envelope supernova (Ic-BL SN). The high peak luminosity ($L_{\mathrm{bol}} \gtrsim 3 \times 10^{44}$ erg $\mathrm{sec}^{-1}$), the short rise time ($t_{\mathrm{rise}}= 3$ days in $g$-band), and the blue colors at peak ($g-r\sim-0.4$) all resemble the high-redshift Ic-BL iPTF16asu, as well as several other unclassified fast transients. The early discovery of SN2018gep (within an hour of shock breakout) enabled an intensive spectroscopic campaign, including the highest-temperature ($T_{\mathrm{eff}}\gtrsim40,000$ K) spectra of a stripped-envelope SN. A retrospective search revealed luminous ($M_g \sim M_r \approx -14\,$mag) emission in the days to weeks before explosion, the first definitive detection of precursor emission for a Ic-BL. We find a limit on the isotropic gamma-ray energy release $E_\mathrm{γ,iso}<4.9 \times 10^{48}$ erg, a limit on X-ray emission $L_{\mathrm{X}} < 10^{40}\,$erg sec$^{-1}$, and a limit on radio emission $νL_ν\lesssim 10^{37}\,$erg sec$^{-1}$. Taken together, we find that the early ($<10\,$days) data are best explained by shock breakout in a massive shell of dense circumstellar material ($0.02\,M_\odot$) at large radii ($3 \times 10^{14}\,$cm) that was ejected in eruptive pre-explosion mass-loss episodes. The late-time ($>10$ days) light curve requires an additional energy source, which could be the radioactive decay of Ni-56.

preprint2019arXiv

GROWTH on S190425z: Searching thousands of square degrees to identify an optical or infrared counterpart to a binary neutron star merger with the Zwicky Transient Facility and Palomar Gattini IR

The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg$^2$ of the initial skymap over the next two nights, corresponding to 46\% integrated probability, ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% integrated probability based on the initial sky map. The revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient &#34;alerts&#34; over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical lightcurves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled-out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible.

preprint2019arXiv

Pre-discovery Activity of New Interstellar Comet 2I/Borisov Beyond 5 AU

Comet 2I/Borisov, the first unambiguous interstellar comet ever found, was discovered in August 2019 at $\sim3$ au from the Sun on its inbound leg. No pre-discovery detection beyond 3 au has yet been reported, mostly due to the comet&#39;s proximity to the Sun as seen from the Earth. Here we present a search for pre-discovery detections of comet Borisov using images taken by the Catalina Sky Survey (CSS), Pan-STARRS and Zwicky Transient Facility (ZTF), with a further comprehensive follow-up campaign being presented in \citet{Bolin2019}. We identified comet Borisov in ZTF images taken in May 2019 and use these data to update its orbit. This allowed us to identify the comet in images acquired as far back as December 2018, when it was 7.8 au from the Sun. The comet was not detected in November 2018 when it was 8.6 au from the Sun, possibly implying an onset of activity around this time. This suggests that the activity of the comet is either driven by a more volatile species other than H$_2$O, such as CO or CO$_2$, or by exothermic crystallization of amorphous ice. We derive the radius of the nucleus to be $<7$ km using the non-detection in November 2018, and estimate an area of $\sim0.5$---$10 \mathrm{km^2}$ has been active between December 2018 and September 2019, though this number is model-dependent and is highly uncertain. The behavior of comet Borisov during its inbound leg is observationally consistent with dynamically new comets observed in our solar system, suggesting some similarities between the two.

preprint2019arXiv

ZTF Early Observations of Type Ia Supernovae I: Properties of the 2018 Sample

Early-time observations of Type Ia supernovae (SNe Ia) are essential to constrain their progenitor properties. In this paper, we present high-quality light curves of 127 SNe Ia discovered by the Zwicky Transient Facility (ZTF) in 2018. We describe our method to perform forced point spread function (PSF) photometry, which can be applied to other types of extragalactic transients. With a planned cadence of six observations per night ($3g+3r$), all of the 127 SNe Ia are detected in both $g$ and $r$ band more than 10\,d (in the rest frame) prior to the epoch of $g$-band maximum light. The redshifts of these objects range from $z=0.0181$ to 0.165; the median redshift is 0.074. Among the 127 SNe, 50 are detected at least 14\,d prior to maximum light (in the rest frame), with a subset of 9 objects being detected more than 17\,d before $g$-band peak. This is the largest sample of young SNe Ia collected to date; it can be used to study the shape and color evolution of the rising light curves in unprecedented detail. We discuss six peculiar events in this sample, including one 02cx-like event ZTF18abclfee (SN\,2018crl), one Ia-CSM SN ZTF18aaykjei (SN\,2018cxk), and four objects with possible super-Chandrasekhar mass progenitors: ZTF18abhpgje (SN\,2018eul), ZTF18abdpvnd (SN\,2018dvf), ZTF18aawpcel (SN\,2018cir) and ZTF18abddmrf (SN\,2018dsx).