Using nonlinear, fully relativistic, simulations we investigate the dynamics and gravitational-wave signature associated with instabilities in neutron star post-merger remnants. For simplified models of the remnant, we establish the presence of instability in stars with moderate $T/|W|$, the ratio between the kinetic and the gravitational potential energies. Detailed analysis of the density oscillation pattern reveals a \emph{local} instability in the inner region of the more realistic differential rotation profile. We apply Rayleigh's inflection theorem and Fj{\o}rtoft's theorem to analyze the stability criteria concluding that this \emph{inner local} instability originates from a shear instability close to the peak of the angular velocity profile and that it later evolves into a fast-rotating $m=2$ oscillation pattern. We discuss the importance of the presence of a co-rotation point in the fluid, its connection with the shear instability, and comparisons to the Rossby Wave and Papaloizou Pringle Instabilities considered in the wider literature.

Implementing a family of differential rotation laws inspired by binary neutron-star merger remnants, we consider the impact of the rotation profile on the low-T/W instability. We use time evolutions of the linearised dynamical equations, in Newtonian gravity, to study non-axisymmetric oscillations and identify the unstable modes. The presence and evolution of the low-T/W instability is monitored with the canonical energy and angular momentum, while the growth time is extracted from the evolved kinetic energy. The results for the new rotation laws highlight similarities with the commonly considered j-constant law. The instability sets in when an oscillation mode co-rotates with the star (i.e. whenever there is a point where the mode's pattern speed matches the bulk angular velocity) and grows faster deep inside the co-rotation region. However, the new profiles add features, like an additional co-rotation point to the problem, which affect the onset of instability. The rotation laws influence more drastically the oscillation frequencies of the l=m=2 f-mode in fast rotating models, but affect the instability growth time at any rotation rate. We also identify models where the low-T/W instability appears to be triggered by inertial modes. We discuss to what extent the inferred qualitative behaviour is likely to be of observational relevance.

Gravitational waves from core-collapse supernovae are produced by the excitation of different oscillation modes in the proto-neutron star (PNS) and its surroundings, including the shock. In this work we study the relationship between the post-bounce oscillation spectrum of the PNS-shock system and the characteristic frequencies observed in gravitational-wave signals from core-collapse simulations. This is a fundamental first step in order to develop a procedure to infer astrophysical parameters of the PNS formed in core-collapse supernovae. Our method combines information from the oscillation spectrum of the PNS, obtained through linear-perturbation analysis in general relativity of a background physical system, with information from the gravitational-wave spectrum of the corresponding non-linear, core-collapse simulation. Using results from the simulation of the collapse of a 35 $M_{\odot}$ presupernova progenitor we show that both types of spectra are indeed related and we are able to identify the modes of oscillation of the PNS, namely g-modes, p-modes, hybrid modes, and standing-accretion-shock-instability (SASI) modes, obtaining a remarkably close correspondence with the time-frequency distribution of the gravitational-wave modes. The analysis presented in this paper provides a proof-of-concept that asteroseismology is indeed possible in the core-collapse scenario, and it may serve as a basis for future work on PNS parameter inference based on gravitational-wave observations.

Improvements in ground-based, advanced gravitational wave (GW) detectors may allow in the near future to observe the GW signal of a nearby core-collapse supernova. For the most common type of progenitors, likely with slowly rotating cores, the dominant GW emission mechanisms are the post-bounce oscillations of the proto-neutron star (PNS) before the explosion. We present a new procedure to compute the eigenmodes of the system formed by the PNS and the stalled accretion shock in general relativity including spacetime perturbations. The new method improves on previous results by accounting for perturbations of both the lapse function and the conformal factor. We apply our analysis to two numerical core-collapse simulations and show that our improved method is able to obtain eigenfrequencies that accurately match the features observed in the GW signal and to predict the qualitative behaviour of quasi-radial oscillations. Our analysis is possible thanks to a newly developed algorithm to classify the computed eigenmodes in different classes (f-, p-, and g-modes), improving previous results which suffered from misclassification issues. We find that most of the GW energy is stored in the lowest order eigenmodes, in particular in the $^2g_1$ mode and in the $^2f$ mode. Our results also suggest that a low-frequency component of the GW signal attributed in previous works to the characteristic frequency of the Standing Accretion Shock Instability should be identified as the fundamental quadrupolar f-mode. We also develop a formalism to estimate the contribution of quasi-radial ($l=0$) modes to the quadrupolar component of the GW emission in the case of a deformed background, with application to rapidly rotating cores. This work provides further support for asteroseismology of core-collapse supernovae and the inference of PNS properties based on GW observations.

We revisit the various approximations employed to study the long-term evolution of the magnetic field in neutron star cores and discuss their limitations and possible improvements. A recent controversy on the correct form of the induction equation and the relevant evolution timescale in superconducting neutron star cores is addressed and clarified. We show that this ambiguity in the estimation of timescales arises as a consequence of nominally large terms that appear in the induction equation, but which are, in fact, mostly irrotational. This subtlety leads to a discrepancy by many orders of magnitude when velocity fields are absent or ignored. Even when internal velocity fields are accounted for, only the solenoidal part of the electric field contributes to the induction equation, which can be substantially smaller than the irrotational part. We also argue that stationary velocity fields must be incorporated in the slow evolution of the magnetic field as the next level of approximation.

We study ambipolar diffusion in strongly magnetised neutron stars, with special focus on the effects of neutrino reaction rates and the impact of a superfluid/superconducting transition in the neutron star core. For axisymmetric magnetic field configurations, we determine the deviation from $\beta-$equilibrium induced by the magnetic force and calculate the velocity of the slow, quasi-stationary, ambipolar drift. We study the temperature dependence of the velocity pattern and clearly identify the transition to a predominantly solenoidal flow. For stars without superconducting/superfluid constituents and with a mixed poloidal-toroidal magnetic field of typical magnetar strength, we find that ambipolar diffusion proceeds fast enough to have a significant impact on the magnetic field evolution only at low core temperatures, $T \lesssim 1-2\times10^8$ K. The ambipolar diffusion timescale becomes appreciably shorter when fast neutrino reactions are present, because the possibility to balance part of the magnetic force with pressure gradients is reduced. We also find short ambipolar diffusion timescales in the case of superconducting cores for $T \lesssim 10^9$ K, due to the reduced interaction between protons and neutrons. In the most favourable scenario, with fast neutrino reactions and superconducting cores, ambipolar diffusion results in advection velocities of several km/kyr. This velocity can substantially reorganize magnetic fields in magnetar cores, in a way that can only be confirmed by dynamical simulations.

We study the torsional magneto-elastic oscillations of relativistic superfluid magnetars and explore the effects of a phase transition in the crust-core interface (nuclear pasta) which results in a weaker elastic response. Exploring various models with different extension of nuclear pasta phases, we find that the differences in the oscillation spectrum present in purely elastic modes (weak magnetic field), are smeared out with increasing strength of the magnetic field. For magnetar conditions, the main characteristic and features of models without nuclear pasta are preserved. We find in general two classes of magneto-elastic oscillations which exhibit a different oscillation pattern. For B_p < 4 x 10^{14} G, the spectrum is characterised by the turning points and edges of the continuum which are confined into the star's core, and have no constant phase. Increasing the magnetic field, we find, in addition, several magneto-elastic oscillations which reach the surface and have an angular structure similar to crustal modes. These global magneto-elastic oscillations show a constant phase and become dominant when B_p > 5 x 10^{14} G. We do not find any evidence of fundamental pure crustal modes in the low frequency range (below 200 Hz) for B_p >= 10^{14} G.

We consider the local dynamics of a realistic neutron star core, including composition gradients, superfluidity and thermal effects. The main focus is on the gravity g-modes, which are supported by composition stratification and thermal gradients. We derive the equations that govern this problem in full detail, paying particular attention to the input that needs to be provided through the equation of state and distinguishing between normal and superfluid regions. The analysis highlights a number of key issues that should be kept in mind whenever equation of state data is compiled from nuclear physics for use in neutron star calculations. We provide explicit results for a particular stellar model and a specific nucleonic equation of state, making use of cooling simulations to show how the local wave spectrum evolves as the star ages. Our results show that the composition gradient is effectively dominated by the muons whenever they are present. When the star cools below the superfluid transition, the support for g-modes at lower densities (where there are no muons) is entirely thermal. We confirm the recent suggestion that the g-modes in this region may be unstable, but our results indicate that this instability will be weak and would only be present for a brief period of the star's life. Our analysis accounts for the presence of thermal excitations encoded in entrainment between the entropy and the superfluid component. Finally, we discuss the complete spectrum, including the normal sound waves and, in superfluid regions, the second sound.

The motion of a test particle in the gravitational field of a non-spherical source endowed with both mass and mass quadrupole moment is investigated when a test radiation field is also present. The background is described by the Erez-Rosen solution, which is a static spacetime belonging to the Weyl class of solutions to the vacuum Einstein's field equations, and reduces to the familiar Schwarzschild solution when the quadrupole parameter vanishes. The radiation flux has a fixed but arbitrary (non-zero) angular momentum. The interaction with the radiation field is assumed to be Thomson-like, i.e., the particles absorb and re-emit radiation, thus suffering for a friction-like drag force. Such an additional force is responsible for the Poynting-Robertson effect, which is well established in the framework of Newtonian gravity and has been recently extended to the general theory of relativity. The balance between gravitational attraction, centrifugal force and radiation drag leads to the occurrence of equilibrium circular orbits which are attractors for the surrounding matter for every fixed value of the interaction strength. The presence of the quadrupolar structure of the source introduces a further degree of freedom: there exists a whole family of equilibrium orbits parametrized by the quadrupole parameter, generalizing previous works. This scenario is expected to play a role in the context of accretion matter around compact objects.

We study the low T/W instability associated with the f-mode of differentially rotating stars. Our stellar models are described by a polytropic equation of state and the rotation profile is given by the standard j-constant law. The properties of the relevant oscillation modes, including the instability growth time, are determined from time evolutions of the linearised dynamical equations in Newtonian gravity. In order to analyse the instability we monitor also the canonical energy and angular momentum. Our results demonstrate that the l=m=2 f-mode becomes unstable as soon as a co-rotation point develops inside the star (i.e. whenever there is a point where the mode's pattern speed matches the bulk angular velocity). Considering various degrees of differential rotation, we show that the instability grows faster deep inside the co-rotation region and deduce an empirical relation that correlates the mode frequency and the star's parameters, which captures the main features of the l=m=2 f-mode growth time. This function is proportional to the product of the kinetic to gravitational energy ratio and the gradient of the star's spin, strengthening further the relationship between the co-rotation point and the low T/W instability. We briefly consider also the l=m=2 r-mode and demonstrate that it never moves far inside the co-rotation region even for significant differential rotation.

We study the time-evolution of axisymmetric oscillations of superfluid magnetars with a poloidal magnetic field and an elastic crust, working in Newtonian gravity. Extending earlier models, we study the effects of composition gradients and entrainment on the magneto-elastic wave spectrum and on the potential identification of the observed quasi periodic oscillations (QPOs). We use two-fluid polytropic equations of state to construct our stellar models, which mimic realistic composition gradient configurations. The basic features of the axial axisymmetric spectrum of normal fluid stars are reproduced by our results and in addition we find several magneto-elastic waves with a mixed character. In the core, these oscillations mimic the shear mode pattern of the crust as a result of the strong dynamical coupling between these two regions. Incorporating the most recent entrainment configurations in our models, we find that they have a double effect on the spectrum: the magnetic oscillations of the core have a frequency enhancement, while the mixed magneto-elastic waves originating in the crust are moved towards the frequencies of the single-fluid case. The distribution of lower-frequency magneto-elastic oscillations for our models is qualitatively similar to the observed magnetar QPOs with $\nu < 155 $Hz. \emph{In particular, some of these QPOs could represent mixed magneto-elastic oscillations with frequencies not greatly different from the crustal modes of an unmagnetised star.} We find that many QPOs could even be accounted for using a model with a relatively weak polar field of $B_{p} \simeq 3\cdot 10^{14}$G, because of the superfluid enhancement of magnetic oscillations. Finally, we discuss the possible identification of 625 and 1837~Hz QPOs either with non-axisymmetric modes or with high-frequency axisymmetric QPOs excited by crustal mode overtones.

We study the dynamical evolution of the gravitational-wave driven instability of the f-mode in rapidly rotating relativistic stars. With an approach based on linear perturbation theory we describe the evolution of the mode amplitude and follow the trajectory of a newborn neutron star through its instability window. The influence on the f-mode instability of the magnetic field and the presence of an unstable r-mode is also considered. Two different configurations are studied in more detail; an N = 1 polytrope with a typical mass and radius and a more massive polytropic N = 0.62 model with gravitational mass M = 1.98 Msun. We study several evolutions with different initial rotation rates and temperature and determine the gravitational waves radiated during the instability. In more massive models, an unstable f-mode with a saturation energy of about 1e-6 Msun c^2 may generate a gravitational-wave signal which can be detected by the Advanced LIGO/Virgo detector from the Virgo cluster. The magnetic field affects the evolution and then the detectability of the gravitational radiation when its strength is higher than 1e12 G, while the effects of an unstable r-mode become dominant when this mode reaches the maximum saturation value allowed by non-linear mode couplings. However, the relative saturation amplitude of the f- and r-modes must be known more accurately in order to provide a definitive answer to this issue. From the thermal evolution we find also that the heat generated by shear viscosity during the saturation phase completely balances the neutrinos' cooling and prevents the star from entering the regime of mutual friction. The evolution time of the instability is therefore longer and the star loses significantly larger amounts of angular momentum via gravitational waves.

The violent giant flares of magnetars excite QPOs which persist for hundreds of seconds, as seen in the X-ray tail following the initial burst. Recent studies, based on single-fluid barotropic magnetar models, have suggested that the lower-frequency QPOs correspond to magneto-elastic oscillations of the star. The higher frequencies, however, in particular the strong 625 Hz peak, have proved harder to explain, except as high mode multipoles. In this work we study the time evolutions of non-axisymmetric oscillations of two-fluid Newtonian magnetars with no crust. We consider models with superfluid neutrons and normal protons, and poloidal and toroidal background field configurations. We show that multi-fluid physics (composition-gradient stratification, entrainment) tends to increase Alfven mode frequencies significantly from their values in a single-fluid barotropic model. The higher-frequency magnetar QPOs may then be naturally interpreted as Alfven oscillations of the multi-fluid stellar core. The lower-frequency QPOs are less easily explained within our purely fluid core model, but we discuss the possibility that these are crustal modes.

We study the damping of the gravitational radiation-driven f-mode instability in rotating neutron stars by nonlinear bulk viscosity in the so-called supra-thermal regime. In this regime the dissipative action of bulk viscosity is known to be enhanced as a result of nonlinear contributions with respect to the oscillation amplitude. Our analysis of the f-mode instability is based on a time-domain code that evolves linear perturbations of rapidly rotating polytropic neutron star models. The extracted mode frequency and eigenfunctions are subsequently used in standard energy integrals for the gravitational wave growth and viscous damping. We find that nonlinear bulk viscosity has a moderate impact on the size of the f-mode instability window, becoming an important factor and saturating the mode's growth at a relatively large oscillation amplitude. We show similarly that nonlinear bulk viscosity leads to a rather high saturation amplitude even for the r-mode instability. In addition, we show that the action of bulk viscosity can be significantly mitigated by the presence of superfluidity in neutron star matter. Apart from revising the f-mode's instability window we provide results on the mode's gravitational wave detectability. Considering an f-mode-unstable neutron star located in the Virgo cluster and assuming a mode amplitude at the level allowed by bulk viscosity, we find that the emitted gravitational wave signal could be detectable by advanced ground-based detectors such as Advanced LIGO/Virgo and the Einstein Telescope.

We study the effects of an elastic crust on the oscillation spectrum of superfluid neutron stars. Within the two fluid formalism, we consider Newtonian stellar models that include the relevant constituents of a mature neutron stars. The core is formed by a mixture of superfluid neutrons and a conglomerate of charged particles, while the inner crust is described by a lattice of nuclei permeated by superfluid neutrons. We linearise the Poisson and the conservation equations of nonrotating superfluid stars and study the effects of elasticity, entrainment and composition stratification on the shear and acoustic modes. In both the core and the crust, the entrainment is derived from recent results for the nucleon effective mass. Solving the perturbation equations as an eigenvalue problem, we find that the presence of superfluid neutrons in the crust and their large effective mass may have significant impact on the star's oscillation spectrum.

In this paper we consider a simple two-fluid model for pulsar glitches. We derive the basic equations that govern the spin evolution of the system from two-fluid hydrodynamics, accounting for the vortex mediated mutual friction force that determines the glitch rise. This leads to a simple "bulk" model that can be used to describe the main properties of a glitch event resulting from vortex unpinning. In order to model the long term relaxation following the glitch our model would require additional assumptions regarding the repinning of vortices, an issue that we only touch upon briefly. Instead, we focus on comparing the phenomenological model to results obtained from time-evolutions of the linearised two-fluid equations, i.e. a "hydrodynamic" model for glitches. This allows us to study, for the first time, dynamics that was "averaged" in the bulk model, i.e. consider the various neutron star oscillation modes that are excited during a glitch. The hydro-results are of some relevance for efforts to detect gravitational waves from glitching pulsars, although the conclusions drawn from our rather simple model are pessimistic as far as the detectability of these events is concerned.

We study time evolutions of superfluid neutron stars, focussing on the nature of the oscillation spectrum, the effect of mutual friction force on the oscillations and the hydrodynamical spin-up phase of pulsar glitches. We linearise the dynamical equations of a Newtonian two-fluid model for rapidly rotating backgrounds. In the axisymmetric equilibrium configurations, the two fluid components corotate and are in beta-equilibrium. We use analytical equations of state that generate stratified and non-stratified stellar models, which enable us to study the coupling between the dynamical degrees of freedom of the system. By means of time evolutions of the linearised dynamical equations, we determine the spectrum of axisymmetric and non-axisymmetric oscillation modes, accounting for the contribution of the gravitational potential perturbations, i.e. without adopting the Cowling approximation. We study the mutual friction damping of the superfluid oscillations and consider the effects of the non-dissipative part of the mutual friction force on the mode frequencies. We also provide technical details and relevant tests for the hydrodynamical model of pulsar glitches discussed by Sidery, Passamonti and Andersson (2010). In particular, we describe the method used to generate the initial data that mimic the pre-glitch state, and derive the equations that are used to extract the gravitational-wave signal.

We investigate the oscillation spectrum of rotating Newtonian neutron stars endowed with purely toroidal magnetic fields, using a time evolution code to evolve linear perturbations in the Cowling approximation. The background star is generated by numerically solving the MHD equilibrium equations and may be nonspherical by virtue of both rotation and magnetic effects; hence our perturbations and background are fully consistent. Whilst the background field is purely toroidal, the perturbed field is mixed poloidal-toroidal. From Fourier analysis of the perturbations we are able to identify a number of magnetically-restored Alfv\'en (or $a$-) modes. We show that in a rotating star pure inertial and $a$-modes are replaced by hybrid magneto-inertial modes, which reduce to $a$-modes in the nonrotating limit and inertial modes in the nonmagnetic limit. We show that the $r$-mode instability is suppressed by magnetic fields in sufficiently slowly rotating stars. In addition, we determine magnetic frequency shifts in the $f$-mode. We discuss the astrophysical relevance of our results, in particular for magnetar oscillations.

Using time evolutions of the relevant linearised equations we study non-axisymmetric oscillations of rapidly rotating and superfluid neutron stars. We consider perturbations of Newtonian axisymmetric background configurations and account for the presence of superfluid components via the standard two-fluid model. Within the Cowling approximation, we are able to carry out evolutions for uniformly rotating stars up to the mass-shedding limit. This leads to the first detailed analysis of superfluid neutron star oscillations in the fast rotation regime, where the star is significantly deformed by the centrifugal force. For simplicity, we focus on background models where the two fluids (superfluid neutrons and protons) co-rotate, are in beta-equilibrium and coexist throughout the volume of the star. We construct sequences of rotating stars for two analytical model equations of state. These models represent relatively simple generalisations of single fluid, polytropic stars. We study the effects of entrainment, rotation and symmetry energy on non-radial oscillations of these models. Our results show that entrainment and symmetry energy can have a significant effect on the rotational splitting of non-axisymmetric modes. In particular, the symmetry energy modifies the inertial mode frequencies considerably in the regime of fast rotation.

We use time-evolutions of the linear perturbation equations to study the oscillations of rapidly rotating neutrons stars. Our models account for the buoyancy due to composition gradients and we study, for the first time, the nature of the resultant g-modes in a fast spinning star. We provide detailed comparisons of non-stratified and stratified models. This leads to an improved understanding of the relationship between the inertial modes of a non-stratified star and the g-modes of a stratified system. In particular, we demonstrate that each g-mode becomes rotation-dominated, i.e. approaches a particular inertial mode, as the rotation rate of the star is increased. We also discuss issues relating to the gravitational-wave driven instability of the various classes of oscillation modes.

We develop a new perturbative framework for studying the r-modes of rotating superfluid neutron stars. Our analysis accounts for the centrifugal deformation of the star, and considers the two-fluid dynamics at linear order in the perturbed velocities. Our main focus is on a simple model system where the total density profile is that of an $n=1$ polytrope. We derive a partially analytic solution for the superfluid analogue of the classical r-mode. This solution is used to analyse the relevance of the vortex mediated mutual friction damping, confirming that this dissipation mechanism is unlikely to suppress the gravitational-wave driven instability in rapidly spinning superfluid neutron stars. Our calculation of the superfluid r-modes is significantly simpler than previous approaches, because it decouples the r-mode from all other inertial modes of the system. This leads to the results being clearer, but it also means that we cannot comment on the relevance of potential avoided crossings (and associated "resonances") that may occur for particular parameter values. Our analysis of the mutual friction damping differs from previous studies in two important ways. Firstly, we incorporate realistic pairing gaps which means that the regions of superfluidity in the star's core vary with temperature. Secondly, we allow the mutual friction parameters to take the whole range of permissible values rather than focussing on a particular mechanism. Thus, we consider not only the weak drag regime, but also the strong drag regime where the fluid dynamics is significantly different.

Non-axisymmetric oscillations of differentially rotating stars are studied using both slow rotation and Cowling approximation. The equilibrium stellar models are relativistic polytropes where differential rotation is described by the relativistic j-constant rotation law. The oscillation spectrum is studied versus three main parameters: the stellar compactness $M/R$, the degree of differential rotation $A$ and the number of maximun couplings $\ell_{\rm max}$. It is shown that the rotational splitting of the non-axisymmetric modes is strongly enhached by increasing the compactness of the star and the degree of differential rotation. Finally, we investigate the relation between the fundamental quadrupole mode and the corotation band of differentially rotating stars.

The post-bounce oscillations of newly-born relativistic stars are expected to lead to gravitational-wave emission through the excitation of nonradial oscillation modes. At the same time, the star is oscillating in its radial modes, with a central density variation that can reach several percent. Nonlinear couplings between radial oscillations and polar nonradial modes lead to the appearance of combination frequencies (sums and differences of the linear mode frequencies). We study such combination frequencies using a gauge-invariant perturbative formalism, which includes bilinear coupling terms between different oscillation modes. For typical values of the energy stored in each mode we find that gravitational waves emitted at combination frequencies could become detectable in galactic core-collapse supernovae with advanced interferometric or wide-band resonant detectors.

The equations describing nonradial adiabatic oscillations of differentially rotating relativistic stars are derived in relativistic slow rotation approximation. The differentially rotating configuration is described by a perturbative version of the relativistic j-constant rotation law. Focusing on the oscillation properties of the stellar fluid, the adiabatic nonradial perturbations are studied in the Cowling approximation with a system of five partial differential equations. In these equations, differential rotation introduces new coupling terms between the perturbative quantites with respect to the uniformly rotating stars. In particular, we investigate the axisymmetric and barotropic oscillations and compare their spectral properties with those obtained in nonlinear hydrodynamical studies. The perturbative description of the differentially rotating background and the oscillation spectrum agree within a few percent with those of the nonlinear studies.

This thesis investigates in the time domain a particular class of second order perturbations of a perfect fluid non-rotating compact star: those arising from the coupling between first order radial and non-radial perturbations. This problem has been treated by developing a gauge invariant formalism based on the 2-parameter perturbation theory (Sopuerta, Bruni and Gualtieri, 2004) where the radial and non-radial perturbations have been separately parameterized. The non-linear perturbations obey inhomogeneous partial differential equations, where the structure of the differential operator is given by the previous perturbative orders and the source terms are quadratic in the first order perturbations. In the exterior spacetime the sources vanish, thus the gravitational wave properties are completely described by the second order Zerilli and Regge-Wheeler functions. As main initial configuration we have considered a first order differentially rotating and radially pulsating star. Although at first perturbative order this configuration does not exhibit any gravitational radiation, we have found a new interesting gravitational signal at non-linear order, in which the radial normal modes are precisely mirrored. In addition, a resonance effect is present when the frequencies of the radial pulsations are close to the first axial w-mode. Finally, we have roughly estimated the damping times of the radial pulsations due to the non-linear gravitational emission. The coupling near the resonance results to be a very effective mechanism for extracting energy from the radial oscillations.

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.

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.