Particle Acceleration Processes Research Articles

Rising Near-ultraviolet Spectra in Stellar Megaflares

Flares from M dwarf stars can attain energies up to 104 times larger than solar flares but are generally thought to result from similar processes of magnetic energy release and particle acceleration. Larger heating rates in the low atmosphere are needed to reproduce the shape and strength of the observed continua in stellar flares, which are often simplified to a blackbody model from the optical to the far-ultraviolet (FUV). The near-ultraviolet (NUV) has been woefully undersampled in spectral observations despite this being where the blackbody radiation should peak. We present Hubble Space Telescope NUV spectra in the impulsive phase of a flare with E TESS ≈ 7.5 × 1033 erg and a flare with E TESS ≈ 1035 erg and the largest NUV flare luminosity observed to date from an M star. The composite NUV spectra are not well represented by a single blackbody that is commonly assumed in the literature. Rather, continuum flux rises toward shorter wavelengths into the FUV, and we calculate that an optical T = 104 K blackbody underestimates the short-wavelength NUV flux by a factor of ≈6. We show that rising NUV continuum spectra can be reproduced by collisionally heating the lower atmosphere with beams of E ≳ 10 MeV protons or E ≳ 500 keV electrons and flux densities of 1013 erg cm−2 s−1. These are much larger than the canonical values describing accelerated particles in solar flares.

Fermi-LAT Detection of the Low-luminosity Radio Galaxy NGC 4278 during the LHAASO Campaign

Abstract We present a study of the high-energy properties of the compact symmetric object NGC 4278, recently associated with a TeV source by the Large High Altitude Air Shower Observatory (LHAASO) collaboration. We conducted a dedicated analysis of a Fermi Large Area Telescope (LAT) region around NGC 4278, limited to the LHAASO campaign conducted from 2021 March to 2022 October. A statistically significant emission (∼4.3σ) was revealed, spatially consistent with the radio position of NGC 4278 and the LHAASO source. The Fermi-LAT source is detected above 8 GeV, exhibiting a hard spectrum (photon index Γ = 1.3 ± 0.3) and a γ-ray luminosity of L >100 MeV ≃ 4 × 1041 erg s–1. A serendipitous Swift X-Ray Telescope (XRT) observation of NGC 4278 during the TeV campaign reveals the source in a high state, with a flux F 0.5 − 8 ke V = 5 − 2 + 3 × 1 0 − 12 erg s − 1 cm − 2 , compatible with the highest luminosity level observed in previous Chandra pointings. The high-energy spectral energy distribution of the source and the intense flux variation observed in the X-ray band support a jet origin for the observed radiation. We suggest that the significant enhancement of the high-energy flux observed during the LHAASO campaign is due to a transient, highly energetic perturbation in the jet. The detection of NGC 4278 at both high and very high energies opens new frontiers in studying particle acceleration processes. It reveals that even compact, low-power radio galaxies can exceed the sensitivity thresholds of GeV and TeV instruments, becoming promising targets for the upcoming Cherenkov Telescope Array Observatory.

The MeerKAT view on Galactic supernova remnants

The integrated radio spectrum of supernova remnants (SNRs) and the spatial variation of the spectral indices across these extended sources are powerful tools for studying the shocks and particle acceleration processes occurring in different SNR regions. Characterization of these processes requires sensitive flux density measurements and high-resolution images, which are not always available due to observing difficulties. We want to show the potentiality of the high-resolution SARAO MeerKAT legacy Galactic Plane Survey (SMGPS) images regarding the morphological and spectral characterization of 29 known galactic SNRs. We used the SMGPS data at $1.284$ GHz coupled with data from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey ($0.072-0.231$ GHz) to characterize the integrated spectrum of each source and search for spatial spectral variation through analysis of sensitive spectral index maps. We were able to redefine the exact morphology of four SNRs (G024.7-00.6, G051.4+00.7, G348.7+0.3, and G351.9+00.1), distinguishing them from unrelated sources or identifying new emission regions associated with them and never observed before. In many other cases, we identified in the SMGPS images several overlaid with the remnants, and we were able to estimate their spectral contribution through inspection of the spatial variation of the spectral indices across the remnants. The integrated spectral indices show a more uniform distribution with respect to what is obtained by considering the values reported in the literature. We show that new sensitive and high-resolution data are crucial to firmly constraining both the integrated and spatially resolved spectrum of SNRs, especially for the less studied objects of the southern hemisphere. The comparison of our SMGPS-GLEAM spectral index maps with IR, molecular, and gamma -ray images allowed us to investigate the nature of the peculiar remnant regions.

Novel scaling laws to derive spatially resolved flare and CME parameters from sun-as-a-star observables

Coronal mass ejections (CMEs) are often associated with X-ray (SXR) flares powered by magnetic reconnection in the low corona, while the CME shocks in the upper corona and interplanetary (IP) space accelerate electrons often producing the type II radio bursts. The CME and the reconnection event are part of the same energy release process as highlighted by the correlation between reconnection flux (ϕrec) that quantifies the strength of the released magnetic free energy during the SXR flare, and the CME kinetic energy that drives the IP shocks leading to type II bursts. Unlike the Sun, these physical parameters cannot be directly inferred in stellar observations. Hence, scaling laws between unresolved sun-as-a-star observables, namely SXR luminosity (LX) and type II luminosity (LR), and the physical properties of the associated dynamical events are crucial. Such scaling laws also provide insights into the interconnections between the particle acceleration processes across low-corona to IP space during solar-stellar “flare-CME-type II” events. Using long-term solar data in the SXR to radio waveband, we derived a scaling law between two novel power metrics for the flare and CME-associated processes. The metrics of “flare power” (Pflare = √(LXϕrec)) and “CME power” (PCME = √(LRVCME2)), where VCME is the CME speed, scale as Pflare ∝ PCME0.76 ± 0.04. In addition, LX and ϕrec show power-law trends with PCME with indices of 1.12 ± 0.05 and 0.61 ± 0.05, respectively. These power laws help infer the spatially resolved physical parameters, VCME and ϕrec, from disk-averaged observables, LX and LR during solar-stellar flare-CME-type II events.

The remarkable microquasar S26: A super-Eddington PeVatron

Context. S26 is an extragalactic microquasar with the most powerful jets ever discovered. They have a kinetic luminosity of Lj ∼ 5 × 1040 erg s−1. This implies that the accretion power to the black hole should be super-Eddington, of the order of Lacc ∼ Lj. However, the observed X-ray flux of this system indicates an apparent very sub-Eddington accretion luminosity of LX ≈ 1037 erg s−1. Aims. We aim to characterize the nature of S26, explain the system emission, and study the feasibility of super-Eddington microquasars as potential PeVatron sources. Methods. We first analyze multi-epoch X-ray observations of S26 obtained with XMM-Newton and model the super-Eddington disk and its wind. We then develop a jet model and study the particle acceleration and radiative processes that occur in shocks generated near the base of the jet and in its terminal region. Results. We find that the discrepancy between the jet and the apparent disk luminosities in S26 is caused by the complete absorption of the disk radiation by the wind ejected from the super-Eddington disk. The nonthermal X-rays are produced near the base of the jet, and the thermal X-rays are emitted in the terminal regions. The radio emission observed with the Australia Telescope Compact Array can be explained as synchrotron radiation produced at the reverse shock in the lobes. We also find that S26 can accelerate protons to PeV energies in both the inner jet and the lobes. The ultra-high energy protons accelerated in the lobes of S26 are injected into the interstellar medium with a total power of ∼1036 erg s−1. Conclusions. We conclude that S26 is a super-Eddington microquasar with a dense disk-driven wind that obscures the X-ray emission from the inner disk, and that the supercritical nature of the system allows the acceleration of cosmic rays to PeV energies.

Bow Shock Ripples and Their Modulation of Whistler Wave Packets: MMS Observations

AbstractThe terrestrial bow shock is the boundary where supersonic solar wind slows down abruptly near the magnetopause. The shock front geometry could be modulated by surface waves to form rippled structures, which impact the acceleration process of the solar wind particles. However, the rippled structures are hard to be identified unambiguously due to the similar signatures in single‐spacecraft observations between rippled shocks and reforming shocks. Here, we utilize the four‐spacecraft observations from the MMS mission to investigate an event of quasi‐perpendicular bow shock crossing. The periodic oscillations of shock normal directions and normal velocities support the scenario of surface wave propagation in the tangential direction. We also reconstruct the shock profile along the normal direction, and its monotonic shape further excludes the occurrence of shock reformation. These ripples are found to modulate the reflected ions and whistler wave packets, which adds to the complexity of the bow shock plasma environments.

Multifrequency polarimetry of high-synchrotron peaked blazars probes the shape of their jets

Multifrequency polarimetry is emerging as a powerful probe of blazar jets, especially with the advent of the Imaging X-ray Polarimetry Explorer (IXPE) space observatory. We studied the polarization of high-synchrotron peaked (HSP) blazars, for which both optical and X-ray emission can be attributed to synchrotron radiation from a population of nonthermal electrons. We adopted an axisymmetric stationary force-free jet model in which the electromagnetic fields are determined by the jet shape. When the jet is nearly parabolic, the X-ray polarization degree is ΠX ∼ 15–50%, and the optical polarization degree is ΠO ∼ 5–25%. The polarization degree is strongly chromatic: ΠX/ΠO ∼ 2–9. This chromaticity is due to the softening of the electron distribution at high energies, and is much stronger than for a uniform magnetic field. The electric vector position angle (EVPA) is aligned with the projection of the jet axis on the plane of the sky. These results compare very well with multifrequency polarimetric observations of HSP blazars. When the jet is instead nearly cylindrical, the polarization degree is large and weakly chromatic (we find ΠX ∼ 70% and ΠO ∼ 60%, close to the expected values for a uniform magnetic field). The EVPA is perpendicular to the projection of the jet axis on the plane of the sky. A cylindrical geometry is therefore practically ruled out by current observations. The polarization degree and the EVPA may be less sensitive to the specific particle acceleration process (e.g., magnetic reconnection or shocks) than previously thought.

Effect of temperature on the wave breaking amplitude of nonlinear relativistic strong plasma waves

AbstractThe wave‐breaking amplitude of a strong nonlinear plasma wave has been determined in a warm plasma. Pseudopotential technique has been adopted to describe the wave‐breaking phenomena in such plasma. The solution is obtained with the consideration that the phase velocity of the plasma wave is equal to the velocity of light in vacuum. This assumption is justified since the wave which is excited usually by the relativistic charged particle bunches or laser pulses has phase speed very close to the speed of light in vacuum. The investigation shows that the breaking amplitude decreases with the increase of the electron temperature. Furthermore, the wavelength of the plasma wave is seen to decrease as we increase the electron temperature. The obtained results have relevance in the astrophysical situation as well as in the field of plasma based charged particle acceleration process.

Three-stage Acceleration of Solar Energetic Particles Detected by Parker Solar Probe

Coronal mass ejections (CMEs) drive powerful shocks and thereby accelerate solar energetic particles (SEPs) as they propagate from the corona into interplanetary space. Here we present the processes of three-stage particle acceleration by a CME-driven shock detected by the in situ spacecraft—Parker Solar Probe (PSP) on 2022 August 27. The onset of SEPs is produced by a fast CME with a speed of 1284 km s−1 when it propagates to ∼2.85 R ⊙. The second stage of particle acceleration occurs when the fast CME catches up and interacts with a preceding slow one in interplanetary space at ∼40 R ⊙ (∼0.19 au). The CME interaction is accompanied by an intense interplanetary type II radio enhancement. Such direct measurement of particle acceleration during interplanetary CME interaction/radio enhancement is rarely recorded in previous studies. The third stage of energetic storm particles is associated with the CME-driven shock passage of the PSP at ∼0.38 au. Obviously, harder particle spectra are found in the latter two stages than the first one, which can arise from a stronger shock produced by the CME interaction and the enriched seed particles inside the preceding CME.

Energetic particle activity in AD Leo: Detection of a solar-like type-IV burst

Context. AD Leo is a young and active M dwarf with high flaring rates across the X-ray-to-radio bands. Flares accelerate particles in the outer coronal layers and often impact exo-space weather. Wide-band radio dynamic spectra let us explore the evolution of particle acceleration activity across the corona. Identifying the emission features and modelling the mechanisms can provide insights into the possible physical scenarios driving the particle acceleration processes.Aims. We performed an 8 h monitoring of AD Leo across the 550-850 MHz band using upgraded-Giant Metrewave Radio Telescope (uGMRT). The possible flare and post-flare emission mechanisms are explored based on the evolution of flux density and polarisation.Methods. The python-based module, Visibility Averaged Dynamic spectrum (VISAD), was developed to obtain the visibility-averaged wide-band dynamic spectra. Direct imaging was also performed with different frequency-time averaging. Based on existing observational results on AD Leo and on solar active region models, radial profiles of electron density and magnetic fields were derived. Applying these models, we explored the possible emission mechanisms and magnetic field profile of the flaring active region.Results. The star displayed a high brightness temperature (TB≈ 1010−1011K) throughout the observation. The emission was also nearly 100% left circularly polarised during bursts. The post-flare phase was characterised by a highly polarised (60–80%) solar-like type IV burst confined above 700 MHz.Conclusions. The flare emission favours a Z-mode or a higher harmonic X-mode electron cyclotron maser emission mechanism. The >700 MHz post-flare activity is consistent with a type-IV radio burst from flare-accelerated particles trapped in magnetic loops, which could be a coronal mass ejection (CME) signature. This is the first solar-like type-IV burst reported on a young active M dwarf belonging to a different age-related activity population (‘C’ branch) compared to the Sun (‘I’ branch). We also find that a multipole expansion model of the active region magnetic field better accounts for the observed radio emission than a solar-like active region profile.

Radio surface fluctuations in radio relics

Recent observations have revealed detailed structures of radio relics across a wide range of frequencies. In this work, we performed three-dimensional magnetohydrodynamical (3D MHD) simulations of merger shocks propagating through a turbulent magnetized intracluster medium. We then employed on-the-fly Lagrangian particles to explore the physical processes behind the origination of radio substructures and their appearance in high and low-frequency observations. We employed two cosmic-ray (CR) electron acceleration models, with a fresh injection of electrons from the thermal pool and the re-acceleration of mildly relativistic electrons. We used the relative surface brightness fluctuations, δSν, to define a “degree of patchiness.” First, we found that patchiness is produced if the shock’s surface has a distribution of Mach numbers, rather than a single Mach number. Second, radio relics appear patchier if the Mach number distribution consists of a large percentage of low Mach numbers (ℳ ≲ 2.5). Furthermore, as the frequency increases, the patchiness also becomes larger. Nevertheless, if radio relics are patchy at high frequencies (e.g., 18.6 GHz), they necessarily will also be patchy at low frequencies (e.g., 150 MHz). Then, to produce noticeable differences in the patchiness at low and high frequencies, the shock front should have a Mach number spread of σℳ ≳ 0.3 − 0.4. Finally, the extent of the patchiness depends on the Mach number distribution as well as the CR acceleration model. We propose δSν as a potential tool for extracting merger shock properties and information on particle acceleration processes at shocks in radio observations.

Properties of Electrons Accelerated by the Ganymede‐Magnetosphere Interaction: Survey of Juno High‐Latitude Observations

AbstractThe encounter between the Jovian co‐rotating plasma and Ganymede gives rise to electromagnetic waves that propagate along the magnetic field lines and accelerate particles by resonant or non‐resonant wave‐particle interaction. They ultimately precipitate into Jupiter's atmosphere and trigger auroral emissions. In this study, we use Juno/JADE, Juno/UVS data, and magnetic field line tracing to characterize the properties of electrons accelerated by the Ganymede‐magnetosphere interaction in the far‐field region. We show that the precipitating energy flux exhibits an exponential decay as a function of downtail distance from the moon, with an e‐folding value of 29°, consistent with previous UV observations from the Hubble Space Telescope (HST). We characterize the electron energy distributions and show that two distributions exist. Electrons creating the Main Alfvén Wing (MAW) spot and the auroral tail always have broadband distribution and a mean characteristic energy of 2.2 keV while in the region connected to the Transhemispheric Electron Beam (TEB) spot the electrons are distributed non‐monotonically, with a higher characteristic energy above 10 keV. Based on the observation of bidirectional electron beams, we suggest that Juno was located within the acceleration region during the 11 observations reported. We thus estimate that the acceleration region is extended, at least, between an altitude of 0.5 and 1.3 Jupiter radius above the 1‐bar surface. Finally, we estimate the size of the interaction region in the Ganymede orbital plane using far‐field measurements. These observations provide important insights for the study of particle acceleration processes involved in moon‐magnetosphere interactions.

Revealing the Variation Mechanism of ON 231 via the Two-component Shock-in-jet Model

The variation mechanism of blazars is a long-standing unresolved problem. In this work, we present a scenario to explain diverse variation phenomena for ON 231, where the jet emissions are composed of the flaring and the less variable components (most probably from the post-flaring blobs), and the variation is dominated by shock-in-jet instead of the Doppler effect. We perform correlation analysis for the multiwavelength light curves and find no significant correlations. For the optical band, ON 231 exhibits a harder when brighter (HWB) trend, and the trend seems to shift at different periods. Correspondingly, the correlation between the degree of polarization and flux exhibits a V-shaped behavior, and a similar translation relation during different periods is also found. These phenomena could be understood via the superposition of the flaring component and slowly varying background component. We also find that the slopes of the HWB trend become smaller at higher flux levels, which indicates the energy-dependent acceleration processes of the radiative particles. For the X-ray band, we discover a trend transition from HWB to softer when brighter (SWB) to HWB. We consider that the X-ray emission is composed of both the synchrotron tail and the synchrotron self-Compton components, which could be described by two log-parabolic functions. By varying the peak frequency, we reproduce the observed trend transition in a quantitative manner. For the γ-ray band, we find the SWB trend, which could be explained naturally if a very-high-energy γ-ray background component exists. Our study elucidates the variation mechanism of intermediate synchrotron-peaked BL Lac objects.

CFD-Based Study on the Flow and Kinetic Energy Characteristics of a Supercritical Suspended Abrasive Water Jet in the Deep-Sea Environment

A supercritical suspended abrasive water jet with dual inputs of pressure and heat is proposed to improve the cutting performance of the conventional suspended abrasive water jet in deep-sea environments. The paper studies the flow and kinetic characteristics of the supercritical suspended abrasive water jet. The CFD simulation method is proposed to investigate these characteristics by integrating a programmed database of supercritical water material properties with Ansys Fluent. The simulation and comparison show that abrasive particle density, abrasive particle size, inlet pressure, and water temperature affect the acceleration process of the abrasive particles. At the nozzle outlet, the velocity of the abrasive particles reaches over 95% of the supercritical water velocity. With the proposed supercritical abrasive water jet, the jet velocity is increased by 192.2% to 402.40 m/s compared to the conventional suspended abrasive water jet, reducing the amount of water used by 67.7% at a specified temperature of 773.15 K. Correspondingly, the medium kinetic energy is increased by 177.7% and the medium kinetic energy ratio is 2.78. The particle kinetic energy is increased by 723.2% and the particle kinetic energy ratio is 8.23.

The Slipping Magnetic Reconnection and Damped Quasiperiodic Pulsations in a Circular Ribbon Flare

The study of circular ribbon (CR) flares is important to understand the three-dimensional magnetic reconnection in the solar atmosphere. We investigate the slipping brightenings and damped quasiperiodic pulsations in a CR flare by multiwavelength observations. During the flaring process, two extreme ultraviolet brightenings (SP1 and SP2) slip synchronously along the ribbon in a counterclockwise direction. The ribbon and fans between them show synchronous enhancement with the microwave and hard X-ray (HXR) CR source. In the magnetohydrostatic extrapolation results and observations, the dome and outer spine display an evident counterclockwise twisting feature. We propose the slipping reconnection occurs between the fan and outer spine in the null point, which covers the region from SP1 to SP2. The fan of SP1 shows the strongest twist and produces the most efficient reconnection. The ribbon after SP1 becomes weak due to the destruction of the fan configuration. The fan of SP2 is in the front of the slipping motion, which initiates new reconnection and brightens the local ribbon. The twisting of the dome continuously promotes new reconnection in the null point, which brightens the ribbon in sequence to display a counterclockwise slipping feature. Thus, the twist of the dome may trigger and dominate the slipping reconnection, and the rotation of the central positive pole could be one possible cause of the twist. After the peak, the microwave and HXR emission shows damped oscillations at a period of 15 s. The collapse of the fan–spine structure may lead to the standing kink oscillations of the fan to modulate the reconnection and particle acceleration process.

Quasi-perpendicular shocks of galaxy clusters in hybrid kinetic simulations

Context. Cosmic ray acceleration in galaxy clusters is still an ongoing puzzle, with relativistic electrons forming radio relics at merger shocks and emitting synchrotron radiation. These shocks are also potential sources of ultra-high-energy cosmic rays, gamma rays, and neutrinos. Our recent work focuses on electron acceleration at low Mach number merger shocks in the hot intracluster medium which is characterized by high plasma beta. Using particle-in-cell (PIC) simulations, we previously showed that electrons are energized through the stochastic shock-drift acceleration process, which is facilitated by multi-scale turbulence, including ion-scale shock surface rippling. For the present work, we performed hybrid-kinetic simulations in a range of various quasi-perpendicular foreshock conditions, including plasma beta, magnetic obliquity, and the shock Mach number. Aims. We study the ion kinetic physics, which is responsible for the shock structure and wave turbulence, that in turn affects the particle acceleration processes. We cover the spatial and temporal scales, which allow the development of large-scale ion turbulence modes in the system. Methods. We applied a recently developed generalized fluid-particle hybrid numerical code that can combine fluid modeling for both electrons and ions with an arbitrary number of kinetic species. We limited this model to a standard hybrid simulation configuration with kinetic ions and fluid electrons. The model utilizes the exact form of the generalized Ohm’s law, allowing for an arbitrary choice of mass and energy densities, as well as the charge-to-mass ratio of the kinetic species. Results. We show that the properties of ion-driven multi-scale magnetic turbulence in merger shocks are in agreement with the ion structures observed in PIC simulations. In typical shocks with the sonic Mach number Ms = 3, the magnetic structures and shock front density ripples grow and saturate at wavelengths reaching approximately four ion Larmor radii. Only shocks with Ms ≳ 2.3 develop ripples. At very weak shocks with Ms ≲ 2.3, weak turbulence is formed downstream of the shock. We observed a moderate dependence of the strength of magnetic field fluctuations on the quasi-perpendicular magnetic field obliquity. However, as the field obliquity decreases, the shock front ripples exhibit longer wavelengths. Finally, we note that the steady-state structure of Ms = 3 shocks in high-beta plasmas shows evidence that there is little difference between 2D and 3D simulations. The turbulence near the shock front seems to be a 2D-like structure in 3D simulations.

Detecting non-thermal emission in a solar microflare using nested sampling

ABSTRACT Microflares are energetically smaller versions of solar flares, demonstrating the same processes of plasma heating and particle acceleration. However, it remains unclear down to what energy scales this impulsive energy release continues, which has implications for how the solar atmosphere is heated. The heating and particle acceleration in microflares can be studied through their X-ray emission, finding predominantly thermal emission at lower energies; however, at higher energies it can be difficult to distinguish whether the emission is due to hotter plasma and/or accelerated electrons. We present the first application of nested sampling to solar flare X-ray spectra, an approach that provides a quantitative degree of confidence for one model over another. We analyse Nuclear Spectroscopic Telescope Array X-ray observations of a small active region microflare (A0.02 GOES/XRS class equivalent) that occurred on 2021 November 17, with a new python package for spectral fitting, sunkit-spex, to compute the parameter posterior distributions and the evidence of different models representing the higher energy emission as due to thermal or non-thermal sources. Calculating the Bayes factor, we show that there is significantly stronger evidence for the higher energy microflare emission to be produced by non-thermal emission from flare-accelerated electrons than by an additional hot thermal source. Qualitative confirmation of this non-thermal source is provided by the lack of hotter (10 MK) emission in Solar Dynamic Observatory’s Atmospheric Imaging Assembly’s extreme ultraviolet data. The nested sampling approach used in this paper has provided clear support for non-thermal emission at the level of 3 × 1024 erg s−1 in this tiny microflare.

Determining the Origin of Very-high-energy Gamma Rays from Galactic Sources by Future Neutrino Observations

Recently, the Large High Altitude Air Shower Observatory (LHAASO) identified 12 gamma-ray sources emitting gamma rays with energies above 100 TeV, making them potential PeV cosmic-ray accelerators (PeVatrons). Neutrino observations are crucial in determining whether the gamma-ray radiation process is of hadronic or leptonic origin. In this paper, we study three detected sources, LHAASO J1908+0621, LHAASO J2018+3651, and LHAASO J2032+4102, which are also the most promising Galactic high-energy neutrino candidate sources with the lowest pretrial p-value based on the stacking searches testing for excess neutrino emission by the IceCube Neutrino Observatory. We study the lepto-hadronic scenario for the observed multiband spectra of these LHAASO sources considering the possible counterpart source of the LHAASO sources. The very-high-energy gamma rays are entirely attributed to the hadronic contribution; therefore, the most optimistic neutrino flux can be derived. Then, we evaluate the statistical significance (p-value) as a function of the observation time of IceCube and the next-generation IceCube-Gen2 neutrino observatory, respectively. Our results tend to disfavor that all gamma rays above 100 GeV from LHAASO J1908+0621 are of purely hadronic origin based on current IceCube observations, but the purely hadronic origin of gamma rays above 100 TeV is still possible. By IceCube-Gen2, the origin of gamma rays above 100 TeV from LHAASO J1908+0621 can be further determined at a 5σ significance level within a running time of ∼3 yr. For LHAASO J2018+3651 and LHAASO J2032+4102, the required running time of IceCube-Gen2 is ∼10 yr (3σ) and ∼10 yr (5σ), respectively. Future observations by the next-generation neutrino telescope will be crucial to understanding the particle acceleration and radiation processes inside the sources.

Probing the non-thermal physics of stellar bow shocks using radio observations

Context. Massive runaway stars produce bow shocks in the interstellar medium. Recent observations revealed radio emission from a few of these objects, but the origin of this radiation remains poorly understood. Aims. We aim to interpret this radio emission and assess under which conditions it could be either thermal (free–free) or non-thermal (synchrotron), and how to use the observational data to infer physical properties of the bow shocks. Methods. We used an extended non-thermal emission model for stellar bow shocks for which we incorporated a consistent calculation of the thermal emission from the forward shock. We fitted this model to the available radio data (spectral and intensity maps), including largely unexplored data at low frequencies. In addition, we used a simplified one-zone model to estimate the gamma-ray emission from particles escaping the bow shocks. Results. We can only explain the radio data from the best sampled systems (BD+43°3654 and BD+60°2522) assuming a hard electron energy distribution below ∼1 GeV, a high efficiency of conversion of (shocked) wind kinetic power into relativistic electrons (∼1 − 5%), and a relatively high magnetic-to-thermal pressure ratio of ηB ∼ 0.2. In the other systems, the interpretation of the observed flux density is more ambiguous, although a non-thermal scenario is also favoured. We also show how complementary observations at other frequencies can allow us to place stronger constraints in the model. We also estimated the gamma-ray fluxes from the HII regions around the bow shocks of BD+43°3654 and BD+60°2522, and obtained luminosities at GeV energies of ∼1033 erg s−1 and 1032 erg s−1, respectively, under reasonable assumptions. Conclusions. Stellar bow shocks can potentially be very efficient particle accelerators. This work provides multi-wavelength predictions of their emission and demonstrates the key role of low-frequency radio observations in unveiling particle acceleration processes. The prospects of detections with next-generation observatories such as SKA and ngVLA are very promising. Finally, BD+43°3654 may be detected in GeV in the near future, while bow shocks in general may turn out to be non-negligible sources of (at least leptonic) low-energy cosmic rays.

Unraveling parameter degeneracy in GRB data analysis

ABSTRACT Gamma-ray burst (GRB) afterglow light curves and spectra provide information about the density of the environment, the energy of the explosion, the properties of the particle acceleration process, and the structure of the decelerating jet. Due to the large number of parameters involved, the model can present a certain degree of parameter degeneracy. In this paper, we generated synthetic photometric data points using a standard GRB afterglow model and fit them using the Markov chain Monte Carlo (MCMC) method. This method has emerged as the preferred approach for analysing and interpreting data in astronomy. We show that, depending on the choice of priors, the parameter degeneracy can go unnoticed by the MCMC method. Furthermore, we apply the MCMC method to analyse the GRB 170817A afterglow. We find that there is a complete degeneracy between the energy of the explosion E, the density of the environment n, and the microphysical parameters describing the particle acceleration process (e.g. ϵe and ϵB), which cannot be determined by the afterglow light curve alone. Our results emphasize the importance of gaining a deep understanding of the degeneracy properties which can be present in GRB afterglows models, as well as the limitations of the MCMC method. In the case of GRB 170817, we get the following values for the physical parameters: E = 8 × 1050–1 × 1053 erg, n = 7 × 10−5–9 × 10−3, ϵe = 10−3–0.3, ϵB = 10−10–0.3.