scispace - formally typeset
Search or ask a question
Journal ArticleDOI

Year 1 of the ZTF high-cadence Galactic plane survey: strategy, goals, and early results on new single-mode hot subdwarf B-star pulsatos

TL;DR: The first results of the high-cadence Galactic plane survey using the Zwicky Transient Facility (ZTF) were presented in this paper, which revealed the Galactic population of short-period variable stars, including short period binaries and stellar pulsators with periods less than a few hours.
Abstract: We present the goals, strategy and first results of the high-cadence Galactic plane survey using the Zwicky Transient Facility (ZTF). The goal of the survey is to unveil the Galactic population of short-period variable stars, including short period binaries and stellar pulsators with periods less than a few hours. Between June 2018 and January 2019, we observed 64 ZTF fields resulting in 2990 deg$^2$ of high stellar density in ZTF-$r$ band along the Galactic Plane. Each field was observed continuously for 1.5 to 6 hrs with a cadence of 40 sec. Most fields have between 200 and 400 observations obtained over 2-3 continuous nights. As part of this survey we extract a total of $\approx$230 million individual objects with at least 80 epochs obtained during the high-cadence Galactic Plane survey reaching an average depth of ZTF-$r$ $\approx$20.5 mag. For four selected fields with 2 million to 10 million individual objects per field we calculate different variability statistics and find that $\approx$1-2% of the objects are astrophysically variable over the observed period. We present a progress report on recent discoveries, including a new class of compact pulsators, the first members of a new class of Roche Lobe filling hot subdwarf binaries as well as new ultracompact double white dwarfs and flaring stars. Finally we present a sample of 12 new single-mode hot subdwarf B-star pulsators with pulsation amplitudes between ZTF-$r$ = 20-76 mmag and pulsation periods between $P$ = 5.8-16 min with a strong cluster of systems with periods $\approx$ 6 min. All of the data have now been released in either ZTF Data Release 3 or data release 4.

Summary (2 min read)

Introduction

  • Key words: surveys – binaries (including multiple): close – stars: oscillations (including pulsations) – white dwarfs.
  • Recently, a new class of short-period pulsating hot stars known as Blue Large-Amplitude Pulsators was discovered by Pietrukowicz et al. (2017).
  • Table 1. Overview of the high-cadence Galactic plane survey.
  • Because more time was available each night most fields were observed for ≈3 h.

2.1 Field selection

  • The main science driver for the survey is to find and study UCBs consisting of fully degenerate or semi-degenerate stars.
  • Hence, the authors selected ZTF fields based on the density of objects residing well below the main sequence, including mostly white dwarfs and hot subdwarf stars.
  • To achieve this, the authors extracted objects with absolute magnitudes which placed them below the main-sequence based on Gaia data release 2 (Gaia Collaboration et al. 2016, 2018).
  • Only objects with declinations >−30◦ were selected.
  • For each object, the ZTF field was calculated and the fields with the largest number of individual objects were selected for their survey.

2.2 Data-processing and light-curve extraction

  • Data-processing and light-curve generation follows the standard procedure for ZTF and occurs at the Infrared Processing and Analysis Center, Caltech.
  • Each reference image is constructed using a minimum of 15 and maximum of 40 good-quality epochal images, yielding depths of ≈2–2.5 mag deeper than the single-epoch images.
  • The grey shaded region corresponds to the underlying Hertzsprung– Russell diagram showing the position of the main sequence and the red giant branch.
  • Further details of the data processing, PSF-fit photometry software, light-curve generation, formats, and overall performance on photometric accuracy are described in Masci et al. (2019).
  • Fig. 3 shows the sky density of the extracted objects overplotted with the fields observed as part of the high-cadence Galactic plane survey.

3 R ESULTS FROM REPRESENTATIVE FIELDS

  • The authors select four representative fields with different stellar densities and calculate statistics for each individual light curve, including median magnitude, reduced χ2, interquartile range (IQR), skewness, inverse von-Neumann, Stetson J, and Stetson K statistics.
  • The two white square areas in FieldID 331 are individual quadrants, which were not processed.
  • Histogram of the IQR extracted at a ZTF-r median magnitude between 15.5 and 16 with the best Gaussian fit and the limit above which an object is called variable.
  • The grey shaded region corresponds to the underlying distribution of all objects in that field.
  • The authors find that out of the selected light curves only 15 per cent for FieldID = 331 and 40 per cent for FieldID = 309 show real astrophysical variability which shows that in the higher density fields the contamination of false positives is significantly larger.

4 PRO G R E S S R E P O RT

  • The primary science driver for the survey is the detection and study of ultracompact binaries and rapid pulsators.
  • Because the data volume is too large to inspect every single light curve, the authors applied different cuts to pre-select light curves in early analysis.
  • As part of their cross-match with the hot subdwarf catalogue, the authors discovered ZTF J2130, the most compact sdB+WD binary and the first member of systems where the sdB fills its.
  • The typical diurnal mean sampling rate of ZTF is too short to resolve the pulsations for many pulsating white dwarfs (ZZ Ceti pulsation modes range from 100–1500 sec Mukadam et al. 2013; Bognár et al. 2020).
  • Using the first data release from the 14 August fields, more than 1500 stellar flare events have been detected (Klein et al., in preparation).

5 EARLY R ESULTS – N EW SINGLE-MODE HOT SUB DWARF PULSATORS

  • Subdwarf B stars (sdBs) are hot stars of spectral type B with luminosities below the main sequence.
  • The authors discovered 12 new sdB pulsators that show only a single pulsation mode as part of this survey.
  • The authors have presented an overview and first results of their dedicated high-cadence Galactic plane survey for short-period variable objects, carried out as part of ZTF during its first year of operation.
  • The authors calculate light curve statistics for each field and use the IQR to estimate the number of variable objects for each field.
  • This research was funded by the Gordon and Betty Moore Foundation through Grant GBMF5076.

Did you find this useful? Give us your feedback

Content maybe subject to copyright    Report

MNRAS 505, 1254–1267 (2021) https://doi.org/10.1093/mnras/stab1344
Advance Access publication 2021 May 12
Year 1 of the ZTF high-cadence Galactic plane survey: strategy, goals, and
early results on new single-mode hot subdwarf B-star pulsatos
Thomas Kupfer,
1,2
Thomas A. Prince ,
3
Jan van Roestel,
3
Eric C. Bellm,
4
Lars Bildsten,
5,1
Michael W. Coughlin ,
6
Andrew J. Drake,
3
Matthew J. Graham ,
3
Courtney Klein ,
7
Shrinivas R. Kulkarni,
3
Frank J. Masci,
8
Richard Walters,
9
Igor Andreoni,
3
Rahul Biswas ,
10
Corey Bradshaw,
2
Dmitry A. Duev ,
3
Richard Dekany,
9
Joseph A. Guidry,
11
J. J. Hermes,
12
Russ R. Laher
8
and Reed Riddle
9
1
Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106, USA
2
Department of Physics and Astronomy, Texas Tech University, PO Box 41051, Lubbock, TX 79409, USA
3
Division of Physics, Mathematics and Astronomy, California Institute of Technology, Pasadena, CA 91125, USA
4
DIRAC Institute, Department of Astronomy, University of Washington, 3910 15th Avenue NE, Seattle, WA 98195, USA
5
Department of Physics, University of California, Santa Barbara, CA 93106, USA
6
School of Physics and Astronomy, University of Minnesota, Minneapolis, MN 55455, USA
7
Department of Physics and Astronomy, University of California, Irvine, CA 92697, USA
8
IPAC, California Institute of Technology, 1200 E. California Blvd, Pasadena, CA 91125, USA
9
Caltech Optical Observatories, California Institute of Technology, Pasadena, CA 91125, USA
10
Department of Physics, The Oskar Klein Center, S tockholm University, AlbaNova, SE-10691 Stockholm, Sweden
11
Department of Astronomy, The University of Texas at Austin, Austin, TX 78712, USA
12
Department of Astronomy, Boston University, 725 Commonwealth Ave., Boston, MA 02215, USA
Accepted 2021 May 5. Received 2021 May 5; in original form 2021 January 14
ABSTRACT
We present the goals, strategy, and first results of the high-cadence Galactic plane survey using the Zwicky Transient Facility
(ZTF). The goal of the survey is to unveil the Galactic population of short-period variable stars, including short-period binaries,
and stellar pulsators with periods less than a few hours. Between 2018 June and 2019 January, we observed 64 ZTF fields resulting
in 2990 deg
2
of high stellar density in the ZTF-r band along the Galactic plane. Each field was observed continuously for 1.5 to
6 h with a cadence of 40 sec. Most fields have between 200 and 400 observations obtained over 2–3 continuous nights. As part
of this survey, we extract a total of 230 million individual objects with at least 80 epochs obtained during the high-cadence
Galactic plane survey reaching an average depth of ZTF–r 20.5 mag. For four selected fields with 2–10 million individual
objects per field, we calculate different variability statistics and find that 1–2 per cent of the objects are astrophysically variable
over the observed period. We present a progress report on recent discoveries, including a new class of compact pulsators, the first
members of a new class of Roche lobe filling hot subdwarf binaries as well as new ultracompact double white dwarfs and flaring
stars. Finally, we present a sample of 12 new single-mode hot subdwarf B-star pulsators with pulsation amplitudes between
ZTF–r = 20–76 mmag and pulsation periods between P = 5.8–16 min with a strong cluster of systems with periods 6 min. All
of the data have now been released in either ZTF Data Release 3 or Data Release 4.
Key words: surveys binaries (including multiple): close stars: oscillations (including pulsations) white dwarfs.
1 INTRODUCTION
Large-scale optical time-domain surveys have opened a new window
to study the variable sky providing hundreds to thousands of epochs
across the whole sky. Starting with the Sloan Digital Sky Survey
(SDSS; York et al. 2000), a new generation of wide-field optical
surveys has exploited new affordable CCD detectors to open the
frontier of data-intensive astronomy. This allows the study of stellar
variability on different time-scales across the full magnitude range.
In particular, surveys covering also low Galactic latitudes are well
suited to study the Galactic distribution of photometric variable stars.
E-mail: tkupfer@ttu.edu
Ground-based surveys include the Optical Gravitational Lensing
Experiment (OGLE; e.g. Soszy
´
nski et al. 2015),PTF(Lawetal.
2009), the Vista Variables in the Via Lactea (; Saito et al. 2012),
ASAS-SN, (Shappee et al. 2014; Jayasinghe et al. 2018), and most
recently ATLAS (Heinze et al. 2018; Tonry et al. 2018). In particular,
the fastest stellar variabilities on time-scales of minutes to hours are
of great interest for a large number of scientific questions. This
includes ultracompact binaries (UCBs), compact pulsators as well as
fast flaring stars.
UCBs are a class of binary stars with orbital periods less than
about 60 min, consisting of a neutron star (NS)/white dwarf (WD)
primary and a Helium-star (He-star)/WD/NS secondary. These UCBs
are sources of low-frequency gravitational wave signals as probed
C
2021 The Author(s)
Published by Oxford University Press on behalf of Royal Astronomical Society
Downloaded from https://academic.oup.com/mnras/article/505/1/1254/6274690 by California Institute of Technology user on 21 July 2021

ZTF high-cadence Galactic plane survey 1255
by the Laser Interferometer Space Antenna (LISA) and are crucial to
our understanding of compact binary evolution and offer pathways
towards Type Ia and other thermonuclear supernovae. Systems with
orbital periods <20 min will be the strongest Galactic LISA sources
and will be detected by LISA within weeks after its operation begins
(Nelemans, Yungelson & Portegies Zwart 2004;Str
¨
oer & Vecchio
2006;Nissankeetal.2012; Littenberg et al. 2013; Korol et al. 2017;
Kremer et al. 2017; Kupfer et al. 2018; Lamberts et al. 2019; Burdge
et al. 2019, 2020a) and as such are ideal multimessenger sources
(Shah, van der Sluys & Nelemans 2012; Shah & Nelemans 2014;
Baker et al. 2019; Littenberg et al. 2019; Kupfer et al. 2019a)
Short time-scale photometric variations can also originate from
astrophysical changes within the internal structure or atmosphere
of the star. Such sources are flaring stars, pulsating stars, or white
dwarfs. Prominent examples of rapidly pulsating stars are δ Scuti
stars (Breger 2000), SX Phoenicis variables (Nemec & Mateo 1990),
or Ap or Am stars (Renson, Gerbaldi & Catalano 1991), and pulsating
white dwarfs (ZZ Ceti variables; Fontaine & Brassard 2008) as well
as pulsating extremely low-mass white dwarfs (Hermes et al. 2013).
The stars exhibit pulsation amplitudes of less than 1 percent up to
several tens of percent on time-scales of few to tens of minutes.
Recently, a new class of short-period pulsating hot stars known
as Blue Large-Amplitude Pulsators (BLAPs) was discovered by
Pietrukowicz et al. (2017). Their pulsation periods are typically
between 20–40 min. Romero et al. (2018) and Byrne & Jeffery (2018)
proposed that the BLAPs are hot pre-helium white dwarfs that are
cooling and contracting, with masses in the range 0.3–0.35 M
.A
recent study by Meng et al. (2020) suggests that BLAPs could be
the surviving companions of type Ia supernovae. In this scenario
mass transfer from the BLAP progenitor leads to the detonation of a
white dwarf companion, leaving behind a BLAP as single star. Their
pulsation properties can best be explained by fundamental radial
mode pulsators. However, the known sample of BLAPs is small and
in-homogeneous, discovered only in the OGLE survey.
A number of fast cadence ground-based surveys have been
executed to study the variable sky down to a few minutes period. The
main goal of these high-cadence surveys is the discovery of rapid
brightness variations seen in UCBs, compact pulsators as well as fast
flaring stars. The first survey at low Galactic latitudes targeting short-
period systems was the Rapid Temporal Survey (RATS; Ramsay
& Hakala 2005; Barclay et al. 2011) covering a total of 46 deg
2
.
Another more recent survey is the OmegaWhite (OW) survey, which
covers a total of 400 deg
2
at low Galactic latitudes (|b| < 10
)aswell
as in the Galactic Bulge using high-cadence optical observations.
Two neighbouring 1 deg
2
fields are alternatingly observed in 39-s
exposures over an observing duration of 2 h, with an observational
median cadence of 2.7 min per field (Macfarlane et al. 2015,
2017a,b; Toma et al. 2016; Kupfer et al. 2017)
As part of the Zwicky Transient Facility (ZTF), the Palomar 48-
inch (P48) telescope images the sky every clear night conducting
several surveys including a Northern Sky Survey with a 3-d cadence,
as well as smaller surveys such as a 1-d Galactic plane survey and
simultaneous observations of the Northern TESS sectors (Graham
et al. 2019; van Roestel et al. 2019; Bellm et al. 2019a,b). As part of
the partnership share of ZTF, we conducted a dedicated high-cadence
Galactic plane survey with a cadence of 40 sec at low Galactic
latitudes aiming to find UCBs and compact pulsators. During that
dedicated survey, we either observed one field or alternated between
two adjacent fields continuously for 1.5–3 h on two to three
consecutive nights in the ZTF-r band. Here, we present an overview
of the ZTF high-cadence Galactic plane survey executed in ZTF year
1. We present the observing strategy as well as some survey statistics.
Tab le 1 . Overview of the high-cadence Galactic plane survey.
Period No. of fields Sky coverage Filter
(deg
2
)
06-15-2018–07-31-2018 16 750 ZTF-r
08-03-2018–08-18-2018 14 650 ZTF-r
11-15-2018–01-15-2019 34 1590 ZTF-r
We show a progress report and finalize with some early results from
the survey. In Section 2, we discuss the observing strategy as well
as the field selection and the data processing. In Section 3, we
present results for four representative fields and in Section 4 we give
a progress report of already published results. In Section 5, we show
some new early results from the high-cadence Galactic plane survey
and summarize and conclude in Section 6.
2 DESIGN OF THE HIGH-CADENCE
GALACTIC PLANE SURVEY
ZTF uses a 47 deg
2
camera consisting of 16 individual CCDs each 6k
× 6k covering the full focal plane of the P48 telescope. The ZTF high-
cadence Galactic plane survey covered a total of 2990 deg
2
split
in 64 individual ZTF fields observed over three observing blocks;
mid-June–July 2018, two weeks in 2018 August and November 15
to 2019 January 15. All observations of the high-cadence Galactic
plane survey were obtained in ZTF-r band. Each field has 200–
400 epochs. The images were processed using ZTF data-processing
pipeline described in full detail in Masci et al. (2019).
During June and July, we observed 16 fields covering 750 deg
2
mostly at low Galactic longitudes. In June, we observed every field
continuously for 1 h15 min and in July for 1 h 25 min. All fields
were observed over 2 or 3 consecutive nights. In August, 14 fields
covering 650 deg
2
were observed over two weeks (see Table 1).
We alternated between two adjacent fields continuously for 2 h
40min each night. The same fields were repeated the following night.
The observations in June/July and August were done under stable
conditions with an average seeing of 2 arcsec. We lost only a total of
five nights due to weather during June/July and August observations.
Between 2018 November 15 and 2019 January 15 high-cadence
Galactic plane observations were scheduled for every night. We
observed an additional 33 fields covering 1550 deg
2
with stable
weather conditions (see Table 1). The overall strategy varied during
those two months due to unstable weather. Most fields were observed
alternating with adjacent fields but in particular low declination
fields were observed continuously. Because more time was available
each night most fields were observed for 3 h. However, due to
the unstable weather conditions about half of the fields were only
observed in a single night and not repeated in subsequent nights. All
other fields were observed over 2 or 3 nights. Although the seeing
varies strongly between 1.7 and 4 arcsec, each field has a limiting
magnitude of >19.5 mag. A detailed overview of the fields observed
in stable weather conditions is given in Tables A1A3.
2.1 Field selection
The main science driver for the survey is to find and study UCBs
consisting of fully degenerate or semi-degenerate stars. Hence, we
selected ZTF fields based on the density of objects residing well
below the main sequence, including mostly white dwarfs and hot
subdwarf stars.
MNRAS 505, 1254–1267 (2021)
Downloaded from https://academic.oup.com/mnras/article/505/1/1254/6274690 by California Institute of Technology user on 21 July 2021

1256 T. Kupfer et al.
Figure 1. Sky density of candidate white dwarfs and hot subdwarfs selected from Gaia DR2 with the selected fields of the ZTF high-cadence Galactic plane
survey observed in ZTF year-1. The squares show individual ZTF fields which have been observed in high-cadence Galactic plane observations. The white lines
correspond to the Galactic equator and |b|=15
o
.
To achieve this, we extracted objects with absolute magnitudes
which placed them below the main-sequence based on Gaia data
release 2 (Gaia Collaboration et al. 2016, 2018). Only objects with de-
clinations >30
were selected. As the goal was not to extract a clean
sample we used a very relaxed tolerance for the parallax precision
( /σ
> 3) and did not include any additional quality cuts, resulting
in 350 000 individual objects (Fig. 2). For each object, the ZTF field
was calculated and the fields with the largest number of individual
objects were selected for our survey. Selected fields with the highest
density have 1000 selected objects per deg
2
. The final field selection
was mainly driven by stellar density but also included other aspects
like visibility and weather. Fig. 1 shows the sky density of the selected
sources overplotted with the fields observed as part of the high-
cadence Galactic plane survey. As the stellar density is highest at low
Galactic latitudes most fields are located at Galactic latitudes 15
o
.
2.2 Data-processing and light-curve extraction
Data-processing and light-curve generation follows the standard
procedure for ZTF and occurs at the Infrared Processing and Analysis
Center, Caltech. The raw camera image data are first instrumentally
calibrated and astrometric solutions are derived using the Gaia DR1
catalogue. Sources are detected and fluxes measured using both aper-
ture (Bertin & Arnouts 1996) and PSF-fit photometry (Stetson 1987).
Sources are photometrically calibrated using the Pan-STARRS1 DR1
catalogue. The epochal image data are then co-added within their
respective survey fields and camera readout channels to construct
reference images. Each reference image is constructed using a
minimum of 15 and maximum of 40 good-quality epochal images,
yielding depths of 2–2.5 mag deeper than the single-epoch images.
The reference image co-adds are archived for use in other down-
stream processing: image differencing and light curve construction.
To support light curve generation, sources are first detected and
Figure 2. Hertzsprung–Russell diagram of the objects used for the field
selection. The grey shaded region corresponds to the underlying Hertzsprung–
Russell diagram showing the position of the main sequence and the red giant
branch. The colour–coded region corresponds to the objects which were
selected for the field selection. The clump around M
G
5 corresponds to the
hot subdwarfs whereas the region below corresponds to the white dwarfs.
extracted from each reference image using PSF-fit photometry
by running the DAOPhot utility (Stetson 1987). These sources
provide the seeds to facilitate positional cross-matching across all
the single-epoch-based PSF extraction catalogues going back to
the beginning of the survey. PSF-fitting on the single-epoch images
is also performed using DAOPhot. These single-epoch images are
direct single exposure images, not difference images. A fixed source-
match radius of 1.5 arcsec is used for the positional matching. Only
the PSF-fit-derived positions and photometry, along with selected
PSF-fit metrics are retained during the source-matching process.
Aperture photometry is not propagated to the light-curve metadata.
MNRAS 505, 1254–1267 (2021)
Downloaded from https://academic.oup.com/mnras/article/505/1/1254/6274690 by California Institute of Technology user on 21 July 2021

ZTF high-cadence Galactic plane survey 1257
Figure 3. Sky density of individual objects with at least 80 epochs from the ZTF high-cadence Galactic plane survey observed in ZTF year-1. The squares show
individual ZTF fields, which have been observed in high-cadence Galactic plane observations. The white lines correspond to the Galactic equator and |b|=15
o
.
Further details of the data processing, PSF-fit photometry soft-
ware, light-curve generation, formats, and overall performance on
photometric accuracy are described in Masci et al. (2019). We
extracted the light curves from each object which has at least 80
epochs obtained as part of the high-cadence Galactic plane survey
from the ZTF light curve data base ‘Kowalski’ and extracted a
total of 230 million light curves across all observed fields. Fig. 3
shows the sky density of the extracted objects overplotted with
the fields observed as part of the high-cadence Galactic plane
survey.
3 RESULTS FROM REPRESENTATIVE FIELDS
We select four representative fields with different stellar densities and
calculate statistics for each individual light curve, including median
magnitude, reduced χ
2
, interquartile range (IQR), skewness, inverse
von-Neumann, Stetson J, and Stetson K statistics. We refer to table 2
in Coughlin et al. (2020) for a detailed description of the individual
statistics. These statistics are used to evaluate light curve variability
and identify variable objects in the high-cadence Galactic plane data.
The selected fields are FieldID = 331 (l=8.
59, b=8.
69),
FieldID = 538 (l = 44.
64, b = 3.
09), FieldID = 769
((l = 91.
87, b = 2.
38)l = 91.87 deg, b =−2.38 deg) and FieldID =
309 (l = 230.
68, b = 2.
98). Fig. 4 shows the stellar density of the
four selected fields. The different colouring correspond to different
densities. Some structure of lower and higher density regions can
be seen in the fields. The white circles are masked regions due to
saturated stars and the horizontal and vertical gaps are chip gaps
between the 16 CCDs. The two white square areas in FieldID
331 are individual quadrants, which were not processed. FieldID
= 331 was observed in three consecutive nights for 1.5 h each
night. FieldID = 538 was observed for two consecutive nights for
1.5 and 2.5 h. FieldID = 769 was observed for two consecutive
nights for 3 h per night and FieldID = 309 was observed for two
consecutive nights for 1.7 and 3 h. See Tables A1A3 for more
details. The four selected fields represent a large range of stellar
densities from 10 million individual light curves for FieldID =
331, 5.8 million for FieldID = 538, 3.8 million for FieldID =
769, and 2.3 million for FieldID = 309.
Different statistics can be used to evaluate whether an object
is considered variable or not. We decided to use the IQR which
corresponds to the difference between the upper and lower quartile
values. The main advantage of IQR is that it is a robust statistic not
affected by extreme outliers due to individual bad photometry data
points and therefore a better estimate of the general spread around the
median magnitude. When there are outliers in a sample, the median
and IQR are best used to summarize a typical value and the variability
in the sample, respectively. Large variability in the light curve leads
to a larger IQR.
To estimate whether an object is variable in the high-cadence
Galactic plane data, we used the following approach. We calculate the
median ZTF-r band magnitude for each object. The full magnitude
range of all objects in a ZTF field is binned into several bins and each
object is added to a bin based on its median magnitude. We used bin
sizes of 0.5 mag between 12 and 16.5 mag where the IQR is constant
at 0.02 mag for non-variable sources. Above 16.5 mag, we used
bin sizes of 0.2 mag because the IQR for non-variable sources is
increasing with increasing magnitudes to 0.2 mag at the faint end
at ZTF-r = 20.5 mag because the underlying noise of the sources
is increasing for fainter sources. Each bin contains at least a few
thousand objects.
For each magnitude bin, we then calculate a histogram of IQR
values. If there are only constant sources with a similar noise pattern,
a normal distribution in an IQR histogram is expected. An excess
of variable light curves will result in a deviation from a normal
distribution towards more objects with larger IQR values. To estimate
the excess and number of variable objects, we fit a Gaussian to each
IQR histogram and define an object as variable if, in a given histogram
bin, less than 5 per cent of the sources in that magnitude bin are
MNRAS 505, 1254–1267 (2021)
Downloaded from https://academic.oup.com/mnras/article/505/1/1254/6274690 by California Institute of Technology user on 21 July 2021

1258 T. Kupfer et al.
Tab le 2 . Photometric properties of the high-gravity-BLAPs.
Object RA (J2000) Dec. (J2000) P A
b
ZTF-r
c
g
a
g r
a
FieldID
(
h
:
min
:
sec
)(
:
:

) (s) (mmag) (mas) (mag) (mag)
ZTF-sdBV1 18:26:36.09 +10:00:22.3 347.287 ± 0.02 19.9 0.5640 ± 0.0372 14.90 0.41 538
ZTF-sdBV2 19:13:54.77 11:05:19.4 360.576 ± 0.02 19.6 0.4568 ± 0.0888 16.24 0.33 385
ZTF-sdBV3 21:03:08.50 +44:43 14.9 365.822 ± 0.01 53.3 0.3550 ± 0.0769 16.07 0.31 768
ZTF-sdBV4 08:06:07.05 20:08 19.1 370.481 ± 0.02 26.0 0.4021 ± 0.0531 16.37 0.32 311
ZTF-sdBV5 21:49:45.56 +45:58 36.5 375.690 ± 0.02 36.1 0.3112 ± 0.0683 16.55 0.35 770
ZTF-sdBV6 07:40:00.71 14:42 02.5 375.782 ± 0.02 21.1 0.5979 ± 0.0394 15.10 0.36 310
ZTF-sdBV7 19:24:08.62 12:58 30.9 377.468 ± 0.02 19.3 0.7399 ± 0.0355 14.64 0.41 385
ZTF-sdBV8 06:54:48.55 25:22 08.3 435.459 ± 0.02 31.4 0.5458 ± 0.0318 13.93 0.44 259
ZTF-sdBV9 18:52:49.07 18:46 08.4 524.829 ± 0.02 73.4 0.1427 ± 0.1538 17.12 0.38 334
ZTF-sdBV10 19:49:59.34 +08:31 06.1 646.306 ± 0.02 32.8 0.3884 ± 0.0642 16.23 0.38 541
ZTF-sdBV11 07:18:43.69 02:29 31.2 688.485 ± 0.02 30.2 0.1571 ± 0.1329 17.52 0.33 411
ZTF-sdBV12 17:07:41.34
15:22 42.9 995.479 ± 0.02 37.3 0.4376 ± 0.0346 14.11 0.35 330
a
From the Pan-STARRS release 1 (PS1) survey (Chambers et al. 2016) reddening corrected using Green et al. (2019).
b
Amplitude in ZTF-r from the ZTF light curves.
c
Parallaxes from Gaia EDR3 (Gaia Collaboration et al. 2016, 2020).
Figure 4. Sky density of ZTF objects with at least 80 epochs obtained during the high-cadence Galactic plane survey for four selected fields. The white circles
are masked regions due to saturated stars and the horizontal and vertical gaps are chip gaps between the 16 CCDs. The two white square areas in FieldID 331
are individual quadrants, which were not processed.
considered constant based on the Gaussian fit. Fig. 5 shows two
examples of IQR histograms with Gaussian fits of FieldID = 309
and 331 for the magnitude range 15.5–16 mag. The orange curve
corresponds to the Gaussian fit of the underlying histogram and the
red-dashed line corresponds to the IQR limit above which we call an
object a variable.
Using this method, we find 15 per cent variable sources for
FieldID = 331, 9 per cent variable sources for FieldID = 769, and
FieldID = 538 and 5 per cent variable sources for FieldID = 309 (see
Fig. 6). The number of variable sources per field decreases with the
stellar density of the field and it is likely that blending could produce
false positives where the object appears variable due to a close-by
MNRAS 505, 1254–1267 (2021)
Downloaded from https://academic.oup.com/mnras/article/505/1/1254/6274690 by California Institute of Technology user on 21 July 2021

Citations
More filters
Journal ArticleDOI
TL;DR: In this paper , a review of the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals.
Abstract: The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe.

56 citations

Journal ArticleDOI
TL;DR: In this article , the authors employ a model developed by Maoz, Hallakoun and Badenes (2018) for the statistical analysis of the local DWD population using two complementary large, multi-epoch, spectroscopic samples: the Sloan Digital Sky Survey (SDSS) and the Supernova Ia Progenitor surveY (SPY).
Abstract: Realistic models of the Galactic double white dwarf (DWD) population are crucial for testing and quantitatively defining the science objectives of the Laser Interferometer Space Antenna (LISA), a future European Space Agency's gravitational-wave observatory. In addition to numerous individually detectable DWDs, LISA will also detect an unresolved confusion foreground produced by the underlying Galactic population, which will affect the detectability of all LISA sources at frequencies below a few mHz. So far, the modelling of the DWD population for LISA has been based on binary population synthesis (BPS) techniques. The aim of this study is to construct an observationally driven population. To achieve this, we employ a model developed by Maoz, Hallakoun and Badenes (2018) for the statistical analysis of the local DWD population using two complementary large, multi-epoch, spectroscopic samples: the Sloan Digital Sky Survey (SDSS), and the Supernova Ia Progenitor surveY (SPY). We calculate the number of LISA-detectable DWDs and the Galactic confusion foreground, based on their assumptions and results. We find that the observationally driven estimates yield 1) 2 - 5 times more individually detectable DWDs than various BPS forecasts, and 2) a significantly different shape of the DWD confusion foreground. Both results have important implications for the LISA mission. A comparison between several variations to our underlying assumptions shows that our observationally driven model is robust, and that the uncertainty on the total number of LISA-detectable DWDs is in the order of 20 per cent.

15 citations

Journal ArticleDOI
TL;DR: In this paper , the discovery of a 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 was presented.
Abstract: 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.

14 citations

Journal ArticleDOI
TL;DR: In this paper, the authors reported the discovery of five new eclipsing AM CVn systems with orbital periods of 61.5, 55.3, 37.4, and 35.4 minutes by searching for deep eclipses in the Zwicky Transient Facility (ZTF) lightcurves of white dwarfs.
Abstract: AM CVn systems are ultra-compact, helium-rich, accreting binaries with degenerate or semi-degenerate donors. We report the discovery of five new eclipsing AM CVn systems with orbital periods of 61.5, 55.5, 53.3, 37.4, and 35.4 minutes. These systems were discovered by searching for deep eclipses in the Zwicky Transient Facility (ZTF) lightcurves of white dwarfs selected using Gaia parallaxes. We obtained phase-resolved spectroscopy to confirm that all systems are AM CVn binaries, and we obtained high-speed photometry to confirm the eclipse and characterize the systems. The spectra of two long-period systems (61.5 and 53.3 minutes) show many emission and absorption lines, indicating the presence of N, O, Na, Mg, Si, and Ca, and also the K and Zn, elements which have never been detected in AM CVn systems before. By modelling the high-speed photometry, we measured the mass and radius of the donor star, potentially constraining the evolutionary channel that formed these AM CVn systems. We determined that the average mass of the accreting white dwarf is $\approx0.8$$\mathrm{M_{\odot}}$, and that the white dwarfs in long-period systems are hotter than predicted by recently updated theoretical models. The donors have a high entropy and are a factor of $\approx$ 2 more massive compared to zero-entropy donors at the same orbital period. The large donor radius is most consistent with He-star progenitors, although the observed spectral features seem to contradict this. The discovery of 5 new eclipsing AM~CVn systems is consistent with the known observed AM CVn space density and estimated ZTF recovery efficiency. Based on this estimate, we expect to find another 1--4 eclipsing AM CVn systems as ZTF continues to obtain data. This will further increase our understanding of the population, but will require high precision data to better characterize these 5 systems and any new discoveries.

14 citations

Journal ArticleDOI
TL;DR: In this paper , the authors correct the Gaia G-band apparent magnitude and GBP − GRP colours of 89.6 million sources in the Galactic plane with good quality photometry combined with supplementary all-sky data totalling 162.3 million sources.
Abstract: Blue Large Amplitude Pulsators (BLAPs) are hot, subluminous stars undergoing rapid variability with periods of under 60 mins. They have been linked with the early stages of pre-white dwarfs and hot subdwarfs. They are a rare class of variable star due to their evolutionary history within interacting binary systems and the short timescales relative to their lifetime in which they are pulsationally unstable. All currently known BLAPs are relatively faint (15 − 19 mag) and are located in the Galactic plane. These stars have intrinsically blue colours but the large interstellar extinction in the Galactic plane prevents them from swift identification using colour-based selection criteria. In this paper, we correct the Gaia G-band apparent magnitude and GBP − GRP colours of 89.6 million sources brighter than 19 mag in the Galactic plane with good quality photometry combined with supplementary all-sky data totalling 162.3 million sources. Selecting sources with colours consistent with the known population of BLAPs and performing a cross-match with the Zwicky Transient Facility (ZTF) DR3, we identify 98 short period candidate variables. Manual inspection of the period-folded light curves reveals 22 candidate BLAPs. Of these targets, 6 are consistent with the observed periods and light curves of the known BLAPs, 10 are within the theoretical period range of BLAPs and 6 are candidate high-gravity BLAPs. We present follow-up spectra of 21 of these candidate sources and propose to classify 1 of them as a BLAP, and tentatively assign an additional 8 of them as BLAPs for future population studies.

9 citations

References
More filters
Journal ArticleDOI
TL;DR: A previously unknown class of variable stars has been found in OGLE survey data, characterized by periodic brightness variations on ~30min timescales, amplitudes of ~0.3 mag and temperatures of ~30,000 K as mentioned in this paper.
Abstract: A previously unidentified class of variable stars has been found in OGLE survey data, characterized by periodic brightness variations on ~30-min timescales, amplitudes of ~0.3 mag and temperatures of ~30,000 K. They are potentially evolved low-mass stars.

52 citations

Journal ArticleDOI
TL;DR: In this article, the authors reported the discovery of the first short-period binary in which a hot subdwarf star (sdOB) filled its Roche lobe and started mass transfer to its companion.
Abstract: We report the discovery of the first short-period binary in which a hot subdwarf star (sdOB) filled 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 = 39.3401(1) minutes, 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_(eff) = 42,400 ± 300 K and a surface gravity of log(g) = 5.77 ± 0.05. A high signal-to-noise ratio 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_(sdOB) = 0.337 ± 0.015 M⊙ and a white dwarf accretor with a mass M_(WD) = 0.545 ± 0.020 M⊙. Theoretical binary modeling indicates the hot subdwarf formed during a common envelope phase when a 2.5–2.8 M⊙ star lost its envelope when crossing the Hertzsprung gap. To match its current P_(orb), T_(eff), log(g), and masses, we estimate a post–common envelope period of P_(orb) ≈ 150 minutes and find that 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 ≈0.1 M⊙, merge with its carbon/oxygen white dwarf companion after ≈17 Myr, and presumably explode as a thermonuclear supernova or form an R CrB star.

45 citations

Journal ArticleDOI
TL;DR: In this article, the evolutionary rate of cooling of the pulsating hydrogen atmosphere (DA) white dwarf ZZ Ceti (Ross-548), as reflected by the drift rate of the 213.13260694s period, was determined.
Abstract: We have finally measured the evolutionary rate of cooling of the pulsating hydrogen atmosphere (DA) white dwarf ZZ Ceti (Ross 548), as reflected by the drift rate of the 213.13260694 s period. Using 41 yr of time-series photometry from 1970 November to 2012 January, we determine the rate of change of this period with time to be dP/dt = (5.2 ± 1.4) × 10–15 s s–1 employing the O – C method and (5.45 ± 0.79) × 10–15 s s–1 using a direct nonlinear least squares fit to the entire lightcurve. We adopt the dP/dt obtained from the nonlinear least squares program as our final determination, but augment the corresponding uncertainty to a more realistic value, ultimately arriving at the measurement of dP/dt = (5.5 ± 1.0) × 10–15 s s–1. After correcting for proper motion, the evolutionary rate of cooling of ZZ Ceti is computed to be (3.3 ± 1.1) × 10–15 s s–1. This value is consistent within uncertainties with the measurement of (4.19 ± 0.73) × 10–15 s s–1 for another similar pulsating DA white dwarf, G 117-B15A. Measuring the cooling rate of ZZ Ceti helps us refine our stellar structure and evolutionary models, as cooling depends mainly on the core composition and stellar mass. Calibrating white dwarf cooling curves with this measurement will reduce the theoretical uncertainties involved in white dwarf cosmochronometry. Should the 213.13 s period be trapped in the hydrogen envelope, then our determination of its drift rate compared to the expected evolutionary rate suggests an additional source of stellar cooling. Attributing the excess cooling to the emission of axions imposes a constraint on the mass of the hypothetical axion particle.

41 citations

Journal ArticleDOI
TL;DR: In this article, a first-stage study of using knowledge from electromagnetic (EM) observations in the gravitational wave (GW) data analysis of Galactic binaries that are predicted to be observed by the new Laser Interferometer Space Antenna in the lowfrequency range, 10 4 Hz < f < 1Hz.
Abstract: We present a first-stage study of the e ect of using knowledge from electromagnetic (EM) observations in the gravitational wave (GW) data analysis of Galactic binaries that are predicted to be observed by the new Laser Interferometer Space Antenna in the lowfrequency range, 10 4 Hz < f < 1Hz. In particular, we examine the extent to which the accuracy of GW parameter estimation improves if we use available information from EM data. We do this by investigating whether correlations exist between the GW parameters that describe these binaries and whether some of these parameters are also available from EM observations. We used verification binaries, which are known as the guaranteed sources for eLISA and will test the functioning of the instrument. We find that of the seven parameters that characterise such a binary, only a few are correlated. The most useful result is the strong correlation between amplitude and inclination, which can be used to constrain the parameter uncertainty in amplitude by making use of the constraint of inclination from EM measurements. The improvement can be up to a factor of 6:5, but depends on the signal-to-noise ratio of the source data. Moreover, we find that this strong correlation depends on the inclination. For mildly face-on binaries ( . 45 ), EM data on inclination can improve the estimate of the GW amplitude by a significant factor. However, for edge-on binaries ( 90 ), the inclination can be determined accurately from GW data alone, thus GW data can be used to select systems that will likely be eclipsing binaries for EM follow-up.

40 citations

Journal ArticleDOI
TL;DR: In this paper, the authors reported the discovery of pulsations in the light curve of Balloon 090100001, a bright sdB star from the balloon-borne survey of Bixler et al.
Abstract: We report the discovery of pulsations in the light curve of Balloon 090100001, a bright sdB star () from the balloon-borne survey of Bixler et al. (1991). Our time series observations made with the IAC80 telescope reveal a powerful short period oscillation with a period of 356 s and an amplitude (~59 mmag) only matched in all of the sdB pulsators by the main pulsation of PG 1605+072. The first harmonic (at 178 s) and another independent period (at 264 s) are also present in the amplitude spectrum of this new EC 14026 star. New spectroscopic data from the William Herschel Telescope provide the parameters K and dex. These place Balloon 090100001 in the low temperature and low gravity corner of the sdBV instability region, making it an interesting target for study of possible frequency changes due to stellar evolution.

38 citations

Frequently Asked Questions (1)
Q1. What are the contributions in "Year 1 of the ztf high-cadence galactic plane survey: strategy, goals, and early results on new single-mode hot subdwarf b-star pulsatos" ?

The authors present the goals, strategy, and first results of the high-cadence Galactic plane survey using the Zwicky Transient Facility ( ZTF ). As part of this survey, the authors extract a total of ≈230 million individual objects with at least 80 epochs obtained during the high-cadence Galactic plane survey reaching an average depth of ZTF–r ≈ 20. The authors present a progress report on recent discoveries, including a new class of compact pulsators, the first members of a new class of Roche lobe filling hot subdwarf binaries as well as new ultracompact double white dwarfs and flaring stars. Finally, the authors present a sample of 12 new single-mode hot subdwarf B-star pulsators with pulsation amplitudes between ZTF–r = 20–76 mmag and pulsation periods between P = 5.