Numerical relativity simulation of GW150914 beyond general relativity

We produce the first astrophysically-relevant numerical binary black hole gravitational waveform in a higher-curvature theory of gravity beyond general relativity. We simulate a system with parameters consistent with GW150914, the first LIGO detection, in order-reduced dynamical Chern-Simons gravity, a theory with motivations in string theory and loop quantum gravity. We present results for the leading-order corrections to the merger and ringdown waveforms, as well as the ringdown quasi-normal mode spectrum. We estimate that such corrections may be discriminated in detections with signal to noise ratio $\gtrsim 180-240$, with the precise value depending on the dimension of the GR waveform family used in data analysis.Comment: 7 pages + appendices, 8 figures, Updated to match Phys. D. Rev articl

An analysis of the mechanisms governing species replacements in crayfish

We investigated mechanisms governing replacement of the native crayfish Orconectes sanborni by an invading crayfish, Orconectes rusticus. The two species had similar life histories, habitat preferences, and feeding patterns in allopatric and sympatric stream areas. Orconectes rusticus young-of-year (YOY) grew faster than O. sanborni YOY in the field. Adult O. rusticus were larger and, hence, dominant over adult O. sanborni; YOY were non-aggressive. In laboratory experiments, adult crayfish (about 28 mm carapace length or larger) were not susceptable to predation by largemouth bass (Micropterus salmoides, 30 cm total length) and did not alter shelter use when fish were present. Orconectes rusticus YOY were less susceptible to predation than O. sanborni YOY. Orconectes rusticus YOY reduced their vulnerability to largemouth bass by occupying shelters more often than YOY O. sanborni. In mixed-species mate-selection experiments, male O. rusticus and male O. sanborni preferentially mated with O. rusticus females. Inappropriate mate selection in sympatry may have caused the 90% reduction in recruitment for both species in 1982. Orconectes rusticus probably maintains greater population growth than O. sanborni, because (1) more gravid O. rusticus females occurred in sympatry, (2) O. rusticus produced more young than O. sanborni, and (3) O. rusticus young grew faster. Reproductive interference, acting synergistically with differences in aggressive dominance and young-of-year susceptibility to predation, appears to serve as the major mechanisms regulating replacement of O. sanborni by O. rusticus in Ohio streams.Funding for this project was provided by the Ohio Department of Natural Resources - Division of Natural Areas and Preserves, a Sigma Xi Grant-in-aid of Research, and an Ohio State University Grant-in-aid of Research

Numerical binary black hole collisions in dynamical Chern-Simons gravity

We produce the first numerical relativity binary black hole gravitational waveforms in a higher-curvature theory beyond general relativity. In particular, we study head-on collisions of binary black holes in order-reduced dynamical Chern-Simons gravity. This is a precursor to producing beyond-general-relativity waveforms for inspiraling binary black hole systems that are useful for gravitational wave detection. Head-on collisions are interesting in their own right, however, as they cleanly probe the quasinormal mode spectrum of the final black hole. We thus compute the leading-order dynamical Chern-Simons modifications to the complex frequencies of the postmerger gravitational radiation. We consider equal-mass systems, with equal spins oriented along the axis of collision, resulting in remnant black holes with spin. We find modifications to the complex frequencies of the quasinormal mode spectrum that behave as a power law with the spin of the remnant, and that are not degenerate with the frequencies associated with a Kerr black hole of any mass and spin. We discuss these results in the context of testing general relativity with gravitational wave observations

Surrogate models for precessing binary black hole simulations with unequal masses

Only numerical relativity simulations can capture the full complexities of binary black hole mergers. These simulations, however, are prohibitively expensive for direct data analysis applications such as parameter estimation. We present two new fast and accurate surrogate models for the outputs of these simulations: the first model, NRSur7dq4, predicts the gravitational waveform and the second model, \RemnantModel, predicts the properties of the remnant black hole. These models extend previous 7-dimensional, non-eccentric precessing models to higher mass ratios, and have been trained against 1528 simulations with mass ratios $q\leq4$ and spin magnitudes $\chi_1,\chi_2 \leq 0.8$, with generic spin directions. The waveform model, NRSur7dq4, which begins about 20 orbits before merger, includes all $\ell \leq 4$ spin-weighted spherical harmonic modes, as well as the precession frame dynamics and spin evolution of the black holes. The final black hole model, \RemnantModel, models the mass, spin, and recoil kick velocity of the remnant black hole. In their training parameter range, both models are shown to be more accurate than existing models by at least an order of magnitude, with errors comparable to the estimated errors in the numerical relativity simulations. We also show that the surrogate models work well even when extrapolated outside their training parameter space range, up to mass ratios $q=6$.Comment: Matches published version. Models publicly available at https://zenodo.org/record/3455886#.XZ9s1-dKjBI and https://pypi.org/project/surfinB

Incorporating Uncertainties in Atomic Data Into the Analysis of Solar and Stellar Observations: A Case Study in Fe XIII

Information about the physical properties of astrophysical objects cannot be measured directly but is inferred by interpreting spectroscopic observations in the context of atomic physics calculations. Ratios of emission lines, for example, can be used to infer the electron density of the emitting plasma. Similarly, the relative intensities of emission lines formed over a wide range of temperatures yield information on the temperature structure. A critical component of this analysis is understanding how uncertainties in the underlying atomic physics propagates to the uncertainties in the inferred plasma parameters. At present, however, atomic physics databases do not include uncertainties on the atomic parameters and there is no established methodology for using them even if they did. In this paper we develop simple models for the uncertainties in the collision strengths and decay rates for Fe XIII and apply them to the interpretation of density sensitive lines observed with the EUV Imagining spectrometer (EIS) on Hinode. We incorporate these uncertainties in a Bayesian framework. We consider both a pragmatic Bayesian method where the atomic physics information is unaffected by the observed data, and a fully Bayesian method where the data can be used to probe the physics. The former generally increases the uncertainty in the inferred density by about a factor of 5 compared with models that incorporate only statistical uncertainties. The latter reduces the uncertainties on the inferred densities, but identifies areas of possible systematic problems with either the atomic physics or the observed intensities.Comment: in press at Ap