Normalized to: Stergioulas, N.
[1]
oai:arXiv.org:2005.14135 [pdf] - 2103546
Astrophysical Implications of Neutron Star Inspiral and Coalescence
Submitted: 2020-05-28
The first inspiral of two neutron stars observed in gravitational waves was
remarkably close, allowing the kind of simultaneous gravitational wave and
electromagnetic observation that had not been expected for several years. Their
merger, followed by a gamma-ray burst and a kilonova, was observed across the
spectral bands of electromagnetic telescopes. These GW and electromagnetic
observations have led to dramatic advances in understanding short gamma-ray
bursts; determining the origin of the heaviest elements; and determining the
maximum mass of neutron stars. From the imprint of tides on the gravitational
waveforms and from observations of X-ray binaries, one can extract the radius
and deformability of inspiraling neutron stars. Together, the radius, maximum
mass, and causality constrain the neutron-star equation of state, and future
constraints can come from observations of post-merger oscillations. We
selectively review these results, filling in some of the physics with
derivations and estimates.
[2]
oai:arXiv.org:2004.00846 [pdf] - 2074532
Equation of state constraints from the threshold binary mass for prompt
collapse of neutron star mergers
Bauswein, Andreas;
Blacker, Sebastian;
Vijayan, Vimal;
Stergioulas, Nikolaos;
Chatziioannou, Katerina;
Clark, James A.;
Bastian, Niels-Uwe F.;
Blaschke, David B.;
Cierniak, Mateusz;
Fischer, Tobias
Submitted: 2020-04-02
Using hydrodynamical simulations for a large set of high-density matter
equations of state (EoSs) we systematically determine the threshold mass
M_thres for prompt black-hole formation in equal-mass and asymmetric neutron
star (NS) mergers. We devise the so far most direct, general and accurate
method to determine the unknown maximum mass of nonrotating NSs from merger
observations revealing M_thres. Considering hybrid EoSs with hadron-quark phase
transition, we identify a new, observable signature of quark matter in NS
mergers. Furthermore, our findings have direct applications in gravitational
wave searches, kilonova interpretations and multi-messenger constraints on NS
properties.
[3]
oai:arXiv.org:1910.04036 [pdf] - 2072571
Numerical-relativity simulations of long-lived remnants of binary
neutron star mergers
Submitted: 2019-10-09, last modified: 2020-03-21
We analyze the properties of the gravitational wave signal emitted after the
merger of a binary neutron star system when the remnant survives for more than
a 80 ms (and up to 140ms). We employ four different piecewise polytropic
equations of state supplemented by an ideal fluid thermal component. We find
that the post-merger phase can be subdivided into three phases: an early
post-merger phase (where the quadrupole mode and a few subdominant features are
active), the intermediate post-merger phase (where only the quadrupole mode is
active) and the late post-merger phase (where convective instabilities trigger
inertial modes). The inertial modes have frequencies somewhat smaller than the
quadrupole modes. In one model, we find an interesting association of a
corotation of the quadrupole mode in parts of the star with a revival of its
amplitude. The gravitational wave emission of inertial modes in the late
post-merger phase is concentrated in a narrow frequency region and is
potentially detectable by the planned third-generation detectors. This allows
for the possibility of probing not only the cold part of the equation of state,
but also its dependence on finite temperature. In view of these results, it
will be important to investigate the impact of various type of viscosities on
the potential excitation of inertial modes in binary neutron star merger
remnants.
[4]
oai:arXiv.org:1910.10856 [pdf] - 2085050
Empirical relations for gravitational-wave asteroseismology of binary
neutron star mergers
Submitted: 2019-10-23
We construct new, multivariate empirical relations for measuring neutron star
radii and tidal deformabilities from the dominant gravitational wave frequency
in the post-merger phase of binary neutron star mergers. The relations
determine neutron star radii and tidal deformabilities for specific neutron
star masses with consistent accuracy and depend only on two observables: the
post-merger peak frequency $f_{\rm peak}$ and the chirp mass $M_{\rm chirp}$.
The former could be measured with good accuracy from gravitational waves
emitted in the post-merger phase using next-generation detectors, whereas the
latter is already obtained with good accuracy from the inspiral phase with
present-day detectors. Our main data set consists of a gravitational wave
catalogue obtained with CFC/SPH simulations. We also extract the $f_{\rm peak}$
frequency from the publicly available CoRe data set, obtained through
grid-based GRHD simulations and find good agreement between the extracted
frequencies of the two data sets. As a result, we can construct empirical
relations for the combined data sets. Furthermore, we investigate empirical
relations for two secondary peaks, $f_{2-0}$ and $f_{\rm spiral}$, and show
that these relations are distinct in the whole parameter space, in agreement
with a previously introduced spectral classification scheme. Finally, we show
that the spectral classification scheme can be reproduced using
machine-learning techniques.
[5]
oai:arXiv.org:1908.11390 [pdf] - 1951580
Probing the Nature of Black Holes: Deep in the mHz Gravitational-Wave
Sky
Baibhav, Vishal;
Barack, Leor;
Berti, Emanuele;
Bonga, Béatrice;
Brito, Richard;
Cardoso, Vitor;
Compère, Geoffrey;
Das, Saurya;
Doneva, Daniela;
Garcia-Bellido, Juan;
Heisenberg, Lavinia;
Hughes, Scott A.;
Isi, Maximiliano;
Jani, Karan;
Kavanagh, Chris;
Lukes-Gerakopoulos, Georgios;
Mueller, Guido;
Pani, Paolo;
Petiteau, Antoine;
Rajendran, Surjeet;
Sotiriou, Thomas P.;
Stergioulas, Nikolaos;
Taylor, Alasdair;
Vagenas, Elias;
van de Meent, Maarten;
Warburton, Niels;
Wardell, Barry;
Witzany, Vojtěch;
Zimmerman, Aaron
Submitted: 2019-08-29
Black holes are unique among astrophysical sources: they are the simplest
macroscopic objects in the Universe, and they are extraordinary in terms of
their ability to convert energy into electromagnetic and gravitational
radiation. Our capacity to probe their nature is limited by the sensitivity of
our detectors. The LIGO/Virgo interferometers are the gravitational-wave
equivalent of Galileo's telescope. The first few detections represent the
beginning of a long journey of exploration. At the current pace of
technological progress, it is reasonable to expect that the gravitational-wave
detectors available in the 2035-2050s will be formidable tools to explore these
fascinating objects in the cosmos, and space-based detectors with peak
sensitivities in the mHz band represent one class of such tools. These
detectors have a staggering discovery potential, and they will address
fundamental open questions in physics and astronomy. Are astrophysical black
holes adequately described by general relativity? Do we have empirical evidence
for event horizons? Can black holes provide a glimpse into quantum gravity, or
reveal a classical breakdown of Einstein's gravity? How and when did black
holes form, and how do they grow? Are there new long-range interactions or
fields in our universe, potentially related to dark matter and dark energy or a
more fundamental description of gravitation? Precision tests of black hole
spacetimes with mHz-band gravitational-wave detectors will probe general
relativity and fundamental physics in previously inaccessible regimes, and
allow us to address some of these fundamental issues in our current
understanding of nature.
[6]
oai:arXiv.org:1904.01306 [pdf] - 1953461
Equation-of-state Constraints and the QCD Phase Transition in the Era of
Gravitational-Wave Astronomy
Bauswein, Andreas;
Bastian, Niels-Uwe Friedrich;
Blaschke, David;
Chatziioannou, Katerina;
Clark, James Alexander;
Fischer, Tobias;
Janka, Hans-Thomas;
Just, Oliver;
Oertel, Micaela;
Stergioulas, Nikolaos
Submitted: 2019-04-02
We describe a multi-messenger interpretation of GW170817, which yields a
robust lower limit on NS radii. This excludes NSs with radii smaller than about
10.7 km and thus rules out very soft nuclear matter. We stress the potential of
this type of constraints when future detections become available. A very
similar argumentation may yield an upper bound on the maximum mass of
nonrotating NSs. We also discuss simulations of NS mergers, which undergo a
first-order phase transition to quark matter. We point out a different
dynamical behavior. Considering the gravitational-wave signal, we identify an
unambiguous signature of the QCD phase transition in NS mergers. The occurrence
of quark matter through a strong first-order phase transition during merging
leads to a characteristic shift of the dominant postmerger frequency. The
frequency shift is indicative for a phase transition if it is compared to the
postmerger frequency which is expected for purely hadronic EoS models. A very
strong deviation of several 100 Hz is observed for hybrid EoSs in an otherwise
tight relation between the tidal deformability and the postmerger frequency. We
address the potential impact of a first-order phase transition on the
electromagnetic counterpart of NS mergers. Our simulations suggest that there
would be no significant qualitative differences between a system undergoing a
phase transition to quark matter and purely hadronic mergers. The quantitative
differences are within the spread which is found between different hadronic EoS
models. This implies on the one hand that GW170817 is compatible with a
possible transition to quark matter. On the other hand these considerations
show that it may not be easy to identify quantitative differences between
purely hadronic mergers and events in which quark matter occurs considering
solely their electromagnetic counterpart or their nucleosynthesis products.
(abridged)
[7]
oai:arXiv.org:1903.05519 [pdf] - 1979546
Rotating Black Holes in Higher Order Gravity
Submitted: 2019-03-13
We develop a new technique for finding black hole solutions in modified
gravity that have "stealth" hair, i.e., hair whose only gravitational effect is
to tune the cosmological constant. We consider scalar-tensor theories in which
gravitational waves propagate at the speed of light, and show that Einstein
metrics can be painted with stealth hair provided there exists a family of
geodesics always normal to spacelike surfaces. We also present a novel
scalar-dressed rotating black hole that has finite scalar field at both the
black hole and cosmological event horizons.
[8]
oai:arXiv.org:1806.05195 [pdf] - 1916735
Black holes, gravitational waves and fundamental physics: a roadmap
Barack, Leor;
Cardoso, Vitor;
Nissanke, Samaya;
Sotiriou, Thomas P.;
Askar, Abbas;
Belczynski, Krzysztof;
Bertone, Gianfranco;
Bon, Edi;
Blas, Diego;
Brito, Richard;
Bulik, Tomasz;
Burrage, Clare;
Byrnes, Christian T.;
Caprini, Chiara;
Chernyakova, Masha;
Chrusciel, Piotr;
Colpi, Monica;
Ferrari, Valeria;
Gaggero, Daniele;
Gair, Jonathan;
Garcia-Bellido, Juan;
Hassan, S. F.;
Heisenberg, Lavinia;
Hendry, Martin;
Heng, Ik Siong;
Herdeiro, Carlos;
Hinderer, Tanja;
Horesh, Assaf;
Kavanagh, Bradley J.;
Kocsis, Bence;
Kramer, Michael;
Tiec, Alexandre Le;
Mingarelli, Chiara;
Nardini, Germano;
Nelemans, Gijs;
Palenzuela, Carlos;
Pani, Paolo;
Perego, Albino;
Porter, Edward K.;
Rossi, Elena M.;
Schmidt, Patricia;
Sesana, Alberto;
Sperhake, Ulrich;
Stamerra, Antonio;
Stein, Leo C.;
Tamanini, Nicola;
Tauris, Thomas M.;
Urena-Lopez, L. Arturo;
Vincent, Frederic;
Volonteri, Marta;
Wardell, Barry;
Wex, Norbert;
Yagi, Kent;
Abdelsalhin, Tiziano;
Aloy, Miguel Angel;
Amaro-Seoane, Pau;
Annulli, Lorenzo;
Arca-Sedda, Manuel;
Bah, Ibrahima;
Barausse, Enrico;
Barakovic, Elvis;
Benkel, Robert;
Bennett, Charles L.;
Bernard, Laura;
Bernuzzi, Sebastiano;
Berry, Christopher P. L.;
Berti, Emanuele;
Bezares, Miguel;
Blanco-Pillado, Jose Juan;
Blazquez-Salcedo, Jose Luis;
Bonetti, Matteo;
Boskovic, Mateja;
Bosnjak, Zeljka;
Bricman, Katja;
Bruegmann, Bernd;
Capelo, Pedro R.;
Carloni, Sante;
Cerda-Duran, Pablo;
Charmousis, Christos;
Chaty, Sylvain;
Clerici, Aurora;
Coates, Andrew;
Colleoni, Marta;
Collodel, Lucas G.;
Compere, Geoffrey;
Cook, William;
Cordero-Carrion, Isabel;
Correia, Miguel;
de la Cruz-Dombriz, Alvaro;
Czinner, Viktor G.;
Destounis, Kyriakos;
Dialektopoulos, Kostas;
Doneva, Daniela;
Dotti, Massimo;
Drew, Amelia;
Eckner, Christopher;
Edholm, James;
Emparan, Roberto;
Erdem, Recai;
Ferreira, Miguel;
Ferreira, Pedro G.;
Finch, Andrew;
Font, Jose A.;
Franchini, Nicola;
Fransen, Kwinten;
Gal'tsov, Dmitry;
Ganguly, Apratim;
Gerosa, Davide;
Glampedakis, Kostas;
Gomboc, Andreja;
Goobar, Ariel;
Gualtieri, Leonardo;
Guendelman, Eduardo;
Haardt, Francesco;
Harmark, Troels;
Hejda, Filip;
Hertog, Thomas;
Hopper, Seth;
Husa, Sascha;
Ihanec, Nada;
Ikeda, Taishi;
Jaodand, Amruta;
Jimenez-Forteza, Philippe Jetzer Xisco;
Kamionkowski, Marc;
Kaplan, David E.;
Kazantzidis, Stelios;
Kimura, Masashi;
Kobayashi, Shiho;
Kokkotas, Kostas;
Krolik, Julian;
Kunz, Jutta;
Lammerzahl, Claus;
Lasky, Paul;
Lemos, Jose P. S.;
Said, Jackson Levi;
Liberati, Stefano;
Lopes, Jorge;
Luna, Raimon;
Ma, Yin-Zhe;
Maggio, Elisa;
Montero, Marina Martinez;
Maselli, Andrea;
Mayer, Lucio;
Mazumdar, Anupam;
Messenger, Christopher;
Menard, Brice;
Minamitsuji, Masato;
Moore, Christopher J.;
Mota, David;
Nampalliwar, Sourabh;
Nerozzi, Andrea;
Nichols, David;
Nissimov, Emil;
Obergaulinger, Martin;
Obers, Niels A.;
Oliveri, Roberto;
Pappas, George;
Pasic, Vedad;
Peiris, Hiranya;
Petrushevska, Tanja;
Pollney, Denis;
Pratten, Geraint;
Rakic, Nemanja;
Racz, Istvan;
Radia, Miren;
Ramazanouglu, Fethi M.;
Ramos-Buades, Antoni;
Raposo, Guilherme;
Rosca-Mead, Roxana;
Rogatko, Marek;
Rosinska, Dorota;
Rosswog, Stephan;
Morales, Ester Ruiz;
Sakellariadou, Mairi;
Sanchis-Gual, Nicolas;
Salafia, Om Sharan;
Samajdar, Anuradha;
Sintes, Alicia;
Smole, Majda;
Sopuerta, Carlos;
Souza-Lima, Rafael;
Stalevski, Marko;
Stergioulas, Nikolaos;
Stevens, Chris;
Tamfal, Tomas;
Torres-Forne, Alejandro;
Tsygankov, Sergey;
Unluturk, Kivanc;
Valiante, Rosa;
van de Meent, Maarten;
Velhinho, Jose;
Verbin, Yosef;
Vercnocke, Bert;
Vernieri, Daniele;
Vicente, Rodrigo;
Vitagliano, Vincenzo;
Weltman, Amanda;
Whiting, Bernard;
Williamson, Andrew;
Witek, Helvi;
Wojnar, Aneta;
Yakut, Kadri;
Yan, Haopeng;
Yazadjiev, Stoycho;
Zaharijas, Gabrijela;
Zilhao, Miguel
Submitted: 2018-06-13, last modified: 2019-02-01
The grand challenges of contemporary fundamental physics---dark matter, dark
energy, vacuum energy, inflation and early universe cosmology, singularities
and the hierarchy problem---all involve gravity as a key component. And of all
gravitational phenomena, black holes stand out in their elegant simplicity,
while harbouring some of the most remarkable predictions of General Relativity:
event horizons, singularities and ergoregions. The hitherto invisible landscape
of the gravitational Universe is being unveiled before our eyes: the historical
direct detection of gravitational waves by the LIGO-Virgo collaboration marks
the dawn of a new era of scientific exploration. Gravitational-wave astronomy
will allow us to test models of black hole formation, growth and evolution, as
well as models of gravitational-wave generation and propagation. It will
provide evidence for event horizons and ergoregions, test the theory of General
Relativity itself, and may reveal the existence of new fundamental fields. The
synthesis of these results has the potential to radically reshape our
understanding of the cosmos and of the laws of Nature. The purpose of this work
is to present a concise, yet comprehensive overview of the state of the art in
the relevant fields of research, summarize important open problems, and lay out
a roadmap for future progress.
[9]
oai:arXiv.org:1807.05449 [pdf] - 1790835
Differentially rotating neutron stars in scalar-tensor theories of
gravity
Submitted: 2018-07-14
We present the first numerical models of differentially rotating stars in
alternative theories of gravity. We chose a particular class of scalar-tensor
theories of gravity that is indistinguishable from GR in the weak field regime
but can lead to significant deviations when strong fields are considered. We
show that the maximum mass that a differentially rotating neutron star can
sustain increases significantly for scalarized solutions and such stars can
reach larger angular momenta. In addition, the presence of a nontrivial scalar
field has the effect of increasing the required axis ratio for reaching a given
value of angular momentum, when compared to a corresponding model of same rest
mass in general relativity. We find that the scalar field also makes rapidly
rotating models less quasi-toroidal than their general-relativistic
counterparts. For large values of the angular momentum and values of the
coupling parameter that are in agreement with the observations, we find a
second turning point for scalarized models along constant angular momentum
sequences, which could have interesting implications for the stability of
remnants created in a binary neutron star merger.
[10]
oai:arXiv.org:1802.03288 [pdf] - 1709342
Convective Excitation of Inertial Modes in Binary Neutron Star Mergers
Submitted: 2018-02-09, last modified: 2018-07-03
We present the first very long-term simulations (extending up to ~140 ms
after merger) of binary neutron star mergers with piecewise polytropic
equations of state and in full general relativity. Our simulations reveal that
at a time of 30-50 ms after merger, parts of the star become convectively
unstable, which triggers the excitation of inertial modes. The excited inertial
modes are sustained up to several tens of milliseconds and are potentially
observable by the planned third-generation gravitational-wave detectors at
frequencies of a few kilohertz. Since inertial modes depend on the rotation
rate of the star and they are triggered by a convective instability in the
postmerger remnant, their detection in gravitational waves will provide a
unique opportunity to probe the rotational and thermal state of the merger
remnant. In addition, our findings have implications for the long-term
evolution and stability of binary neutron star remnants
[11]
oai:arXiv.org:1709.10067 [pdf] - 1614829
Universal and approximate relations for the gravitational-wave damping
timescale of $f$-modes in neutron stars
Submitted: 2017-09-28, last modified: 2018-01-09
Existing estimates of the gravitational-wave damping timescale of the
dominant quadrupole oscillation mode in the case of rapidly rotating stars are
based on using a Newtonian estimate for the energy of the mode, in combination
with the lowest-order post-Newtonian quadrupole formula for estimating the
gravitational-wave luminosity. We investigate a number of other choices for
estimating the gravitational-wave damping timescale in the nonrotating limit
and construct a highly accurate, empirically corrected formula that has a
maximum relative error of only 3% with respect to the perturbative result in
full general relativity. The expressions involved are sufficiently general to
be extended to the case of rapidly rotating stars. We also present a new
higher-order empirical relation for the gravitational-wave damping timescale of
quadrupole oscillations that is accurate in the whole range of expected values
for the compactness of neutron stars, without the need for involving the moment
of inertia.
[12]
oai:arXiv.org:1612.03050 [pdf] - 1595706
Rotating Stars in Relativity
Submitted: 2016-12-09, last modified: 2017-11-30
Rotating relativistic stars have been studied extensively in recent years,
both theoretically and observationally, because of the information they might
yield about the equation of state of matter at extremely high densities and
because they are considered to be promising sources of gravitational waves. The
latest theoretical understanding of rotating stars in relativity is reviewed in
this updated article. The sections on equilibrium properties and on
nonaxisymmetric oscillations and instabilities in $f$-modes and $r$-modes have
been updated. Several new sections have been added on equilibria in modified
theories of gravity, approximate universal relationships, the one-arm spiral
instability, on analytic solutions for the exterior spacetime, rotating stars
in LMXBs, rotating strange stars, and on rotating stars in numerical relativity
including both hydrodynamic and magnetohydrodynamic studies of these objects.
[13]
oai:arXiv.org:1710.06843 [pdf] - 1602710
Neutron-star radius constraints from GW170817 and future detections
Submitted: 2017-10-18, last modified: 2017-11-09
We introduce a new, powerful method to constrain properties of neutron stars
(NSs). We show that the total mass of GW170817 provides a reliable constraint
on the stellar radius if the merger did not result in a prompt collapse as
suggested by the interpretation of associated electromagnetic emission. The
radius R_1.6 of nonrotating NSs with a mass of 1.6 M_sun can be constrained to
be larger than 10.68_{-0.04}^{+0.15} km, and the radius R_max of the
nonrotating maximum mass configuration must be larger than 9.60_{-0.03}^{+0.14}
km. We point out that detections of future events will further improve these
constraints. Moreover, we show that a future event with a signature of a prompt
collapse of the merger remnant will establish even stronger constraints on the
NS radius from above and the maximum mass M_max of NSs from above. These
constraints are particularly robust because they only require a measurement of
the chirp mass and a distinction between prompt and delayed collapse of the
merger remnant, which may be inferred from the electromagnetic signal or even
from the presence/absence of a ringdown gravitational-wave (GW) signal. This
prospect strengthens the case of our novel method of constraining NS
properties, which is directly applicable to future GW events with accompanying
electromagnetic counterpart observations. We emphasize that this procedure is a
new way of constraining NS radii from GW detections independent of existing
efforts to infer radius information from the late inspiral phase or postmerger
oscillations, and it does not require particularly loud GW events.
[14]
oai:arXiv.org:1710.02334 [pdf] - 1652340
Constraining properties of high-density matter in neutron stars with
magneto-elastic oscillations
Submitted: 2017-10-06
We discuss torsional oscillations of highly magnetised neutron stars
(magnetars) using two-dimensional, magneto-elastic-hydrodynamical simulations.
Our model is able to explain both the low- and high-frequency quasi-periodic
oscillations (QPOs) observed in magnetars. The analysis of these oscillations
provides constraints on the breakout magnetic-field strength, on the
fundamental QPO frequency, and on the frequency of a particularly excited
overtone. More importantly, we show how to use this information to generically
constraint properties of high-density matter in neutron stars, employing
Bayesian analysis. In spite of current uncertainties and computational
approximations, our model-dependent Bayesian posterior estimates for SGR
1806-20 yield a magnetic-field strength $\bar B\sim
2.1^{+1.3}_{-1.0}\times10^{15}\,$G and a crust thickness of $\Delta r =
1.6^{+0.7}_{-0.6}$ km, which are both in remarkable agreement with
observational and theoretical expectations, respectively (1-$\sigma$ error bars
are indicated). Our posteriors also favour the presence of a superfluid phase
in the core, a relatively low stellar compactness, $M/R<0.19$, indicating a
relatively stiff equation of state and/or low mass neutron star, and high shear
speeds at the base of the crust, $c_s>1.4\times10^8\,$cm/s. Although the
procedure laid out here still has large uncertainties, these constraints could
become tighter when additional observations become available.
[15]
oai:arXiv.org:1709.06806 [pdf] - 1614819
Magneto-elastic oscillations modulating the emission of magnetars
Submitted: 2017-09-20
Magneto-elastic oscillations of neutron stars are believed to explain
observed quasi-periodic oscillations (QPOs) in the decaying tail of the giant
flares of highly magnetized neutron stars (magnetars). Strong efforts of the
theoretical modelling from different groups have increased our understanding of
this phenomenon significantly. Here, we discuss some constraints on the matter
in neutron stars that arise if the interpretation of the observations in terms
of superfluid, magneto-elastic oscillations is correct. To explain the observed
modulation of the light curve of the giant flare, we describe a model that
allows the QPOs to couple to the stellar exterior through the magnetic field.
In this magnetosphere, the shaking magnetic field induces currents that provide
scattering targets for resonant cyclotron scattering of photons, which is
calculated with a Monte-Carlo approach and coupled to a code that calculates
the momentum distribution of the charge carriers as a one-dimensional
accelerator problem. We show first results of a simplified, but self-consistent
momentum distribution, i.e. a waterbag distribution, and of the corresponding
spectra.
[16]
oai:arXiv.org:1709.02787 [pdf] - 1608446
Universal relations for differentially rotating relativistic stars at
the threshold to collapse
Submitted: 2017-09-08
A binary neutron star merger produces a rapidly and differentially rotating
compact remnant whose lifespan heavily affects the electromagnetic and
gravitational emissions. Its stability depends on both the equation of state
(EOS) and the rotation law and it is usually investigated through numerical
simulations. Nevertheless, by means of a sufficient criterion for secular
instability, equilibrium sequences can be used as a computationally inexpensive
way to estimate the onset of dynamical instability, which, in general, is close
to the secular one. This method works well for uniform rotation and relies on
the location of turning points: stellar models that are stationary points in a
sequence of equilibrium solutions with constant rest mass or angular momentum.
Here, we investigate differentially rotating models (using a large number of
equations of state and different rotation laws) and find that several universal
relations between properly scaled gravitational mass, rest mass and angular
momentum of the turning-point models that are valid for uniform rotation, are
insensitive to the degree of differential rotation, to high accuracy.
[17]
oai:arXiv.org:1702.02567 [pdf] - 1581434
Semi-analytic derivation of the threshold mass for prompt collapse in
binary neutron star mergers
Submitted: 2017-02-08, last modified: 2017-08-03
The threshold mass for prompt collapse in binary neutron star mergers was
empirically found to depend on the stellar properties of the maximum-mass
non-rotating neutron star model. Here we present a semi-analytic derivation of
this empirical relation which suggests that it is rather insensitive to thermal
effects, to deviations from axisymmetry and to the exact rotation law in merger
remnants. We utilize axisymmetric, cold equilibrium models with differential
rotation and determine the threshold mass for collapse from the comparison
between an empirical relation that describes the angular momentum in the
remnant for a given total binary mass and the sequence of rotating equilibrium
models at the threshold to collapse (the latter assumed to be near the turning
point of fixed-angular-momentum sequences). In spite of the various simplifying
assumptions, the empirical relation for prompt collapse is reproduced with good
accuracy, which demonstrates its robustness. We discuss implications of our
methodology and results for understanding other empirical relations satisfied
by neutron-star merger remnants that have been discovered by numerical
simulations and that play a key role in constraining the high-density equation
of state through gravitational-wave observations.
[19]
oai:arXiv.org:1605.02629 [pdf] - 1442772
On the dynamics of tilted black hole-torus systems
Submitted: 2016-05-09, last modified: 2016-06-17
We present results from three-dimensional, numerical relativity simulations
of a {\it tilted} black hole-thick accretion disc system. The simulations are
analysed using tracer particles in the disc which are advected with the flow.
Such tracers, which we employ in these new simulations for the first time,
provide a powerful means to analyse in detail the complex dynamics of tilted
black hole-torus systems. We show how its use helps to gain insight in the
overall dynamics of the system, discussing the origin of the observed black
hole precession and the development of a global non-axisymmetric $m=1$ mode in
the disc. Our three-dimensional simulations show the presence of quasi-periodic
oscillations (QPOs) in the instantaneous accretion rate, with frequencies in a
range compatible with those observed in low mass X-ray binaries with either a
black hole or a neutron star component. The frequency ratio of the dominant low
frequency peak and the first overtone is $o_1/f \sim 1.9$, a frequency ratio
not attainable when modelling the QPOs as $p$-mode oscillations in axisymmetric
tori.
[20]
oai:arXiv.org:1605.07638 [pdf] - 1418932
Coherent magneto-elastic oscillations in superfluid magnetars
Submitted: 2016-05-24
We study the effect of superfluidity on torsional oscillations of highly
magnetised neutron stars (magnetars) with a microphysical equation of state by
means of two-dimensional, magnetohydrodynamical- elastic simulations. The
superfluid properties of the neutrons in the neutron star core are treated in a
parametric way in which we effectively decouple part of the core matter from
the oscillations. Our simulations confirm the existence of two groups of
oscillations, namely continuum oscillations that are confined to the neutron
star core and are of Alfv\'enic character, and global oscillations with
constant phase and that are of mixed magneto-elastic type. The latter might
explain the quasi-periodic oscillations observed in magnetar giant flares,
since they do not suffer from the additional damping mechanism due to phase
mixing, contrary to what happens for continuum oscillations. However, we cannot
prove rigorously that the coherent oscillations with constant phase are normal
modes. Moreover, we find no crustal shear modes for the magnetic field
strengths typical for magnetars.We provide fits to our numerical simulations
that give the oscillation frequencies as functions of magnetic field strength
and proton fraction in the core.
[21]
oai:arXiv.org:1506.04056 [pdf] - 1439692
Numerical relativity simulations of thick accretion disks around tilted
Kerr black holes
Submitted: 2015-06-12, last modified: 2016-02-27
In this work we present 3D numerical relativity simulations of thick
accretion disks around tilted Kerr BH. We investigate the evolution of three
different initial disk models with a range of initial black hole spin
magnitudes and tilt angles. For all the disk-to-black hole mass ratios
considered (0.044-0.16) we observe significant black hole precession and
nutation during the evolution. This indicates that for such mass ratios,
neglecting the self-gravity of the disks by evolving them in a fixed background
black hole spacetime is not justified. We find that the two more massive models
are unstable against the Papaloizou-Pringle (PP) instability and that those
PP-unstable models remain unstable for all initial spins and tilt angles
considered, showing that the development of the instability is a very robust
feature of such PP-unstable disks. Our lightest model, which is the most
astrophysically favorable outcome of mergers of binary compact objects, is
stable. The tilt between the black hole spin and the disk is strongly modulated
during the growth of the PP instability, causing a partial global realignment
of black hole spin and disk angular momentum in the most massive model with
constant specific angular momentum l. For the model with non-constant l-profile
we observe a long-lived m=1 non-axisymmetric structure which shows strong
oscillations of the tilt angle in the inner regions of the disk. This effect
might be connected to the development of Kozai-Lidov oscillations. Our
simulations also confirm earlier findings that the development of the PP
instability causes the long-term emission of large amplitude gravitational
waves, predominantly for the l=m=2 multipole mode. The imprint of the BH
precession on the gravitational waves from tilted BH-torus systems remains an
interesting open issue that would require significantly longer simulations than
those presented in this work.
[22]
oai:arXiv.org:1602.00950 [pdf] - 1351772
Dynamics and gravitational-wave emission of neutron-star merger remnants
Submitted: 2016-02-01
The coalescence of a neutron-star binary is likely to result in the formation
of a neutron-star merger remnant for a large range of binary mass
configurations. The massive merger remnant shows strong oscillations, which are
excited by the merging process, and emits gravitational waves. Here we discuss
possibilities and prospects of inferring unknown stellar properties of neutron
stars by the detection of postmerger gravitational-wave emission, which thus
leads to constraints of the equation of state of high-density matter. In
particular, the dominant oscillation frequency of the postmerger remnant
provides tight limits to neutron-star radii. We mention first steps towards a
practical implementation of future gravitational-wave searches for the
postmerger emission. Moreover, we outline possibilities to estimate the unknown
maximum mass of nonrotating neutron stars from such types of measurements.
Finally, we review the origin and scientific implications of secondary peaks in
the gravitational-wave spectrum of neutron-star mergers and differences in the
dynamical behavior of the postmerger remnant depending on the binary
configuration. These considerations lead to a unified picture of the
post-merger gravitational-wave emission and the post-merger dynamics.
[23]
oai:arXiv.org:1508.05493 [pdf] - 1392691
Exploring properties of high-density matter through remnants of
neutron-star mergers
Submitted: 2015-08-22, last modified: 2015-12-15
Remnants of neutron-star mergers are essentially massive, hot, differentially
rotating neutron stars, which are initially strongly oscillating. They
represent a unique probe for high-density matter because the oscillations are
detectable via gravitational-wave measurements and are strongly dependent on
the equation of state. The impact of the equation of state is apparent in the
frequency of the dominant oscillation mode of the remnant. For a fixed total
binary mass a tight relation between the dominant postmerger frequency and the
radii of nonrotating neutron stars exists. Inferring observationally the
dominant postmerger frequency thus determines neutron star radii with high
accuracy of the order of a few hundred meters. By considering symmetric and
asymmetric binaries of the same chirp mass, we show that the knowledge of the
binary mass ratio is not critical for this kind of radius measurements. We
summarize different possibilities to deduce the maximum mass of nonrotating
neutron stars. We clarify the nature of the three most prominent features of
the postmerger gravitational-wave spectrum and argue that the merger remnant
can be considered to be a single, isolated, self-gravitating object that can be
described by concepts of asteroseismology. The understanding of the different
mechanisms shaping the gravitational-wave signal yields a physically motivated
analytic model of the gravitational-wave emission, which may form the basis for
template-based gravitational-wave data analysis. We explore the observational
consequences of a scenario of two families of compact stars including hadronic
and quark stars. We find that this scenario leaves a distinctive imprint on the
postmerger gravitational-wave signal. In particular, a strong discontinuity in
the dominant postmerger frequency as function of the total mass will be a
strong indication for two families of compact stars. (abridged)
[24]
oai:arXiv.org:1509.08522 [pdf] - 1388837
Observing Gravitational Waves From The Post-Merger Phase Of Binary
Neutron Star Coalescence
Submitted: 2015-09-28
We present an effective, low-dimensionality frequency-domain template for the
gravitational wave signal from the stellar remnants from binary neutron star
coalescence. A principal component decomposition of a suite of numerical
simulations of binary neutron star mergers is used to construct orthogonal
basis functions for the amplitude and phase spectra of the waveforms for a
variety of neutron star equations of state and binary mass configurations. We
review the phenomenology of late merger / post-merger gravitational wave
emission in binary neutron star coalescence and demonstrate how an
understanding of the dynamics during and after the merger leads to the
construction of a universal spectrum. We also provide a discussion of the
prospects for detecting the post-merger signal in future gravitational wave
detectors as a potential contribution to the science case for third generation
instruments. The template derived in our analysis achieves $>90\%$ match across
a wide variety of merger waveforms and strain sensitivity spectra for current
and potential gravitational wave detectors. A Fisher matrix analysis yields a
preliminary estimate of the typical uncertainty in the determination of the
dominant post-merger oscillation frequency $f_{\mathrm{peak}}$ as $\delta
f_{\mathrm{peak}} \sim 50$Hz. Using recently derived correlations between
$f_{\mathrm{peak}}$ and the neutron star radii, this suggests potential
constraints on the radius of a fiducial neutron star of $\sim 220$\,m. Such
measurements would only be possible for nearby ($\sim 30$Mpc) sources with
advanced LIGO but become more feasible for planned upgrades to advanced LIGO
and other future instruments, leading to constraints on the high density
neutron star equation of state which are independent and complementary to those
inferred from the pre-merger inspiral gravitational wave signal.
[26]
oai:arXiv.org:1503.00972 [pdf] - 1262932
Magnetized neutron stars with superconducting cores: Effect of
entrainment
Submitted: 2015-03-03, last modified: 2015-08-16
We construct equilibrium configurations of magnetized, two-fluid neutron
stars using an iterative numerical method. Working in Newtonian framework we
assume that the neutron star has two regions: the core, which is modelled as a
two-component fluid consisting of type-II superconducting protons and
superfluid neutrons, and the crust, a region composed of normal matter. Taking
a new step towards more complete equilibrium models, we include the effect of
entrainment, which implies that a magnetic force acts on neutrons, too. We
consider purely poloidal field cases and present improvements to an earlier
numerical scheme for solving equilibrium equations, by introducing new
convergence criteria. We find that entrainment results in qualitative
differences in the structure of field lines along the magnetic axis.
[27]
oai:arXiv.org:1502.03176 [pdf] - 1095587
A unified picture of the post-merger dynamics and gravitational wave
emission in neutron-star mergers
Submitted: 2015-02-10, last modified: 2015-06-02
We introduce a classification scheme of the post-merger dynamics and
gravitational-wave emission in binary neutron star mergers, after identifying a
new mechanism by which a secondary peak in the gravitational-wave spectrum is
produced. It is caused by a spiral deformation, the pattern of which rotates
slower with respect to the double-core structure in center of the remnant. This
secondary peak is typically well separated in frequency from the secondary peak
produced by a nonlinear interaction between a quadrupole and a quasi-radial
oscillation. The new mechanism allows for an explanation of low-frequency
modulations seen in a number of physical characteristics of the remnant, such
as the central lapse function, the maximum density and the separation between
the two cores. We find empirical relations for both types of secondary peaks
between their gravitational-wave frequency and the compactness of nonrotating
individual neutron stars, that exist for fixed total binary masses. These
findings are derived for equal-mass binaries without intrinsic neutron-star
spin analyzing hydrodynamical simulations without magnetic field effects. Our
classification scheme may form the basis for the construction of detailed
gravitational-wave templates of the post-merger phase. We find that the
quasi-radial oscillation frequency of the remnant decreases with the total
binary mass. For a given merger event our classification scheme may allow to
determine the proximity of the measured total binary mass to the threshold mass
for prompt black hole formation, which can, in turn, yield an estimate of the
maximum neutron-star mass.
[28]
oai:arXiv.org:1503.08769 [pdf] - 972708
Inferring neutron-star properties from gravitational-wave signals of
binary mergers
Submitted: 2015-03-30
The oscillations of a merger remnant forming after the coalescence of two
neutron stars are very characteristic for the high-density equation of state.
The dominant oscillation frequency occurs as a pronounced peak in the kHz range
of the gravitational-wave spectrum. We describe how the dominant oscillation
frequency of the remnant can be employed to infer the radii of non-rotating
neutron stars.
[29]
oai:arXiv.org:1503.06636 [pdf] - 969464
Unusual behavior of epicyclic frequencies around rapidly rotating
compact stars
Submitted: 2015-03-23
We report on numerical calculations of orbital and epicyclic frequencies in
nearly circular orbits around rotating neutron stars and strange quark stars.
The FPS equation of state was used to describe the neutron star structure while
the MIT bag model was used to model the equation of state of strange quark
stars. The uniformly rotating stellar configurations were computed in full
general relativity. We find that the vertical epicyclic frequency is very
sensitive to the oblateness of the rotating star. For models of rotating
neutron stars of moderate mass, as well as for strange quark star models, the
sense of the nodal precession of test particle orbits close to the star changes
at a certain stellar rotation rate. These findings may have implications for
models of kHz QPOs.
[30]
oai:arXiv.org:1406.7375 [pdf] - 870103
On the accuracy of the IWM-CFC approximation in differentially rotating
relativistic stars
Submitted: 2014-06-28, last modified: 2014-09-10
We determine the accuracy of the conformal flatness (IWM-CFC) approximation
for the case of single, but strongly differentially rotating relativistic
stars. We find that for the fastest rotating and most relativistic polytropic
models, the deviation from full general relativity is below 5% for integrated
quantities and below 10% for local quantities, such as the angular velocity.
Furthermore, we study the deviation of the IWM-CFC approximation from full
general relativity by evaluating and comparing different error indicators. We
find that for the models that are not near the maximum mass, a simple error
indicator constructed from local values of the metric potentials is more
indicative of the accuracy of the IWM-CFC approximation than an error indicator
that is based on the Cotton-York tensor. Furthermore, we construct a simple,
linear empirical relation that allows for the estimation of the error made by
the IWM-CFC approximation and which only involves the flattening of the star
due to rotation and the minimum value of the lapse function. Thus, in any
numerical simulation involving rotating relativistic stars, one can readily
know the deviations from full general relativity due to the IWM-CFC
approximation.
[31]
oai:arXiv.org:1408.6526 [pdf] - 864022
The Large Observatory For x-ray Timing
Feroci, M.;
Herder, J. W. den;
Bozzo, E.;
Barret, D.;
Brandt, S.;
Hernanz, M.;
van der Klis, M.;
Pohl, M.;
Santangelo, A.;
Stella, L.;
Watts, A.;
Wilms, J.;
Zane, S.;
Ahangarianabhari, M.;
Albertus, C.;
Alford, M.;
Alpar, A.;
Altamirano, D.;
Alvarez, L.;
Amati, L.;
Amoros, C.;
Andersson, N.;
Antonelli, A.;
Argan, A.;
Artigue, R.;
Artigues, B.;
Atteia, J. -L.;
Azzarello, P.;
Bakala, P.;
Baldazzi, G.;
Balman, S.;
Barbera, M.;
van Baren, C.;
Bhattacharyya, S.;
Baykal, A.;
Belloni, T.;
Bernardini, F.;
Bertuccio, G.;
Bianchi, S.;
Bianchini, A.;
Binko, P.;
Blay, P.;
Bocchino, F.;
Bodin, P.;
Bombaci, I.;
Bidaud, J. -M. Bonnet;
Boutloukos, S.;
Bradley, L.;
Braga, J.;
Brown, E.;
Bucciantini, N.;
Burderi, L.;
Burgay, M.;
Bursa, M.;
Budtz-Jørgensen, C.;
Cackett, E.;
Cadoux, F. R.;
Cais, P.;
Caliandro, G. A.;
Campana, R.;
Campana, S.;
Capitanio, F.;
Casares, J.;
Casella, P.;
Castro-Tirado, A. J.;
Cavazzuti, E.;
Cerda-Duran, P.;
Chakrabarty, D.;
Château, F.;
Chenevez, J.;
Coker, J.;
Cole, R.;
Collura, A.;
Cornelisse, R.;
Courvoisier, T.;
Cros, A.;
Cumming, A.;
Cusumano, G.;
D'Aì, A.;
D'Elia, V.;
Del Monte, E.;
De Luca, A.;
De Martino, D.;
Dercksen, J. P. C.;
De Pasquale, M.;
De Rosa, A.;
Del Santo, M.;
Di Cosimo, S.;
Diebold, S.;
Di Salvo, T.;
Donnarumma, I.;
Drago, A.;
Durant, M.;
Emmanoulopoulos, D.;
Erkut, M. H.;
Esposito, P.;
Evangelista, Y.;
Fabian, A.;
Falanga, M.;
Favre, Y.;
Feldman, C.;
Ferrari, V.;
Ferrigno, C.;
Finger, M.;
Finger, M. H.;
Fraser, G. W.;
Frericks, M.;
Fuschino, F.;
Gabler, M.;
Galloway, D. K.;
Sanchez, J. L. Galvez;
Garcia-Berro, E.;
Gendre, B.;
Gezari, S.;
Giles, A. B.;
Gilfanov, M.;
Giommi, P.;
Giovannini, G.;
Giroletti, M.;
Gogus, E.;
Goldwurm, A.;
Goluchová, K.;
Götz, D.;
Gouiffes, C.;
Grassi, M.;
Groot, P.;
Gschwender, M.;
Gualtieri, L.;
Guidorzi, C.;
Guy, L.;
Haas, D.;
Haensel, P.;
Hailey, M.;
Hansen, F.;
Hartmann, D. H.;
Haswell, C. A.;
Hebeler, K.;
Heger, A.;
Hermsen, W.;
Homan, J.;
Hornstrup, A.;
Hudec, R.;
Huovelin, J.;
Ingram, A.;
Zand, J. J. M. in't;
Israel, G.;
Iwasawa, K.;
Izzo, L.;
Jacobs, H. M.;
Jetter, F.;
Johannsen, T.;
Jacobs, H. M.;
Jonker, P.;
Josè, J.;
Kaaret, P.;
Kanbach, G.;
Karas, V.;
Karelin, D.;
Kataria, D.;
Keek, L.;
Kennedy, T.;
Klochkov, D.;
Kluzniak, W.;
Kokkotas, K.;
Korpela, S.;
Kouveliotou, C.;
Kreykenbohm, I.;
Kuiper, L. M.;
Kuvvetli, I.;
Labanti, C.;
Lai, D.;
Lamb, F. K.;
Laubert, P. P.;
Lebrun, F.;
Lin, D.;
Linder, D.;
Lodato, G.;
Longo, F.;
Lund, N.;
Maccarone, T. J.;
Macera, D.;
Maestre, S.;
Mahmoodifar, S.;
Maier, D.;
Malcovati, P.;
Mandel, I.;
Mangano, V.;
Manousakis, A.;
Marisaldi, M.;
Markowitz, A.;
Martindale, A.;
Matt, G.;
McHardy, I. M.;
Melatos, A.;
Mendez, M.;
Mereghetti, S.;
Michalska, M.;
Migliari, S.;
Mignani, R.;
Miller, M. C.;
Miller, J. M.;
Mineo, T.;
Miniutti, G.;
Morsink, S.;
Motch, C.;
Motta, S.;
Mouchet, M.;
Mouret, G.;
Mulačová, J.;
Muleri, F.;
Muñoz-Darias, T.;
Negueruela, I.;
Neilsen, J.;
Norton, A. J.;
Nowak, M.;
O'Brien, P.;
Olsen, P. E. H.;
Orienti, M.;
Orio, M.;
Orlandini, M.;
Orleanski, P.;
Osborne, J. P.;
Osten, R.;
Ozel, F.;
Pacciani, L.;
Paolillo, M.;
Papitto, A.;
Paredes, J. M.;
Patruno, A.;
Paul, B.;
Perinati, E.;
Pellizzoni, A.;
Penacchioni, A. V.;
Perez, M. A.;
Petracek, V.;
Pittori, C.;
Pons, J.;
Portell, J.;
Possenti, A.;
Poutanen, J.;
Prakash, M.;
Provost, P. Le;
Psaltis, D.;
Rambaud, D.;
Ramon, P.;
Ramsay, G.;
Rapisarda, M.;
Rachevski, A.;
Rashevskaya, I.;
Ray, P. S.;
Rea, N.;
Reddy, S.;
Reig, P.;
Aranda, M. Reina;
Remillard, R.;
Reynolds, C.;
Rezzolla, L.;
Ribo, M.;
de la Rie, R.;
Riggio, A.;
Rios, A.;
Gil, P. Rodríguez-;
Rodriguez, J.;
Rohlfs, R.;
Romano, P.;
Rossi, E. M. R.;
Rozanska, A.;
Rousseau, A.;
Ryde, F.;
Sabau-Graziati, L.;
Sala, G.;
Salvaterra, R.;
Sanna, A.;
Sandberg, J.;
Scaringi, S.;
Schanne, S.;
Schee, J.;
Schmid, C.;
Shore, S.;
Schneider, R.;
Schwenk, A.;
Schwope, A. D.;
Seyler, J. -Y.;
Shearer, A.;
Smith, A.;
Smith, D. M.;
Smith, P. J.;
Sochora, V.;
Soffitta, P.;
Soleri, P.;
Spencer, A.;
Stappers, B.;
Steiner, A. W.;
Stergioulas, N.;
Stratta, G.;
Strohmayer, T. E.;
Stuchlik, Z.;
Suchy, S.;
Sulemainov, V.;
Takahashi, T.;
Tamburini, F.;
Tauris, T.;
Tenzer, C.;
Tolos, L.;
Tombesi, F.;
Tomsick, J.;
Torok, G.;
Torrejon, J. M.;
Torres, D. F.;
Tramacere, A.;
Trois, A.;
Turolla, R.;
Turriziani, S.;
Uter, P.;
Uttley, P.;
Vacchi, A.;
Varniere, P.;
Vaughan, S.;
Vercellone, S.;
Vrba, V.;
Walton, D.;
Watanabe, S.;
Wawrzaszek, R.;
Webb, N.;
Weinberg, N.;
Wende, H.;
Wheatley, P.;
Wijers, R.;
Wijnands, R.;
Wille, M.;
Wilson-Hodge, C. A.;
Winter, B.;
Wood, K.;
Zampa, G.;
Zampa, N.;
Zampieri, L.;
Zdunik, L.;
Zdziarski, A.;
Zhang, B.;
Zwart, F.;
Ayre, M.;
Boenke, T.;
van Damme, C. Corral;
Kuulkers, E.;
Lumb, D.
Submitted: 2014-08-27, last modified: 2014-08-29
The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3
Cosmic Vision framework and participated in the final down-selection for a
launch slot in 2022-2024. Thanks to the unprecedented combination of effective
area and spectral resolution of its main instrument, LOFT will study the
behaviour of matter under extreme conditions, such as the strong gravitational
field in the innermost regions of accretion flows close to black holes and
neutron stars, and the supra-nuclear densities in the interior of neutron
stars. The science payload is based on a Large Area Detector (LAD, 10 m 2
effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of
view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1
arcmin source location accuracy, 300 eV spectral resolution). The WFM is
equipped with an on-board system for bright events (e.g. GRB) localization. The
trigger time and position of these events are broadcast to the ground within 30
s from discovery. In this paper we present the status of the mission at the end
of its Phase A study.
[32]
oai:arXiv.org:1407.7672 [pdf] - 1215977
Modulating the magnetosphere of magnetars by internal magneto-elastic
oscillations
Submitted: 2014-07-29
We couple internal torsional, magneto-elastic oscillations of highly
magnetized neutron stars (magnetars) to their magnetospheres. The corresponding
axisymmetric perturbations of the external magnetic field configuration evolve
as a sequence of linear, force-free equilibria that are completely determined
by the background magnetic field configuration and by the perturbations of the
magnetic field at the surface. The perturbations are obtained from simulations
of magneto-elastic oscillations in the interior of the magnetar. While such
oscillations can excite travelling Alfv\'en waves in the exterior of the star
only in a very limited region close to the poles, they still modulate the near
magnetosphere by inducing a time-dependent twist between the foot-points of
closed magnetic field lines that exit the star at a polar angle $\gtrsim
0.19\,$rad. Moreover, we find that for a dipole-like background magnetic field
configuration the magnetic field modulations in the magnetosphere, driven by
internal oscillations, can only be symmetric with respect to the equator. This
is in agreement with our previous findings, where we interpreted the observed
quasi-periodic oscillations in the X-ray tail of magnetar bursts as driven by
the family of internal magneto-elastic oscillations with symmetric magnetic
field perturbations.
[33]
oai:arXiv.org:1406.5444 [pdf] - 1215111
Prospects For High Frequency Burst Searches Following Binary Neutron
Star Coalescence With Advanced Gravitational Wave Detectors
Submitted: 2014-06-20
The equation of state plays a critical role in the physics of the merger of
two neutron stars. Recent numerical simulations with microphysical equation of
state suggest the outcome of such events depends on the mass of the neutron
stars. For less massive systems, simulations favor the formation of a
hypermassive, quasi-stable neutron star, whose oscillations produce a short,
high frequency burst of gravitational radiation. Its dominant frequency content
is tightly correlated with the radius of the neutron star, and its measurement
can be used to constrain the supranuclear equation of state. In contrast, the
merger of higher mass systems results in prompt gravitational collapse to a
black hole. We have developed an algorithm which combines waveform
reconstruction from a morphology-independent search for gravitational wave
transients with Bayesian model selection, to discriminate between post-merger
scenarios and accurately measure the dominant oscillation frequency. We
demonstrate the efficacy of the method using a catalogue of simulated binary
merger signals in data from LIGO and Virgo, and we discuss the prospects for
this analysis in advanced ground-based gravitational wave detectors. From the
waveforms considered in this work and assuming an optimally oriented source, we
find that the post-merger neutron star signal may be detectable by this
technique to $\sim 10\text{--}25$\,Mpc. We also find that we successfully
discriminate between the post-merger scenarios with $\sim 95\%$ accuracy and
determine the dominant oscillation frequency of surviving post-merger neutron
stars to within $\sim 10$\,Hz, averaged over all detected signals. This leads
to an uncertainty in the estimated radius of a non-rotating 1.6\,M$_{\odot}$
reference neutron star of $\sim 100\,$m.
[34]
oai:arXiv.org:1403.5301 [pdf] - 1208519
Revealing the high-density equation of state through binary neutron star
mergers
Submitted: 2014-03-20, last modified: 2014-06-19
We present a novel method for revealing the equation of state of high-density
neutron star matter through gravitational waves emitted during the postmerger
phase of a binary neutron star system. The method relies on a small number of
detections of the peak frequency in the postmerger phase for binaries of
different (relatively low) masses, in the most likely range of expected
detections. From such observations, one can construct the derivative of the
peak frequency versus the binary mass, in this mass range. Through a detailed
study of binary neutron star mergers for a large sample of equations of state,
we show that one can extrapolate the above information to the highest possible
mass (the threshold mass for black hole formation in a binary neutron star
merger). In turn, this allows for an empirical determination of the maximum
mass of cold, nonrotating neutron stars to within 0.1 M_sun, while the
corresponding radius is determined to within a few percent. Combining this with
the determination of the radius of cold, nonrotating neutron stars of 1.6 M_sun
(to within a few percent, as was demonstrated in Bauswein et al., PRD, 86,
063001, 2012), allows for a clear distinction of a particular candidate
equation of state among a large set of other candidates. Our method is
particularly appealing because it reveals simultaneously the moderate and very
high-density parts of the equation of state, enabling the distinction of
mass-radius relations even if they are similar at typical neutron star masses.
Furthermore, our method also allows to deduce the maximum central energy
density and maximum central rest-mass density of cold, nonrotating neutron
stars with an accuracy of a few per cent.
[35]
oai:arXiv.org:1405.6976 [pdf] - 860766
Orbital and epicyclic frequencies around rapidly rotating compact stars
in scalar-tensor theories of gravity
Submitted: 2014-05-27
We study the orbital and epicyclic frequencies of particles orbiting around
rapidly rotating neutron stars and strange stars in a particular scalar-tensor
theory of gravity. We find very large deviations of these frequencies, when
compared to their corresponding values in general relativity, for the
maximum-mass rotating models. In contrast, for models rotating with spin
frequency of 700Hz (approximately the largest known rotation rate of neutron
stars), the deviations are generally small. Nevertheless, for a very stiff
equation of state and a high mass the deviation of one of the epicyclic
frequencies from its GR value is appreciable even at a spin frequency of 700Hz.
In principle, such a deviation could become important in models of
quasi-periodic oscillations in low-mass x-ray binaries and could serve as a
test of strong gravity (if other parameters are well constraint). Even though
the present paper is concentrated mainly on orbital and epicyclic frequencies,
we present here for the first time rapidly rotating, scalarized equilibrium
compact stars with realistic hadronic equations of state and strange matter
equation of state. We also provide analytical expressions for the exterior
spacetime of scalarized neutron stars and their epicyclic frequencies in the
nonrotating limit.
[36]
oai:arXiv.org:1403.1129 [pdf] - 1208193
Epicyclic frequencies for rotating strange quark stars: the importance
of stellar oblateness
Submitted: 2014-03-05
Kilohertz QPOs can be used as a probe of the inner regions of accretion disks
in compact stars and hence also of the properties of the central object. Most
models of kHz QPOs involve epicyclic frequencies to explain their origin. We
compute the epicyclic frequencies of nearly circular orbits around rotating
strange quark stars. The MIT bag model is used to model the equation of state
of quark matter and the uniformly rotating stellar configurations are computed
in full general relativity. The vertical epicyclic frequency and the related
nodal precession rate of inclined orbits are very sensitive to the oblateness
of the rotating star. For slowly rotating stellar models of moderate and high
mass strange stars, the sense of the nodal precession changes at a certain
rotation rate. At lower stellar rotation rates the orbital nodal precession is
prograde, as it is in the Kerr metric, while at higher rotation rates the
precession is retrograde, as it is for Maclaurin spheroids. Thus,
qualitatively, the orbits around rapidly rotating strange quark stars are
affected more strongly by the effects of stellar oblateness than by the effects
of general relativity. We show that epicyclic and orbital frequencies
calculated numerically for small mass strange stars are in very good agreement
with analytical formulae for Maclaurin spheroids.
[37]
oai:arXiv.org:1309.0605 [pdf] - 787674
Rapidly rotating neutron stars in scalar-tensor theories of gravity
Submitted: 2013-09-03, last modified: 2014-02-24
We present the field equations governing the equilibrium of rapidly rotating
neutron stars in scalar-tensor theories of gravity, as well as representative
numerical solutions. The conditions for the presence of a nontrivial scalar
field and the deviations from the general relativistic solutions are studied.
Two examples of scalar-tensor theories are examined - one case that is
equivalent to the Brans-Dicke theory and a second case, that is perturbatively
equivalent to Einstein's General Relativity in the weak field regime, but can
differ significantly for strong fields. Our numerical results show that rapidly
rotating neutron star models with a nontrivial scalar field exist in both cases
and that the effect of the scalar field is stronger for rapid rotation. If we
consider values of the coupling parameters in accordance with current
observations, only the second example of scalar-tensor theories has significant
influence on the neutron star structure. We show that scalarized, rapidly
rotating neutron stars exist for a larger range of the parameters than in the
static case, since a nontrivial scalar field is present even for values of the
coupling constant $\beta>-4.35$, and that these solutions are energetically
more favorable than the general relativistic ones. In addition, the deviations
of the rapidly rotating scalar-tensor neutron stars from the
general-relativistic solutions is significantly larger than in the static case.
[38]
oai:arXiv.org:1310.7436 [pdf] - 787690
Breakdown of I-Love-Q universality in rapidly rotating relativistic
stars
Submitted: 2013-10-28, last modified: 2014-02-23
It was shown recently, that normalized relations between the moment of
inertia (I), the quadrupole moment (Q) and the tidal deformability (Love
number) exist and for slowly-rotating neutron stars they are almost independent
of the equation of state (EOS). We extend the computation of the I-Q relation
to models rotating up to the mass-shedding limit and show that the universality
of the relations is lost. With increasing rotation rate, the normalized I-Q
relation departs significantly from its slow-rotation limit deviating up to 40%
for neutron stars and up to 75% for strange stars. The deviation is also EOS
dependent and for a broad set of hadronic and strange matter EOS the spread due
to rotation is comparable to the spread due to the EOS, if one considers
sequences with fixed rotational frequency. Still, for a restricted sample of
modern realistic EOS one can parameterize the deviations from universality as a
function of rotation only. The previously proposed I-Love-Q relations should
thus be used with care, because they lose their universality in astrophysical
situations involving compact objects rotating faster than a few hundred Hz.
[40]
oai:arXiv.org:1210.1214 [pdf] - 657312
The runaway instability in general relativistic accretion disks
Submitted: 2012-10-03, last modified: 2013-04-27
When an accretion disk falls prey to the runaway instability, a large portion
of its mass is devoured by the black hole within a few dynamical times. Despite
decades of effort, it is still unclear under what conditions such an
instability can occur. The technically most advanced relativistic simulations
to date were unable to find a clear sign for the onset of the instability. In
this work, we present three-dimensional relativistic hydrodynamics simulations
of accretion disks around black holes in dynamical space-time. We focus on the
configurations that are expected to be particularly prone to the development of
this instability. We demonstrate, for the first time, that the fully
self-consistent general relativistic evolution does indeed produce a runaway
instability.
[41]
oai:arXiv.org:1304.3566 [pdf] - 1165929
Imprints of superfluidity on magneto-elastic QPOs of SGRs
Submitted: 2013-04-12
Our numerical simulations show that axisymmetric, torsional, magneto-elastic
oscillations of magnetars with a superfluid core can explain the whole range of
observed quasi-periodic oscillations (QPOs) in the giant flares of soft
gamma-ray repeaters. There exist constant phase, magneto-elastic QPOs at both
low (f<150 Hz) and high frequencies (f>500 Hz), in full agreement with
observations. The range of magnetic field strengths required to match the
observed QPO frequencies agrees with that from spin-down estimates. These
results strongly suggest that neutrons in magnetar cores are superfluid.
[42]
oai:arXiv.org:1208.6443 [pdf] - 1151015
Magneto-elastic oscillations of neutron stars: exploring different
magnetic field configurations
Submitted: 2012-08-31, last modified: 2013-02-27
We study magneto-elastic oscillations of highly magnetized neutron stars
(magnetars) which have been proposed as an explanation for the quasi-periodic
oscillations (QPOs) appearing in the decaying tail of the giant flares of soft
gamma-ray repeaters (SGRs). We extend previous studies by investigating various
magnetic field configurations, computing the Alfv\'en spectrum in each case and
performing magneto-elastic simulations for a selected number of models. By
identifying the observed frequencies of 28 Hz (SGR 1900+14) and 30 Hz (SGR
1806-20) with the fundamental Alfv\'en QPOs, we estimate the required surface
magnetic field strength. For the magnetic field configurations investigated
(dipole-like poloidal, mixed toroidal-poloidal with a dipole-like poloidal
component and a toroidal field confined to the region of field lines closing
inside the star, and for poloidal fields with an additional quadrupole-like
component) the estimated dipole spin-down magnetic fields are between 8x10^14 G
and 4x10^15 G, in broad agreement with spin-down estimates for the SGR sources
producing giant flares. A number of these models exhibit a rich Alfv\'en
continuum revealing new turning points which can produce QPOs. This allows one
to explain most of the observed QPO frequencies as associated with
magneto-elastic QPOs. In particular, we construct a possible configuration with
two turning points in the spectrum which can explain all observed QPOs of SGR
1900+14. Finally, we find that magnetic field configurations which are entirely
confined in the crust (if the core is assumed to be a type I superconductor)
are not favoured, due to difficulties in explaining the lowest observed QPO
frequencies (f<30 Hz).
[43]
oai:arXiv.org:1209.1497 [pdf] - 1151206
LOFT: the Large Observatory For X-ray Timing
Feroci, M.;
Herder, J. W. den;
Bozzo, E.;
Barret, D.;
Brandt, S.;
Hernanz, M.;
van der Klis, M.;
Pohl, M.;
Santangelo, A.;
Stella, L.;
Watts, A.;
Wilms, J.;
Zane, S.;
Ahangarianabhari, M.;
Alpar, A.;
Altamirano, D.;
Alvarez, L.;
Amati, L.;
Amoros, C.;
Andersson, N.;
Antonelli, A.;
Argan, A.;
Artigue, R.;
Azzarello, P.;
Baldazzi, G.;
Balman, S.;
Barbera, M.;
Belloni, T.;
Bertuccio, G.;
Bianchi, S.;
Bianchini, A.;
Bodin, P.;
Bidaud, J. -M. Bonnet;
Boutloukos, S.;
Braga, J.;
Brown, E.;
Bucciantini, N.;
Burderi, L.;
Bursa, M.;
Budtz-Jørgensen, C.;
Cackett, E.;
Cadoux, F. R.;
Cais, P.;
Caliandro, G. A.;
Campana, R.;
Campana, S.;
Casella, P.;
Chakrabarty, D.;
Chenevez, J.;
Coker, J.;
Cole, R.;
Collura, A.;
Courvoisier, T.;
Cros, A.;
Cumming, A.;
Cusumano, G.;
D'Aì, A.;
D'Elia, V.;
Del Monte, E.;
De Martino, D.;
De Rosa, A.;
Di Cosimo, S.;
Diebold, S.;
Di Salvo, T.;
Donnarumma, I.;
Drago, A.;
Durant, M.;
Emmanoulopoulos, D.;
Evangelista, Y.;
Fabian, A.;
Falanga, M.;
Favre, Y.;
Feldman, C.;
Ferrigno, C.;
Finger, M. H.;
Fraser, G. W.;
Fuschino, F.;
Galloway, D. K.;
Sanchez, J. L. Galvez;
Garcia-Berro, E.;
Gendre, B.;
Gezari, S.;
Giles, A. B.;
Gilfanov, M.;
Giommi, P.;
Giovannini, G.;
Giroletti, M.;
Goldwurm, A.;
Götz, D.;
Gouiffes, C.;
Grassi, M.;
Guidorzi, P. Groot C.;
Haas, D.;
Hansen, F.;
Hartmann, D. H.;
Haswe, C. A.;
Heger, A.;
Homan, J.;
Hornstrup, A.;
Hudec, R.;
Huovelin, J.;
Ingram, A.;
Zand, J. J. M. in't;
Isern, J.;
Israe, G.;
Izzo, L.;
Jonker, P.;
Kaaret, P.;
Karas, V.;
Karelin, D.;
Kataria, D.;
Keek, L.;
Kennedy, T.;
Klochkov, D.;
Kluzniak, W.;
Kokkotas, K.;
Korpela, S.;
Kouveliotou, C.;
Kreykenbohm, I.;
Kuiper, L. M.;
Kuvvetli, I.;
Labanti, C.;
Lai, D.;
Lamb, F. K.;
Lebrun, F.;
Lin, D.;
Linder, D.;
Lodato, G.;
Longo, F.;
Lund, N.;
Maccarone, T. J.;
Macera, D.;
Maier, D.;
Malcovati, P.;
Mangano, V.;
Manousakis, A.;
Marisaldi, M.;
Markowitz, A.;
Martindale, A.;
Matt, G.;
McHardy, I. M.;
Melatos, A.;
Mendez, M.;
Migliari, S.;
Mignani, R.;
Miller, M. C.;
Miller, J. M.;
Mineo, T.;
Miniutti, G.;
Morsink, S.;
Motch, C.;
Motta, S.;
Mouchet, M.;
Muleri, F.;
Norton, A. J.;
Nowak, M.;
O'Brien, P.;
Orienti, M.;
Orio, M.;
Orlandini, M.;
Orleanski, P.;
Osborne, J. P.;
Osten, R.;
Ozel, F.;
Pacciani, L.;
Papitto, A.;
Paul, B.;
Perinati, E.;
Petracek, V.;
Portell, J.;
Poutanen, J.;
Psaltis, D.;
Rambaud, D.;
Ramsay, G.;
Rapisarda, M.;
Rachevski, A.;
Ray, P. S.;
Rea, N.;
Reddy, S.;
Reig, P.;
Aranda, M. Reina;
Remillard, R.;
Reynolds, C.;
Rodríguez-Gil, P.;
Rodriguez, J.;
Romano, P.;
Rossi, E. M. R.;
Ryde, F.;
Sabau-Graziati, L.;
Sala, G.;
Salvaterra, R.;
Sanna, A.;
Schanne, S.;
Schee, J.;
Schmid, C.;
Schwenk, A.;
Schwope, A. D.;
Seyler, J. -Y.;
Shearer, A.;
Smith, A.;
Smith, D. M.;
Smith, P. J.;
Sochora, V.;
Soffitta, P.;
Soleri, P.;
Stappers, B.;
Stelzer, B.;
Stergioulas, N.;
Stratta, G.;
Strohmayer, T. E.;
Stuchlik, Z.;
Suchy, S.;
Sulemainov, V.;
Takahashi, T.;
Tamburini, F.;
Tenzer, C.;
Tolos, L.;
Torok, G.;
Torrejon, J. M.;
Torres, D. F.;
Tramacere, A.;
Trois, A.;
Turriziani, S.;
Uter, P.;
Uttley, P.;
Vacchi, A.;
Varniere, P.;
Vaughan, S.;
Vercellone, S.;
Vrba, V.;
Walton, D.;
Watanabe, S.;
Wawrzaszek, R.;
Webb, N.;
Weinberg, N.;
Wende, H.;
Wheatley, P.;
Wijers, R.;
Wijnands, R.;
Wille, M.;
Wilson-Hodge, C. A.;
Winter, B.;
Wood, K.;
Zampa, G.;
Zampa, N.;
Zampieri, L.;
Zdziarski, A.;
Zhang, B.
Submitted: 2012-09-07
The LOFT mission concept is one of four candidates selected by ESA for the M3
launch opportunity as Medium Size missions of the Cosmic Vision programme. The
launch window is currently planned for between 2022 and 2024. LOFT is designed
to exploit the diagnostics of rapid X-ray flux and spectral variability that
directly probe the motion of matter down to distances very close to black holes
and neutron stars, as well as the physical state of ultra-dense matter. These
primary science goals will be addressed by a payload composed of a Large Area
Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1
degree field of view) experiment operating in the energy range 2-50 keV, with a
10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM
will operate in the same energy range as the LAD, enabling simultaneous
monitoring of a few-steradian wide field of view, with an angular resolution of
<5 arcmin. The LAD and WFM experiments will allow us to investigate variability
from submillisecond QPO's to year-long transient outbursts. In this paper we
report the current status of the project.
[44]
oai:arXiv.org:1109.6233 [pdf] - 499282
Magneto-elastic oscillations of neutron stars with dipolar magnetic
fields
Submitted: 2011-09-28, last modified: 2012-04-12
By means of two dimensional, general-relativistic, magneto-hydrodynamical
simulations we investigate the oscillations of magnetized neutron star models
(magnetars) including the description of an extended solid crust. The aim of
this study is to understand the origin of the QPOs observed in the giant flares
of SGRs. We confirm the existence of three different regimes: (a) a weak
magnetic field regime B<5 x 10^13 G, where crustal shear modes dominate the
evolution; (b) a regime of intermediate magnetic fields 5 x 10^13 G<B< 10^15 G,
where Alfv\'en QPOs are mainly confined to the core of the neutron star and the
crustal shear modes are damped very efficiently; and (c) a strong field regime
B>10^15 G, where magneto-elastic oscillations reach the surface and approach
the behavior of purely Alfv\'en QPOs. When the Alfv\'en QPOs are confined to
the core of the neutron star, we find qualitatively similar QPOs as in the
absence of a crust. The lower QPOs associated with the closed field lines of
the dipolar magnetic field configuration are reproduced as in our previous
simulations without crust, while the upper QPOs connected to the open field
lines are displaced from the polar axis. Additionally, we observe a family of
edge QPOs. Our results do not leave much room for a crustal-mode interpretation
of observed QPOs in SGR giant flares, but can accommodate an interpretation of
these observations as originating from Alfv\'en-like, global, turning-point
QPOs in models with dipolar magnetic field strengths in the narrow range of 5 x
10^15 G < B < 1.4 x 10^16 G. This range is somewhat larger than estimates for
magnetic field strengths in known magnetars. The discrepancy may be resolved in
models including a more complicated magnetic field structure or with models
taking superfluidity of the neutrons and superconductivity of the protons in
the core into account.
[45]
oai:arXiv.org:1105.0122 [pdf] - 966780
Collapse of differentially rotating neutron stars and cosmic censorship
Submitted: 2011-04-30, last modified: 2011-07-07
We present new results on the dynamics and gravitational-wave emission from
the collapse of differentially rotating neutron stars. We have considered a
number of polytropic stellar models having different values of the
dimensionless angular momentum J/M^2, where J and M are the asymptotic angular
momentum and mass of the star, respectively. For neutron stars with J/M^2<1,
i.e., "sub-Kerr" models, we were able to find models that are dynamically
unstable and that collapse promptly to a rotating black hole. Both the dynamics
of the collapse and the consequent emission of gravitational waves resemble the
one seen for uniformly rotating stars, although with an overall decrease in the
efficiency of gravitational-wave emission. For stellar models with J/M^2>1,
i.e., "supra-Kerr" models, on the other hand, we were not able to find models
that are dynamically unstable and all of the computed supra-Kerr models were
found to be far from the stability threshold. For these models a gravitational
collapse is possible only after a very severe and artificial reduction of the
pressure, which then leads to a torus developing nonaxisymmetric instabilities
and eventually contracting to a stable axisymmetric stellar configuration.
While this does not exclude the possibility that a naked singularity can be
produced by the collapse of a differentially rotating star, it also suggests
that cosmic censorship is not violated and that generic conditions for a
supra-Kerr progenitor do not lead to a naked singularity.
[46]
oai:arXiv.org:1105.0368 [pdf] - 1076353
Gravitational waves and nonaxisymmetric oscillation modes in mergers of
compact object binaries
Submitted: 2011-05-02
We study the excitation of nonaxisymmetric modes in the post-merger phase of
binary compact object mergers and the associated gravitational wave emission.
Our analysis is based on general-relativistic simulations, in the spatial
conformal flatness approximation, using smoothed-particle-hydrodynamics for the
evolution of matter, and we use a set of equal and unequal mass models,
described by two nonzero-temperature hadronic equations of state and by one
strange star equation of state. Through Fourier transforms of the evolution of
matter variables, we can identify a number of oscillation modes, as well as
several nonlinear components (combination frequencies). We focus on the
dominant m=2 mode, which forms a triplet with two nonlinear components that are
the result of coupling to the quasiradial mode. A corresponding triplet of
frequencies is identified in the gravitational wave spectrum, when the
individual masses of the compact objects are in the most likely range of 1.2 to
1.35 $M_\odot$. We can thus associate, through direct analysis of the dynamics
of the fluid, a specific frequency peak in the gravitational wave spectrum with
the nonlinear component resulting from the difference between the m=2 mode and
the quasiradial mode. Once such observations become available, both the m=2 and
quasiradial mode frequencies could be extracted, allowing for the application
of gravitational-wave asteroseismology to the post-merger remnant and leading
to tight constraints on the equation of state of high-density matter.
[47]
oai:arXiv.org:1011.5406 [pdf] - 335933
Magneto-elastic torsional oscillations of magnetars
Submitted: 2010-11-24
We extend a general-relativistic ideal magneto-hydrodynamical code to include
the effects of elasticity. Using this numerical tool we analyse the
magneto-elastic oscillations of highly magnetised neutron stars (magnetars). In
simulations without magnetic field we are able to recover the purely crustal
shear oscillations within an accuracy of about a few per cent. For dipole
magnetic fields between 5 x 10^13 and 10^15 G the Alfv\'en oscillations become
modified substantially by the presence of the crust. Those quasi-periodic
oscillations (QPOs) split into three families: Lower QPOs near the equator,
Edge QPOs related to the last open field line and Upper QPOs at larger distance
from the equator. Edge QPOs are called so because they are related to an edge
in the corresponding Alfv\'en continuum. The Upper QPOs are of the same kind,
while the Lower QPOs are turning-point QPOs, related to a turning point in the
continuous spectrum.
[48]
oai:arXiv.org:1007.0856 [pdf] - 1033516
Magneto-elastic oscillations and the damping of crustal shear modes in
magnetars
Submitted: 2010-07-06, last modified: 2010-11-24
In a realistic model of magneto-elastic oscillations in magnetars, we find
that crustal shear oscillations, often invoked as an explanation of
quasi-periodic oscillations (QPOs) seen after giant flares in soft gamma-ray
repeaters (SGRs), are damped by resonant absorption on timescales of at most
0.2s, for a lower limit on the dipole magnetic field strength of 5 x 10^13 G.
At higher magnetic field strengths (typical in magnetars) the damping timescale
is even shorter, as anticipated by earlier toy-models. We have investigated a
range of equations of state and masses and if magnetars are dominated by a
dipole magnetic field, our findings exclude torsional shear oscillations of the
crust from explaining the observed low-frequency QPOs. In contrast, we find
that the Alfv\'en QPO model is a viable explanation of observed QPOs, if the
dipole magnetic field strength exceeds a minimum strength of about several
times 10^14 G to 10^15 G. Then, Alfv\'en QPOs are no longer confined to the
fluid core, but completely dominate in the crust region and have a maximum
amplitude at the surface of the star.
[49]
oai:arXiv.org:1011.5475 [pdf] - 394308
Stationary, Axisymmetric Neutron Stars with Meridional Circulation in
General Relativity
Submitted: 2010-11-24
We present the first stationary, axisymmetric neutron star models with
meridional circulation in general relativity. For that purpose, we developed
GRNS, a new code based on a fixed point iteration. We find a two-dimensional
set of meridional circulation modes, which differ by the number of vortices in
the stream lines of the neutron star fluid. For expected maximal meridional
circulation velocities of about 1000 km/s, the vortices cause surface
deformations of about a percent. The deformations depend on the shape of the
vortices close to the surface and increase with the meridional circulation
velocity. We also computed models of rotating neutron stars with meridional
circulation, where neither the surface rotates nor does the rotation velocity
exceed the circulation velocity.
[50]
oai:arXiv.org:1011.3010 [pdf] - 321206
Stability of general-relativistic accretion disks
Submitted: 2010-11-12
Self-gravitating relativistic disks around black holes can form as transient
structures in a number of astrophysical scenarios such as binary neutron star
and black hole-neutron star coalescences, as well as the core-collapse of
massive stars. We explore the stability of such disks against runaway and
non-axisymmetric instabilities using three-dimensional hydrodynamics
simulations in full general relativity using the THOR code. We model the disk
matter using the ideal fluid approximation with a $\Gamma$-law equation of
state with $\Gamma=4/3$. We explore three disk models around non-rotating black
holes with disk-to-black hole mass ratios of 0.24, 0.17 and 0.11. Due to metric
blending in our initial data, all of our initial models contain an initial
axisymmetric perturbation which induces radial disk oscillations. Despite these
oscillations, our models do not develop the runaway instability during the
first several orbital periods. Instead, all of the models develop unstable
non-axisymmetric modes on a dynamical timescale. We observe two distinct types
of instabilities: the Papaloizou-Pringle and the so-called intermediate type
instabilities. The development of the non-axisymmetric mode with azimuthal
number m = 1 is accompanied by an outspiraling motion of the black hole, which
significantly amplifies the growth rate of the m = 1 mode in some cases.
Overall, our simulations show that the properties of the unstable
non-axisymmetric modes in our disk models are qualitatively similar to those in
Newtonian theory.
[51]
oai:arXiv.org:1003.0779 [pdf] - 262280
On the frequency band of the f-mode CFS instability
Submitted: 2010-03-03
Rapidly rotating neutron stars can be unstable to the
gravitational-wave-driven CFS mechanism if they have a neutral point in the
spectrum of nonaxisymmetric f-modes. We investigate the frequencies of these
modes in two sequences of uniformly rotating polytropes using nonlinear
simulations in full general relativity, determine the approximate locations of
the neutral points, and derive limits on the observable frequency band
available to the instability in these sequences. We find that general
relativity enhances the detectability of a CFS-unstable neutron star
substantially, both by widening the instability window and enlarging the band
into the optimal range for interferometric detectors like LIGO, VIRGO, and
GEO-600.
[52]
oai:arXiv.org:0902.1472 [pdf] - 1001435
Alfven QPOs in magnetars in the anelastic approximation
Submitted: 2009-02-09, last modified: 2009-06-19
We perform two-dimensional simulations of Alfven oscillations in magnetars,
modeled as relativistic stars with a dipolar magnetic field. We use the
anelastic approximation to general relativistic magnetohydrodynamics, which
allows for an effective suppression of fluid modes and an accurate description
of Alfven waves. In addition, we compute Alfven oscillation frequencies along
individual magnetic field lines with a semi-analytic approach, employing a
short-wavelength approximation. Our main findings are as follows: a) we confirm
the existence of two families of quasi-periodic oscillations (QPOs), with
harmonics at integer multiples of the fundamental frequency, as was found in
the linear study of Sotani, Kokkotas & Stergioulas; b) the QPOs appearing near
the magnetic axis are split into two groups, depending on their symmetry across
the equatorial plane. The antisymmetric QPOs have only odd integer-multiple
harmonics; c) the continuum obtained with our semi-analytic approach agrees
remarkably well with QPOs obtained via the two-dimensional simulations,
allowing for a clear interpretation of the QPOs as corresponding to turning
points of the continuum. This agreement will allow for a comprehensive study of
Alfven QPOs for a larger number of different models, without the need for
time-consuming simulations. Finally, we construct empirical relations for the
QPO frequencies and compare them to observations of known Soft Gamma Repeaters.
We find that, under the assumptions of our model and if the magnetic field of
magnetars is characterized by a strong dipolar component, and QPOs are produced
near the magnetic pole, then one can place an upper limit to the mean surface
strength of the magnetic field of about 3-8 10^15 G.
[53]
oai:arXiv.org:0710.1113 [pdf] - 5709
Alfv\'en QPOs in Magnetars
Submitted: 2007-10-04, last modified: 2007-11-19
We investigate torsional Alfv\'en oscillations of relativistic stars with a
global dipole magnetic field, via two-dimensional numerical simulations. We
find that a) there exist two families of quasi-periodic oscillations (QPOs)
with harmonics at integer multiples of the fundamental frequency, b) the
lower-frequency QPO is related to the region of closed field lines, near the
equator, while the higher-frequency QPO is generated near the magnetic axis, c)
the QPOs are long-lived, d) for the chosen form of dipolar magnetic field, the
frequency ratio of the lower to upper fundamental QPOs is ~0.6, independent of
the equilibrium model or of the strength of the magnetic field, and e) within a
representative sample of equations of state and of various magnetar masses, the
Alfv\'en QPO frequencies are given by accurate empirical relations that depend
only on the compactness of the star and on the magnetic field strength. The
lower and upper QPOs can be interpreted as corresponding to the edges or
turning points of an Alfv\'en continuum, according to the model proposed by
Levin (2007). Several of the low-frequency QPOs observed in the X-ray tail of
SGR 1806-20 can readily be identified with the Alfv\'en QPOs we compute. In
particular, one could identify the 18Hz and 30Hz observed frequencies with the
fundamental lower and upper QPOs, correspondingly, while the observed
frequencies of 92Hz and 150Hz are then integer multiples of the fundamental
upper QPO frequency (three times and five times, correspondingly). With this
identification, we obtain an upper limit on the strength of magnetic field of
SGR 1806-20 (if is dominated by a dipolar component) between ~3 and $7\times
10^{15}$G. Furthermore, we discuss the implications for the high-density EOS of
compact stars. (Abridged)
[54]
oai:arXiv.org:astro-ph/0608626 [pdf] - 84503
Torsional Oscillations of Relativistic Stars with Dipole Magnetic Fields
Submitted: 2006-08-29, last modified: 2007-10-11
We present the formalism and numerical results for torsional oscillations of
relativistic stars endowed with a strong dipole magnetic field. We do a
systematic search of parameter space by computing torsional mode frequencies
for various values of the harmonic index $\ell$ and for various overtones,
using an extended sample of models of compact stars, varying in mass,
high-density equation of state and crust model. We show that torsional mode
frequencies are sensitive to the crust model if the high-density equation of
state is very stiff. In addition, torsional mode frequencies are drastically
affected by a dipole magnetic field, if the latter has a strength exceeding
roughly $10^{15}$G and we find that the magnetic field effects are sensitive to
the adopted crust model. Using our extended numerical results we derive
empirical relations for the effect of the magnetic field on torsional modes as
well as for the crust thickness. We compare our numerical results to observed
frequencies in SGRs and find that certain high-density EoS and mass values are
favored over others in the non-magnetized limit. On the other hand, if the
magnetic field is strong, then its effect has to be taken into account in
attempts to formulate a theory of asteroseismology for magnetars.
[55]
oai:arXiv.org:0704.0431 [pdf] - 1513418
Fragmentation of general relativistic quasi-toroidal polytropes
Submitted: 2007-04-03
We investigate the role of rotational instabilities in the context of black
hole formation in relativistic stars. In addition to the standard scenario - an
axially symmetric dynamical instability forming a horizon at the star's center
- the recently found low-$T/|W|$ instabilities are shown to lead to
fragmentation and off-center horizon formation in differentially rotating
stars. This process might be an alternative pathway to produce SMBHs from
supermassive stars with inefficient angular momentum transport.
[56]
oai:arXiv.org:gr-qc/0702099 [pdf] - 112510
Gravitational Waves from Nonlinear Couplings of Radial and Polar
Nonradial Modes in Relativistic Stars
Submitted: 2007-02-17
The post-bounce oscillations of newly-born relativistic stars are expected to
lead to gravitational-wave emission through the excitation of nonradial
oscillation modes. At the same time, the star is oscillating in its radial
modes, with a central density variation that can reach several percent.
Nonlinear couplings between radial oscillations and polar nonradial modes lead
to the appearance of combination frequencies (sums and differences of the
linear mode frequencies). We study such combination frequencies using a
gauge-invariant perturbative formalism, which includes bilinear coupling terms
between different oscillation modes. For typical values of the energy stored in
each mode we find that gravitational waves emitted at combination frequencies
could become detectable in galactic core-collapse supernovae with advanced
interferometric or wide-band resonant detectors.
[57]
oai:arXiv.org:astro-ph/0611666 [pdf] - 87052
Torsional Oscillations of Relativistic Stars with Dipole Magnetic Fields
II. Global Alfv\'en Modes
Submitted: 2006-11-21
We investigate torsional Alfv\'{e}n modes of relativistic stars with a global
dipole magnetic field. It has been noted recently (Glampedakis et al. 2006)
that such oscillation modes could serve as as an alternative explanation (in
contrast to torsional crustal modes) for the SGR phenomenon, if the magnetic
field is not confined to the crust. We compute global Alfv\'{e}n modes for a
representative sample of equations of state and magnetar masses, in the ideal
MHD approximation and ignoring $\ell \pm 2$ terms in the eigenfunction. We find
that the presence of a realistic crust has a negligible effect on Alfv\'{e}n
modes for $B > 4\times 10^{15}$ G. Furthermore, we find strong avoided
crossings between torsional Alfv\'{e}n modes and torsional crust modes. For
magnetar-like magnetic field strengths, the spacing between consecutive
Alfv\'{e}n modes is of the same order as the gap of avoided crossings. As a
result, it is not possible to identify modes of predominantly crustal character
and all oscillations are predominantly Alfv\'{e}n-like. Interestingly, we find
excellent agreement between our computed frequencies and observed frequencies
in two SGRs, for a maximum magnetic field strenght in the range of
(0.8--1.2)$\times 10^{16}$ G.
[58]
oai:arXiv.org:astro-ph/0611601 [pdf] - 86987
Non-axisymmetric instability and fragmentation of general relativistic
quasi-toroidal stars
Submitted: 2006-11-19
In a recent publication, we have demonstrated that differentially rotating
stars admit new channels of black hole formation via fragmentation
instabilities. Since a higher order instability of this kind could potentially
transform a differentially rotating supermassive star into a multiple black
hole system embedded in a massive accretion disk, we investigate the dependence
of the instability on parameters of the equilibrium model. We find that many of
the models constructed exhibit non-axisymmetric instabilities with corotation
points, even for low values of T/|W|, which lead to a fission of the stars into
one, two or three fragments, depending on the initial perturbation.
At least in the models selected here, an m=1 mode becomes unstable at lower
values of T/|W|, which would seem to favor a scenario where one black hole with
a massive accretion disk forms. In this case, we have gained evidence that low
values of compactness of the initial model can lead to a stabilization of the
resulting fragment, thus preventing black hole formation in this scenario.
[59]
oai:arXiv.org:astro-ph/0511394 [pdf] - 77752
Non-linear axisymmetric pulsations of rotating relativistic stars in the
conformal flatness approximation
Submitted: 2005-11-14, last modified: 2006-03-30
We study non-linear axisymmetric pulsations of rotating relativistic stars
using a general relativistic hydrodynamics code under the assumption of a
conformal flatness. We compare our results to previous simulations where the
spacetime dynamics was neglected. The pulsations are studied along various
sequences of both uniformly and differentially rotating relativistic polytropes
with index N = 1. We identify several modes, including the lowest-order l = 0,
2, and 4 axisymmetric modes, as well as several axisymmetric inertial modes.
Differential rotation significantly lowers mode frequencies, increasing
prospects for detection by current gravitational wave interferometers. We
observe an extended avoided crossing between the l = 0 and l = 4 first
overtones, which is important for correctly identifying mode frequencies in
case of detection. For uniformly rotating stars near the mass-shedding limit,
we confirm the existence of the mass-shedding-induced damping of pulsations,
though the effect is not as strong as in the Cowling approximation. We also
investigate non-linear harmonics of the linear modes and notice that rotation
changes the pulsation frequencies in a way that would allow for various
parametric instabilities between two or three modes to take place. We assess
the detectability of each obtained mode by current gravitational wave detectors
and outline how the empirical relations that have been constructed for
gravitational wave asteroseismology could be extended to include the effects of
rotation.
[60]
oai:arXiv.org:gr-qc/0506083 [pdf] - 1592124
Gravitational Waves from Compact Sources
Submitted: 2005-06-15, last modified: 2005-07-06
We review sources of high-frequency gravitational waves, summarizing our
current understanding of emission mechanisms, expected amplitudes and event
rates. The most promising sources are gravitational collapse (formation of
black holes or neutron stars) and subsequent ringing of the compact star,
secular or dynamical rotational instabilities and high-mass compact objects
formed through the merger of binary neutron stars. Significant and unique
information for the various stages of the collapse, the structure of
protoneutron stars and the high density equation of state of compact objects
can be drawn from careful study of gravitational wave signals.
[61]
oai:arXiv.org:gr-qc/0501080 [pdf] - 262567
Black hole formation through fragmentation of toroidal polytropes
Submitted: 2005-01-26, last modified: 2005-02-10
We investigate new paths to black hole formation by considering the general
relativistic evolution of a differentially rotating polytrope with toroidal
shape. We find that this polytrope is unstable to nonaxisymmetric modes, which
leads to a fragmentation into self-gravitating, collapsing components. In the
case of one such fragment, we apply a simplified adaptive mesh refinement
technique to follow the evolution to the formation of an apparent horizon
centered on the fragment. This is the first study of the one-armed instability
in full general relativity.
[62]
oai:arXiv.org:gr-qc/0403029 [pdf] - 419237
Three-dimensional relativistic simulations of rotating neutron-star
collapse to a Kerr black hole
Submitted: 2004-03-06, last modified: 2005-01-30
We present a new three-dimensional fully general-relativistic hydrodynamics
code using high-resolution shock-capturing techniques and a conformal traceless
formulation of the Einstein equations. Besides presenting a thorough set of
tests which the code passes with very high accuracy, we discuss its application
to the study of the gravitational collapse of uniformly rotating neutron stars
to Kerr black holes. The initial stellar models are modelled as relativistic
polytropes which are either secularly or dynamically unstable and with angular
velocities which range from slow rotation to the mass-shedding limit. We
investigate the gravitational collapse by carefully studying not only the
dynamics of the matter, but also that of the trapped surfaces, i.e. of both the
apparent and event horizons formed during the collapse. The use of these
surfaces, together with the dynamical horizon framework, allows for a precise
measurement of the black-hole mass and spin. The ability to successfully
perform these simulations for sufficiently long times relies on excising a
region of the computational domain which includes the singularity and is within
the apparent horizon. The dynamics of the collapsing matter is strongly
influenced by the initial amount of angular momentum in the progenitor star
and, for initial models with sufficiently high angular velocities, the collapse
can lead to the formation of an unstable disc in differential rotation. All the
simulations performed with uniformly rotating initial data and a polytropic or
ideal-fluid equation of state show no evidence of shocks or of the presence of
matter on stable orbits outside the black hole.
[63]
oai:arXiv.org:astro-ph/0312648 [pdf] - 61860
Nonlinear Pulsations in Differentially Rotating Neutron Stars:
Mass-Shedding-Induced Damping and Splitting of the Fundamental Mode
Submitted: 2003-12-30, last modified: 2004-05-04
We study small-amplitude, nonlinear pulsations of uniformly and
differentially rotating neutron stars employing a two-dimensional evolution
code for general-relativistic hydrodynamics. Using Fourier transforms at
several points inside the star, both the eigenfrequencies and two-dimensional
eigenfunctions of pulsations are extracted. The centrifugal forces and the
degree of differential rotation have significant effects on the
mode-eigenfunction. We find that near the mass-shedding limit, the pulsations
are damped due to shocks forming at the surface of the star. This new damping
mechanism may set a small saturation amplitude for modes that are unstable to
the emission of gravitational waves. After correcting for the assumption of the
Cowling approximation (used in our numerical code), we construct empirical
relations that predict the range of gravitational-wave frequencies from
quasi-periodic post-bounce oscillations in the core collapse of massive stars.
We also find that the fundamental quasi-radial mode is split, at least in the
Cowling approximation and mainly in differentially rotating stars, into two
different sequences.
[64]
oai:arXiv.org:gr-qc/0310061 [pdf] - 111360
Approximate Matching of Analytic and Numerical Solutions for Rapidly
Rotating Neutron Stars
Submitted: 2003-10-13, last modified: 2004-02-16
We investigate the properties of a closed-form analytic solution recently
found by Manko et al. (2000) for the exterior spacetime of rapidly rotating
neutron stars. For selected equations of state we numerically solve the full
Einstein equations to determine the neutron star spacetime along constant rest
mass sequences. The analytic solution is then matched to the numerical
solutions by imposing the condition that the quadrupole moment of the numerical
and analytic spacetimes be the same. For the analytic solution we consider,
such a matching condition can be satisfied only for very rapidly rotating
stars. When solutions to the matching condition exist, they belong to one of
two branches. For one branch the current octupole moment of the analytic
solution is very close to the current octupole moment of the numerical
spacetime; the other branch is more similar to the Kerr solution. We present an
extensive comparison of the radii of innermost stable circular orbits (ISCOs)
obtained with a) the analytic solution, b) the Kerr metric, c) an analytic
series expansion derived by Shibata and Sasaki (1998) and d) a highly accurate
numerical code. In most cases where a corotating ISCO exists, the analytic
solution has an accuracy consistently better than the Shibata-Sasaki expansion.
The numerical code is used for tabulating the mass-quadrupole and
current-octupole moments for several sequences of constant rest mass.
[65]
oai:arXiv.org:astro-ph/0308035 [pdf] - 58314
Non-linear resonance in the accretion disk of a millisecond pulsar
Submitted: 2003-08-02
Twin quasi-periodic millisecond modulations of the X-ray flux (kHz QPOs) have
recently been reported from an accreting 2.5 ms X-ray pulsar. We identify modes
of disk oscillations whose frequencies are in agreement with the observed ones
when the rotating neutron star is modeled with realistic equations of state.
The frequency difference of the twin QPOs, equal to about one half of the
neutron-star spin rate, clearly indicates that resonant oscillations of the
accretion disk have been observed. Similar non-linear resonances may also be
spontaneously excited in the accretion disk. The two QPO frequencies in the
pulsar system are close to a 5:7 ratio and this suggests a link with the QPOs
in black hole systems, where frequency ratios of 2:3 and 3:5 have been
reported.
[66]
oai:arXiv.org:gr-qc/0302034 [pdf] - 1237415
Rotating Stars in Relativity
Submitted: 2003-02-10
Rotating relativistic stars have been studied extensively in recent years,
both theoretically and observationally, because of the information one could
obtain about the equation of state of matter at extremely high densities and
because they are considered to be promising sources of gravitational waves. The
latest theoretical understanding of rotating stars in relativity is reviewed in
this updated article. The sections on the equilibrium properties and on the
nonaxisymmetric instabilities in f-modes and r-modes have been updated and
several new sections have been added on analytic solutions for the exterior
spacetime, rotating stars in LMXBs, rotating strange stars, and on rotating
stars in numerical relativity.
[67]
oai:arXiv.org:astro-ph/0204380 [pdf] - 48925
Spin-down of Relativistic Stars with Phase Transitions and PSR
J0537-6910
Submitted: 2002-04-23, last modified: 2003-02-07
Using a highly accurate numerical code, we study the spin down of rotating
relativistic stars, undergoing a quark deconfinement phase transition. Such
phase transitions have been suggested to yield an observable signal in the
braking index of spinning-down pulsars, which is based on a ``backbending''
behaviour of the moment of inertia. We focus on a particular equation of state
that has been used before to study this behaviour, and find that for the
population of normal pulsars the moment of inertia does not exhibit a
backbending behaviour. In contrast, for supramassive millisecond pulsars a very
strong backbending behaviour is found. Essentially, once a quark core appears
in a spinning-down supramassive millisecond pulsar, the star spins up and
continues to do so until it reaches the instability to collapse. This strong
spin-up behaviour makes it easier to distinguish a phase transition in such
pulsars: a negative first time-derivative of the rotational period suffices and
one does not have to measure the braking index. In the spin-up era, the usually
adopted spin-down power law fails to describe the evolution of the angular
velocity. We adopt a general-relativistic spin-down power law and derive the
equations that describe the angular velocity and braking index evolution in
rapidly rotating pulsars.
[68]
oai:arXiv.org:astro-ph/0206470 [pdf] - 50101
An unusually low mass of some "neutron" stars?
Submitted: 2002-06-27
The X-ray emission of RXJ1856.5-3754 has been found to coincide to
unprecedented accuracy with that of a blackbody, of radius 5.8+-0.9 km for the
measured parallax distance of 140 pc (Burwitz et al. 2001, Drake et al. 2002).
If the emission is uniform over the whole surface of a non-rotating star, the
mass of the star cannot exceed 0.75+-0.12 solar mass regardless of its
composition. If the compact object is a quark star described by the MIT-bag
equation of state (a ``strange star''), the mass is no more than 0.3 solar
mass. Comparably small masses are also obtained for the X-ray bursters Aql X-1
and KS1731-260 for some fits to their spectra.
[69]
oai:arXiv.org:astro-ph/0112495 [pdf] - 46868
Time evolution of the linear perturbations of a rotating Newtonian
polytrope
Submitted: 2001-12-20
We present the results of numerical time evolutions of the linearised
perturbations of rapidly and rigidly rotating Newtonian polytropes while making
the Cowling approximation. The evolution code runs stably for hundreds of
stellar rotations, allowing us to compare our results with previously published
eigenmode calculations, for instance the f-mode calculations of Ipser &
Lindblom, and the r-mode calculations of Karino et al. The mode frequencies
were found to be in agreement within the expected accuracy. We have also
examined the inertial modes recently computed by Lockitch & Friedman, and we
were able to extend their slow-rotation results into the rapid rotation regime.
In the longer term, this code will provide a platform for studying a number of
poorly understood problems in stellar oscillation theory, such as the effect of
differential rotation and gravitational radiation reaction on normal mode
oscillations and, with suitable modifications, mode-mode coupling in the mildly
non-linear regime.
[70]
oai:arXiv.org:gr-qc/0110047 [pdf] - 110857
Three-dimensional general relativistic hydrodynamics II: long-term
dynamics of single relativistic stars
Submitted: 2001-10-09
This is the second in a series of papers on the construction and validation
of a three-dimensional code for the solution of the coupled system of the
Einstein equations and of the general relativistic hydrodynamic equations, and
on the application of this code to problems in general relativistic
astrophysics. In particular, we report on the accuracy of our code in the
long-term dynamical evolution of relativistic stars and on some new physics
results obtained in the process of code testing. The tests involve single
non-rotating stars in stable equilibrium, non-rotating stars undergoing radial
and quadrupolar oscillations, non-rotating stars on the unstable branch of the
equilibrium configurations migrating to the stable branch, non-rotating stars
undergoing gravitational collapse to a black hole, and rapidly rotating stars
in stable equilibrium and undergoing quasi-radial oscillations. The numerical
evolutions have been carried out in full general relativity using different
types of polytropic equations of state using either the rest-mass density only,
or the rest-mass density and the internal energy as independent variables. New
variants of the spacetime evolution and new high resolution shock capturing
(HRSC) treatments based on Riemann solvers and slope limiters have been
implemented and the results compared with those obtained from previous methods.
Finally, we have obtained the first eigenfrequencies of rotating stars in full
general relativity and rapid rotation. A long standing problem, such
frequencies have not been obtained by other methods. Overall, and to the best
of our knowledge, the results presented in this paper represent the most
accurate long-term three-dimensional evolutions of relativistic stars available
to date.
[71]
oai:arXiv.org:astro-ph/0110209 [pdf] - 358195
Lower limits on the maximum orbital frequency around rotating strange
stars
Submitted: 2001-10-09
Observations of kHz quasi-periodic oscillations (QPOs) in the X-ray fluxes of
low-mass X-ray binaries (LMXBs) have been used in attempts to constrain the
external metric of the compact members of these binaries, as well as their
masses and the equation of state of matter at supranuclear denisties. We
compute the maximum orbital frequency of stable circular motion around
uniformly rotating strange stars described by the MIT bag model. The
calculations are performed for both normal and supramassive constant baryon
mass sequences of strange stars rotating at all possible rates. We find the
lower limits on the maximum orbital frequency and discuss them for a range of
masses and for all rotational frequencies allowed in the model considered. We
show that for slowly and moderately rotating strange stars the maximum value of
orbital frequency can be a good indicator of the mass of the compact object.
However, for rapidly rotating strange stars the same value of orbital frequency
in the innermost stable circular orbit is obtained for stars with masses
ranging from that of a planetoid to about three solar masses. At sufficiently
high rotation rates of the strange star, the rotational period alone constrains
the stellar mass to a surprisingly narrow range.
[72]
oai:arXiv.org:astro-ph/0105175 [pdf] - 42404
Torsional Oscillations of Magnetized Relativistic Stars
Submitted: 2001-05-10, last modified: 2001-09-19
Strong magnetic fields in relativistic stars can be a cause of crust
fracturing, resulting in the excitation of global torsional oscillations. Such
oscillations could become observable in gravitational waves or in high-energy
radiation, thus becoming a tool for probing the equation of state of
relativistic stars. As the eigenfrequency of torsional oscillation modes is
affected by the presence of a strong magnetic field, we study torsional modes
in magnetized relativistic stars. We derive the linearized perturbation
equations that govern torsional oscillations coupled to the oscillations of a
magnetic field, when variations in the metric are neglected (Cowling
approximation). The oscillations are described by a single two-dimensional wave
equation, which can be solved as a boundary value problem to obtain
eigenfrequencies. We find that in the non-magnetized case, typical oscillation
periods of the fundamental l=2 torsional modes can be nearly a factor of two
larger for relativistic stars than previously computed in the Newtonian limit.
For magnetized stars, we show that the influence of the magnetic field is
highly dependent on the assumed magnetic field configuration and simple
estimates obtained previously in the literature cannot be used for identifying
normal modes observationally.
[73]
oai:arXiv.org:astro-ph/0107043 [pdf] - 880529
Fast Magnetosonic Waves Driven by Gravitational Waves
Submitted: 2001-07-03
The propagation of a gravitational wave (GW) through a magnetized plasma is
considered. In particular, we study the excitation of fast magnetosonic waves
(MSW) by a gravitational wave, using the linearized general-relativistic
hydromagnetic equations. We derive the dispersion relation for the plasma,
treating the gravitational wave as a perturbation in a Minkowski background
space-time. We show that the presence of gravitational waves will drive
magnetosonic waves in the plasma and discuss the potential astrophysical
implications.
[74]
oai:arXiv.org:astro-ph/0012477 [pdf] - 40066
Axisymmetric Modes of Rotating Relativistic Stars in the Cowling
Approximation
Submitted: 2000-12-22, last modified: 2001-04-06
Axisymmetric pulsations of rotating neutron stars can be excited in several
scenarios, such as core-collapse, crust and core-quakes and binary mergers and
could become detectable either in gravitational waves or high-energy radiation.
Here, we present a comprehensive study of all low-order axisymmetric modes of
uniformly and rapidly rotating relativistic stars. Initial stationary
configurations are appropriately perturbed and are numerically evolved using an
axisymmetric, nonlinear relativistic hydrodynamics code, assuming
time-independence of the gravitational field (Cowling approximation). The
simulations are performed using a high-resolution shock-capturing
finite-difference scheme accurate enough to maintain the initial rotation law
for a large number of rotational periods, even for stars at the mass-shedding
limit. Through Fourier transforms of the time evolution of selected fluid
variables, we compute the frequencies of quasi-radial and non-radial modes with
spherical harmonic indices l=0,1,2 and 3, for a sequence of rotating stars from
the non-rotating limit to the mass-shedding limit. The frequencies of the
axisymmetric modes are affected significantly by rotation only when the
rotation rate exceeds about 50% of the maximum allowed. As expected, at large
rotation rates, apparent mode crossings between different modes appear. In
addition to the above modes, several axisymmetric inertial modes are also
excited in our numerical evolutions.
[75]
oai:arXiv.org:gr-qc/0007086 [pdf] - 898900
Nonlinear r-modes in Rapidly Rotating Relativistic Stars
Submitted: 2000-07-31, last modified: 2000-12-08
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.
[76]
oai:arXiv.org:astro-ph/0011083 [pdf] - 39090
Nonlinear evolution of r-modes in rotating relativistic stars
Submitted: 2000-11-03
A numerical study of nonlinear r-modes in isentropic, rapidly rotating
relativistic stars, via 3-D general-relativistic hydrodynamical evolutions, is
presented. On dynamical timescales, we find no evidence for strong coupling of
r-modes to other modes at amplitudes of order one or larger. Therefore, unless
nonlinear saturation sets in on longer timescales, the maximum r-mode amplitude
is of order one. An absolute upper limit on the amplitude is set by causality.
Our simulations also show that r-modes and inertial modes in isentropic stars
are discrete modes, with no evidence for the existence of a continuous part in
the spectrum.
[77]
oai:arXiv.org:astro-ph/0009244 [pdf] - 38086
Evolutionary Self-Energy-Loss Effects in Compact Binary Systems:
Importance of Rapid Rotation and of Equation of State
Submitted: 2000-09-15
The spin-down of a millisecond pulsar in a compact binary system leads to
self-energy losses, that cause the binary's orbital period to increase (the
effect being of first post-Newtonian (1PN) order). If the pulsar's
period-derivative is not exceedingly small, this effect can become measurable
in long-term high-precision timing measurements of the binary's orbital motion.
We use rotating compressible spheroids to obtain an approximate, explicit
expression for the orbital period derivative, due to spin-down, valid to 1PN
order. This expression can be used to observationally constrain the pulsar's
moment of inertia, and, combined with other observations, the high-density
equation of state of compact stars. We apply our expression to representative
models of millisecond pulsars in binary systems and demonstrate the importance
of including higher-order rotational effects as well as the importance of the
choice of equation of state, which both have been neglected, so far, in the
literature. The computed increase in orbital period is larger, by as much as
32%, when the higher-order rotational effects are included, while it can change
by more than a factor of two between different plausible equations of state.
[78]
oai:arXiv.org:gr-qc/0003071 [pdf] - 110515
Towards a Stable Numerical Evolution of Strongly Gravitating Systems in
General Relativity: The Conformal Treatments
Submitted: 2000-03-17, last modified: 2000-05-11
We study the stability of three-dimensional numerical evolutions of the
Einstein equations, comparing the standard ADM formulation to variations on a
family of formulations that separate out the conformal and traceless parts of
the system. We develop an implementation of the conformal-traceless (CT)
approach that has improved stability properties in evolving weak and strong
gravitational fields, and for both vacuum and spacetimes with active coupling
to matter sources. Cases studied include weak and strong gravitational wave
packets, black holes, boson stars and neutron stars. We show under what
conditions the CT approach gives better results in 3D numerical evolutions
compared to the ADM formulation. In particular, we show that our implementation
of the CT approach gives more long term stable evolutions than ADM in all the
cases studied, but is less accurate in the short term for the range of
resolutions used in our 3D simulations.
[79]
oai:arXiv.org:astro-ph/0002114 [pdf] - 34457
R-mode runaway and rapidly rotating neutron stars
Submitted: 2000-02-04, last modified: 2000-03-14
We present a simple spin evolution model that predicts that rapidly rotating
accreting neutron stars will mainly be confined to a narrow range of
spin-frequencies; P= 1.5-5 ms. This is in agreement with current observations
of both neutron stars in the Low-Mass X-ray Binaries and millisecond
radiopulsars. The main ingredients in the model are; i) the instability of
r-modes above a critical spin rate, ii) thermal runaway due to heat released as
viscous damping mechanisms counteract the r-mode growth, and iii) a revised
estimate of the strength of dissipation due to the presence of a viscous
boundary layer at the base of the crust in an old and relatively cold neutron
star. We discuss the gravitational waves that are radiated during the brief
r-mode driven spin-down phase. We also briefly touch on how the new estimates
affects the predicted initial spin periods of hot young neutron stars.
[80]
oai:arXiv.org:astro-ph/9909152 [pdf] - 1469823
Keplerian frequencies and innermost stable circular orbits of rapidly
rotating strange stars
Submitted: 1999-09-09
It has been suggested that the frequency in the co-rotating innermost stable
circular orbit (ISCO) about a compact stellar remnant can be determined through
X-ray observations of low-mass X-ray binaries, and that its value can be used
to constrain the equation of state of ultradense matter.
Upon constructing numerical models of rapidly rotating strange (quark) stars
in general relativity, we find that for stars rotating at the equatorial
mass-shedding limit, the ISCO is indeed above the stellar surface, for a wide
range of central energy densities at a height equal to 11% of the
circumferential stellar radius, which scales inversely with the square root of
the energy density, of self-bound quark matter at zero presure. In contrast to
static stars, the ISCO frequencies of rapidly rotating strange stars can be as
low as 0.9 kHz for a 1.3 solar mass strange star. Hence, the presence of
strange stars in low-mass X-ray binaries cannot be excluded on the basis of the
currently observed frequencies of kHz QPOs, such as the cut-off frequency of
1066 Hz in 4U 1820-30.
[81]
oai:arXiv.org:gr-qc/9908010 [pdf] - 113228
Nonlinear hydrodynamical evolution of rotating relativistic stars:
Numerical methods and code tests
Submitted: 1999-08-03
We present numerical hydrodynamical evolutions of rapidly rotating
relativistic stars, using an axisymmetric, nonlinear relativistic hydrodynamics
code. We use four different high-resolution shock-capturing (HRSC)
finite-difference schemes (based on approximate Riemann solvers) and compare
their accuracy in preserving uniformly rotating stationary initial
configurations in long-term evolutions. Among these four schemes, we find that
the third-order PPM scheme is superior in maintaining the initial rotation law
in long-term evolutions, especially near the surface of the star. It is further
shown that HRSC schemes are suitable for the evolution of perturbed neutron
stars and for the accurate identification (via Fourier transforms) of normal
modes of oscillation. This is demonstrated for radial and quadrupolar
pulsations in the nonrotating limit, where we find good agreement with
frequencies obtained with a linear perturbation code. The code can be used for
studying small-amplitude or nonlinear pulsations of differentially rotating
neutron stars, while our present results serve as testbed computations for
three-dimensional general-relativistic evolution codes.
[82]
oai:arXiv.org:gr-qc/9904009 [pdf] - 113157
Non-linear Evolution of Rotating Relativistic Stars
Submitted: 1999-04-05
We present first results of the non-linear evolution of rotating relativistic
stars obtained with an axisymmetric relativistic hydrodynamics code in a fixed
spacetime. As initial data we use stationary axisymmetric and perturbed
configurations. We find that, in order to prevent (numerical) angular momentum
loss at the surface layers of the star a high-resolution grid (or a numerical
scheme that retains high order at local extrema) is needed. For non-rotating
stars, we compute frequencies of radial and non-radial small-amplitude
oscillations, which are in excellent agreement with linear normal mode
frequencies computed in the Cowling approximation. As a first application of
our code, quasi-radial modes of rapidly rotating relativistic stars are
computed. By generalizing our numerical code to 3-D, we plan to study the
evolution and non-linear dynamics of toroidal oscillations (r-modes) of rapidly
rotating neutron stars, which are a promising source of gravitational waves.
[83]
oai:arXiv.org:astro-ph/9806089 [pdf] - 101682
On the relevance of the r-mode instability for accreting neutron stars
and white dwarfs
Submitted: 1998-06-05, last modified: 1998-12-02
We present a case study for the relevance of the r-mode instability for
accreting compact stars. Our estimates are based on approximations that
facilitate back-of-the-envelope calculations. We discuss two different cases:
1) For recycled millisecond pulsars we argue that the r-mode instability may
be active at rotation periods longer than the Kepler period (which provides the
dynamical limit on rotation) as long as the core temperature is larger than
about $2\times10^5$ K. Our estimates suggest that the instability may have
played a role in the evolution of the fastest spinning pulsars, and that it may
be presently active in the recently discovered 2.49 ms X-ray pulsar SAX
J1808.4-3658 as well as the rapidly spinning neutron stars observed in low mass
X-ray binaries. This provides a new explanation for the remarkably similar
rotation periods inferred from kHz quasi-periodic oscillations in the LMXBs.
The possibility that the rotation of recycled pulsars may be
gravitational-radiation limited is interesting because the gravitational waves
from a neutron star rotating at the instability limit may well be detectable
with the new generation of interferometric detectors. 2) We also consider white
dwarfs, and find that the r-mode instability may possibly be active in short
period white dwarfs. Our order of magnitude estimates (for a white dwarf of
$M=M_\odot$ and $R=0.01 R_\odot$ composed of $C^{12}$) show that the
instability could be operating for rotational periods shorter than $P\approx
27-33$ s. This number is in interesting agreement with the observed periods
($>28$ s) of the rapidly spinning DQ Herculis stars. However, we find that the
instability grows too slowly to affect the rotation of these stars
significantly .
[84]
oai:arXiv.org:gr-qc/9805012 [pdf] - 1498651
Rotating Stars in Relativity
Submitted: 1998-05-05, last modified: 1998-07-07
Rotating relativistic stars are receiving significant attention in recent
years, because of the information they can yield about the equation of state of
matter at extremely high densities and because they are one of the more
possible sources of detectable gravitational waves. We review the latest
theoretical and numerical methods for modeling rotating relativistic stars,
including stars with a strong magnetic field and hot proto-neutron stars. We
also review nonaxisymmetric oscillations and instabilities in rotating stars
and summarize the latest developments regarding the gravitational wave-driven
(CFS) instability in both polar and axial quasi-normal modes.
[85]
oai:arXiv.org:gr-qc/9806008 [pdf] - 898989
Quasi-normal modes of rotating relativistic stars - neutral modes for
realistic equations of state
Submitted: 1998-06-01
We compute zero-frequency (neutral) quasi-normal f-modes of fully
relativistic and rapidly rotating neutron stars, using several realistic
equations of state (EOSs) for neutron star matter. The zero-frequency modes
signal the onset of the gravitational radiation-driven instability. We find
that the l=m=2 (bar) f-mode is unstable for stars with gravitational mass as
low as 1.0 - 1.2 M_\odot, depending on the EOS. For 1.4 M_\odot neutron stars,
the bar mode becomes unstable at 83 % - 93 % of the maximum allowed rotation
rate. For a wide range of EOSs, the bar mode becomes unstable at a ratio of
rotational to gravitational energies T/W \sim 0.07-0.09 for 1.4 M_\odot stars
and T/W \sim 0.06 for maximum mass stars. This is to be contrasted with the
Newtonian value of T/W \sim 0.14. We construct the following empirical formula
for the critical value of T/W for the bar mode,
(T/W)_2 = 0.115 - 0.048 M / M_{max}^{sph}, which is insensitive to the EOS to
within 4 - 6 %. This formula yields an estimate for the neutral mode sequence
of the bar mode as a function only of the star's mass, M, given the maximum
allowed mass, M_{max}^{sph}, of a nonrotating neutron star. The recent
discovery of the fast millisecond pulsar in the supernova remnant N157B,
supports the suggestion that a fraction of proto-neutron stars are born in a
supernova collapse with very large initial angular momentum. Thus, in a
fraction of newly born neutron stars the instability is a promising source of
continuous gravitational waves. It could also play a major role in the
rotational evolution (through the emission of angular momentum) of merged
binary neutron stars, if their post-merger angular momentum exceeds the maximum
allowed to form a Kerr black hole.
[86]
oai:arXiv.org:astro-ph/9805297 [pdf] - 363213
Analytic description of the r-mode instability in uniform density stars
Submitted: 1998-05-22
We present an analytic description of the $r$-mode instability in newly-born
neutron stars, using the approximation of uniform density. Our computation is
consistently accurate to second order in the angular velocity of the star. We
obtain formulae for the growth-time of the instability due to
gravitational-wave emission, for both current and mass multipole radiation and
for the damping timescale, due to viscosity. The $l=m=2$ current-multipole
radiation dominates the timescale of the instability. We estimate the deviation
of the second order accurate results from the lowest order approximation and
show that the uncertainty in the equation of state has only a small effect on
the onset of the $r$-mode instability. The viscosity coefficients and the
cooling process in newly-born neutron stars are, at present, uncertain and our
analytic formaulae enables a quick check of such effects on the development of
the instability.
[87]
oai:arXiv.org:gr-qc/9804048 [pdf] - 363341
Construction of Highly Accurate Models of Rotating Neutron Stars -
Comparison of Three Different Numerical Schemes
Submitted: 1998-04-21
We conduct a direct comparison of three different representative numerical
codes for constructing models of rapidly rotating neutron stars in general
relativity. Our aim is to evaluate the accuracy of the codes and to investigate
how the accuracy is affected by the choice of interpolation, domain of
integration and equation of state. In all three codes, the same physical
parameters, equations of state and interpolation method are used. We construct
25 selected models for polytropic equations of state and 22 models with
realistic neutron star matter equations of state. The three codes agree well
with each other (typical agreement is better than 0.1% to 0.01%) for most
models, except for the extreme assumption of uniform density stars. We conclude
that the codes can be used for the construction of highly accurate initial data
configurations for polytropes of index N > 0.5 (which typically correspond to
realistic neutron stars), when the domain of integration includes all space and
for realistic equations with no phase transitions. With the exception of the
uniform density case, the obtained values of physical parameters for the models
considered in this paper can be regarded as ``standard'' and we display them in
detail for all models.
[88]
oai:arXiv.org:gr-qc/9705056 [pdf] - 112848
Nonaxisymmetric Neutral Modes in Rotating Relativistic Stars
Submitted: 1997-05-21
We study nonaxisymmetric perturbations of rotating relativistic stars.
modeled as perfect-fluid equilibria. Instability to a mode with angular
dependence $\exp(im\phi)$ sets in when the frequency of the mode vanishes. The
locations of these zero-frequency modes along sequences of rotating stars are
computed in the framework of general relativity. We consider models of
uniformly rotating stars with polytropic equations of state, finding that the
relativistic models are unstable to nonaxisymmetric modes at significantly
smaller values of rotation than in the Newtonian limit. Most strikingly, the
m=2 bar mode can become unstable even for soft polytropes of index $N \leq
1.3$, while in Newtonian theory it becomes unstable only for stiff polytropes
of index $N \leq 0.808$. If rapidly rotating neutron stars are formed by the
accretion-induced collapse of white dwarfs, instability associated with these
nonaxisymmetric, gravitational-wave driven modes may set an upper limit on
neutron-star rotation. Consideration is restricted to perturbations that
correspond to polar perturbations of a spherical star. A study of axial
perturbations is in progress.
[89]
oai:arXiv.org:astro-ph/9608179 [pdf] - 95297
Upper Limit Set by Causality on the Rotation and Mass of Uniformly
Rotating Relativistic Stars
Submitted: 1996-08-27
Causality alone suffices to set a lower bound on the period of rotation of
relativistic stars as a function of their maximum observed mass. That is, by
assuming a one-parameter equation of state (EOS) that satisfies v_sound < c and
that allows stars with masses as large as the largest observed neutron-star
mass, M_sph^max, we find P[ms] > 0.282 + 0.196 ( M_sph^max/M_odot-1.442). The
limit does not assume that the EOS agrees with a known low-density form for
ordinary matter, but if one adds that assumption, the minimum period is raised
by a few percent. Thus the current minimum period of uniformly rotating stars,
set by causality, is 0.28ms (0.29ms for stars with normal crust). The
minimizing EOS yields models with a maximally soft exterior supported by a
maximally stiff core. An analogous upper limit set by causality on the maximum
mass of rotating neutron stars requires a low-density match and the limit
depends on the matching density, epsilon_m. We recompute it, obtaining a
slightly revised value, M_rot^max \sim 6.1( 2 * 10^14 g/cm^3 epsilon_m )^1/2
M_odot.
[90]
oai:arXiv.org:astro-ph/9411032 [pdf] - 142697
Comparing Models of Rapidly Rotating Relativistic Stars Constructed by
Two Numerical Methods
Submitted: 1994-11-08
We present the first direct comparison of codes based on two different
numerical methods for constructing rapidly rotating relativistic stars. A code
based on the Komatsu-Eriguchi-Hachisu (KEH) method (Komatsu et al. 1989),
written by Stergioulas, is compared to the Butterworth-Ipser code (BI), as
modified by Friedman, Ipser and Parker. We compare models obtained by each
method and evaluate the accuracy and efficiency of the two codes. The agreement
is surprisingly good. A relatively large discrepancy recently reported
(Eriguchi et al. 1994) is found to arise from the use of two different versions
of the equation of state. We find, for a given equation of state, that
equilibrium models with maximum values of mass, baryon mass, and angular
momentum are (generically) all distinct and either all unstable to collapse or
are all stable. Our implementation of the KEH method will be available as a
public domain program for interested users.
