skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Accurate simulations of the dynamical bar-mode instability in full general relativity

Abstract

We present accurate simulations of the dynamical bar-mode instability in full general relativity focusing on two aspects which have not been investigated in detail in the past, namely, on the persistence of the bar deformation once the instability has reached its saturation and on the precise determination of the threshold for the onset of the instability in terms of the parameter {beta}=T/|W|. We find that generic nonlinear mode-coupling effects appear during the development of the instability and these can severely limit the persistence of the bar deformation and eventually suppress the instability. In addition, we observe the dynamics of the instability to be strongly influenced by the value {beta} and on its separation from the critical value {beta}{sub c} marking the onset of the instability. We discuss the impact these results have on the detection of gravitational waves from this process and provide evidence that the classical perturbative analysis of the bar-mode instability for Newtonian and incompressible Maclaurin spheroids remains qualitatively valid and accurate also in full general relativity.

Authors:
 [1]; ;  [2];  [1];  [3];  [4]
  1. Max-Planck-Institut fuer Gravitationsphysik, Albert-Einstein-Institut, Golm (Germany)
  2. Dipartimento di Fisica, Universita di Parma and INFN, Parma (Italy)
  3. (Italy)
  4. (United States)
Publication Date:
OSTI Identifier:
21011088
Resource Type:
Journal Article
Resource Relation:
Journal Name: Physical Review. D, Particles Fields; Journal Volume: 75; Journal Issue: 4; Other Information: DOI: 10.1103/PhysRevD.75.044023; (c) 2007 The American Physical Society; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS; COMPUTERIZED SIMULATION; COSMOLOGY; COUPLING; DEFORMATION; GENERAL RELATIVITY THEORY; GRAVITATIONAL WAVES; NONLINEAR PROBLEMS

Citation Formats

Baiotti, Luca, Pietri, Roberto De, Manca, Gian Mario, Rezzolla, Luciano, SISSA, International School for Advanced Studies and INFN, Trieste, and Department of Physics, Louisiana State University, Baton Rouge. Accurate simulations of the dynamical bar-mode instability in full general relativity. United States: N. p., 2007. Web. doi:10.1103/PHYSREVD.75.044023.
Baiotti, Luca, Pietri, Roberto De, Manca, Gian Mario, Rezzolla, Luciano, SISSA, International School for Advanced Studies and INFN, Trieste, & Department of Physics, Louisiana State University, Baton Rouge. Accurate simulations of the dynamical bar-mode instability in full general relativity. United States. doi:10.1103/PHYSREVD.75.044023.
Baiotti, Luca, Pietri, Roberto De, Manca, Gian Mario, Rezzolla, Luciano, SISSA, International School for Advanced Studies and INFN, Trieste, and Department of Physics, Louisiana State University, Baton Rouge. Thu . "Accurate simulations of the dynamical bar-mode instability in full general relativity". United States. doi:10.1103/PHYSREVD.75.044023.
@article{osti_21011088,
title = {Accurate simulations of the dynamical bar-mode instability in full general relativity},
author = {Baiotti, Luca and Pietri, Roberto De and Manca, Gian Mario and Rezzolla, Luciano and SISSA, International School for Advanced Studies and INFN, Trieste and Department of Physics, Louisiana State University, Baton Rouge},
abstractNote = {We present accurate simulations of the dynamical bar-mode instability in full general relativity focusing on two aspects which have not been investigated in detail in the past, namely, on the persistence of the bar deformation once the instability has reached its saturation and on the precise determination of the threshold for the onset of the instability in terms of the parameter {beta}=T/|W|. We find that generic nonlinear mode-coupling effects appear during the development of the instability and these can severely limit the persistence of the bar deformation and eventually suppress the instability. In addition, we observe the dynamics of the instability to be strongly influenced by the value {beta} and on its separation from the critical value {beta}{sub c} marking the onset of the instability. We discuss the impact these results have on the detection of gravitational waves from this process and provide evidence that the classical perturbative analysis of the bar-mode instability for Newtonian and incompressible Maclaurin spheroids remains qualitatively valid and accurate also in full general relativity.},
doi = {10.1103/PHYSREVD.75.044023},
journal = {Physical Review. D, Particles Fields},
number = 4,
volume = 75,
place = {United States},
year = {Thu Feb 15 00:00:00 EST 2007},
month = {Thu Feb 15 00:00:00 EST 2007}
}
  • We perform fully general relativistic simulations of rotating stellar core collapse in three spatial dimensions. The hydrodynamic equations are solved using a high-resolution shock-capturing scheme. A parametric equation of state is adopted to model collapsing stellar cores and neutron stars following Dimmelmeier et al. The early stage of the collapse is followed by an axisymmetric code. When the stellar core becomes compact enough, we start a three-dimensional simulation adding a bar-mode nonaxisymmetric density perturbation. The axisymmetric simulations are performed for a wide variety of initial conditions changing the rotational velocity profile, parameters of the equations of state, and the totalmore » mass. It is clarified that the maximum density, the maximum value of the compactness, and the maximum value of the ratio of the kinetic energy T to the gravitational potential energy W ({beta}{identical_to}T/W) achieved during the stellar collapse and bounce depend sensitively on the velocity profile and the total mass of the initial core and equations of state. It is also found that for all the models with a high degree of differential rotation, a funnel structure is formed around the rotational axis after the formation of neutron stars. For selected models in which the maximum value of {beta} is larger than {approx}0.27, three-dimensional numerical simulations are performed. It is found that the bar-mode dynamical instability sets in for the case that the following conditions are satisfied: (i) the progenitor of the stellar core collapse should be rapidly rotating with the initial value of 0.01 < or approx. {beta} < or approx. 0.02 (ii) the degree of differential rotation for the velocity profile of the initial condition should be sufficiently high, and (iii) a depletion factor of pressure in an early stage of collapse should be large enough to induce a significant contraction to form a compact stellar core for which an efficient spin-up can be achieved surmounting the strong centrifugal force. As a result of the onset of the bar-mode dynamical instabilities, the amplitude of gravitational waves can be by a factor of {approx}10 larger than that in the axisymmetric collapse. It is found that a dynamical instability with the m=1 mode is also induced for the dynamically unstable cases against the bar mode, but the perturbation does not grow significantly and, hence, it does not contribute to an outstanding amplification of gravitational waves. No evidence for fragmentation of the protoneutron stars is found in the first few 10 msec after the bounce.« less
  • Numerical-relativity simulation is performed for rapidly spinning black holes (BHs) in a higher-dimensional spacetime of special symmetries for the dimensionality 6{<=}d{<=}8. We find that higher-dimensional BHs, spinning rapidly enough, are dynamically unstable against nonaxisymmetric bar-mode deformation and spontaneously emit gravitational waves, irrespective of d as in the case d=5[M. Shibata and H. Yoshino, Phys. Rev. D 81, 021501(R) (2010).]. The critical values of a nondimensional spin parameter for the onset of the instability are q:=a/{mu}{sup 1/(d-3){approx_equal}}0.74 for d=6, {approx_equal}0.73 for d=7, and {approx_equal}0.77 for d=8 where {mu} and a are mass and spin parameters. Black holes with a spin smallermore » than these critical values (q{sub crit}) appear to be dynamically stable for any perturbation. Long-term simulations for the unstable BHs are also performed for d=6 and 7. We find that they spin down as a result of gravitational-wave emission and subsequently settle to a stable stationary BH of a spin smaller than q{sub crit}. For more rapidly spinning unstable BHs, the time scale, for which the new state is reached, is shorter and fraction of the spin-down is larger. Our findings imply that a highly rapidly spinning BH with q>q{sub crit} cannot be a stationary product in the particle accelerators, even if it would be formed as a consequence of a TeV-gravity hypothesis. Its implications for the phenomenology of a mini BH are discussed.« less
  • We perform axisymmetric simulations for gravitational collapse of a massive iron core to a black hole in full general relativity. The iron cores are modeled by {gamma}=4/3 equilibrium polytrope for simplicity. The hydrodynamic equations are solved using a high-resolution shock-capturing scheme with a parametric equation of state. The Cartoon method is adopted for solving the Einstein equations. Simulations are performed for a wide variety of initial conditions changing the mass ({approx_equal}2.0-3.0M{sub {center_dot}}), the angular momentum, the rotational velocity profile of the core, and the parameters of the equations of state which are chosen so that the maximum mass of themore » cold spherical polytrope is {approx_equal}1.6M{sub {center_dot}}. Then, the criterion for the prompt black hole formation is clarified in terms of the mass and the angular momentum for several rotational velocity profile of the core and equations of state. It is found that (i) with the increase of the thermal energy generated by shocks, the threshold mass for the prompt black hole formation is increased by 20-40%, (ii) the rotational centrifugal force increases the threshold mass by < or approx. 25%, (iii) with the increase of the degree of differential rotation, the threshold mass is also increased, and (iv) the amplification factors shown in the results (i)-(iii) depend sensitively on the equation of state. We also find that the collapse dynamics and the structure of the shock formed at the bounce depend strongly on the stiffness of the adopted equation of state. In particular, as a new feature, a strong bipolar explosion is observed for the collapse of rapidly rotating iron cores with an equation of state which is stiff in subnuclear density and soft in supranuclear density. Gravitational waves are computed in terms of a quadrupole formula. It is also found that the waveform depends sensitively on the equations of state.« less
  • A recent paper by Lucas-Serrano et al. [A. Lucas-Serrano, J. A. Font, J. M. Ibanez, and J. M. Marti, Astron. Astrophys. 428, 703 (2004)] indicates that a high-resolution central (HRC) scheme is robust enough to yield accurate hydrodynamical simulations of special relativistic flows in the presence of ultrarelativistic speeds and strong shock waves. In this paper we apply this scheme in full general relativity (involving dynamical spacetimes), and assess its suitability by performing test simulations for oscillations of rapidly rotating neutron stars and merger of binary neutron stars. It is demonstrated that this HRC scheme can yield results as accuratemore » as those by the so-called high-resolution shock-capturing (HRSC) schemes based upon Riemann solvers. Furthermore, the adopted HRC scheme has increased computational efficiency as it avoids the costly solution of Riemann problems and has practical advantages in the modeling of neutron star spacetimes. Namely, it allows simulations with stiff equations of state by successfully dealing with very low-density unphysical atmospheres. These facts not only suggest that such a HRC scheme may be a desirable tool for hydrodynamical simulations in general relativity, but also open the possibility to perform accurate magnetohydrodynamical simulations in curved dynamic spacetimes.« less
  • We study the effects of magnetic fields on the evolution of differentially rotating neutron stars, which can be formed in stellar core collapse or binary neutron star coalescence. Magnetic braking and the magnetorotational instability (MRI) both act on differentially rotating stars to redistribute angular momentum. Simulations of these stars are carried out in axisymmetry using our recently developed codes which integrate the coupled Einstein-Maxwell-MHD equations. We consider stars with two different equations of state (EOS), a gamma-law EOS with {gamma}=2, and a more realistic hybrid EOS, and we evolve them adiabatically. Our simulations show that the fate of the starmore » depends on its mass and spin. For initial data, we consider three categories of differentially rotating, equilibrium configurations, which we label normal, hypermassive and ultraspinning. Normal configurations have rest masses below the maximum achievable with uniform rotation, and angular momentum below the maximum for uniform rotation at the same rest mass. Hypermassive stars have rest masses exceeding the mass limit for uniform rotation. Ultraspinning stars are not hypermassive, but have angular momentum exceeding the maximum for uniform rotation at the same rest mass. We show that a normal star will evolve to a uniformly rotating equilibrium configuration. An ultraspinning star evolves to an equilibrium state consisting of a nearly uniformly rotating central core, surrounded by a differentially rotating torus with constant angular velocity along magnetic field lines, so that differential rotation ceases to wind the magnetic field. In addition, the final state is stable against the MRI, although it has differential rotation. For a hypermassive neutron star, the MHD-driven angular momentum transport leads to catastrophic collapse of the core. The resulting rotating black hole is surrounded by a hot, massive, magnetized torus undergoing quasistationary accretion, and a magnetic field collimated along the spin axis--a promising candidate for the central engine of a short gamma-ray burst.« less