The Hyper Suprime-Cam SSP Transient Survey in COSMOS: Overview
Naoki Yasuda, Masaomi Tanaka, Nozomu Tominaga, Ji-an Jiang, Takashi J., Moriya, Tomoki Morokuma, Nao Suzuki, Ichiro Takahashi, Masaki S. Yamaguchi,, Keiichi Maeda, Masao Sako, Shiro Ikeda, Akisato Kimura, Mikio Morii, Naonori, Ueda, Naoki Yoshida, Chien-Hsiu Lee, Sherry H. Suyu

TL;DR
This paper provides an overview of a deep transient survey in the COSMOS field using Subaru HSC, identifying thousands of supernova candidates and enabling studies of high-redshift supernovae and other variable objects.
Contribution
It presents a comprehensive survey dataset that doubles the known number of high-redshift Type Ia supernovae and offers valuable data for time-domain astrophysics.
Findings
Identified 1,824 supernova candidates.
Classified 433 as Type Ia supernovae, with 129 having redshift data.
Doubled the number of known Type Ia SNe at z > 1.
Abstract
We present an overview of a deep transient survey of the COSMOS field with the Subaru Hyper Suprime-Cam (HSC). The survey was performed for the 1.77 deg ultra-deep layer and 5.78 deg deep layer in the Subaru Strategic Program over 6- and 4-month periods from 2016 to 2017, respectively. The ultra-deep layer shows a median depth per epoch of 26.4, 26.3, 26.0, 25.6, and 24.6 mag in , , , , and bands, respectively; the deep layer is mag shallower. In total, 1,824 supernova candidates were identified. Based on light curve fitting and derived light curve shape parameter, we classified 433 objects as Type Ia supernovae (SNe); among these candidates, 129 objects have spectroscopic or COSMOS2015 photometric redshifts and 58 objects are located at . Our unique dataset doubles the number of Type Ia SNe at and enables various time-domain analyses ofā¦
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
\jyear
2019 \Receivedreception date \Acceptedaccept date \Publishedpublication date
\KeyWords
supernovae: general ā cosmology: observations ā surveys
The Hyper Suprime-Cam SSP Transient Survey in COSMOS: Overview
Naoki Yasuda11affiliation: Kavli Institute for the Physics and Mathematics of the Universe (WPI), The University of Tokyo Institutes for Advanced Study, The University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8583, Japan **affiliationmark:
āā
Masaomi Tanaka22affiliation: Astronomical Institute, Tohoku University, Aoba, Sendai 980-8578, Japan 11affiliationmark:
āā
Nozomu Tominaga33affiliation: Department of Physics, Faculty of Science and Engineering, Konan University, 8-9-1 Okamoto, Kobe, Hyogo 658-8501, Japan 11affiliationmark:
āā
Ji-an Jiang44affiliation: Institute of Astronomy, Graduate School of Science, The University of Tokyo, 2-21-1 Osawa, Mitaka, Tokyo 181-0015, Japan
āā
Takashi J. Moriya55affiliation: National Astronomical Observatory of Japan, National Institutes of Natural Sciences, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan
āā
Tomoki Morokuma44affiliationmark: 11affiliationmark:
āā
Nao Suzuki11affiliationmark:
āā
Ichiro Takahashi11affiliationmark: 66affiliation: CREST, JST, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
āā
Masaki S. Yamaguchi33affiliationmark:
āā
Keiichi Maeda77affiliation: Department of Astronomy, Kyoto University, Kitashirakawa-Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan 11affiliationmark:
āā
Masao Sako88affiliation: Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA
āā
Shiro Ikeda99affiliation: Research Center for Statistical Machine Learning, The Institute of Statistical Mathematics, 10-3 Midori-cho, Tachikawa, Tokyo 190-8562, Japan 11affiliationmark:
āā
Akisato Kimura1010affiliation: NTT Communication Science Laboratories, 2-4 Hikaridai, Seika-cho, Keihanna Science City, Kyoto 619-0237, Japan
āā
Mikio Morii99affiliationmark:
āā
Naonori Ueda1010affiliationmark:
āā
Naoki Yoshida1111affiliation: Department of Physics, Graduate School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 11affiliationmark: 66affiliationmark:
āā
Chien-Hsiu Lee1212affiliation: Subaru Telescope, National Astronomical Observatory of Japan, 650 N Aāohoku Pl., Hilo, HI 96720, USA
āā
Sherry H. Suyu1313affiliation: Max-Planck-Institut für Astrophysik, Karl-Schwarzschild-Str. 1, 85748 Garching, Germany 1414affiliation: Physik-Department, Technische UniversitƤt München, James-Franck-StraĆe 1, 85748 Garching, Germany
āā
Yutaka Komiyama55affiliationmark: 1515affiliation: Graduate University for Advanced Studies (SOKENDAI), 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan
āā
Nicolas Regnault1616affiliation: Sorbonne UniversitƩ, UniversitƩ Paris Diderot, CNRS/IN2P3, Laboratoire de Physique NuclƩaire et de Hautes Energies, LPNHE, 4 Place Jussieu, 75252 Paris, France
āā
and David Rubin1717affiliation: Physics Division, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA 1818affiliation: Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA
Abstract
We present an overview of a deep transient survey of the COSMOS field with the Subaru Hyper Suprime-Cam (HSC). The survey was performed for the 1.77 deg2 ultra-deep layer and 5.78 deg2 deep layer in the Subaru Strategic Program over 6- and 4-month periods from 2016 to 2017, respectively. The ultra-deep layer shows a median depth per epoch of 26.4, 26.3, 26.0, 25.6, and 24.6 mag in , , , , and bands, respectively; the deep layer is mag shallower. In total, 1,824 supernova candidates were identified. Based on light curve fitting and derived light curve shape parameter, we classified 433 objects as Type Ia supernovae (SNe); among these candidates, 129 objects have spectroscopic or COSMOS2015 photometric redshifts and 58 objects are located at . Our unique dataset doubles the number of Type Ia SNe at and enables various time-domain analyses of Type II SNe, high redshift superluminous SNe, variable stars, and active galactic nuclei.
1 Introduction
Time-domain astronomy, now a major field of astronomy, continues to reveal fascinating variable and transient phenomena in the universe. Inspired by the discovery of dark energy through Type Ia supernovae (SNe Ia, [Riess etĀ al. (1998), Perlmutter etĀ al. (1999)]), many large and systematic surveys have been conducted over the last decade to test the cold dark matter (CDM) model with high precision (Astier etĀ al., 2006; Frieman etĀ al., 2008; Freedman etĀ al., 2009). These surveys demonstrate that the time-domain dataset is very rich, providing information on SNĀ Ia cosmology and other recent discoveries, such as the āsuper-Chandrasekharā supernova (SN) (Howell etĀ al., 2006), superluminous supernova (SLSN, Quimby etĀ al. (2007)), gravitationally lensed SN (Goobar etĀ al., 2017), and rapidly evolving transients (Drout etĀ al., 2014; Pursiainen etĀ al., 2018).
For SNĀ Ia cosmology, the latest Panoramic Survey Telescope and Rapid Response System 1 (Pan-STARRS1) survey () reports that the cumulative number of spectroscopically confirmed SNeĀ Ia is now 1,049 (Scolnic etĀ al., 2017), and the ongoing Dark Energy Survey (Bernstein etĀ al., 2012; Dark Energy Survey Collaboration etĀ al., 2016) () is about to add a few thousand SNeĀ Ia to the Hubble Diagram (DES Collaboration etĀ al., 2018). However, the number of SNeĀ Ia in high redshift () is still limited since the deep surveys have been possible only with the Hubble Space Telescope (HST), whose field of view is very small (Suzuki etĀ al., 2012; Riess etĀ al., 2018). We aim to probe the high redshift universe and trace the expansion history of the universe from the deceleration epoch through the acceleration epoch, to determine whether dark energy is time-variable or not. In this paper, we describe our transient survey using the 8-m Subaru telescope and highlight the expected scientific outcomes.
The Hyper-Suprime Cam (HSC) (Miyazaki etĀ al., 2018; Komiyama etĀ al., 2018; Furusawa etĀ al., 2018; Kawanomoto etĀ al., 2018) on the Subaru Telescope is the only instrument mounted on the prime focus among the large (8-10 m) telescopes. It is unique in its wide field of view (1.77 deg2), with 104 charge-coupled devices (CCDs; 4k2k pixels) providing a 0.168ā/pixel scale. The focal ratio of 1.83 makes the HSC the fastest camera among the large telescopes. An international team, composed of the astronomical communities of Japan, Taiwan and Princeton University, is in the process of completing a 300-night, 5-year HSC Subaru Strategic Program survey (HSC-SSP, 2014ā2019, Aihara etĀ al. (2018a)). Time-domain science is one of the main objectives of this SSP. In this paper, we provide an overview of the HSC transient survey for the COSMOS field.
The HSC transient survey is unique in its depth (26 mag) and volume (1.77 deg2 FoV), as it explores the deepest transients sky for the area deg2 (see Figure 1 for comparison with other transient surveys). The HSC-SSP survey has two ultra-deep fields: COSMOS (Cosmic Evolution Survey, Scoville etĀ al. (2007)) and SXDS (Subaru/XMM-Newton Deep Survey, Furusawa etĀ al. (2008)). A cadenced observation will be conducted on these two fields to achieve various transient science cases. In this paper, we will introduce some examples such as Type Ia and Type II-P SN cosmology and SLSNe.
This paper is organized as follows. In Section 2, we describe the survey strategy, data reduction, and transient findings. We show transient samples in Section 3 and present science highlights in Section 4. Finally, we give a summary in Section 5.
2 Overview of the Transient Survey
2.1 Observations
Observations were conducted as part of the HSC-SSP (Aihara etĀ al., 2018a) from November 2016 to April 2017 on the COSMOS field. As shown in Figure 2, there was one pointing of the ultra-deep layer (solid circle) and four pointings of the deep layer (dashed circles) surrounding the ultra-deep pointing, with significant overlap. The basic observing strategy was to obtain two epochs separated by 7ā10 days in all five broad bands (, , , , and -bands) during each monthly observation run for the ultra-deep layer. In total, we acquired 12 epochs over 6 months for each of five bands. For the deep layer, six epochs were obtained over 4 months, with shorter exposure times than the ultra-deep run, as the total planned exposure time was limited. During this transient survey, wide-layer observation around the COSMOS field was also conducted and the data obtained were included in the analysis. Detailed observation dates and typical exposure times are summarized in Table 2.1. Each image was taken with five-point dithering, as described in Aihara etĀ al. (2018a), to fill the gaps between the HSC CCDs. There are overlaps between ultra-deep and deep pointings, as mentioned earlier; additionally, the exposure time and other statistics vary from one position to the next. In Table 2.1, values around the center of the pointings are presented. Seeing values measured on coadded images are listed in Table 2.1 and are shown in Figure 3. {longtable}lllārrrārrr Observation log. Ultra-deep layer (1.77 deg2) Deep layer (5.78 deg2)
Obs. date MJD filter exptime seeing lim. mag exptime seeing lim. mag
(sec) (arcsec) (mag) (sec) (arcsec) (mag)
\endfirsthead Ultra-deep Deep
Obs. date MJD filter exptime seeing depth exptime seeing depth
(sec) (arcsec) (mag) (sec) (arcsec) (mag)
\endhead\endfoot\endlastfootref 57177.75 HSC-G 7500 0.81 900 0.83
ref 57087.75 HSC-R 8460 0.63 1260 0.66
ref 57048.44 HSC-I 11730 0.60 2700 0.57
ref 57150.73 HSC-Z 28260 0.59 3510 0.59
ref 57150.19 HSC-Y 12960 0.63 1890 0.62
2016-11-23 57715.54 HSC-Z 4200 0.71 25.64
2016-11-23 57715.62 HSC-Y 4500 0.72 25.24
2016-11-25 57717.57 HSC-G 1800 1.09 25.66
2016-11-25 57717.62 HSC-I2 3000 0.80 26.01
2016-11-28 57720.60 HSC-R2 1800 0.76 26.66
2016-11-29 57721.55 HSC-I2 2400 1.15 25.81
2016-11-29 57721.60 HSC-Z 3000 1.04 25.47
2016-12-23 57745.56 HSC-Z 5440 1.05 25.32 2070 1.04 24.83
2016-12-25 57747.53 HSC-R2 1560 1.12 25.88 540 1.14 25.45
2016-12-25 57747.62 HSC-I2 3840 1.23 25.64 810 1.20 25.01
2016-12-26 57748.53 HSC-Y 3540 1.48 22.88 810 1.39 22.36
2017-01-02 57755.45 HSC-Z 3300 0.73 25.59
2017-01-02 57755.51 HSC-I2 3000 0.68 26.51 600 0.67 25.64
2017-01-02 57755.61 HSC-G 1710 0.69 26.75 840 0.68 26.21
2017-01-04 57757.52 HSC-Y 4640 0.78 24.69 1610 0.69 24.27
2017-01-21 57774.50 HSC-Z 5240 0.52 26.31 2650 0.54 25.68
2017-01-23 57776.41 HSC-R2 2280 0.83 26.44 1440 0.84 25.98
2017-01-23 57776.54 HSC-I2 3470 0.70 26.43 1510 0.70 25.84
2017-01-25 57778.45 HSC-G 6480 1.77 26.13 2520 1.70 25.60
2017-01-25 57778.62 HSC-Y 3570 1.13 23.86 810 1.19 23.20
2017-01-26 57779.52 HSC-Z 400 0.73 24.36 200 0.74 24.28
2017-01-30 57783.43 HSC-I2 2670 0.74 26.03 810 0.75 25.31
2017-01-30 57783.55 HSC-Z 5250 0.65 25.85 1350 0.64 25.22
2017-02-01 57785.39 HSC-G 3540 0.66 26.38 1440 0.66 25.77
2017-02-02 57786.45 HSC-R2 1380 0.65 26.60 540 0.63 25.97
2017-02-02 57786.59 HSC-I2 800 0.49 25.83 600 0.48 25.65
2017-02-03 57787.47 HSC-Y 10710 1.20 24.71 2430 1.25 23.96
2017-02-21 57805.37 HSC-Z 3840 0.64 25.69 1350 0.69 25.04
2017-02-23 57807.37 HSC-G 3660 1.40 26.30 1440 1.31 25.86
2017-02-23 57807.48 HSC-R2 1680 0.91 26.33 720 0.95 25.77
2017-02-25 57809.40 HSC-I2 10350 0.75 25.85 3240 0.64 25.21
2017-02-27 57811.40 HSC-Y 1800 0.79 24.44 270 0.88 23.70
2017-03-04 57816.31 HSC-Z 3840 0.64 25.73 1080 0.68 24.81
2017-03-04 57816.47 HSC-I2 3210 0.67 26.33 810 0.65 25.68
2017-03-06 57818.51 HSC-R2 1740 0.73 26.47 540 0.74 25.89
2017-03-07 57819.47 HSC-Y 5110 0.51 25.21 1810 0.59 24.45
2017-03-21 57833.38 HSC-Y 4110 0.60 24.97 810 0.55 24.21
2017-03-22 57834.32 HSC-G 2490 0.84 26.74 1215 0.80 26.22
2017-03-22 57834.43 HSC-Z 4380 0.56 25.82 1350 0.59 25.12
2017-03-23 57835.26 HSC-I2 2340 0.67 25.89 540 0.65 25.30
2017-03-25 57837.27 HSC-R2 2220 0.98 26.11 720 0.98 25.58
2017-03-25 57837.49 HSC-Y 1500 1.80 22.86
2017-03-29 57841.29 HSC-G 1740 0.92 26.53 540 1.01 25.73
2017-03-29 57841.41 HSC-Z 3300 0.74 25.60
2017-03-30 57842.27 HSC-I2 3600 0.98 26.02
2017-03-30 57842.34 HSC-Y 4500 1.02 24.64
2017-04-01 57844.33 HSC-R2 3300 1.18 26.13
2017-04-20 57863.28 HSC-Y 4200 0.62 24.42
2017-04-23 57866.25 HSC-R2 2400 0.94 26.10
2017-04-23 57866.36 HSC-Z 3600 0.81 25.32
2017-04-26 57869.27 HSC-I2 4800 1.25 25.70
2017-04-26 57869.33 HSC-G 1200 0.88 26.45
2017-04-27 57870.35 HSC-I2 1800 0.55 26.09
2017-04-29 57872.26 HSC-Z 3300 0.74 25.30
2017-06-20 57924.28 HSC-Z 900 1.15 23.95
2.2 Data reduction
HSC pipeline (Bosch etĀ al., 2018) version 4.0.5, with default configuration parameters, was used for the data reduction. We applied the same reduction procedure with the HSC-SSP data release (Aihara etĀ al., 2018b) for the standard image reduction, which included bias, dark, flat, and fringe corrections, as well as astrometric and photometric calibrations against the PS1 catalog (Magnier etĀ al., 2013). Typical photometric accuracy is 1%-2% (Aihara etĀ al., 2018b). Based on the astrometric solutions, the images were warped into a predefined sky grid, which we refer to as a warped image. Data were processed on a nightly basis.
For the image subtraction, the difference imaging method of Alard & Lupton (Alard & Lupton, 1998; Alard, 2000) was applied using reference deep, coadded images created from data taken during March 2014 and April 2016. Images with seeing better than 0.7 arcsec were used as reference images, with the exception of the -band. Table 2.1 also includes exposure time and seeing of reference images. Difference imaging was conducted for each warped image, and warped difference images were coadded to create deep difference images for each filter and epoch. With this method, we can avoid subtraction error caused by a discrete change in the point spread function (PSF) at the CCD gaps in coadded images.
Note that - and -band filters were replaced with filters having a more spatially uniform response across the focal plane in 2016 July and February, respectively (Aihara etĀ al., 2018a). Thus, reference images were observed with the old filters, and search images were observed with the new filters. The non-uniformity in the old filters can result in 4-5% magnitude offset for high-redshift SNe Ia, whereas the new filters reduced the offset to less than 1%. A detailed software patch is currently under development to match the precision required for SN cosmology. However, the offset described did not affect any of the results presented in this paper.
Once coadded difference images are created, we can detect and measure sources in difference images. Based on these sources, transient sources can be identified (see section 2.4 for details). Forced photometry was performed at the location of transients for images of all filters and epochs. Here, the location of a transient was measured by taking the mean of positions in each coadded difference image, in which a transient was detected.
2.3 Limiting magnitude/detection efficiency
To estimate the limiting magnitude of each epoch image, we injected artificial stars with magnitudes ranging from 21 to 28 mag in processed CCD images, before warping to the predefined grid. The number density of artificial stars was 20,000 / deg2, which corresponds to about 400 objects per CCD. The corresponding CCD images were processed in the same way as the real data. The source catalog from difference coadded images was then compared with the input artificial star catalog. Figure 4 shows the magnitudes at detection efficiency for each filter as a function of the observing epoch. Figure 5 shows the spatial distribution of the limiting magnitude of the -band on a specific night. Image subtraction examples for artificial stars embedded onto galaxies are shown in FigureĀ 6; the image subtraction technique worked well for them also.
For transient findings, we imposed at least two detections at the same position (see section 2.4). Figure 7 shows the detection rate of constant brightness objects as a function of magnitude. Different lines denote how many times the same object was detected at the same position.
2.4 Transient finding
Sources detected on coadded difference images have been classified as real or bogus by machine learning techniques. We adopted a convolutional neural network (CNN) with a combination of convolution, pooling, and dropout layers, trained by 100,000 artificial stars as a real sample and 100,000 objects within actual observational images as a bogus sample. The trained CNN was validated with 10,000 artificial stars and 10,000 bogus samples. Note that here, a bogus sample includes ārealā transients, as they are taken from actual observational data. The CNN showed a false-positive rate of 4.3% and 6.0% at the true-positive rate of 90% for objects with a signal-to-noise ratio better than 10 and 7.5, respectively. These values are not necessarily better than those cited in a previous study for the HSC data (Morii etĀ al., 2016) using various machine learning methods for measured parameters. This is because our bogus sample was constructed from actual observational data and included a large number of ārealā transients, as the reference images are taken well before the search observation. The actual performance is expected to be better than that indicated by these values. We also applied the area under the curve (AUC) boosting method (Morii etĀ al., 2016) and partial AUC optimization (Ueda & Fujino, 2018) for real/bogus classification. The results are discussed in section 3.1.
After the CNN screening, if the same source was identified at the same place (within arcsec), then that source was registered as a transient. In total, 65,387 transient candidates were identified. For the registered transients, the closest object in reference images was tagged as the host object. Sometimes, the host object identification is not optimal (matched with very faint noise-like objects or matched with deblended children of big galaxies); in such cases, any clear misidentifications were corrected by visual inspection. The host objects were then matched with a compilation of public redshift catalogs, COSMOS2015 (Laigle etĀ al., 2016), and HSC photo-z catalog (Tanaka etĀ al., 2018) objects with a search radius of arcsec. Public redshift catalogs included SDSS DR12 (Alam etĀ al., 2015), PRIMUS DR1 (Coil etĀ al., 2011; Cool etĀ al., 2013), VVDS (Le FĆØvre etĀ al., 2013), zCOSMOS DR3 (Lilly etĀ al., 2009), and 3D-HST (Skelton etĀ al., 2014; Momcheva etĀ al., 2016).
3 Sample of Transients
3.1 Classification of transients
We mainly focus on extragalactic transients and SNe in this paper. To identify SN candidates, 65,387 candidates were classified according to their properties after CNN screening. The flowchart of the classification is shown in FigureĀ 8. First, we excluded objects having light curves dominated by negative PSF fluxes. When a light curve of a certain filter had more negative points than positive points, it was judged to be a negative light curve. If the number of filters with a negative light curve exceeded that with a positive light curve, the object was excluded as a āNegativeā candidate. Second, we checked host objects in reference images. If the host object was a point source, we excluded it, as it was most likely to be a variable star. Note that some active galactic nuclei (AGN) are classified as āNegativeā and āPoint source in referenceā. Finally, visual inspection was performed on the remaining 26,988 candidates by nine experts. The light curve shapes were also visually checked, along with the time series of coadded difference images.
In the visual inspection, we first excluded any bogus detections, which accounted for a large proportion of the remaining candidates (Figure 9). Some of the bogus detections are caused by imperfect image subtraction near bright stars (see HSC17dsge in Figure 9); this occurs because they are not classified as a āPoint source in referenceā if the offset from the center of the star is relatively large. A large portion of these objects was detected in only one or two bands, with nearly constant fluxes (see HSC17bknx); this may occur due to imperfect subtraction. Other bogus detection examples include artifacts in the reference image (HSC17dshg) and too few detections with low signal-to-noise ratios (HSC17brko).
We then classified the clean objects into SN or AGN. When a candidate had a clear offset ( of centroid error) with a host object, it was classified as a SN. AGN can be identified by their association with an X-ray source using the COSMOS2015 catalog. However, because some SNe can occur in X-ray bright galaxies and low-luminosity AGN can elude X-ray detection, we mainly classified AGN based on the light curve shapes (e.g., multiple peaks or a very long duration; see Figure 10). Ultimately, 1,824 objects were classified as SN, and 1,534 as AGN. Marginal cases (139 objects) were flagged as āSN/AGNā. Figure 11 shows the distributions of the apparent distance to the host object. AGN were almost exclusively located within arcsec from the host objects. On the other hand, SNe were distributed more widely. Note that our AGN samples are not complete. Complete AGN samples including candidates classified as āNegativeā and āPoint source in referenceā will be presented in forthcoming papers.
Although visual inspection was performed on all of the objects following CNN screening and the exclusion of āNegativeā and āPoint source in referenceā objects, we also applied the real/bogus classification process with AUC boosting (Morii etĀ al., 2016) and partial AUC optimization, using a deep neural network (Ueda & Fujino, 2018). In Table 2.3, we show the results only for classified objects to compare the reduction factors of SNe, AGN, and bogus objects. By applying AUC boosting and partial AUC optimization, in addition to CNN, the number of bogus objects was reduced to 56% and 70%, while more than 94% and 94% of real objects (SNe and AGN) were preserved, respectively. When we applied both methods, the bogus objects were reduced to 50%, with more than 91 % of the real objects retained. This highlights the usefulness of multiple real/bogus classifiers, as demonstrated in Morii etĀ al. (2016).
3.2 Supernova candidates
In this section, we examine the properties of objects classified as SN. Figure 12 shows the distribution of observed -band magnitudes at the peak. Here āpeakā refers to the brightest magnitude within the observation and not the result of light curve fitting. For the āSNe Iaā sample described below, the difference between this āpeakā magnitude and the peak magnitude obtained from the light curve fit shows a -mag offset (the āpeakā magnitude is fainter) and -mag scatter. A majority of our samples exhibit ā mag at the peak, with a tail to ā mag. Thus, the HSC-SSP transient survey is among the deepest transient surveys (see also Figure 1), detecting a larger number of SNe than SN surveys with the HST (Dawson etĀ al., 2009; Rodney etĀ al., 2014; Graur etĀ al., 2014).
With deep depth, these transients are located at high redshifts. By virtue of the rich dataset in the COSMOS field, 207 and 371 SN candidates have spectroscopic redshifts (hereafter, spec-z) and COSMOS2015 photometric redshifts (COSMOS photo-z), respectively, by identifying potential host galaxies of the SN candidates. Figure 13 shows the redshift distribution of these 578 objects. The distribution has a median of ; 187 objects (32%) are located at . Images of example SNe in each redshift range are displayed in Figure 14. For 141 of the SN candidates, we were unable to identify a clear host galaxy. A detailed analysis of hostless samples will be described in a separate paper.
For classification of SNe Ia and other types of SNe, we applied the SALT2 light curve fitter (Guy etĀ al., 2007) to our SN sample. SALT2 is an empirical model of SN Iaās spectro-photometric evolution over time. The model is constructed using a large data set that includes light curves and spectra of both nearby and distant SNe, up to redshift 1. Available ultraviolet (UV) spectra from the International Ultraviolet Explorer (IUE) are also included. The model is valid in the spectral range of Ć . The public version of the SALT2 package111http://supernovae.in2p3.fr/salt/ provides light curve fitting algorithms based on this model. Light curve fitters can estimate (time of peak brightness in the -band), redshift, (the SN color parameter of the model), (the flux scale or luminosity distance), and (the light curve shape parameter of the model) to fit the observed multi-band photometric data by minimizing . A commonly used fitting tool, snfit, fixes the redshift to a given value, allowing the other four parameters to be determined. Another fitter, snphotoz, fits all five parameters. At the first stage of snphotoz, both and are fixed to 0 (mean values of SN Ia), and then the redshift is scanned to find the minimum value. This redshift is then used as an initial value for a full five-parameter fit. One may question the validity of using the SALT2 model at a redshift beyond 1. Balland etĀ al. (2018) found no evidence of redshift evolution from Very Large Telescope (VLT) spectra of SNe Ia below a redshift of 1.0. There has been no observational study clearly showing the evolution of SN Ia beyond redshift 1, mainly because the observations are limited by the faintness of the objects. In this paper, we simply assumed that the SALT2 model is valid beyond redshift 1.
For all SNe, we ran snfit by constraining the redshift to the best available redshift. For SNe that do not associate with a clear host galaxy (hostless), or which have only photometric redshifts from HSC broad-band photometry, we also ran snphotoz to estimate both the redshift and light curve parameters. For both snfit and snphotoz, we added an option ā-w 2500 8000ā to use a wider wavelength range. When the reduced of snphotoz was less than 70% of the reduced of snfit, we adopted the result of snphotoz. Note that the -band images suffer from scattered light, and some of the objects are affected by imperfect correction (Aihara etĀ al., 2018b). Therefore, we have not used -band data for light curve fits, to ensure a similar classification process for all SN candidates.
With regard to defining āSNe Iaā samples in this paper, we selected SNe with the following characteristics: (1) light curve parameters, color (), and shape () within the range of Scolnic & Kessler (2016)ās āAll G10ā distribution, (2) a brighter than mag, (3) a reduced of less than 10, and (4) the number of degrees of freedom is greater than or equal to 5. In total, 433 SNe were classified as āSNe Iaā. Among them, 129 SNe have spec-z or COSMOS photo-z; the relatively small fraction of these with respect to the total number of SN is mainly due to the large area outside of the original COSMOS field, as shown in Figure 2. For 57 of the SNe with either incorrect or unavailable redshifts, the redshift was recovered by snphotoz. Figure 15 shows representative light curves of āSNe Iaā at different redshifts.
As described above, our āSN Iaā sample was selected solely from photometric information; no SNe spectroscopic information was used. Thus, this sample may include contamination from other SNe types or may be missing genuine SNe Ia. In this paper, we do not focus on detailed SN classification, as more detailed SN classification will be presented in other papers. The number of āSN Iaā (433) is relatively small compared with the entire sample number (1,824). This is due to our conservative criteria. Infact, there are 232 SNe with a small number of available data points, which do not satisfy condition (4) above. If we loosen condition (1), corresponding to the light curve parameters condition of being within to being within , as well as condition (3), in which the condition is reduced from less than 10 to less than 20, 104 additional SNe can be assigned to āSN Iaā. This results in as the SNe Ia fraction. Note that this value is still lower than the general expectation of about . The detailed classification algorithms differ; photometric identifications of SNe Ia in SDSS (Sako etĀ al., 2011) and Pan-STARRS1 (Jones etĀ al., 2017) SN surveys show a similarly low SN Ia fraction.
Distributions of absolute magnitudes are shown in Figure 16; absolute magnitudes as a function of redshift are shown in Figure 17. In these figures, we only show SNe with spec-z and COSMOS photo-z. The absolute magnitudes of āSNe Iaā are clustered at to mag. The redshift distribution of āSNe Iaā has a median of , and 58 objects (45%) are located at .
Note that some SNe, which are not classified as āSNe Iaā, have absolute magnitudes brighter than mag. Although some of these objects are real (such as bright Type IIn SNe or SLSNe, see Section 4), it should be noted that most of non-āSNe Iaā at have only photometric redshifts (Figure 17), and the magnitude distributions are affected by the uncertainties of the photometric redshifts. In fact, if we limit ourselves to using only spectroscopic redshifts, the absolute magnitudes of non-āSNe Iaā have a steep cutoff around mag; there is no object with an absolute magnitude of mag.
4 Science highlights
We discuss some of the latest findings obtained using the HSC transient survey data presented in this paper.
4.1 Type Ia SN cosmology
Two decades have passed since the discovery of dark energy (Perlmutter etĀ al., 1999; Riess etĀ al., 1998), yet its nature still remains as one of the biggest mysteries in modern physics. Today, baryon acoustic oscillation augments the measurement of dark energy, but SNĀ Ia still leads with respect to measurement precision. Now we would like to know whether dark energy changes in time or not. The HSC transient survey is designed to probe redshifts in where time variability becomes sensitive.
We identified 433 āSNeĀ Iaā, based on the light curve and photometric redshifts from host galaxies. In particular, we found 58 āSNe Iaā beyond redshift , with reliable spec-z or COSMOS photo-z. In the past, only the HST could reach this redshift range, and only two dozen SNĀ Ia have been measured (Suzuki etĀ al., 2012; Riess etĀ al., 2018). HSC is the only instrument that can probe SNe from the ground having the necessary photometric accuracy for cosmological analysis. By doubling the number of high-redshift SNeĀ Ia, we expect to impose a tight constraint on the nature of dark energy, in which the cosmological parameter becomes sensitive with high redshifts.
We conducted a spectroscopic follow-up campaign for live SNe with large telescopes: Subaru/Faint Object Camera and Spectrograph (FOCAS), Keck/Low Resolution Imaging Spectrometer (LRIS), VLT/visual and near-UV FOcal Reducer and low dispersion Spectrograph (FORS), Gemini/Gemini Multi-Object Spectrographs (GMOS), and Grand Telescopio CANARIAS (GTC)/Optical System for Imaging and low Intermediate Resolution Spectroscopy (OSIRIS). During the 2017 season, in collaboration with the COSMOS Lyman-Alpha Mapping And Tomography Observations (CLAMATO) project team (Lee etĀ al., 2014), we placed a few slits on live SNe on the Keck/LRIS mask, while other slits were used for Lyman break galaxies. We observed 17 live SN spectra; the details will be reported in a forthcoming paper.
We also performed host galaxy spectroscopic follow-up observations. First, we used the COSMOS2015 catalog to determine if spec-z was known; if so, we adopted that redshift (Laigle etĀ al., 2016). When slit mask observation was conducted, we placed slits on potential host galaxies in the field of view for Keck/LRIS, Subaru/FOCAS, and VLT/FORS. The most efficient observation was conducted by the 3.9-m Anglo-Australian Telescope (AAT)/AAOmega spectrograph, which has 400 fibers in the two-square-degree field of view. In the 2018 February run, we collected 257 host galaxy spectra, with the goal of completing the collection in the upcoming semesters.
Although the HSC can detect SNĀ Ia beyond , the peak flux of SNĀ Ia goes into the infrared (IR), and HSC loses its sensitivity. We conducted an IR imaging follow-up observation via the HST. For cosmological analysis, we introduced light curve width (stretch) correction and color correction (Tripp, 1998), and it is essential to reduce error propagation from colors. Together with HSC data, we designed our HST IR follow-up observation to reduce color-associated error to less than 3%. High-z SNĀ Ia candidates were identified from HSC observations. An observing request was sent to the HST. We successfully observed 26 SNeĀ Ia candidates () with the HSC and HST, but have yet to observe the HST reference images. We hope to collect spectroscopic redshifts of the host galaxies for cosmological analyses in the near future.
4.2 Type II SNe
Type II SNe (SNe II) constitute the most common class among core-collapse SNe (Li etĀ al., 2011), tracing the most typical evolutionary path of massive stars. Characterized by hydrogen features, they are robustly interpreted as an explosion of a red supergiant. Analysis of their light curve properties, which thus far have been limited mostly to the local sample, provides a clue as to the nature of the progenitor evolution and actively debated explosion mechanism (Anderson etĀ al., 2014). In addition, there has been increasing interest in using SNe II for cosmological applications (Hamuy & Pinto, 2002; Nugent etĀ al., 2006; Olivares E. etĀ al., 2010; de Jaeger etĀ al., 2017a)
The SNe light curves are characterized by a rapid increase to the peak value ( days), followed by a slow decline, frequently showing a plateau, for days (GonzÔlez-GaitÔn et al., 2015). As such, the COSMOS HSC transient survey is suited to SNe II discovery and characterization of the entire evolution of their multi-color light curves. The peak in the redshift distribution is expected to occur at , with extension to , far beyond the well-established available SN II samples.
In de Jaeger etĀ al. (2017b), we used one SN II discovered by the COSMOS HSC transient survey to extend the SN II Hubble diagram to . We applied the so-called standard candle method, which uses a correlation between the intrinsic luminosity and the expansion velocity from the spectra (Hamuy & Pinto, 2002; Nugent etĀ al., 2006). Followed by the rapid rise and a plateau-like evolution of SN 2016jhj, as our best high-redshift SN II candidate with the COSMOS transient survey, spectroscopy was performed with Keck/LRIS; the LRIS spectrum confirmed it to be an SN II. Cross correlation to SN II template spectra was performed to measure the expansion velocity, in an attempt to place this object within the luminosityā expansion velocity relationship (Poznanski etĀ al., 2009). The standardized magnitude, as well as the color correction (due to extinction), were then modeled with the cosmological parameters using a Monte Carlo Markov Chain simulation, to derive the best-fit Hubble diagram and the probability distribution of the parameters. The derived cosmological parameters are consistent with the -CDM model. The resulting dispersion in the standardized magnitude was 0.27 mag (i.e., 12-13% in distance). This work represents a proof-of-concept of the capability of high-redshift SN II cosmology ().
We plan to extend the analysis to the photometric color method (de Jaeger etĀ al., 2015, 2017a), which uses only photometric information to derive the distances to SNe II. Minimal requirements of this method include light curve information from two bands, the slope of the plateau in a given band pass, and a color term. This methodology fits well to the HSC sample, which consists of a number of (apparently) faint SNe II for which the spectroscopic follow-up is difficult to coordinate. Given the quality of the multi-band light curves in the HSC sample, both in terms of photometric accuracy and coverage over the entire survey duration, the application of this method is straightforward. This approach will make the best use of the large sampling of HSC-discovered high-redshift SNe II, to derive the Hubble diagram up to .
4.3 Superluminous SNe
We searched for high-redshift superluminous SNe (SLSNe) in this survey. High-redshift SLSN candidates were selected based on the photometric redshifts of host galaxies. We mainly used the photometric redshifts provided in the COSMOS database (Laigle etĀ al., 2016).
The COSMOS HSC transient survey has led to the discovery of nearly ten high-redshift SLSN candidates. Among them, we have thus far successfully confirmed redshifts in three high-redshift SLSNe (Curtin etĀ al., 2018). The identified redshifts are , , and . Unfortunately, the spectra are not good enough to identify spectroscopic type. There are several candidates with higher host photometric redshifts, and , whose spectra have not been obtained (Moriya etĀ al., 2018).
Based on the three SLSNe at , Moriya etĀ al. (2018) estimated SLSN rates at of . This rate, based solely on spectroscopically confirmed SLSNe, is already comparable to the total SLSN rate at , estimated by extrapolating the local SLSN rate (Quimby etĀ al., 2013) using the cosmic star-formation history (Madau & Dickinson, 2014). The SLSN rate at may be higher if we take SLSN candidates without spectroscopic confirmation into account.
5 Summary
We performed a deep transient survey with HSC of the Subaru telescope. The ultra-deep layer, the central 1.77 deg2 of the COSMOS field with one HSC pointing, was observed repeatedly for 6 months in 2016 and 2017 with , , , , and filters, while the deep layer, 5.78 deg2 with four HSC pointings, was observed for 4 months. For each month, data were taken at two epochs per filter. The limiting magnitudes per epoch for the ultra-deep layer are as follows: , , , , and mag in the -, -, -, -, and -band, respectively; Deep-layer values are roughly 0.6 mag shallower. The data set obtained is one of the deepest wide-field transient surveys attempted to date.
In the dataset, 1,824 SN candidates were identified. Among our samples, 207 and 371 objects have spec-z and COSMOS photo-z , respectively. The median redshift is , and 187 objects (32%) are located at . This is among the largest high-redshift SN samples. By light curve fitting using SALT2, 433 (129 with spec-z or COSMOS photo-z) candidates were classified as āSNe Iaā. In particular, 58 objects are located at a redshift beyond . More dedicated photometric classification will be presented in a forthcoming paper.
Our dataset doubles the number of Type Ia SNe candidates at for Type Ia SN cosmology, and the great depth enables a search for SLSNe at even higher redshifts (Moriya etĀ al., 2018; Curtin etĀ al., 2018). The survey also provides Type IIP SNe at medium redshift, which has been demonstrated by the highest-redshift Type IIP SNe for cosmological use (de Jaeger etĀ al., 2017b). In addition to the transient science, deep time-series images can also be used for studies of variable stars and AGN. This survey of the COSMOS field is the first half of the HSC-SSP transient survey. A similar transient survey for the SXDS field will be conducted over a period of 6 months, starting in September 2019. More details regarding the science topics to be covered, as well as the results from the forthcoming SXDS survey, will be presented in separate papers.
{ack}
The Hyper Suprime-Cam (HSC) collaboration includes the astronomical communities of Japan, Taiwan, and Princeton University. The HSC instrumentation and software were developed by the National Astronomical Observatory of Japan (NAOJ), the Kavli Institute for the Physics and Mathematics of the Universe (Kavli IPMU), the University of Tokyo, the High Energy Accelerator Research Organization (KEK), the Academia Sinica Institute for Astronomy and Astrophysics in Taiwan (ASIAA), and Princeton University. Funding was contributed by the FIRST program from the Japanese Cabinet Office, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), the Japan Society for the Promotion of Science (JSPS), the Japan Science and Technology Agency (JST), the Toray Science Foundation, NAOJ, Kavli IPMU, KEK, ASIAA, and Princeton University.
The Pan-STARRS1 Surveys (PS1) have been made possible through contributions by the Institute for Astronomy, the University of Hawaii, the Pan-STARRS Project Office, the Max-Planck Society and its participating institutes, the Max Planck Institute for Astronomy, Heidelberg and the Max Planck Institute for Extraterrestrial Physics, Garching, Johns Hopkins University, Durham University, the University of Edinburgh, Queenās University Belfast, the Harvard-Smithsonian Center for Astrophysics, the Las Cumbres Observatory Global Telescope Network Incorporated, the National Central University of Taiwan, the Space Telescope Science Institute, the National Aeronautics and Space Administration (under Grant No. NNX08AR22G issued through the Planetary Science Division of the NASA Science Mission Directorate), the National Science Foundation (under Grant No. AST-1238877), the University of Maryland, and Eotvos Lorand University (ELTE).
This study used software developed for the Large Synoptic Survey Telescope (LSST). We thank the LSST Project for making their code available as free software at http://dm.lsst.org.
IT and NY acknowledge financial support from JST CREST (JPMHCR1414). MT is supported by an Inoue Science Research Award from the Inoue Foundation for Science and the Grant-in-Aid for Scientific Research programs of JSPS (15H02075, 16H02183) and MEXT (17H06363).
This research is based in part on data collected at the Subaru Telescope and retrieved from the HSC data archive system, which is operated by the Subaru Telescope and Astronomy Data Center at NAOJ.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Aihara et al. (2018 a) Aihara, H., et al. 2018 a, PASJ, 70, S 4
- 2Aihara et al. (2018 b) ā. 2018 b, PASJ, 70, S 8
- 3Alam et al. (2015) Alam, S., et al. 2015, Ap JS, 219, 12
- 4Alard (2000) Alard, C. 2000, A&AS, 144, 363
- 5Alard & Lupton (1998) Alard, C., & Lupton, R. H. 1998, Ap J, 503, 325
- 6Anderson et al. (2014) Anderson, J. P., et al. 2014, Ap J, 786, 67
- 7Astier et al. (2006) Astier, P., et al. 2006, A&A, 447, 31
- 8Balland et al. (2018) Balland, C., et al. 2018, A&A, 614, A 134
