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.

Atomic scale lattice distortions and domain wall profiles

We present an atomic scale theory of lattice distortions using strain related variables and their constraint equations. Our approach connects constrained {\it atomic length} scale variations to {\it continuum} elasticity and describes elasticity at several length scales. We apply the approach to a two-dimensional square lattice with a monatomic basis, and find the elastic deformations and hierarchical atomic relaxations in the vicinity of a domain wall between two different homogeneous strain states. We clarify the microscopic origin of gradient terms, some of which are included phenomenologically in Ginzburg-Landau theory, by showing that they are anisotropic.Comment: 6 figure

Two Steps Block Method for the Solution of General Second Order Initial Value Problems of Ordinary Differential Equation

In this paper, an implicit block linear multistep method for the solution of ordinary differential equation was extended to the general form of differential equation. This method is self starting and does not need a predictor to solve for the unknown in the corrector. The method can also be extended to boundary value problems without additional cost. The method was found to be efficient after being tested with numerical problems of second order

Electron-Phonon Driven Spin Frustration in Multi-Band Hubbard Models: MX Chains and Oxide Superconductors

We discuss the consequences of both electron-phonon and electron-electron couplings in 1D and 2D multi-band (Peierls-Hubbard) models. After briefly discussing various analytic limits, we focus on (Hartree-Fock and exact) numerical studies in the intermediate regime for both couplings, where unusual spin-Peierls as well as long-period, frustrated ground states are found. Doping into such phases or near the phase boundaries can lead to further interesting phenomena such as separation of spin and charge, a dopant-induced phase transition of the global (parent) phase, or real-space (``bipolaronic'') pairing. We discuss possible experimentally observable consequences of this rich phase diagram for halogen-bridged, transition metal, linear chain complexes (MX chains) in 1D and the oxide superconductors in 2D.Comment: 6 pages, four postscript figures (appended), in regular Te

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

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