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.

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.

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.

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

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.

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.

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.

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.

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.