Scalar field induced oscillations of neutron stars and gravitational collapse

We study the interaction of massless scalar fields with self-gravitating neutron stars by means of fully dynamic numerical simulations of the Einstein-Klein-Gordon perfect fluid system. Our investigation is restricted to spherical symmetry and the neutron stars are approximated by relativistic polytropes. Studying the nonlinear dynamics of isolated neutron stars is very effectively performed within the characteristic formulation of general relativity, in which the spacetime is foliated by a family of outgoing light cones. We are able to compactify the entire spacetime on a computational grid and simultaneously impose natural radiative boundary conditions and extract accurate radiative signals. We study the transfer of energy from the scalar field to the fluid star. We find, in particular, that depending on the compactness of the neutron star model, the scalar wave forces the neutron star either to oscillate in its radial modes of pulsation or to undergo gravitational collapse to a black hole on a dynamical timescale. The radiative signal, read off at future null infinity, shows quasi-normal oscillations before the setting of a late time power-law tail.Comment: 12 pages, 13 figures, submitted to Phys. Rev.

Nonlinear r-modes in Rapidly Rotating Relativistic Stars

The r-mode instability in rotating relativistic stars has been shown recently to have important astrophysical implications (including the emission of detectable gravitational radiation, the explanation of the initial spins of young neutron stars and the spin-distribution of millisecond pulsars and the explanation of one type of gamma-ray bursts), provided that r-modes are not saturated at low amplitudes by nonlinear effects or by dissipative mechanisms. Here, we present the first study of nonlinear r-modes in isentropic, rapidly rotating relativistic stars, via 3-D general-relativistic hydrodynamical evolutions. Our numerical simulations show that (1) on dynamical timescales, there is no strong nonlinear coupling of r-modes to other modes at amplitudes of order one -- unless nonlinear saturation occurs on longer timescales, the maximum r-mode amplitude is of order unity (i.e., the velocity perturbation is of the same order as the rotational velocity at the equator). An absolute upper limit on the amplitude (relevant, perhaps, for the most rapidly rotating stars) is set by causality. (2) r-modes and inertial modes in isentropic stars are predominantly discrete modes and possible associated continuous parts were not identified in our simulations. (3) In addition, the kinematical drift associated with r-modes, recently found by Rezzolla, Lamb and Shapiro (2000), appears to be present in our simulations, but an unambiguous confirmation requires more precise initial data. We discuss the implications of our findings for the detectability of gravitational waves from the r-mode instability.Comment: 4 pages, 4 eps figures, accepted in Physical Review Letter

The Post-Quasistatic Approximation as a test bed for Numerical Relativity

It is shown that observers in the standard ADM 3+1 treatment of matter are the same as the observers used in the matter treatment of Bondi: they are comoving and local Minkowskian. Bondi's observers are the basis of the post--quasitatic approximation (PQSA) to study a contracting distribution of matter. This correspondence suggests the possibility of using the PQSA as a test bed for Numerical Relativity. The treatment of matter by the PQSA and its connection with the ADM 3+1 treatment are presented, for its practical use as a calibration tool and as a test bed for numerical relativistic hydrodynamic codes.Comment: 4 pages; to appear as a Brief Report in Physical Review

Relativistic hydrodynamics on spacelike and null surfaces: Formalism and computations of spherically symmetric spacetimes

We introduce a formulation of Eulerian general relativistic hydrodynamics which is applicable for (perfect) fluid data prescribed on either spacelike or null hypersurfaces. Simple explicit expressions for the characteristic speeds and fields are derived in the general case. A complete implementation of the formalism is developed in the case of spherical symmetry. The algorithm is tested in a number of different situations, predisposing for a range of possible applications. We consider the Riemann problem for a polytropic gas, with initial data given on a retarded/advanced time slice of Minkowski spacetime. We compute perfect fluid accretion onto a Schwarzschild black hole spacetime using ingoing null Eddington-Finkelstein coordinates. Tests of fluid evolution on dynamic background include constant density and TOV stars sliced along the radial null cones. Finally, we consider the accretion of self-gravitating matter onto a central black hole and the ensuing increase in the mass of the black hole horizon.Comment: 23 pages, 13 figures, submitted to Phys. Rev.

Matter flows around black holes and gravitational radiation

We develop and calibrate a new method for estimating the gravitational radiation emitted by complex motions of matter sources in the vicinity of black holes. We compute numerically the linearized curvature perturbations induced by matter fields evolving in fixed black hole backgrounds, whose evolution we obtain using the equations of relativistic hydrodynamics. The current implementation of the proposal concerns non-rotating holes and axisymmetric hydrodynamical motions. As first applications we study i) dust shells falling onto the black hole isotropically from finite distance, ii) initially spherical layers of material falling onto a moving black hole, and iii) anisotropic collapse of shells. We focus on the dependence of the total gravitational wave energy emission on the flow parameters, in particular shell thickness, velocity and degree of anisotropy. The gradual excitation of the black hole quasi-normal mode frequency by sufficiently compact shells is demonstrated and discussed. A new prescription for generating physically reasonable initial data is discussed, along with a range of technical issues relevant to numerical relativity.Comment: 27 pages, 12 encapsulated figures, revtex, amsfonts, submitted to Phys. Rev.

Type II critical phenomena of neutron star collapse

We investigate spherically-symmetric, general relativistic systems of collapsing perfect fluid distributions. We consider neutron star models that are driven to collapse by the addition of an initially "in-going" velocity profile to the nominally static star solution. The neutron star models we use are Tolman-Oppenheimer-Volkoff solutions with an initially isentropic, gamma-law equation of state. The initial values of 1) the amplitude of the velocity profile, and 2) the central density of the star, span a parameter space, and we focus only on that region that gives rise to Type II critical behavior, wherein black holes of arbitrarily small mass can be formed. In contrast to previously published work, we find that--for a specific value of the adiabatic index (Gamma = 2)--the observed Type II critical solution has approximately the same scaling exponent as that calculated for an ultrarelativistic fluid of the same index. Further, we find that the critical solution computed using the ideal-gas equations of state asymptotes to the ultrarelativistic critical solution.Comment: 24 pages, 22 figures, RevTeX 4, submitted to Phys. Rev.

General Relativistic Radiant Shock Waves in the Post-Quasistatic Approximation

An evolution of radiant shock wave front is considered in the framework of a recently presented method to study self-gravitating relativistic spheres, whose rationale becomes intelligible and finds full justification within the context of a suitable definition of the post-quasistatic approximation. The spherical matter configuration is divided into two regions by the shock and each side of the interface having a different equation of state and anisotropic phase. In order to simulate dissipation effects due to the transfer of photons and/or neutrinos within the matter configuration, we introduce the flux factor, the variable Eddington factor and a closure relation between them. As we expected the strength of the shock increases the speed of the fluid to relativistic values and for some critical ones is larger than light speed. In addition, we find that energy conditions are very sensible to the anisotropy, specially the strong one. As a special feature of the model, we find that the contribution of the matter and radiation to the radial pressure are the same order of magnitude as in the mant as in the core, moreover, in the core radiation pressure is larger than matter pressure.Comment: To appear in Journal of Physics:Conference Series:"XXIX Spanish Relativity Meeting (ERE 2006): Einstein's Legacy: From the Theoretical Paradise to Astrophysical Observations

Relativistic Hydrodynamics around Black Holes and Horizon Adapted Coordinate Systems

Despite the fact that the Schwarzschild and Kerr solutions for the Einstein equations, when written in standard Schwarzschild and Boyer-Lindquist coordinates, present coordinate singularities, all numerical studies of accretion flows onto collapsed objects have been widely using them over the years. This approach introduces conceptual and practical complications in places where a smooth solution should be guaranteed, i.e., at the gravitational radius. In the present paper, we propose an alternative way of solving the general relativistic hydrodynamic equations in background (fixed) black hole spacetimes. We identify classes of coordinates in which the (possibly rotating) black hole metric is free of coordinate singularities at the horizon, independent of time, and admits a spacelike decomposition. In the spherically symmetric, non-rotating case, we re-derive exact solutions for dust and perfect fluid accretion in Eddington-Finkelstein coordinates, and compare with numerical hydrodynamic integrations. We perform representative axisymmetric computations. These demonstrations suggest that the use of those coordinate systems carries significant improvements over the standard approach, especially for higher dimensional studies.Comment: 10 pages, 4 postscript figures, accepted for publication in Phys. Rev.

Gravitational waves in dynamical spacetimes with matter content in the Fully Constrained Formulation

The Fully Constrained Formulation (FCF) of General Relativity is a novel framework introduced as an alternative to the hyperbolic formulations traditionally used in numerical relativity. The FCF equations form a hybrid elliptic-hyperbolic system of equations including explicitly the constraints. We present an implicit-explicit numerical algorithm to solve the hyperbolic part, whereas the elliptic sector shares the form and properties with the well known Conformally Flat Condition (CFC) approximation. We show the stability andconvergence properties of the numerical scheme with numerical simulations of vacuum solutions. We have performed the first numerical evolutions of the coupled system of hydrodynamics and Einstein equations within FCF. As a proof of principle of the viability of the formalism, we present 2D axisymmetric simulations of an oscillating neutron star. In order to simplify the analysis we have neglected the back-reaction of the gravitational waves into the dynamics, which is small (<2 %) for the system considered in this work. We use spherical coordinates grids which are well adapted for simulations of stars and allow for extended grids that marginally reach the wave zone. We have extracted the gravitational wave signature and compared to the Newtonian quadrupole and hexadecapole formulae. Both extraction methods show agreement within the numerical errors and the approximations used (~30 %).Comment: 17 pages, 9 figures, 2 tables, accepted for publication in PR

Solubility of nickel in slags equilibrated with Ni-S melt

To provide thermodynamic data for converting the nickel matte to liquid nickel, an experimental study was conducted in the phase equilibrium between the Ni-S melt and FeOX-SiO2, FeOX-CaO or CaO-Al2O3 based slag melted in a magnesia crucible at 1773 and 1873 K. pSO2 was controlled at 10.1 kPa while pO2 and pS2 ranged between those where NiO precipitated and Ni3S2 formed, respectively. The nickel content in the slag and the sulfur content in the metal at given pO2 and pS2 were smallest for the CaO-Al2O3 based slag. Both decreased with increasing temperature. At 1873 K, the content of nickel in the CaO-Al2O3 based slag at pO2 of 10 Pa (near the precipitation of NiO) was 4%, while the content of sulfur in alloy is 0.4 mass %. Thus, the CaO-Al2O3 base slag at 1873 K would be suitable for direct converting of Ni3S2 to metallic nickel. The distribution behavior of nickel between the slag and the Ni-S melt was discussed based on the concept of oxidic and sulfidic dissolution