Multi-messenger astrophysics is a fast-growing, interdisciplinary field that combines data, which vary in volume and speed of data processing, from many different instruments that probe the Universe using different cosmic messengers: electromagnetic waves, cosmic rays, gravitational waves and neutrinos. In this Expert Recommendation, we review the key challenges of real-time observations of gravitational wave sources and their electromagnetic and astroparticle counterparts, and make a number of recommendations to maximize their potential for scientific discovery. These recommendations refer to the design of scalable and computationally efficient machine learning algorithms; the cyber-infrastructure to numerically simulate astrophysical sources, and to process and interpret multi-messenger astrophysics data; the management of gravitational wave detections to trigger real-time alerts for electromagnetic and astroparticle follow-ups; a vision to harness future developments of machine learning and cyber-infrastructure resources to cope with the big-data requirements; and the need to build a community of experts to realize the goals of multi-messenger astrophysics.

The centaur (2060) Chiron has exhibited outgassing behaviour and possibly hosts a ring system. On 2011 November 29, Chiron occulted a fairly bright star (R approximately 15 mag) as observed from the 3-m NASA Infrared Telescope Facility (IRTF) on Mauna Kea and the 2-m Faulkes Telescope North (FTN) at Haleakala. Data were taken as visible wavelength images and simultaneous, low-resolution, near-infrared (NIR) spectra. Here, we present a detailed examination of the light-curve features in the optical data and an analysis of the near-infrared spectra. We place a lower limit on the diameter of Chiron's nucleus of 160.2+/-1.3 km. Sharp, narrow dips were observed between 280-360 km from the centre (depending on event geometry). For a central chord and assumed ring plane, the separated features are 298.5 to 302 and 308 to 310.5 km from the nucleus, with normal optical depth approximately 0.5-0.9, and a gap of 9.1+/-1.3 km. These features are similar in equivalent depth to Chariklo's inner ring. The absence of absorbing or scattering material near the nucleus suggests that these sharp dips are more likely to be planar rings than a shell of material. The region of relatively-increased transmission is within the 1:2 spin-orbit resonance, which is consistent with the proposed clearing pattern for a non-axisymmetric nucleus. Characteristics of additional, azimuthally incomplete features are presented, which are likely to be transient, as well as detection of an extended shell or diffuse ring from approximately 900-1500 km. There are no significant features in the NIR light curves, nor any correlation between optical features and NIR spectral slope.

Fundamental coding and software development skills are increasingly necessary for success in nearly every aspect of astronomical and astrophysical research as large surveys and high resolution simulations become the norm. However, professional training in these skills is inaccessible or impractical for many members of our community. Students and professionals alike have been expected to acquire these skills on their own, apart from formal classroom curriculum or on-the-job training. Despite the recognized importance of these skills, there is little opportunity to develop them - even for interested researchers. To ensure a workforce capable of taking advantage of the computational resources and the large volumes of data coming in the next decade, we must identify and support ways to make software development training widely accessible to community members, regardless of affiliation or career level. To develop and sustain a technology capable astronomical and astrophysical workforce, we recommend that agencies make funding and other resources available in order to encourage, support and, in some cases, require progress on necessary training, infrastructure and policies. In this white paper, we focus on recommendations for how funding agencies can lead in the promotion of activities to support the astronomy and astrophysical workforce in the 2020s.

We describe an "Urban Observatory" facility designed for the study of complex urban systems via persistent, synoptic, and granular imaging of dynamical processes in cities. An initial deployment of the facility has been demonstrated in New York City and consists of a suite of imaging systems - both broadband and hyperspectral - sensitive to wavelengths from the visible (~400 nm) to the infrared (~13 micron) operating at cadences of ~0.01 - 30 Hz (characteristically ~0.1 Hz). Much like an astronomical survey, the facility generates a large imaging catalog from which we have extracted observables (e.g., time-dependent brightnesses, spectra, temperatures, chemical species, etc.), collecting them in a parallel source catalog. We have demonstrated that, in addition to the urban science of cities as systems, these data are applicable to a myriad of domain-specific scientific inquiries related to urban functioning including energy consumption and end use, environmental impacts of cities, and patterns of life and public health. We show that an Urban Observatory facility of this type has the potential to improve both a city's operations and the quality of life of its inhabitants.

Time domain and multi-messenger astrophysics are growing and important modes of observational astronomy that will help define astrophysics in the 2020s. Significant effort is being put into developing the components of a follow-up system for dynamically turning survey alerts into data. This system consists of: 1) brokers that will aggregate, classify, and filter alerts; 2) Target Observation Managers (TOMs) for prioritizing targets and managing observations and data; and 3) observatory interfaces, schedulers, and facilities along with data reduction software and science archives. These efforts need continued community support and funding in order to complete and maintain them. Many of the efforts can be community open-source software projects but they will benefit from the leadership of professional software developers. The coordination should be done by institutions that are involved in the follow-up system such as the national observatories (e.g. LSST/Gemini/NOAO Mid-scale/Community Science and Data Center) or a new MMA institute. These tools will help the community to produce the most science from new facilities and will provide new capabilities for all users of the facilities that adopt them.

In the current era of time-domain astronomy, it is increasingly important to have rigorous, data driven models for classifying transients, including supernovae. We present the first application of Principal Component Analysis to the spectra of stripped-envelope core-collapse supernovae. We use one of the largest compiled optical datasets of stripped-envelope supernovae, containing 160 SNe and 1551 spectra. We find that the first 5 principal components capture 79\% of the variance of our spectral sample, which contains the main families of stripped supernovae: Ib, IIb, Ic and broad-lined Ic. We develop a quantitative, data-driven classification method using a support vector machine, and explore stripped-envelope supernovae classification as a function of phase relative to V-band maximum light. Our classification method naturally identifies "transition" supernovae and supernovae with contested labels, which we discuss in detail. We find that the stripped-envelope supernovae types are most distinguishable in the later phase ranges of $10\pm5$ days and $15\pm5$ days relative to V-band maximum, and we discuss the implications of our findings for current and future surveys such as ZTF and LSST.

We present results from three world-wide campaigns that resulted in the detections of two single-chord and one multi-chord stellar occultations by the Plutino object (84922) 2003~VS$_2$. From the single-chord occultations in 2013 and 2014 we obtained accurate astrometric positions for the object, while from the multi-chord occultation on November 7th, 2014, we obtained the parameters of the best-fitting ellipse to the limb of the body at the time of occultation. We also obtained short-term photometry data for the body in order to derive its rotational phase during the occultation. The rotational light curve present a peak-to-peak amplitude of 0.141 $\pm$ 0.009 mag. This allows us to reconstruct the three-dimensional shape of the body, with principal semi-axes $a = 313.8 \pm 7.1$ km, $b = 265.5^{+8.8}_{-9.8}$ km, and $c = 247.3^{+26.6}_{-43.6}$ km, which is not consistent with a Jacobi triaxial equilibrium figure. The derived spherical volume equivalent diameter of $548.3 ^{+29.5}_{-44.6}$ km is about 5\% larger than the radiometric diameter of 2003~VS$_2$ derived from Herschel data of $523 \pm 35$ km, but still compatible with it within error bars. From those results we can also derive the geometric albedo ($0.123 ^{+0.015}_{-0.014}$) and, under the assumption that the object is a Maclaurin spheroid, the density $\rho = 1400^{+1000}_{-300}$ for the plutino. The disappearances and reappearances of the star during the occultations do not show any compelling evidence for a global atmosphere considering a pressure upper limit of about 1 microbar for a pure nitrogen atmosphere, nor secondary features (e.g. rings or satellite) around the main body.

The Large Synoptic Survey Telescope (LSST) can advance scientific frontiers beyond its groundbreaking 10-year survey. Here we explore opportunities for extended operations with proposal-based observing strategies, new filters, or transformed instrumentation. We recommend the development of a mid-decade community- and science-driven process to define next-generation LSST capabilities.

Through the lens of the LSST Science Collaborations' experience, this paper advocates for new and improved ways to fund large, complex collaborations at the interface of data science and astrophysics as they work in preparation for and on peta-scale, complex surveys, of which LSST is a prime example. We advocate for the establishment of programs to support both research and infrastructure development that enables innovative collaborative research on such scales.

We identify minimal observing cadence requirements that enable photometric astronomical surveys to detect and recognize fast and explosive transients and fast transient features. Observations in two different filters within a short time window (e.g., g-and-i, or r-and-z, within < 0.5 hr) and a repeat of one of those filters with a longer time window (e.g., > 1.5 hr) are desirable for this purpose. Such an observing strategy delivers both the color and light curve evolution of transients on the same night. This allows the identification and initial characterization of fast transient -- or fast features of longer timescale transients -- such as rapidly declining supernovae, kilonovae, and the signatures of SN ejecta interacting with binary companion stars or circumstellar material. Some of these extragalactic transients are intrinsically rare and generally all hard to find, thus upcoming surveys like the Large Synoptic Survey Telescope (LSST) could dramatically improve our understanding of their origin and properties. We colloquially refer to such a strategy implementation for the LSST as the Presto-Color strategy (rapid-color). This cadence's minimal requirements allow for overall optimization of a survey for other science goals.

We present a cadence optimization strategy to unveil a large population of kilonovae using optical imaging alone. These transients are generated during binary neutron star and potentially neutron star-black hole mergers and are electromagnetic counterparts to gravitational-wave signals detectable in nearby events with Advanced LIGO, Advanced Virgo, and other interferometers that will come online in the near future. Discovering a large population of kilonovae will allow us to determine how heavy element production varies with the intrinsic parameters of the merger and across cosmic time. The rate of binary neutron star mergers is still uncertain, but only few (less than 15) events with associated kilonovae may be detectable per year within the horizon of next-generation ground-based interferometers. The rapid evolution (hours to days) at optical/infrared wavelengths, relatively low luminosity, and the low volumetric rate of kilonovae makes their discovery difficult, especially during blind surveys of the sky. We propose future large surveys to adopt a rolling cadence in which g-i observations are taken nightly for blocks of 10 consecutive nights. With the current baseline2018a cadence designed for the Large Synoptic Survey Telescope (LSST), less than 7.5 poorly-sampled kilonovae are expected to be detected in both the Wide Fast Deep (WFD) and Deep Drilling Fields (DDF) surveys per year, under optimistic assumptions on their rate, duration, and luminosity. We estimate the proposed strategy to return up to about 272 GW170817-like kilonovae throughout the LSST WFD survey, discovered independently from gravitational-wave triggers.

Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.)

The identification of the electromagnetic counterpart of the gravitational wave event, GW170817, and discovery of neutrinos and gamma-rays from TXS 0506+056 heralded the new era of multi-messenger astrophysics. As the number of multi-messenger events rapidly grow over the next decade, the cyberinfrastructure requirements to handle the increase in data rates, data volume, need for event follow up, and analysis across the different messengers will also explosively grow. The cyberinfrastructure requirements to enhance multi-messenger astrophysics will both be a major challenge and opportunity for astronomers, physicists, computer scientists and cyberinfrastructure specialists. Here we outline some of these requirements and argue for a distributed cyberinfrastructure institute for multi-messenger astrophysics to meet these challenges.

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.

Context. Pluto's tenuous nitrogen (N2) atmosphere undergoes strong seasonal effects due to high obliquity and orbital eccentricity, and has been recently (July 2015) observed by the New Horizons spacecraft. Goals are (i) construct a well calibrated record of the seasonal evolution of surface pressure on Pluto and (ii) constrain the structure of the lower atmosphere using a central flash observed in 2015. Method: eleven stellar occultations by Pluto observed between 2002 and 2016 are used to retrieve atmospheric profiles (density, pressure, temperature) between $\sim$5 km and $\sim$380 km altitude levels (i.e. pressures from about 10 microbar to 10 nanobar). Results: (i) Pressure has suffered a monotonic increase from 1988 to 2016, that is compared to a seasonal volatile transport model, from which tight constraints on a combination of albedo and emissivity of N2 ice are derived; (ii) A central flash observed on 2015 June 29 is consistent with New Horizons REX profiles, provided that (a) large diurnal temperature variations (not expected by current models) occur over Sputnik Planitia and/or (b) hazes with tangential optical depth of about 0.3 are present at 4-7 km altitude levels and/or (c) the nominal REX density values are overestimated by an implausibly large factor of about 20% and/or (d) higher terrains block part of the flash in the Charon facing hemisphere.

This report provides an overview of recent work that harnesses the Big Data Revolution and Large Scale Computing to address grand computational challenges in Multi-Messenger Astrophysics, with a particular emphasis on real-time discovery campaigns. Acknowledging the transdisciplinary nature of Multi-Messenger Astrophysics, this document has been prepared by members of the physics, astronomy, computer science, data science, software and cyberinfrastructure communities who attended the NSF-, DOE- and NVIDIA-funded "Deep Learning for Multi-Messenger Astrophysics: Real-time Discovery at Scale" workshop, hosted at the National Center for Supercomputing Applications, October 17-19, 2018. Highlights of this report include unanimous agreement that it is critical to accelerate the development and deployment of novel, signal-processing algorithms that use the synergy between artificial intelligence (AI) and high performance computing to maximize the potential for scientific discovery with Multi-Messenger Astrophysics. We discuss key aspects to realize this endeavor, namely (i) the design and exploitation of scalable and computationally efficient AI algorithms for Multi-Messenger Astrophysics; (ii) cyberinfrastructure requirements to numerically simulate astrophysical sources, and to process and interpret Multi-Messenger Astrophysics data; (iii) management of gravitational wave detections and triggers to enable electromagnetic and astro-particle follow-ups; (iv) a vision to harness future developments of machine and deep learning and cyberinfrastructure resources to cope with the scale of discovery in the Big Data Era; (v) and the need to build a community that brings domain experts together with data scientists on equal footing to maximize and accelerate discovery in the nascent field of Multi-Messenger Astrophysics.

Stars exhibit a bewildering variety of variable behaviors ranging from explosive magnetic flares to stochastically changing accretion to periodic pulsations or rotations. The principal LSST surveys will have cadences too sparse and irregular to capture most of these phenomena. A novel idea is proposed here to observe a single Galactic field, rich in unobscured stars, in a continuous sequence of $\sim 15$ second exposures for one long winter night in a single photometric band. The result will be a unique dataset of $\sim 1$ million regularly spaced stellar lightcurves. The lightcurves will gives a particularly comprehensive collection of dM star variability. A powerful array of statistical procedures can be applied to the ensemble of lightcurves from the long-standing fields of time series analysis, signal processing and econometrics. Dozens of `features' describing the variability can be extracted and subject to machine learning classification, giving a unique authoritative objective classification of rapidly variable stars. The most effective features can then inform the wider LSST community on the best approaches to variable star identification and classification from the sparse, irregular cadences that dominate the LSST project.

Unlike the ordinary supernovae (SNe) some of which are hydrogen and helium deficient (called Type Ic SNe), broad-lined Type Ic SNe (SNe Ic-bl) are very energetic events, and all SNe coincident with bona fide long duration gamma-ray bursts (LGRBs) are of Type Ic-bl. Understanding the progenitors and the mechanism driving SN Ic-bl explosions vs those of their SNe Ic cousins is key to understanding the SN-GRB relationship and jet production in massive stars. Here we present the largest set of host-galaxy spectra of 28 SNe Ic and 14 SN Ic-bl, all discovered before 2013 by the same untargeted survey, namely the Palomar Transient Factory (PTF). We carefully measure their gas-phase metallicities, stellar masses (M*s) and star-formation rates (SFRs) by taking into account recent progress in the metallicity field and propagating uncertainties correctly. We further re-analyze the hosts of 10 literature SN-GRBs using the same methods and compare them to our PTF SN hosts with the goal of constraining their progenitors from their local environments by conducting a thorough statistical comparison, including upper limits. We find that the metallicities, SFRs and M*s of our PTF SN Ic-bl hosts are statistically comparable to those of SN-GRBs, but significantly lower than those of the PTF SNe Ic. The mass-metallicity relations as defined by the SNe Ic-bl and SN-GRBs are not significantly different from the same relations as defined by the SDSS galaxies, in contrast to claims by earlier works. Our findings point towards low metallicity as a crucial ingredient for SN Ic-bl and SN-GRB production since we are able to break the degeneracy between high SFR and low metallicity. We suggest that the PTF SNe Ic-bl may have produced jets that were choked inside the star or were able break out of the star as unseen low-luminosity or off-axis GRBs.

The discovery of the electromagnetic counterparts to the binary neutron star merger GW170817 has opened the era of GW+EM multi-messenger astronomy. Exploiting this breakthrough requires increasing samples to explore the diversity of kilonova behaviour and provide more stringent constraints on the Hubble constant, and tests of fundamental physics. LSST can play a key role in this field in the 2020s, when the gravitational wave detector network is expected to detect higher rates of merger events involving neutron stars ($\sim$10s per year) out to distances of several hundred Mpc. Here we propose comprehensive target-of-opportunity (ToOs) strategies for follow-up of gravitational-wave sources that will make LSST the premiere machine for discovery and early characterization for neutron star mergers and other gravitational-wave sources.

We present a unique method that allows the LSST to scan the sky for stellar variability on short timescales. The operational component of the strategy requires LSST to take star trail images. The image processing component uses deep learning to sift for transient events on timescales down to 10 ms. We advocate for enabling this observing mode with LSST, as coupling this capability with the LSST's tremendous 319.5 m$^2$deg$^2$ etendue will produce the first wide area optical survey of the universe on these timescales. We explain how these data will advance both planned lines of investigation and enable new research in the areas of stellar flares, cataclysmic variables, active galactic nuclei, Kuiper Belt objects, gamma-ray bursts, and fast radio bursts.

The Wide-Field Infrared Survey Telescope (WFIRST) is expected to launch in the mid-2020s. With its wide-field near-infrared (NIR) camera, it will survey the sky to unprecedented detail. As part of normal operations and as the result of multiple expected dedicated surveys, WFIRST will produce several relatively wide-field (tens of square degrees) deep (limiting magnitude of 28 or fainter) fields. In particular, a planned supernova survey is expected to image 3 deep fields in the LSST footprint roughly every 5 days over 2 years. Stacking all data, this survey will produce, over all WFIRST supernova fields in the LSST footprint, ~12-25 deg^2 and ~5-15 deg^2 regions to depths of ~28 mag and ~29 mag, respectively. We suggest LSST undertake mini-surveys that will match the WFIRST cadence and simultaneously observe the supernova survey fields during the 2-year WFIRST supernova survey, achieving a stacked depth similar to that of the WFIRST data. We also suggest additional observations of these same regions throughout the LSST survey to get deep images earlier, have long-term monitoring in the fields, and produce deeper images overall. These fields will provide a legacy for cosmology, extragalactic, and transient/variable science.

Empirical models of supernova (SN) spectral energy distributions (SEDs) are widely used for SN survey simulations and photometric classifications. The existing library of SED models has excellent optical templates but limited, poorly constrained coverage of ultraviolet (UV) and infrared (IR) wavelengths. However, both regimes are critical for the design and operation of future SN surveys, particularly at IR wavelengths that will be accessible with the James Webb Space Telescope (JWST) and the Wide-Field Infrared Survey Telescope (WFIRST). We create a public repository of improved empirical SED templates using a sampling of Type Ia and core-collapse (CC) photometric light curves to extend the Type Ia parameterized SALT2 model and a set of SN Ib, SN Ic, and SN II SED templates into the UV and near-IR. We apply this new repository of extrapolated SN SED models to examine how future surveys can discriminate between CC and Type Ia SNe at UV and IR wavelengths, and present an open-source software package written in Python, SNSEDextend, that enables a user to generate their own extrapolated SEDs.

We present multi-epoch photometry and spectroscopy of a light echo from eta Carinae's 19th century Great Eruption. This echo shows a steady decline over a decade, sampling the 1850s plateau of the eruption. Spectra show the bulk outflow speed increasing from 150 km/s at early times, up to 600 km/s in the plateau. Later phases also develop remarkably broad emission wings indicating mass accelerated to more than 10,000 km/s. Together with other clues, this provides direct evidence for an explosive ejection. This is accompanied by a transition from narrow absorption lines to emission lines, often with broad P Cygni profiles. These changes imply that the pre-1845 luminosity spikes are distinct from the 1850s plateau. The key reason for this change may be that shock interaction dominates the plateau. The spectral evolution of eta Car closely resembles that of UGC2773-OT, which had clear signatures of shock interaction. We propose a 2-stage scenario for eta Car's eruption: (1) a slow outflow in the decades before the eruption, driven by binary interaction that produced a dense equatorial outflow, followed by (2) explosive energy injection that drove CSM interaction, powering the plateau and sweeping slower CSM into a fast shell that became the Homunculus. We discuss how this sequence could arise from a stellar merger in a triple system, leaving the eccentric binary seen today. This gives a self-consistent scenario that may explain interacting transients across a wide range of initial mass.

In our ongoing study of eta Carinae's light echoes, there is a relatively bright echo that has been fading slowly, reflecting the 1845-1858 plateau of the eruption. A separate paper discusses its detailed evolution, but here we highlight one important result: the H-alpha line shows extremely broad emission wings that reach -10,000km/s to the blue and +20,000km/s to the red. The line profile shape is inconsistent with electron scattering wings, indicating high-velocity outflowing material. These are the fastest outflow speeds ever seen in a non-terminal massive star eruption. The broad wings are absent in early phases of the eruption, but strengthen in the 1850s. These speeds are two orders of magnitude faster than the escape speed from a warm supergiant, and 5-10 times faster than winds from O-type or Wolf-Rayet stars. Instead, they are reminiscent of fast supernova ejecta or outflows from accreting compact objects, profoundly impacting our understanding of eta Car and related transients. This echo views eta Car from latitudes near the equator, so the high speed does not trace a collimated polar jet aligned with the Homunculus. Combined with fast material in the Outer Ejecta, it indicates a wide-angle explosive outflow. The fast material may constitute a small fraction of the total outflowing mass, most of which expands at 600 km/s. This is reminiscent of fast material revealed by broad absorption during the presupernova eruptions of SN2009ip.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5-1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift $\sim$10, signatures of Pop III stars, sources and physics of re-ionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA).

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

Unusually high velocities (< ~0.1c) and correspondingly high kinetic energies have been observed in a subset of Type Ic supernovae (so-called "broad-lined Ic" supernovae; SNe Ic-BL), prompting a search for a central engine model capable of generating such energetic explosions. A clue to the explosion mechanism may lie in the fact that all supernovae that accompany long-duration gamma-ray bursts belong to the SN Ic-BL class. Using a combination of two-dimensional relativistic hydrodynamics and radiation transport calculations, we demonstrate that the central engine responsible for long gamma-ray bursts can also trigger a SN Ic-BL. We find that a reasonable gamma-ray burst engine injected into a stripped Wolf-Rayet progenitor produces a relativistic jet with energy ~10^51 ergs, as well as a SN whose synthetic light curves and spectra are fully consistent with observed SNe Ic-BL during the photospheric phase. As a result of the jet's asymmetric energy injection, the SN spectra and light curves depend on viewing angle. The impact of viewing angle on the spectrum is particularly pronounced at early times, while the viewing angle dependence for the light curves (~10% variation in bolometric luminosity) persists throughout the photospheric phase.

Super-luminous supernovae (SLSNe) are tremendously luminous explosions whose power sources and progenitors are highly debated. Broad-lined SNe Ic (SNe Ic-bl) are the only type of SNe that are connected with long-duration gamma ray bursts (GRBs). Studying the spectral similarity and difference between the populations of hydrogen-poor SLSNe (SLSNe Ic) and of hydrogen-poor stripped-envelope core-collapse SNe, in particular SNe Ic and SNe Ic-bl, can provide crucial observations to test predictions of theories based on various power source models and progenitor models. In this paper, we collected all of the published optical spectra of 32 SLSNe Ic, 21 SNe Ic-bl, as well as 17 SNe Ic, quantified their spectral features, constructed average spectra, and compared them in a systematic way using new tools we have developed. We find that SLSNe Ic and SNe Ic-bl, including those connected with GRBs, have comparable widths for their spectral features and average absorption velocities at all phases. Thus, our findings strengthen the connection between SLSNe Ic and GRBs. In particular, SLSNe Ic have average \FeII absorption velocities of $-15,000 +/- 2,600$ \kms~at 10 days after peak, which are higher than those of SNe Ic by $\sim7,000$ \kms~on average. SLSNe Ic also have significantly broader \FeII lines than SNe Ic. Moreover, we find that such high absorption and width velocities of SLSNe Ic may be hard to explain by the interaction model, and none of 13 SLSNe Ic with measured absorption velocities spanning over 10 days has a convincing flat velocity-evolution, which is inconsistent with the magnetar model in one dimension. Lastly, we compare SN 2011kl, the first SN connected with an ultra-long GRB, with the mean spectrum of SLSNe Ic and of SNe Ic-bl.

In Paper I of this series, we showed that the ratio between stripped-envelope (SE) supernova (SN) and Type II SN rates reveals a significant SE SN deficiency in galaxies with stellar masses $\lesssim 10^{10}~{\rm M}_\odot$. Here, we test this result by splitting the volume-limited subsample of the Lick Observatory Supernova Search (LOSS) SN sample into low- and high-mass galaxies and comparing the relative rates of various SN types found in them. The LOSS volume-limited sample contains 180 SNe and SN impostors and is complete for SNe Ia out to 80 Mpc and core-collapse SNe out to 60 Mpc. All of these transients were recently reclassified by us in Shivvers et al. (2017). We find that the relative rates of some types of SNe differ between low- and high-mass galaxies: SNe Ib and Ic are underrepresented by a factor of ~3 in low-mass galaxies. These galaxies also contain the only examples of SN 1987A-like SNe in the sample and host about 9 times as many SN impostors. Normal SNe Ia seem to be ~30\% more common in low-mass galaxies, making these galaxies better sources for homogeneous SN Ia cosmology samples. The relative rates of SNe IIb are consistent in both low- and high-mass galaxies. The same is true for broad-line SNe Ic, although our sample includes only two such objects. The results presented here are in tension with a similar analysis from the Palomar Transient Factory, especially as regards SNe IIb.

Most types of supernovae (SNe) have yet to be connected with their progenitor stellar systems. Here, we reanalyze the 10-year SN sample collected during 1998-2008 by the Lick Observatory Supernova Search (LOSS) in order to constrain the progenitors of SNe Ia and stripped-envelope SNe (SE SNe, i.e., SNe~IIb, Ib, Ic, and broad-lined Ic). We matched the LOSS galaxy sample with spectroscopy from the Sloan Digital Sky Survey and measured SN rates as a function of galaxy stellar mass, specific star formation rate, and oxygen abundance (metallicity). We find significant correlations between the SN rates and all three galaxy properties. The SN Ia correlations are consistent with other measurements, as well as with our previous explanation of these measurements in the form of a combination of the SN Ia delay-time distribution and the correlation between galaxy mass and age. The ratio between the SE SN and SN II rates declines significantly in low-mass galaxies. This rules out single stars as SE~SN progenitors, and is consistent with predictions from binary-system progenitor models. Using well-known galaxy scaling relations, any correlation between the rates and one of the galaxy properties examined here can be expressed as a correlation with the other two. These redundant correlations preclude us from establishing causality - that is, from ascertaining which of the galaxy properties (or their combination) is the physical driver for the difference between the SE SN and SN II rates. We outline several methods that have the potential to overcome this problem in future works.

Hypertemporal visible imaging of an urban lightscape can reveal the phase of the electrical grid granular to individual housing units. In contrast to in-situ monitoring or metering, this method offers broad, persistent, real-time, and non-permissive coverage through a single camera sited at an urban vantage point. Rapid changes in the phase of individual housing units signal changes in load (e.g., appliances turning on and off), while slower building- or neighborhood-level changes can indicate the health of distribution transformers. We demonstrate the concept by observing the 120 Hz flicker of lights across a NYC skyline. A liquid crystal shutter driven at 119.75 Hz down-converts the flicker to 0.25 Hz, which is imaged at a 4 Hz cadence by an inexpensive CCD camera; the grid phase of each source is determined by analysis of its sinusoidal light curve over an imaging "burst" of some 25 seconds. Analysis of bursts taken at ~15 minute cadence over several hours demonstrates both the stability and variation of phases of halogen, incandescent, and some fluorescent lights. Correlation of such results with ground-truth data will validate a method that could be applied to better monitor electricity consumption and distribution in both developed and developing cities.

We compare the diversity of spectral line velocities in a large sample of type IIb supernovae (SNe IIb) with the expected asphericity in the explosion, as measured from the light echoes of Cassiopeia A (Cas A), which was a historical galactic SN IIb. We revisit the results of Rest et al. (2011a), who used light echoes to observe Cas A from multiple lines of sight and hence determine its asphericity, as seen in the velocity of three spectral lines (He I $\lambda$5876, H$\alpha$ and the Ca II NIR triplet). We confirm and improve on this measurement by reproducing the effect of the light echoes in the spectra of several extragalactic SNe IIb found in the literature as well as mean SN IIb spectra recently created by Liu et al. (2016), and comparing these to the observed light echo spectra of Cas A, including their associated uncertainties. In order to quantify the accuracy of this comparison we smooth the light echo spectra of Cas A using Gaussian processes and use a Monte Carlo method to measure the absorption velocities of these three features in the spectra. We then test the hypothesis that the diversity of ejecta velocities seen in SNe IIb can be explained by asphericity. We do this by comparing the range of velocities seen in the different light echoes, and hence different lines of sight, of Cas A to that seen in the population of SNe IIb. We conclude that these two ranges are of the same order and thus asphericity could be enough to explain the diversity in the expansion velocity alone.

Using the largest spectroscopic dataset of stripped-envelope core-collapse supernovae (stripped SNe), we present a systematic investigation of spectral properties of Type IIb SNe (SNe IIb), Type Ib SNe (SNe Ib), and Type Ic SNe (SNe Ic). Prior studies have been based on individual objects or small samples. Here, we analyze 227 spectra of 14 SNe IIb, 258 spectra of 21 SNe Ib, and 207 spectra of 17 SNe Ic based on the stripped SN dataset of Modjaz et al. (2014) and other published spectra of individual SNe. Each SN in our sample has a secure spectroscopic ID, a date of $V$-band maximum light, and multiple spectra at different phases. We analyze these spectra as a function of subtype and phase in order to improve the SN identification scheme and constrain the progenitors of different kinds of stripped SNe. By comparing spectra of SNe IIb with those of SNe Ib, we find that the strength of H$\alpha$ can be used to quantitatively differentiate between these two subtypes at all epochs. Moreover, we find a continuum in observational properties between SNe IIb and Ib. We address the question of hidden He in SNe Ic by comparing our observations with predictions from various models that either include hidden He or in which He has been burnt. Our results favor the He-free progenitor models for SNe Ic. Finally, we construct continuum-divided average spectra as a function of subtype and phase to quantify the spectral diversity of the different types of stripped SNe.

We present the first systematic investigation of spectral properties of 17 Type Ic Supernovae (SNe Ic), 10 broad-lined SNe Ic (SNe Ic-bl) without observed Gamma-Ray Bursts (GRBs) and 11 SNe Ic-bl with GRBs (SN-GRBs) as a function of time in order to probe their explosion conditions and progenitors. We analyze a total of 407 spectra, which were drawn from published spectra of individual SNe as well as from the densely time-sampled spectra data of Modjaz et al. (2014). In order to quantify the diversity of the SN spectra as a function of SN subtype, we construct average spectra of SNe Ic, SNe Ic-bl without GRBs and SNe Ic-bl with GRBs. We find that SN 1994I is not a typical SN Ic, in contrast to common belief, while the spectra of SN 1998bw/GRB 980425 are representative of mean spectra of SNe Ic-bl. We measure the ejecta absorption and width velocities using a new method described here and find that SNe Ic-bl with GRBs, on average, have quantifiably higher absorption velocities, as well as broader line widths than SNe without observed GRBs. In addition, we search for correlations between SN-GRB spectral properties and the energies of their accompanying GRBs. Finally, we show that the absence of clear He lines in optical spectra of SNe Ic-bl, and in particular of SN-GRBs, is not due to them being too smeared out due to the high velocities present in the ejecta. This implies that the progenitor stars of SN-GRBs are probably He-free, in addition to being H-free, which puts strong constraints on the stellar evolutionary paths needed to produce such SN-GRB progenitors at the observed low metallicities.

How black holes accrete surrounding matter is a fundamental, yet unsolved question in astrophysics. It is generally believed that matter is absorbed into black holes via accretion disks, the state of which depends primarily on the mass-accretion rate. When this rate approaches the critical rate (the Eddington limit), thermal instability is supposed to occur in the inner disc, causing repetitive patterns of large-amplitude X-ray variability (oscillations) on timescales of minutes to hours. In fact, such oscillations have been observed only in sources with a high mass accretion rate, such as GRS 1915+105. These large-amplitude, relatively slow timescale, phenomena are thought to have physical origins distinct from X-ray or optical variations with small amplitudes and fast ($\lesssim$10 sec) timescales often observed in other black hole binaries (e.g., XTE J1118+480 and GX 339-4). Here we report an extensive multi-colour optical photometric data set of V404 Cygni, an X-ray transient source containing a black hole of nine solar masses (and a conpanion star) at a distance of 2.4 kiloparsecs. Our data show that optical oscillations on timescales of 100 seconds to 2.5 hours can occur at mass-accretion rates more than ten times lower than previously thought. This suggests that the accretion rate is not the critical parameter for inducing inner-disc instabilities. Instead, we propose that a long orbital period is a key condition for these large-amplitude oscillations, because the outer part of the large disc in binaries with long orbital periods will have surface densities too low to maintain sustained mass accretion to the inner part of the disc. The lack of sustained accretion -- not the actual rate -- would then be the critical factor causing large-amplitude oscillations in long-period systems.

We present the open-source Python code pyMCZ that determines oxygen abundance and its distribution from strong emission lines in the standard metallicity calibrators, based on the original IDL code of Kewley & Dopita (2002) with updates from Kewley & Ellison (2008), and expanded to include more recently developed calibrators. The standard strong-line diagnostics have been used to estimate the oxygen abundance in the interstellar medium through various emission line ratios in many areas of astrophysics, including galaxy evolution and supernova host galaxy studies. We introduce a Python implementation of these methods that, through Monte Carlo sampling, better characterizes the statistical oxygen abundance confidence region including the effect due to the propagation of observational uncertainties. These uncertainties are likely to dominate the error budget in the case of distant galaxies, hosts of cosmic explosions. Given line flux measurements and their uncertainties, our code produces synthetic distributions for the oxygen abundance in up to 15 metallicity calibrators simultaneously, as well as for E(B-V), and estimates their median values and their 68% confidence regions. We test our code on emission line measurements from a sample of nearby supernova host galaxies (z < 0.15) and compare our metallicity results with those from previous methods. Our metallicity estimates are consistent with previous methods but yield smaller statistical uncertainties. Systematic uncertainties are not taken into account. We offer visualization tools to assess the spread of the oxygen abundance in the different calibrators, as well as the shape of the estimated oxygen abundance distribution in each calibrator, and develop robust metrics for determining the appropriate Monte Carlo sample size. The code is open access and open source and can be found at https://github.com/nyusngroup/pyMCZ (Abridged)

SN 2012dn is a super-Chandrasekhar mass candidate in a purportedly normal spiral (SAcd) galaxy, and poses a challenge for theories of type Ia supernova diversity. Here we utilize the fast and highly parameterized spectrum synthesis tool, SYNAPPS, to estimate relative expansion velocities of species inferred from optical spectra obtained with six facilities. As with previous studies of normal SN Ia, we find that both unburned carbon and intermediate mass elements are spatially coincident within the ejecta near and below 14,000 km/s. Although the upper limit on SN 2012dn's peak luminosity is comparable to some of the most luminous normal SN Ia, we find a progenitor mass exceeding ~1.6 Msun is not strongly favored by leading merger models since these models do not accurately predict spectroscopic observations of SN 2012dn and more normal events. In addition, a comparison of light curves and host-galaxy masses for a sample of literature and Palomar Transient Factory SN Ia reveals a diverse distribution of SN Ia subtypes where carbon-rich material remains unburned in some instances. Such events include SN 1991T, 1997br, and 1999aa where trace signatures of C III at optical wavelengths are presumably detected.

The progenitors of stripped-envelope supernovae (SNe Ibc) remain to be conclsuively identified, but correlations between SN rates and host-galaxy properties can constrain progenitor models. Here, we present one result from a re-analysis of the rates from the Lick Observatory Supernova Search. Galaxies with stellar masses less than $\sim 10^{10}~{\rm M_\odot}$ are less efficient at producing SNe Ibc than more massive galaxies. Any progenitor scenario must seek to explain this new observation.

Using a method to discover and classify supernovae (SNe) in galaxy spectra, we detect 91 Type Ia SNe (SNe Ia) and 16 Type II SNe (SNe II) among 740,000 galaxies of all types and 215,000 star-forming galaxies without active galactic nuclei, respectively, in Data Release 9 of the Sloan Digital Sky Survey. Of these SNe, 15 SNe Ia and 8 SNe II are new discoveries reported here for the first time. We use our SN samples to measure SN rates per unit mass as a function of galaxy stellar mass, star-formation rate (SFR), and specific SFR (sSFR), as derived by the MPA-JHU Galspec pipeline. We show that correlations between SN Ia and SN II rates per unit mass and galaxy stellar mass, SFR, and sSFR can be explained by a combination of the respective SN delay-time distributions (the distributions of times that elapse between the formation of a stellar population and all ensuing SNe), the ages of the surveyed galaxies, the redshifts at which they are observed, and their star-formation histories. This model was first suggested by Kistler et al. for the SN Ia rate-mass correlation, but is expanded here to SNe II and to correlations with galaxy SFR and sSFR. Finally, we measure a volumetric SN II rate at redshift 0.075 of $R_{\rm II,V} = 0.621^{+0.197}_{-0.154}({\rm stat})^{+0.024}_{-0.063}({\rm sys}) \times 10^{-4}$ yr$^{-1}$ Mpc$^{-3}$. Assuming that SNe IIP and IIL account for 60 per cent of all core-collapse (CC) SNe, the CC SN rate is $R_{\rm CC,V} = 1.04^{+0.33}_{-0.26}({\rm stat})^{+0.04}_{-0.11}({\rm sys}) \times 10^{-4}$ yr$^{-1}$ Mpc$^{-3}$.

For over a century, light echoes have been observed around variable stars and transients. The discovery of centuries-old light echoes from supernovae in the Large Magellanic Cloud has allowed the spectroscopic characterization of these events using modern instrumentation, even in the complete absence of any visual record of those events. Here we review the pivotal role the Blanco 4m telescope played in these discoveries.

We present a densely-sampled, homogeneous set of light curves of 64 low redshift (z < 0.05) stripped-envelope supernovae (SN of type IIb, Ib, Ic and Ic-bl). These data were obtained between 2001 and 2009 at the Fred L. Whipple Observatory (FLWO) on Mt. Hopkins in Arizona, with the optical FLWO 1.2-m and the near-infrared PAIRITEL 1.3-m telescopes. Our dataset consists of 4543 optical photometric measurements on 61 SN, including a combination of UBVRI, UBVr'i', and u'BVr'i', and 2142 JHKs near-infrared measurements on 25 SN. This sample constitutes the most extensive multi-color data set of stripped-envelope SN to date. Our photometry is based on template-subtracted images to eliminate any potential host galaxy light contamination. This work presents these photometric data, compares them with data in the literature, and estimates basic statistical quantities: date of maximum, color, and photometric properties. We identify promising color trends that may permit the identification of stripped-envelope SN subtypes from their photometry alone. Many of these SN were observed spectroscopically by the CfA SN group, and the spectra are presented in a companion paper (Modjaz et al. 2014). A thorough exploration that combines the CfA photometry and spectroscopy of stripped-envelope core-collapse SN will be presented in a follow-up paper.

We present 645 optical spectra of 73 supernovae (SNe) of Types IIb, Ib, Ic, and broad-lined Ic. All of these types are attributed to the core collapse of massive stars, with varying degrees of intact H and He envelopes before explosion. The SNe in our sample have a mean redshift <cz> = 4200 km/s. Most of these spectra were gathered at the Harvard-Smithsonian Center for Astrophysics (CfA) between 2004 and 2009. For 53 SNe, these are the first published spectra. The data coverage range from mere identification (1-3 spectra) for a few SNe to extensive series of observations (10-30 spectra) that trace the spectral evolution for others, with an average of 9 spectra per SN. For 44 SNe of the 73 SNe presented here, we have well-determined dates of maximum light to determine the phase of each spectrum. Our sample constitutes the most extensive spectral library of stripped-envelope SNe to date. We provide very early coverage (as early as 30 days before V-band max) for photospheric spectra, as well as late-time nebular coverage when the innermost regions of the SNe are visible (as late as 2 years after explosion, while for SN1993J, we have data as late as 11.6 years). This data set has homogeneous observations and reductions that allow us to study the spectroscopic diversity of these classes of stripped SNe and to compare these to SNe associated with gamma-ray bursts. We undertake these matters in follow-up papers.

Light echoes, light from a variable source scattered off dust, have been observed for over a century. The recent discovery of light echoes around centuries-old supernovae in the Milky Way and the Large Magellanic Cloud have allowed the spectroscopic characterization of these events, even without contemporaneous photometry and spectroscopy using modern instrumentation. Here we review the recent scientific advances using light echoes of ancient and historic transients, and focus on our latest work on SN 1987A's and Eta Carinae's light echoes.

We present follow-up optical imaging and spectroscopy of one of the light echoes of $\eta$ Carinae's 19th-century Great Eruption discovered by Rest et al. (2012). By obtaining images and spectra at the same light echo position between 2011 and 2014, we follow the evolution of the Great Eruption on a three-year timescale. We find remarkable changes in the photometric and spectroscopic evolution of the echo light. The $i$-band light curve shows a decline of $\sim 0.9$ mag in $\sim 1$ year after the peak observed in early 2011 and a flattening at later times. The spectra show a pure-absorption early G-type stellar spectrum at peak, but a few months after peak the lines of the [Ca II] triplet develop strong P-Cygni profiles and we see the appearance of [Ca II] 7291,7324 doublet in emission. These emission features and their evolution in time resemble those observed in the spectra of some Type IIn supernovae and supernova impostors. Most surprisingly, starting $\sim 300$ days after peak brightness, the spectra show strong molecular transitions of CN at $\gtrsim 6800$ \AA. The appearance of these CN features can be explained if the ejecta are strongly Nitrogen enhanced, as it is observed in modern spectroscopic studies of the bipolar Homunculus nebula. Given the spectroscopic evolution of the light echo, velocities of the main features, and detection of strong CN, we are likely seeing ejecta that contributes directly to the Homunculus nebula.

We present time series photometric and spectroscopic data for the transient SN 2009ip from the start of its outburst in September 2012 until November 2013. This data was collected primarily with the new robotic capabilities of the Las Cumbres Observatory Global Telescope Network, a specialized facility for time domain astrophysics, and includes supporting high-resolution spectroscopy from the Southern Astrophysical Research Telescope, Kitt Peak National Observatory, and Gemini Observatory. Based on our nightly photometric monitoring, we interpret the strength and timing of fluctuations in the light curve as interactions between fast-moving ejecta and an inhomogeneous CSM produced by past eruptions of this massive luminous blue variable (LBV) star. Our time series of spectroscopy in 2012 reveals that, as the continuum and narrow H-alpha flux from CSM interactions declines, the broad component of H-alpha persists with SN-like velocities that are not typically seen in LBVs or SN Impostor events. At late times we find that SN 2009ip continues to decline slowly, at <0.01 magnitudes per day, with small fluctuations in slope similar to Type IIn supernovae (SNe IIn) or SN impostors, but no further LBV-like activity. The late time spectrum features broad calcium lines similar to both late time SNe and SN Impostors. In general, we find that the photometric and spectroscopic evolution of SN 2009ip is more similar to SNe IIn than either continued eruptions of an LBV star or SN Impostors. In this context, we discuss the implications for episodic mass loss during the late stages of massive star evolution.

Combining stellar occultation observations probing Pluto's atmosphere from 1988 to 2013 and models of energy balance between Pluto's surface and atmosphere, we conclude that Pluto's atmosphere does not collapse at any point in its 248-year orbit. The occultation results show an increasing atmospheric pressure with time in the current epoch, a trend present only in models with a high thermal inertia and a permanent N2 ice cap at Pluto's north rotational pole.

Las Cumbres Observatory Global Telescope (LCOGT) is a young organization dedicated to time-domain observations at optical and (potentially) near-IR wavelengths. To this end, LCOGT is constructing a world-wide network of telescopes, including the two 2m Faulkes telescopes, as many as 17 x 1m telescopes, and as many as 23 x 40cm telescopes. These telescopes initially will be outfitted for imaging and (excepting the 40cm telescopes) spectroscopy at wavelengths between the atmospheric UV cutoff and the roughly 1-micron limit of silicon detectors. Since the first of LCOGT's 1m telescopes are now being deployed, we lay out here LCOGT's scientific goals and the requirements that these goals place on network architecture and performance, we summarize the network's present and projected level of development, and we describe our expected schedule for completing it. In the bulk of the paper, we describe in detail the technical approaches that we have adopted to attain the desired performance. In particular, we discuss our choices for the number and location of network sites, for the number and sizes of telescopes, for the specifications of the first generation of instruments, for the software that will schedule and control the network's telescopes and reduce and archive its data, and for the structure of the scientific and educational programs for which the network will provide observations.

Here we present observations of 7 large Kuiper Belt Objects. From these observations, we extract a point source catalog with $\sim0.01"$ precision, and astrometry of our target Kuiper Belt Objects with $0.04-0.08"$ precision within that catalog. We have developed a new technique to predict the future occurrence of stellar occultations by Kuiper Belt Objects. The technique makes use of a maximum likelihood approach which determines the best-fit adjustment to cataloged orbital elements of an object. Using simulations of a theoretical object, we discuss the merits and weaknesses of this technique compared to the commonly adopted ephemeris offset approach. We demonstrate that both methods suffer from separate weaknesses, and thus, together provide a fair assessment of the true uncertainty in a particular prediction. We present occultation predictions made by both methods for the 7 tracked objects, with dates as late as 2015. Finally, we discuss observations of three separate close passages of Quaoar to field stars, which reveal the accuracy of the element adjustment approach, and which also demonstrate the necessity of considering the uncertainty in stellar position when assessing potential occultations.

The 2012 explosion of SN2009ip raises questions about our understanding of the late stages of massive star evolution. Here we present a comprehensive study of SN2009ip during its remarkable re-brightening(s). High-cadence photometric and spectroscopic observations from the GeV to the radio band obtained from a variety of ground-based and space facilities (including the VLA, Swift, Fermi, HST and XMM) constrain SN2009ip to be a low energy (E~ 10^50 erg for an ejecta mass ~ 0.5 Msun) and likely asymmetric explosion in a complex medium shaped by multiple eruptions of the restless progenitor star. Most of the energy is radiated as a result of the shock breaking out through a dense shell of material located at 5x10^14 cm with M~0.1 Msun, ejected by the precursor outburst ~40 days before the major explosion. We interpret the NIR excess of emission as signature of dust vaporization of material located further out (R>4x 10^15 cm), the origin of which has to be connected with documented mass loss episodes in the previous years. Our modeling predicts bright neutrino emission associated with the shock break-out if the cosmic ray energy is comparable to the radiated energy. We connect this phenomenology with the explosive ejection of the outer layers of the massive progenitor star, that later interacted with material deposited in the surroundings by previous eruptions. Future observations will reveal if the luminous blue variable (LBV) progenitor star survived. Irrespective of whether the explosion was terminal, SN2009ip brought to light the existence of new channels for sustained episodic mass-loss, the physical origin of which has yet to be identified.

The Taiwanese-American Occultation Survey (TAOS) aims to detect serendipitous occultations of stars by small (about 1 km diameter) objects in the Kuiper Belt and beyond. Such events are very rare (<0.001 events per star per year) and short in duration (about 200 ms), so many stars must be monitored at a high readout cadence. TAOS monitors typically around 500 stars simultaneously at a 5 Hz readout cadence with four telescopes located at Lulin Observatory in central Taiwan. In this paper, we report the results of the search for small Kuiper Belt Objects (KBOs) in seven years of data. No occultation events were found, resulting in a 95% c.l. upper limit on the slope of the faint end of the KBO size distribution of q = 3.34 to 3.82, depending on the surface density at the break in the size distribution at a diameter of about 90 km.

We present an analysis of the maximum light, near ultraviolet (NUV; 2900-5500 A) spectra of 32 low redshift (0.001<z<0.08) Type Ia supernovae (SNe Ia), obtained with the Hubble Space Telescope (HST). We combine this spectroscopic sample with high-quality gri light curves obtained with robotic telescopes to measure photometric parameters, such as stretch, optical colour, and brightness. By comparing our data to a comparable sample of SNe Ia at intermediate-z (0.4<z<0.9), we detect modest spectral evolution (3-sigma), in the sense that our mean low-z NUV spectrum has a depressed flux compared to its intermediate-z counterpart. We also see a strongly increased dispersion about the mean with decreasing wavelength, confirming the results of earlier surveys. These trends are consistent with changes in metallicity as predicted by contemporary SN Ia spectral models. We also examine the properties of various NUV spectral diagnostics in the individual spectra. We find a general correlation between stretch and the velocity (or position) of many NUV spectral features. In particular, we observe that higher stretch SNe have larger Ca II H&K velocities, that also correlate with host galaxy stellar mass. This latter trend is probably driven by the well-established correlation between stretch and stellar mass. We find no trends between UV spectral features and optical colour. Mean spectra constructed according to whether the SN has a positive or negative Hubble residual show very little difference at NUV wavelengths, indicating that the NUV evolution and variation we identify do not directly correlate with Hubble residuals. Our work confirms and strengthens earlier conclusions regarding the complex behaviour of SNe Ia in the NUV spectral region, but suggests the correlations we find are more useful in constraining progenitor models than improving the use of SNe Ia as cosmological probes.

The nearby Type Ia supernova SN 2011fe in M101 (cz=241 km s^-1) provides a unique opportunity to study the early evolution of a "normal" Type Ia supernova, its compositional structure, and its elusive progenitor system. We present 18 high signal-to-noise spectra of SN 2011fe during its first month beginning 1.2 days post-explosion and with an average cadence of 1.8 days. This gives a clear picture of how various line-forming species are distributed within the outer layers of the ejecta, including that of unburned material (C+O). We follow the evolution of C II absorption features until they diminish near maximum light, showing overlapping regions of burned and unburned material between ejection velocities of 10,000 and 16,000 km s^-1. This supports the notion that incomplete burning, in addition to progenitor scenarios, is a relevant source of spectroscopic diversity among SNe Ia. The observed evolution of the highly Doppler-shifted O I 7774 absorption features detected within five days post-explosion indicate the presence of O I with expansion velocities from 11,500 to 21,000 km s^-1. The fact that some O I is present above C II suggests that SN 2011fe may have had an appreciable amount of unburned oxygen within the outer layers of the ejecta.

We present the study on GRB 071112C X-ray and optical light curves. In these two wavelength ranges, we have found different temporal properties. The R-band light curve showed an initial rise followed by a single power-law decay, while the X-ray light curve was described by a single power-law decay plus a flare-like feature. Our analysis shows that the observed temporal evolution cannot be described by the external shock model in which the X-ray and optical emission are produced by the same emission mechanism. No significant color changes in multi-band light curves and a reasonable value of the initial Lorentz factor ({\Gamma}0 = 275 \pm 20) in a uniform ISM support the afterglow onset scenario as the correct interpretation for the early R-band rise. The result suggests the optical flux is dominated by afterglow. Our further investigations show that the X-ray flux could be created by an additional feature related to energy injection and X-ray afterglow. Different theoretical interpretations indicate the additional feature in X-ray can be explained by either late internal dissipation or local inverse-Compton scattering in the external shock.

Eta Carinae (Eta Car) is one of the most massive binary stars in the Milky Way. It became the second-brightest star in the sky during its mid-19th century "Great Eruption," but then faded from view (with only naked-eye estimates of brightness). Its eruption is unique among known astronomical transients in that it exceeded the Eddington luminosity limit for 10 years. Because it is only 2.3 kpc away, spatially resolved studies of the nebula have constrained the ejected mass and velocity, indicating that in its 19th century eruption, Eta Car ejected more than 10 M_solar in an event that had 10% of the energy of a typical core-collapse supernova without destroying the star. Here we report the discovery of light echoes of Eta Carinae which appear to be from the 1838-1858 Great Eruption. Spectra of these light echoes show only absorption lines, which are blueshifted by -210 km/s, in good agreement with predicted expansion speeds. The light-echo spectra correlate best with those of G2-G5 supergiant spectra, which have effective temperatures of ~5000 K. In contrast to the class of extragalactic outbursts assumed to be analogs of Eta Car's Great Eruption, the effective temperature of its outburst is significantly cooler than allowed by standard opaque wind models. This indicates that other physical mechanisms like an energetic blast wave may have triggered and influenced the eruption.

Type Ia supernovae (SNe Ia) have been used empirically as standardized candles to reveal the accelerating universe even though fundamental details, such as the nature of the progenitor system and how the star explodes, remained a mystery. There is consensus that a white dwarf star explodes after accreting matter in a binary system, but the secondary could be anything from a main sequence star to a red giant, or even another white dwarf. The uncertainty stems from the fact that no recent SN Ia has been discovered close enough to detect the stars before explosion. Here we report early observations of SN 2011fe (PTF11kly) in M101 at a distance of 6.4 Mpc, the closest SN Ia in the past 25 years. We find that the exploding star was likely a carbon-oxygen white dwarf, and from the lack of an early shock we conclude that the companion was most likely a main sequence star. Early spectroscopy shows high-velocity oxygen that varies on a time scale of hours and extensive mixing of newly synthesized intermediate mass elements in the outermost layers of the supernova. A companion paper uses pre-explosion images to rule out luminous red giants and most helium stars as companions.

We calculate the average stretch-corrected rise-time of type Ia supernovae (SNe Ia) in the Supernova Legacy Survey. We use the aggregate lightcurves of spectroscopic and photometrically identified SNe Ia to fit the rising part of the lightcurve with a simple quadratic model. We obtain a lightcurve shape corrected, i .e. stretch-corrected, fiducial rise-time of 17.02^{+0.18}_{-0.28} (stat) days. The measured rise-time differs from an earlier finding by the SNLS (Conley et al. 2006) due to the use of different SN Ia templates. We compare it to nearby samples using the same methods and find no evolution in the early part of the lightcurve of SNe Ia up to z=1. We search for variations among different populations, particularly subluminous objects, by dividing the sample in stretch. Bright and slow decliners (s>1.0) have consistent stretch-corrected rise-times compared to fainter and faster decliners (0.8<s<1.0); they are shorter by 0.57^{+0.47}_{-0.50} (stat) days. Subluminous SNe Ia (here defined as objects with s<0.8), although less constrained, are also consistent, with a rise-time of 18.03^{+0.81}_{-1.37} (stat) days. We study several systematic biases and find that the use of different fiducial templates may affect the average rise-time but not the intrinsic differences between populations. Based on our results, we estimate that subluminous SNe Ia are powered by 0.05-0.35 solar masses of radioactive nickel synthesized in the explosion. Our conclusions are the same for the single-stretch and two-stretch parameterizations of the lightcurve.

While it is generally accepted that Type Ia supernovae are the result of the explosion of a carbon-oxygen White Dwarf accreting mass in a binary system, the details of their genesis still elude us, and the nature of the binary companion is uncertain. Kasen (2010) points out that the presence of a non-degenerate companion in the progenitor system could leave an observable trace: a flux excess in the early rise portion of the lightcurve caused by the ejecta impact with the companion itself. This excess would be observable only under favorable viewing angles, and its intensity depends on the nature of the companion. We searched for the signature of a non-degenerate companion in three years of Supernova Legacy Survey data by generating synthetic lightcurves accounting for the effects of shocking and comparing true and synthetic time series with Kolmogorov-Smirnov tests. Our most constraining result comes from noting that the shocking effect is more prominent in rest-frame B than V band: we rule out a contribution from white dwarf-red giant binary systems to Type Ia supernova explosions greater than 10% at 2 sigma, and than 20% at 3 sigma level.

Motivated by the increasing need for observational resources for the study of time varying astronomy, the Las Cumbres Observatory Global Telescope (LCOGT) is a private foundation, whose goal is to build a global network of robotic telescopes for scientific research and education. Once completed, the network will become a unique tool, capable of continuous monitoring from both the Northern and Southern Hemispheres. The network currently includes 2 x 2.0 m telescopes, already making an impact in the field of exoplanet research. In the next few years they will be joined by at least 12 x 1.0 m and 20 x 0.4 m telescopes. The increasing amount of LCOGT observational resources in the coming years will be of great service to the astronomical community in general, and the exoplanet community in particular.

The Taiwanese-American Occultation Survey (TAOS) project has collected more than a billion photometric measurements since 2005 January. These sky survey data-covering timescales from a fraction of a second to a few hundred days-are a useful source to study stellar variability. A total of 167 star fields, mostly along the ecliptic plane, have been selected for photometric monitoring with the TAOS telescopes. This paper presents our initial analysis of a search for periodic variable stars from the time-series TAOS data on one particular TAOS field, No. 151 (RA = 17$^{\rm h}30^{\rm m}6\fs$67, Dec = 27\degr17\arcmin 30\arcsec, J2000), which had been observed over 47 epochs in 2005. A total of 81 candidate variables are identified in the 3 square degree field, with magnitudes in the range 8 < R < 16. On the basis of the periodicity and shape of the lightcurves, 29 variables, 15 of which were previously unknown, are classified as RR Lyrae, Cepheid, delta Scuti, SX Phonencis, semi-regular and eclipsing binaries.

The Taiwanese-American Occultation Survey (TAOS) monitors fields of up to ~1000 stars at 5 Hz simultaneously with four small telescopes to detect occultation events from small (~1 km) Kuiper Belt Objects (KBOs). The survey presents a number of challenges, in particular the fact that the occultation events we are searching for are extremely rare and are typically manifested as slight flux drops for only one or two consecutive time series measurements. We have developed a statistical analysis technique to search the multi-telescope data set for simultaneous flux drops which provides a robust false positive rejection and calculation of event significance. In this paper, we describe in detail this statistical technique and its application to the TAOS data set.

We have analyzed the first 3.75 years of data from TAOS, the Taiwanese American Occultation Survey. TAOS monitors bright stars to search for occultations by Kuiper Belt Objects (KBOs). This dataset comprises 5e5 star-hours of multi-telescope photometric data taken at 4 or 5 Hz. No events consistent with KBO occultations were found in this dataset. We compute the number of events expected for the Kuiper Belt formation and evolution models of Pan & Sari (2005), Kenyon & Bromley (2004), Benavidez & Campo Bagatin (2009), and Fraser (2009). A comparison with the upper limits we derive from our data constrains the parameter space of these models. This is the first detailed comparison of models of the KBO size distribution with data from an occultation survey. Our results suggest that the KBO population is comprised of objects with low internal strength and that planetary migration played a role in the shaping of the size distribution.

We analyzed data accumulated during 2005 and 2006 by the Taiwan-American Occultation Survey (TAOS) in order to detect short-period variable stars (periods of <~ 1 hour) such as delta Scuti. TAOS is designed for the detection of stellar occultation by small-size Kuiper Belt Objects (KBOs) and is operating four 50cm telescopes at an effective cadence of 5Hz. The four telescopes simultaneously monitor the same patch of the sky in order to reduce false positives. To detect short-period variables, we used the Fast Fourier Transform algorithm (FFT) inasmuch as the data points in TAOS light-curves are evenly spaced. Using FFT, we found 41 short-period variables with amplitudes smaller than a few hundredths of a magnitude and periods of about an hour, which suggest that they are low-amplitude delta Scuti stars (LADS). The light-curves of TAOS delta Scuti stars are accessible online at the Time Series Center website (http://timemachine.iic.harvard.edu)

We present the results of a search for occultation events by objects at distances between 100 and 1000 AU in lightcurves from the Taiwanese-American Occultation Survey (TAOS). We searched for consecutive, shallow flux reductions in the stellar lightcurves obtained by our survey between 7 February 2005 and 31 December 2006 with a total of $\sim4.5\times10^{9}$ three-telescope simultaneous photometric measurements. No events were detected, allowing us to set upper limits on the number density as a function of size and distance of objects in Sedna-like orbits, using simple models.

We present a de-trending algorithm for the removal of trends in time series. Trends in time series could be caused by various systematic and random noise sources such as cloud passages, changes of airmass, telescope vibration or CCD noise. Those trends undermine the intrinsic signals of stars and should be removed. We determine the trends from subsets of stars that are highly correlated among themselves. These subsets are selected based on a hierarchical tree clustering algorithm. A bottom-up merging algorithm based on the departure from normal distribution in the correlation is developed to identify subsets, which we call clusters. After identification of clusters, we determine a trend per cluster by weighted sum of normalized light-curves. We then use quadratic programming to de-trend all individual light-curves based on these determined trends. Experimental results with synthetic light-curves containing artificial trends and events are presented. Results from other de-trending methods are also compared. The developed algorithm can be applied to time series for trend removal in both narrow and wide field astronomy.

We conducted a search for occultations of bright stars by Kuiper Belt Objects (KBOs) to estimate the density of sub-km KBOs in the sky. We report here the first results of this occultation survey of the outer solar system conducted in June 2007 and June/July 2008 at the MMT Observatory using Megacam, the large MMT optical imager. We used Megacam in a novel shutterless continuous--readout mode to achieve high precision photometry at 200 Hz. We present an analysis of 220 star hours at signal-to-noise ratio of 25 or greater. The survey efficiency is greater than 10% for occultations by KBOs of diameter d>=0.7 km, and we report no detections in our dataset. We set a new 95% confidence level upper limit for the surface density \Sigma_N(d) of KBOs larger than 1 km: \Sigma_N(d>=1 km) <= 2.0e8 deg^-2, and for KBOs larger than 0.7 km \Sigma_N(d>= 0.7 km) <= 4.8e8 deg^-2.

The Taiwanese-American Occultation Survey (TAOS) operates four fully automatic telescopes to search for occultations of stars by Kuiper Belt Objects. It is a versatile facility that is also useful for the study of initial optical GRB afterglows. This paper provides a detailed description of the TAOS multi-telescope system, control software, and high-speed imaging.

The star BD+29 1748 was resolved to be a close binary from its occultation by the asteroid 87 Sylvia on 2006 December 18 UT. Four telescopes were used to observe this event at two sites separated by some 80 km apart. Two flux drops were observed at one site, whereas only one flux drop was detected at the other. From the long-term variation of Sylvia, we inferred the probable shape of the shadow during the occultation, and this in turn constrains the binary parameters: the two components of BD+29 1748 have a projected separation of 0.097" to 0.110" on the sky with a position angle 104 deg to 107 deg. The asteroid was clearly resolved with a size scale ranging from 130 to 290 km, as projected onto the occultation direction. No occultation was detected for either of the two known moonlets of 87 Sylvia.

Results from the first two years of data from the Taiwanese-American Occultation Survey (TAOS) are presented. Stars have been monitored photometrically at 4 Hz or 5 Hz to search for occultations by small (~3 km) Kuiper Belt Objects (KBOs). No statistically significant events were found, allowing us to present an upper bound to the size distribution of KBOs with diameters 0.5 km < D < 28 km.

The serendipitous detection of stellar occultations by Outer Solar System objects is a powerful method for ascertaining the small end ($r \lesssim 15$ km) of the size distribution of Kuiper Belt Objects and may potentially allow the exploration of objects as far out as the Oort Cloud. The design and implementation of an occultation survey is aided by a detailed understanding of how diffraction and observational parameters affect the detection of occultation events. In this study, stellar occultations are simulated, accounting for diffraction effects, finite source sizes, finite bandwidths, stellar spectra, sampling, and signal-to-noise. Finally, the possibility of detecting small Outer Solar System objects from the Kuiper Belt all the way out to the Oort Cloud is explored for three photometric systems: a proposed space telescope, Whipple (Kaplan et al. 2003), the Taiwanese-American Occultation Survey (Lehner et al. 2006), and the Multi Mirror Telescope (Bianco 2007).

The Taiwan-America Occultation Survey (TAOS) aims to determine the number of small icy bodies in the outer reach of the Solar System by means of stellar occultation. An array of 4 robotic small (D=0.5 m), wide-field (f/1.9) telescopes have been installed at Lulin Observatory in Taiwan to simultaneously monitor some thousand of stars for such rare occultation events. Because a typical occultation event by a TNO a few km across will last for only a fraction of a second, fast photometry is necessary. A special CCD readout scheme has been devised to allow for stellar photometry taken a few times per second. Effective analysis pipelines have been developed to process stellar light curves and to correlate any possible flux changes among all telescopes. A few billion photometric measurements have been collected since the routine survey began in early 2005. Our preliminary result of a very low detection rate suggests a deficit of small TNOs down to a few km size, consistent with the extrapolation of some recent studies of larger (30--100 km) TNOs.