Modeling the Black Hole Excision Problem

We analyze the excision strategy for simulating black holes. The problem is modeled by the propagation of quasi-linear waves in a 1-dimensional spatial region with timelike outer boundary, spacelike inner boundary and a horizon in between. Proofs of well-posed evolution and boundary algorithms for a second differential order treatment of the system are given for the separate pieces underlying the finite difference problem. These are implemented in a numerical code which gives accurate long term simulations of the quasi-linear excision problem. Excitation of long wavelength exponential modes, which are latent in the problem, are suppressed using conservation laws for the discretized system. The techniques are designed to apply directly to recent codes for the Einstein equations based upon the harmonic formulation.Comment: 21 pages, 14 postscript figures, minor contents updat

Observation of coherent π0\pi^0 electroproduction on deuterons at large momentum transfer

The first experimental results for coherent $\pi^0$-electroproduction on the deuteron, $e+d\to e+d +\pi^0$, at large momentum transfer, are reported. The experiment was performed at Jefferson Laboratory at an incident electron energy of 4.05 GeV. A large pion production yield has been observed in a kinematical region for 1.1$<Q^2<$1.8 GeV$^2$, from threshold to 200 MeV excitation energy in the $d\pi^0$ system. The $Q^2$-dependence is compared with theoretical predictions.Comment: 26 page

Oscillations and instabilities of fast and differentially rotating relativistic stars

We study non-axisymmetric oscillations of rapidly and differentially rotating relativistic stars in the Cowling approximation. Our equilibrium models are sequences of relativistic polytropes, where the differential rotation is described by the relativistic $j$-constant law. We show that a small degree of differential rotation raises the critical rotation value for which the quadrupolar f-mode becomes prone to the CFS instability, while the critical value of $T/|W|$ at the mass-shedding limit is raised even more. For softer equations of state these effects are even more pronounced. When increasing differential rotation further to a high degree, the neutral point of the CFS instability first reaches a local maximum and is lowered afterwards. For stars with a rather high compactness we find that for a high degree of differential rotation the absolute value of the critical $T/|W|$ is below the corresponding value for rigid rotation. We conclude that the parameter space where the CFS instability is able to drive the neutron star unstable is increased for a small degree of differential rotation and for a large degree at least in stars with a higher compactness.Comment: 16 pages, 11 figures; paper accepted for publication in Phys. Rev. D (81.084019

Mira's wind explored in scattering infrared CO lines

We have observed the intermediate regions of the circumstellar envelope of Mira (o Ceti) in photospheric light scattered by three vibration-rotation transitions of the fundamental band of CO, from low-excited rotational levels of the ground vibrational state, at an angular distance of beta = 2"-7" away from the star. The data were obtained with the Phoenix spectrometer mounted on the 4 m Mayall telescope at Kitt Peak. The spatial resolution is approximately 0.5" and seeing limited. Our observations provide absolute fluxes, leading to an independent new estimate of the mass-loss rate of approximately 3e-7 Msun/yr, as derived from a simple analytic wind model. We find that the scattered intensity from the wind of Mira for 2" < beta < 7" decreases as beta^-3, which suggests a time constant mass-loss rate, when averaged over 100 years, over the past 1200 years.Comment: accepted for publication in the Astrophysical Journa

Constraint preserving boundary conditions for the Z4c formulation of general relativity

We discuss high order absorbing constraint preserving boundary conditions for the Z4c formulation of general relativity coupled to the moving puncture family of gauges. We are primarily concerned with the constraint preservation and absorption properties of these conditions. In the frozen coefficient approximation, with an appropriate first order pseudo-differential reduction, we show that the constraint subsystem is boundary stable on a four dimensional compact manifold. We analyze the remainder of the initial boundary value problem for a spherical reduction of the Z4c formulation with a particular choice of the puncture gauge. Numerical evidence for the efficacy of the conditions is presented in spherical symmetry.Comment: 18 pages, 8 figure

Determination of the reaction plane in ultrarelativistic nuclear collisions

In the particles produced in a nuclear collision undergo collective flow, the reaction plane can in principle be determined through a global event analysis. We show here that collective flow can be identified by evaluating the reaction plane independently in two separate rapidity intervals, and studying the correlation between the two results. We give an analytical expression for the correlation function between the two planes as a function of their relative angle. We also discuss how this correlation function is related to the anisotropy of the transverse momentum distribution. Email contact: [email protected]: Saclay-T93/026 Email: [email protected]

Constraint-preserving Sommerfeld conditions for the harmonic Einstein equations

The principle part of Einstein equations in the harmonic gauge consists of a constrained system of 10 curved space wave equations for the components of the space-time metric. A new formulation of constraint-preserving boundary conditions of the Sommerfeld type for such systems has recently been proposed. We implement these boundary conditions in a nonlinear 3D evolution code and test their accuracy.Comment: 16 pages, 17 figures, submitted to Phys. Rev.