Schwarzschild Tests of the Wahlquist-Estabrook-Buchman-Bardeen Tetrad Formulation for Numerical Relativity

A first order symmetric hyperbolic tetrad formulation of the Einstein equations developed by Estabrook and Wahlquist and put into a form suitable for numerical relativity by Buchman and Bardeen (the WEBB formulation) is adapted to explicit spherical symmetry and tested for accuracy and stability in the evolution of spherically symmetric black holes (the Schwarzschild geometry). The lapse and shift which specify the evolution of the coordinates relative to the tetrad congruence are reset at frequent time intervals to keep the constant-time hypersurfaces nearly orthogonal to the tetrad congruence and the spatial coordinate satisfying a kind of minimal rate of strain condition. By arranging through initial conditions that the constant-time hypersurfaces are asymptotically hyperbolic, we simplify the boundary value problem and improve stability of the evolution. Results are obtained for both tetrad gauges (``Nester'' and ``Lorentz'') of the WEBB formalism using finite difference numerical methods. We are able to obtain stable unconstrained evolution with the Nester gauge for certain initial conditions, but not with the Lorentz gauge.Comment: (accepted by Phys. Rev. D) minor changes; typos correcte

Adjusted ADM systems and their expected stability properties: constraint propagation analysis in Schwarzschild spacetime

In order to find a way to have a better formulation for numerical evolution of the Einstein equations, we study the propagation equations of the constraints based on the Arnowitt-Deser-Misner formulation. By adjusting constraint terms in the evolution equations, we try to construct an "asymptotically constrained system" which is expected to be robust against violation of the constraints, and to enable a long-term stable and accurate numerical simulation. We first provide useful expressions for analyzing constraint propagation in a general spacetime, then apply it to Schwarzschild spacetime. We search when and where the negative real or non-zero imaginary eigenvalues of the homogenized constraint propagation matrix appear, and how they depend on the choice of coordinate system and adjustments. Our analysis includes the proposal of Detweiler (1987), which is still the best one according to our conjecture but has a growing mode of error near the horizon. Some examples are snapshots of a maximally sliced Schwarzschild black hole. The predictions here may help the community to make further improvements.Comment: 23 pages, RevTeX4, many figures. Revised version. Added subtitle, reduced figures, rephrased introduction, and a native checked. :-

On free evolution of self gravitating, spherically symmetric waves

We perform a numerical free evolution of a selfgravitating, spherically symmetric scalar field satisfying the wave equation. The evolution equations can be written in a very simple form and are symmetric hyperbolic in Eddington-Finkelstein coordinates. The simplicity of the system allow to display and deal with the typical gauge instability present in these coordinates. The numerical evolution is performed with a standard method of lines fourth order in space and time. The time algorithm is Runge-Kutta while the space discrete derivative is symmetric (non-dissipative). The constraints are preserved under evolution (within numerical errors) and we are able to reproduce several known results.Comment: 15 pages, 15 figure

Exploiting gauge and constraint freedom in hyperbolic formulations of Einstein's equations

We present new many-parameter families of strongly and symmetric hyperbolic formulations of Einstein's equations that include quite general algebraic and live gauge conditions for the lapse. The first system that we present has 30 variables and incorporates an algebraic relationship between the lapse and the determinant of the three metric that generalizes the densitized lapse prescription. The second system has 34 variables and uses a family of live gauges that generalizes the Bona-Masso slicing conditions. These systems have free parameters even after imposing hyperbolicity and are expected to be useful in 3D numerical evolutions. We discuss under what conditions there are no superluminal characteristic speeds

Symmetric hyperbolic system in the Ashtekar formulation

We present a first-order symmetric hyperbolic system in the Ashtekar formulation of general relativity for vacuum spacetime. We add terms from constraint equations to the evolution equations with appropriate combinations, which is the same technique used by Iriondo, Leguizam\'on and Reula [Phys. Rev. Lett. 79, 4732 (1997)]. However our system is different from theirs in the points that we primarily use Hermiticity of a characteristic matrix of the system to characterize our system "symmetric", discuss the consistency of this system with reality condition, and show the characteristic speeds of the system.Comment: 4 pages, RevTeX, to appear in Phys. Rev. Lett., Comments added, refs update

Extending the lifetime of 3D black hole computations with a new hyperbolic system of evolution equations

We present a new many-parameter family of hyperbolic representations of Einstein's equations, which we obtain by a straightforward generalization of previously known systems. We solve the resulting evolution equations numerically for a Schwarzschild black hole in three spatial dimensions, and find that the stability of the simulation is strongly dependent on the form of the equations (i.e. the choice of parameters of the hyperbolic system), independent of the numerics. For an appropriate range of parameters we can evolve a single 3D black hole to $t \simeq 600 M$ -- $1300 M$, and are apparently limited by constraint-violating solutions of the evolution equations. We expect that our method should result in comparable times for evolutions of a binary black hole system.Comment: 11 pages, 2 figures, submitted to PR

Energy Norms and the Stability of the Einstein Evolution Equations

The Einstein evolution equations may be written in a variety of equivalent analytical forms, but numerical solutions of these different formulations display a wide range of growth rates for constraint violations. For symmetric hyperbolic formulations of the equations, an exact expression for the growth rate is derived using an energy norm. This expression agrees with the growth rate determined by numerical solution of the equations. An approximate method for estimating the growth rate is also derived. This estimate can be evaluated algebraically from the initial data, and is shown to exhibit qualitatively the same dependence as the numerically-determined rate on the parameters that specify the formulation of the equations. This simple rate estimate therefore provides a useful tool for finding the most well-behaved forms of the evolution equations.Comment: Corrected typos; to appear in Physical Review

CWRML: representing crop wild relative conservation and use data in XML

Background Crop wild relatives are wild species that are closely related to crops. They are valuable as potential gene donors for crop improvement and may help to ensure food security for the future. However, they are becoming increasingly threatened in the wild and are inadequately conserved, both in situ and ex situ. Information about the conservation status and utilisation potential of crop wild relatives is diverse and dispersed, and no single agreed standard exists for representing such information; yet, this information is vital to ensure these species are effectively conserved and utilised. The European Community-funded project, European Crop Wild Relative Diversity Assessment and Conservation Forum, determined the minimum information requirements for the conservation and utilisation of crop wild relatives and created the Crop Wild Relative Information System, incorporating an eXtensible Markup Language (XML) schema to aid data sharing and exchange. Results Crop Wild Relative Markup Language (CWRML) was developed to represent the data necessary for crop wild relative conservation and ensure that they can be effectively utilised for crop improvement. The schema partitions data into taxon-, site-, and population-specific elements, to allow for integration with other more general conservation biology schemata which may emerge as accepted standards in the future. These elements are composed of sub-elements, which are structured in order to facilitate the use of the schema in a variety of crop wild relative conservation and use contexts. Pre-existing standards for data representation in conservation biology were reviewed and incorporated into the schema as restrictions on element data contents, where appropriate. Conclusion CWRML provides a flexible data communication format for representing in situ and ex situ conservation status of individual taxa as well as their utilisation potential. The development of the schema highlights a number of instances where additional standards-development may be valuable, particularly with regard to the representation of population-specific data and utilisation potential. As crop wild relatives are intrinsically no different to other wild plant species there is potential for the inclusion of CWRML data elements in the emerging standards for representation of biodiversity data

Symmetric Hyperbolic System in the Self-dual Teleparallel Gravity

In order to discuss the well-posed initial value formulation of the teleparallel gravity and apply it to numerical relativity a symmetric hyperbolic system in the self-dual teleparallel gravity which is equivalent to the Ashtekar formulation is posed. This system is different from the ones in other works by that the reality condition of the spatial metric is included in the symmetric hyperbolicity and then is no longer an independent condition. In addition the constraint equations of this system are rather simpler than the ones in other works.Comment: 8 pages, no figure

The Cauchy Problem for the Einstein Equations

Various aspects of the Cauchy problem for the Einstein equations are surveyed, with the emphasis on local solutions of the evolution equations. Particular attention is payed to giving a clear explanation of conceptual issues which arise in this context. The question of producing reduced systems of equations which are hyperbolic is examined in detail and some new results on that subject are presented. Relevant background from the theory of partial differential equations is also explained at some lengthComment: 98 page