Discovery of a Pulsar Wind Nebula Candidate Associated with the Galactic PeVatron 1LHAASO J0343+5254u

The Fermi orbiter, mapping the entire gamma-ray sky every three hours, monitors the cosmos for high-energy phenomena both fleeting and enduring.

Context.Vela X is a pulsar wind nebula in which two relativistic particle populations with distinct spatial and spectral distributions dominate the emission at different wavelengths. An extended 2° × 3° nebula is seen in radio and GeV gamma rays. An elongated cocoon prevails in X-rays and TeV gamma rays.Aims.We use ~9.5 yr of data from theFermiLarge Area Telescope (LAT) to disentangle gamma-ray emission from the two components in the energy range from 10 GeV to 2 TeV, bridging the gap between previous measurements at GeV and TeV energies.Methods.We determine the morphology of emission associated to Vela X separately at energies <100 and >100 GeV, and compare it to the morphology seen at other wavelengths. Then, we derive the spectral energy distribution of the two gamma-ray components over the full energy range.Results.The best overall fit to the LAT data is provided by the combination of the two components derived at energies <100 and >100 GeV. The first component has a soft spectrum, spectral index 2.19 ± 0.16−0.22+0.05, and extends over a region of radius 1.°36±0.°04, consistent with the size of the radio nebula. The second component has a harder spectrum, spectral index 0.9 ± 0.3−0.1+0.3, and is concentrated over an area of radius 0.°63±0.°03, coincident with the X-ray cocoon that had already been established as accounting for the bulk of the emission at TeV energies.Conclusions.The spectrum measured for the low-energy component corroborates previous evidence for a roll-over of the electron spectrum in the extended radio nebula at energies of a few tens of GeV possibly due to diffusive escape. The high-energy component has a very hard spectrum: if the emission is produced by electrons with a power-law spectrum, the electrons must be uncooled, and there is a hint that their spectrum may be harder than predictions by standard models ofFermiacceleration at relativistic shocks.

When a star dies, it leaves a mark on its surrounding environment. The energy from the supernova explosion forms an expanding shock wave that interacts with interstellar and circumstellar material, creating what we know as a supernova remnant (SNR). If the original star has a mass that is greater than or equal to 8 solar masses, this can also lead to the formation of a rapidly rotating neutron star called a pulsar. As these objects evolve, they interact with the surrounding environment, producing non-thermal and thermal emission. For an SNR, its non-thermal emission arises from a population of relativistic particles being accelerated at the shock front of the SNR, while its thermal emission arises from the shock front heating ejecta and and swept-up interstellar medium to X-ray emitting temperatures. For pulsars, their non-thermal emission arises from relativistic particles being accelerated at the termination shock of a pulsar wind. These particles interact with surrounding magnetic fields and ambient photon fields producing synchrotron and inverse Compton emission which we observe as a pulsar wind nebula (PWN), while its thermal emission arises from the surface of the neutron star. These properties of SNRs and pulsars provide a unique window into studying the acceleration, injection, propagation and interaction of highly energetic particles called cosmic rays with the interstellar medium. In addition, they providing information about the evolution, and dynamics of these objects; properties of the shock fronts; details about the original progenitor star; and the impact that these objects have on their surroundings. The research presented here focuses on analysing the intimate connection between cosmic rays, the non-thermal emission arising from SNRs interacting with molecular clouds, and pulsar wind nebulae; as well as analysing the observational and evolutionary properties of these objects. In this thesis we model the propagation of cosmic rays through the Galaxy in an attempt to characterise a standard cosmic ray background with uncertainties, to reveal the origin of the cosmic ray electron positron anomaly. Furthermore, we analyse the gamma-ray emission from SNRs Kes 79 and MSH 11-61A, which are known to be interacting with molecular clouds, as well as the non-thermal X-ray emission arising from the PWN of PSR J1741-2054. We find that the emission from both SNRs most likely arises from the decay of neutral pions that resulted from the interaction of relativistic ions which are accelerated at the shock-front of a SNR, with ambient material. For PSR J1741-2054, we characterise the properties, minimum magnetic field and minimum energy of the particle population that produces the observed diffuse synchrotron emission that surrounds and trails the pulsar. In addition, we characterise the X-ray emission arising from Kes 79, MSH 11-61A and PSR J1741-2054, in an attempt to shed light on the origin and nature of these objects and their emission. Using X-ray data from XMM-Newton and Suzaku respectively, we probe the temperature, ionisation state, and elemental abundance of the shocked gas of each SNR. This allows us to determine their evolutionary properties, properties of the shock, and mass of the original progenitor; and constrain the density of the X-ray emitting plasma. Using Chandra, we determined the temperature of PSR J1741-2054, as well as characterised its proper motion, velocity, direction of motion, and presence of small scale structure immediately surrounding the pulsar.

The Tibet ASγ experiment has been observing cosmic gamma rays and cosmic rays in the energy range from teraelectron volts to several tens of petaelectron volts with a surface detector array since 1990. The derivation of cosmic gamma-ray flux is made by finding the excess distribution of the arrival direction of air showers above background cosmic rays. In 2014, the underground water Cherenkov muon detector (MD) was added to separate cosmic gamma rays from the background on the basis of the muon-less feature of the air showers of gamma-ray origin; hybrid observations using these two detectors were started at this time. In the present study, we developed methods to separate gamma-ray-induced air showers and hadronic cosmic-ray-induced ones using the measured particle number density distribution to improve the sensitivity of cosmic gamma-ray measurements using the Tibet air shower array data alone before the installation of the MD. We tested two approaches based on neural networks. The first method used feature values representing the lateral spread of the secondary particles, and the second method used the shower image data. To compare the separation performance of each method, we analyzed Monte Carlo air shower events in the vertically incident direction with mono-initial-energy gamma rays and protons. When discriminated by a single feature, the feature with the highest separation performance has an area under the curve (AUC) value of 0.701 for a gamma-ray energy of 10 TeV and 0.808 for 100 TeV. A separation method with a multilayer perceptron (MLP) based on multiple features has AUC values of 0.761 for a gamma-ray energy of 10 TeV and 0.854 for 100 TeV, which represents an improvement of approximately 5% in the AUC value compared with the single-feature case. We also found that the feature values that effectively contribute to the separation vary depending on the energy. A separation method with a convolutional neural network (CNN) using the shower image data has AUC values of 0.781 for a gamma-ray energy of 10 TeV and 0.901 for 100 TeV, which are approximately 5% higher than those of the MLP method. We applied the CNN separation method to Monte Carlo gamma-ray and cosmic-ray events from the Crab Nebula in the energy range 10–100 TeV. The AUC values range from 0.753 to 0.879, and the significance of the observed gamma-ray excess is improved by 1.3 to 1.8 times compared with the case without the separation procedure.

Three different types of gamma-ray sensors (CsI, CLYC, and CeBr3) were flown on balloon flights as hosted payloads. Two CsI sensors were flown from a September 2014 flight from Fort Sumner, New Mexico for 18 h; CLYC and CeBr3 sensors were flown from Antarctica in December 2016 for 22 days. The data from these flights were used to test and characterize the operation of these sensors in a near-space environment. All sensors returned excellent data. Gamma rays, neutrons, and energetic galactic cosmic rays (GCRs) were measured. Expected atmospheric features, such as the Regener–Pfotzer maximum, were observed, and gamma-ray line emission from materials near the sensors, as well as atmospheric oxygen and nitrogen, were detected. The measured data were compared to simulations of energetic protons, neutrons, and 0.511 MeV gamma rays produced by GCR interactions with the atmosphere. While the simulated protons and neutrons generally matched the data, there were fewer simulated 0.511 MeV gamma rays than measured with the data. This mismatch is likely due to additional 0.511 MeV gamma rays produced in material near the sensors that were not taken into account in the simulations. Discussion is provided for how these types of measurements in space-based missions can be used to characterize upper atmospheric densities at planets with dense atmospheres like the Earth.

In recent papers the authors have followed up the common hypothesis that cosmic ray particles above about 1019 eV are mainly of extragalactic origin by evaluating the consequences for the associated cosmic gamma rays and neutrinos. Here they derive the expected fluxes of gamma rays and neutrinos, with energy of order 1020 eV, arising from the interaction of the protons with microwave background photons. The derived ratio of gamma rays to protons at 1020 eV is approximately 0.03, increasing to 0.15 at 5*1020 eV, for the situation where the cosmic ray sources are distributed nearly uniformly throughout the Universe; the corresponding figures for the neutrino to proton are factors of 5 and 24, respectively. An alternative, or perhaps additive, view is that most of the protons detected come via slow diffusion from the local Virgo cluster of galaxies. Here, significant fluxes of very energetic gamma rays and neutrinos are expected from the direction of this cluster. Specifically, for one steradian central on Virgo, the predicted gamma to proton ratio is approximately 0.2 at 1020 eV increasing to approximately 1.3 at 5*1020 eV. For the same energies, the neutrino to proton ratios are approximately=12 and 75.

We observed diffuse cosmic and atmospheric gamma rays at balloon altitudes with the Sub-MeV gamma-ray Imaging Loaded-on-balloon Experiment I (SMILE-I) as the first step toward a future all-sky survey with a high sensitivity. SMILE-I employed an electron-tracking Compton camera comprised of a gaseous electron tracker as a Compton-scattering target and a scintillation camera as an absorber. The balloon carrying the SMILE-I detector was launched from the Sanriku Balloon Center of the Institute of Space and Astronomical Science/Japan Space Exploration Agency on September 1, 2006, and the flight lasted for 6.8 hr, including level flight for 4.1 hr at an altitude of 32-35 km. During the level flight, we successfully detected 420 downward gamma rays between 100 keV and 1 MeV at zenith angles below 60 degrees. To obtain the flux of diffuse cosmic gamma rays, we first simulated their scattering in the atmosphere using Geant4, and for gamma rays detected at an atmospheric depth of 7.0 g cm-2, we found that 50% and 21% of the gamma rays at energies of 150 keV and 1 MeV, respectively, were scattered in the atmosphere prior to reaching the detector. Moreover, by using Geant4 simulations and the QinetiQ atmospheric radiation model, we estimated that the detected events consisted of diffuse cosmic and atmospheric gamma rays (79%), secondary photons produced in the instrument through the interaction between cosmic rays and materials surrounding the detector (19%), and other particles (2%). The obtained growth curve was comparable to Ling's model, and the fluxes of diffuse cosmic and atmospheric gamma rays were consistent with the results of previous experiments. The expected detection sensitivity of a future SMILE experiment measuring gamma rays between 150 keV and 20 MeV was estimated from our SMILE-I results and was found to be ten times better than that of other experiments at around 1 MeV.

The observation of very-high-energy (VHE, E > 100 GeV) gamma rays is an excellent tool to study the most energetic and violent environments in the Galaxy. This energy range is only accessible with ground-based instruments such as Imaging Atmospheric Cherenkov Telescopes (IACTs) that reconstruct the energy and direction of the primary gamma ray by observing the Cherenkov light from the induced extended air showers in Earths atmosphere. The main goals of Galactic VHE gamma-ray science are the identification of individual sources of cosmic rays (CRs), such as supernova remnants (SNRs), and the study of other extreme astrophysical objects at the highest energies, such as gamma-ray binaries and pulsar wind nebulae (PWNe).One of the main challenges is the discrimination between leptonic and hadronic gamma-ray production channels. To that end, the gamma-ray signal from each individual source needs to be brought into context with the multi-wavelength environment of the astrophysical object in question, particularly with observations tracing the density of the surrounding interstellar medium, or synchrotron radiation from relativistic electrons.In this review presented at the European Cosmic Ray Symposium 2014 (ECRS2014), the most recent developments in the field of Galactic VHE gamma-ray science are highlighted, with particular emphasis on SNRs and PWNe.

The Cygnus region of our Galaxy consists of an active star forming region and a wealth of various astrophysical sources such as pulsar wind nebulae (PWN), supernova remnants (SNRs), and massive star clusters. Massive stellar clusters and associations have been postulated as possible sources of cosmic rays (CRs) in our Galaxy. One example of a gamma-ray source associated with a stellar association lies in the Cygnus region known as the "Cygnus Cocoon". It is an extended region of gamma-ray emission in the Cygnus X region and attributed to a possible superbubble with freshly accelerated CRs which are hypothesized to produce gamma rays via interaction with the ambient gas nuclei. The emission region is an environment of lower particle density and is surrounded by ionization fronts like a carved-out cavity or a cocoon. CRs in the Cocoon could have originated in the OB2 association and been accelerated at the interaction sites of stellar winds of massive type O stars. So far, there is no clear association at TeV energies. In the study presented in this thesis, I used data collected by the HAWC Observatory over 1038 days to disentangle the TeV gamma-ray emission from 2HWC J2031+415, a source which was previously reported in the 2nd HAWC catalog and is collocated with the Cygnus Cocoon, into two components: a pulsar wind nebula and the Cygnus superbubble. The contribution from the Cygnus superbubble is detected at a significance level of ~ 12 sigma with maximum photon energies above 100 TeV, the highest measured yet. Based on the spectrum and morphology of gamma-rays across six decades of energy, and the non-detection of radio and X-ray photons from this region, the gamma-rays are plausibly of a hadronic origin. There is a spectral softening above 1 TeV, which can be explained by two hadronic scenarios. Either there is a leakage of CRs from the superbubble resulting in a spectral break from GeV to TeV or the spectral softening is due to cut-off energy, an upper limit to the particle acceleration by the stellar winds.

An excess of cosmic-ray positrons above 10 GeV with respect to the spallation reaction of cosmic rays (CRs) with the interstellar medium has been measured by Pamela, Fermi-LAT and with unprecedented precision by AMS-02. Various interpretations have been invoked to interpret this excess, such as production from supernova remnants, pulsar wind nebulae (PWNe) and dark matter. A dominant contribution from dark matter is ruled out by the bounds found in gamma rays and other indirect searches. Models where supernova remnants produce secondary CRs struggle to explain the observed CR fluxes by AMS-02. Finally, severe constraints for a significant PWN contribution come from the detection of very high-energy gamma-ray emission from Monogem and Geminga PWNe by Milagro and HAWC experiments. In this contribution we present a detailed study of the GeV gamma-ray halo around Geminga and Monogem, and show the constraints found for the contribution of these PWNe to the cosmic-ray positron excess, combining Milagro and HAWC data with measurements from the Fermi-LAT for the first time. We report the detection of a significant emission from Geminga PWN, derived by including the proper motion of its pulsar. We demonstrate that using gamma-ray data from the LAT is of central importance to provide a precise estimate for a PWN contribution to the cosmic positron flux

Energetic galactic cosmic ray (GCR) particles, arriving within the solar system, are modulated by the overall interplanetary field carried in the solar wind. Localized disturbances related to solar activity cause further reduction in intensity, the largest being Forbush decreases in which fluxes can fall ∼20% over a few days. Understanding Forbush decreases leads to a better understanding of the magnetic field structure related to shock waves and fast streams originating at the Sun since the propagation characteristics of the GCR probe much larger regions of space than do individual spacecraft instruments. We examined the temporal history of the integral GCR fluence (≥100 MeV) measured by the high‐sensitivity telescope (HIST) aboard the Polar spacecraft, along with the solar wind magnetic field and plasma data from the ACE spacecraft during a 40‐day period encompassing the 25 September 1998 Forbush decrease. We also examined the Forbush and (energetic storm particles) ESP event on 28 October 2003. It is the use of HIST in a high‐counting‐rate integral mode that allows previously poorly seen, short‐scale depressions in the GCR fluxes to be observed, adding crucial information on the origin of GCR modulation. Variability on time scales within the frequency range 0.001–1.0 mHz is detected. This paper concentrates on investigating four simple models for explaining short‐term reductions in the GCR intensity of both small and large amplitude. Specifically, these models are a local increase in magnetic scattering power, the passage of a shock discontinuity, and the passage of a tangential discontinuity or magnetic rope in the solar wind plasma. Analysis of the short‐scale GCR depressions during a test period in September through October 1998 shows that they are not correlated with changes in magnetic scattering power or fluctuations in solar wind speed or plasma density. However, magnetic field and plasma data during the test period of Forbush decrease strongly suggest the presence of an interplanetary coronal mass ejection (ICME). Use of a non‐force‐free magnetic rope model in conjunction with the energetic particle data allows modeling of the geometry of the ICME in terms of a magnetic cloud topology. It is only this cloud configuration that allows a satisfactory explanation of the magnitude of the Forbush event of 25 September 1998. Calculations made during the test period point to short‐scale GCR depressions being caused by either small‐scale magnetic flux rope structures or possibly tangential discontinuities in the solar wind.

Pulsars energise particles into lighthouse pencil beams and create extended relativistic outflows, pulsar wind nebulae (PWNe). In the very-high-energy (VHE) gamma-ray wave band, these PWNe represent to date the most populous class of Galactic sources. Nevertheless, the details of the energy conversion mechanisms in the vicinity of pulsars are not well understood, nor is it known which pulsars are able to drive PWNe and emit high-energy radiation. Due to its large field of view and unprecedented sensitivity, H.E.S.S. is the first instrument to allow for deep surveys of the Galactic plane in VHE gamma rays. This work presents the first ever systematic investigation of the connection of VHE gamma-ray sources and PWNe. Besides presenting two new candidate PWNe detected in this search, it is shown that pulsars with large spin-down energy flux are indeed with high probability associated with VHE gamma-ray sources, implying the existence of an efficient mechanism by which a large fraction of pulsar spin-down energy is converted into kinetic energy of particles. The results presented here make it very likely that future more sensitive VHE gamma-ray instruments will detect a rapidly increasing number of lower-luminosity PWNe.

We present observations of the Vela X region at 31 GHz using the Cosmic Background Imager (CBI). We find a strong compact radio source (59 × 41, FWHM) about the Vela pulsar, which we associate with the Vela pulsar wind nebula (PWN) recently discovered at lower radio frequencies. The CBI's 4' resolution for a 45' field of view allows the PWN to be studied in the large-scale context of Vela X. Filamentary structure in Vela X, which stands out in lower frequency maps, is very low level at 31 GHz. By combining the 10 CBI channels, which cover 26-36 GHz, and 8.4 GHz archive data, we study the spectral energy distribution (SED) of the PWN and the brightest filaments. Our results show that the spectral index α (Fν ∝ να) of the PWN is flat, or even marginally positive, with a value of α = 0.10 ± 0.06, while the Vela X filamentary structure has a negative spectral index of α = -0.28 ± 0.09. The SED inhomogeneity observed in Vela X suggests different excitation processes between the PWN and the filaments. We investigate whether the PWN's flat spectrum is a consequence of variability or truly reflects the SED of the object. We also investigate the nature of the Vela X filamentary structure. A faint filament crosses the PWN with its tangent sharing the same position angle as the PWN major axis, suggesting that it might be an extension of the PWN itself. The SED and bulk morphology of Vela X are similar to those of other well-studied plerions, suggesting that it might be powered by the pulsar. The peak of the PWN at 31 GHz is 80'' ± 20'' southwest of the peak at 8.4 GHz. This shift is confirmed by comparing the 31 GHz CBI image with higher resolution 5 GHz Australia Telescope Compact Array observations and is likely to be due to SED variations within the PWN.

We propose that the high energy Cosmic Ray particles up to the upturn commonly called the ankle, from around the spectral turn-down commonly called the knee, mostly come from Blue Supergiant star explosions. At the upturn, i.e., the ankle, Cosmic Rays probably switch to another source class, most likely extragalactic sources. To show this we recently compiled a set of Radio Supernova data where we compute the magnetic field, shock speed and shock radius. This list included both Blue and Red Supergiant star explosions; both data show the same magnetic field strength for these two classes of stars despite very different wind densities and velocities. Using particle acceleration theory at shocks, those numbers can be transformed into characteristic ankle and knee energies. Without adjusting any free parameters both of these observed energies are directly indicated by the supernova data. In the next step in the argument, we use the Supernova Remnant data of the starburst galaxy M82. We apply this analysis to Blue Supergiant star explosions: The shock will race to their outer edge with a magnetic field that is observed to follow over several orders of magnitude B ( r ) × r ∼ c o n s t . , with in fact the same magnetic field strength for such stellar explosions in our Galaxy, and other galaxies including M82. The speed is observed to be ∼0.1 c out to about 10 16 cm radius in the plasma wind. The Supernova shock can run through the entire magnetic plasma wind region at full speed all the way out to the wind-shell, which is of order parsec scale in M82. We compare and identify the Cosmic Ray spectrum in other galaxies, in the starburst galaxy M82 and in our Galaxy with each other; we suggest how Blue Supergiant star explosions can provide the Cosmic Ray particles across the knee and up to the ankle energy range. The data from the ISS-CREAM (Cosmic Ray Energetics and Mass Experiment at the International Space Station) mission will test this cosmic ray concept which is reasonably well grounded in two independent radio supernova data sets. The next step in developing our understanding will be to obtain future more accurate Cosmic Ray data near to the knee, and to use unstable isotopes of Cosmic Ray nuclei at high energy to probe the “piston” driving the explosion. We plan to incorporate these data with the physics of the budding black hole which is probably forming in each of these stars.

The supernova remnant (SNR) 0540–69.3, twin of the Crab Nebula, offers an excellent opportunity to study the continuum emission from a young pulsar and pulsar wind nebula (PWN). We present observations taken with the Very Large Telescope instruments MUSE and X-shooter in the wavelength range 3000–25000 Å, which allow us to study spatial variations of the optical spectra, along with the first near-infrared (NIR) spectrum of the source. We model the optical spectra with a power law (PL) F ν ∝ ν −α and find clear spatial variations (including a torus–jet structure) in the spectral index across the PWN. Generally, we find spectral hardening toward the outer parts, from α ∼ 1.1 to ∼0.1, which may indicate particle reacceleration by the PWN shock at the inner edge of the ejecta or alternatively time variability of the pulsar wind. The optical–NIR spectrum of the PWN is best described by a broken PL, confirming that several breaks are needed to model the full spectral energy distribution of the PWN, and suggesting the presence of more than one particle population. Finally, subtracting the PWN contribution from the pulsar spectrum we find that the spectrum is best described with a broken-PL model with a flat and a positive spectral index, in contrast to the Crab pulsar that has a negative spectral index and no break in the optical. This might imply that pulsar differences propagate to the PWN spectra.