Researcher profile

Nikolaos Stergioulas

Nikolaos Stergioulas contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
39works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

39 published item(s)

preprint2026arXiv

ESO Expanding Horizon White Paper: Revealing the properties of matter at supranuclear densities with gravitational waves

Understanding dense matter under extreme conditions is one of the most fundamental puzzles in modern physics. Complex interactions give rise to emergent, collective phenomena. While nuclear experiments and Earth - based colliders provide valuable insights, much of the quantum chromodynamics phase diagram at high density and low temperature remains accessible only through astrophysical observations of neutron stars, neutron star mergers, and stellar collapse. Astronomical observations thus offer a direct window to the physics on subatomic scales with gravitational waves presenting an especially clean channel. Next-generation gravitational - wave observatories, such as the Einstein Telescope, would serve as unparalleled instruments to transform our understanding of neutron star matter. They will enable the detection of up to tens of thousands of binary neutron star and neutron star - black hole mergers per year, a dramatic increase over the few events accessible with current detectors. They will provide an unprecedented precision in probing cold, dense matter during the binary inspiral, exceeding by at least an order of magnitude what current facilities can achieve. Moreover, these observatories will allow us to explore uncharted regimes of dense matter at finite temperatures produced in a subset of neutron star mergers, areas that remain entirely inaccessible to current instruments. Together with multimessenger observations, these measurements will significantly deepen our knowledge of dense nuclear matter.

preprint2025arXiv

Robust and fast parameter estimation for gravitational waves from binary neutron star merger remnants

We present a robust and efficient methodology for parameter estimation of gravitational waves generated during the post-merger phase of binary neutron star mergers. Our approach leverages an analytic waveform model combined with empirical relations to predict prior ranges for the post-merger frequencies based on measurements of the chirp mass and effective tidal deformability in the inspiral phase. This enables robust inference of the main features of the post-merger spectrum, avoiding possible multi-modality induced by wide priors. Using waveforms derived from numerical relativity, we systematically validate our model across a broad spectrum of neutron star equations of state and mass configurations, demonstrating high fitting factors. Our method can be applied in future detections of gravitational waves from the post-merger phase with third-generation gravitational wave observatories. Furthermore, by integrating the Preconditioned Monte Carlo sampling method within the pocoMC framework, we achieve substantial computational acceleration compared to conventional Bayesian techniques.

preprint2022arXiv

Analytic models of the spectral properties of gravitational waves from neutron star merger remnants

We present a new analytic model describing gravitational wave emission in the post-merger phase of binary neutron star mergers. The model is described by a number of physical parameters that are related to various oscillation modes, quasi-linear combination tones or non-linear features that appear in the post-merger phase. The time evolution of the main post-merger frequency peak is taken into account and it is described by a two-segment linear expression. The effectiveness of the model, in terms of the fitting factor or, equivalently, the reduction in the detection rate, is evaluated along a sequence of equal-mass simulations of varying mass. We find that all parameters of the analytic model correlate with the total binary mass of the system. For high masses, we identify new spectral features originating from the non-linear coupling between the quasi-radial oscillation and the antipodal tidal deformation, the inclusion of which significantly improves the fitting factors achieved by the model. We can thus model the post-merger gravitational-wave emission with an analytic model that achieves high fitting factors for a wide range of total binary masses. Our model can be used for the detection and parameter estimation of the post-merger phase in upcoming searches with upgraded second-generation detectors, such as aLIGO+ and aVirgo+, with future, third-generation detectors (Einstein Telescope and Cosmic Explorer) or with dedicated, high-frequency detectors.

preprint2022arXiv

Astrophysical constraints on compact objects in 4D Einstein-Gauss-Bonnet gravity

We study the properties of compact objects in a particular 4D Horndeski theory originating from higher dimensional Einstein-Gauss-Bonnet gravity. Remarkably, an exact vacuum solution is known. This compact object differs from general relativity mostly in the strong field regime. We discuss some properties of black holes in this framework and investigate in detail the properties of neutron stars, both static and in slow rotation. We find that for relatively modest deviations from general relativity, the secondary object in GW190814 is compatible with being a slowly-rotating neutron star, without resorting to very stiff or exotic equations of state. Remarkably, the equilibrium sequence of neutron stars matches asymptotically to the black hole limit, completetly closing the mass gap between neutron stars and black holes of same radius, although the stability of equilibrium solutions has yet to be determined. As a consequence, there exists a universal endpoint for the neutron star sequence, independent of the equation of state. In light of our results and of current observational constraints, we discuss specific constraints on the coupling constant that parametrizes deviations from general relativity in this theory.

preprint2022arXiv

Astrophysical Implications of Neutron Star Inspiral and Coalescence

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.

preprint2022arXiv

Models of binary neutron star remnants with tabulated equations of state

The emergence of novel differential rotation laws that can reproduce the rotational profile of binary neutron star merger remnants has opened the way for the construction of equilibrium models with properties that resemble those of remnants in numerical simulations. We construct models of merger remnants, using a recently introduced 4-parameter differential rotation law and three tabulated, zero-temperature equations of state. The models have angular momenta that are determined by empirical relations, constructed through numerical simulations. After a systematic exploration of the parameter space of merger remnant equilibrium sequences, which includes the determination of turning points along constant angular momentum sequences, we find that a particular rotation law can reproduce the threshold mass to prompt collapse to a black hole with a relative difference of only $\sim 1\%$ with respect to numerical simulations, in all cases considered. Furthermore, our results indicate a possible correlation between the compactness of equilibrium models of remnants at the threshold mass and the compactness of maximum-mass nonrotating models. Another key prediction of binary neutron star merger simulations is a relatively slowly rotating inner region, where the angular velocity $Ω$ (as measured by an observer at infinity) is mostly due to the frame dragging angular velocity $ω$. In our investigation of the parameter space of the adopted differential rotation law, we naturally find quasi-spherical (Type A) remnant models with this property. Our investigation clarifies the impact of the differential rotation law and of the equation of state on key properties of binary neutron star remnants and lays the groundwork for including thermal effects in future studies.

preprint2022arXiv

New Horizons for Fundamental Physics with LISA

The Laser Interferometer Space Antenna (LISA) has the potential to reveal wonders about the fundamental theory of nature at play in the extreme gravity regime, where the gravitational interaction is both strong and dynamical. In this white paper, the Fundamental Physics Working Group of the LISA Consortium summarizes the current topics in fundamental physics where LISA observations of GWs can be expected to provide key input. We provide the briefest of reviews to then delineate avenues for future research directions and to discuss connections between this working group, other working groups and the consortium work package teams. These connections must be developed for LISA to live up to its science potential in these areas.

preprint2022arXiv

Nonlinear dynamics of oscillating neutron stars in scalar-tensor gravity

The spectrum of oscillating compact objects can be considerably altered in alternative theories of gravity. In particular, it may be enriched by modes with no counterpart in general relativity, tied to the dynamics of additional degrees of freedom generically present in these theories. Detection of these modes, e.g. in the gravitational-wave signal from a binary compact object coalescence, could provide a powerful tool to probe the underlying theory of gravity. To access the potential of such a detection, it is crucial to understand the linear and nonlinear spectral features of dynamically formed, oscillating compact objects in alternative theories of gravity. As a step towards that goal, in this work we present a suite of 1+1 numerical relativity simulations of neutron stars in scalar-tensor theories, we carefully analyze the spectrum of stellar pulsations, and we compare results with expectations from linear perturbation theory. This allows us to build intuition for the case of binary neutron star mergers. Additionally, the models investigated in this work are representatives of two broad classes, in which the scalar field couples either strongly or weakly with the fluid. The distinct phenomenology of the nonlinear dynamics that we identify for each class of models, may find counterparts also in other alternative theories of gravity.

preprint2020arXiv

Numerical-relativity simulations of long-lived remnants of binary neutron star mergers

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.

preprint2020arXiv

Spectral classification of gravitational-wave emission and equation of state constraints in binary neutron star mergers

The first detection of gravitational waves from a binary neutron star merger (GW170817) and the accompanying electromagnetic emission has impressively advanced our understanding of the merger process and has set some first constraints on the macroscopic properties of neutron stars, with direct implications for the high-density equation of state. We discuss work on neutron star mergers focusing on the postmerger gravitational-wave emission. These studies are based on numerical simulations of the merger and survey a large sample of candidate equations of state for neutron star matter. The goal is to connect observables with the underlying physics questions. This offers a way to constrain the properties of high-density matter through the determination of neutron star radii, as inferred by an empirical relation connecting the dominant gravitational wave frequency peak in the postmerger phase to the radius of nonrotating neutron stars of a certain mass. We clarify the physical origin of secondary peaks and discuss a spectral classification scheme, based on their relative strength. Observational prospects for the dominant and the secondary peaks are also discussed. The threshold mass to black hole collapse is connected by another empirical relation to the maximum mass and compactness of nonrotating neutron stars, which can be derived semi-analytically. The observation of GW170817 then sets an absolute minimum radius for neutron stars of typical masses, based only on a minimal number of assumptions. We discuss future prospects, in light of the planned upgrades of the current gravitational wave detectors.

preprint2019arXiv

Empirical relations for gravitational-wave asteroseismology of binary neutron star mergers

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.

preprint2016arXiv

Coherent magneto-elastic oscillations in superfluid magnetars

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énic 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.

preprint2015arXiv

Exploring properties of high-density matter through remnants of neutron-star mergers

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)

preprint2015arXiv

Observing Gravitational Waves From The Post-Merger Phase Of Binary Neutron Star Coalescence

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

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

Unusual behavior of epicyclic frequencies around rapidly rotating compact stars

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.

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

Epicyclic frequencies for rotating strange quark stars: the importance of stellar oblateness

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.

preprint2014arXiv

Modulating the magnetosphere of magnetars by internal magneto-elastic oscillations

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én 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.

preprint2014arXiv

On the accuracy of the IWM-CFC approximation in differentially rotating relativistic stars

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.

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.

preprint2013arXiv

Imprints of superfluidity on magneto-elastic QPOs of SGRs

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.

preprint2013arXiv

Magneto-elastic oscillations of neutron stars: exploring different magnetic field configurations

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én 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én 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én 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).

preprint2012arXiv

Magneto-elastic oscillations of neutron stars with dipolar magnetic fields

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én 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én QPOs. When the Alfvén 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én-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.

preprint2011arXiv

An improved method for constructing models of self-gravitating tori around black holes

General-relativistic models of self-gravitating tori around black holes are constructed with a self-consistent-field method in compactified coordinates. The numerical code is highly accurate and robust, allowing for the construction of models that exactly fill their Roche lobe, when a cusp exists. As a first application, we focus on self-consistent models with cusp, having different values of constant specific angular momentum. Scaling all results with the mass of the black hole, we find evidence that models with constant specific angular momentum that can fill their Roche lobe are still limited by $l<4M_{\rm BH}$ (as is the case for models constructed in a fixed background metric) even for heavy tori.

preprint2011arXiv

Collapse of differentially rotating neutron stars and cosmic censorship

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., &#34;sub-Kerr&#34; 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., &#34;supra-Kerr&#34; 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.

preprint2011arXiv

Gravitational waves and nonaxisymmetric oscillation modes in mergers of compact object binaries

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.

preprint2010arXiv

Magneto-elastic torsional oscillations of magnetars

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én 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én 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.

preprint2010arXiv

On the frequency band of the f-mode CFS instability

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.

preprint2010arXiv

Stability of general-relativistic accretion disks

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 $Γ$-law equation of state with $Γ=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.

preprint2010arXiv

Stationary, Axisymmetric Neutron Stars with Meridional Circulation in General Relativity

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.

preprint2009arXiv

Alfven QPOs in magnetars in the anelastic approximation

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.

preprint2005arXiv

Black hole formation through fragmentation of toroidal polytropes

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.

preprint2003arXiv

Rotating Stars in Relativity

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.

preprint2000arXiv

Nonlinear r-modes in Rapidly Rotating Relativistic Stars

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

preprint1998arXiv

Construction of Highly Accurate Models of Rotating Neutron Stars - Comparison of Three Different Numerical Schemes

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&#39;&#39; and we display them in detail for all models.

preprint1998arXiv

Quasi-normal modes of rotating relativistic stars - neutral modes for realistic equations of state

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&#39;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.

preprint1994arXiv

Comparing Models of Rapidly Rotating Relativistic Stars Constructed by Two Numerical Methods

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.