Evolution of 3D boson stars with waveform extraction

Numerical results from a study of boson stars under nonspherical perturbations using a fully general relativistic 3D code are presented together with the analysis of emitted gravitational radiation. We have constructed a simulation code suitable for the study of scalar fields in space-times of general symmetry by bringing together components for addressing the initial value problem, the full evolution system and the detection and analysis of gravitational waves. Within a series of numerical simulations, we explicitly extract the Zerilli and Newman-Penrose scalar Ψ₄ gravitational waveforms when the stars are subjected to different types of perturbations. Boson star systems have rapidly decaying nonradial quasinormal modes and thus the complete gravitational waveform could be extracted for all configurations studied. The gravitational waves emitted from stable, critical and unstable boson star configurations are analysed and the numerically observed quasinormal mode frequencies are compared with known linear perturbation results. The superposition of the high frequency nonspherical modes on the lower frequency spherical modes was observed in the metric oscillations when perturbations with radial and nonradial components were applied. The collapse of unstable boson stars to black holes was simulated. The apparent horizons were observed to be slightly nonspherical when initially detected and became spherical as the system evolved. The application of nonradial perturbations proportional to spherical harmonics is observed not to affect the collapse time. An unstable star subjected to a large perturbation was observed to migrate to a stable configuration.