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

Quasiequilibrium sequences of black-hole--neutron-star binaries in general relativity

We construct quasiequilibrium sequences of black hole-neutron star binaries for arbitrary mass ratios by solving the constraint equations of general relativity in the conformal thin-sandwich decomposition. We model the neutron star as a stationary polytrope satisfying the relativistic equations of hydrodynamics, and account for the black hole by imposing equilibrium boundary conditions on the surface of an excised sphere (the apparent horizon). In this paper we focus on irrotational configurations, meaning that both the neutron star and the black hole are approximately nonspinning in an inertial frame. We present results for a binary with polytropic index n=1, mass ratio M_{irr}^{BH}/M_{B}^{NS}=5 and neutron star compaction M_{ADM,0}^{NS}/R_0=0.0879, where M_{irr}^{BH} is the irreducible mass of the black hole, M_{B}^{NS} the neutron star baryon rest-mass, and M_{ADM,0}^{NS} and R_0 the neutron star Arnowitt-Deser-Misner mass and areal radius in isolation, respectively. Our models represent valid solutions to Einstein's constraint equations and may therefore be employed as initial data for dynamical simulations of black hole-neutron star binaries.Comment: 5 pages, 1 figure, revtex4, published in Phys.Rev.

Quasiequilibrium black hole-neutron star binaries in general relativity

We construct quasiequilibrium sequences of black hole-neutron star binaries in general relativity. We solve Einstein's constraint equations in the conformal thin-sandwich formalism, subject to black hole boundary conditions imposed on the surface of an excised sphere, together with the relativistic equations of hydrostatic equilibrium. In contrast to our previous calculations we adopt a flat spatial background geometry and do not assume extreme mass ratios. We adopt a Gamma=2 polytropic equation of state and focus on irrotational neutron star configurations as well as approximately nonspinning black holes. We present numerical results for ratios of the black hole's irreducible mass to the neutron star's ADM mass in isolation of M_{irr}^{BH}/M_{ADM,0}^{NS} = 1, 2, 3, 5, and 10. We consider neutron stars of baryon rest mass M_B^{NS}/M_B^{max} = 83% and 56%, where M_B^{max} is the maximum allowed rest mass of a spherical star in isolation for our equation of state. For these sequences, we locate the onset of tidal disruption and, in cases with sufficiently large mass ratios and neutron star compactions, the innermost stable circular orbit. We compare with previous results for black hole-neutron star binaries and find excellent agreement with third-order post-Newtonian results, especially for large binary separations. We also use our results to estimate the energy spectrum of the outgoing gravitational radiation emitted during the inspiral phase for these binaries.Comment: 17 pages, 15 figures, published in Phys. Rev.

Deep learning cardiac motion analysis for human survival prediction

Motion analysis is used in computer vision to understand the behaviour of moving objects in sequences of images. Optimising the interpretation of dynamic biological systems requires accurate and precise motion tracking as well as efficient representations of high-dimensional motion trajectories so that these can be used for prediction tasks. Here we use image sequences of the heart, acquired using cardiac magnetic resonance imaging, to create time-resolved three-dimensional segmentations using a fully convolutional network trained on anatomical shape priors. This dense motion model formed the input to a supervised denoising autoencoder (4Dsurvival), which is a hybrid network consisting of an autoencoder that learns a task-specific latent code representation trained on observed outcome data, yielding a latent representation optimised for survival prediction. To handle right-censored survival outcomes, our network used a Cox partial likelihood loss function. In a study of 302 patients the predictive accuracy (quantified by Harrell's C-index) was significantly higher (p < .0001) for our model C=0.73 (95$\%$ CI: 0.68 - 0.78) than the human benchmark of C=0.59 (95$\%$ CI: 0.53 - 0.65). This work demonstrates how a complex computer vision task using high-dimensional medical image data can efficiently predict human survival

Solving Einstein's equations for rotating spacetimes: Evolution of relativistic star clusters

A numerical relativity code designed to evolve rotating axisymmetric spacetimes is constructed. Both polarization states of gravitational radiation can be tracked. The source of the gravitational field is chosen to be a configuration of collisionless particles. The code is used to evaluate the stability of polytropic and toroidal star clusters. The formation of Kerr black holes by the collapse of unstable clusters is demonstrated. Unstable clusters with J/(M^2) 1 collapse to new equilibrium configurations

Treating instabilities in a hyperbolic formulation of Einstein's equations

We have recently constructed a numerical code that evolves a spherically symmetric spacetime using a hyperbolic formulation of Einstein's equations. For the case of a Schwarzschild black hole, this code works well at early times, but quickly becomes inaccurate on a time scale of 10-100 M, where M is the mass of the hole. We present an analytic method that facilitates the detection of instabilities. Using this method, we identify a term in the evolution equations that leads to a rapidly-growing mode in the solution. After eliminating this term from the evolution equations by means of algebraic constraints, we can achieve free evolution for times exceeding 10000M. We discuss the implications for three-dimensional simulations.Comment: 13 pages, 9 figures. To appear in Phys. Rev.

Implementing an apparent-horizon finder in three dimensions

Locating apparent horizons is not only important for a complete understanding of numerically generated spacetimes, but it may also be a crucial component of the technique for evolving black-hole spacetimes accurately. A scheme proposed by Libson et al., based on expanding the location of the apparent horizon in terms of symmetric trace-free tensors, seems very promising for use with three-dimensional numerical data sets. In this paper, we generalize this scheme and perform a number of code tests to fully calibrate its behavior in black-hole spacetimes similar to those we expect to encounter in solving the binary black-hole coalescence problem. An important aspect of the generalization is that we can compute the symmetric trace-free tensor expansion to any order. This enables us to determine how far we must carry the expansion to achieve results of a desired accuracy. To accomplish this generalization, we describe a new and very convenient set of recurrence relations which apply to symmetric trace-free tensors.Comment: 14 pages (RevTeX 3.0 with 3 figures

Numerical Evolution of Black Holes with a Hyperbolic Formulation of General Relativity

We describe a numerical code that solves Einstein's equations for a Schwarzschild black hole in spherical symmetry, using a hyperbolic formulation introduced by Choquet-Bruhat and York. This is the first time this formulation has been used to evolve a numerical spacetime containing a black hole. We excise the hole from the computational grid in order to avoid the central singularity. We describe in detail a causal differencing method that should allow one to stably evolve a hyperbolic system of equations in three spatial dimensions with an arbitrary shift vector, to second-order accuracy in both space and time. We demonstrate the success of this method in the spherically symmetric case.Comment: 23 pages RevTeX plus 7 PostScript figures. Submitted to Phys. Rev.

Shoot growth of woody trees and shrubs is predicted by maximum plant height and associated traits

1. The rate of elongation and thickening of individual branches (shoots) varies across plant species. This variation is important for the outcome of competition and other plant-plant interactions. Here we compared rates of shoot growth across 44 species from tropical, warm temperate, and cool temperate forests of eastern Australia.2. Shoot growth rate was found to correlate with a suite of traits including the potential height of the species, xylem-specific conductivity, leaf size, leaf area per xylem cross-section, twig diameter (at 40 cm length), wood density and modulus of elasticity.3. Within this suite of traits, maximum plant height was the clearest correlate of growth rates, explaining 50 to 67% of the variation in growth overall (p p 4. Growth rates were not strongly correlated with leaf nitrogen or leaf mass per unit leaf area.5. Correlations between growth and maximum height arose both across latitude (47%, p p p p < 0.0001), reflecting intrinsic differences across species and sites