Non-thermal Emission Research Articles

Impacts of nonthermal emission on the images of a black hole shadow and extended jets in two-temperature GRMHD simulations

The recent 230\,GHz observations from the Event Horizon Telescope (EHT) collaboration are able to image the innermost structure of the M87 galaxy showing the shadow of the black hole, a photon ring, and a ring-like structure that agrees with thermal synchrotron emission from the accretion disc. However, at lower frequencies, M87 is characterized by a large-scale jet with clear signatures of nonthermal emission. It is necessary to explore the impacts of nonthermal emission on black hole shadow images and extended jets, especially at lower frequencies. In this study, we aim to compare models with different electron heating prescriptions to one another and to investigate how these prescriptions and nonthermal electron distributions may affect black hole shadow images and the broadband spectrum energy distribution (SED) function. We performed general relativistic radiative transfer (GRRT) calculations in various two-temperature general relativistic magnetohydrodynamic (GRMHD) models utilizing different black hole spins and different electron heating prescriptions coupled with different electron distribution functions (eDFs). Through a comparison with GRRT images and SEDs, we found that when considering a variable kappa eDF, the parameterized prescription of the $R- electron temperature model with h $ = 1 is similar to the model with electron heating in the morphology of images, and the SEDs at a high frequency. This is consistent with previous studies using the thermal eDF. However, the nuance between them could be differentiated through the diffuse extended structure seen in GRRT images, especially at a lower frequency, and the behavior of SEDs at low frequency. The emission from the nearside jet region is enhanced in the case of electron heating provided by magnetic reconnection and it will increase if the contribution from the regions with stronger magnetization is included or if the magnetic energy contribution to kappa eDF mainly in the magnetized regions is considered. Compared with the thermal eDF, the peaks of the SEDs shift to a lower frequency when we consider nonthermal eDF.

Abell 746: A Highly Disturbed Cluster Undergoing Multiple Mergers

We present deep XMM-Newton, Karl G. Jansky Very Large Array, and upgraded Giant Metrewave Radio Telescope observations of Abell 746, a cluster that hosts a plethora of diffuse emission sources that provide evidence for the acceleration of relativistic particles. Our new XMM-Newton images reveal a complex morphology of the thermal gas with several substructures. We observe an asymmetric temperature distribution across the cluster: the southern regions exhibit higher temperatures, reaching ∼9 keV, while the northern regions have lower temperatures (≤4 keV), likely due to a complex merger. We find evidence of three surface brightness edges and one candidate edge, of which three are merger-driven shock fronts. Combining our new data with published LOw-Frequency ARray observations has unveiled the nature of diffuse sources in this system. The bright NW relic shows thin filaments and a high degree of polarization with aligned magnetic field vectors. We detect a density jump, aligned with the fainter relic to the north. To the south, we detect high-temperature regions, consistent with the shock-heated regions and a density jump coincident with the northern tip of the southern radio structure. Its integrated spectrum shows a high-frequency steepening. Lastly, we find that the cluster hosts large-scale radio halo emission. A comparison of the thermal and nonthermal emission reveals an anticorrelation between the bright radio and X-ray features at the center. Our findings suggest that Abell 746 is a complex system that involves multiple mergers.

A comprehensive X-ray analysis of the massive O-type binary HD 93250 over two decades

Massive star winds are known to be responsible for X-ray emission arising from wind plasma heated by the strong shocks up to temperatures of 10$^6$--10$^7$ K in the case of colliding wind binaries. The investigation of X-ray emission from massive stars thus constitutes a valuable tool for identifying binaries, which is otherwise a difficult task using classical techniques. We investigated thermal and nonthermal X-ray emission from the massive O-type star HD\,93250 to unveil its binary orbital parameters independently. To meet our goal, we analyzed X-ray data obtained with European Photon Imaging Camera on board XMM-Newton spanning over sim 19 years. Additionally, we analyzed NuSTAR observations of HD\,93250 taken at various epochs. We determine the variability timescale of the X-ray emission to be 193.8pm 1.3\,d, in full agreement with the 194.3pm 0.4\,d period derived from the astrometric orbit. The X-ray spectrum of HD\,93250 is well explained by a three-temperature thermal plasma emission model with temperatures of 0.26, 1.0, and 3.3 keV. The resulting X-ray flux varies in compliance with the typical colliding wind emission from eccentric massive binaries where it enhances near periastron passage and decreases gradually close to apastron in proportion to the inverse of the binary separation. The periastron-to-apastron X-ray emission ratio points to an eccentricity range of 0.20-0.25, once again in agreement with the previously determined astrometric orbit. Finally, we do not detect any hard X-ray emission attributable to nonthermal emission above 10 keV. Given the derived plasma temperature, the strong phase-locked variability, and the significant over-luminosity in X-rays, we establish that the X-ray emission from HD\,93250 is dominated by the colliding-wind region. Our results lend support to the idea that X-ray time analysis of massive stars constitutes a relevant tool for investigating their multiplicity and for extracting relevant information on their basic orbital parameters ---such as period and eccentricity--- independently of any orbital solution derived from classical techniques.

Thermal X-ray emission in the western half of the LMC superbubble 30 Dor C

ABSTRACT While 30 Dor C is a unique superbubble in the Large Magellanic Cloud for its luminous non-thermal X-ray emission, the thermal X-ray emission it emanates has not yet been thoroughly investigated and well constrained. Based on the separate ∼1 Ms deep XMM–Newton and Chandra observations, we report the discovery of the thermally emitting plasma in some portions of the western half of 30 Dor C. The thermal emission can be reproduced by a collisional-ionization-equilibrium plasma model with an average electron temperature of ∼0.4 keV. We find a significant overabundance of the intermediate-mass elements such as O, Ne, Mg, and Si, which may be indicative of a recent supernova explosion in 30 Dor C. Dynamical properties in combination with the information of the OB association LH 90 suggest that the internal post-main-sequence stars dominate the power of the superbubble and blow it out in the past ∼1 Myr.

CONGRuENTS (COsmic-ray, Neutrino, Gamma-ray and Radio Non-Thermal Spectra) – II. Population-level correlations between galactic infrared, radio, and γ-ray emission

ABSTRACT Galaxies obey a number of empirical correlations between their radio, γ-ray, and infrared emission, but the physical origins of these correlations remain uncertain. Here, we use the CONGRuENTS model for broad-band non-thermal emission from star-forming galaxies, which self-consistently calculates energy-dependent transport and non-thermal emission from cosmic ray hadrons and leptons, to predict radio and γ-ray emission for a synthetic galaxy population with properties drawn from a large deep-field survey. We show that our synthetic galaxies reproduce observed relations such as the far infrared (FIR)–radio correlation, the FIR–γ correlation, and the distribution of radio spectral indices, and we use the model to explain the physical origins of these relations. Our results show that the FIR–radio correlation arises because the amount of cosmic ray electron power ultimately radiated as synchrotron emission varies only weakly with galaxy star formation rate as a result of the constraints imposed on gas properties by hydrostatic balance and turbulent dynamo action; the same physics dictates the extent of proton calorimetry in different galaxies, and thus sets the FIR–γ–ray correlation. We further show that galactic radio spectral indices result primarily from competition between thermal free–free emission and energy-dependent loss of cosmic ray electrons to bremsstrahlung and escape into galactic haloes, with shaping of the spectrum by inverse Compton, synchrotron, and ionization processes typically playing a subdominant role. In addition to explaining existing observations, we use our analysis to predict a heretofore unseen correlation between the curvature of galaxies’ radio spectra and their pion-driven γ-ray emission, a prediction that will be testable with upcoming facilities.

Indication of a fast ejecta fragment in the atomic cloud interacting with the southwestern limb of SN 1006

Context. Supernova remnants interacting with molecular and atomic clouds are interesting X-ray sources for studies of broadband nonthermal emission. X-ray line emission in these systems can be produced by different processes, such as low-energy cosmic rays (LECRs) interacting with the cloud and fast ejecta fragments moving in the cloud. Aims. This paper is aimed at studying the origin of the non-thermal X-ray emission of the southwestern limb of SN 1006 beyond the main shock to determine whether the emission is due to LECRs diffusing in the cloud or to ejecta knots moving into the cloud. Methods. We analyzed the X-ray emission of the southwestern limb of SN 1006, where the remnant interacts with an atomic cloud, using three different X-ray telescopes: NuSTAR, Chandra, and XMM-Newton. We also performed a combined spectro-imaging analysis of this region. Results. Our analysis of the nonthermal X-ray emission of the southwestern limb of SN 1006 interacting with an atomic cloud has led to the detection of an extended X-ray source in the atomic cloud, approximately 2 pc upstream of the shock front. The source is characterized by a hard continuum (described by a power law with photon index Γ ∼ 1.4) and by Ne, Si, and Fe emission lines. The observed flux suggests that the origin of the X-ray emission is not associated with LECRs interacting with the cloud. On the other hand, the spectral properties of the source, together with the detection of an IR counterpart visible with Spitzer-MIPS at 24 μm, are in good agreement with the general expectations for a fast ejecta fragment moving within the atomic cloud. Conclusions. We detected X-ray and IR emission from a possible ejecta fragment, with an approximate radius of 1 × 1017 cm and approximate mass of 10−3 M⊙ at about 2 pc out of the shell of SN 1006, in the interaction region between the southwestern limb of the remnant and the atomic cloud.

Relative Yield of Thermal and Nonthermal Emission during Weak Flares Observed by STIX during 2021 September 20–25

The disparate nature of the thermal–nonthermal energy partition during flares, particularly during weak flares, is still an open issue. Following the Neupert effect, quantifying the relative yield of X-ray emission in different energy bands can enable the inferring of the underlying energy release mechanism. During 2021 September 20–25, the Solar Orbiter mission—being closer to the Sun (∼0.6 au) and having a moderate separation angle (<40°) from the Sun–Earth line—offered a unique opportunity to analyze multiwavelength emission from ∼200 (mostly weak) flares, commonly observed by the Spectrometer Telescope for Imaging X-rays (STIX), STEREO-A, GOES, and the Solar Dynamics Observatory. Associating the quotient (q f ) of hard X-ray fluence (12–20 keV) and soft X-ray flux (4–10 keV) with the peak soft X-ray flux enabled us to identify strongly nonthermal flares. A multiwavelength investigation of spectral and imaging-mode observations of the 20 strongly nonthermal weak flares reveals an inverse relationship of q f with the emission measure (and density), and a positive relationship with the flare plasma temperature. This indicates that the plasma in tenuous loops attains higher temperatures compared to that in the denser loops, in response to nonthermal energy deposition. This is in agreement with the plasma parameters of the coronal loops, as derived by applying the one-dimensional Palermo–Harvard hydrodynamical code to coronal loop plasma with different initial coronal loop base pressures when subjected to similar heating inputs. Our investigation, therefore, indicates that the plasma parameters of the flaring loop in the initial phase have a decisive role in the thermal–nonthermal energy partitioning.

Cosmic-ray induced ionization rates and non-thermal emissions from nuclei of starburst galaxies

ABSTRACT Cosmic rays are the only agent capable of ionizing the interior of dense molecular clouds and, thus, they are believed to play an essential role in determining the physical and chemical evolution of star-forming regions. In this work, we aim to study cosmic-ray induced ionization rates in starburst environments using non-thermal emissions of cosmic rays from starburst nuclei. To this end, we first revisit cosmic-ray models, which could explain data of non-thermal emissions from radio to X-ray and gamma-ray from nuclei of three prototypical starburst galaxies NGC 253, M82, and Arp 220. These models are then applied to predict ionization rates in starburst environments, which gives values around 10−14 s−1. Such a high value of the ionization rate, which is 2 to 3 orders of magnitude higher than the typical values found in the Milky Way, is probably due to relatively high rates of supernova explosions occurring within the nuclei of these starburst galaxies. We also discuss in more detail the case of NGC 253, where our predicted ionization rate is found to be, in most cases, a few times smaller than the values inferred from molecular line observations of clouds in the starburst nucleus. The general framework provided in this work illustrates how the use of non-thermal emission data could help to provide more insights into ionization rates or, more generally, cosmic-ray impact in starburst environments.

Efficiency of solar microflares in accelerating electrons when rooted in a sunspot

Context. The spectral shape of the X-ray emission in solar flares varies with the event size, with small flares generally exhibiting softer spectra than large events, indicative of a relatively lower number of accelerated electrons at higher energies. Aims. We investigate two microflares of GOES classes A9 and C1 (after background subtraction) observed by STIX onboard Solar Orbiter with exceptionally strong nonthermal emission. We complement the hard X-ray imaging and spectral analysis by STIX with co-temporal observations in the (E)UV and visual range by AIA and HMI to investigate what makes these microflares so efficient in high-energy particle acceleration. Methods. We made a preselection of events in the STIX flare catalog based on the ratio of the thermal to nonthermal quicklook X-ray emission. The STIX spectrogram science data were used to perform spectral fitting to identify the non-thermal and thermal components. The STIX X-ray images were reconstructed to analyze the spatial distribution of the precipitating electrons and the hard X-ray emission they produce. The EUV images from SDO/AIA and SDO/HMI LOS magnetograms were analyzed to better understand the magnetic environment and the chromospheric and coronal response. For the A9 event, EOVSA microwave observations were available, allowing for image reconstruction in the radio domain. Results. We performed case studies of two microflares observed by STIX on October 11, 2021 and November 10, 2022, which showed unusually hard microflare X-ray spectra with power-law indices of the electron flux distributions of δ = (2.98 ± 0.25) and δ = (4.08 ± 0.23), during their non-thermal peaks and photon energies up to 76 keV and 50 keV, respectively. For both events under study, we found that one footpoint is located within a sunspot covering areas with mean magnetic flux densities in excess of 1500 G, suggesting that the hard electron spectra are caused by the strong magnetic fields the flare loops are rooted in. Additionally, we revisited a previously published unusually hard RHESSI microflare and found that in this event, there was also one flare kernel located within a sunspot, which corroborates the result from the two hard STIX microflares under study in this work. Conclusions. The characteristics of the strong photospheric magnetic fields inside the sunspot umbrae and penumbrae where flare loops are rooted play an important role in the generation of exceptionally hard X-ray spectra in these microflares.

Application of neural networks to synchro-Compton blazar emission models

Context. Jets from supermassive black holes at the centers of active galaxies are the most powerful and persistent sources of electromagnetic radiation in the Universe. To infer the physical conditions in the otherwise out-of-reach regions of extragalactic jets, we usually rely on fitting their spectral energy distributions (SEDs). The calculation of radiative models for the jet’s non-thermal emission usually relies on numerical solvers of coupled partial differential equations. Aims. In this work, we use machine learning to tackle the problem of high computational complexity to significantly reduce the SED model evaluation time, which is necessary for SED fittings carried out with Bayesian inference methods. Methods. We computed the SEDs based on the synchrotron self-Compton model for blazar emission using the radiation code ATHEvA. We used them to train neural networks (NNs) to explore whether they can replace the original code, which is computationally expensive. Results. We find that a NN with gated recurrent unit neurons (GRUN) can effectively replace the ATHEvA leptonic code for this application, while it can be efficiently coupled with Markov chain Monte Carlo (MCMC) and nested sampling algorithms for fitting purposes. We demonstrate this approach through an application to simulated data sets, as well as a subsequent application to observational data. Conclusions. We present a proof-of-concept application of NNs to blazar science as the first step in a list of future applications involving hadronic processes and even larger parameter spaces. We offer this tool to the community through a public repository.

Owens Valley Radio Observatory monitoring of LS I +61°303 completes three cycles of the super-orbital modulation

Context. The high-mass X-ray binary LS I +61°303 is composed of a Be-type star and a compact object in an eccentric orbit. The emission from this source is variable and periodic across the electromagnetic spectrum, from radio to very high-energy γ rays. The orbital period has been determined as P1 ≈ 26.5 d, and the source also features a super-orbital period with a value of Plong ≈ 4.6 years. Long-term monitoring of the binary by the Owens Valley Radio Observatory (OVRO) at 15 GHz has now completed 13.8 years, which corresponds to three full cycles of the super-orbital period. This is exactly one long-term cycle more than in the previous publication about OVRO observations of this source. Aims. Our aim is to investigate the presence and the stability of periodic signals in the radio data and to test if they are in agreement with previous results. This will contribute to the understanding of the physical processes behind the non-thermal emission from this source. Methods. We performed a timing analysis of the OVRO radio light curve and made use of the generalized Lomb-Scargle periodogram. We also combined the OVRO data with the full archive of previous radio observations and computed the discrete autocorrelation function. Results. The most powerful features in the periodogram of the OVRO data are two peaks at P1 = 26.49 ± 0.05 d and P2 = 26.93 ± 0.05 d, which are well separated from each other and clearly stand out above the very low noise level. The previously detected long-term period is still present in these new radio data, and our measurement is Plong = 1698 ± 196 d. Dividing the OVRO data into three segments of equal length showed that the two periods, P1 and P2, are present in the periodogram of each of the consecutive long-term cycles. Our analysis of the full radio archive resulted in the detection of the same three periods, and the autocorrelation function showed a regular pattern, proving the continuity of the decades-spanning stability of the super-orbital modulation. In addition, we report a possible systematic modulation of the radio flux density with a timescale of approximately 40 years that has so far remained unnoticed. Conclusions. The physical model of a relativistic jet whose mass loading is modulated with the orbital period P1 and is precessing with the slightly larger period P2, giving rise to a beating with period Plong, had previously been able to reproduce the radio and gigaelectron volt emission from this source. The ongoing presence and the stability of the periodic signals imply that this model is still the most plausible explanation for the physical processes at work in this source.

Evidence for flare-accelerated particles in large scale loops in the behind-the-limb gamma-ray solar flare of September 29, 2022

We report on the detection of the gamma-ray emission above 100 MeV from the solar flare of September 29, 2022, by Fermi LAT with simultaneous coverage in HXR by Solar Orbiter STIX. The Solar Orbiter-Earth separation was 178° at the time of the flare as seen from Earth, with Solar Orbiter observing the east limb. Based on STIX imaging, the flare was located 16° behind the eastern limb as seen from Earth. The STIX and GBM non-thermal emission and the LAT emission above 100 MeV all show similarly shaped time profiles, and the Fermi profiles peaked only 20 s after the STIX signal from the main flare site, setting this flare apart from all the other occulted flares observed by Fermi LAT. The radio spectral imaging based on the Nançay Radioheliograph and ORFEES spectrograph reveal geometries consistent with a magnetic structure that connects the parent active region behind the limb to the visible disk. We studied the basic characteristics of the gamma-ray time profile, in particular, the rise and decay times and the time delay between the gamma-ray and HXR peak fluxes. We compared the characteristics of this event with those of four Fermi LAT behind-the-limb flares and with an on-disk event and found that this event is strikingly similar to the impulsive on-disk flare. Based on multiwavelength observations, we find that the gamma-ray emission above 100 MeV originated from ions accelerated in the parent active region behind the limb and was transported to the visible disk via a large magnetic structure connected to the parent active region behind the limb. Our results strongly suggest that the source of the emission above 100 MeV from the September 29, 2022 flare cannot be the CME-driven shock.

Time evolution of X-ray bright points observed with NuSTAR

ABSTRACT We present analysis on two X-ray bright points observed over several hours during the recent solar minimum (2020 February 21 and 2020 September 12–13) with the Nuclear Spectroscopic Telescope Array (NuSTAR), a sensitive hard X-ray imaging spectrometer. This is so far the most detailed study of bright points in hard X-rays, emission which can be used to search for faint hot and/or non-thermal sources. We investigate the bright points’ time evolution with NuSTAR, and in extreme ultraviolet (EUV) and soft X-rays with Solar Dynamic Observatory/Atmospheric Imaging Assembly (SDO/AIA) and Hinode/X-Ray Telescope. The variability in the X-ray and EUV time profiles is generally not well matched, with NuSTAR detecting spikes that do not appear in EUV. We find that, for the 2020 February bright point, the increased X-ray emission during these spikes is due to material heated to ∼ 4.2–4.4 MK (found from fitting the X-ray spectrum). The 2020 September bright point also shows spikes in the NuSTAR data with no corresponding EUV signature seen by SDO/AIA, though in this case, it was due to an increase in emission measure of material at ∼ 2.6 MK and not a significant temperature change. So, in both cases, the discrepancy is likely due to the different temperature sensitivity of the instruments, with the X-ray variability difficult to detect in EUV due to cooler ambient bright point emission dominating. No non-thermal emission is detected, so we determine upper limits finding that only a steep non-thermal component between 3 and 4 keV could provide the required heating whilst being consistent with a null detection in NuSTAR.

The Non-Thermal Radio Emissions of the Solar Transition Region and the Proposal of an Observational Regime

The transition region is a very thin but most peculiar layer in the solar atmosphere located between the solar chromosphere and the corona. It is a key region for understanding coronal heating, solar eruption triggers, and the origin of solar winds. Here, almost all physical parameters (density, temperature, and magnetic fields) have the maximum gradient. Therefore, this region should be highly dynamic, including fast energy releasing and transporting, plasma heating, and particle accelerating. The physical processes can be categorized into two classes: thermal and non-thermal processes. Thermal processes can be observed at ultraviolet (UV) and extreme ultraviolet (EUV) wavelengths via multi-wavelength images. Non-thermal processes accelerate non-thermal electrons and produce radio emissions via the gyrosynchrotron mechanism resulting from the interaction between the non-thermal electrons and magnetic fields. The frequency range spans from several GHz to beyond 100 GHz, in great number of bursts with narrowband, millisecond lifetime, rapid frequency drifting rates, and being referred to as transition region small-scale microwave bursts (TR-SMBs). This work proposes a new type of Solar Ultra-wide Broadband Millimeter-wave Spectrometer (SUBMS) that can be used to observe TR-SMBs. From SUBMS observations, we can derive rich dynamic information about the transition region, such as information about non-thermal energy release and propagation, the flows of plasma and energetic particles, the magnetic fields and their variations, the generation and transportation of various waves, and the formation and evolution of the source regions of solar eruptions. Such an instrument can actually detect the non-thermal signals in the transition region during no flare as well as the eruptive high-energy processes during solar flares.

Discovery of an old supernova remnant candidate through carbon monoxide line emission

Most old supernova remnants (SNRs) in the Milky Way have not yet been identified. Considering their high potential number and the sufficient momentum-energy transfer to the interstellar medium (ISM), they are a key part of our understanding of the overall role of SNRs in the ISM. Here, we report our discovery of an expanding molecular shell identified by CO line observations, namely It covers a known SNR, specifically G16.0$-$0.5, and is larger in size (i.e.,\ 0.56$^ over 0.20$^ Based on its spatial and kinematic structures, weak nonthermal radio-continuum emission, and derived physical properties, we suggest that it is an old SNR. At a systemic velocity of +41.3 the best estimated kinematic distance of is sim 3.2 kpc, implying its radius of about 15.6 pc. The age of is estimated to be greater than sim $ yr, and, in a dense molecular environment, it has formed dense and thin shell layers. The kinetic energy of the expanding molecular gas of is about $6.4 $ erg, accounting for approximately six percent of the initial SN explosion energy. Although old SNRs have essentially become cold and hard to detect, our discovery suggests that they can be found by searching for CO line emissions.

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.

Cosmic ray acceleration and non-thermal emission from fast luminous optical transient sources

Fast blue optical transients (FBOTs) represent a new class of highly energetic sources observed from radio to X-rays. High luminosity, light curves and spectra of the sources can be understood if they are associated with supernova-like or tidal disruption events. Radio observations of the transient sources revealed a mildly relativistic expansion of some of the remnants. The high power and mildly relativistic shock velocities are providing favorable conditions for very high energy particle acceleration. In this paper we present a model of particle acceleration in mildly relativistic magnetohydrodynamic (MHD) outflow of the transient source. To construct the non-thermal radiation and cosmic ray spectra in a broad range of energies we combined the microscopic particle-in-cell simulations of electron and proton injection at mildly relativistic shock with Monte Carlo technique for high energy particle transport and acceleration. The model allows describing the observed non-thermal radio and X-ray emission of CSS161010, and predicts that such sources can accelerate cosmic rays to energies above 10 PeV, and possible upper limit of maximum energy is 100 PeV. With the expected event rate of FBOTs they can contribute to the very high energy cosmic rays population in galaxies.

Probing diffuse radio emission in bridges between galaxy clusters with uGMRT

Aims. Recent X-ray and Sunyaev–Zeldovich (SZ) observations have detected thermal emission between early-stage merging galaxy clusters. The main purpose of this work is to investigate the properties of the non-thermal emission in the interacting cluster pairs Abell 0399-Abell 0401 and Abell 21-PSZ2 G114.9. Methods. These two unique cluster pairs have been found in an interacting state. In both cases, their connection along a filament is supported by an SZ effect detected by the Planck satellite and, in the special case of Abell 0399-Abell 0401, the presence of a radio bridge has been already confirmed by LOFAR observations at 140 MHz. Here, we analyse new high-sensitivity, wideband (250–500 MHz) uGMRT data of these two systems and describe an injection procedure to place limits on the spectrum of Abell 0399-Abell 0401 and on the radio emission between Abell 21-PSZ2 G114.9. Results. In both cases, the low-surface-brightness diffuse emission is not detected in Band 3 (250–500 MHz). For the A399-A401 pair, we are able to constrain the steep spectral index of the bridge emission to be α &gt; 2.2 with a 95% confidence level between 140 MHz and 400 MHz. We also detect a small patch of the bridge with a flatter spectral index, which may suggest a variable spectral index distribution across the bridge area. For the A21-PSZ2 G114.9 pair, we are able to place an upper limit on the flux density of the bridge emission with two different methods, finding at the central frequency of 383 MHz a conservative value of fu1 &lt; 260 mJy at a 95% confidence level, and a lower value of fu2 &lt; 125 mJy at an 80% confidence level, based on visual inspection and a morphological criterion. Conclusions. Our work provides a constraint on the spectrum in the bridge A399-A401 that disfavours shock acceleration as the main mechanism for the radio emission. The methods that we propose for the limits on the radio emission in the A21-PSZ2 G114.9 system represent a first step towards a systematic study of these sources.

EROSITA studies of the Carina nebula

Context. During the first four all-sky surveys eRASS:4, which was carried out from December 2019 to 2021, the extended Roentgen Survey with an Imaging Telescope Array (eROSITA) on board the Spektrum-Roentgen-Gamma (Spektr-RG, SRG) observed the Galactic H II region, the Carina nebula. Aims. We analysed the eRASS:4 data to study the distribution and spectral properties of the hot interstellar plasma and the bright stellar sources in the Carina nebula. Methods. The spectral extraction regions of the diffuse emission were defined based on the X-ray spectral morphology and multi-wavelength data. The spectra were fit with a combination of thermal and non-thermal emission models. The X-ray bright point sources in the Carina nebula are the colliding wind binary η Car, several O stars, and Wolf–Rayet (WR) stars. We extracted the spectra of the brightest stellar sources, which can be well fit with a multi-component thermal plasma model. Results. The spectra of the diffuse emission in the brighter parts of the Carina nebula are well reproduced by two thermal models, a lower-temperature component (~0.2 keV) and a higher-temperature component (0.6–0.8 keV). An additional non-thermal component dominates the emission above ~1 keV in the Central region around η Car and the other massive stars. Significant orbital variation in the X-ray flux was measured for η Car, WR 22, and WR 25. η Car requires an additional time-variable thermal component in the spectral model, which is associated with the wind-wind collision zone. Conclusions. Properties such as temperature, pressure, and luminosity of the X-ray emitting plasma in the Carina nebula derived from the eROSITA data are consistent with theoretical calculations of emission from superbubbles. This confirms that the X-ray emission is caused by the hot plasma inside the Carina nebula that has been shocked-heated by the stellar winds of the massive stars, in particular, of η Car.

ZWCL 1856.8: A Rare Double Radio Relic System Captured within NuSTAR and Chandra Field of View

Observations of galaxy cluster mergers provide insights on the particle acceleration and heating mechanisms taking place within the intracluster medium. Mergers form shocks that propagate through the plasma, which result in shock/cold fronts in the X-ray, and radio halos and/or relics in the radio regime. The connection between these tracers and the mechanisms driving non-thermal processes, such as inverse Compton, are not well understood. ZWCL 1856.8 is one of the few known double radio relic systems that originate from nearly head-on collisions observed close to the plane of the sky. For the first time, we study NuSTAR and Chandra observations of this system, which contains both relics within their field of view. The spectro-imaging analyses results of the system suggest weak shock fronts with M numbers within 2σ of the radio-derived values, and provide evidence of inverse Compton emission at both relic sites. Our findings have great uncertainties due to the shallow exposure times available. Deeper NuSTAR and Chandra data are crucial for studying the connection of the radio and X-ray emission features and for constraining the thermal versus non-thermal emission contributions in this system. We also present methods and approaches on how to investigate X-ray properties of double relic systems by taking full advantage of the complementary properties of NuSTAR and Chandra missions.