S E M I N A R IDepartament d'Astronomia i Astrofísica
Universitat de València
Stability of Iterative Algorithms for Rotating and Binary Neutron Stars
SFB/TR7 Gravitational Wave Astronomy, Theoretical Physics Institute, University of Jena, Alemanya
Seminari del departament d'Astronomia i AstrofÃsica
Dijous 31 maig 2012
Similar methods have been used to construct models of rapidly rotating or binary stars, in Newtonian and relativistic contexts. The choice of method has been based on numerical experiments, which indicate that particular methods converge quickly to a solution, while others diverge. The theory underlying these differences, however, has not been understood. In an attempt to provide a better theoretical understanding, we analytically examine the behavior of different iterative schemes near an exact static solution. We find the spectrum of the linearized iteration operator and show for self-consistent field methods that iterative instability corresponds to unstable modes of this operator. Minimizing the maximum eigenvalue accelerates convergence and allows computation of highly compact configurations that were previously inaccessible via self-consistent field methods.