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.

The black hole uniqueness and the no-hair theorems imply that the quasinormal spectrum of any astrophysical black hole is determined solely by its mass and spin. The countably infinite number of quasinormal modes of a Kerr black hole are thus related to each other and any deviations from these relations provide a strong hint for physics beyond the general theory of relativity. To test the no-hair theorem using ringdown signals, it is necessary to detect at least two quasinormal modes. In particular, one can detect the fundamental mode along with a subdominant overtone or with another angular mode, depending on the mass ratio and the spins of the progenitor binary. Also in the light of the recent discovery of GW190412, studying how the mass ratio affects the prospect of black hole spectroscopy using overtones or angular modes is pertinent, and this is the major focus of our study. First, we provide ready-to-use fits for the amplitudes and phases of both the angular modes and overtones as a function of mass ratio $q\in[0,10]$. Using these fits we estimate the minimum signal-to-noise ratio for detectability, resolvability, and measurability of subdominant modes/tones. We find that performing black-hole spectroscopy with angular modes is preferable when the binary mass ratio is larger than $q\approx 1.2$ (provided that the source is not located at a particularly disfavoured inclination angle). For nonspinning, equal-mass binary black holes, the overtones seem to be the only viable option to perform a spectroscopy test of the no-hair theorem. However this would require a large ringdown signal-to-noise ratio ($\approx 100$ for a $5\%$ accuracy test with two overtones) and the inclusion of more than one overtone to reduce modelling errors, making black-hole spectroscopy with overtones impractical in the near future.

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.

In viable models of minicharged dark matter, astrophysical black holes might be charged under a hidden $U(1)$ symmetry and are formally described by the same Kerr-Newman solution of Einstein-Maxwell theory. These objects are unique probes of minicharged dark matter and dark photons. We show that the recent gravitational-wave detection of a binary black-hole coalescence by aLIGO provides various observational bounds on the black hole's charge, regardless of its nature. The pre-merger inspiral phase can be used to constrain the dipolar emission of (ordinary and dark) photons, whereas the detection of the quasinormal modes set an upper limit on the final black hole's charge. By using a toy model of a point charge plunging into a Reissner-Nordstrom black hole, we also show that in dynamical processes the (hidden) electromagnetic quasinormal modes of the final object are excited to considerable amplitude in the gravitational-wave spectrum only when the black hole is nearly extremal. The coalescence produces a burst of low-frequency dark photons which might provide a possible electromagnetic counterpart to black-hole mergers in these scenarios.

Validating the black-hole no-hair theorem with gravitational-wave observations of compact binary coalescences provides a compelling argument that the remnant object is indeed a black hole as described by the general theory of relativity. This requires performing a spectroscopic analysis of the post-merger signal and resolving the frequencies of either different angular modes or overtones (of the same angular mode). For a nearly-equal mass binary black-hole system, only the dominant angular mode ($l=m=2$) is sufficiently excited and the overtones are instrumental to perform this test. Here we investigate the robustness of modelling the post-merger signal of a binary black hole coalescence as a superposition of overtones. Further, we study the bias expected in the recovered frequencies as a function of the start time of a spectroscopic analysis and provide a computationally cheap procedure to choose it based on the interplay between the expected statistical error due to the detector noise and the systematic errors due to waveform modelling. Moreover, since the overtone frequencies are closely spaced, we find that resolving the overtones is particularly challenging and requires a loud ringdown signal. Rayleigh's resolvability criterion suggests that in an optimistic scenario a ringdown signal-to-noise ratio larger than $\sim 30$ (achievable possibly with LIGO at design sensitivity and routinely with future interferometers such as Einstein Telescope, Cosmic Explorer, and LISA) is necessary to resolve the overtone frequencies. We then conclude by discussing some conceptual issues associated with black-hole spectroscopy with overtones.

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.

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

In this paper we describe the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to accretion flows in the strong field gravity regime around both stellar-mass and supermassive black-holes. eXTP has the unique capability of using advanced 'spectral-timing-polarimetry' techniques to analyze the rapid variations with three orthogonal diagnostics of the flow and its geometry, yielding unprecedented insight into the inner accreting regions, the effects of strong field gravity on the material within them and the powerful outflows which are driven by the accretion process.

A tentative detection of gravitational-wave echoes in the post-merger signal of GW170817 has been recently claimed at $4.2\sigma$ significance level. It has been speculated that the signal might provide evidence for near-horizon quantum structures in the remnant exotic object. We point out that if the remnant object is an ultracompact neutron star, echoes are expected for objects with radius only slightly smaller than that of an ordinary neutron star. The reported echoes at ~72 Hz are compatible with a toy model of incompressible star with mass approximately $M\in(2,3) M_\odot$ and radius close to the Buchdahl limit, R~9GM/(4c^2). If confirmed, low-frequency gravitational-wave echoes would be in tension with all current neutron-star models and would have dramatic implications for nuclear physics and gravity.

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.

Gravitational-wave astronomy can give us access to the structure of black holes, potentially probing microscopic or even Planckian corrections at the horizon scale, as those predicted by some quantum-gravity models of exotic compact objects. A generic feature of these models is the replacement of the horizon by a reflective surface. Objects with these properties are prone to the so-called ergoregion instability when they spin sufficiently fast. We investigate in detail a simple model consisting of scalar perturbations of a Kerr geometry with a reflective surface near the horizon. The instability depends on the spin, on the compactness, and on the reflectivity at the surface. The instability time scale increases only logarithmically in the black-hole limit and, for a perfectly reflecting object, this is not enough to prevent the instability from occurring on dynamical time scales. However, we find that an absorption rate at the surface as small as 0.4% (reflectivity coefficient as large as $|{\cal R}|^2=0.996$) is sufficient to quench the instability completely. Our results suggest that exotic compact objects are not necessarily ruled out by the ergoregion instability.

In a core-collapse supernova, a huge amount of energy is released in the Kelvin-Helmholtz phase subsequent to the explosion, when the proto-neutron star cools and deleptonizes as it loses neutrinos. Most of this energy is emitted through neutrinos, but a fraction of it can be released through gravitational waves. We model the evolution of a proto-neutron star in the Kelvin-Helmholtz phase using a general relativistic numerical code, and a recently proposed finite temperature, many-body equation of state; from this we consistently compute the diffusion coefficients driving the evolution. To include the many-body equation of state, we develop a new fitting formula for the high density baryon free energy at finite temperature and intermediate proton fraction. We estimate the emitted neutrino signal, assessing its detectability by present terrestrial detectors, and we determine the frequencies and damping times of the quasi-normal modes which would characterize the gravitational wave signal emitted in this stage.

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.

Gravitational waves emitted by distorted black holes---such as those arising from the coalescence of two neutron stars or black holes---carry not only information about the corresponding spacetime but also about the underlying theory of gravity. Although general relativity remains the simplest, most elegant and viable theory of gravitation, there are generic and robust arguments indicating that it is not the ultimate description of the gravitational universe. Here we focus on a particularly appealing extension of general relativity, which corrects Einstein's theory through the addition of terms which are second order in curvature: the topological Gauss-Bonnet invariant coupled to a dilaton. We study gravitational-wave emission from black holes in this theory, and {\bf(i)} find strong evidence that black holes are linearly (mode) stable against both axial and polar perturbations; {\bf(ii)} discuss how the quasinormal modes of black holes can be excited during collisions involving black holes, and finally {\bf(iii)} show that future ringdown detections with large signal-to-noise ratio would improve current constraints on the coupling parameter of the theory.

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.

We study the evolution of the rotation rate of a proto-neutron star, born in a core-collapse supernova, in the first seconds of its life. During this phase, the star evolution can be described as a sequence of stationary configurations, which we determine by solving the neutrino transport and the stellar structure equations in general relativity. We include in our model the angular momentum loss due to neutrino emission. We find that the requirement of a rotation rate not exceeding the mass-shedding limit at the beginning of the evolution implies a strict bound on the rotation rate at later times. Moreover, assuming that the proto-neutron star is born with a finite ellipticity, we determine the emitted gravitational wave signal, and estimate its detectability by present and future ground-based interferometric detectors.

By extending our recent framework to describe the tidal deformations of a spinning compact object, we compute for the first time the tidal Love numbers of a spinning neutron star to linear order in the angular momentum. The spin of the object introduces couplings between electric and magnetic distortions and new classes of spin-induced ("rotational") tidal Love numbers emerge. We focus on stationary tidal fields, which induce axisymmetric perturbations. We present the perturbation equations for both electric-led and magnetic-led rotational Love numbers for generic multipoles and explicitly solve them for various tabulated equations of state and for a tidal field with an electric (even parity) and magnetic (odd parity) component with $\ell=2,3,4$. For a binary system close to the merger, various components of the tidal field become relevant. In this case we find that an octupolar magnetic tidal field can significantly modify the mass quadrupole moment of a neutron star. Preliminary estimates, assuming a spin parameter $\chi\approx0.05$, show modifications $\gtrsim10\%$ relative to the static case, at an orbital distance of five stellar radii. Furthermore, the rotational Love numbers as functions of the moment of inertia are much more sensitive to the equation of state than in the static case, where approximate universal relations at the percent level exist. For a neutron-star binary approaching the merger, we estimate that the approximate universality of the induced mass quadrupole moment deteriorates from $1\%$ in the static case to roughly $6\%$ when $\chi\approx0.05$. Our results suggest that spin-tidal couplings can introduce important corrections to the gravitational waveforms of spinning neutron-star binaries approaching the merger.

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.

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.

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

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.

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.

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.

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.

In this paper we perform a detailed analysis of the effect of various approximations which have been used in the literature to compute the detection rates of compact binary coalescences for interferometric gravitational wave detectors. We evaluate the detection rates for the coalescence of BH-BH, NS-NS, and BH-NS binaries taking into account their specific statistical properties obtained from population synthesis models, the cosmic star formation rate history, and the effects of redshift on the emitted gravitational wave signals. The results are compared with those obtained with procedures that are based on different levels of approximations, such as adopting averaged values for the total mass and symmetric mass ratio for all the systems of a binary population, using these to compute the horizon distance for individual detectors, or estimating the coalescence rate density within this distance by its local value. We show that these approximations introduce a bias which depends not only on the approximation which is used, but also on the statistical properties of the considered binary population.

The advanced interferometer network will herald a new era in observational astronomy. There is a very strong science case to go beyond the advanced detector network and build detectors that operate in a frequency range from 1 Hz-10 kHz, with sensitivity a factor ten better in amplitude. Such detectors will be able to probe a range of topics in nuclear physics, astronomy, cosmology and fundamental physics, providing insights into many unsolved problems in these areas.

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.

Physics in curved spacetime describes a multitude of phenomena, ranging from astrophysics to high energy physics. The last few years have witnessed further progress on several fronts, including the accurate numerical evolution of the gravitational field equations, which now allows highly nonlinear phenomena to be tamed. Numerical relativity simulations, originally developed to understand strong field astrophysical processes, could prove extremely useful to understand high-energy physics processes like trans-Planckian scattering and gauge-gravity dualities. We present a concise and comprehensive overview of the state-of-the-art and important open problems in the field(s), along with guidelines for the next years. This writeup is a summary of the "NR/HEP Workshop" held in Madeira, Portugal from August 31st to September 3rd 2011.

We compute the gravitational wave background (GWB) generated by a cosmological population of (BH-BH) binaries using hybrid waveforms recently produced by numerical simulations of (BH-BH) coalescence, which include the inspiral, merger and ring-down contributions. A large sample of binary systems is simulated using the population synthesis code SeBa, and we extract fundamental statistical information on (BH-BH) physical parameters (primary and secondary BH masses, orbital separations and eccentricities, formation and merger timescales). We then derive the binary birth and merger rates using the theoretical cosmic star formation history obtained from a numerical study which reproduces the available observational data at redshifts $z < 8$. We evaluate the contributions of the inspiral, merger and ring-down signals to the GWB, and discuss how these depend on the parameters which critically affect the number of coalescing (BH-BH) systems. We find that Advanced LIGO/Virgo have a chance to detect the GWB signal from the inspiral phase with a $(S/N)=10$ only for the most optimistic model, which predicts the highest local merger rate of 0.85 Mpc$^{-3}$ Myr$^{-1}$. Third generation detectors, such as ET, could reveal the GWB from the inspiral phase predicted by any of the considered models. In addition, ET could sample the merger phase of the evolution at least for models which predict local merger rates between $[0.053 - 0.85]$ Mpc$^{-3}$ Myr$^{-1}$, which are more than a factor 2 lower the the upper limit inferred from the analysis of the LIGO S5 run\cite{Abadieetal2011}.

This paper has been withdrawn by the authors because it had major revisions.

We study how the frequencies and damping times of oscillations of a newly born, hot proto-neutron star depend on the physical quantities which characterize the star quasi-stationary evolution which follows the bounce. Stellar configurations are modeled using a microscopic equation of state obtained within the Brueckner-Hartree-Fock, nuclear many-body approach, extended to the finite-temperature regime. We discuss the mode frequency behaviour as function of the lepton composition, and of the entropy gradients which prevail in the interior of the star. We find that, in the very early stages, gravitational wave emission efficiently competes with neutrino processes in dissipating the star mechanical energy residual of the gravitational collapse.

Neutron stars can have, in some phases of their life, extremely strong magnetic fields, up to 10^15-10^16 G. These objects, named magnetars, could be powerful sources of gravitational waves, since their magnetic field could determine large deformations. We discuss the structure of the magnetic field of magnetars, and the deformation induced by this field. Finally, we discuss the perspective of detection of the gravitational waves emitted by these stars.

We construct general relativistic models of stationary, strongly magnetized neutron stars. The magnetic field configuration, obtained by solving the relativistic Grad-Shafranov equation, is a generalization of the twisted torus model recently proposed in the literature; the stellar deformations induced by the magnetic field are computed by solving the perturbed Einstein's equations; stellar matter is modeled using realistic equations of state. We find that in these configurations the poloidal field dominates over the toroidal field and that, if the magnetic field is sufficiently strong during the first phases of the stellar life, it can produce large deformations.

We find general relativistic solutions of equilibrium magnetic field configurations in magnetars, extending previous results of Colaiuda et al. (2008). Our method is based on the solution of the relativistic Grad-Shafranov equation, to which Maxwell's equations can be reduced in some limit. We obtain equilibrium solutions with the toroidal magnetic field component confined into a finite region inside the star, and the poloidal component extending to the exterior. These so-called twisted-torus configurations have been found to be the final outcome of dynamical simulations in the framework of Newtonian gravity, and appear to be more stable than other configurations. The solutions include higher order multipoles, which are coupled to the dominant dipolar field. We use arguments of minimal energy to constrain the ratio of the toroidal to the poloidal field.

Two classes of high energy sources in our galaxy are believed to host magnetars, neutron stars whose emission results from the dissipation of their magnetic field. The extremely high magnetic field of magnetars distorts their shape, and causes the emission of a conspicuous gravitational waves signal if rotation is fast and takes place around a different axis than the symmetry axis of the magnetic distortion. Based on a numerical model of the cosmic star formation history, we derive the cosmological background of gravitational waves produced by magnetars, when they are very young and fast spinning. We adopt different models for the configuration and strength of the internal magnetic field (which determines the distortion) as well as different values of the external dipole field strength (which governs the spin evolution of magnetars over a wide range of parameters). We find that the expected gravitational wave background differs considerably from one model to another. The strongest signals are generated for magnetars with very intense toroidal internal fields ($\sim 10^{16}$ G range) and external dipole fields of $\sim 10^{14}$, as envisaged in models aimed at explaining the properties of the Dec 2004 giant flare from SGR 1806-20. Such signals should be easily detectable with third generation ground based interferometers such as the Einstein Telescope.

Astrophysical sources emit gravitational waves in a large variety of processes occurred since the beginning of star and galaxy formation. These waves permeate our high redshift Universe, and form a background which is the result of the superposition of different components, each associated to a specific astrophysical process. Each component has different spectral properties and features that it is important to investigate in view of a possible, future detection. In this contribution, we will review recent theoretical predictions for backgrounds produced by extragalactic sources and discuss their detectability with current and future gravitational wave observatories.

We study the tidal effects of a Kerr black hole on a neutron star in black hole-neutron star binary systems using a semi-analytical approach which describes the neutron star as a deformable ellipsoid. Relativistic effects on the neutron star self-gravity are taken into account by employing a scalar potential resulting from relativistic stellar structure equations. We calculate quasi-equilibrium sequences of black hole-neutron star binaries, and the critical orbital separation at which the star is disrupted by the black hole tidal field: the latter quantity is of particular interest because when it is greater than the radius of the innermost stable circular orbit, a short gamma-ray burst scenario may develop.

We study the tidal disruption of neutron stars in black hole-neutron star coalescing binaries. We calculate the critical orbital separation at which the star is disrupted by the black hole tidal field for several equations of state describing the matter inside the neutron star, and for a large set of the binary parameters. When the disruption occurs before the star reaches the innermost stable circular orbit, the gravitational wave (GW) signal emitted by the system is expected to exhibit a cutoff frequency nu_GWtide, which is a distinctive feature of the waveform. We evaluate nu_GWtide and show that, if this frequency will be found in a detected gravitational wave, it will allow to determine the neutron star radius with an error of a few percent, providing valuable information on the behaviour of matter in the stellar core.

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

Using the results of a numerical simulation which follows the evolution, metal enrichment and energy deposition of both Population III and Population II stars, we predict the redshift dependence of the formation rate of black hole remnants of Population III stars with masses 100- 500Msun and of neutron stars(black holes) remnants of Population II stars with masses 8-20Msun (20-40Msun). We describe the gravitational wave spectrum produced by Population III and Population II sources adopting the most appropriate signals available in the literature and we compute the stochastic backgrounds resulting from the cumulative emission of these sources throughout the history of the Universe. With the aim of assessing whether these backgrounds might act as foregrounds for signals generated in the Inflationary epoch, we compare their amplitudes with the sensitivity of currently planned and future ground/space-based interferometers.

We review the main results obtained in the literature on quasi-normal modes of compact stars and black holes, in the light of recent exciting developments of gravitational wave detectors. Quasi-normal modes are a fundamental feature of the gravitational signal emitted by compact objects in many astrophysical processes; we will show that their eigenfrequencies encode interesting information on the nature and on the inner structure of the emitting source and we will discuss whether we are ready for a gravitational wave asteroseismology.

We find numerical solutions of the coupled system of Einstein-Maxwell's equations with a linear approach, in which the magnetic field acts as a perturbation of a spherical neutron star. In our study, magnetic fields having both poloidal and toroidal components are considered, and higher order multipoles are also included. We evaluate the deformations induced by different field configurations, paying special attention to those for which the star has a prolate shape. We also explore the dependence of the stellar deformation on the particular choice of the equation of state and on the mass of the star. Our results show that, for neutron stars with mass M = 1.4 Msun and surface magnetic fields of the order of 10^15 G, a quadrupole ellipticity of the order of 10^(-6) - 10^(-5) should be expected. Low mass neutron stars are in principle subject to larger deformations (quadrupole ellipticities up to 10^(-3) in the most extreme case). The effect of quadrupolar magnetic fields is comparable to that of dipolar components. A magnetic field permeating the whole star is normally needed to obtain negative quadrupole ellipticities, while fields confined to the crust typically produce positive quadrupole ellipticities.

We propose a new method to study the quasi-normal modes of rotating relativistic stars. Oscillations are treated as perturbations in the frequency domain of the stationary, axisymmetric background describing a rotating star. The perturbed quantities are expanded in circular harmonics, and the resulting 2D-equations they satisfy are integrated using spectral methods in the (r,theta)-plane. The asymptotic conditions at infinity, needed to find the mode frequencies, are implemented by generalizing the standing wave boundary condition commonly used in the non rotating case. As a test, the method is applied to find the quasi-normal mode frequencies of a slowly rotating star.

In this article we study the possibility that, due to non-linear couplings, unstable g-modes associated to convective motions excite stable oscillating g-modes. This problem is of particular interest, since gravitational waves emitted by a newly born proto-neutron star pulsating in its stable g-modes would be in the bandwidth of VIRGO and LIGO. Our results indicate that nonlinear saturation of unstable modes occurs at relatively low amplitudes, and therefore, even if there exists a coupling between stable and unstable modes, it does not seem to be sufficiently effective to explain, alone, the excitation of the oscillating g-modes found in hydrodynamical simulations.

We discuss the possibility that the detection of gravitational waves emitted by compact stars may allow to constrain the MIT bag model of quark matter equation of state. Our results show that the combined knowledge of the frequency of the emitted gravitational wave and of the mass, or the radiation radius, of the source allows one to discriminate between strange stars and neutron stars and set stringent bounds on the bag constants.

We derive the relativistic equations for stellar perturbations, including in a consistent way shear viscosity in the stress-energy tensor, and we numerically integrate our equations in the case of large viscosity. We consider the slow rotation approximation, and we neglect the coupling between polar and axial perturbations. In our approach, the frequency and damping time of the emitted gravitational radiation are directly obtained. We find that, approaching the inviscid limit from the finite viscosity case, the continuous spectrum is regularized. Constant density stars, polytropic stars, and stars with realistic equations of state are considered. In the case of constant density stars and polytropic stars, our results for the viscous damping times agree, within a factor two, with the usual estimates obtained by using the eigenfunctions of the inviscid limit. For realistic neutron stars, our numerical results give viscous damping times with the same dependence on mass and radius as previously estimated, but systematically larger of about 60%.

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

In this paper we compute the gravitational signal emitted when a white dwarf moves around a black hole on a closed or open orbit using the affine model approach. We compare the orbital and the tidal contributions to the signal, assuming that the star moves in a safe region where, although very close to the black hole, the strength of the tidal interaction is insufficient to provoque the stellar disruption. We show that for all considered orbits the tidal signal presents sharp peaks corresponding to the excitation of the star non radial oscillation modes, the amplitude of which depends on how deep the star penetrates the black hole tidal radius and on the type of orbit. Further structure is added to the emitted signal by the coupling between the orbital and the tidal motion.

We construct models of rotating stars using the perturbative approach introduced by J. Hartle in 1967, and a set of equations of state proposed to model hadronic interactions in the inner core of neutron stars. We integrate the equations of stellar structure to third order in the angular velocity and show, comparing our results to those obtained with fully non linear codes, to what extent third order corrections are needed to accurately reproduce the moment of inertia of a star which rotates at rates comparable to that of the fastest isolated pulsars.

Within a fully relativistic framework, we derive and solve numerically the perturbation equations of relativistic stars, including the stresses produced by a non-vanishing shear viscosity in the stress-energy tensor. With this approach, the real and imaginary parts of the frequency of the modes are consistently obtained. We find that, approaching the inviscid limit from the finite viscosity case, the continuous spectrum is regularized and we can calculate the quasi-normal modes for stellar models that do not admit solutions at first order in perturbation theory when the coupling between the polar and axial perturbations is neglected. The viscous damping time is found to agree within factor 2 with the usual estimate obtained by using the eigenfunctions of the inviscid limit and some approximation for the energy dissipation integrals. We find that the frequencies and viscous damping times for relativistic $r-$modes lie between the Newtonian and Cowling results. We compare the results obtained with homogeneous, polytropic and realistic equations of state and find that the frequencies depend only on the rotation rate and on the compactness parameter (M/R), being almost independent of the equation of state. Our numerical results for realistic neutron stars give viscous damping times with the same dependence on mass and radius as previously estimated, but systematically larger of about 60%.

We study the effects of rotation on the quasi normal modes (QNMs) of a newly born proto neutron star (PNS) at different evolutionary stages, until it becomes a cold neutron star (NS). We use the Cowling approximation, neglecting spacetime perturbations, and consider different models of evolving PNS. The frequencies of the modes of a PNS are considerably lower than those of a cold NS, and are further lowered by rotation; consequently, if QNMs were excited in a sufficiently energetic process, they would radiate waves that could be more easily detectable by resonant-mass and interferometric detectors than those emitted by a cold NS. We find that for high rotation rates, some of the g-modes become unstable via the CFS instability; however, this instability is likely to be suppressed by competing mechanisms before emitting a significant amount of gravitational waves.

The frequencies and damping times of the non radial oscillations of neutron stars are computed for a set of recently proposed equations of state (EOS) which describe matter at supranuclear densites. These EOS are obtained within two different approaches, the nonrelativistic nuclear many-body theory and the relativistic mean field theory, that model hadronic interactions in different ways leading to different composition and dynamics. Being the non radial oscillations associated to the emission of gravitational waves, we fit the eigenfrequencies of the fundamental mode and of the first pressure and gravitational-wave mode (polar and axial) with appropriate functions of the mass and radius of the star, comparing the fits, when available, with those obtained by Andersson and Kokkotas in 1998. We show that the identification in the spectrum of a detected gravitational signal of a sharp pulse corresponding to the excitation of the fundamental mode or of the first p-mode, combined with the knowledge of the mass of the star - the only observable on which we may have reliable information - would allow to gain interesting information on the composition of the inner core. We further discuss the detectability of these signals by gravitational detectors.

In this paper we study the effects of rotation on the frequencies of the quasi-normal modes of a proto-neutron star (PNS) born in a gravitational collapse during the first minute of life. Our analysis starts a few tenths of seconds after the PNS formation, when the stellar evolution can be described by a sequence of equilibrium configurations. We use the evolutionary models developed by Pons et al. (1999; 2001) that describe how a non rotating star cools down and contracts while neutrino diffusion and thermalization processes dominate the stellar dynamics. For assigned values of the evolution time, we set the star into slow rotation and integrate the equations of stellar perturbations in the Cowling approximation, both in the time domain and in the frequency domain, to find the quasi-normal mode frequencies. We study the secular instability of the g-modes, that are present in the oscillation spectrum due to the intense entropy and composition gradients that develop in the stellar interior, and we provide an estimate of the growth time of the unstable modes based on a post-Newtonian formula.

We study the gravitational radiation associated to the non--radial oscillations of newly born, hot neutron stars. The frequencies and damping times of the relevant quasi--normal modes are computed for two different models of proto--neutron stars, at different times of evolution, from its birth until it settles down as a cold neutron star. We find that the oscillation properties of proto--neutron stars are remarkably different from those of their cold, old descendants and that this affects the characteristic features of the gravitational signal emitted during the post-collapse evolution. The consequences on the observability of these signals by resonant--mass and interferometric detectors are analyzed. We find that gravitational waves from the pulsations of a newborn proto--neutron star in the galaxy could be detected with a signal to noise ratio of 5 by the first generation interferometers, if the energy stored in the modes is greater than $\sim 10^{-8} M_\odot c^2$, or by a resonant antenna if it is greater than $\sim 10^{-4} M_\odot c^2$. In addition since at early times the frequency of the spacetime modes is much lower than that of a cold neutron star, they would be also detectable with the same signal to noise ratio if a comparable amount of energy is radiated into these modes.

We compute the overlap function between Post-Newtonian (PN) templates and gravitational signals emitted by binary systems composed of one neutron star and one point mass, obtained by a perturbative approach. The calculations are performed for different stellar models and for different detectors, to estimate how effectual and faithful the PN templates are, and to establish whether effects related to the internal structure of neutron stars may possibly be extracted by the matched filtering technique.

A phase transition occurring in the inner core of a neutron star could be associated to a density discontinuity that would affect the frequency spectrum of the non-radial oscillation modes in two ways. Firstly, it would produce a softening of the equation of state, leading to more compact equilibrium configurations and changing the frequency of the fundamental and pressure modes of the neutron star. Secondly, a new non-zero frequency g-- mode would appear, associated to each discontinuity. These discontinuity g--modes have typical frequencies larger than those of g--modes previously studied in the literature (thermal, core g-- modes, or g--modes due to chemical inhomogeneities in the outer layers), and smaller than that of the fundamental mode; therefore they should be distinguishable from the other modes of non radial oscillation. In this paper we investigate how high density discontinuities change the frequency spectrum of the non-radial oscillations, in the framework of the general relativistic theory of stellar perturbations. Our purpose is to understand whether a gravitational signal, emitted at the frequencies of the quasi normal modes, may give some clear information on the equation of state of the neutron star and, in particular, on the parameters that characterize the density discontinuity. We discuss some astrophysical processes that may be associated to the excitation of these modes, and estimate how much gravitational energy should the modes convey to produce a signal detectable by high frequency gravitational detectors.

The effects that the structure of a neutron star would have on the gravitational emission of a binary system are studied in a perturbative regime, and in the frequency domain. Assuming that a neutron star is perturbed by a point mass moving on a close, circular orbit, we solve the equations of stellar perturbations in general relativity to evaluate the energy lost by the system in gravitational waves. We compare the energy output obtained for different stellar models with that found by assuming that the perturbed object is a black hole with the same mass, and we discuss the role played by the excitation of the stellar modes. Ouresults indicate that the stellar structure begins to affect the emitted power when the orbital velocity is v >0.2c (about 185 Hz for a binary system composed of two canonical neutron stars). We show that the differences between different stellar models and a black hole are due mainly to the excitation of the quasinormal modes of the star. Finally, we discuss to what extent and up to which distance the perturbative approach can be used to describe the interaction of a star and a pointlike massive body.

Gravitational waves transport very detailed information on the structure and evolution of astrophysical sources. For instance a binary system in the early stages of its evolution emits a wavetrain at specific frequencies that depend on the characteristics of the obital motion; as the orbit shrinks and circularize, due to radiation reaction effects, the orbiting bodies get closer and tidally interact. This interaction may result in the excitation of the proper modes of oscillation of the stars, and in the emission of gravitational signals that carry information on the mode frequencies, and consequently on the equation of state in the stellar interior. These phenomena may occur either in solar type stars with orbiting planets and in compact binaries, and in this lecture we will discuss different approaches that can be used to study these processes in the framework of General Relativity.

We compute the energy spectra of the gravitational signals emitted when a pointlike mass moves on a closed orbit around a non rotating neutron star, inducing a perturbation of its gravitational field and its internal structure. The Einstein equations and the hydrodynamical equations are perturbed and numerically integrated in the frequency domain. The results are compared with the energy spectra computed by the quadrupole formalism which assumes that both masses are pointlike, and accounts only for the radiation emitted because the orbital motion produces a time dependent quadrupole moment. The results of our perturbative approach show that, in general, the quadrupole formalism overestimates the amount of emitted radiation, especially when the two masses are close. However, if the pointlike mass is allowed to move on an orbit so tight that the keplerian orbital frequency resonates with the frequency of the fundamental quasi-normal mode of the star (2w_K=w_f), this mode can be excited and the emitted radiation can be considerably larger than that computed by the quadrupole approach.

The possibility of exciting the g-modes of a solar-type star as a consequence of the gravitational interaction with a close companion (a planet or a brown dwarf) is studied by a perturbative approach. The amplitude of the emitted gravitational wave is computed and compared with the quadrupole emission of the system, showing that in some cases it can be considerably larger. The effects of radiation reaction are considered to evaluate the timescale of the emission process, and a Roche lobe analysis is used to establish the region where the companion can orbit without being disrupted by tidal interactions with the star.

In this paper we consider the Extra-solar Planetary Systems recently discovered in our Galaxy as potential sources of gravitational waves. We estimate the frequency and characteristic amplitude of the radiation they emit due to the orbital motion, using the quadrupole formalism. In addition, we check whether the conditions needed for the resonant excitation of the f- and g-modes of the central star can be fulfilled. By a Roche-lobe analysis, we show that there could exist systems in which the low-order g-modes could be excited, although this does not happen in the systems discovered up to now.

The excitation of the axial quasi-normal modes of a relativistic star by scattered particles is studied by evolving the time dependent perturbation equations. This work is the first step towards the understanding of more complicated perturbative processes, like the capture or the scattering of particles by rotating stars. In addition, it may serve as a test for the results of the full nonlinear evolution of binary systems.

In this paper we study the perturbations of the charged, dilaton black hole, described by the solution of the low energy limit of the superstring action found by Garfinkle, Horowitz and Strominger. We compute the complex frequencies of the quasi-normal modes of this black hole, and compare the results with those obtained for a Reissner-Nordstr\"{o}m and a Schwarzschild black hole. The most remarkable feature which emerges from this study is that the presence of the dilaton breaks the \emph{isospectrality} of axial and polar perturbations, which characterizes both Schwarzschild and Reissner-Nordstr\"{o}m black holes.

We consider gravitational waves emitted by various populations of compact binaries at cosmological distances. We use population synthesis models to characterize the properties of double neutron stars, double black holes and double white dwarf binaries as well as white dwarf-neutron star, white dwarf-black hole and black hole-neutron star systems. We use the observationally determined cosmic star formation history to reconstruct the redshift distribution of these sources and their merging rate evolution. The gravitational signals emitted by each source during its early-inspiral phase add randomly to produce a stochastic background in the low frequency band with spectral strain amplitude between 10^{-18} Hz^{-1/2} and 5 10^{-17} Hz^{-1/2} at frequencies in the interval [5 10^{-6}-5 10^{-5}] Hz. The overall signal which, at frequencies above 10^{-4}Hz, is largely dominated by double white dwarf systems, might be detectable with LISA in the frequency range [1-10] mHz and acts like a confusion limited noise component which might limit the LISA sensitivity at frequencies above 1 mHz.

Stars and black holes are sources of gravitational radiation in many phases of their life, and the signals they emit exhibit features that are characteristic of the generating process. Emitted since the beginning of star formation, these signals also contribute to create a stochastic background of gravitational waves. We shall show how the spectral properties of this background can be estimated in terms of the energy spectrum of each single event and of the star formation rate history, which is now deducible from astronomical observations. We shall further discuss the process of scattering of masses by stars and black holes, showing that, unlike black holes, stars emit signals that carry a clear signature of the nature of the source.

We study the gravitational wave emission from the first stars which are assumed to be Very Massive Objects (VMOs). We take into account various feedback (both radiative and stellar) effects regulating the collapse of objects in the early universe and thus derive the VMO initial mass function and formation rate. If the final fate of VMOs is to collapse, leaving very massive black hole remnants, then the gravitational waves emitted during each collapse would be seen as a stochastic background. The predicted spectral strain amplitude in a critical density Cold Dark Matter universe peaks in the frequency range \approx 5 \times 10^{-4}-5 \times 10^{-3} Hz where it has a value in the range \approx 10^{-20}-10^{-19} Hz^{-1/2}, and might be detected by LISA. The expected emission rate is roughly 4000 events/yr, resulting in a stationary, discrete sequence of bursts, i.e. a shot--noise signal.

Astrophysical sources of gravitational radiation are likely to have been formed since the beginning of star formation. Realistic source rates of formation throughout the Universe have been estimated from an observation-based determination of the star formation rate density evolution. Both the radiation emitted during the collapse to black holes and the spin-down radiation, induced by the r-mode instability, emitted by hot, young rapidly rotating neutron stars have been considered. We have investigated the overall signal produced by the ensemble of sources exploring the parameter space and discussing its possible detectability.

We estimate the spectral properties of the stochastic background of gravitational radiation emitted by a cosmological population of hot, young, rapidly rotating neutron stars. Their formation rate as a function of redshift is deduced from an observation-based determination of the star formation history in the Universe, and the gravitational energy is assumed to be radiated during the spin-down phase associated to the newly discovered r-mode instability. We calculate the overall signal produced by the ensemble of such neutron stars, assuming various cosmological backgrounds. We find that the spectral strain amplitude has a maximum $\approx (2-4)\times 10^{-26} {Hz}^{-1/2}$, at frequencies $\approx (30-60)$ Hz, while the corresponding closure density, $h^2 \Omega_{GW}$, has a maximum amplitude plateau of $\approx (2.2-3.3) \times 10^{-8}$ in the frequency range $(500-1700)$ Hz. We compare our results with a preliminary analysis done by Owen et al. (1998), and discuss the detectability of this background.

We analyse the stochastic background of gravitational radiation emitted by a cosmological population of core-collapse supernovae. The supernova rate as a function of redshift is deduced from an observation-based determination of the star formation rate density evolution. We then restrict our analysis to the range of progenitor masses leading to black hole collapse. In this case, the main features of the gravitational-wave emission spectra have been shown to be, to some extent, independent of the initial conditions and of the equation of state of the collapsing star, and to depend only on the black hole mass and angular momentum. We calculate the overall signal produced by the ensemble of black-hole collapses throughout the Universe, assuming a flat cosmology with vanishing cosmological constant. Within a wide range of parameter values, we find that the spectral strain amplitude has a maximum at a few hundred Hz with an amplitude between $10^{-28}$ and $10^{-27} Hz^{-1/2}$; the corresponding closure density, $\Omega_{GW}$, has a maximum amplitude ranging between $10^{-11}$ and $10^{-10}$ in the frequency interval $\sim 1.5-2.5$ kHz. Contrary to previous claims, our observation-based determination leads to a duty cycle of order 0.01, making our stochastic backgound a non-continuous one. Although the amplitude of our background is comparable to the sensitivity that can be reached by a pair of advanced LIGO detectors, the characteristic shot-noise structure of the predicted signal might be in principle exploited to design specific detection strategies.