No Arabic abstract
Previously we found that large amplitude $r$-modes could decay catastrophically due to nonlinear hydrodynamic effects. In this paper we found the particular coupling mechanism responsible for this catastrophic decay, and identified the fluid modes involved. We find that for a neutron star described by a polytropic equation of state with polytropic index $Gamma=2$, the coupling strength of the particular three-mode interaction causing the decay is strong enough that the usual picture of the $r$-mode instability with a flow pattern dominated by that of an $r$-mode can only be valid for the dimensionless $r$-mode amplitude less than $O(10^{-2})$.
We study the dynamical evolution of a large amplitude r-mode by numerical simulations. R-modes in neutron stars are unstable growing modes, driven by gravitational radiation reaction. In these simulations, r-modes of amplitude unity or above are destroyed by a catastrophic decay: A large amplitude r-mode gradually leaks energy into other fluid modes, which in turn act nonlinearly with the r-mode, leading to the onset of the rapid decay. As a result the r-mode suddenly breaks down into a differentially rotating configuration. The catastrophic decay does not appear to be related to shock waves at the stars surface. The limit it imposes on the r-mode amplitude is significantly smaller than that suggested by previous fully nonlinear numerical simulations.
In $f(R)$ gravity and Brans-Dicke theory with scalar potentials, we study the structure of neutron stars on a spherically symmetric and static background for two equations of state: SLy and FPS. In massless BD theory, the presence of a scalar coupling $Q$ with matter works to change the star radius in comparison to General Relativity, while the maximum allowed mass of neutron stars is hardly modified for both SLy and FPS equations of state. In Brans-Dicke theory with the massive potential $V(phi)=m^2 phi^2/2$, where $m^2$ is a positive constant, we show the difficulty of realizing neutron star solutions with a stable field profile due to the existence of an exponentially growing mode outside the star. As in $f(R)$ gravity with the $R^2$ term, this property is related to the requirement of extra boundary conditions of the field at the surface of star. For the self-coupling potential $V(phi)=lambda phi^4/4$, this problem can be circumvented by the fact that the second derivative $V_{,phi phi}=3lambdaphi^2$ approaches 0 at spatial infinity. In this case, we numerically show the existence of neutron star solutions for both SLy and FPS equations of state and discuss how the mass-radius relation is modified as compared to General Relativity.
We model the nonlinear saturation of the r-mode instability via three-mode couplings and the effects of the instability on the spin evolution of young neutron stars. We include one mode triplet consisting of the r-mode and two near resonant inertial modes that couple to it. We find that the spectrum of evolutions is more diverse than previously thought. The evolution of the star is dynamic and initially dominated by fast neutrino cooling. Nonlinear effects become important when the r-mode amplitude grows above its first parametric instability threshold. The balance between neutrino cooling and viscous heating plays an important role in the evolution. Depending on the initial r-mode amplitude, and on the strength of the viscosity and of the cooling this balance can occur at different temperatures. If thermal equilibrium occurs on the r-mode stability curve, where gravitational driving equals viscous damping, the evolution may be adequately described by a one-mode model. Otherwise, nonlinear effects are important and lead to various more complicated scenarios. Once thermal balance occurs, the star spins-down oscillating between thermal equilibrium states until the instability is no longer active. For lower viscosity we observe runaway behavior in which the r-mode amplitude passes several parametric instability thresholds. In this case more modes need to be included to model the evolution accurately. In the most optimistic case, we find that gravitational radiation from the r-mode instability in a very young, fast spinning neutron star within about 1 Mpc of Earth may be detectable by advanced LIGO for years, and perhaps decades, after formation. Details regarding the amplitude and duration of the emission depend on the internal dissipation of the modes of the star, which would be probed by such detections.
In the context of f(R)=R + alpha R^2 gravity, we study the existence of neutron and quark stars with no intermediate approximations in the generalised system of Tolman-Oppenheimer-Volkov equations. Analysis shows that for positive alphas the scalar curvature does not drop to zero at the star surface (as in General Relativity) but exponentially decreases with distance. Also the stellar mass bounded by star surface decreases when the value alpha increases. Nonetheless distant observers would observe a gravitational mass due to appearance of a so-called gravitational sphere around the star. The non-zero curvature contribution to the gravitational mass eventually is shown to compensate the stellar mass decrease for growing alphas. We perform our analysis for several equations of state including purely hadronic configurations as well as hyperons and quark stars. In all cases, we assess that the relation between the parameter $alpha$ and the gravitational mass weakly depend upon the chosen equation of state. Another interesting feature is the increase of the star radius in comparison to General Relativity for stars with masses close to maximal, whereas for intermediate masses around 1.4-1.6 solar masses, the radius of star depends upon alpha very weakly. Also the decrease in the mass bounded by star surface may cause the surface redshift to decrease in R^2-gravity when compared to Einsteinian predictions. This effect is shown to hardly depend upon the observed gravitational mass. Finally, for negative values of alpha our analysis shows that outside the star the scalar curvature has damped oscillations but the contribution of the gravitational sphere into the gravitational mass increases indefinitely with radial distance putting into question the very existence of such relativistic stars.
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%.