Showing 1–50 of 114 results for author: Larson, S

LISA Definition Study Report

Authors: Monica Colpi, Karsten Danzmann, Martin Hewitson, Kelly Holley-Bockelmann, Philippe Jetzer, Gijs Nelemans, Antoine Petiteau, David Shoemaker, Carlos Sopuerta, Robin Stebbins, Nial Tanvir, Henry Ward, William Joseph Weber, Ira Thorpe, Anna Daurskikh, Atul Deep, Ignacio Fernández Núñez, César García Marirrodriga, Martin Gehler, Jean-Philippe Halain, Oliver Jennrich, Uwe Lammers, Jonan Larrañaga, Maike Lieser, Nora Lützgendorf , et al. (86 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) is the first scientific endeavour to detect and study gravitational waves from space. LISA will survey the sky for Gravitational Waves in the 0.1 mHz to 1 Hz frequency band which will enable the study of a vast number of objects ranging from Galactic binaries and stellar mass black holes in the Milky Way, to distant massive black-hole mergers and the e… ▽ More The Laser Interferometer Space Antenna (LISA) is the first scientific endeavour to detect and study gravitational waves from space. LISA will survey the sky for Gravitational Waves in the 0.1 mHz to 1 Hz frequency band which will enable the study of a vast number of objects ranging from Galactic binaries and stellar mass black holes in the Milky Way, to distant massive black-hole mergers and the expansion of the Universe. This definition study report, or Red Book, presents a summary of the very large body of work that has been undertaken on the LISA mission over the LISA definition phase. △ Less

Submitted 12 February, 2024; originally announced February 2024.

Comments: 155 pages, with executive summary and table of contents

Understanding Binary Systems -- a Comparison between COSMIC and MESA

Authors: Lailani Kenoly, A. Ken Luu, Celia Toral, Roseanne M. Cheng, Nicole M. Lloyd-Ronning, Shane L. Larson, Gabriel O. Casabona

Abstract: We compare the evolution of binary systems evolved in the MESA stellar evolution code to those in the COSMIC population synthesis code. Our aim is to convey the robustness of the equations that model binary evolution in the COSMIC code, particularly for the cases of high mass stars with closely orbiting compact object companions. Our larger goal is to accurately model the rates of these systems, a… ▽ More We compare the evolution of binary systems evolved in the MESA stellar evolution code to those in the COSMIC population synthesis code. Our aim is to convey the robustness of the equations that model binary evolution in the COSMIC code, particularly for the cases of high mass stars with closely orbiting compact object companions. Our larger goal is to accurately model the rates of these systems, as they are promising candidates for the progenitor systems behind energetic, longer lasting, radio bright GRB jets. These systems also may be key contributors to the rates of binary black hole mergers throughout our universe. △ Less

Submitted 17 May, 2023; originally announced May 2023.

Comments: to appear in RNAAS

Report number: LA-UR-23-24423

Sublimation Origin of Active Asteroid P/2018 P3

Authors: Yoonyoung Kim, Jessica Agarwal, David Jewitt, Max Mutchler, Stephen Larson, Harold Weaver, Michael Mommert

Abstract: Active asteroids show (typically transient) cometary activity, driven by a range of processes. A sub-set, sometimes called main-belt comets, may be driven by sublimation and so could be useful for tracing the present-day distribution of asteroid ice. Object P/2018 P3 has a Tisserand parameter 3.096 but a high eccentricity 0.415, placing it within the dynamical boundary between asteroids and comets… ▽ More Active asteroids show (typically transient) cometary activity, driven by a range of processes. A sub-set, sometimes called main-belt comets, may be driven by sublimation and so could be useful for tracing the present-day distribution of asteroid ice. Object P/2018 P3 has a Tisserand parameter 3.096 but a high eccentricity 0.415, placing it within the dynamical boundary between asteroids and comets. We aim to determine the cause of activity (sublimation or something else) and assess the dynamical stability of P3, in order to better constrain the intrinsic ice content in the main belt. We obtained Hubble Space Telescope images of P3 at the highest angular resolution. We compared the observations with a Monte Carlo model of dust dynamics. We identified and analyzed archival CFHT (2013) and NEOWISE (2018) data. In addition, we numerically integrated the orbits of P3 clones for 100 Myr. P3 has been recurrently active near two successive perihelia (at 1.76 AU), indicative of a sublimation origin. The absence of 4.6 um band excess indicates zero or negligible CO or CO2 gas production from P3. The properties of the ejected dust are remarkably consistent with those found in other main-belt comets (continuous emission of ~0.05-5 mm particles at 0.3-3 m/s speeds), with mass-loss rates of >~2 kg/s. The orbit of P3 is unstable on timescales ~10 Myr. We speculate that P3 has recently arrived from a more stable source (either the Kuiper Belt or elsewhere in the main belt) and has been physically aged at its current location, finally becoming indistinguishable from a weakly sublimating asteroid in terms of its dust properties. Whatever the source of P3, given the dynamical instability of its current orbit, P3 should not be used to trace the native distribution of asteroid ice. △ Less

Submitted 16 August, 2022; originally announced August 2022.

Comments: 8 pages, 12 figures, A&A in press

Journal ref: A&A 666, A163 (2022)

Optical trapping of high-aspect-ratio NaYF hexagonal prisms for kHz-MHz gravitational wave detectors

Authors: George Winstone, Zhiyuan Wang, Shelby Klomp, Greg Felsted, Andrew Laeuger, Daniel Grass, Nancy Aggarwal, Jacob Sprague, Peter J. Pauzauskie, Shane L. Larson, Vicky Kalogera, Andrew A. Geraci

Abstract: We present experimental results on optical trapping of Yb-doped $β-$NaYF sub-wavelength-thickness high-aspect-ratio hexagonal prisms with a micron-scale radius. The prisms are trapped in vacuum using an optical standing wave, oriented with the normal vector to their face along the beam propagation direction, and exhibit characteristic modes of three translational and two torsional degrees of freed… ▽ More We present experimental results on optical trapping of Yb-doped $β-$NaYF sub-wavelength-thickness high-aspect-ratio hexagonal prisms with a micron-scale radius. The prisms are trapped in vacuum using an optical standing wave, oriented with the normal vector to their face along the beam propagation direction, and exhibit characteristic modes of three translational and two torsional degrees of freedom. The measured motional spectra are compared with numerical simulations. This plate-like geometry simultaneously enables trapping with low photon-recoil-heating, high mass, and high trap frequency, potentially leading to advances in high frequency gravitational wave searches in the Levitated Sensor Detector (LSD), currently under construction. The material used here has previously been shown to exhibit internal cooling via laser refrigeration when optically trapped and illuminated with light of suitable wavelength. Employing such laser refrigeration methods in the context of our work may enable higher trapping intensity thus and higher trap frequencies for gravitational wave searches approaching the several hundred kHz range. △ Less

Submitted 25 April, 2022; v1 submitted 22 April, 2022; originally announced April 2022.

Comments: Title and Fig. 1 updated, minor changes, prepared for journal submission

Astrophysics with the Laser Interferometer Space Antenna

Authors: Pau Amaro Seoane, Jeff Andrews, Manuel Arca Sedda, Abbas Askar, Quentin Baghi, Razvan Balasov, Imre Bartos, Simone S. Bavera, Jillian Bellovary, Christopher P. L. Berry, Emanuele Berti, Stefano Bianchi, Laura Blecha, Stephane Blondin, Tamara Bogdanović, Samuel Boissier, Matteo Bonetti, Silvia Bonoli, Elisa Bortolas, Katelyn Breivik, Pedro R. Capelo, Laurentiu Caramete, Federico Cattorini, Maria Charisi, Sylvain Chaty , et al. (134 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery… ▽ More The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe. △ Less

Submitted 25 May, 2023; v1 submitted 11 March, 2022; originally announced March 2022.

Journal ref: Living Reviews in Relativity, Volume 26, Article number: 2 (2023)

Deep Co-Added Sky from Catalina Sky Survey Images

Authors: Akshat Singhal, Varun Bhalerao, Ashish A. Mahabal, Kaustubh Vaghmare, Santosh Jagade, Sumeet Kulkarni, Ajay Vibhute, Ajit K. Kembhavi, Andrew J. Drake, S George Djorgovski, Matthew J. Graham, Ciro Donalek, Eric Christensen, Stephen Larson, Edward C. Beshore

Abstract: A number of synoptic sky surveys are underway or being planned. Typically they are done with small telescopes and relatively short exposure times. A search for transient or variable sources involves comparison with deeper baseline images, ideally obtained through the same telescope and camera. With that in mind we have stacked images from the 0.68~m Schmidt telescope on Mt. Bigelow taken over ten… ▽ More A number of synoptic sky surveys are underway or being planned. Typically they are done with small telescopes and relatively short exposure times. A search for transient or variable sources involves comparison with deeper baseline images, ideally obtained through the same telescope and camera. With that in mind we have stacked images from the 0.68~m Schmidt telescope on Mt. Bigelow taken over ten years as part of the Catalina Sky Survey. In order to generate deep reference images for the Catalina Real-time Transient Survey, close to 0.8 million images over 8000 fields and covering over 27000~sq.~deg. have gone into the deep stack that goes up to 3 magnitudes deeper than individual images. CRTS system does not use a filter in imaging, hence there is no standard passband in which the optical magnitude is measured. We estimate depth by comparing these wide-band unfiltered co-added images with images in the $g$-band and find that the image depth ranges from 22.0--24.2 across the sky, with a 200-image stack attaining an equivalent AB magnitude sensitivity of 22.8. We compared various state-of-the-art software packages for co-adding astronomical images and have used SWarp for the stacking. We describe here the details of the process adopted. This methodology may be useful in other panoramic imaging applications, and to other surveys as well. The stacked images are available through a server at Inter-University Centre for Astronomy and Astrophysics (IUCAA). △ Less

Submitted 30 July, 2021; originally announced August 2021.

Comments: 16 pages, 15 figures

Getting Ready for LISA: The Data, Support and Preparation Needed to Maximize US Participation in Space-Based Gravitational Wave Science

Authors: Kelly Holley-Bockelmann, :, Jillian Bellovary, Peter Bender, Emanuele Berti, Warren Brown, Robert Caldwell, Neil Cornish, Jeremy Darling, Matthew Digman, Mike Eracleous, Kayhan Gultekin, Zoltan Haiman, Kelly Holley-Bockelmann, Joey Key, Shane Larson, Xin Liu, Sean McWilliams, Priyamvada Natarajan, David Shoemaker, Deirdre Shoemaker, Krista Lynne Smith, Marcelle Soares-Santos, Robin Stebbins

Abstract: The NASA LISA Study Team was tasked to study how NASA might support US scientists to participate and maximize the science return from the Laser Interferometer Space Antenna (LISA) mission. LISA is gravitational wave observatory led by ESA with NASA as a junior partner, and is scheduled to launch in 2034. Among our findings: LISA science productivity is greatly enhanced by a full-featured US scienc… ▽ More The NASA LISA Study Team was tasked to study how NASA might support US scientists to participate and maximize the science return from the Laser Interferometer Space Antenna (LISA) mission. LISA is gravitational wave observatory led by ESA with NASA as a junior partner, and is scheduled to launch in 2034. Among our findings: LISA science productivity is greatly enhanced by a full-featured US science center and an open access data model. As other major missions have demonstrated, a science center acts as both a locus and an amplifier of research innovation, data analysis, user support, user training and user interaction. In its most basic function, a US Science Center could facilitate entry into LISA science by hosting a Data Processing Center and a portal for the US community to access LISA data products. However, an enhanced LISA Science Center could: support one of the parallel independent processing pipelines required for data product validation; stimulate the high level of research on data analysis that LISA demands; support users unfamiliar with a novel observatory; facilitate astrophysics and fundamental research; provide an interface into the subtleties of the instrument to validate extraordinary discoveries; train new users; and expand the research community through guest investigator, postdoc and student programs. Establishing a US LISA Science Center well before launch can have a beneficial impact on the participation of the broader astronomical community by providing training, hosting topical workshops, disseminating mock catalogs, software pipelines, and documentation. Past experience indicates that successful science centers are established several years before launch; this early adoption model may be especially relevant for a pioneering mission like LISA. △ Less

Submitted 4 December, 2020; originally announced December 2020.

Comments: 93 pages with a lovely cover page thanks to Bernard Kelly and Elizabeth Ferrara

Component properties and mutual orbit of binary main-belt comet 288P/(300163) 2006 VW139

Authors: Jessica Agarwal, Yoonyoung Kim, David Jewitt, Max Mutchler, Harold Weaver, Stephen Larson

Abstract: The binary asteroid 288P/(300163) is unusual both for its combination of wide-separation and high mass ratio and for its comet-like activity. It is not currently known whether there is a causal connection between the activity and the unusual orbit or if instead the activity helped to overcome a strong detection bias against such sub-arcsecond systems. We aim to find observational constraints discr… ▽ More The binary asteroid 288P/(300163) is unusual both for its combination of wide-separation and high mass ratio and for its comet-like activity. It is not currently known whether there is a causal connection between the activity and the unusual orbit or if instead the activity helped to overcome a strong detection bias against such sub-arcsecond systems. We aim to find observational constraints discriminating between possible formation scenarios and to characterise the physical properties of the system components. We measured the component separation and brightness using point spread function fitting to high-resolution Hubble Space Telescope/Wide Field Camera 3 images from 25 epochs between 2011 and 2020. We constrained component sizes and shapes from the photometry, and we fitted a Keplerian orbit to the separation as a function of time. Approximating the components A and B as prolate spheroids with semi-axis lengths a$<$b and assuming a geometric albedo of 0.07, we find $a_A \leq$ 0.6 km, $b_A \geq$ 1.4 km, $a_B \leq$ 0.5 km, and $b_B \geq$ 0.8 km. We find indications that the dust production may have concentrated around B and that the mutual orbital period may have changed by 1-2 days during the 2016 perihelion passage. Orbit solutions have semi-major axes in the range of (105-109) km, eccentricities between 0.41 and 0.51, and periods of (117.3-117.5) days pre-perihelion and (118.5-119.5) days post-perihelion, corresponding to system masses in the range of (6.67-7.23) $\times$ 10$^{12}$ kg. The mutual and heliocentric orbit planes are roughly aligned. Based on the orbit alignment, we infer that spin-up of the precursor by the YORP effect led to the formation of the binary system. We disfavour (but cannot exclude) a scenario of very recent formation where activity was directly triggered by the break-up, because our data support a scenario with a single active component. △ Less

Submitted 29 September, 2020; originally announced September 2020.

Comments: 39 pages, 7 figures, 2 tables, 3 appendices with 28 additional figures, accepted for publication in A&A

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

GPU-Accelerated Periodic Source Identification in Large-Scale Surveys: Measuring $P$ and $\dot{P}$

Authors: Michael L. Katz, Olivia R. Cooper, Michael W. Coughlin, Kevin B. Burdge, Katelyn Breivik, Shane L. Larson

Abstract: Many inspiraling and merging stellar remnants emit both gravitational and electromagnetic radiation as they orbit or collide. These gravitational wave events together with their associated electromagnetic counterparts provide insight about the nature of the merger, allowing us to further constrain properties of the binary. With the future launch of the Laser Interferometer Space Antenna (LISA), fo… ▽ More Many inspiraling and merging stellar remnants emit both gravitational and electromagnetic radiation as they orbit or collide. These gravitational wave events together with their associated electromagnetic counterparts provide insight about the nature of the merger, allowing us to further constrain properties of the binary. With the future launch of the Laser Interferometer Space Antenna (LISA), follow up observations and models are needed of ultracompact binary (UCB) systems. Current and upcoming long baseline time domain surveys will observe many of these UCBs. We present a new fast periodic object search tool capable of searching for generic periodic signals based on the Conditional Entropy algorithm. This new implementation allows for a grid search over both the period ($P$) and the time derivative of the period ($\dot{P}$). To demonstrate the usage of this tool, we use a small, hand-picked subset of a UCB population generated from the population synthesis code \cosmic, as well as a custom catalog for varying periods at fixed intrinsic parameters. We simulate light curves as likely to be observed by future time domain surveys by using an existing eclipsing binary light curve model accounting for the change in orbital period due to gravitational radiation. We find that a search with $\dot{P}$ values is necessary for detecting binaries at orbital periods less than $\sim$10 min. We also show it is useful in finding and characterizing binaries with longer periods, but at a higher computational cost. Our code is called gce (GPU-Accelerated Conditional Entropy). It is available at https://github.com/mikekatz04/gce. △ Less

Submitted 17 February, 2021; v1 submitted 11 June, 2020; originally announced June 2020.

Comments: 15 pages, 7 figures

Spurious Acceleration Noise on the LISA Spacecraft due to Solar Activity

Authors: Barrett M. Frank, Brandon Piotrzkowski, Brett Bolen, Marco Cavaglià, Shane L. Larson

Abstract: One source of noise for the Laser Interferometer Space Antenna (LISA) will be time-varying changes of the space environment in the form of solar wind particles and photon pressure from fluctuating solar irradiance. The approximate magnitude of these effects can be estimated from the average properties of the solar wind and the solar irradiance. We use data taken by the ACE (Advanced Compton Explor… ▽ More One source of noise for the Laser Interferometer Space Antenna (LISA) will be time-varying changes of the space environment in the form of solar wind particles and photon pressure from fluctuating solar irradiance. The approximate magnitude of these effects can be estimated from the average properties of the solar wind and the solar irradiance. We use data taken by the ACE (Advanced Compton Explorer) satellite and the VIRGO (Variability of solar IRradiance and Gravity Oscillations) instrument on the SOHO satellite over an entire solar cycle to calculate the forces due to solar wind and photon pressure irradiance on the LISA spacecraft. We produce a realistic model of the effects of these environmental noise sources and their variation over the expected course of the LISA mission. △ Less

Submitted 3 July, 2020; v1 submitted 16 December, 2019; originally announced December 2019.

Comments: 10 pages, 7 figures

Building A Field: The Future of Astronomy with Gravitational Waves, A State of The Profession Consideration for Astro2020

Authors: Kelly Holley-Bockelmann, Joey Shapiro Key, Brittany Kamai, Robert Caldwell, Warren Brown, Bill Gabella, Karan Jani, Quentin Baghi, John Baker, Jillian Bellovary, Pete Bender, Emanuele Berti, T. J. Brandt, Curt Cutler, John W. Conklin, Michael Eracleous, Elizabeth C. Ferrara, Bernard J. Kelly, Shane L. Larson, Jeff Livas, Maura McLaughlin, Sean T. McWilliams, Guido Mueller, Priyamvada Natarajan, Norman Rioux , et al. (6 additional authors not shown)

Abstract: Harnessing the sheer discovery potential of gravitational wave astronomy will require bold, deliberate, and sustained efforts to train and develop the requisite workforce. The next decade requires a strategic plan to build -- from the ground up -- a robust, open, and well-connected gravitational wave astronomy community with deep participation from traditional astronomers, physicists, data scienti… ▽ More Harnessing the sheer discovery potential of gravitational wave astronomy will require bold, deliberate, and sustained efforts to train and develop the requisite workforce. The next decade requires a strategic plan to build -- from the ground up -- a robust, open, and well-connected gravitational wave astronomy community with deep participation from traditional astronomers, physicists, data scientists, and instrumentalists. This basic infrastructure is sorely needed as an enabling foundation for research. We outline a set of recommendations for funding agencies, universities, and professional societies to help build a thriving, diverse, and inclusive new field. △ Less

Submitted 16 December, 2019; originally announced December 2019.

Eclipses of continuous gravitational waves as a probe of stellar structure

Authors: Pablo Marchant, Katelyn Breivik, Christopher P. L. Berry, Ilya Mandel, Shane L. Larson

Abstract: Although gravitational waves only interact weakly with matter, their propagation is affected by a gravitational potential. If a gravitational wave source is eclipsed by a star, measuring these perturbations provides a way to directly measure the distribution of mass throughout the stellar interior. We compute the expected Shapiro time delay, amplification and deflection during an eclipse, and show… ▽ More Although gravitational waves only interact weakly with matter, their propagation is affected by a gravitational potential. If a gravitational wave source is eclipsed by a star, measuring these perturbations provides a way to directly measure the distribution of mass throughout the stellar interior. We compute the expected Shapiro time delay, amplification and deflection during an eclipse, and show how this can be used to infer the mass distribution of the eclipsing body. We identify continuous gravitational waves from neutron stars as the best candidates to detect this effect. When the Sun eclipses a far-away source, depending on the depth of the eclipse the time-delay can change by up to $\sim 0.034$ ms, the gravitational-wave strain amplitude can increase by $\sim 4$%, and the apparent position of the source in the sky can vary by $4''$. Accreting neutron stars with Roche-lobe filling companion stars have a high probability of exhibiting eclipses, producing similar time delays but undetectable changes in amplitude and sky location. Even for the most rapidly rotating neutron stars, this time delay only corresponds to a few percent of the phase of the gravitational wave, making it an extremely challenging measurement. However, if sources of continuous gravitational waves exist just below the limit of detection of current observatories, next-generation instruments will be able to observe them with enough precision to measure the signal of an eclipsing star. Detecting this effect would provide a new direct probe to the interior of stars, complementing asteroseismology and the detection of solar neutrinos. △ Less

Submitted 9 December, 2019; originally announced December 2019.

Comments: Accepted for publication in Physical Review D

Constraining Galactic Structure with the LISA White Dwarf Foreground

Authors: Katelyn Breivik, Chiara M. F. Mingarelli, Shane L. Larson

Abstract: White dwarfs comprise 95% of all stellar remnants, and are thus an excellent tracer of old stellar populations in the Milky Way. Current and planned telescopes are not able to directly probe the white dwarf population in its entirety due to its inherently low luminosity. However, the Galactic population of double white dwarf binaries gives rise to a millihertz gravitational-wave foreground detecta… ▽ More White dwarfs comprise 95% of all stellar remnants, and are thus an excellent tracer of old stellar populations in the Milky Way. Current and planned telescopes are not able to directly probe the white dwarf population in its entirety due to its inherently low luminosity. However, the Galactic population of double white dwarf binaries gives rise to a millihertz gravitational-wave foreground detectable by the Laser Interferometer Space Antenna (LISA). Here we show how characterizing this foreground's angular power spectrum will enable us to probe the Galactic structure in a novel way and measure the vertical scale height of the Galaxy's oldest stellar populations. We do this using a binary population synthesis study that incorporates different Galactic spatial distributions for the double white dwarf population. We find that the level of anisotropy in the white dwarf foreground's angular power spectrum is strongly dependent on the vertical scale height of the population. Finally, we show that LISA can probe the vertical scale height of the Galactic double white dwarf population with an accuracy of 50 pc-200pc, depending on angular resolution limits, using the angular power spectrum of the white dwarf foreground. △ Less

Submitted 10 September, 2020; v1 submitted 4 December, 2019; originally announced December 2019.

Comments: 10 pages, 7 figures; reflects manuscript accepted in ApJ

COSMIC Variance in Binary Population Synthesis

Authors: Katelyn Breivik, Scott Coughlin, Michael Zevin, Carl L. Rodriguez, Kyle Kremer, Claire S. Ye, Jeff J. Andrews, Michael Kurkowski, Matthew C. Digman, Shane L. Larson, Frederic A. Rasio

Abstract: The formation and evolution of binary stars is a critical component of several fields in astronomy. The most numerous sources for gravitational wave observatories are inspiraling and/or merging compact binaries, while binary stars are present in nearly every electromagnetic survey regardless of the target population. Simulations of large binary populations serve to both predict and inform observat… ▽ More The formation and evolution of binary stars is a critical component of several fields in astronomy. The most numerous sources for gravitational wave observatories are inspiraling and/or merging compact binaries, while binary stars are present in nearly every electromagnetic survey regardless of the target population. Simulations of large binary populations serve to both predict and inform observations of electromagnetic and gravitational wave sources. Binary population synthesis is a tool that balances physical modeling with simulation speed to produce large binary populations on timescales of days. We present a community-developed binary population synthesis suite: COSMIC which is designed to simulate compact-object binary populations and their progenitors. As a proof of concept, we simulate the Galactic population of compact binaries and their gravitational wave signal observable by the Laser Interferometer Space Antenna (LISA). We find that $\sim10^8$ compact binaries reside in the Milky Way today, while $\sim10^4$ of them may be resolvable by LISA. △ Less

Submitted 10 September, 2020; v1 submitted 3 November, 2019; originally announced November 2019.

Comments: reflects the version submitted to ApJ; 17 pages, 2 Tables, 4 Figures; corrects typo in SNR calculation

A guide to LIGO-Virgo detector noise and extraction of transient gravitational-wave signals

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, C. Adams, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Alford, G. Allen, A. Allocca, M. A. Aloy, P. A. Altin, A. Amato, A. Ananyeva , et al. (1113 additional authors not shown)

Abstract: The LIGO Scientific Collaboration and the Virgo Collaboration have cataloged eleven confidently detected gravitational-wave events during the first two observing runs of the advanced detector era. All eleven events were consistent with being from well-modeled mergers between compact stellar-mass objects: black holes or neutron stars. The data around the time of each of these events have been made… ▽ More The LIGO Scientific Collaboration and the Virgo Collaboration have cataloged eleven confidently detected gravitational-wave events during the first two observing runs of the advanced detector era. All eleven events were consistent with being from well-modeled mergers between compact stellar-mass objects: black holes or neutron stars. The data around the time of each of these events have been made publicly available through the gravitational-wave open science center. The entirety of the gravitational-wave strain data from the first and second observing runs have also now been made publicly available. There is considerable interest among the broad scientific community in understanding the data and methods used in the analyses. In this paper, we provide an overview of the detector noise properties and the data analysis techniques used to detect gravitational-wave signals and infer the source properties. We describe some of the checks that are performed to validate the analyses and results from the observations of gravitational-wave events. We also address concerns that have been raised about various properties of LIGO-Virgo detector noise and the correctness of our analyses as applied to the resulting data. △ Less

Submitted 10 February, 2020; v1 submitted 29 August, 2019; originally announced August 2019.

Journal ref: B P Abbott et al 2020 Class. Quantum Grav. 37 055002

Probing Massive Black Hole Binary Populations with LISA

Authors: Michael L. Katz, Luke Zoltan Kelley, Fani Dosopoulou, Samantha Berry, Laura Blecha, Shane L. Larson

Abstract: ESA and NASA are moving forward with plans to launch LISA around 2034. With data from the Illustris cosmological simulation, we provide analysis of LISA detection rates accompanied by characterization of the merging massive black hole population. Massive black holes of total mass $\sim10^5-10^{10} M_\odot$ are the focus of this study. We evolve Illustris massive black hole mergers, which form at s… ▽ More ESA and NASA are moving forward with plans to launch LISA around 2034. With data from the Illustris cosmological simulation, we provide analysis of LISA detection rates accompanied by characterization of the merging massive black hole population. Massive black holes of total mass $\sim10^5-10^{10} M_\odot$ are the focus of this study. We evolve Illustris massive black hole mergers, which form at separations on the order of the simulation resolution ($\sim$kpc scales), through coalescence with two different treatments for the binary massive black hole evolutionary process. The coalescence times of the population, as well as physical properties of the black holes, form a statistical basis for each evolutionary treatment. From these bases, we Monte Carlo synthesize many realizations of the merging massive black hole population to build mock LISA detection catalogs. We analyze how our massive black hole binary evolutionary models affect detection rates and the associated parameter distributions measured by LISA. With our models, we find massive black hole binary detection rates with LISA of $\sim0.5-1$ yr$^{-1}$ for massive black holes with masses greater than $10^5M_\odot$. This should be treated as a lower limit primarily because our massive black hole sample does not include masses below $10^5M_\odot$, which may significantly add to the observed rate. We suggest reasons why we predict lower detection rates compared to much of the literature. △ Less

Submitted 1 November, 2019; v1 submitted 15 August, 2019; originally announced August 2019.

Comments: 18 pages, 9 figures, 3 tables

Space Based Gravitational Wave Astronomy Beyond LISA

Authors: John Baker, Simon F. Barke, Peter L. Bender, Emanuele Berti, Robert Caldwell, John W. Conklin, Neil Cornish, Elizabeth C. Ferrara, Kelly Holley-Bockelmann, Brittany Kamai, Shane L. Larson, Jeff Livas, Sean T. McWilliams, Guido Mueller, Priyamvada Natarajan, Norman Rioux, Shannon R Sankar, Jeremy Schnittman, Deirdre Shoemaker, Jacob Slutsky, Robin Stebbins, Ira Thorpe, John Ziemer

Abstract: The Laser Interferometer Space Antenna (LISA) will open three decades of gravitational wave (GW) spectrum between 0.1 and 100 mHz, the mHz band. This band is expected to be the richest part of the GW spectrum, in types of sources, numbers of sources, signal-to-noise ratios and discovery potential. When LISA opens the low-frequency window of the gravitational wave spectrum, around 2034, the surge o… ▽ More The Laser Interferometer Space Antenna (LISA) will open three decades of gravitational wave (GW) spectrum between 0.1 and 100 mHz, the mHz band. This band is expected to be the richest part of the GW spectrum, in types of sources, numbers of sources, signal-to-noise ratios and discovery potential. When LISA opens the low-frequency window of the gravitational wave spectrum, around 2034, the surge of gravitational-wave astronomy will strongly compel a subsequent mission to further explore the frequency bands of the GW spectrum that can only be accessed from space. The 2020s is the time to start developing technology and studying mission concepts for a large-scale mission to be launched in the 2040s. The mission concept would then be proposed to Astro2030. Only space based missions can access the GW spectrum between 10 nHz and 1 Hz because of the Earths seismic noise. This white paper surveys the science in this band and mission concepts that could accomplish that science. The proposed small scale activity is a technology development program that would support a range of concepts and a mission concept study to choose a specific mission concept for Astro2030. In this white paper, we will refer to a generic GW mission beyond LISA as bLISA. △ Less

Submitted 25 July, 2019; originally announced July 2019.

Comments: Astro2020 APC White Paper Primary Thematic Science Area: Multi-Messenger Astronomy and Astrophysics Secondary Areas: Cosmology and Fundamental Physics, Galaxy Evolution, Formation and Evolution of Compact Objects

The Laser Interferometer Space Antenna: Unveiling the Millihertz Gravitational Wave Sky

Authors: John Baker, Jillian Bellovary, Peter L. Bender, Emanuele Berti, Robert Caldwell, Jordan Camp, John W. Conklin, Neil Cornish, Curt Cutler, Ryan DeRosa, Michael Eracleous, Elizabeth C. Ferrara, Samuel Francis, Martin Hewitson, Kelly Holley-Bockelmann, Ann Hornschemeier, Craig Hogan, Brittany Kamai, Bernard J. Kelly, Joey Shapiro Key, Shane L. Larson, Jeff Livas, Sridhar Manthripragada, Kirk McKenzie, Sean T. McWilliams , et al. (17 additional authors not shown)

Abstract: The first terrestrial gravitational wave interferometers have dramatically underscored the scientific value of observing the Universe through an entirely different window, and of folding this new channel of information with traditional astronomical data for a multimessenger view. The Laser Interferometer Space Antenna (LISA) will broaden the reach of gravitational wave astronomy by conducting the… ▽ More The first terrestrial gravitational wave interferometers have dramatically underscored the scientific value of observing the Universe through an entirely different window, and of folding this new channel of information with traditional astronomical data for a multimessenger view. The Laser Interferometer Space Antenna (LISA) will broaden the reach of gravitational wave astronomy by conducting the first survey of the millihertz gravitational wave sky, detecting tens of thousands of individual astrophysical sources ranging from white-dwarf binaries in our own galaxy to mergers of massive black holes at redshifts extending beyond the epoch of reionization. These observations will inform - and transform - our understanding of the end state of stellar evolution, massive black hole birth, and the co-evolution of galaxies and black holes through cosmic time. LISA also has the potential to detect gravitational wave emission from elusive astrophysical sources such as intermediate-mass black holes as well as exotic cosmological sources such as inflationary fields and cosmic string cusps. △ Less

Submitted 26 July, 2019; v1 submitted 15 July, 2019; originally announced July 2019.

Comments: White Paper submitted to Astro2020 (2020 Decadal Survey on Astronomy and Astrophysics). v2: fixed a reference

Understanding extreme quasar optical variability with CRTS: II. Changing-state quasars

Authors: Matthew J. Graham, Nicholas P. Ross, Daniel Stern, rew J. Drake, Barry McKernan, K. E. Saavik Ford, S. G. Djorgovski, Ashish Mahabal, Eilat Glikman, Steve Larson, Eric Christensen

Abstract: We present the results of a systematic search for quasars in the Catalina Real-time Transient Survey exhibiting both strong photometric and spectroscopic variability over a decadal baseline. We identify 73 sources with specific patterns of optical and mid-IR photometric behavior and a defined spectroscopic change. These "Changing-State" quasars (CSQs) form a higher luminosity sample to complement… ▽ More We present the results of a systematic search for quasars in the Catalina Real-time Transient Survey exhibiting both strong photometric and spectroscopic variability over a decadal baseline. We identify 73 sources with specific patterns of optical and mid-IR photometric behavior and a defined spectroscopic change. These "Changing-State" quasars (CSQs) form a higher luminosity sample to complement existing sets of "Changing-Look" AGN and quasars in the literature. The CSQs (by selection) exhibit larger photometric variability than the CLQs. The spectroscopic variability is marginally stronger in the CSQs than CLQs as defined by the change in H$β$/[OIII] ratio. We find 36 sources with declining H$β$ flux, 37 sources with increasing H$β$ flux and discover seven sources with $z > 0.8$, further extending the redshift arm. Our CSQ sample compares to the literature CLQ objects in similar distributions of H$β$ flux ratios and differential Eddington ratios between high (bright) and low (dim) states. Taken as a whole, we find that this population of extreme varying quasars is associated with changes in the Eddington ratio and the timescales imply cooling/heating fronts propagating through the disk. △ Less

Submitted 6 May, 2019; originally announced May 2019.

Comments: 43 pages, 22 figures, submitted

Astro2020 Science White Paper: Populations of Black Holes in Binaries

Authors: Thomas J. Maccarone, Laura Chomiuk, James Miller-Jones, Eric Bellm, Katelyn Breivik, Chris L. Fryer, Vicky Kalogera, Shane Larson, Jerome Orosz, James F. Steiner, Jay Strader, John A. Tomsick

Abstract: Black holes in binary star systems are vital for understanding the process of pr oducing gravitational wave sources, understanding how supernovae work, and for p roviding fossil evidence for the high mass stars from earlier in the Universe. At the present time, sample sizes of these objects, and especially of black hole s in binaries, are quite limited. Furthermore, more precise measurements of th… ▽ More Black holes in binary star systems are vital for understanding the process of pr oducing gravitational wave sources, understanding how supernovae work, and for p roviding fossil evidence for the high mass stars from earlier in the Universe. At the present time, sample sizes of these objects, and especially of black hole s in binaries, are quite limited. Furthermore, more precise measurements of the binary parameters are needed, as well. With improvements primarily in X-ray an d radio astronomy capabilities, it should be possible to build much larger sampl es of much better measured black hole binaries. △ Less

Submitted 26 April, 2019; originally announced April 2019.

Comments: white paper submitted to Astro2020 Decadal Survey

The Discovery Potential of Space-Based Gravitational Wave Astronomy

Authors: Neil J. Cornish, Emanuele Berti, Kelly Holley-Bockelmann, Shane Larson, Sean McWilliams, Guido Mueller, Priya Natarajan, Michele Vallisneri

Abstract: A space-based interferometer operating in the previously unexplored mHz gravitational band has tremendous discovery potential. If history is any guide, the opening of a new spectral band will lead to the discovery of entirely new sources and phenomena. The mHz band is ideally suited to exploring beyond standard model processes in the early universe, and with the sensitivities that can be reached w… ▽ More A space-based interferometer operating in the previously unexplored mHz gravitational band has tremendous discovery potential. If history is any guide, the opening of a new spectral band will lead to the discovery of entirely new sources and phenomena. The mHz band is ideally suited to exploring beyond standard model processes in the early universe, and with the sensitivities that can be reached with current technologies, the discovery space for exotic astrophysical systems is vast. △ Less

Submitted 3 April, 2019; v1 submitted 2 April, 2019; originally announced April 2019.

Comments: White paper submitted to the Astro2020 decadal survey

Multi-Messenger Astrophysics with Pulsar Timing Arrays

Authors: Luke Zoltan Kelley, Maria Charisi, Sarah Burke-Spolaor, Joseph Simon, Laura Blecha, Tamara Bogdanovic, Monica Colpi, Julie Comerford, Daniel J. D'Orazio, Massimo Dotti, Michael Eracleous, Matthew Graham, Jenny E. Greene, Zoltán Haiman, Kelly Holley-Bockelmann, Erin Kara, Bernard Kelly, S. Komossa, Shane L. Larson, Xin Liu, Chung-Pei Ma, Scott Noble, Vasileios Paschalidis, Roman R. Rafikov, Vikram Ravi , et al. (7 additional authors not shown)

Abstract: Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will… ▽ More Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts. △ Less

Submitted 18 March, 2019; originally announced March 2019.

Comments: Submitted to the Astro2020 decadal review. This is one of five core white papers written by members of the NANOGrav Collaboration, see also: J.Cordes et al., S.R.Taylor et al., X.Siemens et al., and E.Fonseca et al

Astro2020 science white paper: The gravitational wave view of massive black holes

Authors: Monica Colpi, Kelly Holley-Bockelmann, Tamara Bogdanovic, Priya Natarajan, Jillian Bellovary, Alberto Sesana, Michael Tremmel, Jeremy Schnittman, Julia Comerford, Enrico Barausse, Emanuele Berti, Marta Volonteri, Fazeel Khan, Sean McWilliams, Sarah Burke-Spolaor, Jeff Hazboun, John Conklin, Guido Mueller, Shane Larson

Abstract: Coalescing, massive black-hole (MBH) binaries are the most powerful sources of gravitational waves (GWs) in the Universe, which makes MBH science a prime focus for ongoing and upcoming GW observatories. The Laser Interferometer Space Antenna (LISA) -- a gigameter scale space-based GW observatory -- will grant us access to an immense cosmological volume, revealing MBHs merging when the first cosmic… ▽ More Coalescing, massive black-hole (MBH) binaries are the most powerful sources of gravitational waves (GWs) in the Universe, which makes MBH science a prime focus for ongoing and upcoming GW observatories. The Laser Interferometer Space Antenna (LISA) -- a gigameter scale space-based GW observatory -- will grant us access to an immense cosmological volume, revealing MBHs merging when the first cosmic structures assembled in the Dark Ages. LISA will unveil the yet unknown origin of the first quasars, and detect the teeming population of MBHs of $10^4 - 10^7$ solar masses. forming within protogalactic halos. The Pulsar Timing Array, a galactic-scale GW survey, can access the largest MBHs the Universe, detecting the cosmic GW foreground from inspiraling MBH binaries of about 10^9 solar masses. LISA can measure MBH spins and masses with precision far exceeding that from electromagnetic (EM) probes, and together, both GW observatories will provide the first full census of binary MBHs, and their orbital dynamics, across cosmic time. Detecting the loud gravitational signal of these MBH binaries will also trigger alerts for EM counterpart searches, from decades (PTAs) to hours (LISA) prior to the final merger. By witnessing both the GW and EM signals of MBH mergers, precious information will be gathered about the rich and complex environment in the aftermath of a galaxy collision. The unique GW characterization of MBHs will shed light on the deep link between MBHs of $10^4-10^{10}$ solar masses and the grand design of galaxy assembly, as well as on the complex dynamics that drive MBHs to coalescence. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: 5 pages, 2 figures submitted to the Astro2020 Decadal

Astro2020 Decadal Science White Paper: Gravitational Wave Survey of Galactic Ultra Compact Binaries

Authors: Tyson B. Littenberg, Katelyn Breivik, Warren R. Brown, Michael Eracleous, J. J. Hermes, Kelly Holley-Bockelmann, Kyle Kremer, Thomas Kupfer, Shane L. Larson

Abstract: Ultra-compact binaries (UCBs) are systems containing compact or degenerate stars with orbital periods less than one hour. Tens of millions of UCBs are predicted to exist within theGalaxy emitting gravitational waves (GWs) at mHz frequencies. Combining GW searches with electromagnetic (EM) surveys like Gaia and LSST will yield a comprehensive, multimessenger catalog of UCBs in the galaxy. Joint EM… ▽ More Ultra-compact binaries (UCBs) are systems containing compact or degenerate stars with orbital periods less than one hour. Tens of millions of UCBs are predicted to exist within theGalaxy emitting gravitational waves (GWs) at mHz frequencies. Combining GW searches with electromagnetic (EM) surveys like Gaia and LSST will yield a comprehensive, multimessenger catalog of UCBs in the galaxy. Joint EM and GW observations enable measurements of masses, radii, and orbital dynamics far beyond what can be achieved by independent EM or GW studies. GW+EM surveys of UCBs in the galaxy will yield a trove of unique insights into the nature of white dwarfs, the formation of compact objects, dynamical interactions in binaries, and energetic, accretion-driven phenomena like Type Ia superonovae. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: 5 pages, 1 figure, submitted to the Astro2020 call for science white papers

Astro2020 Decadal Science White Paper: The state of gravitational-wave astrophysics in 2020

Authors: Sean T. McWilliams, Robert Caldwell, Kelly Holley-Bockelmann, Shane L. Larson, Michele Vallisneri

Abstract: While still in its infancy, the budding field of gravitational-wave astronomy has so far exceeded most expectations, and the achievements that have already been made bode well for the decade to come. While the discoveries made possible by LIGO have captured the imagination of experts and nonexperts alike, it is important when looking ahead to consider those discoveries in the context of the field… ▽ More While still in its infancy, the budding field of gravitational-wave astronomy has so far exceeded most expectations, and the achievements that have already been made bode well for the decade to come. While the discoveries made possible by LIGO have captured the imagination of experts and nonexperts alike, it is important when looking ahead to consider those discoveries in the context of the field as a whole. Just as radio, optical, and x-ray radiation probe different physical phenomena occurring on a range of length and energy scales, the future of gravitational-wave astrophysics depends on our ability to open up the entire spectrum. We will describe the scientific prospects for the field of gravitational-wave astronomy as a whole as we enter the coming decade, and we will place the specific contributions from a future space-based gravitational-wave observatory within this context. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: submission to the 2020-2030 Astronomy and Astrophysics Decadal Survey (Astro2020)

Multimessenger science opportunities with mHz gravitational waves

Authors: John Baker, Zoltán Haiman, Elena Maria Rossi, Edo Berger, Niel Brandt, Elmé Breedt, Katelyn Breivik, Maria Charisi, Andrea Derdzinski, Daniel J. D'Orazio, Saavik Ford, Jenny E. Greene, J. Colin Hill, Kelly Holley-Bockelmann, Joey Shapiro Key, Bence Kocsis, Thomas Kupfer, Shane Larson, Piero Madau, Thomas Marsh, Barry McKernan, Sean T. McWilliams, Priyamvada Natarajan, Samaya Nissanke, Scott Noble , et al. (10 additional authors not shown)

Abstract: LISA will open the mHz band of gravitational waves (GWs) to the astronomy community. The strong gravity which powers the variety of GW sources in this band is also crucial in a number of important astrophysical processes at the current frontiers of astronomy. These range from the beginning of structure formation in the early universe, through the origin and cosmic evolution of massive black holes… ▽ More LISA will open the mHz band of gravitational waves (GWs) to the astronomy community. The strong gravity which powers the variety of GW sources in this band is also crucial in a number of important astrophysical processes at the current frontiers of astronomy. These range from the beginning of structure formation in the early universe, through the origin and cosmic evolution of massive black holes in concert with their galactic environments, to the evolution of stellar remnant binaries in the Milky Way and in nearby galaxies. These processes and their associated populations also drive current and future observations across the electromagnetic (EM) spectrum. We review opportunities for science breakthroughs, involving either direct coincident EM+GW observations, or indirect multimessenger studies. We argue that for the US community to fully capitalize on the opportunities from the LISA mission, the US efforts should be accompanied by a coordinated and sustained program of multi-disciplinary science investment, following the GW data through to its impact on broad areas of astrophysics. Support for LISA-related multimessenger observers and theorists should be sized appropriately for a flagship observatory and may be coordinated through a dedicated mHz GW research center. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: Submitted to the Astro2020 Decadal Survey call for science white papers

What we can learn from multi-band observations of black hole binaries

Authors: Curt Cutler, Emanuele Berti, Karan Jani, Ely D. Kovetz, Lisa Randall, Salvatore Vitale, Kaze W. K. Wong, Kelly Holley-Bockelmann, Shane L. Larson, Tyson Littenberg, Sean T. McWilliams, Guido Mueller, Jeremy D. Schnittman, David H. Shoemaker, Michele Vallisneri

Abstract: The LIGO/Virgo gravitational-wave (GW) interferometers have to-date detected ten merging black hole (BH) binaries, some with masses considerably larger than had been anticipated. Stellar-mass BH binaries at the high end of the observed mass range (with "chirp mass" ${\cal M} \gtrsim 25 M_{\odot}$) should be detectable by a space-based GW observatory years before those binaries become visible to gr… ▽ More The LIGO/Virgo gravitational-wave (GW) interferometers have to-date detected ten merging black hole (BH) binaries, some with masses considerably larger than had been anticipated. Stellar-mass BH binaries at the high end of the observed mass range (with "chirp mass" ${\cal M} \gtrsim 25 M_{\odot}$) should be detectable by a space-based GW observatory years before those binaries become visible to ground-based GW detectors. This white paper discusses some of the synergies that result when the same binaries are observed by instruments in space and on the ground. We consider intermediate-mass black hole binaries (with total mass $M \sim 10^2 -10^4 M_{\odot}$) as well as stellar-mass black hole binaries. We illustrate how combining space-based and ground-based data sets can break degeneracies and thereby improve our understanding of the binary's physical parameters. While early work focused on how space-based observatories can forecast precisely when some mergers will be observed on the ground, the reverse is also important: ground-based detections will allow us to "dig deeper" into archived, space-based data to confidently identify black hole inspirals whose signal-to-noise ratios were originally sub-threshold, increasing the number of binaries observed in both bands by a factor of $\sim 4 - 7$. △ Less

Submitted 10 March, 2019; originally announced March 2019.

Comments: White Paper submitted to Astro2020 (2020 Decadal Survey on Astronomy and Astrophysics)

Detecting circumbinary exoplanets and hierarchical stellar triples with the LISA gravitational radiation mission

Authors: Jason H. Steffen, Dong-Hong Wu, Shane L. Larson

Abstract: We investigate the possibility of detecting planetary or stellar companions orbiting white dwarf binaries using the LISA gravitational radiation detector. Specifically, we consider the acceleration of the barycenter of the white dwarf binary due to the orbiting third body as well as the effect of changes in the tidal field across the binary due to the perturber's eccentric orbit. We find that the… ▽ More We investigate the possibility of detecting planetary or stellar companions orbiting white dwarf binaries using the LISA gravitational radiation detector. Specifically, we consider the acceleration of the barycenter of the white dwarf binary due to the orbiting third body as well as the effect of changes in the tidal field across the binary due to the perturber's eccentric orbit. We find that the movement of the barycenter is detectable for both stellar and planetary mass objects. If circumbinary planets occur with frequencies similar to gas giant planets around isolated main sequence stars, then we expect to find of order 10 such planets in four years of LISA observations. For a longer, ten-year mission the accessible parameter space for planetary mass, orbital period, and binary orbital period grows and LISA's associated yield increases to ~100 expected detections. △ Less

Submitted 9 December, 2018; originally announced December 2018.

Comments: Submitted to AAS Journals

Active Asteroid P/2017 S5 (ATLAS)

Authors: David Jewitt, Yoonyoung Kim, Jayadev Rajagopal, Susan Ridgway, Ralf Kotulla, Wilson Liu, Max Mutchler, Jing Li, Harold Weaver, Stephen Larson

Abstract: Observations of active asteroid P/2017 S5 when near perihelion reveal the ejection of large (0.1 to 10 mm) particles at 0.2 to 2 m/s speeds, with estimated mass-loss rates of a few kg/s. The protracted nature of the mass loss (continuous over 150 days) is compatible with a sublimation origin, meaning that this object is likely an ice-bearing main-belt comet. Equilibrium sublimation of exposed wate… ▽ More Observations of active asteroid P/2017 S5 when near perihelion reveal the ejection of large (0.1 to 10 mm) particles at 0.2 to 2 m/s speeds, with estimated mass-loss rates of a few kg/s. The protracted nature of the mass loss (continuous over 150 days) is compatible with a sublimation origin, meaning that this object is likely an ice-bearing main-belt comet. Equilibrium sublimation of exposed water ice covering as little as 0.1 sq. km can match the data. Observations a year after perihelion show the object in an inactive state from which we deduce a nucleus effective radius 450(+100/-60) m (albedo 0.06+/-0.02 assumed). The gravitational escape speed from a body of this size is just 0.3 m/s, comparable to the inferred ejection speed of the dust. Time-series photometry provides tentative evidence for rapid rotation (lightcurve period 1.4 hour) that may also play a role in the loss of mass and which, if real, is a likely consequence of spin-up by sublimation torques. P/2017 S5 shares both physical and orbital similarities with the split active asteroid pair P/2016 J1-A and J1-B, and all three objects are likely members of the 7 Myr old, collisionally produced, Theobalda family. △ Less

Submitted 30 November, 2018; originally announced December 2018.

Comments: 15 pages, 6 figures

Post-Newtonian Dynamics in Dense Star Clusters: Binary Black Holes in the LISA Band

Authors: Kyle Kremer, Carl L. Rodriguez, Pau Amaro-Seoane, Katelyn Breivik, Sourav Chatterjee, Michael L. Katz, Shane L. Larson, Frederic A. Rasio, Johan Samsing, Claire S. Ye, Michael Zevin

Abstract: The dynamical processing of black holes in the dense cores of globular clusters (GCs), makes them efficient factories for producing binary black holes (BBHs). Here we explore the population of BBHs that form dynamically in GCs and may be observable at mHz frequencies or higher with LISA. We use our Monte Carlo stellar dynamics code, which includes gravitational radiation reaction effects for all B… ▽ More The dynamical processing of black holes in the dense cores of globular clusters (GCs), makes them efficient factories for producing binary black holes (BBHs). Here we explore the population of BBHs that form dynamically in GCs and may be observable at mHz frequencies or higher with LISA. We use our Monte Carlo stellar dynamics code, which includes gravitational radiation reaction effects for all BH encounters. By creating a representative local universe of GCs, we show that up to dozens of these systems may be resolvable by LISA with signal-to-noise ratios of at least 5. Approximately one third of these binaries will have measurable eccentricities ($e > 10^{-3}$) in the LISA band and a small number ($\lesssim 5$) may evolve from the LISA band to the LIGO band during the LISA mission. △ Less

Submitted 26 February, 2019; v1 submitted 28 November, 2018; originally announced November 2018.

Comments: Accepted for publication in Physical Review D

Journal ref: Phys. Rev. D 99, 063003 (2019)

Physical properties of near-Earth asteroids with a low delta-${\it v}$: Survey of target candidates for the Hayabusa2 mission

Authors: Sunao Hasegawa, Daisuke Kuroda, Kohei Kitazato, Toshihiro Kasuga, Tomohiko Sekiguchi, Naruhisa Takato, Kentaro Aoki, Akira Arai, Young-Jun Choi, Tetsuharu Fuse, Hidekazu Hanayama, Takashi Hattori, Hsiang-Yao Hsiao, Nobunari Kashikawa, Nobuyuki Kawai, Kyoko Kawakami, Daisuke Kinoshita, Steve Larson, Chi-Sheng Lin, Seidai Miyasaka, Naoya Miura, Shogo Nagayama, Yu Nagumo, Setsuko Nishihara, Yohei Ohba , et al. (18 additional authors not shown)

Abstract: Sample return from the near-Earth asteroid known as 25143 Itokawa was conducted as part of the Hayabusa mission, with a large number of scientific findings being derived from the returned samples. Following the Hayabusa mission, Hayabusa2 was planned, targeting sample return from a primitive asteroid. The primary target body of Hayabusa2 was asteroid 162173 Ryugu; however, it was also necessary to… ▽ More Sample return from the near-Earth asteroid known as 25143 Itokawa was conducted as part of the Hayabusa mission, with a large number of scientific findings being derived from the returned samples. Following the Hayabusa mission, Hayabusa2 was planned, targeting sample return from a primitive asteroid. The primary target body of Hayabusa2 was asteroid 162173 Ryugu; however, it was also necessary to gather physical information for backup target selection. Therefore, we examined five asteroids spectroscopically, 43 asteroids spectrophotometrically, and 41 asteroids through periodic analysis. Hence, the physical properties of 74 near-Earth asteroids were obtained, which helped the Hayabusa2 backup target search and, also, furthered understanding of the physical properties of individual asteroids and their origins. △ Less

Submitted 9 October, 2018; v1 submitted 8 October, 2018; originally announced October 2018.

Comments: 31 pages, 2 figures, and 4 tables, accepted for publication in PASJ

Results of a systematic search for outburst events in 1.4 million galaxies

Authors: A. J. Drake, S. G. Djorgovski, M. J. Graham, D. Stern, A. A. Mahabal, M. Catelan, E. Christensen, S. Larson

Abstract: We present an analysis of nine years of Catalina Surveys optical photometry for 1.4 million spectroscopically confirmed SDSS galaxies. We find 717 outburst events that were not reported by ongoing transient surveys. These events have timescales ranging from weeks to years. More than two thirds of these new events are found in starforming galaxies, while such galaxies only constitute ~20% of our sa… ▽ More We present an analysis of nine years of Catalina Surveys optical photometry for 1.4 million spectroscopically confirmed SDSS galaxies. We find 717 outburst events that were not reported by ongoing transient surveys. These events have timescales ranging from weeks to years. More than two thirds of these new events are found in starforming galaxies, while such galaxies only constitute ~20% of our sample. Based on the properties of the hosts and events, we find that almost all of the new events are likely to be associated with regular supernovae. However, a small number of long-timescale events are found among the galaxies containing AGN. These events have similar properties to those recently found in the analyses of light curves of large samples of AGN. Given the lack of such events among the more than a million passive galaxies in the sample, we suggest that the long outbursts are associated with super-massive black holes or their environments. △ Less

Submitted 2 October, 2018; originally announced October 2018.

Comments: 21 pages, 18 figures. MNRAS accepted

The Nucleus of Active Asteroid 311P/(2013 P5) PANSTARRS

Authors: David Jewitt, Harold Weaver, Max Mutchler, Jing Li, Jessica Agarwal, Stephen Larson

Abstract: The unique inner-belt asteroid 311P/PANSTARRS (formerly P/2013 P5) is notable for its sporadic, comet-like ejection of dust in nine distinct epochs spread over $\sim$250 days in 2013. This curious behavior has been interpreted as the product of localized, equator-ward landsliding from the surface of an asteroid rotating at the brink of instability. We obtained new Hubble Space Telescope observatio… ▽ More The unique inner-belt asteroid 311P/PANSTARRS (formerly P/2013 P5) is notable for its sporadic, comet-like ejection of dust in nine distinct epochs spread over $\sim$250 days in 2013. This curious behavior has been interpreted as the product of localized, equator-ward landsliding from the surface of an asteroid rotating at the brink of instability. We obtained new Hubble Space Telescope observations to directly measure the nucleus and to search for evidence of its rapid rotation. However, instead of providing photometric evidence for rapid nucleus rotation, our data set a lower limit to the lightcurve period, $P \ge$ 5.4 hour. The dominant feature of the lightcurve is a V-shaped minimum, $\sim$0.3 magnitudes deep, that is suggestive of an eclipsing binary. Under this interpretation, the time-series data are consistent with a secondary/primary mass ratio, $m_s/m_p \sim$ 1:6, a ratio of separation/primary radius, $r/r_p \sim$ 4 and an orbit period $\sim$0.8 days. These properties lie within the range of other asteroid binaries that are thought to be formed by rotational breakup. While the lightcurve period is long, centripetal dust ejection is still possible if one or both components rotates rapidly ($\lesssim$ 2 hour) and has a small lightcurve variation because of azimuthal symmetry. Indeed, radar observations of asteroids in critical rotation reveal "muffin-shaped" morphologies which are closely azimuthally symmetric and which show minimal lightcurves. Our data are consistent with 311P being a close binary in which one or both components rotates near the centripetal limit. The mass loss in 2013 suggests that breakup occurred recently and could even be on-going. A search for fragments that might have been recently ejected beyond the Hill sphere reveals none larger than effective radius $r_e \sim$ 10 m. △ Less

Submitted 16 April, 2018; originally announced April 2018.

Comments: 37 pages, 9 figures, Astronomical Journal, in press

LISA sources in Milky Way globular clusters

Authors: Kyle Kremer, Sourav Chatterjee, Katelyn Breivik, Carl L. Rodriguez, Shane L. Larson, Frederic A. Rasio

Abstract: We explore the formation of double-compact-object binaries in Milky Way (MW) globular clusters (GCs) that may be detectable by the Laser Interferometer Space Antenna (LISA). We use a set of 137 fully evolved GC models that, overall, effectively match the properties of the observed GCs in the MW. We estimate that, in total, the MW GCs contain $\sim21$ sources that will be detectable by LISA. These… ▽ More We explore the formation of double-compact-object binaries in Milky Way (MW) globular clusters (GCs) that may be detectable by the Laser Interferometer Space Antenna (LISA). We use a set of 137 fully evolved GC models that, overall, effectively match the properties of the observed GCs in the MW. We estimate that, in total, the MW GCs contain $\sim21$ sources that will be detectable by LISA. These detectable sources contain all combinations of black hole (BH), neutron star, and white dwarf components. We predict $\sim7$ of these sources will be BH-BH binaries. Furthermore, we show that some of these BH-BH binaries can have signal-to-noise ratios large enough to be detectable at the distance of the Andromeda galaxy or even the Virgo cluster. △ Less

Submitted 29 May, 2018; v1 submitted 15 February, 2018; originally announced February 2018.

Comments: 7 pages, 4 figures, 2 appendices. As published in Physical Review Letters

Journal ref: Phys. Rev. Lett. 120, 191103 (2018)

Characterizing accreting double white dwarf binaries with the Laser Interferometer Space Antenna and Gaia

Authors: Katelyn Breivik, Kyle Kremer, Michael Bueno, Shane L. Larson, Scott Coughlin, Vassiliki Kalogera

Abstract: We demonstrate a method to fully characterize mass-transferring double white dwarf (DWD) systems with a helium-rich (He) WD donor based on the mass--radius relationship for He WDs. Using a simulated Galactic population of DWDs, we show that donor and accretor masses can be inferred for up to $\sim\, 60$ systems observed by both Laser Interferometer Space Antenna (LISA) and Gaia. Half of these syst… ▽ More We demonstrate a method to fully characterize mass-transferring double white dwarf (DWD) systems with a helium-rich (He) WD donor based on the mass--radius relationship for He WDs. Using a simulated Galactic population of DWDs, we show that donor and accretor masses can be inferred for up to $\sim\, 60$ systems observed by both Laser Interferometer Space Antenna (LISA) and Gaia. Half of these systems will have mass constraints $Δ\,M_{\rm{D}}\lesssim0.2M_{\odot}$ and $Δ\,M_{\rm{A}}\lesssim2.3\,M_{\odot}$. We also show how the orbital frequency evolution due to astrophysical processes and gravitational radiation can be decoupled from the total orbital frequency evolution for up to $\sim 50$ of these systems. △ Less

Submitted 30 April, 2018; v1 submitted 23 October, 2017; originally announced October 2017.

Comments: Updated to match publication in ApJ Letters; 8 pages, 4 figures

Journal ref: 2018ApJ...854L...1B

On the Progenitor of Binary Neutron Star Merger GW170817

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca , et al. (1073 additional authors not shown)

Abstract: On 2017 August 17 the merger of two compact objects with masses consistent with two neutron stars was discovered through gravitational-wave (GW170817), gamma-ray (GRB 170817A), and optical (SSS17a/AT 2017gfo) observations. The optical source was associated with the early-type galaxy NGC 4993 at a distance of just $\sim$40 Mpc, consistent with the gravitational-wave measurement, and the merger was… ▽ More On 2017 August 17 the merger of two compact objects with masses consistent with two neutron stars was discovered through gravitational-wave (GW170817), gamma-ray (GRB 170817A), and optical (SSS17a/AT 2017gfo) observations. The optical source was associated with the early-type galaxy NGC 4993 at a distance of just $\sim$40 Mpc, consistent with the gravitational-wave measurement, and the merger was localized to be at a projected distance of $\sim$2 kpc away from the galaxy's center. We use this minimal set of facts and the mass posteriors of the two neutron stars to derive the first constraints on the progenitor of GW170817 at the time of the second supernova (SN). We generate simulated progenitor populations and follow the three-dimensional kinematic evolution from the binary neutron star (BNS) birth to the merger time, accounting for pre-SN galactic motion, for considerably different input distributions of the progenitor mass, pre-SN semimajor axis, and SN-kick velocity. Though not considerably tight, we find these constraints to be comparable to those for Galactic BNS progenitors. The derived constraints are very strongly influenced by the requirement of keeping the binary bound after the second SN and having the merger occur relatively close to the center of the galaxy. These constraints are insensitive to the galaxy's star formation history, provided the stellar populations are older than 1 Gyr. △ Less

Submitted 1 December, 2017; v1 submitted 16 October, 2017; originally announced October 2017.

Comments: 18 pages, 8 figures, 2 tables

Report number: LIGO-P1700264

Journal ref: ApJL, 850, L40, 2017

Estimating the Contribution of Dynamical Ejecta in the Kilonova Associated with GW170817

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca , et al. (1078 additional authors not shown)

Abstract: The source of the gravitational-wave signal GW170817, very likely a binary neutron star merger, was also observed electromagnetically, providing the first multi-messenger observations of this type. The two week long electromagnetic counterpart had a signature indicative of an r-process-induced optical transient known as a kilonova. This Letter examines how the mass of the dynamical ejecta can be e… ▽ More The source of the gravitational-wave signal GW170817, very likely a binary neutron star merger, was also observed electromagnetically, providing the first multi-messenger observations of this type. The two week long electromagnetic counterpart had a signature indicative of an r-process-induced optical transient known as a kilonova. This Letter examines how the mass of the dynamical ejecta can be estimated without a direct electromagnetic observation of the kilonova, using gravitational-wave measurements and a phenomenological model calibrated to numerical simulations of mergers with dynamical ejecta. Specifically, we apply the model to the binary masses inferred from the gravitational-wave measurements, and use the resulting mass of the dynamical ejecta to estimate its contribution (without the effects of wind ejecta) to the corresponding kilonova light curves from various models. The distributions of dynamical ejecta mass range between $M_{ej} = 10^{-3} - 10^{-2} M_{\odot}$ for various equations of state, assuming the neutron stars are rotating slowly. In addition, we use our estimates of the dynamical ejecta mass and the neutron star merger rates inferred from GW170817 to constrain the contribution of events like this to the r-process element abundance in the Galaxy when ejecta mass from post-merger winds is neglected. We find that if $\gtrsim10\%$ of the matter dynamically ejected from BNS mergers is converted to r-process elements, GW170817-like BNS mergers could fully account for the amount of r-process material observed in the Milky Way. △ Less

Submitted 8 January, 2018; v1 submitted 16 October, 2017; originally announced October 2017.

Comments: 18 pages, 5 figures, accepted for publication in ApJ Letters

Report number: LIGO-P1700309

Revealing black holes with Gaia

Authors: Katelyn Breivik, Sourav Chatterjee, Shane L. Larson

Abstract: We estimate the population of black holes with luminous stellar companions (BH-LCs) in the Milky Way (MW) observable by Gaia. We evolve a realistic distribution of BH-LC progenitors from zero-age to the current epoch taking into account relevant physics, including binary stellar evolution, BH-formation physics, and star formation rate, to estimate the BH-LC population in the MW today. We predict t… ▽ More We estimate the population of black holes with luminous stellar companions (BH-LCs) in the Milky Way (MW) observable by Gaia. We evolve a realistic distribution of BH-LC progenitors from zero-age to the current epoch taking into account relevant physics, including binary stellar evolution, BH-formation physics, and star formation rate, to estimate the BH-LC population in the MW today. We predict that Gaia will discover between 3800 and 12,000 BH-LCs by the end of its 5 yr mission, depending on BH natal kick strength and observability constraints. We find that the overall yield, and distributions of eccentricities and masses of observed BH-LCs can provide important constraints on the strength of BH natal kicks. Gaia-detected BH-LCs are expected to have very different orbital properties compared to those detectable via radio, X-ray, or gravitational wave observations. △ Less

Submitted 14 November, 2017; v1 submitted 12 October, 2017; originally announced October 2017.

Comments: Accepted for publication in the Astrophysical Journal Letters; 8 pages, 4 figures, 1 table; Comments welcome

A binary main belt comet

Authors: Jessica Agarwal, David Jewitt, Max Mutchler, Harold Weaver, Stephen Larson

Abstract: The asteroids are primitive solar system bodies which evolve both collisionally and through disruptions due to rapid rotation [1]. These processes can lead to the formation of binary asteroids [2-4] and to the release of dust [5], both directly and, in some cases, through uncovering frozen volatiles. In a sub-set of the asteroids called main-belt comets (MBCs), the sublimation of excavated volatil… ▽ More The asteroids are primitive solar system bodies which evolve both collisionally and through disruptions due to rapid rotation [1]. These processes can lead to the formation of binary asteroids [2-4] and to the release of dust [5], both directly and, in some cases, through uncovering frozen volatiles. In a sub-set of the asteroids called main-belt comets (MBCs), the sublimation of excavated volatiles causes transient comet-like activity [6-8]. Torques exerted by sublimation measurably influence the spin rates of active comets [9] and might lead to the splitting of bilobate comet nuclei [10]. The kilometer-sized main-belt asteroid 288P (300163) showed activity for several months around its perihelion 2011 [11], suspected to be sustained by the sublimation of water ice [12] and supported by rapid rotation [13], while at least one component rotates slowly with a period of 16 hours [14]. 288P is part of a young family of at least 11 asteroids that formed from a ~10km diameter precursor during a shattering collision 7.5 million years ago [15]. Here we report that 288P is a binary main-belt comet. It is different from the known asteroid binaries for its combination of wide separation, near-equal component size, high eccentricity, and comet-like activity. The observations also provide strong support for sublimation as the driver of activity in 288P and show that sublimation torques may play a significant role in binary orbit evolution. △ Less

Submitted 10 October, 2017; originally announced October 2017.

Comments: 18 pages, 8 figures, 2 tables

Journal ref: Nature 549, pp. 357-359 (2017)

Accreting Double white dwarf binaries: Implications for LISA

Authors: Kyle Kremer, Katelyn Breivik, Shane L. Larson, Vassiliki Kalogera

Abstract: We explore the long-term evolution of mass-transferring white dwarf binaries undergoing both direct-impact and disk accretion and explore implications of such systems to gravitational wave astronomy. We cover a broad range of initial component masses and show that these systems, the majority of which lie within the LISA sensitivity range, exhibit prominent negative orbital frequency evolution (chi… ▽ More We explore the long-term evolution of mass-transferring white dwarf binaries undergoing both direct-impact and disk accretion and explore implications of such systems to gravitational wave astronomy. We cover a broad range of initial component masses and show that these systems, the majority of which lie within the LISA sensitivity range, exhibit prominent negative orbital frequency evolution (chirp) for a significant fraction of their lifetimes. Using a galactic population synthesis, we predict ~$2700$ double white dwarfs will be observable by LISA with negative chirps less than $-0.1\ yr^{-2}$. We also show that detections of mass-transferring double white dwarf systems by LISA may provide astronomers with unique ways of probing the physics governing close compact object binaries. △ Less

Submitted 12 September, 2017; v1 submitted 4 July, 2017; originally announced July 2017.

Comments: Accepted in The Astrophysical Journal

Journal ref: The Astrophysical Journal, 846:95 (11pp), 2017 September 10

Understanding extreme quasar optical variability with CRTS: I. Major AGN flares

Authors: Matthew J. Graham, S. G. Djorgovski, Andrew J. Drake, Daniel Stern, Ashish A. Mahabal, Eilat Glikman, Steve Larson, Eric Christensen

Abstract: There is a large degree of variety in the optical variability of quasars and it is unclear whether this is all attributable to a single (set of) physical mechanism(s). We present the results of a systematic search for major flares in AGN in the Catalina Real-time Transient Survey as part of a broader study into extreme quasar variability. Such flares are defined in a quantitative manner as being a… ▽ More There is a large degree of variety in the optical variability of quasars and it is unclear whether this is all attributable to a single (set of) physical mechanism(s). We present the results of a systematic search for major flares in AGN in the Catalina Real-time Transient Survey as part of a broader study into extreme quasar variability. Such flares are defined in a quantitative manner as being atop of the normal, stochastic variability of quasars. We have identified 51 events from over 900,000 known quasars and high probability quasar candidates, typically lasting 900 days and with a median peak amplitude of $Δm = 1.25$ mag. Characterizing the flare profile with a Weibull distribution, we find that nine of the sources are well described by a single-point single-lens model. This supports the proposal by Lawrence et al. (2016) that microlensing is a plausible physical mechanism for extreme variability. However, we attribute the majority of our events to explosive stellar-related activity in the accretion disk: superluminous supernovae, tidal disruption events, and mergers of stellar mass black holes. △ Less

Submitted 9 June, 2017; originally announced June 2017.

Comments: 25 pages, 18 figures, accepted for publication by MNRAS

Anatomy of an Asteroid Break-Up: The Case of P/2013 R3

Authors: Dave Jewitt, Jessica Agarwal, Jing Li, Harold Weaver, Max Mutchler, Stephen Larson

Abstract: We present an analysis of new and published data on P/2013 R3, the first asteroid detected while disintegrating. Thirteen discrete components are measured in the interval between UT 2013 October 01 and 2014 February 13. We determine a mean, pair-wise velocity dispersion amongst these components of $Δv = 0.33\pm0.03$ m s$^{-1}$ and find that their separation times are staggered over an interval of… ▽ More We present an analysis of new and published data on P/2013 R3, the first asteroid detected while disintegrating. Thirteen discrete components are measured in the interval between UT 2013 October 01 and 2014 February 13. We determine a mean, pair-wise velocity dispersion amongst these components of $Δv = 0.33\pm0.03$ m s$^{-1}$ and find that their separation times are staggered over an interval of $\sim$5 months. Dust enveloping the system has, in the first observations, a cross-section $\sim$30 km$^2$ but fades monotonically at a rate consistent with the action of radiation pressure sweeping. The individual components exhibit comet-like morphologies and also fade except where secondary fragmentation is accompanied by the release of additional dust. We find only upper limits to the radii of any embedded solid nuclei, typically $\sim$100 to 200 m (geometric albedo 0.05 assumed). Combined, the components of P/2013 R3 would form a single spherical body with radius $\lesssim$400 m, which is our best estimate of the size of the precursor object. The observations are consistent with rotational disruption of a weak (cohesive strength $\sim$50 to 100 N m$^{-2}$) parent body, $\sim$400 m in radius. Estimated radiation (YORP) spin-up times of this parent are $\lesssim$1 Myr, shorter than the collisional lifetime. If present, water ice sublimating at as little as 10$^{-3}$ kg s$^{-1}$ could generate a torque on the parent body rivaling the YORP torque. Under conservative assumptions about the frequency of similar disruptions, the inferred asteroid debris production rate is $\gtrsim$10$^3$ kg s$^{-1}$, which is at least 4% of the rate needed to maintain the Zodiacal Cloud. △ Less

Submitted 28 March, 2017; originally announced March 2017.

Comments: 44 pages, 13 figures, accepted by Astronomical Journal

Laser Interferometer Space Antenna

Authors: Pau Amaro-Seoane, Heather Audley, Stanislav Babak, John Baker, Enrico Barausse, Peter Bender, Emanuele Berti, Pierre Binetruy, Michael Born, Daniele Bortoluzzi, Jordan Camp, Chiara Caprini, Vitor Cardoso, Monica Colpi, John Conklin, Neil Cornish, Curt Cutler, Karsten Danzmann, Rita Dolesi, Luigi Ferraioli, Valerio Ferroni, Ewan Fitzsimons, Jonathan Gair, Lluis Gesa Bote, Domenico Giardini , et al. (59 additional authors not shown)

Abstract: Following the selection of The Gravitational Universe by ESA, and the successful flight of LISA Pathfinder, the LISA Consortium now proposes a 4 year mission in response to ESA's call for missions for L3. The observatory will be based on three arms with six active laser links, between three identical spacecraft in a triangular formation separated by 2.5 million km. LISA is an all-sky monitor and… ▽ More Following the selection of The Gravitational Universe by ESA, and the successful flight of LISA Pathfinder, the LISA Consortium now proposes a 4 year mission in response to ESA's call for missions for L3. The observatory will be based on three arms with six active laser links, between three identical spacecraft in a triangular formation separated by 2.5 million km. LISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using Gravitational Waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the infant Universe at TeV energy scales, has known sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales near the horizons of black holes, all the way to cosmological scales. The LISA mission will scan the entire sky as it follows behind the Earth in its orbit, obtaining both polarisations of the Gravitational Waves simultaneously, and will measure source parameters with astrophysically relevant sensitivity in a band from below $10^{-4}\,$Hz to above $10^{-1}\,$Hz. △ Less

Submitted 23 February, 2017; v1 submitted 2 February, 2017; originally announced February 2017.

Comments: Submitted to ESA on January 13th in response to the call for missions for the L3 slot in the Cosmic Vision Programme

MHz Gravitational Wave Constraints with Decameter Michelson Interferometers

Authors: Aaron S. Chou, Richard Gustafson, Craig Hogan, Brittany Kamai, Ohkyung Kwon, Robert Lanza, Shane L. Larson, Lee McCuller, Stephan S. Meyer, Jonathan Richardson, Chris Stoughton, Raymond Tomlin, Rainer Weiss

Abstract: A new detector, the Fermilab Holometer, consists of separate yet identical 39-meter Michelson interferometers. Strain sensitivity achieved is better than $10^{-21} /{\sqrt{\rm{Hz}}}$ between 1 to 13 MHz from a 130-hr dataset. This measurement exceeds the sensitivity and frequency range made from previous high frequency gravitational wave experiments by many orders of magnitude. Constraints are pla… ▽ More A new detector, the Fermilab Holometer, consists of separate yet identical 39-meter Michelson interferometers. Strain sensitivity achieved is better than $10^{-21} /{\sqrt{\rm{Hz}}}$ between 1 to 13 MHz from a 130-hr dataset. This measurement exceeds the sensitivity and frequency range made from previous high frequency gravitational wave experiments by many orders of magnitude. Constraints are placed on a stochastic background at 382 Hz resolution. The 3$σ$ upper limit on $Ω_{\rm{GW}}$, the gravitational wave energy density normalized to the closure density, ranges from $5.6 \times 10^{12}$ at 1 MHz to $8.4 \times 10^{15}$ at 13 MHz. Another result from the same dataset is a search for nearby primordial black hole binaries (PBHB). There are no detectable monochromatic PBHBs in the mass range $0.83$ - $3.5 \times 10^{21}$g between the earth and the moon. Projections for a chirp search with the same dataset increases the mass range to $0.59 - 2.5 \times 10^{25}$g and distances out to Jupiter. This result presents a new method for placing limits on a poorly constrained mass range of primordial black holes. Additionally, solar system searches for PBHBs place limits on their contribution to the total dark matter fraction. △ Less

Submitted 9 June, 2017; v1 submitted 16 November, 2016; originally announced November 2016.

Journal ref: Phys. Rev. D 95, 063002 (2017)

Gravity Spy: Integrating Advanced LIGO Detector Characterization, Machine Learning, and Citizen Science

Authors: Michael Zevin, Scott Coughlin, Sara Bahaadini, Emre Besler, Neda Rohani, Sarah Allen, Miriam Cabero, Kevin Crowston, Aggelos K Katsaggelos, Shane L Larson, Tae Kyoung Lee, Chris Lintott, Tyson B Littenberg, Andrew Lundgren, Carsten Oesterlund, Joshua R Smith, Laura Trouille, Vicky Kalogera

Abstract: (abridged for arXiv) With the first direct detection of gravitational waves, the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) has initiated a new field of astronomy by providing an alternate means of sensing the universe. The extreme sensitivity required to make such detections is achieved through exquisite isolation of all sensitive components of LIGO from non-gravitational… ▽ More (abridged for arXiv) With the first direct detection of gravitational waves, the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) has initiated a new field of astronomy by providing an alternate means of sensing the universe. The extreme sensitivity required to make such detections is achieved through exquisite isolation of all sensitive components of LIGO from non-gravitational-wave disturbances. Nonetheless, LIGO is still susceptible to a variety of instrumental and environmental sources of noise that contaminate the data. Of particular concern are noise features known as glitches, which are transient and non-Gaussian in their nature, and occur at a high enough rate so that accidental coincidence between the two LIGO detectors is non-negligible. In this paper we describe an innovative project that combines crowdsourcing with machine learning to aid in the challenging task of categorizing all of the glitches recorded by the LIGO detectors. Through the Zooniverse platform, we engage and recruit volunteers from the public to categorize images of glitches into pre-identified morphological classes and to discover new classes that appear as the detectors evolve. In addition, machine learning algorithms are used to categorize images after being trained on human-classified examples of the morphological classes. Leveraging the strengths of both classification methods, we create a combined method with the aim of improving the efficiency and accuracy of each individual classifier. The resulting classification and characterization should help LIGO scientists to identify causes of glitches and subsequently eliminate them from the data or the detector entirely, thereby improving the rate and accuracy of gravitational-wave observations. We demonstrate these methods using a small subset of data from LIGO's first observing run. △ Less

Submitted 28 February, 2017; v1 submitted 14 November, 2016; originally announced November 2016.

Comments: 27 pages, 8 figures, 1 table

Report number: LIGO-P1600303

Journal ref: Class. Quantum Grav. 34 (2017) 064003 (22pp)

Null-stream pointing with pulsar timing arrays

Authors: Jeffrey S. Hazboun, Shane L. Larson

Abstract: Locating sources on the sky is one of the largest challenges in gravitational wave astronomy, owing to the omni-directional nature of gravitational wave detection techniques, and the often intrinsically weak signals being observed. Ground-based detectors can address the pointing problem by observing with a network of detectors, effectively triangulating signal locations by observing the arrival ti… ▽ More Locating sources on the sky is one of the largest challenges in gravitational wave astronomy, owing to the omni-directional nature of gravitational wave detection techniques, and the often intrinsically weak signals being observed. Ground-based detectors can address the pointing problem by observing with a network of detectors, effectively triangulating signal locations by observing the arrival times across the network. Space-based detectors will observe long-lived sources that persist while the detector moves relative to their location on the sky, using Doppler shifts of the signal to locate the sky position. While these methods improve the pointing capability of a detector or network, the angular resolution is still coarse compared to the standards one expects from electromagnetic astronomy. Another technique that can be used for sky localization is null-stream pointing. In the case where multiple independent data streams exist, a single astrophysical source of gravitational waves will appear in each of the data streams. Taking the signals from multiple detectors in linear combination with each other, one finds there is a two parameter family of coefficients that effectively null the gravitational wave signal; those two parameters are the angles that define the sky location of the source. This technique has been demonstrated for a network of ground-based interferometric observatories, and for 6-link space interferometers. This paper derives and extends the null-stream pointing method to the unique case of pulsar timing residuals. The basic method is derived and demonstrated, and the necessity of using the method with multiple sub-arrays of pulsars in the pulsar timing array network is considered. △ Less

Submitted 12 July, 2016; originally announced July 2016.

Comments: 8 pages, 15 figures

Distinguishing Between Formation Channels for Binary Black Holes with LISA

Authors: Katelyn Breivik, Carl L. Rodriguez, Shane L. Larson, Vassiliki Kalogera, Frederic A. Rasio

Abstract: The recent detections of GW150914 and GW151226 imply an abundance of stellar-mass binary-black-hole mergers in the local universe. While ground-based gravitational-wave detectors are limited to observing the final moments before a binary merges, space-based detectors, such as the Laser Interferometer Space Antenna (LISA), can observe binaries at lower orbital frequencies where such systems may sti… ▽ More The recent detections of GW150914 and GW151226 imply an abundance of stellar-mass binary-black-hole mergers in the local universe. While ground-based gravitational-wave detectors are limited to observing the final moments before a binary merges, space-based detectors, such as the Laser Interferometer Space Antenna (LISA), can observe binaries at lower orbital frequencies where such systems may still encode information about their formation histories. In particular, the orbital eccentricity and mass of binary black holes in the LISA frequency band can be used together to discriminate between binaries formed in isolation in galactic fields and those formed in dense stellar environments such as globular clusters. In this letter, we explore the orbital eccentricity and mass of binary-black-hole populations as they evolve through the LISA frequency band. Overall we find that there are two distinct populations discernible by LISA. We show that up to ~90% of binaries formed either dynamically or in isolation have eccentricities measurable by LISA. Finally, we note how measured eccentricities of low-mass binary black holes evolved in isolation could provide detailed constraints on the physics of black-hole natal kicks and common-envelope evolution. △ Less

Submitted 22 September, 2016; v1 submitted 30 June, 2016; originally announced June 2016.

Comments: 7 pages, 4 figures, accepted for publication in ApJL

Hubble Space Telescope Observations of Active Asteroid 324P/La Sagra

Authors: David Jewitt, Jessica Agarwal, Harold Weaver, Max Mutchler, Jing Li, Stephen Larson

Abstract: Hubble Space Telescope observations of active asteroid 324P/La Sagra near perihelion show continued mass loss consistent with the sublimation of near-surface ice. Isophotes of the coma measured from a vantage point below the orbital plane are best matched by steady emission of particles having a nominal size $a \sim$ 100 $μ$m. The inferred rate of mass loss, $dM_d/dt \sim$0.2 kg s$^{-1}$, can be s… ▽ More Hubble Space Telescope observations of active asteroid 324P/La Sagra near perihelion show continued mass loss consistent with the sublimation of near-surface ice. Isophotes of the coma measured from a vantage point below the orbital plane are best matched by steady emission of particles having a nominal size $a \sim$ 100 $μ$m. The inferred rate of mass loss, $dM_d/dt \sim$0.2 kg s$^{-1}$, can be supplied by sublimation of water ice in thermal equilibrium with sunlight from an area as small as 930 m$^2$, corresponding to about 0.2\% of the nucleus surface. Observations taken from a vantage point only 0.6\degr~from the orbital plane of 324P set a limit to the velocity of ejection of dust in the direction perpendicular to the plane, $V_{\perp} <$ 1 m s$^{-1}$. Short-term photometric variations of the near-nucleus region, if related to rotation of the underlying nucleus, rule out periods $\le$ 3.8 hr and suggest that rotation probably does not play a central role in driving the observed mass loss. We estimate that, in the previous orbit, 324P lost about 4$\times$10$^7$ kg in dust particles, corresponding to 6$\times$10$^{-5}$ of the mass of a 550 m spherical nucleus of assumed density $ρ$ = 1000 kg m$^{-3}$. If continued, mass loss at this rate would limit the lifetime of 324P to $\sim$1.6$\times$10$^4$ orbits (about 10$^5$ yr). To survive for the 100 Myr to 400 Myr timescales corresponding, respectively, to dynamical and collisional stability requires a duty cycle $2\times 10^{-4} \le f_d \le 8\times 10^{-4}$. Unless its time in orbit is over-estimated by many orders of magnitude, 324P is revealed as a briefly-active member of a vast population of otherwise dormant ice-containing asteroids. △ Less

Submitted 27 June, 2016; originally announced June 2016.

Comments: 31 pages, 6 figures, 4 tables, Accepted by The Astronomical Journal

Hubble and Keck Telescope Observations of Active Asteroid 288P/300163 (2006 VW139)

Authors: Jessica Agarwal, David Jewitt, Harold Weaver, Max Mutchler, Stephen Larson

Abstract: We present Hubble Space Telescope and Keck 10 meter telescope observations of active asteroid 288P/300163 (2006 VW139) taken to examine ejected dust. The nucleus is a C-type object with absolute magnitude $H_V$ = 17.0$\pm$0.1 and estimated diameter $\sim$2.6 km (for assumed visual geometric albedo $p_V$ = 0.04). Variations in the brightness of the nucleus at the 10% to 15% level are significant in… ▽ More We present Hubble Space Telescope and Keck 10 meter telescope observations of active asteroid 288P/300163 (2006 VW139) taken to examine ejected dust. The nucleus is a C-type object with absolute magnitude $H_V$ = 17.0$\pm$0.1 and estimated diameter $\sim$2.6 km (for assumed visual geometric albedo $p_V$ = 0.04). Variations in the brightness of the nucleus at the 10% to 15% level are significant in both 2011 December and 2012 October but we possess too few data to distinguish variations caused by activity from those caused by rotation. The dust scattering cross-section in 2011 December is $\sim$40 km$^2$, corresponding to a dust mass $\sim$9$\times$10$^6$ kg (88 $μ$m mean particle radius assumed). The full width at half maximum of the debris sheet varies from $\sim$100 km near the nucleus to $\sim$1000 km 30arcsec (40,000 km) east of it. Dust dynamical models indicate ejection speeds between 0.06 and 0.3 m s$^{-1}$, particle sizes between 10 and 300 $μ$m and an inverse square-root relation between particle size and velocity. Overall, the data are most simply explained by prolonged, low velocity ejection of dust, starting in or before 2011 July and continuing until at least 2011 October. These properties are consistent with the sublimation of near-surface ice aided by centrifugal forces. The high spatial resolution of our HST images (52 km per pixel) reveals details that remained hidden in previous ground-based observations, such as the extraordinarily small vertical extent of the dust sheet, ejection speeds well below the nucleus escape speed, and the possibility of a binary nucleus. △ Less

Submitted 26 October, 2015; originally announced October 2015.

Comments: Accepted for publication by AJ

Nucleus and Mass Loss from Active Asteroid 313P/Gibbs

Authors: David Jewitt, Jing Li, Jessica Agarwal, Harold Weaver, Max Mutchler, Stephen Larson

Abstract: We present Hubble Space Telescope observations of active asteroid 313P/Gibbs (formerly P/2014 S4) taken over the five month interval from 2014 October to 2015 March. This object has been recurrently active near perihelion (at 2.4 AU) in two different orbits, a property that is naturally explained by the sublimation of near surface ice but which is difficult to reconcile with other activity mechani… ▽ More We present Hubble Space Telescope observations of active asteroid 313P/Gibbs (formerly P/2014 S4) taken over the five month interval from 2014 October to 2015 March. This object has been recurrently active near perihelion (at 2.4 AU) in two different orbits, a property that is naturally explained by the sublimation of near surface ice but which is difficult to reconcile with other activity mechanisms. We find that the mass loss peaks near 1 kg s$^{-1}$ in October and then declines over the subsequent months by about a factor of five, at nearly constant heliocentric distance. This decrease is too large to be caused by the change in heliocentric distance during the period of observation. However, it is consistent with sublimation from an ice patch shadowed by local topography, for example in a pit like those observed on the nuclei of short-period comet 67P/Churyumov-Gerasimenko. While no unique interpretation is possible, a simple self shadowing model shows that sublimation from a pit with depth to diameter ratio near 1/2 matches the observed rate of decline of the activity, while deeper and shallower pits do not. We estimate the nucleus radius to be 700$\pm$100 m (geometric albedo 0.05 assumed). Measurements of the spatial distribution of the dust were obtained from different viewing geometries. They show that dust was ejected continuously not impulsively, that the effective particle size is large, $\sim$50 $μm$, and that the ejection speed is $\sim$2.5 m s$^{-1}$. The total dust mass ejected is $\sim$10$^7$ kg, corresponding to $\sim$10$^{-5}$ of the nucleus mass. The observations are consistent with partially shadowed sublimation from $\sim$10$^4$ m$^2$ of ice, corresponding to $\sim$0.2\% of the nucleus surface. For ice to survive in 313P for billion-year timescales requires that the duty cycle for sublimation be $\lesssim$10$^{-3}$. △ Less

Submitted 7 July, 2015; originally announced July 2015.

Comments: 34 pages, 11 figures, 4 tables; Astronomical Journal: in press