Title: Dynamical evolution of unstable selfgravitating scalar solitons'
1Dynamical evolution of unstable self-gravitating
scalar solitons.
José A. González
In collaboration with M.Alcubierre, M.Salgado
and D. Sudarsky
2Overview
I. Spherically symmetric code
- Introduction
- Regularity conditions
- Regularization algorithm
- Examples
II. Scalar solitons
- Theories with scalar hair
- Numerical evolutions
3Introduction
When developing spherically symmetric codes in
numerical relativity, the coordinate singularity
at the origin can be a source of serious
instabilities caused by the lack of regularity of
the geometric variables there. The problem
arises because of the presence of terms in the
evolution equations that go as 1/r near the
origin. Regularity of the metric guarantees the
exact cancellations of such terms at the origin,
thus ensuring well-behaved solutions.
4This exact cancellation however, though certainly
true for analytical solutions, usually fails to
hold for numerical solutions. One the finds that
the 1/r terms do not cancel and the numerical
solution becomes ill-behaved near r0. The usual
way to deal with this problem is to use the areal
gauge, where the radial coordinate r is chosen in
such a way that the proper area of spheres of
constant r is always . If, moreover, one
also choses a vanishing shift vector one ends up
in the standar polar/areal gauge.
5In this gauge the exact cancellation of the 1/r
terms is now reduced to imposing the boundary
condition at ,which can be
easily done if one solves for from the
hamiltonian constraint and ignores its evolution
equation. The main drawback of the standard
approach is that the gauge choice has been
completely exhausted. In particular, the
polar/areal gauge can not penetrate apparent
horizons, since inside an apparent horizon it is
impossible to keep the areas of spheres
fixed without a non-trivial shift vector.
631 formalism
7Constraints
Evolution equations
8Regularity Conditions
9Evolution equations
10(No Transcript)
11Constraints
Hamiltonian
Momentum
12with
perhaps functions of time
13Regulatization Algorithm
Arbona and Bona (1999) develop a
regularization technique for the spherically
symmetric version of the Bona-Masso evolution
system. Their technique is based on redefining
the auxiliary Dynamical varible that is part
of the standard BM formulation. Our
regularization method is similar in spirit, if
not in detail, to this one. The main difference
being that the approach of AB was tied to the use
of the BM evolution system, while our algorithm
is much more general.
14We define an auxiliary variable
The behavoir at the origin can easily be
imposed numerically using a grid that staggers
the origin, and asking for to be odd
across We need the evolution equation for .
This can be obtained directly from its definition.
15The last evolution equation clearly has
the dangerous
term, but this can be removed with the help of
the momentum constraint to find
16Consider some arbitrary first order formulation
of the Einstein evolution equations in spherical
symmetry that has the generic form
where
The formulation might be hyperbolic or not,
depending on the characteristic structure of the
matrix M. We will assume that one has arrived at
such a formulation by adding multiples of the
constraints
17Using the evolution equation of
If we define
then
18The final step is to substitute the spatial
derivative of
for that of
with
Using that the spatial derivative of is
19we finally find
This last equation is now regular, and has
precisely the same characteristic structure as
the original system.
20Example
Using the Hamiltonian and Momentum constraints to
remove the terms
and
in the evolution equations of and
respectively. We introduced the variables
21Bona-Masso slicing condition (harmonic)
Initial data
with
22Lapse not regularized
23Metric not regularized
24Lapse regularized
25Metric regularized
26Lapse regularized
Lapse not regularized
27No hair conjecture
- Stationary black holes are completely specified
by the mass, charge and angular momentum. - Einstein-Yang-Mills
- Einstein-Skyrme
- Einstein-Yang-Mills-dilaton
- Einstein-Yang-Mills-Higgs
- Einstein-non-Abelian-Procca
Counterexamples
28- One place where it seemed for a while that there
was hope for a restricted form of the conjecture
was the scalar field arena. - No hair theorems
- Bekenstein ? Convex potentials
- Sudarsky ? Arbitrary positive potentials
29Torii, Maeda y Narita. PRD 64 044007 (2001)
Scalar fields minimally coupled to gravity in
spherically simmetric asymptotically anti-de
Sitter spacetimes
30U. Nucamendi y M. Salgado. PRD 68 044026 (2003)
31Asymptotically flat scalar hairy black holes
and solitons
Unstable under radial perturbations
32Numerical evolution
M.Alcubierre, J.A. González y M. Salgado (2004)
Gauge conditions
Collapse
Explosion
33Initial data (soliton)
Solve the constraint and the Klein-Gordon
equation in the static case a perturbation
in
Collapse
Explosion
34Evolution of the scalar field (Collapse)
35Evolution of the lapse function (Collapse)
36 Evolution of the metric function A (Collapse)
37Aparent horizon (Collapse)
38Convergence of the hamiltonian constraint
(Collapse)
39Evolution of the scalar field (Explosion)
40Evolution of the lapse function (Explosion)
41Evolution of the metric function A (Explosion)
42Convergence of the Hamiltonian constraint
(Explosion)
43Summary
- Spherically symmetric code
- Regularization
- Can handle with aparent horizons
- Numerical evolution of the scalar soliton
Collapse and Explosion