Linearized solutions of the Einstein equations within a Bondi-Sachs framework, and implications for boundary conditions in numerical simulations

We linearize the Einstein equations when the metric is Bondi-Sachs, when the background is Schwarzschild or Minkowski, and when there is a matter source in the form of a thin shell whose density varies with time and angular position. By performing an eigenfunction decomposition, we reduce the problem to a system of linear ordinary differential equations which we are able to solve. The solutions are relevant to the characteristic formulation of numerical relativity: (a) as exact solutions against which computations of gravitational radiation can be compared; and (b) in formulating boundary conditions on the $r=2M$ Schwarzschild horizon.Comment: Revised following referee comment

High-powered Gravitational News

We describe the computation of the Bondi news for gravitational radiation. We have implemented a computer code for this problem. We discuss the theory behind it as well as the results of validation tests. Our approach uses the compactified null cone formalism, with the computational domain extending to future null infinity and with a worldtube as inner boundary. We calculate the appropriate full Einstein equations in computational eth form in (a) the interior of the computational domain and (b) on the inner boundary. At future null infinity, we transform the computed data into standard Bondi coordinates and so are able to express the news in terms of its standard $N_{+}$ and $N_{\times}$ polarization components. The resulting code is stable and second-order convergent. It runs successfully even in the highly nonlinear case, and has been tested with the news as high as 400, which represents a gravitational radiation power of about $10^{13}M_{\odot}/sec$.Comment: 24 pages, 4 figures. To appear in Phys. Rev.

Spectral signatures of the Luttinger liquid to charge-density-wave transition

Electron- and phonon spectral functions of the one-dimensional, spinless-fermion Holstein model at half filling are calculated in the four distinct regimes of the phase diagram, corresponding to an attractive or repulsive Luttinger liquid at weak electron-phonon coupling, and a band- or polaronic insulator at strong coupling. The results obtained by means of kernel polynomial and systematic cluster approaches reveal substantially different physics in these regimes and further indicate that the size of the phonon frequency significantly affects the nature of the quantum Peierls phase transition.Comment: 5 pages, 4 figures; final version, accepted for publication in Physical Review

Ill-posedness in the Einstein equations

It is shown that the formulation of the Einstein equations widely in use in numerical relativity, namely, the standard ADM form, as well as some of its variations (including the most recent conformally-decomposed version), suffers from a certain but standard type of ill-posedness. Specifically, the norm of the solution is not bounded by the norm of the initial data irrespective of the data. A long-running numerical experiment is performed as well, showing that the type of ill-posedness observed may not be serious in specific practical applications, as is known from many numerical simulations.Comment: 13 pages, 3 figures, accepted for publication in Journal of Mathematical Physics (to appear August 2000

Gravitational waveforms with controlled accuracy

A partially first-order form of the characteristic formulation is introduced to control the accuracy in the computation of gravitational waveforms produced by highly distorted single black hole spacetimes. Our approach is to reduce the system of equations to first-order differential form on the angular derivatives, while retaining the proven radial and time integration schemes of the standard characteristic formulation. This results in significantly improved accuracy over the standard mixed-order approach in the extremely nonlinear post-merger regime of binary black hole collisions.Comment: Revised version, published in Phys. Rev. D, RevTeX, 16 pages, 4 figure

Cauchy boundaries in linearized gravitational theory

We investigate the numerical stability of Cauchy evolution of linearized gravitational theory in a 3-dimensional bounded domain. Criteria of robust stability are proposed, developed into a testbed and used to study various evolution-boundary algorithms. We construct a standard explicit finite difference code which solves the unconstrained linearized Einstein equations in the 3+1 formulation and measure its stability properties under Dirichlet, Neumann and Sommerfeld boundary conditions. We demonstrate the robust stability of a specific evolution-boundary algorithm under random constraint violating initial data and random boundary data.Comment: 23 pages including 3 figures and 2 tables, revte

Accuracy of numerical relativity waveforms from binary neutron star mergers and their comparison with post-Newtonian waveforms

We present numerical relativity simulations of nine-orbit equal-mass binary neutron star covering the quasicircular late inspiral and merger. The extracted gravitational waveforms are analyzed for convergence and accuracy. Second order convergence is observed up to contact, i.e. about 3-4 cycles to merger; error estimates can be made up to this point. The uncertainties on the phase and the amplitude are dominated by truncation errors and can be minimized to 0.13 rad and less then 1%, respectively, by using several simulations and extrapolating in resolution. In the latter case finite-radius extraction uncertainties become a source of error of the same order and have to be taken into account. The waveforms are tested against accuracy standards for data analysis. The uncertainties on the waveforms are such that accuracy standards are generically not met for signal-to-noise ratios relevant for detection, except for some best cases using extrapolation from several runs. A detailed analysis of the errors is thus imperative for the use of numerical relativity waveforms from binary neutron stars in quantitative studies. The waveforms are compared with the post-Newtonian Taylor T4 approximants both for point-particle and including the analytically known tidal corrections. The T4 approximants accumulate significant phase differences of 2 rad at contact and 4 rad at merger, underestimating the influence of finite size effects. Tidal signatures in the waveforms are thus important at least during the last six orbits of the merger process.Comment: Physical Review D (Vol.85, No.10) 201

Oscillating elastic defects: competition and frustration

We consider a dynamical generalization of the Eshelby problem: the strain profile due to an inclusion or "defect" in an isotropic elastic medium. We show that the higher the oscillation frequency of the defect, the more localized is the strain field around the defect. We then demonstrate that the qualitative nature of the interaction between two defects is strongly dependent on separation, frequency and direction, changing from "ferromagnetic" to "antiferromagnetic" like behavior. We generalize to a finite density of defects and show that the interactions in assemblies of defects can be mapped to XY spin-like models, and describe implications for frustration and frequency-driven pattern transitions.Comment: 4 pages, 5 figure

The ac-Driven Motion of Dislocations in a Weakly Damped Frenkel-Kontorova Lattice

By means of numerical simulations, we demonstrate that ac field can support stably moving collective nonlinear excitations in the form of dislocations (topological solitons, or kinks) in the Frenkel-Kontorova (FK) lattice with weak friction, which was qualitatively predicted by Bonilla and Malomed [Phys. Rev. B{\bf 43}, 11539 (1991)]. Direct generation of the moving dislocations turns out to be virtually impossible; however, they can be generated initially in the lattice subject to an auxiliary spatial modulation of the on-site potential strength. Gradually relaxing the modulation, we are able to get the stable moving dislocations in the uniform FK lattice with the periodic boundary conditions, provided that the driving frequency is close to the gap frequency of the linear excitations in the uniform lattice. The excitations have a large and noninteger index of commensurability with the lattice (suggesting that its actual value is irrational). The simulations reveal two different types of the moving dislocations: broad ones, that extend, roughly, to half the full length of the periodic lattice (in that sense, they cannot be called solitons), and localized soliton-like dislocations, that can be found in an excited state, demonstrating strong persistent internal vibrations. The minimum (threshold) amplitude of the driving force necessary to support the traveling excitation is found as a function of the friction coefficient. Its extrapolation suggests that the threshold does not vanish at the zero friction, which may be explained by radiation losses. The moving dislocation can be observed experimentally in an array of coupled small Josephson junctions in the form of an {\it inverse Josephson effect}, i.e., a dc-voltage response to the uniformly applied ac bias current.Comment: Plain Latex, 13 pages + 9 PostScript figures. to appear on Journal of Physics: condensed matte