95 research outputs found
Numerical Models of Spin-Orbital Coupling in Neutron Star Binaries
We present a new numerical scheme for solving the initial value problem for
quasiequilibrium binary neutron stars allowing for arbitrary spins. We
construct sequences of circular-orbit binaries of varying separation, keeping
the rest mass and circulation constant along each sequence. The spin angular
frequency of the stars is shown to vary along the sequence, a result that can
be derived analytically in the PPN limit. This spin effect, in addition to
leaving an imprint on the gravitational waveform emitted during binary
inspiral, is measurable in the electromagnetic signal if one of the stars is a
pulsar visible from Earth.Comment: 4 pages, 3 figures. Submitted to the Proceedings of the "X Marcel
Grossmann Meeting on General Relativity" in Rio de Janeiro, Brazil, July
20-26 (2003
Solving the Initial Value Problem of two Black Holes
We solve the elliptic equations associated with the Hamiltonian and momentum
constraints, corresponding to a system composed of two black holes with
arbitrary linear and angular momentum. These new solutions are based on a
Kerr-Schild spacetime slicing which provides more physically realistic
solutions than the initial data based on conformally flat metric/maximal
slicing methods. The singularity/inner boundary problems are circumvented by a
new technique that allows the use of an elliptic solver on a Cartesian grid
where no points are excised, simplifying enormously the numerical problem.Comment: 4 pages, 3 figures. Minor corrections, some points clarified, and one
reference added. To appear in Phys. Rev. Let
Binary black holes on a budget: Simulations using workstations
Binary black hole simulations have traditionally been computationally very
expensive: current simulations are performed in supercomputers involving dozens
if not hundreds of processors, thus systematic studies of the parameter space
of binary black hole encounters still seem prohibitive with current technology.
Here we show how the multi-layered refinement level code BAM can be used on
dual processor workstations to simulate certain binary black hole systems. BAM,
based on the moving punctures method, provides grid structures composed of
boxes of increasing resolution near the center of the grid. In the case of
binaries, the highest resolution boxes are placed around each black hole and
they track them in their orbits until the final merger when a single set of
levels surrounds the black hole remnant. This is particularly useful when
simulating spinning black holes since the gravitational fields gradients are
larger. We present simulations of binaries with equal mass black holes with
spins parallel to the binary axis and intrinsic magnitude of S/m^2= 0.75. Our
results compare favorably to those of previous simulations of this particular
system. We show that the moving punctures method produces stable simulations at
maximum spatial resolutions up to M/160 and for durations of up to the
equivalent of 20 orbital periods.Comment: 20 pages, 8 figures. Final version, to appear in a special issue of
Class. Quantum Grav. based on the New Frontiers in Numerical Relativity
Conference, Golm, July 200
Accurate Evolutions of Orbiting Black-Hole Binaries Without Excision
We present a new algorithm for evolving orbiting black-hole binaries that
does not require excision or a corotating shift. Our algorithm is based on a
novel technique to handle the singular puncture conformal factor. This system,
based on the BSSN formulation of Einstein's equations, when used with a
`pre-collapsed' initial lapse, is non-singular at the start of the evolution,
and remains non-singular and stable provided that a good choice is made for the
gauge. As a test case, we use this technique to fully evolve orbiting
black-hole binaries from near the Innermost Stable Circular Orbit (ISCO)
regime. We show fourth order convergence of waveforms and compute the radiated
gravitational energy and angular momentum from the plunge. These results are in
good agreement with those predicted by the Lazarus approach.Comment: 4 pages, revtex4, 3 figs, references added, typos fixe
Relativistic Models for Binary Neutron Stars with Arbitrary Spins
We introduce a new numerical scheme for solving the initial value problem for
quasiequilibrium binary neutron stars allowing for arbitrary spins. The coupled
Einstein field equations and equations of relativistic hydrodynamics are solved
in the Wilson-Mathews conformal thin sandwich formalism. We construct sequences
of circular-orbit binaries of varying separation, keeping the rest mass and
circulation constant along each sequence. Solutions are presented for
configurations obeying an n=1 polytropic equation of state and spinning
parallel and antiparallel to the orbital angular momentum. We treat stars with
moderate compaction ((m/R) = 0.14) and high compaction ((m/R) = 0.19). For all
but the highest circulation sequences, the spins of the neutron stars increase
as the binary separation decreases. Our zero-circulation cases approximate
irrotational sequences, for which the spin angular frequencies of the stars
increases by 13% (11%) of the orbital frequency for (m/R) = 0.14 ((m/R) = 0.19)
by the time the innermost circular orbit is reached. In addition to leaving an
imprint on the inspiral gravitational waveform, this spin effect is measurable
in the electromagnetic signal if one of the stars is a pulsar visible from
Earth.Comment: 21 pages, 14 figures. A few explanatory sentences added and some
typos corrected. Accepted for publication in Phys. Rev.
Grazing Collisions of Black Holes via the Excision of Singularities
We present the first simulations of non-headon (grazing) collisions of binary
black holes in which the black hole singularities have been excised from the
computational domain. Initially two equal mass black holes are separated a
distance and with impact parameter . Initial data are
based on superposed, boosted (velocity ) solutions of single black
holes in Kerr-Schild coordinates. Both rotating and non-rotating black holes
are considered. The excised regions containing the singularities are specified
by following the dynamics of apparent horizons. Evolutions of up to are obtained in which two initially separate apparent horizons are present
for . At that time a single enveloping apparent horizon forms,
indicating that the holes have merged. Apparent horizon area estimates suggest
gravitational radiation of about 2.6% of the total mass. The evolutions end
after a moderate amount of time because of instabilities.Comment: 2 References corrected, reference to figure update
Momentum constraint relaxation
Full relativistic simulations in three dimensions invariably develop runaway
modes that grow exponentially and are accompanied by violations of the
Hamiltonian and momentum constraints. Recently, we introduced a numerical
method (Hamiltonian relaxation) that greatly reduces the Hamiltonian constraint
violation and helps improve the quality of the numerical model. We present here
a method that controls the violation of the momentum constraint. The method is
based on the addition of a longitudinal component to the traceless extrinsic
curvature generated by a vector potential w_i, as outlined by York. The
components of w_i are relaxed to solve approximately the momentum constraint
equations, pushing slowly the evolution toward the space of solutions of the
constraint equations. We test this method with simulations of binary neutron
stars in circular orbits and show that effectively controls the growth of the
aforementioned violations. We also show that a full numerical enforcement of
the constraints, as opposed to the gentle correction of the momentum relaxation
scheme, results in the development of instabilities that stop the runs shortly.Comment: 17 pages, 10 figures. New numerical tests and references added. More
detailed description of the algorithms are provided. Final published versio
- âŠ