BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Matthew J. Graham

Matthew J. Graham contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEastro-ph.GAastro-ph.SRastro-ph.IMastro-ph.EPgr-qcastro-ph.CO

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
38works
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

38 published item(s)

preprint2023arXiv

Science with the Ultraviolet Explorer (UVEX)

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Autonomous real-time science-driven follow-up of survey transients

Astronomical surveys continue to provide unprecedented insights into the time-variable Universe and will remain the source of groundbreaking discoveries for years to come. However, their data throughput has overwhelmed the ability to manually synthesize alerts for devising and coordinating necessary follow-up with limited resources. The advent of Rubin Observatory, with alert volumes an order of magnitude higher at otherwise sparse cadence, presents an urgent need to overhaul existing human-centered protocols in favor of machine-directed infrastructure for conducting science inference and optimally planning expensive follow-up observations. We present the first implementation of autonomous real-time science-driven follow-up using value iteration to perform sequential experiment design. We demonstrate it for strategizing photometric augmentation of Zwicky Transient Facility Type Ia supernova light-curves given the goal of minimizing SALT2 parameter uncertainties. We find a median improvement of 2-6% for SALT2 parameters and 3-11% for photometric redshift with 2-7 additional data points in g, r and/or i compared to random augmentation. The augmentations are automatically strategized to complete gaps and for resolving phases with high constraining power (e.g. around peaks). We suggest that such a technique can deliver higher impact during the era of Rubin Observatory for precision cosmology at high redshift and can serve as the foundation for the development of general-purpose resource allocation systems.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Confirming new changing--look AGNs discovered through optical variability using a random-forest based light curve classifier

Determining the frequency and duration of changing--look (CL) active galactic nuclei (AGNs) phenomena, where the optical broad emission lines appear or disappear, is crucial to understand the evolution of the accretion flow around supermassive black holes. We present a strategy to select new CL candidates starting from a spectroscopic type 2 AGNs sample and searching for current type 1 photometric variability. We use the publicly available Zwicky Transient Facility (ZTF) alert stream and the Automatic Learning for the Rapid Classification of Events (ALeRCE) light curve classifier to produce a list of CL candidates with a highly automated algorithm, resulting in 60 candidates. Visual inspection reduced the sample to 30. We performed new spectroscopic observations of six candidates of our clean sample, without further refinement, finding the appearance of clear broad Balmer lines in four of them and tentative evidence of type changes in the remaining two, which suggests a promising success rate of $\geq66$ per cent for this CL selection method.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Examining AGN UV/optical Variability Beyond the Simple Damped Random Walk

We present damped harmonic oscillator (DHO) light-curve modeling for a sample of 12,714 spectroscopically confirmed quasars in the Sloan Digital Sky Survey Stripe 82 region. DHO is a second-order continuous-time autoregressive moving-average (CARMA) process, which can be fully described using four independent parameters: a natural oscillation frequency ($ω_{0}$), a damping ratio ($ξ$), a characteristic perturbation timescale ($τ_{\mathrm{perturb}}$), and an amplitude for the perturbing white noise ($σ_{\mathrmε}$). The asymptotic variability amplitude of a DHO process is quantified by $σ_{\mathrm{DHO}}$ -- a function of $ω_{0}$, $ξ$, $τ_{\mathrm{perturb}}$, and $σ_{\mathrmε}$. We find that both $τ_{\mathrm{perturb}}$ and $σ_{\mathrmε}$ follow different dependencies with rest-frame wavelength ($λ_{\mathrm{RF}}$) on either side of 2500 Å, whereas $σ_{\mathrm{DHO}}$ follows a single power-law relation with $λ_{\mathrm{RF}}$. After correcting for wavelength dependence, $σ_{\mathrm{DHO}}$ exhibits anti-correlations with both the Eddington ratio and the black hole mass, while $τ_{\mathrm{perturb}}$ -- with a typical value of days in the rest-frame -- shows an anti-correlation with the bolometric luminosity. Modeling AGN variability as a DHO offers more insight into the workings of accretion disks close to the supermassive black holes (SMBHs) at the center of AGN. The newly discovered short-term variability (characterized by $τ_{\mathrm{perturb}}$ and $σ_{\mathrmε}$) and its correlation with bolometric luminosity pave the way for new algorithms that will derive fundamental properties (e.g., Eddington ratio) of AGN using photometric data alone.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

From Data to Software to Science with the Rubin Observatory LSST

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) dataset will dramatically alter our understanding of the Universe, from the origins of the Solar System to the nature of dark matter and dark energy. Much of this research will depend on the existence of robust, tested, and scalable algorithms, software, and services. Identifying and developing such tools ahead of time has the potential to significantly accelerate the delivery of early science from LSST. Developing these collaboratively, and making them broadly available, can enable more inclusive and equitable collaboration on LSST science. To facilitate such opportunities, a community workshop entitled "From Data to Software to Science with the Rubin Observatory LSST" was organized by the LSST Interdisciplinary Network for Collaboration and Computing (LINCC) and partners, and held at the Flatiron Institute in New York, March 28-30th 2022. The workshop included over 50 in-person attendees invited from over 300 applications. It identified seven key software areas of need: (i) scalable cross-matching and distributed joining of catalogs, (ii) robust photometric redshift determination, (iii) software for determination of selection functions, (iv) frameworks for scalable time-series analyses, (v) services for image access and reprocessing at scale, (vi) object image access (cutouts) and analysis at scale, and (vii) scalable job execution systems. This white paper summarizes the discussions of this workshop. It considers the motivating science use cases, identified cross-cutting algorithms, software, and services, their high-level technical specifications, and the principles of inclusive collaborations needed to develop them. We provide it as a useful roadmap of needs, as well as to spur action and collaboration between groups and individuals looking to develop reusable software for early LSST science.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Impact of the SpaceX Starlink Satellites on the Zwicky Transient Facility Survey Observations

There is a growing concern about an impact of low-Earth-orbit (LEO) satellite constellations on ground-based astronomical observations, in particular, on wide-field surveys in the optical and infrared. The Zwicky Transient Facility (ZTF), thanks to the large field of view of its camera, provides an ideal setup to study the effects of LEO megaconstellations - such as SpaceX's Starlink - on astronomical surveys. Here, we analyze the archival ZTF observations collected between 2019 November and 2021 September and find 5301 satellite streaks that can be attributed to Starlink satellites. We find that the number of affected images is increasing with time as SpaceX deploys more and more satellites. Twilight observations are particularly affected - a fraction of streaked images taken during twilight has increased from less than 0.5% in late 2019 to 18% in 2021 August. We estimate that once the size of the Starlink constellation reaches 10,000, essentially all ZTF images taken during twilight may be affected. However, despite the increase in satellite streaks observed during the analyzed period, the current science operations of ZTF are not yet strongly affected. We also find that redesigning Starlink satellites (by installing visors intended to block sunlight from reaching the satellite antennas to prevent reflection) reduces their brightness by a factor of 4.6 +/- 0.1 with respect to the original design in g, r, and i bands.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

SNIascore: Deep Learning Classification of Low-Resolution Supernova Spectra

We present SNIascore, a deep-learning based method for spectroscopic classification of thermonuclear supernovae (SNe Ia) based on very low-resolution (R $\sim100$) data. The goal of SNIascore is fully automated classification of SNe Ia with a very low false-positive rate (FPR) so that human intervention can be greatly reduced in large-scale SN classification efforts, such as that undertaken by the public Zwicky Transient Facility (ZTF) Bright Transient Survey (BTS). We utilize a recurrent neural network (RNN) architecture with a combination of bidirectional long short-term memory and gated recurrent unit layers. SNIascore achieves a $<0.6\%$ FPR while classifying up to $90\%$ of the low-resolution SN Ia spectra obtained by the BTS. SNIascore simultaneously performs binary classification and predicts the redshifts of secure SNe Ia via regression (with a typical uncertainty of $<0.005$ in the range from $z = 0.01$ to $z = 0.12$). For the magnitude-limited ZTF BTS survey ($\approx70\%$ SNe Ia), deploying SNIascore reduces the amount of spectra in need of human classification or confirmation by $\approx60\%$. Furthermore, SNIascore allows SN Ia classifications to be automatically announced in real-time to the public immediately following a finished observation during the night.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Zwicky Transient Facility phase I sample of hydrogen-rich superluminous supernovae without strong narrow emission lines

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Zwicky Transient Facility and Globular Clusters: The RR Lyrae gri-Band Period-Luminosity-Metallicity and Period-Wesenheit-Metallicity Relations

Based on time-series observations collected from Zwicky Transient Facility (ZTF), we derived period-luminosity-metallicity (PLZ) and period-Wesenheit-metallicity (PWZ) relations for RR Lyrae located in globular clusters. We have applied various selection criteria to exclude RR Lyrae with problematic or spurious light curves. These selection criteria utilized information on the number of data points per light curve, amplitudes, colors, and residuals on the period-luminosity and/or period-Wesenheit relations. Due to blending, a number of RR Lyrae in globular clusters were found to be anomalously bright and have small amplitudes of their ZTF light curves. We used our final sample of ~750 RR Lyrae in 46 globular clusters covering a wide metallicity range (-2.36 < [Fe/H] < -0.54 dex) to derive PLZ and PWZ relations in gri bands. In addition, we have also derived the period-color-metallicity (PCZ) and for the first time, the PQZ relations where the Q-index is extinction-free by construction. We have compared our various relations to empirical and theoretical relations available in literature, and found a good agreement with most studies. Finally, we applied our derived PLZ relation to a dwarf galaxy, Crater II, and found its true distance modulus should be larger than the most recent determination.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Gaia GraL: Gaia DR2 Gravitational Lens Systems. VII. XMM-Newton Observations of Lensed Quasars

We present XMM-Newton X-ray observations of nine confirmed lensed quasars at $1 \lesssim z \lesssim 3$ identified by the Gaia Gravitational Lens program. Eight systems are strongly detected, with 0.3--8.0 keV fluxes $F_{0.3-8.0} \gtrsim 5 \times 10^{-14}\ {\rm erg}\ {\rm cm}^{-2}\ {\rm s}^{-1}$. Modeling the X-ray spectra with an absorbed power law, we derive power law photon indices and 2--10 keV luminosities for the eight detected quasars. In addition to presenting sample properties for larger quasar population studies and for use in planning for future caustic crossing events, we also identify three quasars of interest: a quasar that shows evidence of flux variability from previous ROSAT observations, the most closely-separated individual lensed sources resolved by XMM-Newton, and one of the X-ray brightest quasars known at $z>3$. These sources represent the tip of discovery that will be enabled by SRG/eROSITA.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Inferring kilonova population properties with a hierarchical Bayesian framework I : Non-detection methodology and single-event analyses

We present ${\tt nimbus}$ : a hierarchical Bayesian framework to infer the intrinsic luminosity parameters of kilonovae (KNe) associated with gravitational-wave (GW) events, based purely on non-detections. This framework makes use of GW 3-D distance information and electromagnetic upper limits from multiple surveys for multiple events, and self-consistently accounts for finite sky-coverage and probability of astrophysical origin. The framework is agnostic to the brightness evolution assumed and can account for multiple electromagnetic passbands simultaneously. Our analyses highlight the importance of accounting for model selection effects, especially in the context of non-detections. We show our methodology using a simple, two-parameter linear brightness model, taking the follow-up of GW190425 with the Zwicky Transient Facility (ZTF) as a single-event test case for two different prior choices of model parameters -- (i) uniform/uninformative priors and (ii) astrophysical priors based on surrogate models of Monte Carlo radiative transfer simulations of KNe. We present results under the assumption that the KN is within the searched region to demonstrate functionality and the importance of prior choice. Our results show consistency with ${\tt simsurvey}$ -- an astronomical survey simulation tool used previously in the literature to constrain the population of KNe. While our results based on uniform priors strongly constrain the parameter space, those based on astrophysical priors are largely uninformative, highlighting the need for deeper constraints. Future studies with multiple events having electromagnetic follow-up from multiple surveys should make it possible to constrain the KN population further.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Search for Extra-Tidal RR Lyrae in the Globular Cluster NGC 5024 and NGC 5053

Recently, Kundu et al (2019) reported that the globular cluster NGC 5024 (M53) possesses five extra-tidal RR Lyrae. In fact, four of them were instead known members of a nearby globular cluster NGC 5053. The status of the remaining extra-tidal RR Lyrae is controversial depending on the adopted tidal radius of NGC 5024. We have also searched for additional extra-tidal RR Lyrae within an area of $\sim8$~deg$^2$ covering both globular clusters. This includes other known RR Lyrae within the search area, as well as stars that fall within the expected range of magnitudes and colors for RR Lyrae (and yet outside the cutoff of 2/3 of the tidal radii of each globular clusters for something to be called &#34;extra-tidal&#34;) if they were extra-tidal RR Lyrae candidates for NGC 5024 or NGC 5053. Based on the the proper motion information and their locations on the color-magnitude diagram, none of the known RR Lyrae belong to the extra-tidal RR Lyrae of either globular clusters. In the cases where the stars satisfied the magnitude and color ranges of RR Lyrae, analysis of time series data taken from the Zwicky Transient Facility do not reveal periodicities, suggesting that none of these stars are RR Lyrae. We conclude that there are no extra-tidal RR Lyrae associated with either NGC 5024 or NGC 5053 located within our search area.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Chandra Observations of Candidate Sub-Parsec Binary Supermassive Black Holes

We present analysis of Chandra X-ray observations of seven quasars that were identified as candidate sub-parsec binary supermassive black hole (SMBH) systems in the Catalina Real-Time Transient Survey (CRTS) based on apparent periodicity in their optical light curves. Simulations predict close-separation accreting SMBH binaries will have different X-ray spectra than single accreting SMBHs, including harder or softer X-ray spectra, ripple-like profiles in the Fe K-$α$ line, and distinct peaks in the spectrum due to the separation of the accretion disk into a circumbinary disk and mini-disks around each SMBH. We obtained Chandra observations to test these models and assess whether these quasars could contain binary SMBHs. We instead find that the quasar spectra are all well fit by simple absorbed power law models, with the rest frame 2-10 keV photon indices, $Γ$, and the X-ray-to-optical power slopes, $α_{\rm OX}$, indistinguishable from the larger quasar population. This may indicate that these seven quasars are not truly sub-parsec binary SMBH systems, or it may simply reflect that our sample size was too small to robustly detect any differences. Alternatively, the X-ray spectral changes might only be evident at higher energies than probed by Chandra. Given the available models and current data, no firm conclusions are drawn. These observations will help motivate and direct further work on theoretical models of binary SMBH systems, such as modeling systems with thinner accretion disks and larger binary separations.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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)

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Testing the relativistic Doppler boost hypothesis for the binary candidate quasar PG1302-102 with multi-band Swift data

The bright quasar PG1302-102 has been identified as a candidate supermassive black hole binary from its near-sinusoidal optical variability. While the significance of its optical periodicity has been debated due to the stochastic variability of quasars, its multi-wavelength variability in the ultraviolet (UV) and optical bands is consistent with relativistic Doppler boost caused by the orbital motion in a binary. However, this conclusion was based previously on sparse UV data which were not taken simultaneously with the optical data. Here we report simultaneous follow-up observations of PG1302-102 with the Ultraviolet Optical Telescope on the Neil Gehrels Swift Observatory in six optical + UV bands. The additional nine Swift observations produce light curves roughly consistent with the trend under the Doppler boost hypothesis, which predicts that UV variability should track the optical, but with a ~2.2 times higher amplitude. We perform a statistical analysis to quantitatively test this hypothesis. We find that the data are consistent with the Doppler boost hypothesis when we compare the the amplitudes in optical B-band and UV light curves. However, the ratio of UV to V-band variability is larger than expected and is consistent with the Doppler model, only if either the UV/optical spectral slopes vary, the stochastic variability makes a large contribution in the UV, or the sparse new optical data underestimate the true optical variability. We have evidence for the latter from comparison with the optical light curve from ASAS-SN. Additionally, the simultaneous analysis of all four bands strongly disfavors the Doppler boost model whenever Swift V-band is involved. Additional, simultaneous optical + UV observations tracing out another cycle of the 5.2-year proposed periodicity should lead to a definitive conclusion.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An Unusual Mid-Infrared Flare in a Type 2 AGN: An Obscured Turning-on AGN or Tidal Disruption Event?

We report the discovery of an exceptional MIR flare in a Type 2 AGN, SDSS J165726.81+234528.1, at $z=0.059$. This object brightened by 3 mag in the Wide-field Infrared Survey Explorer (WISE) $W1$ and $W2$ bands between 2015 and 2017 (and is fading since 2018), without significant changes ($\lesssim$ 0.2 mag) in the optical over the same period of time. Based on the WISE light curves and near-IR imaging, the flare is more significant at longer wavelengths, suggesting an origin of hot dust emission. The estimated black hole mass ($\sim 10^{6.5}\,M_\odot$) from different methods places its peak bolometric luminosity around the Eddington limit. The high luminosity of the MIR flare and its multi-year timescale suggest that it most likely originated from reprocessed dust radiation in an extended torus surrounding the AGN, instead of from stellar explosions. The MIR color variability is consistent with known changing-look AGN and tidal disruption events (TDEs), but inconsistent with normal supernovae. We suggest that it is a turning-on Type 2 AGN or TDE, where the optical variability is obscured by the dust torus during the transition. This MIR flare event reveals a population of dramatic nuclear transients that are missed in the optical.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The first high-redshift changing-look quasars

We report on three redshift $z>2$ quasars with dramatic changes in their C IV emission lines, the first sample of changing-look quasars (CLQs) at high redshift. This is also the first time the changing-look behaviour has been seen in a high-ionisation emission line. SDSS J1205+3422, J1638+2827, and J2228+2201 show interesting behaviour in their observed optical light curves, and subsequent spectroscopy shows significant changes in the C IV broad emission line, with both line collapse and emergence being displayed on rest-frame timescales of $\sim$240-1640 days. These are rapid changes, especially when considering virial black hole mass estimates of $M_{\rm BH} > 10^{9} M_{\odot}$ for all three quasars. Continuum and emission line measurements from the three quasars show changes in the continuum-equivalent width plane with the CLQs seen to be on the edge of the full population distribution, and showing indications of an intrinsic Baldwin effect. We put these observations in context with recent state-change models, and note that even in their observed low-state, the C IV CLQs are generally above $\sim$5\% in Eddington luminosity.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Understanding extreme quasar optical variability with CRTS: II. Changing-state quasars

We present the results of a systematic search for quasars in the Catalina Real-time Transient Survey exhibiting both strong photometric and spectroscopic variability over a decadal baseline. We identify 73 sources with specific patterns of optical and mid-IR photometric behavior and a defined spectroscopic change. These &#34;Changing-State&#34; quasars (CSQs) form a higher luminosity sample to complement existing sets of &#34;Changing-Look&#34; AGN and quasars in the literature. The CSQs (by selection) exhibit larger photometric variability than the CLQs. The spectroscopic variability is marginally stronger in the CSQs than CLQs as defined by the change in H$β$/[OIII] ratio. We find 36 sources with declining H$β$ flux, 37 sources with increasing H$β$ flux and discover seven sources with $z > 0.8$, further extending the redshift arm. Our CSQ sample compares to the literature CLQ objects in similar distributions of H$β$ flux ratios and differential Eddington ratios between high (bright) and low (dim) states. Taken as a whole, we find that this population of extreme varying quasars is associated with changes in the Eddington ratio and the timescales imply cooling/heating fronts propagating through the disk.

0 citations0 reviewsNo code linkedOpen access
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).

0 citations0 reviewsNo code linkedOpen access