Normalized to: Rezzolla, L.

**187** article(s) in total. **1169** co-authors, from **1** to **18** common article(s). Median position in authors list is **3,0**.

[1] oai:arXiv.org:2006.14601 [pdf] - 2122171

A lower bound on the maximum mass if the secondary in GW190814 was once
a rapidly spinning neutron star

Comments: 6 pages, 3 figures, Comments welcome

Submitted: **2020-06-25**

The recent detection of GW190814 featured the merger of a binary with a
primary having a mass of $\sim 23\,M_{\odot}$ and a secondary with a mass of
$\sim 2.6\,M_{\odot}$. While the primary was most likely a black hole, the
secondary could be interpreted as either the lightest black hole or the most
massive neutron star ever observed, but also as the indication of a novel class
of exotic compact objects. We here argue that the secondary in GW190814 needs
not be an ab-initio black hole nor an exotic object; rather, based on our
current understanding of the nuclear-matter equation of state, it can be a
rapidly rotating neutron star that collapsed to a rotating black hole at some
point before merger. Using universal relations connecting the masses and spins
of uniformly rotating neutron stars, we estimate the spin, $0.49 \lesssim \chi
\lesssim 0.68$, of the secondary -- a quantity not constrained so far by the
detection -- and a novel strict lower bound on the maximum mass, $M_{\rm TOV} >
2.08^{+0.04}_{-0.04}\, \,M_{\odot}$, of nonrotating neutron stars, consistent
with recent observations of a very massive pulsar. The new lower bound also
remains valid even in the less likely scenario in which the secondary neutron
star never collapsed to a black hole.

[2] oai:arXiv.org:2003.13580 [pdf] - 2097380

Two-moment scheme for general-relativistic radiation hydrodynamics: a
systematic description and new applications

Comments: 22 pages, 15 figures, matches version accepted to MNRAS

Submitted: **2020-03-30**, last modified: **2020-05-19**

We provide a systematic description of the steps necessary -- and of the
potential pitfalls to be encountered -- when implementing a two-moment scheme
within an Implicit-Explicit (IMEX) scheme to include radiative-transfer
contributions in numerical simulations of general-relativistic
(magneto-)hydrodynamics. We make use of the M1 closure, which provides an exact
solution for the optically thin and thick limit, and an interpolation between
these limits. Special attention is paid to the efficient solution of the
emerging set of implicit conservation equations. In particular, we present an
efficient method for solving these equations via the inversion of a $4\times
4$-matrix within an IMEX scheme. While this method relies on a few
approximations, it offers a very good compromise between accuracy and
computational efficiency. After a large number of tests in special relativity,
we couple our new radiation code, \texttt{FRAC}, with the general-relativistic
magnetohydrodynamics code \texttt{BHAC} to investigate the radiative Michel
solution, namely, the problem of spherical accretion onto a black hole in the
presence of a radiative field. By performing the most extensive exploration of
the parameter space for this problem, we find that the accretion's efficiency
can be expressed in terms of physical quantities such as temperature, $T$,
luminosity, $L$, and black-hole mass, $M$, via the expression
$\varepsilon=(L/L_{\rm Edd})/(\dot{M}/\dot{M}_{\rm Edd})= 7.41\times
10^{-7}\left(T/10^6\,\mathrm{K}\right)^{0.22} \left(L/L_\odot\right)^{0.48}
\left(M/M_\odot\right)^{0.48}$, where $L_{\mathrm{Edd}}$ and
$\dot{M}_{\mathrm{Edd}}$ are the Eddington luminosity and accretion rate,
respectively. Finally, we also consider the accretion problem away from
spherical symmetry, finding that the solution is stable under perturbations in
the radiation field.

[3] oai:arXiv.org:2003.10391 [pdf] - 2095621

The heavier the better: how to constrain mass ratios and spins of
high-mass neutron-star mergers

Comments: 6 pages, 4 figures, Accepted version

Submitted: **2020-03-23**, last modified: **2020-05-11**

The first binary neutron-star merger event, GW170817, and its bright
electromagnetic counterpart have provided a remarkable amount of information.
By contrast, the second event, GW190425, with $M_{\rm
tot}=3.4^{+0.3}_{-0.1}\,M_{\odot}$ and the lack of an electromagnetic
counterpart, has hardly improved our understanding of neutron-star physics.
While GW190425 is compatible with a scenario in which the merger has lead to a
prompt collapse to a black hole and little ejected matter to power a
counterpart, determining the mass ratio and the effective spin $\tilde{\chi}$
of the binary remains difficult. This is because gravitational waveforms cannot
yet well constrain the component spins of the binary. However, since the mass
of GW190425 is significantly larger than the maximum mass for nonrotating
neutron stars, $M_{_{\rm TOV}}$, the mass ratio $q$ cannot be too small, as the
heavier star would not be gravitationally stable. Making use of universal
relations and a large number of equations of state, we provide limits in the
$(\tilde{\chi},q)$ plane for GW190425, namely: $q_{\rm min} \geq 0.38$ and $
\tilde{\chi}_{\rm max}\leq 0.20$, assuming $M_\mathrm{tot} \simeq 3.4\,
M_\odot$. Finally, we show how future observations of high-mass binaries can
provide a lower bound on $M_{_{\rm TOV}}$.

[4] oai:arXiv.org:2005.02677 [pdf] - 2090729

Tidal deformability and gravitational-wave phase evolution of magnetised
compact-star binaries

Comments: 21 pages, 6 figures

Submitted: **2020-05-06**

The evolution of the gravitational-wave phase in the signal produced by
inspiralling binaries of compact stars is modified by the nonzero deformability
of the two stars. Hence, the measurement of these corrections has the potential
of providing important information on the equation of state of nuclear matter.
Extensive work has been carried out over the last decade to quantify these
corrections, but it has so far been restricted to stars with zero intrinsic
magnetic fields. While the corrections introduced by the magnetic tension and
magnetic pressure are expected to be subdominant, it is nevertheless useful to
determine the precise conditions under which these corrections become
important. To address this question, we have carried out a second-order
perturbative analysis of the tidal deformability of magnetised compact stars
under a variety of magnetic-field strengths and equations of state describing
either neutron stars or quark stars. Overall, we find that magnetically induced
corrections to the tidal deformability will produce changes in the
gravitational-wave phase evolution that are unlikely to be detected for
realistic magnetic field i.e., $B\sim 10^{10} - 10^{12}\,{\rm G}$. At the same
time, if the magnetic field is unrealistically large, i.e., $B\sim
10^{16}\,{\rm G}$, these corrections would produce a sizeable contribution to
the phase evolution, especially for quark stars. In the latter case, the
induced phase differences would represent a unique tool to measure the
properties of the magnetic fields, providing information that is otherwise hard
to quantify.

[5] oai:arXiv.org:2005.01339 [pdf] - 2088338

On the prospects of imaging Sagittarius A* from space

Comments: 4 pages, 4 figures, published in Proceedings IAU Symposium No. 342,
2018

Submitted: **2020-05-04**

Very Long Baseline Interferometry (VLBI) at sub-millimeter waves has the
potential to image the shadow of the black hole in the Galactic Center,
Sagittarius A* (Sgr A*), and thereby test basic predictions of the theory of
general relativity. We investigate the imaging prospects of a new Space VLBI
mission concept. The setup consists of two satellites in polar or equatorial
circular Medium-Earth Orbits with slightly different radii, resulting in a
dense spiral-shaped uv-coverage with long baselines, allowing for extremely
high-resolution and high-fidelity imaging of radio sources. We simulate
observations of a general relativistic magnetohydrodynamics model of Sgr A* for
this configuration with noise calculated from model system parameters. After
gridding the $uv$-plane and averaging visibilities accumulated over multiple
months of integration, images of Sgr A* with a resolution of up to 4 $\mu$as
could be reconstructed, allowing for stronger tests of general relativity and
accretion models than with ground-based VLBI.

[6] oai:arXiv.org:1912.09340 [pdf] - 2086911

Postmerger Gravitational-Wave Signatures of Phase Transitions in Binary
Mergers

Comments: 6 pages, 4 figures, includes supplemental material, matches version
accepted to PRL

Submitted: **2019-12-19**, last modified: **2020-04-30**

With the first detection of gravitational waves from a binary system of
neutron stars, GW170817, a new window was opened to study the properties of
matter at and above nuclear-saturation density. Reaching densities a few times
that of nuclear matter and temperatures up to $100\,\rm{MeV}$, such mergers
also represent potential sites for a phase transition (PT) from confined
hadronic matter to deconfined quark matter. While the lack of a postmerger
signal in GW170817 has prevented us from assessing experimentally this
scenario, two theoretical studies have explored the postmerger
gravitational-wave signatures of PTs in mergers of binary systems of neutron
stars. We here extend and complete the picture by presenting a novel signature
of the occurrence of a PT. More specifically, using fully general-relativistic
hydrodynamic simulations and employing a suitably constructed equation of state
that includes a PT, we present the occurrence of a "delayed PT", i.e. a PT that
develops only some time after the merger and produces a metastable object with
a quark-matter core, i.e. a hypermassive hybrid star. Because in this scenario,
the postmerger signal exhibits two distinct fundamental gravitational-wave
frequencies -- before and after the PT -- the associated signature promises to
be the strongest and cleanest among those considered so far, and one of the
best signatures of the production of quark matter in the present Universe.

[7] oai:arXiv.org:2004.13782 [pdf] - 2101497

Common-envelope Dynamics of a Stellar-mass Black Hole: General
Relativistic Simulations

Comments: 16 pages, 27 figures. Accepted for publication in ApJ

Submitted: **2020-04-28**

With the goal of providing more accurate and realistic estimates of the
secular behavior of the mass accretion and drag rates in the "common-envelope"
scenario encountered when a black hole or a neutron star moves in the stellar
envelope of a red supergiant star, we have carried out the first general
relativistic simulations of the accretion flow onto a nonrotating black hole
moving supersonically in a medium with regular but different density gradients.
The simulations reveal that the supersonic motion always rapidly reaches a
stationary state and it produces a shock cone in the downstream part of the
flow. In the absence of density gradients we recover the phenomenology already
observed in the well-known Bondi--Hoyle--Lyttleton accretion problem, with
super-Eddington mass accretion rate and a shock cone whose axis is stably
aligned with the direction of motion. However, as the density gradient is made
stronger, the accretion rate also increases and the shock cone is progressively
and stably dragged toward the direction of motion. With sufficiently large
gradients, the shock-cone axis can become orthogonal to the direction, or even
move in the upstream region of the flow in the case of the largest density
gradient. Together with the phenomenological aspects of the accretion flow, we
have also quantified the rates of accretion of mass and momentum onto the black
hole. Simple analytic expressions have been found for the rates of accretion of
mass, momentum, drag force, and bremsstrahlung luminosity, all of which have
been employed in the astrophysical modelling of the secular evolution of a
binary system experiencing a common-envelope evolution...

[8] oai:arXiv.org:2004.01161 [pdf] - 2074562

SYMBA: An end-to-end VLBI synthetic data generation pipeline

Comments: 20 pages, 15 figures, accepted for publication in A&A

Submitted: **2020-04-02**

Realistic synthetic observations of theoretical source models are essential
for our understanding of real observational data. In using synthetic data, one
can verify the extent to which source parameters can be recovered and evaluate
how various data corruption effects can be calibrated. These studies are
important when proposing observations of new sources, in the characterization
of the capabilities of new or upgraded instruments, and when verifying
model-based theoretical predictions in a comparison with observational data. We
present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a
novel synthetic data generation pipeline for Very Long Baseline Interferometry
(VLBI) observations. SYMBA takes into account several realistic atmospheric,
instrumental, and calibration effects. We used SYMBA to create synthetic
observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has
recently captured the first image of a black hole shadow. After testing SYMBA
with simple source and corruption models, we study the importance of including
all corruption and calibration effects. Based on two example general
relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed
case studies to assess the attainable image quality with the current and future
EHT array for different weather conditions. The results show that the effects
of atmospheric and instrumental corruptions on the measured visibilities are
significant. Despite these effects, we demonstrate how the overall structure of
the input models can be recovered robustly after performing calibration steps.
With the planned addition of new stations to the EHT array, images could be
reconstructed with higher angular resolution and dynamic range. In our case
study, these improvements allowed for a distinction between a thermal and a
non-thermal GRMHD model based on salient features in reconstructed images.

[9] oai:arXiv.org:1910.07974 [pdf] - 2068944

Detection of pulses from the Vela pulsar at millimeter wavelengths with
phased ALMA

Comments: 8 pages, 6 figures, published in ApJ letter

Submitted: **2019-10-17**, last modified: **2020-03-24**

We report on the first detection of pulsed radio emission from a radio pulsar
with the ALMA telescope. The detection was made in the Band-3 frequency range
(85-101 GHz) using ALMA in the phased-array mode developed for VLBI
observations. A software pipeline has been implemented to enable a regular
pulsar observing mode in the future. We describe the pipeline and demonstrate
the capability of ALMA to perform pulsar timing and searching. We also measure
the flux density and polarization properties of the Vela pulsar (PSR
J0835$-$4510) at mm-wavelengths, providing the first polarimetric study of any
ordinary pulsar at frequencies above 32 GHz. Finally, we discuss the lessons
learned from the Vela observations for future pulsar studies with ALMA,
particularly for searches near the supermassive black hole in the Galactic
Center, and the potential of using pulsars for polarization calibration of
ALMA.

[10] oai:arXiv.org:1907.10328 [pdf] - 2063742

Beyond second-order convergence in simulations of magnetised binary
neutron stars with realistic microphysics

Comments: 13 pages, 9 figures matches accepted version

Submitted: **2019-07-24**, last modified: **2020-03-13**

We investigate the impact of using high-order numerical methods to study the
merger of magnetised neutron stars with finite-temperature microphysics and
neutrino cooling in full general relativity. By implementing a fourth-order
accurate conservative finite-difference scheme we model the inspiral together
with the early post-merger and highlight the differences to traditional
second-order approaches at the various stages of the simulation. We find that
even for finite-temperature equations of state, convergence orders higher than
second order can be achieved in the inspiral and post-merger for the
gravitational-wave phase. We further demonstrate that the second-order scheme
overestimates the amount of proton-rich shock-heated ejecta, which can have an
impact on the modelling of the dynamical part of the kilonova emission.
Finally, we show that already at low resolution the growth rate of the magnetic
energy is consistently resolved by using a fourth-order scheme.

[11] oai:arXiv.org:2003.03597 [pdf] - 2105646

On the opening angle of magnetised jets from neutron-star mergers: the
case of GRB170817A

Comments: 9 pages, 5 figures, submitted

Submitted: **2020-03-07**

The observations of GW170817/GRB170817A have confirmed that the coalescence
of a neutron-star binary is the progenitor of a short gamma-ray burst. In the
standard picture of a short gamma-ray burst, a collimated highly relativistic
outflow is launched after merger and it successfully breaks out from the
surrounding ejected matter. Using initial conditions inspired from
numerical-relativity binary neutron-star merger simulations, we have performed
general-relativistic hydrodynamic (HD) and magnetohydrodynamic (MHD)
simulations in which the jet is launched and propagates self-consistently. The
complete set of simulations suggests that: (i) MHD jets have an intrinsic
energy and velocity polar structure with a ``hollow core'' subtending an angle
$\theta_{\rm core}\approx4^{\circ}-5^{\circ}$ and an opening angle of
$\theta_{\rm jet}\gtrsim10^{\circ}$; (ii) MHD jets eject significant amounts of
matter and two orders of magnitude more than HD jets; (iii) the energy
stratification in MHD jets naturally yields the power-law energy scaling
$E(>\Gamma\beta)\propto(\Gamma\beta)^{-4.5}$; (iv) MHD jets provide fits to the
afterglow data from GRB170817A that are comparatively better than those of the
HD jets and without free parameters; (v) finally, both of the best-fit HD/MHD
models suggest an observation angle $\theta_{\rm obs} \simeq 21^{\circ}$ for
GRB170817A.

[12] oai:arXiv.org:2002.01777 [pdf] - 2105617

Plasmoid formation in global GRMHD simulations and AGN flares

Comments: 17 pages, MNRAS submitted

Submitted: **2020-02-05**

One of the main dissipation processes acting on all scales in relativistic
jets is thought to be governed by magnetic reconnection. Such dissipation
processes have been studied in idealized environments, such as reconnection
layers, which evolve in merging islands and lead to the production of
plasmoids, ultimately resulting in efficient particle acceleration. In
accretion flows onto black holes, reconnection layers can be developed and
destroyed rapidly during the turbulent evolution of the flow. We present a
series of two-dimensional general-relativistic magnetohydrodynamic simulations
of tori accreting onto rotating black holes focusing our attention on the
formation and evolution of current sheets. Initially, the tori are endowed with
a poloidal magnetic field having a multi-loop structure along the radial
direction and with an alternating polarity. During reconnection processes,
plasmoids and plasmoid chains are developed leading to a flaring activity and
hence to a variable electromagnetic luminosity. We describe the methods
developed to track automatically the plasmoids that are generated and ejected
during the simulation, contrasting the behaviour of multi-loop initial data
with that encountered in typical simulations of accreting black holes having
initial dipolar field composed of one loop only. Finally, we discuss the
implications that our results have on the variability to be expected in
accreting supermassive black holes.

[13] oai:arXiv.org:1909.01411 [pdf] - 1996114

Studying Black Holes on Horizon Scales with VLBI Ground Arrays

Comments: Astro2020 APC White Paper, 10 pages, 7 figures

Submitted: **2019-09-03**, last modified: **2019-11-08**

High-resolution imaging of supermassive black holes is now possible, with new
applications to testing general relativity and horizon-scale accretion and
relativistic jet formation processes. Over the coming decade, the EHT will
propose to add new strategically placed VLBI elements operating at 1.3mm and
0.87mm wavelength. In parallel, development of next-generation backend
instrumentation, coupled with high throughput correlation architectures, will
boost sensitivity, allowing the new stations to be of modest collecting area
while still improving imaging fidelity and angular resolution. The goal of
these efforts is to move from imaging static horizon scale structure to dynamic
reconstructions that capture the processes of accretion and jet launching in
near real time.

[14] oai:arXiv.org:1910.13893 [pdf] - 2063759

On the Deconfinement Phase Transition in Neutron-Star Mergers

Comments: 10 pages, 6 figures

Submitted: **2019-10-30**

We study in detail the nuclear aspects of a neutron-star merger in which
deconfinement to quark matter takes place. For this purpose, we make use of the
Chiral Mean Field (CMF) model, an effective relativistic model that includes
self-consistent chiral symmetry restoration and deconfinement to quark matter
and, for this reason, predicts the existence of different degrees of freedom
depending on the local density/chemical potential and temperature. We then use
the out-of-chemical-equilibrium finite-temperature CMF equation of state in
full general-relativistic simulations to analyze which regions of different QCD
phase diagrams are probed and which conditions, such as strangeness and
entropy, are generated when a strong first-order phase transition appears. We
also investigate the amount of electrons present in different stages of the
merger and discuss how far from chemical equilibrium they can be and, finally,
draw some comparisons with matter created in supernova explosions and heavy-ion
collisions.

[15] oai:arXiv.org:1910.10193 [pdf] - 1985370

First M87 Event Horizon Telescope Results and the Role of ALMA

Comments: 11 pages + cover page, 6 figures

Submitted: **2019-10-22**

In April 2019, the Event Horizon Telescope (EHT) collaboration revealed the
first image of the candidate super-massive black hole (SMBH) at the centre of
the giant elliptical galaxy Messier 87 (M87). This event-horizon-scale image
shows a ring of glowing plasma with a dark patch at the centre, which is
interpreted as the shadow of the black hole. This breakthrough result, which
represents a powerful confirmation of Einstein's theory of gravity, or general
relativity, was made possible by assembling a global network of radio
telescopes operating at millimetre wavelengths that for the first time included
the Atacama Large Millimeter/ submillimeter Array (ALMA). The addition of ALMA
as an anchor station has enabled a giant leap forward by increasing the
sensitivity limits of the EHT by an order of magnitude, effectively turning it
into an imaging array. The published image demonstrates that it is now possible
to directly study the event horizon shadows of SMBHs via electromagnetic
radiation, thereby transforming this elusive frontier from a mathematical
concept into an astrophysical reality. The expansion of the array over the next
few years will include new stations on different continents - and eventually
satellites in space. This will provide progressively sharper and
higher-fidelity images of SMBH candidates, and potentially even movies of the
hot plasma orbiting around SMBHs. These improvements will shed light on the
processes of black hole accretion and jet formation on event-horizon scales,
thereby enabling more precise tests of general relativity in the truly strong
field regime.

[16] oai:arXiv.org:1904.04220 [pdf] - 1978585

Impact of high spins on the ejection of mass in GW170817

Comments: matches accepted version

Submitted: **2019-04-08**, last modified: **2019-10-14**

Following the detection of GW170817 and the accompanying kilonova AT2017gfo,
it has become crucial to model and understand the various channels through
which mass is ejected in neutron-star binary mergers. We discuss the impact
that high stellar spins prior to merger have on the ejection of mass focussing,
in particular, on the dynamically ejected mass by performing
general-relativistic magnetohydrodynamic simulations employing
finite-temperature equations of state and neutrino-cooling effects. Using eight
different models with dimensionless spins ranging from $\chi\simeq-0.14$ to
$\chi\simeq0.29$ we discuss how the presence of different spins affects the
angular distribution and composition of the ejected matter. Most importantly,
we find that the dynamical component of the ejected mass can be strongly
suppressed in the case of high spins aligned with the orbital angular momentum.
In this case, in fact, the merger remnant has an excess angular momentum
yielding a more extended and "colder" object, with reduced ability to shed mass
dynamically. We discuss how this result impacts the analysis of the recent
merger event GW170817 and its kilonova afterglow.

[17] oai:arXiv.org:1812.04526 [pdf] - 1977341

Magnetically inspired explosive outflows from neutron-star mergers

Comments: matched published version ApJ Letters

Submitted: **2018-12-11**, last modified: **2019-10-10**

Binary neutron-star mergers have long been associated with short-duration
gamma-ray bursts (GRBs). This connection was confirmed with the first
coincident detection of gravitational waves together with electromagnetic
radiation from GW170817. The basic paradigm for short-duration GRBs includes an
ultra-relativistic jet, but the low-luminosity prompt emission together with
follow-up radio and X-ray observations have hinted that this picture may be
different in the case of GW170817. In particular, it has been proposed that
large amounts of the magnetic energy that is amplified after the merger, can be
released when the remnant collapses to a black hole, giving rise to a
quasi-spherical explosion impacting on the merger ejecta. Through numerical
simulations we investigate this scenario for a range of viewing angles,
injected energies and matter densities at the time of the collapse. Depending
on the magnitude of the energy injection and the remnant density, we find two
types of outflows: one with a narrow relativistic core and one with a
wide-angle, but mildly relativistic outflow. Furthermore, very wide outflows
are possible, but require energy releases in excess of 10^52 erg.

[18] oai:arXiv.org:1908.10767 [pdf] - 1950665

TeraHertz Exploration and Zooming-in for Astrophysics (THEZA): ESA
Voyage 2050 White Paper

Comments: White Paper submitted in response to the ESA Call Voyage 2050; 25
pages, 12 figures

Submitted: **2019-08-28**

The astrophysical agenda of the 21st century requires a very sharp view of
celestial objects. High angular resolution studies are essential for
fundamental studies of a broad variety of astrophysical phenomena ranging from
relativistic physics of black holes, their gravitational and electromagnetic
imprints, violent transient processes, including those producing detectable
gravitational waves, birth and evolution of planetary systems. Over the past
decades, radio astronomy made huge leap in achieving ground-breaking angular
resolution measured in tens of microarcseconds (one tenth of nanoradian and
better). Recently a global Event Horizon Telescope (EHT) collaboration obtained
first direct images of the shadow of a super-massive black hole in the nucleus
of the active galaxy M87. These observations were conducted at 230 GHz. The two
first generation Space Very Long Baseline Interferometry (VLBI) missions,
VSOP/HALCA led by the Japan Aerospace Exploration Agency (JAXA) and RadioAstron
led by the Russia Roscosmos State Corporation and Russia Academy of Sciences,
achieved the highest angular resolution at frequencies from 0.3 to 22 GHz in
observations conducted in the period 1997 - 2019. The next step in advancing
high angular resolution radio astronomy is in combining high frequency
(millimeter and sub-millimeter wavelengths) and interferometric baselines
exceeding the Earth diameter. The present THEZA White Paper describes a
combination which would unify technology developments in giga-/tera-hertz
instrumentation and space-borne radio astronomy. The current preprint version
of the THEZA White Paper is slightly re-formatted and edited comparing to the
official submitted version.

[19] oai:arXiv.org:1905.04900 [pdf] - 1943278

Optimal neutron-star mass ranges to constrain the equation of state of
nuclear matter with electromagnetic and gravitational-wave observations

Comments: 10 pages, 6 figures; matches accepted version; data available via
https://doi.org/10.5281/zenodo.3260991 animated Fig. 1 available via
https://doi.org/10.5281/zenodo.3363615

Submitted: **2019-05-13**, last modified: **2019-08-08**

Exploiting a very large library of physically plausible equations of state
(EOSs) containing more than $10^{7}$ members and yielding more than $10^{9}$
stellar models, we conduct a survey of the impact that a neutron-star radius
measurement via electromagnetic observations can have on the EOS of nuclear
matter. Such measurements are soon to be expected from the ongoing
\textit{Neutron Star Interior Composition Explorer} mission and will complement
the constraints on the EOS from gravitational-wave detections. Thanks to the
large statistical range of our EOS library, we can obtain a first quantitative
estimate of the commonly made assumption that the high-density part of the EOS
is best constrained when measuring the radius of the most massive, albeit rare,
neutron stars with masses $M\gtrsim2.1\,M_\odot$. At the same time, we find
that radius measurements of neutron stars with masses
$M\simeq1.7-1.85\,M_\odot$ can provide the strongest constraints on the
low-density part of the EOS. Finally, we quantify how radius measurements by
future missions can further improve our understanding of the EOS of matter at
nuclear densities.

[20] oai:arXiv.org:1907.07197 [pdf] - 1966857

General relativistic resistive magnetohydrodynamics with robust
primitive variable recovery for accretion disk simulations

Comments: Matching accepted version in ApJS

Submitted: **2019-07-16**, last modified: **2019-08-06**

Recent advances in black hole astrophysics, particularly the first visual
evidence of a supermassive black hole at the center of the galaxy M87 by the
Event Horizon Telescope (EHT), and the detection of an orbiting "hot spot"
nearby the event horizon of Sgr A* in the Galactic center by the Gravity
Collaboration, require the development of novel numerical methods to understand
the underlying plasma microphysics. Non-thermal emission related to such hot
spots is conjectured to originate from plasmoids that form due to magnetic
reconnection in thin current layers in the innermost accretion zone.
Resistivity plays a crucial role in current sheet formation, magnetic
reconnection, and plasmoid growth in black hole accretion disks and jets. We
included resistivity in the three-dimensional general-relativistic
magnetohydrodynamics (GRMHD) code BHAC and present the implementation of an
Implicit-Explicit scheme to treat the stiff resistive source terms of the GRMHD
equations. The algorithm is tested in combination with adaptive mesh refinement
to resolve the resistive scales and a constrained transport method to keep the
magnetic field solenoidal. Several novel methods for primitive variable
recovery, a key part in relativistic magnetohydrodynamics codes, are presented
and compared for accuracy, robustness, and efficiency. We propose a new
inversion strategy that allows for resistive-GRMHD simulations of low
gas-to-magnetic pressure ratio and highly magnetized regimes as applicable for
black hole accretion disks, jets, and neutron star magnetospheres. We apply the
new scheme to study the effect of resistivity on accreting black holes,
accounting for dissipative effects as reconnection.

[21] oai:arXiv.org:1904.04923 [pdf] - 1929749

The Event Horizon General Relativistic Magnetohydrodynamic Code
Comparison Project

Comments: Accepted version for publication in ApJS (May 28th, 2019)

Submitted: **2019-04-09**, last modified: **2019-08-05**

Recent developments in compact object astrophysics, especially the discovery
of merging neutron stars by LIGO, the imaging of the black hole in M87 by the
Event Horizon Telescope (EHT) and high precision astrometry of the Galactic
Center at close to the event horizon scale by the GRAVITY experiment motivate
the development of numerical source models that solve the equations of general
relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for
the evolution of a magnetized accretion flow where turbulence is promoted by
the magnetorotational instability from a set of nine GRMHD codes: Athena++,
BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL.
Agreement between the codes improves as resolution increases, as measured by a
consistently applied, specially developed set of code performance metrics. We
conclude that the community of GRMHD codes is mature, capable, and consistent
on these test problems.

[22] oai:arXiv.org:1906.10795 [pdf] - 1958169

Constrained transport and adaptive mesh refinement in the Black Hole
Accretion Code

Comments: 21 pages, 16 figures, accepted for publication in Astronomy &
Astrophysics. Link to www.bhac.science will soon be active

Submitted: **2019-06-25**

Worldwide very long baseline radio interferometry arrays are expected to
obtain horizon-scale images of supermassive black hole candidates as well as of
relativistic jets in several nearby active galactic nuclei. This motivates the
development of models for magnetohydrodynamic flows in strong gravitational
fields. The Black Hole Accretion Code (BHAC) intends to aid with the modelling
of such sources by means of general relativistic magnetohydrodynamical (GRMHD)
simulations in arbitrary stationary spacetimes. New additions were required to
guarantee an accurate evolution of the magnetic field when small and large
scales are captured simultaneously. We discuss the adaptive mesh refinement
(AMR) techniques employed in BHAC, essential to keep several problems
computationally tractable, as well as staggered-mesh-based constrained
transport (CT) algorithms to preserve the divergence-free constraint of the
magnetic field, including a general class of prolongation operators for
face-allocated variables compatible with them. Through several standard tests,
we show that the choice of divergence-control method can produce qualitative
differences in simulations of scientifically relevant accretion problems. We
demonstrate the ability of AMR to reduce the computational costs of accretion
simulations while sufficiently resolving turbulence from the magnetorotational
instability. In particular, we describe a simulation of an accreting Kerr black
hole in Cartesian coordinates using AMR to follow the propagation of a
relativistic jet while self-consistently including the jet engine, a problem
set up-for which the new AMR implementation is particularly advantageous. The
CT methods and AMR strategies discussed here are being employed in the
simulations performed with BHAC used in the generation of theoretical models
for the Event Horizon Telescope Collaboration.

[23] oai:arXiv.org:1906.10065 [pdf] - 2005432

Modeling non-thermal emission from the jet-launching region of M 87 with
adaptive mesh refinement

Comments: 15 pages, 10 figures

Submitted: **2019-06-24**

The galaxy M 87 harbors a kiloparsec-scale relativistic jet, whose origin
coincides with a supermassive black hole. Observational mm-VLBI campaigns are
capable of resolving the jet-launching region at the scale of the event
horizon. In order to provide a context for interpreting these observations,
realistic general-relativistic magnetohydrodynamical (GRMHD) models of the
accretion flow are constructed. The characteristics of the observed
spectral-energy distribution (SED) depend on the shape of the electrons'
energy-distribution function (eDF). The dependency on the eDF is omitted in the
modeling of the first Event Horizon Telescope results. In this work, we aim to
model the M 87 SED from radio up to NIR/optical frequencies using a
thermal-relativistic Maxwell- J\"uttner distribution, as well as a relativistic
$\kappa$-distribution function. The electrons are injected based on sub-grid,
particle-in-cell parametrizations for sub-relativistic reconnection. A GRMHD
simulation in Cartesian-Kerr-Schild coordinates, using eight levels of adaptive
mesh refinement (AMR), forms the basis of our model. To obtain spectra and
images, the GRMHD data is post-processed with the ray-tracing code RAPTOR,
which is capable of ray tracing through AMR GRMHD simulation data. We obtain
radio spectra in both the thermal-jet and $\kappa$-jet models consistent with
radio observations. Additionally, the $\kappa$-jet models also recover the
NIR/optical emission. The models recover the observed source sizes and core
shifts and obtain a jet power of $\approx 10^{43}$ ergs/s. In the $\kappa$-jet
models, both the accretion rates and jet powers are approximately two times
lower than the thermal-jet model. The frequency cut-off observed at $\nu
\approx 10^{15}$ Hz is recovered when the accelerator size is $10^6$ - $10^8$
cm, this could potentially point to an upper limit for plasmoid sizes in the
jet of M 87.

[24] oai:arXiv.org:1811.10929 [pdf] - 1897874

Constraining twin stars with GW170817

Comments: 17 pages, 12 figures, 2 tables

Submitted: **2018-11-27**, last modified: **2019-06-11**

If a phase transition is allowed to take place in the core of a compact star,
a new stable branch of equilibrium configurations can appear, providing
solutions with the same mass as the purely hadronic branch and hence giving
rise to twin-star configurations. We perform an extensive analysis of the
features of the phase transition leading twin-star configurations and, at the
same time, fulfilling the constraints coming from the maximum mass of
$2M_\odot$ and the information following gravitational-wave event GW170817. In
particular, we use a general equation of state for the neutron-star matter that
parametrizes the hadron-quark phase transition between the model describing the
hadronic phase and a constant speed of sound for the quark phase. We find that
the largest number of twin-star solutions has masses in the neutron-star branch
in the range $1-2M_\odot$ and twin-branch masses $\gtrsim 2M_\odot$. The
analysis of the masses, radii and tidal deformabilities also reveals that when
twin stars appear, the tidal deformability shows two distinct branches with the
same mass, thus differing considerably from the behaviour expected for neutron
stars. In addition, we find that the data from GW170817 is compatible with the
existence of hybrid stars and could also be interpreted as produced by the
merger of a binary system of hybrid stars or of a hybrid star with a neutron
star. The presence of a hybrid star in the inspiral phase can be established
clearly if future gravitational-wave detections measure chirp masses
$\mathcal{M}\lesssim 1.2M_\odot$ and tidal deformabilities of
$\Lambda_{1.4}\lesssim 400$ for $1.4M_\odot$ stars. Finally, combining all
observational information available, we set constraints on the parameters that
characterise the phase transition, the maximum masses, and the radii of
$1.4M_\odot$ stars described by equations of state leading to twin-star
configurations.

[25] oai:arXiv.org:1905.12658 [pdf] - 1893519

Neutron-Star-Merger Equation of State

Comments: Contribution to the Special Issue "Compact Stars in the QCD Phase
Diagram and in the Multi-Messenger Era of Astronomy" dedicated to the
conference: Compact Stars in the QCD Phase Diagram VII

Submitted: **2019-05-29**, last modified: **2019-06-03**

In this work, we discuss the dense matter equation of state (EOS) for the
extreme range of conditions encountered in neutron stars and their mergers. The
calculation of the properties of such an EOS involves modeling different
degrees of freedom (such as nuclei, nucleons, hyperons, and quarks), taking
into account different symmetries, and including finite density and temperature
effects in a thermodynamically consistent manner. We begin by addressing
subnuclear matter consisting of nucleons and a small admixture of light nuclei
in the context of the excluded volume approach. We then turn our attention to
supranuclear homogeneous matter as described by the Chiral Mean Field (CMF)
formalism. Finally, we present results from realistic neutron-star-merger
simulations performed using the CMF model that predict signatures for
deconfinement to quark matter in gravitational wave signals.

[26] oai:arXiv.org:1901.04138 [pdf] - 1890332

When Did the Remnant of GW170817 Collapse to a Black Hole?

Comments: 20 pages, 7 figures. Incorporated comments from the referee. Accepted
for publication in ApJ

Submitted: **2019-01-14**, last modified: **2019-04-10**

The main hard pulse of prompt gamma-ray emission in GRB$\,$170817A had a
duration of $\sim0.5\,{\rm s}$ and its onset was delayed with respect to the
gravitational-wave chirp signal by $t_{\rm del} \approx 1.74\,{\rm s}$.
Detailed follow-up of the subsequent broadband kilonova emission revealed a
two-component ejecta -- a lanthanide-poor ejecta with mass $M_{\rm
ej,blue}\approx0.025\,M_\odot$ that powered the early but rapidly fading blue
emission and a lanthanide-rich ejecta with mass $M_{\rm ej,red}\approx
0.04\,M_\odot$ that powered the longer lasting redder emission. Both the prompt
gamma-ray onset delay and the existence of the blue ejecta with modest electron
fraction, $0.2\lesssim Y_e\lesssim0.3$, can be explained if the collapse to a
black hole was delayed by the formation of a hypermassive neutron star (HMNS).
Here, we determine the survival time of the merger remnant by combining two
different constraints, namely, the time needed to produce the requisite
blue-ejecta mass and that necessary for the relativistic jet to bore its way
out of the expanding ejecta. In this way, we determine that the remnant of
GW170817 must have collapsed to a black hole after $t_{\rm
coll}=0.98_{-0.26}^{+0.31}\,{\rm s}$. We also discuss how future detections and
the delays between the gravitational and electromagnetic emissions can be used
to constrain the properties of the merged object.

[27] oai:arXiv.org:1904.04934 [pdf] - 1890444

Simulations of imaging the event horizon of Sagittarius A* from space

Comments: 19 pages, 14 figures, accepted for publication in A&A

Submitted: **2019-04-09**

It has been proposed that Very Long Baseline Interferometry (VLBI) at
sub-millimeter waves will allow us to image the shadow of the black hole in the
center of our Milky Way, Sagittarius A* (Sgr A*), and thereby test basic
predictions of general relativity. This paper presents imaging simulations of a
new Space VLBI mission concept. An initial design study of the concept has been
presented as the Event Horizon Imager (EHI). The EHI may be suitable for
imaging Sgr A* at high frequencies (up to ~690 GHz), which has significant
advantages over performing ground-based VLBI at 230 GHz. The concept EHI design
consists of two or three satellites in polar or equatorial circular
Medium-Earth Orbits with slightly different radii. Due to the relative drift of
the satellites along the individual orbits, this setup will result in a dense
spiral-shaped uv-coverage with long baselines (up to ~60 Glambda), allowing for
extremely high-resolution and high-fidelity imaging of radio sources. We
simulate EHI observations of general relativistic magnetohydrodynamics models
of Sgr A* and calculate the expected noise based on preliminary system
parameters. On long baselines, where the signal-to-noise ratio may be low,
fringes could be detected if the system is sufficiently phase stable and the
satellite orbits can be reconstructed with sufficient accuracy. Averaging
visibilities accumulated over multiple epochs of observations could then help
improving the image quality. With three satellites, closure phases could be
used for imaging. Our simulations show that the EHI could be capable of imaging
the black hole shadow of Sgr A* with a resolution of 4 uas (about 8% of the
shadow diameter) within several months of observing time. The EHI concept could
thus be used to measure black hole shadows much more precisely than with
ground-based VLBI, allowing for stronger tests of general relativity and
accretion models.

[28] oai:arXiv.org:1809.08682 [pdf] - 1862874

How to tell an accreting boson star from a black hole

Comments: 11 pages, 7 figures. Improved reconstruction technique, adding more
information in the form of appendices, and a comparison with a Schwarzschild
black hole. The results do not change

Submitted: **2018-09-23**, last modified: **2019-04-07**

Radio-astronomical observations of the supermassive black-hole candidate in
the galactic center will soon offer the possibility to study gravity in its
strongest regimes and to test different models for these compact objects.
Studies based on semi-analytic models and strong-field images of stationary
plasma configurations around boson stars have stressed the difficulty to
distinguish them from black holes. For the first time, we present
general-relativistic magnetohydrodynamic simulations followed by consistent
general-relativistic radiative-transfer calculations in a boson-star spacetime
to revisit the appearance of these objects in a realistic accretion scenario.
We find that the absence of an event horizon in a boson star leads to important
differences in the dynamics of the accretion and results in both the formation
of a small torus in the interior of the boson star and in the absence of an
evacuated high-magnetization funnel in the polar regions. Synthetic
reconstructed images considering realistic astronomical observing conditions
show that differences in the appearance of the two compact objects are large
enough to be detectable. These results, which also apply to other horizonless
compact objects, strengthen confidence in the ability to determine the presence
of an event horizon via radio observations and highlight the importance of
self-consistent multidimensional simulations to study the compact object at the
galactic center.

[29] oai:arXiv.org:1904.00106 [pdf] - 1949339

Using evolutionary algorithms to model relativistic jets: Application to
NGC 1052

Comments: submitted to A&A, 24 pages, 25 figures

Submitted: **2019-03-29**

High-resolution Very-Long-Baseline Interferometry observations of NGC 1052
show a two sided jet with several regions of enhanced emission and a clear
emission gap between the two jets.This gap shrinks with increasing frequency
and vanishes around $\nu\sim43$ GHz. The observed structures are due to both
the macroscopic fluid dynamics interacting with the surrounding ambient medium
including an obscuring torus and the radiation microphysics. In this paper we
investigate the possible physical conditions in relativistic jets of NGC 1052
by directly modelling the observed emission and spectra via state-of-the-art
special-relativistic hydrodynamic (SRHD) simulations and radiative transfer
calculations. To investigate the physical conditions in the relativistic jet we
coupled our radiative transfer code to evolutionary algorithms and performed
simultaneous modelling of the observed jet structure and the broadband radio
spectrum. During the calculation of the radiation we consider both thermal and
non-thermal emission. In order to compare our model to VLBI observations we
take into account the sparse sampling of the u-v plane, the array properties
and the imaging algorithm. We present for the first time an end-to-end pipeline
for fitting numerical simulations to VLBI observations of relativistic jets
taking into account the macrophysics including fluid dynamics and ambient
medium configurations together with thermal/non-thermal emission and the
properties of the observing array. The detailed analysis of our simulations
shows that the structure and properties of the observed relativistic jets in
NGC 1052 can be reconstructed by a slightly over-pressured jet ($d_k\sim1.5$)
embedded in a decreasing pressure ambient medium

[30] oai:arXiv.org:1901.09977 [pdf] - 1841523

A General-relativistic Determination of the Threshold Mass to Prompt
Collapse in Binary Neutron Star Mergers

Comments: 6 pages, 4 figures, published on ApJL

Submitted: **2019-01-28**, last modified: **2019-03-03**

We study the lifetimes of the remnant produced by the merger of two neutron
stars and revisit the determination of the threshold mass to prompt collapse,
$M_{\rm th}$. Using a fully general-relativistic numerical approach and a novel
method for a rigorous determination of $M_{\rm th}$, we show that a nonlinear
universal relation exists between the threshold mass and the maximum
compactness. For the temperature-dependent equations of state considered here,
our results improve a similar linear relation found recently with methods that
are less accurate but yield quantitatively similar results. Furthermore,
exploiting the information from GW170817, we use the universal relation to set
lower limits on the stellar radii for any mass.

[31] oai:arXiv.org:1807.03684 [pdf] - 1834012

Signatures of quark-hadron phase transitions in general-relativistic
neutron-star mergers

Comments: Version accepted by PRL

Submitted: **2018-07-10**, last modified: **2019-02-19**

Merging binaries of neutron stars are not only strong sources of
gravitational waves, but also have the potential of revealing states of matter
at densities and temperatures not accessible in laboratories. A crucial and
long-standing question in this context is whether quarks are deconfined as a
result of the dramatic increase in density and temperature following the
merger. We present the first fully general-relativistic simulations of merging
neutron stars including quarks at finite temperatures that can be switched off
consistently in the equation of state. Within our approach, we can determine
clearly what signatures a quark-hadron phase transition would leave in the
gravitational-wave signal. In particular, we show that if the conditions are
met for a phase transition to take place at several times nuclear saturation
density, they would lead to a post-merger signal considerably different from
the one expected from the inspiral, that can only probe the hadronic part of
the equations of state, and to an anticipated collapse of the merged object. We
also show that the phase transition leads to a very hot and dense quark core
that, when it collapses to a black hole, produces a ringdown signal different
from the hadronic one. Finally, in analogy with what is done in heavy-ion
collisions, we use the evolution of the temperature and density in the merger
remnant to illustrate the properties of the phase transition in a QCD phase
diagram.

[32] oai:arXiv.org:1807.03795 [pdf] - 1793560

Dynamical ejecta and nucleosynthetic yields from eccentric binary
neutron-star mergers

Comments: 24 pages, 14 figures, matches version accepted by PRD

Submitted: **2018-07-10**, last modified: **2018-12-05**

With the recent advent of multi-messenger gravitational-wave astronomy and in
anticipation of more sensitive, next-generation gravitational-wave detectors,
we investigate the dynamics, gravitational-wave emission, and nucleosynthetic
yields of numerous eccentric binary neutron-star mergers having different
equations of state. For each equation of state we vary the orbital properties
around the threshold of immediate merger, as well as the binary mass ratio. In
addition to a study of the gravitational-wave emission including $f$-mode
oscillations before and after merger, we couple the dynamical ejecta output
from the simulations to the nuclear-reaction network code \texttt{SkyNet} to
compute nucleosynthetic yields and compare to the corresponding results in the
case of a quasi-circular merger. We find that the amount and velocity of
dynamically ejected material is always much larger than in the quasi-circular
case, reaching maximal values of $M_{\rm ej, max} \sim 0.1 \, M_{\odot}$ and
$v_{\rm max}/c \sim 0.75$. At the same time, the properties of this material
are rather insensitive to the details of the orbit, such as pericenter distance
or post-encounter apoastron distance. Furthermore, while the composition of the
ejected matter depends on the orbital parameters and on the equation of state,
the relative nucleosynthetic yields do not, thus indicating that kilonova
signatures could provide information on the orbital properties of dynamically
captured neutron-star binaries.

[33] oai:arXiv.org:1811.11183 [pdf] - 1938325

Computing Fast and Reliable Gravitational Waveforms of Binary Neutron
Star Merger Remnants

Comments: Submitted to PRL. 6 pages, 4 figures

Submitted: **2018-11-27**

Gravitational waves have been detected from the inspiral of a binary
neutron-star, GW170817, which allowed constraints to be placed on the neutron
star equation of state. The equation of state can be further constrained if
gravitational waves from a post-merger remnant are detected. Post-merger
waveforms are currently generated by numerical-relativity simulations, which
are computationally expensive. Here we introduce a hierarchical model trained
on numerical-relativity simulations, which can generate reliable post-merger
spectra in a fraction of a second. Our spectra have mean fitting factors of
0.95, which compares to fitting factors of 0.76 and 0.85 between different
numerical-relativity codes that simulate the same physical system. This method
is the first step towards generating large template banks of spectra for use in
post-merger detection and parameter estimation.

[34] oai:arXiv.org:1803.09215 [pdf] - 1782727

Electrical Resistivity and Hall Effect in Binary Neutron-Star Mergers

Comments: 10 pages, 4 figures, v2: minor changes, matches published version;
v1: 9 page, 4 figures

Submitted: **2018-03-25**, last modified: **2018-11-09**

We examine the range of rest-mass densities, temperatures and magnetic fields
involved in simulations of binary neutron-star mergers and identify the
conditions under which the ideal-magnetohydrodynamics approximation breaks down
and hence the magnetic-field decay should be accounted for. We use recent
calculations of the conductivities of warm correlated plasma in envelopes of
compact stars and find that the magnetic-field decay timescales are much larger
than the characteristic timescales of the merger process for lengthscales down
to a meter. Because these are smaller than the currently available resolution
in numerical simulations, the ideal-magnetohydrodynamics approximation is
effectively valid for all realistic simulations. At the same time, we find that
the Hall effect can be important at low densities and low temperatures, where
it can induce a non-dissipative rearrangement of the magnetic field. Finally,
we mark the region in temperature and density where the hydrodynamic
description breaks down.

[35] oai:arXiv.org:1801.05705 [pdf] - 1746425

Electromagnetic emission from blitzars and its impact on non-repeating
fast radio bursts

Comments: Matches accepted version

Submitted: **2018-01-17**, last modified: **2018-09-07**

It has been suggested that a non-repeating fast radio burst (FRB) represents
the final signal of a magnetized neutron star collapsing to a black hole. In
this model, a supramassive neutron star supported by rapid rotation, will
collapse to a black hole several thousand to million years after its birth as a
result of spin down. The collapse violently snaps the magnetic-field lines
anchored on the stellar surface, thus producing an electromagnetic pulse that
will propagate outwards and accelerate electrons producing a massive radio
burst, i.e. a "blitzar". We present a systematic study of the gravitational
collapse of rotating and magnetised neutron stars with special attention to
far-field evolution at late times after the collapse. By considering a series
of neutron stars with rotation ranging from zero to millisecond periods and
different magnetic-field strengths, we show that the blitzar emission is very
robust and always characterised by a series sub-millisecond pulses decaying
exponentially in amplitude. The luminosity and energy released when the
magnetosphere is destroyed are well reproduced by a simple expression in terms
of the stellar magnetic field and radius. Finally, we assess the occurrence of
pair production during a blitzar scenario, concluding that for typical
magnetic-field strengths of $10^{12}\,{\rm G}$ and spin frequencies of a few
Hz, pair production is suppressed. Overall, the very good match between the
results of the simulations and the luminosities normally observed for FRBs
lends credibility to the blitzar model as a simple and yet plausible
explanation for the phenomenology of non-repeating FRBs.

[36] oai:arXiv.org:1803.00549 [pdf] - 1709367

New constraints on radii and tidal deformabilities of neutron stars from
GW170817

Comments: Added discussion on EOSs with phase transitions; no changes for
purely hadronic EOSs; matches version accepted by PRL

Submitted: **2018-03-01**, last modified: **2018-05-31**

We explore in a parameterized manner a very large range of physically
plausible equations of state (EOSs) for compact stars for matter that is either
purely hadronic or that exhibits a phase transition. In particular, we produce
two classes of EOSs with and without phase transitions, each containing one
million EOSs. We then impose constraints on the maximum mass, ($M < 2.16
M_{\odot}$), and on the dimensionless tidal deformability ($\tilde{\Lambda}
<800$) deduced from GW170817, together with recent suggestions of lower limits
on $\tilde{\Lambda}$. Exploiting more than $10^9$ equilibrium models for each
class of EOSs, we produce distribution functions of all the stellar properties
and determine, among other quantities, the radius that is statistically most
probable for any value of the stellar mass. In this way, we deduce that the
radius of a purely hadronic neutron star with a representative mass of
$1.4\,M_{\odot}$ is constrained to be $12.00\!<\!R_{1.4}/{\rm km}\!<\!13.45$ at
a $2$-$\sigma$ confidence level, with a most likely value of
$\bar{R}_{1.4}=12.39\,{\rm km}$; similarly, the smallest dimensionless tidal
deformability is $\tilde{\Lambda}_{1.4}\!>\!375$, again at a $2$-$\sigma$
level. On the other hand, because EOSs with a phase transition allow for very
compact stars on the so-called `twin-star' branch, small radii are possible
with such EOSs although not probable, i.e. $8.53\!<\!R_{1.4}/{\rm
km}\!<\!13.74$ and $\bar{R}_{1.4}=13.06\,{\rm km}$ at a $2$-$\sigma$ level,
with $\tilde{\Lambda}_{1.4}\!>\!35.5$ at a $3$-$\sigma$ level. Finally, since
these EOSs exhibit upper limits on $\tilde{\Lambda}$, the detection of a binary
with total mass of $3.4\,M_{\odot}$ and $\tilde{\Lambda}_{1.7}\!>\!461$ can
rule out twin-star solutions.

[37] oai:arXiv.org:1801.02839 [pdf] - 1689150

ADER discontinuous Galerkin schemes for general-relativistic ideal
magnetohydrodynamics

Comments: 23 pages, 14 figures, 6 tables

Submitted: **2018-01-09**, last modified: **2018-05-25**

We present a new class of high-order accurate numerical algorithms for
solving the equations of general-relativistic ideal magnetohydrodynamics in
curved spacetimes. In this paper we assume the background spacetime to be given
and static, i.e., we make use of the Cowling approximation. The governing
partial differential equations are solved via a new family of fully-discrete
and arbitrary high-order accurate path-conservative discontinuous Galerkin (DG)
finite-element methods combined with adaptive mesh refinement and time accurate
local timestepping. In order to deal with shock waves and other
discontinuities, the highorder DG schemes are supplemented with a novel
a-posteriori subcell finite-volume limiter, which makes the new algorithms as
robust as classical second-order total-variation diminishing finite-volume
methods at shocks and discontinuities, but also as accurate as unlimited
high-order DG schemes in smooth regions of the flow. We show the advantages of
this new approach by means of various classical two- and three-dimensional
benchmark problems on fixed spacetimes. Finally, we present a performance and
accuracy comparisons between Runge-Kutta DG schemes and ADER high-order
finite-volume schemes, showing the higher efficiency of DG schemes.

[38] oai:arXiv.org:1707.09910 [pdf] - 1678450

Conformal and covariant Z4 formulation of the Einstein equations:
strongly hyperbolic first-order reduction and solution with discontinuous
Galerkin schemes

Comments: 27 pages, 8 figures; v2 matches version published in PRD. Discussion
improved and proof of strong hyperbolicity is now given for standard gauges

Submitted: **2017-07-31**, last modified: **2018-04-28**

We present a strongly hyperbolic first-order formulation of the Einstein
equations based on the conformal and covariant Z4 system (CCZ4) with
constraint-violation damping, which we refer to as FO-CCZ4. As CCZ4, this
formulation combines the advantages of a conformal and traceless formulation,
with the suppression of constraint violations given by the damping terms, but
being first order in time and space, it is particularly suited for a
discontinuous Galerkin (DG) implementation. The strongly hyperbolic first-order
formulation has been obtained by making careful use of first and second-order
ordering constraints. A proof of strong hyperbolicity is given for a selected
choice of standard gauges via an analytical computation of the entire
eigenstructure of the FO-CCZ4 system. The resulting governing partial
differential equations system is written in non-conservative form and requires
the evolution of 58 unknowns. A key feature of our formulation is that the
first-order CCZ4 system decouples into a set of pure ordinary differential
equations and a reduced hyperbolic system of partial differential equations
that contains only linearly degenerate fields. We implement FO-CCZ4 in a
high-order path-conservative arbitrary-high-order-method-using-derivatives
(ADER)-DG scheme with adaptive mesh refinement and local time-stepping,
supplemented with a third-order ADER-WENO subcell finite-volume limiter in
order to deal with singularities arising with black holes. We validate the
correctness of the formulation through a series of standard tests in vacuum,
performed in one, two and three spatial dimensions, and also present
preliminary results on the evolution of binary black-hole systems. To the best
of our knowledge, these are the first successful three-dimensional simulations
of moving punctures carried out with high-order DG schemes using a first-order
formulation of the Einstein equations.

[39] oai:arXiv.org:1804.05812 [pdf] - 1745387

The Current Ability to Test Theories of Gravity with Black Hole Shadows

Comments: Published in Nature Astronomy on 16.04.18 (including supplementary
information); simulations at https://blackholecam.org/telling_bhs_apart/

Submitted: **2018-04-16**

Our Galactic Center, Sagittarius A* (Sgr A*), is believed to harbour a
supermassive black hole (BH), as suggested by observations tracking individual
orbiting stars. Upcoming sub-millimetre very-long-baseline-interferometry
(VLBI) images of Sgr A* carried out by the Event-Horizon-Telescope
Collaboration (EHTC) are expected to provide critical evidence for the
existence of this supermassive BH. We assess our present ability to use EHTC
images to determine if they correspond to a Kerr BH as predicted by Einstein's
theory of general relativity (GR) or to a BH in alternative theories of
gravity. To this end, we perform general-relativistic magnetohydrodynamical
(GRMHD) simulations and use general-relativistic radiative transfer (GRRT)
calculations to generate synthetic shadow images of a magnetised accretion flow
onto a Kerr BH. In addition, and for the first time, we perform GRMHD
simulations and GRRT calculations for a dilaton BH, which we take as a
representative solution of an alternative theory of gravity. Adopting the VLBI
configuration from the 2017 EHTC campaign, we find that it could be extremely
difficult to distinguish between BHs from different theories of gravity, thus
highlighting that great caution is needed when interpreting BH images as tests
of GR.

[40] oai:arXiv.org:1710.04638 [pdf] - 1697940

The THESEUS space mission concept: science case, design and expected
performances

Comments: Accepted for publication in Advances in Space Research. Partly based
on the proposal submitted on October 2016 in response to the ESA Call for
next M5 mission, with expanded and updated science sections

Submitted: **2017-10-12**, last modified: **2018-03-27**

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for
investigating the early Universe and at providing a substantial advancement of
multi-messenger and time-domain astrophysics. These goals will be achieved
through a unique combination of instruments allowing GRB and X-ray transient
detection over a broad field of view (more than 1sr) with 0.5-1 arcmin
localization, an energy band extending from several MeV down to 0.3 keV and
high sensitivity to transient sources in the soft X-ray domain, as well as
on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with
both imaging and spectroscopic capabilities. THESEUS will be perfectly suited
for addressing the main open issues in cosmology such as, e.g., star formation
rate and metallicity evolution of the inter-stellar and intra-galactic medium
up to redshift $\sim$10, signatures of Pop III stars, sources and physics of
re-ionization, and the faint end of the galaxy luminosity function. In
addition, it will provide unprecedented capability to monitor the X-ray
variable sky, thus detecting, localizing, and identifying the electromagnetic
counterparts to sources of gravitational radiation, which may be routinely
detected in the late '20s / early '30s by next generation facilities like
aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide
powerful synergies with the next generation of multi-wavelength observatories
(e.g., LSST, ELT, SKA, CTA, ATHENA).

[41] oai:arXiv.org:1712.08153 [pdf] - 1709303

THESEUS: a key space mission concept for Multi-Messenger Astrophysics

Comments: 25 pages, 13 figures, accepted to Advances in Space Research with
minor revisions. Details on the THESEUS instrumentation, science case and
expected performances can also be found in Amati et al. 2017
(arXiv:1710.04638) and in the presentations of the THESEUS Workshop 2017
(http://www.isdc.unige.ch/theseus/workshop2017-programme.html). v2 few typos
corrected

Submitted: **2017-12-21**, last modified: **2018-03-27**

The recent discovery of the electromagnetic counterpart of the gravitational
wave source GW170817, has demonstrated the huge informative power of
multi-messenger observations. During the next decade the nascent field of
multi-messenger astronomy will mature significantly. Around 2030, third
generation gravitational wave detectors will be roughly ten times more
sensitive than the current ones. At the same time, neutrino detectors currently
upgrading to multi km^3 telescopes, will include a 10 km^3 facility in the
Southern hemisphere that is expected to be operational around 2030. In this
review, we describe the most promising high frequency gravitational wave and
neutrino sources that will be detected in the next two decades. In this
context, we show the important role of the Transient High Energy Sky and Early
Universe Surveyor (THESEUS), a mission concept proposed to ESA by a large
international collaboration in response to the call for the Cosmic Vision
Programme M5 missions. THESEUS aims at providing a substantial advancement in
early Universe science as well as playing a fundamental role in multi-messenger
and time-domain astrophysics, operating in strong synergy with future
gravitational wave and neutrino detectors as well as major ground- and
space-based telescopes. This review is an extension of the THESEUS white paper
(Amati et al. 2017), also in light of the discovery of GW170817/GRB170817A that
was announced on October 16th, 2017.

[42] oai:arXiv.org:1707.09419 [pdf] - 1643579

Limits on Magnetic Field Amplification from the r-Mode Instability

Comments: This is essentially the published version

Submitted: **2017-07-28**, last modified: **2018-03-02**

At second order in perturbation theory, the unstable r-mode of a rotating
star includes growing differential rotation whose form and growth rate are
determined by gravitational-radiation reaction. With no magnetic field, the
angular velocity of a fluid element grows exponentially until the mode reaches
its nonlinear saturation amplitude and remains nonzero after saturation. With a
background magnetic field, the differential rotation winds up and amplifies the
field, and previous work where large mode amplitudes were considered suggests
that the amplification may damp out the instability. A background magnetic
field, however, turns the saturated time-independent perturbations
corresponding to adding differential rotation into perturbations whose
characteristic frequencies are of order the Alfv\'en frequency. As found in
previous studies, we argue that magnetic- field growth is sharply limited by
the saturation amplitude of an unstable mode. In contrast to previous work,
however, we show that if the amplitude is small, i.e., of order 10^(-4), then
the limit on the magnetic-field growth is stringent enough to prevent the loss
of energy to the magnetic field from damping or significantly altering an
unstable r-mode in nascent neutron stars with normal interiors and in cold
stars whose interiors are type II superconductors. We show this result first
for a toy model, and we then obtain an analogous upper limit on magnetic field
growth using a more realistic model of a rotating neutron star. Our analysis
depends on the assumption that there are no marginally unstable perturbations,
and this may not hold when differential rotation leads to a magnetorotational
instability.

[43] oai:arXiv.org:1801.10452 [pdf] - 1745375

RAPTOR I: Time-dependent radiative transfer in arbitrary spacetimes

Comments: 18 pages, 14 figures, 5 tables

Submitted: **2018-01-31**

Observational efforts to image the immediate environment of a black hole at
the scale of the event horizon benefit from the development of efficient
imaging codes that are capable of producing synthetic data, which may be
compared with observational data. We aim to present RAPTOR, a new public code
that produces accurate images, animations, and spectra of relativistic plasmas
in strong gravity by numerically integrating the equations of motion of light
rays and performing time-dependent radiative transfer calculations along the
rays. The code is compatible with any analytical or numerical spacetime. It is
hardware-agnostic and may be compiled and run both on GPUs and CPUs. We
describe the algorithms used in RAPTOR and test the code's performance. We have
performed a detailed comparison of RAPTOR output with that of other
radiative-transfer codes and demonstrate convergence of the results. We then
applied RAPTOR to study accretion models of supermassive black holes,
performing time-dependent radiative transfer through general relativistic
magneto-hydrodynamical (GRMHD) simulations and investigating the expected
observational differences between the so-called fast-light and slow-light
paradigms. Using RAPTOR to produce synthetic images and light curves of a GRMHD
model of an accreting black hole, we find that the relative difference between
fast-light and slow-light light curves is less than 5%. Using two distinct
radiative-transfer codes to process the same data, we find integrated flux
densities with a relative difference less than 0.01%. For two-dimensional GRMHD
models, such as those examined in this paper, the fast-light approximation
suffices as long as errors of a few percent are acceptable. The convergence of
the results of two different codes demonstrates that they are, at a minimum,
consistent.

[44] oai:arXiv.org:1711.00314 [pdf] - 1745365

Using gravitational-wave observations and quasi-universal relations to
constrain the maximum mass of neutron stars

Comments: 5 pages, 3 figures, matches accepted version in ApJL

Submitted: **2017-11-01**, last modified: **2018-01-03**

Combining the GW observations of merging systems of binary neutron stars and
quasi-universal relations, we set constraints on the maximum mass that can be
attained by nonrotating stellar models of neutron stars. More specifically,
exploiting the recent observation of the GW event GW 170817 and drawing from
basic arguments on kilonova modeling of GRB 170817A, together with the
quasi-universal relation between the maximum mass of nonrotating stellar models
$M_{\rm TOV}$ and the maximum mass supported through uniform rotation $M_{\rm
max}=\left(1.20^{+0.02}_{-0.05}\right) M_{\rm TOV}$ we set limits for the
maximum mass to be $ 2.01^{+0.04}_{-0.04}\leq M_{\rm TOV}/M_{\odot}\lesssim
2.16^{+0.17}_{-0.15}$, where the lower limit in this range comes from pulsar
observations. Our estimate, which follows a very simple line of arguments and
does not rely on the modeling of the electromagnetic signal in terms of
numerical simulations, can be further refined as new detections become
available. We briefly discuss the impact that our conclusions have on the
equation of state of nuclear matter.

[45] oai:arXiv.org:1712.09224 [pdf] - 1618390

Magnetic field amplification by the r-mode instability

Comments: 5 pages; Contribution to the proceedings of the Physics of Neutron
Stars - 2017 (Saint Petersburg, July 10-14, 2017), associated with a talk by
A.I. Chugunov

Submitted: **2017-12-26**

We discuss magnetic field enhancement by unstable r-modes (driven by the
gravitational radiation reaction force) in rotating stars. In the absence of a
magnetic field, gravitational radiation exponentially increases the r-mode
amplitude $\alpha$, and accelerates differential rotation (secular motion of
fluid elements). For a magnetized star, differential rotation enhances the
magnetic field energy. Rezzolla et al. (2000--2001) argued that if the magnetic
energy grows faster than the gravitational radiation reaction force pumps
energy into the r-modes, then the r-mode instability is suppressed. Chugunov
(2015) demonstrated that without gravitational radiation, differential rotation
can be treated as a degree of freedom decoupled from the r-modes and controlled
by the back reaction of the magnetic field. In particular, the magnetic field
windup does not damp r-modes. Here we discuss the effect of the back reaction
of the magnetic field on differential rotation of unstable r-modes, and show
that it limits the generated magnetic field and the magnetic energy growth rate
preventing suppression of the r-mode instability by magnetic windup at low
saturation amplitudes, $\alpha \ll 1$, predicted by current models.

[46] oai:arXiv.org:1711.00198 [pdf] - 1635348

Uniformly rotating, axisymmetric and triaxial quark stars in general
relativity

Comments: 14 pages, 9 figures; accepted for publication in PRD

Submitted: **2017-11-01**, last modified: **2017-12-20**

Quasi-equilibrium models of uniformly rotating axisymmetric and triaxial
quark stars are computed in general relativistic gravity scenario. The
Isenberg-Wilson-Mathews (IWM) formulation is employed and the Compact Object
CALculator (COCAL) code is extended to treat rotating stars with finite surface
density and new equations of state (EOSs). Besides the MIT bag model for quark
matter which is composed of de-confined quarks, we examine a new EOS proposed
by Lai and Xu that is based on quark clustering and results in a stiff EOS that
can support masses up to $3.3M_\odot$ in the case we considered. We perform
convergence tests for our new code to evaluate the effect of finite surface
density in the accuracy of our solutions and construct sequences of solutions
for both small and high compactness. The onset of secular instability due to
viscous dissipation is identified and possible implications are discussed. An
estimate of the gravitational wave amplitude and luminosity based on quadrupole
formulas is presented and comparison with neutron stars is discussed.

[47] oai:arXiv.org:1712.00265 [pdf] - 1757226

Test-particle dynamics in general spherically symmetric black hole
spacetimes

Comments: 15 pages, 5 figures, 2 Tables

Submitted: **2017-12-01**

To date, the most precise tests of general relativity have been achieved
through pulsar timing, albeit in the weak-field regime. Since pulsars are some
of the most precise and stable "clocks" in the Universe, present observational
efforts are focused on detecting pulsars in the vicinity of supermassive black
holes (most notably in our Galactic Centre), enabling pulsar timing to be used
as an extremely precise probe of strong-field gravity. In this paper a
mathematical framework to describe test-particle dynamics in general black hole
spacetimes is presented, and subsequently used to study a binary system
comprising a pulsar orbiting a black hole. In particular, taking into account
the parameterization of a general spherically symmetric black hole metric,
general analytic expressions for both the advance of the periastron and for the
orbital period of a massive test particle are derived. Furthermore, these
expressions are applied to four representative cases of solutions arising in
both general relativity and in alternative theories of gravity. Finally, this
framework is applied to the Galactic Centre $S$-stars and four distinct pulsar
toy models. It is shown that by adopting a fully general-relativistic
description of test-particle motion which is independent of any particular
theory of gravity, observations of pulsars can help impose better constraints
on alternative theories of gravity than is presently possible.

[48] oai:arXiv.org:1709.06058 [pdf] - 1616565

On the stability and maximum mass of differentially rotating
relativistic stars

Comments: 5 pages, 7 figures, matches accepted version in MNRAS

Submitted: **2017-09-18**, last modified: **2017-11-06**

The stability properties of rotating relativistic stars against prompt
gravitational collapse to a black hole are rather well understood for uniformly
rotating models. This is not the case for differentially rotating neutron
stars, which are expected to be produced in catastrophic events such as the
merger of binary system of neutron stars or the collapse of a massive stellar
core. We consider sequences of differentially rotating equilibrium models using
the $j$-constant law and by combining them with their dynamical evolution, we
show that a sufficient stability criterion for differentially rotating neutron
stars exists similar to the one of their uniformly rotating counterparts.
Namely: along a sequence of constant angular momentum, a dynamical instability
sets in for central rest-mass densities slightly below the one of the
equilibrium solution at the turning point. In addition, following Breu &
Rezzolla (2016), we show that "quasi-universal" relations can be found when
calculating the turning-point mass. In turn, this allows us to compute the
maximum mass allowed by differential rotation, $M_{\rm
max,dr}$, in terms of the maximum mass of the nonrotating configuration,
$M_{_{\rm TOV}}$, finding that $M_{\rm max, dr} \simeq \left(1.54 \pm
0.05\right) M_{_{\rm TOV}}$ for all the equations of state we have considered.

[49] oai:arXiv.org:1709.09630 [pdf] - 1602654

On r-process nucleosynthesis from matter ejected in binary neutron
starmergers

Comments: 24 pages, 18 figures, submitted to PRD

Submitted: **2017-09-27**

When binary systems of neutron stars merge, a very small fraction of their
rest mass is ejected, either dynamically or secularly. This material is
neutron-rich and its nucleosynthesis could provide the astrophysical site for
the production of heavy elements in the universe, together with a kilonova
signal confirming neutron-star mergers as the origin of short gamma-ray bursts.
We perform full general-relativistic simulations of binary neutron-star mergers
employing three different nuclear-physics EOSs, considering both equal- and
unequal-mass configurations, and adopting a leakage scheme to account for
neutrino radiative losses. Using a combination of techniques, we carry out an
extensive and systematic study of the hydrodynamical, thermodynamical, and
geometrical properties of the matter ejected dynamically, employing the WinNet
nuclear-reaction network to recover the relative abundances of heavy elements
produced by each configurations. Among the results obtained, three are
particularly important. First, we find that both the properties of the
dynamical ejecta and the nucleosynthesis yields are robust against variations
of the EOS and masses, and match very well the observed chemical abundances.
Second, using a conservative but robust criterion for unbound matter, we find
that the amount of ejected mass is $\lesssim 10^{-3}\,M_{\odot}$, hence at
least one order of magnitude smaller than what normally assumed in modelling
kilonova signals. Finally, using a simplified and gray-opacity model we assess
the observability of the infrared kilonova emission finding, that for all
binaries the luminosity peaks around $\sim1/2$ day in the $H$-band, reaching a
maximum magnitude of $-13$, and decreasing rapidly after one day. These rather
low luminosities make the prospects for detecting kilonovae less promising than
what assumed so far.

[50] oai:arXiv.org:1709.07706 [pdf] - 1588731

The kilo Hertz quasi-periodic oscillations in neutron star low-mass
X-ray binaries as tori oscillation modes. I

Comments: 9 pages, 5 figures

Submitted: **2017-09-22**

There have been many efforts to explain the dynamical mechanisms behind the
phenomenology of quasi-periodic oscillations (QPOs) seen in the X-ray light
curves of low-mass X-ray binaries. Up to now, none of the models can
successfully explain all the frequencies observed in the power density spectrum
of the light curve. Here we perform several hydrodynamic simulations of
non-self-gravitating relativistic axisymmetric thick tori applied to the
neutron star in the low-mass X-ray binary 4U 1636-53 and show how the observed
oscillation modes triggered by different velocity perturbations give rise to a
set of variability features similar to what we see in the observational X-ray
data. When we match pairs of frequencies from the simulations of constant
angular momentum ($l$) tori with the observed kilo Hertz QPOs, we find that
certain combinations of frequencies lie on top of the observed relation,
provided we assume a smaller mass for the neutron star than is generally
assumed. However, constant-$l$ tori cannot match the entire range of
frequencies observed for 4U 1636-53 due to physical constraints set by the
torus size. We show that our model is consistent with the observed shift in QPO
frequency that accompany state transitions of the accretion disk.

[51] oai:arXiv.org:1707.09475 [pdf] - 1626276

On the importance of viscous dissipation and heat conduction in binary
neutron-star mergers

Comments: 6 pages, 4 figures

Submitted: **2017-07-29**

Inferring the properties of dense matter is one of the most exciting
prospects from the measurement of gravitational waves from neutron star
mergers. However, it will require reliable numerical simulations that
incorporate viscous dissipation and energy transport if these can play a
significant role within the survival time of the post-merger object. We
calculate timescales for typical forms of dissipation and find that thermal
transport and shear viscosity will not be important unless neutrino trapping
occurs, which requires temperatures above about 10 MeV and gradients over
lengthscales of 0.1 km or less. On the other hand, if direct-Urca processes
remain suppressed, leaving modified-Urca processes to establish flavor
equilibrium, then bulk viscous dissipation could provide significant damping to
density oscillations observed right after the merger. When comparing with data
from a state-of-the-art merger simulation we find that the bulk viscosity takes
values close to its resonant maximum in a typical neutron-star merger,
motivating a more careful assessment of the role of bulk viscous dissipation in
the gravitational-wave signal from merging neutron stars.

[52] oai:arXiv.org:1611.07152 [pdf] - 1580675

Rotational properties of hypermassive neutron stars from binary mergers

Comments: 25 pages, 20 figures

Submitted: **2016-11-22**, last modified: **2017-07-28**

Determining the differential-rotation law of compact stellar objects produced
in binary neutron stars mergers or core-collapse supernovae is an old problem
in relativistic astrophysics. Addressing this problem is important because it
impacts directly on the maximum mass these objects can attain and hence on the
threshold to black-hole formation under realistic conditions. Using the results
from a large number of numerical simulations in full general relativity of
binary neutron star mergers described with various equations of state and
masses, we study the rotational properties of the resulting hypermassive
neutron stars. We find that the angular-velocity distribution shows only a
modest dependence on the equation of state, thus exhibiting the traits of
"quasi-universality" found in other aspects of compact stars, both isolated and
in binary systems. The distributions are characterized by an almost uniformly
rotating core and a "disk". Such a configuration is significantly different
from the $j-{\rm constant}$ differential-rotation law that is commonly adopted
in equilibrium models of differentially rotating stars. Furthermore, the
rest-mass contained in such a disk can be quite large, ranging from $\simeq
0.03\,M_{\odot}$ in the case of high-mass binaries with stiff equations of
state, up to $\simeq 0.2\,M_{\odot}$ for low-mass binaries with soft equations
of state. We comment on the astrophysical implications of our findings and on
the long-term evolutionary scenarios that can be conjectured on the basis of
our simulations.

[53] oai:arXiv.org:1611.09720 [pdf] - 1580729

The Black Hole Accretion Code

Comments: revised version accepted for publication by Computational
Astrophysics and Cosmology. Some figures have been down-sampled

Submitted: **2016-11-29**, last modified: **2017-06-18**

We present the black hole accretion code (BHAC), a new multidimensional
general-relativistic magnetohydrodynamics module for the MPI-AMRVAC framework.
BHAC has been designed to solve the equations of ideal general-relativistic
magnetohydrodynamics in arbitrary spacetimes and exploits adaptive mesh
refinement techniques with an efficient block-based approach. Several
spacetimes have already been implemented and tested. We demonstrate the
validity of BHAC by means of various one-, two-, and three-dimensional test
problems, as well as through a close comparison with the HARM3D code in the
case of a torus accreting onto a black hole. The convergence of a turbulent
accretion scenario is investigated with several diagnostics and we find
accretion rates and horizon-penetrating fluxes to be convergent to within a few
percent when the problem is run in three dimensions. Our analysis also involves
the study of the corresponding thermal synchrotron emission, which is performed
by means of a new general-relativistic radiative transfer code, BHOSS. The
resulting synthetic intensity maps of accretion onto black holes are found to
be convergent with increasing resolution and are anticipated to play a crucial
role in the interpretation of horizon-scale images resulting from upcoming
radio observations of the source at the Galactic Center.

[54] oai:arXiv.org:1705.10850 [pdf] - 1622211

Neutron-star Radius from a Population of Binary Neutron Star Mergers

Comments: 9 pages and 7 figures

Submitted: **2017-05-30**

We show how gravitational-wave observations with advanced detectors of tens
to several tens of neutron-star binaries can measure the neutron-star radius
with an accuracy of several to a few percent, for mass and spatial
distributions that are realistic, and with none of the sources located within
100 Mpc. We achieve such an accuracy by combining measurements of the total
mass from the inspiral phase with those of the compactness from the postmerger
oscillation frequencies. For estimating the measurement errors of these
frequencies we utilize analytical fits to postmerger numerical-relativity
waveforms in the time domain, obtained here for the first time, for four
nuclear-physics equations of state and a couple of values for the mass. We
further exploit quasi-universal relations to derive errors in compactness from
those frequencies. Measuring the average radius to well within 10% is possible
for a sample of 100 binaries distributed uniformly in volume between 100 and
300 Mpc, so long as the equation of state is not too soft or the binaries are
not too heavy.

[55] oai:arXiv.org:1705.01300 [pdf] - 1618274

Jet-torus connection in radio galaxies: Relativistic hydrodynamics and
synthetic emission

Comments: 14 pages, 17 figures, submitted to A&A

Submitted: **2017-05-03**

High-resolution Very-Long-Baseline Interferometry observations of active
galactic nuclei have revealed asymmetric structures in the jets of radio
galaxies. These asymmetric structures may be due to internal asymmetries in the
jet, could be induced by the different conditions in the surrounding ambient
medium including the obscuring torus, or a combination of the two. In this
paper we investigate the influence of the ambient medium (including the
obscuring torus) on the observed properties of jets from radio galaxies. We
performed special-relativistic hydrodynamic (RHD) simulations of over-pressured
and pressure-matched jets using the special-relativistic hydrodynamics code
\texttt{Ratpenat}, which is based on a second-order accurate finite-volume
method and an approximate Riemann solver. Using a newly developed emission code
to compute the electromagnetic emission, we have investigated the influence of
different ambient medium and torus configurations on the jet structure and
subsequently computed the non-thermal emission produced by the jet and the
thermal absorption due to the torus. To better compare the emission simulations
with observations we produced synthetic radio maps, taking into account the
properties of the observatory. The detailed analysis of our simulations shows
that the observed asymmetries can be produced by the interaction of the jet
with the ambient medium and by the absorption properties of the obscuring
torus.

[56] oai:arXiv.org:1607.03540 [pdf] - 1637447

Binary neutron-star mergers: a review of Einstein's richest laboratory

Comments: 106 pages; invited review from Reports on Progress of Physics; small
changes, matches version to be published

Submitted: **2016-07-12**, last modified: **2017-03-22**

The merger of binary neutron-stars systems combines in a single process:
extreme gravity, copious emission of gravitational waves, complex microphysics,
and electromagnetic processes that can lead to astrophysical signatures
observable at the largest redshifts. We review here the recent progress in
understanding what could be considered Einstein's richest laboratory,
highlighting in particular the numerous significant advances of the last
decade. Although special attention is paid to the status of models, techniques,
and results for fully general-relativistic dynamical simulations, a review is
also offered on initial data and advanced simulations with approximate
treatments of gravity. Finally, we review the considerable amount of work
carried out on the post-merger phase, including: black-hole formation, torus
accretion onto the merged compact object, connection with gamma-ray burst
engines, ejected material, and its nucleosynthesis.

[57] oai:arXiv.org:1703.03223 [pdf] - 1574672

Gravitational collapse to a Kerr-Newman black hole

Comments: MNRAS accepted

Submitted: **2017-03-09**

We present the first systematic study of the gravitational collapse of
rotating and magnetised neutron stars to charged and rotating (Kerr-Newman)
black holes. In particular, we consider the collapse of magnetised and rotating
neutron stars assuming that no pair-creation takes place and that the charge
density in the magnetosphere is so low that the stellar exterior can be
described as an electrovacuum. Under these assumptions, which are rather
reasonable for a pulsar that has crossed the 'death line', we show that when
the star is rotating, it acquires a net initial electrical charge, which is
then trapped inside the apparent horizon of the newly formed back hole. We
analyse a number of different quantities to validate that the black hole
produced is indeed a Kerr-Newman one and show that, in the absence of rotation
or magnetic field, the end result of the collapse is a Schwarzschild or Kerr
black hole, respectively.

[58] oai:arXiv.org:1606.08879 [pdf] - 1530906

BlackHoleCam: fundamental physics of the Galactic center

Comments: review paper, 36 pages, 12 figures (v2 is the published version)

Submitted: **2016-06-28**, last modified: **2017-02-07**

Einstein's General Theory of Relativity (GR) successfully describes gravity.
The most fundamental predictions of GR are black holes (BHs), but in spite of
many convincing BH candidates in the Universe, there is no conclusive
experimental proof of their existence using astronomical observations in the
electromagnetic spectrum. Are BHs real astrophysical objects? Does GR hold in
its most extreme limit or are alternatives needed? The prime target to address
these fundamental questions is in the center of our own Galaxy, which hosts the
closest and best-constrained supermassive BH candidate in the Universe,
Sagittarius A* (Sgr A*). Three different types of experiments hold the promise
to test GR in a strong-field regime using observations of Sgr A* with
new-generation instruments. The first experiment aims to image the relativistic
plasma emission which surrounds the event horizon and forms a "shadow" cast
against the background, whose predicted size (~50 microarcseconds) can now be
resolved by upcoming VLBI experiments at mm-waves such as the Event Horizon
Telescope (EHT). The second experiment aims to monitor stars orbiting Sgr A*
with the upcoming near-infrared interferometer GRAVITY at the Very Large
Telescope (VLT). The third experiment aims to time a radio pulsar in tight
orbit about Sgr A* using radio telescopes (including the Atacama Large
Millimeter Array or ALMA). The BlackHoleCam project exploits the synergy
between these three different techniques and aims to measure the main BH
parameters with sufficient precision to provide fundamental tests of GR and
probe the spacetime around a BH in any metric theory of gravity. Here, we
review our current knowledge of the physical properties of Sgr A* as well as
the current status of such experimental efforts towards imaging the event
horizon, measuring stellar orbits, and timing pulsars around Sgr A*.

[59] oai:arXiv.org:1606.08192 [pdf] - 1530897

Simulations of recoiling black holes: adaptive mesh refinement and
radiative transfer

Comments: 20 pages, 20 figures, accepted for publication in A&A

Submitted: **2016-06-27**, last modified: **2016-11-10**

(Abridged) We here continue our effort to model the behaviour of matter when
orbiting or accreting onto a generic black hole by developing a new numerical
code employing advanced techniques geared solve the equations of in
general-relativistic hydrodynamics. The new code employs a number of
high-resolution shock-capturing Riemann-solvers and reconstruction algorithms,
exploiting the enhanced accuracy and the reduced computational cost of AMR
techniques. In addition, the code makes use of sophisticated ray-tracing
libraries that, coupled with general-relativistic radiation-transfer
calculations, allow us to compute accurately the electromagnetic emissions from
such accretion flows. We validate the new code by presenting an extensive
series of stationary accretion flows either in spherical or axial symmetry and
performed either in 2D or 3D. In addition, we consider the highly nonlinear
scenario of a recoiling black hole produced in the merger of a supermassive
black hole binary interacting with the surrounding circumbinary disc. In this
way we can present, for the first time, ray-traced images of the shocked fluid
and the light-curve resulting from consistent general-relativistic
radiation-transport calculations from this process. The work presented here
lays the ground for the development of a generic computational infrastructure
employing AMR techniques to deal accurately and self-consistently with
accretion flows onto compact objects. In addition to the accurate handling of
the matter, we provide a self-consistent electromagnetic emission from these
scenarios by solving the associated radiative-transfer problem. While magnetic
fields are presently excluded from our analysis, the tools presented here can
have a number of applications to study accretion flows onto black holes or
neutron stars.

[60] oai:arXiv.org:1607.05767 [pdf] - 1502008

A new method for shadow calculations: application to parameterised
axisymmetric black holes

Comments: 14 pages, 5 figures, 4 tables; Accepted for publication in PRD

Submitted: **2016-07-19**, last modified: **2016-10-15**

Collaborative international efforts under the name of the Event Horizon
Telescope project, using sub- mm very long baseline interferometry, are soon
expected to provide the first images of the shadow cast by the candidate
supermassive black hole in our Galactic center, Sagittarius A*. Observations of
this shadow would provide direct evidence of the existence of astrophysical
black holes. Although it is expected that astrophysical black holes are
described by the axisymmetric Kerr solution, there also exist many other black
hole solutions, both in general relativity and in other theories of gravity,
which cannot presently be ruled out. To this end, we present calculations of
black hole shadow images from various metric theories of gravity as described
by our recent work on a general parameterisation of axisymmetric black holes
[R. Konoplya, L. Rezzolla and A. Zhidenko, Phys. Rev. D 93, 064015 (2016)]. An
algorithm to perform general ray-tracing calculations for any metric theory of
gravity is first outlined and then employed to demonstrate that even for
extremal metric deformation parameters of various black hole spacetimes, this
parameterisation is both robust and rapidly convergent to the correct solution.

[61] oai:arXiv.org:1602.08759 [pdf] - 1498264

Did GW150914 produce a rotating gravastar?

Comments: 6 pages, 4 figures. Extended and improved discussion of the ringdown.
Conclusions unchanged and strengthened; matches version accepted for
publication in PRD

Submitted: **2016-02-28**, last modified: **2016-09-21**

The interferometric LIGO detectors have recently measured the first direct
gravitational-wave signal from what has been interpreted as the inspiral,
merger and ringdown of a binary system of black holes. The signal-to-noise
ratio of the measured signal is large enough to leave little doubt that it does
refer to the inspiral of two massive and ultracompact objects, whose merger
yields a rotating black hole. Yet, the quality of the data is such that some
room is left for alternative interpretations that do not involve black holes,
but other objects that, within classical general relativity, can be equally
massive and compact, namely, gravastars. We here consider the hypothesis that
the merging objects were indeed gravastars and explore whether the merged
object could therefore be not a black hole but a rotating gravastar. After
comparing the real and imaginary parts of the ringdown signal of GW150914 with
the corresponding quantities for a variety of gravastars, and notwithstanding
the very limited knowledge of the perturbative response of rotating gravastars,
we conclude it is not possible to model the measured ringdown of GW150914 as
due to a rotating gravastar.

[62] oai:arXiv.org:1609.04828 [pdf] - 1531614

The Properties of Short gamma-ray burst Jets Triggered by neutron star
mergers

Comments: 6 pages, 5 figures, Submitted to ApJL

Submitted: **2016-09-15**

The most popular model for short gamma-ray bursts (sGRBs) involves the
coalescence of binary neutron stars. Because the progenitor is actually hidden
from view, we must consider under which circumstances such merging systems are
capable of producing a successful sGRB. Soon after coalescence, winds are
launched from the merger remnant. In this paper, we use realistic wind profiles
derived from global merger simulations in order to investigate the interaction
of sGRB jets with these winds using numerical simulations. We analyze the
conditions for which these axisymmetric winds permit relativistic jets to
breakout and produce a sGRB. We find that jets with luminosities comparable to
those observed in sGRBs are only successful when their half-opening angles are
below ~20{\deg}. This jet collimation mechanism leads to a simple physical
interpretation of the luminosities and opening angles inferred for sGRBs. If
wide, low luminosity jets are observed, they might be indicative of a different
progenitor avenue such as the merger of a neutron star with a black hole. We
also use the observed durations of sGRB to place constraints on the lifetime of
the wind phase, which is determined by the time it takes the jet to breakout.
In all cases we find that the derived limits argue against completely stable
remnants for binary neutron star mergers that produce sGRBs.

[63] oai:arXiv.org:1605.01938 [pdf] - 1436376

The final spin from binary black holes in quasi-circular orbits

Comments: 6 pages, 3 figures. Matches ApJL accepted version. Coefficients
$\epsilon_\alpha$, $\epsilon_\beta$ and $\epsilon_\gamma$ slightly changed to
reflect dataset filtered of imprecise NR data; contains comparison with newly
published NR catalog

Submitted: **2016-05-06**, last modified: **2016-06-07**

We revisit the problem of predicting the spin magnitude and direction of the
black hole resulting from the merger of two black holes with arbitrary masses
and spins inspiralling in quasi-circular orbits. We do this by analyzing a
catalog of 619 recent numerical-relativity simulations collected from the
literature and spanning a large variety of initial conditions. By combining
information from the post-Newtonian approximation, the extreme mass-ratio limit
and perturbative calculations, we improve our previously proposed
phenomenological formulae for the final remnant spin. In contrast with
alternative suggestions in the literature, and in analogy with our previous
expressions, the new formula is a simple algebraic function of the initial
system parameters and is not restricted to binaries with spins
aligned/anti-aligned with the orbital angular momentum, but can be employed for
fully generic binaries. The accuracy of the new expression is significantly
improved, especially for almost extremal progenitor spins and for small mass
ratios, yielding a root-mean-square error $\sigma\approx0.002$ for
aligned/anti-aligned binaries and $\sigma\approx0.006$ for generic binaries.
Our new formula is suitable for cosmological applications and can be employed
robustly in the analysis of the gravitational waveforms from advanced
interferometric detectors.

[64] oai:arXiv.org:1604.00246 [pdf] - 1429465

Gravitational-wave signal from binary neutron stars: a systematic
analysis of the spectral properties

Comments: 20 pages, 11 figures, 2 tables; Minor revisions, matches accepted
version on PRD

Submitted: **2016-04-01**, last modified: **2016-06-01**

A number of works have shown that important information on the equation of
state of matter at nuclear density can be extracted from the gravitational
waves emitted by merging neutron-star binaries. We present a comprehensive
analysis of the gravitational-wave signal emitted during the inspiral, merger
and post-merger of 56 neutron-star binaries. This sample of binaries, arguably
the largest studied to date with realistic equations of state, spans across six
different nuclear-physics equations of state and ten masses, allowing us to
sharpen a number of results recently obtained on the spectral properties of the
gravitational-wave signal. Overall we find that: (i) for binaries with masses
differing no more than $20\%$, the frequency at gravitational-wave amplitude's
maximum is related quasi-universally with the tidal deformability of the two
stars; (ii) the spectral properties vary during the post-merger phase, with a
transient phase lasting a few millisecond after the merger and followed by a
quasi-stationary phase; (iii) when distinguishing the spectral peaks between
these two phases, a number of ambiguities in the identification of the peaks
disappear, leaving a simple and robust picture; (iv) using properly identified
frequencies, quasi-universal relations are found between the spectral features
and the properties of the neutron stars; (v) for the most salient peaks
analytic fitting functions can be obtained in terms of the stellar tidal
deformability or compactness. Altogether, these results support the idea that
the equation of state of nuclear matter can be constrained tightly when a
signal in gravitational waves from binary neutron stars is detected.

[65] oai:arXiv.org:1601.02426 [pdf] - 1425302

Dynamical Mass Ejection from Binary Neutron Star Mergers

Comments: 19 pages, 10 figures. We corrected a problem in the formulation of
the neutrino heating scheme and re-ran all of the affected models. The main
conclusions are unchanged. This version also contains one more figure and a
number of improvements on the text

Submitted: **2016-01-11**, last modified: **2016-05-12**

We present fully general-relativistic simulations of binary neutron star
mergers with a temperature and composition dependent nuclear equation of state.
We study the dynamical mass ejection from both quasi-circular and
dynamical-capture eccentric mergers. We systematically vary the level of our
treatment of the microphysics to isolate the effects of neutrino cooling and
heating and we compute the nucleosynthetic yields of the ejecta. We find that
eccentric binaries can eject significantly more material than quasi-circular
binaries and generate bright infrared and radio emission. In all our
simulations the outflow is composed of a combination of tidally- and
shock-driven ejecta, mostly distributed over a broad $\sim 60^\circ$ angle from
the orbital plane, and, to a lesser extent, by thermally driven winds at high
latitudes. Ejecta from eccentric mergers are typically more neutron rich than
those of quasi-circular mergers. We find neutrino cooling and heating to
affect, quantitatively and qualitatively, composition, morphology, and total
mass of the outflows. This is also reflected in the infrared and radio
signatures of the binary. The final nucleosynthetic yields of the ejecta are
robust and insensitive to input physics or merger type in the regions of the
second and third r-process peaks. The yields for elements on the first peak
vary between our simulations, but none of our models is able to explain the
Solar abundances of first-peak elements without invoking additional first-peak
contributions from either neutrino and viscously-driven winds operating on
longer timescales after the mergers, or from core-collapse supernovae.

[66] oai:arXiv.org:1605.01709 [pdf] - 1416837

Electromagnetic fields in the exterior of an oscillating relativistic
star -- II. Electromagnetic damping

Comments: 21 pages, 1 figure; for paper I see arXiv:gr-qc/0406018; MNRAS in
press

Submitted: **2016-05-05**

An important issue in the asteroseismology of compact and magnetized stars is
the determination of the dissipation mechanism which is most efficient in
damping the oscillations when these are produced. In a linear regime and for
low-multipolarity modes these mechanisms are confined to either
gravitational-wave or electromagnetic losses. We here consider the latter and
compute the energy losses in the form of Poynting fluxes, Joule heating and
Ohmic dissipation in a relativistic oscillating spherical star with a dipolar
magnetic field in vacuum. While this approach is not particularly realistic for
rapidly rotating stars, it has the advantage that it is fully analytic and that
it provides expressions for the electric and magnetic fields produced by the
most common modes of oscillation both in the vicinity of the star and far away
from it. In this way we revisit and extend to a relativistic context the
classical estimates of McDermott et al. Overall, we find that
general-relativistic corrections lead to electromagnetic damping time-scales
that are at least one order of magnitude smaller than in Newtonian gravity.
Furthermore, with the only exception of $g$ (gravity) modes, we find that $f$
(fundamental), $p$ (pressure), $i$ (interface) and $s$ (shear) modes are
suppressed more efficiently by gravitational losses than by electromagnetic
ones.

[67] oai:arXiv.org:1602.04953 [pdf] - 1396871

Bjorken flow in one-dimensional relativistic magnetohydrodynamics with
magnetization

Comments: 12 pages, 9 figures; minor revisions, two references added, journal
version accepted by PRD

Submitted: **2016-02-16**, last modified: **2016-04-13**

We study the one-dimensional, longitudinally boost-invariant motion of an
ideal fluid with infinite conductivity in the presence of a transverse magnetic
field, i.e., in the ideal transverse magnetohydrodynamical limit. In an
extension of our previous work Roy et al., [Phys. Lett. B 750, 45 (2015)], we
consider the fluid to have a non-zero magnetization. First, we assume a
constant magnetic susceptibility $\chi_{m}$ and consider an ultrarelativistic
ideal gas equation of state. For a paramagnetic fluid (i.e., with
$\chi_{m}>0$), the decay of the energy density slows down since the fluid gains
energy from the magnetic field. For a diamagnetic fluid (i.e., with
$\chi_{m}<0$), the energy density decays faster because it feeds energy into
the magnetic field. Furthermore, when the magnetic field is taken to be
external and to decay in proper time $\tau$ with a power law $\sim\tau^{-a}$,
two distinct solutions can be found depending on the values of $a$ and
$\chi_m$. Finally, we also solve the ideal magnetohydrodynamical equations for
one-dimensional Bjorken flow with a temperature-dependent magnetic
susceptibility and a realistic equation of state given by lattice-QCD data. We
find that the temperature and energy density decay more slowly because of the
non-vanishing magnetization. For values of the magnetic field typical for
heavy-ion collisions, this effect is, however, rather small. Only for magnetic
fields which are about an order of magnitude larger than expected for heavy-ion
collisions, the system is substantially reheated and the lifetime of the quark
phase might be extended.

[68] oai:arXiv.org:1601.06083 [pdf] - 1388167

Maximum mass, moment of inertia and compactness of relativistic stars

Comments: v2 matches version published on MNRAS

Submitted: **2016-01-22**, last modified: **2016-04-11**

A number of recent works have highlighted that it is possible to express the
properties of general-relativistic stellar equilibrium configurations in terms
of functions that do not depend on the specific equation of state employed to
describe matter at nuclear densities. These functions are normally referred to
as "universal relations" and have been found to apply, within limits, both to
static or stationary isolated stars, as well as to fully dynamical and merging
binary systems. Further extending the idea that universal relations can be
valid also away from stability, we show that a universal relation is exhibited
also by equilibrium solutions that are not stable. In particular, the mass of
rotating configurations on the turning-point line shows a universal behaviour
when expressed in terms of the normalised Keplerian angular momentum. In turn,
this allows us to compute the maximum mass allowed by uniform rotation,
M_{max}, simply in terms of the maximum mass of the nonrotating configuration,
M_{TOV}, finding that M_{max} ~ (1.203 +- 0.022) M_{TOV} for all the equations
of state we have considered. We further show that a universal relation can be
found between the dimensionless moment of inertia and the stellar compactness.
Although this relation is not surprising as it involves two quantities that
have been shown to exhibit universal behaviour with other stellar properties,
our parameterisation represents a refinement over a similar relation by
Lattimer and Schutz (2005), where a different normalisation was used, and could
provide an accurate tool to constrain the equation of state of nuclear matter
when measurements of the moment of inertia become available.

[69] oai:arXiv.org:1602.02378 [pdf] - 1372043

General parametrization of axisymmetric black holes in metric theories
of gravity

Comments: 17 pages, to be published in PRD

Submitted: **2016-02-07**

Following previous work of ours in spherical symmetry, we here propose a new
parametric framework to describe the spacetime of axisymmetric black holes in
generic metric theories of gravity. In this case, the metric components are
functions of both the radial and the polar angular coordinates, forcing a
double expansion to obtain a generic axisymmetric metric expression. In
particular, we use a continued-fraction expansion in terms of a compactified
radial coordinate to express the radial dependence, while we exploit a Taylor
expansion in terms of the cosine of the polar angle for the polar dependence.
These choices lead to a superior convergence in the radial direction and to an
exact limit on the equatorial plane. As a validation of our approach, we build
parametrized representations of Kerr, rotating dilaton, and
Einstein-dilaton-Gauss-Bonnet black holes. The match is already very good at
lowest order in the expansion and improves as new orders are added. We expect a
similar behavior for any stationary and axisymmetric black-hole metric.

[70] oai:arXiv.org:1502.02021 [pdf] - 1309537

General-relativistic resistive-magnetohydrodynamic simulations of binary
neutron stars

Comments: 26 pages, 17 figures; animations available at
http://www.southampton.ac.uk/~kd10g13/movies/index.shtml

Submitted: **2015-02-06**, last modified: **2015-11-09**

We have studied the dynamics of an equal-mass magnetized neutron-star binary
within a resistive magnetohydrodynamic (RMHD) approach in which the highly
conducting stellar interior is matched to an electrovacuum exterior. Because
our analysis is aimed at assessing the modifications introduced by resistive
effects on the dynamics of the binary after the merger and through to collapse,
we have carried out a close comparison with an equivalent simulation performed
within the traditional ideal magnetohydrodynamic approximation. We have found
that there are many similarities between the two evolutions but also one
important difference: the survival time of the hyper massive neutron star
increases in a RMHD simulation. This difference is due to a less efficient
magnetic-braking mechanism in the resistive regime, in which matter can move
across magnetic-field lines, thus reducing the outward transport of angular
momentum. Both the RMHD and the ideal magnetohydrodynamic simulations carried
here have been performed at higher resolutions and with a different grid
structure than those in previous work of ours [L. Rezzolla, B. Giacomazzo, L.
Baiotti, J. Granot, C. Kouveliotou, and M. A. Aloy, Astrophys. J. Letters 732,
L6 (2011)], but confirm the formation of a low-density funnel with an ordered
magnetic field produced by the black hole--torus system. In both regimes the
magnetic field is predominantly toroidal in the highly conducting torus and
predominantly poloidal in the nearly evacuated funnel. Reconnection processes
or neutrino annihilation occurring in the funnel, none of which we model, could
potentially increase the internal energy in the funnel and launch a
relativistic outflow, which, however, is not produced in these simulations.

[71] oai:arXiv.org:1503.09054 [pdf] - 1416744

A coordinate-independent characterization of a black hole shadow

Comments: 12 pages, 8 figures, accepted for publication on MNRAS

Submitted: **2015-03-30**, last modified: **2015-09-29**

A large international effort is under way to assess the presence of a shadow
in the radio emission from the compact source at the centre of our Galaxy,
Sagittarius A$^*$ (Sgr A$^*$). If detected, this shadow would provide the first
direct evidence of the existence of black holes and that Sgr A$^*$ is a
supermassive black hole. In addition, the shape of the shadow could be used to
learn about extreme gravity near the event horizon and to determine which
theory of gravity better describes the observations. The mathematical
description of the shadow has so far used a number of simplifying assumptions
that are unlikely to be met by the real observational data. We here provide a
general formalism to describe the shadow as an arbitrary polar curve expressed
in terms of a Legendre expansion. Our formalism does not presume any knowledge
of the properties of the shadow, e.g. the location of its centre, and offers a
number of routes to characterize the distortions of the curve with respect to
reference circles. These distortions can be implemented in a coordinate
independent manner by different teams analysing the same data. We show that the
new formalism provides an accurate and robust description of noisy
observational data, with smaller error variances when compared to previous
approaches for the measurement of the distortion.

[72] oai:arXiv.org:1506.06620 [pdf] - 1416750

Analytic Bjorken flow in one-dimensional relativistic
magnetohydrodynamics

Comments: 9 pages, 4 figures, corrected typos and symbols, accepted for
publication in Phys.Lett.B

Submitted: **2015-06-22**, last modified: **2015-08-22**

In the initial stage of relativistic heavy-ion collisions, strong magnetic
fields appear due to the large velocity of the colliding charges. The evolution
of these fields appears as a novel and intriguing feature in the
fluid-dynamical description of heavy-ion collisions. In this work, we study
analytically the one-dimensional, longitudinally boost-invariant motion of an
ideal fluid in the presence of a transverse magnetic field. Interestingly, we
find that, in the limit of ideal magnetohydrodynamics, i.e., for infinite
conductivity, and irrespective of the strength of the initial magnetization,
the decay of the fluid energy density $e$ with proper time $\tau$ is the same
as for the time-honored "Bjorken flow" without magnetic field. Furthermore,
when the magnetic field is assumed to decay $\sim \tau^{-a}$, where $a$ is an
arbitrary number, two classes of analytic solutions can be found depending on
whether $a$ is larger or smaller than one. In summary, the analytic solutions
presented here highlight that the Bjorken flow is far more general than
formerly thought. These solutions can serve both to gain insight on the
dynamics of heavy-ion collisions in the presence of strong magnetic fields and
as testbeds for numerical codes.

[73] oai:arXiv.org:1505.00933 [pdf] - 1416746

Recollimation Shocks in Magnetized Relativistic Jets

Comments: 14 pages, 15 figures and 1 table, accepted for publication in ApJ

Submitted: **2015-05-05**, last modified: **2015-06-19**

We have performed two-dimensional special-relativistic magnetohydrodynamic
simulations of non-equilibrium over-pressured relativistic jets in cylindrical
geometry. Multiple stationary recollimation shock and rarefaction structures
are produced along the jet by the nonlinear interaction of shocks and
rarefaction waves excited at the interface between the jet and the surrounding
ambient medium. Although initially the jet is kinematically dominated, we have
considered axial, toroidal and helical magnetic fields to investigate the
effects of different magnetic-field topologies and strengths on the
recollimation structures. We find that an axial field introduces a larger
effective gas-pressure and leads to stronger recollimation shocks and
rarefactions, resulting in larger flow variations. The jet boost grows
quadratically with the initial magnetic field. On the other hand, a toroidal
field leads to weaker recollimation shocks and rarefactions, modifying
significantly the jet structure after the first recollimation rarefaction and
shock. The jet boost decreases systematically. For a helical field, instead,
the behaviour depends on the magnetic pitch, with a phenomenology that ranges
between the one seen for axial and toroidal magnetic fields, respectively. In
general, however, a helical magnetic field yields a more complex shock and
rarefaction substructure close to the inlet that significantly modifies the jet
structure. The differences in shock structure resulting from different field
configurations and strengths may have observable consequences for disturbances
propagating through a stationary recollimation shock.

[74] oai:arXiv.org:1105.5615 [pdf] - 946058

General relativistic radiation hydrodynamics of accretion flows. I:
Bondi-Hoyle accretion

Comments: With updated bibliographyc information

Submitted: **2011-05-27**, last modified: **2015-03-10**

We present a new code for performing general-relativistic
radiation-hydrodynamics simulations of accretion flows onto black holes. The
radiation field is treated in the optically-thick approximation, with the
opacity contributed by Thomson scattering and thermal bremsstrahlung. Our
analysis is concentrated on a detailed numerical investigation of hot
two-dimensional, Bondi-Hoyle accretion flows with various Mach numbers. We find
significant differences with respect to purely hydrodynamical evolutions. In
particular, once the system relaxes to a radiation-pressure dominated regime,
the accretion rates become about two orders of magnitude smaller than in the
purely hydrodynamical case, remaining however super-Eddington as are the
luminosities. Furthermore, when increasing the Mach number of the inflowing
gas, the accretion rates become smaller because of the smaller cross section of
the black hole, but the luminosities increase as a result a stronger emission
in the shocked regions. Overall, our approach provides the first
self-consistent calculation of the Bondi-Hoyle luminosity, most of which is
emitted within r~100 M from the black hole, with typical values L/L_Edd ~ 1-7,
and corresponding energy efficiencies eta_BH ~ 0.09-0.5. The possibility of
computing luminosities self-consistently has also allowed us to compare with
the bremsstrahlung luminosity often used in modelling the electromagnetic
counterparts to supermassive black-hole binaries, to find that in the
optically-thick regime these more crude estimates are about 20 times larger
than our radiation-hydrodynamics results.

[75] oai:arXiv.org:1412.3240 [pdf] - 944817

Spectral properties of the post-merger gravitational-wave signal from
binary neutron stars

Comments: 25 pages, 19 figures, 4 tables; matches version on PRD

Submitted: **2014-12-10**, last modified: **2015-03-06**

Extending previous work by a number of authors, we have recently presented a
new approach in which the detection of gravitational waves from merging neutron
star binaries can be used to determine the equation of state of matter at
nuclear density and hence the structure of neutron stars. In particular, after
performing a large number of numerical-relativity simulations of binaries with
nuclear equations of state, we have found that the post-merger emission is
characterized by two distinct and robust spectral features. While the
high-frequency peak was already shown to be associated with the oscillations of
the hypermassive neutron star produced by the merger and to depend on the
equation of state, we have highlighted that the low-frequency peak is related
to the merger process and to the total compactness of the stars in the binary.
This relation is essentially universal and provides a powerful tool to set
tight constraints on the equation of state. We here provide additional
information on the extensive analysis performed, illustrating the methods used,
the tests considered, as well as the robustness of the results. We also discuss
additional relations that can be deduced when exploring the data and how these
correlate with various properties of the binary. Finally, we present a simple
mechanical toy model that explains the main spectral features of the
post-merger signal and can even reproduce analytically the complex waveforms
emitted right after the merger.

[76] oai:arXiv.org:1502.00551 [pdf] - 929320

High-Order Numerical-Relativity Simulations of Binary Neutron Stars

Comments: Submitted for the ASTRONUM-2014 proceedings. Includes a previously
unpublished high-resolution simulation

Submitted: **2015-02-02**

We report simulations of the inspiral and merger of binary neutron stars
performed with \texttt{WhiskyTHC}, the first of a new generation of numerical
relativity codes employing higher than second-order methods for both the
spacetime and the hydrodynamic evolution. We find that the use of higher-order
schemes improves substantially the quality of the gravitational waveforms
extracted from the simulations when compared to those computed using
traditional second-order schemes. The reduced de-phasing and the faster
convergence rate allow us to estimate the phase evolution of the gravitational
waves emitted, as well as the magnitude of finite-resolution effects, without
the need of phase- or time-alignments or rescalings of the waves, as sometimes
done in other works. Furthermore, by using an additional unpublished simulation
at very high resolution, we confirm the robustness of our high convergence
order of $3.2$.

[77] oai:arXiv.org:1410.8560 [pdf] - 974826

A novel paradigm for short gamma-ray bursts with extended X-ray emission

Comments: 6 pages, 2 figures; clarified several small points; matches version
accepted on ApJ

Submitted: **2014-10-30**, last modified: **2015-01-21**

The merger of a binary of neutron stars provides natural explanations for
many of the features of short gamma-ray bursts (SGRBs), such as the generation
of a hot torus orbiting a rapidly rotating black hole, which can then build a
magnetic jet and provide the energy reservoir to launch a relativistic outflow.
Yet, this scenario has problems explaining the recently discovered long-term
and sustained X-ray emission associated with the afterglows of a subclass of
SGRBs. We propose a new model that explains how an X-ray afterglow can be
sustained by the product of the merger and how the X-ray emission is produced
before the corresponding emission in the gamma-band, although it is observed to
follow it. Overall, our paradigm combines in a novel manner a number of
well-established features of the emission in SGRBs and results from
simulations. Because it involves the propagation of an ultra-relativistic
outflow and its interaction with a confining medium, the paradigm also
highlights a unifying phenomenology between short and long GRBs.

[78] oai:arXiv.org:1402.3258 [pdf] - 1203340

Distortion of neutron stars with a toroidal magnetic field

Comments: 3 pages, 2 figures, to appear in the Proceedings of the 13th Marcel
Grossmann Meeting held in Stockholm 1-7 July, 2012, based on arXiv:1207.4035.
Updated to final published version

Submitted: **2014-02-13**, last modified: **2014-09-12**

Models of rotating relativistic stars with a toroidal magnetic field have
been computed for a sample of eight equations of state of cold dense matter.
Non-rotating models admit important levels of magnetization and quadrupole
distortion accompanied by a seemingly unlimited growth in size. Rotating models
reach the mass-shedding limit at smaller angular velocities than in the
non-magnetized case according to the larger circumferential equatorial radius
induced by the magnetic field. Moreover, they can be classified as
prolate-prolate, oblate-prolate, or oblate-oblate with respect to surface
deformation and quadrupole distortion. Simple expressions for surface and
quadrupole deformation are provided that are valid up to magnetar field
strengths and rapid rotation.

[79] oai:arXiv.org:1408.6526 [pdf] - 864022

The Large Observatory For x-ray Timing

Comments: Proc. SPIE 9144, Space Telescopes and Instrumentation 2014:
Ultraviolet to Gamma Ray, 91442T

Submitted: **2014-08-27**, last modified: **2014-08-29**

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3
Cosmic Vision framework and participated in the final down-selection for a
launch slot in 2022-2024. Thanks to the unprecedented combination of effective
area and spectral resolution of its main instrument, LOFT will study the
behaviour of matter under extreme conditions, such as the strong gravitational
field in the innermost regions of accretion flows close to black holes and
neutron stars, and the supra-nuclear densities in the interior of neutron
stars. The science payload is based on a Large Area Detector (LAD, 10 m 2
effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of
view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1
arcmin source location accuracy, 300 eV spectral resolution). The WFM is
equipped with an on-board system for bright events (e.g. GRB) localization. The
trigger time and position of these events are broadcast to the ground within 30
s from discovery. In this paper we present the status of the mission at the end
of its Phase A study.

[80] oai:arXiv.org:1403.5672 [pdf] - 862593

Constraining the Equation of State of Neutron Stars from Binary Mergers

Comments: 5 pages, 4 figures; matches version on PRL

Submitted: **2014-03-22**, last modified: **2014-08-21**

Determining the equation of state of matter at nuclear density and hence the
structure of neutron stars has been a riddle for decades. We show how the
imminent detection of gravitational waves from merging neutron star binaries
can be used to solve this riddle. Using a large number of accurate
numerical-relativity simulations of binaries with nuclear equations of state,
we find that the postmerger emission is characterized by two distinct and
robust spectral features. While the high-frequency peak has already been
associated with the oscillations of the hypermassive neutron star produced by
the merger and depends on the equation of state, a new correlation emerges
between the low-frequency peak, related to the merger process, and the total
compactness of the stars in the binary. More importantly, such a correlation is
essentially universal, thus providing a powerful tool to set tight constraints
on the equation of state. If the mass of the binary is known from the inspiral
signal, the combined use of the two frequency peaks sets four simultaneous
constraints to be satisfied. Ideally, even a single detection would be
sufficient to select one equation of state over the others. We test our
approach with simulated data and verify it works well for all the equations of
state considered.

[81] oai:arXiv.org:1406.4650 [pdf] - 844928

Future mmVLBI Research with ALMA: A European vision

Comments: Replaced figures 2 and 3: corrected position SRT. Corrected minor
typo in 5.1

Submitted: **2014-06-18**, last modified: **2014-07-01**

Very long baseline interferometry at millimetre/submillimetre wavelengths
(mmVLBI) offers the highest achievable spatial resolution at any wavelength in
astronomy. The anticipated inclusion of ALMA as a phased array into a global
VLBI network will bring unprecedented sensitivity and a transformational leap
in capabilities for mmVLBI. Building on years of pioneering efforts in the US
and Europe the ongoing ALMA Phasing Project (APP), a US-led international
collaboration with MPIfR-led European contributions, is expected to deliver a
beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al.
2013, arXiv:1309.3519).
This report focuses on the future use of mmVLBI by the international users
community from a European viewpoint. Firstly, it highlights the intense science
interest in Europe in future mmVLBI observations as compiled from the responses
to a general call to the European community for future research projects. A
wide range of research is presented that includes, amongst others:
- Imaging the event horizon of the black hole at the centre of the Galaxy
- Testing the theory of General Relativity an/or searching for alternative
theories
- Studying the origin of AGN jets and jet formation
- Cosmological evolution of galaxies and BHs, AGN feedback
- Masers in the Milky Way (in stars and star-forming regions)
- Extragalactic emission lines and astro-chemistry
- Redshifted absorption lines in distant galaxies and study of the ISM and
circumnuclear gas
- Pulsars, neutron stars, X-ray binaries
- Testing cosmology
- Testing fundamental physical constants

[82] oai:arXiv.org:1206.3803 [pdf] - 810962

On the mass radiated by coalescing black-hole binaries

Comments: 9 pages (emulateapj), 4 figures. Matches version in ApJ but includes
slight changes to fig 4 described in Barausse, et al ApJ 786, 76 (2014)
(doi:10.1088/0004-637X/786/1/76), see also
http://www2.iap.fr/users/barausse/erratum_mass_formula.pdf

Submitted: **2012-06-17**, last modified: **2014-04-16**

We derive an analytic phenomenological expression that predicts the final
mass of the black-hole remnant resulting from the merger of a generic binary
system of black holes on quasi-circular orbits. Besides recovering the correct
test-particle limit for extreme mass-ratio binaries, our formula reproduces
well the results of all the numerical-relativity simulations published so far,
both when applied at separations of a few gravitational radii, and when applied
at separations of tens of thousands of gravitational radii. These validations
make our formula a useful tool in a variety of contexts ranging from
gravitational-wave physics to cosmology. As representative examples, we first
illustrate how it can be used to decrease the phase error of the
effective-one-body waveforms during the ringdown phase. Second, we show that,
when combined with the recently computed self-force correction to the binding
energy of nonspinning black-hole binaries, it provides an estimate of the
energy emitted during the merger and ringdown. Finally, we use it to calculate
the energy radiated in gravitational waves by massive black-hole binaries as a
function of redshift, using different models for the seeds of the black-hole
population.

[83] oai:arXiv.org:1401.4544 [pdf] - 1202840

Magnetically driven winds from differentially rotating neutron stars and
X-ray afterglows of short gamma-ray bursts

Comments: 6 pages, 5 figures. Updated to match published article

Submitted: **2014-01-18**, last modified: **2014-04-10**

Besides being among the most promising sources of gravitational waves,
merging neutron star binaries also represent a leading scenario to explain the
phenomenology of short gamma-ray bursts (SGRBs). Recent observations have
revealed a large subclass of SGRBs with roughly constant luminosity in their
X-ray afterglows, lasting $10\!-\!10^4$ s. These features are generally taken
as evidence of a long-lived central engine powered by the magnetic spin-down of
a uniformly rotating, magnetized object. We propose a different scenario in
which the central engine powering the X-ray emission is a differentially
rotating hypermassive neutron star (HMNS) that launches a quasi-isotropic and
baryon-loaded wind driven by the magnetic field, which is built-up through
differential rotation. Our model is supported by long-term, three-dimensional,
general-relativistic, and ideal magnetohydrodynamic simulations, showing that
this isotropic emission is a very robust feature. For a given HMNS, the
presence of a collimated component depends sensitively on the initial magnetic
field geometry, while the stationary electromagnetic luminosity depends only on
the magnetic energy initially stored in the system. We show that our model is
compatible with the observed timescales and luminosities and express the latter
in terms of a simple scaling relation.

[84] oai:arXiv.org:1312.5004 [pdf] - 794535

High-Order Fully General-Relativistic Hydrodynamics: new Approaches and
Tests

Comments: 34 pages, 16 figures. Version accepted on CQG

Submitted: **2013-12-17**, last modified: **2014-02-20**

We present a new approach for achieving high-order convergence in fully
general-relativistic hydrodynamic simulations. The approach is implemented in
WhiskyTHC, a new code that makes use of state-of-the-art numerical schemes and
was key in achieving, for the first time, higher than second-order convergence
in the calculation of the gravitational radiation from inspiraling binary
neutron stars Radice et al. (2013). Here, we give a detailed description of the
algorithms employed and present results obtained for a series of classical
tests involving isolated neutron stars. In addition, using the
gravitational-wave emission from the late inspiral and merger of binary neutron
stars, we make a detailed comparison between the results obtained with the new
code and those obtained when using standard second-order schemes commonly
employed for matter simulations in numerical relativity. We find that even at
moderate resolutions and for binaries with large compactness, the phase
accuracy is improved by a factor 50 or more.

[85] oai:arXiv.org:1207.4035 [pdf] - 783732

Equilibrium models of relativistic stars with a toroidal magnetic field

Comments: Published version with additional correction of remaining typos

Submitted: **2012-07-17**, last modified: **2014-02-13**

We have computed models of rotating relativistic stars with a toroidal
magnetic field and investigated the combined effects of magnetic field and
rotation on the apparent shape (i.e. the surface deformation), which could be
relevant for the electromagnetic emission, and on the internal matter
distribution (i.e. the quadrupole distortion), which could be relevant for the
emission of gravitational waves. Using a sample of eight different cold nuclear
physics equations of state, we have computed models of maximum field strength,
as well as the distortion coefficients for the surface and the quadrupolar
deformations. Surprisingly, we find that non-rotating models admit arbitrary
levels of magnetization, accompanied by a growth of size and quadrupole
distortion to which we could not find a limit. Rotating models, on the other
hand, are subject to a mass-shedding limit at frequencies well below the
corresponding ones for unmagnetized stars. Overall, the space of solutions can
be split into three distinct classes for which the surface deformation and the
quadrupole distortion are either prolate and prolate, oblate and prolate, or
oblate and oblate, respectively. We also derive a simple formula expressing the
relativistic distortion coefficients, which allows one to compute the surface
deformation and the quadrupole distortion up to significant levels of rotation
and magnetization, essentially covering all known magnetars. Such a formula
replaces Newtonian equivalent expressions that overestimate the magnetic
quadrupole distortion by about a factor of 6 and are inadequate for strongly
relativistic objects like neutron stars.

[86] oai:arXiv.org:1309.3885 [pdf] - 843048

On the universality of I-Love-Q relations in magnetized neutron stars

Comments: Published as MNRAS Letters 438, L71-L75 (2014)

Submitted: **2013-09-16**, last modified: **2014-01-28**

Recently, general relations among the quadrupole moment (Q), the moment of
inertia (I), and the tidal deformability (Love number) of a neutron star were
shown to exist. They are nearly independent of the nuclear matter equation of
state and would be of great aid in extracting parameters from observed
gravitational-waves and in testing general relativity. These relations,
however, do not account for strong magnetic fields. We consider this problem by
studying the effect of a strong magnetic field on slowly rotating relativistic
neutron stars and show that, for simple magnetic field configurations that are
purely poloidal or purely toroidal, the relation between Q and I is again
nearly universal. However, different magnetic field geometries lead to
different I-Q relations, and, in the case of a more realistic twisted-torus
magnetic field configuration, the relation depends significantly on the
equation of state, losing its universality. I-Love-Q relations must thus be
used with very great care, since universality is lost for stars with long spin
periods, i.e. P > 10 s, and strong magnetic fields, i.e. B > 10^12 G.

[87] oai:arXiv.org:1307.1409 [pdf] - 809761

Fast radio bursts: the last sign of supramassive neutron stars

Comments: Astronomy & Astrophysics, in press - 6 pages, 1 figure, revised after
referee comments, slightly expanded, corrected numerical error in eq. 3,
normalization of parameter f changed, abstract abridged

Submitted: **2013-07-04**, last modified: **2014-01-21**

Several fast radio bursts have been discovered recently, showing a bright,
highly dispersed millisecond radio pulse. The pulses do not repeat and are not
associated with a known pulsar or gamma-ray burst. The high dispersion suggests
sources at cosmological distances, hence implying an extremely high radio
luminosity, far larger than the power of single pulses from a pulsar. We
suggest that a fast radio burst represents the final signal of a supramassive
rotating neutron star that collapses to a black hole due to magnetic braking.
The neutron star is initially above the critical mass for non-rotating models
and is supported by rapid rotation. As magnetic braking constantly reduces the
spin, the neutron star will suddenly collapse to a black hole several thousand
to million years after its birth. We discuss several formation scenarios for
supramassive neutron stars and estimate the possible observational signatures
{making use of the results of recent numerical general-relativistic
calculations. While the collapse will hide the stellar surface behind an event
horizon, the magnetic-field lines will snap violently. This can turn an almost
ordinary pulsar into a bright radio "blitzar": Accelerated electrons from the
travelling magnetic shock dissipate a significant fraction of the magnetosphere
and produce a massive radio burst that is observable out to z>0.7. Only a few
percent of the neutron stars needs to be supramassive in order to explain the
observed rate. We suggest that fast radio bursts might trace the solitary
formation of stellar mass black holes at high redshifts. These bursts could be
an electromagnetic complement to gravitational-wave emission and reveal a new
formation and evolutionary channel for black holes that are not seen as
gamma-ray bursts. Radio observations of these bursts could trace the
core-collapse supernova rate throughout the universe.

[88] oai:arXiv.org:1312.1862 [pdf] - 881873

Host redshifts from gravitational-wave observations of binary neutron
star mergers

Comments: 10 pages, 4 figures; same as the version before except for
acknowledgments

Submitted: **2013-12-06**, last modified: **2013-12-28**

Inspiralling compact binaries as standard sirens will soon become an
invaluable tool for cosmology when advanced interferometric gravitational-wave
detectors begin their observations in the coming years. However, a degeneracy
in the information carried by gravitational waves between the total rest-frame
mass $M$ and the redshift $z$ of the source implies that neither can be
directly extracted from the signal, but only the combination $M(1+z)$, the
redshifted mass. Recent work has shown that for binary neutron star systems, a
tidal correction to the gravitational-wave phase in the late-inspiral signal
that depends on the rest-frame source mass could be used to break the
mass-redshift degeneracy. We propose here to use the signature encoded in the
post-merger signal to deduce the redshift to the source. This will allow an
accurate extraction of the intrinsic rest-frame mass of the source, in turn
permitting the determination of source redshift and luminosity distance solely
from gravitational-wave observations. This will herald a new era in precision
cosmography and astrophysics. Using numerical simulations of binary neutron
star mergers of very slightly different mass, we model gravitational-wave
signals at different redshifts and use Bayesian parameter estimation to
determine the accuracy with which the redshift can be extracted for a source of
known mass. We find that the Einstein Telescope can determine the source
redshift to $sim 10$--$20%$ at redshifts of $z<0.04$.

[89] oai:arXiv.org:1306.6052 [pdf] - 758965

Beyond second-order convergence in simulations of binary neutron stars
in full general-relativity

Comments: Improved the alignment procedure with PN. Version accepted as a
Letter on MNRAS. 5 pages, 3 figures

Submitted: **2013-06-25**, last modified: **2013-10-23**

Despite the recent rapid progress in numerical relativity, a convergence
order less than the second has so far plagued codes solving the Einstein-Euler
system of equations. We report simulations of the inspiral of binary neutron
stars in quasi-circular orbits computed with a new code employing high-order,
high-resolution shock-capturing, finite-differencing schemes that, for the
first time, go beyond the second-order barrier. In particular, without any
tuning or alignment, we measure a convergence order above three both in the
phase and in the amplitude of the gravitational waves. Because the new code is
able to calculate waveforms with very small phase errors already at modest
resolutions, we are able to obtain accurate estimates of tidal effects in the
inspiral that are essentially free from the large numerical viscosity typical
of lower-order methods, and even for the challenging large compactness and
small-deformability binary considered here. We find a remarkable agreement
between our Richardson-extrapolated waveform and the one from the tidally
corrected post-Newtonian (PN) Taylor-T4 model, with a de-phasing smaller than
0.2 radians during the seven orbits of the inspiral and up to the contact
point. Because our results can be used reliably to assess the validity of the
PN or other approximations at frequencies significantly larger than those
considered so far in the literature, they seem to exclude at these
compactnesses significant tidal amplifications from
next-to-next-to-leading--order terms in the PN expansion.

[90] oai:arXiv.org:1306.4953 [pdf] - 727474

Implementation of a simplified approach to radiative transfer in general
relativity

Comments: 32 pages, 16 figures, updated to version published on PRD. Corrected
an error in table II

Submitted: **2013-06-20**, last modified: **2013-10-02**

We describe in detail the implementation of a simplified approach to
radiative transfer in general relativity by means of the well-known neutrino
leakage scheme (NLS). In particular, we carry out an extensive investigation of
the properties and limitations of the NLS for isolated relativistic stars to a
level of detail that has not been discussed before in a general-relativistic
context. Although the numerous tests considered here are rather idealized, they
provide a well-controlled environment in which to understand the relationship
between the matter dynamics and the neutrino emission, which is important in
order to model the neutrino signals from more complicated scenarios, such as
binary neutron-star mergers. When considering nonrotating hot neutron stars we
confirm earlier results of one-dimensional simulations, but also present novel
results about the equilibrium properties and on how the cooling affects the
stability of these configurations. In our idealized but controlled setup, we
can then show that deviations from the thermal and weak-interaction equilibrium
affect the stability of these models to radial perturbations, leading models
that are stable in the absence of radiative losses, to a gravitational collapse
to a black hole when neutrinos are instead radiated.

[91] oai:arXiv.org:1309.6549 [pdf] - 1179488

Bar-mode instability suppression in magnetized relativistic stars

Comments: 6 pages, 4 figures, to appear on the proceedings of the 4th YRM
(Trieste 2013)

Submitted: **2013-09-25**

We show that magnetic fields stronger than about $10^{15}$ G are able to
suppress the development of the hydrodynamical bar-mode instability in
relativistic stars. The suppression is due to a change in the rest-mass density
and angular velocity profiles due to the formation and to the linear growth of
a toroidal component that rapidly overcomes the original poloidal one, leading
to an amplification of the total magnetic energy. The study is carried out
performing three-dimensional ideal-magnetohydrodynamics simulations in full
general relativity, superimposing to the initial (matter) equilibrium
configurations a purely poloidal magnetic field in the range $10^{14}-10^{16}$
G. When the seed field is a few parts in $10^{15}$ G or above, all the evolved
models show the formation of a low-density envelope surrounding the star. For
much weaker fields, no effect on the matter evolution is observed, while
magnetic fields which are just below the suppression threshold are observed to
slow down the growth-rate of the instability.

[92] oai:arXiv.org:1306.2803 [pdf] - 721537

Twisted-torus configurations with large toroidal magnetic fields in
relativistic stars

Comments: 5 pages, 3 figures, 1 table

Submitted: **2013-06-12**, last modified: **2013-09-19**

Understanding the properties of the internal magnetic field of neutron stars
remains a theoretical challenge. Over the last years, twisted-torus geometries
have been considered both in Newtonian and general-relativistic equilibrium
models, as they represent a potentially good description of neutron star
interiors. All of these works have found an apparent intrinsic limitation to
geometries that are poloidal-field-dominated, with a toroidal-to-poloidal
energy ratio inside the star that are <10%, unless surface currents are
included and magnetic fields are allowed to be discontinuous. This limitation
is in stark contrast with the general expectation that much higher toroidal
fields should be present in the stellar interior and casts doubt about the
stability and hence realism of these configurations. We here discuss how to
overcome this limitation by adopting a new prescription for the azimuthal
currents that leads to magnetized equilibria where the toroidal-to-total
magnetic-field energy ratio can be as high as 90%, thus including geometries
that are toroidal-field-dominated. Moreover, our results show that for a fixed
exterior magnetic-field strength, a higher toroidal-field energy implies a much
higher total magnetic energy stored in the star, with a potentially strong
impact on the expected electromagnetic and gravitational-wave emission from
highly magnetized neutron stars.

[93] oai:arXiv.org:1308.3989 [pdf] - 1173472

Dynamical bar-mode instability in rotating and magnetized relativistic
stars

Comments: 18 pages, 13 figures

Submitted: **2013-08-19**

We present three-dimensional simulations of the dynamical bar-mode
instability in magnetized and differentially rotating stars in full general
relativity. Our focus is on the effects that magnetic fields have on the
dynamics and the onset of the instability. In particular, we perform
ideal-magnetohydrodynamics simulations of neutron stars that are known to be
either stable or unstable against the purely hydrodynamical instability, but to
which a poloidal magnetic field in the range of $10^{14}$--$10^{16}$ G is
superimposed initially. As expected, the differential rotation is responsible
for the shearing of the poloidal field and the consequent linear growth in time
of the toroidal magnetic field. The latter rapidly exceeds in strength the
original poloidal one, leading to a magnetic-field amplification in the the
stars. Weak initial magnetic fields, i.e. $ \lesssim 10^{15}$ G, have
negligible effects on the development of the dynamical bar-mode instability,
simply braking the stellar configuration via magnetic-field shearing, and over
a timescale for which we derived a simple algebraic expression. On the other
hand, strong magnetic fields, i.e. $\gtrsim 10^{16}$ G, can suppress the
instability completely, with the precise threshold being dependent also on the
amount of rotation. As a result, it is unlikely that very highly magnetized
neutron stars can be considered as sources of gravitational waves via the
dynamical bar-mode instability.

[94] oai:arXiv.org:1302.4368 [pdf] - 694536

Magnetorotational instability in relativistic hypermassive neutron stars

Comments: 5 pages, 4 figures. Updated to match published version

Submitted: **2013-02-18**, last modified: **2013-07-17**

A differentially rotating hypermassive neutron star (HMNS) is a metastable
object which can be formed in the merger of neutron-star binaries. The eventual
collapse of the HMNS into a black hole is a key element in generating the
physical conditions expected to accompany the launch of a short gamma-ray
burst. We investigate the influence of magnetic fields on HMNSs by performing
three-dimensional simulations in general-relativistic magnetohydrodynamics. In
particular, we provide direct evidence for the occurrence of the
magnetorotational instability (MRI) in HMNS interiors. For the first time in
simulations of these systems, rapidly-growing and spatially-periodic structures
are observed to form with features like those of the channel flows produced by
the MRI in other systems. Moreover, the growth time and wavelength of the
fastest-growing mode are extracted and compared successfully with analytical
predictions. The MRI emerges as an important mechanism to amplify magnetic
fields over the lifetime of the HMNS, whose collapse to a black hole is
accelerated. The evidence provided here that the MRI can actually develop in
HMNSs could have a profound impact on the outcome of the merger of neutron-star
binaries and on its connection to short gamma-ray bursts.

[95] oai:arXiv.org:1306.4065 [pdf] - 711317

Matter effects on binary neutron star waveforms

Comments: 21 pages, 13 figures

Submitted: **2013-06-17**

Using an extended set of equations of state and a multiple-group
multiple-code collaborative effort to generate waveforms, we improve
numerical-relativity-based data-analysis estimates of the measurability of
matter effects in neutron-star binaries. We vary two parameters of a
parameterized piecewise-polytropic equation of state (EOS) to analyze the
measurability of EOS properties, via a parameter {\Lambda} that characterizes
the quadrupole deformability of an isolated neutron star. We find that, to
within the accuracy of the simulations, the departure of the waveform from
point-particle (or spinless double black-hole binary) inspiral increases
monotonically with {\Lambda}, and changes in the EOS that did not change
{\Lambda} are not measurable. We estimate with two methods the minimal and
expected measurability of {\Lambda} in second- and third- generation
gravitational-wave detectors. The first estimate, using numerical waveforms
alone, shows two EOS which vary in radius by 1.3km are distinguishable in
mergers at 100Mpc. The second estimate relies on the construction of hybrid
waveforms by matching to post-Newtonian inspiral, and estimates that the same
EOS are distinguishable in mergers at 300Mpc. We calculate systematic errors
arising from numerical uncertainties and hybrid construction, and we estimate
the frequency at which such effects would interfere with template-based
searches.

[96] oai:arXiv.org:1209.1634 [pdf] - 651809

A New Spherical Harmonics Scheme for Multi-Dimensional Radiation
Transport I: Static Matter Configurations

Comments: 29 pages, 13 figures. Version matching the one in Journal of
Computational Physics

Submitted: **2012-09-07**, last modified: **2013-04-15**

Recent work by McClarren & Hauck [29] suggests that the filtered spherical
harmonics method represents an efficient, robust, and accurate method for
radiation transport, at least in the two-dimensional (2D) case. We extend their
work to the three-dimensional (3D) case and find that all of the advantages of
the filtering approach identified in 2D are present also in the 3D case. We
reformulate the filter operation in a way that is independent of the timestep
and of the spatial discretization. We also explore different second- and
fourth-order filters and find that the second-order ones yield significantly
better results. Overall, our findings suggest that the filtered spherical
harmonics approach represents a very promising method for 3D radiation
transport calculations.

[97] oai:arXiv.org:1303.6464 [pdf] - 1165538

Three little pieces for computer and relativity

Comments: 42 pages, 11 figures. Plenary talk at "Relativity and Gravitation:
100 Years after Einstein in Prague", June 25 - 29, 2012, Prague, Czech
Republic. To appear in the Proceedings (Edition Open Access). Collects
results appeared in journal articles [72,73, 122-124]

Submitted: **2013-03-26**

Numerical relativity has made big strides over the last decade. A number of
problems that have plagued the field for years have now been mostly solved.
This progress has transformed numerical relativity into a powerful tool to
explore fundamental problems in physics and astrophysics, and I present here
three representative examples. These "three little pieces" reflect a personal
choice and describe work that I am particularly familiar with. However, many
more examples could be made.

[98] oai:arXiv.org:1209.2936 [pdf] - 635499

Universality and intermittency in relativistic turbulent flows of a hot
plasma

Comments: 5 pages, 4 figures. Minor changes to match the version accepted on
ApJL

Submitted: **2012-09-13**, last modified: **2013-03-04**

With the aim of determining the statistical properties of relativistic
turbulence and unveiling novel and non-classical features, we resent the
results of direct numerical simulations of driven turbulence in an
ultrarelativistic hot plasma using high-order numerical schemes. We study the
statistical properties of flows with average Mach number ranging from $\sim
0.4$ to $\sim 1.7$ and with average Lorentz factors up to $\sim 1.7$. We find
that flow quantities, such as the energy density or the local Lorentz factor,
show large spatial variance even in the subsonic case as compressibility is
enhanced by relativistic effects. The velocity field is highly intermittent,
but its power-spectrum is found to be in good agreement with the predictions of
the classical theory of Kolmogorov. Overall, our results indicate that
relativistic effects are able to significantly enhance the intermittency of the
flow and affect the high-order statistics of the velocity field, while leaving
unchanged the low-order statistics, which instead appear to be universal and in
good agreement with the classical Kolmogorov theory. To the best of our
knowledge, these are the most accurate simulations of driven relativistic
turbulence to date.

[99] oai:arXiv.org:1210.8152 [pdf] - 603598

Compact Binary Progenitors of Short Gamma-Ray Bursts

Comments: 6 pages, 3 figures. Extended the discussion on jet efficiency.
Accepted for publication on The Astrophysical Journal Letters

Submitted: **2012-10-30**, last modified: **2012-12-11**

In recent years, detailed observations and accurate numerical simulations
have provided support to the idea that mergers of compact binaries containing
either two neutron stars (NSs) or an NS and a black hole (BH) may constitute
the central engine of short gamma-ray bursts (SGRBs). The merger of such
compact binaries is expected to lead to the production of a spinning BH
surrounded by an accreting torus. Several mechanisms can extract energy from
this system and power the SGRBs. Here we connect observations and numerical
simulations of compact binary mergers, and use the current sample of SGRBs with
measured energies to constrain the mass of their powering tori. By comparing
the masses of the tori with the results of fully general-relativistic
simulations, we are able to infer the properties of the binary progenitors
which yield SGRBs. By assuming a constant efficiency in converting torus mass
into jet energy, epsilon_{jet}=10%, we find that most of the tori have masses
smaller than 0.01M_{sun}, favoring "high-mass" binary NSs mergers, i.e.,
binaries with total masses >~1.5 the maximum mass of an isolated NS. This has
important consequences for the gravitational-wave signals that may be detected
in association with SGRBs, since "high-mass" systems do not form a long-lived
hypermassive NS after the merger. While NS-BH systems cannot be excluded to be
the engine of at least some of the SGRBs, the BH would need to have an initial
spin of ~0.9, or higher.

[100] oai:arXiv.org:1209.6138 [pdf] - 602254

Black-hole production from ultrarelativistic collisions

Comments: Accepted as FTC on CQG; no discussion of LHC (seen as too speculative
but available in v1); expanded considerations on hoop conjecture

Submitted: **2012-09-27**, last modified: **2012-11-29**

Determining the conditions under which a black hole can be produced is a
long-standing and fundamental problem in general relativity. We use numerical
simulations of colliding selfgravitating fluid objects to study the conditions
of black-hole formation when the objects are boosted to ultrarelativistic
speeds. Expanding on previous work, we show that the collision is characterized
by a type-I critical behaviour, with a black hole being produced for masses
above a critical value, M_c, and a partially bound object for masses below the
critical one. More importantly, we show for the first time that the critical
mass varies with the initial effective Lorentz factor <\gamma> following a
simple scaling of the type M_c ~ K <\gamma>^{-1.0}, thus indicating that a
black hole of infinitesimal mass is produced in the limit of a diverging
Lorentz factor. Furthermore, because a scaling is present also in terms of the
initial stellar compactness, we provide a condition for black-hole formation in
the spirit of the hoop conjecture.

[101] oai:arXiv.org:1206.6604 [pdf] - 1124459

Poloidal-Field Instability in Magnetized Relativistic Stars

Comments: 13 pages, 11 figures, updated to match the published version

Submitted: **2012-06-28**, last modified: **2012-11-01**

We investigate the instability of purely poloidal magnetic fields in
nonrotating neutron stars by means of three-dimensional general-relativistic
magnetohydrodynamics simulations, extending the work presented in Ciolfi et al.
(2011). Our aim is to draw a clear picture of the dynamics associated with the
instability and to study the final configuration reached by the system, thus
obtaining indications on possible equilibria in a magnetized neutron star.
Furthermore, since the internal rearrangement of magnetic fields is a highly
dynamical process, which has been suggested to be behind magnetar giant flares,
our simulations can provide a realistic estimate of the electromagnetic and
gravitational-wave emission which should accompany the flare event. Our main
findings are the following: (i) the initial development of the instability
meets all the expectations of perturbative studies in terms of the location of
the seed of the instability, the timescale for its growth and the generation of
a toroidal component; (ii) in the subsequent nonlinear reorganization of the
system, ~90% of magnetic energy is lost in few Alfven timescales mainly through
electromagnetic emission, and further decreases on a much longer timescale;
(iii) all stellar models tend to achieve a significant amount of magnetic
helicity and the equipartition of energy between poloidal and toroidal magnetic
fields, and evolve to a new configuration which does not show a subsequent
instability on dynamical or Alfven timescales; (iv) the electromagnetic
emission matches the duration of the initial burst in luminosity observed in
giant flares, giving support to the internal rearrangement scenario; (v) only a
small fraction of the energy released during the process is converted into
f-mode oscillations and in the consequent gravitational-wave emission, thus
resulting in very low chances of detecting this signal with present and..

[102] oai:arXiv.org:1210.7081 [pdf] - 584063

Universality and intermittency in relativistic turbulent flows of a hot
gas

Comments: Talk given at the ASTRONUM2012 conference on the 25th of June 2012

Submitted: **2012-10-26**

With the aim of determining the statistical properties of relativistic
turbulence and unveiling novel and non-classical features, we present the
results of direct numerical simulations of driven turbulence in an
ultrarelativistic hot plasma using high-order numerical schemes. We study the
statistical properties of flows with average Mach number ranging from $\sim
0.4$ to $\sim 1.7$ and with average Lorentz factors up to $\sim 1.7$. We find
that flow quantities, such as the energy density or the local Lorentz factor,
show large spatial variance even in the subsonic case as compressibility is
enhanced by relativistic effects. The velocity field is highly intermittent,
but its power-spectrum is found to be in good agreement with the predictions of
the classical theory of Kolmogorov.

[103] oai:arXiv.org:1206.6502 [pdf] - 579331

THC: a new high-order finite-difference high-resolution shock-capturing
code for special-relativistic hydrodynamics

Comments: Updated to match the published version

Submitted: **2012-06-27**, last modified: **2012-10-22**

We present THC: a new high-order flux-vector-splitting code for Newtonian and
special-relativistic hydrodynamics designed for direct numerical simulations of
turbulent flows. Our code implements a variety of different reconstruction
algorithms, such as the popular weighted essentially non oscillatory and
monotonicity-preserving schemes, or the more specialised bandwidth-optimised
WENO scheme that has been specifically designed for the study of compressible
turbulence. We show the first systematic comparison of these schemes in
Newtonian physics as well as for special-relativistic flows. In particular we
will present the results obtained in simulations of grid-aligned and oblique
shock waves and nonlinear, large-amplitude, smooth adiabatic waves. We will
also discuss the results obtained in classical benchmarks such as the
double-Mach shock reflection test in Newtonian physics or the linear and
nonlinear development of the relativistic Kelvin-Helmholtz instability in two
and three dimensions. Finally, we study the turbulent flow induced by the
Kelvin-Helmholtz instability and we show that our code is able to obtain
well-converged velocity spectra, from which we benchmark the effective
resolution of the different schemes.

[104] oai:arXiv.org:1108.1423 [pdf] - 520241

Scientific Potential of Einstein Telescope

Comments: Conforms to conference proceedings, several author names corrected

Submitted: **2011-08-05**, last modified: **2012-06-01**

Einstein gravitational-wave Telescope (ET) is a design study funded by the
European Commission to explore the technological challenges of and scientific
benefits from building a third generation gravitational wave detector. The
three-year study, which concluded earlier this year, has formulated the
conceptual design of an observatory that can support the implementation of new
technology for the next two to three decades. The goal of this talk is to
introduce the audience to the overall aims and objectives of the project and to
enumerate ET's potential to influence our understanding of fundamental physics,
astrophysics and cosmology.

[105] oai:arXiv.org:1109.1177 [pdf] - 1083861

On the detectability of dual jets from binary black holes

Comments: 4 pages, 3 figures

Submitted: **2011-09-06**, last modified: **2012-04-03**

We revisit the suggestion that dual jets can be produced during the inspiral
and merger of supermassive black holes when these are immersed in a force-free
plasma threaded by a uniform magnetic field. By performing independent
calculations and by computing the electromagnetic emission in a way which is
consistent with estimates using the Poynting flux, we show that a dual-jet
structure is present but energetically subdominant with respect to a
non-collimated and predominantly quadrupolar emission, which is similar to the
one computed when the binary is in electrovacuum. While our findings set
serious restrictions on the detectability of dual jets from coalescing
binaries, they also increase the chances of detecting an EM counterpart from
these systems.

[106] oai:arXiv.org:1105.0122 [pdf] - 966780

Collapse of differentially rotating neutron stars and cosmic censorship

Comments: 15 pages, 15 figures. Minor changes to the text and to the
references. In press on Phys. Rev. D

Submitted: **2011-04-30**, last modified: **2011-07-07**

We present new results on the dynamics and gravitational-wave emission from
the collapse of differentially rotating neutron stars. We have considered a
number of polytropic stellar models having different values of the
dimensionless angular momentum J/M^2, where J and M are the asymptotic angular
momentum and mass of the star, respectively. For neutron stars with J/M^2<1,
i.e., "sub-Kerr" models, we were able to find models that are dynamically
unstable and that collapse promptly to a rotating black hole. Both the dynamics
of the collapse and the consequent emission of gravitational waves resemble the
one seen for uniformly rotating stars, although with an overall decrease in the
efficiency of gravitational-wave emission. For stellar models with J/M^2>1,
i.e., "supra-Kerr" models, on the other hand, we were not able to find models
that are dynamically unstable and all of the computed supra-Kerr models were
found to be far from the stability threshold. For these models a gravitational
collapse is possible only after a very severe and artificial reduction of the
pressure, which then leads to a torus developing nonaxisymmetric instabilities
and eventually contracting to a stable axisymmetric stellar configuration.
While this does not exclude the possibility that a naked singularity can be
produced by the collapse of a differentially rotating star, it also suggests
that cosmic censorship is not violated and that generic conditions for a
supra-Kerr progenitor do not lead to a naked singularity.

[107] oai:arXiv.org:1105.3069 [pdf] - 400355

A quasi-radial stability criterion for rotating relativistic stars

Comments: 5 pages, 4 figures; accepted to MNRAS Letters

Submitted: **2011-05-16**, last modified: **2011-06-21**

The stability properties of relativistic stars against gravitational collapse
to black hole is a classical problem in general relativity. A sufficient
criterion for secular instability was established by Friedman, Ipser and Sorkin
(1988), who proved that a sequence of uniformly rotating barotropic stars is
secularly unstable on one side of a turning point and then argued that a
stronger result should hold: that the sequence should be stable on the opposite
side, with the turning point marking the onset of secular instability. We show
here that this expectation is not met. By computing in full general relativity
the $F$-mode frequency for a large number of rotating stars, we show that the
neutral-stability point, i.e., where the frequency becomes zero, differs from
the turning point for rotating stars. Using numerical simulations we validate
that the new criterion can be used to assess the dynamical stability of
relativistic rotating stars.

[108] oai:arXiv.org:1103.2426 [pdf] - 395275

Discontinuous Galerkin methods for general-relativistic hydrodynamics:
formulation and application to spherically symmetric spacetimes

Comments: 24 pages, 19 figures. Small changes; matches version to appear in PRD

Submitted: **2011-03-12**, last modified: **2011-06-14**

We have developed the formalism necessary to employ the
discontinuous-Galerkin approach in general-relativistic hydrodynamics. The
formalism is firstly presented in a general 4-dimensional setting and then
specialized to the case of spherical symmetry within a 3+1 splitting of
spacetime. As a direct application, we have constructed a one-dimensional code,
EDGES, which has been used to asses the viability of these methods via a series
of tests involving highly relativistic flows in strong gravity. Our results
show that discontinuous Galerkin methods are able not only to handle strong
relativistic shock waves but, at the same time, to attain very high orders of
accuracy and exponential convergence rates in smooth regions of the flow. Given
these promising prospects and their affinity with a pseudospectral solution of
the Einstein equations, discontinuous Galerkin methods could represent a new
paradigm for the accurate numerical modelling in relativistic astrophysics.

[109] oai:arXiv.org:1105.3971 [pdf] - 966836

Instability-driven evolution of poloidal magnetic fields in relativistic
stars

Comments: 5 pages, 4 figures; small changes; accepted for publication in ApJL

Submitted: **2011-05-19**, last modified: **2011-06-12**

The problem of the stability of magnetic fields in stars has a long history
and has been investigated in detail in perturbation theory. Here we consider
the nonlinear evolution of a non-rotating neutron star with a purely poloidal
magnetic field, in general relativity. We find that an instability develops in
the region of the closed magnetic field lines and over an Alfven timescale, as
predicted by perturbation theory. After the initial unstable growth, our
evolutions show that a toroidal magnetic field component is generated, which
increases until it is locally comparable in strength with the poloidal one. On
longer timescales the system relaxes to a new non-axisymmetric configuration
with a reorganization of the stellar structure and large-amplitude
oscillations, mostly in the fundamental mode. We discuss the energies involved
in the instability and the impact they may have on the phenomenology of
magnetar flares and on their detectability through gravitational-wave emission.

[110] oai:arXiv.org:1010.1739 [pdf] - 955482

On the development of QPOs in Bondi-Hoyle accretion flows

Comments: 11 pages, 7 figures

Submitted: **2010-10-08**, last modified: **2011-04-11**

The numerical investigation of Bondi-Hoyle accretion onto a moving black hole
has a long history, both in Newtonian and in general-relativistic physics. By
performing new two-dimensional and general-relativistic simulations onto a
rotating black hole, we point out a novel feature, namely, that quasi-periodic
oscillations (QPOs) are naturally produced in the shock cone that develops in
the downstream part of the flow. Because the shock cone in the downstream part
of the flow acts as a cavity trapping pressure perturbations, modes with
frequencies in the integer ratios 2:1 and 3:1 are easily produced. The
frequencies of these modes depend on the black-hole spin and on the properties
of the flow, and scale linearly with the inverse of the black-hole mass. Our
results may be relevant for explaining the detection of QPOs in Sagittarius A*,
once such detection is confirmed by further observations. Finally, we report on
the development of the flip-flop instability, which can affect the shock cone
under suitable conditions; such an instability has been discussed before in
Newtonian simulations but was never found in a relativistic regime.

[111] oai:arXiv.org:1101.4298 [pdf] - 1051614

The missing link: Merging neutron stars naturally produce jet-like
structures and can power short Gamma-Ray Bursts

Comments: minor changes; matches published version in ApJL; for high-resolution
images and animations see
http://numrel.aei.mpg.de/images/relativistic-binary-neutron-star-inspirals

Submitted: **2011-01-22**, last modified: **2011-03-30**

Short Gamma-Ray Bursts (SGRBs) are among the most luminous explosions in the
universe, releasing in less than one second the energy emitted by our Galaxy
over one year. Despite decades of observations, the nature of their
"central-engine" remains unknown. Considering a binary of magnetized neutron
stars and solving Einstein equations, we show that their merger results in a
rapidly spinning black hole surrounded by a hot and highly magnetized torus.
Lasting over 35 ms and much longer than previous simulations, our study reveals
that magnetohydrodynamical instabilities amplify an initially turbulent
magnetic field of ~ 10^{12} G to produce an ordered poloidal field of ~ 10^{15}
G along the black-hole spin-axis, within a half-opening angle of ~ 30 deg,
which may naturally launch a relativistic jet. The broad consistency of our
ab-initio calculations with SGRB observations shows that the merger of
magnetized neutron stars can provide the basic physical conditions for the
central-engine of SGRBs.

[112] oai:arXiv.org:1103.3526 [pdf] - 437696

Will black hole-neutron star binary inspirals tell us about the neutron
star equation of state?

Comments: 4 pages, 4 figures

Submitted: **2011-03-17**

The strong tidal forces that arise during the last stages of the life of a
black hole-neutron star binary may severely distort, and possibly disrupt, the
star. Both phenomena will imprint signatures about the stellar structure in the
emitted gravitational radiation. The information from the disruption, however,
is confined to very high frequencies, where detectors are not very sensitive.
We thus assess whether the lack of tidal distortion corrections in
data-analysis pipelines will affect the detection of the inspiral part of the
signal and whether these may yield information on the equation of state of
matter at nuclear densities. Using recent post-Newtonian expressions and
realistic equations of state to model these scenarios, we find that
point-particle templates are sufficient for the detection of black hole-neutron
star inspiralling binaries, with a loss of signals below 1% for both second and
third-generation detectors. Such detections may be able to constrain
particularly stiff equations of state, but will be unable to reveal the
presence of a neutron star with a soft equation of state.

[113] oai:arXiv.org:1009.2468 [pdf] - 319850

Accurate evolutions of inspiralling and magnetized neutron-stars:
equal-mass binaries

Comments: 18 pages, 10 figures. Added two new figures (figures 1 and 7). Small
modifications to the text to match the version published on Phys. Rev. D

Submitted: **2010-09-13**, last modified: **2011-02-10**

By performing new, long and numerically accurate general-relativistic
simulations of magnetized, equal-mass neutron-star binaries, we investigate the
role that realistic magnetic fields may have in the evolution of these systems.
In particular, we study the evolution of the magnetic fields and show that they
can influence the survival of the hypermassive-neutron star produced at the
merger by accelerating its collapse to a black hole. We also provide evidence
that even if purely poloidal initially, the magnetic fields produced in the
tori surrounding the black hole have toroidal and poloidal components of
equivalent strength. When estimating the possibility that magnetic fields could
have an impact on the gravitational-wave signals emitted by these systems
either during the inspiral or after the merger we conclude that for realistic
magnetic-field strengths B<~1e12 G such effects could be detected, but only
marginally, by detectors such as advanced LIGO or advanced Virgo. However,
magnetically induced modifications could become detectable in the case of
small-mass binaries and with the development of gravitational-wave detectors,
such as the Einstein Telescope, with much higher sensitivities at frequencies
larger than ~2 kHz.

[114] oai:arXiv.org:1007.4160 [pdf] - 1033870

Black hole-neutron star mergers and short GRBs: a relativistic toy model
to estimate the mass of the torus

Comments: Added new figure and new table to confirm agreement with simulations;
matches version accepted for publication in ApJ

Submitted: **2010-07-23**, last modified: **2010-11-05**

The merger of a binary system composed of a black hole and a neutron star may
leave behind a torus of hot, dense matter orbiting around the black hole. While
numerical-relativity simulations are necessary to simulate this process
accurately, they are also computationally expensive and unable at present to
cover the large space of possible parameters, which include the relative mass
ratio, the stellar compactness, and the black hole spin. To mitigate this and
provide a first reasonable coverage of the space of parameters, we have
developed a method for estimating the mass of the remnant torus from black
hole-neutron star mergers. The toy model makes use of an improved relativistic
affine model to describe the tidal deformations of an extended tri-axial
ellipsoid orbiting around a Kerr black hole and measures the mass of the
remnant torus by considering which of the fluid particles composing the star
are on bound orbits at the time of the tidal disruption. We tune the toy model
by using the results of fully general-relativistic simulations obtaining
relative precisions of a few percent and use it to extensively investigate the
space of parameters. In this way we find that the torus mass is largest for
systems with highly spinning black holes, small stellar compactnesses, and
large mass ratios. As an example, tori as massive as ~1.33 solar masses can be
produced for a very extended star with compactness of ~0.1 inspiralling around
a black hole with dimensionless spin equal to 0.85 and mass ratio of about 0.3.
However, for a more astrophysically reasonable mass ratio of ~0.14 and a
canonical value of the stellar compactness of ~0.145, the toy model sets a
considerably smaller upper limit to the torus mass of less than ~0.34 solar
masses.

[115] oai:arXiv.org:1009.0019 [pdf] - 296842

Constraining scalar fields with stellar kinematics and collisional dark
matter

Comments: 23 pages, 8 figures; accepted for publication by JCAP after minor
changes

Submitted: **2010-08-31**, last modified: **2010-10-02**

The existence and detection of scalar fields could provide solutions to
long-standing puzzles about the nature of dark matter, the dark compact objects
at the centre of most galaxies, and other phenomena. Yet, self-interacting
scalar fields are very poorly constrained by astronomical observations, leading
to great uncertainties in estimates of the mass $m_\phi$ and the
self-interacting coupling constant $\lambda$ of these fields. To counter this,
we have systematically employed available astronomical observations to develop
new constraints, considerably restricting this parameter space. In particular,
by exploiting precise observations of stellar dynamics at the centre of our
Galaxy and assuming that these dynamics can be explained by a single boson
star, we determine an upper limit for the boson star compactness and impose
significant limits on the values of the properties of possible scalar fields.
Requiring the scalar field particle to follow a collisional dark matter model
further narrows these constraints. Most importantly, we find that if a scalar
dark matter particle does exist, then it cannot account for both the
dark-matter halos and the existence of dark compact objects in galactic nuclei

[116] oai:arXiv.org:1002.4185 [pdf] - 281135

EM counterparts of recoiling black holes: general relativistic
simulations of non-Keplerian discs

Comments: 17 pages, 11 figures, accepted by A&A, movies available at
http://numrel.aei.mpg.de/Visualisations/Archive/BinaryBlackHoles/EMCounterparts/EMCounterparts.html

Submitted: **2010-02-22**, last modified: **2010-09-12**

We investigate the dynamics of a circumbinary disc that responds to the loss
of mass and to the recoil velocity of the black hole produced by the merger of
a binary system of supermassive black holes. We perform the first
two-dimensional general relativistic hydrodynamics simulations of
\textit{extended} non-Keplerian discs and employ a new technique to construct a
"shock detector", thus determining the precise location of the shocks produced
in the accreting disc by the recoiling black hole. In this way we can study how
the properties of the system, such as the spin, mass and recoil velocity of the
black hole, affect the mass accretion rate and are imprinted on the
electromagnetic emission from these sources. We argue that the estimates of the
bremsstrahlung luminosity computed without properly taking into account the
radiation transfer yield cooling times that are unrealistically short. At the
same time we show, through an approximation based on the relativistic
isothermal evolution, that the luminosity produced can reach a peak value above
$L \simeq 10^{43} \ {\rm erg/s} $ at about $\sim 30\,{\rm d}$ after the merger
of a binary with total mass $M\simeq 10^6 M_\odot$ and persist for several days
at values which are a factor of a few smaller. If confirmed by more
sophisticated calculations such a signal could indeed lead to an
electromagnetic counterpart of the merger of binary black-hole system.

[117] oai:arXiv.org:1003.0873 [pdf] - 195351

Understanding the "anti-kick" in the merger of binary black holes

Comments: 4 pages; small changes to match the published version

Submitted: **2010-03-03**, last modified: **2010-06-28**

The generation of a large recoil velocity from the inspiral and merger of
binary black holes represents one of the most exciting results of
numerical-relativity calculations. While many aspects of this process have been
investigated and explained, the "antikick", namely the sudden deceleration
after the merger, has not yet found a simple explanation. We show that the
antikick can be understood in terms of the radiation from a deformed black hole
where the anisotropic curvature distribution on the horizon correlates with the
direction and intensity of the recoil. Our analysis is focussed on
Robinson-Trautman spacetimes and allows us to measure both the energies and
momenta radiated in a gauge-invariant manner. At the same time, this simpler
setup provides the qualitative and quantitative features of merging black
holes, opening the way to a deeper understanding of the nonlinear dynamics of
black-hole spacetimes.

[118] oai:arXiv.org:1001.5281 [pdf] - 174337

On the Shear Instability in Relativistic Neutron Stars

Comments: 20 pages; accepted to the Classical and Quantum Gravity special issue
for MICRA2009

Submitted: **2010-01-28**, last modified: **2010-05-31**

We present new results on instabilities in rapidly and differentially
rotating neutron stars. We model the stars in full general relativity and
describe the stellar matter adopting a cold realistic equation of state based
on the unified SLy prescription. We provide evidence that rapidly and
differentially rotating stars that are below the expected threshold for the
dynamical bar-mode instability, beta_c = T/|W| ~ 0.25, do nevertheless develop
a shear instability on a dynamical timescale and for a wide range of values of
beta. This class of instability, which has so far been found only for small
values of beta and with very small growth rates, is therefore more generic than
previously found and potentially more effective in producing strong sources of
gravitational waves. Overall, our findings support the phenomenological
predictions made by Watts, Andersson and Jones on the nature of the low-T/|W|.

[119] oai:arXiv.org:1001.3074 [pdf] - 1024677

Accurate evolutions of unequal-mass neutron-star binaries: properties of
the torus and short GRB engines

Comments: 35 pages; small changes to match the published version

Submitted: **2010-01-15**, last modified: **2010-05-13**

We present new results from accurate and fully general-relativistic
simulations of the coalescence of unmagnetized binary neutron stars with
various mass ratios. The evolution of the stars is followed through the
inspiral phase, the merger and prompt collapse to a black hole, up until the
appearance of a thick accretion disk, which is studied as it enters and remains
in a regime of quasi-steady accretion. Although a simple ideal-fluid equation
of state with \Gamma=2 is used, this work presents a systematic study within a
fully general relativistic framework of the properties of the resulting
black-hole--torus system produced by the merger of unequal-mass binaries. More
specifically, we show that: (1) The mass of the torus increases considerably
with the mass asymmetry and equal-mass binaries do not produce significant tori
if they have a total baryonic mass M_tot >~ 3.7 M_sun; (2) Tori with masses
M_tor ~ 0.2 M_sun are measured for binaries with M_tot ~ 3.4 M_sun and mass
ratios q ~ 0.75-0.85; (3) The mass of the torus can be estimated by the simple
expression M_tor(q, M_tot) = [c_1 (1-q) + c_2](M_max-M_tot), involving the
maximum mass for the binaries and coefficients constrained from the
simulations, and suggesting that the tori can have masses as large as M_tor ~
0.35 M_sun for M_tot ~ 2.8 M_sun and q ~ 0.75-0.85; (4) Using a novel technique
to analyze the evolution of the tori we find no evidence for the onset of
non-axisymmetric instabilities and that very little, if any, of their mass is
unbound; (5) Finally, for all the binaries considered we compute the complete
gravitational waveforms and the recoils imparted to the black holes, discussing
the prospects of detection of these sources for a number of present and future
detectors.

[120] oai:arXiv.org:0907.0462 [pdf] - 25933

Gravitational-wave detectability of equal-mass black-hole binaries with
aligned spins

Comments: 18 pages, 11 figures, matches published version

Submitted: **2009-07-02**, last modified: **2010-01-08**

Binary black-hole systems with spins aligned or anti-aligned to the orbital
angular momentum provide the natural ground to start detailed studies of the
influence of strong-field spin effects on gravitational wave observations of
coalescing binaries. Furthermore, such systems may be the preferred end-state
of the inspiral of generic supermassive binary black-hole systems. In view of
this, we have computed the inspiral and merger of a large set of binary systems
of equal-mass black holes with spins parallel to the orbital angular momentum
but otherwise arbitrary. Our attention is particularly focused on the
gravitational-wave emission so as to quantify how much spin effects contribute
to the signal-to-noise ratio, to the horizon distances, and to the relative
event rates for the representative ranges in masses and detectors. As expected,
the signal-to-noise ratio increases with the projection of the total black hole
spin in the direction of the orbital momentum. We find that equal-spin binaries
with maximum spin aligned with the orbital angular momentum are more than
"three times as loud" as the corresponding binaries with anti-aligned spins,
thus corresponding to event rates up to 30 times larger. We also consider the
waveform mismatch between the different spinning configurations and find that,
within our numerical accuracy, binaries with opposite spins S_1=-S_2 cannot be
distinguished whereas binaries with spin S_1=S_2 have clearly distinct
gravitational-wave emissions. Finally, we derive a simple expression for the
energy radiated in gravitational waves and find that the binaries always have
efficiencies E_rad/M > 3.6%, which can become as large as E_rad/M = 10% for
maximally spinning binaries with spins aligned with the orbital angular
momentum.

[121] oai:arXiv.org:0912.2330 [pdf] - 142291

Vacuum Electromagnetic Counterparts of Binary Black-Hole Mergers

Comments:

Submitted: **2009-12-11**

As one step towards a systematic modeling of the electromagnetic (EM)
emission from an inspiralling black hole binary we consider a simple scenario
in which the binary moves in a uniform magnetic field anchored to a distant
circumbinary disc. We study this system by solving the Einstein-Maxwell
equations in which the EM fields are chosen with astrophysically consistent
strengths. We consider binaries with spins aligned or anti-aligned with the
orbital angular momentum and study the dependence of gravitational and EM
signals with these spin configurations. Overall we find that the EM radiation
in the lowest l=2, m=2 multipole accurately reflects the gravitational one,
with identical phase evolutions and amplitudes that differ only by a scaling
factor. We also compute the efficiency of the energy emission in EM waves and
find that it is given by E^rad_EM/M ~ 10^-15 (M/10^8 M_Sun)^2 (B/10^4 G)^2,
hence 13 orders of magnitude smaller than the gravitational energy for
realistic magnetic fields. The corresponding luminosity is much smaller than
the accretion luminosity if the system is accreting at near the Eddington rate.
Most importantly, this EM emission is at frequencies of 10^-4 (10^8 M_Sun/M)
Hz, well outside those accessible to astronomical radio observations. As a
result, it is unlikely that the EM emission discussed here can be detected
directly and simultaneously with the gravitational-wave one. However, indirect
processes, driven by changes in the EM fields behavior could yield observable
events. In particular if the accretion rate of the circumbinary disc is small
and sufficiently stable over the timescale of the final inspiral, then the EM
emission may be observable indirectly as it will alter the accretion rate
through the magnetic torques exerted by the distorted magnetic field lines.

[122] oai:arXiv.org:0912.0384 [pdf] - 297258

Gravitational waves from neutron stars: Promises and challenges

Comments: review article for Einstein Telescope special issue in Gen Rel Grav,
39 pages, 7 figures

Submitted: **2009-12-02**

We discuss different ways that neutron stars can generate gravitational
waves, describe recent improvements in modelling the relevant scenarios in the
context of improving detector sensitivity, and show how observations are
beginning to test our understanding of fundamental physics. The main purpose of
the discussion is to establish promising science goals for third-generation
ground-based detectors, like the Einstein Telescope, and identify the various
challenges that need to be met if we want to use gravitational-wave data to
probe neutron star physics.

[123] oai:arXiv.org:0910.2703 [pdf] - 142242

Axisymmetric General Relativistic Simulations of the Accretion-Induced
Collapse of White Dwarfs

Comments: Submitted to PRD, 31 pages, 23 figures

Submitted: **2009-10-14**

(Abridged.) The accretion-induced collapse (AIC) of a white dwarf (WD) may
lead to the formation of a protoneutron star and a collapse-driven supernova
explosion. This process represents a path alternative to thermonuclear
disruption of accreting white dwarfs in Type Ia supernovae. Neutrino and
gravitational-wave (GW) observations may provide crucial information necessary
to reveal a potential AIC. Motivated by the need for systematic predictions of
the GW signature of AIC, we present results from an extensive set of
general-relativistic AIC simulations using a microphysical finite-temperature
equation of state and an approximate treatment of deleptonization during
collapse. Investigating a set of 114 progenitor models in rotational
equilibrium, with a wide range of rotational configurations, temperatures and
central densities, we extend previous Newtonian studies and find that the GW
signal has a generic shape akin to what is known as a "Type III" signal in the
literature. We discuss the detectability of the emitted GWs, showing that the
signal-to-noise ratio for current or next-generation interferometer detectors
could be high enough to detect such events in our Galaxy. Some of our AIC
models form massive quasi-Keplerian accretion disks after bounce. In rapidly
differentially rotating models, the disk mass can be as large as ~0.8-Msun.
Slowly and/or uniformly rotating models produce much smaller disks. Finally, we
find that the postbounce cores of rapidly spinning white dwarfs can reach
sufficiently rapid rotation to develop a nonaxisymmetric rotational
instability.

[124] oai:arXiv.org:0904.2577 [pdf] - 363591

Predicting the direction of the final spin from the coalescence of two
black holes

Comments: 5 pages, 3 figures. Significant changes to text and figures to match
version accepted for publication in ApJ Letters

Submitted: **2009-04-16**, last modified: **2009-09-14**

Knowledge of the spin of the black hole resulting from the merger of a
generic black-hole binary is of great importance for studying the cosmological
evolution of supermassive black holes. Several attempts have been made to model
the spin via simple expressions exploiting the results of numerical-relativity
simulations. While these expressions are in reasonable agreement with the
simulations, they neglect the precession of the binary's orbital plane, and
cannot therefore be applied directly -- i.e., without evolving the system to
small separations using post-Newtonian theory -- to binaries with separations
larger than a few hundred gravitational radii. While not a problem in
principle, this may be impractical if the formulas are employed in cosmological
merger-trees or N-body simulations, which provide the spins and angular
momentum of the two black holes when their separation is of hundreds or
thousands of gravitational radii. The formula that we propose is instead built
on improved assumptions and gives, for any separation, a very accurate
prediction both for the norm of the final spin and for its direction. By
comparing with the numerical data, we also show that the final-spin direction
is very accurately aligned with the binary's total angular momentum at large
separation. Hence, observations of the final-spin direction (e.g. via a jet)
can provide information on the binary's orbital plane at large separations and
could be relevant, for instance, for studying X-shaped radio sources.

[125] oai:arXiv.org:0901.2722 [pdf] - 20453

Can magnetic fields be detected during the inspiral of binary neutron
stars?

Comments: 5 pages, 4 figures. Small changes to text and figures. Matches
version to appear on MNRAS Letters

Submitted: **2009-01-18**, last modified: **2009-09-07**

Using accurate and fully general-relativistic simulations we assess the
effect that magnetic fields have on the gravitational-wave emission produced
during the inspiral and merger of magnetized neutron stars. In particular, we
show that magnetic fields have an impact after the merger, because they are
amplified by a Kelvin-Helmholtz instability, but also during the inspiral, most
likely because the magnetic tension reduces the stellar tidal deformation for
extremely large initial magnetic fields, B_0>~10^{17}G. We quantify the
influence of magnetic fields by computing the overlap, O, between the waveforms
produced during the inspiral by magnetized and unmagnetized binaries. We find
that for any realistic magnetic field strength B_0<~10^{14}G the overlap during
the inspiral is O>~0.999 and is quite insensitive to the mass of the neutron
stars. Only for unrealistically large magnetic fields like B_0~10^{17}G the
overlap does decrease noticeably, becoming at our resolutions O<~0.76/0.67 for
stars with baryon masses M_b~1.4/1.6 Msun, respectively. Because neutron stars
are expected to merge with magnetic fields ~10^{8}-10^{10}G and because present
detectors are sensitive to O<~0.995, we conclude that it is very unlikely that
the present detectors will be able to discern the presence of magnetic fields
during the inspiral of neutron stars.

[126] oai:arXiv.org:0901.4955 [pdf] - 20870

Accurate evolutions of inspiralling neutron-star binaries: assessment of
the truncation error

Comments: 13 pages, 5 figures. Minor changes to match published version. Added
figure 5 right panel

Submitted: **2009-01-30**, last modified: **2009-05-21**

We have recently presented an investigation in full general relativity of the
dynamics and gravitational-wave emission from binary neutron stars which
inspiral and merge, producing a black hole surrounded by a torus (see
arXiv:0804.0594). We here discuss in more detail the convergence properties of
the results presented in arXiv:0804.0594 and, in particular, the deterioration
of the convergence rate at the merger and during the survival of the merged
object, when strong shocks are formed and turbulence develops. We also show
that physically reasonable and numerically convergent results obtained at
low-resolution suffer however from large truncation errors and hence are of
little physical use. We summarize our findings in an "error budget", which
includes the different sources of possible inaccuracies we have investigated
and provides a first quantitative assessment of the precision in the modelling
of compact fluid binaries.

[127] oai:arXiv.org:0710.2335 [pdf] - 900212

A template bank for gravitational waveforms from coalescing binary black
holes: non-spinning binaries

Comments: A typo fixed in Eq.(B11)

Submitted: **2007-10-11**, last modified: **2009-05-11**

Gravitational waveforms from the inspiral and ring-down stages of the binary
black hole coalescences can be modelled accurately by
approximation/perturbation techniques in general relativity. Recent progress in
numerical relativity has enabled us to model also the non-perturbative merger
phase of the binary black-hole coalescence problem. This enables us to
\emph{coherently} search for all three stages of the coalescence of
non-spinning binary black holes using a single template bank. Taking our
motivation from these results, we propose a family of template waveforms which
can model the inspiral, merger, and ring-down stages of the coalescence of
non-spinning binary black holes that follow quasi-circular inspiral. This
two-dimensional template family is explicitly parametrized by the physical
parameters of the binary. We show that the template family is not only
\emph{effectual} in detecting the signals from black hole coalescences, but
also \emph{faithful} in estimating the parameters of the binary. We compare the
sensitivity of a search (in the context of different ground-based
interferometers) using all three stages of the black hole coalescence with
other template-based searches which look for individual stages separately. We
find that the proposed search is significantly more sensitive than other
template-based searches for a substantial mass-range, potentially bringing
about remarkable improvement in the event-rate of ground-based interferometers.
As part of this work, we also prescribe a general procedure to construct
interpolated template banks using non-spinning black hole waveforms produced by
numerical relativity.

[128] oai:arXiv.org:0806.1700 [pdf] - 254898

Relativistic simulations of the phase-transition-induced collapse of
neutron stars

Comments: 28 pages, 27 figures. Minor changes to be consistent with published
version

Submitted: **2008-06-10**, last modified: **2009-05-05**

An increase in the central density of a neutron star may trigger a phase
transition from hadronic matter to deconfined quark matter in the core, causing
it to collapse to a more compact hybrid-star configuration. We present a study
of this, building on previous work by Lin et al. (2006). We follow them in
considering a supersonic phase transition and using a simplified equation of
state, but our calculations are general relativistic (using 2D simulations in
the conformally flat approximation) as compared with their 3D Newtonian
treatment. We also improved the treatment of the initial phase transformation,
avoiding the introduction of artificial convection. As before, we find that the
emitted gravitational-wave spectrum is dominated by the fundamental
quasi-radial and quadrupolar pulsation modes but the strain amplitudes are much
smaller than suggested previously, which is disappointing for the detection
prospects. However, we see significantly smaller damping and observe a
nonlinear mode resonance which substantially enhances the emission in some
cases. We explain the damping mechanisms operating, giving a different view
from the previous work. Finally, we discuss the detectability of the
gravitational waves, showing that the signal-to-noise ratio for current or
second generation interferometers could be high enough to detect such events in
our Galaxy, although third generation detectors would be needed to observe them
out to the Virgo cluster, which would be necessary for having a reasonable
event rate.

[129] oai:arXiv.org:0812.2325 [pdf] - 19434

Modelling the final state from binary black-hole coalescences

Comments: Submitted to CQG, LISA-7 Special Issue

Submitted: **2008-12-12**

Over the last few years enormous progress has been made in the numerical
description of the inspiral and merger of binary black holes. A particular
effort has gone into the modelling of the physical properties of the final
black hole, namely its spin and recoil velocity, as these quantities have
direct impact in astrophysics, cosmology and, of course, general relativity. As
numerical-relativity calculations still remain computationally very expensive
and cannot be used to investigate the complete space of possible parameters,
semi-analytic approaches have been developed and shown to reproduce with very
high precision the numerical results. I here collect and review these efforts,
pointing out the relative strengths and weaknesses, and discuss which
directions are more promising to further improve them.

[130] oai:arXiv.org:0804.0594 [pdf] - 11477

Accurate evolutions of inspiralling neutron-star binaries: prompt and
delayed collapse to black hole

Comments: 35 pages, 29 figures, corrected few typos to match the published
version. High-resolution figures and animations can be found at
http://numrel.aei.mpg.de/Visualisations/Archive/BinaryNeutronStars/Relativistic_Meudon/index.html

Submitted: **2008-04-03**, last modified: **2008-10-24**

Binary neutron-star (BNS) systems represent primary sources for the
gravitational-wave (GW) detectors. We present a systematic investigation in
full GR of the dynamics and GW emission from BNS which inspiral and merge,
producing a black hole (BH) surrounded by a torus. Our results represent the
state of the art from several points of view: (i) We use HRSC methods for the
hydrodynamics equations and high-order finite-differencing techniques for the
Einstein equations; (ii) We employ AMR techniques with "moving boxes"; (iii) We
use as initial data BNSs in irrotational quasi-circular orbits; (iv) We exploit
the isolated-horizon formalism to measure the properties of the BHs produced in
the merger; (v) Finally, we use two approaches, based either on gauge-invariant
perturbations or on Weyl scalars, to calculate the GWs. These techniques allow
us to perform accurate evolutions on timescales never reported before (ie ~30
ms) and to provide the first complete description of the inspiral and merger of
a BNS leading to the prompt or delayed formation of a BH and to its ringdown.
We consider either a polytropic or an ideal fluid EOS and show that already
with this idealized EOSs a very interesting phenomenology emerges. In
particular, we show that while high-mass binaries lead to the prompt formation
of a rapidly rotating BH surrounded by a dense torus, lower-mass binaries give
rise to a differentially rotating NS, which undergoes large oscillations and
emits large amounts of GWs. Eventually, also the NS collapses to a rotating BH
surrounded by a torus. Finally, we also show that the use of a non-isentropic
EOS leads to significantly different evolutions, giving rise to a delayed
collapse also with high-mass binaries, as well as to a more intense emission of
GWs and to a geometrically thicker torus.

[131] oai:arXiv.org:0810.1838 [pdf] - 254901

Beyond ideal MHD: towards a more realistic modeling of relativistic
astrophysical plasmas

Comments: submitted to MNRAS

Submitted: **2008-10-10**

Many astrophysical processes involving magnetic fields and quasi-stationary
processes are well described when assuming the fluid as a perfect conductor.
For these systems, the ideal-magnetohydrodynamics (MHD) description captures
the dynamics effectively and a number of well-tested techniques exist for its
numerical solution. Yet, there are several astrophysical processes involving
magnetic fields which are highly dynamical and for which resistive effects can
play an important role. The numerical modeling of such non-ideal MHD flows is
significantly more challenging as the resistivity is expected to change of
several orders of magnitude across the flow and the equations are then either
of hyperbolic-parabolic nature or hyperbolic with stiff terms. We here present
a novel approach for the solution of these relativistic resistive MHD equations
exploiting the properties of implicit-explicit (IMEX) Runge Kutta methods. By
examining a number of tests we illustrate the accuracy of our approach under a
variety of conditions and highlight its robustness when compared with
alternative methods, such as the Strang-splitting. Most importantly, we show
that our approach allows one to treat, within a unified framework, both those
regions of the flow which are fluid-pressure dominated (such as in the interior
of compact objects) and those which are instead magnetic-pressure dominated
(such as in their magnetospheres)

[132] oai:arXiv.org:0712.3541 [pdf] - 8388

On the final spin from the coalescence of two black holes

Comments: Extended discussion of physical assumptions to make them less
cryptic; matches published version

Submitted: **2007-12-20**, last modified: **2008-07-25**

We provide a compact analytic formula to compute the spin of the black hole
produced by the coalescence of two black holes. The expression, which uses an
analytic fit of numerical-relativity data and relies on four assumptions, aims
at modelling generic initial spin configurations and mass ratios. A comparison
with numerical-relativity simulations already shows very accurate agreements
with all of the numerical data available to date, but we also suggest a number
of ways in which our predictions can be further improved.

[133] oai:arXiv.org:0711.4558 [pdf] - 7495

The influence of the hydrodynamic drag from an accretion torus on
extreme mass-ratio inspirals

Comments: 21 pages, 8 figures. Phys. Rev. D, in press. Typos fixed and minor
changes; added discussion on the gravitational pull exerted by the torus

Submitted: **2007-11-28**, last modified: **2008-05-16**

We have studied extreme mass-ratio inspirals (EMRIs) in spacetimes containing
a rotating black hole and a non self-gravitating torus with constant specific
angular momentum. We have found that the effect of the hydrodynamic drag
exerted by the torus on the satellite is much smaller than the corresponding
one due to radiation reaction, for systems such as those generically expected
in AGNs and at distances from the SMBH which can be probed with LISA. However,
given the uncertainty on the parameters of these systems, there exist
configurations in which the effect of the hydrodynamic drag can be comparable
to the radiation-reaction one in phases of the inspiral which are detectable by
LISA. This is the case, for instance, for a 10^6 M_sun SMBH surrounded by a
corotating torus of comparable mass and with radius of 10^3-10^4 gravitational
radii, or for a 10^5 M_sun SMBH surrounded by a corotating 10^4 M_sun torus
with radius of 10^5 gravitational radii. Should these conditions be met in
astrophysical systems, EMRI-gravitational waves could provide a characteristic
signature of the presence of the torus. In fact, while radiation reaction
always increases the inclination of the orbit with respect to the equatorial
plane, the hydrodynamic drag from a torus corotating with the SMBH always
decreases it. However, even when initially dominating over radiation reaction,
the influence of the hydrodynamic drag decays very rapidly as the satellite
moves into the very strong-field region of the SMBH (i.e., p <~ 5M). Although
our results have been obtained for a specific class of tori, we argue that they
will be qualitatively valid also for more generic distributions of the specific
angular momentum.

[134] oai:arXiv.org:0711.1498 [pdf] - 6873

Magnetic field amplification in proto-neutron stars -- The role of the
neutron-finger instability for dynamo excitation

Comments: Accepted for publication in A&A - 10 pages; corrections after
language editing included

Submitted: **2007-11-09**, last modified: **2008-02-25**

During the first 40 s after their birth, proto-neutron stars are expected to
be subject to at least two types of instability: the convective instability and
the neutron-finger one. Both instabilities involve convective motions and hence
can trigger dynamo actions which may be responsible for the large magnetic
fields in neutron stars and magnetars. We have solved the mean-field induction
equation in a simplified one-dimensional model of both the convective and the
neutron-finger instability zones. Although very idealized, the model includes
the nonlinearities introduced by the feedback processes which tend to saturate
the growth of the magnetic field (alpha-quenching) and suppress its turbulent
diffusion (eta-quenching). The possibility of a dynamo action is studied within
a dynamical model of turbulent diffusivity where the boundary of the unstable
zone is allowed to move. We show that the dynamo action can be operative and
that the amplification of the magnetic field can still be very effective.
Furthermore, we confirm the existence of a critical spin-period, below which
the dynamo is always excited independently of the degree of differential
rotation, and whose value is related to the size of the neutron-finger
instability zone. Finally we provide a relation for the intensity of the final
field as a function of the spin of the star and of its differential rotation.
Although they were obtained by using a toy model, we expect that our results
are able to capture the qualitative and asymptotic behaviour of a mean-field
dynamo action developing in the neutron-finger instability zone. Overall, we
find that such a dynamo is very efficient in producing magnetic fields well
above equipartition and thus that it could represent a possible explanation for
the large surface magnetic fields observed in neutron stars.

[135] oai:arXiv.org:0708.3999 [pdf] - 319800

Spin Diagrams for Equal-Mass Black-Hole Binaries with Aligned Spins

Comments: 4 pages, 3 figs; small changes matching published version

Submitted: **2007-08-29**, last modified: **2008-02-23**

Binary black-hole systems with spins aligned with the orbital angular
momentum are of special interest as they may be the preferred end-state of the
inspiral of generic supermassive binary black-hole systems. In view of this, we
have computed the inspiral and merger of a large set of binary systems of
equal-mass black holes with spins aligned with the orbital angular momentum but
otherwise arbitrary. By least-square fitting the results of these simulations
we have constructed two "spin diagrams" which provide straightforward
information about the recoil velocity |v_kick| and the final black-hole spin
a_fin in terms of the dimensionless spins a_1 and a_2 of the two initial black
holes. Overall they suggest a maximum recoil velocity of |v_kick|=441.94 km/s,
and minimum and maximum final spins a_fin=0.3471 and a_fin=0.9591,
respectively.

[136] oai:arXiv.org:0710.3345 [pdf] - 6129

The final spin from the coalescence of aligned-spin black-hole binaries

Comments: Version matching the published one; small changes throughout to fit
space constraints; corrects error in vii) about spin-up/down

Submitted: **2007-10-17**, last modified: **2008-02-02**

Determining the final spin of a black-hole (BH) binary is a question of key
importance in astrophysics. Modelling this quantity in general is made
difficult by the fact that it depends on the 7-dimensional space of parameters
characterizing the two initial black holes. However, in special cases, when
symmetries can be exploited, the description can become simpler. For black-hole
binaries with unequal masses but with equal spins which are aligned with the
orbital angular momentum, we show that the use of recent simulations and basic
but exact constraints derived from the extreme mass-ratio limit allow to model
this quantity with a simple analytic expression. Despite the simple dependence,
the expression models very accurately all of the available estimates, with
errors of a couple of percent at most. We also discuss how to use the fit to
predict when a Schwarzschild BH is produced by the merger of two spinning BHs,
when the total angular momentum of the spacetime ``flips'' sign, or under what
conditions the final BH is ``spun-up'' by the merger. Finally, suggest an
extension of the fit to include unequal-spin binaries, thus potentially
providing a complete description of the final spin from the coalescence of
generic black-hole binaries with spins aligned to the orbital angular momentum.

[137] oai:arXiv.org:0704.0138 [pdf] - 11

Circular and non-circular nearly horizon-skimming orbits in Kerr
spacetimes

Comments: 15 pages, 4 figures, 5 tables. Small changes following referee's
suggestion, typos corrected. Accepted for publication in PRD

Submitted: **2007-04-02**, last modified: **2007-07-20**

We have performed a detailed analysis of orbital motion in the vicinity of a
nearly extremal Kerr black hole. For very rapidly rotating black holes (spin
a=J/M>0.9524M) we have found a class of very strong field eccentric orbits
whose angular momentum L_z increases with the orbit's inclination with respect
to the equatorial plane, while keeping latus rectum and eccentricity fixed.
This behavior is in contrast with Newtonian intuition, and is in fact opposite
to the "normal" behavior of black hole orbits. Such behavior was noted
previously for circular orbits; since it only applies to orbits very close to
the black hole, they were named "nearly horizon-skimming orbits". Our analysis
generalizes this result, mapping out the full generic (inclined and eccentric)
family of nearly horizon-skimming orbits. The earlier work on circular orbits
reported that, under gravitational radiation emission, nearly horizon-skimming
orbits tend to evolve to smaller orbit inclination, toward prograde equatorial
configuration. Normal orbits, by contrast, always demonstrate slowly growing
orbit inclination (orbits evolve toward the retrograde equatorial
configuration). Using up-to-date Teukolsky-fluxes, we have concluded that the
earlier result was incorrect: all circular orbits, including nearly
horizon-skimming ones, exhibit growing orbit inclination. Using kludge fluxes
based on a Post-Newtonian expansion corrected with fits to circular and to
equatorial Teukolsky-fluxes, we argue that the inclination grows also for
eccentric nearly horizon-skimming orbits. We also find that the inclination
change is, in any case, very small. As such, we conclude that these orbits are
not likely to have a clear and peculiar imprint on the gravitational waveforms
expected to be measured by the space-based detector LISA.

[138] oai:arXiv.org:0707.2559 [pdf] - 3164

Recoil velocities from equal-mass binary black-hole mergers: a
systematic investigation of spin-orbit aligned configurations

Comments: 24 pages, 15 figures, 5 tables

Submitted: **2007-07-17**

Binary black-hole systems with spins aligned with the orbital angular
momentum are of special interest, as studies indicate that this configuration
is preferred in nature. If the spins of the two bodies differ, there can be a
prominent beaming of the gravitational radiation during the late plunge,
causing a recoil of the final merged black hole. We perform an accurate and
systematic study of recoil velocities from a sequence of equal-mass black holes
whose spins are aligned with the orbital angular momentum, and whose individual
spins range from a = +0.584 to -0.584. In this way we extend and refine the
results of a previous study and arrive at a consistent maximum recoil of 448 +-
5 km/s for anti-aligned models as well as to a phenomenological expression for
the recoil velocity as a function of spin ratio. This relation highlights a
nonlinear behavior, not predicted by the PN estimates, and can be readily
employed in astrophysical studies on the evolution of binary black holes in
massive galaxies. An essential result of our analysis is the identification of
different stages in the waveform, including a transient due to lack of an
initial linear momentum in the initial data. Furthermore we are able to
identify a pair of terms which are largely responsible for the kick, indicating
that an accurate computation can be obtained from modes up to l=3. Finally, we
provide accurate measures of the radiated energy and angular momentum, finding
these to increase linearly with the spin ratio, and derive simple expressions
for the final spin and the radiated angular momentum which can be easily
implemented in N-body simulations of compact stellar systems. Our code is
calibrated with strict convergence tests and we verify the correctness of our
measurements by using multiple independent methods whenever possible.

[139] oai:arXiv.org:0704.3764 [pdf] - 837

Phenomenological template family for black-hole coalescence waveforms

Comments: Minor changes, Submitted to Class. Quantum Grav. (Proc. GWDAW11)

Submitted: **2007-04-30**, last modified: **2007-07-11**

Recent progress in numerical relativity has enabled us to model the
non-perturbative merger phase of the binary black-hole coalescence problem.
Based on these results, we propose a phenomenological family of waveforms which
can model the inspiral, merger, and ring-down stages of black hole coalescence.
We also construct a template bank using this family of waveforms and discuss
its implementation in the search for signatures of gravitational waves produced
by black-hole coalescences in the data of ground-based interferometers. This
template bank might enable us to extend the present inspiral searches to
higher-mass binary black-hole systems, i.e., systems with total mass greater
than about 80 solar masses, thereby increasing the reach of the current
generation of ground-based detectors.

[140] oai:arXiv.org:0706.1513 [pdf] - 2089

How to tell a gravastar from a black hole

Comments: 16 pages, 13 figures, minor improvement

Submitted: **2007-06-11**, last modified: **2007-06-22**

Gravastars have been recently proposed as potential alternatives to explain
the astrophysical phenomenology traditionally associated to black holes,
raising the question of whether the two objects can be distinguished at all.
Leaving aside the debate about the processes that would lead to the formation
of a gravastar and the astronomical evidence in their support, we here address
two basic questions: Is a gravastar stable against generic perturbations? If
stable, can an observer distinguish it from a black hole of the same mass? To
answer these questions we construct a general class of gravastars and determine
the conditions they must satisfy in order to exist as equilibrium solutions of
the Einstein equations. For such models we perform a systematic stability
analysis against axial-perturbations, computing the real and imaginary parts of
the eigenfrequencies. Overall, we find that gravastars are stable to axial
perturbations, but also that their quasi-normal modes differ from those of a
black hole of the same mass and thus can be used to discern, beyond dispute, a
gravastar from a black hole.

[141] oai:arXiv.org:gr-qc/0701109 [pdf] - 112469

WhiskyMHD: a new numerical code for general relativistic
magnetohydrodynamics

Comments: minor changes to match the published version

Submitted: **2007-01-19**, last modified: **2007-06-05**

The accurate modelling of astrophysical scenarios involving compact objects
and magnetic fields, such as the collapse of rotating magnetized stars to black
holes or the phenomenology of gamma-ray bursts, requires the solution of the
Einstein equations together with those of general-relativistic
magnetohydrodynamics. We present a new numerical code developed to solve the
full set of general-relativistic magnetohydrodynamics equations in a dynamical
and arbitrary spacetime with high-resolution shock-capturing techniques on
domains with adaptive mesh refinements. After a discussion of the equations
solved and of the techniques employed, we present a series of testbeds carried
out to validate the code and assess its accuracy. Such tests range from the
solution of relativistic Riemann problems in flat spacetime, over to the
stationary accretion onto a Schwarzschild black hole and up to the evolution of
oscillating magnetized stars in equilibrium and constructed as consistent
solutions of the coupled Einstein-Maxwell equations.

[142] oai:arXiv.org:astro-ph/0609473 [pdf] - 85065

Accurate simulations of the dynamical bar-mode instability in full
General Relativity

Comments: RevTeX4, 23 pages, 19 figures. Version in print

Submitted: **2006-09-18**, last modified: **2007-05-13**

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

[143] oai:arXiv.org:0705.1826 [pdf] - 1224

Dynamical non-axisymmetric instabilities in rotating relativistic stars

Comments: To appear on CQG, NFNR special issue. 16 pages, 5 color figures,
movies from http://www.fis.unipr.it/numrel/BarMode/ResearchBarMode.html

Submitted: **2007-05-13**

We present new results on dynamical instabilities in rapidly rotating
neutron-stars. In particular, using numerical simulations in full General
Relativity, we analyse the effects that the stellar compactness has on the
threshold for the onset of the dynamical bar-mode instability, as well as on
the appearance of other dynamical instabilities. By using an extrapolation
technique developed and tested in our previous study [1], we explicitly
determine the threshold for a wide range of compactnesses using four sequences
of models of constant baryonic mass comprising a total of 59 stellar models.
Our calculation of the threshold is in good agreement with the Newtonian
prediction and improves the previous post-Newtonian estimates. In addition, we
find that for stars with sufficiently large mass and compactness, the m=3
deformation is the fastest growing one. For all of the models considered, the
non-axisymmetric instability is suppressed on a dynamical timescale with an m=1
deformation dominating the final stages of the instability. These results,
together with those presented in [1], suggest that an m=1 deformation
represents a general and late-time feature of non-axisymmetric dynamical
instabilities both in full General Relativity and in Newtonian gravity.

[144] oai:arXiv.org:gr-qc/0612123 [pdf] - 112425

Gravitational waves from Extreme Mass Ratio Inspirals in non-pure Kerr
spacetimes

Comments: 21 pages, 12 figures. Short discussion on the accuracy of the
spacetimes added, typos corrected. Accepted for publication in PRD

Submitted: **2006-12-20**, last modified: **2007-02-26**

To investigate the imprint on the gravitational-wave emission from extreme
mass-ratio inspirals in non-pure Kerr spacetimes, we have studied the
``kludge'' waveforms generated in highly-accurate, numerically-generated
spacetimes containing a black hole and a self-gravitating, homogeneous torus
with comparable mass and spin. In order to maximize their impact on the
produced waveforms, we have considered tori that are compact, massive and close
to the central black hole, investigating under what conditions the LISA
experiment could detect their presence. Our results show that for a large
portion of the space of parameters the waveforms produced by EMRIs in these
black hole-torus systems are indistinguishable from pure-Kerr waveforms. Hence,
a ``confusion problem'' will be present for observations carried out over a
timescale below or comparable to the dephasing time.

[145] oai:arXiv.org:astro-ph/0702485 [pdf] - 89539

Dynamics of magnetized relativistic tori oscillating around black holes

Comments: 10 pages,7 figures, submitted to MNRAS

Submitted: **2007-02-19**

We present a numerical study of the dynamics of magnetized, relativistic,
non-self-gravitating, axisymmetric tori orbiting in the background spacetimes
of Schwarzschild and Kerr black holes. The initial models have a constant
specific angular momentum and are built with a non-zero toroidal magnetic field
component, for which equilibrium configurations have recently been obtained. In
this work we extend our previous investigations which dealt with purely
hydrodynamical thick discs, and study the dynamics of magnetized tori subject
to perturbations which, for the values of the magnetic field strength
considered here, trigger quasi-periodic oscillations lasting for tens of
orbital periods. Overall, we have found that the dynamics of the magnetized
tori analyzed is very similar to that found in the corresponding unmagnetized
models. The spectral distribution of the eigenfrequencies of oscillation shows
the presence of a fundamental p mode and of a series of overtones in a harmonic
ratio 2:3:.... These simulations, therefore, extend the validity of the model
of Rezzolla et al.(2003a) for explaining the high-frequency QPOs observed in
the spectra of LMXBs containing a black-hole candidate also to the case of
magnetized discs with purely toroidal magnetic field distribution. If
sufficiently compact and massive, these oscillations can also lead to the
emission of intense gravitational radiation which is potentially detectable for
sources within the Galaxy.

[146] oai:arXiv.org:gr-qc/0610131 [pdf] - 112353

On the accretion-induced QNM excitation of a Schwarzschild black hole

Comments: 15 pages, 12 figures, revised version. To appear in Phys. Rev. D

Submitted: **2006-10-26**, last modified: **2007-02-09**

By combining the numerical solution of the nonlinear hydrodynamics equations
with the solution of the linear inhomogeneous Zerilli-Moncrief and
Regge-Wheeler equations we investigate the properties of the gravitational
radiation emitted during the axisymmetric accretion of matter onto a
Schwarzschild black hole. The matter models considered include quadrupolar dust
shells and thick accretion disks, permitting us to simulate situations which
may be encountered at the end stages of stellar gravitational collapse or
binary neutron star merger. We focus on the interference pattern appearing in
the energy spectra of the emitted gravitational waves and on the amount of
excitation of the quasi-normal modes of the accreting black hole. We show that,
quite generically in the presence of accretion, the black hole ringdown is not
a simple superposition of quasi-normal modes, although the fundamental mode is
usually present and often dominates the gravitational-wave signal. We interpret
this as due to backscattering of waves off the non-exponentially decaying part
of the black-hole potential and to the finite spatial extension of the
accreting matter. Our results suggest that the black-hole QNM contributions to
the full gravitational-wave signal should be extremely small and possibly not
detectable in generic astrophysical scenarios involving the accretion of
extended distributions of matter.

[147] oai:arXiv.org:gr-qc/0608113 [pdf] - 112289

Challenging the paradigm of singularity excision in gravitational
collapse

Comments: 4 pages, 4 figures, accepted for publication on Phys. Rev. Lett

Submitted: **2006-08-26**, last modified: **2006-10-11**

A paradigm deeply rooted in modern numerical relativity calculations
prescribes the removal of those regions of the computational domain where a
physical singularity may develop. We here challenge this paradigm by performing
three-dimensional simulations of the collapse of uniformly rotating stars to
black holes without excision. We show that this choice, combined with suitable
gauge conditions and the use of minute numerical dissipation, improves
dramatically the long-term stability of the evolutions. In turn, this allows
for the calculation of the waveforms well beyond what previously possible,
providing information on the black-hole ringing and setting a new mark on the
present knowledge of the gravitational-wave emission from the stellar collapse
to a rotating black hole.

[148] oai:arXiv.org:gr-qc/0507102 [pdf] - 111918

The Exact Solution of the Riemann Problem in Relativistic MHD

Comments: 36 pages, 13 figures. Minor changes to match published version

Submitted: **2005-07-25**, last modified: **2006-07-25**

We discuss the procedure for the exact solution of the Riemann problem in
special relativistic magnetohydrodynamics (MHD). We consider both initial
states leading to a set of only three waves analogous to the ones in
relativistic hydrodynamics, as well as generic initial states leading to the
full set of seven MHD waves. Because of its generality, the solution presented
here could serve as an important test for those numerical codes solving the MHD
equations in relativistic regimes.

[149] oai:arXiv.org:gr-qc/0606059 [pdf] - 112222

A hybrid approach to black hole perturbations from extended matter
sources

Comments: 12 pages, 4 figures. Phys. Rev. D, in press

Submitted: **2006-06-13**

We present a new method for the calculation of black hole perturbations
induced by extended sources in which the solution of the nonlinear
hydrodynamics equations is coupled to a perturbative method based on
Regge-Wheeler/Zerilli and Bardeen-Press-Teukolsky equations when these are
solved in the frequency domain. In contrast to alternative methods in the time
domain which may be unstable for rotating black-hole spacetimes, this approach
is expected to be stable as long as an accurate evolution of the matter sources
is possible. Hence, it could be used under generic conditions and also with
sources coming from three-dimensional numerical relativity codes. As an
application of this method we compute the gravitational radiation from an
oscillating high-density torus orbiting around a Schwarzschild black hole and
show that our method is remarkably accurate, capturing both the basic
quadrupolar emission of the torus and the excited emission of the black hole.

[150] oai:arXiv.org:astro-ph/0602437 [pdf] - 80013

A powerful hydrodynamic booster for relativistic jets

Comments: 4 pages, 4 figures (1 in color). ApJ Letters accepted

Submitted: **2006-02-20**

Velocities close to the speed of light are a robust observational property of
the jets observed in microquasars and AGNs, and are expected to be behind much
of the phenomenology of GRBs. Yet, the mechanism boosting relativistic jets to
such large Lorentz factors is still essentially unknown. Building on recent
general-relativistic, multidimensional simulations of progenitors of short
GRBs, we discuss a new effect in relativistic hydrodynamics which can act as an
efficient booster in jets. This effect is purely hydrodynamical and occurs when
large velocities tangential to a discontinuity are present in the flow,
yielding Lorentz factors $\Gamma \sim 10^2-10^3$ or larger in flows with
moderate initial Lorentz factors. Although without a Newtonian counterpart,
this effect can be explained easily through the most elementary hydrodynamical
flow: i.e., a relativistic Riemann problem.

[151] oai:arXiv.org:astro-ph/0506702 [pdf] - 74102

Quasi-periodic Oscillations in the X-ray Light Curves from Relativistic
Tori

Comments: 4 pages, 3 figures, submitted to ApJL

Submitted: **2005-06-28**

We use a relativistic ray-tracing code to analyze the X-ray emission from a
pressure-supported oscillating relativistic torus around a black hole. We show
that a strong correlation exists between the {\it intrinsic} frequencies of the
torus normal modes and the {\it extrinsic} frequencies seen in the observed
light curve power spectrum. This correlation demonstrates the feasibility of
the oscillating-torus model to explain the multiple peaks seen in black hole
high-frequency quasi-periodic oscillations. Using an optically thin,
monochromatic emission model, we also determine how a relativistically
broadened emission line and the amplitude of the X-ray modulations are
dependent on the observer's inclination angle and on the torus oscillation
amplitudes. Observations of these features can provide important information
about the torus as well as the black hole.

[152] oai:arXiv.org:gr-qc/0503016 [pdf] - 111778

Gravitational-Wave Emission from Rotating Gravitational Collapse in
three Dimensions

Comments: 4 pages, Phys. Rev. Lett., in press

Submitted: **2005-03-04**

We present the first calculation of gravitational wave emission produced in
the gravitational collapse of uniformly rotating neutron stars to black holes
in fully three-dimensional simulations. The initial stellar models are
relativistic polytropes which are dynamically unstable and with angular
velocities ranging from slow rotation to the mass-shedding limit. An essential
aspect of these simulations is the use of progressive mesh-refinement
techniques which allow to move the outer boundaries of the computational domain
to regions where gravitational radiation attains its asymptotic form. The
waveforms have been extracted using a gauge-invariant approach in which the
numerical spacetime is matched with the non-spherical perturbations of a
Schwarzschild spacetime. Overall, the results indicate that the waveforms have
features related to the properties of the initial stellar models (in terms of
their w-mode oscillations) and of the newly produced rotating black holes (in
terms of their quasi-normal modes). While our waveforms are in good qualitative
agreement with those computed by Stark and Piran in two-dimensional
simulations, our amplitudes are about one order of magnitude smaller and this
difference is mostly likely due to our less severe pressure reduction. For a
neutron star rotating uniformly near mass-shedding and collapsing at 10 kpc,
the signal-to-noise ratio computed uniquely from the burst is S/N ~ 0.25, but
this grows to be S/N <~ 4 in the case of LIGO II.

[153] oai:arXiv.org:gr-qc/0412063 [pdf] - 111713

Computations of primordial black hole formation

Comments: 29 pages, 5 figures, accepted by Class. Quantum Grav

Submitted: **2004-12-13**, last modified: **2005-01-31**

Results are presented from general relativistic numerical computations of
primordial black-hole formation during the radiation-dominated era of the
universe. Growing-mode perturbations are specified within the linear regime and
their subsequent evolution is followed as they become nonlinear. We use a
spherically symmetric Lagrangian code and study both super-critical
perturbations, which go on to produce black holes, and sub-critical
perturbations, for which the overdensity eventually disperses into the
background medium. For super-critical perturbations, we confirm the results of
previous work concerning scaling-laws but note that the threshold amplitude for
a perturbation to lead to black-hole formation is substantially reduced when
the initial conditions are taken to represent purely growing modes. For
sub-critical cases, where an initial collapse is followed by a subsequent
re-expansion, strong compressions and rarefactions are seen for perturbation
amplitudes near to the threshold. We have also investigated the effect of
including a significant component of vacuum energy and have calculated the
resulting changes in the threshold and in the slope of the scaling law.

[154] oai:arXiv.org:gr-qc/0403029 [pdf] - 419237

Three-dimensional relativistic simulations of rotating neutron-star
collapse to a Kerr black hole

Comments: 28 pages. Animations can be found at
http://www.sissa.it/~rezzolla/Whisky/WhiskyI/

Submitted: **2004-03-06**, last modified: **2005-01-30**

We present a new three-dimensional fully general-relativistic hydrodynamics
code using high-resolution shock-capturing techniques and a conformal traceless
formulation of the Einstein equations. Besides presenting a thorough set of
tests which the code passes with very high accuracy, we discuss its application
to the study of the gravitational collapse of uniformly rotating neutron stars
to Kerr black holes. The initial stellar models are modelled as relativistic
polytropes which are either secularly or dynamically unstable and with angular
velocities which range from slow rotation to the mass-shedding limit. We
investigate the gravitational collapse by carefully studying not only the
dynamics of the matter, but also that of the trapped surfaces, i.e. of both the
apparent and event horizons formed during the collapse. The use of these
surfaces, together with the dynamical horizon framework, allows for a precise
measurement of the black-hole mass and spin. The ability to successfully
perform these simulations for sufficiently long times relies on excising a
region of the computational domain which includes the singularity and is within
the apparent horizon. The dynamics of the collapsing matter is strongly
influenced by the initial amount of angular momentum in the progenitor star
and, for initial models with sufficiently high angular velocities, the collapse
can lead to the formation of an unstable disc in differential rotation. All the
simulations performed with uniformly rotating initial data and a polytropic or
ideal-fluid equation of state show no evidence of shocks or of the presence of
matter on stable orbits outside the black hole.

[155] oai:arXiv.org:astro-ph/0411116 [pdf] - 68715

Dynamics of oscillating relativistic tori around Kerr black holes

Comments: 12 pages, 7 figures, accepted for publication by MNRAS

Submitted: **2004-11-04**

We present a comprehensive numerical study of the dynamics of relativistic
axisymmetric accretion tori with a power-law distribution of specific angular
momentum orbiting in the background spacetime of a Kerr black hole. By
combining general relativistic hydrodynamics simulations with a linear
perturbative approach we investigate the main dynamical properties of these
objects over a large parameter space. The astrophysical implications of our
results extend and improve two interesting results that have been recently
reported in the literature. Firstly, the induced quasi-periodic variation of
the mass quadrupole moment makes relativistic tori of nuclear matter densities,
as those formed during the last stages of binary neutron star mergers,
promising sources of gravitational radiation, potentially detectable by
interferometric instruments. Secondly, $p$-mode oscillations in relativistic
tori of low rest-mass densities could be used to explain high frequency
quasi-periodic oscillations observed in X-ray binaries containing a black hole
candidate under conditions more generic than those considered so far.

[156] oai:arXiv.org:astro-ph/0407642 [pdf] - 66484

Oscillations of vertically integrated relativistic tori -II.
Axisymmetric modes in a Kerr spacetime

Comments: 13 pages, MNRAS, in press

Submitted: **2004-07-30**

This is the second of a series of papers investigating the oscillation
properties of relativistic, non-selfgravitating tori orbiting around black
holes. Extending the work done in a Schwarzschild background, we here consider
the axisymmetric oscillations of vertically integrated tori in a Kerr
spacetime. The tori are modeled with a number of different non-Keplerian
distributions of specific angular momentum and we discuss how the oscillation
properties depend on these and on the rotation of the central black hole. We
first consider a local analysis to highlight the relations between acoustic and
epicyclic oscillations in a Kerr spacetime and subsequently perform a global
eigenmode analysis to compute the axisymmetric p modes. In analogy with what
found in a Schwarzschild background, these modes behave as sound waves that are
modified by rotation and are globally trapped in the torus. For constant
distributions of specific angular momentum, the eigenfrequencies appear in a
sequence 2:3:4:... which is essentially independent of the size of the disc and
of the black hole rotation. For non-constant distributions of angular momentum,
on the other hand, the sequence depends on the properties of the disc and on
the spin of the black hole, becoming harmonic for sufficiently large tori. We
also comment on how p modes could explain the high frequency quasi-periodic
oscillations observed in low-mass X-ray binaries with a black hole candidate
and the properties of an equivalent model in Newtonian physics.

[157] oai:arXiv.org:gr-qc/0406018 [pdf] - 111553

Electromagnetic Fields in the Exterior of an Oscillating Relativistic
Star -- I. General Expressions and application to a rotating magnetic dipole

Comments: MNRAS, in press

Submitted: **2004-06-04**

Relativistic stars are endowed with intense electromagnetic fields but are
also subject to oscillations of various types. We here investigate the impact
that oscillations have on the electric and magnetic fields external to a
relativistic star in vacuum. In particular, modelling the star as a
relativistic polytrope with infinite conductivity, we consider the solution of
the general relativistic Maxwell equations both in the vicinity of the stellar
surface and far from it, once a perturbative velocity field is specified. In
this first paper we present general analytic expressions that are not
specialized to a particular magnetic field topology or velocity field. However,
as a validating example and an astrophysically important application, we
consider a dipolar magnetic field and the velocity field corresponding to the
rotation of the misaligned dipole. Besides providing analytic expressions for
the electromagnetic fields produced by this configuration, we calculate, for
the first time, the general relativistic energy loss through dipolar
electromagnetic radiation. We find that the widely used Newtonian expression
under-estimates this loss by a factor of 2-6 depending on the stellar
compactness. This correction could have important consequences in the study of
the spin evolution of pulsars.

[158] oai:arXiv.org:astro-ph/0011316 [pdf] - 39323

General Relativistic Electromagnetic Fields of a Slowly Rotating
Magnetized Neutron Star. I. Formulation of the equations

Comments: A few typos corrected; matches the versions in MNRAS

Submitted: **2000-11-16**, last modified: **2004-05-07**

We present analytic solutions of Maxwell equations in the internal and
external background spacetime of a slowly rotating magnetized neutron star. The
star is considered isolated and in vacuum, with a dipolar magnetic field not
aligned with the axis of rotation. With respect to a flat spacetime solution,
general relativity introduces corrections related both to the monopolar and the
dipolar parts of the gravitational field. In particular, we show that in the
case of infinite electrical conductivity general relativistic corrections due
to the dragging of reference frames are present, but only in the expression for
the electric field. In the case of finite electrical conductivity, however,
corrections due both to the spacetime curvature and to the dragging of
reference frames are shown to be present in the induction equation. These
corrections could be relevant for the evolution of the magnetic fields of
pulsars and magnetars. The solutions found, while obtained through some
simplifying assumption, reflect a rather general physical configuration and
could therefore be used in a variety of astrophysical situations.

[159] oai:arXiv.org:astro-ph/0401078 [pdf] - 61948

A New, Simple Model for Black Hole High Frequency QPOs

Comments: To appear in Proceedings of "X-ray Timing 2003: Rossi and Beyond",
ed. P. Kaaret, F.K. Lamb, & J.H. Swank

Submitted: **2004-01-07**

Observations of X-ray emissions from binary systems have always been
considered important tools to test the validity of General Relativity in
strong-field regimes. The pairs and triplets of high frequency quasi-periodic
oscillations observed in binaries containing a black hole candidate, in
particular, have been proposed as a means to measure more directly the black
hole properties such as its mass and spin. Numerous models have been suggested
over the years to explain the QPOs and the rich phenomenology accompanying
them. Many of these models rest on a number of assumptions and are at times in
conflict with the most recent observations. We here propose a new, simple model
in which the QPOs result from basic p-mode oscillations of a non-Keplerian disc
of finite size. We show that within this new model all of the key properties of
the QPOs: a) harmonic ratios of frequencies even as the frequencies change; b)
variations in the relative strength of the frequencies with spectral energy
distribution and with photon energy; c) small and systematic changes in the
frequencies, can all be explained simply given a single reasonable assumption.

[160] oai:arXiv.org:gr-qc/0310045 [pdf] - 111354

Dynamics of thick discs around Schwarzschild-de Sitter black holes

Comments: 12 pages; A&A, in press

Submitted: **2003-10-07**

We consider the effects of a cosmological constant on the dynamics of
constant angular momentum discs orbiting Schwarzschild-de Sitter black holes.
The motivation behind this study is to investigate whether the presence of a
radial force contrasting the black hole's gravitational attraction can
influence the occurrence of the runaway instability, a robust feature of the
dynamics of constant angular momentum tori in Schwarzschild and Kerr
spacetimes. In addition to the inner cusp near the black hole horizon through
which matter can accrete onto the black hole, in fact, a positive cosmological
constant introduces also an outer cusp through which matter can leave the torus
without accreting onto the black hole. To assess the impact of this outflow on
the development of the instability we have performed time-dependent and
axisymmetric hydrodynamical simulations of equilibrium initial configurations
in a sequence of background spacetimes of Schwarzschild-de Sitter black holes
with increasing masses. The simulations have been performed with an unrealistic
value for the cosmological constant which, however, yields sufficiently small
discs to be resolved accurately on numerical grids and thus provides a first
qualitative picture of the dynamics. The calculations, carried out for a wide
range of initial conditions, show that the mass-loss from the outer cusp can
have a considerable impact on the instability, with the latter being rapidly
suppressed if the outflow is large enough.

[161] oai:arXiv.org:astro-ph/0309783 [pdf] - 59610

Mean-field dynamo action in protoneutron stars

Comments: A&A Letters, in press

Submitted: **2003-09-29**

We have investigated the turbulent mean-field dynamo action in protoneutron
stars that are subject to convective and neutron finger instabilities. While
the first one develops mostly in the inner regions of the star, the second one
is favoured in the outer regions, where the Rossby number is much smaller and a
mean-field dynamo action is more efficient. By solving the mean-field induction
equation we have computed the critical spin period below which no dynamo action
is possible and found it to be ~1 s for a wide range of stellar models. Because
this critical period is substantially longer than the characteristic spin
period of very young pulsars, we expect that a mean-field dynamo will be
effective for most protoneutron stars.

[162] oai:arXiv.org:astro-ph/0307488 [pdf] - 58201

Oscillations of vertically integrated relativistic tori -- I.
Axisymmetric modes in a Schwarzschild spacetime

Comments: 16 pages; MNRAS, in press

Submitted: **2003-07-28**

This is the first of a series of papers investigating the oscillation
properties of relativistic, non-selfgravitating tori orbiting around a black
hole. In this initial paper we consider the axisymmetric oscillation modes of a
torus constructed in a Schwarzschild spacetime. To simplify the treatment and
make it as analytical as possible, we build our tori with vertically integrated
and vertically averaged quantities. Our investigation progresses by steps. We
first consider a local analysis in Newtonian gravity and determine the
properties of acoustic wave propagation within these objects, as well as the
relations between acoustic and epicyclic oscillations. Next, we extend the
local analysis to a general relativistic framework. Finally, we perform a
global analysis and determine both the eigenfunctions and the eigenfrequencies
of the axisymmetric oscillations corresponding to the p modes of relativistic
tori. These behave as sound waves globally trapped in the torus and possess
eigenfrequencies appearing in the simple sequence 2:3:4:..., independently of
the distribution of angular momentum considered. The properties of the modes
investigated here are in good agreement with those observed in recent numerical
simulations and could have a number of different applications. In X-ray binary
systems containing a black hole candidate, for instance, p-mode oscillations
could be used to explain the harmonic relations in the high frequency
quasi-periodic oscillations observed. In systems comprising a massive torus
orbiting a black hole, on the other hand, p-mode oscillations could be used to
explain the development or the suppression of the runaway instability.

[163] oai:arXiv.org:astro-ph/0307487 [pdf] - 58200

A New Simple Model for High Frequency Quasi Periodic Oscillations in
Black Hole Candidates

Comments: 5 pages; MNRAS Letters, in press

Submitted: **2003-07-28**

Observations of X-ray emissions from binary systems have long since been
considered important tools to test General Relativity in strong-field regimes.
The high frequency quasi-periodic oscillations (HFQPOs) observed in binaries
containing a black hole candidate, in particular, have been proposed as a means
to measure more directly the black hole properties such as its mass and spin.
Numerous models have been suggested to explain the HFQPOs and the rich
phenomenology accompanying them. Many of these models rest on a number of
assumptions and are at times in conflict with the most recent observations. We
here propose a new, simple model in which the HFQPOs result from basic $p$-mode
oscillations of a small accretion torus orbiting close to the black hole. We
show that within this model the key properties of the HFQPOs can be explained
simply, given a single reasonable assumption. We also discuss observational
tests that can falsify the model.

[164] oai:arXiv.org:gr-qc/0303057 [pdf] - 111206

Linear and Nonlinear Dynamics of Relativistic Tori

Comments: Invited Talk at the 15th SIGRAV Conference, Frascati Sept. 2002; IOP
Conference Series, one reference added

Submitted: **2003-03-14**, last modified: **2003-03-21**

I present results of two-dimensional general relativistic hydrodynamical
simulations of constant specific angular momentum tori orbiting around a
Schwarzschild black hole. After introducing axisymmetric perturbations, these
objects either become unstable to the runaway instability or respond with
regular oscillations. The latter, in particular, are responsible for
quasi-periodic bursts of accretion onto the black hole as well as for the
emission of intense gravitational radiation, with signal-to-noise ratios at the
detector which are comparable or even larger than the typical ones expected in
stellar-core collapse.

[165] oai:arXiv.org:gr-qc/0302025 [pdf] - 111167

Gravitational Waves from Perturbed Black Holes and Relativistic Stars

Comments: Lectures given at the Summer School on Astroparticle Physics and
Cosmology, ICTP, July 2002 (47 pages, 9 figs)

Submitted: **2003-02-07**

These lectures aim at providing an introduction to the properties of
gravitational waves and in particular to those gravitational waves that are
expected as a consequence of perturbations of black holes and neutron stars.
Imprinted in the gravitational radiation emitted by these objects is, in fact,
a wealth of physical information. In the case of black holes, a detailed
knowledge of the gravitational radiation emitted as a response to perturbations
will reveal us important details about their mass and spin, but also about the
fundamental properties of the event horizon. In the case of neutron stars, on
the other hand, this information can provide a detailed map of their internal
structure and tell us about the equation of state of matter at very high
density, thus filling-in a gap in energies and densities that cannot be
investigated by experiments in terrestrial laboratories.

[166] oai:arXiv.org:gr-qc/0210018 [pdf] - 111089

Quasi-periodic accretion and gravitational waves from oscillating
"toroidal neutron stars" around a Schwarzschild black hole

Comments: 17 pages, 12 figures. Animations can be found at
http://www.sissa.it/~rezzolla/movies.html

Submitted: **2002-10-05**

We present general relativistic hydrodynamics simulations of constant
specific angular momentum tori orbiting a Schwarzschild black hole. These tori
are expected to form as a result of stellar gravitational collapse, binary
neutron star merger or disruption, can reach very high rest-mass densities and
behave effectively as neutron stars but with a toroidal topology (i.e.
``toroidal neutron stars''). Our attention is here focussed on the dynamical
response of these objects to axisymmetric perturbations. We show that, upon the
introduction of perturbations, these systems either become unstable to the
runaway instability or exhibit a regular oscillatory behaviour resulting in a
quasi-periodic variation of the accretion rate as well as of the mass
quadrupole. The latter, in particular, is responsible for the emission of
intense gravitational radiation whose signal-to-noise ratio at the detector is
comparable or larger than the typical one expected in stellar-core collapse,
making these new sources of gravitational waves potentially detectable. We
discuss a systematic investigation of the parameter space both in the linear
and nonlinear regimes, providing estimates of how the gravitational radiation
emitted depends on the mass of the torus and on the strength of the
perturbation.

[167] oai:arXiv.org:gr-qc/0205034 [pdf] - 110983

An Improved Exact Riemann Solver for Multidimensional Relativistic Flows

Comments: 20 pages, 8 figures

Submitted: **2002-05-08**

We extend our approach for the exact solution of the Riemann problem in
relativistic hydrodynamics to the case in which the fluid velocity has
components tangential to the initial discontinuity. As in one-dimensional
flows, we here show that the wave-pattern produced in a multidimensional
relativistic Riemann problem can be predicted entirely by examining the initial
conditions. Our method is logically very simple and allows for a numerical
implementation of an exact Riemann solver which is both straightforward and
computationally efficient. The simplicity of the approach is also important for
revealing special relativistic effects responsible for a smooth transition from
one wave-pattern to another when the tangential velocities in the initial
states are suitably varied. While the content of this paper is focussed on a
flat spacetime, the local Lorentz invariance allows its use also in fully
general relativistic calculations.

[168] oai:arXiv.org:gr-qc/0204073 [pdf] - 110967

New Relativistic Effects in the Dynamics of Nonlinear Hydrodynamical
Waves

Comments: 4 pages, 5 figures

Submitted: **2002-04-23**

In Newtonian and relativistic hydrodynamics the Riemann problem consists of
calculating the evolution of a fluid which is initially characterized by two
states having different values of uniform rest-mass density, pressure and
velocity. When the fluid is allowed to relax, one of three possible
wave-patterns is produced, corresponding to the propagation in opposite
directions of two nonlinear hydrodynamical waves. New effects emerge in a
special relativistic Riemann problem when velocities tangential to the initial
discontinuity surface are present. We show that a smooth transition from one
wave-pattern to another can be produced by varying the initial tangential
velocities while otherwise maintaining the initial states unmodified. These
special relativistic effects are produced by the coupling through the
relativistic Lorentz factors and do not have a Newtonian counterpart.

[169] oai:arXiv.org:gr-qc/0112017 [pdf] - 110897

Frequencies of f-modes in differentially rotating relativistic stars and
secular stability limits

Comments: 4pages,3figures. To appear in ApJL

Submitted: **2001-12-10**, last modified: **2002-02-15**

We have computed the eigenfrequencies of f-modes for a constant-rest-mass
sequences of rapidly rotating relativistic inviscid stars in differential
rotation. The frequencies have been calculated neglecting the metric
perturbations (the relativistic Cowling approximation) and expressed as a
function of the ratio between the rotational kinetic energy and the absolute
value of the gravitational energy of the stellar model beta=T/|W|. The zeros
and the end-points of these sequences mark respectively the onset of the
secular instability driven by gravitational radiation-reaction and the maximum
value of beta at which an equilibrium model exists. In differentially rotating
stars the secular stability limits appear at a beta larger than those found for
uniformly rotating stars. Differential rotation, on the other hand, also allows
for the existence of equilibrium models at values of beta larger than those for
uniformly rotating stars, moving the end-point of the sequences to larger beta.
As a result, for some degrees of differential rotation, the onset of the
secular instability for f-modes is generally favoured by the presence of
differential rotation.

[170] oai:arXiv.org:gr-qc/0112042 [pdf] - 110899

General Relativistic Rossby-Haurwitz waves of a slowly and
differentially rotating fluid shell

Comments: 4 pages. Comment to Class. Quantum Grav. 18(2001)L87

Submitted: **2001-12-18**, last modified: **2001-12-25**

We show that, at first order in the angular velocity, the general
relativistic description of Rossby-Haurwitz waves (the analogues of r-waves on
a thin shell) can be obtained from the corresponding Newtonian one after a
coordinate transformation. As an application, we show that the results recently
obtained by Rezzolla and Yoshida (2001) in the analysis of Newtonian
Rossby-Haurwitz waves of a slowly and differentially rotating, fluid shell
apply also in General Relativity, at first order in the angular velocity.

[171] oai:arXiv.org:astro-ph/0112032 [pdf] - 46405

General Relativistic Electromagnetic Fields of a Slowly Rotating
Magnetized Neutron Star. II. Solution of the Induction Equations

Comments: 13 pages, 5 figures. Accepted for publication by MNRAS. Replaces
previous version without unnecessary mn.sty

Submitted: **2001-12-03**

We have solved numerically the general relativistic induction equations in
the interior background spacetime of a slowly rotating magnetized neutron star.
The analytic form of these equations was discussed in a recent paper (Rezzolla
et al 2001a), where corrections due both to the spacetime curvature and to the
dragging of reference frames were shown to be present. Through a number of
calculations we have investigated the evolution of the magnetic field with
different rates of stellar rotation, different inclination angles between the
magnetic moment and the rotation axis, as well as different values of the
electrical conductivity. All of these calculations have been performed for a
constant temperature relativistic polytropic star and make use of a consistent
solution of the initial value problem which avoids the use of artificial
analytic functions. Our results show that there exist general relativistic
effects introduced by the rotation of the spacetime which tend to decrease the
decay rate of the magnetic field. The rotation-induced corrections are however
generally hidden by the high electrical conductivity of the neutron star matter
and when realistic values for the electrical conductivity are considered, these
corrections become negligible even for the fastest known pulsar.

[172] oai:arXiv.org:gr-qc/0110047 [pdf] - 110857

Three-dimensional general relativistic hydrodynamics II: long-term
dynamics of single relativistic stars

Comments: 19 pages, 17 figures

Submitted: **2001-10-09**

This is the second in a series of papers on the construction and validation
of a three-dimensional code for the solution of the coupled system of the
Einstein equations and of the general relativistic hydrodynamic equations, and
on the application of this code to problems in general relativistic
astrophysics. In particular, we report on the accuracy of our code in the
long-term dynamical evolution of relativistic stars and on some new physics
results obtained in the process of code testing. The tests involve single
non-rotating stars in stable equilibrium, non-rotating stars undergoing radial
and quadrupolar oscillations, non-rotating stars on the unstable branch of the
equilibrium configurations migrating to the stable branch, non-rotating stars
undergoing gravitational collapse to a black hole, and rapidly rotating stars
in stable equilibrium and undergoing quasi-radial oscillations. The numerical
evolutions have been carried out in full general relativity using different
types of polytropic equations of state using either the rest-mass density only,
or the rest-mass density and the internal energy as independent variables. New
variants of the spacetime evolution and new high resolution shock capturing
(HRSC) treatments based on Riemann solvers and slope limiters have been
implemented and the results compared with those obtained from previous methods.
Finally, we have obtained the first eigenfrequencies of rotating stars in full
general relativity and rapid rotation. A long standing problem, such
frequencies have not been obtained by other methods. Overall, and to the best
of our knowledge, the results presented in this paper represent the most
accurate long-term three-dimensional evolutions of relativistic stars available
to date.

[173] oai:arXiv.org:gr-qc/0108057 [pdf] - 110829

Stationary Electromagnetic Fields of a Slowly Rotating Magnetized
Neutron Star in General Relativity

Comments: 12 pages, no figures

Submitted: **2001-08-23**

Following the general formalism presented by Rezzolla, Ahmedov and Miller
(MNRAS, 322, 723 2001), we here derive analytic solutions of the
electromagnetic fields equations in the internal and external background
spacetime of a slowly rotating highly conducting magnetized neutron star. The
star is assumed to be isolated and in vacuum, with a dipolar magnetic field not
aligned with the axis of rotation. Our results indicate that the
electromagnetic fields of a slowly rotating neutron star are modified by
general relativistic effects arising from both the monopolar and the dipolar
parts of the gravitational field. The results presented here differ from the
ones discussed by Rezzolla, Ahmedov and Miller (MNRAS, 322, 723 2001) mainly in
that we here consider the interior magnetic field to be dipolar with the same
radial dependence as the external one. While this assumption might not be a
realistic one, it should be seen as the application of our formalism to a case
often discussed in the literature.

[174] oai:arXiv.org:gr-qc/0103005 [pdf] - 1578308

An Improved Exact Riemann Solver for Relativistic Hydrodynamics

Comments: Extended the discussion on the computational efficiency with improved
estimates; matches the version in press

Submitted: **2001-03-02**, last modified: **2001-07-19**

A Riemann problem with prescribed initial conditions will produce one of
three possible wave patterns corresponding to the propagation of the different
discontinuities that will be produced once the system is allowed to relax. In
general, when solving the Riemann problem numerically, the determination of the
specific wave pattern produced is obtained through some initial guess which can
be successively discarded or improved. We here discuss a new procedure,
suitable for implementation in an exact Riemann solver in one dimension, which
removes the initial ambiguity in the wave pattern. In particular we focus our
attention on the relativistic velocity jump between the two initial states and
use this to determine, through some analytic conditions, the wave pattern
produced by the decay of the initial discontinuity. The exact Riemann problem
is then solved by means of calculating the root of a nonlinear equation.
Interestingly, in the case of two rarefaction waves, this root can even be
found analytically. Our procedure is straightforward to implement numerically
and improves the efficiency of numerical codes based on exact Riemann solvers.

[175] oai:arXiv.org:gr-qc/0107061 [pdf] - 110814

Properties of r modes in rotating magnetic neutron stars. I. Kinematic
Secular Effects and Magnetic Evolution Equations

Comments: 15 pages, 4 figures. To appear in Phys. Rev. D

Submitted: **2001-07-17**

The instability of r-mode oscillations in rapidly rotating neutron stars has
attracted attention as a potential mechanism for producing high frequency,
almost periodic gravitational waves. The analyses carried so far have shown the
existence of these modes and have considered damping by shear and bulk
viscosity. However, the magnetohydrodynamic coupling of the modes with a
stellar magnetic field and its role in the damping of the instability has not
been fully investigated yet. Following our introductory paper (Rezzolla, Lamb
and Shapiro 2000), we here discuss in more detail the existence of secular
higher-order kinematical effects which will produce toroidal fluid drifts. We
also define the sets of equations that account for the time evolution of the
magnetic fields produced by these secular velocity fields and show that the
magnetic fields produced can reach equipartition in less than a year. The full
numerical calculations as well as the evaluation of the impact of strong
magnetic fields on the onset and evolution of the r-mode instability will be
presented in a companion paper.

[176] oai:arXiv.org:gr-qc/0107062 [pdf] - 110815

Properties of r modes in rotating magnetic neutron stars. II. Evolution
of the r modes and stellar magnetic field

Comments: 16 pages, 12 figures. To appear in Phys. Rev. D

Submitted: **2001-07-17**

The evolution of the r-mode instability is likely to be accompanied by
secular kinematic effects which will produce differential rotation with large
scale drifts of fluid elements, mostly in the azimuthal direction. As first
discussed by Rezzolla, Lamb and Shapiro 2000, the interaction of these secular
velocity fields with a pre-existing neutron star magnetic field could result in
the generation of intense and large scale toroidal fields. Following their
derivation in the companion paper, we here discuss the numerical solution of
the evolution equations for the magnetic field. The values of the magnetic
fields obtained in this way are used to estimate the conditions under which the
r-mode instability might be prevented or suppressed. We also assess the impact
of the generation of large magnetic fields on the gravitational wave
detectability of r-mode unstable neutron stars. Our results indicate that the
signal to noise ratio in the detection of gravitational waves from the r-mode
instability might be considerably decreased if the latter develops in neutron
stars with initial magnetic fields larger than 10^10 G.

[177] oai:arXiv.org:gr-qc/0107027 [pdf] - 142749

Rossby-Haurwitz waves of a slowly and differentially rotating fluid
shell

Comments: 8pages, 3figures. Submitted to CQG

Submitted: **2001-07-09**

Recent studies have raised doubts about the occurrence of r modes in
Newtonian stars with a large degree of differential rotation. To assess the
validity of this conjecture we have solved the eigenvalue problem for
Rossby-Haurwitz waves (the analogues of r waves on a thin-shell) in the
presence of differential rotation. The results obtained indicate that the
eigenvalue problem is never singular and that, at least for the case of a
thin-shell, the analogues of r modes can be found for arbitrarily large degrees
of differential rotation. This work clarifies the puzzling results obtained in
calculations of differentially rotating axi-symmetric Newtonian stars.

[178] oai:arXiv.org:gr-qc/0107019 [pdf] - 110807

Gravitational Wave Emission by Cataclysmic Variables: numerical models
of semi-detached binaries

Comments: 7pages, 4figures. To appear in MNRAS

Submitted: **2001-07-05**

Gravitational wave emission is considered to be the driving force for the
evolution of short-period cataclysmic binary stars, making them a potential
test for the validity of General Relativity. In spite of continuous refinements
of the physical description, a 10% mismatch exists between the theoretical
minimum period ($P_{\rm turn} \simeq 70$ min) and the short-period cut-off
($P_{\rm min} \simeq 80$ min) observed in the period distribution for
cataclysmic variable binaries. A possible explanation for this mismatch was
associated with the use of the Roche model. We here present a systematic
comparison between self-consistent, numerically constructed sequences of
hydrostatic models of binary stars and Roche models of semi-detached binaries.
On the basis of our approach, we also derive a value for the minimum period of
cataclysmic variable binaries. The results obtained through the comparison
indicate that the Roche model is indeed very good, with deviations from the
numerical solution which are of a few percent at most. Our results therefore
suggest that additional sources of angular momentum loss or alternative
explanations need to be considered in order to justify the mismatch.

[179] oai:arXiv.org:astro-ph/9911188 [pdf] - 109316

R-Mode Oscillations in Rotating Magnetic Neutron Stars

Comments: 4 pages, 1 postscript figure, uses emulateapj; submitted to ApJ
Letters 1999 Nov 8; accepted 2000 Jan 25; this version is essentially
identical to the original version except that Figure 2 was deleted in order
to fit within the ApJ Letters page limit

Submitted: **1999-11-10**, last modified: **2000-01-25**

We show that r-mode oscillations distort the magnetic fields of neutron stars
and that their occurrence is likely to be limited by this interaction. If the
field is gtrsim 10^{16} (Omega/Omega_B) G, where Omega and Omega_B are the
angular velocities of the star and at which mass shedding occurs, r-mode
oscillations cannot occur. Much weaker fields will prevent gravitational
radiation from exciting r-mode oscillations or damp them on a relatively short
timescale by extracting energy from the modes faster than gravitational wave
emission can pump energy into them. For example, a 10^{10} G poloidal magnetic
field that threads the star's superconducting core is likely to prevent the
ell=2 mode from being excited unless Omega exceeds 0.35 Omega_B. If Omega is
larger than 0.35 Omega_B initially, the ell=2 mode may be excited but is likely
to decay rapidly once Omega falls below 0.35 Omega_B, which happens in lesssim
15^d if the saturation amplitude is gtrsim 0.1. The r-mode oscillations may
play an important role in determining the structure of neutron star magnetic
fields.

[180] oai:arXiv.org:gr-qc/9907098 [pdf] - 113224

Implementing fully relativistic hydrodynamics in three dimensions

Comments: 5 Pages, 4 Figures, submitted to Proceedings of the 8th Canadian
Conference on General Relativity and Relativistic Astrophysics

Submitted: **1999-07-28**

We report on our numerical implementation of fully relativistic hydrodynamics
coupled to Einstein's field equations in three spatial dimensions. We briefly
review several steps in our code development, including our recasting of
Einstein's equations and several tests which demonstrate its advantages for
numerical integrations. We outline our implementation of relativistic
hydrodynamics, and present numerical results for the evolution of both stable
and unstable Oppenheimer-Volkov equilibrium stars, which represent a very
promising first test of our code.

[181] oai:arXiv.org:gr-qc/9905027 [pdf] - 113170

Constructing a Mass-Current Radiation-Reaction Force For Numerical
Simulations

Comments: 17 pages, Latex2e, emulateapj.sty. To appear in ApJ

Submitted: **1999-05-08**

We present a new set of 3.5 Post-Newtonian equations in which Newtonian
hydrodynamics is coupled to the nonconservative effects of gravitational
radiation emission. Our formalism differs in two significant ways from a
similar 3.5 Post-Newtonian approach proposed by Blanchet (1993, 1997). Firstly
we concentrate only on the radiation-reaction effects produced by a
time-varying mass-current quadrupole $S_{ij}$. Secondly, we adopt a gauge in
which the radiation-reaction force densities depend on the fourth time
derivative of $S_{ij}$, rather than on the fifth, as in Blanchet's approach.
This difference makes our formalism particularly well-suited to numerical
implementation and could prove useful in performing fully numerical simulations
of the recently discovered $r$-mode instability for rotating neutron stars
subject to axial perturbations.

[182] oai:arXiv.org:hep-ph/9610225 [pdf] - 115936

Baryon number segregation at the end of the cosmological quark-hadron
transition

Comments: 33 pages, Latex file, 6 postscript figures included in the text
(psfig.tex). To appear in Phys. Rev. D15

Submitted: **1996-10-02**

One of the most interesting questions regarding a possible first order
cosmological quark--hadron phase transition concerns the final fate of the
baryon number contained within the disconnected quark regions at the end of the
transition. We here present a detailed investigation of the hydrodynamical
evolution of an evaporating quark drop, using a multi-component fluid
description to follow the mechanisms of baryon number segregation. With this
approach, we are able to take account of the simultaneous effects of baryon
number flux suppression at the phase interface, entropy extraction by means of
particles having long mean-free-paths, and baryon number diffusion. A range of
computations has been performed to investigate the permitted parameter-space
and this has shown that significant baryon number concentrations, perhaps even
up to densities above that of nuclear matter, represent an inevitable outcome
within this scenario.

[183] oai:arXiv.org:astro-ph/9605033 [pdf] - 94589

Stability of Cosmological Detonation Fronts

Comments: 39 pages, LaTeX file plus 3 embedded PostScript figures; to be
published in Phys. Rev. D (15 July issue)

Submitted: **1996-05-08**

We present results of a linear stability analysis of relativistic detonation
fronts, which have been considered as representing phase interfaces in
cosmological first order phase transitions. After discussing general stability
conditions for detonation fronts, we concentrate on the properties of the
fronts with respect to corrugation instabilities and discuss separately the
cases of Chapman-Jouguet and strong detonation waves. Contrarily to what
recently claimed, we find that strong detonations are both evolutionary and
stable with respect to corrugations of the front. Moreover, Chapman-Jouguet
detonations appear to be unconditionally linearly stable. The implications of
the stability results for first order cosmological phase transitions are
presented and a discussion of the causal structure of reaction fronts is also
given.

[184] oai:arXiv.org:astro-ph/9510039 [pdf] - 93412

Evaporation of cosmological quark drops and relativistic radiative
transfer

Comments: 32 pages, latex; 12 postscript figures are available in a separate
uuencoded compressed file

Submitted: **1995-10-06**

We discuss the results of a full relativistic treatment of the hydrodynamics
of disconnected quark regions during the final stages of the cosmological
quark-hadron transition. In this study, which represents a further development
of a previous analysis of the evaporation of cosmological quark drops, the
effects of long range energy and momentum transfer via electromagnetically
interacting particles are consistently taken into account. For this purpose, a
set of relativistic Lagrangian equations describing the evolution of the
strongly interacting fluids is coupled to a system of equivalent equations
accounting for the hydrodynamics of the fluid of electromagnetically
interacting particles. The complete set of equations has then been solved
numerically and results are presented from this. The inclusion of relativistic
radiative transfer produces significantly different results, with the formation
of high density regions at the end of the drop evaporation being particularly
relevant. A comparison is made with results obtained for the previous
radiation-free model and the cosmological implications concerning baryon number
concentrations are briefly discussed.

[185] oai:arXiv.org:astro-ph/9502064 [pdf] - 1469187

EVAPORATION OF QUARK DROPS DURING THE COSMOLOGICAL Q-H TRANSITION

Comments: 12 pages Phys. Rev. format, uuencoded postscript file including 12
figures

Submitted: **1995-02-13**

We have carried out a study of the hydrodynamics of disconnected quark
regions during the final stages of the cosmological quark-hadron transition. A
set of relativistic Lagrangian equations is presented for following the
evaporation of a single quark drop and results from the numerical solution of
this are discussed. A self-similar solution is shown to exist and the formation
of baryon number density inhomogeneities at the end of the drop contraction is
discussed.

[186] oai:arXiv.org:astro-ph/9411091 [pdf] - 92073

Hydrodynamics of the Cosmological Quark-Hadron Transition in the
Presence of Long-Range Energy and Momentum Transfer

Comments: 21 pages, Plain Tex, nofigures. uucoded file containing the postcript
version of the paper and 9 figures is available on a separate file

Submitted: **1994-11-22**

Two previous papers in this series have presented a study of the growth of
hadronic bubbles during the cosmological Quark--Hadron transition, treating the
material within each phase as a single perfect fluid. Here, we extend the
analysis to include the effects of long-range energy and momentum transfer by
weakly and electromagnetically interacting particles. After a short review of
the formalism adopted, we discuss the numerical strategies used in the computer
code which has been constructed in order to solve this system of equations.
Results for the growth of single hadronic bubbles are also presented.

[187] oai:arXiv.org:astro-ph/9406055 [pdf] - 142694

Relativistic Radiative Transfer for Spherical Flows

Comments: 28 pages, 9 postscript figs, Plain TeX

Submitted: **1994-06-18**

We present a new complete set of Lagrangian relativistic hydrodynamical
equations describing the transfer of energy and momentum between a standard
fluid and a radiation fluid in a general non-stationary spherical flow. The new
set of equations has been derived for a particular application to the study of
the cosmological Quark--Hadron transition but can also be used in other
contexts.