High Energy Astrophysical Phenomena (astro-ph.HE)

Abstract: The aim of the current paper is to apply the method of Bambi (Bambi, 2015) to a source which contains two or more simultaneous triads of variability components. The joint chi-square variable that can be composed in this case, unlike some previous studies, allows the goodness of the fit to be tested. It appears that a good fit requires one of the observation groups to be disregarded. Even then, the model prediction for the mass of the neutron star in the accreting millisecond pulsar IGR J17511-3057 is way too high to be accepted.

Abstract: Using multi-frequency Very Long Baseline Interferometer (VLBI) observations, we probe the jet size in the optically thick hard state jets of two black hole X-ray binary (BHXRB) systems, MAXI J1820+070 and V404 Cygni. Due to optical depth effects, the phase referenced VLBI core positions move along the jet axis of the BHXRB in a frequency dependent manner. We use this "core shift" to constrain the physical size of the hard state jet. We place an upper limit of $0.3$\,au on the jet size measured between the 15 and 5 GHz emission regions of the jet in MAXI J1820+070, and an upper limit of $1.0$\,au between the $8.4$ and $4.8$\,GHz emission regions of V404 Cygni. Our limit on the jet size in MAXI J1820+070 observed in the low-hard state is a factor of $5$ smaller than the values previously observed in the high-luminosity hard state (using time lags between multi-frequency light curves), thus showing evidence of the BHXRB jet scaling in size with jet luminosity. We also investigate whether motion of the radio-emitting region along the jet axis could affect the measured VLBI parallaxes for the two systems, leading to a mild tension with the parallax measurements of Gaia. Having mitigated the impact of any motion along the jet axis in the measured astrometry, we find the previous VLBI parallax measurements of MAXI J1820+070 and V404 Cygni to be unaffected by jet motion. With a total time baseline of $8$ years, due to having incorporated fourteen new epochs in addition to the previously published ones, our updated parallax measurement of V404 Cygni is $0.450 \pm 0.018$\,mas ($2.226 \pm 0.091$\,kpc).

Abstract: Highly magnetized neutron stars have magnetic fields of order of the critical field and can lead to measurable QED effects. We consider the Goldreich-Julian pulsar model with supercritical magnetic fields, induced subcritical electric fields, and a period of milliseconds. We then study the strong field physics, such as Schwinger pair production and the vacuum birefringence including the wrench effect, whose X-ray polarimetry will be observed in future space missions.

Abstract: As one of the brightest galactic ${\gamma}$-ray sources, the Cygnus Cocoon superbubble has been observed by many detectors, such as $Fermi$-LAT, ARGO, HAWC, and LHAASO. However, the origin of $\gamma$-ray emission for the Cygnus Cocoon and the possible contribution to PeV cosmic rays are still under debate. The recent ultrahigh-energy $\gamma$-ray observations by LHAASO up to 1.4 PeV towards the direction of the Cygnus Cocoon, as well as the neutrino event report of IceCube-201120A coming from the same direction, suggest that the Cygnus Cocoon may be one of the sources of high-energy cosmic rays in the Galaxy. In this work, we propose a dual-zone diffusion model for the Cygnus Cocoon: the cocoon region and surrounding interstellar medium (ISM). This scenario can account for the $\gamma$-ray data from GeV to $\sim$ 50 TeV and agree with the one sub-PeV neutrino event result from IceCube so far. Moreover, it predict a non-negligible contribution $\gamma$-ray emission at hundreds TeV from the ISM surrounding the Cygnus Cocoon. This possible diffuse TeV-PeV gamma-ray features can be resolved by the future LHAASO observations.

Abstract: Acceleration processes that occur in astrophysical plasmas produce cosmic rays that are observed on Earth. To study particle acceleration, fully-kinetic particle-in-cell (PIC) simulations are often used as they can unveil the microphysics of energization processes. Tracing of individual particles in PIC simulations is particularly useful in this regard. However, by-eye inspection of particle trajectories includes a high level of bias and uncertainty in pinpointing specific acceleration mechanisms that affect particles. Here we present a new approach that uses neural networks to aid individual particle data analysis. We demonstrate this approach on the test data that consists of 252,000 electrons which have been traced in a PIC simulation of a non-relativistic high Mach number perpendicular shock, in which we observe the two-stream electrostatic Buneman instability to pre-accelerate a portion of electrons to nonthermal energies. We perform classification, regression and anomaly detection by using a Convolutional Neural Network. We show that regardless of how noisy and imbalanced the datasets are, the regression and classification are able to predict the final energies of particles with high accuracy, whereas anomaly detection is able to discern between energetic and non-energetic particles. The methodology proposed may considerably simplify particle classification in large-scale PIC and also hybrid kinetic simulations.

Abstract: Recent observations of cosmic rays (CRs) have revealed a two-component anomaly in the spectra of primary and secondary particles, as well as their ratios, prompting investigation into their common origin. In this study, we incorporate the identification of slow diffusion zones around sources as a common phenomenon into our calculations, which successfully reproduces all previously described anomalies except for the positron spectrum. Crucially, our research offers a clear physical picture of the origin of CR: while high-energy ($\textrm{>200~GV}$, including the knee) particles are primarily produced by fresh accelerators and are confined to local regions, low energy ($\textrm{<200~GV}$) components come from distant sources and travel through the outer diffusive zone outside of the galactic disk. This scenario can be universally applied in the galactic disk, as evidenced by ultra-high energy diffuse $\rm\gamma$-ray emissions detected by the AS$\rm\gamma$ experiment. Furthermore, our results predict that the spectrum of diffuse $\rm\gamma$-ray is spatial-dependent, resting with local sources, which can be tested by LHAASO experiment.

Abstract: Synchrotron radiation from relativistic electrons is usually invoked as the responsible for the nonthermal emission observed in Supernova Remnants (SNRs). Diffusive shock acceleration (DSA) is the most popular mechanism to explain the process of particles acceleration and within its framework a crucial role is played by the turbulent magnetic-field. However, the standard models commonly used to fit X-ray synchrotron emission do not take into account the effects of turbulence in the shape of the resulting photon spectra. An alternative mechanism that properly includes such effects is the jitter radiation, that provides for an additional power-law beyond the classical synchrotron cutoff. We fitted a jitter spectral model to Chandra, NuSTAR, SWIFT/BAT and INTEGRAL/ISGRI spectra of Cassiopeia A and found that it describes the X-ray soft-to-hard range better than any of the standard cutoff models. The jitter radiation allows us to measure the index of the magnetic turbulence spectrum $\nu_B$ and the minimum scale of the turbulence $\lambda_{\rm{min}}$ across several regions of Cas A, with best-fit values $\nu_B \sim 2-2.4$ and $\lambda_{\rm{min}} \lesssim 100$ km.

Abstract: During substructure formation in magnetized astrophysical plasma, dissipation of magnetic energy facilitated by magnetic reconnection affects the system dynamics by heating and accelerating the ejected plasmoids. Numerical simulations are a crucial tool for investigating such systems. In astrophysical simulations, the energy dissipation, reconnection rate and substructure formation critically depend on the onset of reconnection of numerical or physical origin. In this paper, we hope to assess the reliability of the state-of-the-art numerical codes, PLUTO and KORAL by quantifying and discussing the impact of dimensionality, resolution, and code accuracy on magnetic energy dissipation, reconnection rate, and substructure formation. We quantitatively compare results obtained with relativistic and non-relativistic, resistive and non-resistive, as well as two- and three-dimensional setups performing the Orszag-Tang test problem. We find the sufficient resolution in each model, for which numerical error is negligible and the resolution does not significantly affect the magnetic energy dissipation and reconnection rate. The non-relativistic simulations show that at sufficient resolution, magnetic and kinetic energies convert to internal energy and heat up the plasma. The results show that in the relativistic system, energy components undergo mutual conversion during the simulation time, which leads to a substantial increase in magnetic energy at 20\% and 90\% of the total simulation time of $10$ light-crossing times -- the magnetic field is amplified by a factor of five due to relativistic shocks. We also show that the reconnection rate in all our simulations is higher than $0.1$, indicating plasmoid-mediated regime. It is shown that in KORAL simulations magnetic energy is slightly larger and more substructures are captured than in PLUTO simulations.

Abstract: Recent observations of old warm neutron stars suggest the presence of a heating source in these stars, requiring a paradigm beyond the standard neutron-star cooling theory. In this work, we study the scenario where this heating is caused by the friction associated with the creep motion of neutron superfluid vortex lines in the crust. As it turns out, the heating luminosity in this scenario is proportional to the time derivative of the angular velocity of the pulsar rotation, and the proportional constant $J$ has an approximately universal value for all neutron stars. This $J$ parameter can be determined from the temperature observation of old neutron stars because the heating luminosity is balanced with the photon emission at late times. We study the latest data of neutron star temperature observation and find that these data indeed give similar values of $J$, in favor of the assumption that the frictional motion of vortex lines heats these neutron stars. These values turn out to be consistent with the theoretical calculations of the vortex-nuclear interaction.

Abstract: The composition of relativistic gamma-ray burst (GRB) jets and their emission mechanisms are still debated, and they could be matter or magnetically dominated. One way to distinguish these mechanisms arises because a Poynting flux dominated jet may produce low-frequency radio emission during the energetic prompt phase, through magnetic reconnection at the shock front. We present a search for radio emission coincident with three GRB X-ray flares with the LOw Frequency ARray (LOFAR), in a rapid response mode follow-up of long GRB 210112A (at z~2) with a 2 hour duration, where our observations began 511 seconds after the initial swift-BAT trigger. Using timesliced imaging at 120-168 MHz, we obtain upper limits at 3 sigma confidence of 42 mJy averaging over 320 second snapshot images, and 87 mJy averaging over 60 second snapshot images. LOFAR's fast response time means that all three potential radio counterparts to X-ray flares are observable after accounting for dispersion at the estimated source redshift. Furthermore, the radio pulse in the magnetic wind model was expected to be detectable at our observing frequency and flux density limits which allows us to disfavour a region of parameter space for this GRB. However, we note that stricter constraints on redshift and the fraction of energy in the magnetic field are required to further test jet characteristics across the GRB population.

Abstract: Searches for gravitational waves from compact-binary mergers, which to date have reported nearly 100 observations, have previously ignored binaries whose components are both light ($< 2 M_\odot$) and have high dimensionless spin ($>$ 0.05). While previous searches targeted sources that are representative of observed double neutron star binaries in the galaxy, it is already known that neutron stars can regularly be spun up to a dimensionless spin of $\sim$ 0.4, and in principle reach up to $\sim$ 0.7 before breakup would occur. Furthermore, there may be primordial black hole binaries or exotic formation mechanisms to produce subsolar mass black holes. In these cases, it is possible for the binary constituent to be spun up beyond that achievable by a neutron star. A single detection of this type of source would reveal a novel formation channel for compact binaries. To determine if there is evidence for any such sources, we conduct the first search of LIGO and Virgo data from three observing runs for light compact objects with high spin. Our analysis detects previously known observations e.g. GW170817 and GW200115; however, we report no additional mergers. The most significant candidate, not previously known, is consistent with the noise distribution, and so we constrain the merger rate of spinning light binaries.