Formulation Improvements for Critical Collapse Simulations

Abstract

The precise tuning required to observe critical phenomena in gravitational collapse poses a challenge for most numerical codes. First, threshold estimation searches may be obstructed by the appearance of coordinate singularities, indicating the need for a better gauge choice. Second, the constraint violations to which simulations are susceptible may be too large and force searches to terminate prematurely. This is a particularly serious issue for first order formulations. We want our adaptive pseudospectral code bamps to be a robust tool for the study of critical phenomena so, having encountered both of these difficulties in work on the vacuum setting, we turn here to investigate these issues in the classic context of a spherically symmetric massless scalar field. We suggest two general improvements. We propose a necessary condition for a gauge choice to respect discrete self-similarity (DSS). The condition is not restricted to spherical symmetry and could be verified with any 3+1 formulation. After evaluating common gauge choices against this condition, we suggest a DSS-compatible gauge source function in generalized harmonic gauge (GHG). To control constraint violations, we modify the constraint damping parameters of GHG, adapting them to collapse spacetimes. This allows us to improve our tuning of the critical amplitude for several families of initial data, even going from 6 up to 11 digits. This is the most precise tuning achieved with the first order GHG formulation to date. Consequently, we are able to reproduce the well known critical phenomena as well as competing formulations and methods, clearly observing up to 3 echoes.

arXiv comment19 pages, 10 figures

References in corpus (38)

• On the Numerical Integration of Einstein's Field Equations
• New Formalism for Numerical Relativity
• A New Generalized Harmonic Evolution System
• Constraint damping in the Z4 formulation and harmonic gauge
• Critical phenomena in gravitational collapse
• Conformal and covariant formulation of the Z4 system with constraint-violation damping
• Constraint violation in free evolution schemes: comparing BSSNOK with a conformal decomposition of Z4
• Compact binary evolutions with the Z4c formulation
• Simulations of Binary Black Hole Mergers Using Spectral Methods
• Understanding critical collapse of a scalar field
• Stable radiation-controlling boundary conditions for the generalized harmonic Einstein equations
• Einstein's Equations with Asymptotically Stable Constraint Propagation
• Critical Collapse of the Massless Scalar Field in Axisymmetry
• Hyperbolic slicings of spacetime: singularity avoidance and gauge shocks
• The appearance of coordinate shocks in hyperbolic formalisms of General Realtivity
• An Improved Gauge Driver for the Generalized Harmonic Einstein System
• Scaling of curvature in sub-critical gravitational collapse
• All nonspherical perturbations of the Choptuik spacetime decay
• A Pseudospectral Method for Gravitational Wave Collapse
• Global structure of Choptuik's critical solution in scalar field collapse
• The initial value problem as it relates to numerical relativity
• Symmetry-seeking spacetime coordinates
• Evolutions of centered Brill waves with a pseudospectral method
• Aspherical deformations of the Choptuik spacetime
• Critical Collapse of Scalar Fields Beyond Axisymmetry
• Critical phenomena in the gravitational collapse of electromagnetic waves
• Universality of curvature invariants in critical vacuum gravitational collapse
• Hyperbolic Relaxation Method for Elliptic Equations
• Critical phenomena in the collapse of gravitational waves
• Choptuik's critical spacetime exists
• Critical behavior in 3-d gravitational collapse of massless scalar fields
• Evolution of Brill waves with an adaptive pseudospectral method
• Slicing conditions for axisymmetric gravitational collapse of Brill waves
• Shock-avoiding slicing conditions: tests and calibrations
• Adaptive hp-Refinement for Spectral Elements in Numerical Relativity
• Critical phenomena in the gravitational collapse of electromagnetic dipole and quadrupole waves
• Critical phenomena in the collapse of quadrupolar and hexadecapolar gravitational waves
• An Implementation of DF-GHG with Application to Spherical Black Hole Excision

Cited by in corpus (9)

• Boson star head-on collisions with constraint-violating and constraint-satisfying initial data
• Twist-free axisymmetric critical collapse of a complex scalar field
• Simulations of gravitational collapse in null coordinates: II. Critical collapse of an axisymmetric scalar field
• Robustness of inflation to kinetic inhomogeneities
• Simulating binary black hole mergers using discontinuous Galerkin methods
• Twist and higher modes of a complex scalar field at the threshold of collapse
• Quantum correlations in a gravitational collapse simulation with SpheriCo.jl
• Phase transition analogs in laser collisions with a dark-field setup
• Pseudospectral implementation of the Einstein-Maxwell system