Showing 1–50 of 192 results for author: Stubbs, C W

Design and performance of a Collimated Beam Projector for telescope transmission measurement using a broadband light source

Authors: K. Sommer, J. Cohen-Tanugi, B. Plez, M. Betoule, S. Bongard, L. Le Guillou, J. Neveu, E. Nuss, E. Sepulveda, T. Souverin, M. Moniez, C. W. Stubbs

Abstract: Type Ia supernovae are the most direct cosmological probe to study dark energy in the recent Universe, for which the photometric calibration of astronomical instruments remains one major source of systematic uncertainties. To address this, recent advancements introduce Collimated Beam Projectors (CBP), aiming to enhance calibration by precisely measuring a telescope's throughput as a function of w… ▽ More Type Ia supernovae are the most direct cosmological probe to study dark energy in the recent Universe, for which the photometric calibration of astronomical instruments remains one major source of systematic uncertainties. To address this, recent advancements introduce Collimated Beam Projectors (CBP), aiming to enhance calibration by precisely measuring a telescope's throughput as a function of wavelength. This work describes the performance of a prototype portable CBP. The experimental setup consists of a broadband Xenon light source replacing a more customary but much more demanding high-power laser source, coupled with a monochromator emitting light inside an integrating sphere monitored with a photodiode and a spectrograph. Light is injected at the focus of the CBP telescope projecting a collimated beam onto a solar cell whose quantum efficiency has been obtained by comparison with a NIST-calibrated photodiode. The throughput and signal-to-noise ratio achieved by comparing the photocurrent signal in the CBP photodiode to the one in the solar cell are computed. We prove that the prototype, in its current state of development, is capable of achieving 1.2 per cent and 2.3 per cent precision on the integrated g and r bands of the ZTF photometric filter system respectively, in a reasonable amount of integration time. Central wavelength determination accuracy is kept below $\sim$ {0.91} nm and $\sim$ {0.58} nm for g and r bands. The expected photometric uncertainty caused by filter throughput measurement is approximately 5 mmag on the zero-point magnitude. Several straightforward improvement paths are discussed to upgrade the current setup. △ Less

Submitted 2 March, 2024; v1 submitted 5 December, 2023; originally announced December 2023.

Comments: submitted to RAS Techniques & Instruments (RASTI)

Directly Characterizing Dome Seeing: Differential Image Motion Sensor Using Multisources (DIMSUM)

Authors: Ali Kurmus, Elana Urbach, Christopher W. Stubbs

Abstract: Image degradation impedes our ability to extract information from astronomical observations. One factor contributing to this degradation is ``dome seeing", the reduction in image quality due to variations in the index of refraction within the observatory dome. Addressing this challenge, we introduce a novel setup-DIMSUM (Differential Image Motion Sensor Using Multisources)-which offers a simple in… ▽ More Image degradation impedes our ability to extract information from astronomical observations. One factor contributing to this degradation is ``dome seeing", the reduction in image quality due to variations in the index of refraction within the observatory dome. Addressing this challenge, we introduce a novel setup-DIMSUM (Differential Image Motion Sensor Using Multisources)-which offers a simple installation and provides direct characterization of local index of refraction variations. This is achieved by measuring differential image motion using strobed imaging that effectively``freezes" the atmosphere, aligning our captured images with the timescale of thermal fluctuations, thereby giving a more accurate representation of dome seeing effects. Our apparatus has been situated within the Auxiliary Telescope of the Vera C. Rubin Observatory. Early results from our setup are encouraging. Not only do we observe a correlation between the characteristic differential image motion (DIM) values and local temperature fluctuations (a leading cause of variations in index of refraction), but also hint at the potential of DIM measures to characterize dome seeing with greater precision in subsequent tests. Our preliminary findings underscore the potential of DIMSUM as a powerful tool for enhancing image quality in ground-based astronomical observations. Further refinement and data collection will likely solidify its place as a useful component for managing dome seeing in major observatories like the Vera C. Rubin Observatory. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: 12 pages, 7 figures, submitted to Publications of the Astronomical Society of the Pacific (PASP)

GW190425: Pan-STARRS and ATLAS coverage of the skymap and limits on optical emission associated with FRB190425

Authors: S. J. Smartt, M. Nicholl, S. Srivastav, M. E. Huber, K. C. Chambers, K. W. Smith, D. R. Young, M. D. Fulton, J. L. Tonry, C. W. Stubbs, L. Denneau, A. J. Cooper, A. Aamer, J. P. Anderson, A. Andersson, J. Bulger, T. -W Chen, P. Clark, T. de Boer, H. Gao, J. H. Gillanders, A. Lawrence, C. C. Lin, T. B. Lowe, E. A. Magnier , et al. (10 additional authors not shown)

Abstract: GW190425 is the second of only two binary neutron star (BNS) merger events to be significantly detected by the LIGO-Virgo- Kagra gravitational wave detectors. With a detection only in LIGO Livingston, the skymap containing the source was large and no plausible electromagnetic counterpart was found in real time searching in 2019. Here we summarise our ATLAS and Pan-STARRS wide-field optical coverag… ▽ More GW190425 is the second of only two binary neutron star (BNS) merger events to be significantly detected by the LIGO-Virgo- Kagra gravitational wave detectors. With a detection only in LIGO Livingston, the skymap containing the source was large and no plausible electromagnetic counterpart was found in real time searching in 2019. Here we summarise our ATLAS and Pan-STARRS wide-field optical coverage of the skymap beginning within 1 hour and 3 hours respectively of the GW190425 merger time. More recently, a potential coincidence between GW190425 and a fast radio burst FRB 190425 has been suggested, given their spatial and temporal coincidence. The smaller sky localisation area of FRB 190425 and its dispersion measure have led to the identification of a likely host galaxy, UGC 10667 at a distance of 141 +/- 10 Mpc. Our optical imaging covered the galaxy 6.0 hrs after GW190425 was detected and 3.5 hrs after the FRB 190425. No optical emission was detected and further imaging at +1.2 and +13.2 days also revealed no emission. If the FRB 190425 and GW190425 association were real, we highlight our limits on kilonova emission from a BNS merger in UGC 10667. The model for producing FRB 190425 from a BNS merger involves a supramassive magnetised neutron star spinning down by dipole emission on the timescale of hours. We show that magnetar enhanced kilonova emission is ruled out by optical upper limits. The lack of detected optical emission from a kilonova in UGC 10667 disfavours, but does not disprove, the FRB-GW link for this source. △ Less

Submitted 20 September, 2023; originally announced September 2023.

Comments: Submitted to MNRAS, 20th Sept 2023, 9 pages

All-Sky Faint DA White Dwarf Spectrophotometric Standards for Astrophysical Observatories: The Complete Sample

Authors: Tim Axelrod, Abhijit Saha, Thomas Matheson, Edward W. Olszewski, Ralph C. Bohlin, Annalisa Calamida, Jenna Claver, Susana Deustua, Jay B. Holberg, Ivan Hubeny, John W. Mackenty, Konstantin Malanchev, Gautham Narayan, Sean Points, Armin Rest, Elena Sabbi, Christopher W. Stubbs

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

Submitted 16 May, 2023; v1 submitted 12 May, 2023; originally announced May 2023.

Comments: Accepted for publication in Astrophysical Journal. Corrected error in Table 10 and associated Fig 7 in which RP and BP values for the 3 CALSPEC standards had been transposed

Fifteen years of millimeter accuracy lunar laser ranging with APOLLO: dataset characterization

Authors: James B. R. Battat, Eric Adelberger, Nicholas R. Colmenares, Megan Farrah, Daniel P. Gonzales, C. D. Hoyle, Russett J. McMillan, Thomas W. Murphy Jr., Sanchit Sabhlok, Christopher W. Stubbs

Abstract: We present data from the Apache Point Observatory Lunar Laser-ranging Operation (APOLLO) covering the 15-year span from April 2006 through the end of 2020. APOLLO measures the earth-moon separation by recording the round-trip travel time of photons from the Apache Point Observatory to five retro-reflector arrays on the moon. The APOLLO data set, combined with the 50-year archive of measurements fr… ▽ More We present data from the Apache Point Observatory Lunar Laser-ranging Operation (APOLLO) covering the 15-year span from April 2006 through the end of 2020. APOLLO measures the earth-moon separation by recording the round-trip travel time of photons from the Apache Point Observatory to five retro-reflector arrays on the moon. The APOLLO data set, combined with the 50-year archive of measurements from other lunar laser ranging (LLR) stations, can be used to probe fundamental physics such as gravity and Lorentz symmetry, as well as properties of the moon itself. We show that range measurements performed by APOLLO since 2006 have a median nightly accuracy of 1.7 mm, which is significantly better than other LLR stations. △ Less

Submitted 21 April, 2023; originally announced April 2023.

Comments: 16 pages, 9 figures

Using Elliptical Galaxy Kinematics to Compare of the Strength of Gravity in Cosmological Regions of Differing Gravitational Potential -- A First Look

Authors: Eske M. Pedersen, Christopher W. Stubbs

Abstract: Various models of modified gravity invoke ``screening'' mechanisms that are sensitive to the value of the local gravitational potential. This could have observable consequences for galaxies. These consequences might be seen by comparing two proxies for galaxy mass -- their luminosity and their internal kinematics -- as a function of local galaxy density. Motivated by this prospect, we have compare… ▽ More Various models of modified gravity invoke ``screening'' mechanisms that are sensitive to the value of the local gravitational potential. This could have observable consequences for galaxies. These consequences might be seen by comparing two proxies for galaxy mass -- their luminosity and their internal kinematics -- as a function of local galaxy density. Motivated by this prospect, we have compared the observed properties of luminous red galaxies (LRGs) inside and outside of voids in the cosmic large scale structure. We used archival measurements of line widths, luminosities, redshifts, colors, and positions of galaxies in conjunction with recent void catalogs to construct comparison LRG samples inside and outside of voids. We fitted these two samples to the well-established fundamental plane of elliptical galaxies to constrain any differences between the inferred value of the Newtonian gravitational constant G for the two samples. We obtained a null result, with an upper limit on any fractional difference in G within and outside of cosmological voids to be $α=δ$$ G/G \sim$ 40\%. This upper bound is dominated by the small-number statistics of our N $\sim $ 100 within-void LRG sample. With the caveat that environmental effects could influence various parameters such as star formation, we estimate that a 1\% statistical limit on $α$ could be attained with data from 10${^5}$ elliptical galaxies within voids. This is within the reach of future photometric and spectroscopic surveys, both of which are required to pursue this method. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Measurement of telescope transmission using a Collimated Beam Projector

Authors: Nicholas Mondrik, Michael Coughlin, Marc Betoule, Sébastien Bongard, Joseph P. Rice, Ping-Shine Shaw, Christopher W. Stubbs, John T. Woodward, LSST Dark Energy Science Collaboration

Abstract: With the increasingly large number of type Ia supernova being detected by current-generation survey telescopes, and even more expected with the upcoming Rubin Observatory Legacy Survey of Space and Time, the precision of cosmological measurements will become limited by systematic uncertainties in flux calibration rather than statistical noise. One major source of systematic error in determining SN… ▽ More With the increasingly large number of type Ia supernova being detected by current-generation survey telescopes, and even more expected with the upcoming Rubin Observatory Legacy Survey of Space and Time, the precision of cosmological measurements will become limited by systematic uncertainties in flux calibration rather than statistical noise. One major source of systematic error in determining SNe Ia color evolution (needed for distance estimation) is uncertainty in telescope transmission, both within and between surveys. We introduce here the Collimated Beam Projector (CBP), which is meant to measure a telescope transmission with collimated light. The collimated beam more closely mimics a stellar wavefront as compared to flat-field based instruments, allowing for more precise handling of systematic errors such as those from ghosting and filter angle-of-incidence dependence. As a proof of concept, we present CBP measurements of the StarDICE prototype telescope, achieving a standard (1 sigma) uncertainty of 3 % on average over the full wavelength range measured with a single beam illumination. △ Less

Submitted 22 February, 2023; originally announced February 2023.

Perfecting our set of spectrophotometric standard DA white dwarfs

Authors: A. Calamida, T. Matheson, E. W. Olszewski, A. Saha, Tim Axelrod, C. Shanahan, J. Holberg, S. Points, G. Narayan, K. Malanchev, R. Ridden-Harper, N. Gentile-Fusillo, R. Raddi, R. Bohlin, A. Rest, I. Hubeny, S. Deustua, . J. Mackenty, E. Sabbi, C. W. Stubbs

Abstract: We verified for photometric stability a set of DA white dwarfs with Hubble Space Telescope magnitudes from the near-ultraviolet to the near-infrared and ground-based spectroscopy by using time-spaced observations from the Las Cumbres Observatory network of telescopes. The initial list of 38 stars was whittled to 32 final ones which comprise a high quality set of spectrophotometric standards. These… ▽ More We verified for photometric stability a set of DA white dwarfs with Hubble Space Telescope magnitudes from the near-ultraviolet to the near-infrared and ground-based spectroscopy by using time-spaced observations from the Las Cumbres Observatory network of telescopes. The initial list of 38 stars was whittled to 32 final ones which comprise a high quality set of spectrophotometric standards. These stars are homogeneously distributed around the sky and are all fainter than r ~ 16.5 mag. Their distribution is such that at least two of them would be available to be observed from any observatory on the ground at any time at airmass less than two. Light curves and different variability indices from the Las Cumbres Observatory data were used to determine the stability of the candidate standards. When available, Pan-STARRS1, Zwicky Transient Facility and TESS data were also used to confirm the star classification. Our analysis showed that four DA white dwarfs may exhibit evidence of photometric variability, while a fifth is cooler than our established lower temperature limit, and a sixth star might be a binary. In some instances, due to the presence of faint nearby red sources, care should be used when observing a few of the spectrophotometric standards with ground-based telescopes. Light curves and finding charts for all the stars are provided. △ Less

Submitted 20 September, 2022; originally announced September 2022.

Comments: 48 pages, 50 figures, accepted for publication on ApJ

Snowmass2021 Cosmic Frontier White Paper: Rubin Observatory after LSST

Authors: Bob Blum, Seth W. Digel, Alex Drlica-Wagner, Salman Habib, Katrin Heitmann, Mustapha Ishak, Saurabh W. Jha, Steven M. Kahn, Rachel Mandelbaum, Phil Marshall, Jeffrey A. Newman, Aaron Roodman, Christopher W. Stubbs

Abstract: The Vera C. Rubin Observatory will begin the Legacy Survey of Space and Time (LSST) in 2024, spanning an area of 18,000 square degrees in six bands, with more than 800 observations of each field over ten years. The unprecedented data set will enable great advances in the study of the formation and evolution of structure and exploration of physics of the dark universe. The observations will hold cl… ▽ More The Vera C. Rubin Observatory will begin the Legacy Survey of Space and Time (LSST) in 2024, spanning an area of 18,000 square degrees in six bands, with more than 800 observations of each field over ten years. The unprecedented data set will enable great advances in the study of the formation and evolution of structure and exploration of physics of the dark universe. The observations will hold clues about the cause for the accelerated expansion of the universe and possibly the nature of dark matter. During the next decade, LSST will be able to confirm or dispute if tensions seen today in cosmological data are due to new physics. New and unexpected phenomena could confirm or disrupt our current understanding of the universe. Findings from LSST will guide the path forward post-LSST. The Rubin Observatory will still be a uniquely powerful facility even then, capable of revealing further insights into the physics of the dark universe. These could be obtained via innovative observing strategies, e.g., targeting new probes at shorter timescales than with LSST, or via modest instrumental changes, e.g., new filters, or through an entirely new instrument for the focal plane. This White Paper highlights some of the opportunities in each scenario from Rubin observations after LSST. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021

The Pantheon+ Analysis: Cosmological Constraints

Authors: Dillon Brout, Dan Scolnic, Brodie Popovic, Adam G. Riess, Joe Zuntz, Rick Kessler, Anthony Carr, Tamara M. Davis, Samuel Hinton, David Jones, W. D'Arcy Kenworthy, Erik R. Peterson, Khaled Said, Georgie Taylor, Noor Ali, Patrick Armstrong, Pranav Charvu, Arianna Dwomoh, Antonella Palmese, Helen Qu, Benjamin M. Rose, Christopher W. Stubbs, Maria Vincenzi, Charlotte M. Wood, Peter J. Brown , et al. (21 additional authors not shown)

Abstract: We present constraints on cosmological parameters from the Pantheon+ analysis of 1701 light curves of 1550 distinct Type Ia supernovae (SNe Ia) ranging in redshift from $z=0.001$ to 2.26. This work features an increased sample size, increased redshift span, and improved treatment of systematic uncertainties in comparison to the original Pantheon analysis and results in a factor of two improvement… ▽ More We present constraints on cosmological parameters from the Pantheon+ analysis of 1701 light curves of 1550 distinct Type Ia supernovae (SNe Ia) ranging in redshift from $z=0.001$ to 2.26. This work features an increased sample size, increased redshift span, and improved treatment of systematic uncertainties in comparison to the original Pantheon analysis and results in a factor of two improvement in cosmological constraining power. For a Flat$Λ$CDM model, we find $Ω_M=0.334\pm0.018$ from SNe Ia alone. For a Flat$w_0$CDM model, we measure $w_0=-0.90\pm0.14$ from SNe Ia alone, H$_0=73.5\pm1.1$ km s$^{-1}$ Mpc$^{-1}$ when including the Cepheid host distances and covariance (SH0ES), and $w_0=-0.978^{+0.024}_{-0.031}$ when combining the SN likelihood with constraints from the cosmic microwave background (CMB) and baryon acoustic oscillations (BAO); both $w_0$ values are consistent with a cosmological constant. We also present the most precise measurements to date on the evolution of dark energy in a Flat$w_0w_a$CDM universe, and measure $w_a=-0.1^{+0.9}_{-2.0}$ from Pantheon+ alone, H$_0=73.3\pm1.1$ km s$^{-1}$ Mpc$^{-1}$ when including SH0ES, and $w_a=-0.65^{+0.28}_{-0.32}$ when combining Pantheon+ with CMB and BAO data. Finally, we find that systematic uncertainties in the use of SNe Ia along the distance ladder comprise less than one third of the total uncertainty in the measurement of H$_0$ and cannot explain the present "Hubble tension" between local measurements and early-Universe predictions from the cosmological model. △ Less

Submitted 14 November, 2022; v1 submitted 8 February, 2022; originally announced February 2022.

Comments: 34 Pages, 16 Figures, 7 Tables. Published in ApJ. Comments welcome. Papers and data release here: https://pantheonplussh0es.github.io

Journal ref: ApJ 938 110 (2022)

The Pantheon+ Analysis: Dependence of Cosmological Constraints on Photometric-Zeropoint Uncertainties of Supernova Surveys

Authors: Sasha Brownsberger, Dillon Brout, Daniel Scolnic, Christopher W. Stubbs, Adam G. Riess

Abstract: Type Ia supernovae (SNe Ia) measurements of the Hubble constant, H$_0$, the cosmological mass density, $Ω_M$, and the dark energy equation-of-state parameter, $w$, rely on numerous SNe surveys using distinct photometric systems across three decades of observation. Here, we determine the sensitivities of the upcoming SH0ES+Pantheon+ constraints on H$_0$, $Ω_M$, and $w$ to unknown systematics in the… ▽ More Type Ia supernovae (SNe Ia) measurements of the Hubble constant, H$_0$, the cosmological mass density, $Ω_M$, and the dark energy equation-of-state parameter, $w$, rely on numerous SNe surveys using distinct photometric systems across three decades of observation. Here, we determine the sensitivities of the upcoming SH0ES+Pantheon+ constraints on H$_0$, $Ω_M$, and $w$ to unknown systematics in the relative photometric zeropoint calibration between the 17 surveys that comprise the Pantheon+ supernovae data set. Varying the zeropoints of these surveys simultaneously with the cosmological parameters, we determine that the SH0ES+Pantheon+ measurement of H$_0$ is robust against inter-survey photometric miscalibration, but that the measurements of $Ω_M$ and $w$ are not. Specifically, we find that miscalibrated inter-survey systematics could represent a source of uncertainty in the measured value of H$_0$ that is no larger than $0.2$ km s$^{-1}$ Mpc$^{-1}$. This modest increase in H$_0$ uncertainty could not account for the $7$ km s$^{-1}$ Mpc$^{-1}$ "Hubble Tension" between the SH0ES measurement of H$_0$ and the Planck $Λ$CDM-based inference of H$_0$. However, we find that the SH0ES+Pantheon+ best-fit values of $Ω_M$ and $w$ respectively slip, to first order, by $0.04$ and $-0.17$ per $25$ mmag of inter-survey calibration uncertainty, underscoring the vital role that cross-calibration plays in accurately measuring these parameters. Because the Pantheon+ compendium contains many surveys that share low-$z$ Hubble Flow and Cepheid-paired SNe, the SH0ES+Pantheon+ joint constraint of H$_0$ is robust against inter-survey photometric calibration errors, and such errors do not represent an impediment to jointly using SH0ES+Pantheon+ to measure H$_0$ to 1% accuracy. △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: Submitted to ApJ. Comments welcome

Strobed Imaging as a Method for the Determination and Diagnosis of Local Seeing

Authors: Christopher W. Stubbs

Abstract: The image quality budget of many telescopes can have substantial contributions from local seeing, both``mirror'' and ``dome'', which arise from turbulence and temperature variations that are difficult to quantify, measure directly, and ameliorate. We describe a method to determine the ``local'' seeing degradation due to wavefront perturbations within the final tens of meters of the optical path fr… ▽ More The image quality budget of many telescopes can have substantial contributions from local seeing, both``mirror'' and ``dome'', which arise from turbulence and temperature variations that are difficult to quantify, measure directly, and ameliorate. We describe a method to determine the ``local'' seeing degradation due to wavefront perturbations within the final tens of meters of the optical path from celestial sources to a ground-based telescope, using the primary instrument and along the same path taken by light from celestial sources. The concept involves placing strobed emitters along the light path to produce images on the main focal plane that ``freeze'' different realizations of index perturbations. This method has the advantage of measuring directly the image motion and scintillation imparted by the dynamic spatial and temporal structure of local perturbations in the index of refraction along the light path, with a clean separation from seeing induced in the atmosphere above the dome. The strobed-source approach allows for rapid image motion and scintillation to be measured directly on the focal plane, even for large-aperture telescopes with wide field instruments and slow shutters, such as that being constructed for the Rubin Observatory. A conceptual design is presented that uses the ``guider'' CCDs in the Rubin telescope focal plane to make local-seeing measurements on demand, perhaps even during science exposures. △ Less

Submitted 19 June, 2021; originally announced June 2021.

Comments: 6 pages, 3 figures, submitted to MNRAS

The LSST DESC DC2 Simulated Sky Survey

Authors: LSST Dark Energy Science Collaboration, Bela Abolfathi, David Alonso, Robert Armstrong, Éric Aubourg, Humna Awan, Yadu N. Babuji, Franz Erik Bauer, Rachel Bean, George Beckett, Rahul Biswas, Joanne R. Bogart, Dominique Boutigny, Kyle Chard, James Chiang, Chuck F. Claver, Johann Cohen-Tanugi, Céline Combet, Andrew J. Connolly, Scott F. Daniel, Seth W. Digel, Alex Drlica-Wagner, Richard Dubois, Emmanuel Gangler, Eric Gawiser , et al. (55 additional authors not shown)

Abstract: We describe the simulated sky survey underlying the second data challenge (DC2) carried out in preparation for analysis of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) by the LSST Dark Energy Science Collaboration (LSST DESC). Significant connections across multiple science domains will be a hallmark of LSST; the DC2 program represents a unique modeling effort that stresses… ▽ More We describe the simulated sky survey underlying the second data challenge (DC2) carried out in preparation for analysis of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) by the LSST Dark Energy Science Collaboration (LSST DESC). Significant connections across multiple science domains will be a hallmark of LSST; the DC2 program represents a unique modeling effort that stresses this interconnectivity in a way that has not been attempted before. This effort encompasses a full end-to-end approach: starting from a large N-body simulation, through setting up LSST-like observations including realistic cadences, through image simulations, and finally processing with Rubin's LSST Science Pipelines. This last step ensures that we generate data products resembling those to be delivered by the Rubin Observatory as closely as is currently possible. The simulated DC2 sky survey covers six optical bands in a wide-fast-deep (WFD) area of approximately 300 deg^2 as well as a deep drilling field (DDF) of approximately 1 deg^2. We simulate 5 years of the planned 10-year survey. The DC2 sky survey has multiple purposes. First, the LSST DESC working groups can use the dataset to develop a range of DESC analysis pipelines to prepare for the advent of actual data. Second, it serves as a realistic testbed for the image processing software under development for LSST by the Rubin Observatory. In particular, simulated data provide a controlled way to investigate certain image-level systematic effects. Finally, the DC2 sky survey enables the exploration of new scientific ideas in both static and time-domain cosmology. △ Less

Submitted 26 January, 2021; v1 submitted 12 October, 2020; originally announced October 2020.

Comments: 39 pages, 19 figures, version accepted for publication in ApJS

Design and operation of the ATLAS Transient Science Server

Authors: K. W. Smith, S. J. Smartt, D. R. Young, J. L. Tonry, L. Denneau, H. Flewelling, A. N. Heinze, H. J. Weiland, B. Stalder, A. Rest, C. W. Stubbs, J. P. Anderson, T. -W. Chen, P. Clark, A. Do, F. Förster, M. Fulton, J. Gillanders, O. R. McBrien, D. O'Neill, S. Srivastav, D. E. Wright

Abstract: The Asteroid Terrestrial impact Last Alert System (ATLAS) system consists of two 0.5m Schmidt telescopes with cameras covering 29 square degrees at plate scale of 1.86 arcsec per pixel. Working in tandem, the telescopes routinely survey the whole sky visible from Hawaii (above $δ> -50^{\circ}$) every two nights, exposing four times per night, typically reaching $o < 19$ magnitude per exposure when… ▽ More The Asteroid Terrestrial impact Last Alert System (ATLAS) system consists of two 0.5m Schmidt telescopes with cameras covering 29 square degrees at plate scale of 1.86 arcsec per pixel. Working in tandem, the telescopes routinely survey the whole sky visible from Hawaii (above $δ> -50^{\circ}$) every two nights, exposing four times per night, typically reaching $o < 19$ magnitude per exposure when the moon is illuminated and $c < 19.5$ per exposure in dark skies. Construction is underway of two further units to be sited in Chile and South Africa which will result in an all-sky daily cadence from 2021. Initially designed for detecting potentially hazardous near earth objects, the ATLAS data enable a range of astrophysical time domain science. To extract transients from the data stream requires a computing system to process the data, assimilate detections in time and space and associate them with known astrophysical sources. Here we describe the hardware and software infrastructure to produce a stream of clean, real, astrophysical transients in real time. This involves machine learning and boosted decision tree algorithms to identify extragalactic and Galactic transients. Typically we detect 10-15 supernova candidates per night which we immediately announce publicly. The ATLAS discoveries not only enable rapid follow-up of interesting sources but will provide complete statistical samples within the local volume of 100 Mpc. A simple comparison of the detected supernova rate within 100 Mpc, with no corrections for completeness, is already significantly higher (factor 1.5 to 2) than the current accepted rates. △ Less

Submitted 2 June, 2020; v1 submitted 19 March, 2020; originally announced March 2020.

Comments: 27 pages, 12 figures. Accepted for publication in PASP on 2020 May 15

Observational constraints on the optical and near-infrared emission from the neutron star-black hole binary merger S190814bv

Authors: K. Ackley, L. Amati, C. Barbieri, F. E. Bauer, S. Benetti, M. G. Bernardini, K. Bhirombhakdi, M. T. Botticella, M. Branchesi, E. Brocato, S. H. Bruun, M. Bulla, S. Campana, E. Cappellaro, A. J. Castro-Tirado, K. C. Chambers, S. Chaty, T. -W. Chen, R. Ciolfi, A. Coleiro, C. M. Copperwheat, S. Covino, R. Cutter, F. D'Ammando, P. D'Avanzo , et al. (129 additional authors not shown)

Abstract: On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. Preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope (ENGRAVE) collaboration members carried out an intensive multi-… ▽ More On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. Preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope (ENGRAVE) collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical/near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS and VINROUGE projects also carried out a search on this event. Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN) possibly generated by this NS-BH merger, and for the strategy of future searches. Altogether, our observations allow us to exclude a KN with large ejecta mass $M\gtrsim 0.1\,\mathrm{M_\odot}$ to a high ($>90\%$) confidence, and we can exclude much smaller masses in a subsample of our observations. This disfavours the tidal disruption of the neutron star during the merger. Despite the sensitive instruments involved in the campaign, given the distance of S190814bv we could not reach sufficiently deep limits to constrain a KN comparable in luminosity to AT 2017gfo on a large fraction of the localisation probability. This suggests that future (likely common) events at a few hundreds Mpc will be detected only by large facilities with both high sensitivity and large field of view. Galaxy-targeted observations can reach the needed depth over a relevant portion of the localisation probability with a smaller investment of resources, but the number of galaxies to be targeted in order to get a fairly complete coverage is large, even in the case of a localisation as good as that of this event. △ Less

Submitted 22 June, 2020; v1 submitted 5 February, 2020; originally announced February 2020.

Comments: 52 pages, revised version now accepted for publication in A&A. Abstract abridged to meet arXiv requirements

Journal ref: A&A 643, A113 (2020)

Initial Assessment of Monocrystalline Silicon Solar Cells as Large-Area Sensors for Precise Flux Calibration

Authors: Sasha Brownsberger, Nicholas Mondrik, Christopher W. Stubbs

Abstract: As the precision frontier of large-area survey astrophysics advances towards the one millimagnitude level, flux calibration of astronomical instrumentation remains an ongoing challenge. We describe initial testing of silicon solar cells as large-aperture precise calibration photodiodes. We present measurements of dark current, linearity, frequency response, spatial response uniformity, and noise c… ▽ More As the precision frontier of large-area survey astrophysics advances towards the one millimagnitude level, flux calibration of astronomical instrumentation remains an ongoing challenge. We describe initial testing of silicon solar cells as large-aperture precise calibration photodiodes. We present measurements of dark current, linearity, frequency response, spatial response uniformity, and noise characteristics of the Sunpower C60 solar cells, an interdigitated back-contact 125mm x 125mm monocrystalline solar cell. We find that these devices hold considerable promise as large-area flux calibration sensors and warrant further characterization. △ Less

Submitted 29 September, 2019; originally announced September 2019.

Report on LSST Next-generation Instrumentation Workshop, April 11, 12 2019

Authors: Christopher W. Stubbs, Katrin Heitmann

Abstract: The Large Synoptic Survey Telescope (LSST) is a wide-field imaging system of unprecedented etendue. The initial goal of the project is to carry out a ten year imaging survey in six broad passbands (ugrizy) that cover $350 nm < λ< 1.1 μm$. This document reports on the discussions that occurred at workshop (held April 11-12, 2019 at Argonne National Laboratory) that was convened to explore concepts… ▽ More The Large Synoptic Survey Telescope (LSST) is a wide-field imaging system of unprecedented etendue. The initial goal of the project is to carry out a ten year imaging survey in six broad passbands (ugrizy) that cover $350 nm < λ< 1.1 μm$. This document reports on the discussions that occurred at workshop (held April 11-12, 2019 at Argonne National Laboratory) that was convened to explore concepts for using the LSST system once the initial survey is complete. Participants discussed the tradeoffs in science performance, cost, and uniqueness for i) imaging surveys using the initial wide-field CCD instrument, perhaps supplemented with different filters, ii) replacing the focal plane with some alternative sensor technology, and iii) converting the LSST system to a wide-field multi-object fiber-fed spectrograph. Participants concluded that the fiber spectrograph option would be most effective if the focal plane were to feed upwards of 30,000 fibers. Thermal management and power considerations in the LSST instrument barrel make it difficult to accommodate infrared sensors that require very low operating temperatures, and the current generation of buttable IR sensors that would extend sensitivity to 2 $μ$m are, we concluded, cost-prohibitive. Procuring and using an alternative filter set, on the other hand, is modest in cost, would take full advantage of the LSST image reduction pipeline, and could yield considerable additional information. △ Less

Submitted 12 May, 2019; originally announced May 2019.

Cluster Cosmology Constraints from the 2500 deg$^2$ SPT-SZ Survey: Inclusion of Weak Gravitational Lensing Data from Magellan and the Hubble Space Telescope

Authors: S. Bocquet, J. P. Dietrich, T. Schrabback, L. E. Bleem, M. Klein, S. W. Allen, D. E. Applegate, M. L. N. Ashby, M. Bautz, M. Bayliss, B. A. Benson, M. Brodwin, E. Bulbul, R. E. A. Canning, R. Capasso, J. E. Carlstrom, C. L. Chang, I. Chiu, H-M. Cho, A. Clocchiatti, T. M. Crawford, A. T. Crites, T. de Haan, S. Desai, M. A. Dobbs , et al. (55 additional authors not shown)

Abstract: We derive cosmological constraints using a galaxy cluster sample selected from the 2500~deg$^2$ SPT-SZ survey. The sample spans the redshift range $0.25< z<1.75$ and contains 343 clusters with SZ detection significance $ξ>5$. The sample is supplemented with optical weak gravitational lensing measurements of 32 clusters with $0.29<z<1.13$ (from Magellan and HST) and X-ray measurements of 89 cluster… ▽ More We derive cosmological constraints using a galaxy cluster sample selected from the 2500~deg$^2$ SPT-SZ survey. The sample spans the redshift range $0.25< z<1.75$ and contains 343 clusters with SZ detection significance $ξ>5$. The sample is supplemented with optical weak gravitational lensing measurements of 32 clusters with $0.29<z<1.13$ (from Magellan and HST) and X-ray measurements of 89 clusters with $0.25<z<1.75$ (from Chandra). We rely on minimal modeling assumptions: i) weak lensing provides an accurate means of measuring halo masses, ii) the mean SZ and X-ray observables are related to the true halo mass through power-law relations in mass and dimensionless Hubble parameter $E(z)$ with a-priori unknown parameters, iii) there is (correlated, lognormal) intrinsic scatter and measurement noise relating these observables to their mean relations. We simultaneously fit for these astrophysical modeling parameters and for cosmology. Assuming a flat $νΛ$CDM model, in which the sum of neutrino masses is a free parameter, we measure $Ω_\mathrm{m}=0.276\pm0.047$, $σ_8=0.781\pm0.037$, and $σ_8(Ω_\mathrm{m}/0.3)^{0.2}=0.766\pm0.025$. The redshift evolution of the X-ray $Y_\mathrm{X}$-mass and $M_\mathrm{gas}$-mass relations are both consistent with self-similar evolution to within $1σ$. The mass-slope of the $Y_\mathrm{X}$-mass relation shows a $2.3σ$ deviation from self-similarity. Similarly, the mass-slope of the $M_\mathrm{gas}$-mass relation is steeper than self-similarity at the $2.5σ$ level. In a $νw$CDM cosmology, we measure the dark energy equation of state parameter $w=-1.55\pm0.41$ from the cluster data. We perform a measurement of the growth of structure since redshift $z\sim1.7$ and find no evidence for tension with the prediction from General Relativity. We provide updated redshift and mass estimates for the SPT sample. (abridged) △ Less

Submitted 20 May, 2019; v1 submitted 4 December, 2018; originally announced December 2018.

Comments: Accepted for publication in ApJ (v2 is accepted version), the catalog can be found at https://pole.uchicago.edu/public/data/sptsz-clusters/

Photometry and spectroscopy of faint candidate spectrophotometric standard DA white dwarfs

Authors: A. Calamida, T. Matheson, A. Saha, E. Olszewski, G. Narayan, J. Claver, C. Shanahan, J. Holberg, T. Axelrod, R. Bohlin, C. W. Stubbs, S. Deustua, I. Hubeny, J. Mackenty, S. Points, A. Rest, E. Sabbi

Abstract: We present precise photometry and spectroscopy for 23 candidate spectrophotometric standard white dwarfs. The selected stars are distributed in the Northern hemisphere and around the celestial equators and are all fainter than r ~ 16.5 mag. This network of stars, when established as standards, together with the three Hubble Space Telescope primary CALSPEC white dwarfs, will provide a set of spectr… ▽ More We present precise photometry and spectroscopy for 23 candidate spectrophotometric standard white dwarfs. The selected stars are distributed in the Northern hemisphere and around the celestial equators and are all fainter than r ~ 16.5 mag. This network of stars, when established as standards, together with the three Hubble Space Telescope primary CALSPEC white dwarfs, will provide a set of spectrophotometric standards to directly calibrate data products to better than 1%. These new faint standard white dwarfs will have enough signal-to-noise ratio in future deep photometric surveys and facilities to be measured accurately while still avoiding saturation in such surveys. They will also fall within the dynamic range of large telescopes and their instruments for the foreseeable future. This paper discusses the provenance of the observational data for our candidate standard stars. The comparison with models, reconciliation with reddening, and the consequent derivation of the full spectral energy density distributions for each of them is reserved for a subsequent paper. △ Less

Submitted 25 January, 2019; v1 submitted 30 November, 2018; originally announced December 2018.

Comments: 31 pages, 17 figures, 10 tables, ApJ in press (accepted on December 23rd, 2018)

Sub-percent Photometry: Faint DA White Dwarf Spectophotometric Standards for Astrophysical Observatories

Authors: Gautham Narayan, Thomas Matheson, Abhijit Saha, Tim Axelrod, Annalisa Calamida, Edward Olszewski, Jenna Claver, Kaisey S. Mandel, Ralph C. Bohlin, Jay B. Holberg, Susana Deustua, Armin Rest, Christopher W. Stubbs, Clare E. Shanahan, Amali L. Vaz, Alfredo Zenteno, Giovanni Strampelli, Ivan Hubeny, Sean Points, Elena Sabbi, John Mackenty

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

Submitted 6 February, 2019; v1 submitted 29 November, 2018; originally announced November 2018.

Comments: 46 pages, 23 figures, 8 tables, accepted for publication in ApJS

Constraining Temporal Oscillations of Cosmological Parameters Using Type Ia Supernovae

Authors: Sasha R. Brownsberger, Christopher W. Stubbs, Daniel M. Scolnic

Abstract: The existing set of type Ia supernovae (SNe Ia) is now sufficient to detect oscillatory deviations from the canonical $Λ$CDM cosmology. We determine that the Fourier spectrum of the Pantheon data set of spectroscopically well-observed SNe Ia is consistent with the predictions of $Λ$CDM. We also develop and describe two complementary techniques for using SNe Ia to constrain those alternate cosmolog… ▽ More The existing set of type Ia supernovae (SNe Ia) is now sufficient to detect oscillatory deviations from the canonical $Λ$CDM cosmology. We determine that the Fourier spectrum of the Pantheon data set of spectroscopically well-observed SNe Ia is consistent with the predictions of $Λ$CDM. We also develop and describe two complementary techniques for using SNe Ia to constrain those alternate cosmological models that predict deviations from $Λ$CDM that are oscillatory in conformal time. The first technique uses the reduced $χ^2$ statistic to determine the likelihood that the observed data would result from a given model. The second technique uses bootstrap analysis to determine the likelihood that the Fourier spectrum of a proposed model could result from statistical fluctuations around $Λ$CDM. We constrain three oscillatory alternate cosmological models: one in which the dark energy equation of state parameter oscillates around the canonical value of $w_Λ = -1$, one in which the energy density of dark energy oscillates around its $Λ$CDM value, and one in which gravity derives from a scalar field evolving under an oscillatory coupling. We further determine that any alternate cosmological model that produces distance modulus residuals with a Fourier amplitude of $\simeq 36$ millimags is strongly ruled out, given the existing data, for frequencies between $\simeq 0.08\ \textrm{Gyr}^ {-1} h_{100}$ and $\simeq 80\ \textrm{Gyr}^ {-1} h_{100}$. △ Less

Submitted 19 April, 2019; v1 submitted 20 June, 2018; originally announced June 2018.

Journal ref: ApJ, 875, 34 (2019)

Testing of the LSST's photometric calibration strategy at the CTIO 0.9 meter telescope

Authors: Michael W. Coughlin, Susana Deustua, Augustin Guyonnet, Nicholas Mondrik, Joseph P. Rice, Christopher W. Stubbs, John T. Woodward

Abstract: The calibration hardware system of the Large Synoptic Survey Telescope (LSST) is designed to measure two quantities: a telescope's instrumental response and atmospheric transmission, both as a function of wavelength. First of all, a "collimated beam projector" is designed to measure the instrumental response function by projecting monochromatic light through a mask and a collimating optic onto the… ▽ More The calibration hardware system of the Large Synoptic Survey Telescope (LSST) is designed to measure two quantities: a telescope's instrumental response and atmospheric transmission, both as a function of wavelength. First of all, a "collimated beam projector" is designed to measure the instrumental response function by projecting monochromatic light through a mask and a collimating optic onto the telescope. During the measurement, the light level is monitored with a NIST-traceable photodiode. This method does not suffer from stray light effects or the reflections (known as ghosting) present when using a flat-field screen illumination, which has a systematic source of uncertainty from uncontrolled reflections. It allows for an independent measurement of the throughput of the telescope's optical train as well as each filter's transmission as a function of position on the primary mirror. Second, CALSPEC stars can be used as calibrated light sources to illuminate the atmosphere and measure its transmission. To measure the atmosphere's transfer function, we use the telescope's imager with a Ronchi grating in place of a filter to configure it as a low resolution slitless spectrograph. In this paper, we describe this calibration strategy, focusing on results from a prototype system at the Cerro Tololo Inter-American Observatory (CTIO) 0.9 meter telescope. We compare the instrumental throughput measurements to nominal values measured using a laboratory spectrophotometer, and we describe measurements of the atmosphere made via CALSPEC standard stars during the same run. △ Less

Submitted 6 June, 2018; originally announced June 2018.

Observational implications of lowering the LIGO-Virgo alert threshold

Authors: Ryan Lynch, Michael Coughlin, Salvatore Vitale, Christopher W. Stubbs, Erik Katsavounidis

Abstract: The recent detection of the binary-neutron-star merger associated with GW170817 by both LIGO-Virgo and the network of electromagnetic-spectrum observing facilities around the world has made the multi-messenger detection of gravitational-wave events a reality. These joint detections allow us to probe gravitational-wave sources in greater detail and provide us with the possibility of confidently est… ▽ More The recent detection of the binary-neutron-star merger associated with GW170817 by both LIGO-Virgo and the network of electromagnetic-spectrum observing facilities around the world has made the multi-messenger detection of gravitational-wave events a reality. These joint detections allow us to probe gravitational-wave sources in greater detail and provide us with the possibility of confidently establishing events that would not have been detected in gravitational-wave data alone. In this paper, we explore the prospects of using the electromagnetic follow-up of low-significance gravitational-wave event candidates to increase the sample of confident detections with electromagnetic counterparts. We find that the gravitational-wave alert threshold change that would roughly double the number of detectable astrophysical events would increase the false-alarm rate by more than 5 orders of magnitude from 1 per 100 years to more than 1000 per year. We find that the localization costs of following-up low-significance candidates are marginal, as the same changes to false-alarm rate only increase distance/area localizations by less than a factor of 2 and increase volume localization by less than a factor of 4. We argue that EM follow-up thresholds for low-significance candidates should be set on the basis of alert purity ($P_\text{astro}$) and not false-alarm rate. Ideally, such estimates of $P_\text{astro}$ would be provided by LIGO-Virgo, but in their absence we provide estimates of the average purity of the gravitational-wave candidate alerts issued by LIGO-Virgo as a function of false-alarm rate for various LIGO-Virgo observing epochs. △ Less

Submitted 26 June, 2018; v1 submitted 7 March, 2018; originally announced March 2018.

ATLAS: A High-Cadence All-Sky Survey System

Authors: J. L. Tonry, L. Denneau, A. N. Heinze, B. Stalder, K. W. Smith, S. J. Smartt, C. W. Stubbs, H. J. Weiland, A. Rest

Abstract: Technology has advanced to the point that it is possible to image the entire sky every night and process the data in real time. The sky is hardly static: many interesting phenomena occur, including variable stationary objects such as stars or QSOs, transient stationary objects such as supernovae or M dwarf flares, and moving objects such as asteroids and the stars themselves. Funded by NASA, we ha… ▽ More Technology has advanced to the point that it is possible to image the entire sky every night and process the data in real time. The sky is hardly static: many interesting phenomena occur, including variable stationary objects such as stars or QSOs, transient stationary objects such as supernovae or M dwarf flares, and moving objects such as asteroids and the stars themselves. Funded by NASA, we have designed and built a sky survey system for the purpose of finding dangerous near-Earth asteroids (NEAs). This system, the "Asteroid Terrestrial-impact Last Alert System" (ATLAS), has been optimized to produce the best survey capability per unit cost, and therefore is an efficient and competitive system for finding potentially hazardous asteroids (PHAs) but also for tracking variables and finding transients. While carrying out its NASA mission, ATLAS now discovers more bright ($m < 19$) supernovae candidates than any ground based survey, frequently detecting very young explosions due to its 2 day cadence. ATLAS discovered the afterglow of a gamma-ray burst independent of the high energy trigger and has released a variable star catalogue of 5$\times10^{6}$ sources. This, the first of a series of articles describing ATLAS, is devoted to the design and performance of the ATLAS system. Subsequent articles will describe in more detail the software, the survey strategy, ATLAS-derived NEA population statistics, transient detections, and the first data release of variable stars and transient lightcurves. △ Less

Submitted 28 March, 2018; v1 submitted 2 February, 2018; originally announced February 2018.

Comments: 26 pages, 12 figures, submitted to PASP

A kilonova as the electromagnetic counterpart to a gravitational-wave source

Authors: S. J. Smartt, T. -W. Chen, A. Jerkstrand, M. Coughlin, E. Kankare, S. A. Sim, M. Fraser, C. Inserra, K. Maguire, K. C. Chambers, M. E. Huber, T. Kruhler, G. Leloudas, M. Magee, L. J. Shingles, K. W. Smith, D. R. Young, J. Tonry, R. Kotak, A. Gal-Yam, J. D. Lyman, D. S. Homan, C. Agliozzo, J. P. Anderson, C. R. Angus C. Ashall , et al. (96 additional authors not shown)

Abstract: Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a r… ▽ More Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a relatively well confined sky position and distance estimate. Here we report observations and physical modelling of a rapidly fading electromagnetic transient in the galaxy NGC4993, which is spatially coincident with GW170817 and a weak short gamma-ray burst. The transient has physical parameters broadly matching the theoretical predictions of blue kilonovae from neutron star mergers. The emitted electromagnetic radiation can be explained with an ejected mass of 0.04 +/- 0.01 Msol, with an opacity of kappa <= 0.5 cm2/gm at a velocity of 0.2 +/- 0.1c. The power source is constrained to have a power law slope of beta = -1.2 +/- 0.3, consistent with radioactive powering from r-process nuclides. We identify line features in the spectra that are consistent with light r-process elements (90 < A < 140). As it fades, the transient rapidly becomes red, and emission may have contribution by a higher opacity, lanthanide-rich ejecta component. This indicates that neutron star mergers produce gravitational waves, radioactively powered kilonovae, and are a nucleosynthetic source of the r-process elements. △ Less

Submitted 17 October, 2017; v1 submitted 16 October, 2017; originally announced October 2017.

Comments: Nature, in press, DOI 10.1038/nature24303. Data files will be made available at http://www.pessto.org

The Complete Light-curve Sample of Spectroscopically Confirmed Type Ia Supernovae from Pan-STARRS1 and Cosmological Constraints from The Combined Pantheon Sample

Authors: D. M. Scolnic, D. O. Jones, A. Rest, Y. C. Pan, R. Chornock, R. J. Foley, M. E. Huber, R. Kessler, G. Narayan, A. G. Riess, S. Rodney, E. Berger, D. J. Brout, P. J. Challis, M. Drout, D. Finkbeiner, R. Lunnan, R. P. Kirshner, N. E. Sanders, E. Schlafly, S. Smartt, C. W. Stubbs, J. Tonry, W. M. Wood-Vasey, M. Foley , et al. (15 additional authors not shown)

Abstract: We present optical light curves, redshifts, and classifications for 365 spectroscopically confirmed Type Ia supernovae (SNe Ia) discovered by the Pan-STARRS1 (PS1) Medium Deep Survey. We detail improvements to the PS1 SN photometry, astrometry and calibration that reduce the systematic uncertainties in the PS1 SN Ia distances. We combine the subset of 279 PS1 SN Ia ($0.03 < z < 0.68$) with useful… ▽ More We present optical light curves, redshifts, and classifications for 365 spectroscopically confirmed Type Ia supernovae (SNe Ia) discovered by the Pan-STARRS1 (PS1) Medium Deep Survey. We detail improvements to the PS1 SN photometry, astrometry and calibration that reduce the systematic uncertainties in the PS1 SN Ia distances. We combine the subset of 279 PS1 SN Ia ($0.03 < z < 0.68$) with useful distance estimates of SN Ia from SDSS, SNLS, various low-z and HST samples to form the largest combined sample of SN Ia consisting of a total of 1048 SN Ia ranging from $0.01 < z < 2.3$, which we call the `Pantheon Sample'. When combining Planck 2015 CMB measurements with the Pantheon SN sample, we find $Ω_m=0.307\pm0.012$ and $w = -1.026\pm0.041$ for the wCDM model. When the SN and CMB constraints are combined with constraints from BAO and local H0 measurements, the analysis yields the most precise measurement of dark energy to date: $w0 = -1.007\pm 0.089$ and $wa = -0.222 \pm0.407$ for the w0waCDM model. Tension with a cosmological constant previously seen in an analysis of PS1 and low-z SNe has diminished after an increase of $2\times$ in the statistics of the PS1 sample, improved calibration and photometry, and stricter light-curve quality cuts. We find the systematic uncertainties in our measurements of dark energy are almost as large as the statistical uncertainties, primarily due to limitations of modeling the low-redshift sample. This must be addressed for future progress in using SN Ia to measure dark energy. △ Less

Submitted 25 March, 2018; v1 submitted 2 October, 2017; originally announced October 2017.

Comments: Accepted by ApJ. Data can be found here: http://dx.DOI.org/10.17909/T95Q4X

An absolute calibration system for millimeter-accuracy APOLLO measurements

Authors: E. G. Adelberger, J. B. R. Battat, K. J. Birkmeier, N. R. Colmenares, R. Davis, C. D. Hoyle, L. R. Huang, R. J. McMillan, T. W. Murphy Jr., E. Schlerman, C. Skrobol, C. W. Stubbs, A. Zach

Abstract: Lunar laser ranging provides a number of leading experimental tests of gravitation -- important in our quest to unify General Relativity and the Standard Model of physics. The Apache Point Observatory Lunar Laser-ranging Operation (APOLLO) has for years achieved median range precision at the ~2 mm level. Yet residuals in model-measurement comparisons are an order-of-magnitude larger, raising the q… ▽ More Lunar laser ranging provides a number of leading experimental tests of gravitation -- important in our quest to unify General Relativity and the Standard Model of physics. The Apache Point Observatory Lunar Laser-ranging Operation (APOLLO) has for years achieved median range precision at the ~2 mm level. Yet residuals in model-measurement comparisons are an order-of-magnitude larger, raising the question of whether the ranging data are not nearly as accurate as they are precise, or if the models are incomplete or ill-conditioned. This paper describes a new absolute calibration system (ACS) intended both as a tool for exposing and eliminating sources of systematic error, and also as a means to directly calibrate ranging data in-situ. The system consists of a high-repetition-rate (80 MHz) laser emitting short (< 10 ps) pulses that are locked to a cesium clock. In essence, the ACS delivers photons to the APOLLO detector at exquisitely well-defined time intervals as a "truth" input against which APOLLO's timing performance may be judged and corrected. Preliminary analysis indicates no inaccuracies in APOLLO data beyond the ~3 mm level, suggesting that historical APOLLO data are of high quality and motivating continued work on model capabilities. The ACS provides the means to deliver APOLLO data both accurate and precise below the 2 mm level. △ Less

Submitted 28 June, 2017; originally announced June 2017.

Comments: 21 pages, 10 figures, submitted to Classical and Quantum Gravity

Observations of the GRB afterglow ATLAS17aeu and its possible association with GW170104

Authors: B. Stalder, J. Tonry, S. J. Smartt, M. Coughlin, K. C. Chambers, C. W. Stubbs, T. -W. Chen, E. Kankare, K. W. Smith, L. Denneau, A. Sherstyuk, A. Heinze, H. Weiland, A. Rest, D. R. Young, M. E. Huber, H. Flewelling, T. Lowe, E. A. Magnier, A. S. B. Schultz, C. Waters, R. Wainscoat, M. Willman, D. E. Wright, J. K. Chu , et al. (4 additional authors not shown)

Abstract: We report the discovery and multi-wavelength data analysis of the peculiar optical transient, ATLAS17aeu. This transient was identified in the skymap of the LIGO gravitational wave event GW170104 by our ATLAS and Pan-STARRS coverage. ATLAS17aeu was discovered 23.1hrs after GW170104 and rapidly faded over the next 3 nights, with a spectrum revealing a blue featureless continuum. The transient was a… ▽ More We report the discovery and multi-wavelength data analysis of the peculiar optical transient, ATLAS17aeu. This transient was identified in the skymap of the LIGO gravitational wave event GW170104 by our ATLAS and Pan-STARRS coverage. ATLAS17aeu was discovered 23.1hrs after GW170104 and rapidly faded over the next 3 nights, with a spectrum revealing a blue featureless continuum. The transient was also detected as a fading x-ray source by Swift and in the radio at 6 and 15 GHz. A gamma ray burst GRB170105A was detected by 3 satellites 19.04hrs after GW170104 and 4.10hrs before our first optical detection. We analyse the multi-wavelength fluxes in the context of the known GRB population and discuss the observed sky rates of GRBs and their afterglows. We find it statistically likely that ATLAS17aeu is an afterglow associated with GRB170105A, with a chance coincidence ruled out at the 99\% confidence or 2.6$σ$. A long, soft GRB within a redshift range of $1 \lesssim z \lesssim 2.9$ would be consistent with all the observed multi-wavelength data. The Poisson probability of a chance occurrence of GW170104 and ATLAS17aeu is $p=0.04$. This is the probability of a chance coincidence in 2D sky location and in time. These observations indicate that ATLAS17aeu is plausibly a normal GRB afterglow at significantly higher redshift than the distance constraint for GW170104 and therefore a chance coincidence. However if a redshift of the faint host were to place it within the GW170104 distance range, then physical association with GW170104 should be considered. △ Less

Submitted 7 November, 2017; v1 submitted 1 June, 2017; originally announced June 2017.

Comments: 16 pages, 6 figures, accepted to ApJ

The Pan-STARRS1 Surveys

Authors: K. C. Chambers, E. A. Magnier, N. Metcalfe, H. A. Flewelling, M. E. Huber, C. Z. Waters, L. Denneau, P. W. Draper, D. Farrow, D. P. Finkbeiner, C. Holmberg, J. Koppenhoefer, P. A. Price, A. Rest, R. P. Saglia, E. F. Schlafly, S. J. Smartt, W. Sweeney, R. J. Wainscoat, W. S. Burgett, S. Chastel, T. Grav, J. N. Heasley, K. W. Hodapp, R. Jedicke , et al. (101 additional authors not shown)

Abstract: Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3π$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$σ$ point source limiting sensitivities in the stacked 3$π$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across… ▽ More Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3π$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$σ$ point source limiting sensitivities in the stacked 3$π$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($ Δra, Δdec $) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI. △ Less

Submitted 28 January, 2019; v1 submitted 16 December, 2016; originally announced December 2016.

Comments: 38 pages, 29 figures, 12 tables

Pan-STARRS Pixel Analysis : Source Detection and Characterization

Authors: Eugene A. Magnier, W. E. Sweeney, K. C. Chambers, H. A. Flewelling, M. E. Huber, P. A. Price, C. Z. Waters, L. Denneau, P. Draper, R. Jedicke, K. W. Hodapp, N. Kaiser, R. -P. Kudritzki, N. Metcalfe, C. W. Stubbs, R. J. Wainscoast

Abstract: Over 3 billion astronomical objects have been detected in the more than 22 million orthogonal transfer CCD images obtained as part of the Pan-STARRS1 $3π$ survey. Over 85 billion instances of those objects have been automatically detected and characterized by the Pan-STARRS Image Processing Pipeline photometry software, psphot. This fast, automatic, and reliable software was developed for the Pan-… ▽ More Over 3 billion astronomical objects have been detected in the more than 22 million orthogonal transfer CCD images obtained as part of the Pan-STARRS1 $3π$ survey. Over 85 billion instances of those objects have been automatically detected and characterized by the Pan-STARRS Image Processing Pipeline photometry software, psphot. This fast, automatic, and reliable software was developed for the Pan-STARRS project, but is easily adaptable to images from other telescopes. We describe the analysis of the astronomical objects by psphot in general as well as for the specific case of the 3rd processing version used for the first two public releases of the Pan-STARRS $3π$ survey data, DR1 & DR2. △ Less

Submitted 26 January, 2019; v1 submitted 15 December, 2016; originally announced December 2016.

Comments: Pan-STARRS Public Data Release 2 : Paper IV

Pan-STARRS Photometric and Astrometric Calibration

Authors: Eugene. A. Magnier, Edward. F. Schlafly, Douglas P. Finkbeiner, J. L. Tonry, B. Goldman, S. Röser, E. Schilbach, K. C. Chambers, H. A. Flewelling, M. E. Huber, P. A. Price, W. E. Sweeney, C. Z. Waters, L. Denneau, P. Draper, K. W. Hodapp, R. Jedicke, R. -P. Kudritzki, N. Metcalfe, C. W. Stubbs, R. J. Wainscoast

Abstract: We present the details of the photometric and astrometric calibration of the Pan-STARRS1 $3π$ Survey. The photometric goals were to reduce the systematic effects introduced by the camera and detectors, and to place all of the observations onto a photometric system with consistent zero points over the entire area surveyed, the ~30,000 square degrees north of $δ$ = -30 degrees. The astrometric calib… ▽ More We present the details of the photometric and astrometric calibration of the Pan-STARRS1 $3π$ Survey. The photometric goals were to reduce the systematic effects introduced by the camera and detectors, and to place all of the observations onto a photometric system with consistent zero points over the entire area surveyed, the ~30,000 square degrees north of $δ$ = -30 degrees. The astrometric calibration compensates for similar systematic effects so that positions, proper motions, and parallaxes are reliable as well. The Pan-STARRS Data Release 2 (DR2) astrometry is tied to the Gaia DR1 release. △ Less

Submitted 26 January, 2019; v1 submitted 15 December, 2016; originally announced December 2016.

Comments: Pan-STARRS Public Data Release 2 : Paper V

The Pan-STARRS Data Processing System

Authors: Eugene A. Magnier, K. C. Chambers, H. A. Flewelling, J. C. Hoblitt, M. E. Huber, P. A. Price, W. E. Sweeney, C. Z. Waters, L. Denneau, P. Draper, K. W. Hodapp, R. Jedicke, N. Kaiser, R. -P. Kudritzki, N. Metcalfe, C. W. Stubbs, R. J. Wainscoast

Abstract: The Pan-STARRS Data Processing System is responsible for the steps needed to downloaded, archive, and process all images obtained by the Pan-STARRS telescopes, including real-time detection of transient sources such as supernovae and moving objects including potentially hazardous asteroids. With a nightly data volume of up to 4 terabytes and an archive of over 4 petabytes of raw imagery, Pan-STARR… ▽ More The Pan-STARRS Data Processing System is responsible for the steps needed to downloaded, archive, and process all images obtained by the Pan-STARRS telescopes, including real-time detection of transient sources such as supernovae and moving objects including potentially hazardous asteroids. With a nightly data volume of up to 4 terabytes and an archive of over 4 petabytes of raw imagery, Pan-STARRS is solidly in the realm of Big Data astronomy. The full data processing system consists of several subsystems covering the wide range of necessary capabilities. This article describes the Image Processing Pipeline and its connections to both the summit data systems and the outward-facing systems downstream. The latter include the Moving Object Processing System (MOPS) & the public database: the Published Science Products Subsystem (PSPS). △ Less

Submitted 26 January, 2019; v1 submitted 15 December, 2016; originally announced December 2016.

Comments: Pan-STARRS Public Data Release 2 : Paper II

Velocity Segregation and Systematic Biases In Velocity Dispersion Estimates With the SPT-GMOS Spectroscopic Survey

Authors: Matthew. B. Bayliss, Kyle Zengo, Jonathan Ruel, Bradford A. Benson, Lindsey E. Bleem, Sebastian Bocquet, Esra Bulbul, Mark Brodwin, Raffaella Capasso, I-non Chiu, Michael McDonald, David Rapetti, Alex Saro, Brian Stalder, Antony A. Stark, Veronica Strazzullo, Christopher W. Stubbs, Alfredo Zenteno

Abstract: The velocity distribution of galaxies in clusters is not universal; rather, galaxies are segregated according to their spectral type and relative luminosity. We examine the velocity distributions of different populations of galaxies within 89 Sunyaev Zel'dovich (SZ) selected galaxy clusters spanning $ 0.28 < z < 1.08$. Our sample is primarily draw from the SPT-GMOS spectroscopic survey, supplement… ▽ More The velocity distribution of galaxies in clusters is not universal; rather, galaxies are segregated according to their spectral type and relative luminosity. We examine the velocity distributions of different populations of galaxies within 89 Sunyaev Zel'dovich (SZ) selected galaxy clusters spanning $ 0.28 < z < 1.08$. Our sample is primarily draw from the SPT-GMOS spectroscopic survey, supplemented by additional published spectroscopy, resulting in a final spectroscopic sample of 4148 galaxy spectra---2868 cluster members. The velocity dispersion of star-forming cluster galaxies is $17\pm4$% greater than that of passive cluster galaxies, and the velocity dispersion of bright ($m < m^{*}-0.5$) cluster galaxies is $11\pm4$% lower than the velocity dispersion of our total member population. We find good agreement with simulations regarding the shape of the relationship between the measured velocity dispersion and the fraction of passive vs. star-forming galaxies used to measure it, but we find a small offset between this relationship as measured in data and simulations in which suggests that our dispersions are systematically low by as much as 3\% relative to simulations. We argue that this offset could be interpreted as a measurement of the effective velocity bias that describes the ratio of our observed velocity dispersions and the intrinsic velocity dispersion of dark matter particles in a published simulation result. Measuring velocity bias in this way suggests that large spectroscopic surveys can improve dispersion-based mass-observable scaling relations for cosmology even in the face of velocity biases, by quantifying and ultimately calibrating them out. △ Less

Submitted 8 February, 2017; v1 submitted 8 December, 2016; originally announced December 2016.

Comments: Accepted to ApJ; 21 pages, 11 figures, 5 tables

Cluster Mass Calibration at High Redshift: HST Weak Lensing Analysis of 13 Distant Galaxy Clusters from the South Pole Telescope Sunyaev-Zel'dovich Survey

Authors: T. Schrabback, D. Applegate, J. P. Dietrich, H. Hoekstra, S. Bocquet, A. H. Gonzalez, A. von der Linden, M. McDonald, C. B. Morrison, S. F. Raihan, S. W. Allen, M. Bayliss, B. A. Benson, L. E. Bleem, I. Chiu, S. Desai, R. J. Foley, T. de Haan, F. W. High, S. Hilbert, A. B. Mantz, R. Massey, J. Mohr, C. L. Reichardt, A. Saro , et al. (4 additional authors not shown)

Abstract: We present an HST/ACS weak gravitational lensing analysis of 13 massive high-redshift (z_median=0.88) galaxy clusters discovered in the South Pole Telescope (SPT) Sunyaev-Zel'dovich Survey. This study is part of a larger campaign that aims to robustly calibrate mass-observable scaling relations over a wide range in redshift to enable improved cosmological constraints from the SPT cluster sample. W… ▽ More We present an HST/ACS weak gravitational lensing analysis of 13 massive high-redshift (z_median=0.88) galaxy clusters discovered in the South Pole Telescope (SPT) Sunyaev-Zel'dovich Survey. This study is part of a larger campaign that aims to robustly calibrate mass-observable scaling relations over a wide range in redshift to enable improved cosmological constraints from the SPT cluster sample. We introduce new strategies to ensure that systematics in the lensing analysis do not degrade constraints on cluster scaling relations significantly. First, we efficiently remove cluster members from the source sample by selecting very blue galaxies in V-I colour. Our estimate of the source redshift distribution is based on CANDELS data, where we carefully mimic the source selection criteria of the cluster fields. We apply a statistical correction for systematic photometric redshift errors as derived from Hubble Ultra Deep Field data and verified through spatial cross-correlations. We account for the impact of lensing magnification on the source redshift distribution, finding that this is particularly relevant for shallower surveys. Finally, we account for biases in the mass modelling caused by miscentring and uncertainties in the concentration-mass relation using simulations. In combination with temperature estimates from Chandra we constrain the normalisation of the mass-temperature scaling relation ln(E(z) M_500c/10^14 M_sun)=A+1.5 ln(kT/7.2keV) to A=1.81^{+0.24}_{-0.14}(stat.) +/- 0.09(sys.), consistent with self-similar redshift evolution when compared to lower redshift samples. Additionally, the lensing data constrain the average concentration of the clusters to c_200c=5.6^{+3.7}_{-1.8}. △ Less

Submitted 30 October, 2017; v1 submitted 11 November, 2016; originally announced November 2016.

Comments: 49 pages, 11 tables, 38 figures. Matches the version accepted for publication in MNRAS

SPT-GMOS: A Gemini/GMOS-South Spectroscopic Survey of Galaxy Clusters in the SPT-SZ Survey

Authors: M. B. Bayliss, J. Ruel, C. W. Stubbs, S. W. Allen, D. E. Applegate, M. L. N. Ashby, M. Bautz, B. A. Benson, L. E. Bleem, S. Bocquet, M. Brodwin, R. Capasso, J. E. Carlstrom, C. L. Chang, I. Chiu, H-M. Cho, A. Clocchiatti, T. M. Crawford, A. T. Crites, T. de Haan, S. Desai, J. P. Dietrich, M. A. Dobbs, A. N. Doucouliagos, R. J. Foley , et al. (53 additional authors not shown)

Abstract: We present the results of SPT-GMOS, a spectroscopic survey with the Gemini Multi-Object Spectrograph (GMOS) on Gemini South. The targets of SPT-GMOS are galaxy clusters identified in the SPT-SZ survey, a millimeter-wave survey of 2500 squ. deg. of the southern sky using the South Pole Telescope (SPT). Multi-object spectroscopic observations of 62 SPT-selected galaxy clusters were performed between… ▽ More We present the results of SPT-GMOS, a spectroscopic survey with the Gemini Multi-Object Spectrograph (GMOS) on Gemini South. The targets of SPT-GMOS are galaxy clusters identified in the SPT-SZ survey, a millimeter-wave survey of 2500 squ. deg. of the southern sky using the South Pole Telescope (SPT). Multi-object spectroscopic observations of 62 SPT-selected galaxy clusters were performed between January 2011 and December 2015, yielding spectra with radial velocity measurements for 2595 sources. We identify 2243 of these sources as galaxies, and 352 as stars. Of the galaxies, we identify 1579 as members of SPT-SZ galaxy clusters. The primary goal of these observations was to obtain spectra of cluster member galaxies to estimate cluster redshifts and velocity dispersions. We describe the full spectroscopic dataset and resulting data products, including galaxy redshifts, cluster redshifts and velocity dispersions, and measurements of several well-known spectral indices for each galaxy: the equivalent width, W, of [O II] 3727,3729 and H-delta, and the 4000A break strength, D4000. We use the spectral indices to classify galaxies by spectral type (i.e., passive, post-starburst, star-forming), and we match the spectra against photometric catalogs to characterize spectroscopically-observed cluster members as a function of brightness (relative to m*). Finally, we report several new measurements of redshifts for ten bright, strongly-lensed background galaxies in the cores of eight galaxy clusters. Combining the SPT-GMOS dataset with previous spectroscopic follow-up of SPT-SZ galaxy clusters results in spectroscopic measurements for >100 clusters, or ~20% of the full SPT-SZ sample. △ Less

Submitted 16 September, 2016; originally announced September 2016.

Comments: 24 pages in eapj format. Accepted to the Astrophysical Journal Supplements

Testing for X-ray-SZ Differences and Redshift Evolution in the X-ray Morphology of Galaxy Clusters

Authors: D. Nurgaliev, M. McDonald, B. A. Benson, L. Bleem, S. Bocquet, W. R. Forman, G. P. Garmire, N. Gupta, J. Hlavacek-Larrondo, J. J. Mohr, D. Nagai, D. Rapetti, A. A. Stark, C. W. Stubbs, A. Vikhlinin

Abstract: We present a quantitative study of the X-ray morphology of galaxy clusters, as a function of their detection method and redshift. We analyze two separate samples of galaxy clusters: a sample of 36 clusters at 0.35 < z < 0.9 selected in the X-ray with the ROSAT PSPC 400 deg2 survey, and a sample of 90 clusters at 0.25 < z < 1.2 selected via the Sunyaev-Zel'dovich (SZ) effect with the South Pole Tel… ▽ More We present a quantitative study of the X-ray morphology of galaxy clusters, as a function of their detection method and redshift. We analyze two separate samples of galaxy clusters: a sample of 36 clusters at 0.35 < z < 0.9 selected in the X-ray with the ROSAT PSPC 400 deg2 survey, and a sample of 90 clusters at 0.25 < z < 1.2 selected via the Sunyaev-Zel'dovich (SZ) effect with the South Pole Telescope. Clusters from both samples have similar-quality Chandra observations, which allow us to quantify their X-ray morphologies via two distinct methods: centroid shifts and photon asymmetry. The latter technique provides nearly unbiased morphology estimates for clusters spanning a broad range of redshift and data quality. We further compare the X-ray morphologies of X-ray- and SZ-selected clusters with those of simulated clusters. We do not find a statistically significant difference in the measured X-ray morphology of X-ray and SZ-selected clusters over the redshift range probed by these samples, suggesting that the two are probing similar populations of clusters. We find that the X-ray morphologies of simulated clusters are statistically indistinguishable from those of X-ray- or SZ-selected clusters, implying that the most important physics for dictating the large-scale gas morphology (outside of the core) is well-approximated in these simulations. Finally, we find no statistically significant redshift evolution in the X-ray morphology (both for observed and simulated clusters), over the range z ~ 0.3 to z ~ 1, seemingly in contradiction with the redshift-dependent halo merger rate predicted by simulations. △ Less

Submitted 29 June, 2017; v1 submitted 1 September, 2016; originally announced September 2016.

Comments: 13 pages, 5 figures, 2 tables. Published in ApJ

A search for an optical counterpart to the gravitational wave event GW151226

Authors: S. J. Smartt, K. C. Chambers, K. W. Smith, M. E. Huber, D. R. Young, T. -W. Chen, C. Inserra, D. E. Wright, M. Coughlin, L. Denneau, H. Flewelling, A. Heinze, A. Jerkstrand, E. A. Magnier, K. Maguire, B. Mueller, A. Rest, A. Sherstyuk, B. Stalder, A. S. B. Schultz, C. W. Stubbs, J. Tonry, C. Waters, R. Wainscoat, M. Della Valle , et al. (15 additional authors not shown)

Abstract: We present a search for an electromagnetic counterpart of the gravitational wave source GW151226. Using the Pan-STARRS1 telescope we mapped out 290 square degrees in the optical i_ps filter starting 11.5hr after the LIGO information release and lasting for a further 28 days. The first observations started 49.5hr after the time of the GW151226 detection. We typically reached sensitivity limits of i… ▽ More We present a search for an electromagnetic counterpart of the gravitational wave source GW151226. Using the Pan-STARRS1 telescope we mapped out 290 square degrees in the optical i_ps filter starting 11.5hr after the LIGO information release and lasting for a further 28 days. The first observations started 49.5hr after the time of the GW151226 detection. We typically reached sensitivity limits of i_ps = 20.3-20.8 and covered 26.5% of the LIGO probability skymap. We supplemented this with ATLAS survey data, reaching 31% of the probability region to shallower depths of m~19. We found 49 extragalactic transients (that are not obviously AGN), including a faint transient in a galaxy at 7Mpc (a luminous blue variable outburst) plus a rapidly decaying M-dwarf flare. Spectral classification of 20 other transient events showed them all to be supernovae. We found an unusual transient, PS15dpn, with an explosion date temporally coincident with GW151226 which evolved into a type Ibn supernova. The redshift of the transient is secure at z=0.1747 +/- 0.0001 and we find it unlikely to be linked, since the luminosity distance has a negligible probability of being consistent with that of GW151226. In the 290 square degrees surveyed we therefore do not find a likely counterpart. However we show that our survey strategy would be sensitive to NS-NS mergers producing kilonovae at D < 100 Mpc which is promising for future LIGO/Virgo searches. △ Less

Submitted 8 August, 2016; v1 submitted 15 June, 2016; originally announced June 2016.

Comments: Submitted to ApJL on 6th June 2016. Accepted on 8th August 2016. This is the accepted version after referee report and minor revisions

Supplement: Localization and broadband follow-up of the gravitational-wave transient GW150914

Authors: B. P. Abbott, R. Abbott, T. D. Abbott, M. R. Abernathy, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, A. Allocca, P. A. Altin, S. B. Anderson, W. G. Anderson, K. Arai , et al. (1522 additional authors not shown)

Abstract: This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the dif… ▽ More This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands. △ Less

Submitted 21 July, 2016; v1 submitted 26 April, 2016; originally announced April 2016.

Comments: For the main Letter, see arXiv:1602.08492

Report number: LIGO-P1600137-v2

Journal ref: The Astrophysical Journal Supplement Series, 225:8 (15pp), 2016 July

Maximizing the Probability of Detecting an Electromagnetic Counterpart of Gravitational-wave Events

Authors: Michael W. Coughlin, Christopher W. Stubbs

Abstract: Compact binary coalescences are a promising source of gravitational waves for second-generation interferometric gravitational-wave detectors such as advanced LIGO and advanced Virgo. These are among the most promising sources for joint detection of electromagnetic (EM) and gravitational-wave (GW) emission. To maximize the science performed with these objects, it is essential to undertake a followu… ▽ More Compact binary coalescences are a promising source of gravitational waves for second-generation interferometric gravitational-wave detectors such as advanced LIGO and advanced Virgo. These are among the most promising sources for joint detection of electromagnetic (EM) and gravitational-wave (GW) emission. To maximize the science performed with these objects, it is essential to undertake a followup observing strategy that maximizes the likelihood of detecting the EM counterpart. We present a follow-up strategy that maximizes the counterpart detection probability, given a fixed investment of telescope time. We show how the prior assumption on the luminosity function of the electro-magnetic counterpart impacts the optimized followup strategy. Our results suggest that if the goal is to detect an EM counterpart from among a succession of GW triggers, the optimal strategy is to perform long integrations in the highest likelihood regions, with a time investment that is proportional to the $2/3$ power of the surface density of the GW location probability on the sky. In the future, this analysis framework will benefit significantly from the 3-dimensional localization probability. △ Less

Submitted 18 April, 2016; originally announced April 2016.

Cosmological Constraints from Galaxy Clusters in the 2500 square-degree SPT-SZ Survey

Authors: T. de Haan, B. A. Benson, L. E. Bleem, S. W. Allen, D. E. Applegate, M. L. N. Ashby, M. Bautz, M. Bayliss, S. Bocquet, M. Brodwin, J. E. Carlstrom, C. L. Chang, I. Chiu, H-M. Cho, A. Clocchiatti, T. M. Crawford, A. T. Crites, S. Desai, J. P. Dietrich, M. A. Dobbs, A. N. Doucouliagos, R. J. Foley, W. R. Forman, G. P. Garmire, E. M. George , et al. (52 additional authors not shown)

Abstract: (abridged) We present cosmological constraints obtained from galaxy clusters identified by their Sunyaev-Zel'dovich effect signature in the 2500 square degree South Pole Telescope Sunyaev Zel'dovich survey. We consider the 377 cluster candidates identified at z>0.25 with a detection significance greater than five, corresponding to the 95% purity threshold for the survey. We compute constraints on… ▽ More (abridged) We present cosmological constraints obtained from galaxy clusters identified by their Sunyaev-Zel'dovich effect signature in the 2500 square degree South Pole Telescope Sunyaev Zel'dovich survey. We consider the 377 cluster candidates identified at z>0.25 with a detection significance greater than five, corresponding to the 95% purity threshold for the survey. We compute constraints on cosmological models using the measured cluster abundance as a function of mass and redshift. We include additional constraints from multi-wavelength observations, including Chandra X-ray data for 82 clusters and a weak lensing-based prior on the normalization of the mass-observable scaling relations. Assuming a LCDM cosmology, where the species-summed neutrino mass has the minimum allowed value (mnu = 0.06 eV) from neutrino oscillation experiments, we combine the cluster data with a prior on H0 and find sigma_8 = 0.797+-0.031 and Omega_m = 0.289+-0.042, with the parameter combination sigma_8(Omega_m/0.27)^0.3 = 0.784+-0.039. These results are in good agreement with constraints from the CMB from SPT, WMAP, and Planck, as well as with constraints from other cluster datasets. Adding mnu as a free parameter, we find mnu = 0.14+-0.08 eV when combining the SPT cluster data with Planck CMB data and BAO data, consistent with the minimum allowed value. Finally, we consider a cosmology where mnu and N_eff are fixed to the LCDM values, but the dark energy equation of state parameter w is free. Using the SPT cluster data in combination with an H0 prior, we measure w = -1.28+-0.31, a constraint consistent with the LCDM cosmological model and derived from the combination of growth of structure and geometry. When combined with primarily geometrical constraints from Planck CMB, H0, BAO and SNe, adding the SPT cluster data improves the w constraint from the geometrical data alone by 14%, to w = -1.023+-0.042. △ Less

Submitted 21 March, 2016; originally announced March 2016.

Towards a Network of Faint DA White Dwarfs as High-Precision Spectrophotometric Standards

Authors: Gautham Narayan, Tim Axelrod, Jay B. Holberg, Thomas Matheson, Abhijit Saha, Edward W. Olszewski, Jenna Claver, Christopher W. Stubbs, Ralph C. Bohlin, Susana Deustua, Armin Rest

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

Submitted 11 March, 2016; originally announced March 2016.

Comments: (15 pages, 7 figures, 5 tables, accepted for publication in ApJ)

Light Curves of 213 Type Ia Supernovae from the ESSENCE Survey

Authors: Gautham Narayan, Armin Rest, Brad E. Tucker, Ryan J. Foley, W. Michael Wood-Vasey, Peter Challis, Christopher W. Stubbs, Robert P. Kirshner, Claudio Aguilera, Andrew C. Becker, Stephane Blondin, Alejandro Clocchiatti, Ricardo Covarrubias, Guillermo Damke, Tamara M. Davis, Alexei V. Filippenko, Mohan Ganeshalingam, Arti Garg, Peter M. Garnavich, Malcolm Hicken, Saurabh W. Jha, Kevin Krisciunas, Bruno Leibundgut, Weidong Li, Thomas Matheson , et al. (12 additional authors not shown)

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

Submitted 11 March, 2016; originally announced March 2016.

Comments: (40 pages, 21 figures, accepted for publication in ApJS)

Localization and broadband follow-up of the gravitational-wave transient GW150914

Authors: B. P. Abbott, R. Abbott, T. D. Abbott, M. R. Abernathy, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, A. Allocca, P. A. Altin, S. B. Anderson, W. G. Anderson, K. Arai , et al. (1522 additional authors not shown)

Abstract: A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared wit… ▽ More A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams. △ Less

Submitted 21 July, 2016; v1 submitted 26 February, 2016; originally announced February 2016.

Comments: For Supplement, see https://arxiv.org/abs/1604.07864

Report number: LIGO-P1500227-v12

Journal ref: The Astrophysical Journal Letters, 826:L13 (8pp), 2016 July 20

Pan-STARRS and PESSTO search for an optical counterpart to the LIGO gravitational wave source GW150914

Authors: S. J. Smartt, K. C. Chambers, K. W. Smith, M. E. Huber, D. R. Young, E. Cappellaro, D. E. Wright, M. Coughlin, A. S. B. Schultz, L. Denneau, H. Flewelling, A. Heinze, E. A. Magnier, N. Primak, A. Rest, A. Sherstyuk, B. Stalder, C. W. Stubbs, J. Tonry, C. Waters, M. Willman, J. P. Anderson, C. Baltay, M. T. Botticella, H. Campbell , et al. (26 additional authors not shown)

Abstract: We searched for an optical counterpart to the first gravitational wave source discovered by LIGO (GW150914), using a combination of the Pan-STARRS1 wide-field telescope and the PESSTO spectroscopic follow-up programme. As the final LIGO sky maps changed during analysis, the total probability of the source being spatially coincident with our fields was finally only 4.2 per cent. Therefore we discus… ▽ More We searched for an optical counterpart to the first gravitational wave source discovered by LIGO (GW150914), using a combination of the Pan-STARRS1 wide-field telescope and the PESSTO spectroscopic follow-up programme. As the final LIGO sky maps changed during analysis, the total probability of the source being spatially coincident with our fields was finally only 4.2 per cent. Therefore we discuss our results primarily as a demonstration of the survey capability of Pan-STARRS and spectroscopic capability of PESSTO. We mapped out 442 square degrees of the northern sky region of the initial map. We discovered 56 astrophysical transients over a period of 41 days from the discovery of the source. Of these, 19 were spectroscopically classified and a further 13 have host galaxy redshifts. All transients appear to be fairly normal supernovae and AGN variability and none is obviously linked with GW150914. We illustrate the sensitivity of our survey by defining parameterised lightcurves with timescales of 4, 20 and 40 days and use the sensitivity of the Pan-STARRS1 images to set limits on the luminosities of possible sources. The Pan-STARRS1 images reach limiting magnitudes of i = 19.2, 20.0 and 20.8 respectively for the three timescales. For long timescale parameterised lightcurves (with FWHM=~40d) we set upper limits of M_i <= -17.2 -0.9/+1.4 if the distance to GW150914 is D = 400 +/- 200Mpc. The number of type Ia SN we find in the survey is similar to that expected from the cosmic SN rate, indicating a reasonably complete efficiency in recovering supernova like transients out to D = 400 +/- 200 Mpc. △ Less

Submitted 28 July, 2016; v1 submitted 12 February, 2016; originally announced February 2016.

Comments: Accepted for publication on 28 July 2016 in MNRAS.This is the accepted version after referee report and editor acceptance

Precise Astronomical Flux Calibration and its Impact on Studying the Nature of Dark Energy

Authors: Christopher W. Stubbs, Yorke J. Brown

Abstract: Measurements of the luminosity of type Ia supernovae vs. redshift provided the original evidence for the accelerating expansion of the Universe and the existence of dark energy. Despite substantial improvements in survey methodology, systematic uncertainty in flux calibration dominates the error budget for this technique, exceeding both statistics and other systematic uncertainties. Consequently,… ▽ More Measurements of the luminosity of type Ia supernovae vs. redshift provided the original evidence for the accelerating expansion of the Universe and the existence of dark energy. Despite substantial improvements in survey methodology, systematic uncertainty in flux calibration dominates the error budget for this technique, exceeding both statistics and other systematic uncertainties. Consequently, any further collection of type Ia supernova data will fail to refine the constraints on the nature of dark energy unless we also improve the state of the art in astronomical flux calibration to the order of 1%. We describe how these systematic errors arise from calibration of instrumental sensitivity, atmospheric transmission, and Galactic extinction, and discuss ongoing efforts to meet the 1% precision challenge using white dwarf stars as celestial standards, exquisitely calibrated detectors as fundamental metrologic standards, and real-time atmospheric monitoring. △ Less

Submitted 15 January, 2016; originally announced January 2016.

Comments: 25 pages, 7 figures. Accepted Modern Physics Letters A

Journal ref: Modern Physics Letters A vol 30 number 40 1530030 (2015)

Hypercalibration: A Pan-STARRS1-based recalibration of the Sloan Digital Sky Survey

Authors: Douglas P. Finkbeiner, Edward F. Schlafly, David J. Schlegel, Nikhil Padmanabhan, Mario Juric, William S. Burgett, Kenneth C. Chambers, Larry Denneau, Peter W. Draper, Heather Flewelling, Klaus W. Hodapp, Nick Kaiser, E. A. Magnier, N. Metcalfe, Jeffrey S. Morgan, Paul A. Price, Christopher W. Stubbs, John L. Tonry

Abstract: We present a recalibration of the Sloan Digital Sky Survey (SDSS) photometry with new flat fields and zero points derived from Pan-STARRS1 (PS1). Using PSF photometry of 60 million stars with $16 < r < 20$, we derive a model of amplifier gain and flat-field corrections with per-run RMS residuals of 3 millimagnitudes (mmag) in $griz$ bands and 15 mmag in $u$ band. The new photometric zero points ar… ▽ More We present a recalibration of the Sloan Digital Sky Survey (SDSS) photometry with new flat fields and zero points derived from Pan-STARRS1 (PS1). Using PSF photometry of 60 million stars with $16 < r < 20$, we derive a model of amplifier gain and flat-field corrections with per-run RMS residuals of 3 millimagnitudes (mmag) in $griz$ bands and 15 mmag in $u$ band. The new photometric zero points are adjusted to leave the median in the Galactic North unchanged for compatibility with previous SDSS work. We also identify transient non-photometric periods in SDSS ("contrails") based on photometric deviations co-temporal in SDSS bands. The recalibrated stellar PSF photometry of SDSS and PS1 has an RMS difference of {9,7,7,8} mmag in $griz$, respectively, when averaged over $15'$ regions. △ Less

Submitted 3 December, 2015; originally announced December 2015.

Comments: 13 pages, 7 figures, ApJ, in press. "Hypercalibration" refers to using repeat measurements of many stars from multiple surveys to constrain calibration parameters

Detection of Enhancement in Number Densities of Background Galaxies due to Magnification by Massive Galaxy Clusters

Authors: I. Chiu, J. P. Dietrich, J. Mohr, D. E. Applegate, B. A. Benson, L. E. Bleem, M. B. Bayliss, S. Bocquet, J. E. Carlstrom, R. Capasso, S. Desai, C. Gangkofner, A. H. Gonzalez, N. Gupta, C. Hennig, H. Hoekstra, A. von der Linden, J. Liu, M. McDonald, C. L. Reichardt, A. Saro, T. Schrabback, V. Strazzullo, C. W. Stubbs, A. Zenteno

Abstract: We present a detection of the enhancement in the number densities of background galaxies induced from lensing magnification and use it to test the Sunyaev-Zel'dovich effect (SZE) inferred masses in a sample of 19 galaxy clusters with median redshift $z\simeq0.42$ selected from the South Pole Telescope SPT-SZ survey. Two background galaxy populations are selected for this study through their photom… ▽ More We present a detection of the enhancement in the number densities of background galaxies induced from lensing magnification and use it to test the Sunyaev-Zel'dovich effect (SZE) inferred masses in a sample of 19 galaxy clusters with median redshift $z\simeq0.42$ selected from the South Pole Telescope SPT-SZ survey. Two background galaxy populations are selected for this study through their photometric colours; they have median redshifts ${z}_{\mathrm{median}}\simeq0.9$ (low-$z$ background) and ${z}_{\mathrm{median}}\simeq1.8$ (high-$z$ background). Stacking these populations, we detect the magnification bias effect at $3.3σ$ and $1.3σ$ for the low- and high-$z$ backgrounds, respectively. We fit NFW models simultaneously to all observed magnification bias profiles to estimate the multiplicative factor $η$ that describes the ratio of the weak lensing mass to the mass inferred from the SZE observable-mass relation. We further quantify systematic uncertainties in $η$ resulting from the photometric noise and bias, the cluster galaxy contamination and the estimations of the background properties. The resulting $η$ for the combined background populations with $1σ$ uncertainties is $0.83\pm0.24\mathrm{(stat)}\pm0.074\mathrm{(sys)}$, indicating good consistency between the lensing and the SZE-inferred masses. We use our best-fit $η$ to predict the weak lensing shear profiles and compare these predictions with observations, showing agreement between the magnification and shear mass constraints. This work demonstrates the promise of using the magnification as a complementary method to estimate cluster masses in large surveys. △ Less

Submitted 9 February, 2016; v1 submitted 6 October, 2015; originally announced October 2015.

Comments: 16 pages, 10 figures, accepted for publication in MNRAS

Supercal: Cross-Calibration of Multiple Photometric Systems to Improve Cosmological Measurements with Type Ia Supernovae

Authors: D. Scolnic, S. Casertano, A. G. Riess, A. Rest, E. Schlafly, R. J. Foley, D. Finkbeiner, C. Tang, W. S. Burgett, K. C. Chambers, P. W. Draper, H. Flewelling, K. W. Hodapp, M. E. Huber, N. Kaiser, R. P. Kudritzki, E. A. Magnier, N. Metcalfe, C. W. Stubbs

Abstract: Current cosmological analyses which use Type Ia supernova (SN Ia) observations combine SN samples to expand the redshift range beyond that of a single sample and increase the overall sample size. The inhomogeneous photometric calibration between different SN samples is one of the largest systematic uncertainties of the cosmological parameter estimation. To place these different samples on a single… ▽ More Current cosmological analyses which use Type Ia supernova (SN Ia) observations combine SN samples to expand the redshift range beyond that of a single sample and increase the overall sample size. The inhomogeneous photometric calibration between different SN samples is one of the largest systematic uncertainties of the cosmological parameter estimation. To place these different samples on a single system, analyses currently use observations of a small sample of very bright flux standards on the $HST$ system. We propose a complementary method, called `Supercal', in which we use measurements of secondary standards in each system, compare these to measurements of the same stars in the Pan-STARRS1 (PS1) system, and determine offsets for each system relative to PS1, placing all SN observations on a single, consistent photometric system. PS1 has observed $3π$ of the sky and has a relative calibration of better than 5 mmag (for $\sim15<griz<21$ mag), making it an ideal reference system. We use this process to recalibrate optical observations taken by the following SN samples: PS1, SNLS, SDSS, CSP, and CfA1-4. We measure discrepancies on average of 10 mmag, but up to 35 mmag, in various optical passbands. We find that correcting for these differences changes recovered values for the dark energy equation-of-state parameter, $w$, by on average $2.6\%$. This change is roughly half the size of current statistical constraints on $w$. The size of this effect strongly depends on the error in the $B-V$ calibration of the low-$z$ surveys. The Supercal method will allow future analyses to tie past samples to the best calibrated sample. △ Less

Submitted 11 November, 2015; v1 submitted 17 August, 2015; originally announced August 2015.

Comments: 15 pages, 5 figures. Accepted by ApJ

GALEX Detection of Shock Breakout in Type II-P Supernova PS1-13arp: Implications for the Progenitor Star Wind

Authors: S. Gezari, D. O. Jones, N. E. Sanders, A. M. Soderberg, T. Hung, S. Heinis, S. J. Smartt, A. Rest, D. Scolnic, R. Chornock, E. Berger, R. J. Foley, M. E. Huber, P. Price C. W. Stubbs, A. G. Riess, R. P. Kirshner, K. Smith, W. M. Wood-Vasey, D. Schiminovich, D. C. Martin, W. S. Burgett, K. C. Chambers, H. Flewelling, N. Kaiser, J. L. Tonry , et al. (1 additional authors not shown)

Abstract: We present the GALEX detection of a UV burst at the time of explosion of an optically normal Type II-P supernova (PS1-13arp) from the Pan-STARRS1 survey at z=0.1665. The temperature and luminosity of the UV burst match the theoretical predictions for shock breakout in a red supergiant, but with a duration a factor of ~50 longer than expected. We compare the $NUV$ light curve of PS1-13arp to previo… ▽ More We present the GALEX detection of a UV burst at the time of explosion of an optically normal Type II-P supernova (PS1-13arp) from the Pan-STARRS1 survey at z=0.1665. The temperature and luminosity of the UV burst match the theoretical predictions for shock breakout in a red supergiant, but with a duration a factor of ~50 longer than expected. We compare the $NUV$ light curve of PS1-13arp to previous GALEX detections of Type IIP SNe, and find clear distinctions that indicate that the UV emission is powered by shock breakout, and not by the subsequent cooling envelope emission previously detected in these systems. We interpret the ~ 1 d duration of the UV signal with a shock breakout in the wind of a red supergiant with a pre-explosion mass-loss rate of ~ 10^-3 Msun yr^-1. This mass-loss rate is enough to prolong the duration of the shock breakout signal, but not enough to produce an excess in the optical plateau light curve or narrow emission lines powered by circumstellar interaction. This detection of non-standard, potentially episodic high mass-loss in a RSG SN progenitor has favorable consequences for the prospects of future wide-field UV surveys to detect shock breakout directly in these systems, and provide a sensitive probe of the pre-explosion conditions of SN progenitors. △ Less

Submitted 24 February, 2015; originally announced February 2015.

Comments: Accepted for Publication in ApJ

Selection of Burst-like Transients and Stochastic Variables Using Multi-Band Image Differencing in the Pan-STARRS1 Medium-Deep Survey

Authors: S. Kumar, S. Gezari, S. Heinis, R. Chornock, E. Berger, A. Rest, M. E. Huber, R. J. Foley, G. Narayan, G. H. Marion, D. Scolnic, A. Soderberg, A. Lawrence, C. W. Stubbs, R. P. Kirshner, A. G. Riess, S. J. Smartt, K. Smith, W. M Wood-Vasey, W. S. Burgett, K. C. Chambers, H. Flewelling, N. Kaiser, N. Metcalfe, P. A. Price , et al. (2 additional authors not shown)

Abstract: We present a novel method for the light-curve characterization of Pan-STARRS1 Medium Deep Survey (PS1 MDS) extragalactic sources into stochastic variables (SV) and burst-like (BL) transients, using multi-band image-differencing time-series data. We select detections in difference images associated with galaxy hosts using a star/galaxy catalog extracted from the deep PS1 MDS stacked images, and ado… ▽ More We present a novel method for the light-curve characterization of Pan-STARRS1 Medium Deep Survey (PS1 MDS) extragalactic sources into stochastic variables (SV) and burst-like (BL) transients, using multi-band image-differencing time-series data. We select detections in difference images associated with galaxy hosts using a star/galaxy catalog extracted from the deep PS1 MDS stacked images, and adopt a maximum a posteriori formulation to model their difference-flux time-series in four Pan-STARRS1 photometric bands g,r,i, and z. We use three deterministic light-curve models to fit burst-like transients and one stochastic light curve model, the Ornstein-Uhlenbeck process, in order to fit variability that is characteristic of active galactic nuclei (AGN). We assess the quality of fit of the models band-wise source-wise, using their estimated leave-out-one cross-validation likelihoods and corrected Akaike information criteria. We then apply a K-means clustering algorithm on these statistics, to determine the source classification in each band. The final source classification is derived as a combination of the individual filter classifications. We use our clustering method to characterize 4361 extragalactic image difference detected sources in the first 2.5 years of the PS1 MDS, into 1529 BL, and 2262 SV, with a purity of 95.00% for AGN, and 90.97% for SN based on our verification sets. We combine our light-curve classifications with their nuclear or off-nuclear host galaxy offsets, to define a robust photometric sample of 1233 active galactic nuclei and 812 supernovae. We use these samples to identify simple photometric priors that would enable their real-time identification in future wide-field synoptic surveys. △ Less

Submitted 6 January, 2015; originally announced January 2015.

Comments: Accepted for publication in ApJ