Gravitational-wave (GW) detections of binary neutron star coalescences play a crucial role to constrain the microscopic interaction of matter at ultrahigh density. Similarly, if boson stars exist in the universe their coalescence can be used to constrain the fundamental coupling constants of a scalar field theory. We develop the first coherent waveform model for the inspiral of boson stars with quartic interactions. The waveform includes coherently spin-induced quadrupolar and tidal-deformability contributions in terms of the masses and spins of the binary and of a single coupling constant of the theory. We show that future instruments such as the Einstein Telescope and LISA can provide strong, complementary bounds on bosonic self-interactions, while the constraining power of current detectors is marginal.

We develop a new method to measure neutron star parameters and derive constraints on the equation of state of dense matter by fitting the frequencies of simultaneous Quasi Periodic Oscillation modes observed in the X-ray flux of accreting neutron stars in low mass X-ray binaries. To this aim we calculate the fundamental frequencies of geodesic motion around rotating neutron stars based on an accurate general-relativistic approximation for their external spacetime. Once the fundamental frequencies are related to the observed frequencies through a QPO model, they can be fit to the data to obtain estimates of the three parameters describing the spacetime, namely the neutron star mass, angular momentum and quadrupole moment. From these parameters we derive information on the neutron star structure and equation of state. We present a proof of principle of our method applied to pairs of kHz QPO frequencies observed from three systems (4U1608-52, 4U0614+09 and 4U1728-34). We identify the kHz QPOs with the azimuthal and the periastron precession frequencies of matter orbiting the neutron star, and via our Bayesian inference technique we derive constraints on the neutrons stars' masses and radii. This method is applicable to other geodesic-frequency-based QPO models.

Extreme mass-ratio inspirals detectable by the future Laser Interferometer Space Antenna provide a unique way to test general relativity and fundamental physics. Motivated by this possibility, here we study in detail the EMRI dynamics in the presence of a spinning secondary, collecting and extending various results that appeared in previous work and also providing useful intermediate steps and new relations for the first time. We present the results of a frequency-domain code that computes gravitational-wave fluxes and the adiabatic orbital evolution for the case of circular, equatorial orbits with (anti)aligned spins. The spin of the secondary starts affecting the gravitational-wave phase at the same post-adiabatic order as the leading-order self-force terms and introduces a detectable dephasing, which can be used to measure it at $5-25\%$ level, depending on individual spins. In a companion paper we discuss the implication of this effect for tests of the Kerr bound.

An outstanding prediction of general relativity is the fact that the angular momentum $S$ of an isolated black hole with mass $\mu$ is limited by the Kerr bound, $S\leq G\mu^2/c$. Testing this cornerstone is challenging due to the difficulty in modelling spinning compact objects that violate this bound. We argue that precise, model-independent tests can be achieved by measuring gravitational waves from an extreme mass ratio inspiral around a supermassive object, one of the main targets of the future LISA mission. In the extreme mass ratio limit, the dynamics of the small compact object depends only on its multipole moments, which are free parameters. At variance with the comparable-mass case, accurate waveforms are valid also when the spin of the small object greatly exceeds the Kerr bound. By computing the orbital dephasing and the gravitational-wave signal emitted by a spinning point particle in circular, nonprecessing, equatorial motion around a Kerr black hole, we estimate that LISA will be able to measure the spin of the small compact object at the level of $10\%$. Together with mass measurements, this will allow for theory-agnostic, unprecedented constraints on string-theory inspired objects such as "superspinars", almost in their entire parameter space.

Gravitational waves from the merger of two neutron stars cannot be easily distinguished from those produced by a comparable-mass mixed binary in which one of the companions is a black hole. Low-mass black holes are interesting because they could form in the aftermath of the coalescence of two neutron stars, from the collapse of massive stars, from matter overdensities in the primordial Universe, or as the outcome of the interaction between neutron stars and dark matter. Gravitational waves carry the imprint of the internal composition of neutron stars via the so-called tidal deformability parameter, which depends on the stellar equation of state and is equal to zero for black holes. We present a new data analysis strategy powered by Bayesian inference and machine learning to identify mixed binaries, hence low-mass black holes, using the distribution of the tidal deformability parameter inferred from gravitational-wave observations.

We study Extreme Mass Ratio Inspirals (EMRIs), during which a small body spirals into a supermassive black hole, in gravity theories with additional scalar fields. We first argue that no-hair theorems and the properties of known theories that manage to circumvent them introduce a drastic simplification to the problem: the effects of the scalar on supermassive black holes, if any, are mostly negligible for EMRIs in vast classes of theories. We then exploit this simplification to model the inspiral perturbatively and we demonstrate that the scalar charge of the small body leaves a significant imprint on gravitational wave emission. Although much higher precision is needed for waveform modelling, our results strongly suggest that this imprint is observable with LISA, rendering EMRIs a promising probe of scalar fields.

Asymmetric Dark Stars, i.e., compact objects formed from the collapse of asymmetric dark matter could potentially produce a detectable photon flux if dark matter particles self-interact via dark photons that kinetically mix with ordinary photons. The morphology of the emitted spectrum is significantly different and therefore distinguishable from a typical black-body one. Given the above and the fact that asymmetric dark stars can have masses outside the range of neutron stars, the detection of such a spectrum can be considered as a smoking gun signature for the existence of these exotic stars.

Coalescences of binary white dwarfs represent a copious source of information for gravitational wave interferometers operating in the decihertz band. Moreover, according to the double degenerate scenario, they have been suggested to be possible progenitors of supernovae (SNe) Type Ia events. In this paper we discuss the detectability of gravitational waves emitted by the inspiral of double white dwarfs. We focus on the constraints that can be derived on the source's luminosity distance, and on other binary's parameters, such as the angular momentum orientation. We explore the possibility of coincident detections of gravitational and electromagnetic signals; the latter comes from the observation of the supernova counterpart. Confirmation of the double degenerate scenario would allow one to use distances inferred in the gravitational wave channel to consistently calibrate SNe as standard candles. We find that decihertz gravitational wave interferometers can measure the luminosity distance with relative accuracy better than $1\%$ for binaries at 100 Mpc. We show how multimessenger observations can put strong constraints on the Hubble constant, which are tighter than current bounds at low redshift, and how they can potentially shed new light on the differences with early-universe measurements.

Black-hole spectroscopy is arguably the most promising tool to test gravity in extreme regimes and to probe the ultimate nature of black holes with unparalleled precision. These tests are currently limited by the lack of a ringdown parametrization that is both robust and accurate. We develop an observable-based parametrization of the ringdown of spinning black holes beyond general relativity, which we dub ParSpec (Parametrized Ringdown Spin Expansion Coefficients). This approach is perturbative in the spin, but it can be made arbitrarily precise (at least in principle) through a high-order expansion. It requires O(10) ringdown detections, which should be routinely available with the planned space mission LISA and with third-generation ground-based detectors. We provide a preliminary analysis of the projected bounds on parametrized ringdown parameters with LISA and with the Einstein Telescope, and discuss extensions of our model that can be straightforwardly included in the future.

Using a Bayesian approach, we combine measurements of neutron star macroscopic observables obtained by astrophysical and gravitational observations, to derive joint constraints on the equation of state (EoS) of matter at supranuclear density. In our analysis we use two sets of data: (i) the masses and tidal deformabilities measured in the binary neutron star event GW170817, detected by LIGO and Virgo; (ii) the masses and stellar radii measured from observations of nuclear bursts in accreting low-mass X-ray binaries. Using two different parametrizations of the equation of state, we compute the posteriorprobability distributions of the EoS parameters, using which we infer the posterior distribution for the radius and the mass of the two neutron stars of GW170817. The constraints we set on the radii are tighter than previous bounds.

We study the ability of gravitational-wave detectors to constrain the nature of the environment in which compact binaries evolve. We show that the strong dephasing induced by accretion and dynamical friction can constraint the density of the surrounding medium to orders of magnitude below that of accretion disks. Planned detectors, such as LISA or DECIGO, will be able to probe densities typical of those of dark matter.

Future gravitational-wave observations will enable unprecedented and unique science in extreme gravity and fundamental physics answering questions about the nature of dynamical spacetimes, the nature of dark matter and the nature of compact objects.

Linear perturbations of spherically symmetric spacetimes in general relativity are described by radial wave equations, with potentials that depend on the spin of the perturbing field. In previous work we studied the quasinormal mode spectrum of spacetimes for which the radial potentials are slightly modified from their general relativistic form, writing generic small modifications as a power-series expansion in the radial coordinate. We assumed that the perturbations in the quasinormal frequencies are linear in some perturbative parameter, and that there is no coupling between the perturbation equations. In general, matter fields and modifications to the gravitational field equations lead to coupled wave equations. Here we extend our previous analysis in two important ways: we study second-order corrections in the perturbative parameter, and we address the more complex (and realistic) case of coupled wave equations. We highlight the special nature of coupling-induced corrections when two of the wave equations have degenerate spectra, and we provide a ready-to-use recipe to compute quasinormal modes. We illustrate the power of our parametrization by applying it to various examples, including dynamical Chern-Simons gravity, Horndeski gravity and an effective field theory-inspired model.

Supermassive binaries detectable by the future space gravitational-wave interferometer LISA might allow to distinguish black holes from ultracompact horizonless objects, even when the latter are motivated by quantum-gravity considerations. We show that a measurement of very small tidal Love numbers at the level of $10\%$ accuracy (as achievable with "golden binaries") may also allow to distinguish between different models of these exotic compact objects, even when taking into account an intrinsic uncertainty in the object radius putatively due to quantum mechanics. We argue that there is no conceptual obstacle in performing these measurements, the main challenge remains the detectability of small tidal effects and an accurate waveform modelling. Our analysis uses only coordinate-independent quantities related to the proper radial distance and the total mass of the object.

The tidal deformability of a self-gravitating object leaves an imprint on the gravitational-wave signal of an inspiral which is paramount to measure the internal structure of the binary components. We unveil here a surprisingly unnoticed effect: in the extreme-mass ratio limit the tidal Love number of the central object (i.e. the quadrupole moment induced by the tidal field of its companion) affects the gravitational waveform at the leading order in the mass ratio. This effect acts as a magnifying glass for the tidal deformability of supermassive objects but was so far neglected, probably because the tidal Love numbers of a black hole (the most natural candidate for a compact supermassive object) are identically zero. We argue that extreme-mass ratio inspirals detectable by the future LISA mission might place constraints on the tidal Love numbers of the central object which are roughly 8 orders of magnitude more stringent than current ones on neutron stars, potentially probing all models of black hole mimickers proposed so far.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

The High Energy Modular Ensemble of Satellites (HERMES) project is aimed to realize a modular X/gamma-ray monitor for transient events, to be placed on-board of a CubeSat bus. This expandable platform will achieve a significant impact on Gamma Ray Burst (GRB) science and on the detection of Gravitational Wave (GW) electromagnetic counterparts: the recent LIGO/VIRGO discoveries demonstrated that the high-energy transient sky is still a field of extreme interest. The very complex temporal variability of GRBs (up to the millisecond scale) combined with the spatial and temporal coincidence between GWs and their electromagnetic counterparts suggest that upcoming instruments require sub-ms time resolution combined with a transient localization accuracy lower than a degree. The current phase of the ongoing HERMES project is focused on the realization of a technological pathfinder with a small network (3 units) of nano-satellites to be launched in mid 2020. We will show the potential and prospects for short and medium-term development of the project, demonstrating the disrupting possibilities for scientific investigations provided by the innovative concept of a new "modular astronomy" with nano-satellites (e.g. low developing costs, very short realization time). Finally, we will illustrate the characteristics of the HERMES Technological Pathfinder project, demonstrating how the scientific goals discussed are actually already reachable with the first nano-satellites of this constellation. The detector architecture will be described in detail, showing that the new generation of scintillators (e.g. GAGG:Ce) coupled with very performing Silicon Drift Detectors (SDD) and low noise Front-End-Electronics (FEE) are able to extend down to few keV the sensitivity band of the detector. The technical solutions for FEE, Back-End-Electronics (BEE) and Data Handling will be also described.

We investigated the X-ray emission properties of the powerful radio galaxy 3C 459 revealed by a recent Chandra follow-up observation carried out in October 2014 with a 62 ks exposure. We performed an X-ray spectral analysis from a few selected regions on an image obtained from this observation and also compared the X-ray image with a 4.9 GHz VLA radio map available in the literature. The dominant contribution comes from the radio core but significant X-ray emission is detected at larger angular separations from it, surrounding both radio jets and lobes. According to a scenario in which the extended X-ray emission is due to a plasma collisionally heated by jet-driven shocks and not magnetically dominated, we estimated its temperature to be ~0.8 keV. This hot gas cocoon could be responsible for the radio depolarization observed in 3C 459, as recently proposed also for 3C 171 and 3C 305. On the other hand, our spectral analysis and the presence of an oxygen K edge, blueshifted at 1.23 keV, cannot exclude the possibility that the X-ray radiation originating from the inner regions of the radio galaxy could be intercepted by some outflow of absorbing material intervening along the line of sight, as already found in some BAL quasars.

Effective field theory methods suggest that some rather-general extensions of General Relativity include, or are mimicked by, certain higher-order curvature corrections, with coupling constants expected to be small but otherwise arbitrary. Thus, the tantalizing prospect to test the fundamental nature of gravity with gravitational-wave observations, in a systematic way, emerges naturally. Here, we build black hole solutions in such a framework and study their main properties. Once rotation is included, we find the first purely gravitational example of geometries without $\mathbb{Z}_2$-symmetry. Despite the higher-order operators of the theory, we show that linearized fluctuations of such geometries obey second-order differential equations. We find nonzero tidal Love numbers. We study and compute the quasinormal modes of such geometries. These results are of interest to gravitational-wave science but also potentially relevant for electromagnetic observations of the galactic center or $X$-ray binaries.

Here we investigate how LyC-opaque systems present in the intergalactic medium at $z\approx3$ can distort the spectral shape of a uniform UV background (UVB) through radiative transfer (RT) effects. With this aim in mind, we perform a multi-frequency RT simulation through a cosmic volume of $10h^{-1}$~cMpc scale polluted by metals, and self-consistently derive the ions of all the species. The UVB spatial fluctuations are traced by the ratio of He$\, \rm \scriptstyle II\ $ and H$\, \rm \scriptstyle I\ $ column density, $\eta$, and the ratio of C$\,{\rm {\scriptstyle IV\ }}$ and Si$\,{\rm {\scriptstyle IV\ }}$ optical depths, $\zeta$. We find that: (i) $\eta$ spatially fluctuates through over-dense systems ($\Delta$) with statistically significant deviations $\delta\eta >25$\% in 18\% of the volume ; (ii) same fluctuations in $\zeta$ are also present in $34$\% of the enriched domain (only 8\% of the total volume) and derive from a combination of RT induced effects and in-homogeneous metal enrichment, both effective in systems with $\Delta > 1.5$.

The evolutionary status of the low mass X-ray binary SAX J1808.4-3658 is simulated by following the binary evolution of its possible progenitor system through mass transfer, starting at a period of $\sim$6.6 hr. The evolution includes angular momentum losses via magnetic braking and gravitational radiation. It also takes into account the effects of illumination of the donor by both the X-ray emission and the spin down luminosity of the pulsar. The system goes through stages of mass transfer and stages during which it is detached, where only the rotationally powered pulsar irradiates the donor. We show that the pulsar irradiation is a necessary ingredient to reach SAX J1808.4-3658 orbital period when the donor mass is reduced to 0.04 - 0.06 M$_\odot$. We also show that the models reproduce important properties of the system, including the orbital period derivative, which is shown to be directly linked to the evolution through mass transfer cycles. Moreover we find that the effects of the irradiation on the internal structure of the donor are non negligible, causing the companion star to be non completely convective at the values of mass observed for the system and significantly altering its long term evolution, as the magnetic braking remains active along the whole evolution.

The LIGO/Virgo collaboration has recently announced the direct detection of gravitational waves emitted in the coalescence of a neutron star binary. This discovery allows, for the first time, to set new constraints on the behavior of matter at supranuclear density, complementary with those coming from astrophysical observations in the electromagnetic band. In this paper we demonstrate the feasibility of using gravitational signals to solve the relativistic inverse stellar problem, i.e. to reconstruct the parameters of the equation of state (EoS) from measurements of the stellar mass and tidal Love number. We perform Bayesian inference of mock data, based on different models of the star internal composition, modeled through piecewise polytropes. Our analysis shows that the detection of a small number of sources by a network of advanced interferometers would allow to put accurate bounds on the EoS parameters, and to perform a model selection among the realistic equations of state proposed in the literature.

Several quantum-gravity models of compact objects predict microscopic or even Planckian corrections at the horizon scale. We explore the possibility of measuring two model-independent, smoking-gun effects of these corrections in the gravitational waveform of a compact binary, namely the absence of tidal heating and the presence of tidal deformability. For events detectable by the future space-based interferometer LISA, we show that the effect of tidal heating dominates and allows one to constrain putative corrections down to the Planck scale. The measurement of the tidal Love numbers with LISA is more challenging but, in optimistic scenarios, it allows to constrain the compactness of a supermassive exotic compact object down to the Planck scale. Our analysis suggests that highly-spinning, supermassive binaries at 1-20 Gpc provide unparalleled tests of quantum-gravity effects at the horizon scale.

Future very-large-area X-ray instruments (for which the effective area is larger than $>3$m$^2$) will be able to measure the frequencies of quasi-periodic oscillations~(QPOs) observed in the X-ray flux from accreting compact objects with sub-percent precision. If correctly modeled, QPOs can provide a novel way to test the strong-field regime of gravity. By using the relativistic precession model and a modified version of the epicyclic resonance model, we develop a method to test general relativity against a generic class of theories with quadratic curvature corrections. With the instrumentation being studied for future missions such as eXTP, LOFT, or STROBE-X, a measurement of at least two QPO triplets from a stellar mass black hole can set stringent constraints on the coupling parameters of quadratic gravity.

Light bosonic fields are ubiquitous in extensions of the Standard Model. Even when minimally coupled to gravity, these fields might evade the assumptions of the black-hole no-hair theorems and give rise to spinning black holes which can be drastically different from the Kerr metric. Furthermore, they allow for self-gravitating compact solitons, known as (scalar or Proca) boson stars. The quasi-periodic oscillations (QPOs) observed in the X-ray flux emitted by accreting compact objects carry information about the strong-field region, thus providing a powerful tool to constrain deviations from Kerr's geometry and to search for exotic compact objects. By using the relativistic precession model as a proxy to interpret the QPOs in terms of geodesic frequencies, we investigate how the QPO frequencies could be used to test the no-hair theorem and the existence of light bosonic fields near accreting compact objects. We show that a detection of two QPO triplets with current sensitivity can already constrain these models, and that the future eXTP mission or a LOFT-like mission can set very stringent constraints on black holes with bosonic hair and on (scalar or Proca) boson stars. The peculiar geodesic structure of compact scalar/Proca boson stars implies that these objects can easily be ruled out as alternative models for X-ray source GRO J1655-40.

Theoretical models of self-interacting dark matter represent a promising answer to a series of open problems within the so-called collisionless cold dark matter (CCDM) paradigm. In case of asymmetric dark matter, self-interactions might facilitate gravitational collapse and potentially lead to formation of compact objects predominantly made of dark matter. Considering both fermionic and bosonic equations of state, we construct the equilibrium structure of rotating dark stars, focusing on their bulk properties, and comparing them with baryonic neutron stars. We also show that these dark objects admit the $I$-Love-$Q$ universal relations, which link their moments of inertia, tidal deformabilities, and quadrupole moments. Finally, we prove that stars built with a dark matter equation of state are not compact enough to mimic black holes in general relativity, thus making them distinguishable in potential events of gravitational interferometers.

The tidal Love numbers (TLNs) encode the deformability of a self-gravitating object immersed in a tidal environment and depend significantly both on the object's internal structure and on the dynamics of the gravitational field. An intriguing result in classical general relativity is the vanishing of the TLNs of black holes. We extend this result in three ways, aiming at testing the nature of compact objects: (i) we compute the TLNs of exotic compact objects, including different families of boson stars, gravastars, wormholes, and other toy models for quantum corrections at the horizon scale. In the black-hole limit, we find a universal logarithmic dependence of the TLNs on the location of the surface; (ii) we compute the TLNs of black holes beyond vacuum general relativity, including Einstein-Maxwell, Brans-Dicke and Chern-Simons gravity; (iii) We assess the ability of present and future gravitational-wave detectors to measure the TLNs of these objects, including the first analysis of TLNs with LISA. Both LIGO, ET and LISA can impose interesting constraints on boson stars, while LISA is able to probe even extremely compact objects. We argue that the TLNs provide a smoking gun of new physics at the horizon scale, and that future gravitational-wave measurements of the TLNs in a binary inspiral provide a novel way to test black holes and general relativity in the strong-field regime.

We report the detection of two new gamma-ray sources in the Fermi-LAT sky (Pass 8) at energies higher than 20 GeV, and confirmed at lower energies, using a source detection tool based on the Minimum Spanning Tree algorithm. One of these sources, at a Galactic latitude of about -4{\deg}, is a new discovery, while the other was previously reported above 50 GeV in the 2FHL catalogue. We searched for archival multi-wavelength data of possible counterparts and found interesting candidates. Both objects are radio sources and their WISE infrared colours are typical of blazars. While for the former source no optical spectra are available, for the latter a puzzling optical spectrum corresponding to a white dwarf star is found in the 6dF database. We discuss the spectral energy distributions of both sources and possible interpretations.

We have investigated a group of unassociated radio sources included in the 3CR cat- alogue to increase the multi-frequency information on them and possibly obtain an identification. We have carried out an observational campaign with the Swift satellite to observe with the UVOT and the XRT telescopes the field of view of 21 bright NVSS sources within the positional uncertainty region of the 3CR sources. Furthermore, we have searched in the recent AllWISE Source Catalogue for infrared sources matching the position of these NVSS sources. We have detected significant emission in the soft X-ray band for nine of the investigated NVSS sources. To all of them, and in four cases with no soft X-ray association, we have associated a WISE infrared counterpart. Eight of these infrared candidates have not been proposed earlier in the literature. In the five remaining cases our candidate matches one among a few optical candidates suggested for the same 3CR source in previous studies. No source has been detected in the UVOT filters at the position of the NVSS objects, confirming the scenario that all of them are heavily obscured. With this in mind, a spectroscopic campaign, preferably in the infrared band, will be necessary to establish the nature of the sources that we have finally identified.

We report on the detection of a bright, short, structured X-ray burst coming from the supernova remnant RCW 103 on 2016 June 22 caught by the Swift/BAT monitor, and on the follow-up campaign made with Swift/XRT, Swift/UVOT and the optical/NIR GROND detector. The characteristics of this flash, such as duration, and spectral shape, are consistent with typical short bursts observed from soft gamma repeaters. The BAT error circle at 68 per cent confidence range encloses the point-like X-ray source at the centre of the nebula, 1E161348-5055. Its nature has been long debated due to a periodicity of 6.67 hr in X-rays, which could indicate either an extremely slow pulsating neutron star, or the orbital period of a very compact X-ray binary system. We found that 20 min before the BAT trigger, the soft X-ray emission of 1E161348-5055 was a factor of ~100 higher than measured 2 yr earlier, indicating that an outburst had already started. By comparing the spectral and timing characteristics of the source in the two years before the outburst and after the BAT event, we find that, besides a change in luminosity and spectral shape, also the 6.67 hr pulsed profile has significantly changed with a clear phase shift with respect to its low-flux profile. The UV/optical/NIR observations did not reveal any counterpart at the position of 1E161348-5055. Based on these findings, we associate the BAT burst with 1E161348-5055, we classify it as a magnetar, and pinpoint the 6.67 hr periodicity as the magnetar spin period.

The direct discovery of gravitational waves has finally opened a new observational window on our Universe, suggesting that the population of coalescing binary black holes is larger than previously expected. These sources produce an unresolved background of gravitational waves, potentially observables by ground-based interferometers. In this paper we investigate how modified theories of gravity, modeled using the ppE formalism, affect the expected signal, and analyze the detectability of the resulting stochastic background by current and future ground-based interferometers. We find the constraints that AdLIGO would be able to set on modified theories, showing that they may significantly improve the current bounds obtained from astrophysical observations of binary pulsars.

Horndeski gravity holds a special position as the most general extension of Einstein's theory of general relativity with a single scalar degree of freedom and second-order field equations. Because of these features, Horndeski gravity is an attractive phenomenological playground to investigate the consequences of modifications of general relativity in cosmology and astrophysics. We present a review of the progress made so far in the study of compact objects (black holes and neutron stars) within Horndeski gravity. In particular, we review our recent work on slowly rotating black holes and present some new results on slowly rotating neutron stars.

Horndeski's theory of gravity is the most general scalar-tensor theory with a single scalar whose equations of motion contain at most second-order derivatives. A subsector of Horndeski's theory known as "Fab Four" gravity allows for dynamical self-tuning of the quantum vacuum energy, and therefore it has received particular attention in cosmology as a possible alternative to the $\Lambda$CDM model. Here we study compact stars in Fab Four gravity, which includes as special cases general relativity ("George"), Einstein-dilaton-Gauss-Bonnet gravity ("Ringo"), theories with a nonminimal coupling with the Einstein tensor ("John") and theories involving the double-dual of the Riemann tensor ("Paul"). We generalize and extend previous results in theories of the John class and were not able to find realistic compact stars in theories involving the Paul class.

Astrophysical phenomena involving massive black holes (BHs) in close binaries are expected to leave detectable signatures in the electromagnetic and gravitational-wave spectrum. Such imprints may provide precious information to probe the space-time around rotating BHs and to reveal new insights on the nature of gravity in the strong-field regime. To support this observational window, it is crucial to develop suitable tests to verify the predictions of General Relativity. In this framework, the metric recently proposed by Johannsen and Psaltis parametrizes strong-field deviations from a Kerr space-time in a theory-independent way. In the following, we make use of this approach to describe the tidal field produced by spinning BHs. We compute the gravito-magnetic and gravito-electric tidal tensors for particles moving on equatorial circular geodesics, comparing our results with those obtained in the standard General Relativity scenario. Our calculations show significant differences even for distances far form the last stable orbit, which may affect the evolution of the binary and leave detectable signatures. We test our framework computing quasiequilibrium sequences of BH-white dwarf systems by means of the affine model, for different binary configurations.

We study black hole solutions at first order in the Hartle-Thorne slow-rotation approximation in Horndeski gravity theories. We derive the equations of motion including also cases where the scalar depends linearly on time. In the Hartle-Thorne formalism, all first-order rotational corrections are described by a single frame-dragging function. We show that the frame-dragging function is exactly the same as in general relativity for all known black hole solutions in shift symmetric Horndeski theories, with the exception of theories with a linear coupling to the Gauss-Bonnet invariant. Our results extend previous no-hair theorems for a broad class of Horndeski gravity theories.

We investigate the feasibility of detecting 21cm absorption features in the afterglow spectra of high redshift long Gamma Ray Bursts (GRBs). This is done employing simulations of cosmic reionization, together with the instrumental characteristics of the LOw Frequency ARray (LOFAR). We find that absorption features could be marginally (with a S/N larger than a few) detected by LOFAR at z>7 if the GRB originated from PopIII stars, while the detection would be easier if the noise were reduced by one order of magnitude, i.e. similar to what is expected for the first phase of the Square Kilometer Array (SKA1-low). On the other hand, more standard GRBs are too dim to be detected even with ten times the sensitivity of SKA1-low, and only in the most optimistic case can a S/N larger than a few be reached at z>9.

The deformability of a compact object induced by a perturbing tidal field is encoded in the tidal Love numbers, which depend sensibly on the object's internal structure. These numbers are known only for static, spherically-symmetric objects. As a first step to compute the tidal Love numbers of a spinning compact star, here we extend powerful perturbative techniques to compute the exterior geometry of a spinning object distorted by an axisymmetric tidal field to second order in the angular momentum. The spin of the object introduces couplings between electric and magnetic deformations and new classes of induced Love numbers emerge. For example, a spinning object immersed in a quadrupolar, electric tidal field can acquire some induced mass, spin, quadrupole, octupole and hexadecapole moments to second order in the spin. The deformations are encoded in a set of inhomogeneous differential equations which, remarkably, can be solved analytically in vacuum. We discuss certain subtleties in defining the multipole moments of the central object, which are due to the difficulty in separating the tidal field from the linear response of the object in the solution. By extending the standard procedure to identify the linear response in the static case, we prove analytically that the Love numbers of a Kerr black hole remain zero to second order in the spin. As a by-product, we provide the explicit form for a slowly-rotating, tidally-deformed Kerr black hole to quadratic order in the spin, and discuss its geodesic and geometrical properties.

We present a nonlinear post-Friedmann framework for structure formation, generalizing to cosmology the weak-field (post-Minkowskian) approximation, unifying the treatment of small and large scales. We consider a universe filled with a pressureless fluid and a cosmological constant $\Lambda$, the theory of gravity is Einstein's general relativity and the background is the standard flat $\Lambda$CDM cosmological model. We expand the metric and the energy-momentum tensor in powers of $1/c$, keeping the matter density and peculiar velocity as exact fundamental variables. We assume the Poisson gauge, including scalar and tensor modes up to $1/c^4$ order and vector modes up to $1/c^5$ terms. Through a redefinition of the scalar potentials as a resummation of the metric contributions at different orders, we obtain a complete set of nonlinear equations, providing a unified framework to study structure formation from small to superhorizon scales, from the nonlinear Newtonian to the linear relativistic regime. We explicitly show the validity of our scheme in the two limits: at leading order we recover the fully nonlinear equations of Newtonian cosmology; when linearized, our equations become those for scalar and vector modes of first-order relativistic perturbation theory in the Poisson gauge. Tensor modes are nondynamical at the $1/c^4$ order we consider: they are purely nonlinear and describe a distortion of the spatial slices determined at this order by a constraint, quadratic in the scalar and vector variables. The main results of our analysis are as follows: (a) at leading order a purely Newtonian nonlinear energy current sources a frame-dragging gravitomagnetic vector potential, and (b) in the leading-order Newtonian regime and in the linear relativistic regime the two scalar metric potentials are the same, while the nonlinearity of general relativity makes them different.

Blazars are radio-loud active galactic nuclei well known for their non thermal emission spanning a wide range of frequencies. The Roma-BZCAT is, to date, the most comprehensive list of these sources. We performed the cross-match of several catalogs obtained from recent surveys at different frequencies to search for new blazars. We cross-matched the 1$^{st}$ Swift-XRT Point Source catalog with the spectroscopic sample of the 9$^{th}$ Data Release of the Sloan Digital Sky Survey. Then, we performed further cross-matches with the catalogs corresponding to the Faint Images of the Radio Sky at Twenty cm survey and to the AllWISE Data release, focusing on sources with infrared colors similar to those of confirmed $\gamma$-ray blazars included in the Second Fermi-LAT catalog. As a result, we obtained a preliminary list of objects with all the elements needed for a proper blazar classification according to the prescriptions of the Roma-BZCAT. We carefully investigated additional properties such as their morphology and the slope of their spectral energy distribution in the radio domain, the features shown in their optical spectrum, and the luminosity in the soft X rays to exclude generic active galactic nuclei and focus on authentic blazar-like sources. At the end of our screening we obtained a list of 15 objects with firmly established blazar properties.

The 5th edition of the Roma-BZCAT Multifrequency Catalogue of Blazars is available in a printed version and online at the ASDC website (http://www.asdc.asi.it/bzcat); it is also in the NED database. It presents several relevant changes with respect to the past editions which are briefly described in this paper.

Here we introduce GAMESH, a novel pipeline which implements self-consistent radiative and chemical feedback in a computational model of galaxy formation. By combining the cosmological chemical-evolution model GAMETE with the radiative transfer code CRASH, GAMESH can post process realistic outputs of a N-body simulation describing the redshift evolution of the forming galaxy. After introducing the GAMESH implementation and its features, we apply the code to a low-resolution N-body simulation of the Milky Way formation and we investigate the combined effects of self-consistent radiative and chemical feedback. Many physical properties, which can be directly compared with observations in the Galaxy and its surrounding satellites, are predicted by the code along the merger-tree assembly. The resulting redshift evolution of the Local Group star formation rates, reionisation and metal enrichment along with the predicted Metallicity Distribution Function of halo stars are critically compared with observations. We discuss the merits and limitations of the first release of GAMESH, also opening new directions to a full implementation of feedback processes in galaxy formation models by combining semi-analytic and numerical methods.

Quasi-Periodic Oscillations (QPOs) observed in the X-ray flux emitted by accreting black holes, are associated to phenomena occurring near the horizon. Future very large area X-ray instruments will be able to measure QPO frequencies with very high precision, thus probing this strong-field region. By using the relativistic precession model, we show the way in which QPO frequencies could be used to test general relativity against those alternative theories of gravity which predict deviations from the classical theory in the strong-field regime. We consider one of the best motivated strong-curvature corrections to general relativity, namely the Einstein-Dilaton-Gauss-Bonnet theory, and show that a detection of QPOs with the expected sensitivity of the proposed ESA M-class mission LOFT would set the most stringent constraints on the parameter space of this theory.

Shortly after its birth in a gravitational collapse, a proto-neutron star enters in a phase of quasi-stationary evolution characterized by large gradients of the thermodynamical variables and intense neutrino emission. In few tens of seconds the gradients smooth out while the star contracts and cools down, until it becomes a neutron star. In this paper we study this phase of the proto-neutron star life including rotation, and employing finite temperature equations of state. We model the evolution of the rotation rate, and determine the relevant quantities characterizing the star. Our results show that an isolated neutron star cannot reach, at the end of the evolution, the maximum values of mass and rotation rate allowed by the zero-temperature equation of state. Moreover, a mature neutron star evolved in isolation cannot rotate too rapidly, even if it is born from a proto-neutron star rotating at the mass-shedding limit. We also show that the I-Love-Q relations are violated in the first second of life, but they are satisfied as soon as the entropy gradients smooth out.

Black-hole neutron-star coalescing binaries have been invoked as one of the most suitable scenario to explain the emission of short gamma-ray bursts. Indeed, if the black-hole which forms after the merger, is surrounded by a massive disk, neutrino annihilation processes may produce high-energy and collimated electromagnetic radiation. In this paper, we devise a new procedure, to be used in the search for gravitational waves from black-hole-neutron-star binaries, to assign a probability that a detected gravitational signal is associated to the formation of an accreting disk, massive enough to power gamma-ray bursts. This method is based on two recently proposed semi-analytic fits, one reproducing the mass of the remnant disk surrounding the black hole as a function of some binary parameters, the second relating the neutron star compactness, with its tidal deformability. Our approach can be used in low-latency data analysis to restrict the parameter space searching for gravitational signals associated with short gamma-ray bursts, and to gain information on the dynamics of the coalescing system and on the neutron star equation of state.

We propose a procedure to be used in the search for gravitational waves from black hole-neutron star coalescing binaries, in coincidence with short gamma-ray bursts. It is based on two recently proposed semi-analytic fits, one reproducing the mass of the remnant disk surrounding the black hole which forms after the merging as a function of some binary parameters, the second relating the neutron star compactness, i.e. the ratio of mass and radius, with its tidal deformability. Using a Fisher matrix analysis and the two fits, we assign a probability that the emitted gravitational signal is associated to the formation of an accreting disk massive enough to supply the energy needed to power a short gamma ray burst. This information can be used in low-latency data analysis to restrict the parameter space searching for gravitational wave signals in coincidence with short gamma-ray bursts, and to gain information on the dynamics of the coalescing system and on the internal structure of the components. In addition, when the binary parameters will be measured with high accuracy, it will be possible to use this information to trigger the search for off-axis gamma-ray bursts afterglows.

Long-duration Gamma-Ray Bursts (GRBs) are an extremely rare outcome of the collapse of massive stars, and are typically found in the distant Universe. Because of its intrinsic luminosity ($L\sim 3 \times 10^{53}$ erg s$^{-1}$) and its relative proximity ($z=0.34$), GRB 130427A was a unique event that reached the highest fluence observed in the gamma-ray band. Here we present a comprehensive multiwavelength view of GRB 130427A with Swift, the 2-m Liverpool and Faulkes telescopes and by other ground-based facilities, highlighting the evolution of the burst emission from the prompt to the afterglow phase. The properties of GRB 130427A are similar to those of the most luminous, high-redshift GRBs, suggesting that a common central engine is responsible for producing GRBs in both the contemporary and the early Universe and over the full range of GRB isotropic energies.

We present the analysis of the extraordinarily bright Gamma-Ray Burst (GRB) 130427A under the hypothesis that the GRB central engine is an accretion-powered magnetar. In this framework, initially proposed to explain GRBs with precursor activity, the prompt emission is produced by accretion of matter onto a newly-born magnetar, and the observed power is related to the accretion rate. The emission is eventually halted if the centrifugal forces are able to pause accretion. We show that the X-ray and optical afterglow is well explained as the forward shock emission with a jet break plus a contribution from the spin-down of the magnetar. Our modelling does not require any contribution from the reverse shock, that may still influence the afterglow light curve at radio and mm frequencies, or in the optical at early times. We derive the magnetic field ($B\sim 10^{16}$ G) and the spin period ($P\sim 20$ ms) of the magnetar and obtain an independent estimate of the minimum luminosity for accretion. This minimum luminosity results well below the prompt emission luminosity of GRB 130427A, providing a strong consistency check for the scenario where the entire prompt emission is the result of continuous accretion onto the magnetar. This is in agreement with the relatively long spin period of the magnetar. GRB 130427A was a well monitored GRB showing a very standard behavior and, thus, is a well-suited benchmark to show that an accretion-powered magnetar gives a unique view of the properties of long GRBs.

We study how to extract information on the neutron star equation of state from the gravitational wave signal emitted during the coalescence of a binary system composed of two neutron stars or a neutron star and a black hole. We use post-Newtonian templates which include the tidal deformability parameter and, when tidal disruption occurs before merger, a frequency cut-off. Assuming that this signal is detected by Advanced LIGO/Virgo or ET, we evaluate the uncertainties on these parameters using different data analysis strategies based on the Fisher matrix approach, and on recently obtained analytical fits of the relevant quantities. We find that the tidal deformability is more effective than the stellar compactness to discriminate among different possible equations of state.

We searched for gamma-ray blazar candidates among the 382 unidentified hard X-ray sources of the 3rd Palermo BAT Catalog (3PBC) obtained from the analysis of 66 months of SWIFT-BAT survey data and listing 1586 sources. We adopted a recently developed association method based on the peculiar infrared colors which characterize the gamma-ray blazars included in the second catalog of active galactic nuclei detected by the Fermi Large Area Telescope (2LAC). We used this method exploiting the data of the all-sky survey performed by the Wide-Field Infrared Survey Explorer (WISE) to establish correspondences between unidentified 3PBC sources and WISE gamma-ray blazar candidates located within the BAT positional uncertainty region at 99% confidence level. We obtained a preliminary list of candidates for which we analysed all the available data in the SWIFT archive to complement the information in the literature and in the radio, infrared and optical catalogs with the information on their optical-UV and soft X-ray emission. Requiring the presence of radio and soft X-ray counterparts consistent with the infrared positions of the selected WISE sources, as well as a blazar-like radio morphology, we finally obtained a list of 24 gamma-ray blazar candidates.

Neutron stars are extremely relativistic objects which abound in our universe and yet are poorly understood, due to the high uncertainty on how matter behaves in the extreme conditions which prevail in the stellar core. It has recently been pointed out that the moment of inertia I, the Love number lambda and the spin-induced quadrupole moment Q of an isolated neutron star, are related through functions which are practically independent of the equation of state. These surprising universal I-lambda-Q relations pave the way for a better understanding of neutron stars, most notably via gravitational-wave emission. Gravitational-wave observations will probe highly-dynamical binaries and it is important to understand whether the universality of the I-lambda-Q relations survives strong-field and finite-size effects. We apply a Post-Newtonian-Affine approach to model tidal deformations in compact binaries and show that the I-lambda relation depends on the inspiral frequency, but is insensitive to the equation of state. We provide a fit for the universal relation, which is valid up to a gravitational wave frequency of ~900 Hz and accurate to within a few percent. Our results strengthen the universality of I-lambda-Q relations, and are relevant for gravitational-wave observations with advanced ground-based interferometers. We also discuss the possibility of using the Love-compactness relation to measure the neutron-star radius with an uncertainty of about 10% or smaller from gravitational-wave observations.

Here we introduce CRASH3, the latest release of the 3D radiative transfer code CRASH. In its current implementation CRASH3 integrates into the reference algorithm the code Cloudy to evaluate the ionisation states of metals, self-consistently with the radiative transfer through H and He. The feedback of the heavy elements on the calculation of the gas temperature is also taken into account, making of CRASH3 the first 3D code for cosmological applications which treats self-consistently the radiative transfer through an inhomogeneous distribution of metal enriched gas with an arbitrary number of point sources and/or a background radiation. The code has been tested in idealized configurations, as well as in a more realistic case of multiple sources embedded in a polluted cosmic web. Through these validation tests the new method has been proven to be numerically stable and convergent. We have studied the dependence of the results on a number of physical quantities such as the source characteristics (spectral range and shape, intensity), the metal composition, the gas number density and metallicity.

We discuss the feasibility of the detection of the 21cm forest in the diffuse IGM with the radio telescope LOFAR. The optical depth to the 21cm line has been derived using simulations of reionization which include detailed radiative transfer of ionizing photons. We find that the spectra from reionization models with similar total comoving hydrogen ionizing emissivity but different frequency distribution look remarkably similar. Thus, unless the reionization histories are very different from each other (e.g. a predominance of UV vs. x-ray heating) we do not expect to distinguish them by means of observations of the 21cm forest. Because the presence of a strong x-ray background would make the detection of 21cm line absorption impossible, the lack of absorption could be used as a probe of the presence/intensity of the x-ray background and the thermal history of the universe. Along a random line of sight LOFAR could detect a global suppression of the spectrum from z>12, when the IGM is still mostly neutral and cold, in contrast with the more well-defined, albeit broad, absorption features visible at lower redshift. Sharp, strong absorption features associated with rare, high density pockets of gas could be detected also at z~7 along preferential lines of sight.

Using a semi-analytical approach recently developed to model the tidal deformations of neutron stars in inspiralling compact binaries, we study the dynamical evolution of the tidal tensor, which we explicitly derive at second post-Newtonian order, and of the quadrupole tensor. Since we do not assume a priori that the quadrupole tensor is proportional to the tidal tensor, i.e. the so called "adiabatic approximation", our approach enables us to establish to which extent such approximation is reliable. We find that the ratio between the quadrupole and tidal tensors (i.e., the Love number) increases as the inspiral progresses, but this phenomenon only marginally affects the emitted gravitational waveform. We estimate the frequency range in which the tidal component of the gravitational signal is well described using the stationary phase approximation at next-to-leading post-Newtonian order, comparing different contributions to the tidal phase. We also derive a semi-analytical expression for the Love number, which reproduces within a few percentage points the results obtained so far by numerical integrations of the relativistic equations of stellar perturbations.

We perform a systematic study of how the inhomogeneities in the Inter-Galactic Medium (IGM) affect the observability of Lyman Alpha Emitters (LAEs) around the Epoch of Reionization. We focus on the IGM close to the galaxies as the detailed ionization distribution and velocity fields of this region could significantly influence the scattering of Ly-alpha photons off neutral H atoms as they traverse the IGM after escaping from the galaxy. We simulate the surface brightness (SB) maps and spectra of more than 100 LAEs at z=7.7 as seen by an observer at z=0. To achieve this, we extract the source properties of galaxies and their surrounding IGM from cosmological simulations of box sizes 5-30 Mpc/h and follow the coupled radiative transfer of ionizing and Ly-alpha radiation through the IGM using CRASH-alpha. We find that the simulated SB profiles are extended and their detailed structure is affected by inhomogeneities in the IGM, especially at high neutral fractions. The detectability of LAEs and the fraction of the flux observed depend heavily on the shape of the SB profile and the SB threshold (SB_th) of the observational campaign. Only ultradeep observations (e.g. SB_th ~ 10^-23 ergs/s/cm^2/arcsec^2) would be able to obtain the true underlying mass-luminosity relation and luminosity functions of LAEs. The details of our results depend on whether Ly-alpha photons are significantly shifted in the galaxy to longer wavelengths, the mean ionization fraction in the IGM and the clustering of ionizing sources. These effects can lead to an easier escape of Ly-alpha photons with less scattering in the IGM and a concentrated SB profile similar to the one of a point source. Finally, we show that the SB profiles are steeper at high ionization fraction for the same LAE sample which can potentially be observed from the stacked profile of a large number of LAEs.

We develop a semi-analytical approach, based on the post-Newtonian expansion and on the affine approximation, to model the tidal deformation of neutron stars in the coalescence of black hole-neutron star or neutron star-neutron star binaries. Our equations describe, in a unified framework, both the system orbital evolution, and the neutron star deformations. These are driven by the tidal tensor, which we expand at 1/c^3 post-Newtonian order, including spin terms. We test the theoretical framework by simulating black hole-neutron star coalescence up to the onset of mass shedding, which we determine by comparing the shape of the star with the Roche lobe. We validate our approach by comparing our results with those of fully relativistic, numerical simulations.

We investigate the effect of helium on hydrogen reionisation using a hydrodynamical simulation combined with the cosmological radiative transfer code CRASH. The simulations are run in a 35.12/h comoving Mpc box using a variety of assumptions for the amplitude and power-law extreme-UV (EUV) spectral index, alpha, of the ionising emissivity. We use an empirically motivated prescription for ionising sources which ensures all of the models are consistent with constraints on the Thomson scattering optical depth and the hydrogen photo-ionisation rate at z=6. The inclusion of helium slightly delays reionisation due to the small number of ionising photons which reionise neutral helium instead of hydrogen. However, helium has a significant impact on the thermal state of the IGM. Models with alpha=3 produce IGM temperatures at the mean density at z=6 which are about 20 % higher compared to models without helium photo-heating. Harder EUV indices produce even larger IGM temperature boosts. A comparison to recent observational estimates of the IGM temperature at z=5 - 6 suggests that hydrogen reionisation was primarily driven by pop-II stellar sources with a soft EUV index, alpha<3. We also find that faint, as yet undetected galaxies, characterised by a luminosity function with a steepening faint-end slope and an increasing Lyman continuum escape fraction (fesc=0.5), are required to reproduce the ionising emissivity used in our simulations at z>6. Finally, we note there is some tension between recent observational constraints which indicate the IGM is > 10% neutral by volume z=7, and estimates of the ionising emissivity at z=6 which indicate only between 1 and 3 ionising photons are emitted per hydrogen atom over a Hubble time. This tension may be alleviated by either a lower neutral fraction at z=7 or an IGM which still remains a few % neutral by volume at z=6.

We report on our optical monitoring of the BL Lac object S5 1803+78 from 1996 to 2011. The source showed no clear periodicity, but a time scale of about 1300 days between major flares is possibly present. No systematic trend of the color index with flux variations is evident, at variance with other BL Lacs. In one flare, however, the source was bluer in the rising phase and redder in the falling one. Two Gamma-ray flares were detected by Fermi-GST during our monitoring: on the occasion of only one of them we found simultaneous optical brightening. A one-zone Synchrotron Self Compton (SSC) model appears too simple to explain the source behavior.

We present the results of our analysis devoted to the research of sources emitting in the energy bands surveyed by both the Swift-BAT and the Fermi-LAT telescopes. We cross-correlate the Fermi-LAT 1-year point source catalogue (1FGL) of {\gamma}-ray sources and the second Palermo BAT catalogue (2PBC) of hard X-ray sources, establishing a correspondence between sources when their error boxes overlap. We also extract the significance value in the BAT 15-150 keV map, obtained using a dedicated software for the reduction of BAT data, in the direction of the 1FGL sources and take into account those above the significance threshold {\sigma} = 3. We obtain a sample of common sources emitting in both the hard X- and the {\gamma}-ray energy bands and evaluate its content in galactic and extragalactic objects. We assess the fraction of unidentified sources and describe in greater detail the properties of two of them, 1FGL J0137.8+5814 and 1FGL J2056.7+4938, supporting their classification as blazars after the analysis of their broad-band spectral energy distribution. We discuss the blazar content of the collected 1FGL-2PBC sources: we build its redshift distibution and compare it with that of the whole blazar population as reported in the second edition of the BZCAT blazar catalogue.

We present the new parallel version (pCRASH2) of the cosmological radiative transfer code CRASH2 for distributed memory supercomputing facilities. The code is based on a static domain decomposition strategy inspired by geometric dilution of photons in the optical thin case that ensures a favourable performance speed-up with increasing number of computational cores. Linear speed-up is ensured as long as the number of radiation sources is equal to the number of computational cores or larger. The propagation of rays is segmented and rays are only propagated through one sub-domain per time step to guarantee an optimal balance between communication and computation. We have extensively checked pCRASH2 with a standardised set of test cases to validate the parallelisation scheme. The parallel version of CRASH2 can easily handle the propagation of radiation from a large number of sources and is ready for the extension of the ionisation network to species other than hydrogen and helium.

We present the Second Palermo Swift-BAT hard X-ray catalogue obtained by analysing data acquired in the first 54 months of the Swift mission. Using our software dedicated to the analysis of data from coded mask telescopes, we analysed the BAT survey data in three energy bands (15-30 keV, 15-70 keV, 15-150 keV), obtaining a list of 1256 detections above a significance threshold of 4.8 standard deviations. The identification of the source counterparts is pursued using two strategies: the analysis of field observations of soft X-ray instruments and cross-correlation of our catalogue with source databases.The survey covers 50% of the sky to a 15--150 keV flux limit of 1.0 x 10^-11 erg s^-1 cm^-2 and 9.2 x 10^-12 erg s^-1 cm^-2 for |b|< 10 degrees and |b|> 10 degrees, respectively. The Second Palermo Swift-BAT hard X-ray catalogue includes 1079 (86%) hard X-ray sources with an associated counterpart (26 with a double association and 2 with a triple association) and 177 BAT excesses (14%) that still lack a counterpart. The distribution of the BAT sources among the different object classes consists of 19% Galactic sources, 57% extragalactic sources, and 10% sources with a counterpart at softer energies whose nature has not yet been determined. About half of the BAT associated sources lack a counterpart in the ROSAT catalogues. This suggests that either moderate or strong absorption may be preventing their detection in the ROSAT energy band. The comparison of our BAT catalogue with the Fermi Large Area Telescope First Source Catalogue identifies 59 BAT/Fermi correspondences: 48 blazars, 3 Seyfert galaxies, 1 interacting galaxy, 3 high mass X-ray binaries, and 4 pulsars/supernova remnants. This small number of correspondences indicates that different populations make the sky shine in these two different energy bands.

The 2nd edition of the Roma-BZCAT is available on line at the ASDC website (http://www.asdc.asi.it/bzcat) and in the NED database. In this short paper we describe the major updates from the first edition.

We present the first Lyman Alpha Emitter (LAE) study that combines: (i) cosmological SPH simulations run using GADGET-2, (ii) radiative transfer simulations (CRASH), and (iii) a previously developed LAE model. This complete LAE model accounts for the intrinsic LAE Lyman Alpha/continuum luminosity, dust enrichment and Lyman Alpha transmission through the intergalactic medium (IGM), to quantify the effects of reionization, dust and velocity fields on the Lyman Alpha and UV Luminosity Functions (LF). We find that a model neglecting dust sorely fails to reproduce either the slope or the magnitude of the observed Lyman Alpha and UV LFs. Clumped dust is required to simultaneously fit the observed UV and Lyman Alpha LFs, such that the intrinsic Lyman Alpha-to-continuum luminosity is enhanced by a factor f_alpha/f_c ~ 1.5 (3.7) excluding (including) peculiar velocities. The higher value including velocity fields arises since LAEs reside in large potential wells and inflows decrease their Lyman Alpha transmission. For the first time, a degeneracy is found between the the ionization state of the IGM and the clumping of dust inside high-redshift galaxies. The Lyman Alpha LF at z ~ 5.7 can be well reproduced (to within a 5 sigma error) by a wide range of IGM average neutral hydrogen fraction, 3.4/10^3 to 0.16, provided that the increase in the Lyman Alpha transmission through a more ionized IGM is compensated by a decrease in the Lyman Alpha escape fraction from the galaxy due to dust absorption. The physical properties of LAEs are presented, along with a discussion of the assumptions adopted.

We present the Palermo Swift-BAT hard X-ray catalogue obtained from the analysis of the the data relative to the first 39 months of the Swift mission. We have developed a dedicated software to perform data reduction, mosaicking and source detection on the BAT survey data. We analyzed the BAT dataset in three energy bands (14-150 keV, 14-30 keV, 14-70 keV), obtaining a list of 962 detections above a significance threshold of 4.8 standard deviations. The identification of the source counterparts was pursued using three strategies: cross-correlation with published hard X-ray catalogues, analysis of field observations of soft X-ray instruments, cross-correlation with the SIMBAD databases. The survey covers 90% of the sky down to a flux limit of 2.5x10E-11 erg/cm2/s and 50% of the sky down to a flux limit of 1.8x10E-11 erg/cm2/s in the 14-150 keV band. We derived a catalogue of 754 identified sources, of which ~69% are extragalactic, ~27% are Galactic objects, ~4% are already known X-ray or gamma ray emitters whose nature has not been determined yet. The integrated flux of the extragalactic sample is ~1% of the Cosmic X-ray background in the 14-150 keV range.

We propose the Apparent Shrinking Criterion (ASC) to interpret the spatial extent, R_w, of transmitted flux windows in the absorption spectra of high-z quasars. The ASC can discriminate between the two regimes in which R_w corresponds either to the physical size, R_HII, of the quasar HII region, or to the distance, R^{max}_w, at which the transmitted flux drops to =0.1 and a Gunn-Peterson (GP) trough appears. In the first case (HR regime), one can determine the IGM mean HI fraction, x_HI; in the second (PR regime), the value of R_w allows to measure the local photoionization rate and the local enhancement of the photoionization rate, Gamma_G, due to nearby/intervening galaxies. The ASC has been tested against radiative transfer+SPH numerical simulations, and applied to 15 high-z (z>5.8) quasars sample from Fan et al. (2006). All sample quasars are found to be in the PR regime; hence, their observed spectral properties (inner flux profile, extent of transmission window) cannot reliably constrain the value of x_HI. Four sample quasars show evidence for a local enhancement (up to 50%) in the local photoionization rate possibly produced by a galaxy overdensity. We discuss the possible interpretations and uncertainties of this result.

In this paper we report on the improvements implemented in the cosmological radiative transfer code CRASH. In particular we present a new multi-frequency algorithm for spectra sampling which makes use of colored photon packets: we discuss the need for the multi-frequency approach, describe its implementation and present the improved CRASH performance in reproducing the effects of ionizing radiation with an arbitrary spectrum. We further discuss minor changes in the code implementation which allow for more efficient performance and an increased precision.

We present a new catalogue of blazars based on multi-frequency surveys and on an extensive review of the literature. Blazars are classified as BL Lacertae objects, as flat spectrum radio quasars or as blazars of uncertain/transitional type. Each object is identified by a root name, coded as BZB, BZQ and BZU for these three subclasses respectively, and by its coordinates. This catalogue is being built as a tool useful for the identification of the extragalactic sources that will be detected by present and future experiments for X and gamma-ray astronomy, like Swift, AGILE, Fermi-GLAST and Simbol-X. An electronic version is available from the ASI Science Data Center web site at http://www.asdc.asi.it/bzcat

In this paper we present CRASH_alpha, the first radiative transfer code for cosmological application that follows the parallel propagation of Ly_alpha and ionizing photons. CRASH_alpha is a version of the continuum radiative transfer code CRASH with a new algorithm to follow the propagation of Ly_alpha photons through a gas configuration whose ionization structure is evolving. The implementation introduces the time evolution for Ly_alpha photons (a feature commonly neglected in line radiative transfer codes) and, to reduce the computational time needed to follow each scattering, adopts a statistical approach to the Ly_alpha treatment by making extensive use of pre-compiled tables. With this statistical approach we experience a drastic increase of the computational speed and, at the same time, an excellent agreement with the full Ly_alpha radiative transfer computations of the code MCLy_alpha. We find that the emerging spectra keep memory of the ionization history which generates a given ionization configuration of the gas and, to properly account for this effect, a self-consistent joint evolution of line and ionizing continuum radiation as implemented in CRASH_alpha is necessary. A comparison between the results from our code and from Ly_alpha scattering alone on a fixed HI density field shows that the extent of the difference between the emerging spectra depends on the particular configuration considered, but it can be substantial and can thus affect the physical interpretation of the problem at hand. These differences should furthermore be taken into account when computing the impact of the Ly_alpha radiation on e.g. the observability of the 21 cm line from neutral hydrogen at epochs preceeding complete reionization.

We studied the temporal and spectral evolution of the synchrotron emission from the high energy peaked BL Lac object 1E 1207.9+3945. Two recent observations have been performed by the XMM-Newton and Swift satellites; we carried out X-ray spectral analysis for both of them, and photometry in optical-ultraviolet filters for the Swift one. Combining the results thus obtained with archival data we built the long-term X-ray light curve, spanning a time interval of 26 years, and the Spectral Energy Distribution (SED) of this source. The light curve shows a large flux increasing, about a factor of six, in a time interval of a few years. After reaching its maximum in coincidence with the XMM-Newton pointing in December 2000 the flux decreased in later years, as revealed by Swift. The very good statistics available in the 0.5-10 keV XMM-Newton X-ray spectrum points out a highly significant deviation from a single power law. A log-parabolic model with a best fit curvature parameter of 0.25 and a peak energy at ~1 keV describes well the spectral shape of the synchrotron emission. The simultaneous fit of Swift UVOT and XRT data provides a milder curvature (b~0.1) and a peak at higher energies (~15 keV), suggesting a different state of source activity. In both cases UVOT data support the scenario of a single synchrotron emission component extending from the optical/UV to the X-ray band. New X-ray observations are important to monitor the temporal and spectral evolution of the source; new generation gamma-ray telescopes like AGILE and GLAST could for the first time detect its inverse Compton emission.

We investigate the possibility of constraining the ionization state of the intergalactic medium (IGM) close to the end of reionization (z ~ 6) by measuring the size of the HII regions in high-z quasars spectra. We perform a combination of multiphase smoothed particle hydrodynamics (SPH) and 3D radiative transfer (RT) simulations to reliably predict the properties of typical high-z quasar HII regions, embedded in a partly neutral IGM (x_HI=0.1). In this work we assume a fixed configuration for the quasar lifetime and luminosity, i.e. t_Q=10^7 yr and N_gamma=5.2 10^56 s-1. From the analysis of mock spectra along lines of sight through the simulated QSO environment, we find that the HII region size derived from quasar spectra is on average 30 per cent smaller than the physical one. Additional maximum likelihood analysis shows that this offset induces an overestimate of the neutral hydrogen fraction, x_HI, by a factor of 3. By applying the same statistical method to a sample of observed QSOs our study favors a mostly ionized (x_HI < 0.06) universe at z=6.1.

We present the results of a spectral analysis of 5 Swift XRT and UVOT observations of the BL Lac object PKS 0548-322 carried out over the period April-June 2005. The X-ray flux of this high energy peaked BL Lac (HBL) source was found to be approximately constant at a level of F(2-10 keV) ~ 4x10^-11 erg cm^-2 s^-1, a factor of 2 brighter than when observed by BeppoSAX in 1999 and close to the maximum intensity reported in the Einstein Slew Survey. The very good statistics obtained in the 0.3-10 keV Swift X-ray spectrum allowed us to detect highly significant deviations from a simple power law spectral distribution. A log-parabolic model describes well the X-ray data and gives a best fit curvature parameter of 0.18 and peak energy in the Spectral Energy Distribution of about 2 keV. The UV spectral data from Swift UVOT join well with a power law extrapolation of the soft X-ray data points suggesting that the same component is responsible for the observed emission in the two bands. The combination of synchrotron peak in the X-ray band and high intensity state confirms PKS 0548-322 as a prime target for TeV observations. X-ray monitoring and coordinated TeV campaigns are highly advisable.

Using a Monte Carlo technique, we have developed a 3D lyman-alpha radiation transfer code allowing for prescribed arbitrary hydrogen density, ionisation, temperature structures, and dust distribution, and arbitrary velocity fields and UV photon sources. We have examined the lyman-alpha line profiles predicted for several simple geometrical configurations and their dependence on the main input parameters. Overall, we find line profiles reaching from doubly peaked symmetric emission to symmetric Voigt (absorption) in static configurations with increasing dust content, and asymmetric red-(blue-) shifted emission lines with a blue (red) counterpart ranging from absorption to emission (with increasing line/continuum strength) in expanding (infalling) media. The following results are of interest for the interpretation of lya profiles from galaxies. 1/ Standard lya absorption line fitting of global spectra of galaxies may lead to an underestimate of the true hydrogen column density in certain geometrical conditions. 2/ Normal (inverted) P-Cygni like lya profiles can be obtained in expanding (infalling) media from objects without any intrinsic lya emission, as a natural consequence of radiation transfer effects. 3/ The formation and the detailed shape of lya profiles resulting from expanding shells has been thoroughly revised: for sufficiently large column densities, the position of the main lya emission peak is redshifted by twice the expansion velocity.This is in excellent agreement with the observations of z~3 LBGs. This finding indicates also that large scale, fairly symmetric shell structures must be a good description for the outflows in LBGs.(shortened abstract)

Radiative transfer simulations are now at the forefront of numerical astrophysics. They are becoming crucial for an increasing number of astrophysical and cosmological problems; at the same time their computational cost has come to the reach of currently available computational power. Further progress is retarded by the considerable number of different algorithms (including various flavours of ray-tracing and moment schemes) developed, which makes the selection of the most suitable technique for a given problem a non-trivial task. Assessing the validity ranges, accuracy and performances of these schemes is the main aim of this paper, for which we have compared 11 independent RT codes on 5 test problems: (0) basic physics, (1) isothermal H II region expansion and (2) H II region expansion with evolving temperature, (3) I-front trapping and shadowing by a dense clump, (4) multiple sources in a cosmological density field. The outputs of these tests have been compared and differences analyzed. The agreement between the various codes is satisfactory although not perfect. The main source of discrepancy appears to reside in the multi-frequency treatment approach, resulting in different thicknesses of the ionized-neutral transition regions and different temperature structure. The present results and tests represent the most complete benchmark available for the development of new codes and improvement of existing ones. To this aim all test inputs and outputs are made publicly available in digital form.

We address the topic of the Intra-Night Optical Variability of the BL Lac object S5 0716+714. To this purpose a long term observational campaign was performed, from 1996 to 2003, which allowed the collection of a very large data set, containing 10,675 photometric measurements obtained in 102 nights. The source brightness varied in a range of about 2 mag, although the majority of observations were performed when it was in the range 13.0 < R < 13.75. Variability time scales were estimated from the rates of magnitude variation, which were found to have a distribution function well fitted by an exponential law with a mean value of 0.027 mag/h, corresponding to an e-folding time scale of the flux tau_F = 37.6 h. The highest rates of magnitude variation were around 0.10--0.12 mag/h and lasted less than 2 h. These rates were observed only when the source had an R magnitude < 13.4, but this finding cannot be considered significant because of the low statistical occurrence. The distribution of tau_F has a well defined modal value at 19 h. Assuming the recent estimate of the beaming factor delta about 20, we derived a typical size of the emitting region of about 5 times 10^{16}/(1 + z) cm. The possibility to search for a possible correlation between the mean magnitude variation rate and the long term changes of the velocity of superluminal components in the jet is discussed.

Strong observational evidence for a fluctuating ultraviolet background (UVB) has been accumulating through a number of studies of the HI and HeII Lya forest as well as accurate IGM metallicity measurements. UVB fluctuations could arise both from the inhomogeneous distribution of the ionizing sources and/or from radiative transfer (RT) through the filamentary IGM. In this study we investigate, via numerical simulations, the role of RT effects such as shadowing, self-shielding and filtering of the ionizing radiation, in giving raise to a fluctuating UVB. We focus on possible detectable signatures of these effects on quantities derived from Lya forest spectra, as photoionization rate fluctuations, eta parameter (the HeII to HI column density ratio) distributions and the IGM temperature at redshift about 3. We find that RT induces fluctuations up to 60% in the UVB, which are tightly correlated to the density field. The UVB mean intensity is progressively suppressed toward higher densities and photon energies above 4 Ryd, due to the high HeII opacity. Shielding of overdense regions (Delta > 5) from cosmic HeII ionizing radiation, produces a decreaseing trend of eta with overdensity. Furthermore we find that the mean eta value inferred from HI-HeII Lya forest observations can be explained only by properly accounting for the actual IGM opacity. We outline and discuss several implications of our findings.

Recent observations have shown that the intergalactic medium (IGM) is more transparent to Lyalpha photons close to Lyman Break Galaxies (LBGs) than at large distance from them, ie a proximity effect. Cosmological simulations including winds from LBGs have been so far unable to explain this trend. By coupling such simulations with the radiative transfer code CRASH, we investigate whether the addition of the ionizing radiation emitted by LBGs can increase the transmissivity by decreasing the neutral hydrogen fraction in the inner Mpc of the galaxy halo. The transmissivity as a function of distance is roughly reproduced only if LBGs are identified with dwarf galaxies (with masses < 10^9 solar masses), which are undergoing a vigorous (50 solar masses/yr) burst of star formation. Similar star formation rates in larger galaxies are not sufficient to overwhelm the large recombination rates associated with their denser environment. If so, photoionization partly reconciles theory with observations, although we discuss a number of uncertainties affecting both approaches.

We present a largely improved version of CRASH, a 3-D radiative transfer code that treats the effects of ionizing radiation propagating through a given inhomogeneous H/He cosmological density field, on the physical conditions of the gas. The code, based on a Monte Carlo technique, self-consistently calculates the time evolution of gas temperature and ionization fractions due to an arbitrary number of point/extended sources and/or diffuse background radiation with given spectra. In addition, the effects of diffuse ionizing radiation following recombinations of ionized atoms have been included. After a complete description of the numerical scheme, to demonstrate the performances, accuracy, convergency and robustness of the code we present four different test cases designed to investigate specific aspects of radiative transfer: (i) pure hydrogen isothermal Stromgren sphere; (ii) realistic Stromgren spheres; (iii) multiple overlapping point sources, and (iv) shadowing of background radiation by an intervening optically thick layer. When possible, detailed quantitative comparison of the results against either analytical solutions or 1-D standard photoionization codes has been made showing a good level of agreement. For more complicated tests the code yields physically plausible results, which could be eventually checked only by comparison with other similar codes. Finally, we briefly discuss future possible developments and cosmological applications of the code.

Motivated by the relative lack of neutral hydrogen around Lyman Break Galaxies deduced from recent observations, we investigate the properties of the Lyalpha forest around high redshift galaxies. The study is based on improved numerical SPH simulations implementing, in addition to standard processes, a new scheme for multiphase and outflow physics description. Although on large scales our simulations reproduce a number of statistical properties of the IGM (because of the small filling factor of shock-heated gas), they underpredict the Lyalpha optical depth decrease inside 1 Mpc/h of the galaxies by a factor of ~2. We interpret this result as due to the combined effect of infall occurring along the filaments, which prevents efficient halo gas clearing by the outflow, and the insufficient increase of (collisional) hydrogen ionization produced by the temperature increase inside the hot, outflow-carved bubble. Unless an observational selection bias is present, we speculate that local photoionization could be the only viable explanation to solve the puzzle.