BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

David Hilditch

David Hilditch contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
gr-qchep-th

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

11 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access