THE YOUNG AND BRIGHT TYPE IA SUPERNOVA ASASSN-14lp: DISCOVERY, EARLY-TIME OBSERVATIONS, FIRST-LIGHT TIME, DISTANCE TO NGC 4666, AND PROGENITOR CONSTRAINTS

, , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , and

Published 2016 July 27 © 2016. The American Astronomical Society. All rights reserved.
, , Citation B. J. Shappee et al 2016 ApJ 826 144 DOI 10.3847/0004-637X/826/2/144

0004-637X/826/2/144

ABSTRACT

On 2014 December 9.61, the All-sky Automated Survey for SuperNovae (ASAS-SN or "Assassin") discovered ASASSN-14lp just ∼2 days after first light using a global array of 14 cm diameter telescopes. ASASSN-14lp went on to become a bright supernova (V = 11.94 mag), second only to SN 2014J for the year. We present prediscovery photometry (with a detection less than a day after first light) and ultraviolet through near-infrared photometric and spectroscopic data covering the rise and fall of ASASSN-14lp for more than 100 days. We find that ASASSN-14lp had a broad light curve (${\rm{\Delta }}{m}_{15}(B)=0.80\pm 0.05$), a B-band maximum at 2457015.82 ± 0.03, a rise time of ${16.94}_{-0.10}^{+0.11}$ days, and moderate host-galaxy extinction ($E{(B-V)}_{\mathrm{host}}=0.33\pm 0.06$). Using ASASSN-14lp, we derive a distance modulus for NGC 4666 of $\mu =30.8\pm 0.2$, corresponding to a distance of 14.7 ± 1.5 Mpc. However, adding ASASSN-14lp to the calibrating sample of Type Ia supernovae still requires an independent distance to the host galaxy. Finally, using our early-time photometric and spectroscopic observations, we rule out red giant secondaries and, assuming a favorable viewing angle and explosion time, any nondegenerate companion larger than 0.34 ${R}_{\odot }$.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

Even though Type Ia supernovae (SNe Ia) have been used to discover the accelerating expansion of the universe (Riess et al. 1998; Perlmutter et al. 1999), the nature of their progenitor systems remains unknown (for a review see Maoz et al. 2014). It is generally accepted that SNe Ia involve the thermonuclear explosion of a C/O white dwarf (WD), but basic details of this picture, such as the nature of the binary companion and the sequence of events leading to the SN explosion, are still hotly debated. Recently, significant observational (e.g., Gilfanov & Bogdán 2010; Li et al. 2011; Horesh et al. 2012; Schaefer & Pagnotta 2012; Shappee et al. 2013b) and theoretical progress (e.g., Kasen 2010; Seitenzahl et al. 2013; Shappee et al. 2013a; Piro & Nakar 2014; Scalzo et al. 2014) has been made to constrain SN Ia progenitor systems, but a clear picture has yet to emerge.

Detailed, early-time observations are one of the more promising avenues to further constrain the progenitors of SNe Ia. Statistical studies of the rise times (${t}_{\mathrm{rise}}$) of SNe Ia have a long history (Riess et al. 1999; Conley et al. 2006; Strovink 2007; Hayden et al. 2010a, 2010b; Ganeshalingam et al. 2011; González-Gaitán et al. 2012; Firth et al. 2015), but it has been the in-depth, early-time observations of nearby SNe Ia that have improved our understanding of the progenitor systems in recent years (Nugent et al. 2011b; Bloom et al. 2012; Foley et al. 2012; Silverman et al. 2012; Zheng et al. 2013, 2014; Yamanaka et al. 2014; Goobar et al. 2015; Olling et al. 2015; Siverd et al. 2015). The most notable example is SN 2011fe, which was discovered only hours after first light (${t}_{\mathrm{first}};$ Nugent et al. 2011a) by the Palomar Transient Factory (Law et al. 2009) at 6.4 Mpc (Shappee & Stanek 2011). Bloom et al. (2012) used a nondetection of SN 2011fe just 4 hr after ${t}_{\mathrm{first}}$ to strongly constrain the radius of the primary (${R}_{* }\lesssim 0.02$ ${R}_{\odot }$) and the radius of a companion (${R}_{{\rm{c}}}\lesssim 0.1$ ${R}_{\odot }$). These constraints rule out all commonly considered nondegenerate companions and show directly, for the first time, that the exploding star for SN 2011fe was a degenerate object. However, there might be multiple channels for producing normal SNe Ia (e.g., Maguire et al. 2013), highlighting the need to build a larger sample of well-constrained SNe Ia.

With the addition of ASASSN-14lp there are now four SNe Ia with a clear detection within 1 day of first light. ASASSN-14lp was discovered just ${2.23}_{-0.10}^{+0.11}$ days after first light, and we present multiple prediscovery detections. In Section 2, we describe our discovery and observations. In Section 3, we fit the light curve of ASASSN-14lp and determine the distance to its host galaxy. In Section 4, we use the early-time observations to constrain ${t}_{\mathrm{first}}$ and the explosion time (${t}_{\exp }$). While ${t}_{\exp }$ and ${t}_{\mathrm{first}}$ have sometimes been used interchangeably in the literature, differences arise because of a possible dark phase between the explosion and when the SN first starts to optically brighten (Piro & Nakar 2014). Finally, in Section 5, we combine all these results to constrain the radii of the exploding star and any nondegenerate companion. We summarize our results in Section 6.

2. DISCOVERY AND OBSERVATIONS

The All-sky Automated Survey for SuperNovae (ASAS-SN,26 or "Assassin"; Shappee et al. 2014) scans the entire extragalactic sky in both the northern and southern hemispheres roughly once every 2–3 nights in the V band to depths of 16.5–17.3 mag depending on lunation (Shappee et al. 2016, in preparation). We discovered a new source ($V\sim 14.9$ mag) on 2014 December 9.61 (JD ${t}_{\mathrm{disc}}=\mathrm{2,457,001.112}$) at ${\rm{R.A.}}={12}^{{\rm{h}}}{45}^{{\rm{m}}}09\buildrel{\rm{s}}\over{.} 10$, $\mathrm{decl}.=\,-00^\circ 27^{\prime} 32\buildrel{\prime\prime}\over{.} 5$ (J2000). This location is $7\buildrel{\prime\prime}\over{.} 5$ east and $10\buildrel{\prime\prime}\over{.} 5$ north (physical distance of ∼1 kpc) of the center of NGC 4666 ($z=0.0051;$ Meyer et al. 2004), which is a member of the Virgo Cluster Southern Extension (Karachentsev & Nasonova 2013). There are no pre-explosion images of this site available in the Hubble Space Telescope (HST) Archive.

We designated this new source ASASSN-14lp and released its coordinates less than 5 hr after the discovery images were taken (Holoien et al. 2014), which allowed rapid follow-up at multiple wavelenths (e.g., Foley 2014; Thorstensen et al. 2014). We will present the prediscovery and follow-up photometric observations in Section 2.1 and the spectroscopic follow-up observations in Section 2.2.

2.1. Photometric Observations

We have prediscovery images from two sources. First, ASAS-SN has been observing NGC 4666 since 2012 January and had obtained more than 125 epochs at this location both from the northern unit ("Brutus") on its fourth camera (bd) and from the southern unit ("Cassius") on its second camera (bf) before the discovery of ASASSN-14lp. Brutus obtained the discovery image on 2014 December 9.61, with previous observations obtained 1.25, 7.98, 10.98, 11.97, 15.96, and 18.96 days before discovery (see Table 1 and Figure 1 for details). Images are processed by the fully automatic ASAS-SN pipeline (B. J. Shappee et al. 2016, in preparation) using the ISIS image subtraction package (Alard & Lupton 1998; Alard 2000). The reference image and first detection difference image for both bd and bf are shown in Figure 2. ASASSN-14lp is detectable in the bf image taken 1.25 days prior to its discovery, but this image was taken under poor conditions, which prevented ASASSN-14lp from being flagged by our automated detection pipeline. We performed aperture photometry on the subtracted images using the IRAF apphot package and calibrated the results using the AAVSO Photometric All-sky Survey (APASS; Henden et al. 2015). The ASAS-SN detections and 3σ limits are presented in Table 1.

Figure 1.

Figure 1. Optical–UV photometric observations of ASASSN-14lp spanning ∼105 days from 2014 December 8 through 2015 March 23. Observations from ASAS-SN bd/bf cameras (diamonds/stars), CSPII (circles), K. Itagaki (octagons), LCOGT (squares), LT (pentagons), and Swift (hexagons) are shown with filled symbols. Downward-pointing triangles represent upper limits. Error bars are shown, but they are usually smaller than the points. Black lines show the best-fit SNooPy light-curve models. The vertical orange line shows the best-fit ${t}_{\mathrm{first}}$ (see Section 4.1); its uncertainty is smaller than the width of the line. The vertical turquoise line shows the best-fit ${t}_{\exp }$ and its uncertainty (see Section 4.2). Tick marks at the bottom of the plot denote the epochs of the spectra shown in Figure 3.

Standard image High-resolution image
Figure 2.

Figure 2. Comparison between our ASAS-SN reference image from Cassius (top left), ASAS-SN first detection difference image from Cassius acquired 2014 December 8.35 (bottom left), ASAS-SN reference image from Brutus (top middle), ASAS-SN discovery difference image from Brutus acquired 2014 December 9.61 (bottom middle), the archival Sloan Digital Sky Survey (SDSS) image (top right), and our confirmation image from the Las Cumbres Observatory Global Telescope Network (LCOGT) 1 m at the MacDonald Observatory (bottom right). All images are on the same angular scale, and the position of ASASSN-14lp is marked.

Standard image High-resolution image

Table 1.  Photometric Observations

JD Band Magnitude Telescope
(−2,450,000)      
7001.906 UVW2 17.35(0.09) Swift
7003.772 UVM2 17.83(0.10) Swift
7001.902 UVW1 16.32(0.08) Swift
7002.836 u 15.20(0.02) CSP
7001.239 g 14.91(0.03) LCOGT
6837.541 V >17.32 ASAS-SN/bf
7001.834 r 14.31(0.01) CSP
6998.279 clear >17.91 Itagaki
7001.837 i 14.51(0.01) CSP

Note. The Swift and V-band photometry is calibrated in the Vega magnitude system. The unfiltered photometry is calibrated to r band in the Vega magnitude system. The SDSS $u^{\prime} $-, $g^{\prime} $-, $r^{\prime} $-, and i'-band photometry is calibrated in the AB magnitude system.

Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.

Download table as:  DataTypeset image

The second set of prediscovery images were obtained by K. Itagaki from 2014 November 13 to 2015 January 09 with a 0.30 m f/8 reflector and an unfiltered CCD (KAF-1001E) at the Itagaki Astronomical Observatory, Japan. The KAF-1001E CCD has a quantum efficiency curve27 that reaches half its peak efficiency at ∼4100 and ∼8900 Å. We reduced the data using super-sky flats instead of twilight flats (which are unavailable). Unfortunately, only five images were available from the period prior to our discovery of ASASSN-14lp. For the image subtraction procedure we could use only one of the images for the reference image. The other four either contain the SN, are needed to constrain the progenitor (see Section 4), or were taken in poor conditions. We calibrated the unfiltered CCD photometry using the APASS r-band magnitudes because r band most closely matches the peak of the KAF-1001E detector quantum efficiency curve. We matched 36 stars between the unfiltered images and APASS r-band catalog with magnitudes ranging from 10th to 13th mag. The difference between the APASS magnitudes and instrumental magnitudes had a standard deviation of 0.094 mag (without any outlier rejection). The magnitudes and the 3σ limits are presented in Table 1.

After the detection of ASASSN-14lp, we began a multiwavelength, ground- and space-based follow-up campaign from the X-ray to the near-infrared (NIR) wavelengths. Swift(Gehrels et al. 2004) target-of-opportunity observations were performed on both the UltraViolet/Optical Telescope (UVOT; Roming et al. 2005) and the X-ray Telescope (XRT; Burrows et al. 2005). UVOT observations were obtained in V (5468 Å), B (4392 Å), U (3465 Å), UVW1 (2600 Å), UVM2 (2246 Å), and UVW2 (1928 Å; Poole et al. 2008). These observations span more than 30 days and cover the rise and fall of ASASSN-14lp. We used UVOTSOURCE to extract the source counts inside a $5\buildrel{\prime\prime}\over{.} 0$ radius and from a sky region with a radius of $\sim 40\buildrel{\prime\prime}\over{.} 0$. However, since there is significant host galaxy contamination in the optical wavelengths, we present and analyze only the ultraviolet (UV) bands in this paper. The UVOT Vega magnitudes presented in Table 1 used the most recent UVOT calibrations (Poole et al. 2008; Breeveld et al. 2010).

The XRT was operating in photon-counting mode (Hill 2004) during our observations. The data from all epochs were reduced and combined with the software tasks XRTPIPELINE and XSELECT to create a 0.3–10 keV image at the location of ASASSN-14lp. We found no point source at the location of ASASSN-14lp, and the extended emission from NGC 4666 (XMM-08 in Dietrich et al. 2006) prevents us from placing any useful limits on the X-ray flux of the SN.

Optical ground-based images were obtained with the LCOGT 1 m network of telescopes (gri), the 2 m Liverpool telescope (LT) (gri), and the 40-inch Swope telescope (ugri) as part of the Carnegie Supernova Project II (CSPII; Hsiao et al. 2013). To remove the host-galaxy flux from the LCOGT and LT data, we used the HOTPANTS28 package to subtract the SDSS archival gri images of NGC 4666 from our new images. The CSPII data were reduced in a similar manner to Contreras et al. (2010), and we obtained satisfactory subtractions using the SDSS archival ugri images of NGC 4666 as templates.

2.2. Spectroscopic Observations

We obtained an extensive sample of low- and medium-resolution optical spectra of ASASSN-14lp spanning more than 3 months between 2014 December 10 and 2015 March 16. Table 2 shows a summary of all the spectra, including date (UT+JD), telescope/instrument, wavelength range, spectral resolution, and exposure time.

Table 2.  Spectroscopic Observations

  JD $t-{t}_{\mathrm{first}}$ $t-{t}_{B\max }$   Wavelength Range Exposure
UT Date −2,400,000 (days) Telescope/Instrument (Å) (s)
2014 Dec 10.55 57,002.048 +3.2 −13.8 MDM2.4 m/ModSpec 4200–7500 240.0
2014 Dec 12.48 57,003.976 +5.1 −11.8 MDM2.4 m/ModSpec 4200–7500 600.0
2014 Dec 12.52 57,004.021 +5.1 −11.8 APO 3.5 m/DIS 3500–9600 300.0
2014 Dec 14.52 57,006.023 +7.1 −9.8 FLWO1.5 m/FAST 3500–7400 180.0
2014 Dec 15.53 57,007.025 +8.1 −8.8 FLWO1.5 m/FAST 3500–7400 300.0
2014 Dec 15.56 57,007.057 +8.2 −8.8 MDM2.4 m/ModSpec 4200–7500 480.0
2014 Dec 16.54 57,008.042 +9.2 −7.8 MDM2.4 m/ModSpec 4200–7500 600.0
2014 Dec 19.51 57,011.006 +12.1 −4.8 FLWO1.5 m/FAST 3500–7400 240.0
2014 Dec 20.54 57,012.041 +13.2 −3.8 FLWO1.5 m/FAST 3500–7400 180.0
2014 Dec 21.54 57,013.043 +14.2 −2.8 MDM2.4 m/OSMOS 4000–6850 300.0
2014 Dec 22.23 57,014.730 +15.8 −1.1 NOT2.5 m/ALFOSC 3200–9100 900.0
2015 Jan 4.31 57,026.814 +27.9 +11.0 Baade/IMACS 3600–9800 300.0
2015 Jan 5.22 57,027.724 +28.8 +11.9 NOT2.5 m/ALFOSC 3200–9100 720.0
2015 Jan 21.46 57,043.963 +45.1 +28.1 LBT/MODS 3400–10000 1800.0
2015 Jan 21.49 57,043.993 +45.1 +28.2 FLWO1.5 m/FAST 3500–7400 480.0
2015 Jan 25.29 57,047.793 +48.9 +32.0 du Pont/WFCCD 3700–9200 400.0
2015 Feb 10.50 57,063.997 +65.1 +48.2 MDM2.4 m/OSMOS 4000–6850 1200.0
2015 Feb 14.44 57,067.941 +69.1 +52.1 FLWO1.5 m/FAST 3500–7400 600.0
2015 Feb 16.50 57,070.002 +71.1 +54.2 LBT/MODS 3400–10000 900.0
2015 Feb 17.43 57,070.933 +72.0 +55.1 FLWO1.5 m/FAST 3500–7400 600.0
2015 Feb 24.25 57,077.752 +78.9 +61.9 du Pont/WFCCD 3700–9200 300.0
2015 Mar 10.38 57,091.879 +93.0 +76.1 MDM2.4 m/OSMOS 4000–6850 600.0
2015 Mar 16.34 57,097.837 +98.9 +82.0 FLWO1.5 m/FAST 3500–7400 900.0

Download table as:  ASCIITypeset image

The single-slit spectra from the Modular Spectrograph (Modspec) mounted on the MDM Observatory Hiltner 2.4 m telescope, DIS on the APO 3.5 m telescope, the Ohio State Multi-Object Spectrograph (OSMOS; Martini et al. 2011) on the MDM Observatory Hiltner 2.4 m telescope, the Inamori-Magellan Areal Camera & Spectrograph (IMACS; Dressler et al. 2011) on the Baade-Magellan 6.5 m telescope, the Andalucia Faint Object Spectrograph and Camera (ALFOSC) on the Nordic Optical 2.5 m Telescope (NOT), and the Wide Field Reimaging CCD Camera (WFCCD) on the du Pont 100 inch telescope were all reduced with standard routines in the IRAF twodspec and onedspec packages. The spectra obtained at the F. L. Whipple Observatory (FLWO) 1.5 m Tillinghast telescope using the FAST spectrograph (Fabricant et al. 1998) were reduced using a combination of standard IRAF and custom IDL procedures (Matheson et al. 2005). The spectra from the Multi-object Double Spectrographs (Pogge et al. 2010) on the dual 8.4 m Large Binocular Telescope (LBT) on Mount Graham were reduced using a custom pipeline written in IDL.29 The reductions included bias subtraction, flat-fielding, 1D spectral extraction, wavelength calibration using an arc lamp, and flux calibration using a spectroscopic standard usually taken the same night. We calibrate the spectra by extracting synthetic photometric magnitudes for each filter in the spectral range and then finding the best linear fit to the observed magnitudes interpolated to the spectral epoch. For spectra where only a single filter fit within the spectral range, the spectrum was scaled by a constant. Finally, while no filter is completely inside the OSMOS spectral range, virtually all ($\gt 98 \% $) of the g-band transmission function lies in the observed spectral range, so we calibrate to this filter.

The calibrated spectra are shown in Figure 3 with early-time spectra of SN 2009ig (Foley et al. 2012) shown underneath for reference. Strong host-galaxy Na i D absorption is visible in all ASASSN-14lp spectra, which is promising for time-viability studies to constrain the presence of circumstellar material (e.g., Patat et al. 2007; Simon et al. 2009).

Figure 3.

Figure 3. Spectroscopic observations of ASASSN-14lp spanning from 2014 December 10 to 2015 March 16. Phases relative to our estimate of ${t}_{\mathrm{first}}$ and tBmax are indicated. Details of the observations are given in Table 2. Individual spectra are calibrated using broadband photometry interpolated to the spectral epoch. For the first 15 days after first light, spectra of SN 2009ig interpolated to the same epoch are shown in gray (Foley et al. 2012). Vertical dashed gray lines mark strong telluric features. The vertical apricot line marks host-galaxy Na i D absorption.

Standard image High-resolution image

In Figure 4, we show the pre-maximum evolution of the spectroscopic region near the Si ii λ6355 feature. For comparison, we also show the spectra of the two others, SN 2009ig (Foley et al. 2012) and SN 2011fe (Pereira et al. 2013), that were both observed within 1 day of first light and had early-time spectroscopic observations. These spectra were interpolated in time to the same epoch since ${t}_{\mathrm{first}}$ as ASASSN-14lp. While detailed modeling of these spectra is beyond the scope of this paper, it is interesting to note that a high-velocity Si ii component (∼18,000 $\mathrm{km}\,{{\rm{s}}}^{-1}$) is clearly visible in the ASASSN-14lp spectra at early times. The Si ii feature then becomes dominated by a lower-velocity component (∼13,000 $\mathrm{km}\,{{\rm{s}}}^{-1}$) starting 5 days after ${t}_{\mathrm{first}}$. This behavior is qualitatively similar to SN 2009ig but was not observed in SN 2011fe (Piro & Nakar 2014).

Figure 4.

Figure 4. Pre-maximum evolution of the region near the Si ii λ6355 feature. Phases from our derived ${t}_{\mathrm{first}}$ from Section 4.1 are indicated. The Foley et al. (2012) spectra of SN 2009ig (solid gray) and the Pereira et al. (2013) spectra of SN 2011fe (dotted black), interpolated to the same epochs as ASASSN-14lp, are also shown. Like ASASSN-14lp, SN 2009ig also had a high-velocity Si ii component that gives way to a lower-velocity component starting roughly 5 days after first light.

Standard image High-resolution image

3. LIGHT-CURVE FIT AND DISTANCE TO NGC 4666

We used the CSPII light-curve fitting code SuperNovae in Object Oriented Python (SNooPy; Burns et al. 2011) to fit the light curve of ASASSN-14lp and to derive a distance to NGC 4666. We employed SNooPy's systematics function, the spectral template from Hsiao et al. (2007), and the "Best-observed" calibration (Fit 4 from Table 9 of Folatelli et al. 2010). We adopted a Galactic extinction of $E{(B-V)}_{\mathrm{MW}}=0.021$ mag (Schlegel et al. 1998; Schlafly & Finkbeiner 2011) and an ${R}_{V}=3.1$ reddening law (Cardelli et al. 1989).

The results of the SNooPy fit are shown in Table 3. The light curve of ASASSN-14lp is noteworthy because it is broad (${\rm{\Delta }}{m}_{15}(B)=0.80\pm 0.05$) but not outside the range for normal supernovae (${\rm{\Delta }}{m}_{15}(B)\sim 0.79\mbox{--}1.62;$ see Figure 4 in Stritzinger et al. 2011). Employing the relation in Mazzali et al. (2007), the estimate of ${\rm{\Delta }}{m}_{15}(B)$ yields a 56Ni mass of 0.806 ± 0.013 ${M}_{\odot }$, which is high but not unique among SNe Ia (e.g., Figure 1 in Piro et al. 2014). The time of B-band maximum is JD ${t}_{B\max }=2457015.82\pm 0.03$, 14.728 days after discovery. We see that ASASSN-14lp suffered moderate host galaxy extinction, with $E{(B-V)}_{\mathrm{host}}=0.33\pm 0.06$. Finally, the derived distance modulus to NGC 4666 is $\mu =30.8\pm 0.2$. This distance modulus implies a physical distance to NGC 4666 of 14.7 ± 1.5 Mpc.

Table 3.  Derived ASASSN-14lp Properties

tBmax $E{(B-V)}_{\mathrm{host}}$ ${\rm{\Delta }}{m}_{15}$ μ d
(−2,450,000) (mag) (mag) (mag) (Mpc)
7015.82 ± 0.03 0.33 ± 0.06 0.80 ± 0.05 30.8 ± 0.2 14.7 ± 1.5

Note. The quoted errors include systematic errors as reported by SNooPy's systematics function (Burns et al. 2011).

Download table as:  ASCIITypeset image

4. FIRST LIGHT AND EXPLOSION TIME CONSTRAINTS

In this section we present constraints on ${t}_{\mathrm{first}}$ for ASASSN-14lp from fits to the early-time light curve (Section 4.1), and from the early-time expansion velocity we constrain ${t}_{\exp }$ (Section 4.2).

4.1. Early-time Light-curve Fit

To determine ${t}_{\mathrm{first}}$, we first attempted to model the early-time ASASSN-14lp light curve with an expanding fireball model where flux is proportional to ${(t-{t}_{\mathrm{first}})}^{2}$ . However, this led to unsatisfactory fits. We then added more freedom to our model, fitting an independent power-law index (${\alpha }_{i}$) for each band $\propto {(t-{t}_{\mathrm{first}})}^{{\alpha }_{i}}$, but forcing them all to have the same ${t}_{\mathrm{first}}$. We used the emcee package (Foreman-Mackey et al. 2013), a Python-based implementation of the affine-invariant ensemble sampler for Markov chain Monte Carlo, to perform the fit to each light curve. We excluded the Swift UVM2 band from this analysis because ASASSN-14lp was not clearly detected within 5 days of discovery. The best-fit power laws for each filter and their corresponding 1σ uncertainties are shown in Figure 5. We find that JD ${t}_{\mathrm{first}}={2456998.88}_{-0.10}^{+0.11}$, implying that we discovered ASASSN-14lp just ${2.23}_{-0.10}^{+0.11}$ days after ${t}_{\mathrm{first}}$ and that ${t}_{\mathrm{rise}}={t}_{B\max }-{t}_{\mathrm{first}}={16.94}_{-0.10}^{+0.11}$ days. The best-fit α for each band are presented in Table 4. The g, V, r, clear, and i filters have best-fit $\alpha \sim 1.50\mbox{--}1.75$, but the u filter requires a steeper rise.

Figure 5.

Figure 5. Early-time light curves from 3 days before to 6 days after discovery with power-law fits. The colors and symbols are the same as in Figure 1. Colored solid lines are the best-fit power laws, while dashed black lines show each fit's 1σ uncertainties. The vertical orange line shows the best-fit ${t}_{\mathrm{first}}$ and its uncertainty (see Section 4.1). Top panel: light curves normalized by their flux at 5 days after discovery with an added constant. Bottom panels: residuals from the best-fit power laws for each individual filter.

Standard image High-resolution image

Table 4.  Fit Light Curve Parameters

Band ${t}_{\max }$ ${m}_{\max }$ α
  (−2,450,000) (mag)  
UVW2 7013.12 ± 0.77 15.61 ± 0.03 ${1.45}_{-0.43}^{+0.46}$
UVM2 7016.52 ± 1.03 17.36 ± 0.08
u 7015.14 ± 0.21 12.84 ± 0.03 ${2.38}_{-0.05}^{+0.05}$
g 7017.57 ± 0.20 12.15 ± 0.02 ${1.70}_{-0.04}^{+0.04}$
V 7018.11 ± 0.19 11.94 ± 0.01 ${1.50}_{-0.05}^{+0.05}$
r 7018.01 ± 0.11 11.94 ± 0.01 ${1.57}_{-0.04}^{+0.04}$
clear 7016.67 ± 0.21 12.04 ± 0.02 ${1.75}_{-0.16}^{+0.19}$
i 7015.02 ± 0.22 12.50 ± 0.01 ${1.58}_{-0.04}^{+0.04}$

Download table as:  ASCIITypeset image

The first photometric observations around 1 day after ${t}_{\mathrm{first}}$ from both ASAS-SN and K. Itagaki are brighter than the simple power-law fits by about 3σ and 1σ, respectively. These are weak detections, but they perhaps hint that a broken power-law fit would be a better description at these very early times, similar to SN 2013dy (Zheng et al. 2013). However, with only one photometric point in each band during the first 2 days after ${t}_{\mathrm{first}}$, such a fit is not well constrained. Both of these points are fully accounted for when placing constraints on the progenitor system in Section 5.

In Figure 6, we compare ${t}_{\mathrm{rise}}$, ${\rm{\Delta }}{m}_{15}(B)$, and α for ASASSN-14lp against all available prior measurements: SN 2009ig (Foley et al. 2012), SN 2012cg (Silverman et al. 2012), SN 2014J (Marion et al. 2015; Zheng et al. 2014; Siverd et al. 2015), and the sample (including SN 2011fe) from Firth et al. (2015). We use Equation (5) of Conley et al. (2008) to translate the stretch reported in Table 4 of Firth et al. (2015) into ${\rm{\Delta }}{m}_{15}(B)$. We compare our r-band α to those reported by Firth et al. (2015). ASASSN-14lp is on the edge of, but consistent with, these distributions. In summary, ASASSN-14lp is normal, but it has a post-maximum light curve that is broad, an overall rise time that is short, and an early-time rise that is slow.

Figure 6.

Figure 6. Comparison between ${t}_{\mathrm{rise}}$, ${\rm{\Delta }}{m}_{15}(B)$, and α for ASASSN-14lp (red star), SN 2009ig (blue diamond; Foley et al. 2012), SN 2012cg (blue triangle; Silverman et al. 2012), and the Firth et al. (2015) collection of SNe (gray circles), which includes SN 2011fe (blue circle). Error bars for ASASSN-14lp, SN 2009ig, SN 2012cg, and SN 2011fe are smaller than the points.

Standard image High-resolution image

4.2. Early-time Expansion Velocity Fit

Following Piro & Nakar (2014), we use the early-time spectroscopic observations to measure line velocities, and then by assuming a theoretically motivated, time-dependent velocity evolution model, we estimate ${t}_{\exp }$. We then compare ${t}_{\exp }$, determined by fitting these early-time spectroscopic observations, and ${t}_{\mathrm{first}}$, determined by fitting the early-time light curve, to constrain the existence of a possible dark phase.

First, we measured the velocities of the Ca ii H and K and Si ii λ6355 absorption features for the spectra presented in Table 2 taken within 15 days of ${t}_{\mathrm{first}}$. We used the wavelength of the deepest absorption to determine the velocity. Similar to Piro & Nakar (2014), we assume an error of 500 km s−1 as a rough estimate of the typical measurement error (Parrent et al. 2012). As previously mentioned, the spectra in Figure 4 show that there are two components to the Si ii λ6355 line: a high-velocity component that dominates the first ∼5 days after ${t}_{\mathrm{first}}$ and a low-velocity component that dominates afterward. Measuring the velocity of both Si ii λ6355 components at all epochs would require detailed modeling of the spectra, which is beyond the scope of this paper. Instead, we simply measured the deepest absorption of the dominant Si ii λ6355 component in each spectrum.

We then fit these line velocities by assuming that they have particular power-law evolution in time and varying the explosion time. The expected velocity evolution for the spectral lines of very young SNe Ia is a power law with time ($v\propto {t}^{-0.22}$) if the lines are generated at a constant opacity (Piro & Nakar 2014). The bottom panel of Figure 7 shows the Ca ii H and K and Si ii λ6355 velocities and their best fit assuming $v\propto {t}^{-0.22}$. Piro & Nakar (2014) use the spread in the best-fit explosion time between $v\propto {t}^{-0.20}$ and $v\propto {t}^{-0.24}$ as an estimate of the uncertainty. The top panel of Figure 7 shows the goodness of fit, as a function of explosion time, for these three power laws: $v\propto {t}^{-0.20}$, $v\propto {t}^{-0.22}$, and $v\propto {t}^{-0.24}$. This leads us to JD ${t}_{\exp }={2456996.93}_{-0.84}^{+0.81}$, which is ${1.93}_{-0.85}^{+0.81}$ days earlier than ${t}_{\mathrm{first}}$.

Figure 7.

Figure 7. Constraints on ${t}_{\exp }$ from fitting the velocity evolution of the Ca ii H and K and Si ii λ6355 following Piro & Nakar (2014). Top panel: the dotted, solid, and dashed lines show the ${\chi }^{2}$ for fit power laws with indices −0.20, −0.22, and −0.24, respectively, as a function of assumed explosion date. The orange line and region show the best-fit ${t}_{\mathrm{first}}$ and its 1σ error bars from Section 4.1. Bottom panel: observed low-velocity Si ii λ6355 (squares), high-velocity Si ii λ6355 (circles), and Ca ii H and K (diamonds) with best-fit power-law velocity with $v\propto {t}^{-0.22}$.

Standard image High-resolution image

The offset between ${t}_{\exp }$ and ${t}_{\mathrm{first}}$ is similar to that found by Piro & Nakar (2014) for SN 2009ig, which is spectroscopically (see Section 2.2) and photometrically similar (SN 2009ig also had a broad light curve with ${\rm{\Delta }}{m}_{15}(B)=0.89;$ Foley et al. 2012) to ASASSN-14lp. There are three possible explanations for this difference: (1) The modeling method presented in Piro & Nakar (2014) does not accurately represent any SNe Ia, perhaps because of the assumption that the opacity near the photosphere is roughly constant or the velocity profile of the material is different. Although possible, one would still need to explain why this method seemed to work well for SN 2011fe and SN 2012cg and not for SN 2009ig and ASASSN-14lp. (2) The modeling method presented in Piro & Nakar (2014) is fine in general, but something about the SN 2009ig and ASASSN-14lp explosions breaks its assumptions. Piro & Nakar (2014) suggest that the discrepancy might be due to the explosions being asymmetric, with higher-velocity material coming toward the observer (Maeda et al. 2010). (3) There is a dark time between when the explosion occurs and when optical brightening occurs for some SNe (Piro & Nakar 2013), and fitting a power law to the early-time light curve and extrapolating this to earlier times does not always accurately estimate explosion times. Distinguishing between these three possibilities is beyond the scope of this paper, but it is an important question for early-time studies of SNe Ia because the presence of a dark time between explosion and first light can drastically alter inferences about the structure of the progenitor. To help investigate the underlying mechanism that causes these features, it will be useful to explore whether a larger sample of early-time light curves can be put into groups that share basic properties, as we have seen for SN 2009ig and ASASSN-14lp or SN 2011fe and SN 2012cg.

5. PROGENITOR SYSTEM CONSTRAINTS

To place constraints on the progenitor system of ASASSN-14lp, we first create absolute magnitude light curves from the two filters, ASAS-SN V and K. Itagaki's clear, with very early (within 2 days of ${t}_{\mathrm{first}}$) detections or upper limits. In this section we treat the clear filter as r, based on the calibrations in Section 2.1. One could use the distance modulus and reddening determined in Section 3 to compute absolute magnitude light curves from our observed light curves; however, this would be circular because SNooPy assumed an absolute magnitude of ASASSN-14lp to derive its distance. Instead, we scale our observed V- and r-band light curves to the expected absolute maximum given the observed value of ${\rm{\Delta }}{m}_{15}(B)$ and the "Best-observed" calibration of Folatelli et al. (2010) for these filters (their Fits 7 and 10 from Table 9, respectively). To do this, we determine the peak magnitude (${m}_{\max }$) in each filter by fitting a quadratic to each filter's light curve over ±10 days around tBmax. The time of maximum (${t}_{\max }$) and ${m}_{\max }$ for each band are presented in Table 4. We add the errors induced by the intrinsic scatter of SNe Ia (${\sigma }_{\mathrm{SN}}$), the intercept of the Phillips relation (${M}_{{\rm{X}}}$), the slope of the Phillips relation (${b}_{{\rm{X}}}$), and the estimate of ${m}_{\max }$ in quadrature with the photometric errors. This means that the photometric points in the right panel of Figure 8 have correlated errors. For the upper limits, we take a conservative approach and add the quadrature of these errors to the upper limit. We use the early-time absolute magnitude light curves of ASASSN-14lp to constrain its progenitor system in two ways. In Section 5.1, we constrain the radius of a possible nondegenerate companion by placing limits on its interaction with the SN ejecta. In Section 5.2, we constrain the radius of the exploding star itself by placing limits on flux from its shock-heated envelope.

Figure 8.

Figure 8. Constraints on the companion radius (${R}_{{\rm{c}}}$) assuming Roche-lobe overflow and a 1 ${M}_{\odot }$ companion. Left panel: maximum allowed ${R}_{{\rm{c}}}$ as a function of explosion time. The constraints from the ASAS-SN V-band and K. Itagaki clear light curves are represented by the dashed and solid lines, respectively. The first photometric point detecting ASASSN-14lp (${t}_{\det }$) and the best-fit ${t}_{\mathrm{first}}$ and ${t}_{\exp }$ are indicated. The purple, blue, and red points indicate the ${R}_{{\rm{c}}}$ and explosion time of the companion collision light curves shown in the right panel plotted in the same color. Right panel: the early-time, clear light curve of ASASSN-14lp as compared to theoretical light curves for three choices of the stellar radius.

Standard image High-resolution image

5.1. Constraints on a Nondegenerate Companion's Radius

If the progenitor of an SN Ia is a WD accreting from a nondegenerate companion, then the SN ejecta will interact with the companion and potentially produce an observable signature in the rising light curve at early times. Such a signature is dependent on the viewing angle, with the strongest effect occurring when the companion lies along the line of sight between the observer and the SN. At a fixed viewing angle, this emission scales proportionally with ${R}_{{\rm{c}}}$, which allows early-time detections or upper limits to constrain the properties of the companion.

We apply the analytic models provided by Kasen (2010) to place constraints on the radius of a potential companion to ASASSN-14lp. Since the explosion time is uncertain, we explore a variety of explosion times and their corresponding constraints on ${R}_{{\rm{c}}}$. The general procedure is to choose an explosion time and then find the maximum radius companion that is consistent with the early-time data assuming that the secondary filled its Roche lobe.

This introduces a weak dependence on the mass ratio of the binary,30 and we simply assume primary and companion masses of 1.4 and $1\,{M}_{\odot }$, respectively. The results are summarized in Figure 8, which demonstrates that the constraints on ${R}_{{\rm{c}}}$ are strongly dependent on the explosion time. If the explosion time is near the ${t}_{\mathrm{first}}$ determined in Section 4.1, a maximum ${R}_{{\rm{c}}}$ of $\approx 31\,{R}_{\odot }$ is found. On the other hand, if the explosion time is close to ${t}_{\exp }$, then ${R}_{{\rm{c}}}\lesssim 0.35\,{R}_{\odot }$. These estimates assume a viewing angle of 15°.

Although there are a wide range of potential progenitor scenarios, binary systems with nondegenerate companions can be roughly broken into three main cases: (1) Systems with a red giant companion, similar to the observed binaries RS Oph and T CrB (e.g., Hachisu & Kato 2001), with a large radius of $\gt 100\,{R}_{\odot }$; this scenario is difficult to reconcile with ASASSN-14lp for most viewing angles. (2) Systems with a helium star companion, like the helium-nova system V445 Pup (e.g., Kato et al. 2008), with radii of ∼0.3–6 ${R}_{\odot }$; even if the explosion time is earlier, this scenario is consistent with our constraints. (3) Systems with a main-sequence or subgiant companion, similar to U Sco (e.g., Thoroughgood et al. 2001), with radii of ∼0.4–4 ${R}_{\odot }$; a later explosion time does not provide meaningful constraints on this scenario, but a high-mass main-sequence companion is possible for early explosion times only for unlucky viewing angles. Overall, we rule out a red giant companion for all but the most unfavorable viewing angles, but a helium star, a main sequence, and a subgiant companion are all consistent with our observations.

5.2. Constraints on the Primary's Radius

Assuming that there is no appreciable emission from the interaction of the SN ejecta and a companion, the first light expected from an SN Ia is the cooling of the expanding, shock-heated primary star (Piro & Nakar 2013). This emission will, however, be overtaken by radioactive nickel heating after the first few days. The emission from the cooling of the shock-heated primary is roughly proportional to R* (Piro et al. 2010). This fact was utilized by Bloom et al. (2011) to put tight constraints on the radius of the star that exploded in SN 2011fe. We repeat this exercise for ASASSN-14lp and again consider a variety of explosion times, which we summarize in Figure 9. If the explosion time is near the ${t}_{\mathrm{first}}$ determined in Section 4.1, then a maximum R* of $\approx 6\,{R}_{\odot }$ is found. On the other hand, if the explosion time is close to ${t}_{\exp }$, then ${R}_{* }\lesssim 0.6\,{R}_{\odot }$. Neither of these constraints is surprising given the general agreement that the progenitor of a SN Ia is a WD of much smaller size. Nevertheless, since this emission is isotropic, these constraints are firmer than the viewing-angle-dependent constraints from Section 5.1 on interaction with the companion.

Figure 9.

Figure 9. Constraints on the primary star's radius (R*) from a shock-heated envelope as a function of explosion time. The constraints from the ASAS-SN V-band and K. Itagaki clear light curves are represented by the dashed and solid lines, respectively. The first photometric point detecting ASASSN-14lp (${t}_{\det }$) and the best-fit ${t}_{\mathrm{first}}$ and ${t}_{\exp }$ are indicated.

Standard image High-resolution image

6. DISCUSSION AND SUMMARY

We discovered ASASSN-14lp on 2014 December 9.61, just ${2.23}_{-0.10}^{+0.11}$ days after first light. We announced it to the community less than 5 hr later, and we promptly classified it as an SN Ia (Thorstensen et al. 2014). Here we added prediscovery photometry, with our first detection less than a day after ${t}_{\mathrm{first}}$. We also presented ultraviolet through NIR photometric and spectroscopic follow-up data covering the rise and fall of ASASSN-14lp for more than 100 days. ASASSN-14lp had a broad light curve (${\rm{\Delta }}{m}_{15}(B)=0.80\pm 0.05$), a B-band maximum at JD 2,457,015.82 ± 0.03, a rise time of ${16.94}_{-0.10}^{+0.11}$ days, and moderate host-galaxy extinction ($E{(B-V)}_{\mathrm{host}}=0.33\pm 0.06$). Using the light curve of the SN, we determine the distance modulus to the host galaxy, NGC 4666, to be $\mu =30.8\pm 0.2$, or a distance of 14.7 ± 1.5 Mpc.

In a future work, the RV of the host galaxy reddening will be measured using the high-quality optical–NIR CSPII data once ASASSN-14lp has faded sufficiently to allow template images to be obtained. Such measurements are of interest because one of the largest systematic uncertainties limiting the precision of SNe Ia as a cosmological distance indicator is the uncertainty in the host-galaxy reddening law (Mandel et al. 2011).

While there is a Tully–Fisher estimate of the distance (e.g., Springob et al. 2007; Sorce et al. 2014), a better measurement of the distance to NGC 4666 with HST is needed to add ASASSN-14lp to the surprisingly small calibrating sample (currently nine; Riess et al. 2011 and SN 2011fe) of SNe Ia. However, NGC 4666 is highly inclined, so crowding and extinction will make a Cepheid distance determination difficult. Fortunately, the inclination of NGC 4666 is actually an advantage for the tip of the red giant branch method of determining distance because it is easier to probe the halo population with reduced contamination from disk stars (Madore & Freedman 1998). ASASSN-14lp should be added to the ongoing "Carnegie Hubble Program to Measure H0 to 3% Using Population II" project (CHP II; Freedman 2014) or a similar study.

We used our early-time light curve to constrain JD ${t}_{\mathrm{first}}={2456998.88}_{-0.10}^{+0.11}$, implying that ${t}_{\mathrm{rise}}={16.94}_{-0.10}^{+0.11}$ days and that we discovered ASASSN-14lp at day ${2.23}_{-0.10}^{+0.11}$. We then used the early-time spectroscopic data to determine JD ${t}_{\exp }={2456996.93}_{-0.84}^{+0.81}$ based on the methods of Piro & Nakar (2014). This is ${1.93}_{-0.85}^{+0.81}$ days before the ${t}_{\mathrm{first}}$ determined from the light-curve fit, which is similar to the behavior of SN 2009ig (Piro & Nakar 2014). This difference might indicate that there was a significant dark time between the explosion and when ASASSN-14lp began optically brightening or that the explosion of ASASSN-14lp was asymmetrical, but a more detailed study is needed.

Finally, we used our early-time photometric and spectroscopic data along with our derived light-curve properties to place constraints on the progenitor system. First, we constrained ${R}_{* }\,\lt $ 0.5–6 ${R}_{\odot }$ depending on the explosion time (independent of viewing angle). Second, for a viewing angle of 15°, we rule out the presence of a nondegenerate companion with a radius of 0.34–11 ${R}_{\odot }$ or larger. This range originates from the choice of ${t}_{\exp }$ or ${t}_{\mathrm{first}}$ for the explosion time. While this is a favorable viewing angle, the strong constraints make it difficult to reconcile ASASSN-14lp with the presence of a red giant companion.

Though ASAS-SN is not designed to find bright SNe Ia as early as possible, the discovery of ASASSN-14lp within ∼2 days after first light is quite encouraging. Furthermore, during the previous ASAS-SN epoch, within 1 day after first light, ASASSN-14lp was 1.5 mag above the limiting magnitude of ASAS-SN. Even though we failed to trigger on the transient because of poor weather conditions, this observation demonstrates that sub-day discoveries of bright SNe Ia are possible.

It is the nearby, bright SNe Ia that have recently allowed us to further our understanding of their progenitor systems. This is because many of the most constraining measurements can only be performed for these objects: HST pre-explosion imaging (e.g., Li et al. 2011; Graur & Maoz 2012; Graur et al. 2014; Kelly et al. 2014), constraints on stripped material from a companion (e.g., Mattila et al. 2005; Leonard 2007; Shappee et al. 2013a), X-ray and radio constraints on the presence of circumstellar material (e.g., Chomiuk et al. 2012a, 2012b; Horesh et al. 2012; Pérez-Torres et al. 2014), and constraints on circumstellar material from variable narrow absorption lines (e.g., Patat et al. 2007; Simon et al. 2009). Moreover, the SNe where multiple constraints have been possible (e.g., SN 2011fe, SN 2012cg, and SN 2014J) have proven to be the most enlightening.

This leads us to the question, "What is the earliest that bright SNe Ia can be found with ASAS-SN or other deeper surveys?" To address this question, we use the early-time power-law fits of ASASSN-14lp from Table 4 to determine that surveys with limiting magnitudes of V ∼ 17 mag and r ∼ 20.5 mag would have been capable of detecting ASASSN-14lp 9.5 and 1.9 hr after first light, respectively. In this case, cadence becomes the limiting factor; unless it is more rapid than ∼10 hr, deeper surveys do not have an advantage over ASAS-SN. In addition, ASAS-SN has greater sky coverage.

The primary science aim of the ASAS-SN survey is to systematically find nearby supernovae in the entire extragalactic sky, and, despite its small (14 cm) lenses, ASAS-SN has become the dominant discover of bright ($V\gt 17$ mag) SNe thus far in 2015 (B. J. Shappee et al. 2016, in preparation). Furthermore, the Cassius station in Chile is currently (July 2015) being upgraded from two to four telescopes, to match the Brutus station in Hawaii. Long term, the goal is to expand ASAS-SN to at least two more stations, which will minimize temporal gaps due to detrimental weather or (rare) technical problems at any given ASAS-SN station.

The authors thank Jennifer van Saders, Curtis McCully, Peter Brown, Saurabh Jha, and Louisa Diodato for discussions and encouragement, Toru Yusa for aid in translation, Joe Antognini for reading, and Ryan Foley for providing us with the spectra of SN 2009ig. We thank Neil Gehrels for approving our ToO requests and the Swift science operation team for performing the observations. We thank LCOGT and its staff for their continued support of ASAS-SN. We also thank the anonymous referee for his/her helpful comments and suggestions improving this manuscript.

Development of ASAS-SN has been supported by NSF grant AST-0908816 and CCAPP at the Ohio State University. C.S.K. and K.Z.S. are supported by NSF grants AST-1515876 and AST-1515927. ASAS-SN is supported in part by the Mt. Cuba Astronomical Foundation. Operations of the Cassius ASAS-SN station are partially funded from project IC120009 Millennium Institute of Astrophysics (MAS) of the Millennium Science Initiative, Chilean Ministry of Economy.

B.S. is supported by NASA through Hubble Fellowship grant HF-51348.001 awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS 5-26555. T.W.-S.H. is supported by the DOE Computational Science Graduate Fellowship, grant number DE-FG02-97ER25308. Support for J.L.P. is in part provided by FONDECYT through grant 1151445 and by the Ministry of Economy, Development, and Tourism's Millennium Science Initiative through grant IC120009, awarded to the Millennium Institute of Astrophysics, MAS. J.F.B. is supported by NSF grant PHY-1404311. S.D. is supported by "the Strategic Priority Research Program—The Emergence of Cosmological Structures" of the Chinese Academy of Sciences (grant No. XDB09000000). E.H. and M.S. are supported by the Danish Agency for Science and Technology and Innovation realized through a Sapere Aude Level 2 grant. P.R.W. is supported by the Laboratory Directed Research and Development program at LANL.

CSPII is supported by the NSF under grants AST-0306969, AST-0607438, and AST-1008343.

This paper uses data products produced by the OIR Telescope Data Center, supported by the Smithsonian Astrophysical Observatory.

The Liverpool Telescope is operated on the island of La Palma by Liverpool John Moores University in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias with financial support from the UK Science and Technology Facilities Council.

The Nordic Optical Telescope is operated by the Nordic Optical Telescope Scientific Association at the Observatorio del Roque de los Muchachos, La Palma, Spain, of the Instituto de Astrofisica de Canaries.

This paper is based in part on observations obtained through the CNTAC proposal CN2015A-119.

This paper used data obtained with the MODS spectrographs built with funding from NSF grant AST-9987045 and the NSF Telescope System Instrumentation Program (TSIP), with additional funds from the Ohio Board of Regents and the Ohio State University Office of Research. The LBT is an international collaboration among institutions in the United States, Italy, and Germany.

This paper includes data gathered with the 6.5 m Magellan Telescopes located at Las Campanas Observatory, Chile.

IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation.

Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, the U.S. Department of Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education Funding Council for England. The SDSS Web site is http://www.sdss.org/.

This research has made use of data from the AAVSO Photometric All Sky Survey, whose funding has been provided by the Robert Martin Ayers Sciences Fund.

This research has made use of the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.

This research has also made use of the XRT Data Analysis Software (XRTDAS) developed under the responsibility of the ASI Science center (ASDC), Italy.

Footnotes

Please wait… references are loading.
10.3847/0004-637X/826/2/144