High Energy Astrophysical Phenomena (astro-ph.HE)

Abstract: This paper employs an MHD-PIC method to perform numerical simulations of magnetic reconnection-driven turbulence and turbulent reconnection acceleration of particles. Focusing on the dynamics of the magnetic reconnection, the properties of self-driven turbulence, and the behavior of particle acceleration, we find that: (1) when reaching a statistically steady state of the self-driven turbulence, the magnetic energy is almost released by 50\%, while the kinetic energy of the fluid increases by no more than 15\%. (2) the properties of reconnection-driven turbulence are more complex than the traditional turbulence driven by an external force. (3) the strong magnetic field tends to enhance the turbulent reconnection efficiency to accelerate particles more efficiently, resulting in a hard spectral energy distribution. Our study provides a particular perspective on understanding turbulence properties and turbulent reconnection-accelerated particles.

Abstract: We carry out a stacked search for spatial coincidences between all known radio pulsars and TeV neutrinos from the IceCube 10 year (2008-2018) point source catalog as a followup to our previous work on looking for coincidences with individual pulsars. We consider three different weighting schemes to stack the contribution from individual pulsars. We do not find a statistically significant excess using this method. We report the 95% c.l. differential neutrino flux limit as a function of neutrino energy. We have also made our analysis codes publicly available.

Abstract: Timing noise in pulsars is often modelled with a Fourier-basis Gaussian process that follows a power law with periodic boundary conditions on the observation time, $T_\mathrm{span}$. However the actual noise processes can extend well below $1/T_\mathrm{span}$, and many pulsars are known to exhibit quasi-periodic timing noise. In this paper we investigate several adaptions that try to account for these differences between the observed behaviour and the simple power-law model. Firstly, we propose to include an additional term that models the quasi-periodic spin-down variations known to be present in many pulsars. Secondly, we show that a Fourier basis of $1/2T_\mathrm{span}$ can be more suited for estimating long term timing parameters such as the spin frequency second derivative (F2), and is required when the exponent of the power spectrum is greater than ~4. We also implement a Bayesian version of the generalised least squares `Cholesky' method which has different limitations at low frequency, but find that there is little advantage over Fourier-basis methods. We apply our quasi-periodic spin down model to a sample of pulsars with known spin-down variations and show that this improves parameter estimation of F2 and proper motion for the most pathological cases, but in general the results are consistent with a power-law model. The models are all made available through the run_enterprise software package.

Abstract: X-Ray bursts (XRB) are powerful thermonuclear events on the surface of accreting neutron stars (NS), where nucleosynthesis of intermediate-mass elements occurs. Their predicted and observed luminosities sometimes exceed Eddington's value, thus some of the material may escape by means of a stellar wind. This work seeks to determine the mass-loss and chemical composition of the material ejected through radiation-driven winds and its significance for Galactic abundances. It also reports on the evolution of pysical quantities during the wind phase that could help constrain the mass-radius relation in neutron stars. A non-relativistic radiative wind model was implemented and linked, through a new technique, to a series of XRB hydrodynamic simulations, that include over 300 isotopes. This allows us to construct a quasi-stationary time evolution of the wind during the XRB. The simulations resulted in the first realistic quantification of mass-loss for each isotope synthesized in the XRB. The total mass ejected by the wind was about $6\times10^{19}g$, the average ejected mass per unit time represents 2.6% of the accretion rate, with 0.1% of the envelope mass ejected per burst and ~90% of the ejecta composed by $^{60}$Ni, $^{64}$Zn, $^{68}$Ge and $^{58}$Ni. The ejected material also contained a small fraction ($10^{-4}-10^{-5}$) of some light p-nuclei, but not enough to account for their Galactic abundances. Additionally, the observable magnitudes during the wind phase showed remarkable correlations, some of which involve wind parameters like energy and mass outflows, that are determined by the conditions at the base of the wind envelope. These correlations could be used to link observable magnitudes to the physics of the innermost parts of the envelope, close to its interface with the NS crust. This is a promising result regarding the issue of NS radii determination.

Abstract: We investigate the acceleration of cosmic rays at the termination shock that results from the interaction of the collective wind of star clusters with the surrounding interstellar medium. The solution of the transport equation of accelerated particles in the wind-excavated cavity, including energy losses due to CR interactions with neutral gas in the bubble, shows several interesting properties that are discussed in detail. The issue of the maximum energy of the accelerated particles is discussed with special care, because of its implications for the origin of Galactic cosmic rays. Gamma ray emission is produced in the cavity due to inelastic pp scattering, while accelerated particles are advected downstream of the termination shock and diffuse at the same time. Both the spectrum and the morphology of such emission are discussed, with a comparison of our results with the observations of gamma ray emission from the Cygnus OB2 region.

Abstract: Cosmological simulations predict the presence of warm hot thermal gas in the cosmic filaments that connect galaxy clusters. This gas is thought to constitute an important part of the missing baryons in the Universe. In addition to the thermal gas, cosmic filaments could contain a population of relativistic particles and magnetic fields. A detection of magnetic fields in filaments can constrain early magnetogenesis in the cosmos. So far, the resulting diffuse synchrotron emission has only been indirectly detected. We present our search for thermal and non-thermal diffuse emission from inter-cluster regions of 106 paired galaxy clusters by stacking the $0.6-2.3$~keV X-ray and 144~MHz radio data obtained with the eROSITA telescope on board the Spectrum-Roentgen-Gamma (SRG) observatory and LOw Frequency ARray (LOFAR), respectively. The stacked data do not show the presence of X-ray and radio diffuse emission in the inter-cluster regions. This could be due to the sensitivity of the data sets and/or the limited number of cluster pairs used in this study. Assuming a constant radio emissivity in the filaments, we find that the mean radio emissivity is not higher than $1.2\times10^{-44}\,{\rm erg \, s^{-1} \, cm^{-3} \, Hz^{-1}}$. Under equipartition conditions, our upper limit on the mean emissivity translates to an upper limit of $\sim75\,{\rm nG}$ for the mean magnetic field strength in the filaments, depending on the spectral index and the minimum energy cutoff. We discuss the constraint for the magnetic field strength in the context of the models for the formation of magnetic fields in cosmic filaments.