Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userEquilibrium initial data for moving puncture simulations: The stationary 1+log slicing
181
0
0.0
(
0
)
Added by
Thomas W. Baumgarte
Publication date
2009
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We propose and explore a stationary 1+log slicing condition for the construction of solutions to Einsteins constraint equations. For stationary spacetimes, these initial data will give a stationary foliation when evolved with moving puncture gauge conditions that are often used in black hole evolutions. The resulting slicing is time-independent and agrees with the slicing generated by being dragged along a time-like Killing vector of the spacetime. When these initial data are evolved with moving puncture gauge conditions, numerical errors arising from coordinate evolution are minimized. In the construction of initial data for binary black holes it is often assumed that there exists an approximate helical Killing vector that generates the binarys orbit. We show that, unfortunately, 1+log slices that are stationary with respect to such a helical Killing vector cannot be asymptotically flat, unless the spacetime possesses an additional axial Killing vector.
rate research
Read More
We present a single domain Galerkin-Collocation method to calculate puncture initial data sets for single and binary, either in the trumpet or wormhole geometries. The combination of aspects belonging to the Galerkin and the Collocation methods together with the adoption of spherical coordinates in all cases show to be very effective. We have proposed a unified expression for the conformal factor to describe trumpet and spinning black holes. In particular, for the spinning trumpet black holes, we have exhibited the deformation of the limit surface due to the spin from a sphere to an oblate spheroid. We have also revisited the energy content in the trumpet and wormhole puncture data sets. The algorithm can be extended to describe binary black holes.
We demonstrate that numerical relativity codes based on the moving punctures formalism are capable of evolving nearly maximally spinning black hole binaries. We compare a new evolution of an equal-mass, aligned-spin binary with dimensionless spin chi=0.99 using puncture-based data with recent simulations of the SXS Collaboration. We find that the overlap of our new waveform with the published results of the SXS Collaboration is larger than 0.999. To generate our new waveform, we use the recently introduced HiSpID puncture data, the CCZ4 evolution system, and a modified lapse condition that helps keep the horizon radii reasonably large.
This lecture is devoted to the problem of computing initial data for the Cauchy problem of 3+1 general relativity. The main task is to solve the constraint equations. The conformal technique, introduced by Lichnerowicz and enhanced by York, is presented. Two standard methods, the conformal transverse-traceless one and the conformal thin sandwich, are discussed and illustrated by some simple examples. Finally a short review regarding initial data for binary systems (black holes and neutron stars) is given.
We study the effectiveness of stationary-phase approximated post-Newtonian waveforms currently used by ground-based gravitational-wave detectors to search for the coalescence of binary black holes by comparing them to an accurate waveform obtained from numerical simulation of an equal-mass non-spinning binary black hole inspiral, merger and ringdown. We perform this study for the Initial- and Advanced-LIGO detectors. We find that overlaps between the templates and signal can be improved by integrating the match filter to higher frequencies than used currently. We propose simple analytic frequency cutoffs for both Initial and Advanced LIGO, which achieve nearly optimal matches, and can easily be extended to unequal-mass, spinning systems. We also find that templates that include terms in the phase evolution up to 3.5 pN order are nearly always better, and rarely significantly worse, than 2.0 pN templates currently in use. For Initial LIGO we recommend a strategy using templates that include a recently introduced pseudo-4.0 pN term in the low-mass ($M leq 35 MSun$) region, and 3.5 pN templates allowing unphysical values of the symmetric reduced mass $eta$ above this. This strategy always achieves overlaps within 0.3% of the optimum, for the data used here. For Advanced LIGO we recommend a strategy using 3.5 pN templates up to $M=12 MSun$, 2.0 pN templates up to $M=21 MSun$, pseudo-4.0 pN templates up to $65 MSun$, and 3.5 pN templates with unphysical $eta$ for higher masses. This strategy always achieves overlaps within 0.7% of the optimum for Advanced LIGO.
Binary neutron stars in circular orbits can be modeled as helically symmetric, i.e., stationary in a rotating frame. This symmetry gives rise to a first integral of the Euler equation, often employed for constructing equilibrium solutions via iteration. For eccentric orbits, however, the lack of helical symmetry has prevented the use of this method, and the numerical relativity community has often resorted to constructing initial data by superimposing boosted spherical stars without solving the Euler equation. The spuriously excited neutron star oscillations seen in evolutions of such data arise because such configurations lack the appropriate tidal deformations and are stationary in a linearly comoving---rather than rotating---frame. We consider eccentric configurations at apoapsis that are instantaneously stationary in a rotating frame. We extend the notion of helical symmetry to eccentric orbits, by approximating the elliptical orbit of each companion as instantaneously circular, using the ellipses inscribed circle. The two inscribed helical symmetry vectors give rise to approximate instantaneous first integrals of the Euler equation throughout each companion. We use these integrals as the basis of a self-consistent iteration of the Einstein constraints to construct conformal thin-sandwich initial data for eccentric binaries. We find that the spurious stellar oscillations are reduced by at least an order of magnitude, compared with those found in evolutions of superposed initial data. The tidally induced oscillations, however, are physical and qualitatively similar to earlier evolutions. Finally, we show how to incorporate radial velocity due to radiation reaction in our inscribed helical symmetry vectors, which would allow one to obtain truly non-eccentric initial data when our eccentricity parameter $e$ is set to zero.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
