Source author record

Eric C. Bellm

Eric C. Bellm appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

44works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

44 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'' 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$'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 "dark" 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's massive data throughput will be transformational for many other astrophysics domains and Rubin's data access policy sets the stage for a huge potential users' 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'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's nucleus of $\lesssim$1.4 km. We use the size estimates of 1I/'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'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'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 "kilonova" (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'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 "green valley" 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'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 "alerts" 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'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).

preprint2016arXiv

Volumetric Survey Speed: A Figure of Merit for Transient Surveys

Time-domain surveys can exchange sky coverage for revisit frequency, complicating the comparison of their relative capabilities. By using different revisit intervals, a specific camera may execute surveys optimized for discovery of different classes of transient objects. We propose a new figure of merit, the instantaneous volumetric survey speed, for evaluating transient surveys. This metric defines the trade between cadence interval and snapshot survey volume and so provides a natural means of comparing survey capability. The related metric of areal survey speed imposes a constraint on the range of possible revisit times: we show that many modern time-domain surveys are limited by the amount of fresh sky available each night. We introduce the concept of "spectroscopic accessibility" and discuss its importance for transient science goals requiring followup observing. We present an extension of the control time algorithm for cases where multiple consecutive detections are required. Finally, we explore how survey speed and choice of cadence interval determine the detection rate of transients in the peak absolute magnitude-decay timescale phase space.

preprint2015arXiv

Properties and Evolution of the Redback Millisecond Pulsar Binary PSR J2129-0429

PSR J2129-0429 is a "redback" eclipsing millisecond pulsar binary with an unusually long 15.2 hour orbit. It was discovered by the Green Bank Telescope in a targeted search of unidentified Fermi gamma-ray sources. The pulsar companion is optically bright (mean $m_R = 16.6$ mag), allowing us to construct the longest baseline photometric dataset available for such a system. We present ten years of archival and new photometry of the companion from LINEAR, CRTS, PTF, the Palomar 60-inch, and LCOGT. Radial velocity spectroscopy using the Double-Beam Spectrograph on the Palomar 200-inch indicates that the pulsar is massive: $1.74\pm0.18 M_\odot$. The G-type pulsar companion has mass $0.44\pm0.04 M_\odot$, one of the heaviest known redback companions. It is currently 95\% Roche-lobe filling and only mildly irradiated by the pulsar. We identify a clear 13.1 mmag yr$^{-1}$ secular decline in the mean magnitude of the companion as well as smaller-scale variations in the optical lightcurve shape. This behavior may indicate that the companion is cooling. Binary evolution calculations indicate that PSR J2129-0429 has an orbital period almost exactly at the bifurcation period between systems that converge into tighter orbits as black widows and redbacks and those that diverge into wider pulsar--white dwarf binaries. Its eventual fate may depend on whether it undergoes future episodes of mass transfer and increased irradiation.

preprint2014arXiv

Confirmation of a High Magnetic Field in GRO J1008-57

GRO J1008$-$57 is a high-mass X-ray binary for which several claims of a cyclotron resonance scattering feature near 80 keV have been reported. We use NuSTAR, Suzaku, and Swift data from its giant outburst of November 2012 to confirm the existence of the 80 keV feature and perform the most sensitive search to date for cyclotron scattering features at lower energies. We find evidence for a 78$^{+3}_{-2}$ keV line in the NuSTAR and Suzaku data at $> 4σ$ significance, confirming the detection using Suzaku alone by Yamamoto et al. (2014). A search of both the time-integrated and phase-resolved data rules out a fundamental at lower energies with optical depth larger than 5% of the 78 keV line. These results indicate that GRO J1008$-$57 has a magnetic field of $6.7\times10^{12} (1+z)$ G, the highest of known accreting pulsars.

preprint2014arXiv

NuSTAR results and future plans for magnetar and rotation-powered pulsar observations

The Nuclear Spectroscopic Telescope Array (NuSTAR) is the first focusing hard X-ray mission in orbit and operates in the 3-79 keV range. NuSTAR's sensitivity is roughly two orders of magnitude better than previous missions in this energy band thanks to its superb angular resolution. Since its launch in 2012 June, NuSTAR has performed excellently and observed many interesting sources including four magnetars, two rotation-powered pulsars and the cataclysmic variable AE Aquarii. NuSTAR also discovered 3.76-s pulsations from the transient source SGR J1745-29 recently found by Swift very close to the Galactic Center, clearly identifying the source as a transient magnetar. For magnetar 1E 1841-045, we show that the spectrum is well fit by an absorbed blackbody plus broken power-law model with a hard power-law photon index of ~1.3. This is consistent with previous results by INTEGRAL and RXTE. We also find an interesting double-peaked pulse profile in the 25-35 keV band. For AE Aquarii, we show that the spectrum can be described by a multi-temperature thermal model or a thermal plus non-thermal model; a multi-temperature thermal model without a non-thermal component cannot be ruled out. Furthermore, we do not see a spiky pulse profile in the hard X-ray band, as previously reported based on Suzaku observations. For other magnetars and rotation-powered pulsars observed with NuSTAR, data analysis results will be soon available.

preprint2014arXiv

The Zwicky Transient Facility

The Zwicky Transient Facility (ZTF) is a next-generation optical synoptic survey that builds on the experience and infrastructure of the Palomar Transient Factory (PTF). Using a new 47 deg$^2$ survey camera, ZTF will survey more than an order of magnitude faster than PTF to discover rare transients and variables. I describe the survey and the camera design. Searches for young supernovae, fast transients, counterparts to gravitational-wave detections, and rare variables will benefit from ZTF's high cadence, wide area survey.

preprint2014arXiv

X-ray Spectral Components in the Afterglow of GRB 130925A

We have identified spectral features in the late-time X-ray afterglow of the unusually long, slow-decaying GRB 130925A using NuSTAR, Swift-XRT, and Chandra. A spectral component in addition to an absorbed power-law is required at $>4σ$ significance, and its spectral shape varies between two observation epochs at $2\times10^5$ and $10^6$ seconds after the burst. Several models can fit this additional component, each with very different physical implications. A broad, resolved Gaussian absorption feature of several keV width improves the fit, but it is poorly constrained in the second epoch. An additive black body or second power-law component provide better fits. Both are challenging to interpret: the blackbody radius is near the scale of a compact remnant ($10^8$ cm), while the second powerlaw component requires an unobserved high-energy cutoff in order to be consistent with the non-detection by Fermi-LAT.

preprint2013arXiv

NuSTAR detection of hard X-ray phase lags from the accreting pulsar GS 0834-430

The NuSTAR hard X-ray telescope observed the transient Be/X-ray binary GS 0834-430 during its 2012 outburst - the first active state of this system observed in the past 19 years. We performed timing and spectral analysis, and measured the X-ray spectrum between 3-79keV with high statistical significance. We find the phase-averaged spectrum to be consistent with that observed in many other magnetized accreting pulsars. We fail to detect cyclotron resonance scattering features that would allow us to constrain the pulsar's magnetic field in either phase-averaged or phase-resolved spectra. Timing analysis shows a clearly detected pulse period of ~12.29s in all energy bands. The pulse profiles show a strong, energy-dependent hard phase lag of up to 0.3 cycles in phase, or about 4s. Such dramatic energy-dependent lags in the pulse profile have never before been reported in high-mass X-ray binary (HMXB) pulsars. Previously reported lags have been significantly smaller in phase and restricted to low-energies (E<10keV). We investigate the possible mechanisms that might produce this energy-dependent pulse phase shift. We find the most likely explanation for this effect to be a complex beam geometry.

preprint2013arXiv

The Nuclear Spectroscopic Telescope Array (NuSTAR) Mission

The Nuclear Spectroscopic Telescope Array (NuSTAR) mission, launched on 13 June 2012, is the first focusing high-energy X-ray telescope in orbit. NuSTAR operates in the band from 3 -- 79 keV, extending the sensitivity of focusing far beyond the ~10 keV high-energy cutoff achieved by all previous X-ray satellites. The inherently low-background associated with concentrating the X-ray light enables NuSTAR to probe the hard X-ray sky with a more than one-hundred-fold improvement in sensitivity over the collimated or coded-mask instruments that have operated in this bandpass. Using its unprecedented combination of sensitivity, spatial and spectral resolution, NuSTAR will pursue five primary scientific objectives, and will also undertake a broad program of targeted observations. The observatory consists of two co-aligned grazing-incidence X-ray telescopes pointed at celestial targets by a three-axis stabilized spacecraft. Deployed into a 600 km, near-circular, 6degree inclination orbit, the Observatory has now completed commissioning, and is performing consistent with pre-launch expectations. NuSTAR is now executing its primary science mission, and with an expected orbit lifetime of ten years, we anticipate proposing a guest investigator program, to begin in Fall 2014.

preprint2013arXiv

The smooth cyclotron line in Her X-1 as seen with NuSTAR

Her X-1, one of the brightest and best studied X-ray binaries, shows a cyclotron resonant scattering feature (CRSF) near 37 keV. This makes it an ideal target for detailed study with the Nuclear Spectroscopic Telescope Array (NuSTAR), taking advantage of its excellent hard X-ray spectral resolution. We observed Her X-1 three times, coordinated with Suzaku, during one of the high flux intervals of its 35d super-orbital period. This paper focuses on the shape and evolution of the hard X-ray spectrum. The broad-band spectra can be fitted with a powerlaw with a high-energy cutoff, an iron line, and a CRSF. We find that the CRSF has a very smooth and symmetric shape, in all observations and at all pulse-phases. We compare the residuals of a line with a Gaussian optical depth profile to a Lorentzian optical depth profile and find no significant differences, strongly constraining the very smooth shape of the line. Even though the line energy changes dramatically with pulse phase, we find that its smooth shape does not. Additionally, our data show that the continuum is only changing marginally between the three observations. These changes can be explained with varying amounts of Thomson scattering in the hot corona of the accretion disk. The average, luminosity-corrected CRSF energy is lower than in past observations and follows a secular decline. The excellent data quality of NuSTAR provides the best constraint on the CRSF energy to date.

preprint2013arXiv

Tracing the Orphan Stream to 55 kpc with RR Lyrae Stars

We report positions, velocities and metallicities of 50 ab-type RR Lyrae (RRab) stars observed in the vicinity of the Orphan stellar stream. Using about 30 RRab stars classified as being likely members of the Orphan stream, we study the metallicity and the spatial extent of the stream. We find that RRab stars in the Orphan stream have a wide range of metallicities, from -1.5 dex to -2.7 dex. The average metallicity of the stream is -2.1 dex, identical to the value obtained by Newberg et al. (2010) using blue horizontal branch stars. We find that the most distant parts of the stream (40-50 kpc from the Sun) are about 0.3 dex more metal-poor than the closer parts (within ~30 kpc), suggesting a possible metallicity gradient along the stream's length. We have extended the previous studies and have mapped the stream up to 55 kpc from the Sun. Even after a careful search, we did not identify any more distant RRab stars that could plausibly be members of the Orphan stream. If confirmed with other tracers, this result would indicate a detection of the end of the leading arm of the stream. We have compared the distances of Orphan stream RRab stars with the best-fit orbits obtained by Newberg et al. (2010). We find that model 6 of Newberg et al. (2010) cannot explain the distances of the most remote Orphan stream RRab stars, and conclude that the best fit to distances of Orphan stream RRab stars and to the local circular velocity is provided by potentials where the total mass of the Galaxy within 60 kpc is M_{60}~2.7x10^{11} Msun, or about 60% of the mass found by previous studies. More extensive modelling that would consider non-spherical potentials and the possibility of misalignment between the stream and the orbit, is highly encouraged.

preprint2009arXiv

RHESSI Tests of Quasi-Thermal Gamma-Ray Burst Spectral Models

Prompt gamma-ray burst spectra evolve on short time scales, suggesting that time-resolved spectral fits may help diagnose the still unknown prompt emission mechanism. We use broad-band gamma-ray data from the RHESSI spacecraft to test quasi-thermal models with high signal-to-noise time-resolved spectra of nine bright gamma-ray bursts. In contrast to results reported in more narrow energy bands, the quality of the fits of quasi-thermal models is poor in relation to fits of the phenomenological Band function. Moreover, the best-fit parameters for the simplest quasi-thermal model, a black body plus a nonthermal power law, show significant dependence on the fit band. Models that replace the power law with more complicated nonthermal functions are not robust for the data considered here and decrease the physical relevance of the fit black body.