Normalized to: Galeazzi, F.

**19** article(s) in total. **297** co-authors, from **1** to **10** common article(s). Median position in authors list is **4,0**.

[1] oai:arXiv.org:1709.02643 [pdf] - 1593610

Modeling differential rotations of compact stars in equilibriums

Comments: 7 pages, 4 figures

Submitted: **2017-09-08**, last modified: **2017-11-28**

Outcomes of numerical relativity simulations of massive core collapses or
binary neutron star mergers with moderate masses suggest formations of rapidly
and differentially rotating neutron stars. Subsequent fall back accretion may
also amplify the degree of differential rotations. We propose new formulations
for modeling differential rotations of those compact stars, and present
selected solutions of differentially rotating, stationary, and axisymmetric
compact stars in equilibriums. For the cases when rotating stars reach break-up
velocities, the maximum masses of such rotating models are obtained.

[2] 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.

[3] oai:arXiv.org:1606.04604 [pdf] - 1579985

Do triaxial supramassive compact stars exist?

Comments: 5 pages, 3 figures. Revised version to appear in PRD rapid
communication

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

We study quasiequilibrium solutions of triaxially deformed rotating compact
stars -- a generalization of Jacobi ellipsoids under relativistic gravity and
compressible equations of state (EOS). For relatively stiff (piecewise)
polytropic EOSs, we find supramassive triaxial solutions whose masses exceed
the maximum mass of the spherical solution, but are always lower than those of
axisymmetric equilibriums. The difference in the maximum masses of triaxial and
axisymmetric solutions depends sensitively on the EOS. If the difference turns
out to be only about 10%, it will be strong evidence that the EOS of high
density matter becomes substantially softer in the core of neutron stars. This
finding opens a novel way to probe phase transitions of high density nuclear
matter using detections of gravitational waves from new born neutron stars or
magnetars under fallback accretion.

[4] oai:arXiv.org:1512.06397 [pdf] - 1439719

How loud are neutron star mergers?

Comments: 6 pages, 3 figures. Waveforms available at
https://zenodo.org/record/57844

Submitted: **2015-12-20**, last modified: **2016-07-13**

We present results from the first large parameter study of neutron star
mergers using fully general relativistic simulations with finite-temperature
microphysical equations of state and neutrino cooling. We consider equal and
unequal-mass binaries drawn from the galactic population and simulate each
binary with three different equations of state. Our focus is on the emission of
energy and angular momentum in gravitational waves in the postmerger phase. We
find that the emitted gravitational-wave energy in the first
$\sim$$10\,\mathrm{ms}$ of the life of the resulting hypermassive neutron star
(HMNS) is about twice the energy emitted over the entire inspiral history of
the binary. The total radiated energy per binary mass is comparable to or
larger than that of nonspinning black hole inspiral-mergers. About $0.8-2.5\%$
of the binary mass-energy is emitted at kHz frequencies in the early HMNS
evolution. We find a clear dependence of the postmerger GW emission on binary
configuration and equation of state and show that it can be encoded as a broad
function of the binary tidal coupling constant $\kappa^T_2$. Our results also
demonstrate that the dimensionless spin of black holes resulting from
subsequent HMNS collapse are limited to $\lesssim0.7-0.8$. This may
significantly impact the neutrino pair annihilation mechanism for powering
short gamma-ray bursts (sGRB).

[5] oai:arXiv.org:1605.02629 [pdf] - 1442772

On the dynamics of tilted black hole-torus systems

Comments: 11 pages, 11 figures, accepted for publication in MNRAS

Submitted: **2016-05-09**, last modified: **2016-06-17**

We present results from three-dimensional, numerical relativity simulations
of a {\it tilted} black hole-thick accretion disc system. The simulations are
analysed using tracer particles in the disc which are advected with the flow.
Such tracers, which we employ in these new simulations for the first time,
provide a powerful means to analyse in detail the complex dynamics of tilted
black hole-torus systems. We show how its use helps to gain insight in the
overall dynamics of the system, discussing the origin of the observed black
hole precession and the development of a global non-axisymmetric $m=1$ mode in
the disc. Our three-dimensional simulations show the presence of quasi-periodic
oscillations (QPOs) in the instantaneous accretion rate, with frequencies in a
range compatible with those observed in low mass X-ray binaries with either a
black hole or a neutron star component. The frequency ratio of the dominant low
frequency peak and the first overtone is $o_1/f \sim 1.9$, a frequency ratio
not attainable when modelling the QPOs as $p$-mode oscillations in axisymmetric
tori.

[6] 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.

[7] oai:arXiv.org:1506.04056 [pdf] - 1439692

Numerical relativity simulations of thick accretion disks around tilted
Kerr black holes

Comments: 35 pages, 33 figures, accepted for publication in PRD

Submitted: **2015-06-12**, last modified: **2016-02-27**

In this work we present 3D numerical relativity simulations of thick
accretion disks around tilted Kerr BH. We investigate the evolution of three
different initial disk models with a range of initial black hole spin
magnitudes and tilt angles. For all the disk-to-black hole mass ratios
considered (0.044-0.16) we observe significant black hole precession and
nutation during the evolution. This indicates that for such mass ratios,
neglecting the self-gravity of the disks by evolving them in a fixed background
black hole spacetime is not justified. We find that the two more massive models
are unstable against the Papaloizou-Pringle (PP) instability and that those
PP-unstable models remain unstable for all initial spins and tilt angles
considered, showing that the development of the instability is a very robust
feature of such PP-unstable disks. Our lightest model, which is the most
astrophysically favorable outcome of mergers of binary compact objects, is
stable. The tilt between the black hole spin and the disk is strongly modulated
during the growth of the PP instability, causing a partial global realignment
of black hole spin and disk angular momentum in the most massive model with
constant specific angular momentum l. For the model with non-constant l-profile
we observe a long-lived m=1 non-axisymmetric structure which shows strong
oscillations of the tilt angle in the inner regions of the disk. This effect
might be connected to the development of Kozai-Lidov oscillations. Our
simulations also confirm earlier findings that the development of the PP
instability causes the long-term emission of large amplitude gravitational
waves, predominantly for the l=m=2 multipole mode. The imprint of the BH
precession on the gravitational waves from tilted BH-torus systems remains an
interesting open issue that would require significantly longer simulations than
those presented in this work.

[8] oai:arXiv.org:1411.7975 [pdf] - 977292

Properties of hypermassive neutron stars formed in mergers of spinning
binaries

Comments: 17 pages, 20 figures Updated to version published in PRD

Submitted: **2014-11-28**, last modified: **2015-04-08**

We present numerical simulations of binary neutron star mergers, comparing
irrotational binaries to binaries of NSs rotating aligned to the orbital
angular momentum. For the first time, we study spinning BNSs employing nuclear
physics equations of state, namely the ones of Lattimer and Swesty as well as
Shen, Horowitz, and Teige. We study mainly equal mass systems leading to a
hypermassive neutron star (HMNS), and analyze in detail its structure and
dynamics. In order to exclude gauge artifacts, we introduce a novel coordinate
system used for post-processing. The results for our equal mass models show
that the strong radial oscillations of the HMNS modulate the instantaneous
frequency of the gravitational wave (GW) signal to an extend that leads to
separate peaks in the corresponding Fourier spectrum. In particular, the high
frequency peaks which are often attributed to combination frequencies can also
be caused by the modulation of the m=2 mode frequency in the merger phase. As a
consequence for GW data analysis, the offset of the high frequency peak does
not necessarily carry information about the radial oscillation frequency.
Further, the low frequency peak in our simulations is dominated by the
contribution of the plunge and the first 1-2 bounces. The amplitude of the
radial oscillations depends on the initial NS spin, which therefore has a
complicated influence on the spectrum. Another important result is that HMNSs
can consist of a slowly rotating core with an extended, massive envelope
rotating close to Keplerian velocity, contrary to the common notion that a
rapidly rotating core is necessary to prevent a prompt collapse. Finally, our
estimates on the amount of unbound matter show a dependency on the initial NS
spin, explained by the influence of the latter on the amplitude of radial
oscillations, which in turn cause shock waves.

[9] 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$.

[10] 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.

[11] 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.

[12] 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.

[13] oai:arXiv.org:1305.5720 [pdf] - 667680

The Gravitational Universe

Comments: 20 pages; submitted to the European Space Agency on May 24th, 2013
for the L2/L3 selection of ESA's Cosmic Vision program

Submitted: **2013-05-24**

The last century has seen enormous progress in our understanding of the
Universe. We know the life cycles of stars, the structure of galaxies, the
remnants of the big bang, and have a general understanding of how the Universe
evolved. We have come remarkably far using electromagnetic radiation as our
tool for observing the Universe. However, gravity is the engine behind many of
the processes in the Universe, and much of its action is dark. Opening a
gravitational window on the Universe will let us go further than any
alternative. Gravity has its own messenger: Gravitational waves, ripples in the
fabric of spacetime. They travel essentially undisturbed and let us peer deep
into the formation of the first seed black holes, exploring redshifts as large
as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and
unprecedented measurements of black hole masses and spins will make it possible
to trace the history of black holes across all stages of galaxy evolution, and
at the same time constrain any deviation from the Kerr metric of General
Relativity. eLISA will be the first ever mission to study the entire Universe
with gravitational waves. eLISA is an all-sky monitor and will offer a wide
view of a dynamic cosmos using gravitational waves as new and unique messengers
to unveil The Gravitational Universe. It provides the closest ever view of the
early processes at TeV energies, has guaranteed sources in the form of
verification binaries in the Milky Way, and can probe the entire Universe, from
its smallest scales around singularities and black holes, all the way to
cosmological dimensions.

[14] oai:arXiv.org:1303.1258 [pdf] - 647486

Observation of TeV gamma-rays from the unidentified source HESS
J1841-055 with the ARGO-YBJ experiment

Comments: 17 pages, 4 figures, have been accepted for publication in ApJ

Submitted: **2013-03-06**

We report the observation of a very high energy \gamma-ray source, whose
position is coincident with HESS J1841-055. This source has been observed for
4.5 years by the ARGO-YBJ experiment from November 2007 to July 2012. Its
emission is detected with a statistical significance of 5.3 standard
deviations. Parameterizing the source shape with a two-dimensional Gaussian
function we estimate an extension \sigma=(0.40(+0.32,-0.22}) degree, consistent
with the HESS measurement. The observed energy spectrum is dN/dE =(9.0-+1.6) x
10^{-13}(E/5 TeV)^{-2.32-+0.23} photons cm^{-2} s^{-1} TeV^{-1}, in the energy
range 0.9-50 TeV. The integral \gamma-ray flux above 1 TeV is 1.3-+0.4 Crab
units, which is 3.2-+1.0 times the flux derived by HESS. The differences in the
flux determination between HESS and ARGO-YBJ, and possible counterparts at
other wavelengths are discussed.

[15] oai:arXiv.org:1209.0534 [pdf] - 565311

Long-term Monitoring on Mrk 501 for Its VHE gamma Emission and a Flare
in October 2011

Comments: have been accepted for publication at ApJ

Submitted: **2012-09-04**

As one of the brightest active blazars in both X-ray and very high energy
$\gamma$-ray bands, Mrk 501 is very useful for physics associated with jets
from AGNs. The ARGO-YBJ experiment is monitoring it for $\gamma$-rays above 0.3
TeV since November 2007. Starting from October 2011 the largest flare since
2005 is observed, which lasts to about April 2012. In this paper, a detailed
analysis is reported. During the brightest $\gamma$-rays flaring episodes from
October 17 to November 22, 2011, an excess of the event rate over 6 $\sigma$ is
detected by ARGO-YBJ in the direction of Mrk 501, corresponding to an increase
of the $\gamma$-ray flux above 1 TeV by a factor of 6.6$\pm$2.2 from its steady
emission. In particular, the $\gamma$-ray flux above 8 TeV is detected with a
significance better than 4 $\sigma$. Based on time-dependent synchrotron
self-Compton (SSC) processes, the broad-band energy spectrum is interpreted as
the emission from an electron energy distribution parameterized with a single
power-law function with an exponential cutoff at its high energy end. The
average spectral energy distribution for the steady emission is well described
by this simple one-zone SSC model. However, the detection of $\gamma$-rays
above 8 TeV during the flare challenges this model due to the hardness of the
spectra. Correlations between X-rays and $\gamma$-rays are also investigated.

[16] oai:arXiv.org:1201.1973 [pdf] - 461672

Observation of TeV gamma rays from the Cygnus region with the ARGO-YBJ
experiment

Comments: 14 pages, 4 figures

Submitted: **2012-01-10**

We report the observation of TeV gamma-rays from the Cygnus region using the
ARGO-YBJ data collected from 2007 November to 2011 August. Several TeV sources
are located in this region including the two bright extended MGRO J2019+37 and
MGRO J2031+41. According to the Milagro data set, at 20 TeV MGRO J2019+37 is
the most significant source apart from the Crab Nebula. No signal from MGRO
J2019+37 is detected by the ARGO-YBJ experiment, and the derived flux upper
limits at 90% confidence level for all the events above 600 GeV with medium
energy of 3 TeV are lower than the Milagro flux, implying that the source might
be variable and hard to be identified as a pulsar wind nebula. The only
statistically significant (6.4 standard deviations) gamma-ray signal is found
from MGRO J2031+41, with a flux consistent with the measurement by Milagro.

[17] oai:arXiv.org:1106.0896 [pdf] - 966917

Long-term monitoring of the TeV emission from Mrk 421 with the ARGO-YBJ
experiment

Comments: 30 pages, 8 figures

Submitted: **2011-06-05**

ARGO-YBJ is an air shower detector array with a fully covered layer of
resistive plate chambers. It is operated with a high duty cycle and a large
field of view. It continuously monitors the northern sky at energies above 0.3
TeV. In this paper, we report a long-term monitoring of Mrk 421 over the period
from 2007 November to 2010 February. This source was observed by the
satellite-borne experiments Rossi X-ray Timing Explorer and Swift in the X-ray
band. Mrk 421 was especially active in the first half of 2008. Many flares are
observed in both X-ray and gamma-ray bands simultaneously. The gamma-ray flux
observed by ARGO-YBJ has a clear correlation with the X-ray flux. No lag
between the X-ray and gamma-ray photons longer than 1 day is found. The
evolution of the spectral energy distribution is investigated by measuring
spectral indices at four different flux levels. Hardening of the spectra is
observed in both X-ray and gamma-ray bands. The gamma-ray flux increases
quadratically with the simultaneously measured X-ray flux. All these
observational results strongly favor the synchrotron self-Compton process as
the underlying radiative mechanism.

[18] oai:arXiv.org:1101.4261 [pdf] - 1051613

Mean Interplanetary Magnetic Field Measurement Using the ARGO-YBJ
Experiment

Comments: 6 papges,3 figures

Submitted: **2011-01-21**

The sun blocks cosmic ray particles from outside the solar system, forming a
detectable shadow in the sky map of cosmic rays detected by the ARGO-YBJ
experiment in Tibet. Because the cosmic ray particles are positive charged, the
magnetic field between the sun and the earth deflects them from straight
trajectories and results in a shift of the shadow from the true location of the
sun. Here we show that the shift measures the intensity of the field which is
transported by the solar wind from the sun to the earth.

[19] oai:arXiv.org:1101.2664 [pdf] - 1051412

Differentially-rotating neutron star models with a parametrized rotation
profile

Comments: Submitted to A&A, 6 pages, 3 figures

Submitted: **2011-01-13**

We analyze the impact of the choice rotation law on equilibrium sequences of
relativistic differentially-rotating neutron stars in axisymmetry. The maximum
allowed mass for each model is strongly affected by the distribution of angular
velocity along the radial direction and by the consequent degree of
differential rotation. In order to study the wide parameter space implied by
the choice of rotation law, we introduce a functional form that generalizes the
so called "j-const. law" adopted in all previous work. Using this new rotation
law we reproduce the angular velocity profile of differentially-rotating
remnants from the coalescence of binary neutron stars in various 3-dimensional
dynamical simulations. We compute equilibrium sequences of differentially
rotating stars with a polytropic equation of state starting from the
spherically symmetric static case. By analyzing the sequences at constant
ratio, T/|W|, of rotational kinetic energy to gravitational binding energy, we
find that the parameters that best describe the binary neutron star remnants
cannot produce equilibrium configurations with values of T/|W| that exceed
0.14, the criterion for the onset of the secular instability.