Initial Data for Numerical Relativity

Initial data are the starting point for any numerical simulation. In the case of numerical relativity, Einstein's equations constrain our choices of these initial data. We will examine several of the formalisms used for specifying Cauchy initial data in the 3+1 decomposition of Einstein's equations. We will then explore how these formalisms have been used in constructing initial data for spacetimes containing black holes and neutron stars. In the topics discussed, emphasis is placed on those issues that are important for obtaining astrophysically realistic initial data for compact binary coalescence.Comment: 50 pages, LaTeX(livrev.cls), Review article for "Living Reviews in Relativity" (http://www.livingreviews.org/), July 200

Systematics and distributions of the genera Cyrtusa Erichson, Ecarinosphaerula Hatch, Isoplastus Horn, Liocyrtusa Daffner, Lionothus Brown, and Zeadolopus Broun of the United States and Canada (Coleoptera: Leiodidae: Leiodinae: Leiodini)

The following genera of Leiodini (Coleoptera: Leiodidae: Leiodinae) of the continental United States and Canada are reviewed: Cyrtusa Erichson, with two species; Isoplastus Horn, with two species (one new); Liocyrtusa Daffner, with three species; Lionothus Brown, with five species (three new), and Zeadolopus Broun, with four species (all genera are in the “Cyrtusa genus group”) and Ecarinosphaerula Hatch, with one named species (in the “Leiodes genus group”). The new species are Isoplastus floridanus Peck and Cook of Florida; Lionothus bidentatus Peck and Cook of Texas and Oklahoma, Lionothus exiguus Peck and Cook of Florida and Texas, and Lionothus parvoculus Peck and Cook of Arizona and New Mexico. Bionomic data on the species are given, and complete known distributions are mapped

Testing a Simplified Version of Einstein's Equations for Numerical Relativity

Solving dynamical problems in general relativity requires the full machinery of numerical relativity. Wilson has proposed a simpler but approximate scheme for systems near equilibrium, like binary neutron stars. We test the scheme on isolated, rapidly rotating, relativistic stars. Since these objects are in equilibrium, it is crucial that the approximation work well if we are to believe its predictions for more complicated systems like binaries. Our results are very encouraging.Comment: 9 pages (RevTeX 3.0 with 6 uuencoded figures), CRSR-107

Substructure Discovery Using Minimum Description Length and Background Knowledge

The ability to identify interesting and repetitive substructures is an essential component to discovering knowledge in structural data. We describe a new version of our SUBDUE substructure discovery system based on the minimum description length principle. The SUBDUE system discovers substructures that compress the original data and represent structural concepts in the data. By replacing previously-discovered substructures in the data, multiple passes of SUBDUE produce a hierarchical description of the structural regularities in the data. SUBDUE uses a computationally-bounded inexact graph match that identifies similar, but not identical, instances of a substructure and finds an approximate measure of closeness of two substructures when under computational constraints. In addition to the minimum description length principle, other background knowledge can be used by SUBDUE to guide the search towards more appropriate substructures. Experiments in a variety of domains demonstrate SUBDUE's ability to find substructures capable of compressing the original data and to discover structural concepts important to the domain. Description of Online Appendix: This is a compressed tar file containing the SUBDUE discovery system, written in C. The program accepts as input databases represented in graph form, and will output discovered substructures with their corresponding value.Comment: See http://www.jair.org/ for an online appendix and other files accompanying this articl

Towards a wave-extraction method for numerical relativity. V. Extracting the Weyl scalars in the quasi-Kinnersley tetrad from spatial data

We extract the Weyl scalars $\Psi_0$ and $\Psi_4$ in the quasi-Kinnersley tetrad by finding initially the (gauge--, tetrad--, and background--independent) transverse quasi-Kinnersley frame. This step still leaves two undetermined degrees of freedom: the ratio $|\Psi_0|/|\Psi_4|$, and one of the phases (the product $|\Psi_0|\cdot |\Psi_4|$ and the {\em sum} of the phases are determined by the so-called BB radiation scalar). The residual symmetry ("spin/boost") can be removed by gauge fixing of spin coefficients in two steps: First, we break the boost symmetry by requiring that $\rho$ corresponds to a global constant mass parameter that equals the ADM mass (or, equivalently in perturbation theory, that $\rho$ or $\mu$ equal their values in the no-radiation limits), thus determining the two moduli of the Weyl scalars $|\Psi_0|, |\Psi_4|$, while leaving their phases as yet undetermined. Second, we break the spin symmetry by requiring that the ratio $\pi/\tau$ gives the expected polarization state for the gravitational waves, thus determining the phases. Our method of gauge fixing--specifically its second step--is appropriate for cases for which the Weyl curvature is purely electric. Applying this method to Misner and Brill--Lindquist data, we explicitly find the Weyl scalars $\Psi_0$ and $\Psi_4$ perturbatively in the quasi-Kinnersley tetrad.Comment: 13 page

Quasi-circular Orbits for Spinning Binary Black Holes

Using an effective potential method we examine binary black holes where the individual holes carry spin. We trace out sequences of quasi-circular orbits and locate the innermost stable circular orbit as a function of spin. At large separations, the sequences of quasi-circular orbits match well with post-Newtonian expansions, although a clear signature of the simplifying assumption of conformal flatness is seen. The position of the ISCO is found to be strongly dependent on the magnitude of the spin on each black hole. At close separations of the holes, the effective potential method breaks down. In all cases where an ISCO could be determined, we found that an apparent horizon encompassing both holes forms for separations well inside the ISCO. Nevertheless, we argue that the formation of a common horizon is still associated with the breakdown of the effective potential method.Comment: 13 pages, 10 figures, submitted to PR

Summer Cook Associate Professor of Kinesiology (COLA) travels to Australia

I was partially supported by a CIE International development grant to travel to Perth, Australia for one month in 2016. Dr. Timothy Fairchild, a colleague and friend, invited me to visit the Department of Psychology and Exercise Science at Murdoch University in the hopes of establishing a research relationship with the University of New Hampshire Department of Kinesiology. Over the last several years, I have had students earn Summer Undergraduate Research Fellowships (SURF) from the Hamel Center and have always wanted to give students an opportunity to apply for SURF abroad grants or to study abroad while conducting research within our field. The idea of international collaboration and the opportunity to leave the New England winter to work on the coast of Western Australia made the decision to travel very simple