Source author record

David Hilditch

David Hilditch appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

29works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

29 published item(s)

preprint2022arXiv

Gauge structure of the Einstein field equations in Bondi-like coordinates

The characteristic initial (boundary) value problem has numerous applications in general relativity (GR) involving numerical studies, and is often formulated using Bondi-like coordinates. Recently it was shown that several prototype formulations of this type are only weakly hyperbolic. Presently we examine the root cause of this result. In a linear analysis we identify the gauge, constraint and physical blocks in the principal part of the Einstein field equations in such a gauge, and show that the subsystem related to the gauge variables is only weakly hyperbolic. Weak hyperbolicity of the full system follows as a consequence in many cases. We demonstrate this explicitly in specific examples, and thus argue that Bondi-like gauges result in weakly hyperbolic free evolution systems under quite general conditions. Consequently the characteristic initial (boundary) value problem of GR in these gauges is rendered ill-posed in the simplest norms one would like to employ. The possibility of finding good alternative norms, in which well-posedness is achieved, is discussed. So motivated, we present numerical convergence tests with an implementation of full GR which demonstrate the effect of weak hyperbolicity in practice.

preprint2022arXiv

Quasinormal modes of Proca fields in a Schwarzschild-AdS spacetime

We present new results concerning the Proca massive vector field in a Schwarzschild-AdS black hole geometry. We provide a first principles analysis of Proca vector fields in this geometry using both the vector spherical harmonic (VSH) separation method and the Frolov-Krtouš-Kubizňák-Santos (FKKS) method that separates the relevant equations in spinning geometries. The analysis in the VSH method shows, on one hand, that it is arduous to separate the scalar-type from the vector-type polarizations of the electric sector of the Proca field, and on the other hand, it displays clearly the electric and the magnetic mode sectors. The analysis in the FKKS method is performed by taking the nonrotating limit of the Kerr-AdS spacetime, and shows that the ansatz decouples the polarizations in the electric mode sector even in the nonrotating limit. On the other hand, it captures only two of the three possible polarizations, the magnetic mode sector is missing. The reason for the absence of this polarization is related to the degeneracy of the principal tensor in static spherical symmetric spacetimes. The degrees of freedom and quasinormal modes in both separation methods of the Proca field are found. The frequencies of the quasinormal modes are also computed. For the electric mode sector in the VSH method the frequencies are found through an extension, which substitutes number coefficients by matrix coefficients, of the Horowitz-Hubeny numerical procedure, whereas for the magnetic mode sector in the VSH method and the electric sector of the FKKS method it is shown that a direct use of the procedure can be made. The values of the quasinormal mode frequencies obtained for each method are compared and showed to be in good agreement with each other. This further supports the analytical approaches presented here for the behavior of the Proca field in a Schwarzschild-AdS black hole background.

preprint2022arXiv

Regularizing Dual-Frame Generalized Harmonic Gauge at Null Infinity

The dual-frame formalism leads to an approach to extend numerical relativity simulations in generalized harmonic gauge (GHG) all the way to null infinity. A major setback is that without care, even simple choices of initial data give rise to logarithmically divergent terms that would result in irregular variables and equations on the compactified domain, which would in turn prevent accurate numerical approximation. It has been shown, however, that a suitable choice of gauge and constraint addition can be used to prevent their appearance. Presently we give a first order symmetric hyperbolic reduction of general relativity in GHG on compactified hyperboloidal slices that exploits this knowledge and eradicates these log-terms at leading orders. Because of their effect on the asymptotic solution space, specific formally singular terms are systematically chosen to remain. Such formally singular terms have been successfully treated numerically in toy models and result in a formulation with the desirable property that unphysical radiation content near infinity is suppressed.

preprint2022arXiv

Shock-avoiding slicing conditions: tests and calibrations

While the 1+log slicing condition has been extremely successful in numerous numerical relativity simulations, it is also known to develop "gauge-shocks" in some examples. Alternative "shock-avoiding" slicing conditions suggested by Alcubierre prevent these pathologies in those examples, but have not yet been explored and tested very broadly. In this paper we compare the performance of shock-avoiding slicing conditions with those of 1+log slicing for a number of "text-book" problems, including black holes and relativistic stars. While, in some simulations, the shock-avoiding slicing conditions feature some unusual properties and lead to more "gauge-dynamics" than the 1+log slicing condition, we find that they perform quite similarly in terms of stability and accuracy, and hence provide a very viable alternative to 1+log slicing.

preprint2020arXiv

A semi-linear wave model for critical collapse

In spherical symmetry compelling numerical evidence suggests that in general relativity solutions near the threshold of black hole formation exhibit critical behavior. One aspect of this is that threshold solutions themselves are self-similar and are, in a certain sense, unique. To an extent yet to be fully understood, the same phenomena persist beyond spherical symmetry. It is therefore desirable to construct models that exhibit such symmetry at the threshold of blow-up. Starting with deformations of the wave equation, we discuss models which have discretely self-similar threshold solutions. We study threshold solutions in the past light cone of the blow-up point. In spherical symmetry there is a sense in which a unique critical solution exists. Spherical numerical evolutions are also presented for more general models, and exhibit similar behavior. Away from spherical symmetry threshold solutions attain more freedom. Different topologies of blow-up are possible, and even locally the critical solution needs reinterpretation as a parameterized family.

preprint2020arXiv

Analytical and numerical treatment of perturbed black holes in horizon-penetrating coordinates

The deviations of non-linear perturbations of black holes from the linear case are important in the context of ringdown signals with large signal-to-noise ratio. To facilitate a comparison between the two we derive several results of linear perturbation theory in coordinates which may be adopted in numerical work. Specifically, our results are derived in Kerr-Schild coordinates adjusted by a general height function. In the first part of the paper we address the questions: for an initial configuration of a massless scalar field, what is the amplitude of the excited quasinormal mode (QNM) for any observer outside outside the event horizon, and furthermore what is the resulting tail contribution? This is done by constructing the full Green's function for the problem with exact solutions of the confluent Heun equation satisfying appropriate boundary conditions. In the second part of the paper, we detail new developments to our pseudospectral numerical relativity code bamps to handle scalar fields. In the linear regime we employ precisely the Kerr-Schild coordinates treated by our previous analysis. In particular, we evolve pure QNM type initial data along with several other types of initial data and report on the presence of overtone modes in the signal.

preprint2020arXiv

Hyperbolicity of Divergence Cleaning and Vector Potential Formulations of GRMHD

We examine hyperbolicity of general relativistic magnetohydrodynamics with divergence cleaning, a flux-balance law form of the model not covered by our earlier analysis. The calculations rely again on a dual-frame approach, which allows us to effectively exploit the structure present in the principal part. We find, in contrast to the standard flux-balance law form of the equations, that this formulation is strongly hyperbolic, and thus admits a well-posed initial value problem. Formulations involving the vector potential as an evolved quantity are then considered. Carefully reducing to first order, we find that such formulations can also be made strongly hyperbolic. Despite the unwieldy form of the characteristic variables we therefore conclude that of the free-evolution formulations of general relativistic magnetohydrodynamics presently used in numerical relativity, the divergence cleaning and vector potential formulations are preferred.

preprint2020arXiv

Improved existence for the characteristic initial value problem with the conformal Einstein field equations

We adapt Luk's analysis of the characteristic initial value problem in General Relativity to the asymptotic characteristic problem for the conformal Einstein field equations to demonstrate the local existence of solutions in a neighbourhood of the set on which the data are given. In particular, we obtain existence of solutions along a narrow rectangle along null infinity which, in turn, corresponds to an infinite domain in the asymptotic region of the physical spacetime. This result generalises work by Kánnár on the local existence of solutions to the characteristic initial value problem by means of Rendall's reduction strategy. In analysing the conformal Einstein equations we make use of the Newman-Penrose formalism and a gauge due to J. Stewart.

preprint2019arXiv

Conformally flat slices of asymptotically flat spacetimes

For mathematical convenience initial data sets in numerical relativity are often taken to be conformally flat. Employing the dual-foliation formalism, we investigate the physical consequences of this assumption. Working within a large class of asymptotically flat spacetimes we show that the ADM linear momentum is governed by the leading Lorentz part of a boost even in the presence of supertranslation-like terms. Following up, we find that in spacetimes that are asymptotically flat, and admit spatial slices with vanishing linear momentum that are sufficiently close to conformal flatness, any boosted slice can not be conformally flat. Consequently there are no conformally flat boosted slices of the Schwarzschild spacetime. This confirms the previously anticipated explanation for the presence of junk-radiation in Brandt-Bruegmann puncture data.

preprint2019arXiv

The Hyperboloidal Numerical Evolution of a Good-Bad-Ugly Wave Equation

One method for the numerical treatment of future null-infinity is to decouple coordinates from the tensor basis and choose each in a careful manner. This dual-frame approach is hampered by logarithmically divergent terms that appear in a naive choice of evolved variables. Here we consider a system of wave equations that satisfy the weak-null condition and serve as a model system with similar nonlinearities to those present in the Einstein field equations in generalized harmonic gauge. We show that these equations can be explicitly regularized by a nonlinear change of variables. Working in spherical symmetry, a numerical implementation of this model using compactified hyperboloidal slices is then presented. Clean convergence is found for the regularized system. Although more complicated, it is expected that general relativity can be treated similarly.

preprint2018arXiv

The Weak Null Condition in Free-evolution Schemes for Numerical Relativity: Dual Foliation GHG with Constraint Damping

All strategies for the treatment of future null-infinity in numerical relativity involve some form of regularization of the field equations. In a recent proposal that relies on the dual foliation formalism this is achieved by the use of an asymptotically Minkowskian generalized harmonic tensor basis. For the scheme to work however, derivatives of certain coordinate light-speeds must decay fast enough. Presently, we generalize the method of asymptotic expansions for nonlinear wave equations to treat first order symmetric hyperbolic systems. We then use this heuristic tool to extract the expected rates of decay of the metric near null-infinity in a free-evolution setting. We show, within the asymptotic expansion, that by carefully modifying the non-principal part of the field equations by the addition of constraints, we are able to obtain optimal decay rates even when the constraints are violated. The light-speed condition can hence be satisfied, which paves the way for the explicit numerical treatment of future null-infinity. We then study the behavior of the Trautman-Bondi mass under the decay results predicted by the asymptotic expansion. Naively the mass seems to be unbounded, but we see first that the divergent terms can be replaced with a combination of the constraints and the Einstein field equations, and second that the Bondi mass loss formula is recovered within the framework. Both of the latter results hold in the presence of small constraint violations.

preprint2016arXiv

The evolution of hyperboloidal data with the dual foliation formalism: Mathematical analysis and wave equation tests

A long-standing problem in numerical relativity is the satisfactory treatment of future null-infinity. We propose an approach for the evolution of hyperboloidal initial data in which the outer boundary of the computational domain is placed at infinity. The main idea is to apply the `dual foliation' formalism in combination with hyperboloidal coordinates and the generalized harmonic gauge formulation. The strength of the present approach is that, following the ideas of Zenginoglu, a hyperboloidal layer can be naturally attached to a central region using standard coordinates of numerical relativity applications. Employing a generalization of the standard hyperboloidal slices, developed by Calabrese et. al., we find that all formally singular terms take a trivial limit as we head to null-infinity. A byproduct is a numerical approach for hyperboloidal evolution of nonlinear wave equations violating the null-condition. The height-function method, used often for fixed background spacetimes, is generalized in such a way that the slices can be dynamically `waggled' to maintain the desired outgoing coordinate lightspeed precisely. This is achieved by dynamically solving the eikonal equation. As a first numerical test of the new approach we solve the 3D flat space scalar wave equation. The simulations, performed with the pseudospectral bamps code, show that outgoing waves are cleanly absorbed at null-infinity and that errors converge away rapidly as resolution is increased.

preprint2015arXiv

A Pseudospectral Method for Gravitational Wave Collapse

We present a new pseudospectral code, bamps, for numerical relativity written with the evolution of collapsing gravitational waves in mind. We employ the first order generalized harmonic gauge formulation. The relevant theory is reviewed and the numerical method is critically examined and specialized for the task at hand. In particular we investigate formulation parameters, gauge and constraint preserving boundary conditions well-suited to non-vanishing gauge source functions. Different types of axisymmetric twist-free moment of time symmetry gravitational wave initial data are discussed. A treatment of the axisymmetric apparent horizon condition is presented with careful attention to regularity on axis. Our apparent horizon finder is then evaluated in a number of test cases. Moving on to evolutions, we investigate modifications to the generalized harmonic gauge constraint damping scheme to improve conservation in the strong field regime. We demonstrate strong-scaling of our pseudospectral penalty code. We employ the Cartoon method to efficiently evolve axisymmetric data in our 3+1 dimensional code. We perform test evolutions of Schwarzschild perturbed by gravitational waves and by gauge pulses, both to demonstrate the use of our blackhole excision scheme and for comparison with earlier results. Finally numerical evolutions of supercritical Brill waves are presented to demonstrate durability of the excision scheme for the dynamical formation of a blackhole.

preprint2015arXiv

Dual Foliation Formulations of General Relativity

A dual foliation treatment of General Relativity is presented. The basic idea of the construction is to consider two foliations of a spacetime by spacelike hypersurfaces and relate the two geometries. The treatment is expected to be useful in various situations, and in particular whenever one would like to compare objects represented in different coordinates. Potential examples include the construction of initial data and the study of trapped tubes. It is common for studies in mathematical relativity to employ a double-null gauge. In such studies local well-posedness is treated by referring back, for example, to the generalized harmonic formulation, global properties of solutions being treated in a separate formalism. As a first application of the dual foliation formulation we find that one can in fact obtain local well-posedness in the double-null coordinates directly, which should allow their use in numerical relativity with standard methods. With due care it is expected that practically any coordinates can be used with this approach.

preprint2015arXiv

Spherical symmetry as a test case for unconstrained hyperboloidal evolution

We consider the hyperboloidal initial value problem in numerical relativity, motivated by the goal to evolve radiating compact objects such as black hole binaries with a numerical grid that includes null infinity. Unconstrained evolution schemes promise optimal efficiency, but are difficult to regularize at null infinity, where the compactified Einstein equations are formally singular. In this work we treat the spherically symmetry case, which already poses nontrivial problems and constitutes an important first step. We have carried out stable numerical evolutions with the generalized BSSN and Z4 equations coupled to a scalar field. The crucial ingredients have been to find an appropriate evolution equation for the lapse function and to adapt constraint damping terms to handle null infinity.

preprint2015arXiv

The Initial Value Formulation of Dynamical Chern-Simons Gravity

We derive an initial value formulation for dynamical Chern-Simons gravity, a modification of general relativity involving parity-violating higher derivative terms. We investigate the structure of the resulting system of partial differential equations thinking about linearization around arbitrary backgrounds. This type of consideration is necessary if we are to establish well-posedness of the Cauchy problem. Treating the field equations as an effective field theory we find that weak necessary conditions for hyperbolicity are satisfied. For the full field equations we find that there are states from which subsequent evolution is not determined. Generically the evolution system closes, but the full field equations are in no sense hyperbolic. In a cursory mode analysis we find that the equations of motion contain terms that may cause ill-posedness of the initial value problem.

preprint2014arXiv

Hyperbolicity of High Order Systems of Evolution Equations

We study properties of evolution equations which are first order in time and arbitrary order in space (FTNS). Following Gundlach and Martín-García (2006) we define strong and symmetric hyperbolicity for FTNS systems and examine the relationship between these definitions, and the analogous concepts for first order systems. We demonstrate equivalence of the FTNS definition of strong hyperbolicity with the existence of a strongly hyperbolic first order reduction. We also demonstrate equivalence of the FTNS definition, up to N=4, of symmetric hyperbolicity with the existence of a symmetric hyperbolic first order reduction.

preprint2013arXiv

An Introduction to Well-posedness and Free-evolution

These lecture notes accompany two classes given at the NRHEP2 school. In the first lecture I introduce the basic concepts used for analyzing well-posedness, that is the existence of a unique solution depending continuously on given data, of evolution partial differential equations. I show how strong hyperbolicity guarantees well-posedness of the initial value problem. Symmetric hyperbolic systems are shown to render the initial boundary value problem well-posed with maximally dissipative boundary conditions. I discuss the Laplace-Fourier method for analyzing the initial boundary value problem. Finally I state how these notions extend to systems that are first order in time and second order in space. In the second lecture I discuss the effect that the gauge freedom of electromagnetism has on the PDE status of the initial value problem. I focus on gauge choices, strong-hyperbolicity and the construction of constraint preserving boundary conditions. I show that strongly hyperbolic pure gauges can be used to build strongly hyperbolic formulations. I examine which of these formulations is additionally symmetric hyperbolic and finally demonstrate that the system can be made boundary stable.

preprint2013arXiv

Collapse of Nonlinear Gravitational Waves in Moving-Puncture Coordinates

We study numerical evolutions of nonlinear gravitational waves in moving-puncture coordinates. We adopt two different types of initial data -- Brill and Teukolsky waves -- and evolve them with two independent codes producing consistent results. We find that Brill data fail to produce long-term evolutions for common choices of coordinates and parameters, unless the initial amplitude is small, while Teukolsky wave initial data lead to stable evolutions, at least for amplitudes sufficiently far from criticality. The critical amplitude separates initial data whose evolutions leave behind flat space from those that lead to a black hole. For the latter we follow the interaction of the wave, the formation of a horizon, and the settling down into a time-independent trumpet geometry. We explore the differences between Brill and Teukolsky data and show that for less common choices of the parameters -- in particular negative amplitudes -- Brill data can be evolved with moving-puncture coordinates, and behave similarly to Teukolsky waves.

preprint2013arXiv

Hyperbolicity of Hamiltonian formulations in General Relativity

Well-posedness of the initial (boundary) value problem is an essential property, both of meaningful physical models and of numerical applications. To prove well-posedness of wave-type equations their level of hyperbolicity is an essential ingredient. We develop helpful tools and classify a large class of Hamiltonian versions of Einstein's equations with live gauge conditions with respect to their hyperbolicity. Finally we find a symmetric hyperbolic Hamiltonian formulation that allows for gauge conditions which are similar to the puncture gauge.

preprint2013arXiv

Hyperbolicity of Physical Theories with Application to General Relativity

We consider gauge theories from the free evolution point of view, in which initial data satisfying constraints of a theory are given. Because the constraints are compatible with the field equations they remain so. We study a model constrained Hamiltonian theory and identify a particular structure in the equations of motion which we call the standard gauge freedom. The pure gauge subsystem of this model theory is identified and the manner in which the gauge variables couple to the field equations is presented. We demonstrate that the set of gauge choices that can be coupled to the field equations to obtain a, properly defined, wave-like formulation is exactly the set of wave-like pure gauges. Consequently we analyze a parametrized family of formulations of general relativity. The generalization of the harmonic gauge formulation to a five parameter family of gauge conditions is obtained.

preprint2012arXiv

Compact binary evolutions with the Z4c formulation

Numerical relativity simulations of compact binaries with the Z4c and BSSNOK formulations are compared. The Z4c formulation is advantageous in every case considered. In simulations of non-vacuum spacetimes the constraint violations due to truncation errors are between one and three orders of magnitude lower in the Z4c evolutions. Improvements are also found in the accuracy of the computed gravitational radiation. For equal-mass irrotational binary neutron star evolutions we find that the absolute errors in phase and amplitude of the waveforms can be up to a factor of four smaller. The quality of the Z4c numerical data is also demonstrated by a remarkably accurate computation of the ADM mass from surface integrals. For equal-mass non-spinning binary puncture black hole evolutions we find that the absolute errors in phase and amplitude of the waveforms can be up to a factor of two smaller. In the same evolutions we find that away from the punctures the Hamiltonian constraint violation is reduced by between one and two orders of magnitude. Furthermore, the utility of gravitational radiation controlling, constraint preserving boundary conditions for the Z4c formulation is demonstrated. The evolution of spacetimes containing a single compact object confirm earlier results in spherical symmetry. The boundary conditions avoid spurious and non-convergent effects present in high resolution runs with either formulation with a more naive boundary treatment. We conclude that Z4c is preferable to BSSNOK for the numerical solution of the 3+1 Einstein equations with the puncture gauge.

preprint2012arXiv

Constraint damping for the Z4c formulation of general relativity

One possibility for avoiding constraint violation in numerical relativity simulations adopting free-evolution schemes is to modify the continuum evolution equations so that constraint violations are damped away. Gundlach et. al. demonstrated that such a scheme damps low amplitude, high frequency constraint violating modes exponentially for the Z4 formulation of General Relativity. Here we analyze the effect of the damping scheme in numerical applications on a conformal decomposition of Z4. After reproducing the theoretically predicted damping rates of constraint violations in the linear regime, we explore numerical solutions not covered by the theoretical analysis. In particular we examine the effect of the damping scheme on low-frequency and on high-amplitude perturbations of flat spacetime as well and on the long-term dynamics of puncture and compact star initial data in the context of spherical symmetry. We find that the damping scheme is effective provided that the constraint violation is resolved on the numerical grid. On grid noise the combination of artificial dissipation and damping helps to suppress constraint violations. We find that care must be taken in choosing the damping parameter in simulations of puncture black holes. Otherwise the damping scheme can cause undesirable growth of the constraints, and even qualitatively incorrect evolutions. In the numerical evolution of a compact static star we find that the choice of the damping parameter is even more delicate, but may lead to a small decrease of constraint violation. For a large range of values it results in unphysical behavior.

preprint2011arXiv

Constraint preserving boundary conditions for the Z4c formulation of general relativity

We discuss high order absorbing constraint preserving boundary conditions for the Z4c formulation of general relativity coupled to the moving puncture family of gauges. We are primarily concerned with the constraint preservation and absorption properties of these conditions. In the frozen coefficient approximation, with an appropriate first order pseudo-differential reduction, we show that the constraint subsystem is boundary stable on a four dimensional compact manifold. We analyze the remainder of the initial boundary value problem for a spherical reduction of the Z4c formulation with a particular choice of the puncture gauge. Numerical evidence for the efficacy of the conditions is presented in spherical symmetry.

preprint2011arXiv

Numerical stability of the Z4c formulation of general relativity

We study numerical stability of different approaches to the discretization of a conformal decomposition of the Z4 formulation of general relativity. We demonstrate that in the linear, constant coefficient regime a novel discretization for tensors is formally numerically stable with a method of lines time-integrator. We then perform a full set of apples with apples tests on the non-linear system, and thus present numerical evidence that both the new and standard discretizations are, in some sense, numerically stable in the non-linear regime. The results of the Z4c numerical tests are compared with those of BSSNOK evolutions. We typically do not employ the Z4c constraint damping scheme and find that in the robust stability and gauge wave tests the Z4c evolutions result in lower constraint violation at the same resolution as the BSSNOK evolutions. In the gauge wave tests we find that the Z4c evolutions maintain the desired convergence factor over many more light-crossing times than the BSSNOK tests. The difference in the remaining tests is marginal.

preprint2011arXiv

Stability of the puncture method with a generalized BSSN formulation

The puncture method for dealing with black holes in the numerical simulation of vacuum spacetimes is remarkably successful when combined with the BSSN formulation of the Einstein equations. We examine a generalized class of formulations modeled along the lines of the Laguna-Shoemaker system, including BSSN as a special case. The formulation is a two parameter generalization of the choice of variables used in standard BSSN evolutions. Numerical stability of the standard finite difference methods is proven for the formulation in the linear regime around flat space, a special case of which is the numerical stability of BSSN. Numerical evolutions are presented and compared with a standard BSSN implementation. We find that a significant portion of the parameter space leads to stable evolutions and that standard BSSN is located near the edge of the stability region. Non-standard parameter choices typically result in smoother behaviour of the evolution variables close to the puncture and thus hold promise for improved accuracy in, e.g., long-term BH binary inspirals, and for overcoming (numerical) stability problems still encountered in some types of black-hole simulations, e.g., in $D \ge 6$ dimensions.

preprint2011arXiv

The trumpet solution from spherical gravitational collapse with puncture gauges

We investigate the stationary end-state obtained by evolving a collapsing spherical star with the gauges routinely adopted to study puncture black holes. We compare the end-state of the collapse with the trumpet solution found in the evolution of a single wormhole slice and show that the two solutions closely agree. We demonstrate that the agreement is caused by the use of the Gamma-driver shift condition, which allows the matter to fall inwards into a region of spacetime that is not resolved by the numerical grid, and which simultaneously finds the stationary coordinates of the trumpet outside the matter.

preprint2010arXiv

Constraint violation in free evolution schemes: comparing BSSNOK with a conformal decomposition of Z4

We compare numerical evolutions performed with the BSSNOK formulation and a conformal decomposition of a Z4-like formulation of General Relativity. The important difference between the two formulations is that the Z4 formulation has a propagating Hamiltonian constraint, whereas BSSNOK has a zero-speed characteristic variable in the constraint subsystem. In spherical symmetry we evolve both puncture and neutron star initial data. We demonstrate that the propagating nature of the Z4 constraints leads to results that compare favorably with BSSNOK evolutions, especially when matter is present in the spacetime. From the point of view of implementation the new system is a simple modification of BSSNOK.

preprint2010arXiv

Hyperbolic formulations of General Relativity with Hamiltonian structure

With the aim of deriving symmetric hyperbolic free-evolution systems for GR that possess Hamiltonian structure and allow for the popular puncture gauge condition we analyze the hyperbolicity of Hamiltonian systems. We develop helpful tools which are applicable to either the first order in time, second order in space or the fully second order form of the equations of motion. For toy models we find that the Hamiltonian structure can simplify the proof of symmetric hyperbolicity. In GR we use a special structure of the principal part to prove symmetric hyperbolicity of a formulation that includes gauge conditions which are very similar to the puncture gauge.