Search

Fainter standard stars are essential for the calibration of larger telescopes. This work adds to the CALSPEC (calibration spectra) database 19 faint white dwarfs (WDs) with all-sky coverage and V magnitudes between 16.5 and 18.7. Included for these stars is new UV (ultraviolet) HST (Hubble Space Telescope) STIS (Space Telescope Imaging Spectrometer) spectrophotometry between 1150 and 3000~\AA\ with a resolution of $\sim$500. Pure hydrogen WD models are fit to these UV spectra and to six-band HST/WFC3 (Wide Field Camera 3) photometry at 0.28 to 1.6~\micron\ to construct predicted model SEDs (spectral energy distributions) covering wavelengths from 900~\AA\ to the JWST (James Webb Space Telescope) limit of 30~\micron\ using well-established CALSPEC procedures for producing flux standards with the goal of 1\% accuracy.

We introduce a pipeline that performs rapid image subtraction and source selection to detect transients, with a focus on identifying gravitational wave optical counterparts using the Dark Energy Camera (DECam). In this work, we present the pipeline steps from processing raw data to identification of astrophysical transients on individual exposures. We process DECam data and build difference images using the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) Science Pipelines software, and we use flags and principal component analysis to select transients on a per-exposure basis, without associating the results from different exposures. Those candidates will be sent to brokers for further classification and alert distribution. We validate our pipeline using archival exposures that cover various types of objects, and the tested targets include a kilonova (GW170817), supernovae, stellar flares, variable stars (in a resolved galaxy or the Milky Way Bulge), and serendipitous objects. Overall, the data processing produces clean light curves that are comparable with published results, demonstrating the photometric quality of our pipeline. Real transients can be well selected by our pipeline when sufficiently bright (S/N $\gtrsim15$). This pipeline is intended to serve as a tool for the broader research community. Although this pipeline is designed for DECam, our method can be easily applied to other instruments and future LSST observations., Comment: 37 pages, 25 figures, 5 tables; revised and accepted for publication in AJ

Light echoes give us a unique perspective on the nature of supernovae and non-terminal stellar explosions. Spectroscopy of light echoes can reveal details on the kinematics of the ejecta, probe asymmetry, and reveal details on its interaction with circumstellar matter, thus expanding our understanding of these transient events. However, the spectral features arise from a complex interplay between the source photons, the reflecting dust geometry, and the instrumental setup and observing conditions. In this work we present an improved method for modeling these effects in light echo spectra, one that relaxes the simplifying assumption of a light curve weighted sum, and instead estimates the true relative contribution of each phase. We discuss our logic, the gains we obtain over light echo analysis method(s) used in the past, and prospects for further improvements. Lastly, we show how the new method improves our analysis of echoes from Tycho's supernova (SN 1572) as an example.

We follow our first paper with an analysis of the ensemble of the extensive pre-explosion ground- and space-based infrared observations of the red supergiant (RSG) progenitor candidate for the nearby core-collapse supernova SN 2023ixf in Messier 101, together with optical data prior to explosion obtained with the Hubble Space Telescope (HST). We have confirmed the association of the progenitor candidate with the SN, as well as constrained the metallicity at the SN site, based on SN observations with instruments at Gemini-North. The internal host extinction to the SN has also been confirmed from a high-resolution Keck spectrum. We fit the observed spectral energy distribution (SED) for the star, accounting for its intrinsic variability, with dust radiative-transfer modeling, which assume a silicate-rich dust shell ahead of the underlying stellar photosphere. The star is heavily dust-obscured, likely the dustiest progenitor candidate yet encountered. We found median estimates of the star's effective temperature and luminosity of 2770 K and 9.0e4 L_Sun, with 68% credible intervals of 2340--3150 K and (7.5--10.9)e4 L_sun. The candidate may have a Galactic RSG analog, IRC -10414, with a strikingly similar SED and luminosity. Via comparison with single-star evolutionary models we have constrained the initial mass of the progenitor candidate from 12 M_sun to as high as 14 M_sun. We have had available to us an extraordinary view of the SN 2023ixf progenitor candidate, which should be further followed up in future years with HST and the James Webb Space Telescope., Comment: 40 pages, substantive modifications relative to the previous, although the overall conclusions remain the same; to appear in AAS Journals

Observational evidence points to a red supergiant (RSG) progenitor for SN 2023ixf. The progenitor candidate has been detected in archival images at wavelengths (>0.6 micron) where RSGs typically emit profusely. This object is distinctly variable in the infrared (IR). We characterize the variability using pre-explosion mid-IR (3.6 and 4.5 micron) Spitzer and ground-based near-IR (JHKs) archival data jointly covering 19 yr. The IR light curves exhibit significant variability with RMS amplitudes in the range of 0.2-0.4 mag, increasing with decreasing wavelength. From a robust period analysis of the more densely sampled Spitzer data, we measure a period of 1091+/-71 days. We demonstrate using Gaussian Process modeling that this periodicity is also present in the near-IR light curves, thus indicating a common physical origin, which is likely pulsational instability. We use a period-luminosity relation for RSGs to derive a value of M_K=-11.58+/-0.31 mag. Assuming a late M spectral type, this corresponds to log(L/L_sun)=5.27+/-0.12 at T_eff=3200 K and to log(L/L_sun)=5.37+/-0.12 at T_eff=3500 K. This gives an independent estimate of the progenitor's luminosity, unaffected by uncertainties in extinction and distance. Assuming the progenitor candidate underwent enhanced dust-driven mass-loss during the time of these archival observations, and using an empirical period-luminosity-based mass-loss prescription, we obtain a mass-loss rate of around (2-4)x10^-4 M_sun/yr. Comparing the above luminosity with stellar evolution models, we infer an initial mass for the progenitor candidate of 20+/-4 M_sun, making this one of the most massive progenitors for a Type II SN detected to-date., Comment: 20 pages, 2 tables, accepted to ApJ

Hot DA white dwarfs have fully radiative pure hydrogen atmospheres that are the least complicated to model. Pulsationally stable, they are fully characterized by their effective temperature Teff, and surface gravity log g, which can be deduced from their optical spectra and used in model atmospheres to predict their spectral energy distribution (SED). Based on this, three bright DAWDs have defined the spectrophotometric flux scale of the CALSPEC system of HST. In this paper we add 32 new fainter (16.5 < V < 19.5) DAWDs spread over the whole sky and within the dynamic range of large telescopes. Using ground based spectra and panchromatic photometry with HST/WFC3, a new hierarchical analysis process demonstrates consistency between model and observed fluxes above the terrestrial atmosphere to < 0.004 mag rms from 2700 {\AA} to 7750 {\AA} and to 0.008 mag rms at 1.6{\mu}m for the total set of 35 DAWDs. These DAWDs are thus established as spectrophotometric standards with unprecedented accuracy from the near ultraviolet to the near-infrared, suitable for both ground and space based observatories. They are embedded in existing surveys like SDSS, PanSTARRS and GAIA, and will be naturally included in the LSST survey by Rubin Observatory. With additional data and analysis to extend the validity of their SEDs further into the IR, these spectrophotometric standard stars could be used for JWST, as well as for the Roman and Euclid observatories., Comment: Accepted for publication in Astrophysical Journal. Corrected error in Table 10 and associated Fig 7 in which RP and BP values for the 3 CALSPEC standards had been transposed

We report early-time ultraviolet (UV) and optical spectroscopy of the young, nearby Type II supernova (SN) 2022wsp obtained by the Hubble Space Telescope (HST)/STIS at about 10 and 20 days after the explosion. The SN 2022wsp UV spectra are compared to those of other well-observed Type II/IIP SNe, including the recently studied Type IIP SN 2021yja. Both SNe exhibit rapid cooling and similar evolution during early phases, indicating a common behavior among SNe II. Radiative-transfer modeling of the spectra of SN 2022wsp with the TARDIS code indicates a steep radial density profile in the outer layer of the ejecta, a supersolar metallicity, and a relatively high total extinction of E(B-V) = 0.35 mag. The early-time evolution of the photospheric velocity and temperature derived from the modeling agree with the behavior observed from other previously studied cases. The strong suppression of hydrogen Balmer lines in the spectra suggests interaction with a pre-existing circumstellar environment could be occurring at early times. In the SN 2022wsp spectra, the absorption component of the Mg II P Cygni profile displays a double-trough feature on day +10 that disappears by day +20. The shape is well reproduced by the model without fine-tuning the parameters, suggesting that the secondary blueward dip is a metal transition that originates in the SN ejecta., Comment: Submitted to ApJ Letters on 4/11/2023

This paper presents a new optical imaging survey of four deep drilling fields (DDFs), two Galactic and two extragalactic, with the Dark Energy Camera (DECam) on the 4 meter Blanco telescope at the Cerro Tololo Inter-American Observatory (CTIO). During the first year of observations in 2021, $>$4000 images covering 21 square degrees (7 DECam pointings), with $\sim$40 epochs (nights) per field and 5 to 6 images per night per filter in $g$, $r$, $i$, and/or $z$, have become publicly available (the proprietary period for this program is waived). We describe the real-time difference-image pipeline and how alerts are distributed to brokers via the same distribution system as the Zwicky Transient Facility (ZTF). In this paper, we focus on the two extragalactic deep fields (COSMOS and ELAIS-S1), characterizing the detected sources and demonstrating that the survey design is effective for probing the discovery space of faint and fast variable and transient sources. We describe and make publicly available 4413 calibrated light curves based on difference-image detection photometry of transients and variables in the extragalactic fields. We also present preliminary scientific analysis regarding Solar System small bodies, stellar flares and variables, Galactic anomaly detection, fast-rising transients and variables, supernovae, and active galactic nuclei., Comment: 22 pages, 17 figures, 2 tables. Accepted to MNRAS

We present three epochs of early-time ultraviolet (UV) and optical HST/STIS spectroscopy of the young, nearby Type IIP supernova (SN) 2021yja. We complement the HST data with two earlier epochs of Swift UVOT spectroscopy. The HST and Swift UVOT spectra are consistent with those of other well-studied Type IIP supernovae (SNe). The UV spectra exhibit rapid cooling at early times, while less dramatic changes are seen in the optical. We also present Lick/KAIT optical photometry up to the late-time-tail phase, showing a very long plateau and shallow decline compared with other SNe IIP. Our modeling of the UV spectrum with the TARDIS radiative-transfer code produces a good fit for a high-velocity explosion, a low total extinction $E(B-V) = 0.07$ mag, and a subsolar metallicity. We do not find a significant contribution to the UV flux from an additional heating source, such as interaction with the circumstellar medium, consistent with the observed flat plateau. Furthermore, the velocity width of the Mg II $\lambda$2798 line is comparable to that of the hydrogen Balmer lines, suggesting that the UV emission is confined to a region close to the photosphere., Comment: Accepted to ApJ

Ongoing large-scale optical time-domain surveys, such as the Zwicky Transient Facility (ZTF), are producing alerts at unprecedented rates. Analysis of transient sources has so far followed two distinct paths: archival analysis of data on transient sources at a time when they are no longer observable and real-time analysis at the time when the sources are first detected. The latter is the realm of alert brokers such as the Arizona-NOIRLab Temporal Analysis and Response to Events System (ANTARES). In this paper, we synthesize the two analysis paths and present a first systematic study of archival alert-broker data, focusing on extragalactic transients with multipeaked light curves identified in the ANTARES archive of ZTF alerts. Our analysis yields a sample of 37 such sources, including core-collapse supernovae (with two analogs of iPTF14hls), thermonuclear supernovae interacting with their surrounding circumstellar medium, tidal disruption events, luminous blue variables, and as yet unclassified objects. A large fraction of the identified sources is currently active, warranting allocation of follow-up resources in the immediate future to further constrain their nature and the physical processes at work., Comment: Accepted for publication in ApJ Letters

We recently discovered an extragalactic transient, AT2020caa, using the community alert broker ANTARES. This transient apparently exhibited two outbursts in a time span of a year (between 2020 and 2021). Based on a decade-long historical light curve of the candidate host galaxy, we rule out an activity from the galaxy nucleus to explain these outbursts. The measured peak magnitudes (assuming the known spectroscopic redshift of the candidate host galaxy) put AT2020caa in the realm of thermonuclear supernovae (SNe) or luminous core-collapse SNe. A handful of the latter are known to show prior outbursts (POs), thought to be linked to mass loss in massive stars. Using Gemini/GMOS, we obtained a spectrum of the current outburst that shows it to be a Type Ia supernova (SNIa). We examine the nature of AT2020caa's PO and conclude that it is likely a separate SN within the same galaxy., Comment: submitted to RNAAS

We describe the Arizona-NOIRLab Temporal Analysis and Response to Events System (ANTARES), a software instrument designed to process large-scale streams of astronomical time-domain alerts. With the advent of large-format CCDs on wide-field imaging telescopes, time-domain surveys now routinely discover tens of thousands of new events each night, more than can be evaluated by astronomers alone. The ANTARES event broker will process alerts, annotating them with catalog associations and filtering them to distinguish customizable subsets of events. We describe the data model of the system, the overall architecture, annotation, implementation of filters, system outputs, provenance tracking, system performance, and the user interface., Comment: 24 Pages, 8 figures, Accepted by AJ

The ongoing Zwicky Transient Facility (ZTF) survey is generating a massive alert rate from a variety of optical transients and variable stars, which are being filtered down to subsets meeting user-specified criteria by broker systems such as ANTARES. In a beta implementation of the algorithm of Soraisam et al. (2020) on ANTARES, we flagged AT 2020iko from the ZTF real-time alert stream as an anomalous source. This source is located close to a red extended SDSS source. In the first few epochs of detection, it exhibited a V-shaped brightness profile, preceded by non-detections both in ZTF and in ASASSN extending to 2014. Its full light curve shows a precursor event, followed by a main superoutburst and at least two rebrightenings. A low-resolution spectrum of this source points to a dwarf nova (DN) nature. Although some of the features of AT 2020iko indicate an SU UMa-type DN, its large amplitude, presence of rebrightenings, and inferred supercycle period of > 6 yr are in favor of AT 2020iko being a new WZ Sge-type dwarf nova candidate, a subset of rare DNe consisting of extreme mass-ratio (< 0.1) binaries with orbital period around the period minimum. AT 2020iko's precursor event brightened by 6.5 mag, while its decay spanned 3-5 mag. We speculate this superoutburst is associated with a less expanded accretion disk than in typical superoutbursts in WZ Sge systems, with the large depth of the precursor decay implying an extremely small mass-ratio. To the best of our knowledge, such a precursor event has not been recorded for any DN. This result serves to demonstrate the efficacy of our real-time anomaly search algorithm., Comment: 11 pages, 4 figures; accepted to AJ

Several thousand core-collapse supernovae (CCSNe) of different flavors have been discovered so far. However, identifying their progenitors has remained an outstanding open question in astrophysics. Studies of SN host galaxies have proven to be powerful in providing constraints on the progenitor populations. In this paper, we present all CCSNe detected between 2009 and 2017 by the Palomar Transient Factory. This sample includes 888 SNe of 12 distinct classes out to redshift $z\approx1$. We present the photometric properties of their host galaxies from the far-ultraviolet to the mid-infrared and model the host-galaxy spectral energy distributions to derive physical properties. The galaxy mass functions of Type Ic, Ib, IIb, II, and IIn SNe ranges from $10^{5}$ to $10^{11.5}~M_\odot$, probing the entire mass range of star-forming galaxies down to the least-massive star-forming galaxies known. Moreover, the galaxy mass distributions are consistent with models of star-formation-weighted mass functions. Regular CCSNe are hence direct tracers of star formation. Small but notable differences exist between some of the SN classes. Type Ib/c SNe prefer galaxies with slightly higher masses (i.e., higher metallicities) and star-formation rates than Type IIb and II SNe. These differences are less pronounced than previously thought. H-poor SLSNe and SNe~Ic-BL are scarce in galaxies above $10^{10}~M_\odot$. Their progenitors require environments with metallicities of $<0.4$ and $<1$ solar, respectively. In addition, the hosts of H-poor SLSNe are dominated by a younger stellar population than all other classes of CCSNe. Our findings corroborate the notion that low-metallicity \textit{and} young age play an important role in the formation of SLSN progenitors., Comment: 24 pages main text, 14 figures, 9 Tables, catalogue available at http://www.github.com/steveschulze/PTF

We analyze 143 Type Ia supernovae (SNeIa) observed in $H$ band (1.6-1.8 $\mu$m) and find SNeIa are intrinsically brighter in $H$-band with increasing host galaxy stellar mass. We find SNeIa in galaxies more massive than $10^{10.43} M_{\odot}$ are $0.13 \pm 0.04$ mag brighter in $H$ than SNeIa in less massive galaxies. The same set of SNeIa observed at optical wavelengths, after width-color-luminosity corrections, exhibit a $0.10 \pm 0.03$ mag offset in the Hubble residuals. We observe an outlier population ($|\Delta H_{\rm max}| > 0.5$ mag) in the $H$ band and show that removing the outlier population moves the mass threshold to $10^{10.65} M_{\odot}$ and reduces the step in $H$ band to $0.08 \pm 0.04$ mag, but the equivalent optical mass step is increased to $0.13 \pm 0.04$ mag. We conclude the outliers do not drive the brightness--host-mass correlation. Less massive galaxies preferentially host more higher-stretch SNeIa, which are intrinsically brighter and bluer. It is only after correction for width-luminosity and color-luminosity relationships that SNeIa have brighter optical Hubble residuals in more massive galaxies. Thus finding SNeIa are intrinsically brighter in $H$ in more massive galaxies is an opposite correlation to the intrinsic (pre-width-luminosity correction) optical brightness. If dust and the treatment of intrinsic color variation were the main driver of the host galaxy mass correlation, we would not expect a correlation of brighter $H$-band SNeIa in more massive galaxies., Comment: 54 pages, 21 Figures, 15 Tables. Accepted to the Astrophysical Journal

With the advent of the Large Synoptic Survey Telescope (LSST), time-domain astronomy will be faced with an unprecedented volume and rate of data. Real-time processing of variables and transients detected by such large-scale surveys is critical to identifying the more unusual events and allocating scarce follow-up resources efficiently. We develop an algorithm to identify these novel events within a given population of variable sources. We determine the distributions of magnitude changes (dm) over time intervals (dt) for a given passband f, pf(dm|dt), and use these distributions to compute the likelihood of a test source being consistent with the population, or an outlier. We demonstrate our algorithm by applying it to the DECam multi-band time-series data of more than 2000 variable stars identified by Saha et al. (2019) in the Galactic Bulge that are largely dominated by long-period variables and pulsating stars. Our algorithm discovers 18 outlier sources in the sample, including a microlensing event, a dwarf nova, and two chromospherically active RS CVn stars, as well as sources in the Blue Horizontal Branch region of the color-magnitude diagram without any known counterparts. We compare the performance of our algorithm for novelty detection with multivariate KDE and Isolation Forest on the simulated PLAsTiCC dataset. We find that our algorithm yields comparable results despite its simplicity. Our method provides an efficient way for flagging the most unusual events in a real-time alert-broker system., Comment: 20 pages + Appendix, accepted to ApJ

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network., Comment: White paper submitted to the Astro2020 Decadal Survey

We have obtained repeated images of 6 fields towards the Galactic bulge in 5 passbands (u, g, r, i, z) with the DECam imager on the Blanco 4m telescope at CTIO. From over 1.6 billion individual photometric measurements in the field centered on Baade's window, we have detected 4877 putative variable stars. 474 of these have been confirmed as fundamental mode RR Lyrae stars, whose colors at minimum light yield line-of-sight reddening determinations as well as a reddening law towards the Galactic Bulge which differs significantly from the standard R_V = 3.1 formulation. Assuming that the stellar mix is invariant over the 3 square-degree field, we are able to derive a line-of-sight reddening map with sub-arcminute resolution, enabling us to obtain de-reddened and extinction corrected color-magnitude diagrams (CMD's) of this bulge field using up to 2.5 million well-measured stars. The corrected CMD's show unprecedented detail and expose sparsely populated sequences: e.g., delineation of the very wide red giant branch, structure within the red giant clump, the full extent of the horizontal branch, and a surprising bright feature which is likely due to stars with ages younger than 1 Gyr. We use the RR Lyrae stars to trace the spatial structure of the ancient stars, and find an exponential decline in density with Galactocentric distance. We discuss ways in which our data products can be used to explore the age and metallicity properties of the bulge, and how our larger list of all variables is useful for learning to interpret future LSST alerts., Comment: accepted for publication in AAS journals; most likely The Astrophysical Journal

We have established a network of 19 faint (16.5 mag $< V < $19 mag) northern and equatorial DA white dwarfs as spectrophotometric standards for present and future wide-field observatories. Our analysis infers SED models for the stars that are tied to the three CALSPEC primary standards. Our SED models are consistent with panchromatic Hubble Space Telescope ($HST$) photometry to better than 1%. The excellent agreement between observations and models validates the use of non-local-thermodynamic-equilibrium (NLTE) DA white dwarf atmospheres extinguished by interstellar dust as accurate spectrophotometric references. Our standards are accessible from both hemispheres and suitable for ground and space-based observatories covering the ultraviolet to the near infrared. The high-precision of these faint sources make our network of standards ideally suited for any experiment that has very stringent requirements on flux calibration, such as studies of dark energy using the Large Synoptic Survey Telescope (LSST) and the Wide-Field Infrared Survey Telescope ($WFIRST$)., Comment: 46 pages, 23 figures, 8 tables, accepted for publication in ApJS

We present the complete sample of stripped-envelope supernova (SN) spectra observed by the Lick Observatory Supernova Search (LOSS) collaboration over the last three decades: 888 spectra of 302 SNe, 652 published here for the first time, with 384 spectra (of 92 SNe) having photometrically-determined phases. After correcting for redshift and Milky Way dust reddening and reevaluating the spectroscopic classifications for each SN, we construct mean spectra of the three major spectral subtypes (Types IIb, Ib, and Ic) binned by phase. We compare measures of line strengths and widths made from this sample to the results of previous efforts, confirming that O I {\lambda}7774 absorption is stronger and found at higher velocity in Type Ic SNe than in Types Ib or IIb SNe in the first 30 days after peak brightness, though the widths of nebular emission lines are consistent across subtypes. We also highlight newly available observations for a few rare subpopulations of interest., Comment: 13 pages; 14 figures; 3 tables. Accepted for publication in MNRAS

Most massive stars end their lives as Red Supergiants (RSGs), a short-lived evolution 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 (PTF) 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_sun)>4.8) are variable at dm_R>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_sun

We report the discovery and spectroscopic confirmation of a very large star-forming Lyman Break galaxy, G6025, at z_spec=3.721+/-0.003. In the rest-frame ~2100A, G6025 subtends ~24 kpc in physical extent when measured from the 1.5-sigma isophote, in agreement with the parametric size measurements which yield the half-light radius of 4.9+/-0.5 kpc and the semi-major axis of 12.5+/-0.1 kpc. G6025 is also very UV-luminous (~5L*(z~4}) and young (~140+/-60 Myr). Despite its unusual size and luminosity, the stellar population parameters and dust reddening (M_star~M*(z~4)$, and E(B-V)=0.18+/-0.05) estimated from the integrated light, are similar to those of smaller galaxies at comparable redshifts. The ground-based morphology and spectroscopy show two dominant components, both located off-center, embedded in more diffuse emission. We speculate that G6025 may be a scaled-up version of chain galaxies seen in deep HST imaging, or alternatively, a nearly equal-mass merger involving two super-L* galaxies in its early stage. G6025 lies close to but not within a known massive protocluster at z=3.78. We find four companions within 6 Mpc from G6025, two of which lie within 1.6 Mpc. While the limited sensitivity of the existing spectroscopy does not allow us to robustly characterize the local environment of G6025, it likely resides in a locally overdense environment. The luminosity, size, and youth of G6025 make it uniquely suited to study the early formation of massive galaxies in the universe., Comment: accepted for publication

Most Type I superluminous supernovae (SLSNe-I) reported to date have been identified by their high peak luminosities and spectra lacking obvious signs of hydrogen. We demonstrate that these events can be distinguished from normal-luminosity SNe (including Type Ic events) solely from their spectra over a wide range of light-curve phases. We use this distinction to select 19 SLSNe-I and 4 possible SLSNe-I from the Palomar Transient Factory archive (including 7 previously published objects). We present 127 new spectra of these objects and combine these with 39 previously published spectra, and we use these to discuss the average spectral properties of SLSNe-I at different spectral phases. We find that Mn II most probably contributes to the ultraviolet spectral features after maximum light, and we give a detailed study of the O II features that often characterize the early-time optical spectra of SLSNe-I. We discuss the velocity distribution of O II, finding that for some SLSNe-I this can be confined to a narrow range compared to relatively large systematic velocity shifts. Mg II and Fe II favor higher velocities than O II and C II, and we briefly discuss how this may constrain power-source models. We tentatively group objects by how well they match either SN 2011ke or PTF12dam and discuss the possibility that physically distinct events may have been previously grouped together under the SLSN-I label., Comment: 70 pages, 41 figures, 7 tables, accepted for publication in ApJ

The unprecedented volume and rate of transient events that will be discovered by the Large Synoptic Survey Telescope (LSST) demands that the astronomical community update its followup paradigm. Alert-brokers -- automated software system to sift through, characterize, annotate and prioritize events for followup -- will be critical tools for managing alert streams in the LSST era. The Arizona-NOAO Temporal Analysis and Response to Events System (ANTARES) is one such broker. In this work, we develop a machine learning pipeline to characterize and classify variable and transient sources only using the available multiband optical photometry. We describe three illustrative stages of the pipeline, serving the three goals of early, intermediate and retrospective classification of alerts. The first takes the form of variable vs transient categorization, the second, a multi-class typing of the combined variable and transient dataset, and the third, a purity-driven subtyping of a transient class. While several similar algorithms have proven themselves in simulations, we validate their performance on real observations for the first time. We quantitatively evaluate our pipeline on sparse, unevenly sampled, heteroskedastic data from various existing observational campaigns, and demonstrate very competitive classification performance. We describe our progress towards adapting the pipeline developed in this work into a real-time broker working on live alert streams from time-domain surveys., Comment: 33 pages, 14 figures, submitted to ApJS

We characterize the absolute magnitudes and colors of RR Lyrae stars in the globular cluster M5 in the ugriz filter system of the Dark Energy Camera (DECam). We provide empirical Period-Luminosity (P-L) relationships in all 5 bands based on 47 RR Lyrae stars of the type ab and 14 stars of the type c. The P-L relationships were found to be better constrained for the fundamental mode RR Lyrae stars in the riz passbands, with dispersion of 0.03, 0.02 and 0.02 magnitudes, respectively. The dispersion of the color at minimum light was found to be small, supporting the use of this parameter as a means to obtain accurate interstellar extinctions along the line of sight up to the distance of the RR Lyrae star. We found a trend of color at minimum light with pulsational period that, if taken into account, brings the dispersion in color at minimum light to < 0.016 magnitudes for the (r-i), (i-z), and (r-z) colors. These calibrations will be very useful for using RR Lyrae stars from DECam observations as both standard candles for distance determinations and color standards for reddening measurements., Comment: Accepted for publication in AJ

Galaxy-cluster gravitational lenses can magnify background galaxies by a total factor of up to ~50. Here we report an image of an individual star at redshift z=1.49 (dubbed "MACS J1149 Lensed Star 1 (LS1)") magnified by >2000. A separate image, detected briefly 0.26 arcseconds from LS1, is likely a counterimage of the first star demagnified for multiple years by a >~3 solar-mass object in the cluster. For reasonable assumptions about the lensing system, microlensing fluctuations in the stars' light curves can yield evidence about the mass function of intracluster stars and compact objects, including binary fractions and specific stellar evolution and supernova models. Dark-matter subhalos or massive compact objects may help to account for the two images' long-term brightness ratio., Comment: Accepted by Nature Astronomy on February 23, 2018

We present multiband photometry of 60 spectroscopically-confirmed supernovae (SN): 39 SN II/IIP, 19 IIn, one IIb and one that was originally classified as a IIn but later as a Ibn. Forty-six have only optical photometry, six have only near infrared (NIR) photometry and eight have both optical and NIR. The median redshift of the sample is 0.016. We also present 192 optical spectra for 47 of the 60 SN. All data are publicly available. There are 26 optical and two NIR light curves of SN II/IIP with redshifts z > 0.01, some of which may give rise to useful distances for cosmological applications. All photometry was obtained between 2000 and 2011 at the Fred Lawrence Whipple Observatory (FLWO), via the 1.2m and 1.3m PAIRITEL telescopes for the optical and NIR, respectively. Each SN was observed in a subset of the $u'UBVRIr'i'JHK_s$ bands. There are a total of 2932 optical and 816 NIR light curve points. Optical spectra were obtained using the FLWO 1.5m Tillinghast telescope with the FAST spectrograph and the MMT Telescope with the Blue Channel Spectrograph. Our photometry is in reasonable agreement with other samples from the literature. Comparison with Pan-STARRS shows that two-thirds of our individual star sequences have weighted-mean V offsets within $\pm$0.02 mag. In comparing our standard-system SN light curves with common Carnegie Supernova Project objects using their color terms, we found that roughly three-quarters have average differences within $\pm$0.04 mag. The data from this work and the literature will provide insight into SN II explosions, help with developing methods for photometric SN classification, and contribute to their use as cosmological distance indicators., Comment: Accepted to ApJS. TAR of light curves and star sequences here: https://www.cfa.harvard.edu/supernova/fmalcolm2017/cfa_snII_lightcurvesndstars.june2017.tar ... Spectra can be found here: https://www.cfa.harvard.edu/supernova/fmalcolm2017/cfaspec_snII.tar.gz ... Passbands and plot of spectra can be found here: https://www.cfa.harvard.edu/supernova/SNarchive.html

We present the results of an extensive observational campaign on the nearby Type Ibn SN 2015G, including data from radio through ultraviolet wavelengths. SN 2015G was asymmetric, showing late-time nebular lines redshifted by ~1000 km/s. It shared many features with the prototypical SN In 2006jc, including extremely strong He I emssion lines and a late-time blue pseudocontinuum. The young SN 2015G showed narrow P-Cygni profiles of He I, but never in its evolution did it show any signature of hydrogen - arguing for a dense, ionized, and hydrogen-free circumstellar medium moving outward with a velocity of ~1000 km/s and created by relatively recent mass loss from the progenitor star. Ultraviolet through infrared observations show that the fading SN 2015G (which was probably discovered some 20 days post-peak) had a spectral energy distribution that was well described by a simple, single-component blackbody. Archival HST images provide upper limits on the luminosity of SN 2015G's progenitor, while nondetections of any luminous radio afterglow and optical nondetections of outbursts over the past two decades provide constraints upon its mass-loss history., Comment: 19 pages, 13 figures. Published in MNRAS

SweetSpot is a three-year National Optical Astronomy Observatory (NOAO) Survey program to observe Type Ia supernovae (SNe Ia) in the smooth Hubble flow with the WIYN High-resolution Infrared Camera (WHIRC) on the WIYN 3.5-m telescope. We here present data from the first half of this survey, covering the 2011B-2013B NOAO semesters, and consisting of 493 calibrated images of 74 SNe Ia observed in the rest-frame near-infrared (NIR) from $0.02 < z < 0.09$. Because many observed supernovae require host galaxy subtraction from templates taken in later semesters, this release contains only the 186 NIR ($JHK_s$) data points for the 33 SNe Ia that do not require host-galaxy subtraction. The sample includes 4 objects with coverage beginning before the epoch of B-band maximum and 27 beginning within 20 days of B-band maximum. We also provide photometric calibration between the WIYN+WHIRC and Two-Micron All Sky Survey (2MASS) systems along with light curves for 786 2MASS stars observed alongside the SNe Ia. This work is the first in a planned series of three SweetSpot Data Releases. Future releases will include the full set of images from all 3 years of the survey, including host-galaxy reference images and updated data processing and host-galaxy reference subtraction. SweetSpot will provide a well-calibrated sample that will help improve our ability to standardize distance measurements to SNe Ia, examine the intrinsic optical-NIR colors of SNe Ia at different epochs, explore nature of dust in other galaxies, and act as a stepping stone for more distant, potentially space-based surveys., Comment: Published in AJ. 10 tables. 11 figures. Lightcurve plots included as a figureset and available in source tarball. Data online at http://www.phyast.pitt.edu/~wmwv/SweetSpot/DR1_data/

We compare analytic predictions of supernova light curves with recent high quality data from SN2011fe (Ia), from KSN2011b (Ia), and the Palomar Transient Factory (PTF) and the La Silla-QUEST variability survey (LSQ) (Ia). Because of the steady, fast cadence of observations, KSN2011b provides unique new information on SNe Ia: the smoothness of the light curve, which is consistent with significant large-scale mixing during the explosion, possibly due to 3D effects (e.g., Rayleigh-Taylor instabilities), and provides support for a slowly-varying leakage (mean opacity). For a more complex light curve (SN2008D, SNIb), we separate the luminosity due to multiple causes and indicate the possibility of a radioactive plume. The early rise in luminosity is shown to be affected by the opacity (leakage rate) for thermal and non-thermal radiation. A general derivation of Arnett's rule again shows that it depends upon {\em all} processes heating the plasma, not just radioactive ones, so that SNe Ia will differ from SNe Ibc if the latter have multiple heating processes., Comment: 18 pages, 7 figures; accepted ApJ 2017; expanded comparative discussion of types Ia, b, and c; succinct presentation of light curve algorithm to help users catch implementation errors

The Arizona-NOAO Temporal Analysis and Response to Events System (ANTARES) is a joint effort of NOAO and the Department of Computer Science at the University of Arizona to build prototype software to process alerts from time-domain surveys, especially LSST, to identify those alerts that must be followed up immediately. Value is added by annotating incoming alerts with existing information from previous surveys and compilations across the electromagnetic spectrum and from the history of past alerts. Comparison against a knowledge repository of properties and features of known or predicted kinds of variable phenomena is used for categorization. The architecture and algorithms being employed are described.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

Herein we analyse late-time (post-plateau; 103 < t < 1229 d) optical spectra of low-redshift (z < 0.016), hydrogen-rich Type IIP supernovae (SNe IIP). Our newly constructed sample contains 91 nebular spectra of 38 SNe IIP, which is the largest dataset of its kind ever analysed in one study, and many of the objects have complementary photometric data. We determined the peak and total luminosity, velocity of the peak, HWHM intensity, and profile shape for many emission lines. Temporal evolution of these values and various flux ratios are studied. We also investigate the correlations between these measurements and photometric observables, such as the peak and plateau absolute magnitudes and the late-time light curve decline rates in various optical bands. The strongest and most robust result we find is that the luminosities of all spectral features (except those of helium) tend to be higher in objects with steeper late-time V-band decline rates. A steep late-time V-band slope likely arises from less efficient trapping of gamma-rays and positrons, which could be caused by multidimensional effects such as clumping of the ejecta or asphericity of the explosion itself. Furthermore, if gamma-rays and positrons can escape more easily, then so can photons via the observed emission lines, leading to more luminous spectral features. It is also shown that SNe IIP with larger progenitor stars have ejecta with a more physically extended oxygen layer that is well-mixed with the hydrogen layer. In addition, we find a subset of objects with evidence for asymmetric Ni-56 ejection, likely bipolar in shape. We also compare our observations to theoretical late-time spectral models of SNe IIP from two separate groups and find moderate-to-good agreement with both sets of models. Our SNe IIP spectra are consistent with models of 12-15 M_Sun progenitor stars having relatively low metallicity (Z $\le$ 0.01)., Comment: resubmitted to MNRAS, 50 pages (29 of which are tables), 15 figures, 10 tables

The Large Synoptic Survey Telescope (LSST) will be a discovery machine for the astronomy and physics communities, revealing astrophysical phenomena from the Solar System to the outer reaches of the observable Universe. While many discoveries will be made using LSST data alone, taking full scientific advantage of LSST will require ground-based optical-infrared (OIR) supporting capabilities, e.g., observing time on telescopes, instrumentation, computing resources, and other infrastructure. This community-based study identifies, from a science-driven perspective, capabilities that are needed to maximize LSST science. Expanding on the initial steps taken in the 2015 OIR System Report, the study takes a detailed, quantitative look at the capabilities needed to accomplish six representative LSST-enabled science programs that connect closely with scientific priorities from the 2010 decadal surveys. The study prioritizes the resources needed to accomplish the science programs and highlights ways that existing, planned, and future resources could be positioned to accomplish the science goals., Comment: 174 pages; one chapter of this report was previously published as arXiv:1607.04302

We re-examine the classifications of supernovae (SNe) presented in the Lick Observatory Supernova Search (LOSS) volume-limited sample with a focus on the stripped-envelope SNe. The LOSS volume-limited sample, presented by Leaman et al. (2011) and Li et al. (2011b), was calibrated to provide meaningful measurements of SN rates in the local universe; the results presented therein continue to be used for comparisons to theoretical and modeling efforts. Many of the objects from the LOSS sample were originally classified based upon only a small subset of the data now available, however, and recent studies have both updated some subtype distinctions and improved our ability to perform robust classications, especially for stripped-envelope SNe. We re-examine the spectroscopic classifications of all events in the LOSS volume-limited sample (180 SNe and SN impostors) and update them if necessary. We discuss the populations of rare objects in our sample including broad-lined Type Ic SNe, Ca-rich SNe, SN 1987A-like events (we identify SN 2005io as SN 1987A-like here for the first time), and peculiar subtypes. The relative fractions of Type Ia SNe, Type II SNe, and stripped-envelope SNe in the local universe are not affected, but those of some subtypes are. Most significantly, after discussing the often unclear boundary between SNe Ib and Ic when only noisy spectra are available, we find a higher SN Ib fraction and a lower SN Ic fraction than calculated by Li et al. (2011b): spectroscopically normal SNe Ib occur in the local universe $1.7 \pm 0.9$ times more often than do normal SNe Ic., Comment: Accepted for publication in PASP

We present initial results from a program aimed at establishing a network of hot DA white dwarfs to serve as spectrophotometric standards for present and future wide-field surveys. These stars span the equatorial zone and are faint enough to be conveniently observed throughout the year with large-aperture telescopes. Spectra of these white dwarfs are analyzed to generate a non-local-thermodynamic-equilibrium (NLTE) model atmosphere normalized to HST colors, including adjustments for wavelength-dependent interstellar extinction. Once established, this standard star network will serve ground-based observatories in both hemispheres as well as space-based instrumentation from the UV to the near IR. We demonstrate the effectiveness of this concept and show how two different approaches to the problem using somewhat different assumptions produce equivalent results. We discuss lessons learned and the resulting corrective actions applied to our program., Comment: (15 pages, 7 figures, 5 tables, accepted for publication in ApJ)

The ESSENCE survey discovered 213 Type Ia supernovae at redshifts 0.1 < z < 0.81 between 2002 and 2008. We present their R and I-band photometry, measured from images obtained using the MOSAIC II camera at the CTIO 4 m Blanco telescope, along with rapid-response spectroscopy for each object. We use our spectroscopic follow-up observations to determine an accurate, quantitative classification and a precise redshift. Through an extensive calibration program we have improved the precision of the CTIO Blanco natural photometric system. We use several empirical metrics to measure our internal photometric consistency and our absolute calibration of the survey. We assess the effect of various potential sources of systematic bias on our measured fluxes, and we estimate that the dominant term in the systematic error budget from the photometric calibration on our absolute fluxes is ~1%., Comment: (40 pages, 21 figures, accepted for publication in ApJS)

During the first few days after explosion, Type II supernovae (SNe) are dominated by relatively simple physics. Theoretical predictions regarding early-time SN light curves in the ultraviolet (UV) and optical bands are thus quite robust. We present, for the first time, a sample of $57$ $R$-band Type II SN light curves that are well monitored during their rise, having $>5$ detections during the first 10 days after discovery, and a well-constrained time of explosion to within $1-3$ days. We show that the energy per unit mass ($E/M$) can be deduced to roughly a factor of five by comparing early-time optical data to the model of Rabinak & Waxman (2011), while the progenitor radius cannot be determined based on $R$-band data alone. We find that Type II SN explosion energies span a range of $E/M=(0.2-20)\times 10^{51} \; \rm{erg/(10 M}_\odot$), and have a mean energy per unit mass of $\left\langle E/M \right\rangle = 0.85\times 10^{51} \; \rm{erg/(10 M}_\odot$), corrected for Malmquist bias. Assuming a small spread in progenitor masses, this indicates a large intrinsic diversity in explosion energy. Moreover, $E/M$ is positively correlated with the amount of $^{56}\rm{Ni}$ produced in the explosion, as predicted by some recent models of core-collapse SNe. We further present several empirical correlations. The peak magnitude is correlated with the decline rate ($\Delta m_{15}$), the decline rate is weakly correlated with the rise time, and the rise time is not significantly correlated with the peak magnitude. Faster declining SNe are more luminous and have longer rise times. This limits the possible power sources for such events., Comment: Accepted to ApJ

We present new Hubble Space Telescope (HST) multi-epoch ultraviolet (UV) spectra of the bright Type IIb SN 2013df, and undertake a comprehensive anal- ysis of the set of four Type IIb supernovae for which HST UV spectra are available (SN 1993J, SN 2001ig, SN 2011dh, and SN 2013df). We find strong diversity in both continuum levels and line features among these objects. We use radiative-transfer models that fit the optical part of the spectrum well, and find that in three of these four events we see a UV continuum flux excess, apparently unaffected by line absorption. We hypothesize that this emission originates above the photosphere, and is related to interaction with circumstel- lar material (CSM) located in close proximity to the SN progenitor. In contrast, the spectra of SN 2001ig are well fit by single-temperature models, display weak continuum and strong reverse-fluorescence features, and are similar to spectra of radioactive 56Ni-dominated Type Ia supernovae. A comparison of the early shock-cooling components in the observed light curves with the UV continuum levels which we assume trace the strength of CSM interaction suggests that events with slower cooling have stronger CSM emission. The radio emission from events having a prominent UV excess is perhaps consistent with slower blast-wave velocities, as expected if the explosion shock was slowed down by the CSM that is also responsible for the strong UV, but this connection is currently speculative as it is based on only a few events.

The discovery of variable and transient sources is an essential product of synoptic surveys. The alert stream will require filtering for personalized criteria -- a process managed by a functionality commonly described as a Broker. In order to understand quantitatively the magnitude of the alert generation and Broker tasks, we have undertaken an analysis of the most numerous types of variable targets in the sky -- Galactic stars, QSOs, AGNs and asteroids. It is found that LSST will be capable of discovering ~10^5 high latitude |b| > 20 deg) variable stars per night at the beginning of the survey. (The corresponding number for |b| < 20 deg is orders of magnitude larger, but subject to caveats concerning extinction and crowding.) However, the number of new discoveries may well drop below 100/night within 2 years. The same analysis applied to GAIA clarifies the complementarity of the GAIA and LSST surveys. Discovery of variable galactic nuclei (AGNs) and Quasi Stellar Objects (QSOs) are each predicted to begin at ~3000 per night, and decrease by 50X over 4 years. SNe are expected at ~1100/night, and after several survey years will dominate the new variable discovery rate. LSST asteroid discoveries will start at > 10^5 per night, and if orbital determination has a 50% success rate per epoch, will drop below 1000/night within 2 years., Comment: Published in ApJ, 22 pages, 26 figures, 3 tables; subsequent correction of an important error in interpretation of the Besancon models reduces stellar variable numbers significantly, as seen in figures 6-13

The Arizona-NOAO Temporal Analysis and Response to Events System (ANTARES) is a joint project of the National Optical Astronomy Observatory and the Department of Computer Science at the University of Arizona. The goal is to build the software infrastructure necessary to process and filter alerts produced by time-domain surveys, with the ultimate source of such alerts being the Large Synoptic Survey Telescope (LSST). The ANTARES broker will add value to alerts by annotating them with information from external sources such as previous surveys from across the electromagnetic spectrum. In addition, the temporal history of annotated alerts will provide further annotation for analysis. These alerts will go through a cascade of filters to select interesting candidates. For the prototype, `interesting' is defined as the rarest or most unusual alert, but future systems will accommodate multiple filtering goals. The system is designed to be flexible, allowing users to access the stream at multiple points throughout the process, and to insert custom filters where necessary. We describe the basic architecture of ANTARES and the principles that will guide development and implementation., Comment: 8 pages, 2 figures, submitted to Observatory Operations: Strategies, Processes and Systems V, (Proc. SPIE 9149), June 2014, Montreal

Since the discovery of the unusual prototype SN 2002cx, the eponymous class of low-velocity, hydrogen-poor supernovae has grown to include at most another two dozen members identified from several heterogeneous surveys, in some cases ambiguously. Here we present the results of a systematic study of 1077 hydrogen-poor supernovae discovered by the Palomar Transient Factory, leading to nine new members of this peculiar class. Moreover we find there are two distinct subclasses based on their spectroscopic, photometric, and host galaxy properties: The "SN 2002cx-like" supernovae tend to be in later-type or more irregular hosts, have more varied and generally dimmer luminosities, have longer rise times, and lack a Ti II trough when compared to the "SN 2002es-like" supernovae. None of our objects show helium, and we counter a previous claim of two such events. We also find that these transients comprise 5.6+17-3.7% (90% confidence) of all SNe Ia, lower compared to earlier estimates. Combining our objects with the literature sample, we propose that these subclasses have two distinct physical origins., Comment: 49 pages, 36 figures, submitted to ApJ