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.

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.

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.

Gravitational-wave (GW) observations provide a wealth of information on the nature and properties of black holes. Among these, tidal Love numbers or the multipole moments of the inspiralling and final objects are key to a number of constraints. Here, we consider these observations in the context of higher-dimensional scenarios, with flat large extra dimensions. We show that -- as might be anticipated, but not always appreciated in the literature -- physically motivated set-ups are unconstrained by gravitational-wave data. Dynamical processes that do not excite the Kaluza-Klein (KK) modes lead to a signal identical to that in four-dimensional general relativity in vacuum . In addition, any possible excitation of the KK modes is highly suppressed relative to the dominant quadrupolar term; given existing constraints on the extra dimensions and the masses of the objects seen in gravitational-wave observations, KK modes appear at post-Newtonian order $\sim 10^{11}$. Finally, we re-compute the tidal Love numbers of spherical black holes in higher dimensions. We confirm that these are different from zero, but comparing with previous computations we find a different magnitude and sign.

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.

Scalar-tensor theories of gravity where a new scalar degree of freedom couples to the Gauss-Bonnet invariant can exhibit the phenomenon of spontaneous black hole scalarization. These theories admit both the classic black holes predicted by general relativity as well as novel hairy black hole solutions. The stability of the hairy black holes is strongly dependent on the precise form of the scalar-gravity coupling. A radial stability investigation revealed that all scalarized black hole solutions are unstable when the coupling between the scalar field and the Gauss-Bonnet invariant is quadratic in the scalar, whereas stable solutions exist for exponential couplings. Here we elucidate this behavior. We demonstrate that, while the quadratic term controls the onset of the tachyonic instability that gives rise to the black hole hair, the higher-order coupling terms control the nonlinearities that quench that instability, and hence also control the stability of the hairy black hole solutions.

It has recently been shown that nontrivial couplings between a scalar and the Gauss-Bonnet invariant can give rise to black hole spontaneous scalarization. Theories that exhibit this phenomenon are among the leading candidates for testing gravity with upcoming black hole observations. All models considered so far have focused on specific forms for the coupling, neglecting scalar self-interactions. In this work, we take the first steps towards placing this phenomenon on a more robust theoretical footing by considering the leading-order scalar self-interactions as well as the scalar-Gauss-Bonnet coupling. Our approach is consistent with the principles of effective field theory and yields the simplest and most natural model. We find that a mass term for the scalar alters the threshold for the onset of scalarization, and we study the mass range over which scalarized black hole solutions exist. We also demonstrate that the quartic self-coupling is sufficient to produce scalarized solutions that are stable against radial perturbations, without the need to resort to higher-order terms in the Gauss-Bonnet coupling function. Our model therefore represents a canonical model that can be studied further, with the ultimate aim of developing falsifiable tests of black hole scalarization.

In general relativity, Maxwell's equations are embedded in curved spacetime through the minimal prescription, but this could change if strong-gravity modifications are present. We show that with a nonminimal coupling between gravity and a massless vector field, nonperturbative effects can arise in compact stars. We find solutions describing stars with nontrivial vector field configurations, some of which are associated to an instability, while others are not. The vector field can be interpreted either as the electromagnetic field, or as a hidden vector field weakly coupled with the standard model.

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 brief note, we clarify certain aspects related to the magnetic (i.e., odd parity or axial) tidal Love numbers of a star in general relativity. Magnetic tidal deformations of a compact star had been computed in 2009 independently by Damour and Nagar and by Binnington and Poisson. More recently, Landry and Poisson showed that the magnetic tidal Love numbers depend on the assumptions made on the fluid, in particular they are different (and of opposite sign) if the fluid is assumed to be in static equilibrium or if it is irrotational. We show that the zero-frequency limit of the Regge-Wheeler equation forces the fluid to be irrotational. For this reason, the results of Damour and Nagar are equivalent to those of Landry and Poisson for an irrotational fluid, and are expected to be the most appropriate to describe realistic configurations.

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.

We compute the spin-tidal couplings that affect the dynamics of two orbiting bodies at the leading order in the post-Newtonian (PN) framework and to linear order in the spin. These corrections belong to two classes: (i) terms arising from the coupling between the ordinary tidal terms and the point-particle terms, which depend on the standard tidal Love numbers of order $l$ and affect the gravitational-wave (GW) phase at $(2l+5/2)$PN order and (ii) terms depending on the rotational tidal Love numbers, recently introduced in previous work, that affect the GW phase at $(2l+1/2+\delta_{2l})$PN order. For circular orbits and spins orthogonal to the orbital plane, all leading-order spin-tidal terms enter the GW phase at $1.5$PN order relative to the standard, quadrupolar, tidal deformability term (and, thus, before the standard octupolar tidal deformability terms). We present the GW phase that includes all tidal terms up to $6.5$PN order and to linear order in the spin. We comment on a conceptual issue related to the inclusion of the rotational tidal Love numbers in a Lagrangian formulation and on the relevance of spin-tidal couplings for parameter estimation in coalescing neutron-star binaries and for tests of gravity.

We study the nonlinear dynamics of black holes that carry scalar hair and binaries composed of such black holes. The scalar hair is due to a linear or exponential coupling between the scalar and the Gauss--Bonnet invariant. We work perturbatively in the coupling constant of that interaction but nonperturbatively in the fields. We first consider the dynamical formation of hair for isolated black holes of arbitrary spin and determine the final state. This also allows us to compute for the first time the scalar quasinormal modes of rotating black holes in the presence of this coupling. We then study the evolution of nonspinning black-hole binaries with various mass ratios and produce the first scalar waveform for a coalescence. An estimate of the energy loss in scalar radiation and the effect this has on orbital dynamics and the phase of the GWs (entering at quadratic order in the coupling) show that GW detections can set the most stringent constraint to date on theories that exhibit a coupling between a scalar field and the Gauss--Bonnet invariant.

We identify a class of scalar-tensor theories with coupling between the scalar and the Gauss-Bonnet invariant that exhibit spontaneous scalarization for both black holes and compact stars. In particular, these theories formally admit all of the stationary solutions of general relativity, but these are not dynamically preferred if certain conditions are satisfied. Remarkably, black holes exhibit scalarization if their mass lies within one of many narrow bands. We find evidence that scalarization can occur in neutron stars as well.

With the doors beginning to swing open on the new gravitational wave astronomy, this review provides an up-to-date survey of the most important physical mechanisms that could lead to emission of potentially detectable gravitational radiation from isolated and accreting neutron stars. In particular we discuss the gravitational wave-driven instability and asteroseismology formalism of the f- and r-modes, the different ways that a neutron star could form and sustain a non-axisymmetric quadrupolar "mountain" deformation, the excitation of oscillations during magnetar flares and the possible gravitational wave signature of pulsar glitches. We focus on progress made in the recent years in each topic, make a fresh assessment of the gravitational wave detectability of each mechanism and, finally, highlight key problems and desiderata for future work.

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.

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.

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.

Black holes in General Relativity are very simple objects. This property, that goes under the name of "no-hair," has been refined in the last few decades and admits several versions. The simplicity of black holes makes them ideal testbeds of fundamental physics and of General Relativity itself. Here we discuss the no-hair property of black holes, how it can be measured in the electromagnetic or gravitational window, and what it can possibly tell us about our universe.

We review recent progress in the application of numerical relativity techniques to astrophysics and high-energy physics. We focus on some developments that took place within the "Numerical Relativity and High Energy Physics" network, a Marie Curie IRSES action that we coordinated, namely: spin evolution in black hole binaries, high-energy black hole collisions, compact object solutions in scalar-tensor gravity, superradiant instabilities and hairy black hole solutions in Einstein's gravity coupled to fundamental fields, and the possibility to gain insight into these phenomena using analog gravity models.

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.

One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime.

The demand to obtain answers to highly complex problems within strong-field gravity has been met with significant progress in the numerical solution of Einstein's equations - along with some spectacular results - in various setups. We review techniques for solving Einstein's equations in generic spacetimes, focusing on fully nonlinear evolutions but also on how to benchmark those results with perturbative approaches. The results address problems in high-energy physics, holography, mathematical physics, fundamental physics, astrophysics and cosmology.

Gravitational theories with multiple scalar fields coupled to the metric and each other --- a natural extension of the well studied single-scalar-tensor theories --- are interesting phenomenological frameworks to describe deviations from general relativity in the strong-field regime. In these theories, the $N$-tuple of scalar fields takes values in a coordinate patch of an $N$-dimensional Riemannian target-space manifold whose properties are poorly constrained by weak-field observations. Here we introduce for simplicity a non-trivial model with two scalar fields and a maximally symmetric target-space manifold. Within this model we present a preliminary investigation of spontaneous scalarization for relativistic, perfect fluid stellar models in spherical symmetry. We find that the scalarization threshold is determined by the eigenvalues of a symmetric scalar-matter coupling matrix, and that the properties of strongly scalarized stellar configurations additionally depend on the target-space curvature radius. In preparation for numerical relativity simulations, we also write down the $3+1$ decomposition of the field equations for generic tensor-multi-scalar theories.

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.

We study non-radial oscillations of neutron stars with superfluid baryons, in a general relativistic framework, including finite temperature effects. Using a perturbative approach, we derive the equations describing stellar oscillations, which we solve by numerical integration, employing different models of nucleon superfluidity, and determining frequencies and gravitational damping times of the quasi-normal modes. As expected by previous results, we find two classes of modes, associated to superfluid and non-superfluid degrees of freedom, respectively. We study the temperature dependence of the modes, finding that at specific values of the temperature, the frequencies of the two classes of quasi-normal modes show avoided crossings, and their damping times become comparable. We also show that, when the temperature is not close to the avoided crossings, the frequencies of the modes can be accurately computed by neglecting the coupling between normal and superfluid degrees of freedom. Our results have potential implications on the gravitational wave emission from neutron stars.

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.

Scalar-tensor theories are a compelling alternative to general relativity and one of the most accepted extensions of Einstein's theory. Black holes in these theories have no hair, but could grow "wigs" supported by time-dependent boundary conditions or spatial gradients. Time-dependent or spatially varying fields lead in general to nontrivial black hole dynamics, with potentially interesting experimental consequences. We carry out a numerical investigation of the dynamics of single and binary black holes in the presence of scalar fields. In particular we study gravitational and scalar radiation from black-hole binaries in a constant scalar-field gradient, and we compare our numerical findings to analytical models. In the single black hole case we find that, after a short transient, the scalar field relaxes to static configurations, in agreement with perturbative calculations. Furthermore we predict analytically (and verify numerically) that accelerated black holes in a scalar-field gradient emit scalar radiation. For a quasicircular black-hole binary, our analytical and numerical calculations show that the dominant component of the scalar radiation is emitted at twice the binary's orbital frequency.

We analyze damping of oscillations of general relativistic superfluid neutron stars. To this aim we extend the method of decoupling of superfluid and normal oscillation modes first suggested in [Gusakov & Kantor PRD 83, 081304(R) (2011)]. All calculations are made self-consistently within the finite temperature superfluid hydrodynamics. The general analytic formulas are derived for damping times due to the shear and bulk viscosities. These formulas describe both normal and superfluid neutron stars and are valid for oscillation modes of arbitrary multipolarity. We show that: (i) use of the ordinary one-fluid hydrodynamics is a good approximation, for most of the stellar temperatures, if one is interested in calculation of the damping times of normal f-modes; (ii) for radial and p-modes such an approximation is poor; (iii) the temperature dependence of damping times undergoes a set of rapid changes associated with resonance coupling of neighboring oscillation modes. The latter effect can substantially accelerate viscous damping of normal modes in certain stages of neutron-star thermal evolution.

We discuss a general method to study linear perturbations of slowly rotating black holes which is valid for any perturbation field, and particularly advantageous when the field equations are not separable. As an illustration of the method we investigate massive vector (Proca) perturbations in the Kerr metric, which do not appear to be separable in the standard Teukolsky formalism. Working in a perturbative scheme, we discuss two important effects induced by rotation: a Zeeman-like shift of nonaxisymmetric quasinormal modes and bound states with different azimuthal number m, and the coupling between axial and polar modes with different multipolar index l. We explicitly compute the perturbation equations up to second order in rotation, but in principle the method can be extended to any order. Working at first order in rotation we show that polar and axial Proca modes can be computed by solving two decoupled sets of equations, and we derive a single master equation describing axial perturbations of spin s=0 and s=+-1. By extending the calculation to second order we can study the superradiant regime of Proca perturbations in a self-consistent way. For the first time we show that Proca fields around Kerr black holes exhibit a superradiant instability, which is significantly stronger than for massive scalar fields. Because of this instability, astrophysical observations of spinning black holes provide the tightest upper limit on the mass of the photon: mv<4x10^-20 eV under our most conservative assumptions. Spin measurements for the largest black holes could reduce this bound to mv<10^-22 eV or lower.

Generic extensions of the standard model predict the existence of ultralight bosonic degrees of freedom. Several ongoing experiments are aimed at detecting these particles or constraining their mass range. Here we show that massive vector fields around rotating black holes can give rise to a strong superradiant instability which extracts angular momentum from the hole. The observation of supermassive spinning black holes imposes limits on this mechanism. We show that current supermassive black hole spin estimates provide the tightest upper limits on the mass of the photon (mv<4x10^{-20} eV according to our most conservative estimate), and that spin measurements for the largest known supermassive black holes could further lower this bound to mv<10^{-22} eV. Our analysis relies on a novel framework to study perturbations of rotating Kerr black holes in the slow-rotation regime, that we developed up to second order in rotation, and that can be extended to other spacetime metrics and other theories.

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.

Scalar-tensor theories are among the simplest extensions of general relativity. In theories with light scalars, deviations from Einstein's theory of gravity are determined by the scalar mass m_s and by a Brans-Dicke-like coupling parameter \omega_{BD}. We show that gravitational-wave observations of nonspinning neutron star-black hole binary inspirals can be used to set lower bounds on \omega_{BD} and upper bounds on the combination m_s/\sqrt{\omega_{BD}}$. We estimate via a Fisher matrix analysis that individual observations with signal-to-noise ratio \rho would yield (m_s/\sqrt{\omega_{BD}})(\rho/10)<10^{-15}, 10^{-16} and 10^{-19} eV for Advanced LIGO, ET and eLISA, respectively. A statistical combination of multiple observations may further improve these bounds.

Nonlinear dynamics in cosmological backgrounds has the potential to teach us immensely about our universe, and also to serve as prototype for nonlinear processes in generic curved spacetimes. Here we report on dynamical evolutions of black holes in asymptotically de Sitter spacetimes. We focus on the head-on collision of equal mass binaries and for the first time compare analytical and perturbative methods with full blown nonlinear simulations. Our results include an accurate determination of the merger/scatter transition (consequence of an expanding background) for small mass binaries and a test of the Cosmic Censorship conjecture, for large mass binaries. We observe that, even starting from small separations, black holes in large mass binaries eventually lose causal contact, in agreement with the conjecture.

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 study the coupling of massive scalar fields to matter in orbit around rotating black holes. It is generally expected that orbiting bodies will lose energy in gravitational waves, slowly inspiralling into the black hole. Instead, we show that the coupling of the field to matter leads to a surprising effect: because of superradiance, matter can hover into "floating orbits" for which the net gravitational energy loss at infinity is entirely provided by the black hole's rotational energy. Orbiting bodies remain floating until they extract sufficient angular momentum from the black hole, or until perturbations or nonlinear effects disrupt the orbit. For slowly rotating and nonrotating black holes floating orbits are unlikely to exist, but resonances at orbital frequencies corresponding to quasibound states of the scalar field can speed up the inspiral, so that the orbiting body "sinks". These effects could be a smoking gun of deviations from general relativity.

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.

Dynamical Chern-Simons gravity is an interesting extension of General Relativity, which finds its way in many different contexts, including string theory, cosmological settings and loop quantum gravity. In this theory, the gravitational field is coupled to a scalar field by a parity-violating term, which gives rise to characteristic signatures. Here we investigate how Chern-Simons gravity would affect the quasi-circular inspiralling of a small, stellar-mass object into a large non-rotating supermassive black hole, and the accompanying emission of gravitational and scalar waves. We find the relevant equations describing the perturbation induced by the small object, and we solve them through the use of Green's function techniques. Our results show that for a wide range of coupling parameters, the Chern-Simons coupling gives rise to an increase in total energy flux, which translates into a fewer number of gravitational-wave cycles over a certain bandwidth. For space-based gravitational-wave detectors such as LISA, this effect can be used to constrain the coupling parameter effectively.

Black objects in higher dimensional space-times have a remarkably richer structure than their four dimensional counterparts. They appear in a variety of configurations (e.g. black holes, black branes, black rings, black Saturns), and display complex stability phase diagrams. They might also play a key role in high energy physics: for energies above the fundamental Planck scale, gravity is the dominant interaction which, together with the hoop-conjecture, implies that the trans-Planckian scattering of point particles should be well described by black hole scattering. Higher dimensional scenarios with a fundamental Planck scale of the order of TeV predict, therefore, black hole production at the LHC, as well as in future colliders with yet higher energies. In this setting, accurate predictions for the production cross-section and energy loss (through gravitational radiation) in the formation of black holes in parton-parton collisions is crucial for accurate phenomenological modelling in Monte Carlo event generators. In this paper, we use the formalism and numerical code reported in arXiv:1001.2302 to study the head-on collision of two black holes. For this purpose we provide a detailed treatment of gravitational wave extraction in generic D-dimensional space-times, which uses the Kodama-Ishibashi formalism. For the first time, we present the results of numerical simulations of the head-on collision in five space-time dimensions, together with the relevant physical quantities. We show that the total radiated energy, when two black holes collide from rest at infinity, is approximately (0.089\pm 0.006)% of the centre of mass energy, slightly larger than the 0.055% obtained in the four dimensional case, and that the ringdown signal at late time is in very good agreement with perturbative calculations.

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.

Dynamical Chern-Simons gravity is an extension of General Relativity in which the gravitational field is coupled to a scalar field through a parity-violating Chern-Simons term. In this framework, we study perturbations of spherically symmetric black hole spacetimes, assuming that the background scalar field vanishes. Our results suggest that these spacetimes are stable, and small perturbations die away as a ringdown. However, in contrast to standard General Relativity, the gravitational waveforms are also driven by the scalar field. Thus, the gravitational oscillation modes of black holes carry imprints of the coupling to the scalar field. This is a smoking gun for Chern-Simons theory and could be tested with gravitational-wave detectors, such as LIGO or LISA. For negative values of the coupling constant, ghosts are known to arise, and we explicitly verify their appearance numerically. Our results are validated using both time evolution and frequency domain methods.

The numerical evolution of Einstein's field equations in a generic background has the potential to answer a variety of important questions in physics: from applications to the gauge-gravity duality, to modelling black hole production in TeV gravity scenarios, analysis of the stability of exact solutions and tests of Cosmic Censorship. In order to investigate these questions, we extend numerical relativity to more general space-times than those investigated hitherto, by developing a framework to study the numerical evolution of D dimensional vacuum space-times with an SO(D-2) isometry group for D\ge 5, or SO(D-3) for D\ge 6. Performing a dimensional reduction on a (D-4)-sphere, the D dimensional vacuum Einstein equations are rewritten as a 3+1 dimensional system with source terms, and presented in the Baumgarte, Shapiro, Shibata and Nakamura (BSSN) formulation. This allows the use of existing 3+1 dimensional numerical codes with small adaptations. Brill-Lindquist initial data are constructed in D dimensions and a procedure to match them to our 3+1 dimensional evolution equations is given. We have implemented our framework by adapting the LEAN code and perform a variety of simulations of non-spinning black hole space-times. Specifically, we present a modified moving puncture gauge which facilitates long term stable simulations in D=5. We further demonstrate the internal consistency of the code by studying convergence and comparing numerical versus analytic results in the case of geodesic slicing for D=5,6.

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

It has been suggested that rotating black holes could serve as particle colliders with arbitrarily high center-of-mass energy. Astrophysical limitations on the maximal spin, back-reaction effects and sensitivity to the initial conditions impose severe limits on the likelihood of such collisions.

Dynamical Chern-Simons (DCS) modified gravity is an attractive, yet relatively unexplored, candidate to an alternative theory of gravity. The DCS correction couples a dynamical scalar field to the gravitational field. In this framework, we analyze the perturbation formalism and stability properties of spherically symmetric black holes. Assuming that no background scalar field is present, gravitational perturbations with polar and axial parities decouple. We find no effect of the Chern-Simons coupling on the polar sector, while axial perturbations couple to the Chern-Simons scalar field. The axial sector can develop strong instabilities if the coupling parameter beta, associated to the dynamical coupling of the scalar field, is small enough; this yields a constraint on beta which is much stronger than the constraints previously known in the literature.

We study the transformation of multipolar decompositions of gravitational radiation under rotations and boosts. Rotations to the remnant black hole's frame simplify the waveforms from the merger of generic spinning black hole binaries. Boosts may be important to get an accurate gravitational-wave phasing, especially for configurations leading to large recoil velocities of the remnant. As a test of our formalism we revisit the classic problem of point particles falling into a Schwarzschild black hole. Then we highlight by specific examples the importance of choosing the right frame in numerical simulations of unequal-mass, spinning binary black-hole mergers.

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 study equilibrium shapes, stability and possible bifurcation diagrams of fluids in higher dimensions, held together by either surface tension or self-gravity. We consider the equilibrium shape and stability problem of self-gravitating spheroids, establishing the formalism to generalize the MacLaurin sequence to higher dimensions. We show that such simple models, of interest on their own, also provide accurate descriptions of their general relativistic relatives with event horizons. The examples worked out here hint at some model-independent dynamics, and thus at some universality: smooth objects seem always to be well described by both ``replicas'' (either self-gravity or surface tension). As an example, we exhibit an instability afflicting self-gravitating (Newtonian) fluid cylinders. This instability is the exact analogue, within Newtonian gravity, of the Gregory-Laflamme instability in general relativity. Another example considered is a self-gravitating Newtonian torus made of a homogeneous incompressible fluid. We recover the features of the black ring in general relativity.

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 models with extra dimensions, a black hole evaporates both in the bulk and on the visible brane, where standard model fields live. The exact emissivities of each particle species are needed to determine how the black hole decay proceeds. We compute and discuss the absorption cross-sections, the relative emissivities and the total power output of all known fields in the evaporation phase. Graviton emissivity is highly enhanced as the spacetime dimensionality increases. Therefore, a black hole loses a significant fraction of its mass in the bulk. This result has important consequences for the phenomenology of black holes in models with extra dimensions and black hole detection in particle colliders.

We investigate the non-linear coupling between radial and non-radial oscillations of static spherically symmetric neutron stars as a possible mechanism for the generation of gravitational waves that may lead to observable signatures. In this paper we concentrate on the axial sector of the non-radial perturbations. By using a multi-parameter perturbative framework we introduce a complete description of the non-linear coupling between radial and axial non-radial oscillations; we study the gauge invariant character of the associated perturbative variables and develop a computational scheme to evolve the non-linear coupling perturbations in the time domain. We present results of simulations corresponding to different physical situations and discuss the dynamical behaviour of this non-linear coupling. Of particular interest is the occurrence of signal amplifications in the form of resonance phenomena when a frequency associated with the radial pulsations is close to a frequency associated with one of the axial w-modes of the star. Finally, we mention possible extensions of this work and improvements towards more astrophysically motivated scenarios.

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

Linear perturbation theory is appropriate to describe small oscillations of stars, while a mild non-linearity is still tractable perturbatively but requires to consider mode coupling. It is natural to start to look at this problem by considering the coupling between linear radial and non-radial modes. Radial pulsations of a spherical compact objects do not per se emit gravitational waves but, if the coupling is efficient in driving and possibly amplifying the non-radial modes, gravitational radiation could then be produced to a significant level. In this paper we develop the relativistic formalism to study the coupling of radial and non-radial first order perturbations of a compact spherical star. From a mathematical point of view, it is convenient to treat the two sets of perturbations as separately parametrized, using a 2-parameter perturbative expansion of the metric, the energy-momentum tensor and Einstein equations in which $\lambda$ is associated with the radial modes, $\epsilon$ with the non-radial perturbations, and the $\lambda\epsilon$ terms describe the coupling. This approach provides a well-defined framework to consider the gauge dependence of perturbations, allowing us to use $\epsilon$ order gauge-invariant non-radial variables on the static background and to define new second order $\lambda\epsilon$ gauge-invariant variables describing the non-linear coupling. We present the evolution and constraint equations for our variables outlining the setup for numerical computations, and briefly discuss the surface boundary conditions in terms of the second order $\lambda\epsilon$ Lagrangian pressure perturbation.

This is a progress report on our study of the coupling of first-order radial and non-radial relativistic perturbations of a static spherical star. Our goal is to investigate the effects of this coupling on the gravitational wave signal of neutron stars. In particular, we are looking for the existence of resonances and parametric amplifications, changes in the damping time of non-radial oscillations, etc. To that end, we have developed a formalism that introduces gauge invariant quantities to describe the coupling. Their equations have the same structure as the equations for first-order non-radial perturbations plus some source terms, which makes them very appealing for time domain studies.

We have developed a formalism to study non-adiabatic, non-radial oscillations of non-rotating compact stars in the frequency domain, including the effects of thermal diffusion in the framework of general relativistic perturbation theory. When a general equation of state depending on temperature is used, the perturbations of the fluid result in heat flux which is coupled with the spacetime geometry through the Einstein field equations. Our results show that the frequency of the first pressure (p) and gravity (g) oscillation modes is significantly affected by thermal diffusion, while that of the fundamental (f) mode is basically unaltered due to the global nature of that oscillation. The damping time of the oscillations is generally much smaller than in the adiabatic case (more than two orders of magnitude for the p- and g-modes) reflecting the effect of thermal dissipation. Both the isothermal and adiabatic limits are recovered in our treatment and we study in more detail the intermediate regime. Our formalism finds its natural astrophysical application in the study of the oscillation properties of newly born neutron stars, neutron stars with a deconfined quark core phase, or strange stars which are all promising sources of gravitational waves with frequencies in the band of the first generation and advanced ground-based interferometric detectors.

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.

We introduce N-parameter perturbation theory as a new tool for the study of non-linear relativistic phenomena. The main ingredient in this formulation is the use of the Baker-Campbell-Hausdorff formula. The associated machinery allows us to prove the main results concerning the consistency of the scheme to any perturbative order. Gauge transformations and conditions for gauge invariance at any required order can then be derived from a generating exponential formula via a simple Taylor expansion. We outline the relation between our novel formulation and previous developments.

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.

An implicit fundamental assumption in relativistic perturbation theory is that there exists a parametric family of spacetimes that can be Taylor expanded around a background. The choice of the latter is crucial to obtain a manageable theory, so that it is sometime convenient to construct a perturbative formalism based on two (or more) parameters. The study of perturbations of rotating stars is a good example: in this case one can treat the stationary axisymmetric star using a slow rotation approximation (expansion in the angular velocity Omega), so that the background is spherical. Generic perturbations of the rotating star (say parametrized by lambda) are then built on top of the axisymmetric perturbations in Omega. Clearly, any interesting physics requires non-linear perturbations, as at least terms lambda Omega need to be considered. In this paper we analyse the gauge dependence of non-linear perturbations depending on two parameters, derive explicit higher order gauge transformation rules, and define gauge invariance. The formalism is completely general and can be used in different applications of general relativity or any other spacetime theory.

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.

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.