As a first application of the code, Siebel, Font, Müller, and Papadopoulos  studied axisymmetric pulsations of neutron stars, which were initiated by perturbing the density and -component of velocity of a spherically-symmetric equilibrium configuration. The frequencies measured for the radial and non-radial oscillation modes of the star were found to be in good agreement with the results from linearized perturbation studies. The Bondi news function was computed and its amplitude found to be in rough agreement with the value given by the Einstein quadrupole formula. Both computations involve numerical subtleties: The computation of the news involves large terms, which partially cancel to give a small result, and the quadrupole formula requires computing three time derivatives of the fluid variables. These sources of computational error, coupled with ambiguity in the radiation content in the initial data, prevented any definitive conclusions. The total radiated mass loss was approximately 10–9 of the total mass.
Next, the code was applied to the simulation of axisymmetric supernova core collapse . A hybrid equation of state was used to mimic stiffening at collapse to nuclear densities and shock heating during the bounce. The initial equilibrium state of the core was modeled by a polytrope with index . Collapse was initiated by reducing the polytropic index to 1.3. In order to break spherical symmetry, small perturbations were introduced into the -component of the fluid velocity. During the collapse phase, the central density increased by five orders of magnitude. At this stage the inner core bounced at supra-nuclear densities, producing an expanding shock wave, which heated the outer layers. The collapse phase was well approximated by spherical symmetry but non-spherical oscillations were generated by the bounce. The resulting gravitational waves at null infinity were computed by the compactified code. After the bounce, the Bondi news function went through an oscillatory build up and then decayed in an quadrupole mode. However, a comparison with the results predicted by the Einstein quadrupole formula no longer gave the decent agreement found in the case of neutron star pulsations. This discrepancy was speculated to be due to the relativistic velocities of reached in the core collapse as opposed to for the pulsations. However, gauge effects and numerical errors also make important contributions, which cloud any definitive interpretation. This is the first study of gravitational wave production by the gravitational collapse of a relativistic star carried out with a characteristic code. It is clearly a remarkable piece of work, which offers up a whole new approach to the study of gravitational waves from astrophysical sources.
Living Rev. Relativity 15, (2012), 2
This work is licensed under a Creative Commons License.