Universality and properties of neutron star type I critical collapses

We study the neutron star axisymmetric critical solution previously found in the numerical studies of neutron star mergers. Using neutron star-like initial data and performing similar merger simulations, we demonstrate that the solution is indeed a semi-attractor on the threshold plane separating the basin of a neutron star and the basin of a black hole in the solution space of the Einstein equations. In order to explore the extent of the attraction basin of the neutron star semiattractor, we construct initial data phase spaces for these neutron star-like initial data. From these phase spaces, we also observe several interesting dynamical scenarios where the merged object is supported from prompt collapse. The properties of the critical index of the solution, in particular, its dependence on conserved quantities, are then studied. From the study, it is found that a family of neutron star semi-attractors exist that can be classified by both their rest masses and ADM masses.Comment: 13 pages, 12 figures, 1 new reference adde

Late-time evolution of nonlinear gravitational collapse

We study numerically the fully nonlinear gravitational collapse of a self-gravitating, minimally-coupled, massless scalar field in spherical symmetry. Our numerical code is based on double-null coordinates and on free evolution of the metric functions: The evolution equations are integrated numerically, whereas the constraint equations are only monitored. The numerical code is stable (unlike recent claims) and second-order accurate. We use this code to study the late-time asymptotic behavior at fixed $r$ (outside the black hole), along the event horizon, and along future null infinity. In all three asymptotic regions we find that, after the decay of the quasi-normal modes, the perturbations are dominated by inverse power-law tails. The corresponding power indices agree with the integer values predicted by linearized theory. We also study the case of a charged black hole nonlinearly perturbed by a (neutral) self-gravitating scalar field, and find the same type of behavior---i.e., quasi-normal modes followed by inverse power-law tails, with the same indices as in the uncharged case.Comment: 14 pages, standard LaTeX, 18 Encapsulated PostScript figures. A new convergence test and a determination of QN ringing were added, in addition to correction of typos and update of reference

Critical Collapse of a Complex Scalar Field with Angular Momentum

We report a new critical solution found at the threshold of axisymmetric gravitational collapse of a complex scalar field with angular momentum. To carry angular momentum the scalar field cannot be axisymmetric; however, its azimuthal dependence is defined so that the resulting stress energy tensor and spacetime metric are axisymmetric. The critical solution found is non-spherical, discretely self-similar with an echoing exponent of 0.42 (+- 4%), and exhibits a scaling exponent of 0.11 (+- 10%) in near critical collapse. Our simulations suggest that the solution is universal (within the imposed symmetry class), modulo a family-dependent constant phase in the complex plane.Comment: 4 pages, 4 figure

Critical Collapse of the Massless Scalar Field in Axisymmetry

We present results from a numerical study of critical gravitational collapse of axisymmetric distributions of massless scalar field energy. We find threshold behavior that can be described by the spherically symmetric critical solution with axisymmetric perturbations. However, we see indications of a growing, non-spherical mode about the spherically symmetric critical solution. The effect of this instability is that the small asymmetry present in what would otherwise be a spherically symmetric self-similar solution grows. This growth continues until a bifurcation occurs and two distinct regions form on the axis, each resembling the spherically symmetric self-similar solution. The existence of a non-spherical unstable mode is in conflict with previous perturbative results, and we therefore discuss whether such a mode exists in the continuum limit, or whether we are instead seeing a marginally stable mode that is rendered unstable by numerical approximation.Comment: 11 pages, 8 figure

On free evolution of self gravitating, spherically symmetric waves

We perform a numerical free evolution of a selfgravitating, spherically symmetric scalar field satisfying the wave equation. The evolution equations can be written in a very simple form and are symmetric hyperbolic in Eddington-Finkelstein coordinates. The simplicity of the system allow to display and deal with the typical gauge instability present in these coordinates. The numerical evolution is performed with a standard method of lines fourth order in space and time. The time algorithm is Runge-Kutta while the space discrete derivative is symmetric (non-dissipative). The constraints are preserved under evolution (within numerical errors) and we are able to reproduce several known results.Comment: 15 pages, 15 figure

Asymptotic tails of massive scalar fields in Schwarzschild background

We investigate the asymptotic tail behavior of massive scalar fields in Schwarzschild background. It is shown that the oscillatory tail of the scalar field has the decay rate of $t^{-5/6}$ at asymptotically late times, and the oscillation with the period $2\pi/m$ for the field mass $m$ is modulated by the long-term phase shift. These behaviors are qualitatively similar to those found in nearly extreme Reissner-Nordstr\"{o}m background, which are discussed in terms of a resonant backscattering due to the space-time curvature.Comment: 21 pages, 2 figures, accepted for publication in Phys.Rev.

The footprint of cometary dust analogues: II. Morphology as a tracer of tensile strength and application to dust collection by the Rosetta spacecraft

The structure of cometary dust is a tracer of growth processes in the formation of planetesimals. Instrumentation on board the Rosetta mission to comet 67P/Churyumov- Gerasimenko captured dust particles and analysed them in situ. However, these deposits are a product of a collision within the instrument. We conducted laboratory experiments with cometary dust analogues, simulating the collection process by Rosetta instruments (specifically COSIMA, MIDAS). In Paper I we reported that velocity is a key driver in determining the appearance of deposits. Here in Paper II we use materials with different monomer sizes, and study the effect of tensile strength on the appearance of deposits. We find that mass transfer efficiency increases from $\sim$1 up to $\sim$10% with increasing monomer diameter from 0.3 $\mu$m to 1.5 $\mu$m (i.e. tensile strength decreasing from $\sim$12 to $\sim$3 kPa), and velocities increasing from 0.5 to 6 m/s. Also, the relative abundance of small fragments after impact is higher for material with higher tensile strength. The degeneracy between the effects of velocity and material strength may be lifted by performing a closer study of the deposits. This experimental method makes it possible to estimate the mass transfer efficiency in the COSIMA instrument. Extrapolating these results implies that more than half of the dust collected during the Rosetta mission has not been imaged. We analysed two COSIMA targets containing deposits from single collisions. The collision that occurred closest to perihelion passage led to more small fragments on the target.Comment: 13 pages, 11 figures, accepted for publication in MNRA

Black Hole--Scalar Field Interactions in Spherical Symmetry

We examine the interactions of a black hole with a massless scalar field using a coordinate system which extends ingoing Eddington-Finkelstein coordinates to dynamic spherically symmetric-spacetimes. We avoid problems with the singularity by excising the region of the black hole interior to the apparent horizon. We use a second-order finite difference scheme to solve the equations. The resulting program is stable and convergent and will run forever without problems. We are able to observe quasi-normal ringing and power-law tails as well an interesting nonlinear feature.Comment: 16 pages, 26 figures, RevTex, to appear in Phys. Rev.

Evolution equations for slowly rotating stars

We present a hyperbolic formulation of the evolution equations describing non-radial perturbations of slowly rotating relativistic stars in the Regge--Wheeler gauge. We demonstrate the stability preperties of the new evolution set of equations and compute the polar w-modes for slowly rotating stars.Comment: 27 pages, 2 figure