An Einstein-Bianchi system for Smooth Lattice General Relativity. II. 3+1 vacuum spacetimes

We will present a complete set of equations, in the form of an Einstein-Bianchi system, that describe the evolution of generic smooth lattices in spacetime. All 20 independent Riemann curvatures will be evolved in parallel with the leg-lengths of the lattice. We will show that the evolution equations for the curvatures forms a hyperbolic system and that the associated constraints are preserved. This work is a generalisation of our previous paper arXiv:1101.3171 on the Einstein-Bianchi system for the Schwarzschild spacetime to general 3+1 vacuum spacetimes

An ADM 3+1 formulation for Smooth Lattice General Relativity

A new hybrid scheme for numerical relativity will be presented. The scheme will employ a 3-dimensional spacelike lattice to record the 3-metric while using the standard 3+1 ADM equations to evolve the lattice. Each time step will involve three basic steps. First, the coordinate quantities such as the Riemann and extrinsic curvatures are extracted from the lattice. Second, the 3+1 ADM equations are used to evolve the coordinate data, and finally, the coordinate data is used to update the scalar data on the lattice (such as the leg lengths). The scheme will be presented only for the case of vacuum spacetime though there is no reason why it could not be extended to non-vacuum spacetimes. The scheme allows any choice for the lapse function and shift vectors. An example for the Kasner $T^3$ cosmology will be presented and it will be shown that the method has, for this simple example, zero discretisation error.Comment: 18 pages, plain TeX, 5 epsf figues, gzipped ps file also available at http://newton.maths.monash.edu.au:8000/preprints/3+1-slgr.ps.g

The need for psychiatric treatment in the general population: the Camberwell Needs for Care Survey

Background, This paper presents the first results of a two-stage psychiatric population survey, which uses a new method of directly evaluating needs for specific psychiatric treatment and the extent to which they have been met.Method, The sample was drawn at random from the population of an area of inner south London with high levels of deprivation. Seven hundred and sixty subjects aged 18-65 completed the GHQ-28. All those scoring > 5 and half of the rest were invited to take part in the second stage, comprising measures of mental state (SCAN), social role performance (SRPS), life events and difficulties (LEDS) and a Treatment Inventory. This information was used to rate the community version of the Needs for Care Assessment (NFCAS-C).Results, In all, 408 subjects were interviewed in the second stage. The weighted 1 month prevalence of hierarchically ordered ICD-10 psychiatric disorders was 9.8 %, the 1 year prevalence 12.3 %. The equivalent prevalences for depressive episode were 3.1 % and 5.3 % respectively, while those for anxiety states were both 2.8 %. At interview nearly 10% of the population were identified as having a need for the treatment of a psychiatric condition. This rose to 10.4 % if the whole of the preceding year was assessed. Less than half of all potentially meetable needs were met. There was only partial overlap between diagnosis and an adjudged need for treatment.Conclusion. A majority of people with mental health problems do not have proper treatment; given more resources and greater public and medical awareness, most could be treated by family doctors

Is the Regge Calculus a consistent approximation to General Relativity?

We will ask the question of whether or not the Regge calculus (and two related simplicial formulations) is a consistent approximation to General Relativity. Our criteria will be based on the behaviour of residual errors in the discrete equations when evaluated on solutions of the Einstein equations. We will show that for generic simplicial lattices the residual errors can not be used to distinguish metrics which are solutions of Einstein's equations from those that are not. We will conclude that either the Regge calculus is an inconsistent approximation to General Relativity or that it is incorrect to use residual errors in the discrete equations as a criteria to judge the discrete equations.Comment: 27 pages, plain TeX, very belated update to match journal articl

Long term stable integration of a maximally sliced Schwarzschild black hole using a smooth lattice method

We will present results of a numerical integration of a maximally sliced Schwarzschild black hole using a smooth lattice method. The results show no signs of any instability forming during the evolutions to t=1000m. The principle features of our method are i) the use of a lattice to record the geometry, ii) the use of local Riemann normal coordinates to apply the 1+1 ADM equations to the lattice and iii) the use of the Bianchi identities to assist in the computation of the curvatures. No other special techniques are used. The evolution is unconstrained and the ADM equations are used in their standard form.Comment: 47 pages including 26 figures, plain TeX, also available at http://www.maths.monash.edu.au/~leo/preprint

Fast algorithms for computing defects and their derivatives in the Regge calculus

Any practical attempt to solve the Regge equations, these being a large system of non-linear algebraic equations, will almost certainly employ a Newton-Raphson like scheme. In such cases it is essential that efficient algorithms be used when computing the defect angles and their derivatives with respect to the leg-lengths. The purpose of this paper is to present details of such an algorithm.Comment: 38 pages, 10 figure

A simple expression for the ADM mass

We show by an almost elementary calculation that the ADM mass of an asymptotically flat space can be computed as a limit involving a rate of change of area of a closed 2-surface. The result is essentially the same as that given by Brown and York. We will prove this result in two ways, first by direct calculation from the original formula as given by Arnowitt, Deser and Misner and second as a corollary of an earlier result by Brewin for the case of simplicial spaces.Comment: 9 pages, 1 figur

Regge Calculus as a Fourth Order Method in Numerical Relativity

The convergence properties of numerical Regge calculus as an approximation to continuum vacuum General Relativity is studied, both analytically and numerically. The Regge equations are evaluated on continuum spacetimes by assigning squared geodesic distances in the continuum manifold to the squared edge lengths in the simplicial manifold. It is found analytically that, individually, the Regge equations converge to zero as the second power of the lattice spacing, but that an average over local Regge equations converges to zero as (at the very least) the third power of the lattice spacing. Numerical studies using analytic solutions to the Einstein equations show that these averages actually converge to zero as the fourth power of the lattice spacing.Comment: 14 pages, LaTeX, 8 figures mailed in separate file or email author directl

Slice Stretching Effects for Maximal Slicing of a Schwarzschild Black Hole

Slice stretching effects such as slice sucking and slice wrapping arise when foliating the extended Schwarzschild spacetime with maximal slices. For arbitrary spatial coordinates these effects can be quantified in the context of boundary conditions where the lapse arises as a linear combination of odd and even lapse. Favorable boundary conditions are then derived which make the overall slice stretching occur late in numerical simulations. Allowing the lapse to become negative, this requirement leads to lapse functions which approach at late times the odd lapse corresponding to the static Schwarzschild metric. Demanding in addition that a numerically favorable lapse remains non-negative, as result the average of odd and even lapse is obtained. At late times the lapse with zero gradient at the puncture arising for the puncture evolution is precisely of this form. Finally, analytic arguments are given on how slice stretching effects can be avoided. Here the excision technique and the working mechanism of the shift function are studied in detail.Comment: 16 pages, 4 figures, revised version including a study on how slice stretching can be avoided by using excision and/or shift

A Smooth Lattice construction of the Oppenheimer-Snyder spacetime

We present test results for the smooth lattice method using an Oppenheimer-Snyder spacetime. The results are in excellent agreement with theory and numerical results from other authors.Comment: 60 pages, 28 figure