Fast reconnection in relativistic plasmas: the magnetohydrodynamics tearing instability revisited

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

Contribution to Galactic cosmic rays from young stellar clusters

The origin of Galactic cosmic rays (CR) is still a matter of debate. Diffusive shock acceleration (DSA) applied to supernova remnant (SNR) shocks provides the most reliable explanation. However, within the current understanding of DSA several issues remain unsolved, like the CR maximum energy, the chemical composition and the transition region between Galactic and extra-Galactic CRs. These issues motivate the search for other possible Galactic sources. Recently, several young stellar clusters (YSC) have been detected in gamma rays, suggesting that such objects could be powerful sources of Galactic CRs. The energy input could come from winds of massive stars hosted in the clusters which is a function of the cluster total mass and initial mass function of stars. In this work we evaluate the total CR flux produced by a synthetic population of YSCs assuming that the CR acceleration occurs at the termination shock of the collective wind resulting from the sum of cluster's stellar winds. We show that the spectrum produced by YSC can significantly contribute to energies $\gtrsim 100$ TeV if the diffusion inside the wind-blown bubble is Bohm-like and the spectral slope is harder than the one produced by SNRs.Comment: Proceeding to the International Cosmic Ray Conference, ICRC 2023, Nagoya, Japa

Rarefaction acceleration of ultrarelativistic magnetized jets in gamma-ray burst sources

When a magnetically-dominated super-fast magnetosonic GRB jet leaves the progenitor star the external pressure support may drop and the jet may enter the regime of ballistic expansion during which its magnetic acceleration becomes highly ineffective. However, recent numerical simulations suggested that the transition to this regime is accompanied by a sudden "burst" of acceleration. We confirm this finding and attribute the acceleration to the sideways expansion of the jet - the magnetic energy is converted into the kinetic one in the strong magnetosonic rarefaction wave, which is launched when the jet loses its external support. This type of acceleration, the rarefaction acceleration, is specific to relativistic jets because their energy budget can still be dominated by magnetic energy even in highly super-fast magnetosonic regime. Just like the collimation acceleration of externally confined magnetized jets, it is connected with the geometry of magnetic flux sufaces. In both cases, in the acceleration zone the poloidal field lines diverge faster than in the monopolar configuration. On the other hand, whereas the collimation acceleration keeps the product of jet opening angle and Lorentz factor somewhat below unity, the rarefaction acceleration allows to make it significantly larger, in agreement with the standard model of jet breaks in afterglow light curves.Comment: Submitted to MNRA

The ECHO code for astrophysical plasmas: Special and General Relativistic MHD

The main features of a novel numerical code for astrophysical fluids and magneto-fluids named ECHO (Eulerian conservative high-order) are presented. Here the module for special and general relativistic MHD is discussed and applications of the code to astrophysical problems are summarized

GRB Afterglows with Energy Injection from a spinning down NS

We investigate a model for the shallow decay phases of Gamma-ray Burst (GRB) afterglows discovered by Swift/XRT in the first hours following a GRB event. In the context of the fireball scenario, we consider the possibility that long-lived energy injection from a millisecond spinning, ultramagnetic neutron star (magnetar) powers afterglow emission during this phase. We consider the energy evolution in a relativistic shock subject to both radiative losses and energy injection from a spinning down magnetar in spherical symmetry. We model the energy injection term through magnetic dipole losses and discuss an approximate treatment for the dynamical evolution of the blastwave. We obtain an analytic solution for the energy evolution in the shock and associated lightcurves. To fully illustrate the potential of our solution we calculate lightcurves for a few selected X-ray afterglows observed by Swift and fit them using our theoretical lightcurves. Our solution naturally describes in a single picture the properties of the shallow decay phase and the transition to the so-called normal decay phase. In particular, we obtain remarkably good fits to X-ray afterglows for plausible parameters of the magnetar. Even though approximate, our treatment provides a step forward with respect to previously adopted approximations and provides additional support to the idea that a millisecond spinning (1-3 ms), ultramagnetic (B$\sim 10^{14}-10^{15}$ G) neutron star loosing spin energy through magnetic dipole radiation can explain the luminosity, durations and shapes of X-ray GRB afterglows.Comment: 7 pages, 2 figures, submitted to Astronomy & Astrophysics - referee's comments include

MHD models of Pulsar Wind Nebulae

Pulsar Wind Nebulae (PWNe) are bubbles or relativistic plasma that form when the pulsar wind is confined by the SNR or the ISM. Recent observations have shown a richness of emission features that has driven a renewed interest in the theoretical modeling of these objects. In recent years a MHD paradigm has been developed, capable of reproducing almost all of the observed properties of PWNe, shedding new light on many old issues. Given that PWNe are perhaps the nearest systems where processes related to relativistic dynamics can be investigated with high accuracy, a reliable model of their behavior is paramount for a correct understanding of high energy astrophysics in general. I will review the present status of MHD models: what are the key ingredients, their successes, and open questions that still need further investigation.Comment: 18 pages, 5 figures, Invited Review, Proceedings of the "ICREA Workshop on The High-Energy Emission from Pulsars and their Systems", Sant Cugat, Spain, April 12-16, 201

Magnetic Field Structure of the Crab Pulsar Wind Nebula Revealed with IXPE

We report a detailed study of the magnetic-field structure of the Crab pulsar wind nebula, using the X-ray polarization data in 2--8~keV obtained with the Imaging X-ray Polarimetry Explorer. Contamination of the pulsar emission to the data of the nebula region was removed through application of a stringent pulsation phase-cut, extracting a phase range of 0.7--1.0 only. We found that the electric field vector polarization angle (PA) was about $130^{\circ}$ from north to east with the polarization degree (PD) of about 25\% at the pulsar position, indicating that the direction of the toroidal magnetic field is perpendicular to the pulsar spin axis in the region close to the termination shock. The PA gradually deviated from the angle as an increasing function of the distance from the pulsar. There was a region of a low PD to the west of the X-ray torus. Although such a region is expected to be located at the torus edge, where geometrical depolarization due to a steep spatial variation of the PA is expected, the observed low-PD region positionally deviated from the edge. We found that the region of low PD positionally coincided with a dense filament seen in the optical band, and conjecture that the low-PD region may be produced through deflection of the pulsar wind. By comparing the values of the PD at the pulsar position between the data and a model, in which toroidal and turbulent magnetic fields were considered, we estimated the fractional energy of the turbulent magnetic field to be about $2/3$ of the total. We also evaluated a potential polarization of the northern jet in the nebula and derived the PD and PA to be about $30\%$ and $120^{\circ}$, respectively.Comment: This is a pre-copyedited, author-produced version of an article accepted for publication in PASJ following peer revie

Non-stationary Rayleigh-Taylor instability in supernovae ejecta

The Rayleigh-Taylor instability plays an important role in the dynamics of several astronomical objects, in particular, in supernovae (SN) evolution. In this paper we develop an analytical approach to study the stability analysis of spherical expansion of the SN ejecta by using a special transformation in the co-moving coordinate frame. We first study a non-stationary spherical expansion of a gas shell under the pressure of a central source. Then we analyze its stability with respect to a no radial, non spherically symmetric perturbation of the of the shell. We consider the case where the polytropic constant of the SN shell is $\gamma=5/3$ and we examine the evolution of a arbitrary shell perturbation. The dispersion relation is derived. The growth rate of the perturbation is found and its temporal and spatial evolution is discussed. The stability domain depends on the ejecta shell thickness, its acceleration, and the perturbation wavelength.Comment: 16 page

New Studies of the Pulsar Wind Nebula in the Supernova Remnant CTB 80

We investigated the kinematics of the pulsar wind nebula (PWN) associated with PSR B1951+32 in the old supernova remnant CTB 80 using the Fabry-Perot interferometer of the 6m Special Astrophysical Observatory telescope. In addition to the previously known expansion of the system of bright filaments with a velocity of 100-200km/s, we detected weak high-velocity features in the H-alpha line at least up to velocities of 400-450km/s. We analyzed the morphology of the PWN in the H-alpha, [SII], and [OIII] lines using HST data and discuss its nature. The shape of the central filamentary shell, which is determined by the emission in the [OIII] line and in the radio continuum, is shown to be consistent with the bow-shock model for a significant (about 60 degrees) inclination of the pulsar's velocity vector to the plane of the sky. In this case, the space velocity of the pulsar is twice higher than its tangential velocity, i.e., it reaches ~500 km/s, and PSR B1951+32 is the first pulsar whose line-of-sight velocity (of about 400 km/s) has been estimated from the PWN observations. The shell-like H-alpha-structures outside the bow shock front in the east and the west may be associated with both the pulsar's jets and the pulsar-wind breakthrough due to the layered structure of the extended CTB 80 shell.Comment: to appear in Astronomy Letters, 12 pages, 6 postscript figures, two in colour; for a version with high resolution figures see http://www.sao.ru/hq/grb/team/vkom/CTB80_fine.pd