Variational Integrators for Almost-Integrable Systems

We construct several variational integrators--integrators based on a discrete variational principle--for systems with Lagrangians of the form L = L_A + epsilon L_B, with epsilon << 1, where L_A describes an integrable system. These integrators exploit that epsilon << 1 to increase their accuracy by constructing discrete Lagrangians based on the assumption that the integrator trajectory is close to that of the integrable system. Several of the integrators we present are equivalent to well-known symplectic integrators for the equivalent perturbed Hamiltonian systems, but their construction and error analysis is significantly simpler in the variational framework. One novel method we present, involving a weighted time-averaging of the perturbing terms, removes all errors from the integration at O(epsilon). This last method is implicit, and involves evaluating a potentially expensive time-integral, but for some systems and some error tolerances it can significantly outperform traditional simulation methods.Comment: 14 pages, 4 figures. Version 2: added informative example; as accepted by Celestial Mechanics and Dynamical Astronom

Detection of Supernova Neutrinos by Neutrino-Proton Elastic Scattering

We propose that neutrino-proton elastic scattering, $\nu + p \to \nu + p$, can be used for the detection of supernova neutrinos in scintillator detectors. Though the proton recoil kinetic energy spectrum is soft, with $T_p \simeq 2 E_\nu^2/M_p$, and the scintillation light output from slow, heavily ionizing protons is quenched, the yield above a realistic threshold is nearly as large as that from $\bar{\nu}_e + p \to e^+ + n$. In addition, the measured proton spectrum is related to the incident neutrino spectrum, which solves a long-standing problem of how to separately measure the total energy and temperature of $\nu_\mu$, $\nu_\tau$, $\bar{\nu}_\mu$, and $\bar{\nu}_\tau$. The ability to detect this signal would give detectors like KamLAND and Borexino a crucial and unique role in the quest to detect supernova neutrinos.Comment: 10 pages, 9 figures, revtex

Black Hole Spin via Continuum Fitting and the Role of Spin in Powering Transient Jets

The spins of ten stellar black holes have been measured using the continuum-fitting method. These black holes are located in two distinct classes of X-ray binary systems, one that is persistently X-ray bright and another that is transient. Both the persistent and transient black holes remain for long periods in a state where their spectra are dominated by a thermal accretion disk component. The spin of a black hole of known mass and distance can be measured by fitting this thermal continuum spectrum to the thin-disk model of Novikov and Thorne; the key fit parameter is the radius of the inner edge of the black hole's accretion disk. Strong observational and theoretical evidence links the inner-disk radius to the radius of the innermost stable circular orbit, which is trivially related to the dimensionless spin parameter a_* of the black hole (|a_*| < 1). The ten spins that have so far been measured by this continuum-fitting method range widely from a_* \approx 0 to a_* > 0.95. The robustness of the method is demonstrated by the dozens or hundreds of independent and consistent measurements of spin that have been obtained for several black holes, and through careful consideration of many sources of systematic error. Among the results discussed is a dichotomy between the transient and persistent black holes; the latter have higher spins and larger masses. Also discussed is recently discovered evidence in the transient sources for a correlation between the power of ballistic jets and black hole spin.Comment: 30 pages. Accepted for publication in Space Science Reviews. Also to appear in hard cover in the Space Sciences Series of ISSI "The Physics of Accretion onto Black Holes" (Springer Publisher). Changes to Sections 5.2, 6.1 and 7.4. Section 7.4 responds to Russell et al. 2013 (MNRAS, 431, 405) who find no evidence for a correlation between the power of ballistic jets and black hole spi

All-sky search for long-duration gravitational wave transients with initial LIGO

We present the results of a search for long-duration gravitational wave transients in two sets of data collected by the LIGO Hanford and LIGO Livingston detectors between November 5, 2005 and September 30, 2007, and July 7, 2009 and October 20, 2010, with a total observational time of 283.0 days and 132.9 days, respectively. The search targets gravitational wave transients of duration 10-500 s in a frequency band of 40-1000 Hz, with minimal assumptions about the signal waveform, polarization, source direction, or time of occurrence. All candidate triggers were consistent with the expected background; as a result we set 90% confidence upper limits on the rate of long-duration gravitational wave transients for different types of gravitational wave signals. For signals from black hole accretion disk instabilities, we set upper limits on the source rate density between 3.4×10-5 and 9.4×10-4 Mpc-3 yr-1 at 90% confidence. These are the first results from an all-sky search for unmodeled long-duration transient gravitational waves. © 2016 American Physical Society

All-sky search for long-duration gravitational wave transients with initial LIGO

We present the results of a search for long-duration gravitational wave transients in two sets of data collected by the LIGO Hanford and LIGO Livingston detectors between November 5, 2005 and September 30, 2007, and July 7, 2009 and October 20, 2010, with a total observational time of 283.0 days and 132.9 days, respectively. The search targets gravitational wave transients of duration 10-500 s in a frequency band of 40-1000 Hz, with minimal assumptions about the signal waveform, polarization, source direction, or time of occurrence. All candidate triggers were consistent with the expected background; as a result we set 90% confidence upper limits on the rate of long-duration gravitational wave transients for different types of gravitational wave signals. For signals from black hole accretion disk instabilities, we set upper limits on the source rate density between 3.4×10-5 and 9.4×10-4 Mpc-3 yr-1 at 90% confidence. These are the first results from an all-sky search for unmodeled long-duration transient gravitational waves. © 2016 American Physical Society

Search for gravitational waves from Scorpius X-1 in the second Advanced LIGO observing run with an improved hidden Markov model

We present results from a semicoherent search for continuous gravitational waves from the low-mass x-ray binary Scorpius X-1, using a hidden Markov model (HMM) to track spin wandering. This search improves on previous HMM-based searches of LIGO data by using an improved frequency domain matched filter, the J-statistic, and by analyzing data from Advanced LIGO's second observing run. In the frequency range searched, from 60 to 650 Hz, we find no evidence of gravitational radiation. At 194.6 Hz, the most sensitive search frequency, we report an upper limit on gravitational wave strain (at 95% confidence) of h095%=3.47×10-25 when marginalizing over source inclination angle. This is the most sensitive search for Scorpius X-1, to date, that is specifically designed to be robust in the presence of spin wandering. © 2019 American Physical Society

Search for Tensor, Vector, and Scalar Polarizations in the Stochastic Gravitational-Wave Background

The detection of gravitational waves with Advanced LIGO and Advanced Virgo has enabled novel tests of general relativity, including direct study of the polarization of gravitational waves. While general relativity allows for only two tensor gravitational-wave polarizations, general metric theories can additionally predict two vector and two scalar polarizations. The polarization of gravitational waves is encoded in the spectral shape of the stochastic gravitational-wave background, formed by the superposition of cosmological and individually unresolved astrophysical sources. Using data recorded by Advanced LIGO during its first observing run, we search for a stochastic background of generically polarized gravitational waves. We find no evidence for a background of any polarization, and place the first direct bounds on the contributions of vector and scalar polarizations to the stochastic background. Under log-uniform priors for the energy in each polarization, we limit the energy densities of tensor, vector, and scalar modes at 95% credibility to Ω0T<5.58×10-8, Ω0V<6.35×10-8, and Ω0S<1.08×10-7 at a reference frequency f0=25 Hz. © 2018 American Physical Society

Erratum: "A Gravitational-wave Measurement of the Hubble Constant Following the Second Observing Run of Advanced LIGO and Virgo" (2021, ApJ, 909, 218)

[no abstract available

Search for Gravitational Waves Associated with Gamma-Ray Bursts Detected by Fermi and Swift during the LIGO-Virgo Run O3b

We search for gravitational-wave signals associated with gamma-ray bursts (GRBs) detected by the Fermi and Swift satellites during the second half of the third observing run of Advanced LIGO and Advanced Virgo (2019 November 1 15:00 UTC-2020 March 27 17:00 UTC). We conduct two independent searches: A generic gravitational-wave transients search to analyze 86 GRBs and an analysis to target binary mergers with at least one neutron star as short GRB progenitors for 17 events. We find no significant evidence for gravitational-wave signals associated with any of these GRBs. A weighted binomial test of the combined results finds no evidence for subthreshold gravitational-wave signals associated with this GRB ensemble either. We use several source types and signal morphologies during the searches, resulting in lower bounds on the estimated distance to each GRB. Finally, we constrain the population of low-luminosity short GRBs using results from the first to the third observing runs of Advanced LIGO and Advanced Virgo. The resulting population is in accordance with the local binary neutron star merger rate. © 2022. The Author(s). Published by the American Astronomical Society