3,795 research outputs found
Model waveform accuracy standards for gravitational wave data analysis
Model waveforms are used in gravitational wave data analysis to detect and then to measure the properties of a source by matching the model waveforms to the signal from a detector. This paper derives accuracy standards for model waveforms which are sufficient to ensure that these data analysis applications are capable of extracting the full scientific content of the data, but without demanding excessive accuracy that would place undue burdens on the model waveform simulation community. These accuracy standards are intended primarily for broadband model waveforms produced by numerical simulations, but the standards are quite general and apply equally to such waveforms produced by analytical or hybrid analytical-numerical methods
Nonlinear Couplings of R-modes: Energy Transfer and Saturation Amplitudes at Realistic Timescales
Non-linear interactions among the inertial modes of a rotating fluid can be
described by a network of coupled oscillators. We use such a description for an
incompressible fluid to study the development of the r-mode instability of
rotating neutron stars. A previous hydrodynamical simulation of the r-mode
reported the catastrophic decay of large amplitude r-modes. We explain the
dynamics and timescale of this decay analytically by means of a single three
mode coupling. We argue that at realistic driving and damping rates such large
amplitudes will never actually be reached. By numerically integrating a network
of nearly 5000 coupled modes, we find that the linear growth of the r-mode
ceases before it reaches an amplitude of around 10^(-4). The lowest parametric
instability thresholds for the r-mode are calculated and it is found that the
r-mode becomes unstable to modes with 13<n<15 if modes up to n=30 are included.
Using the network of coupled oscillators, integration times of 10^6 rotational
periods are attainable for realistic values of driving and damping rates.
Complicated dynamics of the modal amplitudes are observed. The initial
development is governed by the three mode coupling with the lowest parametric
instability. Subsequently a large number of modes are excited, which greatly
decreases the linear growth rate of the r-mode.Comment: 3 figures 4 pages Submitted to PR
Solving Einstein's Equations With Dual Coordinate Frames
A method is introduced for solving Einstein's equations using two distinct
coordinate systems. The coordinate basis vectors associated with one system are
used to project out components of the metric and other fields, in analogy with
the way fields are projected onto an orthonormal tetrad basis. These field
components are then determined as functions of a second independent coordinate
system. The transformation to the second coordinate system can be thought of as
a mapping from the original ``inertial'' coordinate system to the computational
domain. This dual-coordinate method is used to perform stable numerical
evolutions of a black-hole spacetime using the generalized harmonic form of
Einstein's equations in coordinates that rotate with respect to the inertial
frame at infinity; such evolutions are found to be generically unstable using a
single rotating coordinate frame. The dual-coordinate method is also used here
to evolve binary black-hole spacetimes for several orbits. The great
flexibility of this method allows comoving coordinates to be adjusted with a
feedback control system that keeps the excision boundaries of the holes within
their respective apparent horizons.Comment: Updated to agree with published versio
Gauge drivers for the generalized harmonic Einstein equations
The generalized harmonic representation of Einstein's equations is manifestly hyperbolic for a large class of gauge conditions. Unfortunately most of the useful gauges developed over the past several decades by the numerical relativity community are incompatible with the hyperbolicity of the equations in this form. This paper presents a new method of imposing gauge conditions that preserves hyperbolicity for a much wider class of conditions, including as special cases many of the standard ones used in numerical relativity: e.g., K freezing, Gamma freezing, Bona-Massó slicing, conformal Gamma drivers, etc. Analytical and numerical results are presented which test the stability and the effectiveness of this new gauge-driver evolution system
Shear viscosity of neutron matter from realistic nucleon-nucleon interactions
The calculation of transport properties of Fermi liquids, based on the
formalism developed by Abrikosov and Khalatnikov, requires the knowledge of the
probability of collisions between quasiparticles in the vicinity of the Fermi
surface. We have carried out a numerical study of the shear viscosity of pure
neutron matter, whose value plays a pivotal role in determining the stability
of rotating neutron stars, in which these processes are described using a
state-of-the-art nucleon-nucleon potential model. Within our approach medium
modifications of the scattering cross section are consistently taken into
account, through an effective interaction obtained from the matrix elements of
the bare interaction between correlated states. Inclusion of medium effects
lead to a large increase of the viscosity at densities larger than
fm^{-3}.Comment: 4 pages, 4 figures. Corrected typo
A New Generalized Harmonic Evolution System
A new representation of the Einstein evolution equations is presented that is
first order, linearly degenerate, and symmetric hyperbolic. This new system
uses the generalized harmonic method to specify the coordinates, and
exponentially suppresses all small short-wavelength constraint violations.
Physical and constraint-preserving boundary conditions are derived for this
system, and numerical tests that demonstrate the effectiveness of the
constraint suppression properties and the constraint-preserving boundary
conditions are presented.Comment: Updated to agree with published versio
R-Modes in Superfluid Neutron Stars
The analogs of r-modes in superfluid neutron stars are studied here. These
modes, which are governed primarily by the Coriolis force, are identical to
their ordinary-fluid counterparts at the lowest order in the small
angular-velocity expansion used here. The equations that determine the next
order terms are derived and solved numerically for fairly realistic superfluid
neutron-star models. The damping of these modes by superfluid ``mutual
friction'' (which vanishes at the lowest order in this expansion) is found to
have a characteristic time-scale of about 10^4 s for the m=2 r-mode in a
``typical'' superfluid neutron-star model. This time-scale is far too long to
allow mutual friction to suppress the recently discovered gravitational
radiation driven instability in the r-modes. However, the strength of the
mutual friction damping depends very sensitively on the details of the
neutron-star core superfluid. A small fraction of the presently acceptable
range of superfluid models have characteristic mutual friction damping times
that are short enough (i.e. shorter than about 5 s) to suppress the
gravitational radiation driven instability completely.Comment: 15 pages, 8 figure
Controlling the growth of constraints in hyperbolic evolution systems
Motivated by the need to control the exponential growth of constraint violations in numerical solutions of the Einstein evolution equations, two methods are studied here for controlling this growth in general hyperbolic evolution systems. The first method adjusts the evolution equations dynamically, by adding multiples of the constraints, in a way designed to minimize this growth. The second method imposes special constraint preserving boundary conditions on the incoming components of the dynamical fields. The efficacy of these methods is tested by using them to control the growth of constraints in fully dynamical 3D numerical solutions of a particular representation of the Maxwell equations that is subject to constraint violations. The constraint preserving boundary conditions are found to be much more effective than active constraint control in the case of this Maxwell system
Mining for Observables: A New Challenge in Numerical Relativity
One of the motivations behind numerical relativity is to provide gravitational wave signals of compact objects to observers using the new gravitational wave detectors. Yet, because of the complexities involved, no dependable signals of binary-black hole coalescences have been established. The work in this proceedings is motivated by how numerical relativity can be used today to predict robust features in gravitational wave signals of binary black-hole coalescence by making approximations to the full problem. To illustrate this, we present results from evolving a Klein-Gordon equation on a frozen background. The background is set by a sequence of initial data in which the binary is in quasi-equilibrium. We probe the data resulting from the evolution for the transition between the linear and non-linear regimes using oscillations of the black holes as our guide. This information is used to motivate a qualitative picture of the gravitational signal of a black-hole coalescence
Reducing orbital eccentricity in binary black hole simulations
Binary black hole simulations starting from quasi-circular (i.e., zero radial
velocity) initial data have orbits with small but non-zero orbital
eccentricities. In this paper the quasi-equilibrium initial-data method is
extended to allow non-zero radial velocities to be specified in binary black
hole initial data. New low-eccentricity initial data are obtained by adjusting
the orbital frequency and radial velocities to minimize the orbital
eccentricity, and the resulting ( orbit) evolutions are compared with
those of quasi-circular initial data. Evolutions of the quasi-circular data
clearly show eccentric orbits, with eccentricity that decays over time. The
precise decay rate depends on the definition of eccentricity; if defined in
terms of variations in the orbital frequency, the decay rate agrees well with
the prediction of Peters (1964). The gravitational waveforms, which contain
cycles in the dominant l=m=2 mode, are largely unaffected by the
eccentricity of the quasi-circular initial data. The overlap between the
dominant mode in the quasi-circular evolution and the same mode in the
low-eccentricity evolution is about 0.99.Comment: 27 pages, 9 figures; various minor clarifications; accepted to the
"New Frontiers" special issue of CQ
- …