Palomar Transient Factory Survey Research Articles

Near-infrared Supernova Ia Distances: Host Galaxy Extinction and Mass-step Corrections Revisited

Abstract We present optical and near-infrared (NIR, Y-, J-, H-band) observations of 42 Type Ia supernovae (SNe Ia) discovered by the untargeted intermediate Palomar Transient Factory survey. This new data set covers a broad range of redshifts and host galaxy stellar masses, compared to previous SN Ia efforts in the NIR. We construct a sample, using also literature data at optical and NIR wavelengths, to examine claimed correlations between the host stellar masses and the Hubble diagram residuals. The SN magnitudes are corrected for host galaxy extinction using either a global total-to-selective extinction ratio, R V = 2.0, for all SNe, or a best-fit R V for each SN individually. Unlike previous studies that were based on a narrower range in host stellar mass, we do not find evidence for a “mass step,” between the color- and stretch-corrected peak J and H magnitudes for galaxies below and above log ( M * / M ⊙ ) = 10 . However, the mass step remains significant (3σ) at optical wavelengths (g, r, i) when using a global R V , but vanishes when each SN is corrected using their individual best-fit R V . Our study confirms the benefits of the NIR SN Ia distance estimates, as these are largely exempted from the empirical corrections dominating the systematic uncertainties in the optical.

From core collapse to superluminous: the rates of massive stellar explosions from the Palomar Transient Factory

ABSTRACT We present measurements of the local core-collapse supernova (CCSN) rate using SN discoveries from the Palomar Transient Factory (PTF). We use a Monte Carlo simulation of hundreds of millions of SN light-curve realizations coupled with the detailed PTF survey detection efficiencies to forward model the SN rates in PTF. Using a sample of 86 CCSNe, including 26 stripped-envelope SNe (SESNe), we show that the overall CCSN volumetric rate is $r^\mathrm{CC}_v=9.10_{-1.27}^{+1.56}\times 10^{-5}\, \text{SNe yr}^{-1}\, \text{Mpc}^{-3}\, h_{70}^{3}$ at 〈z〉 = 0.028, and the SESN volumetric rate is $r^\mathrm{SE}_v=2.41_{-0.64}^{+0.81}\times 10^{-5}\, \text{SNe yr}^{-1}\, \text{Mpc}^{-3}\, h_{70}^{3}$. We further measure a volumetric rate for hydrogen-free superluminous SNe (SLSNe-I) using eight events at z ≤ 0.2 of $r^\mathrm{SLSN-I}_v=35_{-13}^{+25}\, \text{SNe yr}^{-1}\text{Gpc}^{-3}\, h_{70}^{3}$, which represents the most precise SLSN-I rate measurement to date. Using a simple cosmic star formation history to adjust these volumetric rate measurements to the same redshift, we measure a local ratio of SLSN-I to SESN of ${\sim}1/810^{+1500}_{-94}$, and of SLSN-I to all CCSN types of ${\sim}1/3500^{+2800}_{-720}$. However, using host galaxy stellar mass as a proxy for metallicity, we also show that this ratio is strongly metallicity dependent: in low-mass (logM* < 9.5 M⊙) galaxies, which are the only environments that host SLSN-I in our sample, we measure an SLSN-I to SESN fraction of $1/300^{+380}_{-170}$ and $1/1700^{+1800}_{-720}$ for all CCSN. We further investigate the SN rates a function of host galaxy stellar mass, and show that the specific rates of all CCSNe decrease with increasing stellar mass.

From core collapse to superluminous: The rates of massive stellar explosions from the Palomar Transient Factory

We present measurements of the local core collapse supernova (SN) rate using SN discoveries from the Palomar Transient Factory (PTF). We use a Monte Carlo simulation of hundreds of millions of SN light curve realizations coupled with the detailed PTF survey detection efficiencies to forward-model the SN rates in PTF. Using a sample of 86 core collapse SNe, including 26 stripped-envelope SNe (SESNe), we show that the overall core collapse SN volumetric rate is $r^\mathrm{CC}_v=9.10_{-1.27}^{+1.56}\times10^{-5}\,\text{SNe yr}^{-1}\,\text{Mpc}^{-3}\, h_{70}^{3}$ at $ \langle z \rangle = 0.028$, and the SESN volumetric rate is $r^\mathrm{SE}_v=2.41_{-0.64}^{+0.81}\times10^{-5}\, \text{SNe yr}^{-1}\,\text{Mpc}^{-3}\, h_{70}^{3}$. We further measure a volumetric rate for hydrogen-free superluminous SNe (SLSNe-I) using 8 events at $z{\le}0.2$ of $r^\mathrm{SLSN-I}_v=35_{-13}^{+25}\, \text{SNe yr}^{-1}\text{Gpc}^{-3}\, h_{70}^{3}$, which represents the most precise SLSN-I rate measurement to date. Using a simple cosmic star-formation history to adjust these volumetric rate measurements to the same redshift, we measure a local ratio of SLSN-I to SESN of $\sim1/810^{+1500}_{-94}$, and of SLSN-I to all CCSN types of $\sim 1/3500^{+2800}_{-720}$. However, using host galaxy stellar mass as a proxy for metallicity, we also show that this ratio is strongly metallicity dependent: in low-mass ($\mathrm{log} M_{*} < 9.5 \mathrm{M}_\odot$) galaxies, which are the only environments that host SLSN-I in our sample, we measure a SLSN-I to SESN fraction of $1/300^{+380}_{-170}$ and $1/1700^{+1800}_{-720}$ for all CCSN. We further investigate the SN rates a function of host galaxy stellar mass and show that the specific rates of all core collapse SNe decrease with increasing stellar mass.

A Catalog of Short Period Spectroscopic and Eclipsing Binaries Identified from the LAMOST and PTF Surveys

Binaries play key roles in determining stellar parameters and exploring stellar evolution models. We build a catalog of 88 eclipsing binaries with spectroscopic information, taking advantage of observations from both the Large Sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST) and the Palomar Transient Factory (PTF) surveys. A software pipeline is constructed to identify binary candidates by examining their light curves. The orbital periods of binaries are derived from the Lomb-Scargle method. The key distinguishing features of eclipsing binaries are recognized by a new filter \textit{Flat Test}. We classify the eclipsing binaries by applying Fourier analysis on the light curves. Among all the binary stars, 13 binaries are identified as eclipsing binaries for the first time. The catalog contains information: position, primary eclipsing magnitude and time, eclipsing depth, the number of photometry and radial velocity observations, largest radial velocity difference, binary type, the effective temperature of observable star $T_{\rm eff}$, and surface gravity of observable star log \emph{g}. The false-positive probability is calculated by using both a Monte Carlo simulation and real data from the SDSS Stripe 82 Standard Catalog. The binaries in the catalog are mostly with a period of less than one day. The period distribution shows a 0.22-day cut-off which is consistent with the low probability of an eclipsing binary rotating with such a period.

On the Granulation and Irregular Variation of Red Supergiants

Abstract The mechanisms and characteristics of the irregular variations of red supergiants (RSGs) are studied based on the RSG samples in the Small Magellanic Cloud (SMC), Large Magellanic Cloud (LMC), and M31. With the time-series data from the All-Sky Automated Survey for SuperNovae and Intermediate Palomar Transient Factory survey, we use the continuous time autoregressive moving average model to estimate the variability features of the light curves and their power spectral density (PSD). The characteristic evolution timescale and amplitude of granulations are further derived from fitting the posterior PSD with the COR function, which is a Harvey-like granulation model. The consistency of theoretical predictions and results is checked to verify the correctness of the assumption that granulations on RSGs contribute to irregular variation. The relations between granulation and stellar parameters are obtained and compared with the results of red giant branch (RGB) stars and Betelgeuse. It is found that the relations are in agreement with predictions from basic physical process of granulation and fall close to the extrapolated relations of RGB stars. The granulations in most of the RSGs evolve at a timescale of several days to a year with a characteristic amplitude of 10–1000 mmag. The results imply that the irregular variations of RSGs can be attributed to the evolution of granulations. When comparing the results from the SMC, LMC, and M31, the timescale and amplitude of granulation seem to increase with metallicity. The analytical relations of the granulation parameters with stellar parameters are derived for the RSG sample of each galaxy.

A Forward Modeling Approach to AGN Variability--Method Description and Early Applications

Abstract We present a numerical framework for the variability of active galactic nuclei (AGNs), which links the variability of AGNs over a broad range of timescales and luminosities to the observed properties of the AGN population as a whole, and particularly the Eddington ratio distribution function. We implemented our framework on GPU architecture, relying on previously published time-series-generating algorithms. After extensive tests that characterize several intrinsic and numerical aspects of the simulations, we describe some applications used for current and future time-domain surveys and for the study of extremely variable sources (e.g., “changing-look” or flaring AGNs). Specifically, we define a simulation setup that reproduces the AGN variability observed in the (intermediate) Palomar Transient Factory survey and use it to forward model longer light curves of the kind that may be observed within the Large Synoptic Survey Telescope (LSST) main survey. Thanks to our efficient implementations, these simulations are able to cover, for example, over 1 Myr with a roughly weekly cadence. We envision that this framework will become highly valuable to prepare for, and best exploit, data from upcoming time-domain surveys, such as, for example, LSST.

The volumetric rate of normal type Ia supernovae in the local Universe discovered by the Palomar Transient Factory

Abstract We present the volumetric rate of normal type Ia supernovae (SNe Ia) discovered by the Palomar Transient Factory (PTF). Using strict data-quality cuts, and considering only periods when the PTF maintained a regular cadence, PTF discovered 90 SNe Ia at z ≤ 0.09 in a well-controlled sample over three years of operation (2010–2012). We use this to calculate the volumetric rate of SN Ia events by comparing this sample to simulations of hundreds of millions of SN Ia light curves produced in statistically representative realizations of the PTF survey. This quantifies the recovery efficiency of each PTF SN Ia event, and thus the relative weighting of each event. From this, the volumetric SN Ia rate was found to be $r_\mathrm{ v}=2.43\pm 0.29\, \text{(stat)}_{-0.19}^{+0.33}\text{(sys)}\times 10^{-5}\, \rm{SNe\,yr}^{-1}\, \text{Mpc}^{-3}\, h_{70}^{3}$. This represents the most precise local measurement of the SN Ia rate. We fit a simple SN Ia delay-time distribution model, ∝ t−β, to our PTF rate measurement combined with a literature sample of rate measurements from surveys at higher redshifts. We find β ∼ 1, consistent with a progenitor channel governed by the gravitational inspiral of binary white dwarfs.

The Broad Absorption Line Tidal Disruption Event iPTF15af: Optical and Ultraviolet Evolution

Abstract We present multiwavelength observations of the tidal disruption event (TDE) iPTF15af, discovered by the intermediate Palomar Transient Factory survey at redshift z = 0.07897. The optical and ultraviolet (UV) light curves of the transient show a slow decay over 5 months, in agreement with previous optically discovered TDEs. It also has a comparable blackbody peak luminosity of . The inferred temperature from the optical and UV data shows a value of (3–5) × 104 K. The transient is not detected in X-rays up to within the first 5 months after discovery. The optical spectra exhibit two distinct broad emission lines in the He ii region, and at later times also Hα emission. Additionally, emission from [N iii] and [O iii] is detected, likely produced by the Bowen fluorescence effect. UV spectra reveal broad emission and absorption lines associated with high-ionization states of N v, C iv, Si iv, and possibly P v. These features, analogous to those of broad absorption line quasars (BAL QSOs), require an absorber with column densities cm−2. This optically thick gas would also explain the nondetection in soft X-rays. The profile of the absorption lines with the highest column density material at the largest velocity is opposite that of BAL QSOs. We suggest that radiation pressure generated by the TDE flare at early times could have provided the initial acceleration mechanism for this gas. Spectral UV line monitoring of future TDEs could test this proposal.

Light Curves of Hydrogen-poor Superluminous Supernovae from the Palomar Transient Factory

Abstract We investigate the light-curve properties of a sample of 26 spectroscopically confirmed hydrogen-poor superluminous supernovae (SLSNe-I) in the Palomar Transient Factory survey. These events are brighter than SNe Ib/c and SNe Ic-BL, on average, by about 4 and 2 mag, respectively. The peak absolute magnitudes of SLSNe-I in rest-frame g band span −22 ≲ M g ≲ −20 mag, and these peaks are not powered by radioactive 56Ni, unless strong asymmetries are at play. The rise timescales are longer for SLSNe than for normal SNe Ib/c, by roughly 10 days, for events with similar decay times. Thus, SLSNe-I can be considered as a separate population based on photometric properties. After peak, SLSNe-I decay with a wide range of slopes, with no obvious gap between rapidly declining and slowly declining events. The latter events show more irregularities (bumps) in the light curves at all times. At late times, the SLSN-I light curves slow down and cluster around the 56Co radioactive decay rate. Powering the late-time light curves with radioactive decay would require between 1 and 10 M ⊙ of Ni masses. Alternatively, a simple magnetar model can reasonably fit the majority of SLSNe-I light curves, with four exceptions, and can mimic the radioactive decay of 56Co, up to ∼400 days from explosion. The resulting spin values do not correlate with the host-galaxy metallicities. Finally, the analysis of our sample cannot strengthen the case for using SLSNe-I for cosmology.

Variability of Red Supergiants in M31 from the Palomar Transient Factory

Abstract Most massive stars end their lives as red supergiants (RSGs), a short-lived evolutionary phase when they are known to pulsate with varying amplitudes. The RSG period–luminosity (PL) relation has been measured in the Milky Way, the Magellanic Clouds and M33 for about 120 stars in total. Using over 1500 epochs of R-band monitoring from the Palomar Transient Factory survey over a five-year period, we study the variability of 255 spectroscopically cataloged RSGs in M31. We find that all RGSs brighter than M K ≈ −10 mag (log(L/L ⊙) &gt; 4.8) are variable at Δm R &gt; 0.05 mag. Our period analysis finds 63 with significant pulsation periods. Using the periods found and the known values of M K for these stars, we derive the RSG PL relation in M31 and show that it is consistent with those derived earlier in other galaxies of different metallicities. We also detect, for the first time, a sequence of likely first-overtone pulsations. Comparison to stellar evolution models from MESA confirms the first-overtone hypothesis and indicates that the variable stars in this sample have 12 M ⊙ &lt; M &lt; 24 M ⊙. As these RSGs are the immediate progenitors to Type II-P core-collapse supernovae (SNe), we also explore the implication of their variability in the initial-mass estimates for SN progenitors based on archival images of the progenitors. We find that this effect is small compared to the present measurement errors.

The Volumetric Rate of Calcium-rich Transients in the Local Universe

Abstract We present a measurement of the volumetric rate of “calcium-rich” optical transients in the local universe, using a sample of three events from the Palomar Transient Factory (PTF). This measurement builds on a detailed study of the PTF transient detection efficiencies and uses a Monte Carlo simulation of the PTF survey. We measure the volumetric rate of calcium-rich transients to be higher than previous estimates: events yr−1 Mpc−3. This is equivalent to 33%–94% of the local volumetric Type Ia supernova rate. This calcium-rich transient rate is sufficient to reproduce the observed calcium abundances in galaxy clusters, assuming an asymptotic calcium yield per calcium-rich event of ∼0.05 . We also study the PTF detection efficiency of these transients as a function of position within their candidate host galaxies. We confirm as a real physical effect previous results that suggest that calcium-rich transients prefer large physical offsets from their host galaxies.

Searching for Be Stars in the Open Clusters with PTF/iPTF. I. Cluster Sample and Be Star Candidates

Abstract We conducted a search for Be star candidates in open clusters using Hα imaging photometry of the Palomar Transient Factory Survey to investigate some connections among Be star phenomena, cluster environments, and ages. Stellar members of clusters were identified by spatial distributions, near-infrared magnitudes and colors, and by proper motions. Among 104 open clusters, we identified 96 Be star candidates in 32 clusters; 11 of our candidates have been reported in previous studies. We found that the clusters with age 7.5 &lt; log(t(year)) 8.5 tend to have more Be star candidates; there is about a 40% occurrence rate within this age bin. The clusters in this age bin also tend to have a higher Be fraction N(Be)/N(Be+B-type). These results suggest that the environments of young and intermediate clusters are favorable to the formation of Be phenomena. Spatial distribution of Be star candidates with different ages implies that they do not form preferentially in the central regions. Furthermore, we showed that the mid-infrared (MIR) colors of the Be star candidates are similar to known Be stars, which could be caused by free–free emission or bound-free emission. Some Be star candidates might have no circumstellar dust according to their MIR colors. Finally, among 96 Be candidates, we discovered that one Be star candidate FSR 0904-1 exhibits long-term variability on the timescale of ∼2000 days with an amplitude of 0.2–0.3 mag, indicating a long timescale of disk evolution.

Be STARS IN THE OPEN CLUSTER NGC 6830

ABSTRACT We report the discovery of two new Be stars, and re-identify one known Be star in the open cluster NGC 6830. Eleven Hα emitters were discovered using the Hα imaging photometry of the Palomar Transient Factory Survey. Stellar membership of the candidates was verified with photometric and kinematic information using 2MASS data and proper motions. The spectroscopic confirmation was carried out by using the Shane 3 m telescope at the Lick observatory. Based on their spectral types, three Hα emitters were confirmed as Be stars with Hα equivalent widths greater than −10 Å. Two objects were also observed by the new spectrograph spectral energy distribution-machine (SED-machine) on the Palomar 60-inch Telescope. The SED-machine results show strong Hα emission lines, which are consistent with the results of the Lick observations. The high efficiency of the SED-machine can provide rapid observations for Be stars in a comprehensive survey in the future.

SEARCHING FOR Be STARS IN THE OPEN CLUSTER NGC 663

We present Be star candidates in the open cluster NGC\,663, identified by H$\alpha$ imaging photometry with the Palomar Transient Factory Survey, as a pilot program to investigate how the Be star phenomena, the emission spectra, extended circumstellar envelopes, and fast rotation, correlate with massive stellar evolution. Stellar membership of the candidates was verified by 2MASS magnitudes and colors, and by proper motions. We discover 4 new Be stars and exclude one known Be star from being a member due to its inconsistent proper motions. The fraction of Be stars to member stars [N(Be)/N(members)] in NGC\,663 is 3.5\%. The spectral type of the 34 Be stars in NGC\,663 shows bimodal peaks at B0--B2 and B5--B7, which is consistent with the statistics in most star clusters. Additionally, we also discover 23 emission-line stars of different types, including non-member Be stars, dwarfs, and giants.

Main-belt comets in the Palomar Transient Factory survey – I. The search for extendedness

Cometary activity in main-belt asteroids probes the ice content of these objects and provides clues to the history of volatiles in the inner solar system. We search the Palomar Transient Factory (PTF) survey to derive upper limits on the population size of active main-belt comets (MBCs). From data collected March 2009 through July 2012, we extracted 2 million observations of 220 thousand known main-belt objects (40% of the known population, down to 1-km diameter) and discovered 626 new objects in multi-night linked detections. We formally quantify the extendedness of a small-body observation, account for systematic variation in this metric (e.g., due to on-sky motion) and evaluate this method's robustness in identifying cometary activity using observations of 115 comets, including two known candidate MBCs and six newly-discovered non-main-belt comets (two of which were originally designated as asteroids by other surveys). We demonstrate a 66% detection efficiency with respect to the extendedness distribution of the 115 sampled comets, and a 100% detection efficiency with respect to extendedness levels greater than or equal to those we observed in the known candidate MBCs P/2010 R2 (La Sagra) and P/2006 VW139. Using a log-constant prior, we infer 95% confidence upper limits of 33 and 22 active MBCs (per million main-belt asteroids down to 1-km diameter), for detection efficiencies of 66% and 100%, respectively. In a follow-up to this morphological search, we will perform a photometric (disk-integrated brightening) search for MBCs.

PTF 11kx: A Type Ia Supernova with a Symbiotic Nova Progenitor

There is a consensus that type Ia supernovae (SNe Ia) arise from the thermonuclear explosion of white dwarf stars that accrete matter from a binary companion. However, direct observation of SN Ia progenitors is lacking, and the precise nature of the binary companion remains uncertain. A temporal series of high-resolution optical spectra of the SN Ia PTF 11kx reveals a complex circumstellar environment that provides an unprecedentedly detailed view of the progenitor system. Multiple shells of circumstellar material are detected, and the SN ejecta are seen to interact with circumstellar material starting 59 days after the explosion. These features are best described by a symbiotic nova progenitor, similar to RS Ophiuchi.

Luminous Supernovae

Stellar Explosions Stars that are born with masses greater than eight times that of the Sun end their lives in luminous explosions known as supernovae. Over the past decade, access to improved sky surveys has revealed rare types of supernovae that are much more luminous than any of those that were known before. Gal-Yam (p. 927 ) reviews these superluminous events and groups them into three classes that share common observational and physical characteristics. Gamma-ray bursts are another type of extreme explosive events related to the death of massive stars, which occur once per day somewhere in the universe and produce short-lived bursts of gamma-ray light. Gehrels and Mészáros (p. 932 ) review what has been learned about these events since the launch of NASA's Swift (2004) and Fermi (2008) satellites. The current interpretation is that gamma-ray bursts are related to the formation of black holes. Type Ia supernovae are used as cosmological distance indicators. They are thought to be the result of the thermonuclear explosion of white dwarf stars in binary systems, but the nature of the stellar companion to the white dwarf is still debated. Dilday et al. (p. 942 ) report high-resolution spectroscopy of the supernova PTF 11kx, which was detected on 26 January 2011 by the Palomar Transient Factory survey. The data suggest a red giant star companion whose material got transferred to the white dwarf.

Gamma-Ray Bursts

Gamma-ray bursts (GRBs) are bright flashes of gamma rays coming from the cosmos. They occur roughly once per day, typically last for tens of seconds, and are the most luminous events in the universe. More than three decades after their discovery, and after pioneering advances from space and ground experiments, they still remain mysterious. The launch of the Swift and Fermi satellites in 2004 and 2008 brought in a trove of qualitatively new data. In this Review, we survey the interplay between these recent observations and the theoretical models of the prompt GRB emission and the subsequent afterglow.

THE PALOMAR TRANSIENT FACTORY ORION PROJECT: ECLIPSING BINARIES AND YOUNG STELLAR OBJECTS

The Palomar Transient Factory (PTF) Orion project is one of the experiments within the broader PTF survey, a systematic automated exploration of the sky for optical transients. Taking advantage of the wide (35 × 23) field of view available using the PTF camera installed at the Palomar 48 inch telescope, 40 nights were dedicated in 2009 December to 2010 January to perform continuous high-cadence differential photometry on a single field containing the young (7–10 Myr) 25 Ori association. Little is known empirically about the formation of planets at these young ages, and the primary motivation for the project is to search for planets around young stars in this region. The unique data set also provides for much ancillary science. In this first paper, we describe the survey and the data reduction pipeline, and present some initial results from an inspection of the most clearly varying stars relating to two of the ancillary science objectives: detection of eclipsing binaries and young stellar objects. We find 82 new eclipsing binary systems, 9 of which are good candidate 25 Ori or Orion OB1a association members. Of these, two are potential young W UMa type systems. We report on the possible low-mass (M-dwarf primary) eclipsing systems in the sample, which include six of the candidate young systems. Forty-five of the binary systems are close (mainly contact) systems, and one of these shows an orbital period among the shortest known for W UMa binaries, at 0.2156509 ± 0.0000071 days, with flat-bottomed primary eclipses, and a derived distance that appears consistent with membership in the general Orion association. One of the candidate young systems presents an unusual light curve, perhaps representing a semi-detached binary system with an inflated low-mass primary or a star with a warped disk, and may represent an additional young Orion member. Finally, we identify 14 probable new classical T-Tauri stars in our data, along with one previously known (CVSO 35) and one previously reported as a candidate weak-line T-Tauri star (SDSS J052700.12+010136.8).