A Caltech Library Service

Numerical simulations of oscillating soliton stars: Excited states in spherical symmetry and ground state evolutions in 3D

Balakrishna, Jayashree and Bondarescu, Ruxandra and Daues, Gregory and Bondarescu, Mihai (2008) Numerical simulations of oscillating soliton stars: Excited states in spherical symmetry and ground state evolutions in 3D. Physical Review D, 77 (2). Art. No. 024028. ISSN 2470-0010.

See Usage Policy.


Use this Persistent URL to link to this item:


Excited state soliton stars are studied numerically for the first time. The stability of spherically symmetric S-branch excited state oscillatons under radial perturbations is investigated using a 1D code. We find that these stars are inherently unstable either migrating to the ground state or collapsing to black holes. Higher excited state configurations are observed to cascade through intermediate excited states during their migration to the ground state. This is similar to excited state boson stars [J. Balakrishna, E. Seidel, and W.-M. Suen, Phys. Rev. D 58, 104004 (1998).]. Ground state oscillatons are then studied in full 3D numerical relativity. Finding the appropriate gauge condition for the dynamic oscillatons is much more challenging than in the case of boson stars. Different slicing conditions are explored, and a customized gauge condition that approximates polar slicing in spherical symmetry is implemented. Comparisons with 1D results and convergence tests are performed. The behavior of these stars under small axisymmetric perturbations is studied and gravitational waveforms are extracted. We find that the gravitational waves damp out on a short time scale, enabling us to obtain the complete waveform. This work is a starting point for the evolution of real scalar field systems with arbitrary symmetries.

Item Type:Article
Related URLs:
URLURL TypeDescription
Additional Information:©2008 The American Physical Society. (Received 17 October 2007; published 16 January 2008) We make extensive use of the Cactus Computational Toolkit, and its infrastructure for solving Einstein’s equations. The simulations have been performed on the NCSA Tungsten cluster under computer allocation TGMCA02N014. We use the bspline package from the GNU software library [33] for interpolation of initial data and the PETSc solver [34] for solving elliptic equations. We especially want to thank Francisco S. Guzman for advice and encouragement at the beginning of the project, without which this paper would not have been completed. We would also like to acknowledge Francisco for developing the 3D evolution routine for the scalar field. We thank Peter Diener for useful discussions and Ed Seidel and Wai-Mo Suen for allowing us to use their 1D scalar field evolution code. We are very grateful to Doina Costescu and Cornel Costescu for hospitality, support, and discussions during R.B.’s stays in Champaign-Urbana. We have received partial funding from NSF grants PHY-0652952 and AST-0606710, and the Sofja Kovalevskaja Program from the Alexander Von Humboldt Foundation.
Issue or Number:2
Record Number:CaltechAUTHORS:BALprd08
Persistent URL:
Usage Policy:No commercial reproduction, distribution, display or performance rights in this work are provided.
ID Code:9682
Deposited By: Archive Administrator
Deposited On:29 Feb 2008
Last Modified:03 Oct 2019 00:02

Repository Staff Only: item control page