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

Search for continuous gravitational waves from neutron stars in globular cluster NGC 6544

We describe a directed search for continuous gravitational waves in data from the sixth initial LIGO science run. The target was the nearby globular cluster NGC 6544 at a distance of ≈2.7 kpc. The search covered a broad band of frequencies along with first and second frequency derivatives for a fixed sky position. The search coherently integrated data from the two LIGO interferometers over a time span of 9.2 days using the matched-filtering F-statistic. We found no gravitational-wave signals and set 95% confidence upper limits as stringent as 6.0×10-25 on intrinsic strain and 8.5×10-6 on fiducial ellipticity. These values beat the indirect limits from energy conservation for stars with characteristic spin-down ages older than 300 years and are within the range of theoretical predictions for possible neutron-star ellipticities. An important feature of this search was use of a barycentric resampling algorithm which substantially reduced computational cost; this method is used extensively in searches of Advanced LIGO and Virgo detector data. © 2017 American Physical Society

Search for high-energy neutrinos from gravitational wave event GW151226 and candidate LVT151012 with ANTARES and IceCube

The Advanced LIGO observatories detected gravitational waves from two binary black hole mergers during their first observation run (O1). We present a high-energy neutrino follow-up search for the second gravitational wave event, GW151226, as well as for gravitational wave candidate LVT151012. We find two and four neutrino candidates detected by IceCube, and one and zero detected by Antares, within ±500 s around the respective gravitational wave signals, consistent with the expected background rate. None of these neutrino candidates are found to be directionally coincident with GW151226 or LVT151012. We use nondetection to constrain isotropic-equivalent high-energy neutrino emission from GW151226, adopting the GW event's 3D localization, to less than 2×1051-2×1054 erg. © 2017 American Physical Society