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 Hole Mergers in 3d Numerical Relativity

The standard approach to the numerical evolution of black hole data using the ADM formulation with maximal slicing and vanishing shift is extended to non-symmetric black hole data containing black holes with linear momentum and spin by using a time-independent conformal rescaling based on the puncture representation of the black holes. We give an example for a concrete three dimensional numerical implementation. The main result of the simulations is that this approach allows for the first time to evolve through a brief period of the merger phase of the black hole inspiral.Comment: 8 pages, 9 figures, REVTeX; expanded discussion, results unchange

Cosmic Strings from Supersymmetric Flat Directions

Flat directions are a generic feature of the scalar potential in supersymmetric gauge field theories. They can arise, for example, from D-terms associated with an extra abelian gauge symmetry. Even when supersymmetry is broken softly, there often remain directions in the scalar field space along which the potential is almost flat. Upon breaking a gauge symmetry along one of these almost flat directions, cosmic strings may form. Relative to the standard cosmic string picture based on the abelian Higgs model, these flat-direction cosmic strings have the extreme Type-I properties of a thin gauge core surrounded by a much wider scalar field profile. We perform a comprehensive study of the microscopic, macroscopic, and observational characteristics of this class of strings. We find many differences from the standard string scenario, including stable higher winding mode strings, the dynamical formation of higher mode strings from lower ones, and a resultant multi-tension scaling string network in the early universe. These strings are only moderately constrained by current observations, and their gravitational wave signatures may be detectable at future gravity wave detectors. Furthermore, there is the interesting but speculative prospect that the decays of cosmic string loops in the early universe could be a source of ultra-high energy cosmic rays or non-thermal dark matter. We also compare the observational signatures of flat-direction cosmic strings with those of ordinary cosmic strings as well as (p,q) cosmic strings motivated by superstring theory.Comment: 58 pages, 16 figures, v2. accepted to PRD, added comments about baryogenesis and boosted decay products from cusp annihilatio

Acidification of Forest Soils: A Model for Analyzing Impacts of Acidic Deposition in Europe - Version II

Acidification is an unfavorable process in forest soils. Timber logging, natural accumulation of biomass in the ecosystem, and acidic deposition are sources of acidification. Acidification causes a risk of damage to plant roots and a subsequent risk of a decline in ecosystem productivity. A dynamic model is introduced for describing the acidification of forest soils. In one-year time steps the model calculates the soil pH as function of acid stress and the buffer mechanisms of the soil. Acid stress is defined as the hydrogen ion input into the top soil. The buffer mechanisms counteract acidification by providing a sink for hydrogen ions. The concepts buffer rate and buffer capacity are used to quantify the buffer mechanisms. The model compares (i) the rate of the acid stress (annual amount) to the buffer rate, and (ii) the accumulated acid stress (over several years) to the buffer capacity. The comparisons produce an estimate of the soil acidity as the output. Since the first version in May 1984 several changes have been implemented following the advice of the experts. For aluminum and iron buffer ranges an equilibrium approach has been introduced. The pH of the silicate, cation exchange and upper aluminum buffer ranges is now a function of base saturation. In the current version of the model forests are assumed to absorb sulfur compounds more effectively than agricultural lands and, moreover, forests are assumed to grow on poor soil types rather than on the average soil type of a grid. The model system as a whole is now available for analyzing the impact of different emission scenarios. The soil acidification model assumes sulfur deposition estimates from the other submodels as input, and as output it computes the total area of forests in Europe with the estimated soil pH lower than any selected threshold value. Additionally it produces estimates of the acidity of European forest soils in a map format

Generic Tracking of Multiple Apparent Horizons with Level Flow

We report the development of the first apparent horizon locator capable of finding multiple apparent horizons in a ``generic'' numerical black hole spacetime. We use a level-flow method which, starting from a single arbitrary initial trial surface, can undergo topology changes as it flows towards disjoint apparent horizons if they are present. The level flow method has two advantages: 1) The solution is independent of changes in the initial guess and 2) The solution can have multiple components. We illustrate our method of locating apparent horizons by tracking horizon components in a short Kerr-Schild binary black hole grazing collision.Comment: 13 pages including figures, submitted to 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 $m$ are separated a distance $\approx10m$ and with impact parameter $\approx2m$. Initial data are based on superposed, boosted (velocity $\approx0.5c$) 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 $t \approx 35m$ are obtained in which two initially separate apparent horizons are present for $t\approx3.8m$. 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

Radiation from low-momentum zoom-whirl orbits

We study zoom-whirl behaviour of equal mass, non-spinning black hole binaries in full general relativity. The magnitude of the linear momentum of the initial data is fixed to that of a quasi-circular orbit, and its direction is varied. We find a global maximum in radiated energy for a configuration which completes roughly one orbit. The radiated energy in this case exceeds the value of a quasi-circular binary with the same momentum by 15%. The direction parameter only requires minor tuning for the localization of the maximum. There is non-trivial dependence of the energy radiated on eccentricity (several local maxima and minima). Correlations with orbital dynamics shortly before merger are discussed. While being strongly gauge dependent, these findings are intuitive from a physical point of view and support basic ideas about the efficiency of gravitational radiation from a binary system.Comment: 9 pages, 6 figures, Amaldi8 conference proceedings as publishe

Kinematic Constraints on Formation of Bound States of Cosmic Strings - Field Theoretical Approach

Superstring theory predicts the potential formation of string networks with bound states ending in junctions. Kinematic constraints for junction formation have been derived within the Nambu-Goto thin string approximation. Here we test these constraints numerically in the framework of the Abelian-Higgs model in the Type-I regime and report on good agreement with the analytical predictions. We also demonstrate that strings can effectively pass through each other when they meet at speeds slightly above the critical velocity permitting bound state formation. This is due to reconnection effects that are beyond the scope of the Nambu-Goto approximation.Comment: 6 pages, 12 eps figures - matches the published versio