Search

Neutron stars harbour extremely powerful magnetic fields, leading to their shape being deformed. Their magnetic deformation depends both on the geometry - and strength - of their internal magnetic field and on their composition, encoded by the equation of state. However, both the details of the internal magnetic structure and the equation of state of the innermost part of neutron stars are mostly unkown. We performed a study of numerical models of magnetised, static, axisymmetric neutron stars in general relativity and in one of its most promising extensions, scalar-tensor theories. We did so by using several realistic equations of state currently allowed by observational and nuclear physics constraints, considering also those for strange quark stars. We show that it is possible to find simple relations among the magnetic deformation of a neutron star, its Komar mass, and its circumferential radius in the case of purely poloidal and purely toroidal magnetic configurations satisfying the equilibrium criterion in the Bernoulli formalism. These relations are quasi-universal, in the sense that they mostly do not depend on the equation of state. Our results, being formulated in terms of potentially observable quantities, could help to understand the magnetic properties of neutron stars interiors and the detectability of continuous gravitational waves by isolated neutron stars, independently of their equation of state. In the case of scalar-tensor theories, these relations depend also on the scalar charge of the neutron stars, thus potentially providing a new way to set constraints on the theory of gravitation., Comment: 19 pages, 5 figures, submitted as proceeding to the Sixteenth Marcel Grossmann Meeting, this article draws heavily from arXiv:2106.00603

Neutron stars are known to host extremely powerful magnetic fields. Among other effects, one of the consequences of harbouring such fields is the deformation of the neutron star structure, leading, together with rotation, to the emission of continuous gravitational waves. On the one hand, the details of their internal magnetic fields are mostly unknown. Likewise, their internal structure, encoded by the equation of state, is highly uncertain. Here we present a study of axisymmetric models of isolated magnetised neutron stars, for various realistic equations of state considered viable by observations and nuclear physics constraints. We show that it is possible to find simple relations between the magnetic deformation of a neutron star, its Komar mass and its circumferential radius. Such relations are quasi-universal, meaning that they are mostly independent on the equation of state of the neutron star and only slightly dependent on the magnetic field configuration. Being formulated in terms of potentially observable quantities, as we discuss, our results could help to constrain the magnetic properties of the neutron star interior and to better assess the detectability of continuous gravitational waves by isolated neutron stars, without knowing their equation of state. Our results are derived both in general relativity and in scalar-tensor theories - one of the most promising extensions of general relativity - in this case by considering also the scalar charge. We show that even in this case general relations hold that account for deviations from general relativity, that could potentially be used to set constraints on the gravitational theory., Comment: Accepted for publication in A&A; minor language corrections

Scalar-tensor theories are among the most promising alternatives to general relativity that have been developed to account for some long standing issues in our understanding of gravity. Some of these theories predict the existence of a non-linear phenomenon, spontaneous scalarisation, which can lead to the appearance of sizeable modifications to general relativity in the presence of compact matter distributions, namely neutron stars. On the one hand, one of the effects of the scalar field is to modify the emission of gravitational waves, both due to variations in the quadrupolar deformation of the star and to the presence of additional modes of emission. On the other hand, neutron stars are known to harbour extremely powerful magnetic fields which can affect their structure and shape, leading in turn to the emission of gravitational waves, this time due to a magnetic quadrupolar deformation. We investigate here how the presence of spontaneous scalarisation can affect the magnetic deformation of neutron stars and their emission of quadrupolar gravitational waves, both of tensor and scalar nature. We will show that it is possible to provide simple parameterisations of the magnetic deformation and gravitational wave power of neutron stars in terms of their baryonic mass, circumferential radius and scalar charge, and that a universal scaling exists independently of the magnetic field geometry and of the parameters of the scalar-tensor theory. Finally, we comment on the observability of the deviations of the gravitational waves strain from general relativity by current and future observatories., Comment: 9 pages, 2 figures, accepted for publication in A&A

Among the possible extensions of General Relativity that have been put forward in order to address some long standing issues in our understanding of the Universe, Scalar-Tensor Theories have received a lot of attention for their simplicity. Interestingly, some of these predict a potentially observable non-linear phenomenon, known as \textit{spontaneous scalarisation}, in the presence of highly compact matter distributions, like the case of neutron stars. Neutron stars are ideal laboratories to investigate the properties of matter under extreme conditions, and in particular they are known to harbour the strongest magnetic fields in the Universe. Here, for the first time, we present a detailed study of magnetised neutron stars in Scalar-Tensor Theories. First, we show that the formalism developed for the study of magnetised neutron stars in General Relativity, based on the \textit{eXtended Conformally Flat Condition}, can easily be extended in the presence of a non-minimally coupled scalar field, retaining many of its numerical advantages. We then carry out a study of the parameter space considering the two extreme geometries of purely toroidal and purely poloidal magnetic fields, varying both the strength of the magnetic field and the intensity of scalarisation. We compare our results with magnetised general-relativistic solutions and un-magnetised scalarised solutions, showing how the mutual interplay between magnetic and scalar fields affect the magnetic and the scalarisation properties of neutron stars. In particular, we focus our discussion on magnetic deformability, maximum mass and range of scalarisation., Comment: accepted for publication by A&A; minor language corrections; minor typos correction

The extraordinary energetic activity of magnetars is usually explained in terms of dissipation of a huge internal magnetic field of the order of $10^{15-16}$G. How such a strong magnetic field can originate during the formation of a neutron star is still subject of active research. An important role can be played by fast rotation: if magnetars are born as millisecond rotators dynamo mechanisms may efficiently amplify the magnetic field inherited from the progenitor star during the collapse. In this case, the combination of rapid rotation and strong magnetic field determine the right physical condition not only for the development of a powerful jet driven explosion, manifesting as a gamma ray burst, but also for a copious gravitational waves emission. Strong magnetic fields are indeed able to induce substantial quadrupolar deformations in the star. In this paper we analyze the joint effect of rotation and magnetization on the structure of a polytropic and axisymmetric neutron star, within the ideal magneto-hydrodynamic regime. We will consider either purely toroidal or purely poloidal magnetic field geometries. Through the sampling of a large parameter space, we generalize previous results in literature, inferring new quantitative relations that allow for a parametrization of the induced deformation, that takes into account also the effects due to the stellar compactness and the current distribution. Finally, in the case of purely poloidal field, we also discuss how different prescriptions on the surface charge distribution (a gauge freedom) modify the properties of the surrounding electrosphere and its physical implications., Comment: 25 pages, 17 figures, 6 tables, accepted for publication in MNRAS

Magnetic fields play a critical role in the phenomenology of neutron stars. There is virtually no observable aspect which is not governed by them. Despite this, only recently efforts have been done to model magnetic fields in the correct general relativistic regime, characteristic of these compact objects. In this work we present, for the first time a comprehensive and detailed parameter study, in general relativity, of the role that the current distribution, and the related magnetic field structure, have in determining the precise structure of neutron stars. In particular, we show how the presence of localized currents can modify the field strength at the stellar surface, and we look for general trends, both in terms of energetic properties, and magnetic field configurations. Here we verify that, among other things, for a large class of different current distributions the resulting magnetic configurations are always dominated by the poloidal component of the current., Comment: 14 pages, 13 figures, accepted for publication in MNRAS

Soft Gamma-Ray Repeaters and Anomalous X-Ray Pulsars are extreme manifestations of the most magnetized neutron stars: magnetars. The phenomenology of their emission and spectral properties strongly support the idea that the magnetospheres of these astrophysical objects are tightly twisted in the vicinity of the star. Previous studies on equilibrium configurations have so far focused on either the internal or the external magnetic field configuration, without considering a real coupling between the two fields. Here we investigate numerical equilibrium models of magnetized neutron stars endowed with a confined twisted magnetosphere, solving the general relativistic Grad-Shafranov equation both in the interior and in the exterior of the compact object. A comprehensive study of the parameters space is provided to investigate the effects of different current distributions on the overall magnetic field structure., Comment: 15 pages, 14 figures. Accepted for publication in MNRAS

High-energy phenomena in astrophysics involve quite generally a combination of relativistic motions and strong gravity. The simultaneous solution of Einstein equations and General Relativistic MHD equations is thus necessary to model with accuracy such phenomena. The so-called Conformally Flat Condition (CFC) allows a simplified treatment of Einstein equations, that can be particularly efficient in those contexts where gravitational wave emission is negligible, like core-collapse, or the formation/evolution of neutron stars. We have developed a set of codes to model axisymmetric MHD flows, in General Relativity, where the solution of Einstein equations is achieved with a semi-spectral scheme. Here, we will show how this framework is particularly well suited to investigate neutron star equilibrium models in the presence of strong magnetic fields and we will present the XNS code, that has been recently developed and here updated to treat poloidal and mixed configurations., Comment: 6 pages, 2 figures, Proceedings of the 8th International Conference on Numerical Modeling of Space Plasma Flows, 1-5 July 2013, Biarritz, France

Magnetic fields play a crucial role in many astrophysical scenarios and, in particular, are of paramount importance in the emission mechanism and evolution of Neutron Stars (NSs). To understand the role of the magnetic field in compact objects it is important to obtain, as a first step, accurate equilibrium models for magnetized NSs. Using the conformally flat approximation we solve the Einstein's equations together with the GRMHD equations in the case of a static axisymmetryc NS taking into account different types of magnetic configuration. This allows us to investigate the effect of the magnetic field on global properties of NSs such as their deformation., Comment: 4 pages, 6 figures, Proceedings of the 4th International Conference High Energy Phenomena in Relativistic Outflows "HEPRO IV", 23-26 July 2013, Heidelberg, Germany. To be publised in International Journal of Modern Physics: Conference Series

We present an extension of relativistic single-particle distribution function for weakly interacting particles at local thermodynamical equilibrium including spin degrees of freedom, for massive spin 1/2 particles. We infer, on the basis of the global equilibrium case, that at local thermodynamical equilibrium particles acquire a net polarization proportional to the vorticity of the inverse temperature four-vector field. The obtained formula for polarization also implies that a steady gradient of temperature entails a polarization orthogonal to particle momentum. The single-particle distribution function in momentum space extends the so-called Cooper-Frye formula to particles with spin 1/2 and allows to predict their polarization in relativistic heavy ion collisions at the freeze-out., Comment: 15 pages. Final version to appear in Annals of Physics

The powerful high-energy phenomena typically encountered in astrophysics invariably involve physical engines, like neutron stars and black hole accretion disks, characterized by a combination of highly magnetized plasmas, strong gravitational fields, and relativistic motions. In recent years numerical schemes for General Relativistic MHD (GRMHD) have been developed to model the multidimensional dynamics of such systems, including the possibility of an evolving spacetime. Such schemes have been also extended beyond the ideal limit including the effects of resistivity, in an attempt to model dissipative physical processes acting on small scales (sub-grid effects) over the global dynamics. Along the same lines, magnetic fields could be amplified by the presence of turbulent dynamo processes, as often invoked to explain the high values of magnetization required in accretion disks and neutron stars. Here we present, for the first time, a further extension to include the possibility of a mean-field dynamo action within the framework of numerical 3+1 (resistive) GRMHD. A fully covariant dynamo closure is proposed, in analogy with the classical theory, assuming a simple alpha-effect in the comoving frame. Its implementation into a finite-difference scheme for GRMHD in dynamical spacetimes [the X-ECHO code: (Bucciantini and Del Zanna 2011)] is described, and a set of numerical test is presented and compared with analytical solutions wherever possible., Comment: 16 pages, 11 figures, accepted for publication in MNRAS

We present a new numerical code, ECHO, based on an Eulerian Conservative High Order scheme for time dependent three-dimensional general relativistic magnetohydrodynamics (GRMHD) and magnetodynamics (GRMD). ECHO is aimed at providing a shock-capturing conservative method able to work at an arbitrary level of formal accuracy (for smooth flows), where the other existing GRMHD and GRMD schemes yield an overall second order at most. Moreover, our goal is to present a general framework, based on the 3+1 Eulerian formalism, allowing for different sets of equations, different algorithms, and working in a generic space-time metric, so that ECHO may be easily coupled to any solver for Einstein's equations. Various high order reconstruction methods are implemented and a two-wave approximate Riemann solver is used. The induction equation is treated by adopting the Upwind Constrained Transport (UCT) procedures, appropriate to preserve the divergence-free condition of the magnetic field in shock-capturing methods. The limiting case of magnetodynamics (also known as force-free degenerate electrodynamics) is implemented by simply replacing the fluid velocity with the electromagnetic drift velocity and by neglecting the matter contribution to the stress tensor. ECHO is particularly accurate, efficient, versatile, and robust. It has been tested against several astrophysical applications, including a novel test on the propagation of large amplitude circularly polarized Alfven waves. In particular, we show that reconstruction based on a Monotonicity Preserving filter applied to a fixed 5-point stencil gives highly accurate results for smooth solutions, both in flat and curved metric (up to the nominal fifth order), while at the same time providing sharp profiles in tests involving discontinuities., Comment: 20 pages, revised version submitted to A&A