Search

Compact binary mergers detectable in gravitational waves can be accompanied by a kilonova, an electromagnetic transient powered by radioactive decay of newly synthesised r-process elements. A few kilonova candidates have been observed during short gamma-ray burst follow-up, and one found associated with a gravitational wave detection, GW170817. However, robust kilonova candidates are yet to be found in un-triggered, wide-field optical surveys, that is, a search not requiring an initial gravitational wave or gamma-ray burst trigger. Here we present the first observing run for the Kilonova and Transients Program (KNTraP) using the Dark Energy Camera. The first KNTraP run ran for 11 nights, covering 31 fields at a nightly cadence in two filters. The program can detect transients beyond the LIGO/Virgo/KAGRA horizon, be agnostic to the merger orientation, avoid the Sun and/or Galactic plane, and produces high cadence multi-wavelength light curves. The data were processed nightly in real-time for rapid identification of transient candidates, allowing for follow-up of interesting candidates before they faded away. Three fast-rising candidates were identified in real-time, however none had the characteristics of the kilonova AT2017gfo associated with GW170817 or with the expected evolution for kilonovae from our fade-rate models. After the run, the data were reprocessed, then subjected to stringent filtering and model fitting to search for kilonovae offline. Multiple KNTraP runs (3+) are expected to detect kilonovae via this optical-only search method. No kilonovae were detected in this first KNTraP run using our selection criteria, constraining the KN rate to $R < 1.8\times10^{5}$ Gpc$^{-3}$ yr$^{-1}$., Comment: 17 pages, 6 figures, 3 tables, submitted to MNRAS

GRB$\,$220831A is a gamma-ray burst (GRB) with a duration and spectral peak energy that places it at the interface between the distribution of long-soft and short-hard GRBs. In this paper, we present the multi-wavelength follow-up campaign to GRB$\,$220831A and its optical, near-infrared, X-ray and radio counterparts. Our deep optical and near-infrared observations do not reveal an underlying host galaxy, and establish that GRB$\,$220831A is observationally hostless to depth, $m_i\gtrsim26.6$ AB mag. Based on the Amati relation and the non-detection of an accompanying supernova, we find that this GRB is most likely to have originated from a collapsar at $z>2$, but it could also possibly be a compact object merger at $z<0.4$ with a large separation distance from its host galaxy. Regardless of its origin, we show that its optical and near-infrared counterpart departs from the evolution expected from a forward shock dominated synchrotron afterglow, exhibiting a steep post-break temporal powerlaw index of $-3.83^{+0.62}_{-0.79}$, too steep to be the jet-break. By analysing a range of models, we find that the observed steep departure from forward shock closure relations is likely due to an internal process producing either a flare or a plateau., Comment: 17 pages, 14 figures, 7 tables, submitted to MNRAS

The Legacy Survey of Space and Time (LSST) at Vera C. Rubin Observatory is planned to begin in the Fall of 2025. The LSST survey cadence has been designed via a community-driven process regulated by the Survey Cadence Optimization Committee (SCOC), which recommended up to 3% of the observing time to carry out Target of Opportunity (ToO) observations. Experts from the scientific community, Rubin Observatory personnel, and members of the SCOC were brought together to deliver a recommendation for the implementation of the ToO program during a workshop held in March 2024. Four main science cases were identified: gravitational wave multi-messenger astronomy, high energy neutrinos, Galactic supernovae, and small potentially hazardous asteroids possible impactors. Additional science cases were identified and briefly addressed in the documents, including lensed or poorly localized gamma-ray bursts and twilight discoveries. Trigger prioritization, automated response, and detailed strategies were discussed for each science case. This document represents the outcome of the Rubin ToO 2024 workshop, with additional contributions from members of the Rubin Science Collaborations. The implementation of the selection criteria and strategies presented in this document has been endorsed in the SCOC Phase 3 Recommendations document (PSTN-056). Although the ToO program is still to be finalized, this document serves as a baseline plan for ToO observations with the Rubin Observatory.

This paper presents the first public data release (DR1) of the FRB Line-of-sight Ionization Measurement From Lightcone AAOmega Mapping (FLIMFLAM) Survey, a wide field spectroscopic survey targeted on the fields of 10 precisely localized Fast Radio Bursts (FRBs). DR1 encompasses spectroscopic data for 10,468 galaxy redshifts across 10 FRBs fields with z<0.4, covering approximately 26 deg^2 of the sky in total. FLIMFLAM is composed of several layers, encompassing the `Wide' (covering ~ degree or >10 Mpc scales), `Narrow', (several-arcminute or ~ Mpc) and integral field unit (`IFU'; ~ arcminute or ~ 100 kpc ) components. The bulk of the data comprise spectroscopy from the 2dF-AAOmega on the 3.9-meter Anglo-Australian Telescope, while most of the Narrow and IFU data was achieved using an ensemble of 8-10-meter class telescopes. We summarize the information on our selected FRB fields, the criteria for target selection, methodologies employed for data reduction, spectral analysis processes, and an overview of our data products. An evaluation of our data reveals an average spectroscopic completeness of 48.43%, with over 80% of the observed targets having secure redshifts. Additionally, we describe our approach on generating angular masks and calculating the target selection functions, setting the stage for the impending reconstruction of the matter density field., Comment: 23 pages, 8 figures, Submitted to ApJS

The relativistic outflows that produce Long GRBs (LGRBs) can be described by a structured jet model where prompt $\gamma$-ray emission is restricted to a narrow region in the jet's core. Viewing the jet off-axis from the core, a population of afterglows without an associated GRB detection can be predicted. In this work, we conduct an archival search for these `orphan' afterglows (OAs) with minute-cadence, deep ($g\sim23$) data from the Dark Energy Camera (DECam) taken as part of the Deeper, Wider, Faster programme (DWF). We introduce a method to select fast-evolving OA candidates within DWF data that comprises a machine learning model, based on a realistic synthetic population of OAs. Using this classifier, we recover 51 OA candidates. Of these candidates, 42 are likely flare events from M-class stars. The remaining nine possess quiescent, coincident sources in archival data with angular profiles consistent with a star and are inconsistent with the expected population of LGRB host galaxies. We therefore conclude that these are likely Galactic events. We calculate an upper limit on the rate of OAs down to $g<22$ AB mag of 7.46\,deg$^{-2}$yr$^{-1}$ using our criteria and constrain possible jet structures. We also place an upper limit of the characteristic angle between the $\gamma$-ray emitting region and the jet's half opening angle. For a smooth power-law and a power-law with core jet model respectively, these values are $58.3^{\circ}$ and $56.6^{\circ}$, for a power-law index of 0.8 and $75.3^{\circ}$ and $76.8^{\circ}$ for a power-law index of 1.2., Comment: 17 pages, 13 figures, accepted for publication in MNRAS

The dispersion measure of fast radio bursts (FRBs), arising from the interactions of the pulses with free electrons along the propagation path, constitutes a unique probe of the cosmic baryon distribution. Their constraining power is further enhanced in combination with observations of the foreground large-scale structure and intervening galaxies. In this work, we present the first constraints on the partition of the cosmic baryons between the intergalactic medium (IGM) and circumgalactic medium (CGM), inferred from the FLIMFLAM spectroscopic survey. In its first data release, the FLIMFLAM survey targeted galaxies in the foreground of 8 localized FRBs. Using Bayesian techniques, we reconstruct the underlying ~Mpc-scale matter density field that is traced by the IGM gas. Simultaneously, deeper spectroscopy of intervening foreground galaxies (at impact parameters $b_\perp \lesssim r_{200}$) and the FRB host galaxies constrains the contribution from the CGM. Applying Bayesian parameter inference to our data and assuming a fiducial set of priors, we infer the IGM cosmic baryon fraction to be $f_{\rm igm}=0.59^{+0.11}_{-0.10}$, and a CGM gas fraction of $f_{\rm gas} = 0.55^{+0.26}_{-0.29}$ for $10^{10}\,M_\odot \lesssim M_{\rm halo}\lesssim 10^{13}\,M_\odot$ halos. The mean FRB host dispersion measure (rest-frame) in our sample is $\langle \rm{DM_{host}}\rangle = 90^{+29}_{-19}\rm{pc~cm^{-3}}$, of which $\langle{\rm DM_{host}^{unk}}\rangle =69^{+28}_{-19}~\rm{pc~cm^{-3}}$ arises from the host galaxy ISM and/or the FRB progenitor environment. While our current $f_{\rm igm}$ and $f_{\rm gas}$ uncertainties are too broad to constrain most galactic feedback models, this result marks the first measurement of the IGM and CGM baryon fractions, as well as the first systematic separation of the FRB host dispersion measure into two components: arising from the halo and from the inner ISM/FRB engine., Comment: 21 pages, 9 figures, 5 tables; Submitted to ApJ

We report for the first time a relationship between galaxy kinematics and net Lyman-alpha equivalent width (net Lya EW) in star forming galaxies during the epoch of peak cosmic star formation. Building on the previously reported broadband imaging segregation of Lya-emitting and Lya-absorbing Lyman break galaxies (LBGs) at z~2 (Paper I in this series) and previously at z~3, we use the Lya spectral type classification method to study the relationship between net Lya EW and nebular emission-line kinematics determined from IFU spectroscopy. We show that z~2 and z~3 LBGs segregate in colour-magnitude space according to their kinematic properties and Lyman-alpha spectral type, and conclude that LBGs with Lya dominant in absorption are almost exclusively rotation-dominated (presumably disc-like) systems, and LBGs with Lya dominant in emission characteristically have dispersion-dominated kinematics. We quantify the relationship between the strength of rotational dynamic support and net Lya EW, and demonstrate the consistency of our result with other properties that scale with net Lya EW and kinematics. Based on these findings, we suggest a method by which large samples of rotation- and dispersion-dominated galaxies might be selected using broadband imaging in as few as three filters and/or net Lya EW alone. Application of this method will enable an understanding of galaxy kinematic behaviour over large scales in datasets from current and future large-area and all-sky photometric surveys that will select hundreds of millions of LBGs in redshift ranges from z~2-6 across many hundreds to thousands of Mpc. Finally, we speculate that the combination of our result linking net Lya EW and nebular emission-line kinematics with the known large-scale clustering behaviour of Lya-absorbing and Lya-emitting LBGs is evocative of a nascent morphology-density relation at z~2-3., Comment: Published in PASA (Version 2: Typos corrected and minor clarifying edits in S3.1 & S3.2.1)

High-redshift Lyman break galaxies (LBGs) are efficiently selected in deep images using as few as three broadband filters, and have been shown to have multiple intrinsic and small- to large-scale environmental properties related to Lyman-alpha. In this paper we demonstrate a statistical relationship between net Lyman-alpha equivalent width (net Lya EW) and the optical broadband photometric properties of LBGs at z~2. We show that LBGs with the strongest net Lya EW in absorption (aLBGs) and strongest net Lya EW in emission (eLBGs) separate into overlapping but discrete distributions in $(U_n-R)$ colour and $R$-band magnitude space, and use this segregation behaviour to determine photometric criteria by which sub-samples with a desired Lya spectral type can be selected using data from as few as three broadband optical filters. We propose application of our result to current and future large-area and all-sky photometric surveys that will select hundreds of millions of LBGs across many hundreds to thousands of Mpc, and for which spectroscopic follow-up to obtain Lya spectral information is prohibitive. To this end, we use spectrophotometry of composite spectra derived from a sample of 798 LBGs divided into quartiles on the basis of net Lya EW to calculate criteria for the selection of Lya absorbing and Lya emitting populations of z~3 LBGs using $ugri$ broadband photometric data from the Vera Rubin Observatory Legacy Survey of Space and Time (LSST)., Comment: Published in Publications of the Astronomical Society of Australia 16 pages, 6 figures, 8 tables

In recent years, certain luminous extragalactic optical transients have been observed to last only a few days. Their short observed duration implies a different powering mechanism from the most common luminous extragalactic transients (supernovae) whose timescale is weeks. Some short-duration transients, most notably AT2018cow, display blue optical colours and bright radio and X-ray emission. Several AT2018cow-like transients have shown hints of a long-lived embedded energy source, such as X-ray variability, prolonged ultraviolet emission, a tentative X-ray quasiperiodic oscillation, and large energies coupled to fast (but subrelativistic) radio-emitting ejecta. Here we report observations of minutes-duration optical flares in the aftermath of an AT2018cow-like transient, AT2022tsd (the "Tasmanian Devil"). The flares occur over a period of months, are highly energetic, and are likely nonthermal, implying that they arise from a near-relativistic outflow or jet. Our observations confirm that in some AT2018cow-like transients the embedded energy source is a compact object, either a magnetar or an accreting black hole., Comment: 79 pages, 3 figures (main text) + 7 figures (extended data) + 2 figures (supplementary information). Published online in Nature on 15 November 2023

The repeating fast radio burst FRB20190520B is an anomaly of the FRB population thanks to its high dispersion measure (DM$=1205\,$pc/cc) despite its low redshift of $z_\mathrm{frb}=0.241$. This excess has been attributed to a large host contribution of $DM_{host}\approx 900\,$pc/cc, far larger than any other known FRB. In this paper, we describe spectroscopic observations of the FRB20190520B field obtained as part of the FLIMFLAM survey, which yielded 701 galaxy redshifts in the field. We find multiple foreground galaxy groups and clusters, for which we then estimated halo masses by comparing their richness with numerical simulations. We discover two separate $M_{halo} >10^{14}\,M_\odot$ galaxy clusters, at $z=0.1867$ and $z=0.2170$, respectively, that are directly intersected by the FRB sightline within their characteristic halo radius \rvir{}. Subtracting off their estimated DM contributions as well that of the diffuse intergalactic medium, we estimate a host contribution of $DM_{host}=430^{+140}_{-220}\,$pc/cc or $DM_{host}=280^{+140}_{-170}\,$pc/cc (observed frame) depending on whether we assume the halo gas extends to $r_{200}$ or $2\times r_{200}$. This significantly smaller $DM_{host}$ -- no longer the largest known value -- is now consistent with H$\alpha$ emission measures of the host galaxy without invoking unusually high gas temperatures. Combined with the observed FRB scattering timescale, we estimate the turbulent fluctuation and geometric amplification factor of the scattering layer to be $\tilde{F} G\approx4.5 - 11\,(\mathrm{pc^2\;km})^{-1/3}$, suggesting most of the gas is close to the FRB host. This result illustrates the importance of incorporating foreground data for FRB analyses, both for understanding the nature of FRBs and to realize their potential as a cosmological probe., Comment: Accepted by ApJL. Interactive figure included (link in text). Note numerical values have changed from v1

The FLIMFLAM survey is collecting spectroscopic data of field galaxies near fast radio burst (FRB) sightlines to constrain key parameters describing the distribution of matter in the Universe. In this work, we leverage the survey data to determine the source of the excess extragalactic dispersion measure (DM), compared to the Macquart relation estimate of four FRBs: FRB20190714A, FRB20200430A, FRB20200906A, and FRB20210117A. By modeling the gas distribution around the foreground galaxy halos and galaxy groups of the sightlines, we estimate $\rm DM_{halos}$, their contribution to the FRB dispersion measures. The FRB20190714A sightline shows a clear excess of foreground halos which contribute roughly 2/3$^{rd}$ of the observed excess DM, thus implying a baryon-dense sightline. FRB20200906A shows a smaller but non-negligible foreground halo contribution, and further analysis of the IGM is necessary to ascertain the true cosmic contribution to its DM. RB20200430A and FRB20210117A show negligible foreground contributions, implying a large host galaxy excess and/or progenitor environment excess., Comment: 14 pages, 6 figures, 3 tables. Data access available here: https://drive.google.com/drive/folders/1gVVwTAmiHJtFB_GD56hgmDylXHYbOTFS

Superluminous supernovae (SLSNe) have been detected to $z\sim4$ and can be detected to $z\gtrsim15$ using current and upcoming facilities. SLSNe are extremely UV luminous, and hence objects at $z\gtrsim7$ are detected exclusively via their rest-frame UV using optical and infrared facilities. SLSNe have great utility in multiple areas of stellar and galactic evolution. Here, we explore the potential use of SLSNe type-I as high-redshift cosmological distance indicators in their rest-frame UV. Using a SLSNe-I sample in the redshift range $1\lesssim z\lesssim 3$, we investigate correlations between the peak absolute magnitude in a synthetic UV filter centered at 250 nm and rise time, colour and decline rate of SLSNe-I light curves. We observe a linear correlation between $M_0(250)$ and the rise time with an intrinsic scatter of 0.29. Interestingly, this correlation is further tightened ($\sigma_{int} \approx 0.2$) by eliminating those SLSNe which show a pre-peak bump in their light curve. This result hints at the possibility that the "bumpy" SLSNe could belong to a different population. Weak correlations are observed between the peak luminosity and colour indices. No relationship is found between UV peak magnitude and the decline rate in contrast to what is typically found in optical band. The correlations found here are promising, and give encouraging insights for the use of SLSNe as cosmological probes at high redshifts using standardising relations in the UV. We also highlight the importance of early, and consistent, photometric data for constraining the light curve properties., Comment: Accepted for publication in MNRAS

Tidal disruption events (TDEs) are bursts of electromagnetic energy released when supermassive black holes (SMBHs) at the centers of galaxies violently disrupt a star that passes too close. TDEs provide a new window to study accretion onto SMBHs; in some rare cases, this accretion leads to launching of a relativistic jet, but the necessary conditions are not fully understood. The best studied jetted TDE to date is Swift J1644+57, which was discovered in gamma-rays, but was too obscured by dust to be seen at optical wavelengths. Here we report the optical discovery of AT2022cmc, a rapidly fading source at cosmological distance (redshift z=1.19325) whose unique lightcurve transitioned into a luminous plateau within days. Observations of a bright counterpart at other wavelengths, including X-rays, sub-millimeter, and radio, supports the interpretation of AT2022cmc as a jetted TDE containing a synchrotron "afterglow", likely launched by a SMBH with spin $a \gtrsim 0.3$. Using 4 years of Zwicky Transient Facility (ZTF) survey data, we calculate a rate of $0.02 ^{+ 0.04 }_{- 0.01 }$ Gpc$^{-3}$ yr$^{-1}$ for on-axis jetted TDEs based on the luminous, fast-fading red component, thus providing a measurement complementary to the rates derived from X-ray and radio observations. Correcting for the beaming angle effects, this rate confirms that about 1% of TDEs have relativistic jets. Optical surveys can use AT2022cmc as a prototype to unveil a population of jetted TDEs., Comment: Accepted for publication in Nature

Reservoirs of dense atomic gas (primarily hydrogen), contain approximately 90 percent of the neutral gas at a redshift of 3, and contribute to 2-3 percent of the total baryons in the Universe. These damped Lyman-${\alpha}$ systems (so called because they absorb Lyman-${\alpha}$ photons from within and from background sources) have been studied for decades, but only through absorption lines present in the spectra of background quasars and gamma-ray bursts. Such pencil beams do not constrain the physical extent of the systems. Here, we report integral-field spectroscopy of a bright, gravitationally lensed galaxy at a redshift of 2.7 with two foreground damped Lyman-${\alpha}$ systems. These systems are $>$ 238 $kpc^2$ in extent, with column densities of neutral hydrogen varying by more than an order of magnitude on $<$ 3 kpc-scales. The mean column densities are $10^{20.46}$ - $10^{20.84} cm^{-2}$ and the total masses are $> 5.5 \times 10^{8}$ - $1.4 \times 10^{9} M_{\odot}$, showing that they contain the necessary fuel for the next generation of star formation, consistent with relatively massive, low-luminosity primeval galaxies at redshifts $>$ 2., Comment: 31 pages, 4 Figures, 4 Extended Data figures, 2 Extended Data tables, Author's version, Accepted: 4 March 2022, Published online by Nature on May 18, 2022

Star-forming galaxies are the sources likely to have reionized the universe. As we cannot observe them directly due to the opacity of the intergalactic medium at $z\gtrsim5$, we study $z\sim3\text{--}5$ galaxies as proxies to place observational constraints on cosmic reionization. Using new deep \textit{Hubble Space Telescope} rest-frame UV F336W and F435W imaging (30-orbit, $\sim40$~arcmin$^2$, $\sim29\text{--}30$~mag depth at 5$\sigma$), we attempt to identify a sample of Lyman continuum galaxies (LCGs). These are individual sources that emit ionizing flux below the Lyman break ($<912~\text{\AA}$). This population would allow us to constrain cosmic reionization parameters such as the number density and escape fraction ($f_{\rm esc}$) of ionizing sources. We compile a comprehensive parent sample that does not rely on the Lyman-break technique for redshifts. We present three new spectroscopic candidates at $z\sim3.7\text{--}4.4$, and 32 new photometric candidates. The high-resolution multi-band HST imaging and new Keck/Low Resolution Imaging Spectrometer (LRIS) redshifts make these promising spectroscopic LCG candidates. Using both a traditional and probabilistic approach, we find the most likely $f_{\rm esc}$ values for the three spectroscopic LCG candidates are $>100\%$, and therefore not physical. We are unable to confirm the true nature of these sources with the best available imaging and direct blue Keck/LRIS spectroscopy. More spectra, especially from the new class of 30 m telescopes, will be required to build a statistical sample of LCGs to place firm observational constraints on cosmic reionization., Comment: Accepted for publication in ApJ. 38 pages (28 in main body), 14 figures

On 2019 August 14 at 21:10:39 UTC, the LIGO/Virgo Collaboration (LVC) detected a possible neutron star-black hole merger (NSBH), the first ever identified. An extensive search for an optical counterpart of this event, designated GW190814, was undertaken using the Dark Energy Camera (DECam) on the 4m Victor M. Blanco Telescope at the Cerro Tololo Inter-American Observatory. Target of Opportunity interrupts were issued on 8 separate nights to observe 11 candidates using the 4.1m Southern Astrophysical Research (SOAR) telescope's Goodman High Throughput Spectrograph in order to assess whether any of these transients was likely to be an optical counterpart of the possible NSBH merger. Here, we describe the process of observing with SOAR, the analysis of our spectra, our spectroscopic typing methodology, and our resultant conclusion that none of the candidates corresponded to the gravitational wave merger event but were all instead other transients. Finally, we describe the lessons learned from this effort. Application of these lessons will be critical for a successful community spectroscopic follow-up program for LVC observing run 4 (O4) and beyond., Comment: 32 pages, 18 figures, accepted for publication by ApJ

The current data acquisition rate of astronomical transient surveys and the promise for significantly higher rates during in the next decade necessitate the development of novel approaches to analyze astronomical data sets and promptly detect objects of interest. The Deeper, Wider, Faster (DWF) program is a survey focused on the identification of fast evolving transients, such as fast radio bursts, gamma-ray bursts, and supernova shock breakouts. It employs a multi-frequency simultaneous coverage of the same part of the sky over several orders of magnitude. Using the Dark Energy Camera mounted on the 4-meter Blanco telescope, DWF captures a 20 second g-band exposure every minute, at a typical seeing of ~ 1" and an airmass of ~ 1.5. These optical data are collected simultaneously with observations conducted over the entire electromagnetic spectrum - from radio to gamma-rays - as well as cosmic ray observations. In this paper, we present a novel real-time light curve analysis algorithm, designed to detect transients in the DWF optical data; this algorithm functions independently from, or in conjunction with, image subtraction. We present a sample of fast transients detected by our algorithm, as well as a false-positive analysis. Our algorithm is customizable and can be tuned to be sensitive to transients evolving over different timescales and flux ranges., Comment: AJ accepted, 19 pages, 11 figures

The dispersion measures (DM) of fast radio bursts (FRBs) encode the integrated electron density along the line-of-sight, which is typically dominated by the intergalactic medium (IGM) contribution in the case of extragalactic FRBs. In this paper, we show that incorporating wide-field spectroscopic galaxy survey data in the foreground of localized FRBs can significantly improve constraints on the partition of diffuse cosmic baryons. Using mock DMs and realistic lightcone galaxy catalogs derived from the Millennium simulation, we define spectroscopic surveys that can be carried out with 4m and 8m-class wide field spectroscopic facilities. On these simulated surveys, we carry out Bayesian density reconstructions in order to estimate the foreground matter density field. In comparison with the `true' matter density field, we show that these can help reduce the uncertainties in the foreground structures by $\sim 2-3\times$ compared to cosmic variance. We calculate the Fisher matrix to forecast that $N=30\: (96)$ localized FRBs should be able to constrain the diffuse cosmic baryon fraction to $\sim 10\%\: (\sim 5\%) $, and parameters governing the size and baryon fraction of galaxy circumgalactic halos to within $\sim 20-25\%\: (\sim 8-12\%)$. From the Fisher analysis, we show that the foreground data increases the sensitivity of localized FRBs toward our parameters of interest by $\sim 25\times$. We briefly introduce FLIMFLAM, an ongoing galaxy redshift survey that aims to obtain foreground data on $\sim 30$ localized FRB fields., Comment: 27 pages, 12 Figures. Accepted by ApJ

We present our 500 pc distance-limited study of stellar fares using the Dark Energy Camera as part of the Deeper, Wider, Faster Program. The data was collected via continuous 20-second cadence g band imaging and we identify 19,914 sources with precise distances from Gaia DR2 within twelve, ~3 square-degree, fields over a range of Galactic latitudes. An average of ~74 minutes is spent on each field per visit. All light curves were accessed through a novel unsupervised machine learning technique designed for anomaly detection. We identify 96 flare events occurring across 80 stars, the majority of which are M dwarfs. Integrated are energies range from $\sim 10^{31}-10^{37}$ erg, with a proportional relationship existing between increased are energy with increased distance from the Galactic plane, representative of stellar age leading to declining yet more energetic are events. In agreement with previous studies we observe an increase in flaring fraction from M0 -> M6 spectral types. Furthermore, we find a decrease in the flaring fraction of stars as vertical distance from the galactic plane is increased, with a steep decline present around ~100 pc. We find that ~70% of identified flares occur on short timescales of ~8 minutes. Finally we present our associated are rates, finding a volumetric rate of $2.9 \pm 0.3 \times 10^{-6}$ flares pc$^{-3}$ hr$^{-1}$., Comment: 20 pages, 14 figures, 6 tables

We present optical follow-up imaging obtained with the Katzman Automatic Imaging Telescope, Las Cumbres Observatory Global Telescope Network, Nickel Telescope, Swope Telescope, and Thacher Telescope of the LIGO/Virgo gravitational wave (GW) signal from the neutron star-black hole (NSBH) merger GW190814. We searched the GW190814 localization region (19 deg$^{2}$ for the 90th percentile best localization), covering a total of 51 deg$^{2}$ and 94.6% of the two-dimensional localization region. Analyzing the properties of 189 transients that we consider as candidate counterparts to the NSBH merger, including their localizations, discovery times from merger, optical spectra, likely host-galaxy redshifts, and photometric evolution, we conclude that none of these objects are likely to be associated with GW190814. Based on this finding, we consider the likely optical properties of an electromagnetic counterpart to GW190814, including possible kilonovae and short gamma-ray burst afterglows. Using the joint limits from our follow-up imaging, we conclude that a counterpart with an $r$-band decline rate of 0.68 mag day$^{-1}$, similar to the kilonova AT 2017gfo, could peak at an absolute magnitude of at most $-17.8$ mag (50% confidence). Our data are not constraining for ''red'' kilonovae and rule out ''blue'' kilonovae with $M>0.5 M_{\odot}$ (30% confidence). We strongly rule out all known types of short gamma-ray burst afterglows with viewing angles $<$17$^{\circ}$ assuming an initial jet opening angle of $\sim$$5.2^{\circ}$ and explosion energies and circumburst densities similar to afterglows explored in the literature. Finally, we explore the possibility that GW190814 merged in the disk of an active galactic nucleus, of which we find four in the localization region, but we do not find any candidate counterparts among these sources., Comment: 86 pages, 9 figures

Some supernovae such as pair-instability supernovae are predicted to have the duration of more than a year in the observer frame. To constrain the rates of supernovae lasting for more than a year, we conducted a long-term deep transient survey using Hyper Suprime-Cam (HSC) on the 8.2m Subaru telescope. HSC is a wide-field (a 1.75 deg2 field-of-view) camera and it can efficiently conduct transient surveys. We observed the same 1.75 deg2 field repeatedly using the g, r, i, and z band filters with the typical depth of 26 mag for 4 seasons (from late 2016 to early 2020). Using these data, we searched for transients lasting for more than a year. Two supernovae were detected in 2 continuous seasons, one supernova was detected in 3 continuous seasons, but no transients lasted for all 4 seasons searched. The discovery rate of supernovae lasting for more than a year with the typical limiting magnitudes of 26 mag is constrained to be 1.4^{+1.3}_{-0.7}(stat.)^{+0.2}_{-0.3}(sys.) events deg-2 yr-1. All the long-lasting supernovae we found are likely Type IIn supernovae and our results indicate that about 40% of Type IIn supernovae have long-lasting light curves. No plausible pair-instability supernova candidates lasting for more than a year are discovered. By comparing the survey results and survey simulations, we constrain the luminous pair-instability supernova rate up to z ~ 3 should be of the order of 100 Gpc-3 yr-1 at most, which is 0.01 - 0.1 per cent of the core-collapse supernova rate., Comment: 18 pages, 9 figures, 6 tables, accepted by The Astrophysical Journal

Identification of anomalous light curves within time-domain surveys is often challenging. In addition, with the growing number of wide-field surveys and the volume of data produced exceeding astronomers ability for manual evaluation, outlier and anomaly detection is becoming vital for transient science. We present an unsupervised method for transient discovery using a clustering technique and the Astronomaly package. As proof of concept, we evaluate 85553 minute-cadenced light curves collected over two 1.5 hour periods as part of the Deeper, Wider, Faster program, using two different telescope dithering strategies. By combining the clustering technique HDBSCAN with the isolation forest anomaly detection algorithm via the visual interface of Astronomaly, we are able to rapidly isolate anomalous sources for further analysis. We successfully recover the known variable sources, across a range of catalogues from within the fields, and find a further 7 uncatalogued variables and two stellar flare events, including a rarely observed ultra fast flare (5 minute) from a likely M-dwarf., Comment: Accepted 7 Aug 2020, 19 pages, 8 figures

The first binary neutron star merger, GW170817, was accompanied by a radioactivity-powered optical/infrared transient called a kilonova. To date, no compelling kilonova has been found during optical surveys of the sky, independent of gravitational-wave triggers. In this work, we searched the first 23 months of the Zwicky Transient Facility (ZTF) data stream for candidate kilonovae in the form of rapidly evolving transients. We combined ZTF alert queries with forced point-spread-function photometry and nightly flux stacking to increase our sensitivity to faint and fast transients. Automatic queries yielded $>11,200$ candidates, 24 of which passed quality checks and strict selection criteria based on a grid of kilonova models tailored for both binary neutron star and neutron star-black hole mergers. None of the candidates in our sample was deemed a possible kilonova after thorough vetting, catalog cross-matching, and study of their color evolution. The sources that passed our selection criteria are dominated by Galactic cataclysmic variables. In addition, we identified two fast transients at high Galactic latitude, one of which is the confirmed afterglow of long-duration GRB190106A, and the other is a possible cosmological afterglow. Using a survey simulation code, we constrained the kilonova rate for a range of models including top-hat and linearly decaying light curves and synthetic light curves obtained with radiative transfer simulations. For prototypical GW170817-like kilonovae, we constrain the rate to be $R < 1775$ Gpc$^{-3}$ yr$^{-1}$ at 95% confidence level by requiring at least 2 high-significance detections. By assuming a population of kilonovae with the same geometry and composition of GW170817 observed under a uniform viewing angle distribution, we obtained a constraint on the rate of $R < 4029$ Gpc$^{-3}$ yr$^{-1}$., Comment: Submitted for publication in ApJ

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged), Comment: Submitted to ApJ

Synoptic searches for the optical counterpart to a binary neutron star (BNS) or neutron star-black hole (NSBH) merger can pose significant challenges towards the discovery of kilonovae and performing multi-messenger science. In this work, we describe the advantage of a global multi-telescope network towards this end, with a particular focus on the key and complementary role the Dark Energy Camera (DECam) plays in multi-facility follow-up. We describe the Global Relay of Observatories Watching Transients Happen (GROWTH) Target-of-Opportunity (ToO) Marshal, a common web application we built to ingest events, plan observations, search for transient candidates, and retrieve performance summary statistics for all of the telescopes in our network. Our infrastructure enabled us to conduct observations of two events during O3a, S190426c and S190510g. Furthermore, our analysis of deep DECam observations of S190814bv conducted by the DESGW team, and access to a variety of global follow-up facilities allowed us to place meaningful constraints on the parameters of the kilonova and the merging binary. We emphasize the importance of a global telescope network in conjunction with a power telescope like DECam in performing searches for the counterparts to gravitational-wave sources., Comment: 9 pages, 4 figures, submitted to Rev. Mex. A. A. (AstroRob 2019 conference proceedings)

On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to a NSBH merger. To search for potential counterparts, the GROWTH collaboration performed real-time image subtraction on 6 nights of public Dark Energy Camera (DECam) images acquired in the three weeks following the merger, covering $>$98\% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and we placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be $M_\mathrm{ej} < 0.04$~$M_{\odot}$ at polar viewing angles, or $M_\mathrm{ej} < 0.03$~$M_{\odot}$ if the opacity is $\kappa < 2$~cm$^2$g$^{-1}$. Assuming a tidal deformability for the neutron star at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be $ \chi < 0.7$ for mass ratios $Q < 6$, with weaker constraints for more compact neutron stars. We publicly release the photometry from this campaign at http://www.astro.caltech.edu/~danny/static/s190814bv., Comment: Accepted to ApJ with no major changes

The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg$^2$ of the initial skymap over the next two nights, corresponding to 46\% integrated probability, ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% integrated probability based on the initial sky map. The revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient "alerts" over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical lightcurves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled-out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible.

We report our observations of HSC16aayt (SN 2016jiu), which was discovered by the Subaru/Hyper Suprime-Cam (HSC) transient survey conducted as part of Subaru Strategic Program (SSP). It shows very slow photometric evolution and its rise time is more than 100 days. The optical magnitude change in 400 days remains within 0.6 mag. Spectra of HSC16aayt show a strong narrow emission line and we classify it as a Type IIn supernova. The redshift of HSC16aayt is 0.6814 +/- 0.0002 from the spectra. Its host galaxy center is at 5 kpc from the supernova location and HSC16aayt might be another example of isolated Type IIn supernovae, although the possible existence of underlying star forming activity of the host galaxy at the supernova location is not excluded., Comment: 14 pages, 12 figures, 1 table, accepted by The Astrophysical Journal

The first two months of the third Advanced LIGO and Virgo observing run (2019 April-May) showed that distant gravitational wave (GW) events can now be readily detected. Three candidate mergers containing neutron stars (NS) were reported in a span of 15 days, all likely located more than 100 Mpc away. However, distant events such as the three new NS mergers are likely to be coarsely localized, which highlights the importance of facilities and scheduling systems that enable deep observations over hundreds to thousands of square degrees to detect the electromagnetic counterparts. On 2019-05-10 02:59:39.292 UT the GW candidate S190510g was discovered and initially classified as a BNS merger with 98% probability. The GW event was localized within an area of 3462 deg2, later refined to 1166 deg2 (90%) at a distance of 227 +- 92 Mpc. We triggered Target of Opportunity observations with the Dark Energy Camera (DECam), a wide-field optical imager mounted at the prime focus of the 4m Blanco Telescope at CTIO in Chile. This Letter describes our DECam observations and our real-time analysis results, focusing in particular on the design and implementation of the observing strategy. Within 24 hours of the merger time, we observed 65% of the total enclosed probability of the final skymap with an observing efficiency of 94%. We identified and publicly announced 13 candidate counterparts. S190510g was re-classified 1.7 days after the merger, after our observations were completed, with a "binary neutron star merger" probability reduced from 98% to 42% in favor of a "terrestrial" classification., Comment: 4 figures, 3 tables, accepted for publication in ApJL

The discovery of a transient kilonova following the gravitational-wave event GW170817 highlighted the critical need for coordinated rapid and wide-field observations, inference, and follow-up across the electromagnetic spectrum. In the Southern hemisphere, the Dark Energy Camera (DECam) on the Blanco 4-m telescope is well-suited to this task, as it is able to cover wide-fields quickly while still achieving the depths required to find kilonovae like the one accompanying GW170817 to $\sim$500 Mpc, the binary neutron star horizon distance for current generation of LIGO/Virgo collaboration (LVC) interferometers. Here, as part of the multi-facility followup by the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration, we describe the observations and automated data movement, data reduction, candidate discovery, and vetting pipeline of our target-of-opportunity DECam observations of S190426c, the first possible neutron star--black hole merger detected via gravitational waves. Starting 7.5hr after S190426c, over 11.28\,hr of observations, we imaged an area of 525\,deg$^2$ ($r$-band) and 437\,deg$^2$ ($z$-band); this was 16.3\% of the total original localization probability and nearly all of the probability density visible from the Southern hemisphere. The machine-learning based pipeline was optimized for fast turnaround, delivering transient candidates for human vetting within 17 minutes, on average, of shutter closure. We reported nine promising counterpart candidates 2.5 hours before the end of our observations. Our observations yielded no detection of a bona fide counterpart to $m_z = 22.5$ and $m_r = 22.9$ at the 5$\sigma$ level of significance, consistent with the refined LVC positioning. We view these observations and rapid inferencing as an important real-world test for this novel end-to-end wide-field pipeline., Comment: submitted to AAS journals

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

We present the Deeper Wider Faster (DWF) program that coordinates more than 30 multi-wavelength and multi-messenger facilities worldwide and in space to detect and study fast transients (millisecond-to-hours duration). DWF has four main components, (1) simultaneous observations, where about 10 major facilities, from radio to gamma-ray, are coordinated to perform deep, wide-field, fast-cadenced observations of the same field at the same time. Radio telescopes search for fast radio bursts while optical imagers and high-energy instruments search for seconds-to-hours timescale transient events, (2) real-time (seconds to minutes) supercomputer data processing and candidate identification, along with real-time (minutes) human inspection of candidates using sophisticated visualisation technology, (3) rapid-response (minutes) follow-up spectroscopy and imaging and conventional ToO observations, and (4) long-term follow up with a global network of 1-4m-class telescopes. The principal goals of DWF are to discover and study counterparts to fast radio bursts and gravitational wave events, along with millisecond-to-hour duration transients at all wavelengths., Comment: 4 pages, 1 figure, proceedings of the IAU Symposium 339

We report our first discoveries of high-redshift supernovae from the Subaru HIgh-Z sUpernova CAmpaign (SHIZUCA), the transient survey using Subaru/Hyper Suprime-Cam. We report the discovery of three supernovae at the spectroscopically-confirmed redshifts of 2.399 (HSC16adga), 1.965 (HSC17auzg), and 1.851 (HSC17dbpf), and two supernova candidates with the host-galaxy photometric redshifts of 3.2 (HSC16apuo) and 4.2 (HSC17dsid), respectively. In this paper, we present their photometric properties and the spectroscopic properties of the confirmed high-redshift supernovae are presented in the accompanying paper (Curtin et al. 2019). The supernovae with the confirmed redshifts of z ~ 2 have the rest ultraviolet peak magnitudes close to -21 mag and they are likely superluminous supernovae. The discovery of three supernovae at z ~ 2 roughly corresponds to the approximate event rate of ~ 900 +/- 520 Gpc-3 yr-1 with Poisson error, which is consistent with the total superluminous supernova rate estimated by extrapolating the local rate based on the cosmic star-formation history. Adding unconfirmed superluminous supernova candidates would increase the event rate. Our superluminous supernova candidates at the redshifts of around 3 and 4 indicate the approximate superluminous supernova rates of ~ 400 +/- 400 Gpc-3 yr-1 (z ~ 3) and ~ 500 +/- 500 Gpc-3 yr-1 (z ~ 4) with Poisson errors. Our initial results demonstrate the outstanding capability of Hyper Suprime-Cam to discover high-redshift supernovae., Comment: 27 pages, 15 figures, 6 tables, accepted by The Astrophysical Journal Supplements

We present Keck spectroscopic observations of three probable high redshift superluminous supernovae (SLSNe) from the Subaru HIgh-Z sUpernova CAmpaign (SHIZUCA), confirming redshifts of 1.851, 1.965 and 2.399. The host galaxies were selected for transient monitoring from multi-band photometric redshifts. The supernovae are detected during their rise, and the classically scheduled spectra are collected near maximum light. The restframe far-ultraviolet (FUV;$\sim$1000A-2500A) spectra include a significant host galaxy flux contribution and we compare our host galaxy subtracted spectra to UV-luminous SNe from the literature. While the signal-to-noise ratios of the spectra presented here are sufficient for redshift confirmation, supernova spectroscopic type confirmation remains inconclusive. The success of the first SHIZUCA Keck spectroscopic follow-up program demonstrates that campaigns such as SHIZUCA are capable of identifying high redshift SLSNe with sufficient accuracy, speed and depth for rapid, well-cadenced and informative follow-up., Comment: 15 pages, 6 figures, submitted to the Astrophysical Journal

We report on SALT low resolution optical spectroscopy and optical/IR photometry undertaken with other SAAO telescopes (MASTER-SAAO and IRSF) of the kilonova AT 2017gfo (aka SSS17a) in the galaxy NGC4993 during the first 10 days of discovery. This event has been identified as the first ever electromagnetic counterpart of a gravitational wave event, namely GW170817, which was detected by the LIGO and Virgo gravitational wave observatories. The event is likely due to a merger of two neutron stars, resulting in a kilonova explosion. SALT was the third telescope to obtain spectroscopy of AT 2017gfo and the first spectrum, 1.2 d after the merger, is quite blue and shows some broad features, but no identifiable spectral lines and becomes redder over time. We compare the spectral and photometric evolution with recent kilonova simulations and conclude that they are in qualitative agreement for post-merger wind models with proton: nucleon ratios of $Y_e$ = 0.25$-$0.30. The blue colour of the first spectrum is consistent with the lower opacity of the Lathanide-free r-process elements in the ejecta. Differences between the models and observations are likely due to the choice of system parameters combined with the absence of atomic data for more elements in the ejecta models., Comment: 5 pages, 4 figures, 1 table

The ability to quickly detect transient sources in optical images and trigger multi-wavelength follow up is key for the discovery of fast transients. These include events rare and difficult to detect such as kilonovae, supernova shock breakout, and "orphan" Gamma-ray Burst afterglows. We present the Mary pipeline, a (mostly) automated tool to discover transients during high-cadenced observations with the Dark Energy Camera (DECam) at CTIO. The observations are part of the "Deeper Wider Faster" program, a multi-facility, multi-wavelength program designed to discover fast transients, including counterparts to Fast Radio Bursts and gravitational waves. Our tests of the Mary pipeline on DECam images return a false positive rate of ~2.2% and a missed fraction of ~3.4% obtained in less than 2 minutes, which proves the pipeline to be suitable for rapid and high-quality transient searches. The pipeline can be adapted to search for transients in data obtained with imagers other than DECam., Comment: Accepted for publication in PASA; 15 pages, 7 figures, 3 tables

An interesting transient has been detected in one of our three Dark Energy Camera deep fields. Observations of these deep fields take advantage of the high red sensitivity of DECam on the Cerro Tololo Interamerican Observatory Blanco telescope. The survey includes the Y band with rest wavelength 1430{\AA} at z = 6. Survey fields (the Prime field 0555-6130, the 16hr field 1600-75 and the SUDSS New Southern Field) are deeper in Y than other infrared surveys. They are circumpolar, allowing all night to be used efficiently, exploiting the moon tolerance of 1 micron observations to minimize conflict with the Dark Energy Survey. As an i-band dropout (meaning that the flux decrement shortward of Lyman alpha is in the i bandpass), the transient we report here is a supernova candidate with z ~ 6, with a luminosity comparable to the brightest known current epoch superluminous supernova (i.e., ~ 2 x 10^11 solar luminosities)., Comment: Reference added

By applying a display ecology to the {\em Deeper, Wider, Faster} proactive, simultaneous telescope observing campaign, we have shown a dramatic reduction in the time taken to inspect DECam CCD images for potential transient candidates and to produce time-critical triggers to standby telescopes. We also show how facilitating rapid corroboration of potential candidates and the exclusion of non-candidates improves the accuracy of detection; and establish that a practical and enjoyable workspace can improve the experience of an otherwise taxing task for astronomers. We provide a critical road-test of two advanced displays in a research context -- a rare opportunity to demonstrate how they can be used rather than simply discuss how they might be used to accelerate discovery., Comment: 19 pages, 11 figues, 2 tables, accepted to Publications of the Astronomical Society of Australia, March 2017

Superluminous supernovae are beginning to be discovered at redshifts as early as the epoch of reionization. A number of candidate mechanisms is reviewed, together with the discovery programs., Comment: for publication in PASA

The high Antarctic plateau provides exceptional conditions for conducting infrared observations of the cosmos on account of the cold, dry and stable atmosphere above the ice surface. This paper describes the scientific goals behind the first program to examine the time-varying universe in the infrared from Antarctica - the Kunlun Infrared Sky Survey (KISS). This will employ a small (50 cm aperture) telescope to monitor the southern skies in the 2.4um Kdark window from China's Kunlun station at Dome A, on the summit of the Antarctic plateau, through the uninterrupted 4-month period of winter darkness. An earlier paper discussed optimisation of the Kdark filter for the best sensitivity (Li et al 2016). This paper examines the scientific program for KISS. We calculate the sensitivity of the camera for the extrema of observing conditions that will be encountered. We present the parameters for sample surveys that could then be carried out for a range of cadences and sensitivities. We then discuss several science programs that could be conducted with these capabilities, involving star formation, brown dwarfs and hot Jupiters, exoplanets around M dwarfs, the terminal phases of stellar evolution, discovering fast transients as part of multi-wavelength campaigns, embedded supernova searches, reverberation mapping of active galactic nuclei, gamma ray bursts and the detection of the cosmic infrared background. Accepted for publication in PASA, 04/08/16., Comment: 21 pages, 6 Figures

The Thirty Meter Telescope (TMT) first light instrument IRIS (Infrared Imaging Spectrograph) will complete its preliminary design phase in 2016. The IRIS instrument design includes a near-infrared (0.85 - 2.4 micron) integral field spectrograph (IFS) and imager that are able to conduct simultaneous diffraction-limited observations behind the advanced adaptive optics system NFIRAOS. The IRIS science cases have continued to be developed and new science studies have been investigated to aid in technical performance and design requirements. In this development phase, the IRIS science team has paid particular attention to the selection of filters, gratings, sensitivities of the entire system, and science cases that will benefit from the parallel mode of the IFS and imaging camera. We present new science cases for IRIS using the latest end-to-end data simulator on the following topics: Solar System bodies, the Galactic center, active galactic nuclei (AGN), and distant gravitationally-lensed galaxies. We then briefly discuss the necessity of an advanced data management system and data reduction pipeline., Comment: 15 pages, 7 figures, SPIE (2016) 9909-05

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged), Comment: 210 pages, 91 figures. Exposure draft. Appendices to the Detailed Science Case can be found at http://mse.cfht.hawaii.edu/docs/

We present and discuss ultraviolet and optical photometry from the Ultraviolet/Optical Telescope and X-ray limits from the X-Ray Telescope on Swift and imaging polarimetry and ultraviolet/optical spectroscopy with the Hubble Space Telescope of ASASSN-15lh. It has been classified as a hydrogen-poor superluminous supernova (SLSN I) more luminous than any other supernova observed. ASASSN-15lh is not detected in the X-rays in individual or coadded observations. From the polarimetry we determine that the explosion was only mildly asymmetric. We find the flux of ASASSN-15lh to increase strongly into the ultraviolet, with a ultraviolet luminosity a hundred times greater than the hydrogen-rich, ultraviolet-bright SLSN II SN 2008es. We find objects as bright as ASASSN-15lh are easily detectable beyond redshifts of ~4 with the single-visit depths planned for the Large Synoptic Survey Telescope. Deep near-infrared surveys could detect such objects past a redshift of ~20 enabling a probe of the earliest star formation. A late rebrightening -- most prominent at shorter wavelengths -- is seen about two months after the peak brightness, which is itself as bright as a superluminous supernova. The ultraviolet spectra during the rebrightening are dominated by the continuum without the broad absorption or emission lines seen in SLSNe or tidal disruption events and the early optical spectra of ASASSN-15lh. Our spectra show no strong hydrogen emission, showing only LyA absorption near the redshift previously found by optical absorption lines of the presumed host. The properties of ASASSN-15lh are extreme when compared to either SLSNe or tidal disruption events., Comment: submitted to AAS Journals 7 May 2016, accepted for publication in ApJ 16 June 2016, this version includes changes following the referee report and proof checking

The formation and evolution of galaxies require large reservoirs of cold, neutral gas. The damped Lya systems (DLAs), seen in absorption towards distant quasars and gamma-ray bursts, are predicted to be the dominant reservoirs for this gas. Detailed properties of DLAs have been studied extensively for decades with great success. However, their size, fundamental in understanding their nature, has remained elusive, as quasar and gamma-ray-burst sightlines only probe comparatively tiny areas of the foreground DLAs. Here, we introduce a new approach to measure the full extent of DLAs in the sightlines toward extended background sources. We present the discovery of a high-column-density (log N(HI) = 21.1 +/-0.4 cm^-2) DLA at z ~ 2.4 covering 90-100% of the luminous extent of a line-of-sight background galaxy. Estimates of the size of the background galaxy range from a minimum of a few kpc^2, to ~100 kpc^2, and demonstrate that high-column density neutral gas can span continuous areas 10^8 - 10^10 times larger than previously explored in quasar or gamma-ray burst sightlines. The DLA presented here is the first from a sample of DLAs in our pilot survey that searches Lyman break and Lyman continuum galaxies at high redshift. The low luminosities, large sizes, and mass contents (>~10^6 - 10^9 M_solar) implied by this DLA and the early data suggest that DLAs contain the necessary fuel for galaxies, with many systems consistent with relatively massive, low-luminosity primeval galaxies., Comment: 4 pages, 3 figures, 1 table, accepted to ApJ Letters