405 research outputs found
The Initial-Boundary Value Problem in General Relativity
In this article we summarize what is known about the initial-boundary value
problem for general relativity and discuss present problems related to it.Comment: 11 pages, 2 figures. Contribution to a special volume for Mario
Castagnino's seventy fifth birthda
Geometrical optics analysis of the short-time stability properties of the Einstein evolution equations
Many alternative formulations of Einstein's evolution have lately been
examined, in an effort to discover one which yields slow growth of
constraint-violating errors. In this paper, rather than directly search for
well-behaved formulations, we instead develop analytic tools to discover which
formulations are particularly ill-behaved. Specifically, we examine the growth
of approximate (geometric-optics) solutions, studied only in the future domain
of dependence of the initial data slice (e.g. we study transients). By
evaluating the amplification of transients a given formulation will produce, we
may therefore eliminate from consideration the most pathological formulations
(e.g. those with numerically-unacceptable amplification). This technique has
the potential to provide surprisingly tight constraints on the set of
formulations one can safely apply. To illustrate the application of these
techniques to practical examples, we apply our technique to the 2-parameter
family of evolution equations proposed by Kidder, Scheel, and Teukolsky,
focusing in particular on flat space (in Rindler coordinates) and Schwarzchild
(in Painleve-Gullstrand coordinates).Comment: Submitted to Phys. Rev.
Constraint-preserving boundary treatment for a harmonic formulation of the Einstein equations
We present a set of well-posed constraint-preserving boundary conditions for
a first-order in time, second-order in space, harmonic formulation of the
Einstein equations. The boundary conditions are tested using robust stability,
linear and nonlinear waves, and are found to be both less reflective and
constraint preserving than standard Sommerfeld-type boundary conditions.Comment: 18 pages, 7 figures, accepted in CQ
Modeling the source of GW150914 with targeted numerical-relativity simulations
In fall of 2015, the two LIGO detectors measured the gravitational wave
signal GW150914, which originated from a pair of merging black holes. In the
final 0.2 seconds (about 8 gravitational-wave cycles) before the amplitude
reached its maximum, the observed signal swept up in amplitude and frequency,
from 35 Hz to 150 Hz. The theoretical gravitational-wave signal for merging
black holes, as predicted by general relativity, can be computed only by full
numerical relativity, because analytic approximations fail near the time of
merger. Moreover, the nearly-equal masses, moderate spins, and small number of
orbits of GW150914 are especially straightforward and efficient to simulate
with modern numerical-relativity codes. In this paper, we report the modeling
of GW150914 with numerical-relativity simulations, using black-hole masses and
spins consistent with those inferred from LIGO's measurement. In particular, we
employ two independent numerical-relativity codes that use completely different
analytical and numerical methods to model the same merging black holes and to
compute the emitted gravitational waveform; we find excellent agreement between
the waveforms produced by the two independent codes. These results demonstrate
the validity, impact, and potential of current and future studies using
rapid-response, targeted numerical-relativity simulations for better
understanding gravitational-wave observations.Comment: 11 pages, 3 figures, submitted to Classical and Quantum Gravit
Modeling Gravitational Recoil Using Numerical Relativity
We review the developments in modeling gravitational recoil from merging
black-hole binaries and introduce a new set of 20 simulations to test our
previously proposed empirical formula for the recoil. The configurations are
chosen to represent generic binaries with unequal masses and precessing spins.
Results of these simulations indicate that the recoil formula is accurate to
within a few km/s in the similar mass-ratio regime for the out-of-plane recoil.Comment: corrections to text, 11 pages, 1 figur
Black Hole Mergers and Unstable Circular Orbits
We describe recent numerical simulations of the merger of a class of equal
mass, non-spinning, eccentric binary black hole systems in general relativity.
We show that with appropriate fine-tuning of the initial conditions to a region
of parameter space we denote the threshold of immediate merger, the binary
enters a phase of close interaction in a near-circular orbit, stays there for
an amount of time proportional to logarithmic distance from the threshold in
parameter space, then either separates or merges to form a single Kerr black
hole. To gain a better understanding of this phenomena we study an analogous
problem in the evolution of equatorial geodesics about a central Kerr black
hole. A similar threshold of capture exists for appropriate classes of initial
conditions, and tuning to threshold the geodesics approach one of the unstable
circular geodesics of the Kerr spacetime. Remarkably, with a natural mapping of
the parameters of the geodesic to that of the equal mass system, the scaling
exponent describing the whirl phase of each system turns out to be quite
similar. Armed with this lone piece of evidence that an approximate
correspondence might exist between near-threshold evolution of geodesics and
generic binary mergers, we illustrate how this information can be used to
estimate the cross section and energy emitted in the ultra relativistic black
hole scattering problem. This could eventually be of use in providing estimates
for the related problem of parton collisions at the Large Hadron Collider in
extra dimension scenarios where black holes are produced.Comment: 16 pages, 12 figures; updated to coincide with journal versio
The discrete energy method in numerical relativity: Towards long-term stability
The energy method can be used to identify well-posed initial boundary value
problems for quasi-linear, symmetric hyperbolic partial differential equations
with maximally dissipative boundary conditions. A similar analysis of the
discrete system can be used to construct stable finite difference equations for
these problems at the linear level. In this paper we apply these techniques to
some test problems commonly used in numerical relativity and observe that while
we obtain convergent schemes, fast growing modes, or ``artificial
instabilities,'' contaminate the solution. We find that these growing modes can
partially arise from the lack of a Leibnitz rule for discrete derivatives and
discuss ways to limit this spurious growth.Comment: 18 pages, 22 figure
Towards a Realistic Neutron Star Binary Inspiral: Initial Data and Multiple Orbit Evolution in Full General Relativity
This paper reports on our effort in modeling realistic astrophysical neutron
star binaries in general relativity. We analyze under what conditions the
conformally flat quasiequilibrium (CFQE) approach can generate
``astrophysically relevant'' initial data, by developing an analysis that
determines the violation of the CFQE approximation in the evolution of the
binary described by the full Einstein theory. We show that the CFQE assumptions
significantly violate the Einstein field equations for corotating neutron stars
at orbital separations nearly double that of the innermost stable circular
orbit (ISCO) separation, thus calling into question the astrophysical relevance
of the ISCO determined in the CFQE approach. With the need to start numerical
simulations at large orbital separation in mind, we push for stable and long
term integrations of the full Einstein equations for the binary neutron star
system. We demonstrate the stability of our numerical treatment and analyze the
stringent requirements on resolution and size of the computational domain for
an accurate simulation of the system.Comment: 22 pages, 18 figures, accepted to Phys. Rev.
Simulation of Binary Black Hole Spacetimes with a Harmonic Evolution Scheme
A numerical solution scheme for the Einstein field equations based on
generalized harmonic coordinates is described, focusing on details not provided
before in the literature and that are of particular relevance to the binary
black hole problem. This includes demonstrations of the effectiveness of
constraint damping, and how the time slicing can be controlled through the use
of a source function evolution equation. In addition, some results from an
ongoing study of binary black hole coalescence, where the black holes are
formed via scalar field collapse, are shown. Scalar fields offer a convenient
route to exploring certain aspects of black hole interactions, and one
interesting, though tentative suggestion from this early study is that behavior
reminiscent of "zoom-whirl" orbits in particle trajectories is also present in
the merger of equal mass, non-spinning binaries, with appropriately fine-tuned
initial conditions.Comment: 16 pages, 14 figures; replaced with published versio
- …