Source author record

Kostas D. Kokkotas

Kostas D. Kokkotas appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

56works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

56 published item(s)

preprint2022arXiv

Charged black holes in de Sitter space: superradiant amplification of charged scalar waves and resonant hyperradiation

Black holes typically inhabit highly dynamical galactic environments and are frequently permeated by accretion media. The inevitable scattering of scalar, electromagnetic and gravitational waves off rotating and charged black holes provides a remarkable demonstration of the energy extraction and subsequent amplification of scattered waves under the expense of the compact object's rotational and electromagnetic energy; a phenomenon known as superradiance. Certain circumstances tremendously favor the energy extraction process in such extent that linear superradiant instabilities are triggered. We examine the superradiant amplification of monochromatic charged scalar fields impinging Reissner-Nordström-de Sitter black holes which are known to exhibit superradiantly unstable quasinormal modes. We find that even though the frequency range of superradiance is reduced when a positive cosmological constant is incorporated, the amplification factors are significantly elevated with respect to those occurring in Reissner-Nordström spacetime. Intriguingly, we confirm that the long-lived quasinormal resonances that reside in the superradiant regime induce resonant peaks when the wave's frequency matches the real part of the quasinormal mode, regardless of the spacetime's modal stability.

preprint2022arXiv

Dynamical behavior of the $C$-metric: charged scalar fields, quasinormal modes and superradiance

The $C$-metric is a boost-symmetric spacetime solution to the vacuum Einstein field equations which describes black holes that are uniformly accelerated under the tension of a cosmic string. Only recently the thermodynamics of accelerating black holes and their modal stability against neutral scalar perturbations were concisely established. The generalization of accelerating black holes to incorporate an electric charge, namely the charged $C$-metric, possess three distinct families of quasinormal mode frequencies; the complex photon surface quasinormal modes associated with unstable null particles at the equatorial plane of the photon surface, the purely imaginary acceleration modes whose existence solely depends on the acceleration of spacetime and the purely imaginary near-extremal modes which dominate the dynamics of the ringdown at late times when the event and Cauchy horizon approach each other. We extend the quasinormal mode analysis to charged scalar fluctuations and find that the photon surface modes are continuously deformed with respect to their neutral counterpart as the scalar charge is increased. We further find that the acceleration and near-extremal families acquire an oscillation frequency when the scalar charge is introduced. Finally, we study the superradiant amplification of charged scalar monochromatic waves impinging charged accelerating black holes. We find that even though the frequency range of superradiant amplification is lessened due to the acceleration, the amplification factors are considerably elevated with respect to those transpiring in Reissner-Nordström black holes and are maximized when the scalar charge is significantly large.

preprint2022arXiv

GRMHD beyond Kerr: An extension of the HARM code for thin disks to non-Kerr spacetimes

Black hole based tests of general relativity have proliferated in recent times with new and improved detectors and telescopes. Modelling of the black hole neighborhood, where most of the radiation carrying strong-field signature originates, is of utmost importance for robust and accurate constraints on possible violations of general relativity. As a first step, this paper presents the extension of general relativistic magnetohydrodynamic simulations of thin accretion disks to parametrically deformed black holes that generalize the Kerr solution. The extension is based on \textsc{harmpi}, a publicly available member of the \textsc{harm} family of codes, and uses a phenomenological metric to study parametric deviations away from Kerr. The extended model is used to study the disk structure, stability, and radiative efficiency. We also compute the Fe K$α$ profiles in simplified scenarios and present an outlook for the future.

preprint2022arXiv

On Black Hole Area Quantization and Echoes

In this work we argue that black hole area quantization of Bekenstein and Mukhanov should not give rise to measurable effects in terms of so-called gravitational wave echoes during black hole mergers. We outline that the quantum spectrum of a black hole should be washed out during and after black hole mergers, and hence one should not expect echoes in this scenario. The extreme broadening of the spectrum is due to the large particle emission rate during ringdown. Our results question key assumptions being made in recent literature on this topic.

preprint2022arXiv

Quasi-periodic oscillations in precursor flares via seismic aftershocks from resonant shattering

GRB 211211A was accompanied by a kilonova, implying a merger origin for the event. A precursor flare, modulated by quasi-periodic oscillations at $\approx 22\, \mbox{Hz}$, was also seen $\approx 1\, \mbox{s}$ prior to the main emission. It is suggested here that the precursor resulted from the resonant shattering of one star's crust prior to coalescence. Seismic aftershocks and low-frequency torsional modes may emanate from the rupture site in this case, explaining the underlying oscillations. This interpretation is directly analogous to proposals for giant flares from soft gamma repeaters, where oscillations at similar frequencies have been observed, involving starquakes followed by crustal vibrations. Torsional mode properties were computed for sequences of slowly rotating, magnetised neutron stars in general relativity. The $\approx 22\, \mbox{Hz}$ modulations in the precursor matches that of nodeless, $\ell =2$ torsional modes for a wide variety of stellar parameters. Global or discrete Alfvén modes are also viable candidates.

preprint2022arXiv

Slowly-rotating curved acoustic black holes: Quasinormal modes, Hawking-Unruh radiation and quasibound states

Astrophysical black holes are generally surrounded by accretion disks, galactic matter and the omnipresent cosmic microwave background radiation, thus allowing for the concurrent propagation of both gravitational and sound waves. Recently, acoustic black holes were embedded in Schwarzschild spacetime allowing for the coexistence of event and acoustic horizons. Here, we obtain a class of perturbative solutions to the field equations of the relativistic Gross-Pitaevskii and Yang-Mills theories, which describe sound waves propagating on a curved slowly-rotating acoustic black hole, akin to Lense-Thirring spacetime. We investigate the quasinormal mode frequencies, Hawking-Unruh radiation, and quasibound states. Our novel metric mimics the gravitational field of astrophysical compact objects in the limiting case of slow rotation, and therefore could, in principle, shed more light into the underlying classical and quantum physics of black holes through analog acoustic probes.

preprint2021arXiv

Precessing magnetars as central engines in short gamma-ray bursts

Short gamma-ray bursts that are followed by long-duration X-ray plateaus may be powered by the birth, and hydrodynamic evolution, of magnetars from compact binary coalescence events. If the rotation and magnetic axes of the system are not orthogonal to each other, the star will undergo free precession, leading to fluctuations in the luminosity of the source. In some cases, precession-induced modulations in the spin-down power may be discernible in the X-ray flux of the plateau. In this work, 25 X-ray light curves associated with bursts exhibiting a plateau are fitted to luminosity profiles appropriate for precessing, oblique rotators. Based on the Akaike Information Criterion, 16 (64 per cent) of the magnetars within the sample display either moderate or strong evidence for precession. Additionally, since the precession period of the star is directly tied to its quadrupolar ellipticity, the fits allow for an independent measure of the extent to which the star is deformed by internal stresses. Assuming these deformations arise due to a mixed poloidal-toroidal magnetic field, we find that the distribution of magnetic-energy ratios is bimodal, with data points clustering around energetically equal and toroidally dominated partitions. Implications of this result for gravitational-wave emission and dynamo activity in newborn magnetars are discussed.

preprint2020arXiv

Evidence for magnetar precession in X-ray afterglows of gamma-ray bursts

Many gamma-ray bursts are followed by periods of extended emission. At least in some cases, the burst afterglow may be powered by a rapidly rotating, highly-magnetised neutron star, which spins down due to electromagnetic and gravitational wave emission. Such a remnant is likely to strongly precess in the early stages of its life, which would lead to modulations in the X-ray luminosity as the triaxiality of the system evolves over time. Using a radiation profile appropriate for a precessing, oblique rotator, we find that Swift-XRT data of a long (080602) and a short (090510) burst matches the model with significantly higher accuracy (mean-square residuals dropping by $\gtrsim 200 \%$ in the early stages of the extended emission) than for an orthogonal rotator. We interpret this as evidence for precession in newborn magnetars.

preprint2020arXiv

Precursor flares of short gamma-ray bursts from crust yielding due to tidal resonances in coalescing binaries of rotating, magnetized neutron stars

As evidenced by the coincident detections of GW170817 and GRB 170817A, short gamma-ray bursts are likely associated with neutron star-neutron star merger events. Although rare, some bursts display episodes of early emission, with precursor flares being observed up to $\sim 10$ seconds prior to the main burst. As the stars inspiral due to gravitational wave emission, the exertion of mutual tidal forces leads to the excitation of stellar oscillation modes, which may come into resonance with the orbital motion. Mode amplitudes increase substantially during a period of resonance as tidal energy is deposited into the star. The neutron star crust experiences shear stress due to the oscillations and, if the resonant amplitudes are large enough, may become over-strained. This over-straining can lead to fractures or quakes which release energy, thereby fueling precursor activity prior to the merger. Using some simple Maclaurin spheroid models, we investigate the influence of magnetic fields and rapid rotation on tidally-forced $f$- and $r$- modes, and connect the associated eigenfrequencies with the orbital frequencies corresponding to precursor events seen in, for example, GRB 090510.

preprint2020arXiv

Stability of topological neutron stars

Tensor-multi-scalar theories (TMST) are among the most natural generalizations of Einstein's theory, they are mathematically self-consistent and free from pathologies. They pass through all the known observations but contrary to standard scalar-tensor theories, TMST offer extremely rich spectrum of solutions and allow for large deviations from general relativity. One of the most interesting objects in these theories are the topological neutron stars. The main goal in the present paper is to study their radial stability since many branches of solutions can exist even for a fixed value of the topological charge. It turns out that only one of these branches is stable for each value of the topological charge while all the rest are unstable. The stable branch is exactly the one that spans from small to large neutron star masses having as well moderate values of the radii in agreement with the observations. As expected, it becomes unstable at the maximum of the mass. The frequencies of the radial modes are examined and it turns out that in most cases the mode frequencies are larger than the general relativistic ones and they increase with the increase of the topological charge. All this shows that the topological neutron stars are viable astrophysical objects that should be explored further in order to determine their observational manifestations.

preprint2020arXiv

Testing spacetime symmetry through gravitational waves from extreme-mass-ratio inspirals

One of the primary aims of upcoming space-borne gravitational wave detectors is to measure radiation in the mHz range from extreme-mass-ratio inspirals. Such a detection would place strong constraints on hypothetical departures from a Kerr description for astrophysically stable black holes. The Kerr geometry, which is unique in general relativity, admits a higher-order symmetry in the form of a Carter constant, which implies that the equations of motion describing test particle motion in a Kerr background are Liouville-integrable. In this article, we investigate whether the Carter symmetry itself is discernible from a generic deformation of the Kerr metric in the gravitational waveforms for such inspirals. We build on previous studies by constructing a new metric which respects current observational constraints, describes a black hole, and contains two non-Kerr parameters, one of which controls the presence or absence of the Carter symmetry, thereby controlling the existence of chaotic orbits, and another which serves as a generic deformation parameter. We find that these two parameters introduce fundamentally distinct features into the orbital dynamics, and evince themselves in the gravitational waveforms through a significant dephasing. Although only explored in the quadrupole approximation, this, together with a Fisher metric analysis, suggests that gravitational wave data analysis may be able to test, in addition to the governing theory of gravity, the underlying symmetries of spacetime.

preprint2020arXiv

The gravitational Higgs mechanism and resulting smoking gun effects

Recently, a toy model was introduced to demonstrate that screening mechanisms in alternative theories of gravitation can hide additional effects. In this model a scalar field is charged under a $U(1)$ symmetry. In sufficiently compact objects the scalar field spontaneously grows, i.e. the object scalarizes, spontaneously breaking the $U(1)$ symmetry. Exactly as in the $U(1)$ Higgs mechanism this leads to the emergence of a mass for the gauge field. The aim of this paper is to provide an example of the physical consequences if we consider this toy model as a prototype of Weak Equivalence Principle (WEP) violations. We model neutron stars with a dipolar magnetic field to compare the magnetic field behaviour of stars in Einstein-Maxwell theory on the one hand and in scalar-tensor theory with the, so-called, gravitational Higgs mechanism on the other hand.

preprint2016arXiv

Saturation of the f-mode instability in neutron stars: II. Applications and results

We present the first results on the saturation of the f-mode instability in neutron stars, due to nonlinear mode coupling. Emission of gravitational waves drives the f-mode (fundamental mode) unstable in fast-rotating, newborn neutron stars. The initial growth phase of the mode is followed by its saturation, because of energy leaking to other modes of the star. The saturation point determines the strain of the generated gravitational-wave signal, which can then be used to extract information about the neutron star equation of state. The parent (unstable) mode couples via parametric resonances with pairs of daughter modes, with the triplets' evolution exhibiting a rich variety of behaviors. We study both supernova- and merger-derived neutron stars, simply modeled as polytropes in a Newtonian context, and show that the parent may couple to many different daughter pairs during the star's evolution through the instability window, with the saturation amplitude changing by orders of magnitude.

preprint2015arXiv

A gravitational wave afterglow in binary neutron star mergers

We study in detail the f-mode secular instability for rapidly rotating neutron stars, putting emphasis on supermassive models which do not have a stable nonrotating counterpart. Such neutron stars are thought to be the generic outcome of the merger of two standard mass neutron stars. In addition we take into account the effects of strong magnetic field and r-mode instability, that can drain a substantial amount of angular momentum. We find that the gravitational wave signal emitted by supramassive neutron stars can reach above the Advance LIGO sensitivity at distance of about 20Mpc and the detectability is substantially enhanced for the Einstein Telescope. The event rate will be of the same order as the merging rates, while the analysis of the signal will carry information for the equation of state of the post-merging neutron stars and the strength of the magnetic fields.

preprint2015arXiv

Asteroseismology of rapidly rotating neutron stars - an alternative approach

In the present paper we examine gravitational wave asteroseismology relations for f-modes of rapidly rotating neutron stars. An approach different to the previous studies is employed - first, the moment of inertia is used instead of the stellar radius, and second, the normalization of the oscillation frequencies and damping times is different. It was shown that in the non-rotating case this can lead to a much stronger equation of state independence and our goal is to generalize the static relations to the rapidly rotating case and values of the spherical mode number $l\ge2$. We employ realistic equations of state that cover a very large range of stiffness in order to check better the universality of the relations. At the end we explore the inverse problem, i.e. obtain the neutron star parameters from the observed gravitational frequencies and damping times. It turns out that with this new set of relations we can solve the inverse problem with a very good accuracy using three frequencies that was not possible in the previous studies where one needs also the damping times. The asteroseismology relations are also particularly good for the massive rapidly rotating models that are subject to secular instabilities.

preprint2015arXiv

Gravitational wave asteroseismology of neutron and strange stars in $R^2$ gravity

We study the oscillations of neutron and strange stars in $R^2$ gravity. More precisely the nonradial $f$-modes are examined and the differences with pure general relativity are investigated. Using these results we build several gravitational wave asteroseismology relations. Our goal is to determine up to what extend these relations are equation of state independent and whether they deviate enough from general relativity in order to produce an observable effect. The results show that the differences coming from $R^2$ gravity are up to 10\% and that will be difficult to be observed in the near future. On the other hand the small deviations in some of the asteroseismology relations show that they are not only equation of state independent, but they are also quite insensitive to the gravitational theory. That is why solving the inverse problem can give us quite robust estimates of the neutron star parameters.

preprint2015arXiv

R-mode astronomy

Next generation gravitational wave detectors will start taking data in the near future. Here we discuss the chances to detect the continuous emission from r-mode oscillations in compact stars and study which properties of compact stars we can infer from such novel data. In particular we show that the combination of the gravitational wave data with electromagnetic multi-messenger observations could give us detailed insight into compact star properties, ranging from precise mass-radius measurements to the determination of the equation of state and the phase structure of dense matter.

preprint2015arXiv

Radiation from charged particles on eccentric orbits in a dipolar magnetic field around a Schwarzschild black hole

We obtain an approximate solution for the motion of a charged particle around a Schwarzschild black hole immersed in a weak dipolar magnetic field. We focus on eccentric bound orbits in the equatorial plane of the Schwarzschild black hole and derive an analytic expression for the spectral distribution of the electromagnetic emission from a charged particle on such an orbit. Two sets of harmonic contributions appear, with specific frequency spacing. The expression can be written in compact form, if it is truncated up to the lowest order harmonic contributions.

preprint2015arXiv

Rapidly rotating neutron stars in $R$-squared gravity

$f(R)$ theories of gravity are one of the most popular alternative explanations for dark energy and therefore studying the possible astrophysical implications of these theories is an important task. In the present paper we make a substantial advance in this direction by considering rapidly rotating neutron stars in $R^2$ gravity. The results are obtained numerically and the method we use is non-perturbative and self-consistent. The neutron star properties, such as mass, radius and moment of inertia, are studied in detail and the results show that rotation magnifies the deviations from general relativity and the maximum mass and moment of inertia can reach very high values. This observation is similar to previous studies of rapidly rotating neutron stars in other alternative theories of gravity, such as the scalar-tensor theories, and it can potentially lead to strong astrophysical manifestations.

preprint2015arXiv

Saturation of the f-mode instability in neutron stars: I. Theoretical framework

The basic formulation describing quadratic mode coupling in rotating Newtonian stars is presented, focusing on polar modes. Due to the Chandrasekhar-Friedman-Schutz mechanism, the f-mode (fundamental oscillation) is driven unstable by the emission of gravitational waves. If the star falls inside the so-called instability window, the mode's amplitude grows exponentially, until it is halted by nonlinear effects. Quadratic perturbations form three-mode networks inside the star, which evolve as coupled oscillators, exchanging energy. Coupling of the unstable f-mode to other (stable) modes can lead to a parametric resonance and the subsequent saturation of its amplitude, thus suppressing the instability. The saturation point determines the amplitude of the gravitational-wave signal obtained from an individual source, as well as the evolutionary path of the latter inside the instability window.

preprint2015arXiv

The I-Q relations for rapidly rotating neutron stars in $f(R)$ gravity

In the present paper we study the behavior of the normalized $I$-$Q$ relation for neutron stars in a particular class of $f(R)$ theories of gravity, namely the $R^2$ gravity that is one of the most natural and simplest extensions of general relativity in the strong field regime. We study both the slowly and rapidly rotating cases. The results show that the $I$-$Q$ relation remain nearly equation of state independent for fixed values of the normalized rotational parameter, but the deviations from universality can be a little bit larger compared to the general relativistic case. What is the most interesting in our studies, is that the differences with the pure Einstein's theory can be large reaching above 20\%. This is qualitative different from the majority of alternative theories of gravity, where the normalized $I$-$Q$ relations are almost indistinguishable from the general relativistic case, and can lead to observational constraints on the $f(R)$ theories in the future.

preprint2015arXiv

The stochastic background of gravitational waves due to the f-mode instability in neutron stars

This paper presents an estimate for the spectral properties of the stochastic background of gravitational waves emitted by a population of hot, young, rapidly rotating neutron stars throughout the Universe undergoing $f$-mode instabilities, formed through either core-collapse supernova explosions or the merger of binary neutron star systems. Their formation rate, from which the gravitational wave event rate is obtained, is deduced from observation-based determinations of the cosmic star formation rate. The gravitational wave emission occurs during the spin-down phase of the $f$-mode instability. For low magnetized neutron stars and assuming 10\% of supernova events lead to $f$-mode unstable neutron stars, the background from supernova-derived neutron stars peaks at $Ω_{\text{gw}} \sim 10^{-9}$ for the $l=m=2$ $f$-mode, which should be detectable by cross-correlating a pair of second generation interferometers (e.g. Advanced LIGO/Virgo) with an upper estimate for the signal-to-noise ratio of $\approx$ 9.8. The background from supramassive neutron stars formed from binary mergers peaks at $Ω_{\text{gw}} \sim 10^{-10}$ and should not be detectable, even with third generation interferometers (e.g. Einstein Telescope).

preprint2014arXiv

Breakdown of I-Love-Q universality in rapidly rotating relativistic stars

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.

preprint2014arXiv

Electromagnetic waves from neutron stars and black holes driven by polar gravitational perturbations

Neutron stars and black holes are the most compact astrophysical objects we can think of and as a consequence they are the main sources of gravitational waves. There are many astrophysically relevant scenarios in which these objects are immersed in or endowed with strong magnetic fields, in such a way that gravitational perturbations can couple to electromagnetic ones and can potentially trigger synergistic electromagnetic signatures. In a recent paper we derived the main equations for gravito-electromagnetic perturbations and studied in detail the case of polar electromagnetic perturbations driven by axial gravitational perturbations. In this paper we deal with the case of axial electromagnetic perturbations driven by polar black-hole or neutron stars oscillations, in which the energy emitted in case is considerably larger than in the previous case. In the case of neutron stars the phenomenon lasts considerably longer since the fluid acts as an energy reservoir that shakes the magnetic field for a timescale of the order of secs.

preprint2014arXiv

Non-perturbative and self-consistent models of neutron stars in R-squared gravity

In the present paper we investigate non-perturbatively and self-consistently the structure of neutron stars in $R$-squared gravity by simultaneously solving the interior and exterior problem. The mass-radius relations are obtained for several equations of state and for wide range of the $R$-squared gravity parameter $a$. Even though the deviation from general relativity for nonzero values of $a$ can be large, they are still comparable with the variations due to different modern realistic equations of state. That is why the current observations of the neutron star masses and radii alone can not put constraints on the value of the parameter $a$. We also compare our results with those obtained within the perturbative method and we discuss the differences between them.

preprint2014arXiv

Orbital and epicyclic frequencies around rapidly rotating compact stars in scalar-tensor theories of gravity

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.

preprint2014arXiv

Rapidly rotating neutron stars in scalar-tensor theories of gravity

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 $β>-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.

preprint2014arXiv

Slowly rotating neutron and strange stars in $R^2$ gravity

In the present paper we investigate self-consistently slowly rotating neutron and strange stars in R-squared gravity. For this purpose we first derive the equations describing the structure of the slowly rotating compact stars in $f(R)$-gravity and then simultaneously solve the exterior and the interior problem. The structure of the slowly rotating neutron stars is studied for two different hadronic equations of state and a strange matter equation of state. The moment of inertia and its dependence on the stellar mass and the $R$-squared gravity parameter $a$ is also examined in details. We find that the neutron star moment of inertia for large values of the parameter $a$ can be up to $30\%$ larger compared to the corresponding general relativistic models. This is much higher than the change in the maximum mass induced by $R$-squared gravity and is beyond the EOS uncertainty. In this way the future observations of the moment of inertia of compact stars could allow us to distinguish between general relativity and $f(R)$ gravity, and more generally to test the strong field regime of gravity.

preprint2014arXiv

Universal I-Q relations for rapidly rotating neutron and strange stars in scalar-tensor theories

We study how rapid rotation influences the relation between the normalized moment of inertia $\bar{I}$ and quadrupole moment $\bar{Q}$ for scalarized neutron stars. The questions one has to answer are whether the EOS universality is preserved in this regime and what are the deviations from general relativity. Our results show that the $\bar{I}-\bar{Q}$ relation is nearly EOS independent for scalarized rapidly rotating stars, but the differences with pure Einstein's theory increase compared to the slowly rotating case. In general, smaller negative values of the scalar field coupling parameters $β$ lead to larger deviations, but these deviations are below the expected accuracy of the future astrophysical observations if one considers values of $β$ in agreement with the current observational constraint. An important remark is that although the normalized $\bar{I}-\bar{Q}$ relation is quite similar for scalar-tensor theories and general relativity, the unnormalized moment of inertia and quadrupole moment can be very different in the two theories. This demonstrates that although the $\bar{I}-\bar{Q}$ relations are potentially very useful for some purposes, they might not serve us well when trying to distinguish between different theories of gravity.

preprint2013arXiv

Gravitational wave asteroseismology of fast rotating neutron stars with realistic equations of state

In the present paper we study the oscillations of fast rotating neutron stars with realistic equations of state (EoS) within the Cowling approximation. We derive improved empirical relations for gravitational wave asteroseismology with f-modes and for the first time we consider not only quadrupolar oscillations but also modes with higher spherical order (l=|m|=3,4). After performing a systematic comparison with polytropic EoS, it is shown that the empirical relations found in this case approximately also hold for realistic EoS. Even more, we show that these relations will not change significantly even if the Cowling approximation is dropped and the full general relativistic case is considered, although the normalization used here(frequencies and damping times in the nonrotating limit) could differ considerably. We also address the inverse problem, i.e. we investigate in detail what kind of observational data is required in order to determine characteristical neutron star parameters. It is shown that masses, radii and rotation rates can be estimated quite accurately using the derived asteroseismology relations. We also compute the instability window for certain models, i.e. the limiting curve in a T-Ω-plane where the secular Chandrasekhar-Friedman-Schutz(CFS) instability overcomes dissipative effects, and show that some of the modern realistic EoS will lead to a larger instability window compared to all of the polytropic ones presented so far in the literature. Additionally, we calculate the r-mode instability window and compare it with the f-mode-case. The overall results for the instability window suggest that it is vital to take into account oscillations with l=3,4 when considering gravitational wave asteroseismology using the f-mode in rapidly rotating neutron stars, as these modes can become CFS unstable for a much larger range of parameters than pure quadrupolar oscillations.

preprint2013arXiv

Gravitationally Driven Electromagnetic Perturbations of Neutron Stars and Black Holes

Gravitational perturbations of neutron stars and black holes are well known sources of gravitational radiation. If the compact object is immersed in or endowed with a magnetic field, the gravitational perturbations would couple to electromagnetic perturbations and potentially trigger synergistic electromagnetic signatures. We present a detailed analytic calculation of the dynamics of coupled gravitational and electromagnetic perturbations for both neutron stars and black holes. We discuss the prospects for detecting the electromagnetic waves in these scenarios and the potential that these waves have for providing information about their source.

preprint2013arXiv

Investigating the retention of intermediate-mass black holes in star clusters using N-body simulations

Contrary to supermassive and stellar-mass black holes (SBHs), the existence of intermediate-mass black holes (IMBHs) with masses ranging between 10^{2-5} Msun has not yet been confirmed. The main problem in the detection is that the innermost stellar kinematics of globular clusters (GCs) or small galaxies, the possible natural loci to IMBHs, are very difficult to resolve. However, if IMBHs reside in the centre of GCs, a possibility is that they interact dynamically with their environment. A binary formed with the IMBH and a compact object of the GC would naturally lead to a prominent source of gravitational radiation, detectable with future observatories. We use N-body simulations to study the evolution of GCs containing an IMBH and calculate the gravitational radiation emitted from dynamically formed IMBH-SBH binaries and the possibility that the IMBH escapes the GC after an IMBH-SBH merger. We run for the first time direct-summation integrations of GCs with an IMBH including the dynamical evolution of the IMBH with the stellar system and relativistic effects, such as energy loss in gravitational waves (GWs) and periapsis shift, and gravitational recoil. We find in one of our models an intermediate mass-ratio inspiral (IMRI), which leads to a merger with a recoiling velocity higher than the escape velocity of the GC. The GWs emitted fall in the range of frequencies that a LISA-like observatory could detect, like the European eLISA or in mission options considered in the recent preliminary mission study conducted in China. The merger has an impact on the global dynamics of the cluster, as an important heating source is removed when the merged system leaves the GC. The detection of one IMRI would constitute a test of GR, as well as an irrefutable proof of the existence of IMBHs.

preprint2012arXiv

Gravitational Waves and Hydromagnetic Instabilities in Rotating Magnetized Neutron Stars

We perform nonlinear general relativistic ideal magnetohydrodynamic simulations of poloidal magnetic fields in rotating polytropic neutron stars. We have three primary goals: i) to understand the nature of magnetohydrodynamic instabilities inherent to poloidal magnetic fields in non-rotating and rotating neutron stars, ii) to explore the possible space of stable equilibrium configurations and iii) to understand gravitational wave emissions caused by the catastrophic reconfiguration of magnetic fields associated with giant magnetar flares. Our key physical contributions can be summarized as follows: i) gravitational waves from f-modes caused by magnetar flares are unlikely to be detected in the current or near-future generation of gravitational waves observatories, ii) gravitational waves from Alfven waves propagating inside the neutron star are more likely candidates, although this interpretation relies on the unknown damping time of these modes, iii) any magnetic field equilibria derived from our simulations are characterized as non-axisymmetric, with approximately 65% of their magnetic energy in the poloidal field, iv) rotation acts to separate the timescales of different instabilities in our system, with the varicose mode playing a more major role due to a delayed kink instability and v) despite the slowing growth rate of the kink mode, it is always present in our simulations, even for models where the rotational period is of the same order as the Alfven timescale.

preprint2011arXiv

Are gravitational waves from giant magnetar flares observable?

Are giant flares in magnetars viable sources of gravitational radiation? Few theoretical studies have been concerned with this problem, with the small number using either highly idealized models or assuming a magnetic field orders of magnitude beyond what is supported by observations. We perform nonlinear general-relativistic magnetohydrodynamics simulations of large-scale hydromagnetic instabilities in magnetar models. We utilise these models to find gravitational wave emissions over a wide range of energies, from 10^40 to 10^47 erg. This allows us to derive a systematic relationship between the surface field strength and the gravitational wave strain, which we find to be highly nonlinear. In particular, for typical magnetar fields of a few times 10^15 G, we conclude that a direct observation of f-modes excited by global magnetic field reconfigurations is unlikely with present or near-future gravitational wave observatories, though we also discuss the possibility that modes in a low-frequency band up to 100 Hz could be sufficiently excited to be relevant for observation.

preprint2011arXiv

Gravitational wave asteroseismology with fast rotating neutron stars

We investigate damping and growth times of the f-mode for rapidly rotating stars and a variety of different polytropic equations of state in the Cowling approximation. We discuss the differences in the eigenfunctions of co- and counterrotating modes and compute the damping times of the f-mode for several EoS and all rotation rates up to the Kepler-limit. This is the first study of the damping/growth time of this type of oscillations for fast rotating neutron stars in a general relativistic framework. We use these frequencies and damping/growth times to create robust empirical formulae which can be used for gravitational wave asteroseismology. The estimation of the damping/growth time is based on the quadrupole formula and our results agree very well with Newtonian ones in the appropriate limit.

preprint2011arXiv

Hydromagnetic Instabilities in Neutron Stars

We model the non-linear ideal magnetohydrodynamics of poloidal magnetic fields in neutron stars in general relativity assuming a polytropic equation of state. We identify familiar hydromagnetic modes, in particular the 'sausage/varicose' mode and 'kink' instability inherent to poloidal magnetic fields. The evolution is dominated by the kink instability, which causes a cataclysmic reconfiguration of the magnetic field. The system subsequently evolves to new, non-axisymmetric, quasi-equilibrium end-states. The existence of this branch of stable quasi-equilibria may have consequences for magnetar physics, including flare generation mechanisms and interpretations of quasi-periodic oscillations.

preprint2011arXiv

Magnetars oscillations in the presence of a crust

We study axisymmetric perturbations of neutron star endowed with a strong magnetic field (magnetars), considering the coupled oscillations of the fluid core with the solid crust. We recover discrete oscillations based mainly in the crust and a continuum in the core. We also confirm the presence of "discrete Alfvén modes" in the gap between two contiguous continua (see Van Hooven& Levin (2011)) and, in addition, we can resolve some of them also inside the continua. Our results can explain both the lower and the higher observed quasi periodical oscillations (QPOs) in SGR 1806-20 and SGR 1900+14 and put constrains on the mass, radius and crust thickness of the two magnetars.

preprint2011arXiv

On the saturation amplitude of the f-mode instability

We investigate strong nonlinear damping effects which occur during high amplitude oscillations of neutron stars, and the gravitational waves they produce. For this, we use a general relativistic nonlinear hydrodynamics code in conjunction with a fixed spacetime (Cowling approximation) and a polytropic equation of state (EOS). Gravitational waves are estimated using the quadrupole formula. Our main interest are l=m=2 f modes subject to the CFS (Chandrasekhar, Friedman, Schutz) instability, but we also investigate axisymmetric and quasiradial modes. We study various models to determine the influence of rotation rate and EOS. We find that axisymmetric oscillations at high amplitudes are predominantly damped by shock formation, while the nonaxisymmetric f modes are mainly damped by wave breaking and, for rapidly rotating models, coupling to nonaxisymmetric inertial modes. From the observed nonlinear damping, we derive upper limits for the saturation amplitude of CFS-unstable f modes. Finally, we estimate that the corresponding gravitational waves for an oscillation amplitude at the upper limit should be detectable with the advanced LIGO and VIRGO interferometers at distances above 10 MPc. This strongly depends on the stellar model, in particular on the mode frequency.

preprint2011arXiv

The f-mode instability in relativistic neutron stars

Rapidly spinning neutron stars are known to harbour pulsation modes that may become unstable and grow in amplitude by emitting gravitational radiation. Among the various stellar modes, the f-mode is the one typically considered as a promising source of gravitational radiation for ground-based detectors such as LIGO and VIRGO. Improving the existing work in Newtonian stellar models, we present the first calculation of the basic properties of the f-mode instability in rapidly rotating relativistic neutron stars, adopting the Cowling approximation. Using a relativistic polytropic stellar model, we obtain a minimum gravitational growth timescale (for the dominant l=m=4 mode) of the order of 10^3-10^4 s near the Kepler spin frequency Omega_K, which is substantially shorter than the Newtonian value. By accounting for dissipation in neutron star matter, i.e. shear/bulk viscosity and superfluid mutual friction, we calculate the associated f-mode instability window. For our specific stellar model, the instability is active above 0.92 \times Omega_K and for temperatures \sim (10^9 - 2 \times 10^{10}) K, characteristic of newborn neutron stars.

preprint2010arXiv

High-order perturbations of a spherical collapsing star

In Ref. [1, 2] a formalism to deal with high-order perturbations of a general spherical background was developed. In this article, we apply it to the particular case of a perfect fluid background. We have expressed the perturbations of the energy-momentum tensor at any order in terms of the perturbed fluid's pressure, density and velocity. In general, these expressions are not linear and have sources depending on lower order perturbations. For the second-order case we make the explicit decomposition of these sources in tensor spherical harmonics. Then, a general procedure is given to evolve the perturbative equations of motions of the perfect fluid for any value of the harmonic label. Finally, with the problem of a spherical collapsing star in mind, we discuss the high-order perturbative matching conditions across a timelike surface, in particular the surface separating the perfect fluid interior from the exterior vacuum.

preprint2010arXiv

MYRIAD: A new N-body code for simulations of Star Clusters

We present a new C++ code for collisional N-body simulations of star clusters. The code uses the Hermite fourth-order scheme with block time steps, for advancing the particles in time, while the forces and neighboring particles are computed using the GRAPE-6 board. Special treatment is used for close encounters, binary and multiple sub-systems that either form dynamically or exist in the initial configuration. The structure of the code is modular and allows the appropriate treatment of more physical phenomena, such as stellar and binary evolution, stellar collisions and evolution of close black-hole binaries. Moreover, it can be easily modified so that the part of the code that uses GRAPE-6, could be replaced by another module that uses other accelerating-hardware like the Graphics Processing Units (GPUs). Appropriate choice of the free parameters give a good accuracy and speed for simulations of star clusters up to and beyond core collapse. Simulations of Plummer models consisting of equal-mass stars reached core collapse at t~17 half-mass relaxation times, which compares very well with existing results, while the cumulative relative error in the energy remained below 0.001. Also, comparisons with published results of other codes for the time of core collapse for different initial conditions, show excellent agreement. Simulations of King models with an initial mass-function, similar to those found in the literature, reached core collapse at t~0.17, which is slightly smaller than the expected result from previous works. Finally, the code accuracy becomes comparable and even better than the accuracy of existing codes, when a number of close binary systems is dynamically created in a simulation. This is due to the high accuracy of the method that is used for close binary and multiple sub-systems.

preprint2010arXiv

Oscillations and instabilities of fast and differentially rotating relativistic stars

We study non-axisymmetric oscillations of rapidly and differentially rotating relativistic stars in the Cowling approximation. Our equilibrium models are sequences of relativistic polytropes, where the differential rotation is described by the relativistic $j$-constant law. We show that a small degree of differential rotation raises the critical rotation value for which the quadrupolar f-mode becomes prone to the CFS instability, while the critical value of $T/|W|$ at the mass-shedding limit is raised even more. For softer equations of state these effects are even more pronounced. When increasing differential rotation further to a high degree, the neutral point of the CFS instability first reaches a local maximum and is lowered afterwards. For stars with a rather high compactness we find that for a high degree of differential rotation the absolute value of the critical $T/|W|$ is below the corresponding value for rigid rotation. We conclude that the parameter space where the CFS instability is able to drive the neutron star unstable is increased for a small degree of differential rotation and for a large degree at least in stars with a higher compactness.

preprint2009arXiv

Relativistic g-modes in rapidly rotating neutron stars

We study the g-modes of fast rotating stratified neutron stars in the general relativistic Cowling approximation, where we neglect metric perturbations and where the background models take into account the buoyant force due to composition gradients. This is the first paper studying this problem in a general relativistic framework. In a recent paper by Passamonti et al.(2009), a similar study was performed within the Newtonian framework, where the authors presented results about the onset of CFS-unstable g-modes and the close connection between inertial- and gravity-modes for sufficiently high rotation rates and small composition gradients. This correlation arises from the interplay between the buoyant force which is the restoring force for g-modes and the Coriolis force which is responsible for the existence of inertial modes. In our relativistic treatment of the problem, we find an excellent qualitatively agreement with respect to the Newtonian results.

preprint2008arXiv

Gravitational Wave Astronomy

As several large scale interferometers are beginning to take data at sensitivities where astrophysical sources are predicted, the direct detection of gravitational waves may well be imminent. This would open the gravitational-wave window to our Universe, and should lead to a much improved understanding of the most violent processes imaginable; the formation of black holes and neutron stars following core collapse supernovae and the merger of compact objects at the end of binary inspiral.

preprint2003arXiv

High frequency sources of gravitational waves

Sources of high frequency gravitational waves are reviewed. Gravitational collapse, rotational instabilities and oscillations of the remnant compact objects are potentially important sources of gravitational waves. Significant and unique information for the various stages of the collapse, the evolution of protoneutron stars and the details of the equations of state of such objects can be drawn from careful study of the gravitational wave signal.

preprint2003arXiv

Stability of five-dimensional rotating black holes projected on the brane

We study the stability of five-dimensional Myers-Perry black holes with a single angular momentum under linear perturbations, and we compute the quasinormal modes (QNM's) of the black hole metric projected on the brane, using Leaver's continued fraction method. In our numerical search we do not find unstable modes. The damping time of modes having l=m=2 and l=m=1 tends to infinity as the black hole spin tends to the extremal value, showing a behaviour reminiscent of the one observed for ordinary 4-dimensional Kerr black holes.

preprint2001arXiv

On the r-mode spectrum of relativistic stars: Inclusion of the radiation reaction

We consider both mode calculations and time evolutions of axial r-modes for relativistic uniformly rotating non-barotropic neutron stars, using the slow-rotation formalism, in which rotational corrections are considered up to linear order in the angular velocity Ω. We study various stellar models, such as uniform density models, polytropic models with different polytropic indices n, and some models based on realistic equations of state. For weakly relativistic uniform density models, and polytropes with small values of n, we can recover the growth times predicted from Newtonian theory when standard multipole formulae for the gravitational radiation are used. However, for more compact models, we find that relativistic linear perturbation theory predicts a weakening of the instability compared to the Newtonian results. When turning to polytropic equations of state, we find that for certain ranges of the polytropic index n, the r-mode disappears, and instead of a growth, the time evolutions show a rapid decay of the amplitude. This is clearly at variance with the Newtonian predictions. It is, however, fully consistent with our previous results obtained in the low-frequency approximation.

preprint1999arXiv

Construction of Astrophysical Initial Data for Perturbations of Relativistic Stars

We develop a framework for constructing initial data sets for perturbations about spherically symmetric matter distributions. This framework facilitates setting initial data representing astrophysical sources of gravitational radiation involving relativistic stars. The procedure is based on Lichnerowicz-York's conformal approach to solve the constraints in Einstein's equations. The correspondence of these initial data sets in terms of the standard gauge perturbation variables in the Regge-Wheeler perturbation variables is established, and examples of initial data sets of merging neutron stars under the close-limit approximation are presented.

preprint1999arXiv

Quasi-Normal Modes of Stars and Black Holes

Perturbations of stars and black holes have been one of the main topics of relativistic astrophysics for the last few decades. They are of particular importance today, because of their relevance to gravitational wave astronomy. In this review we present the theory of quasi-normal modes of compact objects from both the mathematical and astrophysical points of view. The discussion includes perturbations of black holes (Schwarzschild, Reissner-Nordström, Kerr and Kerr-Newman) and relativistic stars (non-rotating and slowly-rotating). The properties of the various families of quasi-normal modes are described, and numerical techniques for calculating quasi-normal modes reviewed. The successes, as well as the limits, of perturbation theory are presented, and its role in the emerging era of numerical relativity and supercomputers is discussed.

preprint1997arXiv

Gravitational waves from pulsating stars: Evolving the perturbation equations for a relativistic star

We consider the perturbations of a relativistic star as an initial-value problem. Having discussed the formulation of the problem (the perturbation equations and the appropriate boundary conditions at the centre and the surface of the star) in detail we evolve the equations numerically from several different sets of initial data. In all the considered cases we find that the resulting gravitational waves carry the signature of several of the star's pulsation modes. Typically, the fluid $f$-mode, the first two $p$-modes and the slowest damped gravitational $w$-mode are present in the signal. This indicates that the pulsation modes may be an interesting source for detectable gravitational waves from colliding neutron stars or supernovae. We also survey the literature and find several indications of mode presence in numerical simulations of rotating core collapse and coalescing neutron stars. If such mode-signals can be detected by future gravitational-wave antennae one can hope to infer detailed information about neutron stars. Since a perturbation evolution should adequately describe the late time behaviour of a dynamically excited neutron star, the present work can also be used as a bench-mark test for future fully nonlinear simulations.

preprint1996arXiv

Spacetime modes of relativistic stars

The problem of relativistic stellar pulsations is studied in a somewhat ad hoc approximation that ignores all fluid motions. This {\em Inverse Cowling Approximation} (ICA) is motivated by two observations: 1) For highly damped ($w$-mode) oscillations the fluid plays very little role. 2) If the fluid motion is neglected the problem for polar oscillation modes becomes similar to that for axial modes. Using the ICA we find a polar mode spectrum that has all features of the $w$-mode spectrum of the full problem. Moreover, in the limit of superdense stars, we find the ICA spectrum to be qualitatively similar to that of the axial modes. These results clearly show the importance of general relativity for the pulsation modes of compact stars, and that there are modes whose existence do not depend on motions of the fluid at all; pure ``spacetime'' modes.

preprint1995arXiv

A NEW NUMERICAL APPROACH TO THE OSCILLATION MODES OF RELATIVISTIC STARS

The oscillation modes of a simple polytropic stellar model are studied. Using a new numerical approach (based on integration for complex coordinates) to the problem for the stellar exterior we have computed the eigenfrequencies of the highly damped w-modes. The results obtained agree well with recent ones of Leins, Nollert and Soffel (1993) Specifically, we are able to explain why several modes in this regime of the complex frequency plane could not be identified within the WKB approach of Kokkotas and Schutz (1992). Furthermore, we have established that the ``kink'' that was a prominent feature of the spectra of Kokkotas and Schutz, but did not appear in the results of Leins {\em et al.}, was a numerical artefact. Using our new numerical code we are also able to compute, for the first time, several of the slowly damped (p) modes for the considered stellar models. For very compact stars we find, somewhat surprisingly, that the damping of these modes does not decrease monotonically as one proceeds to higher oscillation frequencies. The existence of low-order modes that damp away much faster than anticipated may have implications for questions regarding stellar stability and the lifetime of gravitational-wave sources. The present results illustrate the accuracy and reliability of the complex-coordinate method and indicate that the method could prove to be of great use also in problems involving rotating stars. There is no apparent reason why the complex-coordinate approach should not extend to rotating stars, whereas it is accepted that all previous methods will fail to do so.

preprint1995arXiv

On the detectability of post-Newtonian effects in gravitational-wave emission of a coalescing binary

The effect of the recently obtained 2nd post-Newtonian corrections on the accuracy of estimation of parameters of the gravitational-wave signal from a coalescing binary is investigated. It is shown that addition of this correction degrades considerably the accuracy of determination of individual masses of the members of the binary. However the chirp mass and the time parameter in the sinal is still determined to a very good accuracy. The performance of the Newtonian filter is investigated and it is compared with performance of post-Newtonian search templates introduced recently. It is shown that both search templates can extract accurately useful information about the binary.

preprint1995arXiv

ON THE OSCILLATION SPECTRA OF ULTRA COMPACT STARS

Quasinormal modes of ultra compact stars with uniform energy density have been calculated. For less compact stars, there is only one very slowly damped polar mode (corresponding to the Kelvin f-mode) for each spherical harmonic index $l$. Further long-lived modes become possible for a sufficiently compact star (roughly when $M/R \ge 1/3$). We compare the characteristic frequencies of these resonant polar modes to the axial modes first found by Chandrasekhar and Ferrari [{\em Proc. Roy. Soc. London A} {\bf 434} 449 (1991)]. We find that the two spectra approach each other as the star is made more compact. The oscillation frequencies of the corresponding polar and axial modes agree to within a percent for stars more compact than $M/R = 0.42$. At the same time, the damping times are slightly different. The results illustrate that there is no real difference between the origin of these axial and polar modes: They are essentially spacetime modes.

preprint1995arXiv

TIDAL AND TIDAL-RESONANT EFFECTS IN COALESCING BINARIES

Tidal and tidal-resonant effects in coalescing compact binary systems are investigated by direct numerical integration of the equations of motion. For the stars polytropic models are used. The tidal effects are found to be dominated by the (non-resonant) $f$-modes. The effect of the $g$-mode-tidal resonances is obtained. The tidal interaction is shown to be of interest especially for low-mass binaries. There exists a characteristic final plunge orbit beyond which the system cannot remain stable even if radiation reaction is not taken into account; in agreement with results obtained by Lai et al. \shortcite{Lai93}. The importance of the investigated effects for the observation of gravitational waves on Earth is discussed.