-
Graviton-Photon Oscillations in an Expanding Universe
Authors:
Peter Anninos,
Tony Rothman,
Andrea Palessandro
Abstract:
Through the Gertsenshtein effect, the presence of a large external B-field may allow photons and gravitons to mix in a way that resembles neutrino oscillations and is even more similar to axion-photon mixing. Assuming a background B-field (or any field that behaves like one), we examine the Gertsenshtein mechanism in the context of FLRW expanding universe models, as well as de Sitter space. The co…
▽ More
Through the Gertsenshtein effect, the presence of a large external B-field may allow photons and gravitons to mix in a way that resembles neutrino oscillations and is even more similar to axion-photon mixing. Assuming a background B-field (or any field that behaves like one), we examine the Gertsenshtein mechanism in the context of FLRW expanding universe models, as well as de Sitter space. The conformal invariance of Maxwell's equations and the conformal noninvariance of the Einstein equations preclude the operation of the Gertsenshtein effect at all scales. In general we find for the matter- and radiation-dominated cases, graviton-oscillations are possible only at late conformal times or when the wavelengths are much shorter than the horizon ($kη\gg 1$), but that the time-dependent oscillations eventually damp out in any case. The presence of charged particles additionally damps out the oscillations. For the de Sitter universe, we find that oscillations are possible only at early conformal times ($η\ll H^{-1}$) and for wavelengths short compared to the Hubble radius, but eventually freeze in when wavelengths become longer than the Hubble radius. In principle a Gertsenshtein-like mechanism might influence the balance of particle species in an inflationary phase before freezing in; however, we find that in all our models the mixing length is larger than the Hubble radius. We discuss several possible remedies to this situation.
△ Less
Submitted 1 April, 2024; v1 submitted 16 January, 2024;
originally announced January 2024.
-
The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project
Authors:
Oliver Porth,
Koushik Chatterjee,
Ramesh Narayan,
Charles F. Gammie,
Yosuke Mizuno,
Peter Anninos,
John G. Baker,
Matteo Bugli,
Chi-kwan Chan,
Jordy Davelaar,
Luca Del Zanna,
Zachariah B. Etienne,
P. Chris Fragile,
Bernard J. Kelly,
Matthew Liska,
Sera Markoff,
Jonathan C. McKinney,
Bhupendra Mishra,
Scott C. Noble,
Héctor Olivares,
Ben Prather,
Luciano Rezzolla,
Benjamin R. Ryan,
James M. Stone,
Niccolò Tomei
, et al. (3 additional authors not shown)
Abstract:
Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic…
▽ More
Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL. Agreement between the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems.
△ Less
Submitted 5 August, 2019; v1 submitted 9 April, 2019;
originally announced April 2019.
-
Hydrodynamic Stability of Cosmological Quark-Hadron Phase Transitions
Authors:
P. Chris Fragile,
Peter Anninos
Abstract:
Recent results from linear perturbation theory suggest that first-order cosmological quark-hadron phase transitions occurring as deflagrations may be ``borderline'' unstable, and those occurring as detonations may give rise to growing modes behind the interface boundary. However, since nonlinear effects can play important roles in the development of perturbations, unstable behavior cannot be ass…
▽ More
Recent results from linear perturbation theory suggest that first-order cosmological quark-hadron phase transitions occurring as deflagrations may be ``borderline'' unstable, and those occurring as detonations may give rise to growing modes behind the interface boundary. However, since nonlinear effects can play important roles in the development of perturbations, unstable behavior cannot be asserted entirely by linear analysis, and the uncertainty of these recent studies is compounded further by nonlinearities in the hydrodynamics and self-interaction fields. In this paper we investigate the growth of perturbations and the stability of quark-hadron phase transitions in the early Universe by solving numerically the fully nonlinear relativistic hydrodynamics equations coupled to a scalar field with a quartic self-interaction potential regulating the transitions. We consider single, perturbed, phase transitions propagating either by detonation or deflagration, as well as multiple phase and shock front interactions in 1+2 dimensional spacetimes.
△ Less
Submitted 4 March, 2003;
originally announced March 2003.
-
Non-Oscillatory Central Difference and Artificial Viscosity Schemes for Relativistic Hydrodynamics
Authors:
Peter Anninos,
P. Chris Fragile
Abstract:
High resolution, non-oscillatory, central difference (NOCD) numerical schemes are introduced as alternatives to more traditional artificial viscosity (AV) and Godunov methods for solving the fully general relativistic hydrodynamics equations. These new approaches provide the advantages of Godunov methods in capturing ultra-relativistic flows without the cost and complication of Riemann solvers,…
▽ More
High resolution, non-oscillatory, central difference (NOCD) numerical schemes are introduced as alternatives to more traditional artificial viscosity (AV) and Godunov methods for solving the fully general relativistic hydrodynamics equations. These new approaches provide the advantages of Godunov methods in capturing ultra-relativistic flows without the cost and complication of Riemann solvers, and the advantages of AV methods in their speed, ease of implementation, and general applicability without explicitly using artificial viscosity for shock capturing. Shock tube, wall shock, and dust accretion tests, all with adiabatic equations of state, are presented and compared against equivalent solutions from both AV and Godunov based codes. In the process we address the accuracy of time-explicit NOCD and AV methods over a wide range of Lorentz factors.
△ Less
Submitted 16 September, 2002; v1 submitted 16 June, 2002;
originally announced June 2002.
-
Computational Cosmology: from the Early Universe to the Large Scale Structure
Authors:
Peter Anninos
Abstract:
In order to account for the observable Universe, any comprehensive theory or model of cosmology must draw from many disciplines of physics, including gauge theories of strong and weak interactions, the hydrodynamics and microphysics of baryonic matter, electromagnetic fields, and spacetime curvature, for example. Although it is difficult to incorporate all these physical elements into a single c…
▽ More
In order to account for the observable Universe, any comprehensive theory or model of cosmology must draw from many disciplines of physics, including gauge theories of strong and weak interactions, the hydrodynamics and microphysics of baryonic matter, electromagnetic fields, and spacetime curvature, for example. Although it is difficult to incorporate all these physical elements into a single complete model of our Universe, advances in computing methods and technologies have contributed significantly towards our understanding of cosmological models, the Universe, and astrophysical processes within them. A sample of numerical calculations (and numerical methods) applied to specific issues in cosmology are reviewed in this article: from the Big Bang singularity dynamics to the fundamental interactions of gravitational waves; from the quark-hadron phase transition to the large scale structure of the Universe. The emphasis, although not exclusively, is on those calculations designed to test different models of cosmology against the observed Universe.
△ Less
Submitted 30 August, 2001;
originally announced August 2001.
-
Instability of Extremal Relativistic Charged Spheres
Authors:
Peter Anninos,
Tony Rothman
Abstract:
With the question, ``Can relativistic charged spheres form extremal black holes?" in mind, we investigate the properties of such spheres from a classical point of view. The investigation is carried out numerically by integrating the Oppenheimer-Volkov equation for relativistic charged fluid spheres and finding interior Reissner-Nordström solutions for these objects. We consider both constant den…
▽ More
With the question, ``Can relativistic charged spheres form extremal black holes?" in mind, we investigate the properties of such spheres from a classical point of view. The investigation is carried out numerically by integrating the Oppenheimer-Volkov equation for relativistic charged fluid spheres and finding interior Reissner-Nordström solutions for these objects. We consider both constant density and adiabatic equations of state, as well as several possible charge distributions, and examine stability by both a normal mode and an energy analysis. In all cases, the stability limit for these spheres lies between the extremal ($Q = M$) limit and the black hole limit ($R = R_+$). That is, we find that charged spheres undergo gravitational collapse before they reach $Q = M$, suggesting that extremal Reissner-Nordtröm black holes produced by collapse are ruled out. A general proof of this statement would support a strong form of the cosmic censorship hypothesis, excluding not only stable naked singularities, but stable extremal black holes. The numerical results also indicate that although the interior mass-energy $m(R)$ obeys the usual $m/R < 4/9$ stability limit for the Schwarzschild interior solution, the gravitational mass $M$ does not. Indeed, the stability limit approaches $R_+$ as $Q \to M$. In the Appendix we also argue that Hawking radiation will not lead to an extremal Reissner-Nordström black hole. All our results are consistent with the third law of black hole dynamics, as currently understood.
△ Less
Submitted 30 August, 2001;
originally announced August 2001.
-
Radiation recoil from highly distorted black holes
Authors:
Steven Brandt,
Peter Anninos
Abstract:
We present results from numerical evolutions of single black holes distorted by axisymmetric, but equatorially asymmetric, gravitational (Brill) waves. Net radiated energies, apparent horizon embeddings, and recoil velocities are shown for a range of Brill wave parameters, including both even and odd parity distortions of Schwarzschild black holes. We find that a wave packet initially concentrat…
▽ More
We present results from numerical evolutions of single black holes distorted by axisymmetric, but equatorially asymmetric, gravitational (Brill) waves. Net radiated energies, apparent horizon embeddings, and recoil velocities are shown for a range of Brill wave parameters, including both even and odd parity distortions of Schwarzschild black holes. We find that a wave packet initially concentrated on the black hole throat, a likely model also for highly asymmetric stellar collapse and late stage binary mergers, can generate a maximum recoil velocity of about 150 (23) km/sec for even (odd) parity perturbations, significantly less than that required to eject black holes from galactic cores.
△ Less
Submitted 6 July, 1999;
originally announced July 1999.
-
Head--on Collision of Two Unequal Mass Black Holes
Authors:
Peter Anninos,
Steven Brandt
Abstract:
We present results from the first fully nonlinear numerical calculations of the head--on collision of two unequal mass black holes. Selected waveforms of the most dominant l=2, 3 and 4 quasinormal modes are shown, as are the total radiated energies and recoil velocities for a range of mass ratios and initial separations. Our results validate the close and distant separation limit perturbation st…
▽ More
We present results from the first fully nonlinear numerical calculations of the head--on collision of two unequal mass black holes. Selected waveforms of the most dominant l=2, 3 and 4 quasinormal modes are shown, as are the total radiated energies and recoil velocities for a range of mass ratios and initial separations. Our results validate the close and distant separation limit perturbation studies, and suggest that the head--on collision scenario is not likely to produce an astrophysically significant recoil effect.
△ Less
Submitted 5 June, 1998;
originally announced June 1998.
-
Stable characteristic evolution of generic 3-dimensional single-black-hole spacetimes
Authors:
The Binary Black Hole Grand Challenge Alliance,
:,
R. Gomez,
L. Lehner,
R. Marsa,
J. Winicour,
A. Abrahams,
A. Anderson,
P. Anninos,
T. Baumgarte,
N. Bishop,
S. Brandt J. Browne,
K. Camarda,
M. Choptuik,
R. Correl,
G. Cook,
C. Evans,
L. Finn,
G. Fox,
T. Haupt,
M. Huq,
L. Kidder,
S. Klasky,
P. Laguna,
W. Landry
, et al. (20 additional authors not shown)
Abstract:
We report new results which establish that the accurate 3-dimensional numerical simulation of generic single-black-hole spacetimes has been achieved by characteristic evolution with unlimited long term stability. Our results cover a selection of distorted, moving and spinning single black holes, with evolution times up to 60,000M.
We report new results which establish that the accurate 3-dimensional numerical simulation of generic single-black-hole spacetimes has been achieved by characteristic evolution with unlimited long term stability. Our results cover a selection of distorted, moving and spinning single black holes, with evolution times up to 60,000M.
△ Less
Submitted 20 January, 1998;
originally announced January 1998.
-
New Coordinate Systems for Axisymmetric Black Hole Collisions
Authors:
P. Anninos,
S. R. Brandt,
P. Walker
Abstract:
We describe a grid generation procedure designed to construct new classes of orthogonal coordinate systems for binary black hole spacetimes. The computed coordinates offer an alternative approach to current methods, in addition to providing a framework for potentially more stable and accurate evolutions of colliding black holes. As a particular example, we apply our procedure to generate appropr…
▽ More
We describe a grid generation procedure designed to construct new classes of orthogonal coordinate systems for binary black hole spacetimes. The computed coordinates offer an alternative approach to current methods, in addition to providing a framework for potentially more stable and accurate evolutions of colliding black holes. As a particular example, we apply our procedure to generate appropriate numerical grids to evolve Misner's axisymmetric initial data set representing two equal mass black holes colliding head-on. These new results are compared with previously published calculations, and we find generally good agreement in both the waveform profiles and total radiated energies over the allowable range of separation parameters. Furthermore, because no specialized treatment of the coordinate singularities is required, these new grids are more easily extendible to unequal mass and spinning black hole collisions.
△ Less
Submitted 10 December, 1997;
originally announced December 1997.
-
Boosted three-dimensional black-hole evolutions with singularity excision
Authors:
The Binary Black Hole Grand Challenge Alliance,
:,
G. B. Cook,
M. F. Huq,
S. A. Klasky,
M. A. Scheel,
A. M. Abrahams,
A. Anderson,
P. Anninos,
T. W. Baumgarte,
N. T. Bishop,
S. R. Brandt,
J. C. Browne,
K. Camarda,
M. W. Choptuik,
C. R. Evans,
L. S. Finn,
G. C. Fox,
R. Gomez,
T. Haupt,
L. E. Kidder,
P. Laguna,
W. Landry,
L. Lehner,
J. Lenaghan
, et al. (21 additional authors not shown)
Abstract:
Binary black hole interactions provide potentially the strongest source of gravitational radiation for detectors currently under development. We present some results from the Binary Black Hole Grand Challenge Alliance three- dimensional Cauchy evolution module. These constitute essential steps towards modeling such interactions and predicting gravitational radiation waveforms. We report on singl…
▽ More
Binary black hole interactions provide potentially the strongest source of gravitational radiation for detectors currently under development. We present some results from the Binary Black Hole Grand Challenge Alliance three- dimensional Cauchy evolution module. These constitute essential steps towards modeling such interactions and predicting gravitational radiation waveforms. We report on single black hole evolutions and the first successful demonstration of a black hole moving freely through a three-dimensional computational grid via a Cauchy evolution: a hole moving ~6M at 0.1c during a total evolution of duration ~60M.
△ Less
Submitted 26 November, 1997;
originally announced November 1997.
-
New Coordinate Systems for Axisymmetric Black Hole Collisions
Authors:
S. R. Brandt,
P. Walker,
P. Anninos
Abstract:
We describe a numerical grid generating procedure to construct new classes of orthogonal coordinate systems that are specially adapted to binary black hole spacetimes. The new coordinates offer an alternative approach to the conventional Čadež coordinates, in addition to providing a potentially more stable and flexible platform to extend previous calculations of binary black hole collisions.
We describe a numerical grid generating procedure to construct new classes of orthogonal coordinate systems that are specially adapted to binary black hole spacetimes. The new coordinates offer an alternative approach to the conventional Čadež coordinates, in addition to providing a potentially more stable and flexible platform to extend previous calculations of binary black hole collisions.
△ Less
Submitted 6 November, 1997;
originally announced November 1997.
-
Gravitational wave extraction and outer boundary conditions by perturbative matching
Authors:
The Binary Black Hole Grand Challenge Alliance,
:,
A. M. Abrahams,
L. Rezzolla,
M. E. Rupright,
A. Anderson,
P. Anninos,
T. W. Baumgarte,
N. T. Bishop,
S. R. Brandt,
J. C. Browne,
K. Camarda,
M. W. Choptuik,
G. B. Cook,
C. R. Evans,
L. S. Finn,
G. Fox,
R. Gomez,
T. Haupt,
M. F. Huq,
L. E. Kidder,
S. Klasky,
P. Laguna,
W. Landry,
L. Lehner
, et al. (20 additional authors not shown)
Abstract:
We present a method for extracting gravitational radiation from a three-dimensional numerical relativity simulation and, using the extracted data, to provide outer boundary conditions. The method treats dynamical gravitational variables as nonspherical perturbations of Schwarzschild geometry. We discuss a code which implements this method and present results of tests which have been performed wi…
▽ More
We present a method for extracting gravitational radiation from a three-dimensional numerical relativity simulation and, using the extracted data, to provide outer boundary conditions. The method treats dynamical gravitational variables as nonspherical perturbations of Schwarzschild geometry. We discuss a code which implements this method and present results of tests which have been performed with a three dimensional numerical relativity code.
△ Less
Submitted 30 September, 1997;
originally announced September 1997.
-
A Phase Space Approach to the Gravitational Arrow of Time
Authors:
Tony Rothman,
Peter Anninos
Abstract:
We attempt to find a function that characterizes gravitational clumping and that increases monotonically as inhomogeneity increases. We choose $S = lnΩ$ as the candidate ``gravitational entropy'' function, where $Ω$ is the phase-space volume below the Hamiltonian H of the system under consideration. We compute $Ω$ for transverse electromagnetic waves and for gravitational wave, radiation and den…
▽ More
We attempt to find a function that characterizes gravitational clumping and that increases monotonically as inhomogeneity increases. We choose $S = lnΩ$ as the candidate ``gravitational entropy'' function, where $Ω$ is the phase-space volume below the Hamiltonian H of the system under consideration. We compute $Ω$ for transverse electromagnetic waves and for gravitational wave, radiation and density perturbations in an expanding FLRW universe. These calculations are carried out in the linear regime under the assumption that the phases of the oscillators comprising the system are random. Entropy is thus attributed to the lack of knowledge of the exact field configuration. We find that $Ω$, and hence $lnΩ$ behaves as required. We also carry out calculations for Bianchi IX cosmological models and find that, even in this homogeneous case, the function can be interpreted sensibly. We compare our results with Penrose's C^2 hypothesis. Because S is defined to resemble the fundamental statistical mechanics definition of entropy, we are able to recover the entropy in a variety of familiar circumstances including, evidently, black-hole entropy. The results point to the utility of the relativistic ADM Hamiltonian formalism in establishing a connection between general relativity and statistical mechanics, although fully nonlinear calculations will need to be performed to remove any doubt.
△ Less
Submitted 21 December, 1996;
originally announced December 1996.
-
Dynamics of Gravitational Waves in 3D: Formulations, Methods, and Tests
Authors:
Peter Anninos,
Joan Masso,
Edward Seidel,
Wai-Mo Suen,
Malcolm Tobias
Abstract:
The dynamics of gravitational waves is investigated in full 3+1 dimensional numerical relativity, emphasizing the difficulties that one might encounter in numerical evolutions, particularly those arising from non-linearities and gauge degrees of freedom. Using gravitational waves with amplitudes low enough that one has a good understanding of the physics involved, but large enough to enable non-…
▽ More
The dynamics of gravitational waves is investigated in full 3+1 dimensional numerical relativity, emphasizing the difficulties that one might encounter in numerical evolutions, particularly those arising from non-linearities and gauge degrees of freedom. Using gravitational waves with amplitudes low enough that one has a good understanding of the physics involved, but large enough to enable non-linear effects to emerge, we study the coupling between numerical errors, coordinate effects, and the nonlinearities of the theory. We discuss the various strategies used in identifying specific features of the evolution. We show the importance of the flexibility of being able to use different numerical schemes, different slicing conditions, different formulations of the Einstein equations (standard ADM vs. first order hyperbolic), and different sets of equations (linearized vs. full Einstein equations). A non-linear scalar field equation is presented which captures some properties of the full Einstein equations, and has been useful in our understanding of the coupling between finite differencing errors and non-linearites. We present a set of monitoring devices which have been crucial in our studying of the waves, including Riemann invariants, pseudo-energy momentum tensor, hamiltonian constraint violation, and fourier spectrum analysis.
△ Less
Submitted 24 October, 1996;
originally announced October 1996.
-
Finding Apparent Horizons in Dynamic 3D Numerical Spacetimes
Authors:
P. Anninos,
K. Camarda,
J. Libson,
J. Massó,
E. Seidel,
W. -M. Suen
Abstract:
We have developed a general method for finding apparent horizons in 3D numerical relativity. Instead of solving for the partial differential equation describing the location of the apparent horizons, we expand the closed 2D surfaces in terms of symmetric trace--free tensors and solve for the expansion coefficients using a minimization procedure. Our method is applied to a number of different spa…
▽ More
We have developed a general method for finding apparent horizons in 3D numerical relativity. Instead of solving for the partial differential equation describing the location of the apparent horizons, we expand the closed 2D surfaces in terms of symmetric trace--free tensors and solve for the expansion coefficients using a minimization procedure. Our method is applied to a number of different spacetimes, including numerically constructed spacetimes containing highly distorted axisymmetric black holes in spherical coordinates, and 3D rotating, and colliding black holes in Cartesian coordinates.
△ Less
Submitted 20 May, 1998; v1 submitted 25 September, 1996;
originally announced September 1996.
-
The collision of boosted black holes
Authors:
John Baker,
Andrew Abrahams,
Peter Anninos,
Steve Brandt,
Richard Price,
Jorge Pullin,
Edward Seidel
Abstract:
We study the radiation from a collision of black holes with equal and opposite linear momenta. Results are presented from a full numerical relativity treatment and are compared with the results from a ``close-slow'' approximation. The agreement is remarkable, and suggests several insights about the generation of gravitational radiation in black hole collisions.
We study the radiation from a collision of black holes with equal and opposite linear momenta. Results are presented from a full numerical relativity treatment and are compared with the results from a ``close-slow'' approximation. The agreement is remarkable, and suggests several insights about the generation of gravitational radiation in black hole collisions.
△ Less
Submitted 27 August, 1996;
originally announced August 1996.
-
The Near-Linear Regime of Gravitational Waves in Numerical Relativity
Authors:
Peter Anninos,
Joan Masso,
Edward Seidel,
Wai-Mo Suen,
Malcolm Tobias
Abstract:
We report on a systematic study of the dynamics of gravitational waves in full 3D numerical relativity. We find that there exists an interesting regime in the parameter space of the wave configurations: a near-linear regime in which the amplitude of the wave is low enough that one expects the geometric deviation from flat spacetime to be negligible, but nevertheless where nonlinearities can exci…
▽ More
We report on a systematic study of the dynamics of gravitational waves in full 3D numerical relativity. We find that there exists an interesting regime in the parameter space of the wave configurations: a near-linear regime in which the amplitude of the wave is low enough that one expects the geometric deviation from flat spacetime to be negligible, but nevertheless where nonlinearities can excite unstable modes of the Einstein evolution equations causing the metric functions to evolve out of control. The implications of this for numerical relativity are discussed.
△ Less
Submitted 24 October, 1996; v1 submitted 17 January, 1996;
originally announced January 1996.
-
Head-on collision of two black holes: comparison of different approaches
Authors:
Peter Anninos,
Richard H. Price,
Jorge Pullin,
Ed Seidel,
Wai-Mo Suen
Abstract:
A benchmark problem for numerical relativity has been the head-on collision of two black holes starting from the ``Misner initial data,'' a closed form momentarily stationary solution to the constraint equations with an adjustable closeness parameter $μ_0$. We show here how an eclectic mixture of approximation methods can provide both an efficient means of determining the time development of the…
▽ More
A benchmark problem for numerical relativity has been the head-on collision of two black holes starting from the ``Misner initial data,'' a closed form momentarily stationary solution to the constraint equations with an adjustable closeness parameter $μ_0$. We show here how an eclectic mixture of approximation methods can provide both an efficient means of determining the time development of the initial data and a good understanding of the physics of the problem. When the Misner data is chosen to correspond to holes initially very close together, a common horizon surrounds both holes and the geometry exterior to the horizon can be treated as a non-spherical perturbation of a single Schwarzschild hole. When the holes are initially well separated the problem can be treated with a different approximation scheme, ``the particle-membrane method.'' For all initial separations, numerical relativity is in principle applicable, but is costly and of uncertain accuracy. We present here a comparison of the different approaches. We compare waveforms, for $\ell=2$ and $\ell=4$ radiation, for different values of $μ_0$, from the three different approaches to the problem.
△ Less
Submitted 23 May, 1995;
originally announced May 1995.
-
Three dimensional numerical relativity: the evolution of black holes
Authors:
Peter Anninos,
Karen Camarda,
Joan Masso,
Edward Seidel,
Wai-Mo Suen,
John Towns
Abstract:
We report on a new 3D numerical code designed to solve the Einstein equations for general vacuum spacetimes. This code is based on the standard 3+1 approach using cartesian coordinates. We discuss the numerical techniques used in developing this code, and its performance on massively parallel and vector supercomputers. As a test case, we present evolutions for the first 3D black hole spacetimes. W…
▽ More
We report on a new 3D numerical code designed to solve the Einstein equations for general vacuum spacetimes. This code is based on the standard 3+1 approach using cartesian coordinates. We discuss the numerical techniques used in developing this code, and its performance on massively parallel and vector supercomputers. As a test case, we present evolutions for the first 3D black hole spacetimes. We identify a number of difficulties in evolving 3D black holes and suggest approaches to overcome them. We show how special treatment of the conformal factor can lead to more accurate evolution, and discuss techniques we developed to handle black hole spacetimes in the absence of symmetries. Many different slicing conditions are tested, including geodesic, maximal, and various algebraic conditions on the lapse. With current resolutions, limited by computer memory sizes, we show that with certain lapse conditions we can evolve the black hole to about $t=50M$, where $M$ is the black hole mass. Comparisons are made with results obtained by evolving spherical initial black hole data sets with a 1D spherically symmetric code. We also demonstrate that an ``apparent horizon locking shift'' can be used to prevent the development of large gradients in the metric functions that result from singularity avoiding time slicings. We compute the mass of the apparent horizon in these spacetimes, and find that in many cases it can be conserved to within about 5\% throughout the evolution with our techniques and current resolution.
△ Less
Submitted 14 March, 1995;
originally announced March 1995.
-
Horizon Boundary Condition for Black Hole Spacetimes
Authors:
Peter Anninos,
Greg Daues,
Joan Masso,
Edward Seidel,
Wai-Mo Suen
Abstract:
It was recently shown that spacetime singularities in numerical relativity could be avoided by excising a region inside the apparent horizon in numerical evolutions. In this paper we report on the details of the implementation of this scheme. The scheme is based on using (1)~a horizon locking coordinate which locks the coordinate system to the geometry, and (2)~a finite differencing scheme which…
▽ More
It was recently shown that spacetime singularities in numerical relativity could be avoided by excising a region inside the apparent horizon in numerical evolutions. In this paper we report on the details of the implementation of this scheme. The scheme is based on using (1)~a horizon locking coordinate which locks the coordinate system to the geometry, and (2)~a finite differencing scheme which respects the causal structure of the spacetime. We show that the horizon locking coordinate can be affected by a number of shift conditions, such as a ``distance freezing'' shift, an ``area freezing'' shift, an ``expansion freezing'' shift, or the minimal distortion shift. The causal differencing scheme is illustrated with the evolution of scalar fields, and its use in evolving the Einstein equations is studied. We compare the results of numerical evolutions with and without the use of this horizon boundary condition scheme for spherical black hole spacetimes. With the boundary condition a black hole can be evolved accurately well beyond $t=1000 M$, where $M$ is the black hole mass.
△ Less
Submitted 22 December, 1994;
originally announced December 1994.
-
3D Numerical Relativity at NCSA
Authors:
Peter Anninos,
Karen Camarda,
Joan Masso,
Edward Seidel,
Wai-Mo Suen,
Malcolm Tobias,
John Towns
Abstract:
Numerical relativity is finally coming of age with the development of massively parallel computers. 3D problems, which were completely intractable several years ago due to limited computer power, can now be performed with grid sizes of about $200^3$. We report on several new codes developed for solving the full 3D Einstein equations, and present results of using them to evolve black holes and gr…
▽ More
Numerical relativity is finally coming of age with the development of massively parallel computers. 3D problems, which were completely intractable several years ago due to limited computer power, can now be performed with grid sizes of about $200^3$. We report on several new codes developed for solving the full 3D Einstein equations, and present results of using them to evolve black holes and gravitational waves.
△ Less
Submitted 20 December, 1994;
originally announced December 1994.
-
Event Horizons of Numerical Black Holes
Authors:
Peter Anninos,
David Bernstein,
Steven Brandt,
Joan Masso,
Edward Seidel,
Larry Smarr,
Wai-Mo Suen,
Paul Walker
Abstract:
We have developed a powerful and efficient method for locating the event horizon of a black hole spacetime, making possible the study of the dynamics of event horizons in numerical relativity. We describe the method and apply it to a colliding black hole spacetime.
We have developed a powerful and efficient method for locating the event horizon of a black hole spacetime, making possible the study of the dynamics of event horizons in numerical relativity. We describe the method and apply it to a colliding black hole spacetime.
△ Less
Submitted 20 December, 1994;
originally announced December 1994.
-
Singularity Avoidance in Numerical Black Hole Spacetimes
Authors:
Peter Anninos,
Greg Daues,
Joan Masso,
Edward Seidel,
Wai-Mo Suen
Abstract:
Spacetime singularities in numerical relativity can be avoided by excising a region of the computational domain from inside the apparent horizon. We report on results of such a scheme that is based on using ({\it i}) a horizon locking coordinate which locks the coordinate system to the geometry, and ({\it ii}) a finite differencing scheme which respects the causal structure of the spacetime. Wit…
▽ More
Spacetime singularities in numerical relativity can be avoided by excising a region of the computational domain from inside the apparent horizon. We report on results of such a scheme that is based on using ({\it i}) a horizon locking coordinate which locks the coordinate system to the geometry, and ({\it ii}) a finite differencing scheme which respects the causal structure of the spacetime. With this technique a black hole can be evolved accurately well beyond $t=1000M$, where $M$ is the black hole mass.
△ Less
Submitted 20 December, 1994;
originally announced December 1994.
-
The Head-On Collision of Two Equal Mass Black Holes: Numerical Methods
Authors:
Peter Anninos,
David Hobill,
Edward Seidel,
Larry Smarr,
Wai-Mo Suen
Abstract:
The head-on collision of two nonrotating axisymmetric equal mass black holes is treated numerically. We take as initial data the single parameter family of time-symmetric solutions discovered by Misner which consists of two Einstein-Rosen bridges that can be placed arbitrarily distant from one another. A number of problems associated with previous attempts to evolve these data sets have been ove…
▽ More
The head-on collision of two nonrotating axisymmetric equal mass black holes is treated numerically. We take as initial data the single parameter family of time-symmetric solutions discovered by Misner which consists of two Einstein-Rosen bridges that can be placed arbitrarily distant from one another. A number of problems associated with previous attempts to evolve these data sets have been overcome. In this article, we discuss our choices for coordinate systems, gauges, and the numerical algorithms that we have developed to evolve this system.
△ Less
Submitted 31 August, 1994;
originally announced August 1994.
-
Dynamics of Apparent and Event Horizons
Authors:
Peter Anninos,
David Bernstein,
Steven Brandt,
Joseph Libson,
Joan Masso,
Edward Seidel,
Larry Smarr,
Wai-Mo Suen,
Paul Walker
Abstract:
The dynamics of apparent and event horizons of various black hole spacetimes, including those containing distorted, rotating and colliding black holes, are studied. We have developed a powerful and efficient new method for locating the event horizon, making possible the study of both types of horizons in numerical relativity. We show that both the event and apparent horizons, in all dynamical bl…
▽ More
The dynamics of apparent and event horizons of various black hole spacetimes, including those containing distorted, rotating and colliding black holes, are studied. We have developed a powerful and efficient new method for locating the event horizon, making possible the study of both types of horizons in numerical relativity. We show that both the event and apparent horizons, in all dynamical black hole spacetimes studied, oscillate with the quasinormal frequency.
△ Less
Submitted 3 March, 1994;
originally announced March 1994.
-
The Collision of Two Black Holes
Authors:
Peter Anninos,
David Hobill,
Edward Seidel,
Larry Smarr,
Wai-Mo Suen
Abstract:
We study the head-on collision of two equal mass, nonrotating black holes. We consider a range of cases from holes surrounded by a common horizon to holes initially separated by about $20M$, where $M$ is the mass of each hole. We determine the waveforms and energies radiated for both the $\ell = 2$ and $\ell=4$ waves resulting from the collision. In all cases studied the normal modes of the fina…
▽ More
We study the head-on collision of two equal mass, nonrotating black holes. We consider a range of cases from holes surrounded by a common horizon to holes initially separated by about $20M$, where $M$ is the mass of each hole. We determine the waveforms and energies radiated for both the $\ell = 2$ and $\ell=4$ waves resulting from the collision. In all cases studied the normal modes of the final black hole dominate the spectrum. We also estimate analytically the total gravitational radiation emitted, taking into account the tidal heating of horizons using the membrane paradigm, and other effects. For the first time we are able to compare analytic calculations, black hole perturbation theory, and strong field, nonlinear numerical calculations for this problem, and we find excellent agreement.
△ Less
Submitted 14 September, 1993;
originally announced September 1993.