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

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

Alternatives to standard puncture initial data for binary black hole evolution

Standard puncture initial data have been widely used for numerical binary black hole evolutions despite their shortcomings, most notably the inherent lack of gravitational radiation at the initial time that is later followed by a burst of spurious radiation. We study the evolution of three alternative initial data schemes. Two of the three alternatives are based on post-Newtonian expansions that contain realistic gravitational waves. The first scheme is based on a second-order post-Newtonian expansion in Arnowitt, Deser, and Misner transverse-traceless (ADMTT) gauge that has been resummed to approach standard puncture data at the black holes. The second scheme is based on asymptotic matching of the 4-metrics of two tidally perturbed Schwarzschild solutions to a first-order post-Newtonian expansion in ADMTT gauge away from the black holes. The final alternative is obtained through asymptotic matching of the 4-metrics of two tidally perturbed Schwarzschild solutions to a second-order post-Newtonian expansion in harmonic gauge away from the black holes. When evolved, the second scheme fails to produce quasicircular orbits (and instead leads to a nearly head-on collision). This failure can be traced back to inaccuracies in the extrinsic curvature due to low order matching. More encouraging is that the latter two alternatives lead to quasicircular orbits and show gravitational radiation from the onset of the evolution, as well as a reduction of spurious radiation. Current deficiencies compared to standard punctures data include more eccentric trajectories during the inspiral and larger constraint violations, since the alternative data sets are only approximate solutions of Einstein's equations. The eccentricity problem can be ameliorated by adjusting the initial momentum parameters.Comment: 11 pages, 11 figures, 1 appendix, typos corrected, removed duplicate reference, matches published versio

Evolution of Binary Black Hole Spacetimes

We describe early success in the evolution of binary black hole spacetimes with a numerical code based on a generalization of harmonic coordinates. Indications are that with sufficient resolution this scheme is capable of evolving binary systems for enough time to extract information about the orbit, merger and gravitational waves emitted during the event. As an example we show results from the evolution of a binary composed of two equal mass, non-spinning black holes, through a single plunge-orbit, merger and ring down. The resultant black hole is estimated to be a Kerr black hole with angular momentum parameter a~0.70. At present, lack of resolution far from the binary prevents an accurate estimate of the energy emitted, though a rough calculation suggests on the order of 5% of the initial rest mass of the system is radiated as gravitational waves during the final orbit and ringdown.Comment: 4 pages, 3 figure

Inner boundary conditions for black hole Initial Data derived from Isolated Horizons

We present a set of boundary conditions for solving the elliptic equations in the Initial Data Problem for space-times containing a black hole, together with a number of constraints to be satisfied by the otherwise freely specifiable standard parameters of the Conformal Thin Sandwich formulation. These conditions altogether are sufficient for the construction of a horizon that is instantaneously in equilibrium in the sense of the Isolated Horizons formalism. We then investigate the application of these conditions to the Initial Data Problem of binary black holes and discuss the relation of our analysis with other proposals that exist in the literature.Comment: 13 pages. Major general revision. Section V comparing with previous approaches restructured; discussion on the lapse boundary condition extended. Appendix with some technical details added. Version accepted for publication in Phys.Rev.

Design and development of a water vapor electrolysis unit

Design and development of water vapor electrolysis unit for oxygen productio

Gravitational waves from nonspinning black hole-neutron star binaries: dependence on equations of state

We report results of a numerical-relativity simulation for the merger of a black hole-neutron star binary with a variety of equations of state (EOSs) modeled by piecewise polytropes. We focus in particular on the dependence of the gravitational waveform at the merger stage on the EOSs. The initial conditions are computed in the moving-puncture framework, assuming that the black hole is nonspinning and the neutron star has an irrotational velocity field. For a small mass ratio of the binaries (e.g., MBH/MNS = 2 where MBH and MNS are the masses of the black hole and neutron star, respectively), the neutron star is tidally disrupted before it is swallowed by the black hole irrespective of the EOS. Especially for less-compact neutron stars, the tidal disruption occurs at a more distant orbit. The tidal disruption is reflected in a cutoff frequency of the gravitational-wave spectrum, above which the spectrum amplitude exponentially decreases. A clear relation is found between the cutoff frequency of the gravitational-wave spectrum and the compactness of the neutron star. This relation also depends weakly on the stiffness of the EOS in the core region of the neutron star, suggesting that not only the compactness but also the EOS at high density is reflected in gravitational waveforms. The mass of the disk formed after the merger shows a similar correlation with the EOS, whereas the spin of the remnant black hole depends primarily on the mass ratio of the binary, and only weakly on the EOS. Properties of the remnant disks are also analyzed.Comment: 27pages, 21 figures; erratum is added on Aug 5. 201

Predictions for the fracture toughness of cancellous bone of fracture neck of femur patients

Current protocol in determining if a patient is osteoporotic and their fracture risk is based on dual energy X-ray absorptiometry (DXA). DXA gives an indication of their bone mineral density (BMD) which is the product of both the porosity and density of the mineralized bone tissue; this is usually taken at the hip. The DXA results are assessed using the fracture risk assessment tool as recommended by the World Health Organization. While this provides valuable data on a person’s fracture risk advancements in medical imagining technology enables development of more robust and accurate risk assessment tools. In order to develop such tools in vitro analysis of bone is required to assess the morphological properties of bone osteoporotic bone tissue and how these pertain to the fracture toughness (Kcmax) of the tissue.Support was provided by the EPSRC (EP/K020196: Point-ofCare High Accuracy Fracture Risk Prediction), the UK Department of Transport under the BOSCOS (Bone Scanning for Occupant Safety) project, and approved by Gloucester and Cheltenham NHS Trust hospitals under ethical consent (BOSCOS – Mr. Curwen CI REC ref 01/179G)