Search

Relativistic magnetic reconnection is one of the most fundamental mechanisms considered responsible for the acceleration of relativistic particles in astrophysical jets and magnetospheres of compact objects. Understanding the properties of the dissipation of magnetic fields and the formation of non-ideal electric fields is of paramount importance to quantify the efficiency of reconnection at energizing charged particles. Recent results from particle-in-cell (PIC) simulations suggest that the fundamental properties of how magnetic fields dissipate in a current sheet might be captured by an "effective resistivity" formulation, which would locally enhance the amount of magnetic energy dissipated and favor the onset of fast reconnection. Our goal is to assess this ansatz quantitatively by comparing fluid models of magnetic reconnection with a non-constant magnetic diffusivity and fully-kinetic models. We perform 2D resistive relativistic magnetohydrodynamic (ResRMHD) simulations of magnetic reconnection combined to PIC simulations using the same initial conditions (namely a Harris current sheet). We explore the impact of crucial parameters such as the plasma magnetization, its mass density, the grid resolution, and the characteristic plasma skin depth. Our ResRMHD models with effective resistivity can quantitatively reproduce the dynamics of fully-kinetic models of relativistic magnetic reconnection. In particular, they lead to reconnection rates consistent with PIC simulations, while for constant-resistivity fluid models the reconnection dynamics is generally 10 times slower. Even at modest resolutions the adoption of an effective resistivity can qualitatively capture the properties of kinetic reconnection models and produce reconnection rates compatible with collisionless models, i.e. of the order of $\sim10^{-1}$., Comment: 17 pages, 17 figures; accepted for publication on A&A

Space plasmas are three-dimensional dynamic entities. Except under very special circumstances, their structure in space and their behavior in time are not related in any simple way. Therefore, single spacecraft in situ measurements cannot unambiguously unravel the full space-time structure of the heliospheric plasmas of interest in the inner heliosphere, in the Geospace environment, or the outer heliosphere. This shortcoming leaves numerous central questions incompletely answered. Deficiencies remain in at least two important subjects, Space Weather and fundamental plasma turbulence theory, due to a lack of a more complete understanding of the space-time structure of dynamic plasmas. Only with multispacecraft measurements over suitable spans of spatial separation and temporal duration can these ambiguities be resolved. We note that these characterizations apply to turbulence across a wide range of scales, and also equally well to shocks, flux ropes, magnetic clouds, current sheets, stream interactions, etc. In the following, we will describe the basic requirements for resolving space-time structure in general, using turbulence' as both an example and a principal target or study. Several types of missions are suggested to resolve space-time structure throughout the Heliosphere., Comment: White Paper submitted to: Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033. arXiv admin note: substantial text overlap with arXiv:1903.06890

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

Neutron stars, and magnetars in particular, are known to host the strongest magnetic fields in the Universe. The origin of these strong fields is a matter of controversy. In this preliminary work, via numerical simulations, we study, for the first time in non-ideal general relativistic magnetohydrodynamic (GRMHD) regime, the growth of the magnetic field due to the action of the mean-field dynamo due to sub-scale, unresolved turbulence. The dynamo process, combined with the differential rotation of the (proto-)star, is able to produce an exponential growth of any initial magnetic seed field up to the values required to explain the observations. By varying the dynamo coefficient we obtain different growth rates. We find a quasi-linear dependence of the growth rates on the intensity of the dynamo. Furthermore, the time interval in which exponential growth occurs and the growth rates also seems to depend on the initial configuration of the magnetic field., Comment: 12 pages, 4 figures, accepted for publication on Universe

Pulsars out of their parent SNR directly interact with the ISM producing so called Bow-Shock Pulsar Wind Nebulae, the relativistic equivalents of the heliosphere/heliotail system. These have been directly observed from Radio to X-ray, and are found also associated to TeV halos, with a large variety of morphologies. They offer a unique environment where the pulsar wind can be studied by modelling its interaction with the surrounding ambient medium, in a fashion that is different/complementary from the canonical Plerions. These systems have also been suggested as the possible origin of the positron excess detected by AMS and PAMELA, in contrast to dark matter. I will present results from 3D Relativistic MHD simulations of such nebulae. On top of these simulations we computed the expected emission signatures, the properties of high energy particle escape, the role of current sheets in channeling cosmic rays, the level of turbulence and magnetic amplification, and how they depend on the wind structure and magnetisation., Comment: 10 pages, 3 figures, Proceedings of ASTRONUM-2019, Accepted in Journal of Physics Conference Series

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 recent imaging of the M87 black hole at millimeter wavelengths by the Event Horizon Telescope (EHT) collaboration has triggered a renewed interest in numerical models for the accretion of magnetized plasma in the regime of general relativistic magnetohydrodynamics (GRMHD). Here non-ideal simulations, including both the resistive effects and, above all, the mean-field dynamo action due to sub-scale, unresolved turbulence, are applied for the first time to such systems in the fully nonlinear regime. Combined with the differential rotation of the disk, the dynamo process is able to produce an exponential growth of any initial seed magnetic field up to the values required to explain the observations, when the instability tends to saturate even in the absence of artificial quenching effects. Before reaching the final saturation stage we observe a secondary regime of exponential growing, where the magnetic field increases more slowly due to accretion, which is modifying the underlying equilibrium. By varying the dynamo coefficient we obtain different growth rates, though the field seems to saturate at approximately the same level, at least for the limited range of parameters explored here, providing substantial values for the MAD parameter for magnetized accretion. For reasonable values of the central mass density and the commonly employed recipes for synchrotron emission by relativistically hot electrons, our model is able to reproduce naturally the observed flux of Sgr A*, the next target for EHT., Comment: 15 pages, 13 figures

We describe a novel Godunov-type numerical method for solving the equations of resistive relativistic magnetohydrodynamics. In the proposed approach, the spatial components of both magnetic and electric fields are located at zone interfaces and are evolved using the constrained transport formalism. Direct application of Stokes' theorem to Faraday's and Ampere's laws ensures that the resulting discretization is divergence-free for the magnetic field and charge-conserving for the electric field. Hydrodynamic variables retain, instead, the usual zone-centred representation commonly adopted in finite-volume schemes. Temporal discretization is based on Runge-Kutta implicit-explicit (IMEX) schemes in order to resolve the temporal scale disparity introduced by the stiff source term in Ampere's law. The implicit step is accomplished by means of an improved and more efficient Newton-Broyden multidimensional root-finding algorithm. The explicit step relies on a multidimensional Riemann solver to compute the line-averaged electric and magnetic fields at zone edges and it employs a one-dimensional Riemann solver at zone interfaces to update zone-centred hydrodynamic quantities. For the latter, we introduce a five-wave solver based on the frozen limit of the relaxation system whereby the solution to the Riemann problem can be decomposed into an outer Maxwell solver and an inner hydrodynamic solver. A number of numerical benchmarks demonstrate that our method is superior in stability and robustness to the more popular charge-conserving divergence cleaning approach where both primary electric and magnetic fields are zone-centered. In addition, the employment of a less diffusive Riemann solver noticeably improves the accuracy of the computations., Comment: 25 pages, 14 figures

This paper briefly reviews a number of fundamental measurements that need to be made in order to characterize turbulence in space plasmas such as the solar wind. It has long been known that many of these quantities require simultaneous multipoint measurements to attain a proper characterization that would reveal the fundamental physics of plasma turbulence. The solar wind is an ideal plasma for such an investigation, and it now appears to be technologically feasible to carry out such an investigation, following the pioneering Cluster and MMS missions. Quantities that need to be measured using multipoint measurements include the two-point, two-time second correlation function of velocity, magnetic field and density, and higher order statistical objects such as third and fourth order structure functions. Some details of these requirements are given here, with a eye towards achieving closure on fundamental questions regarding the cascade rate, spectral anisotropy, characteristic coherent structures, intermittency, and dissipation mechanisms that describe plasma turbuelence, as well as its variability with plasma parameters in the solar wind. The motivation for this discussion is the current planning for a proposed Helioswarm mission that would be designed to make these measurements,leading to breakthrough understanding of the physics of space and astrophysical turbulence., Comment: White paper submitted to the PLASMA 2020 Decadal Survey Committee

Neutron stars are among the most fascinating astrophysical sources, being characterized by strong gravity, densities about the nuclear one or even above, and huge magnetic fields. Their observational signatures can be extremely diverse across the electromagnetic spectrum, ranging from the periodic and low-frequency signals of radio pulsars, up to the abrupt high-energy gamma-ray flares of magnetars, where energies of ~10^46 erg are released in a few seconds. Fast-rotating and highly magnetized neutron stars are expected to launch powerful relativistic winds, whose interaction with the supernova remnants gives rise to the non-thermal emission of pulsar wind nebulae, which are known cosmic accelerators of electrons and positrons up to PeV energies. In the extreme cases of proto-magnetars (magnetic fields of ~10^15 G and millisecond periods), a similar mechanism is likely to provide a viable engine for the still mysterious gamma-ray bursts. The key ingredient in all these spectacular manifestations of neutron stars is the presence of strong magnetic fields in their constituent plasma. Here we will present recent updates of a couple of state-of-the-art numerical investigations by the high-energy astrophysics group in Arcetri: a comprehensive modeling of the steady-state axisymmetric structure of rotating magnetized neutron stars in general relativity, and dynamical 3-D MHD simulations of relativistic pulsar winds and their associated nebulae., Comment: EPS 44th Conference on Plasma Physics (June 2017, Belfast), paper accepted for publication on Plasma Physics and Controlled Fusion

Motivated by recent observational searches of sub-structure in starless molecular cloud cores, we investigate the evolution of density perturbations on scales smaller than the Jeans length embedded in contracting isothermal clouds, adopting the same formalism developed for the expanding Universe and the solar wind. We find that initially small amplitude, Jeans-stable perturbations (propagating as sound waves in the absence of a magnetic field), are amplified adiabatically during the contraction, approximately conserving the wave action density, until they either become nonlinear and steepen into shocks at a time $t_{\rm nl}$, or become gravitationally unstable when the Jeans length decreases below the scale of the perturbations at a time $t_{\rm gr}$. We evaluate analytically the time $t_{\rm nl}$ at which the perturbations enter the non-linear stage using a Burgers' equation approach, and we verify numerically that this time marks the beginning of the phase of rapid dissipation of the kinetic energy of the perturbations. We then show that for typical values of the rms Mach number in molecular cloud cores, $t_{\rm nl}$ is smaller than $t_{\rm gr}$, and therefore density perturbations likely dissipate before becoming gravitational unstable. Solenoidal modes grow at a faster rate than compressible modes, and may eventually promote fragmentation through the formation of vortical structures., Comment: 8 pages, 4 figures

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

Pulsar Wind Nebulae (PWNe) constitute an ideal astrophysical environment to test our current understanding of relativistic plasma processes. It is well known that magnetic fields play a crucial role in their dynamics and emission properties. At present, one of the main issues concerns the level of magnetic turbulence present in these systems, which in the absence of space resolved X-ray polarization measures cannot be directly constrained. In this work we investigate, for the first time using simulated synchrotron maps, the effect of a small scale fluctuating component of the magnetic field on the emission properties in X-ray. We illustrate how to include the effects of a turbulent component in standard emission models for PWNe, and which consequences are expected in terms of net emissivity and depolarization, showing that the X-ray surface brightness maps can provide already some rough constraints. We then apply our analysis to the Crab and Vela nebulae and, by comparing our model with Chandra and Vela data, we found that the typical energies in the turbulent component of the magnetic field are about 1.5 to 3 times the one in the ordered field., Comment: 9 pages, 8 figures, accepted for publication in MNRAS

In the last decade, the relativistic magnetohydrodynamic (MHD) modelling of pulsar wind nebulae, and of the Crab nebula in particular, has been highly successful, with many of the observed dynamical and emission properties reproduced down to the finest detail. Here, we critically discuss the results of some of the most recent studies: namely the investigation of the origin of the radio emitting particles and the quest for the acceleration sites of particles of different energies along the termination shock, by using wisps motion as a diagnostic tool; the study of the magnetic dissipation process in high magnetization nebulae by means of new long-term three-dimensional simulations of the pulsar wind nebula evolution; the investigation of the relativistic tearing instability in thinning current sheets, leading to fast reconnection events that might be at the origin of the Crab nebula gamma-ray flares., Comment: 30 pages, 12 figures

Magnetic reconnection in laboratory, space and astrophysical plasmas is often invoked to explain explosive energy release and particle acceleration. However, the timescales involved in classical models within the macroscopic MHD regime are far too slow to match the observations. Here we revisit the tearing instability by performing visco-resistive two-dimensional numerical simulations of the evolution of thin current sheets, for a variety of initial configurations and of values of the Lunquist number $S$, up to $10^7$. Results confirm that when the critical aspect ratio of $S^{1/3}$ is reached in the reconnecting current sheets, the instability proceeds on ideal (Alfv\'enic) macroscopic timescales, as required to explain observations. Moreover, the same scaling is seen to apply also to the local, secondary reconnection events triggered during the nonlinear phase of the tearing instability, thus accelerating the cascading process to increasingly smaller spatial and temporal scales. The process appears to be robust, as the predicted scaling is measured both in inviscid simulations and when using a Prandtl number $P=1$ in the viscous regime., Comment: Accepted for publication in Plasma Physics and Controlled Fusion

We investigate the possible implications of quark deconfinement on the phenomenology of Long Gamma-Ray Bursts focusing, in particular, on the possibility to describe multiple prompt emission phases in the context of the proto-magnetar model. Starting from numerical models of rotating Hadron Stars and Quark Stars in full general relativity we track the electromagnetic spin-down evolution in both the hadronic and quark phase, linking the two families through conservation of baryon number and angular momentum. We give estimates of the timescales and the energetics involved in the spin-down process deriving, in the relevant spin range, the relation between the initial and the final masses and rotational energies, whenever hadron-quark conversion is possible. We show how the results can be used in relevant astrophysical cases such as the double burst GRB 110709B., Comment: 5 pages, 3 figures, 2 tables

Fast reconnection operating in magnetically dominated plasmas is often invoked in models for magnetar giant flares, for magnetic dissipation in pulsar winds, or to explain the gamma-ray flares observed in the Crab nebula, hence its investigation is of paramount importance in high-energy astrophysics. Here we study, by means of two dimensional numerical simulations, the linear phase and the subsequent nonlinear evolution of the tearing instability within the framework of relativistic resistive magnetohydrodynamics, as appropriate in situations where the Alfven velocity approaches the speed of light. It is found that the linear phase of the instability closely matches the analysis in classical MHD, where the growth rate scales with the Lundquist number S as S^-1/2, with the only exception of an enhanced inertial term due to the thermal and magnetic energy contributions. In addition, when thin current sheets of inverse aspect ratio scaling as S^-1/3 are considered, the so-called "ideal" tearing regime is retrieved, with modes growing independently on S and extremely fast, on only a few light crossing times of the sheet length. The overall growth of fluctuations is seen to solely depend on the value of the background Alfven velocity. In the fully nonlinear stage we observe an inverse cascade towards the fundamental mode, with Petschek-type supersonic jets propagating at the external Alfven speed from the X-point, and a fast reconnection rate at the predicted value R~(ln S)^-1., Comment: 14 pages, 9 figures, accepted for publication (MNRAS)

Classical MHD reconnection theories, both the stationary Sweet-Parker model and the tearing instability, are known to provide rates which are too slow to explain the observations. However, a recent analysis has shown that there exists a critical threshold on current sheet's thickness, namely a/L~S^(-1/3), beyond which the tearing modes evolve on fast macroscopic Alfvenic timescales, provided the Lunquist number S is high enough, as invariably found in solar and astrophysical plasmas. Therefore, the classical Sweet-Parker scenario, for which the diffusive region scales as a/L~S^(-1/2) and thus can be up to ~100 times thinner than the critical value, is likely to be never realized in nature, as the current sheet itself disrupts in the elongation process. We present here two-dimensional, compressible, resistive MHD simulations, with S ranging from 10^5 to 10^7, that fully confirm the linear analysis. Moreover, we show that a secondary plasmoid instability always occurs when the same critical scaling is reached on the local, smaller scale, leading to a cascading explosive process, reminiscent of the flaring activity., Comment: Proceedings of ASTRONUM 2015, Avignon, France (11 pages, 6 figures)

The phenomenology of the emission of pulsars and magnetars depends dramatically on the structure and properties of their magnetic field. In particular it is believed that the outbursting and flaring activity observed in AXPs and SRGs is strongly related to their internal magnetic field. Recent observations have moreover shown that charges are present in their magnetospheres supporting the idea that their magnetic field is tightly twisted in the vicinity of the star. In principle these objects offer a unique opportunity to investigate physics in a regime beyond what can be obtained in the laboratory. We will discuss the properties of equilibrium models of magnetized neutron stars, and we will show how internal and external currents can be related. These magnetic field configurations will be discussed considering also their stability, relevant for their origin and possibly connected to events like SNe and GRBs. We will also show what kind of deformations they induce in the star, that could lead to emission of gravitational waves. In the case of a twisted magnetosphere we will show how the amount of twist regulates their general topology. A general formalism based on the simultaneous numerical solution of the general relativistic Grad-Shafranov equation and Einstein equations will be presented., Comment: 9 pages, 3 figures, Proceedings of the 10th International Conference on Numerical Modeling of Space Plasma Flows, 8-12 June 2015, Avignon, France

We study the impact of three stochastic parametrizations in the ocean component of a coupled model, on forecast reliability over seasonal timescales. The relative impacts of these schemes upon the ocean mean state and ensemble spread are analyzed. The oceanic variability induced by the atmospheric forcing of the coupled system is, in most regions, the major source of ensemble spread. The largest impact on spread and bias came from the Stochastically Perturbed Parametrization Tendency (SPPT) scheme - which has proven particularly effective in the atmosphere. The key regions affected are eddy-active regions, namely the western boundary currents and the Southern Ocean. However, unlike its impact in the atmosphere, SPPT in the ocean did not result in a significant decrease in forecast error. Whilst there are good grounds for implementing stochastic schemes in ocean models, our results suggest that they will have to be more sophisticated. Some suggestions for next-generation stochastic schemes are made., Comment: 24 pages, 3 figures